| 5 | 1/1 | 返回列表 |
| 查看: 5442 | 回復: 12 | ||||||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||||||
xtdut銀蟲 (小有名氣)
|
[交流]
【求助】如何計算有機小分子的三重態(tài)能級(T1-S0)?
|
|||||||
|
看了論壇里很多大牛關于類似問題的回復,但是還是不知道在Gaussian中如何具體操作。拿最簡單的吡啶環(huán)為例來說吧: 1.先優(yōu)化基態(tài)的最低能量構型: # opt b3lyp/6-31+g(d) Charge = 0 Multiplicity = 1 2.接著計算三重態(tài)的所謂垂直激發(fā)能: # td=(triplets,nstates=3,root=1) b3lyp/6-31+g(d) Charge = 0 Multiplicity = 1 3.然后優(yōu)化此三重態(tài)的最低能量構型: # opt b3lyp/6-31+g(d) Charge = 0 Multiplicity = 3 問題是第二步計算了三個三重態(tài),此步優(yōu)化的是不是這三個三重態(tài)呢?還是應該在第二步中只計算一個三重態(tài),然后在此步中優(yōu)化此三重態(tài)也就是我們需要的T1態(tài)呢?還有這一步中自旋多重度為3沒有錯誤吧?引申一下,若是在第二步中計算了十個單重態(tài)(也即十個吸收峰,假設它們均存在),在此步中優(yōu)化單重態(tài)的最低能量構型如何計算,自旋多重度也是3對嗎? 4.最關鍵一步,就是計算從最低能量構型的三重態(tài)回落到基態(tài)的過程,不知道如何計算? 我知道對各位大牛來說,這是個最簡單計算發(fā)射光譜(不論是熒光光譜還是磷光光譜)的事情,原理也就是F-C原理搞清楚了,計算應該沒有問題,可是對我們菜鳥來說,這些具體的操作才是問題的關鍵,希望各位大牛能夠耐心指點一下,而不是一句F-C原理了事,謝謝大家了! 另外相關問題是:在計算類似問題時,我們這些非專業(yè)的人也就是用普通帶XP的計算機算算而已,分子中要是有3到7個環(huán)左右(而這在有機合成中是非常常見的)時,可能光優(yōu)化構型就得好幾天,這中途要是碰上什么事情(停電、誤操作等)停機了,就只能從頭開始算,有沒有什么辦法能夠接著上次終止的地方繼續(xù)計算的,如何操作?我知道對各位專業(yè)人士來說這種事情無法想象,可國情如此,只好這么不專業(yè)地蠻干了。 |
分子模擬要術 | 不懂計算 | 量子化學-第一性原理 | Gaussian |
量化學者 |
鐵桿木蟲 (正式寫手)
|
首先你的3個步驟都是正確的。 第二步計算td,gaussian默認三個態(tài),如果是對稱性高的分子建議用 Largest Abelian subgroup 后的Nop 個數(shù)。為什么計算一個態(tài)需要用3個態(tài)或更多的態(tài)呢,那是因為求解這個態(tài),需要用其它的態(tài)來對它就解。第2步,自旋多重度還是寫1,因為垂直激發(fā)是基于基態(tài)的構型,所以電荷和自旋多重度都是基態(tài)的。 第3步,優(yōu)化激發(fā)態(tài)的構型,電荷和自旋多重度當然需要激發(fā)態(tài)的,所以電荷和自旋多重度為0 3. 要想計算有機小分子的三重態(tài)能級(T1-S0),即是算激發(fā)態(tài)的發(fā)射,但是gaussian中沒這么直接的,我們可以轉(zhuǎn)換一下思路:要算發(fā)射,基于F-C原理,這時用優(yōu)化好的激發(fā)態(tài)的構型,但是電荷和自旋多重度寫0 1,那這時的構型相當于垂直時的基態(tài)構型,這時再做td,從輸出文件中找到Excited State 1: Triplet-A 2.6936 eV 460.28 nm f=0.0000 |
銀蟲 (小有名氣)
|
太感謝您了,還是有點不太明白,麻煩您耐心給指導一下: 也就是說第四步應該這樣寫: # td=(singlets,nstates=3,root=1) b3lyp/6-31+g(d) 0 1 還是應該這樣寫: # td=(triplets,nstates=3,root=1) b3lyp/6-31+g(d) 0 1 不太明白此時該計算單重態(tài)的發(fā)射還是三重態(tài)的發(fā)射?另外還想請教您一下,我在第三步優(yōu)化時優(yōu)化了三個三重態(tài),還是只優(yōu)化了第一個三重態(tài)(也即T1態(tài)),如果我想知道第二或第三個三重態(tài)的能級(也即T2或者T3時),是不是只需要在第四步中將root改為等于2或者3即可呢(當然nstates的值也應該相應地改大點)? |
鐵桿木蟲 (正式寫手)
|
用下面這一種 #p td=(triplets,nstates=15,root=1,direct) b3lyp/6-31+g(d) test 0 1 因為是算三重態(tài),所以命令行寫triplets;而用的是基態(tài)的構型,所以電荷和自旋多重度寫0 1 你第3步優(yōu)化的是第一個三重態(tài)(即T1態(tài)),后面說的對的,但要算第二或三激發(fā)態(tài)的td,你應該優(yōu)化好了激發(fā)態(tài)再做td. 如果你研究的分子比較大,建議把nstates寫大一些,要不可能會出錯 [ Last edited by wuy069 on 2009-11-6 at 14:12 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 調(diào)劑310 +12 | 溫柔的晚安 2026-03-25 | 13/650 |
|
|---|---|---|---|---|
|
[考研] 08工科求調(diào)劑286 +3 | tgs_001 2026-03-28 | 3/150 |
|
|
[考研] 334分 一志愿武理 材料求調(diào)劑 +7 | 李李不服輸 2026-03-26 | 7/350 |
|
|
[考研] 305求調(diào)劑 +8 | RuiFairyrui 2026-03-28 | 8/400 |
|
|
[考研] 壓國家一區(qū)線,求導師收留,有恩必謝! +7 | 迷人的哈哈 2026-03-28 | 7/350 |
|
|
[考研] 346求調(diào)劑 一志愿070303有機化學 +3 | 蘿卜燉青菜 2026-03-28 | 3/150 |
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
[考研] 一志愿南昌大學324求調(diào)劑 +7 | hanamiko 2026-03-27 | 7/350 |
|
|
[考研] 085404求調(diào)劑,總分309,本科經(jīng)歷較為豐富 +4 | 來財aa 2026-03-25 | 4/200 |
|
|
[考研] 安徽大學專碩生物與醫(yī)藥專業(yè)(086000)324分,英語已過四六級,六級521,求調(diào)劑 +4 | 美味可樂雞翅 2026-03-26 | 4/200 |
|
|
[考研] 305求調(diào)劑 +5 | 哇盧卡庫 2026-03-26 | 5/250 |
|
|
[考研] 315調(diào)劑 +4 | 0860求調(diào)劑 2026-03-26 | 5/250 |
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
[考研] 調(diào)劑 +4 | 柚柚yoyo 2026-03-26 | 4/200 |
|
|
[考研] 機械學碩310分,數(shù)一英一,一志愿211本科雙非找調(diào)劑信息 +3 | @357 2026-03-25 | 3/150 |
|
|
[考研] 一志愿 南京郵電大學 288分 材料考研 求調(diào)劑 +3 | jl0720 2026-03-26 | 3/150 |
|
|
[考研] 網(wǎng)絡空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-23 | 4/200 |
|
|
[考研] 接收2026碩士調(diào)劑(學碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|