 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1532 | 回復(fù): 5 | |||||
beefly專家顧問 (職業(yè)作家)
地溝油冶煉專家
|
[交流]
GAMESS的gMCQDPT做RASPT2計(jì)算可能存在BUG 已有4人參與
|
|
對于大活性空間的RASPT2型計(jì)算，GAMESS的gMCQDPT模塊（$GMCPT group & MRPT=GMCPT in $MRMP group）可能存在BUG。 最早懷疑程序有問題，是三年前算CuCl2的激發(fā)能�；钚钥臻g是21電子占據(jù)17軌道（Cu:3d4s4d + Cl:3p + Cl:3p）。實(shí)驗(yàn)值，以及Molpro的RASPT2計(jì)算結(jié)果，見 J. Chem. Phys. 130, 154313, 2009，表3。用同樣的活性空間，GAMESS的RASSCF（$GMCPT group & CISTEP=GMCCI in $MCSCF group）可以得到與Molpro非常接近的參考態(tài)絕對能量（未在表3列出），但是加上相關(guān)能的RASPT2（gMCQDPT）結(jié)果非常差。 最近發(fā)現(xiàn)的例子是FeO-。最新一期的JCTC發(fā)表了RASPT2級別的計(jì)算結(jié)果（J. Chem. Theory Comput., 7, 3961, 2011，表9）。其實(shí)早在一年前，我用Molpro得到了類似的RASPT2結(jié)果，但是和GAMESS的gMCQDPT相差太遠(yuǎn)。 程序： GAMESS 12.JAN.2009(R3)； Molpro 2009.1作為參考 二者都測試了主頁提供的二進(jìn)制程序，和自己從源碼編譯的程序，結(jié)果沒有差別 計(jì)算方法： 活性空間：13電子占據(jù)17軌道（Fe:3d4s4d + O:2p3p）。一般來說，O:3p不是必須的，但是對FeO-來說，F(xiàn)e:4d和O:3p的能量非常接近以至于無法分開，只好都放入活性空間。在J. Phys. Chem. A, 113, 10588 2009這篇文章中，就是因?yàn)闆]有加上O:3p，導(dǎo)致Fe:4d的某些分量被O:3p替代，因此對FeO-的強(qiáng)關(guān)聯(lián)描述不夠準(zhǔn)確，結(jié)果和實(shí)驗(yàn)不符。 組態(tài)限制：Fe:4d + O:3p只允許0-2個(gè)電子占據(jù) 基組：DKH2標(biāo)量相對論ANO-RCC（收縮到VQZ） 排除入侵態(tài)： Molpro的RASPT2：ipea=0.25,shift=0.2 GAMESS的gMCQDPT：EDSHFT=0.03 結(jié)果： 這里只列出鍵長＝1.64 Ang.的單點(diǎn)能 RASSCF MOLPRO GAMESS 6Sigma+ -1346.4359457 -1346.4359467 Hartree 4Delta -1346.4300975 -1346.4300976 Hartree RASSCF結(jié)果完全一致。 RASPT2 MOLPRO GAMESS 6Sigma+ -1347.1556336 -1347.1554413 Hartree 4Delta -1347.1547722 -1347.1682683 Hartree Tv 0.023 -0.349 eV MOLPRO遇到入侵態(tài)會自動停止計(jì)算；GAMESS的參考態(tài)權(quán)重都在89%-90%之間。可見以上結(jié)果不存在入侵態(tài)問題。 6Sigma+態(tài)的總能量比較接近，但是4Delta態(tài)總能量相差較大。gMCQDPT模塊嚴(yán)重高估了4Delta態(tài)的相關(guān)能，導(dǎo)致垂直激發(fā)能從0.023eV一下降到-0.349eV�？晒﹨⒖嫉闹禐槲墨I(xiàn)JCTC, 7, 3961, 2011中給出的RASPT2激發(fā)能0.02 eV（用MOLCAS7開發(fā)版做的，鍵長略有不同），和MOLPRO結(jié)果一致。 不同的RASPT2方法得到不同結(jié)果是正常的，但是gMCQDPT誤差太大，連定性結(jié)果都得不到。注意GAMESS的RASSCF計(jì)算也使用gMCQDPT模塊，可見用這個(gè)模塊做RASSCF計(jì)算還是可靠的，但是在BUG排除之前，不建議用于RASPT2級別的計(jì)算。 [ Last edited by beefly on 2011-12-14 at 20:59 ] |
回復(fù)此樓» 本帖已獲得的紅花（最新10朵）
calfwoo» 猜你喜歡
化學(xué)化工學(xué)院 招收化工、化學(xué)、材料等相關(guān)方向研究生（學(xué)碩、專碩都有調(diào)劑名額）
已經(jīng)有1人回復(fù)
-
河海大學(xué) （211.雙一流高校）鈣鈦礦光電實(shí)驗(yàn)室丁勇團(tuán)隊(duì)歡迎有志青年！
已經(jīng)有3人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費(fèi)?
已經(jīng)有251人回復(fù)
-
一志愿南大理論與計(jì)算化學(xué)方向（070300）342分，求調(diào)劑
已經(jīng)有1人回復(fù)
-
歡迎化學(xué)、化工、物理、材料、或計(jì)算機(jī)背景的考研咨詢
已經(jīng)有0人回復(fù)
-
黑龍江省原子能研究院化學(xué)課題組招收2026級化學(xué)專業(yè)或化學(xué)相關(guān)專業(yè)碩士生
已經(jīng)有0人回復(fù)
-
高效微納米儲能材料與產(chǎn)業(yè)化推進(jìn)實(shí)驗(yàn)室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有0人回復(fù)
-
溫州大學(xué)新能源材料微觀溯源課題組招收2026級博士研究生（材料、化學(xué)方向）
已經(jīng)有6人回復(fù)
-
材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流
已經(jīng)有5人回復(fù)
-
26申博
已經(jīng)有2人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對您同樣有幫助:
- 2-氨基芐醇中的羥基做成氯或者溴，有什么條件
已經(jīng)有16人回復(fù)
- 各位大俠，我可以優(yōu)化結(jié)構(gòu)，但只要一計(jì)算頻率就出現(xiàn)下面的問題？
已經(jīng)有4人回復(fù)
- 計(jì)算 Gaussian 03已經(jīng) 弄壞兩臺電腦了~~~~~
已經(jīng)有46人回復(fù)
- gaussian計(jì)算
已經(jīng)有4人回復(fù)
- 求助一個(gè)計(jì)算標(biāo)準(zhǔn)偏差的方法
已經(jīng)有12人回復(fù)
- 在Gaussian 09 計(jì)算中怎么做核凍結(jié)處理？
已經(jīng)有6人回復(fù)
- 重獎(jiǎng)?wù)骷瘒鴥?nèi)國外做理論計(jì)算課題組的信息
已經(jīng)有44人回復(fù)
calfwoo
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1812.1
- 紅花: 2
- 帖子: 171
- 在線: 276.9小時(shí)
- 蟲號: 97355
- 注冊: 2005-11-08
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
2樓2011-12-14 16:22:10
Jasminer
鐵桿木蟲 (著名寫手)
- 應(yīng)助: 52 (初中生)
- 金幣: 7407.8
- 散金: 1845
- 紅花: 27
- 帖子: 1002
- 在線: 728.8小時(shí)
- 蟲號: 620156
- 注冊: 2008-10-08
- 性別: GG
- 專業(yè): 應(yīng)用有機(jī)化學(xué)
3樓2011-12-14 20:46:27
calfwoo
金蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1812.1
- 紅花: 2
- 帖子: 171
- 在線: 276.9小時(shí)
- 蟲號: 97355
- 注冊: 2005-11-08
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
4樓2011-12-15 14:38:44
dongwen176
木蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1831.3
- 紅花: 1
- 帖子: 179
- 在線: 195.6小時(shí)
- 蟲號: 91819
- 注冊: 2005-09-07
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
5樓2011-12-17 10:13:01
ZCY2009
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 200.1
- 帖子: 46
- 在線: 26.3小時(shí)
- 蟲號: 885454
- 注冊: 2009-10-27
- 專業(yè): 理論和計(jì)算化學(xué)
6樓2018-01-23 11:39:09
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 beefly 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料，紡織，生物（0856、0710），化學(xué)招生啦 +3 | Eember. 2026-03-17 | 7/350 |
|
|---|---|---|---|---|
|
[考研] 344求調(diào)劑 +4 | knight344 2026-03-16 | 4/200 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 332求調(diào)劑 +6 | Zz版 2026-03-13 | 6/300 |
|
|
|
[考研] 求調(diào)劑，總分315，考的生物醫(yī)藥，一志愿湖南師范大學(xué)。調(diào)劑到任何專業(yè)都可以 +4 | 小丁想進(jìn)步 2026-03-11 | 5/250 |
|
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學(xué)123 2026-03-15 | 8/400 |
|
|
|
[考研] 0854控制工程 359求調(diào)劑 可跨專業(yè) +3 | 626776879 2026-03-14 | 9/450 |
|
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語 2026-03-15 | 5/250 |
|
|
|
[考研] 0703一志愿211 285分求調(diào)劑 +5 | ly3471z 2026-03-13 | 5/250 |
|
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
|
[考研] 0703化學(xué)調(diào)劑，求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
|
[考研] 求老師收留調(diào)劑 +4 | jiang姜66 2026-03-14 | 5/250 |
|
|
|
[考研] 265求調(diào)劑 +4 | 威化餅07 2026-03-12 | 4/200 |
|
|
|
[考研] 【0703化學(xué)調(diào)劑】-一志愿華中師范大學(xué)-六級475 +5 | Becho359 2026-03-11 | 5/250 |
|
|
|
[基金申請]
有必要更換申報(bào)口嗎
 20+3
|
fannyamoy 2026-03-11 | 3/150 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +6 | 步川酷紫123 2026-03-13 | 6/300 |
|
|
|
[考研] 一志愿西南交大，材料專碩317求調(diào)劑 +5 | lx8568 2026-03-11 | 5/250 |
|
|
|
[考研] 求調(diào)劑 +5 | 一定有學(xué)上- 2026-03-12 | 5/250 |
|
|
|
[考研] 土木第一志愿276求調(diào)劑，科研和技能十分豐富，求新興方向的導(dǎo)師收留 +3 | 土木小天才 2026-03-12 | 3/150 |
|
|
|
[考博] 2026年博士申請 +3 | QwQwQW10 2026-03-11 | 3/150 |
|