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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » gaussian09计算出错，求高手指教
61/1返回列表
查看: 1715  |  回复: 5
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

cuishuai1234

铜虫 (初入文坛)

[求助] gaussian09计算出错，求高手指教

gaissian09优化分子，结果出现如下错误
The electronic state of the initial guess is 7-A'.
Initial guess = 0.0000 = 0.0000 = 3.0000 =12.0000 S= 3.0000
Defaulting to unpruned grid for atomic number  63.
Defaulting to unpruned grid for atomic number  63.
Defaulting to unpruned grid for atomic number  63.
Defaulting to unpruned grid for atomic number  63.
Defaulting to unpruned grid for atomic number  63.
Defaulting to unpruned grid for atomic number  63.
Defaulting to unpruned grid for atomic number  63.
Defaulting to unpruned grid for atomic number  63.
Convergence failure.
Error termination via Lnk1e in /opt/soft/g09/l508.exe at Mon Jul 23 15:14:27 2012.
Job cpu time:  0 days  6 hours 22 minutes 22.6 seconds.
File lengths (MBytes):  RWF=     34 Int=      0 D2E=      0 Chk=      6 Scr=      1


输入文件是：
%ch=6-DMDHOPDA+Eu.chk
%mem=400mw
%nprocshared=8
# opt=maxcycle=1000 bp86/genecp iop(5/13=1) scf=(maxcycle=500,VShift,dsymm,ssd)

Title Card Required

3 7
  C                  0.02824629   -1.50935937   -2.38659573
  C                  0.58832866   -2.25077244   -1.16006815
  H                  1.69362728   -2.24339362   -1.19371637
  H                  0.27797084   -3.31077503   -1.18734229
  C                  0.58832866   -2.25077244    1.16006815
  C                  0.02824629   -1.50935937    2.38659573
  H                  0.27797084   -3.31077503    1.18734229
  H                  1.69362728   -2.24339362    1.19371637
  O                  0.11979949   -1.61844052    0.00000000
  O                 -0.64113220   -0.49537049   -2.26910628
  O                 -0.64113220   -0.49537049    2.26910628
  C                  1.11123016   -3.26015963   -3.82870459
  H                  1.74870742   -3.48026131   -2.97136103
  H                  1.76953707   -3.11913473   -4.69329094
  C                 -0.25220632   -1.42145174   -4.80228175
  C                  0.66874977   -0.37244734   -5.43970730
  H                 -1.19083870   -0.95147709   -4.49967054
  H                 -0.48704693   -2.21409366   -5.52518169
  C                  0.05987591    0.24988516   -6.70206976
  H                  0.86501699    0.40785303   -4.69420754
  H                  1.63869663   -0.82937253   -5.68526895
  C                  0.95827896    1.30982429   -7.35345308
  H                 -0.15537462   -0.54228357   -7.43534819
  H                 -0.91000369    0.70247779   -6.44954140
  H                  1.17485359    2.10011369   -6.61997366
  H                  1.92850575    0.85735503   -7.60795602
  N                  0.32459040   -2.05636186   -3.61372972
  C                  0.34708906    1.93865952   -8.61169706
  C                  1.24443524    2.99699341   -9.26675004
  H                 -0.62156810    2.39349521   -8.35632770
  H                  0.12644389    1.14766882   -9.34431784
  C                  0.62504162    3.62282327  -10.52038291
  H                  1.46579541    3.78635311   -8.53440420
  H                  2.21180932    2.54226568   -9.52423518
  H                  0.42380787    2.86255794  -11.28534712
  H                 -0.32661368    4.11647291  -10.28787963
  C                  1.11123016   -3.26015963    3.82870459
  H                  1.76953707   -3.11913473    4.69329094
  H                  1.74870742   -3.48026131    2.97136103
  C                 -0.25220632   -1.42145174    4.80228175
  C                  0.66874977   -0.37244734    5.43970730
  H                 -0.48704693   -2.21409366    5.52518169
  H                 -1.19083870   -0.95147709    4.49967054
  C                  0.05987591    0.24988516    6.70206976
  H                  1.63869663   -0.82937253    5.68526895
  H                  0.86501699    0.40785303    4.69420754
  C                  0.95827896    1.30982429    7.35345308
  H                 -0.91000369    0.70247779    6.44954140
  H                 -0.15537462   -0.54228357    7.43534819
  H                  1.92850575    0.85735503    7.60795602
  H                  1.17485359    2.10011369    6.61997366
  C                  0.34708906    1.93865952    8.61169706
  C                  1.24443524    2.99699341    9.26675004
  H                  0.12644389    1.14766882    9.34431784
  H                 -0.62156810    2.39349521    8.35632770
  C                  0.62504162    3.62282327   10.52038291
  H                  2.21180932    2.54226568    9.52423518
  H                  1.46579541    3.78635311    8.53440420
  H                 -0.32661368    4.11647291   10.28787963
  H                  0.42380787    2.86255794   11.28534712
  N                  0.32459040   -2.05636186    3.61372972
  H                  0.47707969   -4.13771332   -4.02479650
  H                  0.47707969   -4.13771332    4.02479650
  H                  1.28960930    4.37324418  -10.96369520
  H                  1.28960930    4.37324418   10.96369520
  Eu                -1.37449859    0.39828346    0.00000000

H C N O 0
3-21g
求高手指教，刚刚接触高斯。
回复此楼

» 猜你喜欢

» 本主题相关商家推荐: (我也要在这里推广)

» 本主题相关价值贴推荐，对您同样有帮助:
好好学习
1楼 2012-07-23 15:38:27
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
回帖置顶 ( 共有1个 )

cuishuai1234

铜虫 (初入文坛)
cuishuai1234: 回帖置顶 2012-07-24 14:34:11
引用回帖:
4楼: Originally posted by cxyuan09 at 2012-07-23 16:40:17
我分析问题还是出在scf关键词的使用上
那里面的SSD选项是否适合你当前的体系，我不太清楚，因为以前没接触过稀土，如果不是必须，不妨去掉试试看。
dsymm可以继续用，仍然不收敛时可用symm，有时更管用
如果还不 ...

之前用的xqc计算，不收敛后改为qc，之后提示说要用ssd，这不还是不收敛。这个分子以前用b3lyp方法算出来了，该为bp86后就怎么也不收敛，应该不是自旋多重度的事，请问大侠还有其他的解决办法吗？谢谢啦！
赞一下 回复此楼
好好学习
6楼2012-07-23 16:59:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
普通回帖

cxyuan09

金虫 (正式写手)

【答案】应助回帖

感谢参与，应助指数 +1
如果输入文件只有这么多，那么是你的基组定义有问题
那个重金属没有基组和赝势的定义
赞一下 回复此楼
2楼2012-07-23 15:52:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cuishuai1234

铜虫 (初入文坛)
引用回帖:
2楼: Originally posted by cxyuan09 at 2012-07-23 15:52:54
如果输入文件只有这么多，那么是你的基组定义有问题
那个重金属没有基组和赝势的定义

这是输入文件的一部分
%ch=6-DMDHOPDA+Eu.chk
%mem=400mw
%nprocshared=8
# opt=maxcycle=1000 bp86/genecp iop(5/13=1) scf=(maxcycle=500,VShift,dsymm,ssd)

Title Card Required

3 7
  C                  0.02824629   -1.50935937   -2.38659573
  C                  0.58832866   -2.25077244   -1.16006815
  H                  1.69362728   -2.24339362   -1.19371637
  H                  0.27797084   -3.31077503   -1.18734229
  C                  0.58832866   -2.25077244    1.16006815
  C                  0.02824629   -1.50935937    2.38659573
  H                  0.27797084   -3.31077503    1.18734229
  H                  1.69362728   -2.24339362    1.19371637
  O                  0.11979949   -1.61844052    0.00000000
  O                 -0.64113220   -0.49537049   -2.26910628
  O                 -0.64113220   -0.49537049    2.26910628
  C                  1.11123016   -3.26015963   -3.82870459
  H                  1.74870742   -3.48026131   -2.97136103
  H                  1.76953707   -3.11913473   -4.69329094
  C                 -0.25220632   -1.42145174   -4.80228175
  C                  0.66874977   -0.37244734   -5.43970730
  H                 -1.19083870   -0.95147709   -4.49967054
  H                 -0.48704693   -2.21409366   -5.52518169
  C                  0.05987591    0.24988516   -6.70206976
  H                  0.86501699    0.40785303   -4.69420754
  H                  1.63869663   -0.82937253   -5.68526895
  C                  0.95827896    1.30982429   -7.35345308
  H                 -0.15537462   -0.54228357   -7.43534819
  H                 -0.91000369    0.70247779   -6.44954140
  H                  1.17485359    2.10011369   -6.61997366
  H                  1.92850575    0.85735503   -7.60795602
  N                  0.32459040   -2.05636186   -3.61372972
  C                  0.34708906    1.93865952   -8.61169706
  C                  1.24443524    2.99699341   -9.26675004
  H                 -0.62156810    2.39349521   -8.35632770
  H                  0.12644389    1.14766882   -9.34431784
  C                  0.62504162    3.62282327  -10.52038291
  H                  1.46579541    3.78635311   -8.53440420
  H                  2.21180932    2.54226568   -9.52423518
  H                  0.42380787    2.86255794  -11.28534712
  H                 -0.32661368    4.11647291  -10.28787963
  C                  1.11123016   -3.26015963    3.82870459
  H                  1.76953707   -3.11913473    4.69329094
  H                  1.74870742   -3.48026131    2.97136103
  C                 -0.25220632   -1.42145174    4.80228175
  C                  0.66874977   -0.37244734    5.43970730
  H                 -0.48704693   -2.21409366    5.52518169
  H                 -1.19083870   -0.95147709    4.49967054
  C                  0.05987591    0.24988516    6.70206976
  H                  1.63869663   -0.82937253    5.68526895
  H                  0.86501699    0.40785303    4.69420754
  C                  0.95827896    1.30982429    7.35345308
  H                 -0.91000369    0.70247779    6.44954140
  H                 -0.15537462   -0.54228357    7.43534819
  H                  1.92850575    0.85735503    7.60795602
  H                  1.17485359    2.10011369    6.61997366
  C                  0.34708906    1.93865952    8.61169706
  C                  1.24443524    2.99699341    9.26675004
  H                  0.12644389    1.14766882    9.34431784
  H                 -0.62156810    2.39349521    8.35632770
  C                  0.62504162    3.62282327   10.52038291
  H                  2.21180932    2.54226568    9.52423518
  H                  1.46579541    3.78635311    8.53440420
  H                 -0.32661368    4.11647291   10.28787963
  H                  0.42380787    2.86255794   11.28534712
  N                  0.32459040   -2.05636186    3.61372972
  H                  0.47707969   -4.13771332   -4.02479650
  H                  0.47707969   -4.13771332    4.02479650
  H                  1.28960930    4.37324418  -10.96369520
  H                  1.28960930    4.37324418   10.96369520
  Eu                -1.37449859    0.39828346    0.00000000

H C N O 0
3-21g
****
Eu 0
S 5 1.00
70059.4200000 0.000097
10776.2350000 0.000707
2482.4901000 0.003362
702.1526000 0.010206
216.7926000 0.016597
S 1 1.00
41.4950000 1.
S 1 1.00
29.4969000 1.
S 1 1.00
15.0309000 1.
S 1 1.00
3.7772000 1.
S 1 1.00
1.9219000 1.
S 1 1.00
0.7131000 1.
S 1 1.00
0.3097000 1.
S 1 1.00
0.0549000 1.
S 1 1.00
0.0227000 1.
P 6 1.00
4058.0207000 0.000048
1023.3259000 0.000357
339.3852000 0.001687
122.3485000 0.005693
29.0558000 0.072398
20.6269000 -0.046807
P 1 1.00
14.6179000 1.
P 1 1.00
5.4444000 1.
P 1 1.00
2.7414000 1.
P 1 1.00
1.2283000 1.
P 1 1.00
0.5769000 1.
P 1 1.00
0.2489000 1.
P 1 1.00
0.0800000 1.
D 6 1.00
389.0841000 0.000534
117.7602000 0.004258
44.9869000 0.017160
22.8599000 0.012427
7.8331000 0.270416
4.1134000 0.445272
D 1 1.00
2.0999000 1.
D 1 1.00
1.0043000 1.
D 1 1.00
0.3422000 1.
D 1 1.00
0.1051000 1.
F 5 1.00
128.2341000 0.003424
46.2685000 0.034313
20.7366000 0.123194
9.4479000 0.251824
4.3149000 0.351292
F 1 1.00
1.9121000 1.
F 1 1.00
0.7897000 1.
F 1 1.00
0.2859000 1.
****

Eu 0
ECP28MWB 5 28
H-Komponente
1
2 1.000000 0.000000
S-H
1
2 23.471384 607.659331
P-H
1
2 16.772479 264.385476
D-H
1
2 13.981343 115.381375
F-H
1
2 23.962888 -49.400794
G-H
1
2 21.232458 -26.748273

--Link1--
%ch=6-DMDHOPDA+Eu.chk
%mem=400mw
%nprocshared=8
# opt=maxcycle=1000 bp86/genecp iop(5/13=1) scf=(maxcycle=500,VShift,dsymm,ssd)

Title Card Required

3 7

H C N O 0
6-31g(d)
****
Eu 0
S 5 1.00
70059.4200000 0.000097
10776.2350000 0.000707
2482.4901000 0.003362
702.1526000 0.010206
216.7926000 0.016597
S 1 1.00
41.4950000 1.
S 1 1.00
29.4969000 1.
S 1 1.00
15.0309000 1.
S 1 1.00
3.7772000 1.
S 1 1.00
1.9219000 1.
S 1 1.00
0.7131000 1.
S 1 1.00
0.3097000 1.
S 1 1.00
0.0549000 1.
S 1 1.00
0.0227000 1.
P 6 1.00
4058.0207000 0.000048
1023.3259000 0.000357
339.3852000 0.001687
122.3485000 0.005693
29.0558000 0.072398
20.6269000 -0.046807
P 1 1.00
14.6179000 1.
P 1 1.00
5.4444000 1.
P 1 1.00
2.7414000 1.
P 1 1.00
1.2283000 1.
P 1 1.00
0.5769000 1.
P 1 1.00
0.2489000 1.
P 1 1.00
0.0800000 1.
D 6 1.00
389.0841000 0.000534
117.7602000 0.004258
44.9869000 0.017160
22.8599000 0.012427
7.8331000 0.270416
4.1134000 0.445272
D 1 1.00
2.0999000 1.
D 1 1.00
1.0043000 1.
D 1 1.00
0.3422000 1.
D 1 1.00
0.1051000 1.
F 5 1.00
128.2341000 0.003424
46.2685000 0.034313
20.7366000 0.123194
9.4479000 0.251824
4.3149000 0.351292
F 1 1.00
1.9121000 1.
F 1 1.00
0.7897000 1.
F 1 1.00
0.2859000 1.
****

Eu 0
ECP28MWB 5 28
H-Komponente
1
2 1.000000 0.000000
S-H
1
2 23.471384 607.659331
P-H
1
2 16.772479 264.385476
D-H
1
2 13.981343 115.381375
F-H
1
2 23.962888 -49.400794
G-H
1
2 21.232458 -26.748273



--Link1--
%chk=6-DMDHOPDA+Eu.chk
%mem=400mw
%nprocshared=8
# bp86/genecp guess=read geom=check scf(maxcycle=500,VShift,dsymm,ssd) iop(5/13=1) freq=noraman

Title Card Required

3 7

H C N O 0
6-311g(d,p)
****
Eu 0
S 5 1.00
70059.4200000 0.000097
10776.2350000 0.000707
2482.4901000 0.003362
702.1526000 0.010206
216.7926000 0.016597
S 1 1.00
41.4950000 1.
S 1 1.00
29.4969000 1.
S 1 1.00
15.0309000 1.
S 1 1.00
3.7772000 1.
S 1 1.00
1.9219000 1.
S 1 1.00
0.7131000 1.
S 1 1.00
0.3097000 1.
S 1 1.00
0.0549000 1.
S 1 1.00
0.0227000 1.
P 6 1.00
4058.0207000 0.000048
1023.3259000 0.000357
339.3852000 0.001687
122.3485000 0.005693
29.0558000 0.072398
20.6269000 -0.046807
P 1 1.00
14.6179000 1.
P 1 1.00
5.4444000 1.
P 1 1.00
2.7414000 1.
P 1 1.00
1.2283000 1.
P 1 1.00
0.5769000 1.
P 1 1.00
0.2489000 1.
P 1 1.00
0.0800000 1.
D 6 1.00
389.0841000 0.000534
117.7602000 0.004258
44.9869000 0.017160
22.8599000 0.012427
7.8331000 0.270416
4.1134000 0.445272
D 1 1.00
2.0999000 1.
D 1 1.00
1.0043000 1.
D 1 1.00
0.3422000 1.
D 1 1.00
0.1051000 1.
F 5 1.00
128.2341000 0.003424
46.2685000 0.034313
20.7366000 0.123194
9.4479000 0.251824
4.3149000 0.351292
F 1 1.00
1.9121000 1.
F 1 1.00
0.7897000 1.
F 1 1.00
0.2859000 1.
****

Eu 0
ECP28MWB 5 28
H-Komponente
1
2 1.000000 0.000000
S-H
1
2 23.471384 607.659331
P-H
1
2 16.772479 264.385476
D-H
1
2 13.981343 115.381375
F-H
1
2 23.962888 -49.400794
G-H
1
2 21.232458 -26.748273
赞一下 回复此楼
好好学习
3楼2012-07-23 15:59:14
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cxyuan09

金虫 (正式写手)

【答案】应助回帖

我分析问题还是出在scf关键词的使用上
那里面的SSD选项是否适合你当前的体系，我不太清楚，因为以前没接触过稀土，如果不是必须，不妨去掉试试看。
dsymm可以继续用，仍然不收敛时可用symm，有时更管用
如果还不管用，可在scf关键词的最后加个xqc选项。
赞一下 回复此楼
4楼2012-07-23 16:40:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cxyuan09

金虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★
cuishuai1234: 金币+5, ★有帮助 2012-07-23 20:33:42
另外，你的重金属只有3配位，
如何选择自选多重度也是一个值得研究的问题，
这个选不对，也容易导致scf不收敛
赞一下 回复此楼
5楼2012-07-23 16:48:37
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cuishuai1234 的主题更新
61/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[教师之家] 售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急 +3 XLGfUIdOmEAO 2026-08-27 4/200 2026-08-27 22:26 by 3RAFyYXYSwGS
[基金申请] 为什么资助数各大高校都创新高，自己申请怎么就这么难 +8 Kittylucky 2026-08-27 8/400 2026-08-27 18:33 by raiy123
[教师之家] 导师吐槽：我怎么摊上了这么个极品研究生！ +7 苏东坡二世 2026-08-23 7/350 2026-08-27 18:14 by 瞬息宇宙
[教师之家] 售SCI一区文章，我:8.O.551.O.5.4,科目全,可伽急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 04:33 by Ie9AyIAvGbvs
[硕博家园] 售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急 +3 LIbGuocjEEYw 2026-08-26 4/200 2026-08-27 02:01 by Ie9AyIAvGbvs
[基金申请] 范进中举一文的中心思想 +9 炎黄贵胄 2026-08-22 10/500 2026-08-26 15:40 by semaglutide
[基金申请] 2026国自然函评费到账 +22 羊腰板 2026-08-21 25/1250 2026-08-26 15:00 by zuocuiping
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件，有关系吗？ +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 能否退出参与的面上项目解除限项 +23 koalala 2026-08-24 26/1300 2026-08-26 14:29 by 宝贝虫子
[基金申请] 哪位高人中了，把查询到的截图贴出来让我看看，让我长长见识 +4 yuleib84 2026-08-26 5/250 2026-08-26 13:40 by yuleib84
[基金申请] 系统进不去 +4 yanglien 2026-08-26 5/250 2026-08-26 11:10 by wenfengw83
[基金申请] 项目信息和经费信息在系统里都可以看到了 +6 wittyboy 2026-08-26 14/700 2026-08-26 10:55 by wittyboy
[基金申请] 国合现在查不到了吗？ +10 chengyan1220 2026-08-24 20/1000 2026-08-26 08:57 by peasantsprig
[基金申请] 在坚冰还盖着北海的时候，我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构，以效率低为荣，以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
[基金申请] 如果此刻你正在为国基感到焦虑，不妨来听听这首《基金之外》 +8 scalable 2026-08-24 8/400 2026-08-25 12:52 by jnhyjjm
[基金申请] 没有任何消息-是不是就凉了 +9 图啦图啦 2026-08-24 10/500 2026-08-25 11:59 by 南海小哥
[基金申请] 建议基金发布提前给出明确的时间点 +13 kulium 2026-08-21 16/800 2026-08-24 16:27 by superceng
[基金申请] 看来今天不会放榜了？ +8 chengyan1220 2026-08-21 11/550 2026-08-21 17:52 by dcqxinyang
[基金申请] 基金啊基金 +4 longfie172 2026-08-20 4/200 2026-08-21 08:58 by mark mao
信息提示
关闭
请填处理意见
关闭