第4章：常用力场介绍——OPLS-AA、COMPASS、PCFF、GAFF、TraPPE-UA
力场这东西，说白了就是分子模拟的「游戏规则」。你选对了力场，结果八九不离十；选错了，算出来的东西可能连定性趋势都是反的。我做了这么多年高分子模拟，踩过的坑不少，今天就跟你聊聊几个主流力场该怎么选。
4.1 OPLS-AA：有机小分子和蛋白质的万金油
OPLS-AA 的全称是 Optimized Potentials for Liquid Simulations - All Atom。它把所有原子都显式处理，一个氢原子都不放过。这个力场最早是 Jorgensen 课题组搞出来的，初衷是模拟液体性质。
适用场景：
- 有机小分子溶液（比如苯、乙醇、水混合体系）
- 蛋白质-配体结合自由能计算
- 脂质双分子层模拟
- 高分子稀溶液（浓度别太高）
我个人经验：OPLS-AA 对液态密度和汽化焓的预测非常准。我做过一个聚乙二醇（PEG）水溶液的模拟，用 OPLS-AA 算出来的密度跟实验值只差 0.5%。但要注意，它对长链高分子的构象分布描述一般，尤其是当链长超过 50 个重复单元时。
避坑指南：我曾经用 OPLS-AA 模拟聚乳酸（PLA）的玻璃化转变温度，结果算出来比实验值高了将近 20 K。后来发现是二面角参数对酯基的描述不够准确。如果你做聚酯类材料，建议搭配专门的二面角修正参数。
4.2 COMPASS：凝聚态高分子的老朋友
COMPASS 是 Materials Studio 里的当家力场，全称是 Condensed-phase Optimized Molecular Potentials for Atomistic Simulation Studies。它属于第二代力场，参数拟合时同时考虑了气相和凝聚相数据。
适用场景：
- 高分子本体材料（熔体、玻璃态）
- 高分子共混物相容性预测
- 力学性能模拟（杨氏模量、泊松比）
- 界面体系（高分子-金属、高分子-无机填料）
| 性能 | COMPASS | OPLS-AA |
|---|---|---|
| 密度预测 | 优秀（误差 < 1%） | 良好（误差 1-3%） |
| 内聚能密度 | 优秀 | 良好 |
| 链构象 | 良好 | 一般 |
| 计算速度 | 中等 | 较慢 |
我的建议：如果你用 Materials Studio 做模拟，COMPASS 是首选。但如果你用 LAMMPS 或 GROMACS，COMPASS 的参数文件不太好找。我一般只在 MS 里用 COMPASS，转到开源软件就换别的力场。
4.3 PCFF：从聚合物到材料的全能选手
PCFF 是 Polymer Consistent Force Field 的缩写，它跟 COMPASS 算是「亲戚」——都源于 CFF 系列力场。但 PCFF 更早，参数覆盖范围也更广。
适用场景：
- 多种高分子材料（聚乙烯、聚丙烯、尼龙、聚碳酸酯）
- 高分子-溶剂相互作用
- 热力学性质（比热容、热膨胀系数）
- 交联高分子体系
我记得有一次做环氧树脂的交联模拟，用 PCFF 算出来的交联密度和实验值吻合得很好。但要注意，PCFF 对含氟高分子的描述不太行——氟原子的电负性太强，PCFF 的静电参数处理得不够精细。
4.4 GAFF：药物分子的标准配置
GAFF 是 General Amber Force Field 的简称，它是 Amber 软件包的一部分。这个力场的特点是参数化非常系统——几乎覆盖了所有有机官能团。
适用场景：
- 药物分子与靶点结合模拟
- 有机小分子晶体结构预测
- 核酸、糖类等生物分子
- 高分子短链（分子量 < 2000）
注意：GAFF 不适合长链高分子。我试过用 GAFF 模拟 100 个重复单元的聚乙烯，结果链的均方回转半径比实验值小了 30%。原因很简单——GAFF 的参数主要来自小分子，对长链的构象熵描述不足。
4.5 TraPPE-UA：粗粒化的效率之选
TraPPE-UA 是 Transferable Potentials for Phase Equilibria - United Atom 的缩写。它把 CH₃、CH₂、CH 等基团合并成一个「联合原子」，大大减少了计算量。
适用场景：
- 相平衡计算（气液平衡、液液平衡）
- 高分子熔体流变学
- 大尺度体系（超过 10 万个原子）
- 需要长时间模拟（微秒级）
为什么会有人用联合原子力场？你想想看，一个 100 个重复单元的聚乙烯链，全原子模型有 300 多个原子，而 TraPPE-UA 只有 100 个「珠子」。计算速度能快 3-5 倍。
我曾经踩过的坑：用 TraPPE-UA 模拟聚苯乙烯的玻璃化转变温度，结果比实验值低了 15 K。原因是联合原子模型丢失了苯环上氢原子的细节，导致链间 π-π 堆积作用被低估。如果你关注的是材料的力学性能，建议还是用全原子力场。
4.6 力场选择决策树
说了这么多，到底该怎么选？我画了一张图，帮你快速决策：
4.7 我的最终建议
选力场没有绝对的对错，关键看你的研究问题。我个人的经验是：
- 做生物相关：优先 OPLS-AA 或 GAFF
- 做高分子本体：COMPASS 或 PCFF 更靠谱
- 做相平衡或大体系：TraPPE-UA 是效率之王
- 不确定时：先用小体系跑几个力场对比一下，哪个跟实验吻合好就用哪个
最后说一句：力场参数不是万能的。我见过有人拿着 OPLS-AA 去算金属有机框架（MOF），结果一塌糊涂。记住，每个力场都有自己的「舒适区」，超出这个范围，结果就要打问号了。