MSDGW F: Difference between revisions

From VASP Wiki
(Created page with "{{TAGDEF|MSDGW_F|[real]|-1}} Description: A positive value of {{TAG|MSDGW_F}} triggers the mixed stochastic-deterministic compression algorithm of Altman and co-workers.{{cite|altman:prl:2024}} ---- {{TAG|MSDGW_F}} is the constant energy ratio <math>F</math> of the compression algorithm. If set to a positive value, energies beyond the protected space defined by {{TAG|MSDGW_NP}} are subdivided into energy bins of width <math>\Delta E_i</math> and replaced by other energi...")
 
No edit summary
Line 13: Line 13:
   bands after compression:      240   
   bands after compression:      240   


To test different compression settings, it is recommended perform the GW/CRPA calculation in steps and to set {{TAG|MSDGW_F}} only in the actual GW step. This avoids repeating the expensive exact diagonalization of the Kohn-Sham hamiltonian.
To test different compression settings, it is recommended perform the GW/CRPA calculation in steps and to set {{TAG|MSDGW_F}} only in the actual GW step. This avoids repeating the expensive exact diagonalization of the Kohn-Sham hamiltonian.  


==Caveats==
==Caveats==
If the GW/CRPA calculation is performed in steps in combination with {{FILE|WAVEDER}}, it is recommended to set {{TAG|LPEAD}}=T. Otherwise the long-wave limit is not calculated and the k-point convergence of band-gaps is deteriorated.  
If the GW/CRPA calculation is performed in steps in combination with {{FILE|WAVEDER}}, {{TAG|LPEAD}} must be set in addition to {{TAG|LOPTICS}} to recalculated the long-wave limit.


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

Revision as of 11:39, 18 February 2025

MSDGW_F = [real]
Default: MSDGW_F = -1 

Description: A positive value of MSDGW_F triggers the mixed stochastic-deterministic compression algorithm of Altman and co-workers.[1]


MSDGW_F is the constant energy ratio of the compression algorithm. If set to a positive value, energies beyond the protected space defined by MSDGW_NP are subdivided into energy bins of width and replaced by other energies , such that . The original orbitals are replaced by MSDGW_NXI randomly linear combined orbitals. Larger values of MSDGW_F increase the compression level at expense of accuracy. The same holds true for smaller values of MSDGW_NXI.

This compression algorithm has been developed to reduce the large number of virtual states required for the calculation of the screened interaction in GW calculations. It has been demonstrated[1] thay one can reduce the virtual manifold by more than 50 per cent and speed up the GW step by a factor of 2 or more with an resulting error of only 50 meV or less on the quasi particle band gap. Comparable results can be achieved for Constrained Random Phase Approximation (CRPA) calculations.

Mind: Avaliable as of VASP version > 6.5.0.

Use cases

  • Recommended for ALGO=EVG0W[R|RK]|CRPA[R|RK]:

The compression can be used for all type of GW calculations regardless if the calculation is performed in steps or in the all-in-one mode. Following lines in stdout and OUTCAR indicate that the compression has been performed:

 => avg. energy ratio F (%):  1.00                                                                                                                                                                                                                                                                                          
 bands after compression:      240  

To test different compression settings, it is recommended perform the GW/CRPA calculation in steps and to set MSDGW_F only in the actual GW step. This avoids repeating the expensive exact diagonalization of the Kohn-Sham hamiltonian.

Caveats

If the GW/CRPA calculation is performed in steps in combination with WAVEDER, LPEAD must be set in addition to LOPTICS to recalculated the long-wave limit.

Related tags and articles

MSDGW_NXI, MSDGW_SEED, MSDGW_NP

References