First-principles periodic calculation of four-body spin terms in high-Tc cuprate superconductors

A general mapping between the energy of pertinent magnetic solutions and the diagonal terms of the spin Hamiltonian in a local representation provides the first general framework to extract accurate values for the many body terms of extended spin Hamiltonians from periodic first-principle calculations. Estimates of these terms for La2CuO4, the paradigm of high-Tc superconductor parent compounds, and for the SrCu2O3 ladder compound are reported. For La2CuO4, present results support experimental evidence by Toader et al. [Phys. Rev. Lett. 94, 197202 (2005)]. For SrCu2O3 even larger four-body spin amplitudes are found together with Jl/Jr=1 and non-negligible ferromagnetic interladder exchange.

First-principles periodic calculation of four-body spin terms in high-Tc Cuprate superconductors

A general mapping between the energy of pertinent magnetic solutions and the diagonal terms of the spin Hamiltonian in a local representation provides the first general framework to extract accurate values for the many body terms of extended spin Hamiltonians from periodic first-principle calculations. Estimates of these terms for La2CuO4, the paradigm of high-Tc superconductor parent compounds, and for the SrCu2O3 ladder compound are reported. For La2CuO4, present results support experimental evidence by Toader et al. [Phys. Rev. Lett. 94, 197202 (2005)]. For SrCu2O3 even larger four-body spin amplitudes are found together with Jl/Jr=1 and non-negligible ferromagnetic interladder exchange.

Towards and ab initio description of magnetism in ionic solids

The physical contributions to the KNiF3 magnetic exchange coupling integral have been obtained from specially designed ab initio cluster model calculations. Three important mechanisms have been identified. These are the delocalization of the magnetic orbitals into the anion p band, the variational contribution of the second-order interactions, and the many-body terms hidden in the two-body operator and the Heisenberg Hamiltonian.

X-ray photoelectron spectroscopy of oxygen adsorbates on Al(111): Theory experiment

We present the result of polar angle resolved x¿ray photoemission spectroscopy on Al(111)/O and cluster calculations of the O(1s) binding energy (BE) for various model situations. In the experimental data two O(1s) peaks are observed, separated by 1.3 eV. The angular behavior (depth¿resolution) could indicate that the lower BE peak is associated with an O atom under the surface, and the higher BE peak with an O atom above the surface. Equally, it could indicate oxygen islands on the surface where the perimeter atoms have a higher O(1s) BE than the interior atoms. The cluster calculations show that the former interpretation cannot be correct, since an O ads below the surface has a higher calculated O(1s) BE than one above. Cluster calculations simulating oxygen islands are, however, consistent with the experimental data.

Nonempirical cluster-model study of the chemisorption of atomic hydrogen on the (111) surface of diamondlike crystals

The interaction of atomic hydrogen with C4H9, Si4H9, and Ge4H9 model clusters has been studied using all-electron and pseudopotential ab initio Hartree-Fock computations with basis sets of increasing flexibility. The results show that the effect of polarization functions is important in order to reproduce the experimental findings, but their inclusion only for the atoms directly involved in the chemisorption bond is usually sufficient. For the systems H-C4H9 and H-Si4H9 all-electron and pseudopotential results are in excellent agreement when basis sets of comparable quality are used. Besides, semiempirical modified-neglect-of-differential-overlap computations provide quite reliable results both for diamond and silicon and have been used to investigate larger model clusters. The results confirm the local nature of chemisorption and further justify the use of minimal X4H9 model clusters.

Chemisorption of atomic aluminum on Si(111): Evidence for an adsorbate-induced relaxation based on ab initio cluster-model calculations

Interaction models of atomic Al with Si4H9, Si4H7, and Si6H9 clusters have been studied to simulate Al chemisorption on the Si(111) surface in the atop, fourfold atop, and open sites. Calculations were carried out using nonempirical pseudopotentials in the framework of the ab initio Hartree-Fock procedure. Equilibrium bond distances, binding energies for adsorption, and vibrational frequencies of the adatoms are calculated. Several basis sets were used in order to show the importance of polarization effects, especially in the binding energies. Final results show the importance of considering adatom-induced relaxation effects to specify the order of energy stabilities for the three different sites, the fourfold atop site being the preferred one, in agreement with experimental findings.

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.

Chemisorption of group-III metals on the Si(111) and Ge(111) surfaces: An ab initio study

Chemisorption of group-III metal adatoms on Si(111) and Ge(111) has been studied through the ab initio Hartree-Fock method including nonempirical pseudopotentials and using cluster models to simulate the surface. Three different high-symmetry sites (atop, eclipsed, and open) have been considered by using X4H9, X4H7, and X6H9 (X=Si,Ge) cluster models. In a first step, ideal surface geometries have been used. Metal-induced reconstruction upon chemisorption has also been taken into account. Equilibrium distances, binding energies, and vibrational frequencies have been obtained and compared with available experimental data. From binding-energy considerations, the atop and eclipsed sites seem to be the most favorable ones and thus a coadsorption picture may be suggested. Group-III metals exhibit a similar behavior and the same is true for Si(111) and Ge(111) surfaces when chemisorption is considered.

Evidence for oxygen-island formation on Al(111): Cluster-model theory and x-ray photoelectron spectroscopy

The O 1s x-ray photoelectron spectroscopy spectrum for Al(111)/O at 300 K shows two components whose behavior as a function of time and variation of detection angle are consistent with either (a) a surface species represented by the higher binding-energy (BE) component and a subsurface species represented by the lower BE component, or (b) small close-packed oxygen islands with the interior atoms represented by the lower BE component and the perimeter atoms by the higher BE component. We have modeled both situations using ab initio Hartree-Fock wave functions for clusters of Al and O atoms. For an O atom in a threefold site, it was found that a below-surface position gave a higher O 1s BE than an above-surface position, incompatible with interpretation (a). This change in the O 1s BE could arise because the bond for O to Al may have a more covalent character when the O is below the surface than when it is above the surface. We present evidence consistent with this view. An O adatom island with all the O atoms in threefold sites gives calculated O 1s BE's which are significantly higher for the perimeter O atoms. Further, the results for an isolated O island without the Al substrate present also give higher BE¿s for the perimeter atoms. Both these results are consistent with interpretation (b). Published scanning-tunneling-microscopy data supports the suggestion that the chemisorbed state consists of small, close-packed islands, whereas the presence of two vibrational modes in high-resolution electron-energy-loss spectroscopy data has been interpreted as representing surface and subsurface oxygen atoms. In light of the present results, we suggest that a vibrational interpretation in terms of interior and perimeter adatoms should be considered.

Evidence for two different bonding mechanisms of Al on Si(111)

By means of the ab initio cluster-model approach, we present theoretical evidence for two different mechanisms of bonding of atomic Al to Si(111). On the atop site (T1) the interaction of atomic Al to Si(111) is characteristic of an ionic bond whereas interaction above the threefold eclipsed site (T4) leads to the formation of a typical covalent bond. Moreover, both sites have a similar interaction energy if electronic correlation effects are included. While the conclusions regarding the nature of the chemisorption bond in the two sites do not depend either on the cluster-model size, the kind of embedding hydrogen atoms used, or the quality of the wave function (Hartree-Fock or configuration interaction), the chemisorption energy depends strongly on the wave function used. In fact, inclusion of correlation energy is necessary to properly describe the interaction energies.

Ab initio valence-bond cluster model for ionic solids: Alkaline-earth oxides

A linear M-O-M (M=metal, O=oxygen) cluster embedded in a Madelung field, and also including the quantum effects of the neighboring ions, is used to represent the alkaline-earth oxides. For this model an ab initio wave function is constructed as a linear combination of Slater determinants written in an atomic orbital basis set, i.e., a valence-bond wave function. Each valence-bond determinant (or group of determinants) corresponds to a resonating valence-bond structure. We have obtained ab initio valence-bond cluster-model wave functions for the electronic ground state and the excited states involved in the optical-gap transitions. Numerical results are reasonably close to the experimental values. Moreover, the model contains the ionic model as a limiting case and can be readily extended and improved.

Measures of ionicity of alkaline-earth oxides from the analysis of ab initio cluster wave functions

We present an analysis of the M-O chemical bonding in the binary oxides MgO, CaO, SrO, BaO, and Al2O3 based on ab initio wave functions. The model used to represent the local environment of a metal cation in the bulk oxide is an MO6 cluster which also includes the effect of the lattice Madelung potential. The analysis of the wave functions for these clusters leads to the conclusion that all the alkaline-earth oxides must be regarded as highly ionic oxides; however, the ionic character of the oxides decreases as one goes from MgO, almost perfectly ionic, to BaO. In Al2O3 the ionic character is further reduced; however, even in this case, the departure from the ideal, fully ionic, model of Al3+ is not exceptionally large. These conclusions are based on three measures, a decomposition of the Mq+-Oq- interaction energy, the number of electrons associated to the oxygen ions as obtained from a projection operator technique, and the analysis of the cation core-level binding energies. The increasing covalent character along the series MgO, CaO, SrO, and BaO is discussed in view of the existing theoretical models and experimental data.

Ground and excited states of KNiF3: An ab initio cluster-model approach

Finite cluster models and a variety of ab initio wave functions have been used to study the electronic structure of bulk KNiF3. Several electronic states, including the ground state and some charge-transfer excited states, have been considered. The study of the cluster-model wave functions has permitted an understanding of the nature of the chemical bond in the electronic ground state. This is found to be highly ionic and the different ionic and covalent contributions to the bonding have been identified and quantified. Finally, we have studied the charge-transfer excited states leading to the optical gap and have found that calculated and experimental values are in good agreement. The wave functions corresponding to these excited states have also been analyzed and show that although KNiF3 may be described as a ligand-to-metal charge-transfer insulator there is a strong configuration mixing with the metal-to-metal charge-transfer states.

Extent and limitations of density functional theory in describing magnetic systems

The performance of density-functional theory to solve the exact, nonrelativistic, many-electron problem for magnetic systems has been explored in a new implementation imposing space and spin symmetry constraints, as in ab initio wave function theory. Calculations on selected systems representative of organic diradicals, molecular magnets and antiferromagnetic solids carried out with and without these constraints lead to contradictory results, which provide numerical illustration on this usually obviated problem. It is concluded that the present exchange-correlation functionals provide reasonable numerical results although for the wrong physical reasons, thus evidencing the need for continued search for more accurate expressions.