Magnetic coupling in the weac ferromagnet CuF2

CuF2 is known to be an antiferromagnetic compound with a weak ferromagnetism due to the anisotropy of its monoclinic unit cell (Dzialoshinsky-Moriya mechanism). We investigate the magnetic ordering of this compound by means of ab initio periodic unrestricted Hartree-Fock calculations and by cluster calculations which employ state-of-the-art configuration interaction expansions and modern density functional theory techniques. The combined use of periodic and cluster models permits us to firmly establish that the antiferromagnetic order arises from the coupling of one-dimensional subunits which themselves exhibit a very small ferromagnetic coupling between Cu neighbor cations. This magnetic order could be anticipated from the close correspondence between CuF2 and rutile crystal structures.

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.

Neutral atoms in ionic lattices: Excited states of KCl:Ag(0)

The optical-absorption spectrum of a cationic Ag0 atom in a KCl crystal has been studied theoretically by means of a series of cluster models of increasing size. Excitation energies have been determined by means of a multiconfigurational self-consistent field procedure followed by a second-order perturbation correlation treatment. Moreover results obtained within the density-functional framework are also reported. The calculations confirm the assignment of bands I and IV to transitions of the Ag-5s electron into delocalized states with mainly K-4s,4p character. Bands II and III have been assigned to internal transitions on the Ag atom, which correspond to the atomic Ag-4d to Ag-5s transition. We also determine the lowest charge transfer (CT) excitation energy and confirm the assignment of band VI to such a transition. The study of the variation of the CT excitation energy with the Ag-Cl distance R gives additional support to a large displacement of the Cl ions due to the presence of the Ag0 impurity. Moreover, from the present results, it is predicted that on passing to NaCl:Ag0 the CT onset would be out of the optical range while the 5s-5p transition would undergo a redshift of 0.3 eV. These conclusions, which underline the different character of involved orbitals, are consistent with experimental findings. The existence of a CT transition in the optical range for an atom inside an ionic host is explained by a simple model, which also accounts for the differences with the more common 3d systems. The present study sheds also some light on the R dependence of the s2-sp transitions due to s2 ions like Tl+.

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.

Adsorption energy and spin state of first-row transition metals adsorbed on MgO(100)

Slab and cluster model spin-polarized calculations have been carried out to study various properties of isolated first-row transition metal atoms adsorbed on the anionic sites of the regular MgO(100) surface. The calculated adsorption energies follow the trend of the metal cohesive energies, indicating that the changes in the metal-support and metal-metal interactions along the series are dominated by atomic properties. In all cases, except for Ni at the generalized gradient approximation level, the number of unpaired electron is maintained as in the isolated metal atom. The energy required to change the atomic state from high to low spin has been computed using the PW91 and B3LYP density-functional-theory-based methods. PW91 fails to predict the proper ground state of V and Ni, but the results for the isolated and adsorbed atom are consistent within the method. B3LYP properly predicts the ground state of all first-row transition atom the high- to low-spin transition considered is comparable to experiment. In all cases, the interaction with the surface results in a reduced high- to low-spin transition energy.

Adsoption-induced restructuring of gold nanochains

The chemical properties of single-atomic chains of gold atoms are investigated using density functional calculations. The nanochains are shown to be unusually chemically active with strong chemisorption of oxygen atoms and carbon monoxide. The chemisorption energies vary significantly with the strain/stress conditions for the chain. Oxygen atoms are found to energetically prefer to get incorporated into a chain forming a new type of gold-oxygen nanochain with a conductance of one quantum unit. We suggest that the long bond lengths observed in electron microscopy investigations of gold chains can be due to oxygen incorporation.

Ground- and excited-state properties of M-center oxygen vacancy aggregates in the bulk and the surface of MgO

Aggregates of oxygen vacancies (F centers) represent a particular form of point defects in ionic crystals. In this study we have considered the combination of two oxygen vacancies, the M center, in the bulk and on the surface of MgO by means of cluster model calculations. Both neutral and charged forms of the defect M and M+ have been taken into account. The ground state of the M center is characterized by the presence of two doubly occupied impurity levels in the gap of the material; in M+ centers the highest level is singly occupied. For the ground-state properties we used a gradient corrected density functional theory approach. The dipole-allowed singlet-to-singlet and doublet-to-doublet electronic transitions have been determined by means of explicitly correlated multireference second-order perturbation theory calculations. These have been compared with optical transitions determined with the time-dependent density functional theory formalism. The results show that bulk M and M+ centers give rise to intense absorptions at about 4.4 and 4.0 eV, respectively. Another less intense transition at 1.3 eV has also been found for the M+ center. On the surface the transitions occur at 1.6 eV (M+) and 2 eV (M). The results are compared with recently reported electron energy loss spectroscopy spectra on MgO thin films.

Electronic structure of a neutral oxygen vacancy in SrTiO3

The electronic structure of an isolated oxygen vacancy in SrTiO3 has been investigated with a variety of ab initio quantum mechanical approaches. In particular we compared pure density functional theory (DFT) approaches with the Hartree-Fock method, and with hybrid methods where the exchange term is treated in a mixed way. Both local cluster models and periodic calculations with large supercells containing up to 80 atoms have been performed. Both diamagnetic (singlet state) and paramagnetic (triplet state) solutions have been considered. We found that the formation of an O vacancy is accompanied by the transfer of two electrons to the 3d(z2) orbitals of the two Ti atoms along the Ti-Vac-Ti axis. The two electrons are spin coupled and the ground state is diamagnetic. New states associated with the defect center appear in the gap just below the conduction band edge. The formation energy computed with respect to an isolated oxygen atom in the triplet state is 9.4 eV.

Theoretical modeling of photon and electron-stimulated Na and K desorption from SiO2

The mechanism of generation of atomic Na and K from SiO2 samples has been studied using explicitly correlated wave function and density functional theory cluster calculations. Possible pathways for the photon and electron stimulated desorption of Na and K atoms from silicates are proposed, thus providing new insight on the generation of the tenuous Na and K atmosphere of the Moon.

First-principles study of the neutral molecular metal Ni(tmdt)2.

The electronic structure of the molecular solid Ni(tmdt)2, the only well characterized neutral molecular metal to date, has been studied by means of first-principles density functional calculations. It is shown that these calculations correctly describe the metallic vs semiconducting behavior of molecular conductors of this type. The origin of the band overlap leading to the metallic character and the associated Fermi surfaces has been studied.

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.

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.

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.