24小时热门版块排行榜    

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

ww1987

铁杆木虫 (小有名气)

[求助] 计算出错 这是为什么

请帮我看看 用B3LYP -----6-31+G(d,p)水平算mg金属卟啉 结果一开始的双电子积分就不给我算 这是为什么  为什么啊

Rotational constants (GHZ):      0.2649516      0.2630332      0.1320100
   539 basis functions,   908 primitive gaussians,   539 cartesian basis functions
    86 alpha electrons       86 beta electrons
       nuclear repulsion energy      2110.3944871538 Hartrees.
NAtoms=   37 NActive=   37 NUniq=   37 SFac= 7.50D-01 NAtFMM=   80 NAOKFM=F Big=F
Harris functional with IExCor=  402 diagonalized for initial guess.
ExpMin= 1.46D-02 ExpMax= 1.17D+04 ExpMxC= 1.76D+03 IAcc=3 IRadAn=         0 AccDes= 0.00D+00
HarFok:  IExCor=  402 AccDes= 0.00D+00 IRadAn=         0 IDoV= 1
ScaDFX=  1.000000  1.000000  1.000000  1.000000
Spurious integrated density or basis function:
NE=  172 NElCor=    0 El error=9.99D-04 rel=5.81D-06 Tolerance=1.00D-03
Shell   161     absolute error=1.32D-02              Tolerance=1.20D-02
Shell   161       signed error=1.32D-02              Tolerance=1.00D-01
Inaccurate quadrature in CalDSu.
Error termination via Lnk1e in C:\G09W\l401.exe at Wed May 04 16:02:37 2011.
Job cpu time:  0 days  0 hours  1 minutes  5.0 seconds.
File lengths (MBytes):  RWF=     26 Int=      0 D2E=      0 Chk=      1 Scr=      1
回复此楼
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

ww1987

铁杆木虫 (小有名气)

我算到这一步又不行了  为什么啊
Closed shell SCF:
Requested convergence on RMS density matrix=1.00D-08 within 128 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.
No special actions if energy rises.
Using DIIS extrapolation, IDIIS=  1040.
Integral symmetry usage will be decided dynamically.
       320780 words used for storage of precomputed grid.
IEnd=     1217065 IEndB=     1217065 NGot=   157286400 MDV=   156377036
LenX=   156377036 LenY=   156085976
Fock matrices will be formed incrementally for  20 cycles.

Cycle   1  Pass 1  IDiag  1:
FoFCou: FMM=F IPFlag=           0 FMFlag=      100000 FMFlg1=           0
         NFxFlg=           0 DoJE=F BraDBF=F KetDBF=F FulRan=T
         Omega=  0.000000  0.000000  1.000000  0.000000  0.000000 ICntrl=       0 IOpCl=  0
         NMat0=    1 NMatS0=    1 NMatT0=    0 NMatD0=    1 NMtDS0=    0 NMtDT0=    0
         I1Cent=           0 NGrid=           0.
Symmetry not used in FoFCou.
Spurious integrated density or basis function:
NE=  172 NElCor=    0 El error=1.05D-03 rel=6.09D-06 Tolerance=1.00D-03
Shell   125     absolute error=1.27D-02              Tolerance=1.20D-02
Shell   125       signed error=1.27D-02              Tolerance=1.00D-01
Inaccurate quadrature in CalDSu.
Error termination via Lnk1e in C:\G09W\l502.exe at Wed May 04 22:38:03 2011.
Job cpu time:  0 days  2 hours 12 minutes 11.0 seconds.
File lengths (MBytes):  RWF=    178 Int=      0 D2E=      0 Chk=     17 Scr=      1
5楼2011-05-04 22:48:30
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
查看全部 5 个回答

gmy1990

荣誉版主 (著名写手)

优秀版主优秀版主

【答案】应助回帖

ww1987(金币+5): 2011-05-04 23:01:34
ww1987(金币+5): 2011-05-05 10:07:52
引用回帖:
Originally posted by ww1987 at 2011-05-04 20:14:19:
请帮我看看 用B3LYP -----6-31+G(d,p)水平算mg金属卟啉 结果一开始的双电子积分就不给我算 这是为什么  为什么啊

Rotational constants (GHZ):      0.2649516      0.2630332      0.1320100
   539 basi ...

镁卟啉我也做过,自己试着找原因,一个简单的优化哪有这么多问题,自己尝试自己学,不是我泼冷水
2楼2011-05-04 21:13:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

求学之子

新虫 (小有名气)

【答案】应助回帖

是不是你没有考虑电荷和多重度问题啊
3楼2011-05-04 21:25:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

求学之子

新虫 (小有名气)

【答案】应助回帖

是不是你没有考虑电荷和多重度问题啊
4楼2011-05-04 21:25:17
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
信息提示
请填处理意见