24小时热门版块排行榜    

查看: 918  |  回复: 6

黄金比例

木虫 (正式写手)

[交流] MP2 计算,结果不一致 已有2人参与

今天发现个奇怪的现象。
用#P MP2/cc-pVDZ OPT=CalcFC 算出来的E2 和 EUMP2 与
用#P MP2/cc-pVDZ OPT=CalcFC FREQ=NoRaman 算出来的,频率分析部分算出来的E2 和 EUMP2 不同,这是为什么???
为什么算出来的电子相关能会不一样呢???
计算的时候,到低用哪个???
结果如下:
#P MP2/cc-pVDZ OPT=CalcFC
----------------------------------------
Spin components of T(2) and E(2):
     alpha-alpha T2 =       0.1115891003D-01 E2=     -0.3131962411D-01
     alpha-beta  T2 =       0.8118108037D-01 E2=     -0.2402864251D+00
     beta-beta   T2 =       0.1115891003D-01 E2=     -0.3131962411D-01
The integrals were generated   1 times.
ANorm=    0.1050475559D+01
E2 =    -0.3029256734D+00 EUMP2 =    -0.79537475761735D+02
==========================================================
#P MP2/cc-pVDZ OPT=CalcFC FREQ=NoRaman
---------------------------------------------------------------------
Spin components of T(2) and E(2):
     alpha-alpha T2 =       0.1115891004D-01 E2=     -0.3131962412D-01
     alpha-beta  T2 =       0.8118108046D-01 E2=     -0.2402864253D+00
     beta-beta   T2 =       0.1115891004D-01 E2=     -0.3131962412D-01
ANorm=    0.1485596783D+01
E2 =    -0.3029256735D+00 EUMP2 =    -0.79537475761888D+02
=========================================================
回复此楼

» 猜你喜欢

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

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

fooo

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
本来就不应该一样,如果一样那才怪了呢。EUMP2=E2+EHF
2楼2013-04-28 03:52:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

黄金比例

木虫 (正式写手)

引用回帖:
2楼: Originally posted by fooo at 2013-04-27 07:52:54
本来就不应该一样,如果一样那才怪了呢。EUMP2=E2+EHF

我用mp2进行 opt 和 freq 计算时,优化出来的stationary point 的 EUMP2 能量
和在该stationary point 结构上进行频率分析得出来 EUMP2 能量是不一样的,
用 HF 和 B3LYP 计算时,两部分的能量是一样的,为什么用mp2计算时,就不一样了???

我知道EUMP2=E2+E(HF)
对于同一个点的计算怎么就不一样呢
3楼2013-04-28 07:27:32
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yanrding

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
其实基本上是一样的,属于误差吧。
ANorm这个值查了不少。

另外,
#P MP2/cc-pVDZ OPT=CalcFC FREQ=NoRaman
这个计算对最后一个点应该算了2次,你比较一下这两次是否完全一样。
华北电力大学(北京)数理学院,计算化学,凝聚态物理。欢迎合作!招硕士/博士研究生,招老师。
4楼2013-04-28 08:05:35
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

黄金比例

木虫 (正式写手)

引用回帖:
4楼: Originally posted by yanrding at 2013-04-27 12:05:35
其实基本上是一样的,属于误差吧。
ANorm这个值查了不少。

另外,
#P MP2/cc-pVDZ OPT=CalcFC FREQ=NoRaman
这个计算对最后一个点应该算了2次,你比较一下这两次是否完全一样。

在优化部分,最后两次的EUMP2是不相等的。

而在频率分析部分,之算了两次,两次的EUMP2 是相等的!

这不是误差吧?!!!
5楼2013-04-28 08:29:55
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yanrding

金虫 (正式写手)


小木虫: 金币+0.5, 给个红包,谢谢回帖
最后一个点,怎么算了那么多次?
opt的最后两次当然会不一样,那么不是同一个构型。
freq里面应该只有一次。

opt的最后一个点,在freq里面会再算一次。这样,最后一个点会被算2次。
华北电力大学(北京)数理学院,计算化学,凝聚态物理。欢迎合作!招硕士/博士研究生,招老师。
6楼2013-04-28 08:44:42
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

黄金比例

木虫 (正式写手)

引用回帖:
6楼: Originally posted by yanrding at 2013-04-27 12:44:42
最后一个点,怎么算了那么多次?
opt的最后两次当然会不一样,那么不是同一个构型。
freq里面应该只有一次。

opt的最后一个点,在freq里面会再算一次。这样,最后一个点会被算2次。

真的和你说的不一样。
在freq时,就算了两次,结果一样,
而且与opt的最后一个点的能量不一样!!!
就是搞不懂。
7楼2013-04-28 12:25:54
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 黄金比例 的主题更新
普通表情 高级回复 (可上传附件)
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考博] 东华理工大学化材专业26届硕士博士申请 +6 zlingli 2026-03-13 6/300 2026-03-15 20:00 by ryzcf
[考研] 309求调剂 +4 花与叶@ 2026-03-10 4/200 2026-03-14 21:26 by a不易
[考研] 材料与化工 323 英一+数二+物化,一志愿:哈工大 本人本科双一流 +4 自由的_飞翔 2026-03-13 5/250 2026-03-14 19:39 by hmn_wj
[考研] 一志愿哈工大材料324分求调剂 +5 闫旭东 2026-03-14 5/250 2026-03-14 14:53 by 木瓜膏
[考研] 材料080500调剂求收留 +3 一颗meteor 2026-03-13 3/150 2026-03-14 10:54 by peike
[考研] 化学工程321分求调剂(南京工业,浙江工业) +3 大米饭! 2026-03-09 4/200 2026-03-14 02:34 by JourneyLucky
[考研] 一志愿北京化工大学材料与化工296分求调剂 +16 稻妻小编 2026-03-09 18/900 2026-03-14 02:00 by JourneyLucky
[考研] 调剂 +3 13853210211 2026-03-10 3/150 2026-03-14 00:47 by JourneyLucky
[考研] 312求调剂 +6 陌宸希 2026-03-10 6/300 2026-03-14 00:40 by JourneyLucky
[考研] 材料371求调剂 +9 鳄鱼? 2026-03-11 11/550 2026-03-13 22:53 by JourneyLucky
[考研] 0703,333分求调剂 一志愿郑州大学-物理化学 +3 李魔女斗篷 2026-03-11 3/150 2026-03-13 22:24 by JourneyLucky
[考研] 26调剂/材料/英一数二/总分289/已过A区线 +6 步川酷紫123 2026-03-13 6/300 2026-03-13 21:59 by 星空星月
[考研] 293求调剂 +3 世界首富 2026-03-11 3/150 2026-03-13 16:27 by JourneyLucky
[考研] 0817化学工程与技术考研312分调剂 +3 T123 tt 2026-03-12 3/150 2026-03-13 10:49 by houyaoxu
[考研] 290求调剂 +3 ADT 2026-03-13 3/150 2026-03-13 10:19 by peike
[考研] 283求调剂,材料、化工皆可 +8 苏打水7777 2026-03-11 10/500 2026-03-13 09:06 by Linda Hu
[考研] 调剂 +5 呵唔哦豁 2026-03-10 5/250 2026-03-10 22:00 by 28375m
[考研] 0703化学调剂 +3 三dd. 2026-03-10 3/150 2026-03-10 15:45 by peike
[硕博家园] 木虫好像不热闹了,是不是? +4 偏振片 2026-03-10 4/200 2026-03-10 09:51 by longwave
[考研] 数二英二309分请求调剂 +3 dtdxzxx 2026-03-09 4/200 2026-03-09 19:56 by yuningshan
信息提示
请填处理意见