24小时热门版块排行榜    

查看: 904  |  回复: 2

xufang0122

新虫 (初入文坛)

[求助] 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
谢谢!
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

双双liu1211

新虫 (初入文坛)

楼主,你的反应力场是哪里找的呀?

发自小木虫Android客户端
2楼2018-04-23 14:26:59
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

aiwodejia

木虫 (小有名气)

能和我说一下你得到的数据都是怎么处理的么?

发自小木虫Android客户端
3楼2018-05-19 20:24:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 xufang0122 的主题更新
信息提示
请填处理意见