PROOUT: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
 
(2 intermediate revisions by 2 users not shown)
Line 5: Line 5:
\equiv \langle \beta_{lm}^{\alpha}|S|\psi_{n\mathbf{k}}\rangle
\equiv \langle \beta_{lm}^{\alpha}|S|\psi_{n\mathbf{k}}\rangle
= \underbrace{\langle \beta_{lm}^{\alpha}|\psi_{n\mathbf{k}}\rangle}_{P^{\text{SOFT},\alpha}_{lmn\mathbf{k}}} +  
= \underbrace{\langle \beta_{lm}^{\alpha}|\psi_{n\mathbf{k}}\rangle}_{P^{\text{SOFT},\alpha}_{lmn\mathbf{k}}} +  
\underbrace{\sum_{ij} \langle \beta^\alpha_{lm}|p_i\rangle Q_{ij} \langle p_j | \psi_{n\mathbf{k}}\rangle}_{P^{\text{AUG},\alpha}_{lmn\mathbf{k}}}.
\underbrace{\sum_{ij} \langle \beta^\alpha_{lm}|\tilde{p}_i\rangle Q_{ij} \langle \tilde{p}_j | \psi_{n\mathbf{k}}\rangle}_{P^{\text{AUG},\alpha}_{lmn\mathbf{k}}}.
</math>
</math>


Here, the two terms on the right-hand side are called soft and augmentation part, respectively, and <math>S</math> is the overlap matrix,
Here, the two terms on the right-hand side are called soft and augmentation part, respectively, and <math>S</math> is the overlap matrix,


<math>S = 1+\sum_{ij} |p_i\rangle Q_{ij} \langle p_j|.</math>
<math>S = 1+\sum_{ij} |\tilde{p}_i\rangle Q_{ij} \langle \tilde{p}_j|.</math>


The radial and the angular part of <math>\beta^\alpha_{lm}(\mathbf{r})</math> are described by a linear combination of spherical bessel functions <math>\phi_n(r)</math> and by spherical harmonics <math>Y^\alpha_{lm}(\hat{\mathbf{r}})</math>, respectively:  
The radial and the angular part of <math>\beta^\alpha_{lm}(\mathbf{r})</math> are described by a linear combination of spherical bessel functions <math>\phi_n(r)</math> and by spherical harmonics <math>Y^\alpha_{lm}(\hat{\mathbf{r}})</math>, respectively:  
Line 19: Line 19:
</math>
</math>


The spherical bessel functions <math>\phi_n(r)</math> are parametrized to be non-zero within a radius determined by {{TAG|RWIGS}}. In addition, it so happens that <math>|p_i\rangle</math> has a similar form as <math>|\beta^\alpha_{lm}\rangle</math> which simplifies the computations above.
The spherical bessel functions <math>\phi_n(r)</math> are parametrized to be non-zero within a radius determined by {{TAG|RWIGS}}. In addition, it so happens that <math>|\tilde{p}_i\rangle</math> has a similar form as <math>|\beta^\alpha_{lm}\rangle</math> which simplifies the computations above.


The {{TAG|PROOUT}} file contains similar information as the {{TAG|PROCAR}} file but the following differences exist:
The {{TAG|PROOUT}} file contains similar information as the {{TAG|PROCAR}} file but the following differences exist:
Line 38: Line 38:
*below : augmentation part
*below : augmentation part
*last line: real part of <math>P^{\text{AUG},\alpha}_{lmn\mathbf{k}}</math> for every <math>lm</math>-quantum number (inner loop), ion per type, ion type, band and k point (outer loop).
*last line: real part of <math>P^{\text{AUG},\alpha}_{lmn\mathbf{k}}</math> for every <math>lm</math>-quantum number (inner loop), ion per type, ion type, band and k point (outer loop).
{{NB|warning|This file is not correctly written when {{TAG|LNONCOLLINEAR}} {{=}} .TRUE. for versions of VASP <{{=}} 6.2.1 }}
{{NB|warning|For VASP version <{{=}} 6.2.1, {{FILE|PROOUT}} is not correctly written when {{TAG|LNONCOLLINEAR}} {{=}} True. }}


== Related Tags and Sections ==
== Related Tags and Sections ==

Latest revision as of 15:02, 17 August 2022

This file is written if LORBIT=5 and the RWIGS tag is set in the INCAR file. It contains the projection of the Kohn-Sham orbitals onto a localized orbitals basis which can be written as

Here, the two terms on the right-hand side are called soft and augmentation part, respectively, and is the overlap matrix,

The radial and the angular part of are described by a linear combination of spherical bessel functions and by spherical harmonics , respectively:

The spherical bessel functions are parametrized to be non-zero within a radius determined by RWIGS. In addition, it so happens that has a similar form as which simplifies the computations above.

The PROOUT file contains similar information as the PROCAR file but the following differences exist:

  • The PROOUT file writes the real and imaginary parts of , and the real part of the augmentation part .
  • The PROCAR file contains the information of the squared projection, , whereas the PROOUT file contains .
  • The arrangement of the output is very different in both files.

This information makes it possible to construct, e.g., partial DOS projected onto bonding and anti-bonding molecular orbitals and the so-called crystal-orbital-overlap-population (COOP) function. Depending on the application, users might find it more practical to use the information contained in the PROJCAR and LOCPROJ files. These are controlled by the LOCPROJ tag.

For spin-polarized calculations with ISPIN=2 and noncollinear calculations with LNONCOLLINEAR=.TRUE., PROOUT.1 and PROOUT.2 are written, referring to the Kohn-Sham orbitals' two spin components.

Format

  • line 1: PROOUT
  • line 2: Number of kpoints, bands and ions
  • line 3: Twice the number of types followed by the number of ions for each type
  • line 4: The Fermi weights for each k point (inner loop) and band (outer loop)
  • line 5 ...: Real and imaginary part of for every -quantum number (inner loop), band, ion per type, k point and ion type (outer loop).
  • below : augmentation part
  • last line: real part of for every -quantum number (inner loop), ion per type, ion type, band and k point (outer loop).
Warning: For VASP version <= 6.2.1, PROOUT is not correctly written when LNONCOLLINEAR = True.

Related Tags and Sections

LORBIT, LOCPROJ, PROJCAR, PROCAR