| 5 | 1/1 | 返回列表 |
| 查看: 5445 | 回復: 12 | ||||||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | ||||||||
xtdut銀蟲 (小有名氣)
|
[交流]
【求助】如何計算有機小分子的三重態(tài)能級(T1-S0)?
|
|||||||
|
看了論壇里很多大牛關(guān)于類似問題的回復,但是還是不知道在Gaussian中如何具體操作。拿最簡單的吡啶環(huán)為例來說吧: 1.先優(yōu)化基態(tài)的最低能量構(gòu)型: # 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)的最低能量構(gòu)型: # opt b3lyp/6-31+g(d) Charge = 0 Multiplicity = 3 問題是第二步計算了三個三重態(tài),此步優(yōu)化的是不是這三個三重態(tài)呢?還是應(yīng)該在第二步中只計算一個三重態(tài),然后在此步中優(yōu)化此三重態(tài)也就是我們需要的T1態(tài)呢?還有這一步中自旋多重度為3沒有錯誤吧?引申一下,若是在第二步中計算了十個單重態(tài)(也即十個吸收峰,假設(shè)它們均存在),在此步中優(yōu)化單重態(tài)的最低能量構(gòu)型如何計算,自旋多重度也是3對嗎? 4.最關(guān)鍵一步,就是計算從最低能量構(gòu)型的三重態(tài)回落到基態(tài)的過程,不知道如何計算? 我知道對各位大牛來說,這是個最簡單計算發(fā)射光譜(不論是熒光光譜還是磷光光譜)的事情,原理也就是F-C原理搞清楚了,計算應(yīng)該沒有問題,可是對我們菜鳥來說,這些具體的操作才是問題的關(guān)鍵,希望各位大牛能夠耐心指點一下,而不是一句F-C原理了事,謝謝大家了! 另外相關(guān)問題是:在計算類似問題時,我們這些非專業(yè)的人也就是用普通帶XP的計算機算算而已,分子中要是有3到7個環(huán)左右(而這在有機合成中是非常常見的)時,可能光優(yōu)化構(gòu)型就得好幾天,這中途要是碰上什么事情(停電、誤操作等)停機了,就只能從頭開始算,有沒有什么辦法能夠接著上次終止的地方繼續(xù)計算的,如何操作?我知道對各位專業(yè)人士來說這種事情無法想象,可國情如此,只好這么不專業(yè)地蠻干了。 |
分子模擬要術(shù) | 不懂計算 | 量子化學-第一性原理 | Gaussian |
量化學者 |
銀蟲 (小有名氣)
|
太感謝您了,還是有點不太明白,麻煩您耐心給指導一下: 也就是說第四步應(yīng)該這樣寫: # td=(singlets,nstates=3,root=1) b3lyp/6-31+g(d) 0 1 還是應(yīng)該這樣寫: # 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的值也應(yīng)該相應(yīng)地改大點)? |
鐵桿木蟲 (正式寫手)
|
首先你的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)的構(gòu)型,所以電荷和自旋多重度都是基態(tài)的。 第3步,優(yōu)化激發(fā)態(tài)的構(gòu)型,電荷和自旋多重度當然需要激發(fā)態(tài)的,所以電荷和自旋多重度為0 3. 要想計算有機小分子的三重態(tài)能級(T1-S0),即是算激發(fā)態(tài)的發(fā)射,但是gaussian中沒這么直接的,我們可以轉(zhuǎn)換一下思路:要算發(fā)射,基于F-C原理,這時用優(yōu)化好的激發(fā)態(tài)的構(gòu)型,但是電荷和自旋多重度寫0 1,那這時的構(gòu)型相當于垂直時的基態(tài)構(gòu)型,這時再做td,從輸出文件中找到Excited State 1: Triplet-A 2.6936 eV 460.28 nm f=0.0000 |
鐵桿木蟲 (正式寫手)
|
用下面這一種 #p td=(triplets,nstates=15,root=1,direct) b3lyp/6-31+g(d) test 0 1 因為是算三重態(tài),所以命令行寫triplets;而用的是基態(tài)的構(gòu)型,所以電荷和自旋多重度寫0 1 你第3步優(yōu)化的是第一個三重態(tài)(即T1態(tài)),后面說的對的,但要算第二或三激發(fā)態(tài)的td,你應(yīng)該優(yōu)化好了激發(fā)態(tài)再做td. 如果你研究的分子比較大,建議把nstates寫大一些,要不可能會出錯 [ Last edited by wuy069 on 2009-11-6 at 14:12 ] |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料專碩調(diào)劑 +9 | 椰椰。 2026-03-29 | 9/450 |
|
|---|---|---|---|---|
|
[考研] 279求調(diào)劑 +6 | j的立方 2026-03-29 | 6/300 |
|
|
[考研] 生物技術(shù)與工程 +7 | 1294608413 2026-03-25 | 8/400 |
|
|
[考研] 085601材料工程找調(diào)劑 +16 | oatmealR 2026-03-29 | 17/850 |
|
|
[考研] 282求調(diào)劑 +4 | wcq131415 2026-03-24 | 4/200 |
|
|
[考研] 343求調(diào)劑 +6 | 愛羈絆 2026-03-29 | 6/300 |
|
|
[考研] 332求92調(diào)劑 +8 | 蕉蕉123 2026-03-28 | 8/400 |
|
|
[考研] 學碩274求調(diào)劑 +9 | Li李魚 2026-03-26 | 9/450 |
|
|
[考研] 085701求調(diào)劑初試286分 +4 | secret0328 2026-03-28 | 4/200 |
|
|
[考研] 本科新能源科學與工程,一志愿華理能動285求調(diào)劑 +7 | AZMK 2026-03-28 | 11/550 |
|
|
[考研] 083000學碩274求調(diào)劑 +8 | Li李魚 2026-03-26 | 8/400 |
|
|
[考研] 347求調(diào)劑 +3 | 山頂見α 2026-03-25 | 3/150 |
|
|
[考研] 292求調(diào)劑 +14 | 鵝鵝鵝額額額額?/a> 2026-03-25 | 15/750 |
|
|
[考研] 085701環(huán)境工程求調(diào)劑 +9 | 多久上課 2026-03-27 | 9/450 |
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +4 | 崔wj 2026-03-26 | 4/200 |
|
|
[考研] 321求調(diào)劑 +6 | wasdssaa 2026-03-26 | 6/300 |
|
|
[考研] 中國科學院深圳先進技術(shù)研究院-光纖傳感課題組招生-中國科學院大學、深圳理工大學聯(lián)培 +5 | YangTyu1 2026-03-26 | 5/250 |
|
|
[考研]
|
WWW西西弗斯 2026-03-24 | 8/400 |
|
|
[考研] 07化學303求調(diào)劑 +5 | 睿08 2026-03-25 | 5/250 |
|