LCHIMAG

LCHIMAG = .TRUE. | .FALSE.
Default: LCHIMAG = .FALSE. 

Description: calculate the chemical shifts by means of linear response.

---

For LCHIMAG=.TRUE., VASP calculates the chemical shift tensors.

The chemical shielding tensor is defined as:

$\sigma_{ij}(\mathbf{R}) = - \frac{ \partial B^{\mathrm{in}}_i(\mathbf{R})}{ \partial B^{\mathrm{ext}}_j}$

Here $\mathbf{R}$ denotes the atomic nuclear site, $i$ and $j$ denote cartesian indices, $B^{\mathrm{ext}}$ an applied DC external magnetic field and $B^{\mathrm{in}}(\mathbf{R})$ the induced magnetic field at the nucleus. NMR experiments yield information on the symmetric part of the tensor. Typical NMR experiments yield information on the shielding relative to that of a reference compound:

$\delta_{ij}(\mathbf{R}) = \sigma_{ij}^{\mathrm{ref}} - \sigma_{ij}(\mathbf{R})$

In this (approximate) formula $\sigma_{ij}^{\mathrm{ref}}$ is the isotropic shielding of the nucleus in the reference compound. $\delta_{ij}(\mathbf{R})$ is the chemical shift tensor.

In VASP the chemical "shift" tensor is calculated as:

$\delta_{ij}(\mathbf{R})\mathrm{[VASP]} = \frac{ \partial B^{\mathrm{in}}_i(\mathbf{R})}{ \partial B^{\mathrm{ext}}_j}$

This is minus the shielding tensor. It is not the true chemical shift tensor. To convert it to the real shift tensor one should add the reference shielding:

$\delta_{ij}(\mathbf{R}) = \sigma_{ij}^{\mathrm{ref}} + \delta_{ij}(\mathbf{R})\mathrm{[VASP]}$

VASP calculates chemical "shifts" for non-metallic crystalline systems using the linear response method of Yates, Pickard, and Mauri

Input

An example INCAR is given below:

PREC = Accurate        
ENCUT = 600.0          
EDIFF = 1E-8           
ISMEAR = 0; SIGMA = 0.1 
LREAL = A              

LCHIMAG = .TRUE.       
DQ = 0.001             
ICHIBARE = 1           
LNMR_SYM_RED = .TRUE.  
NLSPLINE = .TRUE.      

Warning: The chemical shifts are calculated from the orbital magnetic response under the assumption that the system is an insulator. Smearing schemes intended for metals can generate nonsense.

The chemical shift is calculated via the induced current (cf. LWRTCUR) . It has contributions from the plane wave grid and one-center contributions (the induced field at the center of a PAW sphere due to the augmentation current inside that sphere). Two-center contributions (induced fields due to augmentation currents in other PAW spheres) are neglected by default. These contributions can be switched on using LLRAUG.

Output

The isotropic chemical shieldings are printed towards the end of the OUTCAR file, after the self-consistent calculation has finished. The chemical shift tensors both before and after space group symmetrization. These are the absolute tensors for the infinite lattice, excluding core contributions. They can be searched for under the UNSYMMETRIZED TENSORS and SYMMETRIZED TENSORS after Absolute Chemical Shift tensors. Additionally, the magnetic susceptibility is printed shortly after and found under ORBITAL MAGNETIC SUSCEPTIBILITY.

To obtain the full absolute tensors requires adding both the $\mathbf{G=0}$ contribution and the contributions due to the core electrons. The latter consist of contributions for each chemical species separately (depending on POTCAR) and a global $\mathbf{G=0}$ susceptibility contribution.

The reference shift experienced by the core is given first:

  Core NMR properties

  typ  El   Core shift (ppm)
 ----------------------------
    1  C     -200.5098801
 ----------------------------

  Core contribution to magnetic susceptibility:     -0.31  10^-6 cm^3/mole
 --------------------------------------------------------------------------

Important: The chemical shieldings calculated are the negative of the chemical shift. Note that the isotropic chemical shift $\delta_{\mathrm{iso}}\mathrm{[VASP]}$ (ISO_SHIFT) is actually minus the isotropic shielding. To make it a real shift one should add the reference shielding.

Next, the tensor is processed and its chemical shielding anisotropy (CSA) characteristics are printed in the OUTCAR. The tensor is symmetrized ($\sigma_{ij} = \sigma_{ji}$ is enforced) and diagonalized. From the three diagonal values the isotropic chemical "shift" $\delta_{\mathrm{iso}}\mathrm{[VASP]}$, span $\Omega$, and skew $\kappa$ are calculated and printed see Ref. for unambiguous definitions. Note that $\kappa$ is ill-defined if $\Omega = 0$. Units are ppm, except for the skew. A typical output is given below:

                                                                                                          
   ---------------------------------------------------------------------------------
    CSA tensor (J. Mason, Solid State Nucl. Magn. Reson. 2, 285 (1993))
   ---------------------------------------------------------------------------------
               EXCLUDING G=0 CONTRIBUTION             INCLUDING G=0 CONTRIBUTION
           -----------------------------------   -----------------------------------
    ATOM    ISO_SHIFT        SPAN        SKEW     ISO_SHIFT        SPAN        SKEW
   ---------------------------------------------------------------------------------
    (absolute, valence only)
       1    4598.8125      0.0000      0.0000     4589.9696      0.0000      0.0000
       2     291.5486      0.0000      0.0000      282.7058      0.0000      0.0000
       3     736.5979    344.8803      1.0000      727.7550    344.8803      1.0000
       4     736.5979    344.8803      1.0000      727.7550    344.8803      1.0000
       5     736.5979    344.8803      1.0000      727.7550    344.8803      1.0000
   ---------------------------------------------------------------------------------
    (absolute, valence and core)
       1   -6536.1417      0.0000      0.0000    -6547.9848      0.0000      0.0000
       2   -5706.3864      0.0000      0.0000    -5718.2296      0.0000      0.0000
       3   -2369.4015    344.8803      1.0000    -2381.2446    344.8803      1.0000
       4   -2369.4015    344.8803      1.0000    -2381.2446    344.8803      1.0000
       5   -2369.4015    344.8803      1.0000    -2381.2446    344.8803      1.0000
   ---------------------------------------------------------------------------------
    IF SPAN.EQ.0, THEN SKEW IS ILL-DEFINED
   ---------------------------------------------------------------------------------

The isotropic chemical shielding for each atom, excluding and including G=0 contributions, as well as the span and skew (descriptions of asymmetry), follow. Finally, core contributions are taken into account for the ISO_SHIFT, SPAN, and SKEW.

Important: The columns excluding the $\mathbf{G=0}$ contribution are useful for supercell calculations on molecules. The columns including the $\mathbf{G=0}$ contribution are for crystals. The upper block gives the shielding due to only the electrons included in the SCF calculation. The lower block has the contributions due to the frozen PAW cores added. These core contributions are rigid . They depend on POTCAR and are isotropic, i.e. affect neither SPAN nor SKEW.

By default the orbital magnetic susceptibility is calculated using the so-called pGv-approximation, i.e. Eqs. 46-48 of Yates et al. . As of vasp.6.4.0 also the vGv-approximation of the susceptibility is calculated. By default, the pGv result is applied for the $\mathbf{G=0}$ contribution to the shifts. With LVGVAPPL one can force VASP to use the vGv result for the $\mathbf{G=0}$ contribution instead. With LVGVCALC one can suppress calculation of the vGv susceptibility. For details see LVGVCALC.

Advice and recommendations

Input parameters

• A larger ENCUT value than usual, generally much higher than the value given by ENMAX in the POTCAR file, e.g. 800 eV for C.
• A small EDIFF is typically required to provide converged chemical shifts, e.g. 1E-8 eV.
• Tighter precision, e.g. PREC = Accurate.

Two additional terms may make a difference depending on your system, which should be tested with and without to determine their importance:

• Non-spherical contributions to the gradient of the density inside PAW spheres, i.e. LASPH = .TRUE.
• Occasionally, e.g. for systems containing H or first-row elements, and short bonds, the two-center contributions are important. In this case, LLRAUG = .TRUE. should be used .

Important: The treatment of the orbital magnetism is non-relativistic. This is fine for light nuclei. The standard POTCARs are scalar-relativistic and account partially for relativistic effects. The accuracy can be improved using LBONE, which restores the small B-component of the wave function inside the PAW spheres. Spin-orbit coupling is not implemented for chemical shift calculations.

PAW pseudopotentials

No special POTCAR files are necessary. The GIPAW is applied using the projectors functions and partial waves that are stored in the regular POTCAR files. A few remarks, however, on accuracy in relation to the different POTCAR flavours:

• Results sensitively depend on the quality, i.e., completeness of the partial wave/projector function set in the energy range needed for good chemical transferability. Result obtained with different POTCAR flavours can differ a few ppm for first and second row sp-bonded elements (except for H).
• Use POTCAR files generated with a consistent exchange-correlation functional. The PAW reconstruction with AE partial waves is crucial as the field on the nucleus needs to be calculated. So avoid, if possible, overriding LEXCH from POTCAR with an explicit GGA-tag in the INCAR.

Insufficient memory

In the case of insufficient memory, there are a few options. Speed has been favored over saving memory. Since the linear response calculation is parallel over k-points, this can be used to economize on memory:

• First, do a regular self-consistent calculation at high accuracy for the full k-point mesh. Save the CHGCAR file.
• Second, do a chemical shift calculation for each k-point in the IBZ separately, starting from CHGCAR, i.e., using ICHARG=11.
• Finally, calculate the shifts as a k-point weighted average of the symmetrized shifts of the individual k-points.

Related tags and articles

DQ, ICHIBARE, LNMR_SYM_RED, NLSPLINE, LLRAUG, LBONE, LVGVCALC, LVGVAPPL

Examples that use this tag