24小时热门版块排行榜    

Znn3bq.jpeg
汕头大学海洋科学接受调剂
查看: 4641  |  回复: 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 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[考研] 本科郑州大学,一志愿华东师范大学282求调剂 +28 熊哥xtk 2026-04-07 33/1650 2026-04-13 17:48 by 邹gv
[考研] 245求调剂 +3 冰糖橘?汽水 2026-04-13 7/350 2026-04-13 16:42 by 冰糖橘?汽水
[考研] 0856专硕求调剂 希望是a区院校 +23 好好休息好不好 2026-04-09 26/1300 2026-04-13 14:22 by 张zhihao
[考研] 290求调剂 +18 柯淮然 2026-04-12 20/1000 2026-04-13 12:56 by cyh—315
[考研] 一志愿华工085600 331分 +7 天下ww 2026-04-09 7/350 2026-04-13 09:01 by lhj2009
[考研] 346分,工科0854求调剂,专硕 +6 moser233 2026-04-12 7/350 2026-04-12 22:11 by fqwang
[考研] 一志愿华中农微生物,288分,三年实验经历 +11 代fish 2026-04-09 11/550 2026-04-12 10:21 by Hayaay
[考研] 一志愿西北工业大学289 085602 +33 yang婷 2026-04-10 34/1700 2026-04-12 08:11 by Art1977
[考研] 332求调剂 +14 蕉蕉123 2026-04-10 14/700 2026-04-12 00:27 by 蓝云思雨
[考研] 359求调剂 +5 胃痉挛累了 2026-04-11 5/250 2026-04-11 19:55 by lbsjt
[考研] 农业管理302分求调剂 +3 xuening1 2026-04-10 3/150 2026-04-11 10:18 by zhq0425
[考研] 求调剂 +13 雪逢冬 2026-04-10 13/650 2026-04-11 09:58 by 猪会飞
[考研] 282,求调剂 +12 jggshjkkm 2026-04-09 14/700 2026-04-11 09:39 by 猪会飞
[考研] 314求调剂 +18 xhhdjdjsjks 2026-04-09 19/950 2026-04-10 18:53 by HPUCZ
[考研] 材料专业344求调剂 +16 hualkop 2026-04-10 21/1050 2026-04-10 17:28 by laoshidan
[考研] 282,电气工程专业,求调剂,不挑专业 +9 jggshjkkm 2026-04-10 9/450 2026-04-10 14:55 by 逆水乘风
[考研] 二次调剂求老师收留 +3 笑笑袁 2026-04-08 3/150 2026-04-08 23:50 by 醉在风里
[考研] 一志愿华南师范大学0702物理学305调剂 +4 念常安 2026-04-07 6/300 2026-04-08 22:53 by bljnqdcc
[考研] 331求调剂 +5 luoxin0706. 2026-04-08 5/250 2026-04-08 22:15 by zhouyuwinner
[考研] 一志愿西南090202求调剂 +4 在线求有学上 2026-04-07 4/200 2026-04-07 19:47 by biomichael
信息提示
请填处理意见