24小时热门版块排行榜

>论坛更新日志 (2879)
>虫友互识 (968)
>休闲灌水 (501)
>导师招生 (62)
>文献求助 (29)
>论文道贺祈福 (24)
>基金申请 (24)
>微米和纳米 (23)
>论文投稿 (22)
>博后之家 (19)
>考博 (19)
>硕博家园 (18)
>考研 (17)
>教师之家 (16)
>公派出国 (13)
>找工作 (11)

返回列表

xufang0122

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 158.3
帖子: 12
在线: 9.4小时
虫号: 1265722
注册: 2011-04-14
专业: 工程热物理相关交叉领域

[求助] reaxfff反应力场的问题

我采用reaxff反应力场进行模拟，在计算一半就发生中断了。之前没有加入fix    2 all reax/c/species 10 100 1000 species.out 这条命令的时候可以正常运行，请问这条命令有什么错误吗？下面是我的In文件，请高手指点。
dimension    3
boundary    p p p
units          real

atom_style charge
read_restart testart.equil

pair_style reax/c lmp_control
pair_coeff * * ffield.reax.1 C O H N S
neighbor 2.0 bin
neigh_modify every 1 delay 0 check yes
fix          1 all qeq/reax 1 0.0 10.0 1e-6 reax/c
fix          2 all reax/c/species 10 100 1000 species.out

fix 3 all nvt temp 300.0 300.0 100.0
timestep       0.25

thermo       5000
thermo_style custom step temp press density pe ke etotal vol

dump          1 all custom 5000 2coal.lammpstrj id type x y z
dump          2 all custom 5000 2.dump id element x y z
dump_modify    2 element C O H N S

run 100000
unfix          3
undump       1
undump       2
write_restart testart.1equil
run          0
谢谢！

回复此楼