计算化学公社

标题: 求助高斯优化离子液体与硝酸铕结合后期出现波动不收敛的问题 [打印本页]
作者
Author: soda123    时间: 2022-9-27 09:32
标题: 求助高斯优化离子液体与硝酸铕结合后期出现波动不收敛的问题
小弟最近刚接触计算离子液体方面，在结构优化一个离子液体与硝酸铕结合时高斯总是不收敛，已经尝试更换了溶剂模型，但仍然在后期优化时上下波动巨大，望论坛内老师指导 下面附上我的输入文件和高斯优化能量图
  1. %NprocShared=40
  2. %mem=160000mb
  3. %chk=sucesseun1888dehp.chk
  4. # opt=(cart,maxstep=5,notrust) freq genecp int=fine scrf=(solvent=n-dodecane,read) pbe1pbe
  5. scf=(maxcycle=520,vshift,novaracc,xqc)

  7. n1888

  9. 0 7
  10. N                 -4.94050400   -0.40444400   -1.56212800
  11. O                 -1.24281700   -1.25458500   -1.39469900
  12. C                 -1.08410800   -2.22837000   -2.20938600
  13. O                 -1.96287300   -2.80669500   -2.84890500
  14. C                  0.35647500   -2.69598800   -2.42007700
  15. O                  1.22127700   -1.93083900   -1.59783700
  16. C                  2.58590600   -2.08912800   -1.87464400
  17. C                  3.32325400   -1.11200500   -0.97223800
  18. O                  2.65880900   -0.35264300   -0.24511100
  19. N                  4.65755300   -1.07814900   -0.99179400
  20. C                 -4.92096100   -1.33610900   -2.72572600
  21. C                 -4.53875200   -1.14358800   -0.29025800
  22. C                 -5.19740500   -2.50920000   -0.04182600
  23. C                 -4.31260600   -3.72572800   -0.33521200
  24. C                 -3.16252600   -3.96387500    0.64590700
  25. C                 -3.59080700   -4.18470700    2.10763400
  26. C                 -3.38213300   -2.97036300    3.01154000
  27. C                 -3.84688800   -3.19430100    4.44689900
  28. C                 -3.86901000    0.67028000   -1.76731000
  29. C                 -4.03834000    1.64359200   -2.93502700
  30. C                 -3.58872900   -1.99111700    5.34611800
  31. C                 -4.84836400    2.93193400   -2.72551400
  32. C                 -4.27764200    3.91887000   -1.70258400
  33. C                 -2.95824700    4.56126400   -2.13848600
  34. C                 -2.53906200    5.82023300   -1.37174200
  35. C                 -2.03165600    5.63832600    0.06287100
  36. C                 -3.09983600    5.39589900    1.12596900
  37. C                 -6.27486200    0.28252200   -1.42768100
  38. C                 -7.48024100   -0.53011200   -0.98603100
  39. C                 -8.73219500    0.35264300   -1.04220500
  40. C                 -9.99700300   -0.32668300   -0.49689300
  41. C                -10.31850900   -0.12167400    0.98801800
  42. C                 -9.32628200   -0.68300400    2.01734000
  43. C                 -8.22753900    0.28459700    2.45860200
  44. C                 -7.28247500   -0.32495900    3.48828400
  45. C                  5.46963900   -1.94487500   -1.84429500
  46. C                  5.84874200   -3.28624800   -1.21155400
  47. C                  6.72196000   -3.21208200    0.04308700
  48. C                  8.12963800   -2.63974600   -0.15187700
  49. C                  9.04824000   -3.51593500   -1.00130900
  50. C                 10.44097400   -2.93855800   -1.26862700
  51. C                 11.30254400   -2.67434200   -0.02283500
  52. C                 11.21948400   -1.24908200    0.51870300
  53. C                  5.32889500   -0.03998200   -0.20273300
  54. C                  5.44701200    1.27449400   -0.97000100
  55. C                  6.12517100    2.37937700   -0.15973000
  56. C                  7.59386400    2.10940200    0.17004500
  57. C                  8.32025800    3.30206500    0.79385700
  58. C                  7.79329000    3.71905700    2.16491200
  59. C                  8.59310300    4.85789600    2.79175700
  60. C                  8.06415400    5.27623300    4.15852600
  61. H                  0.61922200   -2.56140300   -3.47802400
  62. H                  0.43797100   -3.76266700   -2.17540800
  63. H                  2.79207400   -1.85181100   -2.92962700
  64. H                  2.90588000   -3.12392800   -1.68979800
  65. H                 -3.92634200   -1.79057800   -2.79751200
  66. H                 -5.14969200   -0.77420100   -3.63116900
  67. H                 -5.67568200   -2.10910600   -2.58769700
  68. H                 -3.45343000   -1.25668200   -0.36271500
  69. H                 -4.74651600   -0.43839900    0.51942100
  70. H                 -6.12470100   -2.60739000   -0.61015200
  71. H                 -5.50995100   -2.53607600    1.00913500
  72. H                 -3.89763300   -3.64955000   -1.34768600
  73. H                 -4.96225600   -4.61225700   -0.32809400
  74. H                 -2.44079500   -3.13933100    0.59991600
  75. H                 -2.61620400   -4.84463100    0.28538300
  76. H                 -3.02077600   -5.02150600    2.53293100
  77. H                 -4.64573500   -4.49993900    2.14519800
  78. H                 -3.90301500   -2.09264700    2.59827800
  79. H                 -2.31520600   -2.71224200    3.01030900
  80. H                 -4.91966800   -3.43851000    4.45245600
  81. H                 -3.33266000   -4.07568800    4.85655800
  82. H                 -3.80004300    1.20245000   -0.81361900
  83. H                 -2.94195700    0.10621400   -1.88574200
  84. H                 -3.01080400    1.93038400   -3.18837100
  85. H                 -4.40904100    1.11890600   -3.82428700
  86. H                 -3.91640300   -2.17769000    6.37535500
  87. H                 -4.12170700   -1.10377600    4.98164800
  88. H                 -2.52067500   -1.74559700    5.37160500
  89. H                 -5.88925500    2.70799800   -2.46351000
  90. H                 -4.89499100    3.43901200   -3.69976300
  91. H                 -4.16097400    3.42640700   -0.72823200
  92. H                 -5.02471300    4.71056700   -1.54702300
  93. H                 -3.05794500    4.83830400   -3.19958800
  94. H                 -2.14507300    3.82651700   -2.09562400
  95. H                 -3.37191800    6.54136800   -1.37217000
  96. H                 -1.73387200    6.29338300   -1.95027500
  97. H                 -1.48004500    6.54806800    0.33890600
  98. H                 -1.29530400    4.82356100    0.09099500
  99. H                 -2.65732500    5.42465100    2.12761200
  100. H                 -3.88454200    6.16376500    1.07969500
  101. H                 -3.57922700    4.41756800    1.01889300
  102. H                 -6.11123300    1.10262300   -0.72085500
  103. H                 -6.47423100    0.72376600   -2.40753100
  104. H                 -7.33716100   -0.90984000    0.02940500
  105. H                 -7.64026700   -1.39687700   -1.63943000
  106. H                 -8.55585300    1.29245200   -0.49955300
  107. H                 -8.89849400    0.63728600   -2.09025800
  108. H                 -9.94773900   -1.40468000   -0.71091500
  109. H                -10.85781300    0.04740400   -1.06555400
  110. H                -11.29903800   -0.58355100    1.16205800
  111. H                -10.46078600    0.95272900    1.17752200
  112. H                 -8.87801900   -1.61284700    1.63543100
  113. H                 -9.88224400   -0.97828300    2.91776600
  114. H                 -8.69963400    1.18294100    2.88053500
  115. H                 -7.64935000    0.63489400    1.59314900
  116. H                 -6.51539000    0.39027700    3.80544600
  117. H                 -7.82938900   -0.64539100    4.38310500
  118. H                 -6.76926300   -1.20859700    3.08665800
  119. H                  6.36627100   -1.37516600   -2.10883700
  120. H                  4.94125800   -2.12026000   -2.78723400
  121. H                  4.93015200   -3.83434500   -0.96501100
  122. H                  6.35162800   -3.87680400   -1.98836700
  123. H                  6.19326900   -2.63458100    0.81295300
  124. H                  6.81314300   -4.22813800    0.45016000
  125. H                  8.57670700   -2.50607000    0.84130100
  126. H                  8.08146300   -1.63333200   -0.59427300
  127. H                  8.57276300   -3.71794600   -1.97137200
  128. H                  9.15572500   -4.49376600   -0.50832000
  129. H                 10.34537400   -2.00737000   -1.84690900
  130. H                 10.96680100   -3.64392900   -1.92520700
  131. H                 12.35220900   -2.88494900   -0.26455200
  132. H                 11.03116700   -3.38985300    0.76677100
  133. H                 11.55235900   -0.52795200   -0.23781800
  134. H                 11.85887100   -1.12336600    1.40000100
  135. H                 10.19984200   -0.97267200    0.80719000
  136. H                  6.31180700   -0.42789500    0.07259600
  137. H                  4.75493500    0.10779600    0.71551800
  138. H                  6.00629300    1.10516900   -1.90188000
  139. H                  4.43855700    1.59632300   -1.25583800
  140. H                  5.55317200    2.54720800    0.76246900
  141. H                  6.05848700    3.31449100   -0.73237400
  142. H                  8.11572300    1.82089600   -0.75473000
  143. H                  7.67883800    1.24982100    0.85091300
  144. H                  8.26160700    4.15936500    0.10668200
  145. H                  9.38829800    3.05761400    0.88498700
  146. H                  7.80951500    2.84990800    2.84006300
  147. H                  6.74050500    4.02406900    2.08732700
  148. H                  8.58244800    5.72229900    2.11274400
  149. H                  9.64619700    4.55452300    2.88045100
  150. H                  8.09296800    4.43898800    4.86640900
  151. H                  7.02331800    5.61533700    4.09199200
  152. H                  8.65496700    6.09441700    4.58552300
  153. N                  0.03968300    1.79181500   -2.40599400
  154. O                  1.11027200    1.15545600   -2.15526300
  155. O                 -0.06977900    2.49774000   -3.38722800
  156. O                 -0.90426400    1.64009500   -1.57704300
  157. N                  0.32869500   -1.58185000    2.28781000
  158. O                 -0.35517400   -1.90566100    1.26076900
  159. O                  0.29749000   -2.24208300    3.30391700
  160. O                  1.03901500   -0.54495200    2.17129100
  161. N                 -0.68835600    2.06950900    1.54544400
  162. O                  0.45775700    2.13401400    1.01068200
  163. O                 -1.10283200    2.93354100    2.28550700
  164. O                 -1.38460700    1.04070700    1.26470300
  165. Eu                 0.27911000   -0.03000200   -0.14958300

  167. C H O N 0
  168. 6-31g(d)
  169. ****
  170. Eu 0
  171. S 5 1.00
  172. 70059.4200000 0.000097
  173. 10776.2350000 0.000707
  174. 2482.4901000 0.003362
  175. 702.1526000 0.010206
  176. 216.7926000 0.016597
  177. S 1 1.00
  178. 41.4950000 1.0
  179. S 1 1.00
  180. 29.4969000 1.0
  181. S 1 1.00
  182. 15.0309000 1.0
  183. S 1 1.00
  184. 3.7772000 1.0
  185. S 1 1.00
  186. 1.9219000 1.0
  187. S 1 1.00
  188. 0.7131000 1.0
  189. S 1 1.00
  190. 0.3097000 1.0
  191. S 1 1.00
  192. 0.0549000 1.0
  193. S 1 1.00
  194. 0.0227000 1.0
  195. P 6 1.00
  196. 4058.0207000 0.000048
  197. 1023.3259000 0.000357
  198. 339.3852000 0.001687
  199. 122.3485000 0.005693
  200. 29.0558000 0.072398
  201. 20.6269000 -0.046807
  202. P 1 1.00
  203. 14.6179000 1.0
  204. P 1 1.00
  205. 5.4444000 1.0
  206. P 1 1.00
  207. 2.7414000 1.0
  208. P 1 1.00
  209. 1.2283000 1.0
  210. P 1 1.00
  211. 0.5769000 1.0
  212. P 1 1.00
  213. 0.2489000 1.0
  214. P 1 1.00
  215. 0.0800000 1.0
  216. D 6 1.00
  217. 389.0841000 0.000534
  218. 117.7602000 0.004258
  219. 44.9869000 0.017160
  220. 22.8599000 0.012427
  221. 7.8331000 0.270416
  222. 4.1134000 0.445272
  223. D 1 1.00
  224. 2.0999000 1.0
  225. D 1 1.00
  226. 1.0043000 1.0
  227. D 1 1.00
  228. 0.3422000 1.0
  229. D 1 1.00
  230. 0.1051000 1.0
  231. F 5 1.00
  232. 128.2341000 0.003424
  233. 46.2685000 0.034313
  234. 20.7366000 0.123194
  235. 9.4479000 0.251824
  236. 4.3149000 0.351292
  237. F 1 1.00
  238. 1.9121000 1.0
  239. F 1 1.00
  240. 0.7897000 1.0
  241. F 1 1.00
  242. 0.2859000 1.0
  243. G 4 1.00
  244. 20.7366000 0.026516
  245. 9.4479000 0.099066
  246. 4.3149000 0.213378
  247. 1.9121000 0.453713
  248. G 1 1.00
  249. 0.7897000 1.0
  250. G 1 1.00
  251. 0.2859000 1.0
  252. ****

  254. Eu 0
  255. ECP28MWB 5 28
  256. H-Komponente
  257. 1
  258. 2 1.000000 0.000000
  259. S-H
  260. 1
  261. 2 23.471384 607.659331
  262. P-H
  263. 1
  264. 2 16.772479 264.385476
  265. D-H
  266. 1
  267. 2 13.981343 115.381375
  268. F-H
  269. 1
  270. 2 23.962888 -49.400794
  271. G-H
  272. 1
  273. 2 21.232458 -26.748273

  275. surface=ses addsph
复制代码




作者
Author: wzkchem5    时间: 2022-9-27 15:44
本帖最后由 wzkchem5 于 2022-9-27 08:47 编辑

从能量最低的那个结构开始继续优化，此外考虑不要用opt=cart。
另外，这个计算必须加色散校正，因为你的体系明显是色散作用非常重要的体系，如果不加色散的话，配体的烷基“尾巴”摆的角度肯定不对，不加色散的结果是不能发表的。这个可能也和不收敛有一定关系，参见https://mattermodeling.stackexch ... structure/9668#9668
作者
Author: soda123    时间: 2022-9-27 22:12
wzkchem5 发表于 2022-9-27 15:44
从能量最低的那个结构开始继续优化，此外考虑不要用opt=cart。
另外，这个计算必须加色散校正，因为你的体 ...

谢谢您的回答！我再计算试一下
作者
Author: soda123    时间: 2022-9-29 14:59
wzkchem5 发表于 2022-9-27 15:44
从能量最低的那个结构开始继续优化，此外考虑不要用opt=cart。
另外，这个计算必须加色散校正，因为你的体 ...

老师您好，我按照您说的方法计算后在后期又出现了波动的问题，但波动的幅度明显减小，请问我现在是等待任务再算下去看下一步趋势还是再做一些修改再计算
作者
Author: soda123    时间: 2022-9-29 14:59
这是现在计算的图片
作者
Author: wzkchem5    时间: 2022-9-29 15:45
soda123 发表于 2022-9-29 07:59
老师您好，我按照您说的方法计算后在后期又出现了波动的问题，但波动的幅度明显减小，请问我现在是等待任 ...

检查一下能量波动是因为结构突变还是波函数突变
作者
Author: soda123    时间: 2022-9-29 21:54
结构突变老师 在两个结构优化点之间我观察到高斯优化的结构出现了大的波动
作者
Author: soda123    时间: 2022-9-29 21:55
wzkchem5 发表于 2022-9-29 15:45
检查一下能量波动是因为结构突变还是波函数突变

结构突变老师，我观察高斯优化的两个突变点的优化结构出现了大的变化
作者
Author: wzkchem5    时间: 2022-9-29 22:34
soda123 发表于 2022-9-29 14:55
结构突变老师，我观察高斯优化的两个突变点的优化结构出现了大的变化

那可能就是高斯优化这种体系本来就比较低效。如果仍然没有收敛趋势的话，可以把任务kill掉，然后从能量最低的结构继续优化
作者
Author: soda123    时间: 2022-9-29 23:44
wzkchem5 发表于 2022-9-29 22:34
那可能就是高斯优化这种体系本来就比较低效。如果仍然没有收敛趋势的话，可以把任务kill掉，然后从能量最 ...

好的老师，谢谢您的回答！
作者
Author: soda123    时间: 2022-10-14 08:46
wzkchem5 发表于 2022-9-29 22:34
那可能就是高斯优化这种体系本来就比较低效。如果仍然没有收敛趋势的话，可以把任务kill掉，然后从能量最 ...

老师您好，我从能量最低开始算后出现了几何优化不收敛的问题，我也尝试了sob老师文章中的方法，但仍然出现后期震荡的情况，麻烦您解答一下 下面是我的优化步数图和输出文件
作者
Author: soda123    时间: 2022-10-14 08:47
图式如下
作者
Author: soda123    时间: 2022-10-14 08:48
这个是log文件
作者
Author: wzkchem5    时间: 2022-10-14 14:23
soda123 发表于 2022-10-14 01:46
老师您好，我从能量最低开始算后出现了几何优化不收敛的问题，我也尝试了sob老师文章中的方法，但仍然出 ...

这种体系本来就难收敛，因为太柔性了。几个可以尝试的点
（1）输出文件明确说了不建议用GDIIS，因为从这个版本的高斯开始，GDIIS基本不会是最好的选择。试试关掉GDIIS
（2）可能初始Hessian质量不好，试试计算一下初始Hessian
（3）可能格点不够大（考虑到Eu比较重，需要比较大的径向格点，而把f电子积准可能对角向格点也有一定要求），试试增加格点
作者
Author: soda123    时间: 2022-10-18 16:02
wzkchem5 发表于 2022-10-14 14:23
这种体系本来就难收敛，因为太柔性了。几个可以尝试的点
（1）输出文件明确说了不建议用GDIIS，因为从这 ...

好的，谢谢老师！



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