24小时热门版块排行榜    

查看: 3032  |  回复: 4
当前只显示满足指定条件的回帖,点击这里查看本话题的所有回帖

liuzhengjun0427

至尊木虫 (著名写手)

[求助] Excessive mixing of frozen core and valence orbitals.

本人用B3LYP/6-311++G**正常优化完一个分子,然后用MP2/6-311++G**算单点时出现如下错误:
Largest valence mixing into a core orbital is  4.91D-01
Largest core mixing into a valence orbital is  4.87D-01
Excessive mixing of frozen core and valence orbitals.
Error termination via Lnk1e in /home/chem/gaussian/g03/l801.exe at Mon Sep  5 13:18:29 2011.


求高手解答!!!
回复此楼

» 猜你喜欢

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

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   关注TA 给TA发消息 送TA红花 TA的回帖

liuzhengjun0427

至尊木虫 (著名写手)

引用回帖:
2楼: Originally posted by beefly at 2011-09-06 12:48:00:
说明gaussian默认设置的冻芯轨道可能不合理。

在mp2计算中搜索“NFC”(也就是冻结轨道的个数)的数值,然后作一个HF计算,找到轨道输出部分,看NFC和NFC+1对应的轨道是否有较大能隙。如果比较接近,需要找一个 ...

我的输出文件中没搜到您说的“NFC”,以下是输出文件的最后一部分: Cycle  17  Pass 1  IDiag  1:
RMSU=  5.22D-09    CP:  7.79D-01  6.96D-01  8.55D-01  1.11D+00  1.19D+00
                    CP:  1.63D+00  1.41D+00  1.45D+00  1.42D+00  1.64D+00
                    CP:  1.81D+00  1.48D+00  1.64D+00  1.36D+00  1.34D+00
                    CP:  1.51D+00
E= -1354.40044392809     Delta-E=       -0.000000000044 Rises=F Damp=F
DIIS: error= 2.67D-08 at cycle  17 NSaved=  17.
NSaved=17 IEnMin=17 EnMin= -1354.40044392809     IErMin=17 ErrMin= 2.67D-08
ErrMax= 2.67D-08 EMaxC= 1.00D-01 BMatC= 1.01D-13 BMatP= 4.13D-13
IDIUse=1 WtCom= 1.00D+00 WtEn= 0.00D+00
Coeff-Com: -0.100D-06 0.110D-05 0.375D-05 0.101D-04-0.968D-05-0.988D-04
Coeff-Com: -0.817D-04 0.624D-03 0.101D-02-0.291D-02-0.463D-02 0.127D-01
Coeff-Com:  0.339D-01-0.599D-01-0.231D+00 0.224D+00 0.103D+01
Coeff:     -0.100D-06 0.110D-05 0.375D-05 0.101D-04-0.968D-05-0.988D-04
Coeff:     -0.817D-04 0.624D-03 0.101D-02-0.291D-02-0.463D-02 0.127D-01
Coeff:      0.339D-01-0.599D-01-0.231D+00 0.224D+00 0.103D+01
Gap=     0.323 Goal=   None    Shift=    0.000
RMSDP=5.57D-09 MaxDP=3.79D-07 DE=-4.37D-11 OVMax= 4.41D-07

SCF Done:  E(RHF) =  -1354.40044393     A.U. after   17 cycles
             Convg  =    0.5567D-08             -V/T =  2.0009
             S**2   =   0.0000
KE= 1.353223733552D+03 PE=-6.299120414074D+03 EE= 2.032819316259D+03
Leave Link  502 at Mon Sep  5 13:18:23 2011, MaxMem=  262144000 cpu:    6879.0
(Enter /home/chem/gaussian/g03/l801.exe)
ExpMin= 7.60D-03 ExpMax= 9.34D+04 ExpMxC= 3.17D+03 IAcc=3 IRadAn=         5 AccDes= 0.00D+00
HarFok:  IExCor= 205 AccDes= 0.00D+00 IRadAn=         5 IDoV=1
ScaDFX=  1.000000  1.000000  1.000000  1.000000
Largest valence mixing into a core orbital is  4.91D-01
Largest core mixing into a valence orbital is  4.87D-01
Excessive mixing of frozen core and valence orbitals.
Error termination via Lnk1e in /home/chem/gaussian/g03/l801.exe at Mon Sep  5 13:18:29 2011.
Job cpu time:  0 days  1 hours 56 minutes 47.6 seconds.
File lengths (MBytes):  RWF=    152 Int=      0 D2E=      0 Chk=      7 Scr=      1
3楼2011-09-06 17:04:34
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

beefly

专家顾问 (职业作家)

地沟油冶炼专家

说明gaussian默认设置的冻芯轨道可能不合理。

在mp2计算中搜索“NFC”(也就是冻结轨道的个数)的数值,然后作一个HF计算,找到轨道输出部分,看NFC和NFC+1对应的轨道是否有较大能隙。如果比较接近,需要找一个合理的分界线,然后设定mp2(window=n),其中的n=新的NFC+1
beefly《西太平洋大学现代英汉词典》[bi:fli]牛肉一般地
2楼2011-09-06 12:48:00
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

liuzhengjun0427

至尊木虫 (著名写手)

以前重来没遇到过这种情况,没想到算个单点还会出错。
4楼2011-09-06 17:07:09
已阅   关注TA 给TA发消息 送TA红花 TA的回帖

yflchx

金虫 (正式写手)

引用回帖:
1楼: Originally posted by liuzhengjun0427 at 2011-09-06 08:58:14:
本人用B3LYP/6-311++G**正常优化完一个分子,然后用MP2/6-311++G**算单点时出现如下错误:
Largest valence mixing into a core orbital is  4.91D-01
...

用MP2(FULL)吧。
5楼2011-09-06 18:45:28
已阅   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见