| 24小時(shí)熱門版塊排行榜 |
| 26 | 1/1 | 返回列表 |
| 查看: 4945 | 回復(fù): 25 | |||||
[交流]
【求助】G09熒光光譜例子 運(yùn)行出錯(cuò)
|
|||||
|
用gaussian09計(jì)算乙醛分子在乙醇溶液中的熒光光譜。根據(jù)手冊(cè),一共分成7個(gè)步驟,http://www.gaussian.com/g_tech/g_ur/k_scrf.htm 但是計(jì)算到第六步出了問(wèn)題。 輸入文件如下: %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 出錯(cuò)信息: PCMU: Requested matrix not present on file. Error termination via Lnk1e in /public/g09/l601.exe 請(qǐng)教大家,如何解決這個(gè)問(wèn)題呢? 謝謝! |
高斯及相關(guān)軟件的使用 | Gaussian模擬光譜 |
» 搶金幣啦!回帖就可以得到:
+1/279
+1/81
+2/42
+1/36
+1/35
+1/33
+5/20
+1/14
+1/11
+1/7
+1/7
+1/6
+1/5
+1/5
+1/3
+1/3
+1/2
+1/2
+1/1
+1/1
木蟲(chóng) (正式寫手)
木蟲(chóng) (正式寫手)
鐵桿木蟲(chóng) (正式寫手)
|
請(qǐng)將StateSpecific關(guān)鍵詞改為ExternalIteration,這個(gè)是Gaussian 09的一個(gè)BUG,我專門寫信問(wèn)過(guò)Gaussian公司,這個(gè)BUG在B01中已經(jīng)修正,但是在A01和A02中都存在,所以你需要把這個(gè)例子中所有的StateSpecific改為ExternalIteration,這兩個(gè)關(guān)鍵詞作用是一樣的,在B01中使用這兩個(gè)關(guān)鍵詞的結(jié)果是一致的,但是在A01中和A02中使用StateSpecific,即便能得到結(jié)果,比如這個(gè)例子中計(jì)算吸收光譜那一步,其結(jié)果也是錯(cuò)的,它僅僅算了真正的StateSpecific中的第一步而已。 Gaussian公司的回復(fù)郵件: Thank you for your inquiry. The problem with the example shown in the manual is related to the "StateSpecific" option. Unfortunately, the "StateSpecific" option was connected incorrectly in the current version of G09, but this will be fixed in the next minor revision. As a workaround, you can use the "ExternalIteration" option instead of the "StateSpecific" option (for both the state specific vertical excitation and the emission state specific solvation). The "ExternalIteration" option is a generalization of the "StateSpecific" option. So, please, try the example in the manual again with the only change of removing "StateSpecific" and adding "ExternalIteration" instead. For example, %nproc=16 %mem=2000mb %chk=06-ac # B3LYP/6-31+G(d,p) TD=(Read,NStates=6,Root=1) SCRF=(Solvent=Ethanol,ExternalIteration,Read) Geom=Check Guess=Read Acetaldehyde emission state-specific solvation at first excited state optimized geometry 0 1 NonEq=write |
榮譽(yù)版主 (文學(xué)泰斗)
![]() |
專家經(jīng)驗(yàn): +14 |
至尊木蟲(chóng) (著名寫手)
木蟲(chóng) (著名寫手)
木蟲(chóng) (正式寫手)
|
您能解釋一下那個(gè)例子的第6步和第7步嗎? 第6步保存non-equilibrium溶劑效應(yīng),第7步的B3lyp只是計(jì)算能量的,不能像td計(jì)算那樣給出發(fā)射波長(zhǎng),那應(yīng)該怎樣得到發(fā)射波長(zhǎng)呢? 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 |
鐵桿木蟲(chóng) (正式寫手)
|
下面是Gaussian公司關(guān)于7步計(jì)算中的能量的一個(gè)計(jì)算的說(shuō)明,包含了你的問(wèn)題的答案,供參考: Dear Dr. Zhou, As I was saying in my previous email, the "StateSpecific" option in G09 rev. A.02 is connected incorrectly and it is not producing the results that it is supposed to produce, so please do not use the "StateSpecific" option in G09 rev. A.02. For the absorption/emission example in solution, you need to use the "ExternalIteration" option instead, which does produce the correct results. If you follow the whole example given in the manual but correcting the inputs to have "ExternalIteration" instead of "StateSpecific", 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 ... Total Energy, E(TD-HF/TD-KS) = -153.690919717 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 ... Total Energy, E(TD-HF/TD-KS) = -153.705918841 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. Please, let me know if you have any further questions, |
鐵桿木蟲(chóng) (正式寫手)
木蟲(chóng) (小有名氣)
鐵桿木蟲(chóng) (著名寫手)
鐵桿木蟲(chóng) (正式寫手)
鐵桿木蟲(chóng) (正式寫手)
鐵桿木蟲(chóng) (正式寫手)
木蟲(chóng) (正式寫手)




| 26 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 280求調(diào)劑 +3 | 咕嚕曉曉 2026-03-18 | 3/150 |
|
|---|---|---|---|---|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
[考研] 311求調(diào)劑 +10 | 冬十三 2026-03-15 | 11/550 |
|
|
[考研] 281求調(diào)劑(0805) +3 | 煙汐憶海 2026-03-16 | 9/450 |
|
|
[考研] 341求調(diào)劑 +5 | 搗蛋豬豬 2026-03-11 | 7/350 |
|
|
[考研] 328求調(diào)劑,英語(yǔ)六級(jí)551,有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-16 | 8/400 |
|
|
[考研] 材料專碩326求調(diào)劑 +6 | 墨煜姒莘 2026-03-15 | 7/350 |
|
|
[考研] 275求調(diào)劑 +4 | 太陽(yáng)花天天開(kāi)心 2026-03-16 | 4/200 |
|
|
[考研] 267一志愿南京工業(yè)大學(xué)0817化工求調(diào)劑 +6 | SUICHILD 2026-03-12 | 6/300 |
|
|
[考研] 274求調(diào)劑 +5 | 時(shí)間點(diǎn) 2026-03-13 | 5/250 |
|
|
[基金申請(qǐng)] 國(guó)自科面上基金字體 +6 | iwuli 2026-03-12 | 7/350 |
|
|
[考研] 304求調(diào)劑 +5 | 素年祭語(yǔ) 2026-03-15 | 5/250 |
|
|
[考研] 283求調(diào)劑 +10 | 小樓。 2026-03-12 | 14/700 |
|
|
[考研] 070305求調(diào)劑 +3 | mlpqaz03 2026-03-14 | 4/200 |
|
|
[考研] 080500,材料學(xué)碩302分求調(diào)劑學(xué)校 +4 | 初識(shí)可樂(lè) 2026-03-14 | 5/250 |
|
|
[考研] 中科大材料與化工319求調(diào)劑 +3 | 孟鑫材料 2026-03-14 | 3/150 |
|
|
[考研] 297求調(diào)劑 +4 | 學(xué)海漂泊 2026-03-13 | 4/200 |
|
|
[考研] 290求調(diào)劑 +9 | ADT 2026-03-11 | 9/450 |
|
|
[考研] 311求調(diào)劑 +3 | 冬十三 2026-03-13 | 3/150 |
|
|
[考研] 材料與化工085600調(diào)劑求老師收留 +9 | jiaanl 2026-03-11 | 9/450 |
|