计算化学公社

标题: 加入溶剂之后优化时不收敛 [打印本页]

作者
Author: fankui1990129 时间: 2015-6-10 12:29
标题: 加入溶剂之后优化时不收敛
本帖最后由 fankui1990129 于 2015-6-10 12:44 编辑

RT，本来很容易收敛的结构加入溶剂之后完全不收敛了。四个收敛标准忽上忽下的，没有下降的趋势。请问这样正常么？有什么办法处理呢？
（本来想把输出文件放上来的，结果发现重新算的时候覆盖掉了，郁闷啊）
%chk=E:\chanwu(on-water).chk# opt b3lyp/6-311+g(d,p) scrf=(solvent=water) geom=connectivity

on-water product

0 1
C                -6.06121800 0.02404000 0.09609600
C                -5.61867400 -1.20699500 0.40062900
C                -6.48133000 -2.16871800 0.75596100
C                -7.78878200 -1.88048700 0.80319000
C                -8.21675100 -0.64526300 0.49750800
C                -7.37012400 0.34280600 0.13345300
C                -7.88971200 1.55513100 -0.16245700
C                -7.30783800 2.70281500 -0.56134700
C                -6.02246300 2.91142900 -0.76749600
N                -4.89956800 3.12104900 -0.96688700
C                -8.09396600 3.74129500 -0.78752700
N                -8.79332000 4.64059900 -1.00260900
H                -5.28849200 0.75188400 -0.18306100
H                -4.53967900 -1.43430000 0.35964200
H                -6.12134100 -3.18054900 1.00577300
H                -8.51109000 -2.66214200 1.09374700
H                -9.30608200 -0.47291800 0.55545000
H                -8.99133200 1.63341900 -0.06860200
O                -6.28450400 5.06882400 1.09385700
H                -6.68604000 5.91633100 1.17575500
H                -5.36285700 5.25290900 1.03556800

1 2 2.0 6 2.0 13 1.0
2 3 2.0 14 1.0
3 4 2.0 15 1.0
4 5 2.0 16 1.0
5 6 2.0 17 1.0
6 7 2.0
7 8 2.0 18 1.0
8 9 2.0 11 2.0
9 10 3.0
10
11 12 3.0
12
13
14
15
16
17
18
19 20 1.0 21 1.0
20
21

作者
Author: 卡开发发 时间: 2015-6-10 15:47
不知道四个标准离收敛的有多远，如果比较远的话可能得调整结构以及采用粗略点的方法/精度先找找；如果在收敛附近震荡的话opt=GDIIS值得试试，也许还有更好的方法。

作者
Author: fankui1990129 时间: 2015-6-10 16:39
第四个指标非常远，前两个指标基本上可以了。

我调整过结构，并且初步优化过，还是不行

作者
Author: 卡开发发 时间: 2015-6-10 16:51

fankui1990129 发表于 2015-6-10 16:39
第四个指标非常远，前两个指标基本上可以了。

我调整过结构，并且初步优化过，还是不行

先试试看opt=gdiis看看，也许有用。

作者
Author: sobereva 时间: 2015-6-11 00:20
参考此文，特别是文末的那几种常用套路

量子化学计算中帮助几何优化收敛的常用方法
http://sobereva.com/164

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