计算化学公社

标题: 求助在mrcc使用XYG3泛函做几何优化 [打印本页]

作者
Author: precious 时间: 2020-8-11 20:34
标题: 求助在mrcc使用XYG3泛函做几何优化
这是我的输入文件：

# XYG3 calculation for water with the cc-pVTZ basis set.
basis=cc-pVTZ
calc=XYG3
mem=500MB
gopt=full

unit=bohr
geom=xyz
3
H       0.143161836236       0.061417548106       0.000000000000
O       1.951596837673       0.062227994840       0.000000000000
H       2.397408854190       1.814851587354       0.000000000000

报错信息：（附件图片也有，更清晰一点）
Executing prop...
Allocation of 500.0 Mbytes of memory...
Analyitc gradients are not implemented for double hybrid functionals
for which the DFT contribution to the energy is not stationary wrt
the variation of the MO coefficients.

Fatal error in exec prop.
Program will stop.

求助！各位大神要怎么解决？？？

作者
Author: hebrewsnabla 时间: 2020-8-11 20:42
本帖最后由 hebrewsnabla 于 2020-8-11 21:30 编辑

该程序中没有该泛函的解析梯度，不能优化。

据说Q-Chem里面有XYGJ-OS的解析梯度。

作者
Author: 喵星大佬 时间: 2020-8-11 20:47
双杂化泛函没有解析梯度，英文写在上面

作者
Author: precious 时间: 2020-8-12 10:55

hebrewsnabla 发表于 2020-8-11 20:42
该程序中没有该泛函的解析梯度，不能优化。

据说Q-Chem里面有XYGJ-OS的解析梯度。

嗯，Q-chem好像不免费，好像nwchem也支持XYG3.

作者
Author: precious 时间: 2020-8-12 11:01

喵星大佬发表于 2020-8-11 20:47
双杂化泛函没有解析梯度，英文写在上面

谢谢，之前看到一个高斯通过设置Iop运行B2GPPLYP得方法。不晓得这个泛函可不可以通过设置Iop在高斯上面运行。C:\Users\User\Desktop\1.png

作者
Author: hebrewsnabla 时间: 2020-8-12 12:44
本帖最后由 hebrewsnabla 于 2020-8-12 13:00 编辑

precious 发表于 2020-8-12 11:01
谢谢，之前看到一个高斯通过设置Iop运行B2GPPLYP得方法。不晓得这个泛函可不可以通过设置Iop在高斯上面运 ...

这是算单点。梯度没那么简单。

一般来说没必要用双杂化来优化。如果执意要用，高斯、orca支持多种双杂化泛函的优化，如B2PLYP,B2GP-PLYP,DSD-PBEP86等，没必要抱着XYG3不放。

作者
Author: precious 时间: 2020-8-12 13:08

hebrewsnabla 发表于 2020-8-12 12:44
这是算单点。梯度没那么简单。

一般来说没必要用双杂化来优化。如果执意要用，高斯、orca支持多种双杂 ...

也是，主要导师执意想要试一试，谢谢！

作者
Author: hebrewsnabla 时间: 2020-8-12 13:21

precious 发表于 2020-8-12 13:08
也是，主要导师执意想要试一试，谢谢！

这个程序可以做XYG3的梯度： https://github.com/ajz34/Py_xDH/ ，需要结合PySCF。不过还不成熟，你可以问问作者有没有推荐的替代品。

作者
Author: sobereva 时间: 2020-8-17 09:22
用XYG3做优化没什么意思，此文已经体现出XYG3做优化精度比起很多常用的普通泛函比都吃亏：J. Chem. Theory Comput., 12, 459 (2016)，还明显不如用B2PLYP-D3(BJ)。

除前面提到的以外，Firefly也支持XYG3

作者
Author: precious 时间: 2020-8-18 14:11
谢谢各位的意见。

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