Chemistry Labs

Chimie théorique et computationnelle

Dynamique moléculaire (MD)

Comment l’intégration numérique des équations de Newton transforme un modèle d’interaction en trajectoires moléculaires et moyennes d’ensemble.

1 Collège2 Lycée3 Université4 Avancé5 Recherche

IntuitionSuivre les atomes au cours du temps

La dynamique moléculaire traite les atomes comme des particules en mouvement. Un champ de force fournit leurs interactions ; les lois de Newton font évoluer positions et vitesses par petits pas de temps, produisant une trajectoire plutôt qu’une structure statique.

Observez une trajectoire modèle et l’effet de la température et des paramètres d’interaction sur le mouvement atomique. Une trajectoire simulée est un échantillon fini, qui ne prouve pas à elle seule l’équilibre ni la convergence.

ScolaireForces, énergie et température

Définition: Modèle d’énergie potentielle

Un champ de force est une fonction U(r₁,…,rN) qui attribue une énergie potentielle à une configuration atomique. La force sur l’atome i est le gradient opposé, Fᵢ = −∇ᵢU.

mid2ri/dt2=Fi=−∇iUm_i d²r_i/dt² = F_i = −∇_i U

Exemple: Échelle de vitesse thermique

Pour une composante cartésienne de vitesse dans un système classique à température T, l’équipartition donne ½m⟨vₓ²⟩=½kBT. Quel effet une hausse de T a-t-elle sur l’échelle de vitesse quadratique moyenne ?

Solution

La vitesse quadratique moyenne de la composante vaut √(kBT/m) ; elle croît comme √T, et non linéairement avec T.

UniversitaireIntégration et mécanique statistique

Velocity-Verlet est très utilisé car il est réversible en temps et symplectique pour une dynamique conservative. Un pas met à jour la position avec vitesse et accélération actuelles, calcule les nouvelles forces, puis met à jour la vitesse avec l’accélération moyenne.

⟨A⟩NVT=∫Aexp(−βH)drdp/∫exp(−βH)drdp;β=1/(kBT)⟨A⟩_NVT = ∫A exp(−βH) dr dp / ∫exp(−βH) dr dp; β = 1/(k_B T)
EnsembleGrandeurs contrôléesUsage courant
NVEN, volume, énergieDynamique d’un système isolé
NVTN, volume, températurePropriétés à l’équilibre thermique
NPTN, pression, températureDensité à pression imposée

Définition: Échantillonnage et ergodicité

Les moyennes d’ensemble ne peuvent être estimées à partir d’une trajectoire que si la dynamique explore suffisamment l’espace des phases pertinent. Des images corrélées ne sont pas indépendantes ; l’incertitude demande une analyse par blocs ou des répliques indépendantes.

AvancéLimites du modèle et pratique avancée

Les champs de force classiques omettent généralement la réorganisation électronique explicite ; ils ne décrivent donc pas directement rupture de liaison, transfert de charge ou excitation électronique. Les champs réactifs et potentiels appris élargissent le modèle, mais exigent une validation sur des données quantiques adaptées et des contrôles prudents d’extrapolation.

RechercheFrontière : trajectoires plus longues et fidèles

Références