计算化学公社

 找回密码 Forget password
 注册 Register
Views: 10320|回复 Reply: 9
打印 Print 上一主题 Last thread 下一主题 Next thread

[Gaussian/gview] MaxDer out of range in HSE 高斯 l914 报错(已解决)

[复制链接 Copy URL]

6

帖子

0

威望

31

eV
积分
37

Level 2 能力者

关键词:# 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社长的指导!

sob01.png (50.86 KB, 下载次数 Times of downloads: 64)

sob01.png

19

帖子

0

威望

329

eV
积分
348

Level 3 能力者

10#
发表于 Post on 2018-11-19 10:38:23 | 只看该作者 Only view this author
sobereva 发表于 2018-11-19 10:19
没有完全写在明面上。找个小体系一试便知能不能算。CAM-B3LYP、ωB97XD、LC-PBE0、LC-BLYP都没这个问题。 ...

明白了,谢谢老师

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

9#
发表于 Post on 2018-11-19 10:19:27 | 只看该作者 Only view this author
灰灰不啦叽 发表于 2018-11-19 09:55
老师,请问是LC-wPBE泛函本身不支持三阶导数,还是这种调节w的方法导致它不支持三阶导数,还有类似的CAM- ...

没有完全写在明面上。找个小体系一试便知能不能算。CAM-B3LYP、ωB97XD、LC-PBE0、LC-BLYP都没这个问题。
对于G16 B.01来说,wPBE、PBEh、OHSE1PBE、OHSE2PBE都没有三阶导数
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

19

帖子

0

威望

329

eV
积分
348

Level 3 能力者

8#
发表于 Post on 2018-11-19 09:55:07 | 只看该作者 Only view this author
sobereva 发表于 2018-11-19 09:09
泛函的导数和能量的导数不是一码事
算TDDFT能量本身就已经需要泛函的二阶导数(通过矩阵元定义就可知) ...

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

490

帖子

2

威望

5306

eV
积分
5836

Level 6 (一方通行)

7#
发表于 Post on 2018-11-19 09:14:50 | 只看该作者 Only view this author
sobereva 发表于 2018-11-19 09:09
泛函的导数和能量的导数不是一码事
算TDDFT能量本身就已经需要泛函的二阶导数(通过矩阵元定义就可知) ...

原来如此。
谢谢sob~~~

6万

帖子

99

威望

6万

eV
积分
125141

管理员

公社社长

6#
发表于 Post on 2018-11-19 09:09:51 | 只看该作者 Only view this author
让你变成回忆 发表于 2018-11-17 11:31
他的意思就是利用优化w的方法来优化了激发态的结构和计算了频率,但是他提到的是用的LC-BLYP*这个泛函。
...

泛函的导数和能量的导数不是一码事
算TDDFT能量本身就已经需要泛函的二阶导数(通过矩阵元定义就可知),因此再多加一阶能量的导数,不支持三阶导数的泛函肯定就算不了了
北京科音自然科学研究中心http://www.keinsci.com)致力于计算化学的发展和传播,长期开办极高质量的各种计算化学类培训:初级量子化学培训班中级量子化学培训班高级量子化学培训班量子化学波函数分析与Multiwfn程序培训班分子动力学与GROMACS培训班CP2K第一性原理计算培训班,内容介绍以及往届资料购买请点击相应链接查看。这些培训是计算化学从零快速入门以及进一步全面系统性提升研究水平的高速路!培训各种常见问题见《北京科音办的培训班FAQ》
欢迎加入北京科音微信公众号获取北京科音培训的最新消息,并避免错过网上有价值的计算化学文章!
欢迎加入人气极高、专业性特别强的理论与计算化学综合交流群思想家公社QQ群(群号见此链接),合计达一万多人。北京科音培训班的学员在群中可申请VIP头衔,提问将得到群主Sobereva的最优先解答。
思想家公社的门口Blog:http://sobereva.com(发布大量原创计算化学相关博文)
Multiwfn主页:http://sobereva.com/multiwfn(十分强大、极为流行的量子化学波函数分析程序)
Google Scholar:https://scholar.google.com/citations?user=tiKE0qkAAAAJ
ResearchGate:https://www.researchgate.net/profile/Tian_Lu

490

帖子

2

威望

5306

eV
积分
5836

Level 6 (一方通行)

5#
发表于 Post on 2018-11-17 11:31:42 | 只看该作者 Only view this author
他的意思就是利用优化w的方法来优化了激发态的结构和计算了频率,但是他提到的是用的LC-BLYP*这个泛函。
我以前用过LC-wPBE(w进行了调控)来优化激发态的结构,报的错和你类似。
我也在群里看到了sob老师的回答,说了这个泛函没有三阶导数,但是我比较疑惑的是:优化不是需要一阶导数,频率计算不是需要二阶吗?为啥不能做呢?
还是请sob老师来详细解答把!

19

帖子

0

威望

329

eV
积分
348

Level 3 能力者

4#
发表于 Post on 2018-11-17 11:02:21 | 只看该作者 Only view this author
让你变成回忆 发表于 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

490

帖子

2

威望

5306

eV
积分
5836

Level 6 (一方通行)

3#
发表于 Post on 2018-11-17 10:50:25 | 只看该作者 Only view this author
灰灰不啦叽 发表于 2018-11-17 09:32
您好,我也在用lc-wpbe泛函计算激发态,遇到了错误提示为
The selected state is a singlet.
CISGrd: IG ...

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

19

帖子

0

威望

329

eV
积分
348

Level 3 能力者

2#
发表于 Post on 2018-11-17 09:32:46 | 只看该作者 Only view this author
您好,我也在用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 05:31 , Processed in 0.176367 second(s), 24 queries , Gzip On.

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