LRPAFORCE: Difference between revisions
No edit summary |
No edit summary |
||
| (3 intermediate revisions by 2 users not shown) | |||
| Line 6: | Line 6: | ||
This tag can be optionally set in low-scaling [[ACFDT/RPA_calculations#RPAFORCES|RPA calculations]] or [[GW_calculations#RPAFORCES| GW calculations]]. It allows computing the [[Forces#RPA-forces|RPA forces]]{{cite|ramberger:prl:118}} on each ion. | This tag can be optionally set in low-scaling [[ACFDT/RPA_calculations#RPAFORCES|RPA calculations]] or [[GW_calculations#RPAFORCES| GW calculations]]. It allows computing the [[Forces#RPA-forces|RPA forces]]{{cite|ramberger:prl:118}} on each ion. | ||
Setting | |||
{{TAGBL|ALGO}} = RPAR ; LRPAFORCE = .TRUE. | {{TAGBL|ALGO}} = RPAR ; LRPAFORCE = .TRUE. | ||
or equivalently | |||
{{TAGBL|ALGO}} = G0W0R ; LRPAFORCE = .TRUE. | {{TAGBL|ALGO}} = G0W0R ; LRPAFORCE = .TRUE. | ||
determines | determines the RPA total energy with corresponding forces and the quasiparticle energies within the GW approximation. | ||
The {{TAG|LRPAFORCE}} tag can be used in combination with the standard relaxation options {{TAG|IBRION}} and {{TAG|NSW}}. | The {{TAG|LRPAFORCE}} tag can be used in combination with the standard relaxation options {{TAG|IBRION}} and {{TAG|NSW}}. | ||
Generally, the energy calculated by the RPA can be quite noisy as a function of the ionic positions, in particular, if {{TAG|PRECFOCK}} = FAST and {{TAG|NMAXFOCKAE}} = 1 is set | |||
(these are the default values for RPA calculations). Most of the noise is related to the exact exchange energy, and we are working on methods to improve this issue. | |||
Currently, to reduce the noise in the energy and forces, it is sensible to set {{TAG|PRECFOCK}} = Normal (typically doubling the execution time and memory requirement). It is also possible to set {{TAG|LMAXFOCKAE}} = -1 (which implicitly sets {{TAG|NMAXFOCKAE}} = 0). This makes the correlation energies and the related forces less noisy, but technically less accurate (i.e. part of the correlation energy will be missing at high transition energies). | |||
Overall, RPA forces must be used carefully and only after extensive testing of all relevant parameters. | |||
{{NB|mind|The RPA stress tensor is not available.}} | {{NB|mind|The RPA stress tensor is not available.}} | ||
{{NB|warning|Only {{TAG|ISIF}}{{=}}0 is supported.}} | {{NB|warning|Only {{TAG|ISIF}}{{=}}0 is supported.}} | ||
== Related tags and articles == | == Related tags and articles == | ||
| Line 29: | Line 34: | ||
---- | ---- | ||
[[Category:INCAR tag]][[Category:GW | [[Category:INCAR tag]][[Category:GW]][[Category:Molecular dynamics]][[Category:Ionic minimization]][[Category:Forces]][[Category:Low-scaling GW and RPA]] | ||
== References == | == References == | ||
Latest revision as of 11:48, 10 December 2024
LRPAFORCE = .TRUE. | .FALSE.
Default: LRPAFORCE = .FALSE.
Description: LRPAFORCE=.TRUE. calculates forces in the random-phase approximation (RPA).
Available as of VASP.6.1.
This tag can be optionally set in low-scaling RPA calculations or GW calculations. It allows computing the RPA forces[1] on each ion. Setting
ALGO = RPAR ; LRPAFORCE = .TRUE.
or equivalently
ALGO = G0W0R ; LRPAFORCE = .TRUE.
determines the RPA total energy with corresponding forces and the quasiparticle energies within the GW approximation.
The LRPAFORCE tag can be used in combination with the standard relaxation options IBRION and NSW.
Generally, the energy calculated by the RPA can be quite noisy as a function of the ionic positions, in particular, if PRECFOCK = FAST and NMAXFOCKAE = 1 is set (these are the default values for RPA calculations). Most of the noise is related to the exact exchange energy, and we are working on methods to improve this issue. Currently, to reduce the noise in the energy and forces, it is sensible to set PRECFOCK = Normal (typically doubling the execution time and memory requirement). It is also possible to set LMAXFOCKAE = -1 (which implicitly sets NMAXFOCKAE = 0). This makes the correlation energies and the related forces less noisy, but technically less accurate (i.e. part of the correlation energy will be missing at high transition energies). Overall, RPA forces must be used carefully and only after extensive testing of all relevant parameters.
| Mind: The RPA stress tensor is not available. |
| Warning: Only ISIF=0 is supported. |
Related tags and articles
NBANDS, ALGO, IBRION, NSW, POTIM, RPA calculations, GW calculations