量子化学微溶剂化通过基于自由能的溶剂化位点识别

《Journal of Chemical Theory and Computation》：Quantum Chemical Microsolvation via Free Energy–Based Solvation Site Identification

　　

在量子化学对溶液的描述里，溶剂可作为隐式介质只计入本体性质，也可显式地把溶剂分子纳入计算。显式处理虽更准确，但计算代价常过高。作为折中，微溶剂化（microsolvation）模型只加入有限数目溶剂分子以捕捉最关键溶质?溶剂相互作用，但如何确定最优溶剂分子数目、位置与取向仍是难题。现有方法要么按次序逐步搭建溶剂簇，要么仅按几何准则挑选，未直接量化相互作用强度。研究人员提出 Free Energy Based Identification of Solvation Sites 2（FEBISS 2）：一种基于 Grid Inhomogeneous Solvation Theory（GIST，网格非均匀溶剂化理论）对分子动力学（MD）轨迹分析的改进微溶剂化协议。该方法从自由能（Free Energy）贡献中识别热力学有利溶剂化位点，自动放置、取向并依相互作用强度排序溶剂分子；其取向采用统计平均、与放置次序无关，且适用于任意刚性溶剂而非仅水。该方法为微溶剂化簇中溶剂分子数目的选取提供物理依据。研究人员在小有机分子、过渡金属催化剂及多种溶剂体系中展示其适用性，并表明显式微溶剂化即便在传统认为弱溶剂敏感体系中也会改变电子结构，从而凸显 FEBISS 2 可作为量子化学研究里构建物理有意义微溶剂化簇的实用方法。

研究背景方面，许多化学过程在溶液中进行，溶剂在微观层面可直接影响分子相互作用甚至参与反应。量子化学中常用极化连续介质模型（PCM）或导体样屏蔽模型（COSMO）等隐式溶剂化方法，成本低但难以描述强方向性相互作用如氢键；完全显式溶剂又需昂贵的第一性原理分子动力学或量子力学/分子力学（QM/MM）方法。微溶剂化或簇?连续介质方法只把最重要的溶剂分子作显式处理，其余作隐式处理，但“选哪些溶剂分子、放哪里、怎么取向”缺乏基于相互作用强度的客观标准。自下而上方法逐步加分子、存在顺序依赖；自上而下方法从 MD 轨迹取快照，却多按几何距离筛选，溶质?溶剂相互作用强度只在选完之后评估。因此研究人员开展本研究，发展 FEBISS 2，用 GIST 从 MD 轨迹得到溶质周围溶剂化自由能空间分布，并以此驱动溶剂分子放置、取向与排序。

该研究发表在《Journal of Chemical Theory and Computation》。研究人员得出：FEBISS 2 能对任意刚性溶剂给出统计平均取向，可用质心（COM）或指定中心原子放置，配套预参数化十种常用溶剂，为有机与有机金属体系反应机理研究中的微溶剂化设置提供通用框架；显式微溶剂化会改变被认为弱溶剂敏感体系的电子结构。

关键技术方法上，研究人员先做受限溶质经典 MD 模拟；再用 CPPTRAJ 做 GIST 分析，把模拟子体积剖分为 0.5 ? 边长体素，得到每体素的溶质?溶剂能量、熵与自由能及溶剂密度；随后按体素占据率累计放置溶剂中心点或 COM，以占据率加权求溶剂分子自由能；再用四元数处理对称等价取向、由 MD 各帧平均出溶剂取向；最后生成微溶剂化簇并可做 QM/MM 或全量子优化。样本来自文献体系：尿素、(+)-borneol、Brooker 部花青染料、钼亚胺烷基卡宾 Schrock 型烯烃复分解催化剂。

研究结果如下。

Grid Inhomogeneous Solvation Theory：研究人员简述 IST 与 GIST，自由能 ΔG_solv = ΔE_solv ? TΔS_solv，并分解为溶质?溶剂与溶剂?溶剂重组项；GIST 用离散体素近似空间积分，输出每体素热力学量，但要求溶质构象刚性或逐构象分析。

Computational Protocol：研究人员给出流程，含构象采样（可选）、受限 MD、GIST 分析、溶剂放置与自由能赋值、排序选择、平均取向确定、微溶剂化簇生成及后处理。

Classical MD Simulation of the Restrained Solute：研究人员将溶质约束在原点做经典 MD，每运行只能分析一个构象；柔性溶质需先采构象并逐构象跑受限 MD。

GIST Analysis：研究人员在约 30×30×30 ?³ 子体积、60×60×60 体素上算 ΔG_solv,k、溶剂密度与取向；以 P₀ 为中心点，P₁、P₂ 定取向，体素占据率 θ_k 由 P₀ 出现频率得。

Solvent Placement and Free Energy Assignment：研究人员把 P₀ 放最高占据体素中心，向邻壳累计至总占据为 1，自由能按 θ_k 加权求和，随后把已用体素占据率置零避免重复放置。

Ranking and Selection of Solvents：研究人员按自由能排序列出溶剂分子，用交互条形图与溶质?溶剂最近原子距离配色辅助用户选分子数目，距离阈值仅作粗指南。

Determination of Average Solvent Orientations：研究人员用三点定义溶剂取向，对氯仿等 C_3v 对称体系用参考四元数 q_ref 与对称等价四元数映射，再按 Markley 平均法得统计平均取向。

Generation of Microsolvated Clusters：研究人员输出含溶质与选定溶剂的 PDB 簇结构，供后续量子化学计算。

Post-Processing：研究人员指出平均取向簇可经 QM/MM 或 GFN2-xTB、DFT、波函数理论全量子优化；要捕获超出静电的显式溶剂效应须对全簇做量子计算。

Urea in Water：研究人员对尿素?水体系用 COM 与氧原子放置均得到冠状第一水壳，最强水分子 ?12.3 kcal mol^?1，位于两个酰胺氢之间形成分叉氢键；两种放置绝对作用差小于 0.5 kcal mol^?1，位置高度相似，并与 3D-RISM、从头算 MD、自下而上微溶剂化结果一致。

(+)-Borneol：研究人员在 DMSO 中得到仅一个强相互作用 DMSO，其氧与 (+)-borneol 羟基形成氢键；距离着色显示近溶质分子多但自由能跳跃少，说明几何近不等于热力学强。

Brooker’s merocyanine dye：研究人员在甲醇、氯仿、DMSO、四氯化碳中研究该染料，经隐式 C-PCM 下 DFT 优化微溶剂化簇，显示不同溶剂环境对染料结构与溶剂组织方式影响可被 FEBISS 2 捕获。

Molybdenum imido alkylidene N-heterocyclic carbene catalyst：研究人员对二氯甲烷中的 Schrock 型钼烯烃复分解催化剂做构象分析后受限 MD 与 GIST，生成微溶剂化簇并 QM/MM 优化；全簇 DFT 单点算 Mayer 键级，表明显式二氯甲烷改变催化剂关键键级，说明均相催化中显式溶剂不可忽略。

讨论与结论部分，研究人员总结：旧 FEBISS 靠元素密度定取向，对多位点溶剂不唯一；FEBISS 2 用溶剂子结构、四元数与对称处理，适用任意刚性溶剂。微溶剂化数目应由自由能曲线跳跃与近溶质强相互作用决定，而非单纯加更多水分子。显式微溶剂化不仅影响氢键网络，也重构电子结构，即便是传统弱溶剂敏感体系亦然。研究人员结论为：FEBISS 2 是以热力学量为依据的、与放置次序无关的微溶剂化构建协议；它把 MD 得到的溶液相信息转化为可量化、可复现、物理合理的显式溶剂簇，为有机、有机金属与均相催化机理研究中的溶剂效应建模提供实用且可推广的方法论。