Calculating the hyperfine coupling constant
The hyperfine coupling constant (cf. hyperfine splitting) describes the interaction between the nuclear magnetic dipole moment and the magnetic field generated by the electrons (i.e. the nuclear spin-electron spin coupling. The hyperfine coupling constant is calculated using LHYPERFINE [1]. The hyperfine splitting often includes the interaction between the nuclear quadrupole moment and the electric field gradient (EFG), which is calculated separately using LEFG and the description in performing an EFG calculation. The hyperfine splitting can be measured using electron paramagnetic resonance (EPR), also called electron spin resonance (ESR), and in atomic spectroscopy. The theory is already covered in the NMR category page and LHYPERFINE page, so it will not be reiterated here.
No additional files are required besides the four standard POSCAR, POTCAR, INCAR, and KPOINTS, unless specifically mentioned. It is important to have a well-converged structure, as the hyperfine coupling constant can be very sensitive to it. It is calculated post-SCF.
Input
The hyperfine coupling constants are calculated using LHYPERFINE
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= 2.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).