请选择 进入手机版 | 继续访问电脑版

计算化学公社

 找回密码
 现在注册!
查看: 826|回复: 4

[Gaussian/gview] gaussian tddft计算与单点能计算求助

[复制链接]

35

帖子

0

威望

227

eV
积分
262

Level 3 能力者

发表于 2018-6-4 10:12:20 | 显示全部楼层 |阅读模式
首先是tddft的问题,我的输入文件如下(这是log文件里的输入文件部分)
%chk=d2-c6h4cl2-bp86
%nproc=8
Will use up to    8 processors via shared memory.
%mem=16gb
----------------------------------------------------------------------
# bp86/6-311+g* SCRF=(solvent=o-DiChloroBenzene) td(singlets,nstates=2
00) Integral(Grid=Ultrafine,Acc2e=12)

d2-c6h4cl2-bp86
---------------
Symbolic Z-matrix:
Charge =  0 Multiplicity = 1
C                    -0.66791  -6.84038  -6.8591
C                    -3.07753  -6.29051  -5.46838
C                    -1.14854  -5.53167  -6.77118
C                    -1.39908  -7.84232  -6.27723
C                    -2.59879  -7.56906  -5.58664
C                    -2.35816  -5.25483  -6.06994
C                    -0.60105  -4.27378  -7.19604
C                    -2.53836  -3.82968  -6.0711
N                    -1.53094  -3.34703  -6.84527
H                     0.24741  -7.04664  -7.37539
H                    -1.05719  -8.85562  -6.34045
H                    -3.13764  -8.38162  -5.14263

……
但是算了有10天了,还是没有算出来,
以下是现在的输出信息
New state     65 was old state     66
New state     66 was old state     65
New state    134 was old state    135
New state    135 was old state    134
Iteration     4 Dimension  1949 NMult  1795 NNew    154
New state    126 was old state    127
New state    127 was old state    126
Iteration     5 Dimension  2106 NMult  1949 NNew    157
New state    133 was old state    134
New state    134 was old state    133
Iteration     6 Dimension  2249 NMult  2106 NNew    143
New state    109 was old state    110
New state    110 was old state    109
Iteration     7 Dimension  2382 NMult  2249 NNew    133
New state     96 was old state     97
New state     97 was old state     96
Iteration     8 Dimension  2503 NMult  2382 NNew    121
Iteration     9 Dimension  2613 NMult  2503 NNew    110
Iteration    10 Dimension  2712 NMult  2613 NNew     99
Iteration    11 Dimension  2796 NMult  2712 NNew     84
Iteration    12 Dimension  2871 NMult  2796 NNew     75
Iteration    13 Dimension  2932 NMult  2871 NNew     61

请问这正常吗?如果不正常请问应该怎么做呢

第二个问题 单点能计算,我用b3lyp/6-31g*结构优化,用wb97x/6-311+g*算tddft,这些都是在Pcm条件下计算的,计算结果是正常的,但是我用wb97x/6-311+g*算单点能时,也出现了提交了计算很久但是输出信息很少的问题,以下是我的输入文件与输出文件片段
输入文件
%nproc=8
%mem=16GB
# wb97x/6-311+g* iop(3/107=0120000000,3/108=0120000000) scf=tight nosymm scf=(maxcycle=1000)

SubPc-ADP-C60-b-v2

0 1
C    -9.606843     0.864127     4.013825
C    -8.160326     0.764237     3.926137
C    -7.454728     1.598095     3.056540
C    -8.168225     2.564706     2.237198
C    -9.558891     2.657211     2.321500
C   -10.291099     1.793479     3.226907
C   -10.140284    -0.477260     4.180657
C    -9.022668    -1.406109     4.194830
C    -7.799328    -0.639835     4.038450
C    -6.746722    -1.153088     3.278028
C    -6.361909     1.061698     2.263898
C    -7.514605     2.622818     0.942856
C    -8.276899     2.755341    -0.222504
C    -9.722076     2.854705    -0.131645
C   -10.352105     2.813436     1.112571
C   -11.580571     2.053358     1.277483
C   -11.541190     1.419876     2.581053
……

输出文件
Rotational constants (GHZ):      0.0115892      0.0036167      0.0030580
Standard basis: 6-311+G(d) (5D, 7F)
  3272 basis functions,  5358 primitive gaussians,  3415 cartesian basis functions
   460 alpha electrons      460 beta electrons
       nuclear repulsion energy     27888.4086998632 Hartrees.
NAtoms=  185 NActive=  185 NUniq=  185 SFac= 1.00D+00 NAtFMM=   60 NAOKFM=T Big=T
Integral buffers will be    131072 words long.
Raffenetti 2 integral format.
Two-electron integral symmetry is turned off.
One-electron integrals computed using PRISM.
NBasis=  3272 RedAO= T EigKep=  1.03D-06  NBF=  3272
NBsUse=  3181 1.00D-06 EigRej=  9.62D-07 NBFU=  3181
ExpMin= 3.15D-02 ExpMax= 1.14D+04 ExpMxC= 1.72D+03 IAcc=3 IRadAn=         5 AccDes= 0.00D+00
Harris functional with IExCor= 4538 and IRadAn=       5 diagonalized for initial guess.
HarFok:  IExCor= 4538 AccDes= 0.00D+00 IRadAn=         5 IDoV= 1 UseB2=F ITyADJ=14
ICtDFT=  3500011 ScaDFX=  1.000000  1.000000  1.000000  1.000000
FoFCou: FMM=F IPFlag=           0 FMFlag=      100000 FMFlg1=        2001
         NFxFlg=           0 DoJE=T BraDBF=F KetDBF=T FulRan=T
         wScrn=  0.000000 ICntrl=     500 IOpCl=  0 I1Cent=   200000004 NGrid=           0
         NMat0=    1 NMatS0=      1 NMatT0=    0 NMatD0=    1 NMtDS0=    0 NMtDT0=    0
Symmetry not used in FoFCou.
Requested convergence on RMS density matrix=1.00D-08 within1000 cycles.
Requested convergence on MAX density matrix=1.00D-06.
Requested convergence on             energy=1.00D-06.
No special actions if energy rises.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
Restarting incremental Fock formation.
请各位老师指点迷津,谢谢


1万

帖子

25

威望

1万

eV
积分
36401

管理员

公社社长

发表于 2018-6-4 10:25:35 | 显示全部楼层
用纯泛函做TDDFT计算是极为糟糕的,看
乱谈激发态的计算方法
http://sobereva.com/265

你这体系真有必要算200个态?算的越多耗时显著越长。设多大合适看
Gaussian中用TDDFT计算激发态和吸收、荧光、磷光光谱的方法
http://sobereva.com/314

没特殊必要的话,写Ultrafine完全是多余的,徒增大量耗时

第二个问题,仔细看
常见的多余的和被滥用的Gaussian关键词
http://sobereva.com/331
解决SCF不收敛问题的方法
http://sobereva.com/61
记得写#P便于监控收敛过程
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

35

帖子

0

威望

227

eV
积分
262

Level 3 能力者

 楼主| 发表于 2018-6-12 15:22:29 | 显示全部楼层
sobereva 发表于 2018-6-4 10:25
用纯泛函做TDDFT计算是极为糟糕的,看
乱谈激发态的计算方法
http://sobereva.com/265

不好意思没有及时回复,用纯泛函是为了做个比较,和杂化泛函长程矫正泛函做个比较,然后关于scf不收敛,请问#p是写在哪个位置呢

1万

帖子

25

威望

1万

eV
积分
36401

管理员

公社社长

发表于 2018-6-12 15:24:55 | 显示全部楼层
#后头写P
北京科音自然科学研究中心:http://www.keinsci.com  致力于计算化学的发展和传播,不定期开办各层次量子化学、分子动力学、波函数分析与Multiwfn程序等主题的培训。欢迎加入“北京科音”微信公众号获取培训最新消息和计算化学资讯
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(最流行的量子化学波函数分析程序)
计算化学公社论坛:http://bbs.keinsci.com(高水平、高人气、综合性计算化学交流论坛)
思想家公社QQ群1号:18616395,2号:466017436。用于讨论理论、计算化学,两个群讨论范畴相同,可加入任意其一但不可都加入,申请信息必须注明具体研究方向,否则一概不批。研究方向和理论、计算化学无关者勿加,以免浪费宝贵的空位

此账号为诸Sobereva共用
Money and papers are rubbish, get a real life!

35

帖子

0

威望

227

eV
积分
262

Level 3 能力者

 楼主| 发表于 2018-6-12 15:40:36 | 显示全部楼层

好的,谢谢啦
您需要登录后才可以回帖 登录 | 现在注册!

本版积分规则

手机版|北京科音自然科学研究中心|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949-1号 )

GMT+8, 2018-12-19 20:50 , Processed in 0.103557 second(s), 23 queries .

快速回复 返回顶部 返回列表