计算化学公社

 找回密码 Forget password
 注册 Register

MaxDer out of range in HSE 高斯 l914 报错(已解决)

查看数: 10317 | 评论数: 9 | 收藏 Add to favorites 5
关灯 | 提示:支持键盘翻页<-左 右->
    组图打开中,请稍候......
发布时间: 2018-4-24 14:04

正文摘要:

关键词:# td(nstates=20)  6-31g(d,p)  scrf=conver=6  pop=nbo6  density  em=gd3bj  LC-wPBE  IOp(3/107=0196500000,3/108=0196500000) ...

回复 Reply

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

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

没有完全写在明面上。找个小体系一试便知能不能算。CAM-B3LYP、ωB97XD、LC-PBE0、LC-BLYP都没这个问题。
对于G16 B.01来说,wPBE、PBEh、OHSE1PBE、OHSE2PBE都没有三阶导数
灰灰不啦叽 发表于 Post on 2018-11-19 09:55:07
sobereva 发表于 2018-11-19 09:09
泛函的导数和能量的导数不是一码事
算TDDFT能量本身就已经需要泛函的二阶导数(通过矩阵元定义就可知) ...

老师,请问是LC-wPBE泛函本身不支持三阶导数,还是这种调节w的方法导致它不支持三阶导数,还有类似的CAM-B3LYP、ωB97XD、LC-PBE0、LC-BLYP这些范围分离泛函如何判断是否支持激发态优化,老师可以推荐一下查这个性质的资料的方式吗
让你变成回忆 发表于 Post on 2018-11-19 09:14:50
sobereva 发表于 2018-11-19 09:09
泛函的导数和能量的导数不是一码事
算TDDFT能量本身就已经需要泛函的二阶导数(通过矩阵元定义就可知) ...

原来如此。
谢谢sob~~~
sobereva 发表于 Post on 2018-11-19 09:09:51
让你变成回忆 发表于 2018-11-17 11:31
他的意思就是利用优化w的方法来优化了激发态的结构和计算了频率,但是他提到的是用的LC-BLYP*这个泛函。
...

泛函的导数和能量的导数不是一码事
算TDDFT能量本身就已经需要泛函的二阶导数(通过矩阵元定义就可知),因此再多加一阶能量的导数,不支持三阶导数的泛函肯定就算不了了
让你变成回忆 发表于 Post on 2018-11-17 11:31:42
他的意思就是利用优化w的方法来优化了激发态的结构和计算了频率,但是他提到的是用的LC-BLYP*这个泛函。
我以前用过LC-wPBE(w进行了调控)来优化激发态的结构,报的错和你类似。
我也在群里看到了sob老师的回答,说了这个泛函没有三阶导数,但是我比较疑惑的是:优化不是需要一阶导数,频率计算不是需要二阶吗?为啥不能做呢?
还是请sob老师来详细解答把!
灰灰不啦叽 发表于 Post on 2018-11-17 11:02:21
让你变成回忆 发表于 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
让你变成回忆 发表于 Post on 2018-11-17 10:50:25
灰灰不啦叽 发表于 2018-11-17 09:32
您好,我也在用lc-wpbe泛函计算激发态,遇到了错误提示为
The selected state is a singlet.
CISGrd: IG ...

这个泛函好像我以前用过,记得是不能够用于激发态优化。
灰灰不啦叽 发表于 Post on 2018-11-17 09:32:46
您好,我也在用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,但是并不是很理解这是什么意思
希望老师您能帮我解惑,谢谢

手机版 Mobile version|北京科音自然科学研究中心 Beijing Kein Research Center for Natural Sciences|京公网安备 11010502035419号|计算化学公社 — 北京科音旗下高水平计算化学交流论坛 ( 京ICP备14038949号-1 )|网站地图

GMT+8, 2026-2-21 02:29 , Processed in 0.175132 second(s), 25 queries , Gzip On.

快速回复 返回顶部 返回列表 Return to list