计算化学公社

标题: 尝试将GAUSSIAN中的波函数导入ORCA进行计算，SCF还是不收敛，求各位大神指点 [打印本页]

作者
Author: 魔卡少女晗 时间: 2021-5-17 17:40
标题: 尝试将GAUSSIAN中的波函数导入ORCA进行计算，SCF还是不收敛，求各位大神指点
本帖最后由魔卡少女晗于 2021-5-17 20:23 编辑

按照老师的帖子我把GAUSSIAN的波函数文件导入orca，这趟一圈还是不收敛，请各位老师帮忙看看是哪里出了问题

gaussian的输入文件：
%chk=benzene-scf.chk
# PBE1PBE/def2QZVPP stable=opt nosymm

3

0 1
Ag                1.66317300 -0.56261200 0.20677800
Cl                0.55288200 1.77893200 -0.29237600
Cl                -2.74649000 1.33563600 -1.41781100
Cl                -1.83927300 0.07879800 1.61093300
Cl                -1.97844200 3.60090500 1.12790100

得到的out文件：
The wavefunction is stable under the perturbations considered.
The wavefunction is already stable.

ORCA的 inp文件：
%pal nprocs 16 end

!VERYSLOWCONV PWPB95 D3 ma-def2-QZVPP def2/JK def2-QZVPP/C RIJK grid6  noFinalGrid tightscf miniprint nopop pal4
%scf DirectResetFreq 1 end
%maxcore  1500
* xyz 0 1
Ag                1.66317300 -0.56261200 0.20677800
Cl                0.55288200 1.77893200 -0.29237600
Cl                -2.74649000 1.33563600 -1.41781100
Cl                -1.83927300 0.07879800 1.61093300
Cl                -1.97844200 3.60090500 1.12790100
Al                -1.51333100 1.70814600 0.27936700

然后FCHK文件转化为MKL，再转化为GBW文件，和INP文件一起放在目录里：

但是跑出来的结果还是不收敛，太难过了，是哪个步骤出了问题呢123  -3988.9322430035  -0.000127275453 0.00194245  0.00000613  0.0123555 0.9200
124  -3988.9323522213  -0.000109217764 0.00188797  0.00000594  0.0113542 0.9200

            *****************************************************
            *                   ERROR                      *
            *       SCF NOT CONVERGED AFTER 125 CYCLES       *
            *****************************************************

  Last Energy change       ... -1.0922e-04  Tolerance : 1.0000e-08
  Last MAX-Density change ... 1.6360e-03  Tolerance : 1.0000e-07

================================================================================
MOs in valence range (operator=0) identified as core orbitals.
Correct this by swapping MOs
                        MOs 46 and 43
                        MOs 47 and 44
                                             done
================================================================================

         **** DENSITY FILE WAS UPDATED (benzene-scf.scfp) ****
         **** ENERGY FILE WAS UPDATED (benzene-scf.en.tmp) ****

   --------------------------------------------------------------------
                                    WARNING
   The wavefunction IS NOT YET CONVERGED! It shows however signs of
   convergence. Therefore the wavefunction will be stored and can be
   used as input for another calculation.
   DO NOT USE THIS WAVEFUNCTION  FOR ANYHTING ELSE. It is NOT RELIABLE
   --------------------------------------------------------------------

作者
Author: hebrewsnabla 时间: 2021-5-17 17:58
泛函和基组完全不一样，凭什么能收敛呢。

作者
Author: zjxitcc 时间: 2021-5-17 18:11
类似帖子很多，一搜就有了。例如http://bbs.keinsci.com/thread-22911-1-1.html 我在9L中给出了操作示例。

作者
Author: 魔卡少女晗 时间: 2021-5-17 19:19

hebrewsnabla 发表于 2021-5-17 17:58
泛函和基组完全不一样，凭什么能收敛呢。

就是因为gaussian没有 PWPB95才想用ORCA的

，老师，那么只有基组相同可以用这个方法吗？是要给GJF文件加个弥散吗？

作者
Author: 魔卡少女晗 时间: 2021-5-17 19:23

zjxitcc 发表于 2021-5-17 18:11
类似帖子很多，一搜就有了。例如http://bbs.keinsci.com/thread-22911-1-1.html 我在9L中给出了操作示例。

老师，那我现在应该在高斯里面对DEF2QCVPP加个link1，变成MA-DEF2QZVPP之后再用这个波函数吗？泛函不一样，可以用GAUSSIAN输出的波函数吗

作者
Author: hdhxx123 时间: 2021-5-17 19:25

魔卡少女晗发表于 2021-5-17 19:19
就是因为gaussian没有 PWPB95才想用ORCA的，老师，那么只有基组相同可以用这个方法吗？是要给GJF文件 ...

先在小基组收敛再用大基组收敛，看这里http://bbs.keinsci.com/thread-20748-1-1.html

作者
Author: hebrewsnabla 时间: 2021-5-17 19:45

魔卡少女晗发表于 2021-5-17 19:19
就是因为gaussian没有 PWPB95才想用ORCA的，老师，那么只有基组相同可以用这个方法吗？是要给GJF文件 ...

基组可以一样，当然要用一样的。泛函应该选相近的。

饭要一口一口吃。先用相同泛函、相同基组试一下，如果没有问题，再尝试不同泛函、相同基组。
但是这里面注意点很多，不是一次就能成功的，估计你第一步就会出问题。

作者
Author: sobereva 时间: 2021-5-17 20:15
本来Gaussian就没有PWPB95-D3，因此明显不可能保证轨道传到ORCA后100%能收敛

作者
Author: wzkchem5 时间: 2021-5-17 21:20
出现MOs in valence range (operator=0) identified as core orbitals.这个提示，基本可以肯定是因为基组线性相关。存在基组线性相关的情况下，用高斯或者其他软件预收敛波函数是没用的。
解决基组线性相关的方法：（1）减小基组/辅助基组；（2）把积分阈值设严；（3）加大格点。
对于你的情况，你的体系大概率不需要ma-def2-QZVPP这么大的基组，事实上def2-QZVPP已经很弥散了，一般不会遇到def2-QZVPP还要加弥散的情况。另外如果你需要这么高精度的结果的话，可以考虑显式考虑Ag的相对论效应，赝势的误差对于你的需求可能不能忽略了。

作者
Author: 魔卡少女晗 时间: 2021-5-18 17:08

hdhxx123 发表于 2021-5-17 19:25
先在小基组收敛再用大基组收敛，看这里http://bbs.keinsci.com/thread-20748-1-1.html

非常有用，谢谢老师，我挨个试试看

作者
Author: 魔卡少女晗 时间: 2021-5-18 17:08

sobereva 发表于 2021-5-17 20:15
本来Gaussian就没有PWPB95-D3，因此明显不可能保证轨道传到ORCA后100%能收敛

嗯嗯好的，谢谢老师，我用其他方法尝试下

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