Electron density relativistic effective core potentials

Slovenia), using the DFT implementation in the Gaussian03 code. Revision C.02 (8). The orbitals were described by a mixed basis set. A fully uncontracted basis set from LANL2DZ was used for the valence electrons of Re (9), augmented by two / functions Q =1.14 and 0.40) in the full optimization. Re core electrons were treated by the Hay-Wadt relativistic effective core potential (ECP) given by the standard LANL2 parameter set (electron-electron and nucleus-electron). The 6-3IG basis set was used to describe the rest of the system. The B3PW91 density functional was used in all calculations. 

Abstract. Investigation of P,T-parity nonconservation (PNC) phenomena is of fundamental importance for physics. Experiments to search for PNC effects have been performed on TIE and YbF molecules and are in progress for PbO and PbF molecules. For interpretation of molecular PNC experiments it is necessary to calculate those needed molecular properties which cannot be measured. In particular, electronic densities in heavy-atom cores are required for interpretation of the measured data in terms of the P,T-odd properties of elementary particles or P,T-odd interactions between them. Reliable calculations of the core properties (PNC effect, hyperfine structure etc., which are described by the operators heavily concentrated in atomic cores or on nuclei) usually require accurate accounting for both relativistic and correlation effects in heavy-atom systems. In this paper, some basic aspects of the experimental search for PNC effects in heavy-atom molecules and the computational methods used in their electronic structure calculations are discussed. The latter include the generalized relativistic effective core potential (GRECP) approach and the methods of nonvariational and variational one-center restoration of correct shapes of four-component spinors in atomic cores after a two-component GRECP calculation of a molecule. Their efficiency is illustrated with calculations of parameters of the effective P,T-odd spin-rotational Hamiltonians in the molecules PbF, HgF, YbF, BaF, TIF, and PbO. 

We review the Douglas-Kroll-Hess (DKH) approach to relativistic density functional calculations for molecular systems, also in comparison with other two-component approaches and four-component relativistic quantum chemistry methods. The scalar relativistic variant of the DKH method of solving the Dirac-Kohn-Sham problem is an efficient procedure for treating compounds of heavy elements including such complex systems as transition metal clusters, adsorption complexes, and solvated actinide compounds. This method allows routine ad-electron density functional calculations on heavy-element compounds and provides a reliable alternative to the popular approximate strategy based on relativistic effective core potentials. We discuss recent method development aimed at an efficient treatment of spin-orbit interaction in the DKH approach as well as calculations of g tensors. Comparison with results of four-component methods for small molecules reveals that, for many application problems, a two-component treatment of spin-orbit interaction can be competitive with these more precise procedures. 

The relativistic calculations on the electronic structure of actinide compounds were reviewed by Pyykko (1987). He also reviewed relativistic quantum chemistry in 1988, whereas the relativistic calculations were limited to small molecules containing one heavy atom only (Pyykko 1988). Calculations on the uranyl and neptunyl ions were introduced in the review article. The general information on the computational chemistry of heavy elements and relativistic calculation techniques appear in the book written by Balasubramanian (1997). There are several first-principle approaches to the electronic structure of actinide compounds. The relativistic effective core potential (ECP) and relativistic density functional methods are widely used for complex systems containing actinide elements. Pepper and Bursten (1991) reviewed relativistic quantum chemistry, while Schreckenbach et al. (1999) reviewed density functional calculations on actinide compounds in which theoretical background and application to actinide compounds were described in detail. The Encyclopedia of computational chemistry also contains examples including lanthanide and actinide elements (Schleyer et al. 1998). The various methods for the computational approach to the chemistry of transuranium elements are briefly described and summarized below. 

In the linear combination of atomic orbitals (LCAO) approximation of the molecular orbital, the energy of the overlap electron density between the atomic orbitals Xu and Xv due to attraction by the core. See Density Functional Theory (DFT), Hartree-Fock (HF), and the Self-consistent Field Relativistic Effective Core Potential Techniques for Molecules Containing Very Heavy Atoms Transition Metal Chemistry and Transition Metals Applications. 

In this chapter we tried to point out certain general observations and our personal preferences for choosing computational tools that represent a good compromise between accuracy and computational cost. When dealing with the heavy elements, relativistic corrections are imperative and, to a good approximation one can include them via small-core relativistic effective core potentials. All electron calculations with relativistic Hamiltonians are becoming more doable for large molecular systems, but with the exception of certain problems they do not necessarily afford an increase in accuracy worth the extra cost. Density functional theory has been proven to be a useful tool to probe the electronic structures and... 

The most common way of including relativistic effects in a calculation is by using relativisticly parameterized effective core potentials (RECP). These core potentials are included in the calculation as an additional term in the Hamiltonian. Core potentials must be used with the valence basis set that was created for use with that particular core potential. Core potentials are created by htting a potential function to the electron density distribution from an accurate relativistic calculation for the atom. A calculation using core potentials does not have any relativistic terms, but the effect of relativity on the core electrons is included. 

The electronic structure calculations were carried out using the hybrid density functional method B3LYP [15] as implemented in the GAUSSIAN-94 package [16], in conjunction with the Stevens-Basch-Krauss (SBK) [17] effective core potential (ECP) (a relativistic ECP for Zr atom) and the standard 4-31G, CEP-31 and (8s8p6d/4s4p3d) basis sets for the H, (C, P and N), and Zr atoms, respectively. 

With the exception of Ligand Field Theory, where the inclusion of atomic spin-orbit coupling is easy, a complete molecular treatment of relativity is difficult although not impossible. The work of Ellis within the Density Functional Theory DVXa framework is notable in this regard [132]. At a less rigorous level, it is relatively straightforward to develop a partial relativistic treatment. The most popular approach is to modify the potential of the core electrons to mimic the potential appropriate to the relativistically treated atom. This represents a specific use of so-called Effective Core Potentials (ECPs). Using ECPs reduces the numbers of electrons to be included explicitly in the calculation and hence reduces the execution time. Relativistic ECPs within the Hartree-Fock approximation [133] are available for all three transition series. A comparable frozen core approximation [134] scheme has been adopted for... 

In principle, it should also be possible to add a semi-loced potential to the non-relativistic all-electron Hamiltonian to eirrive at a quasi-rela-tivistic all-electron method. One such suggestion has been made by Delley [76], but the resulting method has only been tested for valence properties, which could also have been obtained by valence-only methods. Effective core potential methods have the advantage of a reduced computational effort (compared to all-electron methods) and are a valuable tool as long as one is aware of the limited domain of valence-only methods. Properties for which density variations in the atomic core are important should not be calculated this way. Examples are the electric field gradient at the nucleus or the nuclear magnetic shielding. 

Mitin A (2006) V., C. van Wullen Two-component relativistic density-functional calculations of the dimers of the halogens from bromine through element 117 using effective core potential and all-electron methods. J Chem Phys 124(6) 64305 Peng D, Liu W, Xiao Y, Cheng L (2007) Making four- and two-component relativistic density functional methods fully equivalent based on the idea of from atoms to molecule . J Chem Phys 127(10) 104106... 

---