Electric Field Gradient: Difference between revisions

From VASP Wiki
No edit summary
m (Pyykko quadrupole moments update)
 
(9 intermediate revisions by 3 users not shown)
Line 2: Line 2:


The following flags control the behaviour of VASP (the given value is the default value):
The following flags control the behaviour of VASP (the given value is the default value):
*{{TAG|LEFG}}=.FALSE. This tag switches on the calculation of the electric field gradient tensors ({{TAG|LEFG}}=''.TRUE.''). The EFG tensors are symmetric. The principal components <math>V_{ii}</math> and asymmetry parameter <math>\eta</math> are printed out for each atom. Following convention is made for the principal components <math>V_{ii}</math>:
*{{TAG|LEFG}}=.TRUE. This tag switches on the calculation of the electric field gradient tensors. The EFG tensors are symmetric. The principal components <math>V_{ii}</math> and asymmetry parameter <math>\eta</math> are printed out for each atom. Following convention is made for the principal components <math>V_{ii}</math>:


<math>|V_{zz}| > |V_{xx}| > |V_{yy}|</math>
<math>|V_{zz}| > |V_{xx}| > |V_{yy}|</math>
Line 14: Line 14:
<math>C_q = \frac{e Q V_{zz}}{h}</math>
<math>C_q = \frac{e Q V_{zz}}{h}</math>


The {{TAG|QUAD_EFG}} tag consists of the nuclear quadrupole moment in millibarns for each atomic species, in the same order as in the {{TAG|POTCAR}} file. The output of <math>C_q</math> is in MHz. See reference <ref name="pyykko2008"/> for a compilation of nuclear quadrupole moments.
The {{TAG|QUAD_EFG}} tag consists of the nuclear quadrupole moment in millibarns for each atomic species, in the same order as in the {{TAG|POTCAR}} file. The output of <math>C_q</math> is in MHz. See references <ref name="pyykko2008"/> and <ref name="pyykko2017"/> for a compilation of nuclear quadrupole moments.


Suppose a solid contains Al, C and Si, than the {{TAG|QUAD_EFG}} tag could read:
Suppose a solid contains Al, C and Si, than the {{TAG|QUAD_EFG}} tag could read:


  QUAD_EFG = 146.6 33.27 0
  QUAD_EFG = 146.6 33.27 0
}
<math>^{27}</math>Al is the stable isotope of Al with a natural abundance of 100~\% and $Q = 146.6$.
The stable isotopes $^{12}$C and $^{13}$C are not quadrupolar nuclei, however, the radioactive
$^{11}$C is. It has $Q = 33.27$. For Si it is pointless to calculate a $C_q$: Again all stable isotopes have $I \le 1/2$. No moments are known
for the other isotopes.


Beware: several definitions of $C_q$ are used in the NMR community.
<math>^{27}\mathrm{Al}</math> is the stable isotope of Al with a natural abundance of 100% and <math>Q = 146.6</math>. The stable isotopes <math>^{12}\mathrm{C}</math> and <math>^{13}\mathrm{C}</math> are not quadrupolar nuclei, however, the radioactive <math>^{11}\mathrm{C}</math> is. It has <math>Q = 33.27</math>. For Si it is pointless to calculate a <math>C_q</math> since all stable isotopes have <math>I \le 1/2</math>. No moments are known for the other isotopes.
\end{itemize}


Beware: for heavy nuclei inaccuracies are to be expected because of an incomplete treatement of relativistic effects.
Beware:  
*Several definitions of <math>C_q</math> are used in the NMR community.
*For heavy nuclei inaccuracies are to be expected because of an incomplete treatment of relativistic effects.
 
== Related Tags and Sections ==
{{TAG|LEFG}}
 
{{sc|Electric Field Gradient|Examples|Examples that use this tag}}


== References ==
== References ==
Line 34: Line 35:
<references>
<references>
<ref name="petrilli1998">[http://journals.aps.org/prb/abstract/10.1103/PhysRevB.57.14690 H. M. Petrilli, P. E. Blöchl, P. Blaha, and K. Schwarz, Phys. Rev. B 57, 14690 (1998).]</ref>
<ref name="petrilli1998">[http://journals.aps.org/prb/abstract/10.1103/PhysRevB.57.14690 H. M. Petrilli, P. E. Blöchl, P. Blaha, and K. Schwarz, Phys. Rev. B 57, 14690 (1998).]</ref>
<ref name="pyykko2008">[http://www.tandfonline.com/doi/abs/10.1080/00268970802018367 P. Pyykkö, Mol. Phys. 106, 1965-1974 (2008), online compilation: http://www.chem.helsinki.fi/pyykko/Q2008.pdf.]</ref>
<ref name="pyykko2008">[http://www.tandfonline.com/doi/abs/10.1080/00268970802018367 P. Pyykkö, Mol. Phys. 106, 1965-1974 (2008)]</ref>
<ref name="pyykko2017">[https://www.tandfonline.com/doi/full/10.1080/00268976.2018.1426131 P. Pyykkö, Mol. Phys. 116,  1328-1338 (2018)]</ref>
</references>
</references>


----
----


[[The_VASP_Manual|Contents]]
[[Category:NMR]][[Category:Electric-field gradient]]
 
[[Category:INCAR]]

Latest revision as of 16:11, 25 November 2022

Electric field gradients at the positions of the atomic nuclei can be calculated by VASP using the method described in reference [1].

The following flags control the behaviour of VASP (the given value is the default value):

  • LEFG=.TRUE. This tag switches on the calculation of the electric field gradient tensors. The EFG tensors are symmetric. The principal components and asymmetry parameter are printed out for each atom. Following convention is made for the principal components :

The asymmetry parameter . For so-called "quadrupolar nuclei", i.e. nuclei with nuclear spin , NMR experiments can access and .

Beware: Attaining convergence can require somewhat smaller EDIFF parameters than the default of EDIFF=1.e-4 and a somewhat larger cutoff ENCUT than the default one with PREC=A. Moreover, the calculation of EFGs typically requires high quality PAW data sets. Semi-core electrons can be important (check the POSCAR files with *_pv or *_sv) as well as explicit inclusion of augmentation channels with -projectors.

  • QUAD_EFG=1 This tag allows the conversion by VASP of the values into the often encountered in NMR literature. The conversion formula is given as follows( is the element and isotope specific quadrupole moment):

The QUAD_EFG tag consists of the nuclear quadrupole moment in millibarns for each atomic species, in the same order as in the POTCAR file. The output of is in MHz. See references [2] and [3] for a compilation of nuclear quadrupole moments.

Suppose a solid contains Al, C and Si, than the QUAD_EFG tag could read:

QUAD_EFG = 146.6 33.27 0

is the stable isotope of Al with a natural abundance of 100% and . The stable isotopes and are not quadrupolar nuclei, however, the radioactive is. It has . For Si it is pointless to calculate a since all stable isotopes have . No moments are known for the other isotopes.

Beware:

  • Several definitions of are used in the NMR community.
  • For heavy nuclei inaccuracies are to be expected because of an incomplete treatment of relativistic effects.

Related Tags and Sections

LEFG

Examples that use this tag

References