sobereva 发表于 2025-6-19 07:27
先用B3LYP-D3(BJ)或wB97XD做opt freq，拿无虚频的结构当M06-D3的初猜，结合很小步长上限做opt freq，并且加 ...

谢谢老师，问题这样就解决了

先用B3LYP-D3(BJ)或wB97XD做opt freq，拿无虚频的结构当M06-D3的初猜，结合很小步长上限做opt freq，并且加上CPHF=grid=fine。若还有虚频，上传输出文件

leoliu 发表于 2025-6-18 12:12
老师，这是我orca的输入文件
%maxcore     1000
%pal nprocs   4 end

从https://github.com/wzkchem5/CompoundScripts下载最新版的compound script，搭配最新版ORCA（6.1.0版），再试试。
这个compound script的很早的版本是有bug的

wzkchem5 发表于 2025-6-17 14:33
试试用ORCA的iterative optimization方法自动消虚频，再把收敛的结构放到高斯里优化

老师，这是我orca的输入文件
%maxcore     1000
%pal nprocs   4 end
* xyz  -1   1
C      1.56420700   -0.07587300   -0.00410300
O      1.54179300    1.16386800    0.00165300
O      0.30723000   -0.74777900   -0.01624100
C     -0.93289400   -0.04703500   -0.00015300
C     -1.11055900    0.81226800   -1.24723800
H     -2.12966000    1.23405800   -1.28038100
H     -0.38345300    1.63353900   -1.26128000
H     -0.96434600    0.19908100   -2.15109000
C     -1.97592100   -1.15705100   -0.00499500
H     -1.85598200   -1.80012500    0.88055000
H     -2.99640700   -0.74056800    0.00119200
H     -1.86121400   -1.78650900   -0.90094300
C     -1.09274300    0.78419700    1.26824500
H     -2.11577100    1.19233300    1.33315300
H     -0.91922400    0.15370800    2.15512100
H     -0.37526600    1.61379300    1.28480500
O      2.51207600   -0.86588300   -0.00187100
*

%compound "iterativeOptimization.cmp"
end

但是报错：
                        Compound Error
Routine        : AbortNonExistingVariable
Line (    62)  : OptDone :
Reason of error: Variable 'OptDone' is not defined before usage. Please define the variable.

Hint: In case the reported input line does not seem to be problematic,
          please check if there is a missing ';' above this line.
我是第一次使用orca，网上也没找到有效的解决办法。

wzkchem5 发表于 2025-6-17 14:33
试试用ORCA的iterative optimization方法自动消虚频，再把收敛的结构放到高斯里优化

好的，谢谢老师，我试试。

leoliu 发表于 2025-6-17 12:24
是有机金属体系，导师参考不久前JACS说用M06D3来算，不然我肯定不用这个

试试用ORCA的iterative optimization方法自动消虚频，再把收敛的结构放到高斯里优化

wzkchem5 发表于 2025-6-17 12:14
M06以外的泛函试过吗？
M06、M06-2X对格点都很敏感，格点不是特别大的话很容易出虚频。而且M06对于有机体 ...

是有机金属体系，导师参考不久前JACS说用M06D3来算，不然我肯定不用这个

M06以外的泛函试过吗？
M06、M06-2X对格点都很敏感，格点不是特别大的话很容易出虚频。而且M06对于有机体系也没有特别的精度优势，不如用B3LYP-D3

图上的是OBoc，两个情况很相似。