Simulation de dynamique moléculaire

La dynamique moléculaire intègre les équations du mouvement de Newton pour un système d'atomes en interaction, générant des trajectoires à partir desquelles sont calculées les propriétés structurales, dynamiques et thermodynamiques.

Definition

Une technique de simulation qui propage les positions et les vitesses atomiques dans le temps sous l'influence des forces dérivant d'une fonction d'énergie potentielle, produisant une trajectoire déterministe à travers l'espace des phases.

Scope

Couvre l'intégration numérique des équations classiques du mouvement, les intégrateurs réversibles dans le temps tels que la famille de Verlet, le contrôle de la température et de la pression via des thermostats et des barostats, le traitement des interactions électrostatiques à longue portée et des conditions aux limites périodiques, ainsi que l'extraction des propriétés d'équilibre et de transport à partir des trajectoires.

Core questions

• Comment l'équation du mouvement de Newton est-elle intégrée de manière stable sur de nombreux pas de temps ?
• Comment la température et la pression sont-elles contrôlées pour échantillonner un ensemble thermodynamique choisi ?
• Comment les interactions électrostatiques à longue portée sont-elles gérées efficacement ?
• Comment les propriétés observables sont-elles récupérées à partir d'une trajectoire finie ?

Key theories

Intégration de Verlet

Un schéma simple, réversible dans le temps et symplectique pour l'intégration des équations du mouvement qui conserve bien l'énergie sur de longues simulations et sous-tend la plupart des codes de dynamique moléculaire.

Échantillonnage ergodique

Selon l'hypothèse ergodique, les moyennes temporelles le long d'une trajectoire suffisamment longue sont égales aux moyennes d'ensemble, reliant ainsi la simulation aux observables de la mécanique statistique.

Mechanisms

Chaque étape calcule les forces à partir du gradient de potentiel, fait avancer les positions et les vitesses par l'intégrateur, et applique des corrections de thermostat ou de barostat ; la répétition de ce processus construit une trajectoire dont les moyennes fournissent les propriétés thermodynamiques et dynamiques.

Clinical relevance

La dynamique moléculaire révèle les changements conformationnels, les événements de liaison, la diffusion et les processus adjacents aux réactions dans les protéines, les membranes et les matériaux, offrant une compréhension mécanistique qui complète l'expérimentation en biophysique et en découverte de médicaments.

History

Suite aux simulations de sphères dures d'Alder et Wainwright et à l'étude du potentiel continu de Rahman sur l'argon liquide, les algorithmes de Verlet de 1967 et le développement ultérieur des thermostats, barostats et de l'électrostatique basée sur la méthode d'Ewald ont établi la dynamique moléculaire comme une discipline de simulation mature.

Key figures

• Loup Verlet
• Berni Alder
• Aneesur Rahman
• Herman Berendsen

Seminal works

• verlet1967
• frenkel2002

Frequently asked questions

Qu'est-ce qui limite le pas de temps en dynamique moléculaire ?

Les mouvements les plus rapides, typiquement les vibrations de liaison impliquant l'hydrogène, fixent cette limite ; la contrainte de ces liaisons permet un pas de temps de quelques femtosecondes, tandis que les temps totaux accessibles atteignent les microsecondes ou plus.

La dynamique moléculaire est-elle déterministe ?

L'intégration est déterministe étant donné les conditions initiales, mais les trajectoires sont chaotiques, de sorte que de minuscules différences s'amplifient rapidement ; les résultats significatifs proviennent des moyennes statistiques plutôt que des chemins individuels.