 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 5394 | 回復(fù): 12 | ||||||||
xtdut銀蟲 (小有名氣)
|
[交流]
【求助】如何計算有機小分子的三重態(tài)能級（T1-S0）？
|
|
看了論壇里很多大牛關(guān)于類似問題的回復(fù)，但是還是不知道在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è)地蠻干了。 |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 分子模擬要術(shù) | 不懂計算 | 量子化學(xué)-第一性原理 | Gaussian |
| 量化學(xué)者 |
» 猜你喜歡
26年儲能、電池方向博士申請
已經(jīng)有1人回復(fù)
-
求調(diào)劑，本科吉大，一志愿華工339
已經(jīng)有2人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費?
已經(jīng)有99人回復(fù)
-
意大利博洛尼亞大學(xué)招聘研究助理2名（2D材料/鋰離子與鈉離子電池方向）
已經(jīng)有1人回復(fù)
-
關(guān)于碳量子點抗菌的研究
已經(jīng)有15人回復(fù)
-
哈爾濱工業(yè)大學(xué)（深圳）綠色化學(xué)團隊招收2026年秋季入學(xué)博士生
已經(jīng)有11人回復(fù)
-
南大化學(xué)調(diào)劑
已經(jīng)有11人回復(fù)
-
澳大利亞昆士蘭科技大學(xué)PhD博士獎學(xué)金
已經(jīng)有0人回復(fù)
-
0703調(diào)劑
已經(jīng)有2人回復(fù)
-
1
已經(jīng)有14人回復(fù)
-
0703 調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)商家推薦: (我也要在這里推廣)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 循環(huán)伏安法測有機聚合物 得到聚機物能級 的相關(guān)問題
已經(jīng)有8人回復(fù)
- 雙原子或三原子小分子的激發(fā)態(tài)能級和結(jié)構(gòu)的database
已經(jīng)有3人回復(fù)
- 什么樣的有機染料在受到光激發(fā)后能夠轉(zhuǎn)變成長壽命的三重激發(fā)態(tài)？
已經(jīng)有8人回復(fù)
- 怎么計算或者在哪可以查到有機染料的HOMO能級
已經(jīng)有6人回復(fù)
- 怎么計算或者在哪可以查到有機染料的HOMO能級
已經(jīng)有3人回復(fù)
- 靜態(tài)計算還是態(tài)密度計算得到的femi能級更準(zhǔn)確？？
已經(jīng)有9人回復(fù)
- 【求助】ms中也搜尋過渡態(tài)，那ms能用來搜尋有機小分子催化的過渡態(tài)嗎？
已經(jīng)有5人回復(fù)
- 【求助】VASP怎樣精確計算出費米能級處的態(tài)密度值？
已經(jīng)有22人回復(fù)
- 【求助】問如何選擇基組精確計算有機小分子結(jié)構(gòu)和能量的問題
已經(jīng)有5人回復(fù)
- 【求助】【求助】TNT 的HOMO，LUMO能級
已經(jīng)有5人回復(fù)
- 【求助】電子態(tài)密度計算dos圖中feimi能級的確定
已經(jīng)有7人回復(fù)
- 【求助】環(huán)糊精對有機小分子包合的空間模擬怎么做
已經(jīng)有12人回復(fù)
2樓2009-11-06 11:52:17
wuy069
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 20 (小學(xué)生)
- 金幣: 7921.9
- 散金: 15
- 紅花: 10
- 帖子: 506
- 在線: 549.1小時
- 蟲號: 760723
- 注冊: 2009-04-30
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
★ ★ ★ ★ ★
xtdut(金幣+5,VIP+0):非常感謝您的回復(fù)！ 11-6 13:11
xtdut(金幣+5,VIP+0):非常感謝您的回復(fù)！ 11-6 13:11
|
首先你的3個步驟都是正確的。 第二步計算td，gaussian默認(rèn)三個態(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)型，電荷和自旋多重度當(dāng)然需要激發(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)型相當(dāng)于垂直時的基態(tài)構(gòu)型，這時再做td，從輸出文件中找到Excited State 1: Triplet-A 2.6936 eV 460.28 nm f=0.0000 |
3樓2009-11-06 12:06:19
xtdut
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 936.1
- 帖子: 97
- 在線: 212.6小時
- 蟲號: 379496
- 注冊: 2007-05-22
- 性別: GG
- 專業(yè): 材料光電子
謝謝回復(fù)，金幣已送出！
|
太感謝您了，還是有點不太明白，麻煩您耐心給指導(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 不太明白此時該計算單重態(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即可呢（當(dāng)然nstates的值也應(yīng)該相應(yīng)地改大點）？ |
4樓2009-11-06 13:10:44
wuy069
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 20 (小學(xué)生)
- 金幣: 7921.9
- 散金: 15
- 紅花: 10
- 帖子: 506
- 在線: 549.1小時
- 蟲號: 760723
- 注冊: 2009-04-30
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
★ ★ ★ ★ ★
xtdut(金幣+5,VIP+0):非常感謝！另外關(guān)于停機后接著繼續(xù)計算，你有沒有什么好的辦法？ 11-6 14:21
xtdut(金幣+5,VIP+0):非常感謝！另外關(guān)于停機后接著繼續(xù)計算，你有沒有什么好的辦法？ 11-6 14:21
|
用下面這一種 #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 ] |
5樓2009-11-06 14:08:01
wuy069
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 20 (小學(xué)生)
- 金幣: 7921.9
- 散金: 15
- 紅花: 10
- 帖子: 506
- 在線: 549.1小時
- 蟲號: 760723
- 注冊: 2009-04-30
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
★ ★ ★ ★ ★
xtdut(金幣+5,VIP+0):非常感謝您的回復(fù)！您也說了優(yōu)化可以支持接著算，那能量也可以嗎？也可以用restart嗎？希望能得到您再次耐心地指點。 11-7 08:17
xtdut(金幣+5,VIP+0):非常感謝您的回復(fù)！您也說了優(yōu)化可以支持接著算，那能量也可以嗎？也可以用restart嗎？希望能得到您再次耐心地指點。 11-7 08:17
|
優(yōu)化任務(wù)接著繼續(xù)計算： % chk=Yourcheckfilename # B3LYP/6-31G(d) opt=restart Title 0 1 或者 % chk=Yourcheckfilename # B3LYP/6-31G(d) opt geom=check guess=read Title 0 1 如果是頻率計算，g09前的版本不支持接著算，g09計算時需要chk文件，命令行要寫savenm |
6樓2009-11-06 17:49:37
cjhztt
木蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 3421.6
- 散金: 11
- 紅花: 1
- 帖子: 84
- 在線: 43.3小時
- 蟲號: 417201
- 注冊: 2007-07-01
- 專業(yè): 無機化學(xué)
7樓2009-11-08 20:13:07
jiewei
專家顧問 (著名寫手)
小才
-
專家經(jīng)驗: +52 - 應(yīng)助: 60 (初中生)
- 金幣: 7815.2
- 散金: 1399
- 紅花: 32
- 沙發(fā): 3
- 帖子: 1778
- 在線: 581.4小時
- 蟲號: 702271
- 注冊: 2009-02-16
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
- 管轄: 量子化學(xué)
8樓2009-11-08 23:16:29
wuy069
鐵桿木蟲 (正式寫手)
- 應(yīng)助: 20 (小學(xué)生)
- 金幣: 7921.9
- 散金: 15
- 紅花: 10
- 帖子: 506
- 在線: 549.1小時
- 蟲號: 760723
- 注冊: 2009-04-30
- 性別: GG
- 專業(yè): 理論和計算化學(xué)
9樓2009-11-09 09:38:51
minisun
金蟲 (正式寫手)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1839
- 散金: 6
- 帖子: 517
- 在線: 95.6小時
- 蟲號: 526725
- 注冊: 2008-03-17
- 專業(yè): 理論和計算化學(xué)
|
"要算發(fā)射，基于F-C原理，這時用優(yōu)化好的激發(fā)態(tài)的構(gòu)型，但是電荷和自旋多重度寫0 1，那這時的構(gòu)型相當(dāng)于垂直時的基態(tài)構(gòu)型，這時再做td，從輸出文件中找到Excited State 1: Triplet-A 2.6936 eV 460.28 nm f=0.0000 本文來自: 小木蟲論壇 http://www.gaoyang168.com/bbs/viewthread.php?tid=1648347 請教3樓的：用激發(fā)態(tài)的構(gòu)型為基態(tài)計算發(fā)射光譜沒錯，這時多重度是3啊，為什么還寫1呢？請不吝賜教！ |
10樓2010-12-27 12:31:14
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 xtdut 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 材料與化工求調(diào)劑 +4 | 為學(xué)666 2026-03-16 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 211本，11408一志愿中科院277分，曾在中科院自動化所實習(xí) +6 | Losir 2026-03-12 | 7/350 |
|
|
|
[考研] 085601求調(diào)劑 +3 | Du.11 2026-03-16 | 3/150 |
|
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑，一志愿985 +3 | 張vvvv 2026-03-15 | 5/250 |
|
|
|
[考研] 化學(xué)調(diào)劑0703 +8 | 啊我我的 2026-03-11 | 8/400 |
|
|
|
[基金申請]
今年的國基金是打分制嗎？
 50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
|
[考研] 中科院材料273求調(diào)劑 +4 | yzydy 2026-03-15 | 4/200 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 5/250 |
|
|
|
[考研] 招收0805（材料）調(diào)劑 +3 | 18595523086 2026-03-13 | 3/150 |
|
|
|
[考研] 0703，333分求調(diào)劑 一志愿鄭州大學(xué)-物理化學(xué) +3 | 李魔女斗篷 2026-03-11 | 3/150 |
|
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
|
[考研] 求材料調(diào)劑 +5 | 隔壁陳先生 2026-03-12 | 5/250 |
|
|
|
[考研] ［0860］321分求調(diào)劑，ab區(qū)皆可 +4 | 寶貴熱 2026-03-13 | 4/200 |
|
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
|
[考研] 工科，求調(diào)劑 +3 | 我887 2026-03-11 | 3/150 |
|
|
|
[考研] 314求調(diào)劑 +7 | 無懈可擊的巨人 2026-03-12 | 7/350 |
|
|
|
[考研] 求調(diào)劑 +3 | 程雨杭 2026-03-12 | 3/150 |
|
|
|
[考研] 0817化學(xué)工程與技術(shù)考研312分調(diào)劑 +3 | T123 tt 2026-03-12 | 3/150 |
|
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|