Dirac spherical

The DIRAC package [55,56], devised by Saue and collaborators, rather than exploiting the group theoretical properties of Dirac spherical 4-spinors as in BERTHA, treats each component in terms of a conventional quantum chemical basis of real-valued Cartesian functions. The approach used in DIRAC, building on earlier work by Rosch [80] for semi-empirical models, uses a quaternion matrix representation of one electron operators in a basis of Kramers pairs. The transformation properties of these matrices, analysed in [55], are used to build point group transformation properties into the Fock matrix. 

Campbell s Theorem, 174 Cartwright, M. L., 388 Caywood, T. E., 313 C-coefficients, 404 formulas for, 406 recursion relations, 406 relation to spherical harmonics, 408 tabulations of, 408 Wigner s formula, 408 Central field Dirac equation in, 629 Central force law... 

For a Dirac particle in a central field 4> is spherically symmetrical and A = 0. Setting the potential energy V(r) — e(r), the Dirac Hamiltonian becomes... 

Despite the complication due to the interdependence of orbital and spin angular momenta, the Dirac equation for a central field can be separated in spherical polar coordinates [63]. The energy eigenvalues for the hydrogen atom (V(r) = e2/r, in electrostatic units), are equivalent to the relativistic terms of the old quantum theory [64]... 

There are two approaches to map crystal charge density from the measured structure factors by inverse Fourier transform or by the multipole method [32]. Direct Fourier transform of experimental structure factors was not useful due to the missing reflections in the collected data set, so a multipole refinement is a better approach to map charge density from the measured structure factors. In the multipole method, the crystal charge density is expanded as a sum of non-spherical pseudo-atomic densities. These consist of a spherical-atom (or ion) charge density obtained from multi-configuration Dirac-Fock (MCDF) calculations [33] with variable orbital occupation factors to allow for charge transfer, and a small non-spherical part in which local symmetry-adapted spherical harmonic functions were used. 

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... 

Until now we have described electrons states as obeying the wave equation without specifying which wave equation. In general it is necessary to solve the Dirac equation even for a local potential. If the local potential is spherically symmetric the Dirac equation reduces to ... 

Apart of historical reasons, there are several features of the Dirac-Pauli representation which make its choice rather natural. In particular, it is the only representation in which, in a spherically-symmetric case, large and small components of the wavefunction are eigenfunctions of the orbital angular momentum operator. However, this advantage of the Dirac-Pauli representation is irrelevant if we study non-spherical systems. It appears that the representation of Weyl has several very interesting properties which make attractive its use in variational calculations. Also several other representations seem to be worth of attention. Usefulness of these ideas is illustrated by an example. 

The Dirac-Pauli representation is most commonly used in all applications of the Dirac theory to studies on electronic structure of atoms and molecules. Apart of historical reasons, there are several features of this representation which make its choice quite natural. Probably the most important is a well defined symmetry of and in the case of spherically-symmetric potentials V. The Dirac Hamiltonian... 

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. 

Here pj7m) denotes the spherical-wave free-electron function with the usual notations for Dirac angular quantum numbers. The numbers jlm are fixed by the overlap with the bound-electron wave function a) = njlm) where n is the principal quantum number. Integration over p is interpreted as integration over energies Ep = /p2 + m2. 

Band Structure for the Free-Electron Case. If the electron is free, then the Bloch functions are simple plane waves, because the wavef unctions nk(r) used for the expansion Eq. (8.4.2) are themselves plane waves. For an electron gas with no lattice and no imposed symmetry, Fermi-Dirac statistics apply At 0 K all electrons pair up (spin-up and spin-down), with an occupancy of 2 for every k value from k 0 to the Fermi wavevector kF =1.92/rs = 3.63 a0/rs, and from zero energy up to the Fermi energy eF = h2kF2 /rn 50.1 eV rs/a0) 2, where rs is the radius per conduction electron and a0 is the Bohr radius, and the energy levels are spherically symmetric in k-space. The Fermi surface is a sphere of radius kF. Note that the ratio (rs/a0) varies from 0.2 to 1.0 nm for metals (Table 8.3). 

In these expressions the index i runs over electrons and a runs over nuclei. The Fermi contact term describes the magnetic interaction between the electron spin and nuclear spin magnetic moments when there is electron spin density at the nucleus. This condition is imposed by the presence of the Dirac delta function S(rai) in the expression. The dipole-dipole coupling term describes the classical interaction between the magnetic dipole moments associated with the electron and nuclear spins. It depends on the relative orientations of the two moments described in equation (7.145) and falls off as the inverse cube of the separations of the two dipoles. The cartesian form of the dipole-dipole interaction to some extent masks the simplicity of this term. Using the results of spherical tensor algebra from the previous chapter, we can bring this into the open as... 

Attempts to formulate a causal description of electron spin have not been completely successful. Two approaches were to model the motion on either a rigid sphere with the Pauli equation [102] as basis, or a point particle using Dirac s equation, which is pursued here no further. The methodology is nevertheless of interest and consistent with the spherical rotation model. The basic problem is to formulate a wave function in polar form E = RetS h as a spinor, by expressing each complex component in spinor form... 

The sampling of the points is usually taken as a sum of points centered around the atoms and given as the product of a certain adjustable constant sampling weight t and spherical Fermi-Dirac distribution centered at each atomic site and given by... 

Experimental studies on the K0/Ka x-ray intensity ratio for 3d elements have shown [18-23] that this ratio changes under influence of the chemical environment of the 3d atom. Brunner et al. [22] explained their experimental results due to the change in screening of 3p electrons by 3d valence electrons as well as the polarization effect. Band et al. [24] used the scattered-wave (SW) Xa MO method [25] and calculated the chemical effect on the K0/Ka ratios for 3d elements. They performed the SW-Xa MO calculations for different chemical compounds of Cr and Mn. The spherically averaged self-consistent-field (SCF) potential and the total charge of valence electrons in the central atom, obtained by the MO calculations, were used to solve the Dirac equation for the central atom and the x-ray transition probabilities were calculated. 

In addition, Dirac s theory provides a direct explanation for the fact that the electron magnetic dipole moment is about twice the value expected classically on the basis of a spherical charged particle rotating around one... 

The potential surrouding each atom in a molecule is not the same as that for the free atom, because electron transfer occurs between atoms in the molecule. This means that atomic orbitals in the molecule are distinct from those in the free atom. Accordingly, it is necessary to use atomic orbitals optimized for each atomic potential in the molecule, as basis functions. In the present methods, the molecular wave functions were expressed as linear combinations of atomic orbitals obtained by numerically solving the Dirac-Slater or Hartree-Fock-Slater equations in the atomic-like potential derived from the spherical average of the molecular charge density around the nuclei [15]. Thus the atomic orbitals used as basis functions were automatically optimized for the molecule and thus the minimum size of the present basis set has enough flexibility to form accurate molecular orbitals. 

---