In this article, we investigate molecular dynamics (MD) trajectories of a butane molecule, as obtained using different types of thermostats. Results show that at low temperature. where the harmonic approximation holds, the Nose'-Hoover (NH) thermostat fails to reproduce the statistical mechanical behavior, even using simulation lengths of millions of time steps, whereas the Gaussian isokinetic (IG) thermostat reproduces quite well the expected statistical mechanical values. The Berendsen's coupling (BC) provides good results for basic properties such as the average potential and kinetic energies but fails in reproducing the canonical fluctuations. Moreover, using the speed of divergence of initially nearby trajectories in phase space as a measure of the dynamical chaoticity, we found that the NH thermostat provides very slow divergence for the physical phase space degrees of freedom, concentrating most of its chaoticity in the dynamics of the thermostat virtual degree of freedom. On the contrary, the IG thermostat provides always highly diverging trajectories in phase space, characterized by a high chaoticity of each degree of freedom. Finally, the BC thermostat provides a moderate chaotic behavior for all of the degrees of freedom. Such results suggest that even assuming for both the "rigorous" algorithms (NH and IG) a full ergodic behavior the NH thermostat could require an extremely long time to achieve convergence of the time averaged properties.

Dynamical and statistical mechanical characterization of temperature coupling algorithms / D'Alessandro, Maira; Tenenbaum, Alexander; Andrea, Amadei. - In: JOURNAL OF PHYSICAL CHEMISTRY. B, CONDENSED MATTER, MATERIALS, SURFACES, INTERFACES & BIOPHYSICAL. - ISSN 1520-6106. - STAMPA. - 106:19(2002), pp. 5050-5057. [10.1021/jp013689i]

Dynamical and statistical mechanical characterization of temperature coupling algorithms

D'ALESSANDRO, Maira;TENENBAUM, Alexander;
2002

Abstract

In this article, we investigate molecular dynamics (MD) trajectories of a butane molecule, as obtained using different types of thermostats. Results show that at low temperature. where the harmonic approximation holds, the Nose'-Hoover (NH) thermostat fails to reproduce the statistical mechanical behavior, even using simulation lengths of millions of time steps, whereas the Gaussian isokinetic (IG) thermostat reproduces quite well the expected statistical mechanical values. The Berendsen's coupling (BC) provides good results for basic properties such as the average potential and kinetic energies but fails in reproducing the canonical fluctuations. Moreover, using the speed of divergence of initially nearby trajectories in phase space as a measure of the dynamical chaoticity, we found that the NH thermostat provides very slow divergence for the physical phase space degrees of freedom, concentrating most of its chaoticity in the dynamics of the thermostat virtual degree of freedom. On the contrary, the IG thermostat provides always highly diverging trajectories in phase space, characterized by a high chaoticity of each degree of freedom. Finally, the BC thermostat provides a moderate chaotic behavior for all of the degrees of freedom. Such results suggest that even assuming for both the "rigorous" algorithms (NH and IG) a full ergodic behavior the NH thermostat could require an extremely long time to achieve convergence of the time averaged properties.
2002
01 Pubblicazione su rivista::01a Articolo in rivista
Dynamical and statistical mechanical characterization of temperature coupling algorithms / D'Alessandro, Maira; Tenenbaum, Alexander; Andrea, Amadei. - In: JOURNAL OF PHYSICAL CHEMISTRY. B, CONDENSED MATTER, MATERIALS, SURFACES, INTERFACES & BIOPHYSICAL. - ISSN 1520-6106. - STAMPA. - 106:19(2002), pp. 5050-5057. [10.1021/jp013689i]
File allegati a questo prodotto
Non ci sono file associati a questo prodotto.

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11573/33841
 Attenzione

Attenzione! I dati visualizzati non sono stati sottoposti a validazione da parte dell'ateneo

Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 44
  • ???jsp.display-item.citation.isi??? 42
social impact