使用ORCA的DKH2进行过渡金属配合物几何优化出错

wangzuwei · 发表于 Post on 2024-4-26 16:44:36

尝试了一下用ORCA的DKH2结合BP86给一个过渡金属配合物做几何优化。Rh用的是SARC-DKH-TZVPP，其余原子是DKH-def2-TZVP，辅助基组SARC/J。优化过程中遇到了几个问题；
1. 基组线性相关。尝试了把sthresh设置为 1 e-6，似乎并没有解决问题。
2. 向错误的几何优化。最开始的结构是比较接近于正确结构的，Rh-Rh键应该在2.4A左右。但九轮优化过后，Rh-Rh键跟轴向上的化学键似乎都错误的向内收缩了。
3. lambda方程不收敛导致报错退出。但感觉也不是震荡，从第八轮开始就完全没有变化了。

Warning: the length of the step is outside the trust region - taking restricted step instead
The input lambda is .... -70.904285
iter: 1 x= -85.158887 g= 0.160662 f(x)= 2.290172
iter: 2 x= -105.328773 g= 0.049384 f(x)= 0.996062
iter: 3 x= -131.813012 g= 0.015895 f(x)= 0.420958
iter: 4 x= -161.324510 g= 0.005641 f(x)= 0.166472
iter: 5 x= -184.510253 g= 0.002426 f(x)= 0.056238
iter: 6 x= -193.518021 g= 0.001430 f(x)= 0.012885
iter: 7 x= -194.455684 g= 0.001192 f(x)= 0.001118
iter: 8 x= -194.464407 g= 0.001170 f(x)= 0.000010
iter: 9 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 10 x= -194.464408 g= 0.001170 f(x)= -0.000000
iter: 11 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 12 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 13 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 14 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 15 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 16 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 17 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 18 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 19 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 20 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 21 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 22 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 23 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 24 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 25 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 26 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 27 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 28 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 29 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 30 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 31 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 32 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 33 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 34 x= -194.464408 g= 0.001170 f(x)= 0.000000
-----------------------------中间省略-----------------------------------
iter: 94 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 95 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 96 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 97 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 98 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 99 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 100 x= -194.464408 g= 0.001170 f(x)= 0.000000
iter: 101 x= -194.464408 g= 0.001170 f(x)= 0.000000
[file orca_gstep/orca_gstep.cpp, line 806]: Error (ORCA_GSTEP): The lambda equations have not converged

复制代码

我把输入输出文件也上传上来了，求各位老师指点！

wzkchem5 · 发表于 Post on 2024-4-26 16:59:49

DKH有variational instability问题，一般对于平衡结构的单点计算是没问题的，但是做结构优化有时容易让键长严重偏短。几何优化一般建议用ZORA。
下一版ORCA会支持sf-X2C，比ZORA准，且没有variational instability问题

wangzuwei · 发表于 Post on 2024-4-26 17:07:38

wzkchem5 发表于 2024-4-26 16:59
DKH有variational instability问题，一般对于平衡结构的单点计算是没问题的，但是做结构优化有时容易让键长 ...

多谢解答！我用ZORA再试一下

wangzuwei · 发表于 Post on 2024-4-26 17:08:19

wzkchem5 发表于 2024-4-26 16:59
DKH有variational instability问题，一般对于平衡结构的单点计算是没问题的，但是做结构优化有时容易让键长 ...

顺便打听一下，下一版的ORCA发布有消息了吗（

wzkchem5 · 发表于 Post on 2024-4-26 17:40:54

wangzuwei 发表于 2024-4-26 10:08
顺便打听一下，下一版的ORCA发布有消息了吗（

估计还有两个月左右吧，正在做发布之前的内测

zjxitcc · 发表于 Post on 2024-4-27 00:17:07

体系非常大，建议把非金属原子的基组降为DKH-def2-SVP，Rh的基组降为SARC-DKH-TZVP

		自动登录 Automatic login	找回密码 Forget password
密码 Password			注册 Register

[ORCA] 使用ORCA的DKH2进行过渡金属配合物几何优化出错