Slabs of stabilized jellium: Quantum-size and self-compression effects

We examine thin films of two simple metals (aluminum and lithium) in the stabilized jellium model, a modification of the regular jellium model in which a constant potential is added inside the metal to stabilize the system for a given background density. We investigate quantum-size effects on the surface energy and the work function. For a given film thickness we also evaluate the density yielding energy stability, which is found to be slightly higher than the equilibrium density of the bulk system and to approach this value in the limit of thick slabs. A comparison of our self-consistent calculations with the predictions of the liquid-drop model shows the validity of this model.Comment: 7 pages, 6 figures, to appear in Phys. Rev.

Verification of the mixed layer depth in the OceanMAPS operational forecast model for Austral autumn

The ocean mixed layer depth is an important parameter describing the exchange of fluxes between the atmosphere and ocean. In ocean modelling a key factor in the accurate representation of the mixed layer is the parameterization of vertical mixing. An ideal opportunity to investigate the impact of different mixing schemes was provided when the Australian Bureau of Meteorology upgraded its operational ocean forecasting model, OceanMAPS to version 3.0. In terms of the mixed layer, the main difference between the old and new model versions was a change of vertical mixing scheme from that of Chen et al. (1994) to the General Ocean Turbulence Model.The model estimates of the mixed layer depth were compared with those derived from Argo observations. Both versions of the model exhibited a deep bias in tropical latitudes and a shallow bias in the Southern Ocean, consistent with previous studies. The bias, however, was greatly reduced in version 3.0, and variance between model runs decreased. Additionally, model skill against climatology also improved significantly. Further analysis discounted changes to model resolution outside of the Australian region having a significant impact on these results, leaving the change in vertical\ud mixing scheme as the main factor in the assessed improvements to mixed layer depth representation.</p

Lattice Dynamics and the High Pressure Equation of State of Au

Elastic constants and zone-boundary phonon frequencies of gold are calculated by total energy electronic structure methods to twofold compression. A generalized force constant model is used to interpolate throughout the Brillouin zone and evaluate moments of the phonon distribution. The moments are used to calculate the volume dependence of the Gruneisen parameter in the fcc solid. Using these results with ultrasonic and shock data, we formulate the complete free energy for solid Au. This free energy is given as a set of closed form expressions, which are valid to compressions of at least V/V_0 = 0.65 and temperatures up to melting. Beyond this density, the Hugoniot enters the solid-liquid mixed phase region. Effects of shock melting on the Hugoniot are discussed within an approximate model. We compare with proposed standards for the equation of state to pressures of ~200 GPa. Our result for the room temperature isotherm is in very good agreement with an earlier standard of Heinz and Jeanloz.Comment: 13 pages, 8 figures. Accepted by Phys. Rev.

Block-Diagonalization and f-electron Effects in Tight-Binding Theory

We extend a tight-binding total energy method to include f-electrons, and apply it to the study of the structural and elastic properties of a range of elements from Be to U. We find that the tight-binding parameters are as accurate and transferable for f-electron systems as they are for d-electron systems. In both cases we have found it essential to take great care in constraining the fitting procedure by using a block-diagonalization procedure, which we describe in detail.Comment: 9 pages, 6 figure

Simple DFT-LSDA modeling of the molecular-like aspects of ultra-thin film properties

Ordered ultra-thin films (UTF`s) are atomic n-layers (n = 1,2,3,...) with translational symmetry in-plane and molecular-like inter-planar spacings. Though commonly used (especially at relatively large n-values) as models of crystalline surfaces, they are intrinsically interesting and of growing technological significance as the basic building blocks of multi-layer electronic devices. Predicting the structure and properties of even a simple diatomic 1-layer means addressing aspects of molecular binding (and boundary conditions) in the context of an extended, periodically bounded system. At the level of refinement provided by the local spin density approximation to Density Functional Theory, the baseline standard of today`s predictive, chemically specific solid-state calculations, a number of technical and fundamental issues arise. The authors focus on treatment of the isolated atoms, on basis sets, and on numerical precision, as illustrated by the Fe atom and BN 1- and 2-layer calculations. Computational requirements are illustrated by a brief summary of recently completed calculations on crystalline sapphire, {alpha}-Al{sub 2}O{sub 3}, which used the same code