TEBEG: Difference between revisions
No edit summary |
No edit summary |
||
Line 14: | Line 14: | ||
*{{TAG|Nose-Hoover thermostat}}: | *{{TAG|Nose-Hoover thermostat}}: | ||
In this thermostat the number of degrees of freedom including constraines are already acounted for in the potential energy term. In this this method the center of mass is conserved. This lowers the degrees of freedom by one which is | In this thermostat the number of degrees of freedom including constraines are already acounted for in the potential energy term. In this this method the center of mass is conserved. This lowers the degrees of freedom by one which is also taken into account in the {{TAG|OUTCAR}} file. Hence the temperature output in the {{TAG|OUTCAR}} file is the actual simulation temperature. This means that the real simulation temperature is: T={{TAG|TEBEG}}×N<sub>ions</sub>/(N<sub>ions</sub>-1). | ||
Consequently, the temperature written by VASP (e.g. in the is incorrect and has to be corrected in accordance with the above. Usually the effect is rather small and subtle, but one should correct the error if very precise results are required; in that case the temperature should be specified according to: {{TAG|TEBEG}}=T<sub>requested</sub>×(N<sub>ions</sub>-1)/N<sub>ions</sub>. | Consequently, the temperature written by VASP (e.g. in the is incorrect and has to be corrected in accordance with the above. Usually the effect is rather small and subtle, but one should correct the error if very precise results are required; in that case the temperature should be specified according to: {{TAG|TEBEG}}=T<sub>requested</sub>×(N<sub>ions</sub>-1)/N<sub>ions</sub>. | ||
*{{TAG|Andersen thermostat}}: | |||
Same as for {{TAG|Nose-Hoover thermostat}}. | |||
*{{TAG|Langevin thermostat}}: | *{{TAG|Langevin thermostat}}: | ||
As for the {{TAG|Nose-Hoover thermostat}} in this thermostat the number of degrees of freedom including constraines are already acounted for. The center of mass is not conserved in this method since stochastic forces are acting on the atoms. | As for the {{TAG|Nose-Hoover thermostat}} in this thermostat the number of degrees of freedom including constraines are already acounted for. The center of mass is also subtracted but the centre of mass is not conserved in this method since stochastic forces are acting on the atoms. Consequently, the temperature in the {{TAG|OUTCAR}} file has to be corrected to T={{TAG|TEBEG}}×(N<sub>ions</sub>-1)/N<sub>ions</sub> | ||
Revision as of 10:32, 27 June 2019
TEBEG = [real]
Default: TEBEG = 0
Description: TEBEG sets the start temperature for an ab-initio molecular dynamics run (IBRION=0) and other routines (e.g. Electron-phonon interactions from Monte-Carlo sampling).
If no initial velocities are supplied on the POSCAR file, the velocities are set randomly according to a Maxwell-Boltzmann distribution at the initial temperature TEBEG. Velocities are only used for molecular dynamics (IBRION=0).
Mind: VASP defines the temperature as
This temperature ist written to the OUTCAR file. Depending on the type of thermostat this temperature has to be rescaled to obtain the real simulation temperature.
In this thermostat the number of degrees of freedom including constraines are already acounted for in the potential energy term. In this this method the center of mass is conserved. This lowers the degrees of freedom by one which is also taken into account in the OUTCAR file. Hence the temperature output in the OUTCAR file is the actual simulation temperature. This means that the real simulation temperature is: T=TEBEG×Nions/(Nions-1). Consequently, the temperature written by VASP (e.g. in the is incorrect and has to be corrected in accordance with the above. Usually the effect is rather small and subtle, but one should correct the error if very precise results are required; in that case the temperature should be specified according to: TEBEG=Trequested×(Nions-1)/Nions.
Same as for Nose-Hoover thermostat.
As for the Nose-Hoover thermostat in this thermostat the number of degrees of freedom including constraines are already acounted for. The center of mass is also subtracted but the centre of mass is not conserved in this method since stochastic forces are acting on the atoms. Consequently, the temperature in the OUTCAR file has to be corrected to T=TEBEG×(Nions-1)/Nions