24小时热门版块排行榜    

查看: 2685  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 北科281学硕材料求调剂 +6 tcxiaoxx 2026-03-20 6/300 2026-03-22 20:23 by edmund7
[考研] 招08考数学 +4 laoshidan 2026-03-20 8/400 2026-03-22 19:56 by 小皮蛋酱
[考研] 287求调剂 +8 晨昏线与星海 2026-03-19 9/450 2026-03-22 17:01 by i_cooler
[考研] 寻找调剂 +4 倔强芒? 2026-03-21 4/200 2026-03-22 16:14 by 木托莫露露
[考研] 求调剂院校信息 +6 CX 330 2026-03-21 6/300 2026-03-22 15:25 by 无懈可击111
[考研] 一志愿中南大学化学学硕0703总分337求调剂 +3 niko- 2026-03-22 3/150 2026-03-22 15:15 by 杨杨杨紫
[考研] 354求调剂 +7 Tyoumou 2026-03-18 10/500 2026-03-22 11:11 by 人来盛
[考研] 初试 317 +7 半拉月丙 2026-03-20 7/350 2026-03-21 22:26 by peike
[考研] 【考研调剂】化学专业 281分,一志愿四川大学,诚心求调剂 +11 吃吃吃才有意义 2026-03-19 11/550 2026-03-21 18:23 by 学员8dgXkO
[考研] 297求调剂 +11 戏精丹丹丹 2026-03-17 12/600 2026-03-21 17:47 by ColorlessPI
[考研] 302求调剂 +12 呼呼呼。。。。 2026-03-17 12/600 2026-03-21 17:29 by ColorlessPI
[考研] 299求调剂 +4 某某某某位 2026-03-21 4/200 2026-03-21 16:30 by barlinike
[考研] 22 350 本科985求调剂,求老登收留 +3 李轶男003 2026-03-20 3/150 2026-03-21 13:28 by 搏击518
[考研] 一志愿华南师大 070300(化学)304分求调剂 +3 0703武芊慧雪304 2026-03-18 3/150 2026-03-21 00:48 by JourneyLucky
[考研] 一志愿南昌大学,327分,材料与化工085600 +9 Ncdx123456 2026-03-19 9/450 2026-03-20 23:41 by lovewei0727
[考研] 一志愿 西北大学 ,070300化学学硕,总分287,双非一本,求调剂。 +4 晨昏线与星海 2026-03-19 4/200 2026-03-20 22:15 by JourneyLucky
[考研] 08工学调剂 +5 用户573181 2026-03-20 5/250 2026-03-20 15:47 by xia_2003
[考研] 298-一志愿中国农业大学-求调剂 +9 手机用户 2026-03-17 9/450 2026-03-20 14:24 by 无懈可击111
[考研] 材料与化工求调剂 +7 为学666 2026-03-16 7/350 2026-03-19 14:48 by 尽舜尧1
[考研] 312求调剂 +8 陌宸希 2026-03-16 9/450 2026-03-18 12:39 by Linda Hu
信息提示
请填处理意见