orca SCF迷一样的感受

413

各位大家好我最近优化了一个结构，跑了一天，然后拿优化好的结构重新做个单点算点性质
奇怪的地方是未优化的结构SCF很快就收敛了
orca自己优化好的结构，SCF竟然报错
我也是也将之前的gbw文件拷贝成新文件名的gbw，但是仍未解决问题
请问大家有这种体验吗？怎么解决这种问题吗？或者我哪里做的不对才会出现这种情况？

喵星大佬

看看TRAH启动了没，试试吧这玩意先关了

wzkchem5

报的是什么错？SCF报错不一定是不收敛，也有可能是其他原因，只有看具体报错原因才能知道是怎么回事

wzkchem5

喵星大佬 发表于 2022-3-9 17:53
看看TRAH启动了没，试试吧这玩意先关了

TRAH虽然有时会增加迭代步数，但是几乎从不会使一个本来能收敛的计算变成不能收敛。
当然某些泛函缺少kernel导致TRAH报错之类的情况除外，但是那也是看了输出文件知道是这个原因以后才关TRAH，不建议不分析输出文件就关TRAH来解决SCF报错问题

冰释之川

wzkchem5 发表于 2022-3-10 01:01
TRAH虽然有时会增加迭代步数，但是几乎从不会使一个本来能收敛的计算变成不能收敛。
当然某些泛函缺少ke ...

下一个版本能否让不支持的泛函自动关闭TRAH，不然用起来不太人性化啊。。

413

喵星大佬 发表于 2022-3-10 00:53
看看TRAH启动了没，试试吧这玩意先关了

请问怎么检查这个TRAH启动了没？

413

倒没报错我用了跟优化一样的输入文件，只是去掉了opt关键词
SCF的能量就不太对，我就给停掉了，试了多次，也跑过更多的圈数，也还是一样的结果。
还试过把opt的gbw文件拷过来，也一样，现在再试高斯的chk转化过来

zjxitcc

建议在1L中多展示一些信息，例如ORCA版本，你的体系是否含过渡金属，你用的什么基组

413

wzkchem5 发表于 2022-3-10 00:59
报的是什么错？SCF报错不一定是不收敛，也有可能是其他原因，只有看具体报错原因才能知道是怎么回事

体系为碳氢吸附在金银合金团簇上


------------------------------
INITIAL GUESS: MODEL POTENTIAL
------------------------------
Loading Hartree-Fock densities                     ... done
  calling /public/home/users/orca_4.2.1/orca_4_2_1_linux_x86-64_shared_openmpi314/orca atom47.inp > atom47.out in order to generate an atomic fitting density for atom 0 (Ag with ECP) on-the-fly...
    atom 0 (Ag), assumed electronic state with S=2: 1s2 2s2 2p6 3s2 3p6 4s2 3d10 4p6 5s2 4d9           ... done
  calling /public/home/users/orca_4.2.1/orca_4_2_1_linux_x86-64_shared_openmpi314/orca atom79.inp > atom79.out in order to generate an atomic fitting density for atom 12 (Au with ECP) on-the-fly...
    atom 12 (Au), assumed electronic state with S=2: 1s2 2s2 2p6 3s2 3p6 4s2 3d10 4p6 5s2 4d10 5p6 6s1 4f14 5d10           ... done
Calculating cut-offs                               ... done
Setting up the integral package                    ... done
Initializing the effective Hamiltonian             ... done
Starting the Coulomb interaction                   ... done (   1.4 sec)
Reading the grid                                   ... done
Mapping shells                                     ... done
Starting the XC term evaluation                    ... done (   0.7 sec)
  promolecular density results
     # of electrons  =   3246.018386740
     EX              = -12947.568549971
     EC              =   -185.685677767
     EX+EC           = -13133.254227738
Transforming the Hamiltonian                       ... done (   4.8 sec)
Diagonalizing the Hamiltonian                      ... done (   7.0 sec)
Back transforming the eigenvectors                 ... done (   0.2 sec)
Now organizing SCF variables                       ... done
                      ------------------
                      INITIAL GUESS DONE (  23.6 sec)
                      ------------------
--------------
SCF ITERATIONS
--------------
ITER       Energy         Delta-E        Max-DP      RMS-DP      [F,P]     Damp
               ***  Starting incremental Fock matrix formation  ***
  0  -7600.7457555027   0.0000000000004856.57061663 20.53065389  0.3333300 0.7000
  1  -7556.3272119007  44.4185436019625453.90642201 20.34005235  0.3486989 0.7000
  2  -7361.3936161631 194.93359573764610059.00072807 36.13300185  0.5636771 0.7000

ORCA finished by error termination in SCF

renzhogn424

4.2.1，那没有TRAH的，好像是能量变化过大造成报错

413

renzhogn424 发表于 2022-3-10 10:51
4.2.1，那没有TRAH的，好像是能量变化过大造成报错

想请教下基组这么指定对吗，这样跑倒是没报错。
! pbe angs PRINTBASIS PRINTMOS SARC/J
%maxcore 5000
%scf MaxIter 512 Thresh 1e-6 end
%pal nprocs 28 end
%basis
NewGTO Au "SARC-ZORA-TZVP" end #
DelECP Au
NewGTO Ag "SARC-ZORA-TZVP" end #
DelECP Ag
NewGTO C "ZORA-def2-SVP" end #
DelECP C
NewGTO H "ZORA-def2-SVP" end #
DelECP H
end

zjxitcc

413 发表于 2022-3-10 11:48
想请教下基组这么指定对吗，这样跑倒是没报错。
! pbe angs PRINTBASIS PRINTMOS SARC/J
%maxcore 5000 ...

你没开启相对论哈密顿，却使用了相对论全电子基组。你检查过输出文件里有开启ZORA哈密顿么？

413

zjxitcc 发表于 2022-3-10 12:00
你没开启相对论哈密顿，却使用了相对论全电子基组。你检查过输出文件里有开启ZORA哈密顿么？

请问正确的格式是再加个ZORA吗？

! pbe angs PRINTBASIS PRINTMOS SARC/J ZORA
%maxcore 5000
%scf MaxIter 512 Thresh 1e-6 end
%pal nprocs 28 end
%basis
NewGTO Au "SARC-ZORA-TZVP" end #
DelECP Au
NewGTO Ag "SARC-ZORA-TZVP" end #
DelECP Ag
NewGTO C "ZORA-def2-SVP" end #
DelECP C
NewGTO H "ZORA-def2-SVP" end #
DelECP H
end

hebrewsnabla

在结构优化过程中，后一步会读前一步的波函数，因此收敛比较容易。而取出收敛的结构单独做SCF的话，就没有这一步，是从头生成初猜的，所以不收敛也是可以理解的。

还试过把opt的gbw文件拷过来，也一样，

这不太可能，应该立即收敛才对。检查你的输出文件中有没有显示读入波函数成功，或者有没有拷错文件。

renzhogn424

413 发表于 2022-3-10 11:48
想请教下基组这么指定对吗，这样跑倒是没报错。
! pbe angs PRINTBASIS PRINTMOS SARC/J
%maxcore 5000 ...

这个相对论基组我倒是没有用过，不大清楚