溶液中配体交换反应与超分子结合复合物的色散项恰当处理方法

《Computational and Theoretical Chemistry》：How to treat dispersion terms properly for ligand exchange reactions and supramolecular binding complexes in solution

　　

伦敦色散相互作用（London dispersion interactions）在溶液体系的密度泛函理论（Density Functional Theory, DFT）计算中扮演重要角色。传统DFT-D3方法虽被广泛应用，但因忽略溶质-溶质与溶剂-溶质相互作用的复杂性及差异，倾向于高估色散作用，导致活化自由能能垒普遍偏低。因此，确定适用于溶液体系色散作用的准确处理方法具有重要意义。研究人员针对配体交换反应与超分子结合复合物，开发并验证了溶液体系色散校正的综合处理方案。对于配体交换反应，需同时考虑溶剂-溶质色散与溶质-溶质色散；而对于超分子结合复合物，溶质-溶质相互作用区域不存在溶剂-溶质色散作用。

该研究由北京师范大学化学学院Zhan Su、Dan-Yang Liu与De-Cai Fang完成，发表于《Computational and Theoretical Chemistry》。研究背景聚焦于溶液体系量子化学计算的长期挑战：传统DFT方法因未包含激发行列式，对伦敦色散作用的处理存在固有缺陷，而广泛使用的DFT-D3色散校正方法在溶液体系中常因忽略溶质-溶质与溶剂-溶质色散的差异化特征，导致结果高估，活化自由能能垒偏低，甚至与实验观测矛盾。例如，有机金属配合物配体交换反应的计算中，引入色散校正的泛函反而给出与实验相反的趋势；超分子主客体结合能也被高估。这一问题源于极化连续介质模型（Polarizable Continuum Model, PCM）仅考虑溶剂-溶质静电相互作用，且现有程序对溶剂-溶质色散能的计算未正确引入1/2因子——该因子反映了溶质分子与溶剂空腔达到热平衡时仅持有半数相互作用能的物理本质。此外，溶质-溶质色散在新形成的复合体空腔中同样仅需计入半数贡献。为解决这些问题，研究人员旨在通过验证改进的DFT-D3a方法在不同泛函中的适用性，明确溶液体系色散项的恰当处理规则。

关键技术方法方面，研究人员采用Gaussian 09与ORCA 6.0程序开展计算。选取B3LYP、PBE0、TPSS及ωb97x四种泛函，前三者搭配Grimme的D3(BJ)色散校正，ωb97x搭配Najibi与Goerigk的D3BJ校正，基组为TZP-DKH(-dfg)（氢原子去除d函数）。研究以配体交换反应与超分子结合复合物为案例，通过对比传统DFT-D3与改进DFT-D3a方法的计算结果，结合溶剂-溶质色散能（E^SSD）与溶质-溶质色散能（ΔE_D）的半量修正公式，评估不同场景下的色散处理策略。

研究结果部分，首先在计算细节中明确了泛函、色散校正方案与基组的配置。在配体交换反应（M(L₁)₂+ 2L₂→ M(L₂)₂+ 2L₁）研究中，研究人员选取五个关键反应进行分析。结果表明，传统无分散校正方法所得反应焓偏正于实验值，而标准DFT-D3方法则因高估色散作用导致偏差。采用改进的DFT-D3a方法（即同时计入溶质-溶质与溶剂-溶质色散能的1/2贡献）后，不同泛函的计算结果与实验值吻合度显著提升，证实该方法可有效解决传统方法对大位阻配体结合能的高估问题。在超分子结合复合物研究中，针对富勒烯（C₆₀、C₇₀）与buckycatcher的结合、分子镊子与小芳香分子的包结体系，传统B97-D方法预测的结合能显著过高，而DFT-D3a+IDSCRF/TZP-DKH(-dfg)计算则给出合理结果。这进一步证明，此类体系中溶质-溶质相互作用区域无溶剂-溶质色散贡献，仅需考虑溶质-溶质色散的半量修正，无需叠加溶剂-溶质色散项。

讨论与结论部分，研究人员指出，溶液体系中色散作用的处理需区分两类场景：配体交换反应中，溶剂-溶质色散与溶质-溶质色散均需以半量形式纳入能量计算；超分子结合复合物中，仅溶质-溶质色散需半量修正，溶剂-溶质色散可忽略。该规则在不同泛函中均表现出一致性，解决了传统DFT-D3方法在溶液体系中系统性高估色散作用的缺陷。这一成果为溶液相反应机理的理论研究提供了可靠的色散校正方案，尤其对催化选择性、蛋白质折叠及超分子自组装等领域的能量计算具有重要指导意义。研究同时明确了Gaussian程序中溶剂-溶质色散能计算缺乏1/2因子的局限性，为后续程序优化提供了理论依据。