分子动力学模拟基础：力场、周期性边界条件、溶剂模型
做自由能计算，说白了，你得先学会跑分子动力学模拟。MD模拟就像是给分子拍电影，一帧一帧地记录原子怎么动。我个人习惯把MD叫做“计算机里的烧杯实验”——你在硅片上放一堆原子，定义好它们怎么相互作用，然后让时间往前走。
这章我们要聊的三个东西——力场、周期性边界条件和溶剂模型——是MD模拟的三大支柱。任何一个没搭好，你的自由能计算结果就跟空中楼阁似的，看着漂亮，一碰就塌。
1. 力场：分子世界的“物理规则”
力场是什么？就是一套数学公式加参数，用来描述原子里和原子间的相互作用。你可能觉得：量子化学不也能算相互作用吗？为什么非要用力场？
嗯，这里要注意：量子化学方法算一个几百原子的体系，算个几皮秒的轨迹，就得等上几个礼拜。而自由能计算常常需要跑几十纳秒甚至微秒级别的模拟。所以我们必须牺牲一点精度，换速度。力场就是这个妥协的产物。
典型的力场能量表达式长这样：
E_total = E_bond + E_angle + E_dihedral + E_vdW + E_electrostatic翻译过来就是：
- 键伸缩 (E_bond)：两个成键原子像弹簧一样振动
- 角弯曲 (E_angle)：三个原子构成的键角在平衡位置附近摆动
- 二面角 (E_dihedral)：四个原子绕键旋转的扭转势能
- 范德华 (E_vdW)：用Lennard-Jones势描述的非键相互作用
- 静电 (E_electrostatic)：库仑相互作用
常用的几大类力场，我做个小表格对比一下：
| 力场家族 | 代表性力场 | 特点 | 常见应用场景 |
|---|---|---|---|
| AMBER | ff14SB, ff19SB | 蛋白质/核酸参数覆盖全，水模型配套好 | 生物分子自由能计算 |
| CHARMM | C36, C36m | 脂质参数业界标准，蛋白-膜体系首选 | 膜蛋白模拟、药物膜通透性 |
| OPLS | OPLS-AA, OPLS4 | 有机分子参数丰富，小分子配体表现不错 | 药物-靶点结合自由能 |
| GROMOS | 54A7, 2016H66 | 联合原子模型，计算效率高 | 大体系、粗粒化过渡 |
我个人的建议：做自由能扰动（FEP）计算的时候，尽量保持力场的一致性。比如你配体和受体都用同一套力场参数，否则引入的“力场误差”会被自由能循环放大。我在项目中遇到过，把OPLS配体塞进AMBER受体里，结果算出来的结合自由能偏差了将近2 kcal/mol——这个误差直接导致一个先导化合物被错误筛掉。
2. 周期性边界条件：打破“盒子墙”
你想想看，模拟一个蛋白质在溶液中的行为，你的计算资源只允许你放几万个原子。但真正的溶液有无数个水分子。怎么办？
周期性边界条件（PBC）就是干这个的。你把模拟体系放在一个“盒子”里，然后在这个盒子周围复制无数个镜像。当一个原子从盒子左边跑出去，对应的镜像就从右边跑进来。
核心思想：你只用算一个盒子的原子，但通过镜像复制，模拟出无限大体系的“体相”环境。
常见的盒子形状有立方体、长方体、截角八面体、菱形十二面体等。我个人比较推荐用截角八面体或菱形十二面体——在同等最小镜像距离下，原子数更少，计算更快。
一个典型的PBC配置在输入文件里长这样：
periodic boundary conditions
cubic box side length: 60.0 Å
cutoff: 12.0 Å
neighbor list update: every 10 steps我曾经踩过的坑：
- 截断半径不能超过盒子最小边长的一半。否则一个原子会跟自己镜像里的自己相互作用——这叫“自映像”，物理上根本就是错的。
- 模拟蛋白质-配体体系时，盒子边缘到溶质至少留出10 Å的缓冲区。太紧了溶剂层不够，溶质会跟自己的镜像“聊天”。
3. 溶剂模型：水，不可或缺的配角
药物分子在生物体内工作，周围全是水。水不是惰性背景——氢键竞争、疏水效应、水桥作用，哪样都跟自由能直接相关。
溶剂模型分两大类：
3.1 显式水模型
最实在的方法——把每个水分子都放进去。常见的显式水模型有：
- SPC/E：三点模型，计算快，水密度和介电常数模拟得好
- TIP3P：三点模型，跟AMBER力场配套最熟
- TIP4P-Ew：四点模型，体相水性质更准，适合纯水体系
- OPC：四点模型，目前精度最高的一批
选水模型有个简单的原则——跟你用的力场保持一致。AMBER默认用TIP3P，CHARMM家族配的是modified TIP3P，OPLS用TIP4P或者SPC。
小提示：算自由能差异（比如结合自由能、相对结合自由能）时，水模型的选择对ΔΔG的影响往往比对ΔG小。但如果你非要把不同力场家族混着用——请准备好面对物理化学同行的灵魂拷问。
3.2 隐式水模型
不显式放水分子，用连续介电模型代替溶剂环境。常见的有：
- GB (Generalized Born)：算得快，适合小分子筛选
- PBSA (Poisson-Boltzmann)：精度比GB略好，但计算更贵
优点是你不用花95%的计算资源去跑一堆水分子。但隐式水对水桥、空腔水等细节的捕捉很差。我在做自由能微扰（FEP）时很少用隐式水——误差太大了。
知识体系全景
下面这张图，我手绘了本章的核心逻辑：
看到没，力场、PBC、溶剂模型这三块，各自支撑着MD模拟的一个侧面。你缺了任何一个，自由能计算就是无本之木。
以我的经验，很多新手把注意力全放在自由能方法本身——什么FEP、TI、WHAM背得滚瓜烂熟，结果模拟第一步的平衡都跑不稳。力场参数选错是家常便饭，PBC截断设个25 Å的也见了不少，还有人拿真空模拟来算结合自由能——你没看错，真有人这么干。