 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 4647 | 回復(fù): 12 | |||||||||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | |||||||||
戴世杰木蟲 (著名寫手)
|
[交流]
用高斯09計(jì)算手冊(cè)中乙醛在乙醇溶液中的熒光光譜的問(wèn)題
|
||||||||
|
參照手冊(cè)，把乙醛在乙醇溶液中的熒光光譜計(jì)算了一下，但有很多不清楚的地方，希望大家指點(diǎn)。 第一步：基態(tài)的結(jié)構(gòu)優(yōu)化和頻率（平衡溶劑化）。這是對(duì)基態(tài)包含PCM平衡溶劑化的標(biāo)準(zhǔn)Opt Freq計(jì)算。 第二步：包含線性響應(yīng)溶劑化的垂直激發(fā)。這是垂直激發(fā)的TD-DFT計(jì)算，因此在基態(tài)平衡結(jié)構(gòu)計(jì)算，使用默認(rèn)溶劑，線性響應(yīng)，非平衡溶劑化。我們進(jìn)行單點(diǎn)TD-DFT計(jì)算，這對(duì)非平衡溶劑化是默認(rèn)的。 該任務(wù)的結(jié)果將用于識(shí)別感興趣的態(tài)及其順序。 這些結(jié)果給出激發(fā)態(tài)溶劑化的合理描述，但不如特定態(tài)溶劑化計(jì)算的結(jié)果好。在本例中，我們看到n→π*態(tài)是第一激發(fā)態(tài)。下一步，我們將使用特定態(tài)方法，對(duì)垂直激發(fā)步驟產(chǎn)生更好的描述。 我的問(wèn)題：1.什么是線性響應(yīng)溶劑化的垂直激發(fā)？2.什么是平衡和非平衡溶劑化 第三步：垂直激發(fā)的特定態(tài)溶劑化。這需要兩個(gè)任務(wù)步驟：首先做基態(tài)計(jì)算，為了存儲(chǔ)基態(tài)的非平衡溶劑化信息，需要在PCM輸入部分指定NonEq=write。接下來(lái)做實(shí)際的特定態(tài)的計(jì)算，需要用NonEq=read讀取非平衡溶劑化的必要信息。 我的問(wèn)題：什么是特定態(tài)溶劑化？ 第四步：放松激發(fā)態(tài)結(jié)構(gòu)。接下來(lái)，我們用平衡的線性響應(yīng)溶劑化進(jìn)行TD-DFT幾何優(yōu)化，以找到激發(fā)態(tài)勢(shì)能面上的最小能量點(diǎn)。由于這是TD-DFT優(yōu)化，程序默認(rèn)為平衡溶劑化。作為這類體系的典型情況，基態(tài)具有平面對(duì)稱性，而在激發(fā)態(tài)對(duì)稱性被破壞，因此在優(yōu)化的開始對(duì)基態(tài)結(jié)構(gòu)作了略微的擾動(dòng)以打破對(duì)稱性。 第五步：激發(fā)態(tài)結(jié)構(gòu)的振動(dòng)頻率。現(xiàn)在我們通過(guò)運(yùn)行頻率計(jì)算來(lái)確認(rèn)第四步找到的結(jié)構(gòu)是能量極小點(diǎn)。 第六步：溶液中特定態(tài)的發(fā)射（第一部分）。這一步在激發(fā)態(tài)平衡構(gòu)型計(jì)算特定態(tài)的平衡溶劑化，并通過(guò)輸入中的PCM選項(xiàng)NonEq=write，把溶劑化數(shù)據(jù)的輸出寫到文件，用于下一步計(jì)算。 第七步：到最終基態(tài)的發(fā)射（第二部分）。最后，我們用非平衡溶劑化和來(lái)自激發(fā)態(tài)的靜態(tài)溶劑化，在激發(fā)態(tài)結(jié)構(gòu)計(jì)算基態(tài)能量。 因此我總結(jié)了一下，溶液中熒光光譜的計(jì)算步驟：首先做基態(tài)結(jié)構(gòu)的平衡溶劑化;然后做基態(tài)的垂直激發(fā)，這一步是非平衡溶劑化的，用于識(shí)別感興趣的態(tài)及其順序;并接著做特定態(tài)溶劑化計(jì)算以更好地描述，首先要存儲(chǔ)基態(tài)的非平衡溶劑化信息，因?yàn)檫@個(gè)信息用來(lái)做特定態(tài)溶劑化計(jì)算;緊接著做激發(fā)態(tài)的平衡溶劑化;最后做發(fā)射就行了。之所以比氣相中復(fù)雜，主要是由于基態(tài)的垂直激發(fā)和激發(fā)態(tài)的發(fā)射都受到溶劑化作用。 附件是我計(jì)算的輸入，輸出和點(diǎn)文件 ，有什么不對(duì)的地方請(qǐng)大家指正。 |
回復(fù)此樓» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問(wèn)題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : td.zip
2014-08-17 17:32:39, 13.57 M
» 收錄本帖的淘帖專輯推薦
 下載 | Gaussian模擬光譜 | 經(jīng)驗(yàn) | 高斯相關(guān) |
| 熒光計(jì)算 | 量子化學(xué)基礎(chǔ)問(wèn)題 |
» 猜你喜歡
2026湖北師范大學(xué)化學(xué)化工學(xué)院招收研究生調(diào)劑
已經(jīng)有0人回復(fù)
-
26博士申請(qǐng)
已經(jīng)有0人回復(fù)
-
物理化學(xué)論文潤(rùn)色/翻譯怎么收費(fèi)?
已經(jīng)有85人回復(fù)
-
26年材料/化學(xué)相關(guān)專業(yè)博士研究生招生 溫州大學(xué)化學(xué)與材料工程學(xué)院(26年秋季入學(xué))
已經(jīng)有0人回復(fù)
-
0703調(diào)劑
已經(jīng)有2人回復(fù)
-
2026年博士研究生招生
已經(jīng)有1人回復(fù)
-
304求調(diào)劑
已經(jīng)有6人回復(fù)
-
關(guān)于研究生調(diào)劑的些許建議
已經(jīng)有1人回復(fù)
-
塑料油分析測(cè)試
已經(jīng)有1人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 高斯重復(fù)乙醛的原子光譜程序出錯(cuò)
已經(jīng)有5人回復(fù)
- 如何計(jì)算熒光光譜
已經(jīng)有11人回復(fù)
- 【求助】G09熒光光譜例子 運(yùn)行出錯(cuò)
已經(jīng)有25人回復(fù)
小范范1989
木蟲 (著名寫手)
- 應(yīng)助: 123 (高中生)
- 金幣: 3821.3
- 散金: 968
- 紅花: 243
- 帖子: 1979
- 在線: 452.8小時(shí)
- 蟲號(hào): 2013403
- 注冊(cè): 2012-09-19
- 性別: GG
- 專業(yè): 原子和分子物理
3樓2014-09-22 17:08:47
hairan
木蟲 (著名寫手)
瞌睡木蟲
- QC強(qiáng)帖: 3
- 應(yīng)助: 61 (初中生)
- 貴賓: 1.06
- 金幣: 4974.3
- 散金: 8072
- 紅花: 44
- 帖子: 2901
- 在線: 733.4小時(shí)
- 蟲號(hào): 432069
- 注冊(cè): 2007-08-12
- 性別: GG
- 專業(yè): 理論和計(jì)算化學(xué)
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
gmy1990: 金幣+5 2014-08-18 02:26:33
戴世杰: 金幣+8 2014-08-21 12:40:20
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
gmy1990: 金幣+5 2014-08-18 02:26:33
戴世杰: 金幣+8 2014-08-21 12:40:20
|
我以前回復(fù)過(guò)兩次類似的問(wèn)題，但找不到相關(guān)的帖子了。 下面的英文是我問(wèn)Gaussian官方技術(shù)支持關(guān)于熒光計(jì)算的問(wèn)題，也是基于說(shuō)明書的算例的問(wèn)題。 樓主耐心點(diǎn)看看下面的英文。 應(yīng)該能解決你大部分的問(wèn)題。 *************************************************************************** 1. A common TD job will give the enough information of verticle excitation and emission in gas phase, but a statespecific job need to be run for that in solvation. I want to the function of the statespecific job. （這一句是我問(wèn)的問(wèn)題） The "StateSpecific" approach in this context applies to the excited state, and yes it involves solving self-consistently the "fast" component of the solvent polarization for the target state and non-equilibrium solvation for the "slow" component of the solvent polarization (i.e. the "slow" component of the solvent polarization comes from the origin state). The total polarization is always partitioned into two components, "slow" and "fast". The "slow" part can be regarded as the reorganization of the solvent molecules as a response to a change in the electronic density of the solute. The "fast" part can be regarded as the response of the electrons in the solvent to a change in the electronic density of the solute. For a change in the electronic density of the solute such as a vertical electronic transition, the "slow" component of the polarization is much slower than the timescale of the electronic transition, so the solvent does not have time to respond in this way to the vertical electronic transition. The "fast" component of the polarization, on the other hand, is closer in timescale to the vertical electronic transition on the solute. In an equilibrium solvation calculation, both components are in equilibrium with the solute's density. By default, all ground state calculations assume equilibrium solvation, as well as geometry optimizations of excited states (and also any calculation that involves the computation of the relaxed density of the excited state). Again, in equilibrium solvation processes, both the "slow" and "fast" component of the solvent polarization are in equilibrium with the excited state density. In the case of a TD energy calculation (no excited state density or geometry optimization) for the computation of a vertical electronic excitation, the default is to do a non-equilibrium process. For this case then, the "fast" component of the polarization "responds" to the change in the solute density from ground to excited state, but the "slow" component did not have time to "respond" so it still comes from the one that was in equilibrium with the solute's ground state density. This is the case of both "Step 2" and "Step 3" in the example shown in the manual. In "Step 2", an energy calculation using TD is performed, thus it defaults to non-equilibrium solvation. The solvation effects on the excited states energies are computed by means of a linear response approach. The absorption energies via the linear response approach only are those reported directly in the output of this "Step 2" job. In "Step 3", a step further is taken and a correction of the linear response excitation energy is performed by solving the "fast" component of the solvent polarization self-consistently with the selected excited state density (the "State-Specific" approach). This is generally an improvement over the excitation energies obtained by linear response alone. Note that since this "State-Specific" approach involves the calculation of the excited state density, the program would default to doing an equilibrium solvation calculation on the excited state. However, the goal of "Step 3" is to compute the vertical excitation energy, so as mentioned above, we would like to use the "slow" component of the solvent polarization from the ground state calculation (in "Step 3", the first part does an equilibrium calculation on the ground state saving the solvent reaction field to the checkpoint file) and solving self-consistently the "fast" component with the excited state density (the second part of "Step 3" reads the reaction field from the checkpoint file, the one from the ground state calculation, keeps the "slow" component as is, and solves the "fast" component self-consistently with the excited state density). The absorption energy via the "State-Specific" approach is the energy difference between the excited state energy after all PCM corrections from the non-equilibrium calculation in "Step 3" and the ground state energy resulting from the equilibrium process (either first part of "Step 3" or final, optimized geometry, energy from "Step 1", the two ground state energies should be the same). In a TD geometry optimization of an excited state, since one is looking for the equilibrium geometry, the default is to do equilibrium solvation, so the two components, "slow" and "fast", of the polarization are in equilibrium with the solute's excited state density. All "Step 4", "Step 5" and "Step 6" use equilibrium solvation for the selected excited state. The emission energy (vertical energy of the excited to ground state transition) by means of a linear response approach can be found in the output of "Step 4". For the final (optimized) geometry in "Step 4", the "excitation energy" shown in this output would be equal to the emission energy since it is the result of an equilibrium calculation on the selected excited state. "Step 6" and "Step 7" are analogous to the two parts of "Step 3" but this time for the opposite transition (excited to ground states). Thus, "Step 6" is analogous to the first part of "Step 3", it is an equilibrium calculation on the origin state (now the excited state) in which both "slow" and "fast" components of solvent polarization are solved self-consistently with the excited state density (this calculation can be regarded as a correction of the excited state energy beyond the linear response approach, which was done in "Step 4"  . Now, this "Step 6" saves the solvent reaction field to the checkpoint file. "Step 7" reads this information from the file (just like the second part of "Step 3" and performs a non-equilibrium calculation of the ground state energy, using the "slow" component of the solvent polarization from the excited state calculation ("Step 6" and only doing the "fast" component of the solvent polarization self-consistent with the ground state density. The emission energy via the "State-Specific" approach would be the energy difference between the excited state energy after all PCM corrections from "Step 6" and the ground state energy resulting from the non-equilibrium process in "Step 7".These are a few relevant references that deal with different aspects of PCM, equilibrium vs. non-equilibrium solvation, solvation effects in electronic transition calculations: - A review that covers many of the aspects of continuum solvation models (including a section in non-equilibrium): J. Tomasi, B. Mennucci and R. Cammi, "Quantum Mechanical Continuum Solvation Models," Chem. Rev., 105 (2005) 2999-3093. - A couple of early references dealing with solvation in electronic transitions: M. Cossi and V. Barone, “Solvent effect on vertical electronic transitions by the polarizable continuum model,” J. Chem. Phys., 112 (2000) 2427-35. M. Cossi and V. Barone, “Time-dependent density functional theory for molecules in liquid solutions,” J. Chem. Phys., 115 (2001) 4708-17. - A couple of references related to the State-Specific (or external iteration) treatment of solvation in electronic transitions: R. Improta, V. Barone, G. Scalmani, and M. J. Frisch, “A state-specific polarizable continuum model time dependent density functional method for excited state calculations in solution,” J. Chem. Phys., 125 (2006) 054103. R. Improta, G. Scalmani, M. J. Frisch, and V. Barone, “Toward effective and reliable fluorescence energies in solution by a new State Specific Polarizable Continuum Model Time Dependent Density Functional Theory Approach,” J. Chem. Phys., 127 (2007) 074504. 2. For seven example steps in the guide of scrf keywords, each step need read the data from chk file from the second step on. But the guide of those steps does not specify which step's chk should be used in current step. For example, the third step read the chk of first step or second step? Each step use the chk of former step? （這一段也是我問(wèn)的問(wèn)題） - Step 2: requires a copy of the checkpoint file from Step 1. - Step 3: will work with either a copy of checkpoint file from Step 1 or from Step 2. - Step 4: requires a copy of the checkpoint file from Step 2. - Step 5: requires a copy of the checkpoint file from Step 4. - Step 6: requires a copy of the checkpoint file from Step 4. - Step 7: requires a copy of the checkpoint file from Step 6. Please, note that there was a problem with G09 rev. A.02 because the "StateSpecific" option was connected incorrectly, but this has been fixed in our latest release, G09 rev. B.01. So, if you are using G09 rev. A.02, you would need to use "ExternalIteration" or "SelfConsistent" instead of the "StateSpecific" (shown in the example in the manual). If you are using G09 rev. B.01, then you can use "ExternalIteration", "SelfConsistent" or "StateSpecific", as the three are synonyms. If you follow the whole example given in the manual with G09 rev. A.02 but using "ExternalIteration", you should obtain the following results on each one of the 7 proposed jobs in the example: 1) Energy of ground state optimized geometry in solution: SCF Done: E(RB3LYP) = -153.851761719 A.U. after 1 cycles 2) Vertical excitation (absorption) to first excited state from Linear Response calculation in solution (using ground state optimized geometry from job 1): Excited State 1: Singlet-A" 4.3767 eV 283.28 nm f=0.0000 =0.000 3) Vertical excitation (absorption) to first excited state from a non-equilibrium state-specific calculation (using "ExternalIteration" in solution (using ground state optimized geometry from job 1):SC-PCM: The total energy including the PCM contribution is -153.687679830 a.u. Taking the difference of this energy minus the ground state energy (job 1), the excitation energy is 277.69 nm. 4) Geometry optimization of first excited state in solution (Final energy): Excited State 1: Singlet-A 3.2076 eV 386.53 nm f=0.0014 =0.000 5) Frequency calculation to confirm that the geometry optimized in job 4 is a stationary point, and a minimum in the potential energy surface of the excited state. 6) Energy of first excited state from a equilibrium state-specific calculation (using "ExternalIteration" in solution (using optimized geometry from job 4):SC-PCM: The total energy including the PCM contribution is -153.707151097 a.u. 7) Energy of ground state from a non-equilibrium calculation in solution (using first excited state optimized geometry from job 4 and its reaction field from job 6): SCF Done: E(RB3LYP) = -153.822029934 A.U. after 10 cycles The difference between these last two energies gives the vertical emission energy, which, in this case, is 396.62 nm. |
» 本帖已獲得的紅花（最新10朵）
2樓2014-08-17 22:15:20
小范范1989
木蟲 (著名寫手)
- 應(yīng)助: 123 (高中生)
- 金幣: 3821.3
- 散金: 968
- 紅花: 243
- 帖子: 1979
- 在線: 452.8小時(shí)
- 蟲號(hào): 2013403
- 注冊(cè): 2012-09-19
- 性別: GG
- 專業(yè): 原子和分子物理
★
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
送紅花一朵
小木蟲: 金幣+0.5, 給個(gè)紅包，謝謝回帖
送紅花一朵 |
這么好的帖子為什么沒人頂？我也是計(jì)算熒光的時(shí)候，收到hairan的指點(diǎn)，到了這個(gè)帖子上，我感覺我能理解這些步驟了。比原來(lái)好多了。 再就是，我感覺要定下心來(lái)閱讀以下這個(gè)英文的說(shuō)明，還有后面的解釋，非常非常的好。只要是你想學(xué)，一定會(huì)學(xué)會(huì)的。 我們實(shí)驗(yàn)室的師姐，我問(wèn)他高斯計(jì)算熒光，他說(shuō)就是守則上的沒看明白，我感覺就是沒用心，只要用心，什么事都能學(xué)會(huì)。 最后，謝謝hairan，希望小木蟲成為我們科研的交流平臺(tái)。 |
5樓2014-09-30 21:48:59
zhuyunfy
新蟲 (小有名氣)
- 應(yīng)助: 2 (幼兒園)
- 金幣: 721
- 帖子: 110
- 在線: 37.2小時(shí)
- 蟲號(hào): 2568635
- 注冊(cè): 2013-07-27
- 專業(yè): 結(jié)構(gòu)化學(xué)
6樓2015-06-21 23:22:02
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 070300化學(xué)279求調(diào)劑 +8 | 哈哈哈^_^ 2026-03-31 | 10/500 |
|
|---|---|---|---|---|
|
|
[考研] 材料科學(xué)與工程求調(diào)劑 +13 | 深V宿舍吧 2026-03-29 | 13/650 |
|
|
[考研] 375求調(diào)劑 +7 | 雨夏整夜 2026-03-29 | 7/350 |
|
|
|
[考研] 一志愿a區(qū)211，085601-307分求調(diào)劑 +10 | 黨嘉豪 2026-03-31 | 23/1150 |
|
|
|
[考研] 304求調(diào)劑 +8 | 素年祭語(yǔ) 2026-03-31 | 11/550 |
|
|
|
[考研] 308求調(diào)劑 +11 | 墨墨漠 2026-03-25 | 11/550 |
|
|
|
[考研] 318求調(diào)劑 +10 | 陳晨79 2026-03-30 | 10/500 |
|
|
|
[考研] 353求調(diào)劑 +3 | 江上楓_26 2026-03-28 | 3/150 |
|
|
|
[考研] 學(xué)碩274求調(diào)劑 +17 | Li李魚 2026-03-26 | 17/850 |
|
|
|
[考研] 269求調(diào)劑 +4 | 我想讀研11 2026-03-31 | 4/200 |
|
|
|
[考研] 083000學(xué)碩274求調(diào)劑 +12 | Li李魚 2026-03-26 | 12/600 |
|
|
|
[考研] 282求調(diào)劑 +3 | 呼吸都是減肥 2026-03-31 | 3/150 |
|
|
|
[考研] 323分 食品與營(yíng)養(yǎng)調(diào)劑 +3 | 嘿ooo 2026-03-31 | 3/150 |
|
|
|
[考研] 085601一志愿中山大學(xué)深圳材料工程330求調(diào)劑 +5 | pipiver 2026-03-30 | 5/250 |
|
|
|
[有機(jī)交流]
 甲基亞磺磺酸鈉和甲基磺酸酯反應(yīng)機(jī)理
 10+3
|
kaobao456 2026-03-29 | 4/200 |
|
|
|
[考研] 324求調(diào)劑 +9 | hanamiko 2026-03-26 | 11/550 |
|
|
|
[考研] 環(huán)境科學(xué)與工程334分求調(diào)劑 +6 | 王一一依依 2026-03-30 | 8/400 |
|
|
|
[考研] 085600，材料與化工321分求調(diào)劑 +10 | 大饞小子 2026-03-28 | 10/500 |
|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
|
[有機(jī)交流]
高溫高壓反應(yīng)求助
10+4
|
chibby 2026-03-25 | 4/200 |
|