Verlet-Integration
Der Verlet-Algorithmus und seine Geschwindigkeitsform sind die Standard-Integratoren der Molekulardynamik, geschätzt, weil sie zeitumkehrbar und symplektisch sind und die Energie über die Millionen von Schritten, die eine Simulation erfordert, gut erhalten.
Definition
Die Verlet-Integration ist eine zeitumkehrbare, symplektische Methode zur Integration der Newtonschen Bewegungsgleichungen, die Partikelpositionen unter Verwendung aktueller und vorheriger Positionen sowie der Beschleunigung aktualisiert, was zu stabilen Trajektorien für die Molekulardynamik führt.
Scope
Dieses Thema behandelt die Verlet-Integrator-Familie: das ursprüngliche Positions-Verlet-Schema, die äquivalenten Leapfrog- und Geschwindigkeits-Verlet-Formulierungen, ihre Zeitumkehrbarkeit und symplektische Struktur sowie die daraus resultierende langfristige Energieerhaltung. Es ordnet diese Methoden in die umfassendere Theorie der symplektischen Integration Hamiltonscher Systeme ein.
Core questions
- Wie leitet das Verlet-Schema Positionen und Geschwindigkeiten aus den Kräften ab?
- Warum ist der Verlet-Algorithmus zeitumkehrbar und symplektisch?
- Warum erhält die Verlet-Integration die Energie über sehr lange Simulationen gut?
- Wie stehen die Positions-, Leapfrog- und Geschwindigkeits-Verlet-Formulierungen zueinander in Beziehung?
Key theories
- Symplektische und zeitumkehrbare Struktur
- Die Verlet-Integration bewahrt die symplektische Geometrie des Phasenraums und ist invariant unter Zeitumkehr, was zusammen die systematische Energiedrift verhindert, die nicht-symplektische Integratoren konservativer Systeme plagt.
- Erhaltung des Schatten-Hamiltonians
- Obwohl die diskrete Verlet-Trajektorie die wahre Energie nicht exakt erhält, erhält sie einen eng verwandten Schatten-Hamiltonian nahezu, wodurch der Energiefehler begrenzt und oszillierend bleibt, anstatt zu wachsen.
- Äquivalente Formulierungen
- Die Positions-Verlet-, Leapfrog- und Geschwindigkeits-Verlet-Schemata erzeugen dieselbe Trajektorie, unterscheiden sich jedoch darin, wie und wann Geschwindigkeiten verfügbar sind, wobei die Geschwindigkeits-Verlet-Methode bevorzugt wird, wenn synchronisierte Positionen und Geschwindigkeiten benötigt werden.
Clinical relevance
Die Verlet-Integration ist der standardmäßige Zeitentwicklungsmotor in praktisch allen Molekulardynamik-Codes, von einfachen Lennard-Jones-Fluiden bis hin zu großen biomolekularen Simulationen, und dasselbe symplektische Prinzip wird in der langfristigen Orbitalintegration in der Astronomie verwendet.
History
Das Schema wurde vom Astronomen Carl Stormer Anfang des 20. Jahrhunderts verwendet und von Loup Verlet in seiner Studie über Lennard-Jones-Fluide von 1967 für die molekulare Simulation populär gemacht; spätere Analysen zeigten, dass es sich um einen symplektischen Integrator handelt, was seine ausgezeichnete Langzeitstabilität erklärt.
Key figures
- Loup Verlet
- Carl Stormer
- Ernst Hairer
Related topics
Seminal works
- verlet1967
- hairer1993
Frequently asked questions
- Warum wird Verlet in der Molekulardynamik einem Runge-Kutta-Verfahren höherer Ordnung vorgezogen?
- Obwohl Runge-Kutta pro Schritt genauer sein kann, ist es nicht symplektisch und driftet über lange Läufe langsam in der Energie. Verlet's symplektische, zeitumkehrbare Struktur hält die Energie über Millionen von Schritten begrenzt, was für Gleichgewichtssimulationen weitaus wichtiger ist als die Genauigkeit pro Schritt.
- Erhält Verlet die Energie exakt?
- Nein. Es erhält einen nahegelegenen Schatten-Hamiltonian anstelle der exakten Energie, sodass die gemessene Energie innerhalb eines begrenzten Bereichs oszilliert, anstatt abzudriften, was für die Berechnung stabiler thermodynamischer Mittelwerte ausreicht.