计算化学公社

标题: 求助：溶剂条件下优化特戊酸根离子（-OPiv）以及-OBoc虚频无法消除 [打印本页]

作者
Author: leoliu 时间: 2025-6-17 11:51
标题: 求助：溶剂条件下优化特戊酸根离子（-OPiv）以及-OBoc虚频无法消除

(, 下载次数 Times of downloads: 31)
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也试过, 但基本上都是类似的结果。求各位大佬指导。计算资源相对这两个分子还是够多的，折腾一个多星期了，真折腾不动了。

作者
Author: leoliu 时间: 2025-6-17 11:53
图上的是OBoc，两个情况很相似。

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

作者
Author: leoliu 时间: 2025-6-17 12:24

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

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

作者
Author: wzkchem5 时间: 2025-6-17 14:33

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

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

作者
Author: leoliu 时间: 2025-6-17 14:56

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

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

作者
Author: leoliu 时间: 2025-6-18 12:12

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，网上也没找到有效的解决办法。

作者
Author: wzkchem5 时间: 2025-6-18 12:48

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的

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

作者
Author: leoliu 时间: 2025-6-19 14:42

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

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

欢迎光临计算化学公社 (http://ccc.keinsci.com/)