应《网络安全法》要求，自2017年10月1日起，未进行实名认证将不得使用互联网跟帖服务。为保障您的帐号能够正常使用，请尽快对帐号进行手机号验证，感谢您的理解与支持！

24小时热门版块排行榜    

CyRhmU.jpeg
小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » gaussian 09溶剂化荧光发射的例子第4步运行出错
南方科技大学公共卫生及应急管理学院2026级博士研究生招生报考通知（长期有效）
101/1返回列表
查看: 2485  |  回复: 9
只看楼主@他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

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
回复此楼

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

我的

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼2012-04-08 18:44:49
已阅   回复此楼   关注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朵）

cj4566
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 的主题更新
101/1返回列表
普通表情 龙 兔 虎 猫
信息提示
关闭
请填处理意见
关闭