| 查看: 2052 | 回復(fù): 20 | ||||
清風(fēng)凌云007金蟲 (小有名氣)
|
[求助]
用G09優(yōu)化激發(fā)態(tài)及熒光光譜計(jì)算 已有1人參與
|
|||
|
我要算一個(gè)體系的熒光光譜,看了好多小木蟲的帖子越看越糊涂了,還看了G09的英文說(shuō)明,在SCRF這有說(shuō)明(http://www.gaussian.com/g_tech/g_ur/k_scrf.htm),我只算氣態(tài)的,總結(jié)了下是不是需要step1、2、4、5、6?step7是否需要計(jì)算?希望計(jì)算高手們、前輩們給指點(diǎn)迷津,我這樣算到底對(duì)不對(duì),如果不對(duì)應(yīng)該怎么算呢,現(xiàn)在真是非常著急,都不想放假了,先在此謝過(guò)各位大神。∥业挠(jì)算步驟如下(現(xiàn)在在進(jìn)行第三步), 1、# opt freq b3lyp/6-31g(d) 基態(tài)分子結(jié)構(gòu)優(yōu)化及頻率計(jì)算 2、# td(nstates=20) b3lyp/6-31g(d) 吸收光譜計(jì)算 3、# opt td(nstates=20,root=1,direct) b3lyp/6-31g(d) 單重態(tài)第一激發(fā)態(tài)優(yōu)化(優(yōu)化基礎(chǔ)應(yīng)該用優(yōu)化基態(tài)結(jié)果還是吸收光譜計(jì)算結(jié)果?????) 4、# freq td(nstates=20,root=1,direct) b3lyp/6-31g(d) 單重態(tài)第一激發(fā)態(tài)優(yōu)化結(jié)構(gòu)頻率計(jì)算 5、# td(nstates=20) b3lyp/6-31g(d) 單重態(tài)第一激發(fā)態(tài)發(fā)射 |
經(jīng)驗(yàn) | 自己感覺好的東西 |
木蟲 (著名寫手)
|
手冊(cè)上的例子還有我們?cè)谡搲镉懻摰亩际窃谌軇┲械模枰紤]溶劑的平衡化和非平衡化的作用, 還有在TDDFT和DFT過(guò)程中,溶劑的作用。所以需要1-7,確切的說(shuō)是1.3.4.5.6.7這幾部。 如果說(shuō),你只計(jì)算在空氣中的發(fā)射譜,那么這個(gè)很簡(jiǎn)單,下面的步驟 1:# opt freq b3lyp/6-31g(d) 基態(tài)分子結(jié)構(gòu)優(yōu)化及頻率計(jì)算 2:# opt td(nstates=20,root=1,direct) b3lyp/6-31g(d) 單重態(tài)第一激發(fā)態(tài)優(yōu)化(優(yōu)化基礎(chǔ)應(yīng)該用優(yōu)化基態(tài)結(jié)果還是吸收光譜計(jì)算結(jié)果????? ps:針對(duì)你的問(wèn)題,是這樣的,這是激發(fā)態(tài)的優(yōu)化,如果說(shuō)計(jì)算發(fā)射,我們要看從s1-s0的發(fā)射,基于s0 優(yōu)化的結(jié)構(gòu)優(yōu)化s1 在進(jìn)行完上述步驟之后,在優(yōu)化s1的輸出文件中,會(huì)有如下的數(shù)據(jù): Excitation energies and oscillator strengths: Excited State 1: Singlet-A 1.7417 eV 711.84 nm f=0.0001 <S**2>=0.000 167 ->168 -0.70045 This state for optimization and/or second-order correction. Total Energy, E(TD-HF/TD-KS) = -2006.02553586 Copying the excited state density for this state as the 1-particle RhoCI density. Excited State 2: Singlet-A 2.5732 eV 481.82 nm f=0.1422 <S**2>=0.000 167 ->169 -0.69630 Excited State 3: Singlet-A 2.5993 eV 476.99 nm f=0.0028 <S**2>=0.000 165 ->168 0.44525 166 ->168 0.52944 默認(rèn)的是3個(gè),因?yàn)槟J(rèn)的nstates=3.看第一個(gè)td的波長(zhǎng)就是發(fā)射的熒光的波長(zhǎng) ps:很多人,還在優(yōu)化了s1的結(jié)構(gòu)后,又做了一個(gè)td,其實(shí)沒有必要,當(dāng)然你以可以進(jìn)行,計(jì)算出來(lái)的結(jié)果,和我們優(yōu)化s1出來(lái)的結(jié)構(gòu)的第一td是一樣的。 綜上所述,計(jì)算空氣中的熒光,只需要2步,優(yōu)化s0,基于s0優(yōu)化s1,然后從s1輸出文件找第一個(gè)td波長(zhǎng)。 over |

木蟲 (著名寫手)

金蟲 (小有名氣)
金蟲 (小有名氣)
鐵桿木蟲 (著名寫手)
木蟲 (著名寫手)
|
仔細(xì)看一看我的這個(gè)帖子。 http://www.gaoyang168.com/bbs/viewthread.php?tid=7928054 還有相關(guān)的鏈接。 http://www.gaoyang168.com/bbs/viewthread.php?tid=5546124&fpage=1&target=self&page=1 |

金蟲 (小有名氣)
送紅花一朵 |
多謝大神指教,我看了兩個(gè)帖子,還是有一點(diǎn)不明白,就是怎么做發(fā)射光譜圖呢,要是作圖是要在優(yōu)化激發(fā)態(tài)基礎(chǔ)上做TD 嗎?也就是我的最后一步 5、# td(nstates=20) b3lyp/6-31g(d) 單重態(tài)第一激發(fā)態(tài)發(fā)射 還有就是畫圖要用chk文件嗎?只有out或log文件應(yīng)該怎么畫圖呢?在我的輸入文件(如下)中輸入%chk=2.gif.chk就是生成了chk文件吧,但是我打不開,看了好多論壇上的帖子有的說(shuō)要轉(zhuǎn)化成fchk,但是我按照說(shuō)的轉(zhuǎn)化不了 %chk=2.chk %mem=100MW %nproc=1 # opt b3lyp/6-31+g |
金蟲 (小有名氣)
鐵桿木蟲 (著名寫手)
金蟲 (小有名氣)
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 +5 | 十三加油 2026-03-21 | 5/250 |
|
|---|---|---|---|---|
|
[考研] 材料學(xué)碩333求調(diào)劑 +3 | 北道巷 2026-03-18 | 3/150 |
|
|
[考研] 306求0703調(diào)劑一志愿華中師范 +5 | 紙魚ly 2026-03-21 | 5/250 |
|
|
[考研] 268求調(diào)劑 +9 | 簡(jiǎn)單點(diǎn)0 2026-03-17 | 9/450 |
|
|
[考研] 材料與化工(0856)304求 B區(qū) 調(diào)劑 +3 | 邱gl 2026-03-21 | 3/150 |
|
|
[考研] 材料工程(專)一志愿985 初試335求調(diào)劑 +3 | hiloiy 2026-03-17 | 4/200 |
|
|
[考研] 材料 336 求調(diào)劑 +3 | An@. 2026-03-18 | 4/200 |
|
|
[考研] 一志愿武漢理工材料工程專碩調(diào)劑 +9 | Doleres 2026-03-19 | 9/450 |
|
|
[考研] 一志愿蘇州大學(xué)材料求調(diào)劑,總分315(英一) +5 | sbdksD 2026-03-19 | 5/250 |
|
|
[考研] 一志愿北京化工大學(xué)0703化學(xué)318分,有科研經(jīng)歷,求調(diào)劑 +4 | 一瓶苯甲酸 2026-03-14 | 4/200 |
|
|
[考研] 295材料求調(diào)劑,一志愿武漢理工085601專碩 +5 | Charlieyq 2026-03-19 | 5/250 |
|
|
[考研] 086500 325 求調(diào)劑 +3 | 領(lǐng)帶小熊 2026-03-19 | 3/150 |
|
|
[論文投稿]
申請(qǐng)回稿延期一個(gè)月,編輯同意了。但系統(tǒng)上的時(shí)間沒變,給編輯又寫郵件了,沒回復(fù)
10+3
|
wangf9518 2026-03-17 | 4/200 |
|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +4 | 生物工程調(diào)劑 2026-03-16 | 12/600 |
|
|
[考研] 344求調(diào)劑 +6 | knight344 2026-03-16 | 7/350 |
|
|
[考研] 085601專碩,總分342求調(diào)劑,地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
[考研] 277調(diào)劑 +5 | 自由煎餅果子 2026-03-16 | 6/300 |
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
[考研] 一志愿,福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|