First principles LDA+U and GGA+U study of cerium oxides: Dependence on the effective U-parameter

The electronic structure and properties of cerium oxides (CeO2 and Ce2O3) have been studied in the framework of the LDA+U and GGA(PW91)+U implementations of density functional theory. The dependence of selected observables of these materials on the effective U parameter has been investigated in detail. The examined properties include lattice constants, bulk moduli, density of states, and formation energies of CeO2 and Ce2O3. For CeO2, the LDA+U results are in better agreement with experiment than the GGA+U results whereas for the computationally more demanding Ce2O3 both approaches give comparable accuracy. Furthermore, as expected, Ce2O3 is much more sensitive to the choice of the U value. Generally, the PW91 functional provides an optimal agreement with experiment at lower U energies than LDA does. In order to achieve a balanced description of both kinds of materials, and also of nonstoichiometric CeO2¿x phases, an appropriate choice of U is suggested for LDA+U and GGA+U schemes. Nevertheless, an optimum value appears to be property dependent, especially for Ce2O3. Optimum U values are found to be, in general, larger than values determined previously in a self-consistent way.

Theoretical approach to the magnetostructural correlations in spin-Peierls compound CuGeO3

A theoretical density-functional study has been carried out to analyze the exchange coupling in the chains of CuGeO3 using discrete models. The results show a good agreement with the experimental exchange coupling constant (J) together with a strong dependence of J with the Cu-O-Cu angle. The calculation of the J values for a distorted model indicates a larger degree of dimerization than those reported previously.

Controlling the spin of metal atoms adsorbed on oxide surfaces: Ni on regular and defective sites of the MgO(001) surface

We have analyzed the relative energy of nonmagnetic and magnetic low-lying electronic states of Ni atoms adsorbed on regular and defective sites of the MgO(001) surface. To this end cluster and periodic surface models are used within density functional theory. For Ni atoms adsorbed on oxygen vacancies at low coverage, the interaction energy between the metal and the support is much larger than on regular sites. Strong bonding results in a diamagnetic adsorbed species and the energy required to reach the high-spin state increases. Moreover, a correlation appears between the low-spin to high-spin energy difference and the interaction energy hypothesizing that it is possible to prepare the surface to tune the high-spin to low-spin energy difference. Magnetic properties of adsorbed thin films obtained upon increasing coverage are more difficult to interpret. This is because the metallic bond is readily formed and dominates over the effect of the atoms directly bound to the vacancy.

First principles analysis of the stability and diffusion of oxygen vacancies in meetal oxides

Oxygen vacancies in metal oxides are known to determine their chemistry and physics. The properties of neutral oxygen vacancies in metal oxides of increasing complexity (MgO, CaO, alpha-Al2O3, and ZnO) have been studied using density functional theory. Vacancy formation energies, vacancy-vacancy interaction, and the barriers for vacancy migration are determined and rationalized in terms of the ionicity, the Madelung potential, and lattice relaxation. It is found that the Madelung potential controls the oxygen vacancy properties of highly ionic oxides whereas a more complex picture arises for covalent ZnO.

MgO/Ag(001) interface structure and STM images from first principles

The interface of MgO/Ag(001) has been studied with density functional theory applied to slabs. We have found that regular MgO films show a small adhesion to the silver substrate, the binding can be increased in off-stoichiometric regimes, either by the presence of O vacancies at the oxide film or by a small excess of O atoms at the interface between the ceramic to the metal. By means of theoretical methods, the scanning tunneling microscopy signatures of these films is also analyzed in some detail. For defect free deposits containing 1 or 2 ML and at low voltages, tunnelling takes place from the surface Ag substrate, and at large positive voltages Mg atoms are imaged. If defects, oxygen vacancies, are present on the surface of the oxide they introduce much easier channels for tunnelling resulting in big protrusions and controlling the shape of the image, the extra O stored at the interface can also be detected for very thin films.

When Langmuir is too simple: H2 dissociation on Pd(111) at high coverage

Recent experiments of H2 adsorption on Pd(111) [T. Mitsui et al., Nature (London) 422, 705 (2003)] have questioned the classical Langmuir picture of second order adsorption kinetics at high surface coverage requiring pairs of empty sites for the dissociative chemisorption. Experiments find that at least three empty sites are needed. Through density functional theory, we find that H2 dissociation is favored on ensembles of sites that involve a Pd atom with no direct interaction with adsorbed hydrogen. Such active sites are formed by aggregation of at least 3 H-free sites revealing the complex structure of the "active sites."

Effective t-J model Hamiltonian parameters of monolayered cuprate superconductors from ab initio electronic structure calculations

The magnetic coupling constant of selected cuprate superconductor parent compounds has been determined by means of embedded cluster model and periodic calculations carried out at the same level of theory. The agreement between both approaches validates the cluster model. This model is subsequently employed in state-of-the-art configuration interaction calculations aimed to obtain accurate values of the magnetic coupling constant and hopping integral for a series of superconducting cuprates. Likewise, a systematic study of the performance of different ab initio explicitly correlated wave function methods and of several density functional approaches is presented. The accurate determination of the parameters of the t-J Hamiltonian has several consequences. First, it suggests that the appearance of high-Tc superconductivity in existing monolayered cuprates occurs with J/t in the 0.20¿0.35 regime. Second, J/t=0.20 is predicted to be the threshold for the existence of superconductivity and, third, a simple and accurate relationship between the critical temperatures at optimum doping and these parameters is found. However, this quantitative electronic structure versus Tc relationship is only found when both J and t are obtained at the most accurate level of theory.

Structure and energetics of mixed 4He-3He drops

Using a finite-range density functional, we have investigated the energetics and structural features of mixed helium clusters. The possibility of doping the cluster with a molecule of sulfur hexafluoride is also considered. It is seen that the repulsion introduced by the impurity strongly modifies the properties of the smallest drops. Although only a qualitative comparison is possible, the gross features displayed by our calculations are in agreement with recent experimental findings.

Multipole surface-plasmon modes on simple metals

The average multipole surface-plasmon energy for simple metals, as well as that of ordinary surface and bulk plasmons, is obtained using energy-weighted moments of the electronic response to sufficiently general external perturbations. A local approximation of exchange and correlation effects is used within a jellium model. Band-structure effects are incorporated through an effective electronic mass. Taking advantage of the transparency of the method, we analyze under what circumstances such modes might be observable. It is shown that due to an interplay between Coulomb and kinetic energies, the multipole modes become unobservable for increasing values of the transferred momentum (q) parallel to the surface. The value of q at which the multipole mode becomes unobservable is much smaller than the cutoff value for Landau damping. The effect of the electronic surface diffuseness is also analyzed. We compare our results with previous density-functional calculations and with recent experimental data for Na, K, and Cs.

Magnetic-field dependence of hole levels in self-assembled InGaAs quantum dots

Recent magnetotransport experiments of holes in InGaAs quantum dots [D. Reuter, P. Kailuweit, A. D. Wieck, U. Zeitler, O. Wibbelhoff, C. Meier, A. Lorke, and J. C. Maan, Phys. Rev. Lett. 94, 026808 (2005)] are interpreted by employing a multiband k¿p Hamiltonian, which considers the interaction between heavy hole and light hole subbands explicitly. No need of invoking an incomplete energy shell filling is required within this model. The crucial role we ascribe to the heavy hole-light hole interaction is further supported by one-band local-spin-density functional calculations, which show that Coulomb interactions do not induce any incomplete hole shell filling and therefore cannot account for the experimental magnetic field dispersion.

Probing vortices in 4He nanodroplets

We present static and dynamical properties of linear vortices in 4He droplets obtained from density functional calculations. By comparing the adsorption properties of different atomic impurities embedded in pure droplets and in droplets where a quantized vortex has been created, we suggest that Ca atoms should be the dopant of choice to detect vortices by means of spectroscopic experiments.

Ab initio cluster-model study of the on-top chemisorption of F and Cl on Si(111) and Ge(111) surfaces

The interaction of atomic F and Cl with Si4H9 and Ge4H9 cluster models has been studied by using ab initio pseudopotentials and basis sets of increasing complexity. The results show that the effect of d orbitals is important in order to reproduce the experimental findings. However, the use of polarization functions in the atoms which are directly involved in the chemisorption bond leads to results which are very close to those obtained using extended basis sets. The local nature of the chemisorption bond is also interpreted by means of a Mulliken population analysis. For F-Si4H9 and Cl-Si4H9 the present results are in good agreement with previous ab initio all-electron calculations, and for the chemisorption of Cl on Si(111) and Ge(111) surfaces, good agreement is found with respect to the available experimental results as well as with previous slab calculations based on the local-density-functional formalism.

Importance of the embedding environment on the strain within small rings in siliceous materials

The effect of the local environment on the energetic strain within small (SiO)N rings (with N=2,3) in silica materials is investigated via periodic model systems employing density functional calculations. Through comparison of the energies of various nonterminated systems containing small rings in strained and relatively unstrained environments, with alpha quartz, we demonstrate how small ring strain is affected by the nature of the embedding environment. We compare our findings with numerous previously reported calculations, often predicting significantly different small-ring strain energies, leading to a critical assessment of methods of calculating accurate localized ring energies. The results have relevance for estimates of the strain-induced response (e.g., chemical, photo, and radio) of small silica rings, and the propensity for them to form in bulk glasses, thin films, and nanoclusters.

Reduced ceria nanofilms from structure prediction

Experimentally, Ce2O3 films are used to study cerium oxide in its fully or partially reduced state, as present in many applications. We have explored the space of low energy Ce2O3 nanofilms using structure prediction and density functional calculations, yielding more than 30 distinct nanofilm structures. First, our results help to rationalize the roles of thermodynamics and kinetics in the preparation of reduced ceria nanofilms with different bulk crystalline structures (e.g. A-type or bixbyite) depending on the support used. Second, we predict a novel, as yet experimentally unresolved, nanofilm which has a structure that does not correspond to any previously reported bulk A2B3 phase and which has an energetic stability between that of A-type and bixbyite. To assist identification and fabrication of this new Ce2O3 nanofilm we calculate some observable properties and propose supports for its epitaxial growth.