24小时热门版块排行榜    

查看: 3140  |  回复: 9

cj4566

木虫 (正式写手)

[求助] gaussian 09溶剂化荧光发射的例子第4步运行出错

我在运行高斯09的溶剂中荧光发射的例子时出错,希望前辈们指点一下啊。按照http://www.gaussian.com/g_tech/g_ur/k_scrf.htm的步骤,在第4步的时候出错,最后出现的信息是:
  Cartesian force constants read from checkpoint file:
  04-ac.chk
  Length of force constants on CHK file is     0 NAT3TT is   231
  Error termination via Lnk1e in /chimie/soft/g09/l103.exe

执行的计算是:
%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
回复此楼

» 收录本帖的淘帖专辑推荐

我的

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)

★ ★
gmy1990: 金币+2 2012-04-09 11:28:35
去掉 “=RCFC”
因为你的前一步没有频率计算。
# B3LYP/6-31+G(d,p) TD=(Read,NStates=6,Root=1) SCRF=(Solvent=Ethanol)
  Geom=Modify Guess=Read Opt

这样应该没问题
2楼2012-04-09 10:24:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
2楼: Originally posted by Illusionist at 2012-04-09 10:24:43:
去掉 “=RCFC”
因为你的前一步没有频率计算。
# B3LYP/6-31+G(d,p) TD=(Read,NStates=6,Root=1) SCRF=(Solvent=Ethanol)
  Geom=Modify Guess=Read Opt

这样应该没问题

太感谢了,我也觉得去掉“=RCFC”就不会报错。但是想不通为什么例子要读那个力常数。还有一个问题,下面的4 1 2 3 10.0
5 1 2 7 -50.0 是改变一下二面角来破坏对称性吧,这里后面的角度是可以随便改的吧,对其他体系也只要稍微改变结构就行?谢谢了!
3楼2012-04-10 11:39:41
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)

引用回帖:
3楼: Originally posted by cj4566 at 2012-04-10 11:39:41:
太感谢了,我也觉得去掉“=RCFC”就不会报错。但是想不通为什么例子要读那个力常数。还有一个问题,下面的4 1 2 3 10.0
5 1 2 7 -50.0 是改变一下二面角来破坏对称性吧,这里后面的角度是可以随便改的吧,对其他 ...

这个我也不太感确定,反正我算的时候就是不改的。
4楼2012-04-10 13:41:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
4楼: Originally posted by Illusionist at 2012-04-10 13:41:31:
这个我也不太感确定,反正我算的时候就是不改的。

那对其他体系呢,需要自己随便改改来破坏对称性吧。
5楼2012-04-11 14:49:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)

引用回帖:
5楼: Originally posted by cj4566 at 2012-04-11 14:49:12:
那对其他体系呢,需要自己随便改改来破坏对称性吧。

我算的时候没改,因为不知道怎么改············我也不是高手····

» 本帖已获得的红花(最新10朵)

6楼2012-04-11 16:29:00
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

送鲜花一朵
引用回帖:
6楼: Originally posted by Illusionist at 2012-04-11 16:29:00:
我算的时候没改,因为不知道怎么改············我也不是高手····

你好,在第四步的时候你没有破坏对称性直接算?也就是不用例子中的Geom=Modify而还是用Geom=Check吧。另外还有一个问题想请教一下,第二步计算垂直激发的目的是找出感兴趣的状态,例子中感兴趣的是root=1,所以在第三步TD(NStates=6,Root=1)。对于自己的体系,我想模拟用紫外最大吸收波长激发的荧光光谱谱图,那么我应该找的是模拟的UV中最大吸收峰对应的激发态吧,也就是振子强度f最大值对应的root?望指点,谢谢了!!!
7楼2012-04-22 15:00:20
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

Illusionist

银虫 (正式写手)

引用回帖:
7楼: Originally posted by cj4566 at 2012-04-22 15:00:20:
你好,在第四步的时候你没有破坏对称性直接算?也就是不用例子中的Geom=Modify而还是用Geom=Check吧。另外还有一个问题想请教一下,第二步计算垂直激发的目的是找出感兴趣的状态,例子中感兴趣的是root=1,所以在 ...

我用的的确是Geom=Check。第二步计算垂直激发的目的是找出感兴趣的状态,例子中感兴趣的是root=1,所以在第三TD(NStates=6,Root=1)。是这样的。
关于振子强度f我也不太了解··不好意思。
8楼2012-04-23 07:28:12
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

引用回帖:
8楼: Originally posted by Illusionist at 2012-04-23 07:28:12:
我用的的确是Geom=Check。第二步计算垂直激发的目的是找出感兴趣的状态,例子中感兴趣的是root=1,所以在第三TD(NStates=6,Root=1)。是这样的。
关于振子强度f我也不太了解··不好意思。

ok,谢谢了!
9楼2012-04-23 09:22:11
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

cj4566

木虫 (正式写手)

关于例子的step4,我跟高斯官方确认了一下,step4应该读取step2的的checkpoint文件,同时去掉“=RCFC”,具体见邮件回复:

Thank you for your inquiry. "Step 4" should read the checkpoint file from "Step 2", since both of them use the linear-response approach to compute the excited state wave functions ("Step 3" follows the "StateSpecific" approach, so it is a little different). Nevertheless, the TD calculation in "Step 2" removed the force constants from the checkpoint file that came from "Step 1". In general, there would not be much harm in just removing the "RCFC" option from "Step 4" since it is a geometry optimization on the excited state but the force constants were from the ground state, and the curvature of the potential energy surfaces for the two states can be different in that region.


If the force constants from the ground state are useful for your calculation on the excited state, then use a copy of "01-ac.chk" for your "Step 4" calculation instead of a copy of "02-ac.chk", in which case, you would then have to remove the "Read" option from the "TD" keyword since "01-ac.chk" does not contain transition densities.
10楼2012-04-24 06:28:53
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 cj4566 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[找工作] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:54 by yCO1Ll7aHtsw
[考博] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +4 ZYXdhzDAy9ZX 2026-09-28 4/200 2026-09-29 16:47 by yCO1Ll7aHtsw
[博后之家] 售SCI-T0P文章,我:8O.5.5.1.O.54,科目齐全,可+急 +3 ZYXdhzDAy9ZX 2026-09-28 3/150 2026-09-29 16:30 by etmYJ6d2rquH
[考博] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:46 by etmYJ6d2rquH
[论文投稿] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +4 tqUTOxClQUMF 2026-09-28 4/200 2026-09-29 15:15 by etmYJ6d2rquH
[博后之家] 售SCI文章,我:8O5.5.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 15:04 by etmYJ6d2rquH
[公派出国] 售SCI一区T0P文章,我:8.O.55.1.O54,科目全,可伽急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-29 14:58 by etmYJ6d2rquH
[硕博家园] 售SCI一区T0P文章,我:8O.55.1.O.5.4,科目齐全,可+急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-29 14:30 by etmYJ6d2rquH
[考研] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 3/150 2026-09-29 14:24 by etmYJ6d2rquH
[找工作] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +6 IDs3scOF0tjC 2026-09-28 6/300 2026-09-29 14:13 by etmYJ6d2rquH
[博后之家] 售SCI一区文章,我:8.O.551.O.5.4,科目全,可伽急 +3 IDs3scOF0tjC 2026-09-28 3/150 2026-09-29 13:56 by etmYJ6d2rquH
[教师之家] 某top大学教授说“能够在市场中兑现的能力才是真能力”无比同意! +8 zju2000 2026-09-26 9/450 2026-09-29 13:47 by ahsoarli
[教师之家] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +4 GDBe8tDZqE8z 2026-09-28 4/200 2026-09-29 09:49 by JzYBbHIXWSrW
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 tqUTOxClQUMF 2026-09-28 3/150 2026-09-28 23:46 by b8fCuTHqckEt
[考研] 售SCI一区文章,我:8.O.55.1.O.54,科目齐全,可伽急 +3 GDBe8tDZqE8z 2026-09-28 3/150 2026-09-28 23:12 by b8fCuTHqckEt
[硕博家园] 售SCI一区T0P文章,我:8.O.55.1.O.54,科目齐全,可+急 +3 CfXuS1rDhLYN 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[公派出国] 售SCI一区文章,我:8O5.5.1.O5.4,科目全,可伽急 +4 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 22:55 by ez6fGg9abYaj
[有机交流] 同一个分子,一条来自文献,一条来自AI——不告诉你答案,你会选哪条? +3 tianxiaoxian 2026-09-28 6/300 2026-09-28 22:44 by tianxiaoxian
[找工作] 售SCI文章,我:8O.5.5.1O.54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 5/250 2026-09-28 22:32 by ez6fGg9abYaj
[考研] 售一区SCI文章T0P,我:8O.551.O54,科目全,可十急 +3 IDs3scOF0tjC 2026-09-28 4/200 2026-09-28 18:51 by 9lS3ad5oOymn
信息提示
请填处理意见