过渡金属掺杂团簇检测产生稳定性波函数报错l508

lulu顺利毕业

各位老师，我在使用bp86/lanl2dz scf=xqc stable=opt对过渡金属掺杂的碳团簇产生稳定波函数时，结果出现L508报错，No lower point found -- run aborted，应该如何尝试解决？下面是我检测产生稳定波函数的输出文件，还请各位老师多多指教。
# bp86/lanl2dz scf=xqc stable=opt
...
>>>>>>>>>> Convergence criterion not met.
SCF Done:  E(UB-P86) =  -3857.38574397     A.U. after   33 cycles
            NFock= 32  Conv=0.51D-02     -V/T= 2.1744
<Sx>= 0.0000 <Sy>= 0.0000 <Sz>= 3.0000 <S**2>=12.0553 S= 3.0079
<L.S>=  0.00000000000   
Annihilation of the first spin contaminant:
S**2 before annihilation    12.0553,   after    12.0009
Defaulting to unpruned grid for atomic number  74.
...
Defaulting to unpruned grid for atomic number  74.
Gradient too large for Newton-Raphson -- use scaled steepest descent instead.
     Accept linear search using points  1 and  2.
...
Gradient too large for Newton-Raphson -- use scaled steepest descent instead.
     Minimum is close to point  5 DX= -6.14D-03 DF= -9.10D-05 DXR=  1.56D-02 DFR=  2.50D-04 which will be used.
Gradient too large for Newton-Raphson -- use scaled steepest descent instead.
     Minimum is close to point  2 DX= -7.93D-03 DF= -2.93D-04 DXR=  4.13D-02 DFR=  1.71D-03 which will be used.
Gradient too large for Newton-Raphson -- use scaled steepest descent instead.
...

Gradient too large for Newton-Raphson -- use scaled steepest descent instead.
LinEq1:  Iter=  0 NonCon=     1 RMS=5.77D-03 Max=2.63D+00 NDo=     1
AX will form     1 AO Fock derivatives at one time.

LinEq1:  Iter=  1 NonCon=     1 RMS=2.30D-02 Max=1.03D+01 NDo=     1
...
Linear equations converged to 6.540D-08 6.540D-07 after   233 iterations.
Gradient too large for Newton-Raphson -- use scaled steepest descent instead.
Gradient too large for Newton-Raphson -- use scaled steepest descent instead.
     Minimum is close to point  2 DX=  1.76D-03 DF= -2.12D-08 DXR=  8.71D-03 DFR=  7.59D-05 which will be used.
Gradient too large for Newton-Raphson -- use scaled steepest descent instead.
Restarting incremental Fock formation.
Gradient too large for Newton-Raphson -- use scaled steepest descent instead.
     Minimum is close to point  2 DX= -1.59D-02 DF= -1.92D-06 DXR=  8.64D-02 DFR=  7.53D-03 which will be used.
...
Linear equations converged to 2.407D-10 2.407D-09 after   225 iterations.
Restarting incremental Fock formation.
Search did not lower the energy significantly.
No lower point found -- try reversing direction.
...
Restarting incremental Fock formation.
Search did not lower the energy significantly.
No lower point found -- run aborted.
Error termination via Lnk1e in /home/gauss16/g16/l508.exe at Tue Nov 21 00:47:21 2023.
Job cpu time:      98 days 18 hours 20 minutes 49.0 seconds.
Elapsed time:       6 days  4 hours 57 minutes 56.1 seconds.
File lengths (MBytes):  RWF=   4876 Int=      0 D2E=      0 Chk=     31 Scr=  

zjxitcc

如此暴力计算过于可怕，一个SCF算了6天4h报错，表明计算不合理，用户不自知。表面上看是SCF难以收敛，实际上有很多问题，因此不能简单当做SCF不收敛处理。
（1）应交待体系电荷、自旋多重度等信息。如果是单重态你这么算，我会痛骂。
（2）检查几何结构是否合理。
（3）LANL2DZ在2023年已经太过时了，如果体系不小，可以给过渡金属原子用def2TZVP，非金属原子用6-31G(d,p)。而且把LANL2DZ写在关键词行时，高斯默认开启6D 10F，极易导致基函数线性依赖，虽然高斯会自动砍去线性依赖的基函数，但SCF收敛仍然不易。
如果体系不涉密、可以让别人看的话，最好上传压缩后的输入输出文件。

ahxb

你的”过渡金属掺杂的碳团簇“体系中，碳团簇有多大？过渡金属怎么掺杂的？

lulu顺利毕业

zjxitcc 发表于 2023-11-23 11:00
如此暴力计算过于可怕，一个SCF算了6天4h报错，表明计算不合理，用户不自知。表面上看是SCF难以收敛，实际 ...

感谢老师的回复。关于您的解答我还有一些问题。
（1）体系电荷数是648，上面的输出文件是自旋多重度7。闭壳层结构，应该是没有磁性，但是就目前，我在自旋多重度3.5下已经通过产生稳定波函数，优化，检测稳定性得到了没有虚频的结构。发现能量要比在自旋多重度1下得到的结构能量低。
（2）如何检查几何结构是否合理，结构并不是我自己构建的。
（3）关于您说的改变泛函和基组的问题。我可以尝试一下，看看是否优化结果会更好一些。
由于我现在没有任何人指导，全都是自己摸索着做，之前已经做了很多无用功，也没找到通往正确的道理。还望您能多说一些，计算方法也好，理论知识也好，再次非常感谢您的回复和解答。

lulu顺利毕业

ahxb 发表于 2023-11-23 15:32
你的”过渡金属掺杂的碳团簇“体系中，碳团簇有多大？过渡金属怎么掺杂的？

Mo12C80中空笼状团簇

zjxitcc

lulu顺利毕业 发表于 2023-11-23 19:42
感谢老师的回复。关于您的解答我还有一些问题。
（1）体系电荷数是648，上面的输出文件是自旋多重度7。 ...

（1）自旋多重度3.5 是什么意思？自旋多重度是2S+1，如果体系是偶数个电子，自旋多重度只有1,3,5,7,...这些可能，没有3.5。
（2）我没有说换泛函，只是建议换基组。
（3）你的单重态怎么算的？波函数稳定吗？

lulu顺利毕业

zjxitcc 发表于 2023-11-23 19:57
（1）自旋多重度3.5 是什么意思？自旋多重度是2S+1，如果体系是偶数个电子，自旋多重度只有1,3,5,7,...这 ...

不好意思，回复不够严谨。
（1）分别在自旋多重度3和自旋多重度5下，进行产生稳定波函数，优化，检测稳定性这三步计算。
（2）如果我要更换的话，是否要在现有优化后稳定的结构基础上进行计算？还是在最开始的初猜结构上进行计算呢？
（3）单重态我在计算时，还不知道要检测波函数的稳定性。什么前提计算都没有，就上来直接对结构进行优化，后来发现优化后的结构，虽然没有虚频率，但是波函数检测不稳定。于是我就在单重态下，现在正在进行第一步计算产生稳定波函数。我说的能量降低是将自旋多重度3和自旋多重度5经过上述3步完成的优化，和自旋多重度1检测波函数不稳定的结构进行的对比。这样的结果比对是不是没有意义的？
（4）对于一个可能有磁性的结构，我目前只知道通过上述三步，比对在不同自旋多重度下的能量，能量越低越稳定。这种计算方法是否可靠？

zjxitcc

lulu顺利毕业 发表于 2023-11-23 20:20
不好意思，回复不够严谨。
（1）分别在自旋多重度3和自旋多重度5下，进行产生稳定波函数，优化，检测稳 ...

（1）可以。
（2）在现有的结构基础上进行。
（3）没有意义，你的单重态还在计算中，结果都不确定。
（4）从计算步骤上说是可靠的，毕竟计算步骤没错就行。从现实看，不能排除万一体系过于复杂，UBP86没法给出正确的相对能量顺序。

ahxb

lulu顺利毕业 发表于 2023-11-23 19:49
Mo12C80中空笼状团簇

你这里是scf不收敛。对scf收敛有帮助的方法很多，可以参考社长的文章http://sobereva.com/61。对你的体系，我推荐使用scf=vshift=100（再取大些也可以），并结合stable=opt在scf收敛后得到稳定波函数，必要时可尝试适当降低scf收敛精度（最低到scf=conver=6）。
也就是说，用你原来使用的泛函和基组时，关键词写为#p BP86/lanl2dz scf=vshift=100 stable=opt。另外如果你想听取其他人的建议混用基组，要注意def2tzvp对Mo是赝势基组，需要参考http://sobereva.com/60第6节的写法定义。

sobereva

根本不是波函数分析问题，发到波函数分析与Multiwfn版块做什么？不要乱找版块发帖，这次给你移动了，下次直接删帖扣分处理

lulu顺利毕业

好的，感谢各位老师的解答，已经认真研读，接下来也会尝试解决。
关于社长的回复，本人下次不会在出错，初次发帖没找好定位。
感谢各位的帮助与回复。

此外，还有一个疑问。优化后的结构，最终一定要是稳定的波函数，才能进行后续性质的研究计算吗？

zjxitcc

lulu顺利毕业 发表于 2023-11-24 09:12
好的，感谢各位老师的解答，已经认真研读，接下来也会尝试解决。
关于社长的回复，本人下次不会在出错，初 ...

是的