计算化学公社

标题: 求助:离子复合物体系隐式溶剂模型下结构优化不收敛 [打印本页]

作者
Author:
wangyueda    时间: 2022-10-17 12:51
标题: 求助:离子复合物体系隐式溶剂模型下结构优化不收敛
本帖最后由 wangyueda 于 2022-10-17 18:00 编辑

请问下各位老师,我一个阴阳离子/分子复合物体系,在气相条件下优化可以收敛,但是一在隐式溶剂模型下opt+freq时总是不收敛,参照了卢老师的相关博文http://sobereva.com/61 http://sobereva.com/164,依次尝试了在DFT方法分步结构优化(1. 去掉弥散和隐式溶剂模型 2. 1优化后结构阴离子加上弥散 3. 2优化后结构加上隐式溶剂模型跑个100来步 4. 3中受力最小的结构在calcall下计算),以及半经验计算(PM6D3),还是无法收敛,这种情况应该咋解决?(做计算以来第一次遇到。。。)附上DFT和半经验优化的结构文件各一个

作者
Author:
chands    时间: 2022-10-17 14:02
不熟悉你的体系。你做锂离子电池的离子液体?
你的(FSO2)2N-是O配位的吗?我觉得是N配位,因为N原子带负电。

作者
Author:
scf    时间: 2022-10-17 14:22
本帖最后由 scf 于 2022-10-17 14:23 编辑

如果溶剂部分里加上非极性部分,会对优化有所改进吗?
http://sobereva.com/327
"Gaussian09/16默认用的是曾经最主流的PCM模型(具体来说用的是IEFPCM第二版形式),但直接使用的话,默认并不计算非极性部分。要计算非极性部分的话,得在scrf里面写read,并且坐标末尾空一行写
Dis
Rep
Cav
这代表将构成溶剂的非极性部分的值都计算出来并纳入到输出的单点能里。"

另外https://manual.q-chem.com/5.2/Ch12.S2.html
”. Q-Chem employs a Switching/Gaussian or “SWIG” implementation of these PCMs.523, 524, 525, 373, 521 This approach resolves a long-standing—though little-publicized—problem with standard PCMs, namely, that the boundary-element methods used to discretize the solute/continuum interface may lead to discontinuities in the potential energy surface for the solute molecule. These discontinuities inhibit convergence of geometry optimizations“ 可能隐式优化在一些程序包里就是比较困难

https://www.orcasoftware.de/tutorials_orca/prop/CPCM.html
"The usual point charge scheme might lead to instabilities in the energy, e.g. if two points end up too close. In ORCA, we use a smeared Gaussian charge to obtain a smoother potential energy surface [Karplus1999] [Neese2020] instead."

http://bbs.keinsci.com/thread-25930-1-1.html
"
wzkchem5

5976
帖子        

0
威望        

3543
eV

积分
    9519

Level 6 (一方通行)

    发消息 Send PM

        
6#
发表于 Post on 2021-10-25 15:45:46 | 只看该作者 Only view this author

    scf 发表于 2021-10-25 04:45
    qchem里有一种SWIG方法,据说收敛会好
    https://manual.q-chem.com/5.2/Ch12.S2.html


orca也有类似方法"
作者
Author:
wangyueda    时间: 2022-10-17 14:53
chands 发表于 2022-10-17 14:02
不熟悉你的体系。你做锂离子电池的离子液体?
你的(FSO2)2N-是O配位的吗?我觉得是N配位,因为N原子带负电 ...

是O配位的,跑了MD算了RDF,Li-N相较于Li-O峰很小
作者
Author:
wangyueda    时间: 2022-10-17 14:56
scf 发表于 2022-10-17 14:22
如果溶剂部分里加上非极性部分,会对优化有所改进吗?
http://sobereva.com/327
"Gaussian09/16默认用的 ...

没试过,之前貌似也没看到这种非极性部分会对结构优化有啥影响的




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