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.

Density-functional calculations of magnetoplasmons in quantum rings

We have studied the structure and dipole charge-density response of nanorings as a function of the magnetic field using local-spin-density-functional theory. Two small rings consisting of 12 and 22 electrons confined by a positively charged background are used to represent the cases of narrow and wide rings. The results are qualitatively compared with experimental data existing on microrings and on antidots. A smaller ring containing five electrons is also analyzed to allow for a closer comparison with a recent experiment on a two-electron quantum ring.

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.

Magnetic structure of Li2CuO2: From ab initio calculations to macroscopic simulations

The magnetic structure of the edge-sharing cuprate compound Li2CuO2 has been investigated with highly correlated ab initio electronic structure calculations. The first- and second-neighbor in-chain magnetic interactions are calculated to be 142 and -22 K, respectively. The ratio between the two parameters is smaller than suggested previously in the literature. The interchain interactions are antiferromagnetic in nature and of the order of a few K only. Monte Carlo simulations using the ab initio parameters to define the spin model Hamiltonian result in a Nel temperature in good agreement with experiment. Spin population analysis situates the magnetic moment on the copper and oxygen ions between the completely localized picture derived from experiment and the more delocalized picture based on local-density calculations.

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.

Magnetic structure of Li2CuO2: From ab initio calculations to macroscopic simulations

The magnetic structure of the edge-sharing cuprate compound Li2CuO2 has been investigated with highly correlated ab initio electronic structure calculations. The first- and second-neighbor in-chain magnetic interactions are calculated to be 142 and -22 K, respectively. The ratio between the two parameters is smaller than suggested previously in the literature. The interchain interactions are antiferromagnetic in nature and of the order of a few K only. Monte Carlo simulations using the ab initio parameters to define the spin model Hamiltonian result in a Nel temperature in good agreement with experiment. Spin population analysis situates the magnetic moment on the copper and oxygen ions between the completely localized picture derived from experiment and the more delocalized picture based on local-density calculations.

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.

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.

A comparative analysis by means of quantum molecular similarity measures of density distributions derived from conventional ab initio and density functional methods

A procedure based on quantum molecular similarity measures (QMSM) has been used to compare electron densities obtained from conventional ab initio and density functional methodologies at their respective optimized geometries. This method has been applied to a series of small molecules which have experimentally known properties and molecular bonds of diverse degrees of ionicity and covalency. Results show that in most cases the electron densities obtained from density functional methodologies are of a similar quality than post-Hartree-Fock generalized densities. For molecules where Hartree-Fock methodology yields erroneous results, the density functional methodology is shown to yield usually more accurate densities than those provided by the second order Møller-Plesset perturbation theory

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.

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.

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.

Ab initio absorption spectrum of NiO combining molecular dynamics with the embedded cluster approach in a discrete reaction field

We developed a procedure that combines three complementary computational methodologies to improve the theoretical description of the electronic structure of nickel oxide. The starting point is a Car-Parrinello molecular dynamics simulation to incorporate vibrorotational degrees of freedom into the material model. By means ofcomplete active space self-consistent field second-order perturbation theory (CASPT2) calculations on embedded clusters extracted from the resulting trajectory, we describe localized spectroscopic phenomena on NiO with an efficient treatment of electron correlation. The inclusion of thermal motion into the theoretical description allowsus to study electronic transitions that, otherwise, would be dipole forbidden in the ideal structure and results in a natural reproduction of the band broadening. Moreover, we improved the embedded cluster model by incorporating self-consistently at the complete active space self-consistent field (CASSCF) level a discrete (or direct) reaction field (DRF) in the cluster surroundings. The DRF approach offers an efficient treatment ofelectric response effects of the crystalline embedding to the electronic transitions localized in the cluster. We offer accurate theoretical estimates of the absorption spectrum and the density of states around the Fermi level of NiO, and a comprehensive explanation of the source of the broadening and the relaxation of the charge transferstates due to the adaptation of the environment