计算化学公社

标题: 求助 ORCA里波函数稳定性测试及后续计算的问题 [打印本页]

作者
Author: wsy151162 时间: 2024-11-25 09:52
标题: 求助 ORCA里波函数稳定性测试及后续计算的问题
本帖最后由 wsy151162 于 2024-11-26 09:29 编辑

非常新新手，想请教各位ORCA6在波函数稳定性测试和后续计算的问题
通过查看论坛内的各类帖子和手册查到orca6可以使用
%scf stabperform true end  做稳定性测试，刚开始显示波函数不稳定，但是后面它自己重新又收敛到了稳定的波函数。
接下来我就用！MORead %moinp ".gbw" 和%scf GuessMode CMatrix end 读取在稳定性测试中得到的波函数做tddft, 得到的T1依然是负值后，我就对上一步得到的波函数再做了一次稳定性测试，第一轮仍然是无法得到稳定的波函数，需要到第二轮，感觉有点鬼打墙了，想请教一下各位要怎么样利用到这稳定的波函数去做下一步呢？

文件包含两个稳定性测试的inp和out文件，非常感谢各位。

Orbitals have been transformed and reconvergence is now attempted

               **************************************************
               * The wavefunction in try  1/ 5 is unstable.    *
               *       Trying to re-converging now          *
               **************************************************

----------------------------------------D-I-I-S--------------------------------------------
Iteration Energy (Eh)          Delta-E RMSDP    MaxDP    DIISErr Damp  Time(sec)
-------------------------------------------------------------------------------------------
            ***  Starting incremental Fock matrix formation  ***
1 -55867.9023760067648254    0.00e+00  3.35e-05  7.02e-03  3.36e-03  0.700 166.1
Warning: op=0 Small HOMO/LUMO gap ( 0.052) - skipping pre-diagonalization
      Will do a full diagonalization
2 -55867.9034939686825965 -1.12e-03  3.68e-05  7.41e-03  3.56e-03  0.700  88.6
                           ***Turning on AO-DIIS***
3 -55867.9052610028884374 -1.77e-03  3.37e-05  6.21e-03  3.67e-03  0.700  78.0
4 -55867.9070216857799096 -1.76e-03  1.15e-04  1.91e-02  3.59e-03  0.000  77.3
5 -55867.9138436696375720 -6.82e-03  9.31e-05  2.59e-02  2.44e-03  0.000  86.6
6 -55867.9153572440409334 -1.51e-03  8.76e-05  2.39e-02  1.19e-03  0.000  82.7
7 -55867.9158106163522461 -4.53e-04  3.54e-05  1.13e-02  8.16e-04  0.000  82.4
            *** Restarting incremental Fock matrix formation ***
                           ****Resetting DIIS****
Warning: op=0 Small HOMO/LUMO gap ( 0.078) - skipping pre-diagonalization
      Will do a full diagonalization
8 -55867.9158250061736908 -1.44e-05  1.06e-04  1.81e-02  4.75e-04  0.000 179.4
9 -55867.9158237872034078    1.22e-06  5.58e-06  1.61e-03  9.33e-04  0.000 119.6
10 -55867.9158276784146437 -3.89e-06  2.98e-05  8.65e-03  1.43e-04  0.000  86.5
11 -55867.9158325425523799 -4.86e-06  1.53e-05  4.64e-03  1.33e-04  0.000  95.9
12 -55867.9158315450040391    9.98e-07  1.09e-05  2.34e-03  1.33e-04  0.000  91.1
               **** Energy Check signals convergence ****

            *****************************************************
            *                   SUCCESS                      *
            *          SCF CONVERGED AFTER  12 CYCLES       *
            *****************************************************

Recomputing exchange energy using gridx3    ... done (  180.685 sec)
Old exchange energy                         : -318.973661956 Eh
New exchange energy                         : -318.973620204 Eh
Exchange energy change after final integration :    0.000041752 Eh
Total energy after final integration          : -55867.915789793 Eh
         **** ENERGY FILE WAS UPDATED (A-I-triF-fromg.en.tmp) ****

-------------------------------------------------------------------------------------------
                           WAVEFUNCTION STABILITY ANALYSIS
-------------------------------------------------------------------------------------------

                     ****Iteration 0****
Lowest Energy       :    0.056442185673
Maximum Energy change  :    0.060887386274 (vector 2)
Maximum residual norm  :    0.005148511405
（省略）……
                     ****Iteration 8****
Lowest Energy       :    0.047129713967
Maximum Energy change  :    0.000051868459 (vector 1)
Maximum residual norm  :    0.000019153692

   *** CONVERGENCE OF RESIDUAL NORM REACHED ***
The eigenvalues of the stability matrix:
E( 0) = 0.04712971 Eh
E( 1) = 0.04759661 Eh
E( 2) = 0.05535287 Eh

The stability analysis shows that the wavefunction is stable

作者
Author: fatpig 时间: 2024-11-25 10:01
用上次的 .gbw 不停重复
%scf stabperform true end
这一步
直到一开始就显示stable为止...

作者
Author: zjxitcc 时间: 2024-11-25 10:38
不要写GuessMode CMatrix。stabperform true不足以真的检验波函数稳定性，要这5行都写
STABPerform true
STABRestartUHFifUnstable true
STABMaxIter 500
STABDTol 1e-5
STABRTol 1e-5

作者
Author: wzkchem5 时间: 2024-11-25 12:17

zjxitcc 发表于 2024-11-25 03:38
不要写GuessMode CMatrix。stabperform true不足以真的检验波函数稳定性，要这5行都写
STABPerform true
...

有没有测试过ORCA 6是否仍然必须写这5行？ORCA 6的SCF代码是从头写的，其中波函数稳定性分析部分也是从头写的，我个人测试过对于某些小体系而言，默认的stability analysis收敛阈值可以把orbital Hessian的最小本征值收敛到1e-6，应当已经足够了

作者
Author: wsy151162 时间: 2024-11-25 20:24

zjxitcc 发表于 2024-11-25 10:38
不要写GuessMode CMatrix。stabperform true不足以真的检验波函数稳定性，要这5行都写
STABPerform true
...

谢谢老师，我加上这5行之后又做了2次，它还是一样的结果，第一轮不是稳定的，第二轮得到稳定的波函数，它们的输出文件中连收敛过程都非常地类似，跟没有加其他几行命令的时候也非常类似。得到能量值也没有明显的变化，想请问一下老师这是正常的吗？我是继续用上次的 .gbw 不停重复吗？非常感谢！

下面是输出文件的节选
以下是第一轮不稳定的数据
stable 3:
The eigenvalues of the stability matrix:
E( 0) = -0.04929669 Eh
E( 1) = 0.03432445 Eh
E( 2) = 0.03758749 Eh

stable 4:
The eigenvalues of the stability matrix:
E( 0) = -0.04929664 Eh
E( 1) = 0.03432367 Eh
E( 2) = 0.03758801 Eh

以下是第二轮稳定后的数据
stable 3:
The eigenvalues of the stability matrix:
E( 0) = 0.04712420 Eh
E( 1) = 0.04757106 Eh
E( 2) = 0.05533941 Eh

stable 4:
The eigenvalues of the stability matrix:
E( 0) = 0.04712808 Eh
E( 1) = 0.04757454 Eh
E( 2) = 0.05533986 Eh

作者
Author: zjxitcc 时间: 2024-11-25 20:30

wsy151162 发表于 2024-11-25 20:24
谢谢老师，我加上这5行之后又做了2次，它还是一样的结果，第一轮不是稳定的，第二轮得到稳定的波函数，它 ...

你有注意到3L回复中“不要写GuessMode CMatrix”吗？如果有注意到，可以上传输入、输出文件看看。

作者
Author: wsy151162 时间: 2024-11-26 09:36

zjxitcc 发表于 2024-11-25 20:30
你有注意到3L回复中“不要写GuessMode CMatrix”吗？如果有注意到，可以上传输入、输出文件看看。

老师您好，我注意到了，在stable3和stable4里面已经删掉了GuessMode CMatrix，我把它们的输入和输出文件放到stable3 4的压缩包里上传进来了，非常感谢老师

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