请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 385|回复: 12

[Gaussian/gview] MP2能量严重异常

[复制链接]

12

帖子

0

威望

42

eV
积分
54

Level 2 能力者

发表于 2018-5-29 12:59:52 | 显示全部楼层 |阅读模式
本帖最后由 keshsong 于 2018-5-30 09:58 编辑

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)]\\@

175

帖子

0

威望

1447

eV
积分
1622

Level 5 (御坂)

发表于 2018-5-29 13:56:14 | 显示全部楼层
你这实际上是做的UHF和UMP2的计算。最好分两步,第一步做UHF,并检查波函数稳定性,必须得到稳定的波函数才能继续;第二步做UMP2。如
#p uhf/basis stable=opt

Title

1 2
xxx


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

12

帖子

0

威望

42

eV
积分
54

Level 2 能力者

 楼主| 发表于 2018-5-29 14:32:13 | 显示全部楼层
zjxitcc 发表于 2018-5-29 13:56
你这实际上是做的UHF和UMP2的计算。最好分两步,第一步做UHF,并检查波函数稳定性,必须得到稳定的波函数才 ...

谢谢指教。

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

175

帖子

0

威望

1447

eV
积分
1622

Level 5 (御坂)

发表于 2018-5-29 14:40:17 | 显示全部楼层
keshsong 发表于 2018-5-29 14:32
谢谢指教。

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

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

12

帖子

0

威望

42

eV
积分
54

Level 2 能力者

 楼主| 发表于 2018-5-29 16:13:25 | 显示全部楼层
本帖最后由 keshsong 于 2018-5-29 16:38 编辑

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

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)]\\@





1万

帖子

25

威望

1万

eV
积分
32295

管理员

公社社长

发表于 2018-5-29 19:00:16 | 显示全部楼层
上传压缩后的输出文件,一看就清楚了
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

175

帖子

0

威望

1447

eV
积分
1622

Level 5 (御坂)

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

151

帖子

0

威望

1149

eV
积分
1300

Level 4 (黑子)

发表于 2018-5-29 20:25:36 | 显示全部楼层
本帖最后由 Jasminer 于 2018-5-30 12:25 编辑

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

12

帖子

0

威望

42

eV
积分
54

Level 2 能力者

 楼主| 发表于 2018-5-30 09:58:05 | 显示全部楼层
sobereva 发表于 2018-5-29 19:00
上传压缩后的输出文件,一看就清楚了

输出文件已上传,请帮忙分析解决,谢谢!

cation-DOX.log

54.41 KB, 下载次数: 7

213

帖子

0

威望

2603

eV
积分
2816

Level 5 (御坂)

发表于 2018-5-30 10:50:57 | 显示全部楼层
keshsong 发表于 2018-5-30 09:58
输出文件已上传,请帮忙分析解决,谢谢!

%nprocshared=28

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

12

帖子

0

威望

42

eV
积分
54

Level 2 能力者

 楼主| 发表于 2018-5-30 11:16:57 | 显示全部楼层
yflchx 发表于 2018-5-30 10:50
%nprocshared=28

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

谢谢,是这个问题。

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

213

帖子

0

威望

2603

eV
积分
2816

Level 5 (御坂)

发表于 2018-5-30 17:20:35 | 显示全部楼层
keshsong 发表于 2018-5-30 11:16
谢谢,是这个问题。

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

设置28的依据是什么?

测试过加速效果吗?

1万

帖子

25

威望

1万

eV
积分
32295

管理员

公社社长

发表于 2018-5-30 18:10:42 | 显示全部楼层
去掉弥散函数再试,当前给出的E2能量巨大E2 =    -0.2499925509D+03,极有可能是数值不稳定性造成的

另外,对开壳层体系用常规UMP2去计算是糟糕透顶的做法,很容易被内行审稿人强烈批评要求用DFT重算。MP2对绝大多数类型的问题早过时了,尤其是算开壳层更是遭人诟病,比DFT贵得多还差得多。你这样的体系用M06-2X/def2-TZVP结果明显强得多得多。
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-7-18 03:10 , Processed in 0.204273 second(s), 25 queries .

快速回复 返回顶部 返回列表