Performance of correlation functionals in ab initio chemisorption cluster-model calculations: Alkali metals on Si(111)

The performance of different correlation functionals has been tested for alkali metals, Li to Cs, interacting with cluster models simulating different active sites of the Si(111) surface. In all cases, the ab initio Hartree-Fock density has been obtained and used as a starting point. The electronic correlation energy is then introduced as an a posteriori correction to the Hartree-Fock energy using different correlation functionals. By making use of the ionic nature of the interaction and of different dissociation limits we have been able to prove that all functionals tested introduce the right correlation energy, although to a different extent. Hence, correlation functionals appear as an effective and easy way to introduce electronic correlation in the ab initio Hartree-Fock description of the chemisorption bond in complex systems where conventional configuration interaction techniques cannot be used. However, the calculated energies may differ by some tens of eV. Therefore, these methods can be employed to get a qualitative idea of how important correlation effects are, but they have some limitations if accurate binding energies are to be obtained.

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.

Ab initio theoretical comparative study of magnetic coupling in KNiF3 and K2NiF4s

The origin of magnetic coupling in KNiF3 and K2 NiF4 is studied by means of an ab initio cluster model approach. By a detailed study of the mapping between eigenstates of the exact nonrelativistic and spin model Hamiltonians it is possible to obtain the magnetic coupling constant J and to compare ab initio cluster-model values with those resulting from ab initio periodic Hartree-Fock calculations. This comparison shows that J is strongly determined by two-body interactions; this is a surprising and unexpected result. The importance of the ligands surrounding the basic metal-ligand-metal interacting unit is reexamined by using two different partitions and the constrained space orbital variation method of analysis. This decomposition enables us to show that this effect is basically environmental. Finally, dynamical electronic correlation effects have found to be critical in determining the final value of the magnetic coupling constant.

Local character of magnetic coupling in ionic solids

Magnetic interactions in ionic solids are studied using parameter-free methods designed to provide accurate energy differences associated with quantum states defining the Heisenberg constant J. For a series of ionic solids including KNiF3, K2NiF4, KCuF3, K2CuF4, and high- Tc parent compound La2CuO4, the J experimental value is quantitatively reproduced. This result has fundamental implications because J values have been calculated from a finite cluster model whereas experiments refer to infinite solids. The present study permits us to firmly establish that in these wide-gap insulators, J is determined from strongly local electronic interactions involving two magnetic centers only thus providing an ab initio support to commonly used model Hamiltonians.

Ab Initio study of magnetic interactions in the KCuF3 and K2CuF4 low-dimensional Systems

The ab initio cluster model approach has been used to study the electronic structure and magnetic coupling of KCuF3 and K2CuF4 in their various ordered polytype crystal forms. Due to a cooperative Jahn-Teller distortion these systems exhibit strong anisotropies. In particular, the magnetic properties strongly differ from those of isomorphic compounds. Hence, KCuF3 is a quasi-one-dimensional (1D) nearest neighbor Heisenberg antiferromagnet whereas K2CuF4 is the only ferromagnet among the K2MF4 series of compounds (M=Mn, Fe, Co, Ni, and Cu) behaving all as quasi-2D nearest neighbor Heisenberg systems. Different ab initio techniques are used to explore the magnetic coupling in these systems. All methods, including unrestricted Hartree-Fock, are able to explain the magnetic ordering. However, quantitative agreement with experiment is reached only when using a state-of-the-art configuration interaction approach. Finally, an analysis of the dependence of the magnetic coupling constant with respect to distortion parameters is presented.

Electronic structure and magnetic interactions of the spin-chain compounds Ca2CuO3 and Sr2CuO3

The results are presented of a combined periodic and cluster model approach to the electronic structure and magnetic interactions in the spin-chain compounds Ca2CuO3 and Sr2CuO3. An extended t-J model is presented that includes in-chain and interchain hopping and magnetic interaction processes with parameters extracted from ab initio calculations. For both compounds, the in-chain magnetic interaction is found to be around -240 meV, larger than in any of the other cuprates reported in the literature. The interchain magnetic coupling is found to be weakly antiferromagnetic, -1 meV. The effective in-chain hopping parameters are estimated to be ~650 meV for both compounds, whereas the value of the interchain hopping parameter is 30 meV for Sr2CuO3 and 40 meV for Ca2CuO3, in line with the larger interchain distance in the former compound. These effective parameters are shown to be consistent with expressions recently suggested for the Néel temperature and the magnetic moments, and with relations that emerge from the t-J model Hamiltonian. Next, we investigate the physical nature of the band gap. Periodic calculations indicate that an interpretation in terms of a charge-transfer insulator is the most appropriate one, in contrast to the suggestion of a covalent correlated insulator recently reported in the literature.

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.

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.

Ionic-covalent transition in titanium oxides

The nature of the chemical bond in three titanium oxides of different crystal structure and different formal oxidation state has been studied by means of the ab initio cluster-model approach. The covalent and ionic contributions to the bond have been measured from different theoretical techniques. All the analysis is consistent with an increasing of covalence in the TiO, Ti2O3, and TiO2 series as expected from chemical intuition. Moreover, the use of the ab initio cluster-model approach combined with different theoretical techniques has permitted us to quantify the degree of ionic character, showing that while TiO can approximately be described as an ionic compound, TiO2 is better viewed as a rather covalent oxide.

Deformed-jellium model for the fission of multiply charged metal clusters

A deformed-jellium model is used to calculate the fission barrier height of positive doubly charged sodium clusters within an extended Thomas-Fermi approximation. The fissioning cluster is continuously deformed from the parent configuration until it splits into two fragments. Although the shape of the fission barrier obviously depends on the parametrization of the fission path, we have found that remarkably, the maximum of the barrier corresponds to a configuration in which the emerging fragments are already formed and rather well apart. The implication of this finding in the calculation of critical numbers for fission is illustrated in the case of multiply charged Na clusters.