使用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

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

泰然汪汪汪

SMD 配合 M052X/6-31G*

风中一点白

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

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