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

24小时热门版块排行榜    

小木虫论坛-学术科研互动平台 » 计算模拟区 » 量子化学 » Gaussian » Gaussian TDDFT 激发态 频率计算 输入文件
91/1返回列表
查看: 4508  |  回复: 8
只看楼主 @他人 存档 新回复提醒 (忽略) 收藏    在APP中查看

shenzhendsh

铁虫 (小有名气)

[求助] Gaussian TDDFT 激发态 频率计算 输入文件

我要计算单重激发态的频率，输入文件如下：
%chk=chk.chk
%nprocshared=8
%mem=3200MB
# freq td=(nstates=30,root=1,singlet) b3lyp  tzvp

title

Cartesians coordination
请问为何一直算不完，而且有很多SCF Done能量？正确的输入命令是什么呢？
对于激发态过渡态的优化，输入命令又是什么呢？
刚接触激发态计算，还望大家不吝赐教！非常感谢！
回复此楼

» 猜你喜欢

» 本主题相关价值贴推荐，对您同样有帮助:
1楼 2013-11-12 15:59:09
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与，应助指数 +1
shenzhendsh: 金币+1 2013-11-12 17:25:10
gmy1990: 金币+3 2013-11-12 17:25:36
一直算不完是什么意思？算了一天？两天？而且你算了30个激发态，如果体系大的话算个十几天都很正常。
SCF的能量不用管它，都会算的，只需关注最后一次循环的输出结果即可。
如果要做激发态过渡态优化，加opt关键词就可以了，但注意如此做的话最好用--Link1--把opt与freq分成两个任务来算，不然很有可能会报错。
赞一下(1人) 回复此楼
2楼2013-11-12 16:21:39
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shenzhendsh

铁虫 (小有名气)
引用回帖:
2楼: Originally posted by 枪下游魂 at 2013-11-12 16:21:39
一直算不完是什么意思？算了一天？两天？而且你算了30个激发态，如果体系大的话算个十几天都很正常。
SCF的能量不用管它，都会算的，只需关注最后一次循环的输出结果即可。
如果要做激发态过渡态优化，加opt关键词 ...

多谢回复！体系含三个V，8个O，一个C，四个H，8个CPU跑了两天了，而且SCF Done的能量几乎不变。
做激发态过渡态优化时，只加opt不行吧，我加了opt（ts,calcfc）,结果报错如下：No analytic 2nd derivatives for this method.
>Error termination via Lnk1e in /home/soft/g09/l1.exe at Tue Nov 12 11:30:43 2013.
优化时应加什么关键词呢？多谢回复！
赞一下 回复此楼
3楼2013-11-12 16:27:56
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

枪下游魂

木虫 (著名写手)
引用回帖:
3楼: Originally posted by shenzhendsh at 2013-11-12 16:27:56
多谢回复！体系含三个V，8个O，一个C，四个H，8个CPU跑了两天了，而且SCF Done的能量几乎不变。
做激发态过渡态优化时，只加opt不行吧，我加了opt（ts,calcfc）,结果报错如下：No analytic 2nd derivatives for t ...

算激发态频率的耗时有时候会比基态频率的多得多，只要高斯还在跑，应该就没问题。
报错意思应是没有解析二阶导数。但我不算过渡态，不知道具体是哪的问题，也许并不是opt本身的问题。这个只能等高人回答了
赞一下(1人) 回复此楼
4楼2013-11-12 16:48:48
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

pwzhou

铁杆木虫 (正式写手)

【答案】应助回帖

★ ★ ★ ★ ★ ★ ★ ★ ★ ★
感谢参与，应助指数 +1
zhou2009: 金币+5 2013-11-14 14:27:24
shenzhendsh: 金币+5 2013-11-17 16:27:42
Gaussian的TDDFT算激发态频率是数值的，所以需要算6N+1次（N是体系的原子数目）TD，所以你会找到6N+1个"SCF Done"，什么时候"SCF Done"的数目到达这个了，什么时候你的激发态频率计算就结束了。2天时间太少，TDDFT算个激发态频率20天到一个月都是正常的，慢慢等吧。没事数数"SCF Done"的数目，大概估算一下需要算多久吧。
赞一下(1人) 回复此楼
5楼2013-11-14 13:53:10
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

shenzhendsh

铁虫 (小有名气)
引用回帖:
5楼: Originally posted by pwzhou at 2013-11-14 13:53:10
Gaussian的TDDFT算激发态频率是数值的，所以需要算6N+1次（N是体系的原子数目）TD，所以你会找到6N+1个"SCF Done"，什么时候"SCF Done"的数目到达这个了，什么时候你的激发态频率计算就结束了。 ...

非常感谢！频率计算确实是那样的。我还想请教一个问题，请问激发态的TS如何优化？ opt(ts,calcfc) 是不行的，因为TD只能数值分析，无法计算力常数。请问opt的命令该如何写呢？谢谢！
赞一下 回复此楼
6楼2013-11-17 16:27:31
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

blueybz

木虫 (正式写手)
引用回帖:
6楼: Originally posted by shenzhendsh at 2013-11-17 16:27:31
非常感谢！频率计算确实是那样的。我还想请教一个问题，请问激发态的TS如何优化？ opt(ts,calcfc) 是不行的，因为TD只能数值分析，无法计算力常数。请问opt的命令该如何写呢？谢谢！...

激发态的TS有这个概念吗
赞一下 回复此楼
7楼2013-11-18 11:19:46
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

小范范1989

木虫 (著名写手)
引用回帖:
5楼: Originally posted by pwzhou at 2013-11-14 13:53:10
Gaussian的TDDFT算激发态频率是数值的，所以需要算6N+1次（N是体系的原子数目）TD，所以你会找到6N+1个"SCF Done"，什么时候"SCF Done"的数目到达这个了，什么时候你的激发态频率计算就结束了。 ...

大神你好，我在这问您个问题，就是我有个分子，大约含有110额原子，之前优化激发态一直出错，使用的方法是b3lyp/6-31g*。假如说，我计算激发态的频率，能使用b3lyp/6-31g*吗？sob给的解说是这样的：Gaussian没有TDDFT解析二阶导数，这是几乎做不动的，耗时我估计是600~1200次激发能计算的耗时。必须要计算频率的话，只能是在CIS下优化并计算频率，因为CIS下有解析二阶导数，速度和TDDFT计算频率耗时完全不在一个数量级。
我想，基态是b3lyp/6-31g*,激发态用CIS能行吗？b3lyp真的做不了我的这个计算？我曾经计算过一个70个分子，很长很长时间，我该怎么办？谢谢pwzhou的指点。谢谢您
赞一下(1人) 回复此楼
It doesn't matter how slow you are, as long as you're determined to get there, you'll get there.
8楼2015-01-14 23:41:04
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

忧凝听雨轩

铁虫 (初入文坛)
你好，请问激发态的频率计算怎么做？为什么我计算老出错
分子的激发态优化后做频率计算
输入文件是
%chk=b.chk
%rwf=b.rwf
%mem=200MW
%nprocshare=12
# freq td=(singlets,nstates=6,root=1) b3lyp/6-31g(d)
计算到下边就停止了
SavETr:  write IOETrn=   770 NScale= 10 NData=  16 NLR=1 LETran=     118.
          IDoAtm=1111111111111111111111111111111111111111
          Differentiating once with respect to electric field.
                with respect to dipole field.
          Differentiating once with respect to nuclear coordinates.
          There are     1 degrees of freedom in the 1st order CPHF.  IDoFFX=0.
LinEq1:  Iter=  0 NonCon=     1 RMS=2.02D-03 Max=2.21D-01
AX will form     1 AO Fock derivatives at one time.
LinEq1:  Iter=  1 NonCon=     1 RMS=4.31D-04 Max=2.25D-02
LinEq1:  Iter=  2 NonCon=     1 RMS=2.91D-04 Max=1.36D-02
LinEq1:  Iter=  3 NonCon=     1 RMS=2.27D-04 Max=1.60D-02
LinEq1:  Iter=  4 NonCon=     1 RMS=1.99D-04 Max=1.26D-02
LinEq1:  Iter=  5 NonCon=     1 RMS=1.04D-04 Max=6.16D-03
LinEq1:  Iter=  6 NonCon=     1 RMS=5.06D-05 Max=1.92D-03
LinEq1:  Iter=  7 NonCon=     1 RMS=2.28D-05 Max=1.17D-03
LinEq1:  Iter=  8 NonCon=     1 RMS=1.49D-05 Max=1.22D-03
LinEq1:  Iter=  9 NonCon=     1 RMS=1.48D-05 Max=1.31D-03
LinEq1:  Iter= 10 NonCon=     1 RMS=7.12D-06 Max=2.68D-04
LinEq1:  Iter= 11 NonCon=     1 RMS=3.05D-06 Max=1.36D-04
LinEq1:  Iter= 12 NonCon=     1 RMS=2.41D-06 Max=1.46D-04
LinEq1:  Iter= 13 NonCon=     1 RMS=1.28D-06 Max=5.58D-05
赞一下 回复此楼
9楼2016-01-19 22:48:51
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 shenzhendsh 的主题更新
91/1返回列表
普通表情 龙 兔 虎 猫
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 289求调剂 +4 yang婷 2026-03-02 4/200 2026-03-02 14:31 by 向上的胖东
[考研] 江苏省农科院招调剂1名 +4 Qwertyuop 2026-03-01 4/200 2026-03-02 14:27 by 升格阿达
[考研] 290分材料工程085601求调剂 数二英一 +3 llx0610 2026-03-02 3/150 2026-03-02 14:15 by yc258
[考研] 0856调剂 +7 刘梦微 2026-02-28 7/350 2026-03-02 14:11 by liyongv
[考研] 292求调剂 +7 yhk_819 2026-02-28 7/350 2026-03-02 12:43 by 无际的草原
[基金申请] 成果系统访问量大，请15分钟后再尝试。由此给您造成的不便，敬请谅解。 +5 xhuama 2026-03-02 5/250 2026-03-02 12:34 by stidwellNK
[考研] 281求调剂 +5 2026计算机_诚心 2026-03-01 8/400 2026-03-02 11:05 by 汪！？！
[考研] 一志愿郑大材料学硕298分，求调剂 +6 wsl111 2026-03-01 6/300 2026-03-02 11:00 by ydudjddnd
[考研] 化工专硕342，一志愿大连理工大学，求调剂 +6 kyf化工 2026-02-28 7/350 2026-03-02 10:56 by 无际的草原
[考研] 275求调剂 +3 L-xin? 2026-03-01 6/300 2026-03-02 10:22 by 热情沙漠
[考研] 0854复试调剂 276 +4 wmm9 2026-03-01 6/300 2026-03-02 09:28 by 热情沙漠
[考研] 材料学调剂 +10 提神豆沙包 2026-02-28 12/600 2026-03-02 09:26 by 李老师！
[考研] 322求调剂 +3 熊境喆 2026-03-01 3/150 2026-03-02 08:44 by houyaoxu
[考研] 高分子化学与物理调剂 +6 好好好1233 2026-02-28 13/650 2026-03-02 07:27 by 好好好1233
[考研] 321求调剂一志愿东北林业大学材料与化工英二数二 +4 虫虫虫虫虫7 2026-03-01 7/350 2026-03-01 16:52 by caszguilin
[基金申请] 刚录用，没有期刊号，但是在线可看的论文可以放为代表作吗 10+3 arang1 2026-03-01 3/150 2026-03-01 16:43 by babero
[考研] 303求调剂 +4 今夏不夏 2026-03-01 4/200 2026-03-01 14:46 by 嘟嘟小浣熊
[考研] 295复试调剂 +3 简木ChuFront 2026-03-01 3/150 2026-03-01 14:27 by zzxw520th
[论文投稿] Optics letters投稿被拒求助 30+3 luckyry 2026-02-26 4/200 2026-03-01 09:06 by babero
[高分子] 求环氧树脂研发1名 +3 孙xc 2026-02-25 11/550 2026-02-28 16:57 by ichall
信息提示
关闭
请填处理意见
关闭