PROCAR: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
Line 24: Line 24:
For {{TAG|ISPIN}}=2 PROCAR contains a second set of projections for the spin down channel.  
For {{TAG|ISPIN}}=2 PROCAR contains a second set of projections for the spin down channel.  


For {{TAG|LNONCOLLINEAR}}=.TRUE. three additional projections for each ion are printed and the output is similar to
ion      s    py    pz    px    dxy    dyz    dz2    dxz  x2-y2    tot 
    1  0.144  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.145
    2  0.291  0.000  0.006  0.000  0.000  0.000  0.000  0.000  0.000  0.298
    3  0.291  0.000  0.006  0.000  0.000  0.000  0.000  0.000  0.000  0.298
tot    0.727  0.000  0.013  0.000  0.000  0.000  0.000  0.000  0.000  0.740
    1 -0.011 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000  0.000 -0.000 -0.011
    2 -0.023 -0.000 -0.000  0.000  0.000 -0.000 -0.000  0.000 -0.000 -0.023
    3 -0.023 -0.000 -0.000  0.000  0.000 -0.000 -0.000  0.000 -0.000 -0.023
tot  -0.057 -0.000 -0.001  0.000  0.000 -0.000 -0.000  0.000 -0.000 -0.058
    1 -0.142 -0.000  0.000  0.000  0.000  0.000 -0.000 -0.000 -0.000 -0.142
    2 -0.286  0.000 -0.006 -0.000 -0.000  0.000 -0.000 -0.000  0.000 -0.293
    3 -0.286  0.000 -0.006 -0.000 -0.000  0.000 -0.000 -0.000  0.000 -0.293
tot  -0.715  0.000 -0.012 -0.000 -0.000  0.000 -0.000 -0.000  0.000 -0.727
    1 -0.024 -0.000  0.000 -0.000 -0.000  0.000 -0.000  0.000 -0.000 -0.024
    2 -0.048  0.000 -0.001  0.000  0.000  0.000 -0.000  0.000  0.000 -0.049
    3 -0.048  0.000 -0.001  0.000  0.000  0.000 -0.000  0.000  0.000 -0.049
tot  -0.119  0.000 -0.002  0.000  0.000  0.000 -0.000  0.000  0.000 -0.121
Here the entries correspond to the projected magnetizations
<math>1/2\sum_{\mu,\nu=1}^2\sigma_{\mu\nu}^j
\langle \chi_{n {\bf k}}^\mu | Y_{lm}^\alpha \rangle
\langle  Y_{lm}^\alpha | \chi_{n {\bf k}}^\nu \rangle
</math> and are calculated for the spinor
of the spinor <math>|\Psi_{n{\bf k}}\rangle=\left(\begin{matrix}\chi_{n{\bf k}}^\uparrow \\\chi_{n{\bf k}}^\downarrow \end{matrix}\right)</math> and the Pauli matrices:
<math>
\sigma^x = \left(\begin{matrix}
  0  &  1 \\
  1  &  0 \\
\end{matrix}\right), \quad
\sigma^y = \left(\begin{matrix}
  0  &  -i \\
  i  &  0 \\
\end{matrix}\right), \quad
\sigma^z = \left(\begin{matrix}
  1  &  0 \\
  0  &  -1 \\
\end{matrix}\right)
</math>
The first set is the total (absolute) magnetization, while the remaining three sets of entries correspond to the three directions <math>j=1,2,3</math>.





Revision as of 11:07, 13 March 2019

For static calculations, the file PROCAR contains the spd- and site projected wave function character of each band. The wave function character is calculated by projecting the wave functions onto spherical harmonics that are non zero within spheres of a radius RWIGS around each ion. LORBIT determines what information is written to PROCAR. For instance RWIGS must be specified in the INCAR file in order to obtain the file for LORBIT<10.

Mind: The spd- and site projected character of each band is not evaluated in the parallel version if NPARFailed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): \neq 1.

# of k-points:    5         # of bands:   26         # of ions:    3

k-point     1 :    0.00000000 0.00000000 0.00000000     weight = 0.06250000

band     1 # energy  -17.37867948 # occ.  1.00000000

ion      s     py     pz     px    dxy    dyz    dz2    dxz  x2-y2    tot  
    1  0.144  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.145
    2  0.291  0.000  0.006  0.000  0.000  0.000  0.000  0.000  0.000  0.298
    3  0.291  0.000  0.006  0.000  0.000  0.000  0.000  0.000  0.000  0.298
tot    0.727  0.000  0.013  0.000  0.000  0.000  0.000  0.000  0.000  0.740

The header contains the information about the number of k-points, bands and how many ions are considered. The next line prints the k-point with the three coordinates in the first Brillouin zone and the corresponding k-point weight for the numerical integration followed by the band number and the energy and occupancy of the state. Each (k-point,band) pair contains the projections for every ion Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): |\langle Y_{{lm}}^{{\alpha }}|\phi _{{n{\mathbf {k}}}}\rangle |^{2} , where Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): Y_{{lm}}^{\alpha } is the spherical harmonic centered at ion index Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): \alpha , Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): l,m the angular moment and magnetic quantum and Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): \phi _{{n{\mathbf {k}}}} the wavefunction. The line and column with "tot" is the corresponding sum of the line and column, respectively.

For ISPIN=2 PROCAR contains a second set of projections for the spin down channel.

For LNONCOLLINEAR=.TRUE. three additional projections for each ion are printed and the output is similar to

ion      s     py     pz     px    dxy    dyz    dz2    dxz  x2-y2    tot  
    1  0.144  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.000  0.145
    2  0.291  0.000  0.006  0.000  0.000  0.000  0.000  0.000  0.000  0.298
    3  0.291  0.000  0.006  0.000  0.000  0.000  0.000  0.000  0.000  0.298
tot    0.727  0.000  0.013  0.000  0.000  0.000  0.000  0.000  0.000  0.740
    1 -0.011 -0.000 -0.000 -0.000 -0.000 -0.000 -0.000  0.000 -0.000 -0.011
    2 -0.023 -0.000 -0.000  0.000  0.000 -0.000 -0.000  0.000 -0.000 -0.023
    3 -0.023 -0.000 -0.000  0.000  0.000 -0.000 -0.000  0.000 -0.000 -0.023
tot   -0.057 -0.000 -0.001  0.000  0.000 -0.000 -0.000  0.000 -0.000 -0.058 
    1 -0.142 -0.000  0.000  0.000  0.000  0.000 -0.000 -0.000 -0.000 -0.142
    2 -0.286  0.000 -0.006 -0.000 -0.000  0.000 -0.000 -0.000  0.000 -0.293
    3 -0.286  0.000 -0.006 -0.000 -0.000  0.000 -0.000 -0.000  0.000 -0.293
tot   -0.715  0.000 -0.012 -0.000 -0.000  0.000 -0.000 -0.000  0.000 -0.727
    1 -0.024 -0.000  0.000 -0.000 -0.000  0.000 -0.000  0.000 -0.000 -0.024
    2 -0.048  0.000 -0.001  0.000  0.000  0.000 -0.000  0.000  0.000 -0.049
    3 -0.048  0.000 -0.001  0.000  0.000  0.000 -0.000  0.000  0.000 -0.049
tot   -0.119  0.000 -0.002  0.000  0.000  0.000 -0.000  0.000  0.000 -0.121

Here the entries correspond to the projected magnetizations Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): 1/2\sum _{{\mu ,\nu =1}}^{2}\sigma _{{\mu \nu }}^{j}\langle \chi _{{n{{\bf {k}}}}}^{\mu }|Y_{{lm}}^{\alpha }\rangle \langle Y_{{lm}}^{\alpha }|\chi _{{n{{\bf {k}}}}}^{\nu }\rangle and are calculated for the spinor of the spinor Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): |\Psi _{{n{{\bf {k}}}}}\rangle =\left({\begin{matrix}\chi _{{n{{\bf {k}}}}}^{\uparrow }\\\chi _{{n{{\bf {k}}}}}^{\downarrow }\end{matrix}}\right) and the Pauli matrices:

Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): \sigma ^{x}=\left({\begin{matrix}0&1\\1&0\\\end{matrix}}\right),\quad \sigma ^{y}=\left({\begin{matrix}0&-i\\i&0\\\end{matrix}}\right),\quad \sigma ^{z}=\left({\begin{matrix}1&0\\0&-1\\\end{matrix}}\right)

The first set is the total (absolute) magnetization, while the remaining three sets of entries correspond to the three directions Failed to parse (SVG (MathML can be enabled via browser plugin): Invalid response ("Math extension cannot connect to Restbase.") from server "https://www.vasp.at/wiki/restbase/vasp.at/v1/":): j=1,2,3 .