The ZORA equation is the zeroth order regular approximation to the Dirac equation. The relativistic and non-relativistic Kohn-Sham DFT equations can be written as:
[ T + VKS] φi = εi φi
The potential VKS is the local Kohn-Sham potential, see also the Kohn-Sham MO model. T is the kinetic energy operator, which is different for each relativistic method:
TDirac = σ⋅p c2/(2c2 + εi – VKS) σ⋅p
TZORA = σ⋅p c2/(2c2 – VKS) σ⋅p
TNR = p2/2
TDirac is the energy dependent Dirac kinetic energy operator, TZORA is the ZORA kinetic energy operator, and TNR is the non-relativistic kinetic energy operator.
The scaled ZORA spinor energy εiscaled is:
εiscaled = εiZORA / < φi | σ⋅p c2/(2c2 – VKS)2 σ⋅p | φi >
The scaled ZORA spinor energies are exactly equal to the Dirac energies for hydrogen like systems,
where VKS = -Z/r.
A magnetic field can be included by minimal substitution for a negatively charged particle:
p → Π = p + A/c
The ZORA kinetic energy can be split in a scalar relativistic part TSR and a spin-orbit coupling TSO term:
TZORA = TSR + TSO = p⋅ c2/(2c2 – VKS) p + c2/(2c2 – VKS)2 σ⋅ (∇ VKS × p )
Spin-orbit coupling results from the interaction of the electron magnetic moments with the magnetic field generated by their own orbital motion. The computation is more involved than for scalar relativistic or non-relativistic calculations because spin-orbit coupling does not have the ordinary molecular symmetry. Instead, it is invariant under the molecular point double group, which affects both spin and spatial coordinates. One of its effects is to split energies of the molecule that would otherwise be identical in energy. This interaction is responsible for many of the details of spectroscopic properties.
- ADF uses double group symmetry and provides symmetry labels
- spin-orbit coupling is part of the ZORA Hamiltonian
- spin-orbit coupling can be included for most spectroscopic properties
- all electron relativistic basis sets for all elements
- spin-orbit coupling matrix elements with perturbative spin-orbit coupling TDDFT
- singlet-triplet contributions for spin-orbit coupled excitations