24小时热门版块排行榜    

查看: 905  |  回复: 0

fw7668172

捐助贵宾 (初入文坛)

[求助] 运行错误H原子种类不对 ,lj势语句错误

angle_coeff        1 2.60185 104.52
bond_coeff        1 10.84 0.9572

#angle_coeff        2  0 109.47
#bond_coeff        2  0  0.142
pair_style        lj/cut/tip4p/long 3 2 1 1 0.125 9 8
pair_coeff        1 1 0.0 0.0
pair_coeff        1 4 0.00324   3.19
pair_coeff        4 4 0.0070575 3.16435
pair_coeff        2 2 0.0 0.0
pair_coeff        1 2 0.0 0.0
pair_coeff        2 3 0.0 0.0
pair_coeff        2 4 0.00324   3.19
pair_coeff        3 4 0.0 0.0
pair_coeff        3 3 0.0 0.0
pair_coeff        1 3 0.0 0.0
kspace_style        pppm/tip4p 1.0e-5
group        g_A type 1
group        g_X type 2
group        g_C union g_A g_X
group        g_0 type 4
group        g_H type 3
group        g_W union g_0 g_H

#variable        field equal 0.04*cos(2*PI*20e12*1e-12*time)
#minimize        1e-5 1e-8 10000 100000
#reset_timestep        0

fix        1 g_C setforce 0.0 0.0 0.0
fix        2 g_W nvt temp 300 300 0.1
fix        3 g_W recenter INIT INIT NULL
fix        4 g_W shake 1e-5 500 0 b 1 a 1
#fix        5 g_W efield 0.0 0.0 v_field

velocity        g_W create 300.0 790316 dist gaussian
run        0
velocity        g_W scale 300.0
velocity        g_C set 0.0 0.0 0.0

timestep        0.001
compute 1 g_W temp
thermo        100
thermo_style        custom step c_1  epair  
#thermo_modify temp 1

#dump        1 g_W custom 200 c-h.lammpstrj type x y z
dump        1 all custom 200 move.lammpstrj mass id type x y z q element

run        200000
write_data data.*

#write_restart        restart.stable_singleC_water

#fix        5 g_W efield 0.0 0.0 0.08
#compute        1 g_W chunk/atom bin/3d x lower 1.5 y lower 1.5 z 0.0 1 nchunk every ids nfreq bound y 48.1425 49.6425 units box
#fix        5 g_W ave/chunk 100 5 10000 1 density/number ave running file angle.profile
#dump        1 g_W custom 500 angle.lammpstrj type x y z
#dump        2 all custom 200 single_C.lammpstrj id mol type x y z q



运行错误ERROR on proc 0:TIP4P hydrogen has incorrect atom type
模型是水滴和石墨炔
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 fw7668172 的主题更新
信息提示
请填处理意见