Theoretical and computational chemistry
Monte Carlo and QM/MM models
Two complementary tools for molecular systems: statistical sampling by Monte Carlo and a quantum/classical partition for chemically active regions.
IntuitionTwo ways to make molecular models useful
Monte Carlo proposes configurations and accepts them with carefully chosen probabilities to sample a target distribution. QM/MM instead divides a system: quantum mechanics treats bonds and electrons where chemistry changes, while molecular mechanics handles the surrounding environment more cheaply.
SchoolAccepting or rejecting a trial
Definition: Metropolis acceptance
For a symmetric proposal in the canonical ensemble, a trial with energy change ΔU is accepted with probability min(1, exp(−βΔU)). Downhill moves are accepted; uphill moves can still occur thermally.
Example: A thermally uphill move
At a temperature where βΔU = ln 4 for an uphill trial, what is the Metropolis acceptance probability?
Solution
Pacc = exp(−ln 4) = 1/4. The move is sometimes accepted, allowing thermal exploration beyond local minima.
UndergraduatePartitioning a chemical system
Definition: QM/MM energy
A common electrostatic-embedding expression is E = E_QM + E_MM + E_QM/MM. The QM region responds to the MM environment’s electrostatic potential; cross-boundary covalent bonds require link atoms or a more sophisticated boundary treatment.
| Choice | Trade-off |
|---|---|
| QM region | Accuracy for chemistry vs computational cost |
| Embedding | Environmental polarization vs complexity |
| MC proposal | Exploration rate vs acceptance and decorrelation |
AdvancedCoupling sampling to quantum chemistry
Monte Carlo can sample conformations or reaction coordinates, while QM/MM energies evaluate chemically important configurations. For detailed balance, the acceptance rule must use the target energy and the proposal probability ratio; asymmetric moves cannot blindly use the simple Metropolis expression.
ResearchFrontier: sampling, embedding, and uncertainty
References
- 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