# Nuclear spin-rotation constants¶

## Introduction¶

In this tutorial we introduce the calculation of nuclear spin-rotation tensors as given in the DIRAC code, based on the theoretical
developments by I. Agustín Aucar *et al.* For details, you are welcome to consult [Aucar_JCP2012].

The nuclear spin-rotation (SR) tensor elements of a nucleus \(N\) are given by

They have three terms: the first of them is independent of the electronic variables, whereas the second and third are given by an expectation value (EV) and a linear response (LR) function, respectively.

In tensorial notation, the nuclear spin-rotation contributions of a nucleus \(N\) are given, in SI units, by:

where \({\bf R}_{MN}\) is the position of nucleus \(M\) with respect to the position of nucleus \(N\); \({\bf R}_{N,CM}\) is the position of nucleus \(N\) with respect to the molecular center of mass; \({\bf I}\) is the inertia tensor of the molecule and \({\bf J}_e = \left({\bf r} - {\bf R}_{CM} \right) \times {\bf p}+{\bf S}_e\) is the electronic total angular momentum.

For molecules at their equilibrium geometry, it is found that the sum of the first two terms, \({\bf M}_N^{NU} (eq) + {\bf M}_N^{EV} (eq)\), is equal to a new tensor which is independent of the electronic dynamics (see Eq. 60 of [Aucar_JCP2012]),

Therefore, the electronic dependence of the nuclear spin-rotation tensor of a nucleus \(N\) in a molecule in equilibrium is completely given by its linear response term, \({\bf M}_N^{LR}\) (see Eq. 59 of [Aucar_JCP2012]).

Note

The current implementation gives by default (.PRINT values up to 3) results only for molecules in equilibrium.

## Application to the HF molecule¶

As an example, we show a calculation of the SR constant of the fluorine nucleus at the Hydrogen fluoride molecule. The input file spinrot.inp is given by

```
**DIRAC
.TITLE
Spin-rotation constant
.WAVE FUNCTION
.PROPERTIES
**HAMILTONIAN
.URKBAL
**INTEGRALS
*READIN
.UNCONTRACT
**WAVE FUNCTIONS
.SCF
*SCF
.EVCCNV
1.0E-7
**PROPERTIES
.SPIN-R
.PRINT
4
*LINEAR RESPONSE
.THRESH
1.0D-7
*END OF
```

whereas the molecular input file HF_cv3z.mol is

```
INTGRL
Hydrogen fluoride. Experimental bond length: 0.917 A
dyall.cv3z basis set
C 2 A .10D-15
9.0 1
F 0.00000000000 0.00000000000 0.00000000000 Isotope=19
LARGE BASIS dyall.cv3z
1.0 1
H 0.00000000000 0.00000000000 0.91700000000 Isotope=1
LARGE BASIS dyall.cv3z
FINISH
```

The calculation is run using:

```
pam --inp=spinrot --mol=HF_cv3z
```

As a result, SR constants are obtained at the coupled Hartree-Fock level. The code also works at the DFT level.

It is also interesting to note that as .URKBAL is requested in the present calculation, the results will be very close to those obtained using the RKB prescription, because the diamagnetic-like (e-p) contributions to \(M_{N,\alpha\beta}^{LR}\) are almost zero (see Eq. 60 of [Aucar_JCP2012]).

## Reading the output file¶

As the .PRINT flag in the input file is set to 4, the results are fully detailed and given in both kHz and ppm units.

The SR constant of the fluorine nucleus, in kHz, will look like:

```
Spin-rotation constants (kHz) for F
-----------------------------------
Nuclear g-value: 5.257736
Total spin-rotation constant (SRC) : -317.17730994
Nuclear contribution to SRC (M^nuc) : 52.47420202
Electronic contribution to SRC (M^LR) : -369.65151197
M^LR-L(e-e) : -381.52309467
M^LR-S(e-e) : 11.85860365
M^LR-L(e-p) : 4.27896875
M^LR-S(e-p) : -4.26598970
******************************************************
******** M^NU ********: 55.00521470
******** M^EV ********: -2.51499774
******** M^LR ********: -369.65151197
******** M^total ********: -317.16129500
******************************************************
```

One should recall that in the case of linear molecules, as the present one, only one tensor element is printed out (the spin-rotation constant) because for these molecules the SR tensor has only two equal and non-zero elements.

As it can be seen, the total SR constant of the fluorine nucleus is given, and then separated in its two terms (nuc and LR). The latter contribution, the linear response function, is further separated in their (e-e) and (e-p) parts, as well as in their \(\mathbf{L}\) and \(\mathbf{S}\) parts. It is seen how the (e-p) contribution to the linear response term of the SR constants is almost zero.

In addition, the results are fully detailed, showing the three terms mentioned above (NU, EV and LR).

Finally, if one is interested in the relationship between the SR constants and .SHIELDING (see [Aucar_JPCL2016] and [AucarChap2019]), which is given, in SI units, by

then we print also the SR constants in ppm, by multiplying each SR tensor element \(M_{N,\alpha\beta}\) (given in kHz) by \(\frac{m_p I_{\beta \gamma}}{g_N}\,\frac{2\, \pi \, 10^9}{m_e \, \hslash}\), where \(m_p\) and \(m_e\) are the proton and electron masses, respectively, \(\hslash\) is the reduced Planck’s constant, and \(g_N\) is the g-factor of nucleus N.

In this way, one obtains the following output:

```
Spin-rotation constants (in ppm) for F
--------------------------------------
Total spin-rotation constant (SRC) : -88.19431889
Nuclear contribution to SRC (M^nuc) : 14.59097597
Electronic contribution to SRC (M^LR) : -102.78529486
M^LR-L(e-e) : -106.08630700
M^LR-S(e-e) : 3.29740318
M^LR-L(e-p) : 1.18981000
M^LR-S(e-p) : -1.18620104
******************************************************
******** M^NU ********: 15.29474932
******** M^EV ********: -0.69932024
******** M^LR ********: -102.78529486
******** M^total ********: -88.18986578
******************************************************
```