NMR Chemical Shifts

NMR Chemical shifts have been implemented [1] [2] [3] [4] [5] in a separate property program NMR. It requires the adf.rkf (TAPE21) result file from an ADF calculation. The NMR module can be combined with the ZORA treatment for relativistic effects and with Spin-Orbit effects, making it suitable for treatment of heavy elements. See also the general review on relativistic computations of NMR parameters [6].

Two GUI tutorials on NMR Chemical shift calculations are currently available:

Important notes

adf.rkf (TAPE21) and TAPE10

NMR requires an ASCII input file and adf.rkf (TAPE21) and TAPE10 result files from an ADF calculation on the molecule to be analyzed. The ADF result files adf.rkf (TAPE21) and TAPE10 must be present with names TAPE21 and TAPE10 in the directory where you execute NMR, or alternatively one can use the keys ADFFile and TAPE10File to point to an adf.rkf file and TAPE10 file. Use the keywords SAVE TAPE10 in the adf calculation in order to obtain a TAPE10 result file.

Recalculation of TAPE10 by NMR

Warning: the NMR property program will not always give the correct result for every SCF potential in the ADF calculation, like for example the SAOP potential, or if one uses COSMO in the ADF calculation, if one lets the NMR recalculate TAPE10. This is due to the GIAO method used in this program, which requires the calculation of the SCF potential, which is not done correctly for potentials, other than the standard LDA and GGA potentials. To obtain correct results one should, in addition to the use of TAPE21, also use TAPE10 that ADF generates, using the keywords SAVE TAPE10, and use it as input for the NMR property program. On TAPE10 the SCF potential is written, which is read in by the NMR program.

Atomic calculation

NMR calculations on 1 atom must have symmetry NOSYM.

Spin-orbit coupling

NMR calculations on systems computed by ADF with Spin Orbit relativistic effects included must have used NOSYM symmetry in the ADF calculation. NMR can also be combined with ADF ZORA calculations. The NMR program reads from adf.rkf (TAPE21) the relativistic option that is used in the ADF calculation, and will use the same relativistic option in the NMR calculations.

Kernel used in spin-orbit coupling

One can include an improved exchange-correlation kernel, as was implemented by J. Autschbach [7], but this is not used by default, see subkey USE. This is an important option that should be seriously considered and has been advocated in Ref. [7].

Bug spin-orbit part ADF2008 - ADF2013

In the ADF2008.01 a bug was introduced in the spin-orbit part of the calculated chemical shielding, which caused the calculated chemical shielding to be gauge dependent. This bug is relevant for spin-orbit coupled calculations for ADF versions ADF2008-ADF2013. In ADF2014 this bug has been fixed.

Unscaled ZORA default ADF2014

There is gauge dependence if the scaled ZORA method is used in the calculation of NMR chemical shieldings. Therefore the default method for NMR chemical shielding calculations is changed in ADF2014 to use the unscaled ZORA method.

SAOP

The use of the model SAOP potential leads to isotropic chemical shifts which are substantially improved over both LDA and GGA functionals, and of similar accuracy as results with a self-interaction-corrected functional (SIC), see [8]. SAOP is computationally expedient and routinely applicable to all systems, requiring virtually the same computational effort as LDA and GGA calculations.

NICS

The Nucleus-Independent Chemical Shift (NICS) can be calculated at any point in the molecule.

Hybrids

Starting from ADF2009.01 Hartree-Fock and the hybrid potentials can used in combination with NMR chemical shielding calculations. see Refs. [9] [10]. Use SAVE TAPE10 in the ADF calculation. The use of frozen cores and hybrids gives gauge dependent results for the NMR chemical shieldings, therefore the NMR program will stop in this case.

Meta-GGA’s and meta-hybrids

Meta-GGA’s and meta-hybrids should not be used in combination with NMR chemical shielding calculations. The results are wrong due to an incorrect inclusion of GIAO terms.

Input options

The input file for NMR uses the block key NMR, with several (optional) subkeys, each having a series of options. For analysis a separate block key can be used.

$AMSBIN/nmr << eor
ADFFile adffile
TAPE10File tape10file
ZSOAO2007
RECALCULATETAPE10
NMR
End
Analysis
End
eor

ADFFile adffile

Path to adf.rkf (TAPE21) file from which nmr reads data and to which nmr writes data. Default TAPE21.

TAPE10File tape10file

Path to TAPE10 file from which nmr reads data. Default TAPE10.

ZSOAO2007 keyword

In the ADF2008.01 a bug was introduced in the spin-orbit part of the calculated chemical shielding, which caused the calculated chemical shielding to be gauge dependent. This bug is relevant for spin-orbit coupled calculations for ADF versions ADF2008-ADF2013. Workaround in ADF versions ADF2008-ADF2013 is to include the keyword ZSOAO2007 in the NMR part of the input, which causes a one-center approximation to be used. The bug has been fixed in ADF2014 by introducing an extra gauge-correction term for the spin-orbit coupled part. One can still get the (slightly) incorrect results in ADF2014 by using the keyword WRONGSOGAUGE, and not including ZSOAO2007.

REACALCULATETAPE10 keyword

If there is no TAPE10 present the NMR program will stop. One can use the key REACALCULATETAPE10 such that TAPE10 will be recalculated by the NMR module. Not recommended to be used. Better use ‘SAVE TAPE10’ in the ADF calculation, and use this TAPE10 as input for NMR.

NMR block key

NMR
  Out  OutOptions
  Calc CalcOptions
  Use UseOptions
  U1K  U1KOptions
  Nuc  NucOptions
  Atoms  AtomsOptions
  Ghosts  GhostsOptions
  Analysis  AnalysisOptions
End
Out OutOptions

The subkey Out controls printed output. Its options specify the details by their (optional) presence. The following OutOptions are recognized (Default ISO TENS):

All
Implies all the other options
ISO
Isotropic shielding constants
Tens
Shielding tensors
Eig
Eigenvectors
U1
The U1 matrix
F1
The first order change in the Fock matrix
S1
The first order change in the Overlap matrix
AOP
The paramagnetic AO matrix (= the matrix in the representation of elementary atomic basis functions)
AOD
The diamagnetic AO matrix
AOF
The Fermi-contact AO matrix
REFS
Literature references
INFO
General information
Calc CalcOptions

The subkey Calc controls what is actually calculated. The following options are available (Default ALL):

All
Implies all of the other options to this key
Para
The paramagnetic part
Dia
The diamagnetic part
FC
The Fermi-contact part in case of the Pauli Hamiltonian
SO
The Fermi-contact part in case of the ZORA Hamiltonian
Use UseOptions

The subkey Use controls some optional options (default none)

FXC
Improves the exchange-correlation kernel used, as was implemented by J. Autschbach [7]. Important only in case of spin-orbit coupled calculations. This may give some (small) gauge dependent results when using this. Important option that should be seriously considered and has been advocated in Ref. [7].
SCALED
Implies the scaled ZORA method, which gives (slightly) gauge dependent results, as was shown in Ref. [7]. Note that in case of the ZORA Hamiltonian default the unscaled ZORA method is used. For chemical shifts, only compare results with the same options.
SO1C
Before ADF2008.01 in the the spin-orbit term a 1-center approximation was used, which does not suffer from gauge dependence. This 1-center approximation can be used with USE SO1C.
U1K U1KOptions

The subkey U1K determines which terms are included in the calculation of the U1 matrix (first order changes in MO coefficients). Options (Default Best):

Best
The best (recommended) options for each relativistic option are included for this subkey. Implies ‘None’ for non-relativistic and scalar relativistic ZORA, ‘SO’ + ‘SOFULL’ for spin-orbit coupled ZORA, and ‘MV + ‘Dar’ for the Pauli Hamiltonian.
None
Implies none of the other options to this key.
All
Implies all the other options to this key.
MV
The mass-velocity term .
Dar
The Darwin term.
ZMAN
The Spin-Zeeman term. Can be included only in case of spin-orbit coupled Pauli Hamiltonian.
SO
ZORA spin-orbit part.
SOFULL
ZORA spin-orbit part.

Note: for chemical shifts, only compare results with the same options. If the subkey U1K is used with the option ALL in the ZORA calculation, then the scaled ZORA orbital energies are used in the making of the U1 matrix, which is not recommended. Recommended is to use ‘U1K Best’ in all cases, which uses plain ZORA orbital energies in the making of the U1 matrix.

NUC NucOptions

This subkey NUC specifies for which nuclei the NMR shielding is calculated. Default all nuclei are calculated, i.e. as for omitting the subkeys ATOMS and NUC. Else you may use this options by simply typing Nuc in the NMR block (without any further data); this means: for no nuclei at all. Alternatively you may type the index of the atom(s) you want to see analyzed. Default all nuclei are calculated, i.e. as for omitting this subkey.

Example:

NUC 2 1

The numbers refer to the internal numbering of the nuclei as it appears somewhere early in the general ADF output. This internal numbering is also the internal NMR numbering, but it is not necessarily the same as the input ordering. Use the subkey ATOMS to specify the nuclei according to this input ordering in the ADF calculation.

Note that the number of nuclei has a significant consequence for the total CPU time.

Atoms AtomsOptions

This subkey ATOMS specifies for which nuclei the NMR shielding is calculated. Default all nuclei are calculated, i.e. as for omitting the subkeys ATOMS and NUC.

Example:

ATOMS 2 1

The numbers refer to the input ordering in the ADF calculation. Use the subkey NUC to specify the nuclei according to the internal NMR numbers of the atoms.

GHOSTS

The subkey GHOSTS is a block type subkey. The format is:

Ghosts
  xx1 yy1 zz1
  xx2 yy2 zz2
  ...
SubEnd

With this key, the user can specify ANY point(s) within the molecule at which the shielding is to be calculated (whatever the physical meaning of this shielding is). One can think of those points as neutrons within the molecule. There is a publication by P. Schleyer et al. using a similar feature (J. Am. Chem. Soc. 118, 6317, 1996). They call it NICS, Nucleus-Independent Chemical Shift. Note that the NICS value is minus 1 times the isotropic part of the shielding tensor that is calculated at these points:

xx1 yy1 zz1

real numbers that specify the Cartesian coordinates of ‘ghost’ 1, etc.

The coordinates have to be specified in Angstrom. The same set of coordinates has to be specified as ‘point charges with charge zero’ using the subsubkey MultipolePotential in the ElectrostaticEmbedding subkey of the key System in the AMS block. This is necessary in order to allow the appropriate distribution of integration points around the ghosts.

E.g., if you want to have two ‘ghosts’ with the coordinates xx1 yy1 zz1 and xx2 yy2 zz2 then you must also have in the input the key MultipolePotential key as follows:

System
  ...
  ElectrostaticEmbedding
    MultipolePotential
      Coordinates
        xx1 yy1 zz1   0.0
        xx2 yy2 zz2   0.0
      End
    End
  End
End

(the last number is the charge at these coordinates - zero).

Eventually, this step should be programmed internally but for now the procedure outlined above works. No check is done to verify whether those ‘point charges’ are taken care of or not, but their omission leads to unpredictable results.

Only Cartesian coordinates are possible for ghosts, even if the atoms were originally specified using internal coordinates. This shouldn’t be a problem, though (e.g., one could start an ADF run of the molecule of interest, and get very soon the Cartesian coordinates of the atoms in the output. This run would then be aborted, and restarted with the ghosts specified as desired.) The ghosts are numbered in the output as NNUC+1, NNUC+2 … where NNUC is the total number of nuclei in this molecule. Default: no ghosts.

Analysis AnalysisOptions
The subkey Analysis controls the MO analysis. After the word (subkey) Analysis you type an integer, which then specifies that the first so many MOs are to be analyzed. Default no Analysis. The value of this analysis subkey in the block key NMR is somewhat limited. The separate ANALYSIS block key can give more analysis of the NMR chemical shielding.

Analysis block key

The NMR shielding tensor, can be analyzed in detail, see Refs. [12] [13] [6]. For the analysis option with the ANALYSIS block key there are some restrictions. In the ADF calculation all electron basis sets should have been used, and SYMMETRY NOSYM. Can not be used in case of non-relativistic calculations. The ADF calculation should use relativistic scalar ZORA or relativistic spinorbit ZORA. In case of scalar relativistic ZORA the keyword FAKESO should be added to the NMR input (outside of the NMR or Analysis block keys). The analysis utilizes the ZORA spin-orbit branch of the NMR code. For scalar ZORA, the NMR analysis contributions will appear in equivalent pairs for spin-orbitals even if the ADF calculation is closed-shell spin restricted. The MO numbering then also reflects this doubling of MOs. In the analysis, canonical MOs number 1 and 2 are the alpha and beta spin ADF MO 1, canonical MOs number 3 and 4 correspond to the alpha and beta spin ADF MO 2, and so on. In case of spinorbit relativistic ZORA the keyword FAKESO should not be included. Curly brackets are not part of the key. For an NBO analysis of NMR, see the section on NBO analysis.

Analysis
  print threshold
  canonical
  {components}
End
{FakeSO}
print threshold
The print keyword selects printout of contributions relative to the total diamagnetic, paramagnetic. For example in case of ‘print 0.01’ only contributions greater than 1% are printed. Set to zero to print ALL contributions.
canonical
It enables an analysis of the shielding in terms of the canonical MOs.
components
The components keyword is optional and enables an analysis not only of the isotropic shielding but also of the diagonal Cartesian components of the tensor XX, YY, and ZZ). In order to analyze the principal shielding tensor components with canonical MOs you can calculate the shielding tensor first with the NMR code, rotate the molecule such that the principal axes system aligns with the Cartesian coordinate system, and then repeat the NMR calculation with the analysis features switched on.

Paramagnetic NMR Chemical Shifts

Knowledge of the g-tensor and hyperfine A-tensor can be used in the prediction and analysis of paramagnetic NMR shifts, see Refs. [11] [14] [15]. Because of the dependence on the g-tensor, prediction of paramagnetic NMR (pNMR) shifts is not straightforward, as the dependence of the pNMR shift on excess \(\alpha\) or \(\beta\) electron spin density at a nucleus is to be combined with the sign and magnitude of the isotropic g-value. Ref. [11] describes in detail how to calculate pNMR contact chemical shifts and pNMR pseudo-contact chemical shifts, using the ADF program.

Of course, like for NMR chemical shielding of closed shell molecules, one also needs to calculate of the so called orbital dependent part of the NMR chemical shielding of the open shell molecule. For open shell molecules in ADF this can only be calculated without taken into account spin-orbit coupling:

$AMSBIN/ams << eor
  ...
  Engine ADF
    SpinPolarization vlaue
    Unrestricted Yes
    Relativity
      Level Scalar
    End
    Basis
      Core None
    End
    SAVE TAPE10
  EndEngine
eor

$AMSBIN/nmr << eor
  ADFFile path/to/adf.rkf
  TAPE10File path/to/TAPE10
  ALLINONE
  nmr
     u1k best
     calc all
     out iso tens
  end
eor

Note that one can only do this at the scalar relativistic ZORA level, one needs to use all electron basis sets, and one needs to include the key ALLINONE in the input for NMR,

References

[1]G. Schreckenbach and T. Ziegler, The calculation of NMR shielding tensors using GIAO’s and modern density functional theory, Journal of Physical Chemistry 99, 606 (1995)
[2]G. Schreckenbach and T. Ziegler, The calculation of NMR shielding tensors based on density functional theory and the frozen-core approximation, International Journal of Quantum Chemistry 60, 753 (1996)
[3]G. Schreckenbach and T. Ziegler, Calculation of NMR shielding tensors based on density functional theory and a scalar relativistic Pauli-type Hamiltonian. The application to transition metal complexes, International Journal of Quantum Chemistry 61, 899 (1997)
[4]S.K. Wolff and T. Ziegler, Calculation of DFT-GIAO NMR shifts with inclusion of spin-orbit coupling, Journal of Chemical Physics 109, 895 (1998)
[5]S.K. Wolff, T. Ziegler, E. van Lenthe and E.J. Baerends, Density functional calculations of nuclear magnetic shieldings using the zeroth-order regular approximation (ZORA) for relativistic effects: ZORA nuclear magnetic resonance , Journal of Chemical Physics 110, 7689 (1999)
[6](1, 2) J. Autschbach and S. Zheng, Relativistic computations of NMR parameters from first principles: Theory and applications, Annual Reports on NMR Spectroscopy 67, 1 (2009)
[7](1, 2, 3, 4, 5) J. Autschbach, The role of the exchange-correlation response kernel and scaling corrections in relativistic density functional nuclear magnetic shielding calculations with the zeroth-order regular approximation, Molecular Physics 111, 2544 (2013)
[8]J. Poater, E. van Lenthe and E.J. Baerends, Nuclear magnetic resonance chemical shifts with the statistical average of orbital-dependent model potentials in Kohn.Sham density functional theory, Journal of Chemical Physics 118, 8584 (2003)
[9]M. Krykunov, T. Ziegler and E. van Lenthe, Hybrid density functional calculations of nuclear magnetic shieldings using Slater-type orbitals and the zeroth-order regular approximation, International Journal of Quantum Chemistry 109, 1676 (2009)
[10]M. Krykunov, T. Ziegler and E. van Lenthe, Implementation of a hybrid DFT method for calculating NMR shieldings using Slater-type orbitals with spin-orbital coupling included. Applications to 187 Os, 195 Pt and 13 C in heavy metal complexes, Journal of Physical Chemistry A 113, 11495 (2009)
[11](1, 2) J. Autschbach, S. Patchkovskii, and B. Pritchard, Calculation of Hyperfine Tensors and Paramagnetic NMR Shifts Using the Relativistic Zeroth-Order Regular Approximation and Density Functional Theory, Journal of Chemical Theory and Computation 7, 2175 (2011)
[12]J. Autschbach, Analyzing NMR shielding tensors calculated with two-component relativistic methods using spin-free localized molecular orbitals, Journal of Chemical Physics 128, 164112 (2008)
[13]J. Autschbach and S. Zheng, Analyzing Pt chemical shifts calculated from relativistic density functional theory using localized orbitals: The role of Pt 5d lone pairs, Magnetic Resonance in Chemistry 46, S45 (2008)
[14]S. Moon, and S. Patchkovskii, First-principles calculations of paramagnetic NMR shifts, in Calculation of NMR and EPR parameters, ISBN13: 9783527307791, M. Kaupp, M. Bühl, V.G. Malkin, Editors, (Wiley, Weinheim, 2004).
[15]P. Hrobárik, Ro. Reviakine, A.V. Arbuznikov, O.L. Malkina, V.G. Malkin, F.H. Köhler, and M. Kaupp, Density functional calculations of NMR shielding tensors for paramagnetic systems with arbitrary spin multiplicity: Validation on 3d metallocenes, Journal of Chemical Physics 126, 024107 (2007)