团簇融化

#原子数量

variable npart  equal 500

#LJ单位温度

variable tempstart   equal 0.5

variable tempstop    equal 0.5

#模型基本设置

units    lj

dimension 2

atom_style  atomic

boundary        p p p

neighbor        0.5   bin

neigh_modify    every 1 delay 0 check yes

#模拟盒子设置

region box block -20 20 -20 20 -0.1 0.1

create_box 1 box

#2维计算设置

fix 2d  all enforce2d

#在box中央设置box2

region box2 block -10 10 -10 10 -0.1 0.1

#在box2中随机生成原子

create_atoms 1 random ${npart} 324523 box2

#摩尔质量

mass   *  1

#LJ力场设置

pair_style lj/cut 2.5

pair_coeff * * 1.0 1.0

#能量最小化

minimize 1e-4 1e-4 1000 1000

reset_timestep 0

#保存原子轨迹

dump        img all  atom  1000 dump.xyz

#nve系综，langevin控温

fix integrator all nve

fix dynamics   all langevin ${tempstart} ${tempstop} 1.0 2352325

#设置屏幕输出热力学信息

thermo_style custom step temp  ke pe

thermo    100

#模拟步长设置

timestep  0.01

#运行20000步

run    20000