| 5 | 1/1 | 返回列表 |
| 查看: 2469 | 回復(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è)試過這兩種方式,得出來的HF能量是不同的,同時(shí),幾何結(jié)構(gòu)也不相同,我想知道一下,三重態(tài)優(yōu)化,我的方法對(duì)不對(duì)。謝謝 |

木蟲 (著名寫手)
|
參考別人說的,我總結(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來解決。 第二:優(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,那么算出來的波函數(shù)是T1的,在做triplets的TDDFT,就會(huì)的到更高的三重態(tài),即T2(root=1) 3: 當(dāng)用三重態(tài)作為參考態(tài)的時(shí)候,triplet選項(xiàng)無效。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。 |

新蟲 (初入文壇)
|
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ì)稱性來使用基態(tài)方法來優(yōu)化其他激發(fā)態(tài)。 |
木蟲 (著名寫手)

新蟲 (初入文壇)
|
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +3 | www_q 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
[考博] 申博26年 +3 | 八6八68 2026-03-19 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 ,六級(jí)已過,有科研經(jīng)歷 +12 | 曦熙兮 2026-03-15 | 12/600 |
|
|
[考研] 復(fù)試調(diào)劑 +4 | z1z2z3879 2026-03-14 | 6/300 |
|
|
[考研] 【考研調(diào)劑】化學(xué)專業(yè) 281分,一志愿四川大學(xué),誠心求調(diào)劑 +5 | 吃吃吃才有意義 2026-03-19 | 5/250 |
|
|
[考研] 一志愿西安交通大學(xué)材料工程專業(yè) 282分求調(diào)劑 +5 | 楓橋ZL 2026-03-18 | 7/350 |
|
|
[考研] 化學(xué)求調(diào)劑 +3 | 臨澤境llllll 2026-03-17 | 4/200 |
|
|
[考研] 0703化學(xué)調(diào)劑 +5 | pupcoco 2026-03-17 | 8/400 |
|
|
[考研] 材料專碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +5 | 梨花珞晚風(fēng) 2026-03-17 | 5/250 |
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
[考研] 0854可跨調(diào)劑,一作一項(xiàng)核心論文五項(xiàng)專利,省、國級(jí)證書40+數(shù)一英一287 +8 | 小李0854 2026-03-16 | 8/400 |
|
|
[考研] 280求調(diào)劑 +6 | 咕嚕曉曉 2026-03-18 | 7/350 |
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
[考研] 085601求調(diào)劑 +4 | Du.11 2026-03-16 | 4/200 |
|
|
[考研] 藥學(xué)383 求調(diào)劑 +3 | 藥學(xué)chy 2026-03-15 | 4/200 |
|
|
[基金申請(qǐng)]
今年的國基金是打分制嗎?
50+3
|
zhanghaozhu 2026-03-14 | 3/150 |
|
|
[考研] 070300化學(xué)學(xué)碩求調(diào)劑 +6 | 太想進(jìn)步了0608 2026-03-16 | 6/300 |
|
|
[考研] 22408總分284求調(diào)劑 +3 | InAspic 2026-03-13 | 3/150 |
|