ENCUTGWSOFT: Difference between revisions

From VASP Wiki
(document the change from squeezed Coulomb kernel to smoothed one for newer VASP versions.)
 
(2 intermediate revisions by one other user not shown)
Line 2: Line 2:
{{DEF|ENCUTGWSOFT|{{TAG|ENCUTGW}}<math>\times 0.8</math>| for {{TAG|ALGO}}{{=}}''ACFDT''|{{TAG|ENCUTGW}}<math>\times 0.8</math>| as of VASP.6.3 | {{TAG|ENCUTGW}}|else}}
{{DEF|ENCUTGWSOFT|{{TAG|ENCUTGW}}<math>\times 0.8</math>| for {{TAG|ALGO}}{{=}}''ACFDT''|{{TAG|ENCUTGW}}<math>\times 0.8</math>| as of VASP.6.3 | {{TAG|ENCUTGW}}|else}}
{{NB|important|For vasp.6.3 and later releases {{TAG|ENCUTGWSOFT}} always defaults to {{TAG|ENCUTGW}}<math>\times 0.8</math>.}}
{{NB|important|For vasp.6.3 and later releases {{TAG|ENCUTGWSOFT}} always defaults to {{TAG|ENCUTGW}}<math>\times 0.8</math>.}}
Descprition: The flag {{TAG|ENCUTGWSOFT}} sets the energy cutoff for response function, such that it allows to truncate the Coulomb kernel slowly between the energy
Descprition: The flag {{TAG|ENCUTGWSOFT}} sets the energy cutoff for the response function, such that it allows to truncate the Coulomb kernel slowly between the energy
specified by {{TAG|ENCUTGWSOFT}} and {{TAG|ENCUTGW}} using a cosine window function.
specified by {{TAG|ENCUTGWSOFT}} and {{TAG|ENCUTGW}} using a cosine window function.
----
----
RPA/ACFDT correlation energies converge very slowly with respect to <math>\mathbf{G}_{\rm max }</math>.
Thus VASP automatically extrapolates to the infinite basis set limit using a linear regression to the equation: {{cite|harl:2008}}{{cite|harl:2010}}{{cite|klimes:2014}}
<math>E_{\mathrm{c}}({\mathbf{G}})=E_{\mathrm{c}}(\infty)+\frac{A}{{\mathbf{G}}^3}</math>.
This  usually leads to much smoother energy-volume curves in {{TAG|ACFDT calculations}} and {{TAG|MP2 calculations}}.
This  usually leads to much smoother energy-volume curves in {{TAG|ACFDT calculations}} and {{TAG|MP2 calculations}}.
The modified Coulomb kernel is in this case:
The modified Coulomb kernel is in this case:
Line 19: Line 24:
\qquad \mbox{for} \quad  \mathrm{ENCUTGWSOFT}=\frac{\hbar^2G_{min}^2}{2m_e}<\frac{\hbar^2 G^2}{2m_e}<\frac{\hbar^2G_{max}^2}{2m_e}=\mathrm{ENCUTGW}</math>  
\qquad \mbox{for} \quad  \mathrm{ENCUTGWSOFT}=\frac{\hbar^2G_{min}^2}{2m_e}<\frac{\hbar^2 G^2}{2m_e}<\frac{\hbar^2G_{max}^2}{2m_e}=\mathrm{ENCUTGW}</math>  


This kernel ''squeezes'' contributions from large wave vectors <math>G>G_{max}</math> into the window given by {{TAGBL|ENCUTGWSOFT}}. Effectively, this extrapolates the RPA correlation energy to the {{TAG|ENCUTGW}} <math>\to \infty</math> limit, assuming that the basis set incompleteness error falls off as <math>1/\mathrm{ENCUTGW}^{3/2}</math>.
This kernel ''squeezes'' contributions from large wave vectors <math>G>G_{max}</math> into the window given by {{TAGBL|ENCUTGWSOFT}}.
For GW type calculations the squeezed Coulomb kernel was the default (when {{TAGBL|ENCUTGWSOFT}} was set in the INCAR file) before version vasp.6.3, but
in newer releases the code always defaults to a smoothed Coulomb kernel (both for GW and RPA type calculations). If one desires to recover the behavior
of vasp.6.2 and older versions, {{TAG|LSCK}}=.TRUE. must be set in the INCAR file for GW type calculations if {{TAGBL|ENCUTGWSOFT}} is set in the INCAR file.
 
 
{{NB|mind|The infinite basis set limit extrapolation for RPA/ACFDT is described in more detail [[ACFDT/RPA_calculations#Basis_set_convergence|here]].}}
== Related tags and articles ==
== Related tags and articles ==
{{TAG|PRECFOCK}},
{{TAG|PRECFOCK}},
Line 25: Line 36:
{{TAG|ENCUTGW}},
{{TAG|ENCUTGW}},
{{TAG|GW calculations}},
{{TAG|GW calculations}},
{{TAG|LSCK}}
{{TAG|LSCK}},
[[ACFDT/RPA_calculations#Basis_set_convergence|RPA/ACFDT basis set convergence]]
[[ACFDT/RPA_calculations#Basis_set_convergence|RPA/ACFDT basis set convergence]],
{{sc|ENCUTGWSOFT|Examples|Examples that use this tag}}
{{sc|ENCUTGWSOFT|Examples|Examples that use this tag}}
----
----


[[Category:INCAR tag]][[Category:Many-body perturbation theory]][[Category:GW]]
[[Category:INCAR tag]][[Category:Many-body perturbation theory]][[Category:GW]]

Latest revision as of 05:42, 29 August 2024

ENCUTGWSOFT = [real] 

Default: ENCUTGWSOFT = ENCUTGW for ALGO=ACFDT
= ENCUTGW as of VASP.6.3
= ENCUTGW else
Important: For vasp.6.3 and later releases ENCUTGWSOFT always defaults to ENCUTGW.

Descprition: The flag ENCUTGWSOFT sets the energy cutoff for the response function, such that it allows to truncate the Coulomb kernel slowly between the energy specified by ENCUTGWSOFT and ENCUTGW using a cosine window function.


RPA/ACFDT correlation energies converge very slowly with respect to . Thus VASP automatically extrapolates to the infinite basis set limit using a linear regression to the equation: [1][2][3]

.

This usually leads to much smoother energy-volume curves in ACFDT calculations and MP2 calculations. The modified Coulomb kernel is in this case:

If LSCK is set to .TRUE., the squeezed Coulomb kernel is used instead of the cosine window:[4]

This kernel squeezes contributions from large wave vectors into the window given by ENCUTGWSOFT. For GW type calculations the squeezed Coulomb kernel was the default (when ENCUTGWSOFT was set in the INCAR file) before version vasp.6.3, but in newer releases the code always defaults to a smoothed Coulomb kernel (both for GW and RPA type calculations). If one desires to recover the behavior of vasp.6.2 and older versions, LSCK=.TRUE. must be set in the INCAR file for GW type calculations if ENCUTGWSOFT is set in the INCAR file.


Mind: The infinite basis set limit extrapolation for RPA/ACFDT is described in more detail here.

Related tags and articles

PRECFOCK, ENCUT, ENCUTGW, GW calculations, LSCK, RPA/ACFDT basis set convergence, Examples that use this tag