LTEMPER: Difference between revisions
No edit summary |
m (References added, equations added) |
||
| Line 5: | Line 5: | ||
---- | ---- | ||
VASP supports various modes where simultaneous calculations for different INCAR, KPOINTS, POTCAR, or POSCAR files are performed. | VASP supports various modes where simultaneous calculations for different INCAR, KPOINTS, POTCAR, or POSCAR files are performed. | ||
The parallel tempering mode is explained in this section, however, please read the section {{TAG|IMAGES}} first. | The parallel tempering mode is explained in this section, however, please read the section {{TAG|IMAGES}} first. Parallel tempering is also known as replica-exchange molecular dynamics.<ref name="sugi99"/> | ||
If the tags {{TAG| IMAGES}}=nn and {{TAG| LTEMPER}}=.TRUE. are set in the {{TAG| INCAR}} file, VASP performs parallel tempering calculations. In this case, it is expedient to supply different {{TAG| INCAR}} and {{TAG| POSCAR}} files in each subdirectory 01, 02, 03, ..., nn. For each subdiretory a different simulation temperature should be supplied using the tags {{TAG| TEBEG}} in the INCAR files 01/INCAR, 02/INCAR ... nn/INCAR. | If the tags {{TAG| IMAGES}}=nn and {{TAG| LTEMPER}}=.TRUE. are set in the {{TAG| INCAR}} file, VASP performs parallel tempering calculations. In this case, it is expedient to supply different {{TAG| INCAR}} and {{TAG| POSCAR}} files in each subdirectory 01, 02, 03, ..., nn. For each subdiretory a different simulation temperature should be supplied using the tags {{TAG| TEBEG}} in the INCAR files 01/INCAR, 02/INCAR ... nn/INCAR. | ||
In the course of the simulations, VASP will attempt to swap the temperatures between the images. Swapping attempts are made every | In the course of the simulations, VASP will attempt to swap the temperatures between the images. Swapping attempts are made every | ||
{{TAG|NTEMPER}} MD steps. Information on the current temperatures for each image can be found in the {{TAG|OUTCAR}} | {{TAG|NTEMPER}} MD steps, and accepted with a likelyhood of | ||
:<math> p = \min \left( 1, \frac{ \exp \left( -\frac{E_j}{kT_i} - \frac{E_i}{kT_j} \right) }{ \exp \left( -\frac{E_i}{kT_i} - \frac{E_j}{kT_j} \right) } \right) = \min \left( 1, e^{(E_i - E_j) \left( \frac{1}{kT_i} - \frac{1}{kT_j} \right)} \right) ,</math> | |||
where <math> E_i, E_j</math> and <math> T_i, T_j </math> are the energies and temperatures of the considered two replicas. Note that VASP swaps the temperatures and not the positions between images. This means that the temperatures in each subdirectory changes as the MD progresses. Information on the current temperatures for each image can be found in the {{TAG|OUTCAR}} files around the lines (all OUTCAR files will show the same information): | |||
parallel tempering new | parallel tempering new | ||
The average acceptance ratios are also written to the {{TAG|OUTCAR}} file. For efficient parallel tempering the acceptance ratio should not fall much below 0. | The average acceptance ratios are also written to the {{TAG|OUTCAR}} file. For efficient parallel tempering the acceptance ratio should not fall much below 0.3. If the acceptance ratio is too small, one usually needs to increase the number of images. | ||
== Related Tags and Sections == | == Related Tags and Sections == | ||
| Line 24: | Line 29: | ||
== References == | == References == | ||
<references> | <references> | ||
<ref name="sugi99">[http://dx.doi.org/10.1016/S0009-2614(99)01123-9 Y. Sugita and Y. Okamoto, Chemical Physics Letters. 314, 141 (1999).]</ref> | |||
</references> | </references> | ||
---- | ---- | ||
[[Category:INCAR]][[Category:Parallel Tempering]] | [[Category:INCAR]][[Category:Parallel Tempering]] | ||
Revision as of 07:49, 18 March 2022
LTEMPER = [logical]
Default: LTEMPER = .FALSE.
Description: LTEMPER specified whether parallel tempering is used. The flag must be used in combination with IMAGES.
VASP supports various modes where simultaneous calculations for different INCAR, KPOINTS, POTCAR, or POSCAR files are performed. The parallel tempering mode is explained in this section, however, please read the section IMAGES first. Parallel tempering is also known as replica-exchange molecular dynamics.[1]
If the tags IMAGES=nn and LTEMPER=.TRUE. are set in the INCAR file, VASP performs parallel tempering calculations. In this case, it is expedient to supply different INCAR and POSCAR files in each subdirectory 01, 02, 03, ..., nn. For each subdiretory a different simulation temperature should be supplied using the tags TEBEG in the INCAR files 01/INCAR, 02/INCAR ... nn/INCAR.
In the course of the simulations, VASP will attempt to swap the temperatures between the images. Swapping attempts are made every NTEMPER MD steps, and accepted with a likelyhood of
where and are the energies and temperatures of the considered two replicas. Note that VASP swaps the temperatures and not the positions between images. This means that the temperatures in each subdirectory changes as the MD progresses. Information on the current temperatures for each image can be found in the OUTCAR files around the lines (all OUTCAR files will show the same information):
parallel tempering new
The average acceptance ratios are also written to the OUTCAR file. For efficient parallel tempering the acceptance ratio should not fall much below 0.3. If the acceptance ratio is too small, one usually needs to increase the number of images.