Berendsen thermostat: Difference between revisions
| Line 7: | Line 7: | ||
Though the thermostat does not generate a correct [[canonical ensemble]] (especially for small systems), for large systems on the order of hundreds or thousands of atoms/molecules, the approximation yields roughly correct results for most calculated properties.<ref>{{cite journal | first=T. | last=Morishita| journal=The Journal of Chemical Physics| title=Fluctuation formulas in molecular-dynamics simulations with the weak coupling heat bath | year=2000 | volume=113 | pages=2976–2982| doi=10.1063/1.1287333 }}</ref> The scheme is widely used due to the efficiency with which it relaxes a system to some target (bath) temperature. In many instances, systems are initially equilibrated using the Berendsen scheme, while properties are calculated using the widely known [[Nosé-Hoover thermostat]], which correctly generates trajectories consistent with a canonical ensemble. |
Though the thermostat does not generate a correct [[canonical ensemble]] (especially for small systems), for large systems on the order of hundreds or thousands of atoms/molecules, the approximation yields roughly correct results for most calculated properties.<ref>{{cite journal | first=T. | last=Morishita| journal=The Journal of Chemical Physics| title=Fluctuation formulas in molecular-dynamics simulations with the weak coupling heat bath | year=2000 | volume=113 | pages=2976–2982| doi=10.1063/1.1287333 }}</ref> The scheme is widely used due to the efficiency with which it relaxes a system to some target (bath) temperature. In many instances, systems are initially equilibrated using the Berendsen scheme, while properties are calculated using the widely known [[Nosé-Hoover thermostat]], which correctly generates trajectories consistent with a canonical ensemble. |
||
==See also== |
|||
* [[Molecular mechanics]] |
|||
* [[List of software for molecular mechanics modeling|Software for molecular mechanics modeling]] |
|||
==References== |
==References== |
||
Revision as of 12:03, 18 January 2015
The Berendsen thermostat[1] is an algorithm to re-scale the velocities of particles in molecular dynamics simulations to control the simulation temperature.
Basic Description
In this scheme, the system is weakly coupled to a heat bath with some temperature. The thermostat suppresses fluctuations of the kinetic energy of the system and therefore cannot produce trajectories consistent with the canonical ensemble. The temperature of the system is corrected such that the deviation exponentially decays with some time constant .
Though the thermostat does not generate a correct canonical ensemble (especially for small systems), for large systems on the order of hundreds or thousands of atoms/molecules, the approximation yields roughly correct results for most calculated properties.[2] The scheme is widely used due to the efficiency with which it relaxes a system to some target (bath) temperature. In many instances, systems are initially equilibrated using the Berendsen scheme, while properties are calculated using the widely known Nosé-Hoover thermostat, which correctly generates trajectories consistent with a canonical ensemble.
See also
References
- ^ Berendsen, H. J. C.; Postma, J. P. M.; van Gunsteren, W. F.; DiNola, A.; Haak, J. R. (1984). "Molecular-Dynamics with Coupling to an External Bath". Journal of Chemical Physics. 81 (8): 3684–3690. Bibcode:1984JChPh..81.3684B. doi:10.1063/1.448118.
- ^ Morishita, T. (2000). "Fluctuation formulas in molecular-dynamics simulations with the weak coupling heat bath". The Journal of Chemical Physics. 113: 2976–2982. doi:10.1063/1.1287333.
 | This physics-related article is a stub. You can help Wikipedia by expanding it. |