Intégration de Verlet

L'algorithme de Verlet et sa forme en vitesses sont les intégrateurs standards de la dynamique moléculaire, appréciés pour leur réversibilité temporelle, leur caractère symplectique et leur excellente conservation de l'énergie sur les millions d'étapes qu'une simulation requiert.

Definition

L'intégration de Verlet est une méthode réversible dans le temps et symplectique pour l'intégration des équations du mouvement de Newton, qui met à jour les positions des particules en utilisant les positions actuelles et précédentes ainsi que l'accélération, produisant des trajectoires stables pour la dynamique moléculaire.

Scope

Ce sujet aborde la famille d'intégrateurs de Verlet : le schéma original de Verlet en positions, les formulations équivalentes de leapfrog et de Verlet en vitesses, leur réversibilité temporelle et leur structure symplectique, ainsi que la conservation de l'énergie à long terme qui en résulte. Il situe ces méthodes dans le cadre plus large de la théorie de l'intégration symplectique des systèmes hamiltoniens.

Core questions

• Comment le schéma de Verlet fait-il progresser les positions et les vitesses à partir des forces ?
• Pourquoi l'algorithme de Verlet est-il réversible dans le temps et symplectique ?
• Pourquoi l'intégration de Verlet conserve-t-elle bien l'énergie sur de très longues simulations ?
• Comment les formulations de Verlet en positions, leapfrog et Verlet en vitesses sont-elles liées les unes aux autres ?

Key theories

Structure symplectique et réversible dans le temps

L'intégration de Verlet préserve la géométrie symplectique de l'espace des phases et est invariante par renversement du temps, ce qui, ensemble, empêche la dérive énergétique systématique qui afflige les intégrateurs non symplectiques des systèmes conservatifs.

Conservation de l'hamiltonien d'ombre

Bien que la trajectoire discrète de Verlet ne conserve pas exactement l'énergie réelle, elle conserve presque un hamiltonien d'ombre étroitement lié, maintenant l'erreur énergétique bornée et oscillatoire plutôt que croissante.

Formulations équivalentes

Les schémas de Verlet en positions, leapfrog et Verlet en vitesses produisent la même trajectoire mais diffèrent dans la manière et le moment où les vitesses sont disponibles, le Verlet en vitesses étant préféré lorsque des positions et des vitesses synchronisées sont nécessaires.

Clinical relevance

L'intégration de Verlet est le moteur de pas de temps par défaut dans pratiquement tous les codes de dynamique moléculaire, des fluides de Lennard-Jones simples aux grandes simulations biomoléculaires, et le même principe symplectique est utilisé dans l'intégration orbitale à long terme en astronomie.

History

Le schéma a été utilisé par l'astronome Carl Stormer au début du XXe siècle et popularisé pour la simulation moléculaire par Loup Verlet dans son étude de 1967 sur les fluides de Lennard-Jones ; une analyse ultérieure a montré qu'il s'agissait d'un intégrateur symplectique, expliquant son excellente stabilité à long terme.

Key figures

• Loup Verlet
• Carl Stormer
• Ernst Hairer

Seminal works

• verlet1967
• hairer1993

Frequently asked questions

Pourquoi Verlet est-il préféré à une méthode de Runge-Kutta d'ordre supérieur en dynamique moléculaire ?

Bien que Runge-Kutta puisse être plus précis par pas, il n'est pas symplectique et dérive lentement en énergie sur de longues simulations. La structure symplectique et réversible dans le temps de Verlet maintient l'énergie bornée sur des millions d'étapes, ce qui est bien plus important que la précision par pas pour les simulations à l'équilibre.

Verlet conserve-t-il l'énergie exactement ?

Non. Il conserve un hamiltonien d'ombre proche plutôt que l'énergie exacte, ainsi l'énergie mesurée oscille dans une bande bornée au lieu de dériver, ce qui est suffisant pour le calcul de moyennes thermodynamiques stables.