计算化学公社

标题: 求助：请问反应力场弛豫时如何关闭键参数以避免发生反应呢 [打印本页]

作者
Author: zd7777777 时间: 2021-12-16 10:13
标题: 求助：请问反应力场弛豫时如何关闭键参数以避免发生反应呢
如题，在反应立场中进行氧原子轰击材料模拟，需要保证材料在氧原子撞击之前成分不变。氧原子持续撞击后进行产物分析

我按照文献温度进行弛豫时，材料本身就发生变化了。所以怎么保证材料在被撞击之前（弛豫后）保持成分不变呢

作者
Author: 菓瀛 时间: 2022-9-14 19:21
请问您解决了么，我目前也遇到这样的问题了

作者
Author: mengxiangidea 时间: 2022-9-15 21:58

菓瀛发表于 2022-9-14 19:21
请问您解决了么，我目前也遇到这样的问题了

别人都去年的贴了，又在这遇见你了，具体你模型什么样，需要组合，还是很混合

作者
Author: 菓瀛 时间: 2022-9-17 10:06
本帖最后由菓瀛于 2022-9-17 10:08 编辑

mengxiangidea 发表于 2022-9-15 21:58
别人都去年的贴了，又在这遇见你了，具体你模型什么样，需要组合，还是很混合

我是在reaxff力场下模拟甲烷、氨气、氧气混合燃烧，在能量最小化和弛豫阶段至少有三分之一氨气分子会分解并与氧气反应。所以我尝试在力场文件中关闭N-O、O-H键，但是就一直出现"fix qeq/reax CG convergence failed after 200 iterations at step 1(src/REAXFF/fix_qeq_reaxff.cpp:771)"的警告，也不计算。
这里红色是我修改的部分
C 1.3742 4.0000  12.0000 1.9684 0.1723 0.8712 1.2385 4.0000
   9.4606 2.1346 4.0000  31.0823  79.5548 5.7254 6.9235 0.0000
   1.2104 0.0000 183.7012 5.7419  33.3951  11.9957 0.8563 0.0000
   -2.8983 2.5000 1.0564 4.0000 2.9663 0.0000 0.0000 0.0000
H 0.6867 1.0000 1.0080 1.3525 0.0616 0.8910  -0.1000 1.0000
   9.3858 5.0013 1.0000 0.0000 121.1250 3.8446  10.0839 1.0000
   -0.1000 0.0000  58.4228 3.8461 3.2540 1.0000 1.0698 0.0000
-15.7683 2.1504 1.0338 1.0000 2.8793 0.0000 0.0000 0.0000
O 1.3142 2.0000  15.9990 1.9741 0.0880 0.8712 1.1139 6.0000
   10.2186 7.7719 4.0000  29.5271 116.0768 8.5000 7.1412 2.0000
   0.9909  14.9473  69.2812 9.1371 1.6258 0.1863 0.9745 0.0000
   -3.5965 2.5000 1.0493 4.0000 2.9225 0.0000 0.0000 0.0000
N 1.2450 3.0000  14.0000 1.9951 0.1088 1.0512 1.1911 5.0000
   9.9303 7.8431 4.0000  32.4758 100.0000 6.7768 6.8035 2.0000
   1.0636 0.1045 128.0119 2.1604 2.9464 2.5181 0.9745 0.0000
   -4.0959 2.0047 1.0183 4.0000 2.8793 0.0000 0.0000 0.0000
10    ! Nr of bonds; Edis1;LPpen;n.u.;pbe1;pbo5;13corr;pbo6
                     pbe2;pbo3;pbo4;Etrip;pbo1;pbo2;ovcorr
  1  1 141.9346 113.4487  67.6027 0.1554  -0.3045 1.0000  30.4515 0.4283
      0.0801  -0.2113 8.5395 1.0000  -0.0933 6.6967 1.0000 0.0000
  1  2 163.6889 0.0000 0.0000  -0.4525 0.0000 1.0000 6.0000 0.5921
      12.1053 1.0000 0.0000 1.0000  -0.0097 8.6351 0.0000 0.0000
  2  2 169.8421 0.0000 0.0000  -0.3591 0.0000 1.0000 6.0000 0.7503
      9.3119 1.0000 0.0000 1.0000  -0.0169 5.9406 0.0000 0.0000
  1  3 164.0476 117.4881  72.1261  -0.6031  -0.1795 1.0000  14.9755 0.5413
      1.2626  -0.3063 7.0000 1.0000  -0.1588 4.5000 0.0000 0.0000
  3  3 110.4748 155.6441  40.0000 0.1150  -0.1054 1.0000  28.5221 0.2000
      0.9590  -0.2635 8.5715 1.0000  -0.1007 6.8548 1.0000 0.0000
  1  4 130.7147 175.2276  97.2523  -0.0368  -0.4942 1.0000  26.7545 0.5133
      0.3296  -0.3653 7.0000 1.0000  -0.1171 5.1025 1.0000 0.0000
3  4  0.0000  0.0000 0.0000 0.0000  0.0000 0.0000 0.0000 0.0000
      0.0000  0.0000 0.0000 0.0000  0.0000 0.0000 0.0000 0.0000
  4  4 157.7518  67.1322 160.9732  -0.5869  -0.1824 1.0000  12.0000 0.7136
      0.8204  -0.1657  10.6490 1.0000  -0.0967 4.5976 1.0000 0.0000
2  3  0.0000 0.0000 0.0000  0.0000 0.0000 0.0000 0.0000 0.0000
      0.0000 0.0000 0.0000 0.0000  0.0000 0.0000 0.0000 0.0000
  2  4  212.1772 0.0000 0.0000  -0.3585 0.0000 1.0000 6.0000 0.3316
      10.4316 1.0000 0.0000 1.0000  -0.0658 6.4545 0.0000 0.0000
下面是我的in文件
units                   real
boundary          p p p

atom_style          charge
read_data          NC150.1.data

pair_style                   reax/c lmp_control
pair_coeff                   * * ffield.reax(1).chon C H N O

fix                      q all qeq/reax 1 0.0 10.0 1e-6 param.qeq
neighbor             2 bin
neigh_modify    every 10 delay 0 check no
timestep             0.25

minimize             1e-8 1e-8 10000 10000
reset_timestep    0

velocity             all create 298.0 9898

thermo                1000
thermo_style       custom step press pe temp ke
dump                1 all xyz 100 NC.xyz
fix                      1 all npt temp 298.0 298.0 100  iso 0.0 0.0 250
run                      30000
unfix                   1
undump             1
reset_timestep    0

thermo             1000
thermo_style    custom step press pe temp ke
dump                1 all xyz 100 NC.xyz
fix                      1 all nvt temp 298.0 298.0 100
run                      30000
unfix                   1
undump             1
reset_timestep    0

pair_style                reax/c lmp_control
pair_coeff                * * ffield.reax(2).chon C H N O

thermo             1000
thermo_style    custom step time press pe temp ke
fix                   1 all nve
fix                   2 all qeq/reax 1 0.0 10.0 1e-6 param.qeq
fix                   3 all temp/berendsen 298.0 2500.0 100
fix                   4 all reaxff/species 10 10 100 species_NC.out  element C H N O
fix                   5 all reaxff/bonds 100 bonds.reaxc
dump             1 all custom 100 dump.reac id type x y z
dump             2 all atom 100 NC.lammpstrj
run                4000000

作者
Author: uranusshi 时间: 2023-2-10 15:36
关闭键是这个意思吗，学习了

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