求助：请问反应力场弛豫时如何关闭键参数以避免发生反应呢

zd7777777

如题，在反应立场中进行氧原子轰击材料模拟，需要保证材料在氧原子撞击之前成分不变。氧原子持续撞击后进行产物分析

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

菓瀛

请问您解决了么，我目前也遇到这样的问题了

mengxiangidea

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

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

菓瀛

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

uranusshi

关闭键是这个意思吗，学习了