计算化学公社

标题: 请问G16优化双配体Pd催化剂出现这类问题如何解决? [打印本页]

作者
Author:
红米饭1234    时间: 2021-10-25 09:44
标题: 请问G16优化双配体Pd催化剂出现这类问题如何解决?
请教下各位大佬,我使用BP86/DEF2SVP几何优化双配体Pd催化剂,在气相中能够正常结束且无虚频。那么我在溶剂化模型及D3BJ情况下下优化结构总是不收敛。以下是我的所有输出文件。请问还有其他办法帮助收敛吗?

作者
Author:
GoldenBaby    时间: 2021-10-25 11:08
你好像没有尝试计算精确的Hessian矩阵吧
作者
Author:
红米饭1234    时间: 2021-10-25 11:21
GoldenBaby 发表于 2021-10-25 11:08
你好像没有尝试计算精确的Hessian矩阵吧

应该如何计算精确的矩阵?
作者
Author:
sobereva    时间: 2021-10-25 11:33
上传一堆输出文件的时候,明确说清楚每个文件各对应什么,别让答疑者一个个去打开查看、去猜,不仅耽误答疑者的时间,还可能造成误解

甲苯的极性非常低,你这个体系又没有局部显著电荷,做opt freq完全没必要带溶剂模型。带了溶剂模型容易因为数值噪音使得个别情况下难收敛。加不加D3和是否能收敛没直接关系。
作者
Author:
scf    时间: 2021-10-25 11:45
qchem里有一种SWIG方法,据说收敛会好
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, introduce serious artifacts in vibrational frequency calculations, and make ab initio molecular dynamics calculations virtually impossible.523, 524 In contrast, Q-Chem’s SWIG PCMs afford potential energy surfaces that are rigorously continuous and smooth. Unlike earlier attempts to obtain smooth PCMs, the SWIG approach largely preserves the properties of the underlying integral-equation solvent models, so that solvation energies and molecular surface areas are hardly affected by the smoothing procedure.

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

orca也有类似方法
作者
Author:
红米饭1234    时间: 2021-10-25 16:23
sobereva 发表于 2021-10-25 11:33
上传一堆输出文件的时候,明确说清楚每个文件各对应什么,别让答疑者一个个去打开查看、去猜,不仅耽误答疑 ...

谢谢sob老师,我在仔细看看,因为源参考文献里面的溶剂是MTBE,想着甲苯和这个相似,因为做的Pd催化有机反应机理,想着优化双配体催化剂研究反应机理。
作者
Author:
红米饭1234    时间: 2021-10-25 16:25
scf 发表于 2021-10-25 11:45
qchem里有一种SWIG方法,据说收敛会好
https://manual.q-chem.com/5.2/Ch12.S2.html

我用高斯研究过渡态啊,不用qchem
作者
Author:
wzkchem5    时间: 2021-10-25 16:47
红米饭1234 发表于 2021-10-25 09:23
谢谢sob老师,我在仔细看看,因为源参考文献里面的溶剂是MTBE,想着甲苯和这个相似,因为做的Pd催化有机反 ...

我有点怀疑MTBE和甲苯能不能看作是足够相似的。假如能的话,那就意味着MTBE的极性也非常低,所以opt freq可以不加溶剂;假如不能的话,那就意味着你这么替代是不成立的,应当去查MTBE的介电常数等信息,然后自定义溶剂
作者
Author:
biogon    时间: 2021-10-25 17:09
wzkchem5 发表于 2021-10-25 16:47
我有点怀疑MTBE和甲苯能不能看作是足够相似的。假如能的话,那就意味着MTBE的极性也非常低,所以opt freq ...

MTBE是醚,介电常数目测和甲苯差了很大了
作者
Author:
wzkchem5    时间: 2021-10-25 17:25
biogon 发表于 2021-10-25 10:09
MTBE是醚,介电常数目测和甲苯差了很大了

我也这么觉得,我感觉用乙醚代替MTBE都比用甲苯代替合适
作者
Author:
喵星大佬    时间: 2021-10-26 11:36
wzkchem5 发表于 2021-10-25 16:47
我有点怀疑MTBE和甲苯能不能看作是足够相似的。假如能的话,那就意味着MTBE的极性也非常低,所以opt freq ...

MTBE的极性其实不小,但这也看极性怎么定义了,要是按照和硅胶的气相吸附吉布斯自由能的话这玩意快接近乙酸乙酯了,但是按照液相介电常数的话又是另一回事
作者
Author:
红米饭1234    时间: 2021-10-27 09:57
biogon 发表于 2021-10-25 17:09
MTBE是醚,介电常数目测和甲苯差了很大了

不是啊

Toluene: ε=2.3741

MTBE :Dielectric Constant @ 20°C 2.6

MTBE的介电常数,第二页
作者
Author:
红米饭1234    时间: 2021-10-27 09:58
wzkchem5 发表于 2021-10-25 16:47
我有点怀疑MTBE和甲苯能不能看作是足够相似的。假如能的话,那就意味着MTBE的极性也非常低,所以opt freq ...

MTBE的介电常数和甲苯差不多呀
作者
Author:
wzkchem5    时间: 2021-10-27 15:29
红米饭1234 发表于 2021-10-27 02:58
MTBE的介电常数和甲苯差不多呀

查了一下,确实MTBE和甲苯比和乙醚更接近。那就说明优化结构时不需要考虑溶剂效应,因为甲苯一般是不用考虑的。当然如果溶质极性特别强另说
作者
Author:
sai77    时间: 2021-10-29 09:00
我怎么感觉你这初始结构不太合理,P上两个叔丁基这么大位阻,反式配位能量更低些吧
作者
Author:
红米饭1234    时间: 2021-10-31 17:13
sai77 发表于 2021-10-29 09:00
我怎么感觉你这初始结构不太合理,P上两个叔丁基这么大位阻,反式配位能量更低些吧

啥叫反式配位啊
作者
Author:
wzkchem5    时间: 2021-10-31 18:27
红米饭1234 发表于 2021-10-31 10:13
啥叫反式配位啊

自己找金属有机化学的教材看看
作者
Author:
sai77    时间: 2021-10-31 19:08
本帖最后由 sai77 于 2021-10-31 19:09 编辑
红米饭1234 发表于 2021-10-31 17:13
啥叫反式配位啊


左边是你画的cis配位,右边trans




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