 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 1625 | 回復(fù): 5 | ||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||||
tephoon78木蟲 (正式寫手)
|
[交流]
【求助】如何計(jì)算在溶劑中的發(fā)射波長(zhǎng)？ 已有3人參與
|
|||
|
（1）在計(jì)算吸收波長(zhǎng)時(shí)， Excited State 1: Singlet-A 3.2076 eV 386.53 nm f=0.0014 12 -> 13 -0.70615 意味著電子從12軌道到13軌道的躍遷。但是在使用td計(jì)算發(fā)射光譜的時(shí)候，計(jì)算結(jié)果也類似，也會(huì)出現(xiàn) Excited State 1: Singlet-A 2.159eV 326.53 nm f=0.0026 12 -> 13 -0.8396 那么這個(gè)12 -> 13 該如何解釋？ （2) 在Gaussian 09 user’s reference (page 249-252)給出的例子，讓人看不明白。例子的解釋如下： 第1步優(yōu)化基態(tài)幾何結(jié)構(gòu)和頻率。 第2步計(jì)算線性垂直激發(fā)能。 第3步分成2步，首先保留基態(tài)的溶劑效應(yīng)，并寫入chk文件；再次，從chk文件中讀溶劑效應(yīng)，然后計(jì)算線性垂直激發(fā)能。 第4步是激發(fā)態(tài)的幾何優(yōu)化。 第5步是計(jì)算頻率。 到這一步，前面的步驟都是可以理解的，下面的兩步讓人看不明 第6步保留激發(fā)的溶劑效應(yīng)，并寫入chk文件。 第7步使用B3LYP方法計(jì)算，而不是TD. 既然第4步是激發(fā)態(tài)幾何構(gòu)型優(yōu)化，那么它的溶劑效應(yīng)就是激發(fā)態(tài)的。在這一步可以直接得到波長(zhǎng)。第6、7兩步純粹多余。 即使是第6步保存了溶劑效應(yīng)，第7步的B3lyp并不能計(jì)算得到波長(zhǎng)，而是計(jì)算能量的。 誰能解釋一下？ 我將這個(gè)例子復(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模擬光譜 |
» 猜你喜歡
澳大利亞昆士蘭科技大學(xué)PhD博士獎(jiǎng)學(xué)金
已經(jīng)有0人回復(fù)
-
山西大學(xué)資源與環(huán)境工程研究所化工分離/環(huán)境功能材料小組招收2026年博士生
已經(jīng)有3人回復(fù)
-
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有137人回復(fù)
-
【211博士招生】環(huán)境化學(xué)、地學(xué)、毒理方向，擅長(zhǎng)污染物遷移轉(zhuǎn)化降解、計(jì)算模擬者優(yōu)先
已經(jīng)有17人回復(fù)
-
求助如何利用vasp計(jì)算磁光克爾效應(yīng)
已經(jīng)有1人回復(fù)
-
070300總分298求調(diào)劑
已經(jīng)有2人回復(fù)
-
一志愿廈門大學(xué)化學(xué)學(xué)碩307求調(diào)劑
已經(jīng)有10人回復(fù)
-
一志愿南師大0703化學(xué) 275求調(diào)劑
已經(jīng)有5人回復(fù)
-
理論與計(jì)算化學(xué)方向調(diào)劑——初試388分（數(shù)學(xué)144）
已經(jīng)有0人回復(fù)
-
A 區(qū)過線就可申！重慶雙非公辦高校 0856調(diào)劑，接受跨專業(yè)！納米材料 / 潤(rùn)滑熱門方向
已經(jīng)有0人回復(fù)
-
歡迎調(diào)劑長(zhǎng)江師范學(xué)院湯老師課題組的碩士研究生
已經(jīng)有0人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 怎么把激發(fā)和發(fā)射波長(zhǎng)弄到一張圖里，一半一半那種
已經(jīng)有9人回復(fù)
- 求助熒光計(jì)算發(fā)射波長(zhǎng)值與實(shí)驗(yàn)差別比較大 有50nm，怎么辦？
已經(jīng)有3人回復(fù)
- 【求助】羅丹明溶液中加入納米金和氯金酸溶液，熒光發(fā)射波長(zhǎng)移至543nm，請(qǐng)教大家
已經(jīng)有3人回復(fù)
- 【求助】熒光檢測(cè)中出現(xiàn)激發(fā)波長(zhǎng)比發(fā)射波長(zhǎng)長(zhǎng)的情況
已經(jīng)有8人回復(fù)
- 【求助】激發(fā)光譜下，發(fā)射波長(zhǎng)一半處出現(xiàn)的強(qiáng)干擾峰。
已經(jīng)有13人回復(fù)
- 【求助】發(fā)射波長(zhǎng)與激發(fā)波長(zhǎng)會(huì)相差300nm嗎
已經(jīng)有7人回復(fù)
tephoon78
木蟲 (正式寫手)
- 應(yīng)助: 4 (幼兒園)
- 金幣: 1682
- 散金: 724
- 帖子: 677
- 在線: 130.3小時(shí)
- 蟲號(hào): 330598
- 注冊(cè): 2007-03-24
- 專業(yè): 理論和計(jì)算化學(xué)
3樓2010-06-04 14:01:05
abbott
金蟲 (著名寫手)
不要用QQ問我東西
- 應(yīng)助: 16 (小學(xué)生)
- 金幣: 1071.6
- 散金: 2787
- 紅花: 10
- 帖子: 1015
- 在線: 105.2小時(shí)
- 蟲號(hào): 452267
- 注冊(cè): 2007-11-05
- 性別: GG
- 專業(yè): 計(jì)算機(jī)硬件技術(shù)
2樓2010-06-03 17:33:31
loovfnd
至尊木蟲 (著名寫手)
- 應(yīng)助: 35 (小學(xué)生)
- 貴賓: 1.111
- 金幣: 14361.4
- 散金: 236
- 紅花: 4
- 帖子: 2790
- 在線: 187小時(shí)
- 蟲號(hào): 219674
- 注冊(cè): 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小時(shí)
- 蟲號(hào): 56225
- 注冊(cè): 2005-01-07
- 性別: GG
- 專業(yè): 化學(xué)動(dòng)力學(xué)
tephoon78(金幣+20): 2011-03-26 16:56:11
|
請(qǐng)看這個(gè)帖子：http://www.gaoyang168.com/bbs/viewthread.php?tid=2826231&fpage=1 我已經(jīng)在回復(fù)中給出了Gaussian公司關(guān)于這個(gè)例子的一些詳細(xì)的解釋。另外，有啥不明白的可以看看Gaussian在這個(gè)部分引用的3篇JCP的文獻(xiàn)，看看State－Specific和line Response到底有什么區(qū)別。 |
6樓2011-03-26 15:46:39
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)354求調(diào)劑 +15 | 101次希望 2026-03-28 | 15/750 |
|
|---|---|---|---|---|
|
[考研] 考研生物與醫(yī)藥調(diào)劑 +6 | 鐵憨憨123425 2026-03-31 | 6/300 |
|
|
|
[考研] 機(jī)械學(xué)碩總分317求調(diào)劑�。。�！ +6 | Acaciad 2026-03-25 | 6/300 |
|
|
|
[考研] 本科211安全工程，初試290分，求調(diào)劑 +3 | 2719846834 2026-03-28 | 3/150 |
|
|
|
[考研] 289求調(diào)劑 +6 | BrightLL 2026-03-29 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +4 | 圖鑒212 2026-03-30 | 4/200 |
|
|
|
[考研] 福建理工大學(xué)材料學(xué)院先進(jìn)合金團(tuán)隊(duì)招收考研調(diào)劑學(xué)生 +3 | 大華金商都 2026-03-30 | 4/200 |
|
|
|
[考研] 327求調(diào)劑 +5 | 小卡不卡. 2026-03-29 | 5/250 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +6 | 丹青奶蓋 2026-03-26 | 8/400 |
|
|
|
[考研] 310求調(diào)劑 +10 | 爭(zhēng)取九點(diǎn)睡 2026-03-30 | 10/500 |
|
|
|
[考研] 考研調(diào)劑 +5 | Sanmu-124 2026-03-26 | 5/250 |
|
|
|
[考研] 【求調(diào)劑】085601材料工程專碩 | 總分272 | +7 | 腳滑的守法公民 2026-03-27 | 7/350 |
|
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +12 | c297914 2026-03-28 | 12/600 |
|
|
|
[考研] 321求調(diào)劑 +7 | 璞玉~~ 2026-03-25 | 8/400 |
|
|
|
[考研] 071000生物學(xué)求調(diào)劑，初試成績(jī)343 +7 | 小小甜面團(tuán) 2026-03-25 | 7/350 |
|
|
|
[考研] 一志愿南京航空航天大學(xué)材料學(xué)碩求調(diào)劑 +3 | @taotao 2026-03-28 | 3/150 |
|
|
|
[考研] 285求調(diào)劑 +4 | AZMK 2026-03-27 | 7/350 |
|
|
|
[考研] 08開頭275求調(diào)劑 +4 | 拉誰不重要 2026-03-26 | 4/200 |
|
|
|
[論文投稿] Journal of Mechanical Science and Technology +3 | Russ_ss 2026-03-25 | 5/250 |
|
|
|
[考研] 一志愿吉大071010，316分求調(diào)劑 +3 | xgbiknn 2026-03-27 | 3/150 |
|