计算化学公社

标题: 求助：lammps在优化分子晶胞时发生断键，如何解决？ [打印本页]

作者
Author: 江湖采花贼 时间: 2021-12-8 21:27
标题: 求助：lammps在优化分子晶胞时发生断键，如何解决？
最近在研究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

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

作者
Author: fineren 时间: 2021-12-9 11:23
别用反应力场

作者
Author: 江湖采花贼 时间: 2021-12-9 17:29

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

能说说为啥不，谢谢

作者
Author: fineren 时间: 2021-12-10 09:08
优化晶胞没必要用反应力场，后者对你的体系不合适的话就会断键，非反应力场不会断键

欢迎光临计算化学公社 (http://bbs.keinsci.com/)