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):

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

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

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