| 5 | 1/1 | 返回列表 |
| 查看: 2468 | 回復(fù): 5 | |||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | |||
小范范1989木蟲 (著名寫手)
|
[交流]
高斯優(yōu)化最低三重態(tài)的問題。 已有1人參與
|
||
|
之前,我對(duì)于優(yōu)化單重態(tài)和三重態(tài)有自己的理解,如下: 1:優(yōu)化s1的結(jié)構(gòu)的時(shí)候, #p b3lyp/6-31g opt td(singlets,root=1) 這種方法應(yīng)該是對(duì)的吧,同時(shí)是不是會(huì)默認(rèn)計(jì)算三個(gè)td?那和下面的計(jì)算方法有什么區(qū)別呢? #p b3lyp/6-31g opt td(nstates=6,singlets,root=1) 這個(gè)事計(jì)算6個(gè)td?和上一個(gè),就這一點(diǎn)區(qū)別? 對(duì)于這兩個(gè)的區(qū)別,我還在測(cè)試。 2:就是優(yōu)化三重態(tài)的問題,這個(gè)是我的主要問題:如下 第一種 #P opt b3lyp/6-31g td(triplets,root=1) 0 1 坐標(biāo) 我的問題是:此處的自旋多重讀為1,但是td中要寫triplets。 同時(shí),坐標(biāo)最好為s0的基態(tài)結(jié)構(gòu)坐標(biāo)、這樣準(zhǔn)確。是吧? 第二種 優(yōu)化t1的第二種方式: #P b3lyp/6-31g opt 0 3 坐標(biāo) 這個(gè)自旋多重讀為3, 我的問題是,我測(cè)試過這兩種方式,得出來(lái)的HF能量是不同的,同時(shí),幾何結(jié)構(gòu)也不相同,我想知道一下,三重態(tài)優(yōu)化,我的方法對(duì)不對(duì)。謝謝 |

木蟲 (著名寫手)

新蟲 (初入文壇)
|
1. NStates=M 計(jì)算M個(gè)激發(fā)態(tài),默認(rèn)M=3 #p b3lyp/6-31g opt td(singlets,root=1) 計(jì)算3個(gè)激發(fā)態(tài) #p b3lyp/6-31g opt td(nstates=6,singlets,root=1) 計(jì)算6個(gè)激發(fā)態(tài) 您說的3個(gè)TD/6個(gè)TD不嚴(yán)謹(jǐn)。高斯只做一個(gè)TD計(jì)算,只是考慮的電子態(tài)數(shù)目不同。 2. 這個(gè)問題比較有意思,這里詳細(xì)講一下。 第一種 #p b3lyp/6-31g opt td(triplets,root=1) Title 0 1 這種方法是以基態(tài)單重態(tài)作為參考態(tài),優(yōu)化第一個(gè)激發(fā)態(tài)(三重態(tài)) 第二種 #p b3lyp/6-31g opt Title 0 3 這種方法是優(yōu)化三重態(tài)“基態(tài)” 兩種方式描述的都是同一個(gè)真實(shí)的電子態(tài)。該電子態(tài)在所有三重態(tài)中能量最低,所以是三重態(tài)的“基態(tài)”;此外,如果把單重態(tài)也考慮進(jìn)去,由于有能量更低的單重態(tài),以其為參考態(tài),我們要研究的該電子態(tài)又是激發(fā)態(tài)。 雖然描述的是同樣的電子態(tài),但兩種方法得到的結(jié)果(Hartree)并不相同!因?yàn)镈FT和TD-DFT都是近似的方法,而且采用了不同的一系列近似,相當(dāng)于使用兩種不同的方法,得到的數(shù)值自然不會(huì)相同。 推薦選擇第二種方法,即“基態(tài)”方法。使用基態(tài)方法可以優(yōu)化所有軌道能級(jí),而TD只能以基態(tài)作為參考態(tài)優(yōu)化特定的電子態(tài)。事實(shí)上有文獻(xiàn)試圖通過限制對(duì)稱性來(lái)使用基態(tài)方法來(lái)優(yōu)化其他激發(fā)態(tài)。 |
新蟲 (初入文壇)
|
木蟲 (著名寫手)
|
參考別人說的,我總結(jié)的,不知道是否正確,還請(qǐng)指正: 激發(fā)態(tài)的優(yōu)化 一:優(yōu)化s1 方法就一種: #p b3lyp/6-31g* opt td(singlets,root=1) 0 1 坐標(biāo) (此處坐標(biāo)為優(yōu)化后s0的結(jié)構(gòu),即基態(tài)) 說明: 1:對(duì)于此計(jì)算,是一種TDDFT的方法,默認(rèn)計(jì)算nstates=3, 與td(nstates=6,singlets,root=1)沒有區(qū)別。 2:計(jì)算可以加入溶劑:scrf=(solvent=n-hexane)?梢约尤胙h(huán)次次和收斂標(biāo)準(zhǔn):scf(maxcyc=400,conver=7) 3:計(jì)算出錯(cuò):You need to solve for more vectors in order to follow this state.可以在td中加入direct來(lái)解決。 第二:優(yōu)化t1 第一種方法:#p b3lyp/6-31g* opt td(triplets,root=1) 0 1 坐標(biāo) (為優(yōu)化后s0的坐標(biāo),即基態(tài)) 1:這種方法是以基態(tài)單重態(tài)作為參考態(tài),優(yōu)化第一個(gè)激發(fā)態(tài)(三重態(tài))。利用的TDDFT的方法。 2:默認(rèn)計(jì)算nstates=3。 第二種方法:#P b3lyp/6-31g* opt 0 3 坐標(biāo) (最好也是s0的坐標(biāo),但自旋多重度不同) 1:這種方法是優(yōu)化三重態(tài)“基態(tài)”。利用的為DFT的方法, 2:DFT和TDDFT計(jì)算的都是同一個(gè)真實(shí)的電子態(tài)。DFT計(jì)算的電子態(tài)在所有三重態(tài)中能量最低,所以是三重態(tài)的“基態(tài)”;此外,如果把單重態(tài)也考慮進(jìn)去,由于有能量更低的單重態(tài),以其為參考態(tài),我們要研究的該電子態(tài)又是激發(fā)態(tài)。 3:雖然描述的是同樣的電子態(tài),但兩種方法得到的結(jié)果(Hartree)并不相同!因?yàn)镈FT和TD-DFT都是近似的方法,而且采用了不同的一系列近似,相當(dāng)于使用兩種不同的方法,得到的數(shù)值自然不會(huì)相同。 4:推薦選擇第二種方法(DFT),即“基態(tài)”方法。使用基態(tài)方法可以優(yōu)化所有軌道能級(jí),而TD只能以基態(tài)作為參考態(tài)優(yōu)化特定的電子態(tài)。 第三種方法:#P b3lyp/6-31g* opt td(triplets,root=1) 0 3 坐標(biāo) 1:這種方式為,利用T1的結(jié)構(gòu)優(yōu)化T2的結(jié)構(gòu)。 2:自旋多重度為0 3,那么算出來(lái)的波函數(shù)是T1的,在做triplets的TDDFT,就會(huì)的到更高的三重態(tài),即T2(root=1) 3: 當(dāng)用三重態(tài)作為參考態(tài)的時(shí)候,triplet選項(xiàng)無(wú)效。Gaussian的U-TDDFT不能計(jì)算自旋反轉(zhuǎn),因此所有激發(fā)態(tài)也都是三重態(tài),這時(shí)候root=1對(duì)應(yīng)T2,參考態(tài)才是T1。另外,你只能得到近似的三重態(tài),因?yàn)樽孕窍拗艱FT有自旋污染。 4: 所以說:TDDFT是以基態(tài)為基礎(chǔ)計(jì)算激發(fā)態(tài)的。 不論激發(fā)態(tài)單線態(tài)還是三線態(tài),對(duì)于基態(tài)是單線態(tài)的分子,都應(yīng)該把自旋多重度設(shè)置為1。 |

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 294求調(diào)劑材料與化工專碩 +14 | 陌の森林 2026-03-18 | 14/700 |
|
|---|---|---|---|---|
|
[考研] 一志愿南京理工大學(xué)085701資源與環(huán)境302分求調(diào)劑 +3 | 葵梓衛(wèi)隊(duì) 2026-03-18 | 5/250 |
|
|
[考研] 288求調(diào)劑,一志愿華南理工大學(xué)071005 +5 | ioodiiij 2026-03-17 | 5/250 |
|
|
[考研] 本人考085602 化學(xué)工程 專碩 +17 | 不知道叫什么! 2026-03-15 | 19/950 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 324分 085600材料化工求調(diào)劑 +3 | llllkkkhh 2026-03-18 | 3/150 |
|
|
[考研] 求調(diào)劑,一志愿:南京航空航天大學(xué)大學(xué) ,080500材料科學(xué)與工程學(xué)碩,總分289分 +3 | @taotao 2026-03-19 | 3/150 |
|
|
[考研] 一志愿南昌大學(xué),327分,材料與化工085600 +3 | Ncdx123456 2026-03-19 | 3/150 |
|
|
[考研] 一志愿天大材料與化工(085600)總分338 +5 | 蔡大美女 2026-03-13 | 5/250 |
|
|
[考研] 0703化學(xué) 305求調(diào)劑 +4 | FY_yy 2026-03-14 | 4/200 |
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
[考研] 材料,紡織,生物(0856、0710),化學(xué)招生啦 +3 | Eember. 2026-03-17 | 9/450 |
|
|
[碩博家園] 湖北工業(yè)大學(xué) 生命科學(xué)與健康學(xué)院-課題組招收2026級(jí)食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|
[考研] 機(jī)械專碩325,尋找調(diào)劑院校 +3 | y9999 2026-03-15 | 5/250 |
|
|
[考研] 333求調(diào)劑 +3 | 文思客 2026-03-16 | 7/350 |
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-14 | 3/150 |
|
|
[考研] 321求調(diào)劑 +5 | 大米飯! 2026-03-15 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷,論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
[考研] 294求調(diào)劑 +3 | Zys010410@ 2026-03-13 | 4/200 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|