Chemistry Labs

理论与计算化学

蒙特卡罗方法与 QM/MM 模型

介绍分子体系的两种互补工具：蒙特卡罗统计采样，以及针对化学活性区域的量子/经典分区。

1 初中2 高中3 大学4 进阶5 科研

直觉让分子模型发挥作用的两种方法

蒙特卡罗提出构型，并按精心选择的概率接受，以采样目标分布。QM/MM 则划分体系：化学变化区域的键与电子由量子力学处理，周围环境由成本较低的分子力学处理。

探索试探移动及量子区域与分子力学环境的划分。图示是概念模型；定量结论需要明确的哈密顿量和收敛采样。

中学阶段接受或拒绝试探构型

定义: Metropolis 接受准则

在正则系综的对称提议下，能量变化 ΔU 的试探以 min(1, exp(−βΔU)) 的概率接受。能量下降的移动必然接受；热涨落也允许能量上升。

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

例题: 热力学不利移动

对于满足 βΔU=ln 4 的上坡试探，Metropolis 接受概率是多少？

解答

Pacc=exp(−ln 4)=1/4。该移动有时会被接受，从而允许热涨落探索越过局部极小值。

大学划分化学体系

定义: QM/MM 能量

常见静电嵌入表达式为 E=E_QM+E_MM+E_QM/MM。QM 区域响应 MM 环境的静电势；跨边界共价键需要连接原子或更精细的边界处理。

EQM/MM=EQM+EMM+EcouplingE_QM/MM = E_QM + E_MM + E_coupling
选择权衡
QM 区域化学精度与计算成本
嵌入环境极化与复杂度
MC 提议探索速率与接受率、去相关

进阶将采样与量子化学结合

蒙特卡罗可采样构象或反应坐标，再以 QM/MM 能量评估化学关键构型。满足细致平衡时，接受规则必须纳入目标能量与提议概率比；非对称移动不能直接套用简单 Metropolis 公式。

科研前沿：采样、嵌入与不确定度

参考文献

  • 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