ENCUTGWSOFT: Difference between revisions
No edit summary |
(document the change from squeezed Coulomb kernel to smoothed one for newer VASP versions.) |
||
(3 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}}. | 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]], | |||
{{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