Atomic spinors

An important advantage of ECP basis sets is their ability to incorporate approximately the physical effects of relativistic core contraction and associated changes in screening on valence orbitals, by suitable adjustments of the radius of the effective core potential. Thus, the ECP valence atomic orbitals can approximately mimic those of a fully relativistic (spinor) atomic calculation, rather than the non-relativistic all-electron orbitals they are nominally serving to replace. The partial inclusion of relativistic effects is an important physical correction for heavier atoms, particularly of the second transition series and beyond. Thus, an ECP-like treatment of heavy atoms is necessary in the non-relativistic framework of standard electronic-structure packages, even if the reduction in number of... 

The twin facts that heavy-atom compounds like BaF, T1F, and YbF contain many electrons and that the behavior of these electrons must be treated relati-vistically introduce severe impediments to theoretical treatments, that is, to the inclusion of sufficient electron correlation in this kind of molecule. Due to this computational complexity, calculations of P,T-odd interaction constants have been carried out with relativistic matching of nonrelativistic wavefunctions (approximate relativistic spinors) [42], relativistic effective core potentials (RECP) [43, 34], or at the all-electron Dirac-Fock (DF) level [35, 44]. For example, the first calculation of P,T-odd interactions in T1F was carried out in 1980 by Hinds and Sandars [42] using approximate relativistic wavefunctions generated from nonrelativistic single particle orbitals. 

Table 8 Second-order many-body perturbation theory corrections to beryllium-like ions using non-relativistic (E ), Dirac-Coulomb (E ) and Dirac-Coulomb-Breit (E ) hamiltonians, obtained using the atomic precursor to BERTHA, known as SWIRLES. Basis sets are even-tempered S-spinors of dimension N= 17, with exponent sets, Xi generated by Xi = abi-i, with a = 0.413, and p = 1.376. Angular momenta in the range 0 < / < 6 have been included in the partial wave expansion of each second-order energy, and the total relativistic correction toE has been collected as Ef. All energies in hartree.

Heavy atoms exhibit large relativistic effects, often too large to be treated perturba-tively. The Schrodinger equation must be supplanted by an appropriate relativistic wave equation such as Dirac-Coulomb or Dirac-Coulomb-Breit. Approximate one-electron solutions to these equations may be obtained by the self-consistent-field procedure. The resulting Dirac-Fock or Dirac-Fock-Breit functions are conceptually similar to the familiar Hartree-Fock functions the Hartree-Fock orbitals are replaced, however, by four-component spinors. Correlation is no less important in the relativistic regime than it is for the lighter elements, and may be included in a similar manner. 

The no-pair DCB Hamiltonian (6) is used as a starting point for variational or many-body relativistic calculations [9], The procedure is similar to the nonrelativistic case, with the Hartree-Fock orbitals replaced by the four-component Dirac-Fock-Breit (DFB) functions. The spherical symmetry of atoms leads to the separation of the one-electron equation into radial and spin-angular parts [10], The radial four-spinor has the so-called large component the upper two places and the small component Q, in the lower two. The quantum number k (with k =j+ 1/2) comes from the spin-angular equation, and n is the principal quantum number, which counts the solutions of the radial equation with the same k. Defining... 

Abstract. BERTHA is a 4-component relativistic molecular structure program based on relativistic Gaussian (G-spinor) basis sets which is intended to make affordable studies of atomic and molecular electronic structure, particularly of systems containing high-Z elements. This paper reviews some of the novel technical features embodied in the code, and assesses its current status, its potential and its prospects. 

Whilst this demonstrates that calculations using the methods of this paper may prove very useful in studies of molecules containing only low-Z atoms, a major objective has been to study systems containing heavier atoms. So far, only a limited number of molecular calculations have been carried out with BERTHA at the DHF level, mainly in connection with studies of hyperfine and PT-odd effects in heavy polar molecules such as YbF [33] and TIF [13]. The reader is referred to the literature for an assessment of these calculations and for technical details on the construction of basis sets which must not only describe molecular bonding properly but also give a good representation of spinors close to the heavy nuclei to handle the short-range electron-nuclear electroweak interactions. 

We demonstrated that by the selection of a representation of the Dirac Hamiltonian in the spinor space one may strongly influence the performance of the variational principle. In a vast majority of implementations the standard Pauli representation has been used. Consequently, computational algorithms developed in relativistic theory of many-electron systems have been constructed so that they are applicable in this representation only. The conditions, under which the results of these implementations are reliable, are very well understood and efficient numerical codes are available for both atomic and molecular calculations (see e.g. [16]). However, the representation of Weyl, if the external potential is non-spherical, or the representation of Biedenharn, in spherically-symmetric cases, seem to be attractive and, so far, hardly explored options. 

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. 

All the above restoration schemes are called nonvariational as compared to the variational one-center restoration (VOCR, see below) procedure proposed in [79, 80]. Proper behavior of the molecular orbitals (four-component spinors) in atomic cores of molecules can be restored in the scope of a variational procedure if the molecular pseudoorbitals (two-component pseudospinors) match correctly the original orbitals (large components of bispinors) in the valence region after the molecular RECP calculation. As is demonstrated in [69, 44], this condition is rather correct when the shape-consistent RECP is involved to the molecular calculation with explicitly... 

In the inner core region, the pseudospinors are smoothed, so that the electronic density with the pseudo-wave function is not correct. When operators describing properties of interest are heavily concentrated near or on nuclei, their mean values are strongly affected by the wave function in the inner region. The four-component molecular spinors must, therefore, be restored in the heavy-atom cores. 

Finally, the atomic two-component pseudospinors in the molecular basis are replaced by equivalent four-component spinors and the expansion coefficients from Eq. (6.3) are preserved ... 

The molecular four-component spinors constructed this way are orthogonal to the inner core spinors of the atom, because the atomic basis functions used in Eq. (6.4) are generated with the inner core shells treated as frozen. 

Variational one-center restoration. In the variational technique of one-center restoration (VOCR) [79, 80], the proper behavior of the four-component molecular spinors in the core regions of heavy atoms can be restored as an expansion in spherical harmonics inside the sphere with a restoration radius, Rvoa, that should not be smaller than the matching radius, Rc, used at the RECP generation. The outer parts of spinors are treated as frozen after the RECP calculation of a considered molecule. This method enables one to combine the advantages of two well-developed approaches, molecular RECP calculation in a gaussian basis set and atomic-type one-center calculation in numerical basis functions, in the most optimal way. This technique is considered theoretically in [80] and some results concerning the efficiency of the one-center reexpansion of orbitals on another atom can be found in [75]. 

The energy spectrum of atoms and ions with j j coupling can be found using the relativistic Hamiltonian of iV-electron atoms (2.1)-(2.7). Its irreducible tensorial form is presented in Chapter 19. The relativistic one-electron wave functions are four-component spinors (2.15). They are the eigenfunctions of the total angular momentum operator for the electron and are used to determine one-electron and two-electron matrix elements of relativistic interaction operators. These matrix elements, in the representation of occupation numbers, are the parameters that enter into the expansions of the operators corresponding to physical quantities (see general expressions (13.22) and (13.23)). 

All electrons, protons and neutrons, the elementary constituents of atoms, are fermions and therefore intrinsically endowed with an amount h/2 of angular momentum, known as spin. Like mass and charge, the other properties of fermions, the nature of spin is poorly understood. In quantum theory spin is treated purely mathematically in terms of operators and spinors, without physical connotation. 

Another four-component fully relativistic DFT program is the Beijing density functional (BDF) method [71]. Four-component numerical atomic spinors obtained by finite-difference atomic calculations are used for cores,... 

To construct the Dirac-Fock equations, it is assumed that the wave function for an atom having N electrons may be expressed as an antisymmetrized product of four-component Dirac spinors of the form shown in Eq. (9). For cases where a single antisymmetrized product is an eigenfunction of the total angular momentum operator J2, the JV-electron atomic wave function may be written... 

Equation (29) is an analog of Eq. (5) for the case of atomic basis functions represented as two-component spinors. 

---