24小时热门版块排行榜    

Znn3bq.jpeg
查看: 2722  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 人工智能320调剂08工类还有机会吗 +12 振—TZ 2026-04-10 12/600 2026-04-12 07:54 by wj165256
[考研] 22408调剂315分 +3 zhuangyan123 2026-04-09 3/150 2026-04-12 00:25 by 蓝云思雨
[考研] 本人女孩 +7 吼吼, 2026-04-10 9/450 2026-04-11 14:45 by ACS Nano——
[考研] 22408 352分求调剂 +4 努力的夏末 2026-04-09 4/200 2026-04-11 10:42 by maddjdld
[考研] 283求调剂 +22 那个噜子 2026-04-09 22/1100 2026-04-11 10:41 by 逆水乘风
[考研] 22408 366分,本科211,一志愿西工大 +4 Rubt 2026-04-09 4/200 2026-04-10 19:51 by chemisry
[考研] 265求调剂 +12 风说她早忘了 2026-04-10 13/650 2026-04-10 18:56 by chemisry
[考研] 284求调剂 +7 让我上岸吧阿西 2026-04-09 7/350 2026-04-09 18:59 by haironglove
[考研] 一志愿厦大生物学332求调剂 +10 池池池池池池 2026-04-08 10/500 2026-04-09 17:10 by 独醉梦孤城
[考研] 化学工程与技术专业一志愿哈工程 291分B区 国家级大创负责人 有一作论文 +13 Emmy~ 2026-04-09 13/650 2026-04-09 14:47 by only周
[考研] 调剂 +12 月@163.com 2026-04-08 12/600 2026-04-09 14:27 by rl1980
[考研] 086004 求调剂 309 +7 Yin DY 2026-04-08 7/350 2026-04-09 13:59 by Delta2012
[考研] 材料工程322 +18 哈哈哈吼吼吼哈 2026-04-07 19/950 2026-04-09 10:44 by cymywx
[考研] 机械工程313分找工科调剂 +3 双一流本科机械 2026-04-08 3/150 2026-04-08 20:41 by 土木硕士招生
[考研] 电子信息346 +4 zuoshaodian 2026-04-08 4/200 2026-04-08 11:54 by zzucheup
[考研] 求调剂 +15 熊二想上岸 2026-04-06 15/750 2026-04-08 04:53 by 无际的草原
[考研] 344求调剂 +11 魏子per 2026-04-07 11/550 2026-04-07 23:01 by JourneyLucky
[考研] 一志愿南科大生物学297分,求调剂推荐 +8 Y-yyusx 2026-04-06 9/450 2026-04-07 19:38 by biomichael
[考研] 22408 318分求调剂 +4 勤奋的小笼包 2026-04-06 6/300 2026-04-07 15:05 by 纸鹤555
[考研] 071000生物学调剂 +7 拉提桃 2026-04-06 7/350 2026-04-06 18:55 by 52305043001
信息提示
请填处理意见