A general property operator in 4-component calculations is generated from linear combinations of the basic form:

\[\hat{P} = f M_{4 \times 4} \Omega\]

with the scalar factor \(f\) and the scalar operator \(\Omega\), and where

\[M_{4 \times 4}\]

is one of the following \(4 \times 4\) matrices:

\[I_{4 \times 4}, \gamma_5,\]\[\alpha_x, \alpha_y, \alpha_z\]\[\Sigma_x, \Sigma_y, \Sigma_z\]\[\beta \Sigma_x, \beta \Sigma_y, \beta \Sigma_z\]\[i \beta \alpha_x, i \beta \alpha_y, i \beta \alpha_z\]

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):

\[\hat{H} = \hat{H}_0 + 0.01 \cdot \hat{z}\]

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