MP2能量严重异常

keshsong

1，4-二氧6环，MP2做正离子的单点能， 结果MP2能量降低100多原子单位。
请教一下是什么问题？
输出文件最后HF=-305.536...\MP2=-441.876...


(Enter /opt/software/g09/l9999.exe)
1\1\GINC-CU23\SP\UMP2-FC\6-311+G(d,p)\C4H8O2(1+,2)\\28-May-2018\
0\\#p mp2/6-311+g(d,p)\\Title Card Required\\1,2\C,0,1.399779,0.000046
,-0.225186\C,0,0.623831,1.23763,0.212373\C,0,-1.33482,-0.000051,0.1789
83\C,0,0.623919,-1.237589,0.212377\H,0,0.604877,1.304713,1.312333\H,0,
1.055971,2.156859,-0.186699\H,0,1.489469,0.000049,-1.315466\H,0,2.4047
38,0.00009,0.212277\H,0,-1.338281,-0.000047,1.284242\H,0,-2.348673,-0.
000082,-0.217651\H,0,0.604962,-1.304682,1.312335\H,0,1.056105,-2.15679
,-0.186705\O,0,-0.71288,1.173227,-0.28037\O,0,-0.712798,-1.173268,-0.2
80374\\Version=ES64L-G09RevD.01\State=2-A\HF=-305.5362184\MP2=-441.876
4175\PUHF=-305.5495274\PMP2-0=-441.8839109\S2=0.797285\S2-1=0.755408\S
2A=0.751062\RMSD=7.439e-09\PG=C01 [X(C4H8O2)]\\@

zjxitcc

你这实际上是做的UHF和UMP2的计算。最好分两步，第一步做UHF，并检查波函数稳定性，必须得到稳定的波函数才能继续；第二步做UMP2。如
#p uhf/basis stable=opt

Title

1 2
xxx


--Link1--
#p ump2/basis guess=read geom=allcheck

keshsong

zjxitcc 发表于 2018-5-29 13:56
你这实际上是做的UHF和UMP2的计算。最好分两步，第一步做UHF，并检查波函数稳定性，必须得到稳定的波函数才 ...

谢谢指教。

正离子的构型是MP2的优化的二氧六环构型，不知道哪里出问题了。

zjxitcc

keshsong 发表于 2018-5-29 14:32
谢谢指教。

正离子的构型是MP2的优化的二氧六环构型，不知道哪里出问题了。

我上面给出了建议，检测UHF波函数稳定性。若不稳定，做UMP2是无意义的。

keshsong

谢谢，根据建议做了测试，还是异常且数值有变化。

Eigenvectors of the stability matrix:

Eigenvector   1:  2.035-A    Eigenvalue= 0.0755579  <S**2>=0.786
      9B -> 20B        0.16536
     10B -> 20B       -0.13886
     11B -> 20B        0.47817
     12B -> 20B       -0.12648
     14B -> 20B        0.16172
     15B -> 20B       -0.48312
     16B -> 20B       -0.39381
     17B -> 20B        0.41968
     18B -> 20B       -0.10484
     19B -> 20B        0.14478
SavETr:  write IOETrn=   770 NScale= 10 NData=  16 NLR=1 NState=    1 LETran=      28.
The wavefunction is stable under the perturbations considered.
Leave Link  914 at Tue May 29 16:08:46 2018, MaxMem=    33554432 cpu:        46.6
(Enter /opt/software/g09/l508.exe)
The wavefunction is already stable.
Leave Link  508 at Tue May 29 16:08:46 2018, MaxMem=    33554432 cpu:         0.0
(Enter /opt/software/g09/l601.exe)
Copying SCF densities to generalized density rwf, IOpCl= 1 IROHF=0.

(Enter /opt/software/g09/l9999.exe)
1\1\GINC-CU38\SP\UMP2-FC\6-311+G(d,p)\C4H8O2(1+,2)\29-May-2018\
0\\#p mp2/6-311+g(d,p)\\Title Card Required\\1,2\C,0,1.160384,0.733355
,-0.199442\C,0,1.160478,-0.733211,0.199456\C,0,-1.160375,-0.733365,0.1
99443\C,0,-1.160486,0.733207,-0.199434\H,0,2.023222,-1.251806,-0.22668
1\H,0,1.192944,0.81418,-1.296571\H,0,2.023061,1.252054,0.226702\H,0,-2
.023066,-1.252057,-0.226691\H,0,-1.192911,-0.814196,1.296569\H,0,-1.19
3065,0.814051,-1.296557\H,0,-2.02323,1.251775,0.226738\H,0,1.193005,-0
.814054,1.296581\O,0,-0.000088,1.389797,0.296657\O,0,0.000093,-1.38978
,-0.296686\\Version=ES64L-G09RevD.01\State=2-A\HF=-305.5631347\MP2=-45
7.3473781\PUHF=-305.5722407\PMP2-0=-457.3524559\S2=0.781191\S2-1=0.753
302\S2A=0.750556\RMSD=4.090e-09\PG=C01 [X(C4H8O2)]\\@

sobereva

上传压缩后的输出文件，一看就清楚了

zjxitcc

UHF结束之后是否有Warning分子轨道系数超级大？
你这个UHF能量确实降低了，说明一开始的UHF能量-305.5362184是内部不稳定的，拿它算UMP2没有意义。后面这个-305.5631347比较低。
试试看提升积分精度int=acc2e=12，看看是否有什么不同。
（为节约时间，可以利用当前UHF的chk文件，读取其中的轨道，新的UHF将会在几圈之内收敛，不必从零算起）

Jasminer

我试了一下，没出现你说的现象
你换个软件版本试试？

keshsong

sobereva 发表于 2018-5-29 19:00
上传压缩后的输出文件，一看就清楚了

输出文件已上传，请帮忙分析解决，谢谢！

yflchx

keshsong 发表于 2018-5-30 09:58
输出文件已上传，请帮忙分析解决，谢谢！

%nprocshared=28

这个调整一下，甚至先拿1试试。

keshsong

yflchx 发表于 2018-5-30 10:50
%nprocshared=28

这个调整一下，甚至先拿1试试。

谢谢，是这个问题。

搞不懂的是为什么这个分子出错，而有些分子不出错。

yflchx

keshsong 发表于 2018-5-30 11:16
谢谢，是这个问题。

搞不懂的是为什么这个分子出错，而有些分子不出错。

设置28的依据是什么？

测试过加速效果吗？

sobereva

去掉弥散函数再试，当前给出的E2能量巨大E2 =    -0.2499925509D+03，极有可能是数值不稳定性造成的

另外，对开壳层体系用常规UMP2去计算是糟糕透顶的做法，很容易被内行审稿人强烈批评要求用DFT重算。MP2对绝大多数类型的问题早过时了，尤其是算开壳层更是遭人诟病，比DFT贵得多还差得多。你这样的体系用M06-2X/def2-TZVP结果明显强得多得多。