24小时热门版块排行榜    

查看: 2684  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿西安交通大学材料工程专业 282分求调剂 +11 枫桥ZL 2026-03-18 13/650 2026-03-22 20:26 by edmund7
[考研] 306求调剂 +6 chuanzhu川烛 2026-03-18 6/300 2026-03-22 19:39 by 公瑾逍遥
[考研] 289求调剂 +7 怀瑾握瑜l 2026-03-20 7/350 2026-03-22 15:57 by ColorlessPI
[考研] 生物学一志愿985,分数349求调剂 +4 zxts12 2026-03-21 7/350 2026-03-22 09:57 by zxts12
[考研] 材料求调剂 +5 @taotao 2026-03-21 5/250 2026-03-21 20:55 by lbsjt
[考研] 材料与化工(0856)304求B区调剂 +3 邱gl 2026-03-20 7/350 2026-03-21 19:05 by 15709483992
[考研] 0805 316求调剂 +3 大雪深藏 2026-03-18 3/150 2026-03-21 18:55 by 学员8dgXkO
[考研] 求调剂 +3 白QF 2026-03-21 3/150 2026-03-21 13:12 by zhukairuo
[考研] 301求调剂 +10 yy要上岸呀 2026-03-17 10/500 2026-03-21 03:14 by JourneyLucky
[考研] 303求调剂 +5 睿08 2026-03-17 7/350 2026-03-21 03:11 by JourneyLucky
[考研] 一志愿武理材料305分求调剂 +6 想上岸的鲤鱼 2026-03-18 7/350 2026-03-21 01:03 by JourneyLucky
[考研] 材料专业求调剂 +6 hanamiko 2026-03-18 6/300 2026-03-21 00:24 by JourneyLucky
[考研] 0817 化学工程 299分求调剂 有科研经历 有二区文章 +22 rare12345 2026-03-18 22/1100 2026-03-20 20:39 by zhukairuo
[考研] 295材料求调剂,一志愿武汉理工085601专硕 +5 Charlieyq 2026-03-19 5/250 2026-03-20 20:35 by JourneyLucky
[考研] 材料学求调剂 +4 Stella_Yao 2026-03-20 4/200 2026-03-20 20:28 by ms629
[考研] 328求调剂,英语六级551,有科研经历 +4 生物工程调剂 2026-03-16 12/600 2026-03-19 11:10 by 生物工程调剂
[考研] 【同济软件】软件(085405)考研求调剂 +3 2026eternal 2026-03-18 3/150 2026-03-18 19:09 by 搏击518
[考研] 0854可跨调剂,一作一项核心论文五项专利,省、国级证书40+数一英一287 +8 小李0854 2026-03-16 8/400 2026-03-18 14:35 by 搏击518
[考研] 308求调剂 +4 是Lupa啊 2026-03-16 4/200 2026-03-17 17:12 by ruiyingmiao
[考研] 一志愿苏州大学材料工程(085601)专硕有科研经历三项国奖两个实用型专利一项省级立项 +6 大火山小火山 2026-03-16 8/400 2026-03-17 15:05 by 无懈可击111
信息提示
请填处理意见