| 5 | 1/1 | 返回列表 |
| 查看: 3586 | 回復: 7 | |||||
| 當前只顯示滿足指定條件的回帖,點擊這里查看本話題的所有回帖 | |||||
thelight金蟲 (小有名氣)
|
[求助]
求熒光計算的完整例子
|
||||
|
用高斯計算熒光,已經(jīng)完全重復了高斯09 中乙醛的熒光例子(共7步),仍然不是很明白。感覺對熒光的計算有很多人在做,但在網(wǎng)上沒有找到完整的、包括數(shù)據(jù)處理的例子。在此提出一些不明白的地方,請高手指點,最好能給一個完整的例子。 1、完成了乙醛例子的7 步后,熒光的峰位置、強度、量子產(chǎn)率是否都能得到,或者說,通過高斯能夠得到哪些信息?從結果文件中如何找到這些信息? 2、實驗測到的熒光強度與激發(fā)波長有關,計算結果是否體現(xiàn)這一點? 3、如果不考慮溶劑,只考慮真空中的分子,是否可以省略某幾步? 4、第三步的第一步讀結構所使用的chk 文件是否是第二步的結果chk文件?第三步的第一步做的是保存基態(tài)的非平衡溶劑信息,是不是應該使用第一步的結果chk文件? 5、第三步的第二步是使用 State-specific 特定態(tài)方法來垂直激發(fā),既然選擇了 Root=1 這個特定的態(tài),關心的就是這個激發(fā)態(tài),為什么還要計算6個激發(fā)態(tài) NStates=6 ? 6、第四步是優(yōu)化激發(fā)態(tài),使用了 NStates=6 Root=1,這究竟是對第一激發(fā)態(tài)進行了優(yōu)化,還是對六個激發(fā)態(tài)進行了優(yōu)化? 7、第四步,例子中用 Geom=Modify 對基態(tài)的結構進行了改變。是否可以不做這個改變,直接用 Geom=Check ? 8、第三、四、五、六步使用了 Root=1,對于自己的體系,如何決定該使用 Root=? 9、第六步,從激發(fā)態(tài)用TD 計算發(fā)射光譜時,我的理解,是計算 Root=1 的激發(fā)態(tài)到基態(tài),為什么還要使用NStates=6 ? [ Last edited by thelight on 2013-9-5 at 01:47 ] |
經(jīng)驗 | 熒光計算 | 量子化學基礎問題 |
金蟲 (小有名氣)
|
非常感謝!明白了一些,還有些問題仍需解惑。 1、做實驗的人最關心的是,為什么會有熒光以及熒光強度與結構中的什么因素有關,如果通過計算只能得到“熒光峰的位置”,豈不是對實驗沒有什么指導、參考、幫助了嗎,畢竟通過實驗已經(jīng)得到了熒光峰的位置。 2、高斯參考手冊列出了一些 SCRF 選項可以接受的溶劑,如果使用混合溶劑,比如說99:1 的乙腈:水,是不是就沒辦法了? 3、如果不考慮溶劑,對于真空分子,我的理解,應該需要 一、二、四、七,其中第一步是基態(tài)優(yōu)化,第二步是從基態(tài)到激發(fā)態(tài),第四步是激發(fā)態(tài)優(yōu)化,若無第七步,如何得到發(fā)射譜信息? 4、第四步中的Geom=Modify,如果我計算的分子并沒有對稱性,是否可以不做這個改變,直接用 Geom=Check |
木蟲 (著名寫手)
瞌睡木蟲
|
最近也接觸了一些熒光計算,試著回答一下你的問題,由于激發(fā)態(tài)計算很復雜,我也還在摸索中,因此我不保證回答的都對。 1、用Gaussian的state-specific算法算完之后,只能得到熒光峰的位置,而強度、量子產(chǎn)率等信息的計算要考慮的因素太多,目前只考慮state-specific算法還無法得到。熒光峰的位置是手工計算得到的,第四步得到的優(yōu)化結構的能量和第七部計算得到的使用激發(fā)態(tài)的優(yōu)化結構及非平衡溶劑計算得到的基態(tài)能量的差,通過轉(zhuǎn)換得到波長。 2、計算結果不能體現(xiàn)這一點; 3、不考慮溶劑時,一般只需要第一、第二和第四步,如果需要確定是否找到了激發(fā)態(tài)勢能面的平衡結構需要第五步,由于Gaussian軟件現(xiàn)在使用數(shù)值法運行激發(fā)態(tài)的頻率計算,對稍大一點的分子,其計算量太大,因此點群為C1的分子就不算了;其余步驟都可以省略; 4、第三步需要的chk文件可以用第一步的,也可以用第二步的;只要Route里面定義的不是激發(fā)態(tài)的單點計算,程序默認都是使用平衡溶劑; 5、激發(fā)態(tài)計算很復雜,有時,不通過計算很難確定哪些態(tài)的能量低,哪些高,Nstates比root多出來的數(shù)值所對應的多算的態(tài)是為了保證能找到正確的對應了root數(shù)值的態(tài); 6、激發(fā)態(tài)優(yōu)化只對root數(shù)值對應的態(tài)進行優(yōu)化; 7、不可以,乙醛這個分子在基態(tài)有對稱性,而根據(jù)原理,它在激發(fā)態(tài)的軌道占據(jù)情況需要打破這個對稱性。如果直接讀chk,則在激發(fā)態(tài)優(yōu)化中還是使用并將繼續(xù)保持基態(tài)的對稱結構,這樣激發(fā)態(tài)的計算結果肯定不準; 8、這個比較復雜,爭議也比較大,我和別人討論時,有人說應該是去計算第一激發(fā)態(tài),因為從能量上來說第一激發(fā)態(tài)所需要的激發(fā)能最小,最容易;也有人說應該根據(jù)計算出來的振子強度(f)來看,從上往下尋找第一個陣子強度比較大的激發(fā)態(tài)作為root,也就是說要確保這個激發(fā)態(tài)有吸收/發(fā)射,這樣才能和實驗結果進行對比; 9、解釋同第五條。 |

木蟲 (著名寫手)
瞌睡木蟲

木蟲 (著名寫手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600 材料與化工 329分求調(diào)劑 +14 | Mr. Z 2026-03-25 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 291求調(diào)劑 +5 | Y-cap 2026-03-29 | 6/300 |
|
|
[考研] 352分-085602-一志愿985 +5 | 海納百川Ly 2026-03-29 | 5/250 |
|
|
[考研] 330分求調(diào)劑 +5 | qzenlc 2026-03-29 | 5/250 |
|
|
[考研] 081200-11408-276學碩求調(diào)劑 +6 | 崔wj 2026-03-26 | 6/300 |
|
|
[考研] 083000學碩274求調(diào)劑 +8 | Li李魚 2026-03-26 | 8/400 |
|
|
[考研] 071000生物學求調(diào)劑,初試成績343 +7 | 小小甜面團 2026-03-25 | 7/350 |
|
|
[考研] 311(085601)求調(diào)劑 +4 | liziyeyeye 2026-03-28 | 4/200 |
|
|
[考研] 312,生物學求調(diào)劑 +3 | 小譯同學abc 2026-03-28 | 3/150 |
|
|
[考研] 299求調(diào)劑 +7 | 嗯嗯嗯嗯2 2026-03-27 | 7/350 |
|
|
[考研] 283求調(diào)劑 +7 | A child 2026-03-28 | 7/350 |
|
|
[考研] 272求調(diào)劑 +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
[考研] 材料292調(diào)劑 +12 | 橘頌思美人 2026-03-23 | 12/600 |
|
|
[考研] 0856調(diào)劑 +5 | 求求讓我有書讀?/a> 2026-03-26 | 6/300 |
|
|
[考研] 321求調(diào)劑 +6 | Ymlll 2026-03-24 | 6/300 |
|
|
[考研] 網(wǎng)絡空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
[考研] 各位老師您好:本人初試372分 +5 | jj涌77 2026-03-25 | 6/300 |
|
|
[考研] 285求調(diào)劑 +3 | AZMK 2026-03-24 | 3/150 |
|
|
[考研] 282求調(diào)劑 +3 | wcq131415 2026-03-24 | 3/150 |
|
|
[考研] 086003食品工程求調(diào)劑 +6 | 淼淼111 2026-03-24 | 6/300 |
|