| 查看: 1918 | 回復(fù): 12 | |||
piaoma木蟲 (小有名氣)
|
[交流]
【討論】用gaussian進行SCS-mp2計算
|
|
最近琢磨怎么用gaussian來做scs-mp2 文獻見S.Grimme J. Chem. Phys 2003,118,9095 需要用到gaussian計算結(jié)果的數(shù)據(jù),如下(為g03的輸出文件) Spin components of T(2) and E(2): alpha-alpha T2 = 0.9386991328D-01 E2= -0.2551138229D+00 alpha-beta T2 = 0.4903096895D+00 E2= -0.1496720103D+01 beta-beta T2 = 0.9386991328D-01 E2= -0.2551138229D+00 ANorm= 0.1295395506D+01 E2 = -0.2006947749D+01 EUMP2 = -0.49555972975718D+03 1)alpha-alpha、alpha-beta 和beta-beta 應(yīng)指電子的自旋情況。請問 T2 和E2指什么? 2)分子能量應(yīng)為EUMP2 那E2又是指什么呢? 請指教. [ Last edited by piaoma on 2010-3-21 at 20:32 ] |
鐵桿木蟲 (職業(yè)作家)
站在計算化學(xué)入門的門檻上

木蟲 (小有名氣)
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
T2應(yīng)當(dāng)是其文中公式8,是構(gòu)建校正能的中間量,這不用考慮。 列出來的三個E2分項是每類電子對兒對總校正能的貢獻,三項相加就是下面的總校正能E2=-0.2006947749D+01 ?傂U蹺2加到HF能量上就是EUMP2。 SCS-MP2校正能應(yīng)當(dāng)為1.2*-0.1496720103D+01 + 1/3*(2*-0.2551138229D+00),加到HF能量上就是SCS-MP2能量。建議用文中的結(jié)構(gòu)和基組計算與原文對照一下。 |
鐵桿木蟲 (職業(yè)作家)
站在計算化學(xué)入門的門檻上

鐵桿木蟲 (職業(yè)作家)
站在計算化學(xué)入門的門檻上
|
根據(jù)您的提示以及JCP, 118,9095中結(jié)論部分的這句話“One important aspect of the idea is to require that the total correlation energy after scaling (roughly) equals the unscaled MP2 energy.“ J. Chem. Phys. 129, 024103 (2008),這篇文獻中的公式5 E(MP2) = E(HF) + Ec = E(HF) + c↑↓Ec↑↓ + 0.5c↑↑Ec↑↑ + 0.5c↓↓Ec↓↓, 應(yīng)該是錯的,我覺得應(yīng)該沒有前邊的系數(shù)1/2了。c↑↓=1.2, c↓↓=0.33 請您指教。理由如下: Ec = c↑↓Ec↑↓ + 0.5c↑↑Ec↑↑ + 0.5c↓↓Ec↓↓=1.2*-0.1496720103D+01 + 0.33* -0.2551138229D+00 = -1.8806 不等于E2 = -0.2006947749D+01 , 而略去二分之一時: Ec = c↑↓Ec↑↓ + 0.5c↑↑Ec↑↑ + 0.5c↓↓Ec↓↓=1.2*-0.1496720103D+01 + 2*0.33* -0.2551138229D+00 = -1.9643 約等于E2 =-0.2006947749D+01 |

木蟲 (小有名氣)
鐵桿木蟲 (職業(yè)作家)
站在計算化學(xué)入門的門檻上

至尊木蟲 (著名寫手)
本人已永久離開小木蟲
|
文中的這句話容易造成誤解: ES and ET are given by contributions from electron pairs with aa, bb, and ab spin as... 他的順序表面上反了,但是并沒有說respectively。而且他的意思是讓讀者看下面的式子,下面的式子含義正確。所以這沒必要追究這句話。正確意思就是ET是自旋平行的貢獻,ES是反平行的貢獻。 文中那個1/2應(yīng)當(dāng)沒有實際意義。文中Ec=ES+ET,高斯輸出中Ec=E2=E2↑↓+ E2↓↓+E2↑↑,顯然E2↓↓+E2↑↑=ET,故高斯輸出的E2↓↓和E2↑↑就對應(yīng)文中式2的兩項(包括了1/2)。 |
木蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 315分,誠求調(diào)劑,材料與化工085600 +3 | 13756423260 2026-03-22 | 3/150 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
[考研] 305分求調(diào)劑(食品工程) +4 | Sxy112 2026-03-21 | 6/300 |
|
|
[考研] 260求調(diào)劑 +3 | 朱芷琳 2026-03-20 | 4/200 |
|
|
[考研] 材料學(xué)碩301分求調(diào)劑 +7 | Liyouyumairs 2026-03-21 | 7/350 |
|
|
[考研] 0805 316求調(diào)劑 +3 | 大雪深藏 2026-03-18 | 3/150 |
|
|
[考研] 278求調(diào)劑 +9 | 煙火先于春 2026-03-17 | 9/450 |
|
|
[考研] 279求調(diào)劑 +5 | 紅衣隱官 2026-03-21 | 5/250 |
|
|
[考研] 279分求調(diào)劑 一志愿211 +14 | chaojifeixia 2026-03-19 | 15/750 |
|
|
[考研] 南昌大學(xué)材料專碩311分求調(diào)劑 +6 | 77chaselx 2026-03-20 | 6/300 |
|
|
[考研] 346求調(diào)劑[0856] +4 | WayneLim327 2026-03-16 | 7/350 |
|
|
[考研] 一志愿 西北大學(xué) ,070300化學(xué)學(xué)碩,總分287,雙非一本,求調(diào)劑。 +3 | 晨昏線與星海 2026-03-18 | 3/150 |
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩,總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
[考研]
|
簡木ChuFront 2026-03-19 | 8/400 |
|
|
[考研] 08工學(xué)調(diào)劑 +5 | 用戶573181 2026-03-20 | 5/250 |
|
|
[考研] 一志愿中國海洋大學(xué),生物學(xué),301分,求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
[考研] 283求調(diào)劑 +3 | 聽風(fēng)就是雨; 2026-03-16 | 3/150 |
|