Calculating the hyperfine coupling constant
Hyperfine coupling
The hyperfine coupling constants are calculated using LHYPERFINE, i.e. the hyperfine splitting [1]. Specifically, the coupling between the nuclear magnetic dipole moment and the magnetic field generated by the electrons (or nuclear spin-electron spin coupling) is referred to. For the interaction between the nuclear quadrupole moment and the electric field gradient, see LEFG and the description in insert link back to efg above when hyperfine made into its own how-to page. The hyperfine splitting can be measured using electron paramagnetic resonance (EPR), also called electron spin resonance (ESR), as well as in atomic spectroscopy.
There is one additional keyword that must be defined:
- NGYROMAG defines the nuclear gyromagnetic ratios for each element in your POSCAR file. The defaults are set to 1, which will return meaningless results. Reasonable values may be found here.
An example INCAR file is given here:
ENCUT = 500 # Plane-wave energy cutoff in eV ISMEAR = 0; SIGMA = 0.01 # Defines the type of smearing; smearing width in eV EDIFF = 1E-8 # Energy cutoff criterion for the SCF loop, in eV PREC = Accurate # Sets the "precision" mode LHYPERFINE = .TRUE. # Turns on calculating the hyperfine coupling tensor NGYROMAG = 10.7084 42.577478461 # Specifies the nuclear gyromagnetic ratios for the ions - C and H in this case ISPIN = 2 # Turns on spin-polarization - noncollinear can also be used
| Important: Make sure to replace the NGYROMAG in the INCAR with the values for the isotopes in your system. |
Output
You can find the output for the hyperfine calculation in the OUTCAR file after the SCF cycle finishes. The total magnetic moment is listed, then the Fermi contact term:
Total magnetic moment S= 1.0000000 Fermi contact (isotropic) hyperfine coupling parameter (MHz) ------------------------------------------------------------- ion A_pw A_1PS A_1AE A_1c A_tot ------------------------------------------------------------- 1 - - - - - 2 - - - - - -------------------------------------------------------------
Note the A_tot does not include the core contribution term A_1c [2]. The dipolar hyperfine coupling parameter comes next and finally the total hyperfine coupling parameter
Dipolar hyperfine coupling parameters (MHz) --------------------------------------------------------------------- ion A_xx A_yy A_zz A_xy A_xz A_yz --------------------------------------------------------------------- 1 - - - - - - 2 - - - - - - --------------------------------------------------------------------- Total hyperfine coupling parameters after diagonalization (MHz) (convention: |A_zz| > |A_xx| > |A_yy|) ---------------------------------------------------------------------- ion A_xx A_yy A_zz asymmetry (A_yy - A_xx)/ A_zz ---------------------------------------------------------------------- 1 - - - - 2 - - - - ---------------------------------------------------------------------
Recommendations and advice
Distinct from the chemical shielding and EFG, the hyperfine constant is less dependent on EDIFF and ENCUT, generally converging relatively quickly with respect to both. However, it is extremely strongly influenced by the method used. HSE06 was found to give values close to experimental values for molecular systems [1].
References
- ↑ a b K. Szasz, T. Hornos, M. Marsman, and A. Gali, Hyperfine coupling of point defects in semiconductors by hybrid density functional calculations: The role of core spin polarization, Phys. Rev. B, 88, 075202 (2013).
- ↑ O. V. Yazyev, I. Tavernelli, L. Helm, and U. R. Roethlisberger, Core spin-polarization correction in pseudopotential-based electronic structure calculations, Phys. Rev. B 71, 115110 (2006).