LHYPERFINE: Difference between revisions
Vaspmaster (talk | contribs) No edit summary |
Vaspmaster (talk | contribs) 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 ρ<sub>s</sub> is the spin density, μ<sub>0</sub> is the magnetic susceptibility of free space, | |||
γ<sub>e</sub> the electron gyromagnetic ratio, γ<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. | |||
δ<sub>T</sub>('''r''') is a smeared out δ 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: