24小时热门版块排行榜    

查看: 4327  |  回复: 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的回帖

忧凝听雨轩

铁虫 (初入文坛)

你好,请问激发态的频率计算怎么做?为什么我计算老出错
分子的激发态优化后做频率计算
输入文件是
%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的回帖
查看全部 9 个回答

枪下游魂

木虫 (著名写手)

【答案】应助回帖

★ ★ ★ ★
感谢参与,应助指数 +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的回帖
信息提示
请填处理意见