| 查看: 4142 | 回復(fù): 14 | ||||
Illusionist銀蟲 (正式寫手)
|
[求助]
激發(fā)態(tài)的homo和lumo
|
|||
|
小弟要計(jì)算激發(fā)態(tài)的homo和lumo 看了論壇的帖子 http://www.gaoyang168.com/bbs/viewthread.php?tid=1459478 但是還是有些不明白 1、基態(tài)構(gòu)型優(yōu)化:# b3lyp/6-31g(d) opt freq 2、垂直激發(fā)計(jì)算:# b3lyp/6-31g(d) td=(singlets,nstates=6,root=1)。坐標(biāo)部分采用1優(yōu)化后的坐標(biāo) 3、激發(fā)態(tài)結(jié)構(gòu)優(yōu)化:# rcis(root=1,singlets) opt freq。坐標(biāo)部分采用2的坐標(biāo) 4、再一次TD的話就可以找出熒光譜和磷光譜了。 不明的地方: 我的做法是: 1、基態(tài)構(gòu)型優(yōu)化:# b3lyp/6-31g(d) opt freq 2.然后用這個優(yōu)化# opt cis=singlets/6-311g geom=connectivity 3.# cis=direct/6-311+g(d,p) td(singlets,nstates=6,root=1) 算的不對 4.再用# cis=(nstates=6,root=1,direct)/6-311+g(d,p) geom=connectivity 也不行 5.我用# b3lyp/6-31+g(d) td(singlets,nstates=6,root=1)結(jié)果差不多 為什么一定要用 “2、垂直激發(fā)計(jì)算:# b3lyp/6-31g(d) td=(singlets,nstates=6,root=1)。坐標(biāo)部分采用1優(yōu)化后的坐標(biāo)” 這個計(jì)算。 不明白為什么cis算得就是不行。是我的設(shè)置還是順序什么的有問題么? |
專家顧問 (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經(jīng)驗(yàn): +458 |

銀蟲 (正式寫手)
專家顧問 (職業(yè)作家)
地溝油冶煉專家
![]() |
專家經(jīng)驗(yàn): +458 |
|
第一步應(yīng)該用dft。這一步確定基態(tài)平衡結(jié)構(gòu) 第二步在基態(tài)平衡結(jié)構(gòu)計(jì)算激發(fā)能,根據(jù)Franck-Condon原理,近似對應(yīng)電子吸收譜的能量 第三步用tddft優(yōu)化第一激發(fā)態(tài),最后得到的激發(fā)能近似對應(yīng)這個電子態(tài)的電子發(fā)射譜能量。這里的的freq是數(shù)值方法(因?yàn)間aussian的tddft目前還沒實(shí)現(xiàn)解析二階導(dǎo)數(shù)),相當(dāng)耗時,如果能肯定是能量極小點(diǎn)(例如,根據(jù)文獻(xiàn),或者沒有使用對稱性而且也沒對結(jié)構(gòu)加上限制),沒必要做freq計(jì)算。 HF和DFT的標(biāo)度遠(yuǎn)比post-HF低得多,因此如果體系不是太大(例如,800基函數(shù)以下),可以從頭至尾都用最大的基組做DFT/TDDFT計(jì)算。這樣可以省去第四步。 |

金蟲 (小有名氣)
很不愿意努力但必須努力的大喵
送紅花一朵 |
樓主,我還是想問一下您,對于激發(fā)態(tài)我確實(shí)有太多的問題,您在 1. # opt b3lyp/6-31+(d,p) geom=connectivity 是對基態(tài)優(yōu)化,其相應(yīng)的軌道計(jì)算也就是: # b3lyp/6-31+(d,p) pop=full geom=connectivity 而在 3. # opt td=(nstates=6,root=1)/6-31+g(d,p) geom=check guess=read 是對激發(fā)態(tài)進(jìn)行優(yōu)化,那么其對應(yīng)的軌道計(jì)算也就是: # td=(nstates=6,root=1)/6-31+g(d,p) pop=full geom=check guess=read 我這樣計(jì)算是對的嗎?還是說激發(fā)態(tài)的前線軌道還有別的方法? |

銀蟲 (正式寫手)
銀蟲 (正式寫手)
銀蟲 (正式寫手)
|
那具體應(yīng)該怎么做啊 只要用cis就可以了么 還是 # td=(nstates=6,root=1) b3lyp/6-311g geom=connectivity 那結(jié)構(gòu)優(yōu)化也要用td=(nstates=6,root=1) b3lyp/6-311g ? 小弟新手,但是有點(diǎn)急。昨晚用 # td=(nstates=6,root=1) b3lyp/6-311++g(d,p) geom=connectivity算激發(fā)態(tài)homo和lumo感覺挺準(zhǔn)的。不過之前也進(jìn)行了cis的優(yōu)化。 迷茫啊···· 選擇應(yīng)助,我給你分。太感謝了 |
木蟲 (正式寫手)

銀蟲 (正式寫手)
|
第一步 優(yōu)化是指優(yōu)化 基態(tài)結(jié)構(gòu)吧,用這個算法應(yīng)該夠了吧 # opt freq hf/6-311g geom=connectivity 第二步td 第二部td的意義是什么。是為了計(jì)算激發(fā)態(tài),通過td得到一個較粗糙的激發(fā)態(tài)模型 那么是否該這么寫 # td=(nstates=6,root=1)/6-31+g(d,p) geom=connectivity 第三步 進(jìn)行更精確的模型 # opt freq td=(nstates=6,root=1) b3lyp/6-311g geom=connectivity 第四步 # td=(nstates=6,root=1) b3lyp/6-311+g(d,p) geom=connectivity |
木蟲 (正式寫手)

銀蟲 (正式寫手)
|
第一步 確定基態(tài)平衡結(jié)構(gòu) # opt freq b3lyp/6-311g geom=connectivity 第二部 在基態(tài)平衡結(jié)構(gòu)計(jì)算激發(fā)能,根據(jù)Franck-Condon原理,近似對應(yīng)電子吸收譜的能量。 # td=(nstates=6,root=1) b3lyp/6-311++g(d,p) geom=connectivity 第三步 用tddft優(yōu)化第一激發(fā)態(tài),最后得到的激發(fā)能近似對應(yīng)這個電子態(tài)的電子發(fā)射譜能量。這里的的freq是數(shù)值方法(因?yàn)間aussian的tddft目前還沒實(shí)現(xiàn)解析二階導(dǎo)數(shù)),相當(dāng)耗時,如果能肯定是能量極小點(diǎn)(例如,根據(jù)文獻(xiàn),或者沒有使用對稱性而且也沒對結(jié)構(gòu)加上限制),沒必要做freq計(jì)算。 # opt freq td=(nstates=6,root=1) b3lyp/6-311++g(d,p) geom=connectivity 問一下 第二步的結(jié)果對第三步的影響有多大? 不使用擴(kuò)散函數(shù)的結(jié)果會不會受影響。即可否使用 # td=(nstates=6,root=1) b3lyp/6-311g geom=connectivity 這樣比較節(jié)約時間 非常謝謝 |
金蟲 (小有名氣)
很不愿意努力但必須努力的大喵

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章 +3 | rare12345 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 299求調(diào)劑 +5 | △小透明* 2026-03-17 | 5/250 |
|
|
[考研] 290求調(diào)劑 +10 | ADT 2026-03-11 | 10/500 |
|
|
[考研] 312求調(diào)劑 +7 | 陌宸希 2026-03-16 | 8/400 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +4 | share_joy 2026-03-16 | 4/200 |
|
|
[考研] 考研化學(xué)學(xué)碩調(diào)劑,一志愿985 +4 | 張vvvv 2026-03-15 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | p asserby. 2026-03-15 | 4/200 |
|
|
[考研] 26考研求調(diào)劑 +6 | 丶宏Sir 2026-03-13 | 6/300 |
|
|
[考研] 08工科 320總分 求調(diào)劑 +4 | 梨花珞晚風(fēng) 2026-03-17 | 4/200 |
|
|
[考研] 東南大學(xué)364求調(diào)劑 +5 | JasonYuiui 2026-03-15 | 5/250 |
|
|
[考研] 286求調(diào)劑 +3 | lemonzzn 2026-03-16 | 5/250 |
|
|
[考研] 0703化學(xué)調(diào)劑,求各位老師收留 +8 | 秋有木北 2026-03-14 | 8/400 |
|
|
[考研] 327求調(diào)劑 +6 | 拾光任染 2026-03-15 | 11/550 |
|
|
[考研]
|
笨笨兔子 2026-03-12 | 3/150 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
[考研] 材料與化工求調(diào)劑一志愿 985 總分 295 +8 | dream…… 2026-03-12 | 8/400 |
|
|
[考研] 求b區(qū)學(xué)校調(diào)劑 +3 | 周56 2026-03-11 | 3/150 |
|
|
[考研] 314求調(diào)劑 +7 | 無懈可擊的巨人 2026-03-12 | 7/350 |
|
|
[考研] 0856化學(xué)工程280分求調(diào)劑 +4 | shenzxsn 2026-03-11 | 4/200 |
|
|
[考研] 333求調(diào)劑 +3 | 152697 2026-03-12 | 4/200 |
|