 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1916 | 回復(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 ] |
回復(fù)此樓» 猜你喜歡
深圳理工大學(xué)招聘AI4Science與分子模擬方向科研人員
已經(jīng)有0人回復(fù)
-
全獎博士 英國利物浦大學(xué) × 臺灣清華大學(xué) 聯(lián)合培養(yǎng)
已經(jīng)有0人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費?
已經(jīng)有58人回復(fù)
-
尋博導(dǎo)
已經(jīng)有13人回復(fù)
-
基金正文30頁指的是報告正文還是整個申請書
已經(jīng)有9人回復(fù)
-
哈爾濱工業(yè)大學(xué)（深圳）綠色化學(xué)團隊招收2026年秋季入學(xué)博士生
已經(jīng)有11人回復(fù)
-
澳大利亞昆士蘭科技大學(xué)PhD博士獎學(xué)金
已經(jīng)有0人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- Gaussian MS計算中四氧化三鐵模型的建立
已經(jīng)有8人回復(fù)
- Gaussian 熱化學(xué)計算 問題請教
已經(jīng)有8人回復(fù)
- Gaussian計算的苯紫外光譜和教材上的差別太大，求助
已經(jīng)有10人回復(fù)
- gaussian可以計算多個分子的紅外光譜嗎
已經(jīng)有5人回復(fù)
- MP2 計算，結(jié)果不一致
已經(jīng)有6人回復(fù)
- 關(guān)于Gaussian 能量的計算
已經(jīng)有7人回復(fù)
- Gaussian 計算錯誤
已經(jīng)有4人回復(fù)
- Gaussian IRC計算是否正常結(jié)束？
已經(jīng)有6人回復(fù)
- gaussian未進行優(yōu)化計算后，分子構(gòu)型發(fā)生變化
已經(jīng)有3人回復(fù)
- 用Gaussian03或09，MP2(full)/aug-cc-pvdz如何優(yōu)化HI？
已經(jīng)有8人回復(fù)
- 在Gaussian 09 計算中怎么做核凍結(jié)處理？
已經(jīng)有6人回復(fù)
- 【beefly個人文集】用Gaussian計算超重元素
已經(jīng)有7人回復(fù)
- 【求助】高斯MP2方法計算能量的HF能量和MP2能量的區(qū)別？？？？
已經(jīng)有4人回復(fù)
- 【求助】怎樣用Gaussian計算分子直徑？？
已經(jīng)有3人回復(fù)
- 【討論】DFT方法屬不屬于從頭算？
已經(jīng)有13人回復(fù)
lihb734
鐵桿木蟲 (職業(yè)作家)
站在計算化學(xué)入門的門檻上
- QC強帖: 2
- 應(yīng)助: 132 (高中生)
- 金幣: 7495.4
- 散金: 9378
- 紅花: 74
- 帖子: 3107
- 在線: 1177.3小時
- 蟲號: 448403
- 注冊: 2007-11-01
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
2樓2010-03-21 21:27:25
piaoma
木蟲 (小有名氣)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 2725.5
- 散金: 157
- 帖子: 162
- 在線: 590.2小時
- 蟲號: 201836
- 注冊: 2006-02-28
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
3樓2010-03-21 22:17:05
sobereva
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
- QC強帖: 17
- 應(yīng)助: 226 (大學(xué)生)
- 貴賓: 5.4
- 金幣: 2180.5
- 紅花: 217
- 帖子: 1241
- 在線: 366.8小時
- 蟲號: 281271
- 注冊: 2006-09-23
- 專業(yè): 宇宙學(xué)
★
piaoma(金幣+3):謝謝參與
piaoma(金幣+50): 2010-03-22 10:54
piaoma(金幣+3):謝謝參與
piaoma(金幣+50): 2010-03-22 10:54
|
T2應(yīng)當(dāng)是其文中公式8，是構(gòu)建校正能的中間量，這不用考慮。 列出來的三個E2分項是每類電子對兒對總校正能的貢獻，三項相加就是下面的總校正能E2=-0.2006947749D+01 �？傂Ｕ蹺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)和基組計算與原文對照一下。 |
4樓2010-03-22 04:28:30
lihb734
鐵桿木蟲 (職業(yè)作家)
站在計算化學(xué)入門的門檻上
- QC強帖: 2
- 應(yīng)助: 132 (高中生)
- 金幣: 7495.4
- 散金: 9378
- 紅花: 74
- 帖子: 3107
- 在線: 1177.3小時
- 蟲號: 448403
- 注冊: 2007-11-01
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
5樓2010-03-22 09:21:51
lihb734
鐵桿木蟲 (職業(yè)作家)
站在計算化學(xué)入門的門檻上
- QC強帖: 2
- 應(yīng)助: 132 (高中生)
- 金幣: 7495.4
- 散金: 9378
- 紅花: 74
- 帖子: 3107
- 在線: 1177.3小時
- 蟲號: 448403
- 注冊: 2007-11-01
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
★ ★ ★
piaoma(金幣+40): 2010-03-22 10:54
yjcmwgk(金幣+3): 2010-03-23 13:37
piaoma(金幣+40): 2010-03-22 10:54
yjcmwgk(金幣+3): 2010-03-23 13:37
|
根據(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 |
6樓2010-03-22 09:39:22
piaoma
木蟲 (小有名氣)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 2725.5
- 散金: 157
- 帖子: 162
- 在線: 590.2小時
- 蟲號: 201836
- 注冊: 2006-02-28
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
7樓2010-03-22 15:50:54
lihb734
鐵桿木蟲 (職業(yè)作家)
站在計算化學(xué)入門的門檻上
- QC強帖: 2
- 應(yīng)助: 132 (高中生)
- 金幣: 7495.4
- 散金: 9378
- 紅花: 74
- 帖子: 3107
- 在線: 1177.3小時
- 蟲號: 448403
- 注冊: 2007-11-01
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
8樓2010-03-22 16:01:26
sobereva
至尊木蟲 (著名寫手)
本人已永久離開小木蟲
- QC強帖: 17
- 應(yīng)助: 226 (大學(xué)生)
- 貴賓: 5.4
- 金幣: 2180.5
- 紅花: 217
- 帖子: 1241
- 在線: 366.8小時
- 蟲號: 281271
- 注冊: 2006-09-23
- 專業(yè): 宇宙學(xué)
★ ★ ★
yjcmwgk(金幣+3): 2010-03-23 13:37
yjcmwgk(金幣+3): 2010-03-23 13:37
|
文中的這句話容易造成誤解： 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）。 |
9樓2010-03-22 16:53:02
piaoma
木蟲 (小有名氣)
- 應(yīng)助: 6 (幼兒園)
- 金幣: 2725.5
- 散金: 157
- 帖子: 162
- 在線: 590.2小時
- 蟲號: 201836
- 注冊: 2006-02-28
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
10樓2010-03-22 16:59:10
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 piaoma 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0856材料專碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 280求調(diào)劑 +11 | 咕嚕曉曉 2026-03-18 | 12/600 |
|
|
|
[考研] 354求調(diào)劑 +6 | Tyoumou 2026-03-18 | 9/450 |
|
|
|
[考研] 302求調(diào)劑 +12 | 呼呼呼。。。。 2026-03-17 | 12/600 |
|
|
|
[考研] 268求調(diào)劑 +9 | 簡單點0 2026-03-17 | 9/450 |
|
|
|
[考研] 材料與化工（0856）304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
|
[考研] 330求調(diào)劑0854 +3 | assdll 2026-03-21 | 3/150 |
|
|
|
[考研] 070300化學(xué)319求調(diào)劑 +7 | 錦鯉0909 2026-03-17 | 7/350 |
|
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
|
[考研] 332求調(diào)劑 +4 | ydfyh 2026-03-17 | 4/200 |
|
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
|
[考研] 294求調(diào)劑材料與化工專碩 +15 | 陌の森林 2026-03-18 | 15/750 |
|
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
|
[考研] 中南大學(xué)化學(xué)學(xué)碩337求調(diào)劑 +3 | niko- 2026-03-19 | 6/300 |
|
|
|
[考研] 289求調(diào)劑 +6 | 懷瑾握瑜l 2026-03-20 | 6/300 |
|
|
|
[考研] 261求B區(qū)調(diào)劑，科研經(jīng)歷豐富 +3 | 牛奶很忙 2026-03-20 | 4/200 |
|
|
|
[考研] 廣西大學(xué)家禽遺傳育種課題組2026年碩士招生（接收計算機專業(yè)調(diào)劑） +3 | 123阿標(biāo) 2026-03-17 | 3/150 |
|
|
|
[考研] 085600材料與化工調(diào)劑 324分 +10 | llllkkkhh 2026-03-18 | 12/600 |
|
|
|
[考博] 26博士申請 +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|