Andersen-Thermostat

Das Andersen-Thermostat ist in der Molekulardynamik eine Methode, mit der ein klassisches Vielteilchensystem stochastisch an ein Wärmebad gekoppelt werden kann, dadurch sind Simulationen bei einer konstanten Temperatur (kanonisches Ensemble) möglich.

Problemstellung und Funktionsweise

Werden in einer Molekulardynamik-Simulation die newtonschen Bewegungsgleichungen ohne Modifikationen gelöst, so ist die Energie exakt erhalten, was dem mikrokanonischen Ensemble (NVE) der statistischen Mechanik entspricht. In den meisten Fällen ist man jedoch daran interessiert ein System bei einer konstanten, vorgegebenen Temperatur zu simulieren, also ein System, welches an ein Wärmebad gekoppelt ist.

Das Andersen-Thermostat führt ein (virtuelles) Wärmebad ein, welches mit den Teilchen stochastisch Stöße ausführt und durch die stochastische Stoßfrequenz $\nu$ charakterisiert wird. Die Zeitdifferenz zweier aufeinander folgender Kollisionen folgt dann einer Poisson-Verteilung:

P(t;\nu )=\nu e^{-\nu t}.

Ein Zeitschritt mit dem Andersen-Thermostat läuft dann folgendermaßen ab:

Numerische Integration der Bewegungsgleichungen für den Zeitschritt $\Delta t$ .
Auswahl der Teilchen, welche an einem Stoß mit dem Wärmebad beteiligt sind; die Auswahlwahrscheinlichkeit beträgt $\nu \Delta t$ .
Kollisionsschritt: für jedes Teilchen, welches mit dem Bad stößt, wird eine neue (zufällige) Geschwindigkeit aus der Maxwell-Boltzmann-Verteilung der Temperatur $T$ gezogen

Die stochastischen Kollisionen verändern die Dynamik des Systems auf eine unphysikalische Weise, dynamische Eigenschaften wie Diffusionskoeffizienten können mit der Methode daher nicht bestimmt werden.

Einzelnachweise

↑ Hans C. Andersen: Molecular dynamics simulations at constant pressure and/or temperature. In: The Journal of Chemical Physics. Band 72, Nr. 4, 15. Februar 1980, ISSN 0021-9606, S. 2384–2393, doi:10.1063/1.439486.
↑ Daan Frenkel, Berend Smit: Understanding molecular simulation: From algorithms to applications (= Computational science series . Band 1). 2. Auflage. Academic Press, San Diego [u. a.] 2002, ISBN 978-0-12-267351-1, doi:10.1016/b978-0-12-267351-1.x5000-7.

This article is issued from Wikipedia. The text is licensed under Creative Commons - Attribution - Sharealike. Additional terms may apply for the media files.

[1] Hans C. Andersen: Molecular dynamics simulations at constant pressure and/or temperature. In: The Journal of Chemical Physics. Band 72, Nr. 4, 15. Februar 1980, ISSN 0021-9606, S. 2384–2393, doi:10.1063/1.439486.

[2] Daan Frenkel, Berend Smit: Understanding molecular simulation: From algorithms to applications (= Computational science series . Band 1). 2. Auflage. Academic Press, San Diego [u. a.] 2002, ISBN 978-0-12-267351-1, doi:10.1016/b978-0-12-267351-1.x5000-7.