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

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

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

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

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

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[碩博家園] 湖北工業(yè)大學 生命科學與健康學院-課題組招收2026級食品/生物方向碩士 +3 | 1喜春8 2026-03-17 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 材料工程專碩274一志愿211求調(diào)劑 +6 | 薛云鵬 2026-03-15 | 6/300 |
|
|
[考研] 302求調(diào)劑 +4 | 小賈同學123 2026-03-15 | 8/400 |
|
|
[考研] 藥學383 求調(diào)劑 +3 | 藥學chy 2026-03-15 | 4/200 |
|
|
[考研] 化學調(diào)劑0703 +8 | 啊我我的 2026-03-11 | 8/400 |
|
|
[考研] 070303一志愿西北大學學碩310找調(diào)劑 +5 | d如愿上岸 2026-03-12 | 8/400 |
|
|
[考研] 0856求調(diào)劑 +3 | 劉夢微 2026-03-15 | 3/150 |
|
|
[考博] 東華理工大學化材專業(yè)26屆碩士博士申請 +6 | zlingli 2026-03-13 | 6/300 |
|
|
[考研] 中科大材料專碩319求調(diào)劑 +3 | 孟鑫材料 2026-03-13 | 3/150 |
|
|
[考研] 復試調(diào)劑 +3 | 呼呼?~+123456 2026-03-14 | 3/150 |
|
|
[考研] 308 085701 四六級已過求調(diào)劑 +7 | 溫喬喬喬喬 2026-03-12 | 14/700 |
|
|
[考研] 求調(diào)劑(材料與化工327) +4 | 愛吃香菜啦 2026-03-11 | 4/200 |
|
|
[考研] 求材料調(diào)劑 085600英一數(shù)二總分302 前三科235 精通機器學習 一志愿哈工大 +4 | 林yaxin 2026-03-12 | 4/200 |
|
|
[考研] 0856材料與化工301求調(diào)劑 +5 | 奕束光 2026-03-13 | 5/250 |
|
|
[考研] 290求調(diào)劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
[碩博家園] 085600 260分求調(diào)劑 +3 | 天空還下雨么 2026-03-13 | 5/250 |
|
|
[考研] 290求調(diào)劑 +7 | ADT 2026-03-12 | 7/350 |
|
|
[論文投稿]
投稿問題
5+4
|
星光燦爛xt 2026-03-12 | 6/300 |
|
|
[考研] 290求調(diào)劑 +3 | ADT 2026-03-13 | 3/150 |
|
|
[考博] 福州大學楊黃浩課題組招收2026年專業(yè)學位博士研究生,2026.03.20截止 +3 | Xiangyu_ou 2026-03-12 | 3/150 |
|