LSELFENERGY: Difference between revisions

From VASP Wiki
No edit summary
No edit summary
Line 5: Line 5:
----
----


If {{TAG|LSELFENERGY}}=''.FALSE.'',  QP shifts are evaluated. This is the default behavior.  
* If {{TAGO|LSELFENERGY|.FALSE.}} (default)quasi-particle (QP) shifts are evaluated.  
QP shifts <psi_nk| G(iteration)W_0 |psi_nk>: iteration 1
for sc-GW calculations column KS-energies equals QP-energies in previous step
and V_xc(KS)=  KS-energies - (<T + V_ion + V_H > + <T+V_H+V_ion>^1  + <V_x>^1)
k-point  1 :      0.0000    0.0000    0.0000
  band No.  KS-energies  QP-energies  sigma(KS)  V_xc(KS)    V^pw_x(r,r')  Z            occupation Imag(sigma)
      1      -7.1627      -8.3040    -14.5626    -12.7276    -21.6682      0.6219      2.0000      1.2037
      2      -2.0901      -3.4347    -15.7660    -14.2799    -21.7439      0.9048      2.0000      0.6914
      3      -2.0901      -3.4347    -15.7660    -14.2799    -21.7439      0.9048      2.0000      0.6914
      4      -2.0901      -3.4347    -15.7660    -14.2799    -21.7439      0.9048      2.0000      0.6914
      5      0.4603      -0.4663    -13.7603    -12.5200    -18.1532      0.7471      2.0000      0.2167
      6      0.4603      -0.4663    -13.7603    -12.5200    -18.1532      0.7471      2.0000      0.2167


If {{TAG|LSELFENERGY}}=''.TRUE.'', the frequency dependent self-energy <math> \langle \psi_{n {\mathbf{k}}} | \Sigma(\omega) |\psi_{n {\mathbf{k}}} \rangle </math> is evaluated and printed to {{FILE|vasprun.xml}}. An example output looks like the following:
* If {{TAGO|LSELFENERGY|.TRUE.}}, the frequency dependent self-energy <math> \langle \psi_{n {\mathbf{k}}} | \Sigma(\omega) |\psi_{n {\mathbf{k}}} \rangle </math> is evaluated and printed to {{FILE|vasprun.xml}}. An example output looks like the following:
  <varray name="selfenergy" >                                                                                                                                                                                                                                               
  <varray name="selfenergy" >                                                                                                                                                                                                                                               
   <v>    -150.00000000    -25.40060536      0.24429448 </v>                                                                                                                                                                                                               
   <v>    -150.00000000    -25.40060536      0.24429448 </v>                                                                                                                                                                                                               
Line 15: Line 28:
   <v>    -148.80000000    -25.42221682      0.25406890 </v>                                                                                                                                                                                                               
   <v>    -148.80000000    -25.42221682      0.25406890 </v>                                                                                                                                                                                                               
   <v>    -148.50000000    -25.42762671      0.25647992 </v>                                                                                                                                                                                                               
   <v>    -148.50000000    -25.42762671      0.25647992 </v>                                                                                                                                                                                                               
   <v>    -148.20000000    -25.43303731      0.25888834 </v>      
   <v>    -148.20000000    -25.43303731      0.25888834 </v>
To print the self-energy is a slight extra computational effort since, within the GW algorithms, the self-energy is usually just evaluated near KS eigenenergies <math>\epsilon_{nk}</math> and not the entire frequency range.     
{{NB|mind|In quartic-scaling GW the self-energy is given on the real-frequency axis, while for low-scaling GW the self-energy is given on the imaginary-frequency axis.}}


If quartic-scaling GW algorithms are selected, e.g. {{TAG|ALGO}}=EVGW0, the first column corresponds to points on the real-frequency axis. The second and third column are the real and imaginary parts of the self-energy at a given k-point and band index. In this case the evaluation of QP shifts is also bypassed and a small set of self-energy points are printed to {{FILE|OUTCAR}} instead, similar to the following output
== Format ==
 
=== Real frequencies ===
 
If '''quartic-scaling GW algorithms''' are selected, e.g. {{TAG|ALGO}}=EVGW0, the first column corresponds to points on the real-frequency axis (in eV). The second and third columns are the real and imaginary parts of the self-energy (in eV) at a given band index and '''k''' point. To identify the band index and '''k''' point, the ordering has to be taken from the {{FILE|OUTCAR}}: Instead of the QP shifts, a small set of self-energy points are printed to {{FILE|OUTCAR}}, similar to the following output
  calculating selfenergy CALC_SELFENERGY_LINEAR between w=-150.00 150.00                                                                                                                                                                                                     
  calculating selfenergy CALC_SELFENERGY_LINEAR between w=-150.00 150.00                                                                                                                                                                                                     
                                                                                                                                                                                                                                                                              
                                                                                                                                                                                                                                                                              
  k-point  1 :      0.0000    0.0000    0.0000                                                                                                                                                                                                                            
  k-point  1 :      0.0000    0.0000    0.0000                                                                                                                                                                                    
   band No.  band energies    occupation                                                                                                                                                                                                                                  
   band No.  band energies    occupation                                                                                                                                                                                          
                                                 
                                                                                                                                                                                                                                   
       1    -11.7125     2.00000  selfenergy along real axis                                                                                                                                                                                                              
       1    -11.4323     2.00000  selfenergy along real axis                                                                                                                                                                      
   -150.0000000  -25.4006054     0.2442945                                                                                                                                                                                                                                 
   -150.0000000  -24.0756124     0.2065066                                                                                                                                                                                         
   -147.0000000  -25.4546849     0.2685030                                                                                                                                                                                                                                 
   -147.0000000  -24.1277845     0.2302741                                                                                                                                                                                         
   -144.0000000  -25.5088867     0.2923888      
   -144.0000000  -24.1803224     0.2537669                                                                                                                                                                                         
...
  ...
  147.0000000  -20.3498375    -2.8348252                                                                                                                                                                                         
  150.0000000  -20.2310127    -2.7491028                                                                                                                                                                                         
                                                                                                                                                                                                                                   
      2      -2.7832      2.00000  selfenergy along real axis                                                                                                                                                                     
  -150.0000000  -13.0060959    0.1938781                                                                                                                                                                                         
  -147.0000000  -13.0530569     0.2231126                                                                                                                                                                                         
  -144.0000000  -13.1030584    0.2520593
  ...
 
Here, the first and second band at the Gamma point are printed.
The line with <code>selfenergy along real axis</code> contains band No., the KS energy, and the occupation of this state.
The output in {{FILE|vasprun.xml}} has the same ordering of bands and '''k''' points, i.e., the band index is always the fastest. The frequency grid cannot be controlled. It is always 1000 points in the range of -150 to 150.
 
=== Imaginary frequencies ===


For [[Practical_guide_to_GW_calculations#Low_scaling_GW_algorithms|low-scaling GW algorithms]], the QP-energies are evaluated and printed to {{FILE|OUTCAR}}. In contrast, the {{FILE|vasprun.xml}} file now contains the self-energy on the imaginary frequency axis for the points selected by the Minimax routines.{{cite|Kaltak:PRB:2020}}
For [[Practical_guide_to_GW_calculations#Low_scaling_GW_algorithms|low-scaling GW algorithms]], the QP shifts are evaluated and printed to {{FILE|OUTCAR}} for both, {{TAG|LSELFENERGY}}=T and F. For {{TAG|LSELFENERGY}}=T, the {{FILE|vasprun.xml}} file additionally contains the self-energy for a given band index and '''k''' point on the imaginary-frequency axis:


  <varray name="selfenergy along imaginary axis" >                                                                                                                                                                                                                           
  <varray name="selfenergy along imaginary axis" >                                                                                                                                                                                                                           
Line 37: Line 70:
   <v>      1.3058951128526406    -23.7949927690568188    -0.1482973553314607</v>         
   <v>      1.3058951128526406    -23.7949927690568188    -0.1482973553314607</v>         
  ...
  ...
Here the number of evaluated frequency points depends on {{TAG|NOMEGA}}.
To identify the band index and '''k''' point, the ordering has to be taken from the QP shifts block in the {{FILE|OUTCAR}}. The band index is faster than the '''k''' point.
The imaginary frequency is selected by the Minimax routines{{cite|Kaltak:PRB:2020}}, and the number of points depends on {{TAG|NOMEGA}}.


== Related tags and articles ==
== Related tags and articles ==

Revision as of 10:14, 20 February 2024

LSELFENERGY = [logical]
Default: LSELFENERGY = .FALSE. 

Description: Controls whether the frequency-dependent self-energy is calculated or not.


  • If LSELFENERGY = .FALSE. (default), quasi-particle (QP) shifts are evaluated.
QP shifts <psi_nk| G(iteration)W_0 |psi_nk>: iteration 1
for sc-GW calculations column KS-energies equals QP-energies in previous step 
and V_xc(KS)=  KS-energies - (<T + V_ion + V_H > + <T+V_H+V_ion>^1  + <V_x>^1)

k-point   1 :       0.0000    0.0000    0.0000
 band No.  KS-energies  QP-energies   sigma(KS)   V_xc(KS)     V^pw_x(r,r')   Z            occupation Imag(sigma)

     1      -7.1627      -8.3040     -14.5626     -12.7276     -21.6682       0.6219       2.0000       1.2037
     2      -2.0901      -3.4347     -15.7660     -14.2799     -21.7439       0.9048       2.0000       0.6914
     3      -2.0901      -3.4347     -15.7660     -14.2799     -21.7439       0.9048       2.0000       0.6914
     4      -2.0901      -3.4347     -15.7660     -14.2799     -21.7439       0.9048       2.0000       0.6914
     5       0.4603      -0.4663     -13.7603     -12.5200     -18.1532       0.7471       2.0000       0.2167
     6       0.4603      -0.4663     -13.7603     -12.5200     -18.1532       0.7471       2.0000       0.2167
  • If LSELFENERGY = .TRUE., the frequency dependent self-energy is evaluated and printed to vasprun.xml. An example output looks like the following:
<varray name="selfenergy" >                                                                                                                                                                                                                                               
 <v>    -150.00000000     -25.40060536       0.24429448 </v>                                                                                                                                                                                                              
 <v>    -149.70000000     -25.40600800       0.24673910 </v>                                                                                                                                                                                                              
 <v>    -149.40000000     -25.41141065       0.24918372 </v>                                                                                                                                                                                                              
 <v>    -149.10000000     -25.41681330       0.25162834 </v>                                                                                                                                                                                                              
 <v>    -148.80000000     -25.42221682       0.25406890 </v>                                                                                                                                                                                                              
 <v>    -148.50000000     -25.42762671       0.25647992 </v>                                                                                                                                                                                                              
 <v>    -148.20000000     -25.43303731       0.25888834 </v>  

To print the self-energy is a slight extra computational effort since, within the GW algorithms, the self-energy is usually just evaluated near KS eigenenergies and not the entire frequency range.

Mind: In quartic-scaling GW the self-energy is given on the real-frequency axis, while for low-scaling GW the self-energy is given on the imaginary-frequency axis.

Format

Real frequencies

If quartic-scaling GW algorithms are selected, e.g. ALGO=EVGW0, the first column corresponds to points on the real-frequency axis (in eV). The second and third columns are the real and imaginary parts of the self-energy (in eV) at a given band index and k point. To identify the band index and k point, the ordering has to be taken from the OUTCAR: Instead of the QP shifts, a small set of self-energy points are printed to OUTCAR, similar to the following output

calculating selfenergy CALC_SELFENERGY_LINEAR between w=-150.00 150.00                                                                                                                                                                                                     
                                                                                                                                                                                                                                                                           
k-point   1 :       0.0000    0.0000    0.0000                                                                                                                                                                                     
 band No.  band energies     occupation                                                                                                                                                                                            
                                                                                                                                                                                                                                   
     1     -11.4323      2.00000  selfenergy along real axis                                                                                                                                                                       
 -150.0000000   -24.0756124     0.2065066                                                                                                                                                                                          
 -147.0000000   -24.1277845     0.2302741                                                                                                                                                                                          
 -144.0000000   -24.1803224     0.2537669                                                                                                                                                                                          
 ...
  147.0000000   -20.3498375    -2.8348252                                                                                                                                                                                          
  150.0000000   -20.2310127    -2.7491028                                                                                                                                                                                          
                                                                                                                                                                                                                                   
     2      -2.7832      2.00000  selfenergy along real axis                                                                                                                                                                       
 -150.0000000   -13.0060959     0.1938781                                                                                                                                                                                          
 -147.0000000   -13.0530569     0.2231126                                                                                                                                                                                          
 -144.0000000   -13.1030584     0.2520593 
 ...

Here, the first and second band at the Gamma point are printed. The line with selfenergy along real axis contains band No., the KS energy, and the occupation of this state. The output in vasprun.xml has the same ordering of bands and k points, i.e., the band index is always the fastest. The frequency grid cannot be controlled. It is always 1000 points in the range of -150 to 150.

Imaginary frequencies

For low-scaling GW algorithms, the QP shifts are evaluated and printed to OUTCAR for both, LSELFENERGY=T and F. For LSELFENERGY=T, the vasprun.xml file additionally contains the self-energy for a given band index and k point on the imaginary-frequency axis:

<varray name="selfenergy along imaginary axis" >                                                                                                                                                                                                                          
  <v>       0.1570801806644298    -23.7971327349915711     -0.0179645094529067</v>                                                                                                                                                                                         
  <v>       0.4718033117773284    -23.7968260788831572     -0.0537856829291801</v>                                                                                                                                                                                         
  <v>       0.8108804505648752    -23.7966749072893577     -0.0932103588740998</v>                                                                                                                                                                                         
  <v>       1.3058951128526406    -23.7949927690568188     -0.1482973553314607</v>         
...

To identify the band index and k point, the ordering has to be taken from the QP shifts block in the OUTCAR. The band index is faster than the k point. The imaginary frequency is selected by the Minimax routines[1], and the number of points depends on NOMEGA.

Related tags and articles

Examples that use this tag