LEFG: Difference between revisions
No edit summary |
No edit summary |
||
Line 19: | Line 19: | ||
To convert the ''V''<sub>zz</sub> values into the ''C''<sub>q</sub> often encountered in NMR literature, one has to specify the nuclear quadrupole moment by means of the {{TAG|QUAD_EFG}}-tag. | To convert the ''V''<sub>zz</sub> values into the ''C''<sub>q</sub> often encountered in NMR literature, one has to specify the nuclear quadrupole moment by means of the {{TAG|QUAD_EFG}}-tag. | ||
==Recommendations and advice== | |||
Tight settings are required for calculating the electric field gradient. | |||
===Input parameters=== | |||
* A larger {{TAG|ENCUT}} value than usual, generally much higher than the value given by ENMAX in the {{FILE|POTCAR}} file, e.g. 800 eV for C in diamond, rather than the standard 400 eV. | |||
* A small {{TAG|EDIFF}} is required to provide converged chemical shifts, e.g. <code>1E-8</code> eV. | |||
* Tighter precision, e.g. {{TAG|PREC}} = Accurate. | |||
* Non-spherical contributions to the gradient of the density inside PAW spheres, i.e. {{TAG|LASPH}} = .TRUE. | |||
===Structure=== | |||
* The structure is extremely important, so using the experimental structure can improve results. | |||
===PAW pseudopotentials=== | |||
* The use of PAW potentials has a strong influence, GW {{FILE|POTCAR}} files often improve values. | |||
* Semi-core electrons can be important (check the {{TAG|POSCAR}} files with ''*_pv'' or ''*_sv'') as well as explicit inclusion of augmentation channels with <math>d</math>-projectors. | |||
{{NB|important|Several definitions of <math>C_q</math> are used in the NMR community, ensure that you are comparing between the same definitions in calculation and experiment.}} | |||
{{NB|important|For heavy nuclei inaccuracies are to be expected because of an incomplete treatment of relativistic effects.}} | |||
== Related tags and articles == | == Related tags and articles == |
Revision as of 16:15, 27 February 2025
LEFG = .TRUE. | .FALSE.
Default: LEFG = .FALSE.
Description: The LEFG computes the Electric Field Gradient at positions of the atomic nuclei.
For LEFG=.TRUE., the electric field gradient tensors at the positions of the atomic nuclei are calculated using the method of Petrilli et al. [1].
The EFG tensors are symmetric. The principal components Vii and asymmetry parameter η are printed for each atom. Following convention the principal components Vii are ordered such that:
The asymmetry parameter is defined as . For so-called "quadrupolar nuclei", i.e., nuclei with nuclear spin I>1/2, NMR experiments can access Vzz and η.
Mind: Attaining convergence can require somewhat smaller EDIFF than the default of 1.e-4
and somewhat larger cutoff ENCUT than default with PREC=A. Moreover, the calculation of EFGs typically requires high quality PAW data sets. Semi-core electrons can be important (check with *_pv or *_sv POTCARs) as well as explicit inclusion of augmentation channel(s) with d-projectors. |
To convert the Vzz values into the Cq often encountered in NMR literature, one has to specify the nuclear quadrupole moment by means of the QUAD_EFG-tag.
Recommendations and advice
Tight settings are required for calculating the electric field gradient.
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 in diamond, rather than the standard 400 eV.
- A small EDIFF is required to provide converged chemical shifts, e.g.
1E-8
eV. - Tighter precision, e.g. PREC = Accurate.
- Non-spherical contributions to the gradient of the density inside PAW spheres, i.e. LASPH = .TRUE.
Structure
- The structure is extremely important, so using the experimental structure can improve results.
PAW pseudopotentials
- The use of PAW potentials has a strong influence, GW POTCAR files often improve values.
- Semi-core electrons can be important (check the POSCAR files with *_pv or *_sv) as well as explicit inclusion of augmentation channels with -projectors.
Important: Several definitions of are used in the NMR community, ensure that you are comparing between the same definitions in calculation and experiment. |
Important: For heavy nuclei inaccuracies are to be expected because of an incomplete treatment of relativistic effects. |