计算化学公社

标题: 使用ORCA的DKH2进行过渡金属配合物几何优化出错 [打印本页]
作者
Author: wangzuwei    时间: 2024-4-26 16:44
标题: 使用ORCA的DKH2进行过渡金属配合物几何优化出错
尝试了一下用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方程不收敛导致报错退出。但感觉也不是震荡，从第八轮开始就完全没有变化了。
  1. Warning: the length of the step is outside the trust region - taking restricted step instead
  2. The input lambda is                     ....   -70.904285
  3.    iter:   1  x=  -85.158887  g=    0.160662 f(x)=     2.290172
  4.    iter:   2  x= -105.328773  g=    0.049384 f(x)=     0.996062
  5.    iter:   3  x= -131.813012  g=    0.015895 f(x)=     0.420958
  6.    iter:   4  x= -161.324510  g=    0.005641 f(x)=     0.166472
  7.    iter:   5  x= -184.510253  g=    0.002426 f(x)=     0.056238
  8.    iter:   6  x= -193.518021  g=    0.001430 f(x)=     0.012885
  9.    iter:   7  x= -194.455684  g=    0.001192 f(x)=     0.001118
  10.    iter:   8  x= -194.464407  g=    0.001170 f(x)=     0.000010
  11.    iter:   9  x= -194.464408  g=    0.001170 f(x)=     0.000000
  12.    iter:  10  x= -194.464408  g=    0.001170 f(x)=    -0.000000
  13.    iter:  11  x= -194.464408  g=    0.001170 f(x)=     0.000000
  14.    iter:  12  x= -194.464408  g=    0.001170 f(x)=     0.000000
  15.    iter:  13  x= -194.464408  g=    0.001170 f(x)=     0.000000
  16.    iter:  14  x= -194.464408  g=    0.001170 f(x)=     0.000000
  17.    iter:  15  x= -194.464408  g=    0.001170 f(x)=     0.000000
  18.    iter:  16  x= -194.464408  g=    0.001170 f(x)=     0.000000
  19.    iter:  17  x= -194.464408  g=    0.001170 f(x)=     0.000000
  20.    iter:  18  x= -194.464408  g=    0.001170 f(x)=     0.000000
  21.    iter:  19  x= -194.464408  g=    0.001170 f(x)=     0.000000
  22.    iter:  20  x= -194.464408  g=    0.001170 f(x)=     0.000000
  23.    iter:  21  x= -194.464408  g=    0.001170 f(x)=     0.000000
  24.    iter:  22  x= -194.464408  g=    0.001170 f(x)=     0.000000
  25.    iter:  23  x= -194.464408  g=    0.001170 f(x)=     0.000000
  26.    iter:  24  x= -194.464408  g=    0.001170 f(x)=     0.000000
  27.    iter:  25  x= -194.464408  g=    0.001170 f(x)=     0.000000
  28.    iter:  26  x= -194.464408  g=    0.001170 f(x)=     0.000000
  29.    iter:  27  x= -194.464408  g=    0.001170 f(x)=     0.000000
  30.    iter:  28  x= -194.464408  g=    0.001170 f(x)=     0.000000
  31.    iter:  29  x= -194.464408  g=    0.001170 f(x)=     0.000000
  32.    iter:  30  x= -194.464408  g=    0.001170 f(x)=     0.000000
  33.    iter:  31  x= -194.464408  g=    0.001170 f(x)=     0.000000
  34.    iter:  32  x= -194.464408  g=    0.001170 f(x)=     0.000000
  35.    iter:  33  x= -194.464408  g=    0.001170 f(x)=     0.000000
  36.    iter:  34  x= -194.464408  g=    0.001170 f(x)=     0.000000
  37.   -----------------------------中间省略-----------------------------------
  38.    iter:  94  x= -194.464408  g=    0.001170 f(x)=     0.000000
  39.    iter:  95  x= -194.464408  g=    0.001170 f(x)=     0.000000
  40.    iter:  96  x= -194.464408  g=    0.001170 f(x)=     0.000000
  41.    iter:  97  x= -194.464408  g=    0.001170 f(x)=     0.000000
  42.    iter:  98  x= -194.464408  g=    0.001170 f(x)=     0.000000
  43.    iter:  99  x= -194.464408  g=    0.001170 f(x)=     0.000000
  44.    iter: 100  x= -194.464408  g=    0.001170 f(x)=     0.000000
  45.    iter: 101  x= -194.464408  g=    0.001170 f(x)=     0.000000
  46. [file orca_gstep/orca_gstep.cpp, line 806]: Error (ORCA_GSTEP): The lambda equations have not converged
复制代码

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

作者
Author: wzkchem5    时间: 2024-4-26 16:59
DKH有variational instability问题，一般对于平衡结构的单点计算是没问题的，但是做结构优化有时容易让键长严重偏短。几何优化一般建议用ZORA。
下一版ORCA会支持sf-X2C，比ZORA准，且没有variational instability问题
作者
Author: wangzuwei    时间: 2024-4-26 17:07
wzkchem5 发表于 2024-4-26 16:59
DKH有variational instability问题，一般对于平衡结构的单点计算是没问题的，但是做结构优化有时容易让键长 ...

多谢解答！我用ZORA再试一下
作者
Author: wangzuwei    时间: 2024-4-26 17:08
wzkchem5 发表于 2024-4-26 16:59
DKH有variational instability问题，一般对于平衡结构的单点计算是没问题的，但是做结构优化有时容易让键长 ...

顺便打听一下，下一版的ORCA发布有消息了吗（
作者
Author: wzkchem5    时间: 2024-4-26 17:40
wangzuwei 发表于 2024-4-26 10:08
顺便打听一下，下一版的ORCA发布有消息了吗（

估计还有两个月左右吧，正在做发布之前的内测
作者
Author: zjxitcc    时间: 2024-4-27 00:17
体系非常大，建议把非金属原子的基组降为DKH-def2-SVP，Rh的基组降为SARC-DKH-TZVP



欢迎光临 计算化学公社 (http://bbs.keinsci.com/) Powered by Discuz! X3.3