| 5 | 1/1 | 返回列表 |
| 查看: 4149 | 回復(fù): 14 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖,點(diǎn)擊這里查看本話題的所有回帖 | ||||
Illusionist銀蟲(chóng) (正式寫(xiě)手)
|
[求助]
激發(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.然后用這個(gè)優(yōu)化# opt cis=singlets/6-311g geom=connectivity 3.# cis=direct/6-311+g(d,p) td(singlets,nstates=6,root=1) 算的不對(duì) 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)” 這個(gè)計(jì)算。 不明白為什么cis算得就是不行。是我的設(shè)置還是順序什么的有問(wèn)題么? |
金蟲(chóng) (小有名氣)
很不愿意努力但必須努力的大喵

專(zhuān)家顧問(wèn) (職業(yè)作家)
地溝油冶煉專(zhuān)家
![]() |
專(zhuān)家經(jīng)驗(yàn): +458 |
|
你看的這份答案太老,現(xiàn)在過(guò)時(shí)了。 當(dāng)時(shí)gaussian程序的tddft沒(méi)有解析梯度,所以四五個(gè)原子以上的分子做激發(fā)態(tài)優(yōu)化比較難,cis速度更快。gaussian09的tddft已經(jīng)有了解析梯度,所以可以從頭至尾用dft和tddft。但是很多人沒(méi)買(mǎi)gaussian09,或者沒(méi)注意gaussian09的這個(gè)新功能,或者屬于知其然不知其所以然,仍然堅(jiān)持tddft/cis混用。 |

銀蟲(chóng) (正式寫(xiě)手)
|
那具體應(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感覺(jué)挺準(zhǔn)的。不過(guò)之前也進(jìn)行了cis的優(yōu)化。 迷茫啊···· 選擇應(yīng)助,我給你分。太感謝了 |
木蟲(chóng) (正式寫(xiě)手)

| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|