计算化学公社

标题: MaxDer out of range in HSE 高斯 l914 报错（已解决） [打印本页]

作者
Author: lyc460360210 时间: 2018-4-24 14:04
标题: MaxDer out of range in HSE 高斯 l914 报错（已解决）
关键词：# td(nstates=20) 6-31g(d,p) scrf=conver=6 pop=nbo6 density em=gd3bj LC-wPBE IOp(3/107=0196500000,3/108=0196500000)

错误提示： MaxDer out of range in HSE.
Error termination via Lnk1e in /home/g16/l914.exe

解决办法：如图。感谢sob社长的指导！

作者
Author: 灰灰不啦叽 时间: 2018-11-17 09:32
您好，我也在用lc-wpbe泛函计算激发态，遇到了错误提示为
The selected state is a singlet.
CISGrd: IGrad=2 NXY=2 DFT=T
CISAX will form    1 AO SS matrices at one time.
NMat=    1 NSing=    1 JSym2X= 0.
MaxDer out of range in HSEd.
Error termination via Lnk1e in /public1/soft/gaussian/G16_x64_AVX2/g16/l914.exe

我计算的命令是
#p td(singlets,root=1,nstates=6) opt  LC-wPBE/def2SVP iop(3/107=0160900000,3/108=0160900000)

我想请问下社长说的lc-wpbe泛函不能产生激发态密度，是通过哪个关键词控制然后解决问题的。
我查阅的高斯手册，找到Acc2E=N 是设定双电子积分的精度参数为N，但是并不是很理解这是什么意思
希望老师您能帮我解惑，谢谢

作者
Author: 让你变成回忆 时间: 2018-11-17 10:50

灰灰不啦叽发表于 2018-11-17 09:32
您好，我也在用lc-wpbe泛函计算激发态，遇到了错误提示为
The selected state is a singlet.
CISGrd: IG ...

这个泛函好像我以前用过，记得是不能够用于激发态优化。

作者
Author: 灰灰不啦叽 时间: 2018-11-17 11:02

让你变成回忆发表于 2018-11-17 10:50
这个泛函好像我以前用过，记得是不能够用于激发态优化。

According to the above methodology, the tuning optDFTw
procedure based on the RS functionals (LC-BLYP and LC-oPBE)
with the 6-31+G(d) basis set is performed to determine the
optimal o values when J2 reaches the minimum;18a,34,35 henceforth,
we denote the optimally-tuned RS functionals as
LC-BLYP* and LC-oPBE*, which were applied throughout this
work. The ground geometries were then reoptimized using the
new o values, for the sake of comparison, and the calculations
were also performed with the widely used non-optimized
CAM-B3LYP functional;31 and the excited state structures were
calculated using the time-dependent TD-LC-BLYP* functional
with the 6-31+G(d) basis set, followed by calculations of
harmonic vibrational frequencies and normal modes to obtain
equilibrium geometries and to calculate the thermal vibration
correlation functions. All calculations were performed using
the Gaussian 09 software.
我是从一篇文献上看到作者这么处理的，自己也想要做一下尝试

文献：Phys. Chem. Chem. Phys., 2018, 20, 6548--6561

作者
Author: 让你变成回忆 时间: 2018-11-17 11:31
他的意思就是利用优化w的方法来优化了激发态的结构和计算了频率，但是他提到的是用的LC-BLYP*这个泛函。
我以前用过LC-wPBE（w进行了调控）来优化激发态的结构，报的错和你类似。
我也在群里看到了sob老师的回答，说了这个泛函没有三阶导数，但是我比较疑惑的是：优化不是需要一阶导数，频率计算不是需要二阶吗？为啥不能做呢？
还是请sob老师来详细解答把！

作者
Author: sobereva 时间: 2018-11-19 09:09

让你变成回忆发表于 2018-11-17 11:31
他的意思就是利用优化w的方法来优化了激发态的结构和计算了频率，但是他提到的是用的LC-BLYP*这个泛函。
...

泛函的导数和能量的导数不是一码事
算TDDFT能量本身就已经需要泛函的二阶导数（通过矩阵元定义就可知），因此再多加一阶能量的导数，不支持三阶导数的泛函肯定就算不了了

作者
Author: 让你变成回忆 时间: 2018-11-19 09:14

sobereva 发表于 2018-11-19 09:09
泛函的导数和能量的导数不是一码事
算TDDFT能量本身就已经需要泛函的二阶导数（通过矩阵元定义就可知） ...

原来如此。
谢谢sob~~~

作者
Author: 灰灰不啦叽 时间: 2018-11-19 09:55

sobereva 发表于 2018-11-19 09:09
泛函的导数和能量的导数不是一码事
算TDDFT能量本身就已经需要泛函的二阶导数（通过矩阵元定义就可知） ...

老师，请问是LC-wPBE泛函本身不支持三阶导数，还是这种调节w的方法导致它不支持三阶导数，还有类似的CAM-B3LYP、ωB97XD、LC-PBE0、LC-BLYP这些范围分离泛函如何判断是否支持激发态优化，老师可以推荐一下查这个性质的资料的方式吗

作者
Author: sobereva 时间: 2018-11-19 10:19

灰灰不啦叽发表于 2018-11-19 09:55
老师，请问是LC-wPBE泛函本身不支持三阶导数，还是这种调节w的方法导致它不支持三阶导数，还有类似的CAM- ...

没有完全写在明面上。找个小体系一试便知能不能算。CAM-B3LYP、ωB97XD、LC-PBE0、LC-BLYP都没这个问题。
对于G16 B.01来说，wPBE、PBEh、OHSE1PBE、OHSE2PBE都没有三阶导数

作者
Author: 灰灰不啦叽 时间: 2018-11-19 10:38

sobereva 发表于 2018-11-19 10:19
没有完全写在明面上。找个小体系一试便知能不能算。CAM-B3LYP、ωB97XD、LC-PBE0、LC-BLYP都没这个问题。 ...

明白了，谢谢老师

欢迎光临计算化学公社 (http://bbs.keinsci.com/)