求助：溶剂条件下优化特戊酸根离子（-OPiv）以及-OBoc虚频无法消除

leoliu

sob老师的帖子我已认真地阅读。尝试过的高精度关键词如下：
1. m06/def2svp em=gd3 opt=(calcall, verytight, notrust) freq scrf=(pcm,solvent=generic, read(此处自定义为1,4-dioxane)) int=superfine
2. m06/def2svpd em=gd3 opt=(calcall, tight, notrust) freq scrf=(pcm,solvent=generic, read(此处自定义为1,4-dioxane)) int=superfine
早期的各种尝试就不展示了，精度是一步步提上来了。
在最初的优化中， -OPiv就远比其他有机负离子慢和难收敛。对于两种离子，优化基本上都会成功，能量和freq最后converge都是yes，但不管怎么尝试，都会有-60的虚频。沿虚频方向调整我也尝试了几次，但结果基本上没有区别。smd/pcm我也筛过， superfine/ultrafine也试过, 但基本上都是类似的结果。求各位大佬指导。计算资源相对这两个分子还是够多的，折腾一个多星期了，真折腾不动了。

leoliu

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

wzkchem5

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

leoliu

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

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

wzkchem5

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

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

leoliu

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

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

leoliu

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

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的

sobereva

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

leoliu

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

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