计算化学公社

标题: 开启机器学习分子动力学计算，却仍一直在进行从头计算 [打印本页]
作者
Author: GRZ    时间: 2024-3-5 22:25
标题: 开启机器学习分子动力学计算，却仍一直在进行从头计算
采用vasp 6.4.2进行MLMD计算，温度设置为400K时进行计算比较顺利，在一定步数的从头计算之后会采用机器学习力场进行计算。
但是将温度设置为300K时，则会一直进行从头计算，直到1500步达到ML_MB的上限报错，但是在400K时，不到100步就开始使用机器学习力场进行预测。
各位大佬，小弟想请教一下为什么会出现更改温度后不进行机器学习力场预测这种问题？

以下是INCAR文件，300K和400K时仅对温度进行了修改。
  1. ALGO = Fast    # Fast
  2. EDIFF = 0.00001
  3. ENCUT = 400
  4. ISMEAR = 1  #0
  5. SIGMA = 0.1  #0.05
  6. ISPIN = 2
  7. MAGMOM = 80*5.0   12*1.0
  8. ISTART = 1
  9. LREAL = Auto
  10. LWAVE = False
  11. NELM = 300
  12. NWRITE=1
  13. PREC = Normal
  14. NPAR = 4

  16. #FOR MD
  17. NSW=60000
  18.    IBRION=0 #FOR MD
  19.    POTIM=0.5
  20.    ISIF = 2
  21.    MDALGO = 4
  22.    NHC_NCHAINS = 3
  23.    NBLOCK=1
  24.    KBLOCK=10
  25.    NELMIN=6
  26.    LCHARG   =.F.
  27.    TEBEG=300 # TEmperature at BEGining
  28.    TEEND=300 #TEmperature at the END
  29.    SMASS=1
  30.    RANDOM_SEED = 12345678   0    0

  32. #Machine learning paramters
  33.    ML_LMLFF = .TRUE.
  34.    ML_ISTART = 0
  35.    ML_LBASIS_DISCARD=.TRUE
复制代码
检查了一下300K时的ML_ABN文件，我已经跑了1000多步的从头计算，然而ML_ABN文件中却只有6个configurations，不知道这是什么原因。
此外对比了一下300K和400K时的ML_LOGFILE文件，在300K的计算在第二步时求解REGR用了50步，第8列的误差反而越来越大，不知道这个地方的不收敛该怎么解决。
  1. --------------------------------------------------------------------------------
  2. STATUS                  0 threshold  2      T      F         0         0
  3. --------------------------------------------------------------------------------
  4. STATUS                  1 critical   4      T      F         0         1
  5. LCONF                   1 Fe         0        80 H          0        12
  6. SPRSC                   1         1         1 Fe        80        25 H         12        12
  7. REGR                    1    1    1   2.03427090E-01   5.58399039E-01   1.00179913E-09   7.32063736E+01
  8. REGR                    1    1    2   6.54871393E-02   5.73960425E-01   3.19867932E-09   7.55316317E+01
  9. REGR                    1    1    3   1.39152913E-02   5.89413461E-01   1.54586841E-08   7.78722005E+01
  10. REGR                    1    1    4   1.40342867E-03   6.01123925E-01   1.56321403E-07   7.97374649E+01
  11. REGR                    1    1    5   1.62097401E-04   6.06333506E-01   1.36514974E-06   8.07505051E+01
  12. REGR                    1    1    6   6.10359679E-05   6.11595603E-01   3.65698596E-06   8.18132444E+01
  13. REGR                    1    1    7   4.05461288E-05   6.17966013E-01   5.56237121E-06   8.28869789E+01
  14. REGR                    1    1    8   3.35324841E-05   6.22468325E-01   6.77479725E-06   8.35985039E+01
  15. REGR                    1    1    9   3.06198798E-05   6.25024261E-01   7.44968943E-06   8.39952588E+01
  16. REGR                    1    1   10   2.93085577E-05   6.26357739E-01   7.79960796E-06   8.42008973E+01
  17. REGRF                   1    1   11   2.86959634E-05   6.27026457E-01   7.97461693E-06   8.43037312E+01    1.66079299E+15   2.74002942E+09
  18. NDESC                   1        24 Fe       452 H        452
  19. NDESC_SIC               1 Fe       226 H        226
  20. STDAB                   1   0.00000000E+00   4.49911633E-02   0.00000000E+00
  21. ERR                     1   1.33160442E-01   1.78565807E-02   1.50997419E+01
  22. CFE                     1   0.00000000E+00   0.00000000E+00   0.00000000E+00
  23. LASTE                   1   1.33160442E-01   1.78565807E-02   1.50997419E+01
  24. BEE                     1   1.00000000E-06   3.46410162E-02   2.00000000E-02   2.00000000E-03   0.00000000E+00   0.00000000E+00
  25. BEEF                    1   9.81859850E-04   1.20314674E-02   3.87184321E-03   2.00000000E-03   1.01482645E+00   6.42993156E-01
  26. --------------------------------------------------------------------------------
  27. STATUS                  2 critical   4      T      T         0         1
  28. LCONF                   2 Fe        25        89 H         12        24
  29. SPRSC                   2         2         2 Fe        89        30 H         24        24
  30. REGR                    2    1    1   2.19264203E-01   6.91866816E-01   2.00000000E-15   1.84764826E+02
  31. REGR                    2    1    2   4.52633698E-01   1.42824151E+00   2.00000000E-15   3.84527513E+02
  32. REGR                    2    1    3   2.25581875E-01   1.42360323E+00   4.00000000E-15   3.82699101E+02
  33. REGR                    2    1    4   3.87340298E-01   2.44442909E+00   4.00000000E-15   6.59668553E+02
  34. REGR                    2    1    5   1.93885338E-01   2.44714513E+00   8.00000000E-15   6.59026436E+02
  35. REGR                    2    1    6   3.32699860E-01   4.19920789E+00   8.00000000E-15   1.13557778E+03
  36. REGR                    2    1    7   1.66767586E-01   4.20975088E+00   1.60000000E-14   1.13576878E+03
  37. REGR                    2    1    8   2.61683214E-01   6.60572695E+00   1.60000000E-14   1.79025671E+03
  38. REGR                    2    1    9   1.31196056E-01   6.62362181E+00   3.20000000E-14   1.79094402E+03
  39. REGR                    2    1   10   1.83874186E-01   9.28315310E+00   3.20000000E-14   2.52142367E+03
  40. REGR                    2    1   11   9.21090040E-02   9.30051143E+00   6.40000000E-14   2.52173861E+03
  41. REGR                    2    1   12   1.19958298E-01   1.21125348E+01   6.40000000E-14   3.29801076E+03
  42. REGR                    2    1   13   6.01651977E-02   1.21501066E+01   1.28000000E-13   3.29797506E+03
  43. REGR                    2    1   14   7.66761443E-02   1.54844223E+01   1.28000000E-13   4.22365545E+03
  44. REGR                    2    1   15   3.84115826E-02   1.55141125E+01   2.56000000E-13   4.22369654E+03
  45. REGR                    2    1   16   4.97467327E-02   2.00922836E+01   2.56000000E-13   5.49335100E+03
  46. REGR                    2    1   17   2.49178403E-02   2.01282089E+01   5.12000000E-13   5.49341864E+03
  47. REGR                    2    1   18   3.30876444E-02   2.67276380E+01   5.12000000E-13   7.32211244E+03
  48. REGR                    2    1   19   1.65655977E-02   2.67628178E+01   1.02400000E-12   7.32205198E+03
  49. REGR                    2    1   20   2.19214242E-02   3.54155095E+01   1.02400000E-12   9.72041545E+03
  50. REGR                    2    1   21   1.09606810E-02   3.54154091E+01   2.04800000E-12   9.72038788E+03
  51. REGR                    2    1   22   1.39033002E-02   4.49234005E+01   2.04800000E-12   1.23554704E+04
  52. REGR                    2    1   23   6.96826942E-03   4.50307989E+01   4.09600000E-12   1.23554620E+04
  53. REGR                    2    1   24   8.44227457E-03   5.45562099E+01   4.09600000E-12   1.50168415E+04
  54. REGR                    2    1   25   4.23019788E-03   5.46733138E+01   8.19200000E-12   1.50168313E+04
  55. REGR                    2    1   26   5.10653651E-03   6.59995776E+01   8.19200000E-12   1.81935431E+04
  56. REGR                    2    1   27   2.55829413E-03   6.61294918E+01   1.63840000E-11   1.81935351E+04
  57. REGR                    2    1   28   3.25320448E-03   8.40922692E+01   1.63840000E-11   2.32222091E+04
  58. REGR                    2    1   29   1.62929344E-03   8.42313992E+01   3.27680000E-11   2.32222127E+04
  59. REGR                    2    1   30   2.28063791E-03   1.17904680E+02   3.27680000E-11   3.26237345E+04
  60. REGR                    2    1   31   1.14251044E-03   1.18131271E+02   6.55360000E-11   3.26237300E+04
  61. REGR                    2    1   32   1.78520365E-03   1.84583326E+02   6.55360000E-11   5.11659008E+04
  62. REGR                    2    1   33   8.93791241E-04   1.84829288E+02   1.31072000E-10   5.11659017E+04
  63. REGR                    2    1   34   1.49237737E-03   3.08612385E+02   1.31072000E-10   8.57242616E+04
  64. REGR                    2    1   35   7.47378398E-04   3.09104433E+02   2.62144000E-10   8.57242658E+04
  65. REGR                    2    1   36   1.21531657E-03   5.02636602E+02   2.62144000E-10   1.39868513E+05
  66. REGR                    2    1   37   6.09666730E-04   5.04297929E+02   5.24288000E-10   1.39868515E+05
  67. REGR                    2    1   38   9.02955017E-04   7.46897153E+02   5.24288000E-10   2.08062133E+05
  68. REGR                    2    1   39   4.51477509E-04   7.46897154E+02   1.04857600E-09   2.08062133E+05
  69. REGR                    2    1   40   6.14337274E-04   1.01632252E+03   1.04857600E-09   2.83827296E+05
  70. REGR                    2    1   41   3.07723880E-04   1.01815965E+03   2.09715200E-09   2.83827296E+05
  71. REGR                    2    1   42   4.03364217E-04   1.33460285E+03   2.09715200E-09   3.73091253E+05
  72. REGR                    2    1   43   2.02063289E-04   1.33712526E+03   4.19430400E-09   3.73091253E+05
  73. REGR                    2    1   44   2.64106839E-04   1.74768969E+03   4.19430400E-09   4.89073304E+05
  74. REGR                    2    1   45   1.32053420E-04   1.74768969E+03   8.38860800E-09   4.89073305E+05
  75. REGR                    2    1   46   1.72394156E-04   2.28158793E+03   8.38860800E-09   6.39555806E+05
  76. REGR                    2    1   47   8.61970781E-05   2.28158793E+03   1.67772160E-08   6.39555806E+05
  77. REGR                    2    1   48   1.13209669E-04   2.99659594E+03   1.67772160E-08   8.40990129E+05
  78. REGR                    2    1   49   5.66048345E-05   2.99659594E+03   3.35544320E-08   8.40990129E+05
  79. REGRF                   2    1   50   7.62449763E-05   4.03632285E+03   3.35544320E-08   1.13401357E+06   -1.34457656E+16   1.57770125E+15
  80. NDESC                   2        24 Fe       452 H        452
  81. NDESC_SIC               2 Fe       226 H        226
  82. STDAB                   2   4.45091206E-06   5.75450163E-02   4.03952569E-02
  83. ERR                     2   1.33160442E-01   1.78565807E-02   1.50997419E+01
  84. CFE                     2   5.46201079E-02   2.28753996E-02   4.01021682E-01
  85. LASTE                   2   1.33160442E-01   1.78565807E-02   1.50997419E+01
  86. BEE                     2   9.81858683E-04   1.19753653E-02   3.87251056E-03   2.00000000E-03   1.01472259E+00   6.43172560E-01
  87. BEEF                    2   1.00000000E-06   3.46410162E-02   2.00000000E-02   2.00000000E-03   0.00000000E+00   0.00000000E+00
  88. --------------------------------------------------------------------------------
  89. STATUS                  3 critical   4      T      F         0         1
  90. LCONF                   3 Fe        30       110 H         24        36
  91. SPRSC                   3         3         3 Fe       110        35 H         36        34
  92. REGR                    3    1    1   5.22685770E-04  -7.54835046E-04  -1.21455889E-13   1.04923914E+02
  93. REGR                    3    1    2   8.42665437E-08   3.36200129E-02   3.35544320E-08   3.25987270E+03
  94. REGR                    3    1    3   8.84127246E-05   7.05484480E+01   6.71088640E-08   2.95607793E+04
  95. REGR                    3    1    4   4.91276053E-05   7.84022058E+01   1.34217728E-07   3.29252129E+04
  96. REGR                    3    1    5   5.35272916E-05   8.54236167E+01   1.34217728E-07   3.59228173E+04
  97. REGRF                   3    1    6   5.35272916E-05   8.54236167E+01   1.34217728E-07   3.59228173E+04   -2.53531470E+16   1.18902997E+13
  98. NDESC                   3        24 Fe       452 H        452
  99. NDESC_SIC               3 Fe       226 H        226
  100. STDAB                   3   6.80049051E-05   7.69205421E-02   4.38224826E-02
  101. ERR                     3   2.61772686E-04   5.65812477E-01   1.70058187E+00
  102. CFE                     3   0.00000000E+00   0.00000000E+00   0.00000000E+00
  103. LASTE                   3   2.84116947E-04   5.70810154E-01   1.73283363E+00
  104. BEE                     3   1.00000000E-06   3.46410162E-02   2.00000000E-02   2.00000000E-03   0.00000000E+00   0.00000000E+00
  105. BEEF                    3   3.95522187E-06   2.29808927E-01   3.38149687E-02   2.00000000E-03   2.43673172E-01   1.71571998E-01
复制代码

在400K时，第二步的误差是逐渐下降的
  1. STATUS                  0 threshold  2      T      F         0         0
  2. --------------------------------------------------------------------------------
  3. STATUS                  1 critical   4      T      F         0         1
  4. LCONF                   1 Fe         0        80 H          0        12
  5. SPRSC                   1         1         1 Fe        80        25 H         12        12
  6. REGR                    1    1    1   2.03427111E-01   5.58398988E-01   1.00179880E-09   7.32063669E+01
  7. REGR                    1    1    2   6.54871649E-02   5.73960373E-01   3.19867735E-09   7.55316247E+01
  8. REGR                    1    1    3   1.39153029E-02   5.89413409E-01   1.54586677E-08   7.78721935E+01
  9. REGR                    1    1    4   1.40343056E-03   6.01123878E-01   1.56321159E-07   7.97374583E+01
  10. REGR                    1    1    5   1.62097575E-04   6.06333461E-01   1.36514799E-06   8.07504988E+01
  11. REGR                    1    1    6   6.10360115E-05   6.11595558E-01   3.65698258E-06   8.18132381E+01
  12. REGR                    1    1    7   4.05461508E-05   6.17965968E-01   5.56236704E-06   8.28869727E+01
  13. REGR                    1    1    8   3.35324992E-05   6.22468282E-01   6.77479283E-06   8.35984979E+01
  14. REGR                    1    1    9   3.06198917E-05   6.25024219E-01   7.44968503E-06   8.39952530E+01
  15. REGR                    1    1   10   2.93085681E-05   6.26357698E-01   7.79960364E-06   8.42008917E+01
  16. REGRF                   1    1   11   2.86959729E-05   6.27026417E-01   7.97461269E-06   8.43037257E+01    1.66074114E+15   2.74002979E+09
  17. NDESC                   1        24 Fe       452 H        452
  18. NDESC_SIC               1 Fe       226 H        226
  19. STDAB                   1   0.00000000E+00   4.49911636E-02   0.00000000E+00
  20. ERR                     1   1.33160414E-01   1.78565807E-02   1.50997415E+01
  21. CFE                     1   0.00000000E+00   0.00000000E+00   0.00000000E+00
  22. LASTE                   1   1.33160414E-01   1.78565807E-02   1.50997415E+01
  23. BEE                     1   1.00000000E-06   3.46410162E-02   2.00000000E-02   2.00000000E-03   0.00000000E+00   0.00000000E+00
  24. BEEF                    1   9.81859806E-04   1.20314673E-02   3.87184317E-03   2.00000000E-03   1.01482644E+00   6.42993152E-01
  25. --------------------------------------------------------------------------------
  26. STATUS                  2 critical   4      T      T         0         1
  27. LCONF                   2 Fe        25        89 H         12        24
  28. SPRSC                   2         2         2 Fe        89        31 H         24        24
  29. REGR                    2    1    1   4.09510137E+00   6.83152741E-01   2.00000000E-15   1.83261423E+02
  30. REGR                    2    1    2   1.78824428E+00   2.98318373E-01   2.00000000E-15   7.94947179E+01
  31. REGRF                   2    1    3   1.79236326E+00   2.99005511E-01   2.00000000E-15   7.96778236E+01   -1.36179849E+16   8.34109683E+13
  32. NDESC                   2        24 Fe       452 H        452
  33. NDESC_SIC               2 Fe       226 H        226
  34. STDAB                   2   1.96774589E-05   6.13894063E-02   4.83205647E-02
  35. ERR                     2   3.33592945E-06   2.83593686E-02   5.47773247E-02
  36. CFE                     2   6.23424347E-02   2.51345630E-02   4.65148645E-01
  37. LASTE                   2   3.33290111E-06   2.83506991E-02   5.78562090E-02
  38. BEE                     2   9.81858332E-04   1.19683389E-02   3.87305125E-03   2.00000000E-03   1.01470925E+00   6.43202425E-01
  39. BEEF                    2   1.00000000E-06   3.08795700E-02   6.41811580E-03   2.00000000E-03   2.73977324E-02   1.84484828E-02
复制代码









欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3