墙的设置

nwall (0)

设置为1时，z=0 处有一面墙，设置为2时，z=z-box 处也有一面墙。

墙只能与 pbc=xy 一起使用。当设置为2时，可以使用压力耦合和 Ewald 求和（通常最好使用 x/y 可压缩性设置为 0 的半各向同性压力耦合，否则表面积将发生变化）。

墙通过可选的 wall-atomtype 与系统的其余部分交互。

自动添加能量组 wall0 和 wall1（对于 nwall=2），以监控能量组与每个墙的相互作用。

将禁用z方向上的质心运动移除。

wall-atomtype

每个墙在力场中的原子类型名称。例如通过在拓扑中定义一种特殊的墙原子类型及其自身的组合规则，将允许独立调整每个原子类型与墙的交互。

此处输入原子类型个数对应nwall，e.g. nwall = 2; wall-atomtype = opls_136 opls_136

wall-type

9-3

对墙后体积进行积分的LJ势: 9–3势

10-4

对墙面进行积分的LJ势: 10–4势

12-6

直接的LJ势, 由与墙的Z距离决定

table

用户定义的势, 根据与墙的Z距离进行索引, 以类似于energygrp-table的选项读入, 其中的第一个名称为“正常”能量组, 第二名称为wall0或wall1, 只使用表中的色散和排斥列.

wall-r-linpot (-1) [nm]

与墙的距离在此值以下时, 势能线性连续, 因此力为常数. 当一些原子超过墙时, 将此选项设置为正值对平衡尤其有用. 当此值 ≤0(对wall-type=table则是<0)时, 原子超过墙后会产生致命错误.

wall-density [nm-3] / [nm-2]

每面墙的原子数密度, 适用于类型为9–3和10–4的墙.

若墙类型一致，此处输入原子数密度个数也对应nwall，e.g. nwall = 2; wall-density = 64 64

wall-ewald-zfac (3)

第三个盒矢量的缩放因子, 仅用于Ewald加和, 最小值为2. Ewald加和只能与nwall=2联用, 并需要使用ewald-geometry=3dc. 盒子中真空层的作用是降低周期性映象之间不合实际的库仑相互作用.

总结

nwall = 2

wall_type = 10-4

wall_density = 5 5

wall_atomtype = opls_136  opls_136

wall_r_linpot = -1

wall_ewald_zfac = 3

ewald_geometry=3dc

pbc=xy