24小时热门版块排行榜    

查看: 1554  |  回复: 0

hjy2940

铜虫 (小有名气)

[求助] lammps求小分子密度

最近刚学lammps, 准备先算个小分子密度来verify,可能是对计算流程不是很熟悉,自己写的code有问题,好像需要再加个fix或者其他修改,大家能给点意见吗

#initialization
units real
atom_style full
dimension 3
boundary p p p


#force field information
pair_style lj/cut/coul/long 12
bond_style harmonic
angle_style harmonic
dihedral_style opls
kspace_style ewald 1e-4

read_data data.mol

velocity all create 1000.0 825577 dist gaussian

# settings
neighbor 2.0 bin
neigh_modify delay 5

#thermo setting
thermo_style multi

timestep 2.0

#minimize
min_style sd
minimize 1.0e-10 1.0e-10 5000 5000

variable mass equal mass(all)/6.022e23
variable vol equal vol*1.0e-24
variable density equal v_mass/v_vol

fix density all ave/time 100.0 1.0 100.0 v_density file density.mol

thermo 100
thermo_style custom step temp press v_density

run 1000
回复此楼
材料
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

智能机器人

Robot (super robot)

我们都爱小木虫

相关版块跳转 我要订阅楼主 hjy2940 的主题更新
信息提示
请填处理意见