TEBEG: Difference between revisions

From VASP Wiki
(Created page with '{{TAGDEF|TEBEG|[real]|0}} Description: {{TAG|TEBEG}} sets the start temperature for an ab-initio molecular dynamics run ({{TAG|IBRION}}=0). ---- If no initial velocities are sup…')
 
No edit summary
 
(29 intermediate revisions by 3 users not shown)
Line 1: Line 1:
{{TAGDEF|TEBEG|[real]|0}}
{{TAGDEF|TEBEG|[real]|0}}


Description: {{TAG|TEBEG}} sets the start temperature for an ab-initio molecular dynamics run ({{TAG|IBRION}}=0).
Description: {{TAG|TEBEG}} sets the starting temperature (in K) for an ab-initio molecular dynamics run ({{TAG|IBRION}}=0) and other routines (e.g. {{TAG|Electron-phonon interactions from Monte-Carlo sampling}}).
----
----
If no initial velocities are supplied on the {{FILE|POSCAR}} file, the velocities are set randomly according to a Maxwell-Boltzmann distribution at the initial temperature {{TAG|TEBEG}}. Velocities are only used for molecular dynamics ({{TAG|IBRION}}=0).
If no initial velocities are supplied on the {{FILE|POSCAR}} file, the velocities are set randomly according to a Maxwell-Boltzmann distribution at the initial temperature {{TAG|TEBEG}}. Velocities are only used for molecular dynamics ({{TAG|IBRION}}=0).


'''Mind''': VASP defines the temperature as  
'''Mind''': If {{TAG|MDALGO}}>0 is used VASP defines the temperature as  
:<math>
:<math>
T= \frac{1}{3 k_B T N_{\rm ions}} \sum_n M_n | \vec v_n |^2.
T= \frac{1}{ k_B T 3 (N_{\rm ions}-N_{\rm constraints})} \sum\limits_{n}^{N_{\rm ions}} M_n | \vec v_n |^2.
</math>
</math>
But, because the center of mass is conserved, there are only 3(N<sub>ions</sub>-1) degrees of freedom (the sum of all velocities is zero, if a random initialization is chosen). This means that the real simulation temperature is: T={{TAG|TEBEG}}&times;N<sub>ions</sub>/(N<sub>ions</sub>-1)


Consequently, the temperature written by VASP (e.g. in the {{FILE|OUTCAR}} file) 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>&times;(N<sub>ions</sub>-1)/N<sub>ions</sub>
This temperature ist written to the {{FILE|OUTCAR}} file.
Depending on the type of thermostat this temperature has to be rescaled to obtain the real simulation temperature.  


== Related Tags and Sections ==
*{{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 also taken into account in the {{TAG|OUTCAR}} file.
 
*{{TAG|Andersen thermostat}}:
Same as for {{TAG|Nose-Hoover thermostat}}.
 
*{{TAG|Langevin thermostat}}:
As for the {{TAG|Nose-Hoover thermostat}} and {{TAG|Andersen 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, hence this method has 3 degrees of freedom more than the {{TAG|Nose-Hoover thermostat}} and {{TAG|Andersen thermostat}}.
 
== Related tags and articles ==
{{TAG|TEEND}},
{{TAG|IBRION}},
{{TAG|IBRION}},
{{TAG|SMASS}}
{{TAG|SMASS}}
{{sc|TEBEG|Examples|Examples that use this tag}}
----
----
[[The_VASP_Manual|Contents]]


[[Category:INCAR]][[Category:Dynamics]]
[[Category:INCAR tag]][[Category:Molecular dynamics]]

Latest revision as of 12:22, 7 September 2022

TEBEG = [real]
Default: TEBEG = 0 

Description: TEBEG sets the starting temperature (in K) 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: If MDALGO>0 is used 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.

Same as for Nose-Hoover thermostat.

As for the Nose-Hoover thermostat and Andersen 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, hence this method has 3 degrees of freedom more than the Nose-Hoover thermostat and Andersen thermostat.

Related tags and articles

TEEND, IBRION, SMASS

Examples that use this tag