计算化学公社

标题: 水分子辅助消除找过渡态报错，不知道原因是什么 [打印本页]

作者
Author: jaywang7 时间: 2023-10-20 09:58
标题: 水分子辅助消除找过渡态报错，不知道原因是什么
一个消除反应的过渡态，用水分子辅助消除，输入和输出如下，这是为什么报错呢？
输入：
%nprocshared=24
%mem=4000MB
%chk=TS2-I-1.chk
#p opt=(calcfc,ts,noeigen) genecp m062x

TS2-I-1

0 1
C                -0.57398989 -1.50014890 -0.09352162
C                0.06067011 -1.49460190 -1.27960762
C                -0.66895789 -1.75933190 -2.57616262
C                -2.13892089 -2.06973690 -2.42125362
C                -2.84264289 -1.41377390 -1.22420662
C                -2.05382089 -1.74609890 0.03802838
H                -0.22459989 -2.59897090 -3.11295262
H                -2.58280449 -0.99224037 1.13199989
O                -2.73065289 -2.82217190 -3.14435262
O                0.08169511 -1.32125290 1.06539838
H                -0.53887389 -1.30986590 1.80072538
H                -2.23092589 -2.80123490 0.27780938
H                -0.57489789 -0.89475190 -3.24293462
C                -4.27274089 -1.93366790 -1.11661162
O                -4.61560589 -2.69398390 -0.25346662
O                -5.05309689 -1.48736690 -2.08633262
C                -6.36634289 -2.05177090 -2.12121462
H                -6.30022189 -3.13051190 -2.25857262
H                -6.86126289 -1.58430490 -2.96661562
H                -6.89377789 -1.83516790 -1.19286562
C                1.51560511 -1.23571190 -1.36235962
O                2.27468211 -1.03889590 -0.45444962
O                1.91596811 -1.24069790 -2.65603162
C                3.30495811 -0.99285390 -2.85205662
H                3.90207511 -1.75174590 -2.34623562
H                3.57617911 -0.01317990 -2.45837462
H                3.46396211 -1.03227690 -3.92625662
I                -2.85933554 1.21117347 -1.60643529
O                -2.46171771 0.16409764 0.88770596
H                -2.10286064 0.36268540 1.75513666
H                -2.02302654 1.01297990 0.12918834

I    0
S 1 1.00
   0.7242000             1.0000000
S 1 1.00
   0.4653000             1.0000000
S 1 1.00
   0.1336000             1.0000000
P 1 1.00
   1.2900000             1.0000000
P 1 1.00
   0.3180000             1.0000000
P 1 1.00
   0.1053000             1.0000000
P 1 1.00
   0.0308                1.0
D 1 1.00
   0.2940                1.0
****
O H C 0
6-311g*
****

I    0
I-ECP    3    46
f potential
  5
0    1.0715702          -0.0747621
1    44.1936028          -30.0811224
2    12.9367609          -75.3722721
2    3.1956412          -22.0563758
2    0.8589806          -1.6979585
s-f potential
  5
0 127.9202670             2.9380036
1    78.6211465          41.2471267
2    36.5146237          287.8680095
2    9.9065681          114.3758506
2    1.9420086          37.6547714
p-f potential
  5
0    13.0035304             2.2222630
1    76.0331404          39.4090831
2    24.1961684          177.4075002
2    6.4053433          77.9889462
2    1.5851786          25.7547641
d-f potential
  5
0    40.4278108             7.0524360
1    28.9084375          33.3041635
2    15.6268936          186.9453875
2    4.1442856          71.9688361
2    0.9377235             9.3630657
输出：

作者
Author: sobereva 时间: 2023-10-21 05:58
过渡态初猜明显不合理时，过渡态优化自然没法往期望的结构走，达到步数上限时没收敛就这样报错
在确保当前模型能用来描述实际反应机理的情况下，尝试修改过渡态初猜

如果溶剂环境发生的，记得带溶剂模型

作者
Author: jaywang7 时间: 2023-10-21 10:54

sobereva 发表于 2023-10-21 05:58
过渡态初猜明显不合理时，过渡态优化自然没法往期望的结构走，达到步数上限时没收敛就这样报错
在确保当前 ...

谢谢sob老师但是我看优化是往期望的结构走的但是优化了150步报错不收敛但是在view里看能量图那些都是降低并且达到平稳了这是为什么呢

作者
Author: sobereva 时间: 2023-10-22 00:07

jaywang7 发表于 2023-10-21 10:54
谢谢sob老师但是我看优化是往期望的结构走的但是优化了150步报错不收敛但是在view里看能量图那些都是 ...

明显跟你的初猜早就偏离甚巨了

作者
Author: jaywang7 时间: 2023-10-22 13:32

sobereva 发表于 2023-10-22 00:07
明显跟你的初猜早就偏离甚巨了

对于这个初猜哪里有问题 sob老师能指点下嘛我用柔性扫描也是报错

作者
Author: jaywang7 时间: 2023-10-23 11:12

sobereva 发表于 2023-10-22 00:07
明显跟你的初猜早就偏离甚巨了

谢谢 sob老师后面调整了下初猜没什么问题了

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