A general property operator in 4-component calculations is generated from linear combinations of the basic form:
with the scalar factor \(f\) and the scalar operator \(\Omega\), and where
is one of the following \(4 \times 4\) matrices:
There are 19 basic operator types used in DIRAC, liste in this Table:
Keyword | Operator form | Nr. factors |
---|---|---|
DIAGONAL | \(f I_{4 \times 4} \Omega\) | 1 |
GAMMA5 | \(f \gamma_5 \Omega\) | 1 |
XALPHA | \(f \alpha_x \Omega\) | 1 |
YALPHA | \(f \alpha_y \Omega\) | 1 |
ZALPHA | \(f \alpha_z \Omega\) | 1 |
XSIGMA | \(f \Sigma_x \Omega\) | 1 |
YSIGMA | \(f \Sigma_y \Omega\) | 1 |
ZSIGMA | \(f \Sigma_z \Omega\) | 1 |
XAVECTOR | \(f_1 \alpha_y \Omega_z - f_2 \alpha_z \Omega_y\) | 2 |
YAVECTOR | \(f_1 \alpha_z \Omega_x - f_2 \alpha_x \Omega_z\) | 2 |
ZAVECTOR | \(f_1 \alpha_x \Omega_y - f_2 \alpha_y \Omega_x\) | 2 |
ALPHADOT | \(f_1 \alpha_x \Omega_x + f_2 \alpha_y \Omega_y + f_3 \alpha_z \Omega_z\) | 3 |
XBETASIG | \(f \beta \Sigma_x \Omega\) | 1 |
YBETASIG | \(f \beta \Sigma_y \Omega\) | 1 |
ZBETASIG | \(f \beta \Sigma_z \Omega\) | 1 |
XiBETAAL | \(f i \beta \alpha_x \Omega\) | 1 |
YiBETAAL | \(f i \beta \alpha_y \Omega\) | 1 |
ZiBETAAL | \(f i \beta \alpha_z \Omega\) | 1 |
BETAGAMMA5 | \(f i \beta \gamma_5 \Omega\) | 1 |
Operators are specified by the keyword .OPERATOR with the following arguments:
.OPERATOR
'operator name'
operator type
labels for each component
FACTORS
factors for each component
COMFACTOR
common factor for all components
Note that the arguments following the keyword .OPERATOR must start with a blank. The arguments are optional, except for the operator label.
Keyword | Description | Symmetry | Components | Operators |
---|---|---|---|---|
MOLFIELD | Nuclear attraction integrals | Symmetric | MOLFIELD | \(\Omega_1 = \sum_K V_{iK}\) |
OVERLAP | Overlap integrals | Symmetric | OVERLAP | \(\Omega_1 = 1\) |
BETAMAT | Overlap integrals, only SS-block | Symmetric | BETAMAT | \(\Omega_1 = 1\) |
DIPLEN | Dipole length integrals | Symmetric | XDIPLEN | \(\Omega_1 = x\) |
YDIPLEN | \(\Omega_2 = y\) | |||
ZDIPLEN | \(\Omega_3 = z\) | |||
DIPVEL | Dipole velocity integrals | Anti-symmetric | XDIPVEL | |
YDIPVEL | ||||
ZDIPVEL | ||||
QUADRUP | Quadrupole moments integrals | Symmetric | XXQUADRU | |
XYQUADRU | ||||
XZQUADRU | ||||
YYQUADRU | ||||
YZQUADRU | ||||
ZZQUADRU | ||||
SPNORB | Spatial spin-orbit integrals | Anti-symmetric | X1SPNORB | |
Y1SPNORB | ||||
Z1SPNORB | ||||
SECMOM | Second moments integrals | Symmetric | XXSECMOM | \(\Omega_1 = xx\) |
XYSECMOM | \(\Omega_2 = xy\) | |||
XZSECMOM | \(\Omega_3 = xz\) | |||
YYSECMOM | \(\Omega_4 = yy\) | |||
YZSECMOM | \(\Omega_5 = yz\) | |||
ZZSECMOM | \(\Omega_6 = zz\) |
Keyword | Description |
---|---|
THETA | Traceless theta quadrupole integrals |
CARMOM | Cartesian moments integrals, symmetric integrals, (l + 1)(l + 2)/2 components ( l = i + j + k) |
SPHMOM | Spherical moments integrals (real combinations), symmetric integrals, (2l + 1) components ( m = +0, -1, +1, ..., +l) |
SOLVENT | Electronic solvent integrals |
FERMI C | One-electron Fermi contact integrals |
PSO | Paramagnetic spin-orbit integrals |
SPIN-DI | Spin-dipole integrals |
DSO | Diamagnetic spin-orbit integrals |
SDFC | Spin-dipole + Fermi contact integrals |
HDO | Half-derivative overlap integrals |
S1MAG | Second order contribution from overlap matrix to magnetic properties |
ANGLON | Angular momentum around the nuclei |
ANGMOM | Electronic angular momentum around the origin |
LONMOM | London orbital contribution to angular momentum |
MAGMOM | One-electron contributions to magnetic moment |
KINENER | Electronic kinetic energy |
DSUSNOL | Diamagnetic susceptibility without London contribution |
DSUSLH | Angular London orbital contribution to diamagnetic susceptibility |
DIASUS | Angular London orbital contribution to diamagnetic susceptibility |
NUCSNLO | Nuclear shielding integrals without London orbital contribution |
NUCSLO | London orbital contribution to nuclear shielding tensor integrals |
NUCSHI | Nuclear shielding tensor integrals |
NEFIELD | Electric field at the individual nuclei |
ELFGRDC | Electric field gradient at the individual nuclei, cartesian |
ELFGRDS | Electric field gradient at the individual nuclei, spherical |
S1MAGL | Bra-differentiation of overlap matrix with respect to magnetic field |
S1MAGR | Ket-differentiation of overlap matrix with respect to magnetic field |
HDOBR | Ket-differentiation of HDO-integrals with respect to magnetic field |
NUCPOT | Potential energy at the nuclei |
HBDO | Half B-differentiated overlap matrix |
SQHDO | Half-derivative overlap integrals not to be antisymmetrized |
DSUSCGO | Diamagnetic susceptibility with common gauge origin |
NSTCGO | Nuclear shielding integrals with common gauge origin |
EXPIKR | Cosine and sine integrals |
MASSVEL | Mass velocity integrals |
DARWIN | Darwin type integrals |
CM1 | First order magnetic field derivatives of electric field |
CM2 | Second order magnetic field derivatives of electric field |
SQHDOR | Half-derivative overlap integrals not to be anti-symmetrized |
SQOVLAP | Second order derivatives overlap integrals |
We give here several concrete examples on how to construct operators for various properties.
The kinetic part of the Dirac Hamiltonian may be specified by:
.OPERATOR
'Kin energy'
ALPHADOT
XDIPVEL
YDIPVEL
ZDIPVEL
COMFACTOR
-68.51799475
where -68.51799475 is \(-c/2\).
The speed of light \(c\) is an important parameter in relativistic theory, but its explicit value in atomic units not necessarily remembered. A simpler way to specify the kinetic energy operator is therefore:
.OPERATOR
'Kin energy'
ALPHADOT
XDIPVEL
YDIPVEL
ZDIPVEL
CMULT
COMFACTOR
-0.5
where the keyword CMULT assures multiplication of the common factor -0.5 by \(c\). This option has the further advantage that CMULT follows any user-specified modification of the speed of light, as provided by .CVALUE.
Another example:
.OPERATOR
'B_x'
XAVECTOR
ZDIPLEN
YDIPLEN
CMULT
COMFACTOR
-0.5
The program will assume all operators to be Hermitian and will therefore insert an imaginary phase i if necessary (applies to antisymmetric scalar operators).
If no other arguments are given, the program assumes the operator to be a diagonal operator and expects the operator name to be the component label, for instance:
.OPERATOR
OVERLAP
Another example is the finite perturbation calculation with the \(\hat{z}\) dipole length operator added to the Hamiltonian (don’t forget to decrease the symmetry of your system):
.OPERATOR
ZDIPLEN
COMFACTOR
0.01
Here is an example where the Fermi-contact (FC) integrals for a certain nucleus are added to the Hamiltonian in a finite-field calculation. Let’s assume you are looking at a PbX dimer (order in the .mol file: 1. Pb, 2. X) and you want to add to the Dirac-Coulomb **HAMILTONIAN the FC integrals for the Pb nucleus as a perturbation with a given field-strength (FACTORS).
Important note: The raw density values obtained after the fit of your finite-field energies need to be scaled by \(\frac{1}{\frac{4 \pi g_{e}}{3}} = \frac{1}{8.3872954891254192}\), a factor that originates from the definition of the operator for calculating the density at the nucleus:
**HAMILTONIAN
.OPERATOR
'Density at nucleus'
DIAGONAL
'FC Pb 01'
FACTORS
-0.000000001
Here is next example of how-to calculate the electron density at the nucleus as an expectation value \(\langle 0 \vert \delta(r-R) \vert 0 \rangle\) for a Dirac-Coulomb HF wave function including a decomposition of the molecular orbital contributions to the density:
**DIRAC
.WAVE FUNCTION
.PROPERTIES
**HAMILTONIAN
**WAVE FUNCTION
.SCF
**PROPERTIES
.RHONUC
*EXPECTATION
.ORBANA
*END OF
In the following example I am calculating a cartesian moment expectation value \(\langle 0 \vert x^1 y^2 z^3 \vert 0 \rangle\) for a Levy-Leblond HF wave function:
**DIRAC
.WAVE FUNCTION
.PROPERTIES
**HAMILTONIAN
.LEVY-LEBLOND
**WAVE FUNCTION
.SCF
**PROPERTIES
*EXPECTATION
.OPERATOR
CM010203
*END OF