气体扩散

#气体粒子数量

variable npart  equal 100

#模拟基本参数设置

units    lj

dimension 3

atom_style  atomic

boundary        p p p

neighbor        0.5   bin

neigh_modify    every 1 delay 0 check yes

#三维box尺寸

region box block -10 10 -10 10 -10 10

#生成box

create_box 1 box

#随机产生100个气体原子

create_atoms 1 random ${npart} 324523 box

#设置摩尔质量1

mass   1  1

#力场设置，soft势

pair_style soft  1.0

pair_coeff * * 10.0

#温度初始化

velocity all create 2.0  34234123 dist gaussian

#能量最小化

minimize 1e-4 1e-4 1000 1000

reset_timestep 0

#保存轨迹文件

dump        img all  atom  1000 gas.xyz

#设置nve系综

fix 1 all nve

#热力学输出

thermo_style custom step temp  ke pe

thermo    100

#模拟步长

timestep  0.001

#运行100000步

run    100000