计算化学公社

标题: 使用Gaussian计算过渡金属溶剂化能的问题 [打印本页]

作者
Author: 风中一点白 时间: 2023-12-14 10:17
标题: 使用Gaussian计算过渡金属溶剂化能的问题

本人根据Sob老师《谈谈隐式溶剂模型下溶解自由能和体系自由能的计算》（http://sobereva.com/327）文章第2节第（8）小节的描述。

计算过渡金属离子溶剂化能的时候，因为单个过渡金属离子的电荷密度太高，先用配体与金属离子配位得到<span]计算过渡金属离子溶剂化能的时候，因为单个过渡金属离子的电荷密度太高，先用配体与金属离子配位得到<span]优化后，用添加SMD隐式溶剂化的金属配位团簇的能量减去未添加SMD的能量。泛函使用b3lyp/6-311+g**，隐式溶剂化方法scrf=(smd,solvent=generic,read)

计算结果表明，能量值和金属离子价态强相关，例如：2价离子6eV左右，3价离子11.9eV左右。该能量值与金属种类关系不大。

例如：用配体A（涉及保密不方便透露）分别对Fe(II)、Fe(III)、Mn(II)、Mn(III)、Ni(II)、Cr(III)进行6配位。

先结构优化，再用b3lyp/6-311+g** scrf=(smd,solvent=generic,read)计算隐式溶剂化能量，再用b3lyp/6-311+g**计算无溶剂化时能量，两者能量相减，得到溶剂化能。

Fe(II)：6.22eV Mn(II)：6.07eV Ni(II)：6.29eV

Fe(III)：11.9eV Mn(III)：11.86eV Cr(III)：11.98eV

想问一下，这样计算溶解自由能的方法是否合理？

作者
Author: 泰然汪汪汪 时间: 2023-12-14 14:09
SMD 配合 M052X/6-31G*

作者
Author: 风中一点白 时间: 2023-12-14 15:59

泰然汪汪汪发表于 2023-12-14 14:09
SMD 配合 M052X/6-31G*

我先试试修改泛函计算。
我疑惑的是溶剂化能数值与金属离子价态强相关，与种类关系不大，这个结果是否合理？

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