| 24小時(shí)熱門(mén)版塊排行榜 |
| 10 | 1/1 | 返回列表 |
| 查看: 3102 | 回復(fù): 9 | ||||
小范范1989木蟲(chóng) (著名寫(xiě)手)
|
[交流]
高斯09優(yōu)化激發(fā)態(tài)并且計(jì)算激發(fā)態(tài)頻率的討論與求教。
|
|
先說(shuō)一下我的背景:點(diǎn)上計(jì)算資源有限,50多個(gè)人用一個(gè)很老的集群,有些老師感覺(jué)太滿,都自己帶電腦,用單機(jī)跑分子。 然后,我要計(jì)算的分子大約含有原子110個(gè),曾經(jīng),我計(jì)算一個(gè)分子(大約80個(gè)原子),激發(fā)態(tài)的頻率計(jì)算了半年,出來(lái)結(jié)果還有一個(gè)虛頻,我當(dāng)時(shí)的心,就和現(xiàn)在只穿一個(gè)背心逛街感覺(jué)是一樣的。 我的計(jì)算的最終目的是:我需要的是基態(tài)的頻率,還有激發(fā)態(tài)的頻率。 其次:說(shuō)一下我遇到的問(wèn)題 1:基態(tài)的優(yōu)化和頻率計(jì)算都是在b3lyp/6-31g*的基礎(chǔ)上完成的,還算很順利,沒(méi)有虛頻。耗時(shí)也不是很久。 2:優(yōu)化激發(fā)態(tài),一直出錯(cuò),采用的是b3lyp/6-31g*.我最近好多帖子都是問(wèn)這個(gè)問(wèn)題,相信很多蟲(chóng)友已經(jīng)見(jiàn)過(guò)了,基本都能看到我的陰魂。我在想,激發(fā)態(tài)優(yōu)化,我都優(yōu)化2個(gè)月,如果計(jì)算激發(fā)態(tài)的頻率,是不是要等我博6畢業(yè)才能出來(lái)?所以我等不起。 3:看了sob的帖子,說(shuō)激發(fā)態(tài)的優(yōu)化,可以采用pbe0,或者是m062x等,我這就不過(guò)多的累述了。給sob一個(gè)郵件,sob的回復(fù)是: 至于頻率計(jì)算,如果你在TDDFT下進(jìn)行,由于Gaussian沒(méi)有TDDFT解析二階導(dǎo)數(shù),這是幾乎做不動(dòng)的,耗時(shí)我估計(jì)是600~1200次激發(fā)能計(jì)算的耗時(shí)。必須要計(jì)算頻率的話,只能是在CIS下優(yōu)化并計(jì)算頻率,因?yàn)镃IS下有解析二階導(dǎo)數(shù),速度和TDDFT計(jì)算頻率耗時(shí)完全不在一個(gè)數(shù)量級(jí)。頻率對(duì)計(jì)算級(jí)別的敏感度遠(yuǎn)小于激發(fā)能,因此CIS下的頻率的準(zhǔn)確度應(yīng)當(dāng)沒(méi)問(wèn)題。 這是的想問(wèn)的主要問(wèn)題:sob建議我激發(fā)態(tài)的優(yōu)化和頻率的計(jì)算使用CIS是吧?我第一感覺(jué),是不是CIS不更耗時(shí)?再就是,我激發(fā)態(tài)用cis,基態(tài)用b3lyp,然后拿這兩個(gè)頻率去進(jìn)行下一步的計(jì)算,(清華的momoap計(jì)算杜欣斯基轉(zhuǎn)動(dòng)矩陣)。能靠譜嗎? 還有一事就是,CIS計(jì)算,就是把我的b3lyp換成CIS就行了是吧?其他的輸入文件的地方不用改? 謝謝指點(diǎn),謝謝。 |
經(jīng)驗(yàn) |

木蟲(chóng) (著名寫(xiě)手)

木蟲(chóng) (著名寫(xiě)手)
至尊木蟲(chóng) (文壇精英)

木蟲(chóng) (著名寫(xiě)手)
|
謝謝槍神的回復(fù),我打算按照sob的指點(diǎn)試試,激發(fā)態(tài)我用cis優(yōu)化,并且算頻率。但是對(duì)于另外的參數(shù),比如說(shuō)絕熱激發(fā)能和垂直激發(fā)能,我是不是應(yīng)該按照我的b3lyp/6-31g*的形式去找?但是,激發(fā)態(tài)6-31g*出不來(lái),3-21g出來(lái),能否用基態(tài)是b3lyp/6-31g*,激發(fā)態(tài)優(yōu)化是3-21g,來(lái)找絕熱激發(fā)能? 可是我感覺(jué),我整個(gè)的計(jì)算的流程,基態(tài)和激發(fā)態(tài)優(yōu)化的方式不一樣,基泰和激發(fā)態(tài)頻率的方式不一樣,在進(jìn)行下一步,借助清華momap去算輻射躍遷,無(wú)輻射躍遷速率,這樣一系列的省略下來(lái),不知道最后的結(jié)果和實(shí)驗(yàn)差不多少。一直擔(dān)心這個(gè)。 謝謝神槍的指點(diǎn),后面的如果運(yùn)到其他我不能單獨(dú)解決的問(wèn)題,還希望您指點(diǎn)一二,thanks |

木蟲(chóng) (著名寫(xiě)手)

木蟲(chóng) (著名寫(xiě)手)
|
那6-31g可以出得來(lái)嗎?垂直激發(fā)能應(yīng)該沒(méi)問(wèn)題,畢竟不涉及結(jié)構(gòu)優(yōu)化。絕熱激發(fā)能是怎么個(gè)情況?我在理論上確實(shí)小白。。。 就像sob神說(shuō)的,如果你單單只需要頻率數(shù)據(jù)的話其精度的影響不會(huì)太大。 另外不好意思我看錯(cuò)了一點(diǎn),你算的應(yīng)該是三重態(tài)的基態(tài),也就是T0,如果是T1的話應(yīng)該用td(triplet)來(lái)算。T1的計(jì)算應(yīng)該也有很多老帖討論過(guò)的 最后cis的寫(xiě)法是替換td,不是b3lyp,如果我沒(méi)記錯(cuò)的話 |
木蟲(chóng) (著名寫(xiě)手)
|
謝謝神槍指點(diǎn),我還有幾個(gè)問(wèn)題: 1:絕熱激發(fā)能,就是s1的穩(wěn)態(tài)能量和s0的穩(wěn)態(tài)能量之差,我當(dāng)時(shí)問(wèn)過(guò),這個(gè)帖子:http://ccc.keinsci.com/forum.php ... =365&extra= 2:那我就采用sob的提議,CIS計(jì)算頻率。 3:難道以前我都理解錯(cuò)了?我以為我的這個(gè)方式就是優(yōu)化T1, sob的回復(fù):又是三重態(tài)又是root=1,直接用0 3計(jì)算就行了,不需要TDDFT,又快又不會(huì)出現(xiàn)這些問(wèn)題 我當(dāng)時(shí)在這個(gè)帖子問(wèn)過(guò): http://ccc.keinsci.com/forum.php ... =369&extra= 再就是,三重態(tài)基態(tài),是T0?T0沒(méi)見(jiàn)過(guò)啊,我以前都認(rèn)為是t1,還請(qǐng)槍神指點(diǎn)。 4:我昨晚上看了看,是不是這樣寫(xiě) #p RCIS(root=1,nstates=40)/6-31g* opt freq |

木蟲(chóng) (著名寫(xiě)手)
木蟲(chóng) (著名寫(xiě)手)

| 10 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 352求調(diào)劑 +3 | 大米飯! 2026-03-22 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 277材料科學(xué)與工程080500求調(diào)劑 +7 | 自由煎餅果子 2026-03-16 | 7/350 |
|
|
[考研] 293求調(diào)劑 +3 | 濤濤Wjt 2026-03-22 | 5/250 |
|
|
[考研] 298求調(diào)劑 +6 | 上岸6666@ 2026-03-20 | 6/300 |
|
|
[考研] 一志愿華中農(nóng)業(yè)071010,總分320求調(diào)劑 +5 | 困困困困坤坤 2026-03-20 | 6/300 |
|
|
[考研] 298求調(diào)劑一志愿211 +3 | 上岸6666@ 2026-03-20 | 3/150 |
|
|
[考研] 291 求調(diào)劑 +3 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 3/150 |
|
|
[考研] 0856材料專碩353求調(diào)劑 +4 | NIFFFfff 2026-03-20 | 4/200 |
|
|
[考研] 材料求調(diào)劑 +5 | @taotao 2026-03-21 | 5/250 |
|
|
[考研] 268求調(diào)劑 +9 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 9/450 |
|
|
[考研] 22 350 本科985求調(diào)劑,求老登收留 +3 | 李軼男003 2026-03-20 | 3/150 |
|
|
[考研] 材料學(xué)學(xué)碩080502 337求調(diào)劑-一志愿華中科技大學(xué) +4 | 順順順mr 2026-03-18 | 5/250 |
|
|
[考研] 求調(diào)劑 +6 | Mqqqqqq 2026-03-19 | 6/300 |
|
|
[考研] 化學(xué)求調(diào)劑 +4 | 臨澤境llllll 2026-03-17 | 5/250 |
|
|
[考研] 一志愿西南交大,求調(diào)劑 +5 | 材化逐夢(mèng)人 2026-03-18 | 5/250 |
|
|
[考研] 求調(diào)劑 +3 | eation27 2026-03-20 | 3/150 |
|
|
[考研] 0703化學(xué)調(diào)劑 +4 | 18889395102 2026-03-18 | 4/200 |
|
|
[考研] 286求調(diào)劑 +6 | lemonzzn 2026-03-16 | 10/500 |
|
|
[考研] 考研調(diào)劑 +3 | 淇ya_~ 2026-03-17 | 5/250 |
|
|
[考研] [導(dǎo)師推薦]西南科技大學(xué)國(guó)防/材料導(dǎo)師推薦 +3 | 尖角小荷 2026-03-16 | 6/300 |
|