 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 6 | 1/1 | 返回列表 |
| 查看: 1610 | 回復(fù): 5 | ||||
tephoon78木蟲 (正式寫手)
|
[交流]
【求助】如何計算在溶劑中的發(fā)射波長？ 已有3人參與
|
|
（1）在計算吸收波長時， Excited State 1: Singlet-A 3.2076 eV 386.53 nm f=0.0014 12 -> 13 -0.70615 意味著電子從12軌道到13軌道的躍遷。但是在使用td計算發(fā)射光譜的時候，計算結(jié)果也類似，也會出現(xiàn) Excited State 1: Singlet-A 2.159eV 326.53 nm f=0.0026 12 -> 13 -0.8396 那么這個12 -> 13 該如何解釋？ （2) 在Gaussian 09 user’s reference (page 249-252)給出的例子，讓人看不明白。例子的解釋如下： 第1步優(yōu)化基態(tài)幾何結(jié)構(gòu)和頻率。 第2步計算線性垂直激發(fā)能。 第3步分成2步，首先保留基態(tài)的溶劑效應(yīng)，并寫入chk文件；再次，從chk文件中讀溶劑效應(yīng)，然后計算線性垂直激發(fā)能。 第4步是激發(fā)態(tài)的幾何優(yōu)化。 第5步是計算頻率。 到這一步，前面的步驟都是可以理解的，下面的兩步讓人看不明 第6步保留激發(fā)的溶劑效應(yīng)，并寫入chk文件。 第7步使用B3LYP方法計算，而不是TD. 既然第4步是激發(fā)態(tài)幾何構(gòu)型優(yōu)化，那么它的溶劑效應(yīng)就是激發(fā)態(tài)的。在這一步可以直接得到波長。第6、7兩步純粹多余。 即使是第6步保存了溶劑效應(yīng)，第7步的B3lyp并不能計算得到波長，而是計算能量的。 誰能解釋一下？ 我將這個例子復(fù)制下來。 Fluoresence example: Emission (Fluorescence) from First Excited State (n→π*) of Acetaldehyde Here we study the cycle: Acetaldehyde Excitation and Emission Cycle The primary process of interest is the emission, but this example shows how to study the complete cycle including the solvent effects. Step 1: Ground state geometry optimization and frequencies (equilibrium solvation). This is a standard Opt Freq calculation on the ground state including PCM equilibrium solvation. %chk=01-ac # B3LYP/6-31+G(d,p) Opt Freq SCRF=(Solvent=Ethanol) Acetaldehyde ground state 0 1 C C,1,RA X,2,1.,1,A O,2,RB,3,A,1,180.,0 X,1,1.,2,90.,3,0.,0 H,1,R1,2,A1,5,0.,0 H,1,R23,2,A23,5,B23,0 H,1,R23,2,A23,5,-B23,0 H,2,R4,1,A4,3,180.,0 RA=1.53643 RB=1.21718 R1=1.08516 R23=1.08688 R4=1.10433 A=62.1511 A1=110.51212 A23=109.88119 A4=114.26114 B23=120.56468 Step 2: Vertical excitation with linear response solvation. This is a TD-DFT calculation of the vertical excitation, therefore at the ground state equilibrium geometry, with the default solvation: linear response, non-equilibrium. We perform a single-point TD-DFT calculation, which defaults to non-equilibrium solvation. The results of this job will be used to identify which state or states are of interest and their ordering. These results give a reasonable description of the solvation of the excited state, but not quite as good as that from a state-specific solvation calculation. In this case, we see that the n->π* state is the first excited state. Next, we will use the state-specific method to produce a better description of the vertical excitation step. %chk=02-ac # B3LYP/6-31+G(d,p) TD=NStates=6 SCRF=(Solvent=Ethanol) Geom=Check Guess=Read Acetaldehyde: linear response vertical excited states 0 1 Step 3: State-specific solvation of the vertical excitation. This will require two job steps: first the ground state calculation is done, specifying NonEq=write in the PCM input section, in order to store the information about non-equilibrium solvation based on the ground state. Second, the actual state-specific calculation is done, reading in the necessary information for non-equilibrium solvation using NonEq=read. %chk=03-ac # B3LYP/6-31+G(d,p) SCRF=(Solvent=Ethanol,Read) Geom=Check Guess=Read Acetaldehyde: prepare for state-specific non-eq solvation by saving the solvent reaction field from the ground state 0 1 NonEq=write --link1-- %chk=03-ac # B3LYP/6-31+G(d,p) TD(NStates=6,Root=1) SCRF=(Solvent=Ethanol,StateSpecific,Read) Geom=Check Guess=Read Acetaldehyde: read non-eq solvation from ground state and compute energy of the first excited with the state-specific method 0 1 NonEq=read Step 4: Relaxation of the excited state geometry. Next, we perform a TD-DFT geometry optimization, with equilibrium, linear response solvation, in order to find the minimum energy point on the excited state potential energy surface. Since this is a TD-DFT optimization, the program defaults to equilibrium solvation. As is typical of such cases, the molecule has a plane of symmetry in the ground state but the symmetry is broken in the excited state, so the ground state geometry is perturbed slightly to break symmetry at the start of the optimization. %chk=04-ac # B3LYP/6-31+G(d,p) TD=(Read,NStates=6,Root=1) SCRF=(Solvent=Ethanol) Geom=Modify Guess=Read Opt=RCFC Acetaldehyde: excited state opt Modify geometry to break Cs symmetry since first excited state is A" 0 1 4 1 2 3 10.0 5 1 2 7 -50.0 Step 5: Vibrational frequencies of the excited state structure. Now we run a frequency calculation to verify that the geometry located in step 4 is a minimum. The results could also be used as part of a Franck-Condon calculation if desired (see below). This is a numerical frequency calculation. %chk=05-ac # B3LYP/6-31+G(d,p) TD=(Read,NStates=6,Root=1) Freq SCRF=(Solvent=Ethanol) Geom=Check Guess=Read Acetaldehyde excited state freq 0 1 Step 6: Emission state-specific solvation (part 1). This step does state-specific equilibrium solvation of the excited state at its equilibrium geometry, writing out the solvation data for the next step via the PCM NonEq=write input. %chk=06-ac # B3LYP/6-31+G(d,p) TD=(Read,NStates=6,Root=1) SCRF=(Solvent=Ethanol,StateSpecific,Read) Geom=Check Guess=Read Acetaldehyde emission state-specific solvation at first excited state optimized geometry 0 1 NonEq=write Step 7: Emission to final ground state (part 2). Finally, we compute the ground state energy with non-equibrium solvation, at the excited state geometry and with the static solvation from the excited state. %chk=07-ac # B3LYP/6-31+G(d,p) SCRF=(Solvent=Ethanol,Read) Geom=Check Guess=Read Acetaldehyde: ground state non-equilibrium at excited state geometry. 0 1 NonEq=read Steps 1, 2, and 4 would be sufficient to compute the excitation and emission energies in the gas-phase (along with step 5 to confirm the nature of stationary point). They are not sufficient when solvent effects are included because the energies computed in step 4 correspond to the ground state solvent reaction field, while the emission takes place in the reaction field created in response to the excited state charge distribution. This is what is accounted for properly in steps 6 and 7. [ Last edited by tephoon78 on 2010-6-6 at 08:10 ] |
回復(fù)此樓» 收錄本帖的淘帖專輯推薦
 Gaussian模擬光譜 |
» 猜你喜歡
高效微納米儲能材料與產(chǎn)業(yè)化推進(jìn)實驗室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有5人回復(fù)
-
高效微納米儲能材料與產(chǎn)業(yè)化推進(jìn)實驗室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有0人回復(fù)
-
物理化學(xué)論文潤色/翻譯怎么收費?
已經(jīng)有298人回復(fù)
-
高效微納米儲能材料與產(chǎn)業(yè)化推進(jìn)實驗室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有2人回復(fù)
-
化學(xué)化工學(xué)院 招收化工、化學(xué)、材料等相關(guān)方向研究生（學(xué)碩、專碩都有調(diào)劑名額）
已經(jīng)有1人回復(fù)
-
高效微納米儲能材料與產(chǎn)業(yè)化推進(jìn)實驗室 蘇州 招收碩士研究生（含調(diào)劑）
已經(jīng)有0人回復(fù)
-
溫州大學(xué)新能源材料微觀溯源課題組招收2026級博士研究生（材料、化學(xué)方向）
已經(jīng)有6人回復(fù)
-
材料與化工 323 英一+數(shù)二+物化，一志愿：哈工大 本人本科雙一流
已經(jīng)有5人回復(fù)
-
26申博
已經(jīng)有2人回復(fù)
-
巴黎西岱大學(xué)（Université Paris Cité）Fayna MAMMERI 導(dǎo)師尋找CSC公派博士生
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 怎么把激發(fā)和發(fā)射波長弄到一張圖里，一半一半那種
已經(jīng)有9人回復(fù)
- 求助熒光計算發(fā)射波長值與實驗差別比較大 有50nm，怎么辦？
已經(jīng)有3人回復(fù)
- 【求助】羅丹明溶液中加入納米金和氯金酸溶液，熒光發(fā)射波長移至543nm，請教大家
已經(jīng)有3人回復(fù)
- 【求助】熒光檢測中出現(xiàn)激發(fā)波長比發(fā)射波長長的情況
已經(jīng)有8人回復(fù)
- 【求助】激發(fā)光譜下，發(fā)射波長一半處出現(xiàn)的強(qiáng)干擾峰。
已經(jīng)有13人回復(fù)
- 【求助】發(fā)射波長與激發(fā)波長會相差300nm嗎
已經(jīng)有7人回復(fù)
abbott
金蟲 (著名寫手)
不要用QQ問我東西
- 應(yīng)助: 16 (小學(xué)生)
- 金幣: 1071.6
- 散金: 2787
- 紅花: 10
- 帖子: 1015
- 在線: 105.2小時
- 蟲號: 452267
- 注冊: 2007-11-05
- 性別: GG
- 專業(yè): 計算機(jī)硬件技術(shù)
2樓2010-06-03 17:33:31
tephoon78
木蟲 (正式寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 1682
- 散金: 724
- 帖子: 677
- 在線: 130.3小時
- 蟲號: 330598
- 注冊: 2007-03-24
- 專業(yè): 理論和計算化學(xué)
3樓2010-06-04 14:01:05
faqianliu
銀蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 1042.7
- 散金: 75
- 紅花: 1
- 帖子: 165
- 在線: 180.6小時
- 蟲號: 263012
- 注冊: 2006-07-01
- 專業(yè): 無機(jī)納米化學(xué)
4樓2010-06-04 19:32:27
loovfnd
至尊木蟲 (著名寫手)
- 應(yīng)助: 35 (小學(xué)生)
- 貴賓: 1.111
- 金幣: 14361.4
- 散金: 236
- 紅花: 4
- 帖子: 2790
- 在線: 187小時
- 蟲號: 219674
- 注冊: 2006-03-14
- 專業(yè): 物理有機(jī)化學(xué)
5樓2011-03-26 15:40:32
pwzhou
鐵桿木蟲 (正式寫手)
- QC強(qiáng)帖: 2
- 應(yīng)助: 51 (初中生)
- 金幣: 7551.4
- 散金: 13
- 紅花: 14
- 帖子: 412
- 在線: 316.9小時
- 蟲號: 56225
- 注冊: 2005-01-07
- 性別: GG
- 專業(yè): 化學(xué)動力學(xué)
tephoon78(金幣+20): 2011-03-26 16:56:11
|
請看這個帖子：http://www.gaoyang168.com/bbs/viewthread.php?tid=2826231&fpage=1 我已經(jīng)在回復(fù)中給出了Gaussian公司關(guān)于這個例子的一些詳細(xì)的解釋。另外，有啥不明白的可以看看Gaussian在這個部分引用的3篇JCP的文獻(xiàn)，看看State－Specific和line Response到底有什么區(qū)別。 |
6樓2011-03-26 15:46:39
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 tephoon78 的主題更新
| 6 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 085600材料與化工306 +4 | z1z2z3879 2026-03-21 | 4/200 |
|
|---|---|---|---|---|
|
[考研] 考研調(diào)劑 +3 | 呼呼？~+123456 2026-03-21 | 3/150 |
|
|
|
[考研] 297求調(diào)劑 +3 | 喜歡還是不甘心 2026-03-20 | 3/150 |
|
|
|
[考研] 求助 +5 | 夢里的無言 2026-03-21 | 6/300 |
|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-18 | 3/150 |
|
|
|
[考研] 機(jī)械專碩299求調(diào)劑至材料 +3 | kkcoco25 2026-03-16 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +3 | wyyyqx 2026-03-17 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +10 | yy要上岸呀 2026-03-17 | 10/500 |
|
|
|
[考研] 一志愿重慶大學(xué)085700資源與環(huán)境專碩，總分308求調(diào)劑 +3 | 墨墨漠 2026-03-18 | 3/150 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-18 | 6/300 |
|
|
|
[考研] 材料專業(yè)求調(diào)劑 +6 | hanamiko 2026-03-18 | 6/300 |
|
|
|
[考研] 308求調(diào)劑 +3 | 阿姐阿姐家啊 2026-03-18 | 3/150 |
|
|
|
[考研] 324求調(diào)劑 +5 | lucky呀呀呀鴨 2026-03-20 | 5/250 |
|
|
|
[考研] 317求調(diào)劑 +5 | 申子申申 2026-03-19 | 9/450 |
|
|
|
[考博] 招收博士1-2人 +3 | QGZDSYS 2026-03-18 | 3/150 |
|
|
|
[考研] 一志愿中國海洋大學(xué)，生物學(xué)，301分，求調(diào)劑 +5 | 1孫悟空 2026-03-17 | 6/300 |
|
|
|
[考研] 一志愿985，本科211，0817化學(xué)工程與技術(shù)319求調(diào)劑 +10 | Liwangman 2026-03-15 | 10/500 |
|
|
|
[考研] 301求調(diào)劑 +4 | A_JiXing 2026-03-16 | 4/200 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 275求調(diào)劑 +4 | 太陽花天天開心 2026-03-16 | 4/200 |
|