应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » 计算含la系的Dy对应富勒烯体系，总是不收敛，大家帮我分析分析
51/1返回列表
查看: 1644  |  回复: 12
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看
当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

zc311213

木虫 (小有名气)

[交流] 计算含la系的Dy对应富勒烯体系，总是不收敛，大家帮我分析分析

情况：计算的是DySc2@C80,已经尝试将近一月，但总不收敛，加赝势用ms和gaussian都算过，下面贴的是gaussian计算的输入和输出相关

问题：交流下解决此体系不收敛问题方案，注：即使是解决思路也行，最好具体点，谢谢了


输入文件如下：

%nprocshared=12
%mem=2000MW
%rwf=1,2000mb,2,2000mb,3,2000mb,4,2000mb,5,2000mb,6,2000mb,7,2000mb,8,2000mb,9,-1
# opt=maxcycle=200 b3lyp/genecp int(NoXCTest) guess=indo geom=connectivity

Title Card Required

0 2
C                  0.19731700   -0.84518300   -4.05654300
C                 -1.13574200   -0.24245700   -4.02219600
C                  1.38615800   -0.04503100   -3.84624000
C                  1.21754800    1.38408200   -3.69612200
C                 -0.11731400    1.98955000   -3.66346100
C                 -1.31053000    1.17577800   -3.78388700
C                  0.09254300   -2.20664900   -3.53160900
C                 -2.06913300   -1.23052900   -3.47800400
C                  2.50645000   -0.66846800   -3.16944900
C                  2.09663300    2.15743900   -2.81716500
C                 -0.06781800    3.13869300   -2.75545400
C                 -1.31106100   -2.44658200   -3.17097200
C                 -2.48875600    1.59167900   -3.04558500
C                  1.29946100    3.24205300   -2.23299200
C                  1.17869700   -2.79642600   -2.77259000
C                  2.40185300   -2.02596500   -2.63552900
C                  3.41116400    0.09945400   -2.30230900
C                 -3.19780500   -0.82237500   -2.66235400
C                  3.17838800    1.52211800   -2.07778700
C                 -3.41980200    0.60454200   -2.48841100
C                 -2.45425800    2.75399900   -2.14469500
C                 -1.22574600    3.52667900   -1.96049300
C                 -1.67417000   -3.30767700   -2.05380000
C                  3.23328500   -2.10416400   -1.42958100
C                  3.89081600   -0.79160200   -1.22354900
C                  0.82849000   -3.70264800   -1.69853600
C                  1.54645400    3.70660700   -0.88483900
C                 -3.94899400    1.14394900   -1.23049900
C                 -3.60109300   -1.71008100   -1.58303200
C                 -3.37380400    2.48912800   -1.01883300
C                  3.48928300    2.05092600   -0.75590600
C                 -0.57582800   -3.95467300   -1.33613600
C                 -0.97504500    4.09600900   -0.63644800
C                 -2.88216200   -2.96816900   -1.29246000
C                  2.66608300    3.11849600   -0.16714200
C                  0.40236300    4.15666600   -0.10801000
C                  1.63898700   -3.77029800   -0.48047200
C                  2.82458900   -2.94035300   -0.31360700
C                  4.15893100   -0.27680000    0.13074100
C                  3.96946700    1.17255000    0.33449900
C                 -4.11691800   -1.16557000   -0.32069300
C                 -4.26935100    0.26879100   -0.11832600
C                 -3.09496700    2.99504200    0.33182600
C                 -1.89649800    3.84546300    0.50165600
C                 -0.63403300   -4.17199300    0.11933700
C                 -2.95224900   -3.19923200    0.17926200
C                  2.60415400    2.89972700    1.27663200
C                  0.73904000   -4.04628200    0.64239600
C                  0.34379800    3.93229600    1.33810800
C                 -3.69743000   -2.05722100    0.76594800
C                  3.12738600   -2.46268600    1.02433800
C                  3.38995800    1.70069000    1.59274200
C                  3.77078900   -1.14984800    1.24926600
C                 -1.79614200   -3.75956300    0.90734700
C                 -1.05924400    3.71382000    1.71343500
C                 -4.03844100    0.78003900    1.21780400
C                 -3.44871100    2.11063000    1.44267500
C                  1.42566000    3.26423000    2.03927500
C                  0.99344500   -3.49282500    1.95598600
C                  2.21152200   -2.72623900    2.13695600
C                  3.20963000   -0.60518800    2.50403400
C                 -3.38543600   -1.52101400    2.08459200
C                  2.98002600    0.82108700    2.67446100
C                 -3.59680600   -0.10056200    2.29932900
C                 -2.62819100    2.03284700    2.65877300
C                 -1.40617600    2.80373300    2.79308900
C                 -1.52685500   -3.24702800    2.24674000
C                  2.26072700   -1.58508400    3.04809900
C                  1.07832700    2.43420500    3.16971700
C                 -0.15690800   -3.12986800    2.76146400
C                 -2.32380000   -2.16348600    2.83898000
C                 -2.72576900    0.66924000    3.18565900
C                  1.83998800    1.22565700    3.47718300
C                 -0.31889100    2.20272800    3.54118400
C                  1.08489500   -1.17308100    3.78993500
C                 -0.10809500   -1.98451100    3.66911500
C                 -1.44294600   -1.38473200    3.71346500
C                 -1.61172300    0.04494600    3.86643400
C                  0.91012200    0.24168300    4.03137300
C                 -0.42316600    0.84430000    4.07303600
N                 -0.51982447    0.03293313    1.37632584
Sc                -1.18406411    1.71528091   -0.24896163
Sc                -1.24751649   -1.68493495   -0.26030860
Dy                 1.82706606   -0.01930809   -0.20426762

1 2 1.0 3 1.0 7 1.0
2 6 1.0 8 1.0
3 4 1.5 9 1.0
4 5 1.0 10 1.0
5 6 1.0 11 1.0
6 13 1.0
7 12 1.0 15 1.0
8 12 1.0 18 1.0
9 16 1.0 17 1.0
10 14 1.0 19 1.0
11 14 1.0 22 1.0
12 23 1.0
13 20 1.0 21 1.0
14 27 1.5
15 16 1.0 26 1.0
16 24 1.0
17 19 1.0 25 1.0
18 20 1.0 29 1.0
19 31 1.0
20 28 1.0
21 22 1.0 30 1.0
22 33 1.0
23 32 1.0 34 1.0
24 25 1.0 38 1.0
25 39 1.0
26 32 1.0 37 1.0
27 35 1.0 36 1.0
28 30 1.0 42 1.0
29 34 1.0 41 1.0
30 43 1.0
31 35 1.0 40 1.0
32 45 1.0
33 36 1.0 44 1.0
34 46 1.0
35 47 1.0
36 49 1.0
37 38 1.0 48 1.0
38 51 1.0
39 40 1.0 53 1.0
40 52 1.0
41 42 1.0 50 1.0
42 56 1.0
43 44 1.0 57 1.0
44 55 1.0
45 48 1.0 54 1.0
46 50 1.0 54 1.0
47 52 1.0 58 1.0
48 59 1.0
49 55 1.0 58 1.0
50 62 1.0
51 53 1.0 60 1.0
52 63 1.0
53 61 1.0
54 67 1.0
55 66 1.0
56 57 1.0 64 1.0
57 65 1.0
58 69 1.5
59 60 1.0 70 1.0
60 68 1.0
61 63 1.0 68 1.0
62 64 1.0 71 1.0
63 73 1.0
64 72 1.0
65 66 1.0 72 1.0
66 74 1.0
67 70 1.0 71 1.0
68 75 1.0
69 73 1.0 74 1.0
70 76 1.0
71 77 1.0
72 78 1.5
73 79 1.0
74 80 1.0
75 76 1.0 79 1.5
76 77 1.0
77 78 1.0
78 80 1.5
79 80 1.0
80
81 82 1.0 83 1.0 84 1.0
82 84 1.0
83 84 1.0
84

Dy     0
S   6   1.00
21841438.1922420            0.0009760        
8610654.2403350            -0.0000320        
3394619.2459490             0.0019940        
1338276.9187260             0.0020150        
527595.2857840              0.0053940        
207996.4032010              0.0093770        
S   1   1.00
  81999.4130160              1.0000000        
S   1   1.00
  32327.0192730              1.0000000        
S   1   1.00
  12744.4348280              1.0000000        
S   1   1.00
   5024.2992620              1.0000000        
S   1   1.00
   1980.7534360              1.0000000        
S   1   1.00
    780.8818640              1.0000000        
S   1   1.00
    307.8507780              1.0000000        
S   1   1.00
    121.3654790              1.0000000        
S   1   1.00
     47.8464910              1.0000000        
S   1   1.00
     18.8627500              1.0000000        
S   1   1.00
      7.4363520              1.0000000        
S   1   1.00
      2.9316680              1.0000000        
S   1   1.00
      1.1557660              1.0000000        
S   1   1.00
      0.4556430              1.0000000        
S   1   1.00
      0.1796300              1.0000000        
S   1   1.00
      0.0708170              1.0000000        
S   1   1.00
      0.0278140              1.0000000        
P   5   1.00
  35772.9954310              0.0021270        
  14694.4030130              0.0018130        
   6035.9910400              0.0100910        
   2479.3921750              0.0277090        
   1018.4550500              0.0935520        
P   1   1.00
    418.3487790              1.0000000        
P   1   1.00
    171.8443050              1.0000000        
P   1   1.00
     70.5881480              1.0000000        
P   1   1.00
     28.9953550              1.0000000        
P   1   1.00
     11.9103650              1.0000000        
P   1   1.00
      4.8923970              1.0000000        
P   1   1.00
      2.0096400              1.0000000        
P   1   1.00
      0.8254960              1.0000000        
P   1   1.00
      0.3390870              1.0000000        
P   1   1.00
      0.1392860              1.0000000        
P   1   1.00
      0.0395080              1.0000000        
D   4   1.00
   1487.0265820              0.0029570        
    643.7378760              0.0097430        
    278.6758880              0.0490390        
    120.6395540              0.1724900        
D   1   1.00
     52.2251930              1.0000000        
D   1   1.00
     22.6084290              1.0000000        
D   1   1.00
      9.7872510              1.0000000        
D   1   1.00
      4.2369280              1.0000000        
D   1   1.00
      1.8341780              1.0000000        
D   1   1.00
      0.7940200              1.0000000        
D   1   1.00
      0.3437340              1.0000000        
D   1   1.00
      0.1124520              1.0000000        
F   4   1.00
     59.9032870              0.0377380        
     20.7195140              0.1791130        
      7.1665230              0.4235130        
      2.4787770              0.4462420        
F   1   1.00
      0.8573660              1.0000000        
F   1   1.00
      0.2866140              1.0000000        
****
C N Sc  0
3-21g
****
Dy     0
S   6   1.00
21841438.1922420            0.0009760        
8610654.2403350            -0.0000320        
3394619.2459490             0.0019940        
1338276.9187260             0.0020150        
527595.2857840              0.0053940        
207996.4032010              0.0093770        
S   1   1.00
  81999.4130160              1.0000000        
S   1   1.00
  32327.0192730              1.0000000        
S   1   1.00
  12744.4348280              1.0000000        
S   1   1.00
   5024.2992620              1.0000000        
S   1   1.00
   1980.7534360              1.0000000        
S   1   1.00
    780.8818640              1.0000000        
S   1   1.00
    307.8507780              1.0000000        
S   1   1.00
    121.3654790              1.0000000        
S   1   1.00
     47.8464910              1.0000000        
S   1   1.00
     18.8627500              1.0000000        
S   1   1.00
      7.4363520              1.0000000        
S   1   1.00
      2.9316680              1.0000000        
S   1   1.00
      1.1557660              1.0000000        
S   1   1.00
      0.4556430              1.0000000        
S   1   1.00
      0.1796300              1.0000000        
S   1   1.00
      0.0708170              1.0000000        
S   1   1.00
      0.0278140              1.0000000        
P   5   1.00
  35772.9954310              0.0021270        
  14694.4030130              0.0018130        
   6035.9910400              0.0100910        
   2479.3921750              0.0277090        
   1018.4550500              0.0935520        
P   1   1.00
    418.3487790              1.0000000        
P   1   1.00
    171.8443050              1.0000000        
P   1   1.00
     70.5881480              1.0000000        
P   1   1.00
     28.9953550              1.0000000        
P   1   1.00
     11.9103650              1.0000000        
P   1   1.00
      4.8923970              1.0000000        
P   1   1.00
      2.0096400              1.0000000        
P   1   1.00
      0.8254960              1.0000000        
P   1   1.00
      0.3390870              1.0000000        
P   1   1.00
      0.1392860              1.0000000        
P   1   1.00
      0.0395080              1.0000000        
D   4   1.00
   1487.0265820              0.0029570        
    643.7378760              0.0097430        
    278.6758880              0.0490390        
    120.6395540              0.1724900        
D   1   1.00
     52.2251930              1.0000000        
D   1   1.00
     22.6084290              1.0000000        
D   1   1.00
      9.7872510              1.0000000        
D   1   1.00
      4.2369280              1.0000000        
D   1   1.00
      1.8341780              1.0000000        
D   1   1.00
      0.7940200              1.0000000        
D   1   1.00
      0.3437340              1.0000000        
D   1   1.00
      0.1124520              1.0000000        
F   4   1.00
     59.9032870              0.0377380        
     20.7195140              0.1791130        
      7.1665230              0.4235130        
      2.4787770              0.4462420        
F   1   1.00
      0.8573660              1.0000000        
F   1   1.00
      0.2866140              1.0000000        
****



输出的不收敛如下：

。。。。。。
Spurious integrated density or basis function:
NE=  595 NElCor=    0 El error=1.43D-03 rel=2.40D-06 Tolerance=1.00D-03
Shell   259     absolute error=7.62D-02              Tolerance=1.20D-02
Shell   259       signed error=7.62D-02              Tolerance=1.00D-01
Spurious integrated density or basis function:
NE=  595 NElCor=    0 El error=1.43D-03 rel=2.40D-06 Tolerance=1.00D-03
Shell   259     absolute error=7.62D-02              Tolerance=1.20D-02
Shell   259       signed error=7.62D-02              Tolerance=1.00D-01
Spurious integrated density or basis function:
NE=  595 NElCor=    0 El error=1.43D-03 rel=2.40D-06 Tolerance=1.00D-03
Shell   259     absolute error=7.62D-02              Tolerance=1.20D-02
Shell   259       signed error=7.62D-02              Tolerance=1.00D-01
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB3LYP) =  -16242.0568259     A.U. after  129 cycles
            NFock=128  Conv=0.74D-05     -V/T= 2.0029
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 0.5000 <S**2>= 2.8652 S= 1.2650
<L.S>= 0.000000000000E+00
Annihilation of the first spin contaminant:
S**2 before annihilation     2.8652,   after     4.6489
Convergence failure -- run terminated.
Error termination via Lnk1e in /home-gg/Soft/Gaussian.D01.Linda/g09/l502.exe at Mon Mar 21 00:28:41 2016.
Job cpu time:       2 days  2 hours 45 minutes  2.3 seconds.
File lengths (MBytes):  RWF=    942 Int=      0 D2E=      0 Chk=     23 Scr=      1
回复此楼

» 猜你喜欢
yes
1楼 2016-04-07 21:16:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

deitydeng

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
引用回帖:
5楼: Originally posted by zc311213 at 2016-04-08 22:51:09
这个现在是已经分出来了，产率还行，我实际上也按照你的想法在gaussian和MS上算过，用La原子换Dy，可以成功收敛，再换回Dy不行！不过还是感谢分享你的想法...

标题就错了，DySc2@C80 分出来最次是jacs或NC级别的，您这是DySc2N@C80 吓死宝宝了，还以为自己落伍了那

发自小木虫IOS客户端
赞一下 回复此楼
7楼2016-04-09 04:17:50
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答

Woodeninsect

银虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
能不能换一个Dy基组和赝势，试试怎么样呢？我用大核算Dy的配合物，能收敛
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

zc311213
2楼2016-04-08 19:10:02
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

deitydeng

木虫 (正式写手)
★
小木虫: 金币+0.5, 给个红包，谢谢回帖
你这个结构不收敛很正常 三金属肯定不稳定 只在质朴上看到过 还没背分出过 技术上建议你用orca 先算三个sc 的 在替换成dy

发自小木虫IOS客户端
赞一下 回复此楼

» 本帖已获得的红花（最新10朵）

zc311213
3楼2016-04-08 20:24:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

zc311213

木虫 (小有名气)
送红花一朵
引用回帖:
2楼: Originally posted by Woodeninsect at 2016-04-08 19:10:02
能不能换一个Dy基组和赝势，试试怎么样呢？我用大核算Dy的配合物，能收敛

也用了sdd赝势，计算不收敛，可否推荐一个你计算成功的基组和赝势
赞一下 回复此楼
yes
4楼2016-04-08 22:46:03
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 13 个回答
普通表情 龙 兔 虎 猫 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 考研调剂 +6 来好运来来来 2026-03-21 7/350 2026-03-25 22:43 by 418490947
[考研] 291 求调剂 +5 化工2026届毕业� 2026-03-21 6/300 2026-03-25 22:28 by 544594351
[考研] 290分调剂求助 +3 吉祥止止陈 2026-03-25 3/150 2026-03-25 19:58 by barlinike
[考研] 求调剂 一志愿 本科 北科大 化学 343 +4 13831862839 2026-03-24 5/250 2026-03-25 09:47 by 无际的草原
[考研] 材料学求调剂 +6 Stella_Yao 2026-03-20 6/300 2026-03-25 00:37 by baoball
[考研] 0854 考研调剂 招生了！AI 方向 +5 pk3725069 2026-03-19 17/850 2026-03-24 17:30 by zhouxuan..
[考研] 292求调剂 +4 鹅鹅鹅额额额额� 2026-03-24 4/200 2026-03-24 16:41 by peike
[考博] 申博26年 +4 八6八68 2026-03-19 4/200 2026-03-24 15:49 by 小Ben呵呵
[考研] 279分求调剂 一志愿211 +18 chaojifeixia 2026-03-19 20/1000 2026-03-24 10:34 by dolphin_ycj
[考研] 344求调剂 +3 desto 2026-03-24 3/150 2026-03-24 10:09 by 搏击518
[考研] 一志愿武理材料工程348求调剂 +6 ￣^￣゜汗 2026-03-19 9/450 2026-03-23 19:53 by pswait
[考研] 276求调剂 +3 YNRYG 2026-03-21 4/200 2026-03-23 08:31 by 醉在风里
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 269专硕求调剂 +6 金恩贝 2026-03-21 6/300 2026-03-22 14:31 by ColorlessPI
[考研] 一志愿深大，0703化学，总分302，求调剂 +4 七月-七七 2026-03-21 4/200 2026-03-21 18:20 by 学员8dgXkO
[考研] 22 350 本科985求调剂，求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 329求调剂 +9 想上学吖吖 2026-03-19 9/450 2026-03-20 22:01 by luoyongfeng
[考研] 求调剂一志愿南京航空航天大学289分 +3 @taotao 2026-03-19 3/150 2026-03-20 21:34 by JourneyLucky
[考研] 招收调剂硕士 +4 lidianxing 2026-03-19 12/600 2026-03-20 12:25 by lidianxing
[考研] 320求调剂0856 +3 不想起名字112 2026-03-19 3/150 2026-03-19 22:53 by 学员8dgXkO
信息提示
关闭
请填处理意见
关闭