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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » gauss计算单点能L508出错
41/1返回列表
查看: 3115  |  回复: 3
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

风雨双龙山

新虫 (正式写手)

[求助] gauss计算单点能L508出错

SCF在小基组优化结构后用大基组求单点能出现L508错误，如何解决
输入文件如下：
%chk=Co-ON-sp-b.chk
%mem=6GB
%NProcShared=8
%rwf=1,19900000mb,2,19900000mb,3,19900000mb,4,19900000mb
#p b3lyp gen geom=connectivity pseudo=read scf=xqc test

Title Card Required

-5 3
O                  0.87774000    1.27069700    0.18546400
O                 -0.79101400   -0.00064000   -1.25088300
O                 -1.03098500   -0.00438500    1.27757800
O                 -1.90800600   -1.40633700    3.38733400
W                 -2.70878500   -1.73899100    1.78608200
O                 -1.26938400   -2.95344000    1.25807100
W                  0.11007900   -3.58714100    0.09027900
O                 -1.01167200   -2.95182900   -1.34545300
W                 -2.37196300   -1.73642900   -1.98509400
O                 -2.96541000   -1.64310600   -0.21275700
O                  0.88515100   -1.26615300    0.18359500
O                 -3.17559000   -0.00940400   -2.35749400
W                 -2.38461400    1.72362100   -1.98125300
O                 -1.18409200    1.32311200   -3.47154200
W                  0.23153600    0.00434600   -3.48519000
O                 -1.17434600   -1.32469600   -3.47427600
O                 -2.97830200    1.62063900   -0.21000300
W                 -2.72357600    1.71248600    1.79023400
O                 -1.91960500    1.38340500    3.39002300
O                 -4.01658400    2.78480600    2.22154900
O                 -1.29298900    2.94117700    1.26528900
W                  0.08137200    3.58785600    0.09976500
O                  1.37394900    3.45913400    1.51696400
W                  2.32758300    1.84982500    2.11981400
O                  2.88766100    0.00893000    1.96056100
W                  2.34144300   -1.83731500    2.11600500
O                  1.40034700   -3.45257200    1.50975300
O                 -3.60570000   -0.01671800    1.83072300
O                 -3.52525100    2.78389100   -2.74154500
O                 -1.03421100    2.94767200   -1.33882300
O                 -0.21481000    5.29352100    0.07384500
O                  1.58471700    3.45230000   -1.12005700
W                  2.70211600    1.88921300   -1.37705000
O                  3.34329000    2.12596600    0.36443000
O                  3.97740600    2.54693900   -2.34813000
O                  1.39604700    1.31146600   -2.71364700
O                  3.15774700    0.01395500   -1.31373700
W                  2.71747800   -1.86366900   -1.38080900
O                  3.35929400   -2.10071700    0.36052200
O                  3.49080600    2.53134300    3.21519100
O                  0.95472300    1.36339100    3.18231900
O                  0.68640000    0.00775000   -5.15756600
O                  1.40558800   -1.29596500   -2.71652300
O                 -3.50544200   -2.80329200   -2.74675800
O                 -3.99329300   -2.82230700    2.21554400
O                  0.96530000   -1.36227800    3.17956500
O                  3.99790400   -2.50967500   -2.35287600
O                  1.61226000   -3.43698500   -1.12794100
O                 -0.17303900   -5.29498800    0.06065900
O                  3.50972100   -2.51157600    3.21036200
P                 -0.01611600   -0.00008200    0.09883600
O                 -0.32837800   -0.00738700    5.43744000
N                 -1.21413700   -0.01268300    6.24219200
Co                -0.49186600   -0.00579400    3.40229900

1 51 1.5
2 51 1.5
3 51 1.5
4 5 1.0
5 6 1.0 10 1.0 28 1.0 45 1.0
6 7 1.0
7 8 1.0 27 1.0 48 1.0 49 1.0
8 9 1.0
9 10 1.0 12 1.0 16 1.0 44 1.0
10
11 51 1.5
12 13 1.0
13 14 1.0 17 1.0 29 1.0 30 1.0
14 15 1.0
15 16 1.0 36 1.0 42 1.0 43 1.0
16
17 18 1.0
18 19 1.0 20 1.0 21 1.0 28 1.0
19
20
21 22 1.0
22 23 1.0 30 1.0 31 1.0 32 1.0
23 24 1.0
24 25 1.0 34 1.0 40 1.0 41 1.0
25 26 1.0
26 27 1.0 39 1.0 46 1.0 50 1.0
27
28
29
30
31
32 33 1.0
33 34 1.0 35 1.0 36 1.0 37 1.0
34
35
36
37 38 1.0
38 39 1.0 43 1.0 47 1.0 48 1.0
39
40
41
42
43
44
45
46
47
48
49
50
51
52 53 2.0
53
54

N O P 0
6-31+g*
****
W Co 0
sdd
****

W Co 0
sdd

结果显示
FoF2E skips out because all densities are zero.
     ILin=16 X=6.104D-05 Y=-4.291579100166D+03 DE= 6.93D-08 F= -7.85D-10
The polynomial fit failed.  Using point  1.
     A contracting polynomial of degree 16 produced  0.0000
Search did not lower the energy significantly.
No lower point found -- run aborted.
Error termination via Lnk1e in /share/apps/g09.D01/g09/l508.exe at Fri Nov 25 02:09:58 2016.
Job cpu time:      48 days 19 hours  9 minutes 28.0 seconds.
File lengths (MBytes):  RWF=   3385 Int=      0 D2E=      0 Chk=     35 Scr=      1
跪求指点
回复此楼

» 猜你喜欢
1楼 2016-12-13 11:01:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

paramecium86

版主 (著名写手)
如果认为结构没有什么大的出入 而且SCF用了qc还不收敛 可以试试几个方法  1用了带diffuse function的基组 可以试试SCF(novaracc,noincfock) int=acc2e=12  改改格点积分的精度。2 包含有过渡金属 那就试试 SCF=vshift=xxx  3 还可以试试放松收敛的标准 SCF=conver=N  默认可能是8  试一试放宽到6
赞一下(1人) 回复此楼
2楼2016-12-14 20:36:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

风雨双龙山

新虫 (正式写手)
引用回帖:
2楼: Originally posted by paramecium86 at 2016-12-14 20:36:46
如果认为结构没有什么大的出入 而且SCF用了qc还不收敛 可以试试几个方法  1用了带diffuse function的基组 可以试试SCF(novaracc,noincfock) int=acc2e=12  改改格点积分的精度。2 包含有过渡金属 那就试试 SCF=vshi ...

多谢指点
回复此楼
3楼2016-12-15 22:04:16
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

791618275

捐助贵宾 (正式写手)
引用回帖:
3楼: Originally posted by 风雨双龙山 at 2016-12-15 02:04:16
多谢指点...

同学，你找到哪错了么
赞一下 回复此楼
4楼2017-09-03 10:36:05
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 风雨双龙山 的主题更新
41/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭