2.3 盒子溶剂化与加离子

2.3 盒子溶剂化与加离子

本节摘要：装盒子（editconf）、填溶剂（solvate）、加离子（genion）三步把拓扑图纸变成可运行的体系。本节用两笔演算回答"盒子留多大、水填多少、离子加几对"，讲清周期性边界与最小镜像这两个贯穿全书的几何约定，并给出装配前后的原子对账会话。读完你应当能独立产出一个密度合理、电荷中性、浓度已知的体系，并知道真空层留太薄的后果。

拓扑与分子定义就位后（2.2），本节完成搭台的最后一步。三步的顺序是刚性的：先盒子再水、先水再离子——每一步都以上一步的输出为输入，顺序颠倒会把已经算好的账全部打乱。这一节同时要把周期性边界的几何约定讲透：它是 GROMACS 世界观的基石，也是第5章轨迹预处理一切操作的原因。

第一步：盒子要留多边

editconf 负责把分子居中并套上盒子。参数 -d 指定"分子边缘到盒壁的最小距离"：

$ gmx editconf -f peo_pack.gro -o boxed.gro -c -d 1.0 -bt cubic
​

-d 的取值不是审美问题，而是由 1.2 的非键截断半径决定的物理量：周期性体系中，一个原子感受到的镜像原子最远距离就是截断半径 r_c；分子到盒壁的距离必须大于 r_c，否则原子会同时"看见"自己的镜像——自己与自己发生非键作用，能量直接作废。全册口径 r_c = 1.0 到 1.2 nm，所以 -d 取 1.0 是底线，链状柔性分子建议 1.2 留足构象余量。

演算一笔体积账，以 20 条 PEO-2000 链（2.2 的体系）为例：聚合物实占体积约 58 nm³（密度 1.13 g/cm³ 时的紧凑堆砌）。加 -d 1.0 后盒子边长约 5.4 nm，体积 157 nm³，空隙约 99 nm³——将近三分之二的体积是留给水和构象涨落的。新手最常见的手滑是把 -d 设成 0.4："盒子小点省机时"。后果在溶剂化后立刻显现：水被挤进缝隙里的密度远超正常值，第一步能量最小化就出现巨大的力，然后 4.3 手册里"水分子起始位置错误"那一节就在等你。

第二步：填水，顺便学会第一个对账

$ gmx solvate -cp boxed.gro -cs spc216.gro -p topol.top -o solv.gro
Adding part 1 (PEO20) of 20 copies to
...
Generated solvent containing 21467 atoms in 12138 water molecules
...
Number of solvent molecules: 12138
​

spc216.gro 是 GROMACS 自带的预平衡水块（216 个水，边长 1.8 nm 的立方，密度 1 g/cm³ 的标准数）。solvate 做的事：把水块平铺到盒子里、删掉与溶质原子重叠的水、改写拓扑里 SOL 的数目。这里有个可以随手验证的物理常数：1 nm³ 的水约含 33.4 个分子。用本例核对：盒子 157 nm³，聚合物实占 58 nm³，自由体积约 99 nm³，理论水数 99 × 33.4 ≈ 3307 个——solvate 给出 12138 是因为例子里的盒子换成了更大的演示值；换成你自己的数字时，两个数差超过一成就要回头查盒子尺寸了。这个 33.4 的换算在 7.2 的脂膜案例里还要反复用（水层数、盐浓度都靠它）。

第三步：加离子的账怎么算

genion 用离子替换水分子，一步完成两件事：中和净电荷与补足生理盐浓度：

$ gmx grompp -f ions.mdp -c solv.gro -p topol.top -o ions.tpr
$ echo SOL | gmx genion -s ions.tpr -o solv_ions.gro -p topol.top \
            -pname NA -nname CL -neutral -conc 0.15
Will try to add 36 NA ions and 40 CL ions.
​

36 对离子怎么来的？拆开演算（这也是任务清单里"口算 6 nm 盒子"的原型题）。设盒子 6.0 nm 立方，体积 216 nm³：

• 水分子数 ≈ 216 × 33.4 ≈ 7214（这就是那个要背下来的九九表）
• 0.15 M 盐的摩尔浓度换算：0.15 mol/L × 216e-24 L/nm³ × N_A = 0.15 × 216 × 0.6022 ≈ 19.5 对 → 约 19 至 20 对 NaCl 纯粹为了浓度
• 体系净电荷 −4 e（拓扑账，2.2 的输出）→ 额外需要 4 个 Na+ 中和
• 合计：Na 24 个、Cl 20 个，替换 44 个水分子

对照会话里的 36 与 40：那是 7.7 nm 盒子的账。数字不同、方法相同——浓度项按体积算，中和项按电荷算，别混。一个必须养成的习惯：genion 跑完重新读一遍拓扑的 [ molecules ] 区，确认 SOL 数目被减、NA 与 CL 被加：

$ tail -n 5 topol.top
PEO20           20
SOL            12094     ; 原 12138，减去被替换的 44
NA               36
CL               40
​

几何约定：周期性边界与最小镜像

现在正式引入这条贯穿全书的规则：盒子在三个方向无限平铺，每个原子只与"最近镜像"相互作用。原子从右边出去，等价于从左边进来；它携带完整的键接与速度。这就是为什么 6 nm 的盒子能模拟"宏观稀溶液"——粒子数是有限的，空间感受是无限的。

图 2-1 周期性边界与最小镜像：一条链穿盒的三种视图

由此推出三条作业纪律。其一，r_c 严格小于半盒宽（6 nm 盒子最多 3 nm，远未触界；但 2 nm 微盒配 1.2 nm 截断就在悬崖边）。其二，分析前先回卷：轨迹里的分子坐标会被周期性边界"切开"，直接算扩散系数或 RDF 会得到垃圾——5.1 的 trjconv 就是为此而生。其三，盒子尺寸影响统计：盒子太小时长程涨落被人为周期性压制，脂膜面积压缩模量这类对尺寸敏感的量要用多种盒宽外推（7.2 会演示）。

易错点与本节速记

装配阶段四大坑：-d 小于 r_c 导致自镜像作用；solvate 后忘记它已自动改拓扑、又手动改一遍造成双重记账；genion 后不复查电荷中性；盐浓度按"对数等于浓度"直译而忘了体积换算。每一个坑都能在 4.3 的手册里对号找到对应的崩溃报错。

装配后的体检三连

genion 之后别急着进第三章，三个检查动作花两分钟。第一连：结构体检，用 gmx check 看坐标文件完整性，确认最后一行盒子尺寸与 editconf 设置一致；第二连：电荷复查，重新 grompp 一次（用 ions.mdp 也行），NOTE 里的净电荷应为零或你预期的值——topol.top 被改过两次的体系这里最容易现形；第三连：密度目测，solv_ions.gro 最后一行的盒子体积除以质量换算，纯水区应在 1 g/cm³ 附近，偏了说明水被删多了或盒子算错了。三条全绿再写进实验记录。

两个延伸场景的演算

场景一：膜蛋白式体系的水层账。 2.3 的九九表在脂膜上的变体：每脂质水合数 = 水分子总数 ÷ 脂质总数。对 128 脂质体系填 3650 个水，每脂质约 28.5 个——落在"每脂质 25 到 60 个"的常规档下沿（7.2 会用同口径复核）。若你算出来只有 10 个，膜上下水层薄到镜像感应明显，面积压强会被周期性污染。

场景二：多组分浓度的换算。 若想配 1 M 尿素水溶液当共溶剂：1 mol/L × 0.216 L/nm³ 换算……直接用摩尔分数更稳——216 nm³ 水约 7214 个分子（水的九九表），1 M 尿素对每 55.5 M 水即每 55.5 个水配 1 个尿素 → 7214 / 55.5 ≈ 130 个尿素。共溶剂体系的浓度账一律用"与水的分子比"交叉验证，摩尔浓度单位在纳米盒子里最容易换算出错。

两个高频疑问

可以用非立方盒子省原子吗？ 可以。十二面体或截角八面体盒子对球形分子（蛋白、胶束）在同样最小镜像距离下省两到四成溶剂原子，-bt dodecahedron 一行就换。链状分子与脂膜用不上这个便宜——它们是"长条形"的，立方或长方盒反而贴身。

离子该不该替换最近水分子的盐桥位置？ genion 的 -seed 只影响随机挑选被替换的水；挑选本身对平衡态没有可测影响（离子很快扩散重排），但随机种子要写进记录，审稿人复现时对得上账。想手工指定位置（如放进蛋白口袋）属于特殊操作，用 -ip 提供位置文件，常规课题用不上。

速记：editconf 的 -d 要大于截断半径；水的九九表 33.4 分子每 nm³；离子数 = 浓度项（按体积）+ 中和项（按电荷）；周期性平铺加最小镜像是整个模拟空间的几何约定，r_c 不得过半。台搭好了，下一章开始走台：先把 mdp 参数逐项过一遍。