刘建川・Gromacs 实战｜如何手动添加力场中不支持的离子参数
大家好，我是刘建川，14 年深耕 GROMACS 分子模拟领域，始终专注分子模拟实战。我将拆解分子模拟实操技巧，帮你理清操作步骤、避开报错误区，真正做到学以致用。
在分子模拟体系中，经常会用到各类离子。如果选用的力场没有内置对应离子参数，GROMACS 会直接抛出报错： Atomtype ** not found
遇到这类问题，需要我们手动把缺失的离子参数补充到力场配置文件中。下面以 amber99 s b 力场为例，完整讲解添加方法，其他主流力场均可参照同样逻辑修改。
一、力场核心配置文件说明
GROMACS 调用任意力场，都会依赖两个关键文件：
ffbonded.itp 存放键合相关参数，包含键长、键角、二面角力常数等。由于离子不存在成键作用，添加离子无需修改此文件。
ffnonbonded.itp 存放非键相互作用参数，核心为Lennard-Jones势参数，是我们添加离子必须修改的文件。
二、离子参数获取与单位适配
Amber 力场采用标准 LJ 势函数形式，需要提前获取离子对应的 sigma 和 epsilon 参数。 单位要求：
sigma 单位：nm epsilon 单位：kJ/mol
不同文献、不同力场的LJ势表达式存在差异，拿到原始参数后，必须根据公式做单位和形式换算，匹配当前所用力场的数学形式，直接套用极易造成模拟报错和能量异常。
三、具体实操修改步骤
1. 编辑 ffnonbonded.itp 文件
按照文件内现有原子类型格式，在 atomtypes 字段下方，新增自定义离子的原子序号、质量、电荷、sigma、epsilon 等完整参数，严格对齐原有排版格式。
2. 编辑 ions.itp 文件
修改完非键参数文件后，还需要在 ions.itp 中补充离子分子定义，这样在用GROMACS命令添加离子时，程序才能正常识别。
按照文件内现有离子格式逐条仿写，重点注意： [ moleculetype ] 下方的离子名称，必须和 ffnonbonded.itp 中自定义的原子类型名称保持完全一致，否则依然会提示原子类型找不到。
掌握这套方法，就能自行给任意力场扩充缺失的金属离子、卤素离子等参数，彻底解决模拟中原子类型缺失的报错问题。
