求助：lammps在优化分子晶胞时发生断键，如何解决？

江湖采花贼 · 发表于 Post on 2021-12-8 21:27:29

最近在研究lammps优化分子晶体的晶胞，但是在优化结束后发现分子断键，分裂成原子核一些片段，下面是优化的命令：units          real
newton       on
dimension    3
boundary       p p p
atom_style  charge
read_data expand.data
mass 1 1.00794  # H
mass 2 12.0107  # C
pair_style    reax/c NULL
pair_coeff    * * ffield.reax C H
neigh_modify every 1 delay 0 check yes
fix  1 all qeq/reax 1 0.0 10.0 1.0e-6 reax/c
min_style cg
minimize       1.0e-8 1.0e-8 15000 30000
dump  1 all xyz 100 out1.xyz
dump_modify 1 append yes
timestep    0.05
velocity all create 300 12345
fix 2 all npt temp 20.0 5.0 10.0 aniso 1.0 1.0 1000.0
run 10000
reset_timestep 0
min_style    cg
dump 2 all xyz 100 out2.xyz
minimize       1.0e-8 1.0e-8 30000 30000

求大神们指点，为何会出现这种情况？如何改正？跪谢！！！！

fineren · 发表于 Post on 2021-12-9 11:23:50

别用反应力场

江湖采花贼 · 发表于 Post on 2021-12-9 17:29:37

fineren 发表于 2021-12-9 11:23
别用反应力场

能说说为啥不，谢谢

fineren · 发表于 Post on 2021-12-10 09:08:28

优化晶胞没必要用反应力场，后者对你的体系不合适的话就会断键，非反应力场不会断键

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[Lammps] 求助：lammps在优化分子晶胞时发生断键，如何解决？