计算化学公社

标题: 请教稀土配合物优化失败如何解决 [打印本页]
作者
Author: bomsaude    时间: 2014-12-17 18:33
标题: 请教稀土配合物优化失败如何解决
本帖最后由 bomsaude 于 2014-12-17 18:59 编辑

接触高斯有半年了，原先按书上和网上的一些小分子进行了了一些演练，现在想练习一些大一点的分子，几何优化就不收敛了。稀土配合物是最常见的Eu(TTA)3phen，C H O N F采用6-31g基组，S采用6-31g(d)基组，Eu采用赝势基组，输入文件如下，不知道有什么不对或不合适的地方？看到文献上写All the molecular structures were fully optimized in the gas-phase without symmetry restrictions by the gradient procedure，那么关键词中用不用加nosymm?

  1. <p>%chk=euttaphen.chk
  2. # opt b3lyp/genecp geom=connectivity</p><p> </p><p>Opt of EuTTAphen</p><p> </p><p>0 1
  3. Eu                 1.46400000    9.32300000    4.44890000
  4. H                  6.51600000   12.70290000    9.83990000
  5. H                  6.66190000   15.49600000    3.37000000
  6. H                  5.66000000   10.54290000   -2.32690000
  7. S                  4.34400000   11.48090000    7.97790000
  8. S                  4.48500000   13.32190000    3.64000000
  9. S                  4.01100000    9.72300000    0.15300000
  10. F                 -2.73490000    9.45400000    7.36390000
  11. F                 -1.72290000    9.15790000    9.21190000
  12. F                 -1.76290000    7.65590000    7.77890000
  13. F                 -1.68790000   10.32790000    0.96400000
  14. F                 -2.36300000   11.83600000    2.13100000
  15. F                 -1.23700000   12.29400000    0.49000000
  16. F                 -0.96600000    5.24890000    3.32690000
  17. F                 -0.68400000    5.14900000    1.20990000
  18. F                 -1.99290000    6.57890000    2.03690000
  19. O                  2.22190000   10.31900000    6.44800000
  20. O                 -0.26500000    9.15490000    6.05600000
  21. O                  2.37990000   11.40190000    3.79000000
  22. O                 -0.19500000   10.50090000    3.25400000
  23. O                  0.23600000    7.58990000    3.45100000
  24. O                  2.34690000    8.94290000    2.29990000
  25. N                  3.96890000    8.71100000    4.65500000
  26. N                  2.04700000    7.08900000    5.64100000
  27. C                  1.83000000   10.49390000    7.61700000
  28. C                  4.33600000    7.48400000    5.11200000
  29. C                  0.52400000   10.10990000    8.05590000
  30. H                  0.27200000   10.32600000    8.92390000
  31. C                 -0.37200000    9.44590000    7.27090000
  32. C                  2.77800000   11.08090000    8.55790000
  33. C                  2.12800000   12.21400000    2.85290000
  34. C                  0.96500000   12.18590000    2.07390000
  35. H                  0.90500000   12.76290000    1.34690000
  36. C                  3.33190000    6.63100000    5.63790000
  37. C                 -1.64590000    8.96100000    7.93600000
  38. C                  1.11400000    7.21490000    1.30600000
  39. H                  1.05200000    6.65590000    0.56490000
  40. C                  4.92100000    9.49600000    4.14100000
  41. H                  4.67300000   10.30700000    3.75990000
  42. C                  1.10890000    6.28500000    6.12500000
  43. H                  0.23300000    6.59690000    6.14400000
  44. C                  3.86600000   11.92090000   10.43500000
  45. H                  3.95100000   12.19200000   11.32090000
  46. C                  2.68100000   11.34690000    9.86400000
  47. H                  1.91300000   11.17390000   10.35990000
  48. C                 -0.11300000   11.33090000    2.33800000
  49. C                  6.62990000    7.94090000    4.60700000
  50. H                  7.52890000    7.70200000    4.59790000
  51. C                  5.98890000    5.75390000    5.53890000
  52. H                  6.86200000    5.43790000    5.48090000
  53. C                  3.67390000    5.36590000    6.08190000
  54. C                  5.68990000    7.04390000    5.08300000
  55. C                  2.05690000    8.21000000    1.29090000
  56. C                  4.82000000   12.00790000    9.51990000
  57. C                  0.26300000    7.00390000    2.34690000
  58. C                  6.26400000    9.14790000    4.16090000
  59. H                  6.90800000    9.75090000    3.86300000
  60. C                  3.13700000   14.27290000    1.69700000
  61. H                  2.49690000   14.38100000    1.03090000
  62. C                  2.66500000    4.56790000    6.59100000
  63. H                  2.87190000    3.72500000    6.92800000
  64. C                  5.04800000    4.97600000    6.05490000
  65. H                  5.29000000    4.15090000    6.40800000
  66. C                  1.35700000    5.00300000    6.60500000
  67. H                  0.67000000    4.46190000    6.92190000
  68. C                 -1.37300000   11.42090000    1.47100000
  69. C                 -0.82200000    5.90800000    2.22500000
  70. C                  3.10990000   13.26200000    2.61600000
  71. C                  2.85400000    8.02890000   -1.14790000
  72. H                  2.23090000    7.38900000   -1.40600000
  73. C                  5.03690000   14.71190000    2.83900000
  74. C                  4.25900000   15.14790000    1.86900000
  75. H                  4.41390000   15.91600000    1.36790000
  76. C                  4.42490000    9.46190000   -1.53290000
  77. C                  3.65390000    8.46100000   -1.93200000
  78. H                  3.72790000    8.10690000   -2.79000000
  79. C                  2.88100000    8.48400000    0.08400000</p><p> </p><p>1 17 1.0 18 1.0 19 1.0 20 1.0 21 1.0 22 1.0 23 1.0 24 1.0
  80. 2 54 1.0
  81. 3 71 1.0
  82. 4 74 1.0
  83. 5 30 1.0 54 1.0
  84. 6 68 1.0 71 1.0
  85. 7 74 1.0 77 1.0
  86. 8 35 1.0
  87. 9 35 1.0
  88. 10 35 1.0
  89. 11 66 1.0
  90. 12 66 1.0
  91. 13 66 1.0
  92. 14 67 1.0
  93. 15 67 1.0
  94. 16 67 1.0
  95. 17 25 1.0
  96. 18 29 1.0
  97. 19 31 1.0
  98. 20 46 1.0
  99. 21 55 1.0
  100. 22 53 1.0
  101. 23 26 1.0 38 1.0
  102. 24 34 1.0 40 1.0
  103. 25 27 1.0 30 1.0
  104. 26 34 1.0 52 1.0
  105. 27 28 1.0 29 1.0
  106. 28
  107. 29 35 1.0
  108. 30 44 1.0
  109. 31 32 1.0 68 1.0
  110. 32 33 1.0 46 1.0
  111. 33
  112. 34 51 1.0
  113. 35
  114. 36 37 1.0 53 1.0 55 1.0
  115. 37
  116. 38 39 1.0 56 1.0
  117. 39
  118. 40 41 1.0 64 1.0
  119. 41
  120. 42 43 1.0 44 1.0 54 1.0
  121. 43
  122. 44 45 1.0
  123. 45
  124. 46 66 1.0
  125. 47 48 1.0 52 1.0 56 1.0
  126. 48
  127. 49 50 1.0 52 1.0 62 1.0
  128. 50
  129. 51 60 1.0 62 1.0
  130. 52
  131. 53 77 1.0
  132. 54
  133. 55 67 1.0
  134. 56 57 1.0
  135. 57
  136. 58 59 1.0 68 1.0 72 1.0
  137. 59
  138. 60 61 1.0 64 1.0
  139. 61
  140. 62 63 1.0
  141. 63
  142. 64 65 1.0
  143. 65
  144. 66
  145. 67
  146. 68
  147. 69 70 1.0 75 1.0 77 1.0
  148. 70
  149. 71 72 1.0
  150. 72 73 1.0
  151. 73
  152. 74 75 1.0
  153. 75 76 1.0
  154. 76
  155. 77</p><p>C H O F N 0
  156. 6-31g
  157. ****
  158. S 0
  159. 6-31g(d)
  160. ****
  161. Eu     0
  162. S   5   1.00
  163.   70059.4200000              0.0000970        
  164.   10776.2350000              0.0007070        
  165.    2482.4901000              0.0033620        
  166.     702.1526000              0.0102060        
  167.     216.7926000              0.0165970        
  168. S   1   1.00
  169.      41.4950000              1.0000000        
  170. S   1   1.00
  171.      29.4969000              1.0000000        
  172. S   1   1.00
  173.      15.0309000              1.0000000        
  174. S   1   1.00
  175.       3.7772000              1.0000000        
  176. S   1   1.00
  177.       1.9219000              1.0000000        
  178. S   1   1.00
  179.       0.7131000              1.0000000        
  180. S   1   1.00
  181.       0.3097000              1.0000000        
  182. S   1   1.00
  183.       0.0549000              1.0000000        
  184. S   1   1.00
  185.       0.0227000              1.0000000        
  186. P   6   1.00
  187.    4058.0207000              0.0000480        
  188.    1023.3259000              0.0003570        
  189.     339.3852000              0.0016870        
  190.     122.3485000              0.0056930        
  191.      29.0558000              0.0723980        
  192.      20.6269000             -0.0468070        
  193. P   1   1.00
  194.      14.6179000              1.0000000        
  195. P   1   1.00
  196.       5.4444000              1.0000000        
  197. P   1   1.00
  198.       2.7414000              1.0000000        
  199. P   1   1.00
  200.       1.2283000              1.0000000        
  201. P   1   1.00
  202.       0.5769000              1.0000000        
  203. P   1   1.00
  204.       0.2489000              1.0000000        
  205. P   1   1.00
  206.       0.0800000              1.0000000        
  207. D   6   1.00
  208.     389.0841000              0.0005340        
  209.     117.7602000              0.0042580        
  210.      44.9869000              0.0171600        
  211.      22.8599000              0.0124270        
  212.       7.8331000              0.2704160        
  213.       4.1134000              0.4452720        
  214. D   1   1.00
  215.       2.0999000              1.0000000        
  216. D   1   1.00
  217.       1.0043000              1.0000000        
  218. D   1   1.00
  219.       0.3422000              1.0000000        
  220. D   1   1.00
  221.       0.1051000              1.0000000        
  222. F   5   1.00
  223.     128.2341000              0.0034240        
  224.      46.2685000              0.0343130        
  225.      20.7366000              0.1231940        
  226.       9.4479000              0.2518240        
  227.       4.3149000              0.3512920        
  228. F   1   1.00
  229.       1.9121000              1.0000000        
  230. F   1   1.00
  231.       0.7897000              1.0000000        
  232. F   1   1.00
  233.       0.2859000              1.0000000        
  234. G   4   1.00
  235.      20.7366000              0.0265160        
  236.       9.4479000              0.0990660        
  237.       4.3149000              0.2133780        
  238.       1.9121000              0.4537130        
  239. G   1   1.00
  240.       0.7897000              1.0000000        
  241. G   1   1.00
  242.       0.2859000              1.0000000        
  243. ****</p><p> </p><p>Eu 0
  244. ECP28MWB 5 28
  245. H-Komponente
  246. 1
  247. 2 1.000000 0.000000
  248. S-H
  249. 1
  250. 2 23.471384 607.659331
  251. P-H
  252. 1
  253. 2 16.772479 264.385476
  254. D-H
  255. 1
  256. 2 13.981343 115.381375
  257. F-H
  258. 1
  259. 2 23.962888 -49.400794
  260. G-H
  261. 1
  262. 2 21.232458 -26.748273</p><p> </p><p> </p>
复制代码


log文件的部分内容如下


  1. ......

  3. There are   685 symmetry adapted basis functions of A   symmetry.
  4. Integral buffers will be    131072 words long.
  5. Raffenetti 2 integral format.
  6. Two-electron integral symmetry is turned on.
  7.    685 basis functions,  1685 primitive gaussians,   723 cartesian basis functions
  8.    231 alpha electrons      231 beta electrons
  9.        nuclear repulsion energy     11089.1497305945 Hartrees.
  10. NAtoms=   77 NActive=   77 NUniq=   77 SFac= 1.00D+00 NAtFMM=   50 NAOKFM=T Big=T
  11. One-electron integrals computed using PRISM.
  12.    1 Symmetry operations used in ECPInt.
  13. ECPInt:  NShTT=   29890 NPrTT=  254997 LenC2=   21948 LenP2D=   80677.
  14. LDataN:  DoStor=T MaxTD1= 9 Len=  602
  15. NBasis=   685 RedAO= T  NBF=   685
  16. NBsUse=   685 1.00D-06 NBFU=   685
  17. Defaulting to unpruned grid for atomic number  63.
  18. Harris functional with IExCor=  402 diagonalized for initial guess.
  19. ExpMin= 2.27D-02 ExpMax= 7.01D+04 ExpMxC= 1.08D+04 IAcc=3 IRadAn=         0 AccDes= 0.00D+00
  20. HarFok:  IExCor=  402 AccDes= 0.00D+00 IRadAn=         0 IDoV= 1
  21. ScaDFX=  1.000000  1.000000  1.000000  1.000000
  22. Defaulting to unpruned grid for atomic number  63.
  23. Defaulting to unpruned grid for atomic number  63.
  24. FoFCou: FMM=F IPFlag=           0 FMFlag=      100000 FMFlg1=        2001
  25.          NFxFlg=           0 DoJE=T BraDBF=F KetDBF=T FulRan=T
  26.          Omega=  0.000000  0.000000  1.000000  0.000000  0.000000 ICntrl=     500 IOpCl=  0
  27.          NMat0=    1 NMatS0=    1 NMatT0=    0 NMatD0=    1 NMtDS0=    0 NMtDT0=    0
  28.          I1Cent=           4 NGrid=           0.

  30. ......

  32. The electronic state of the initial guess is 1-A.
  33. Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
  34. Requested convergence on MAX density matrix=1.00D-06.
  35. Requested convergence on             energy=1.00D-06.
  36. No special actions if energy rises.
  37. Defaulting to unpruned grid for atomic number  63.
  38. Defaulting to unpruned grid for atomic number  63.
  39. EnCoef did     3 forward-backward iterations
  40. EnCoef did   100 forward-backward iterations
  41. EnCoef did   100 forward-backward iterations
  42. EnCoef did   100 forward-backward iterations
  43. EnCoef did     4 forward-backward iterations
  44. Restarting incremental Fock formation.
  45. Matrix for removal  4 Erem= -4537.45091694709     Crem= 0.000D+00
  46. Matrix for removal  3 Erem= -4639.05934792761     Crem= 0.000D+00
  47. Matrix for removal  3 Erem= -4702.55800949193     Crem= 0.000D+00
  48. Matrix for removal  3 Erem= -4723.24313399146     Crem= 0.000D+00
  49. Restarting incremental Fock formation.
  50. Restarting incremental Fock formation.
  51. Restarting incremental Fock formation.
  52. Restarting incremental Fock formation.
  53. Restarting incremental Fock formation.
  54. >>>>>>>>>> Convergence criterion not met.
  55. SCF Done:  E(RB3LYP) =  -4746.74628519     A.U. after  129 cycles
  56.              Convg  =    0.7593D-05             -V/T =  2.0847
  57. Convergence failure -- run terminated.
  58. Error termination via Lnk1e in /home/user/g09/l502.exe at Wed Dec 17 00:33:33 2014.
  59. Job cpu time:  0 days 15 hours 52 minutes 36.1 seconds.
  60. File lengths (MBytes):  RWF=    250 Int=      0 D2E=      0 Chk=      8 Scr=      1</font>
复制代码


出现这种情况，我应该增大循环次数还是修改收敛标准？望有经验的人给予指点，谢谢了！


作者
Author: sobereva    时间: 2014-12-17 20:03
如果初始结构本身就没有对称性就不用加nosymm
解决SCF不收敛没有万能的办法，有可能解决此问题的方法在此贴有汇总：
解决SCF不收敛问题的方法
http://sobereva.com/61
应结合实际情况，将其中的方法挨个试试，或者组合使用。增大循环次数是最没用、还白浪费计算时间的做法。

过渡金属体系不容易收敛是正常的。

另外应当用#P，输出迭代的细节信息，以便于监控收敛到了什么程度，这有助于选择合适的方法帮助收敛。
作者
Author: bomsaude    时间: 2014-12-18 00:02
sobereva 发表于 2014-12-17 20:03
如果初始结构本身就没有对称性就不用加nosymm
解决SCF不收敛没有万能的办法，有可能解决此问题的方法在此 ...

谢谢，我按照您说的方法试试看。



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