⚗️ Surface d'énergie potentielle

C'est une surface d'énergie potentielle analytique bidimensionnelle largement utilisée, construite à partir d'une somme de quatre termes exponentiels de type gaussien, chacun avec sa propre amplitude, position et courbure. Elle a été conçue pour avoir exactement trois minima locaux et deux points-selles de premier ordre disposés de sorte qu'un chemin d'un minimum à un autre doive passer par les deux points-selles, en faisant un cas de test standard pour les algorithmes qui recherchent des chemins réactionnels et des états de transition.

La particule suit une dynamique de type Langevin : à chaque pas de temps, sa vitesse est mise à jour à l'aide du gradient négatif du potentiel (la force la tirant vers le bas), d'un terme de friction proportionnel à sa vitesse actuelle, et d'une secousse thermique aléatoire mise à l'échelle par le curseur de température. Le gradient est calculé avec une approximation par différences finies centrées, et la particule rebondit avec une vitesse réduite si elle atteint le bord du domaine tracé.

La friction amortit la vitesse de la particule à chaque pas, donc une friction plus élevée rend le mouvement plus lent et plus suramorti, se stabilisant dans un minimum sans beaucoup d'oscillation. La température ajoute une force aléatoire mise à l'échelle selon sa valeur à chaque pas, donnant à la particule suffisamment d'énergie cinétique pour parfois remonter au-dessus d'une barrière qu'elle a déjà franchie — exactement analogue à la façon dont l'énergie thermique dans une vraie molécule lui permet de franchir une barrière d'activation même si la réaction est ascendante en énergie potentielle.

L'énergie d'activation Ea est la différence d'énergie entre un état de transition et le puits réactif d'où il part — dans cette simulation, Ea = V(TS1) − V(R), lue en direct depuis le panneau d'affichage. Une énergie d'activation plus élevée signifie que la barrière entre réactif et produit est plus haute, donc à une température fixe, la particule a besoin d'une fluctuation thermique plus grande et plus rare pour la franchir, reflétant exactement pourquoi les vraies réactions à haute énergie d'activation se déroulent plus lentement à une température donnée.

De nombreuses réactions chimiques réelles ne vont pas directement des réactifs aux produits en une seule étape ; elles passent plutôt par un ou plusieurs intermédiaires — des espèces stables mais à plus haute énergie — chacun séparé de ses voisins par son propre état de transition. Les réactions catalysées par des enzymes, les substitutions nucléophiles SN1, et les mécanismes organiques à plusieurs étapes en sont des exemples classiques. Cartographier l'ensemble de la surface d'énergie potentielle, comme le fait cette simulation sous une forme bidimensionnelle simplifiée, est exactement la façon dont les chimistes informaticiens identifient ces intermédiaires et l'état de transition limitant la vitesse qui contrôle la vitesse globale de la réaction.