| 查看: 5611 | 回复: 25 | |||||
[交流]
【求助】G09荧光光谱例子 运行出错
|
|||||
|
用gaussian09计算乙醛分子在乙醇溶液中的荧光光谱。根据手册，一共分成7个步骤，http://www.gaussian.com/g_tech/g_ur/k_scrf.htm 但是计算到第六步出了问题。 输入文件如下： %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 出错信息： PCMU: Requested matrix not present on file. Error termination via Lnk1e in /public/g09/l601.exe 请教大家，如何解决这个问题呢？ 谢谢！ |
回复此楼» 收录本帖的淘帖专辑推荐
 高斯及相关软件的使用 | Gaussian模拟光谱 |
» 猜你喜欢
美国杰克逊州立大学周剑歌教授的化学计算与AI实验室诚聘博士后研究员
已经有4人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有3人回复
-
物理化学论文润色/翻译怎么收费?
已经有170人回复
-
面上再次挂了，太难了，躺也躺不了，倦也卷不过，小学校之殇！
已经有37人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
Tc数据拟合
已经有1人回复
-
韩国釜山大学储能电池与功能材料实验室诚招博后（电池/电化学/器件制备/MOF合成）
已经有187人回复
-
祈福！
已经有79人回复
-
散金求好运
已经有99人回复
-
博士生招生, 材料模拟方向，University of South Carolina
已经有0人回复
-
Postdoc position, materials simulations/machine learning, U. of South Carolina
已经有0人回复
» 本主题相关商家推荐: (我也要在这里推广)
» 本主题相关价值贴推荐，对您同样有帮助:
- 原子荧光光谱测重金属
已经有10人回复
- 普通分光光度计、荧光分光光度计和荧光光谱仪的区别
已经有4人回复
- gaussian 计算得到的荧光光谱选择其最大发射光谱是看f最大对应的吗？
已经有14人回复
- G09计算拉曼光谱出错
已经有7人回复
- 【求助】按照G09说明书计算荧光出错。出在第二步上。
已经有13人回复
- 【求助】荧光光谱的问题
已经有6人回复
- 【求助】G09进行荧光光谱计算的例子 出错
已经有5人回复
- 【求助】g03如何计算吸收光谱和荧光光谱的特征峰，命令如何编写
已经有20人回复
- 【求助】计算荧光光谱的疑惑，求教！
已经有18人回复
- 【求助】关于荧光光谱上的倍频峰
已经有17人回复
- 荧光光谱，磷光光谱，光致发光光谱三者的异同。有效期至2010年6月20日
已经有8人回复
- 【求助】荧光光谱峰是什么样的？
已经有21人回复
» 抢金币啦！回帖就可以得到:
-
告别流水线爽文！108万字完结仙侠《女娲石之灵石奇缘（仙界篇）》值得细品
+1/4110
-
散金祈福
+3/900
-
每年这时候我都再仔细读一遍《范进中举》
+1/486
-
华西生物制药研究院招聘科研助理
+2/420
-
散金求好运
+1/203
-
上海交通大学张航课题组招聘博士后（电化学储能）
+1/93
-
2026年财务管理、经济与信息技术国际会议 (FMEIT 2026)
+1/92
-
山东征女友，坐标济南
+1/91
-
哈尔滨工业大学（校本部）电子与信息工程学院2027级博士招生
+1/68
-
祈祷老公中国基
+1/45
-
武汉，诚征男友，不异地，年龄大了
+1/36
-
坐标北京诚征男友
+1/17
-
pcold 原核表达载体
+1/17
-
新审查指南施行半年，AI方向的专利申请是不是更看"技术味"了？
+1/11
-
2027环境工程博士申请
+1/9
-
中国科学院生态环境研究中心环境工程相关方向招收保送硕士和联合培养博士和硕士研究生
+1/4
-
清华大学环境课题组招聘环境工作专业客座生研究生2-3名（应用导向）
+1/2
-
北理工集成电路杰青团队 | 诚招科助理
+1/2
-
【产品测评分享】体外转录产量不理想、dsRNA 副产物高？T7 RNA 聚合酶实测体验
+1/1
-
清华大学环境课题组招聘环境工作专业客座生研究生2-3名（应用导向）
+1/1
2楼2011-01-24 14:30:26
3楼2011-01-24 16:18:28
4楼2011-01-25 09:24:00
★
tandz(金币+1):谢谢参与
zhou2009: 2011-01-26 16:06:12
tandz(金币+1):谢谢参与
zhou2009: 2011-01-26 16:06:12
|
请将StateSpecific关键词改为ExternalIteration，这个是Gaussian 09的一个BUG，我专门写信问过Gaussian公司，这个BUG在B01中已经修正，但是在A01和A02中都存在，所以你需要把这个例子中所有的StateSpecific改为ExternalIteration，这两个关键词作用是一样的，在B01中使用这两个关键词的结果是一致的，但是在A01中和A02中使用StateSpecific，即便能得到结果，比如这个例子中计算吸收光谱那一步，其结果也是错的，它仅仅算了真正的StateSpecific中的第一步而已。 Gaussian公司的回复邮件： 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 |
5楼2011-01-25 10:53:35
7楼2011-01-25 15:01:39
paperhunter
荣誉版主 (文学泰斗)
-
专家经验: +14 - 应助: 14603 (教授)
- 贵宾: 11.659
- 金币: 345467
- 帖子: 127554
- 在线: 8880.5小时
- 虫号: 971012
8楼2011-01-25 15:18:40
11楼2011-01-25 15:35:05
12楼2011-01-25 15:56:57
|
您能解释一下那个例子的第6步和第7步吗? 第6步保存non-equilibrium溶剂效应，第7步的B3lyp只是计算能量的,不能像td计算那样给出发射波长,那应该怎样得到发射波长呢? 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 |
14楼2011-01-26 14:56:41
★
tandz(金币+1):谢谢参与
tandz(金币+1):谢谢参与
|
下面是Gaussian公司关于7步计算中的能量的一个计算的说明，包含了你的问题的答案，供参考： 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, |
16楼2011-01-26 22:24:34
17楼2011-01-26 22:28:24
18楼2011-01-26 22:33:55
19楼2011-01-27 14:17:07
20楼2011-03-24 17:58:20
21楼2011-03-24 19:14:36
22楼2011-03-24 19:30:21
23楼2011-03-24 19:31:07
24楼2011-03-25 11:21:29
25楼2011-03-25 11:57:40
26楼2011-04-26 09:55:16
简单回复
zyxme6楼
2011-01-25 11:45
tandz(金币+1):谢谢参与
半熟的小苹果9楼
2011-01-25 15:19
tandz(金币+1):谢谢参与
半熟的小苹果10楼
2011-01-25 15:19
tandz(金币+1):谢谢参与
zyxme13楼
2011-01-25 16:13
tandz(金币+1):谢谢参与
sduqianqia15楼
2011-01-26 15:22
tandz(金币+1):谢谢参与
