LHYPERFINE: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
Line 3: Line 3:
Description: compute the hyperfine tensors at the atomic sites (available as of vasp.5.3.2).
Description: compute the hyperfine tensors at the atomic sites (available as of vasp.5.3.2).
----
----
To have VASP compute the hyperfine tensors at the atomic sites, set


LHYPERFINE = .TRUE.
The hyperfine tensor A<sup>I</sup> describes the interaction between a nuclear spin S<sup>I</sup> (located at site '''R'''<sub>I</sub>) and the electronic spin distribution S<sup>e</sup> (in most cases associated with a paramagnetic defect state):
:<math>
E=\sum_{ij} S^e_i A^I_{ij} S^I_j
</math>
In general it is written as the sum of an isotropic part, the so-called Fermi contact term, and an anisotropic (dipolar) part.
The Fermi contact term is given by
:<math>
(A^I_{\mathrm{iso}})_{ij}= \frac{2}{3}\frac{\mu_0\gamma_e\gamma_I}{<S_z>}\delta_{ij}\int \delta_T(\mathbf{r})\rho_s(\mathbf{r}+\mathbf{R}_I)d\mathbf{r}
</math>
where &rho;<sub>s</sub> is the spin density, &mu;<sub>0</sub> is the magnetic susceptibility of free space,
&gamma;<sub>e</sub> the electron gyromagnetic ratio, &gamma;<sub>I</sub> the nuclear gyromagnetic ratio of the nucleus at '''R'''<sub>I</sub>, and \<S<sub>z</sub>\> the expectation value of the ''z''-component of the total electronic spin.
&delta;<sub>T</sub>('''r''') is a smeared out &delta; function, as described in
the Appendix of Ref.<ref name="bloechl:prb:00"/>.
The dipolar contributions to the hyperfine tensor are given by
:<math>
(A^I_{\mathrm{ani}})_{ij}=\frac{\mu_0}{4\pi}\frac{\gamma_e\gamma_I}{<S_z>}
\int \frac{\rho_s(\mathbf{r}+\mathbf{R}_I)}{r^3}\frac{3r_ir_j-\delta_{ij}r^2}{r^2} d\mathbf{r}
</math>
In the equations above ''r''=|'''r'''|, ''r''<sub>i</sub> the i-th component of '''r''', and '''r''' is
taken relative to the position of the nucleus '''R'''<sub>I</sub>.
The nuclear gyromagnetic ratios should be specified by means of the {{TAG|NGYROMAG}}-tag:
NGYROMAG = gamma_1  gamma_2 ... gamma_N
where one should specify one number for each of the ''N'' species on the {{FILE|POSCAR}} file.
If one does not set {{TAG|NGYROMAG}} in the {{FILE|INCAR}} file, VASP assumes a factor of 1 for each species.
As usual, all output is written to the {{FILE|OUTCAR}} file. VASP writes three blocks of data, that look
something like:


== Related Tags and Sections ==
== Related Tags and Sections ==
{{TAG|NGYROMAG}}
{{TAG|NGYROMAG}}


== References ==
<references>
<ref name="bloechl:prb:00">[http://link.aps.org/doi/10.1103/PhysRevB.62.6158 P. E. Blöchl, Phys. Rev. B 62, 6158 (2000).]</ref>
</references>
----
----
[[The_VASP_Manual|Contents]]
[[The_VASP_Manual|Contents]]


[[Category:INCAR]]
[[Category:INCAR]]

Revision as of 18:40, 3 December 2012

LHYPERFINE = .TRUE. | .FALSE.
Default: LHYPERFINE = .FALSE. 

Description: compute the hyperfine tensors at the atomic sites (available as of vasp.5.3.2).


To have VASP compute the hyperfine tensors at the atomic sites, set

LHYPERFINE = .TRUE.

The hyperfine tensor AI describes the interaction between a nuclear spin SI (located at site RI) and the electronic spin distribution Se (in most cases associated with a paramagnetic defect state):

[math]\displaystyle{ E=\sum_{ij} S^e_i A^I_{ij} S^I_j }[/math]

In general it is written as the sum of an isotropic part, the so-called Fermi contact term, and an anisotropic (dipolar) part.

The Fermi contact term is given by

[math]\displaystyle{ (A^I_{\mathrm{iso}})_{ij}= \frac{2}{3}\frac{\mu_0\gamma_e\gamma_I}{\lt S_z\gt }\delta_{ij}\int \delta_T(\mathbf{r})\rho_s(\mathbf{r}+\mathbf{R}_I)d\mathbf{r} }[/math]

where ρs is the spin density, μ0 is the magnetic susceptibility of free space, γe the electron gyromagnetic ratio, γI the nuclear gyromagnetic ratio of the nucleus at RI, and \<Sz\> the expectation value of the z-component of the total electronic spin.

δT(r) is a smeared out δ function, as described in the Appendix of Ref.[1].

The dipolar contributions to the hyperfine tensor are given by

[math]\displaystyle{ (A^I_{\mathrm{ani}})_{ij}=\frac{\mu_0}{4\pi}\frac{\gamma_e\gamma_I}{\lt S_z\gt } \int \frac{\rho_s(\mathbf{r}+\mathbf{R}_I)}{r^3}\frac{3r_ir_j-\delta_{ij}r^2}{r^2} d\mathbf{r} }[/math]

In the equations above r=|r|, ri the i-th component of r, and r is taken relative to the position of the nucleus RI.

The nuclear gyromagnetic ratios should be specified by means of the NGYROMAG-tag:

NGYROMAG = gamma_1  gamma_2 ... gamma_N

where one should specify one number for each of the N species on the POSCAR file. If one does not set NGYROMAG in the INCAR file, VASP assumes a factor of 1 for each species.

As usual, all output is written to the OUTCAR file. VASP writes three blocks of data, that look something like:

Related Tags and Sections

NGYROMAG

References


Contents