求助qchem中XYG3以及XYGJ-OS优化

liurui

老师，由于计算需要需要将不同泛函优化的结构做比较，也需要用双杂化泛函XYG3以及XYGJ-OS泛函来优化，但是看qchem说明书，XYG3以及XYGJ-OS泛函是基于B3LYP优化的结构来进行单点计算，老师想请问是否有软件用XYG3之类的双杂化泛函来进行优化，谢谢老师

wzkchem5

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

hebrewsnabla

wzk老师说的对。

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

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

ionexchangeC

结合特殊的IOp和Gaussian的调用功能，甚至你能在Gaussian里面做XYG3泛函的结构优化

hebrewsnabla

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

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

另外没有RI也很不划算。

liurui

谢谢老师们的回复，但是还有一点点问题，还是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
谢谢各位老师！！！

ionexchangeC

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

另外没有RI也很不划算。

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

hebrewsnabla

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

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

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

hxd_yi

NWChem应该可以。
但是这些个泛函算起来太慢了，还是建议只算单点。

liurui

谢谢各位老师的回复，帮助非常大，学生去试着试一下。