lammps计算聚合物例子_LAMMPS-GCMC学习（2）

fix gcmc

一、 语法

ID, group-ID：fix编号和group的编号

gcmc：不解释

N：每N步来一个这个fix

X：每N步进行约X次GCMC交换（应该是交换原子的意思吧，未确认）

M：每N步进行约M个MC试移动

type：插入原子的类型（gcmc应该会在随机位置插入溶质原子）

seed：随机种子数（正整数）

T：大热浴系统的温度（温度单位）

mu：大系统的化学势（能量单位），看文献应该是0.*eV这么大，自己试

displace：MC移动原子距离的最大值（长度单位）

keyword：关键词看情况添加，以下只介绍模拟金属可能会用到的

a) mcmoves：Patomtrans=MC moves中原子平移的比例

b) rigid：fix-ID=fix rigid/small命令的ID

c) shake：fix-ID=fix shake命令的ID

d) region：region-ID，进行GCMC交换原子和MC移动原子的区域

e) pressure：大系统的压强（压强单位）

f) fugacity_coeff：大系统的逸度系数（无单位）

g) full_energy：进行GCMC交换和MC移动原子时，计算系统的总能量（能量单位）

h) charge：插入原子的电荷量（电荷单位）

i) group：group-ID=要把原子插入到哪组原子中（字符串）

j) grouptype：type（插入原子的类型（正整数）） group_ID（要把原子插入哪个组中）

k) tfac_insert：按该比例提高/降低插入原子的温度（无单位）

l) overlap_cutoff：删除重叠原子的最大cufoff 距离

二、 几个例子

fix 2 gas gcmc 10 1000 1000 2 29494 298.0 -0.5 0.01

fix 3 water gcmc 10 100 100 0 3456 3.0 -2.5 0.1 mol my_one_water maxangle 180 full_energy

fix 4 my_gas gcmc 1 10 10 1 12345 300.0 -12.5 1.0 region disk

三、 GCMC全面描述

GCMC࿰