求助 ORCA里波函数稳定性测试及后续计算的问题

wsy151162

非常新新手，想请教各位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

fatpig

用上次的 .gbw 不停重复
%scf stabperform true end
这一步
直到一开始就显示stable为止...

zjxitcc

不要写GuessMode CMatrix。stabperform true不足以真的检验波函数稳定性，要这5行都写
STABPerform true
STABRestartUHFifUnstable true
STABMaxIter 500
STABDTol 1e-5
STABRTol 1e-5

wzkchem5

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，应当已经足够了

wsy151162

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

zjxitcc

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

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

wsy151162

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

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