| 5 | 1/1 | 返回列表 |
| 查看: 5463 | 回復(fù): 12 | ||||||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||||||
xtdut銀蟲 (小有名氣)
|
[交流]
【求助】如何計(jì)算有機(jī)小分子的三重態(tài)能級(jí)(T1-S0)?
|
|||||||
|
看了論壇里很多大牛關(guān)于類似問題的回復(fù),但是還是不知道在Gaussian中如何具體操作。拿最簡(jiǎn)單的吡啶環(huán)為例來說吧: 1.先優(yōu)化基態(tài)的最低能量構(gòu)型: # opt b3lyp/6-31+g(d) Charge = 0 Multiplicity = 1 2.接著計(jì)算三重態(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 問題是第二步計(jì)算了三個(gè)三重態(tài),此步優(yōu)化的是不是這三個(gè)三重態(tài)呢?還是應(yīng)該在第二步中只計(jì)算一個(gè)三重態(tài),然后在此步中優(yōu)化此三重態(tài)也就是我們需要的T1態(tài)呢?還有這一步中自旋多重度為3沒有錯(cuò)誤吧?引申一下,若是在第二步中計(jì)算了十個(gè)單重態(tài)(也即十個(gè)吸收峰,假設(shè)它們均存在),在此步中優(yōu)化單重態(tài)的最低能量構(gòu)型如何計(jì)算,自旋多重度也是3對(duì)嗎? 4.最關(guān)鍵一步,就是計(jì)算從最低能量構(gòu)型的三重態(tài)回落到基態(tài)的過程,不知道如何計(jì)算? 我知道對(duì)各位大牛來說,這是個(gè)最簡(jiǎn)單計(jì)算發(fā)射光譜(不論是熒光光譜還是磷光光譜)的事情,原理也就是F-C原理搞清楚了,計(jì)算應(yīng)該沒有問題,可是對(duì)我們菜鳥來說,這些具體的操作才是問題的關(guān)鍵,希望各位大牛能夠耐心指點(diǎn)一下,而不是一句F-C原理了事,謝謝大家了! 另外相關(guān)問題是:在計(jì)算類似問題時(shí),我們這些非專業(yè)的人也就是用普通帶XP的計(jì)算機(jī)算算而已,分子中要是有3到7個(gè)環(huán)左右(而這在有機(jī)合成中是非常常見的)時(shí),可能光優(yōu)化構(gòu)型就得好幾天,這中途要是碰上什么事情(停電、誤操作等)停機(jī)了,就只能從頭開始算,有沒有什么辦法能夠接著上次終止的地方繼續(xù)計(jì)算的,如何操作?我知道對(duì)各位專業(yè)人士來說這種事情無法想象,可國(guó)情如此,只好這么不專業(yè)地蠻干了。 |
分子模擬要術(shù) | 不懂計(jì)算 | 量子化學(xué)-第一性原理 | Gaussian |
量化學(xué)者 |
木蟲 (小有名氣)
| 關(guān)于停電的事情,直接讀取文件的最后一個(gè)結(jié)構(gòu)計(jì)算啊。 |
鐵桿木蟲 (正式寫手)
|
首先你的3個(gè)步驟都是正確的。 第二步計(jì)算td,gaussian默認(rèn)三個(gè)態(tài),如果是對(duì)稱性高的分子建議用 Largest Abelian subgroup 后的Nop 個(gè)數(shù)。為什么計(jì)算一個(gè)態(tài)需要用3個(gè)態(tài)或更多的態(tài)呢,那是因?yàn)榍蠼膺@個(gè)態(tài),需要用其它的態(tài)來對(duì)它就解。第2步,自旋多重度還是寫1,因?yàn)榇怪奔ぐl(fā)是基于基態(tài)的構(gòu)型,所以電荷和自旋多重度都是基態(tài)的。 第3步,優(yōu)化激發(fā)態(tài)的構(gòu)型,電荷和自旋多重度當(dāng)然需要激發(fā)態(tài)的,所以電荷和自旋多重度為0 3. 要想計(jì)算有機(jī)小分子的三重態(tài)能級(jí)(T1-S0),即是算激發(fā)態(tài)的發(fā)射,但是gaussian中沒這么直接的,我們可以轉(zhuǎn)換一下思路:要算發(fā)射,基于F-C原理,這時(shí)用優(yōu)化好的激發(fā)態(tài)的構(gòu)型,但是電荷和自旋多重度寫0 1,那這時(shí)的構(gòu)型相當(dāng)于垂直時(shí)的基態(tài)構(gòu)型,這時(shí)再做td,從輸出文件中找到Excited State 1: Triplet-A 2.6936 eV 460.28 nm f=0.0000 |
銀蟲 (小有名氣)
|
太感謝您了,還是有點(diǎn)不太明白,麻煩您耐心給指導(dǎo)一下: 也就是說第四步應(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 不太明白此時(shí)該計(jì)算單重態(tài)的發(fā)射還是三重態(tài)的發(fā)射?另外還想請(qǐng)教您一下,我在第三步優(yōu)化時(shí)優(yōu)化了三個(gè)三重態(tài),還是只優(yōu)化了第一個(gè)三重態(tài)(也即T1態(tài)),如果我想知道第二或第三個(gè)三重態(tài)的能級(jí)(也即T2或者T3時(shí)),是不是只需要在第四步中將root改為等于2或者3即可呢(當(dāng)然nstates的值也應(yīng)該相應(yīng)地改大點(diǎn))? |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑!生物與醫(yī)藥專碩 +4 | 逆轉(zhuǎn)陸先生 2026-04-01 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 319求調(diào)劑 +13 | 太容易1018 2026-04-01 | 13/650 |
|
|
[考研] 材料化工340求調(diào)劑 +4 | jhx777 2026-03-30 | 4/200 |
|
|
[考研] 一志愿鄭大材料工程290求調(diào)劑 +19 | Youth_ 2026-03-30 | 19/950 |
|
|
[考研] 261求B區(qū)調(diào)劑 +5 | 明仔· 2026-04-01 | 7/350 |
|
|
[考研] 一志愿北交大材料工程總分358 +8 | cs0106 2026-04-01 | 9/450 |
|
|
[考研] 085600,材料與化工321分調(diào)劑 +8 | 大饞小子 2026-03-27 | 10/500 |
|
|
[考研] 316求調(diào)劑 +11 | 舟自梗 2026-04-01 | 12/600 |
|
|
[考研] 301求調(diào)劑 +13 | A_JiXing 2026-04-01 | 13/650 |
|
|
[考研] 食品學(xué)碩362求調(diào)劑 +3 | xuanxianxian 2026-04-01 | 3/150 |
|
|
[考研] 350求調(diào)劑 +7 | 阿佳~ 2026-03-31 | 7/350 |
|
|
[考研] 英一數(shù)二生物信息學(xué)287分,本科生物科學(xué),求調(diào)劑 +3 | 碧水xyz 2026-03-29 | 4/200 |
|
|
[考研] 283求調(diào)劑 +9 | A child 2026-03-28 | 9/450 |
|
|
[考研] 085701環(huán)境工程,267求調(diào)劑 +17 | minht 2026-03-26 | 17/850 |
|
|
[考研] 一志愿西電085401數(shù)一英一299求調(diào)劑 六級(jí)521 +4 | 愛吃大鴨梨 2026-03-31 | 4/200 |
|
|
[考研] 287求調(diào)劑 +17 | land xuxu 2026-03-26 | 17/850 |
|
|
[考研] 福建理工大學(xué)材料學(xué)院先進(jìn)合金團(tuán)隊(duì)招收考研調(diào)劑學(xué)生 +3 | 大華金商都 2026-03-30 | 4/200 |
|
|
[考研] 327求調(diào)劑 +5 | 小卡不卡. 2026-03-29 | 5/250 |
|
|
[考研] 數(shù)一英一271專碩(085401)求調(diào)劑,可跨 +7 | 前行必有光 2026-03-28 | 8/400 |
|
|
[考研] 085602 化工專碩 338分 求調(diào)劑 +12 | 路癡小琪 2026-03-27 | 12/600 |
|