Chimie théorique et computationnelle

Monte Carlo et modèles QM/MM

Deux outils complémentaires pour les systèmes moléculaires : l’échantillonnage statistique Monte Carlo et le partage quantique/classique QM/MM des régions réactives.

1 Collège 2 Lycée 3 Université 4 Avancé 5 Recherche

IntuitionDeux façons de rendre les modèles moléculaires utiles

Monte Carlo propose des configurations et les accepte selon des probabilités choisies pour échantillonner une distribution cible. QM/MM partage le système : la mécanique quantique traite les liaisons et électrons dans la région réactive, tandis que la mécanique moléculaire décrit à moindre coût l’environnement.

Explorez les déplacements d’essai et la séparation entre région quantique et environnement de mécanique moléculaire. L’affichage est conceptuel ; les conclusions quantitatives exigent un hamiltonien défini et un échantillonnage convergé.

ScolaireAccepter ou rejeter un essai

Définition: Acceptation de Metropolis

Pour une proposition symétrique dans l’ensemble canonique, un essai de variation d’énergie ΔU est accepté avec la probabilité min(1, exp(−βΔU)). Les descentes énergétiques sont acceptées ; les montées restent possibles thermiquement.

P_acc = min(1, exp(−βΔU)); β = 1/(k_B T)

Exemple: Déplacement défavorable thermiquement

À une température telle que βΔU = ln 4 pour un essai défavorable, quelle est la probabilité d’acceptation de Metropolis ?

Solution

Pacc = exp(−ln 4) = 1/4. Le déplacement est parfois accepté, permettant l’exploration thermique au-delà des minima locaux.

UniversitairePartitionner un système chimique

Définition: Énergie QM/MM

Une expression courante en embedding électrostatique est E=E_QM+E_MM+E_QM/MM. La région QM répond au potentiel électrostatique de l’environnement MM ; les liaisons covalentes à la frontière exigent des atomes de liaison ou un traitement plus élaboré.

E_QM/MM = E_QM + E_MM + E_coupling

Choix	Compromis
Région QM	Précision chimique contre coût de calcul
Embedding	Polarisation de l’environnement contre complexité
Proposition MC	Vitesse d’exploration contre acceptation et décorrélation

AvancéCoupler échantillonnage et chimie quantique

Monte Carlo peut échantillonner conformations ou coordonnées réactionnelles, tandis que les énergies QM/MM évaluent les configurations chimiquement importantes. Le bilan détaillé exige l’énergie cible et le rapport des probabilités de proposition ; les déplacements asymétriques ne peuvent utiliser aveuglément la formule simple de Metropolis.

RechercheFrontière : échantillonnage, embedding et incertitude

Références

Equation of State Calculations by Fast Computing Machines · N. Metropolis, A. W. Rosenbluth, M. N. Rosenbluth, A. H. Teller, E. Teller, 1953
Theoretical studies of enzymic reactions: dielectric, electrostatic and steric stabilization of the carbonium ion in the reaction of lysozyme · A. Warshel, M. Levitt, 1976
A combined quantum mechanical and molecular mechanical potential for molecular dynamics simulations · M. J. Field, P. A. Bash, M. Karplus, 1990