计算化学公社

标题: 求助qchem中XYG3以及XYGJ-OS优化 [打印本页]

作者
Author: liurui 时间: 2022-5-18 22:08
标题: 求助qchem中XYG3以及XYGJ-OS优化
老师，由于计算需要需要将不同泛函优化的结构做比较，也需要用双杂化泛函XYG3以及XYGJ-OS泛函来优化，但是看qchem说明书，XYG3以及XYGJ-OS泛函是基于B3LYP优化的结构来进行单点计算，老师想请问是否有软件用XYG3之类的双杂化泛函来进行优化，谢谢老师

作者
Author: wzkchem5 时间: 2022-5-18 22:25
不是的，XYG3、XYGJ-OS是基于B3LYP的SCF波函数来进行的MP2校正（不管在算单点时还是在做结构优化时都是如此），而不是基于B3LYP的结构算的单点。所以XYG3、XYGJ-OS是可以用于结构优化的

作者
Author: hebrewsnabla 时间: 2022-5-18 22:36
wzk老师说的对。

Qchem应该可以做XYGJ-OS的优化。

除此以外，我们组的基于PySCF的开源程序 github.com/ajz34/dh 可以做所有XYG3型泛函的优化，包括XYG3和XYGJ-OS。

作者
Author: ionexchangeC 时间: 2022-5-18 22:46
结合特殊的IOp和Gaussian的调用功能，甚至你能在Gaussian里面做XYG3泛函的结构优化

作者
Author: hebrewsnabla 时间: 2022-5-18 22:57

ionexchangeC 发表于 2022-5-18 22:46
结合特殊的IOp和Gaussian的调用功能，甚至你能在Gaussian里面做XYG3泛函的结构优化

解析梯度不行吧？数值梯度是可以的。

另外没有RI也很不划算。

作者
Author: liurui 时间: 2022-5-18 23:03
谢谢老师们的回复，但是还有一点点问题，还是XYG3以及XYGJ-OS在qchem中的优化结构，做了B3LYP的优化计算，然后基于此进行优化，qchem说明书写的如下，
$molecule
0 1
N 0.000000 0.000000 0.547775
N 0.000000 0.000000 -0.547775
$end
$rem
EXCHANGE xyg3
BASIS 6-311+G(3df,2p)
$end
按照这个计算，得到的是能量计算的结果，并没有结构优化结果，然后试了一下在命令里加了 JOB_TYPE  =  Optimization命令，就报错不能计算了，命令为
$rem
BASIS  =  aug-cc-pvDZ
GUI  =  2
  JOB_TYPE  =  Optimization
EXCHANGE  = xyg3
CC_PRINT_PREC = 16
SCF_ALGORITHM  =  DIIS
SCF_CONVERGENCE  =  8
CC_MEMORY = 8000 !8 GB
mem_total = 8500
$end
谢谢各位老师！！！

作者
Author: ionexchangeC 时间: 2022-5-18 23:05

hebrewsnabla 发表于 2022-5-18 22:57
解析梯度不行吧？数值梯度是可以的。

另外没有RI也很不划算。

确实，极度不划算，也就理论上可行

作者
Author: hebrewsnabla 时间: 2022-5-18 23:22

liurui 发表于 2022-5-18 23:03
谢谢老师们的回复，但是还有一点点问题，还是XYG3以及XYGJ-OS在qchem中的优化结构，做了B3LYP的优化计算， ...

提问必须说清楚报错信息。

我印象中qchem里面只实现了XYGJ-OS优化。

作者
Author: hxd_yi 时间: 2022-5-19 00:00
NWChem应该可以。
但是这些个泛函算起来太慢了，还是建议只算单点。

作者
Author: liurui 时间: 2022-5-19 10:02
谢谢各位老师的回复，帮助非常大，学生去试着试一下。

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