24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 4640  |  回复: 8

shenzhendsh

铁虫 (小有名气)

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

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

title

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

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

已阅   回复此楼   关注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分成两个任务来算,不然很有可能会报错。
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本身的问题。这个只能等高人回答了
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"的数目,大概估算一下需要算多久吧。
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的指点。谢谢您
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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 一志愿沪9,326求生物学调剂 +7 刘墨墨 2026-04-13 7/350 2026-04-13 17:23 by ss592599143
[考研] 材料085601调剂 +30 何润采123 2026-04-10 32/1600 2026-04-13 09:51 by szhize
[考研] 一志愿中国科学院上海有机所,有机化学356分找调剂 +12 Nadiums 2026-04-09 12/600 2026-04-13 08:50 by lhj2009
[考研] 339求调剂 +8 hanwudada 2026-04-11 9/450 2026-04-12 15:36 by laoshidan
[考研] 307求调剂 +10 tzq94092 2026-04-10 10/500 2026-04-12 08:18 by wise999
[考研] 一志愿西北工业大学289 085602 +33 yang婷 2026-04-10 34/1700 2026-04-12 08:11 by Art1977
[考研] 280求调剂 +13 wzzz王 2026-04-09 13/650 2026-04-12 00:31 by 勇攀高峰0126
[考研] 284求调剂 +11 archer.. 2026-04-09 12/600 2026-04-11 20:23 by 蓝云思雨
[考研] 调剂 +5 文道星台 2026-04-11 5/250 2026-04-11 15:01 by 凯凯要变帅
[考研] 一志愿211,化学310分,本科重点双非,求调剂 +23 努力奋斗112 2026-04-08 23/1150 2026-04-10 23:29 by 314126402
[考研] 263能源动力专硕求调剂 +4 加大号饭盒袋 2026-04-10 4/200 2026-04-10 20:52 by gong120082
[考研] 材料专业344求调剂 +16 hualkop 2026-04-10 21/1050 2026-04-10 17:28 by laoshidan
[考研] 求调剂 材料与工程 324分 专硕 +19 翩翩一书生 2026-04-10 21/1050 2026-04-10 11:41 by wp06
[考研] 0703化学求调剂 +21 不知名的小卅 2026-04-08 21/1050 2026-04-09 18:55 by l_paradox
[考研] 一志愿中科院105500专业总分315求调剂 +6 lallalh 2026-04-09 7/350 2026-04-09 17:51 by lallalh
[考研] 0860004 求调剂 309分 +6 Yin DY 2026-04-09 6/300 2026-04-09 10:19 by 啊李999
[考研] 286求调剂 +19 Faune 2026-04-08 20/1000 2026-04-09 08:36 by 哦哦123
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 一志愿西南090202求调剂 +4 在线求有学上 2026-04-07 4/200 2026-04-07 19:47 by biomichael
[考研] 333求调剂 +6 合乘杨习夕 2026-04-06 6/300 2026-04-07 09:44 by 猪会飞
信息提示
请填处理意见