计算化学公社

标题: oniom计算出错SchOr2 failed for MMProj. [打印本页]

作者
Author: guojiayi 时间: 2017-3-6 08:56
标题: oniom计算出错SchOr2 failed for MMProj.
输出文件是这样的：
Iteration    8 internal iteration    8 dimension is    9
Root    1 not converged, maximum delta is 0.000010153172794
Eigenvalues at current iteration:
Root    1 : -0.423226723016302 Change is -0.000000000924195
DEMax= 9.24D-10
Iteration    9 internal iteration    9 dimension is 10
Root    1 has converged.
Eigenvalues at current iteration:
Root    1 : -0.423226723142998 Change is -0.000000000126696
DEMax= 1.27D-10
Convergence on expansion vectors, NOT on wavefunctions!
H-V products:    10
Lowest eigenvalue=  -0.42323
Eigenvector required to have negative eigenvalue:
QM components of vector:
                        D57    D58    D54    D55    D56
1                   0.32947 0.31501 0.31006 0.29560 0.26367
                        D53    D51    R24    D52    D50
1                   0.24425 0.23501 0.22819 0.22055 0.16920
Largest MM components of vector:
Atom 113 cartesian X -0.00479
Atom 116 cartesian X -0.00359
Atom    80 cartesian Y -0.00258
Atom 116 cartesian Y -0.00145
Atom    60 cartesian X  0.00109
Atom 113 cartesian Y -0.00097
Atom    48 cartesian Y  0.00083
Atom    48 cartesian X  0.00082
Atom    81 cartesian Y  0.00078
Atom 122 cartesian X -0.00077
H-V products w/o non-bonded:    23
SchOr2 failed for MMProj.
Error termination via Lnk1e in /public/software/Gaussian09/g09/l103.exe at Fri Mar  3 18:00:35 2017.
Job cpu time:  1 days  2 hours  4 minutes  7.6 seconds.
File lengths (MBytes):  RWF= 1375 Int=    0 D2E=    0 Chk=    28 Scr=    1

输入文件是这样的：%chk=ONI-HF-B-5-ts-1-new.1.6.2.5.chk
%nprocshared=12
%mem=20GB
#p opt=(calcfc,ts,noeigentest,quadmac) freq(noraman) oniom(b3pw91/6-31g(d,p):uff)  geom=connectivity nosymm

这个是在计算oniom，固定内层、优化外层的过程中出现的，这个怎么解决啊？谢谢大家

作者
Author: sobereva 时间: 2017-3-6 09:02
非资源帖不要发到资源分享版里，给你移动了
如新人必读里写明的，标题要反映清楚帖子内容，给你改了

作者
Author: sobereva 时间: 2017-3-6 09:04
你这不是优化极小点，而是优化过渡态。
过渡态肯定是要在内层用QM级别描述的，你既然固定了内层还怎么可能优化出内层的过渡态结构。明显ONIOM使用不对

作者
Author: guojiayi 时间: 2017-3-6 09:17

sobereva 发表于 2017-3-6 09:04
你这不是优化极小点，而是优化过渡态。
过渡态肯定是要在内层用QM级别描述的，你既然固定了内层还怎么可能 ...

嗯，好哒，谢谢老师

作者
Author: myyang 时间: 2017-4-15 21:48

sobereva 发表于 2017-3-6 09:04
你这不是优化极小点，而是优化过渡态。
过渡态肯定是要在内层用QM级别描述的，你既然固定了内层还怎么可能 ...

sob 老师，你好，我最近在做ONIOM 找过渡态的时候也出现楼主一样的错误：
Convergence on expansion vectors, NOT on wavefunctions!
H-V products:    36
Lowest eigenvalue=  -0.01099
Eigenvector required to have negative eigenvalue:
QM components of vector:
                        R33    R37    R41    R34    D219
1                   0.64559  -0.48505 0.31303  -0.11524  -0.09853
                     D222    A73    R63    A116    D91
1                -0.08886  -0.08423 0.08045 0.07885  -0.07585
Largest MM components of vector:
QM block rotation 1  0.02522
QM block rotation 2 -0.01968
QM block translation Y -0.01342
QM block translation Z  0.00955
QM block rotation 3  0.00765
QM block translation X  0.00754
H-V products w/o non-bonded:    16
SchOr2 failed for MMProj.
Error termination via Lnk1e in /opt/Chem_soft/g09/l103.exe at Sun Apr 16 00:33:31 2017
Job cpu time:  0 days  1 hours 18 minutes 55.8 seconds.
File lengths (MBytes):  RWF= 847 Int=    0 D2E=    0 Chk=    34 Scr=    1

我的关键词行：
#p  oniom(b3lyp/6-31g:amber=(softonly,qeq)) opt=(ts,calcfc,quadmac，noeigen)  freq nosymm scf(conver=6)

在我的体系中：QM 区域为发生化学反应的区域 MM区域为溶剂区域

谢谢Sober老师！

作者
Author: 维维最爱小黄人 时间: 2017-6-10 18:37
我也出现这样的错误，只要用opt=quadmac就会出现这个错误，你是怎么解决的呢？

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