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

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

Ab initio benchmark study for the oxidative addition of CH4 to Pd: importance of basis-set flexibility and polarization

To obtain a state-of-the-art benchmark potential energy surface (PES) for the archetypal oxidative addition of the methane C-H bond to the palladium atom, we have explored this PES using a hierarchical series of ab initio methods (Hartree-Fock, second-order Møller-Plesset perturbation theory, fourth-order Møller-Plesset perturbation theory with single, double and quadruple excitations, coupled cluster theory with single and double excitations (CCSD), and with triple excitations treated perturbatively [CCSD(T)]) and hybrid density functional theory using the B3LYP functional, in combination with a hierarchical series of ten Gaussian-type basis sets, up to g polarization. Relativistic effects are taken into account either through a relativistic effective core potential for palladium or through a full four-component all-electron approach. Counterpoise corrected relative energies of stationary points are converged to within 0.1-0.2 kcal/mol as a function of the basis-set size. Our best estimate of kinetic and thermodynamic parameters is -8.1 (-8.3) kcal/mol for the formation of the reactant complex, 5.8 (3.1) kcal/mol for the activation energy relative to the separate reactants, and 0.8 (-1.2) kcal/mol for the reaction energy (zero-point vibrational energy-corrected values in parentheses). This agrees well with available experimental data. Our work highlights the importance of sufficient higher angular momentum polarization functions, f and g, for correctly describing metal-d-electron correlation and, thus, for obtaining reliable relative energies. We show that standard basis sets, such as LANL2DZ+ 1f for palladium, are not sufficiently polarized for this purpose and lead to erroneous CCSD(T) results. B3LYP is associated with smaller basis set superposition errors and shows faster convergence with basis-set size but yields relative energies (in particular, a reaction barrier) that are ca. 3.5 kcal/mol higher than the corresponding CCSD(T) values

Ab initio study of NOx compounds adsorption on SnO2 surface

An ab initio study of the adsorption processes on NOx compounds on (1 1 0) SnO2 surface is presented with the aim of providing theoretical hints for the development of improved NOx gas sensors. From first principles calculations (DFT¿GGA approximation), the most relevant NO and NO2 adsorption processes are analyzed by means of the estimation of their adsorption energies. The resulting values and the developed model are also corroborated with experimental desorption temperatures for NO and NO2, allowing us to explain the temperature-programmed desorption experiments. The interference of the SO2 poisoning agent on the studied processes is discussed and the adsorption site blocking consequences on sensing response are analyzed.

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.

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.

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.

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.

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.