理论与计算化学
势能面与自由能计算
将分子几何与势能景观联系起来,再用统计力学和增强采样估计自由能差与能垒。
直觉能量景观与可达状态
势能面为每种分子构型赋予能量。斜率表示力,谷底暗示稳定结构,但温度会改变统计占据区域:决定平衡概率的是自由能,而不仅是势能。
中学阶段读懂能量图
定义: 势能面
固定电子态与核坐标 R 时,V(R) 即势能。驻点梯度为零;局部极小在所有非约束方向曲率为正,而一级鞍点有一个下降方向。
例题: 反应坐标上的能垒
在简化的一维坐标上,反应物极小值比过渡态极大值低 12 kJ mol⁻¹。从反应物到过渡态的势能垒是多少?
解答
能垒为 V‡−V反应物=12 kJ mol⁻¹。这是势能差,并不自动等同于有限温度下的活化自由能。
大学从微观态到自由能
定义: 平均力势
对集体变量 ξ,自由能剖面 W(ξ) 由平衡概率密度 p(ξ) 定义:W(ξ)=−kBT ln p(ξ)+常数。未显式保留在 ξ 中的其他坐标都已积分消去。
两个能谷的平衡占据概率决定自由能差,而非仅由各极小点能量决定。熵贡献反映每个能谷占据的构型空间体积。
| 方法 | 策略 | 注意事项 |
|---|---|---|
| 伞形采样 | 对重叠窗口施加偏置 | 窗口重叠与重加权至关重要 |
| 元动力学 | 沿集体变量累积偏置 | 不佳变量会掩盖慢模态 |
| 自由能微扰 | 在状态间重加权样本 | 需要相空间重叠 |
进阶能垒、坐标与不确定度
反应坐标是约化描述;选择不当可能隐藏正交慢重排并产生滞后。自由能方法需证明窗口间重叠、可逆性或收敛,并估计不确定度;仅有视觉平滑的曲线并不足以构成证据。
科研前沿:稀有事件与可靠自由能
参考文献
- Understanding Molecular Simulation: From Algorithms to Applications · D. Frenkel, B. Smit, 2002
- Optimized Monte Carlo data analysis · A. M. Ferrenberg, R. H. Swendsen, 1989
- Calculating free energy differences from computer simulations using perturbation theory · R. W. Zwanzig, 1954