甲基化增强多甲基咪唑与CO2的π–π堆积作用

《Journal of Molecular Spectroscopy》：Methylation strengthens π–π CO2 stacking with polymethylated imidazoles

　　

研究人员采用MP2、CCSD及CCSD(T)等高精度ab initio波函数方法，系统研究了咪唑及其甲基衍生物与CO2形成的复合物，旨在探索多甲基化咪唑作为环境友好型高效CO2捕集体系的可行性。研究发现，体系中最稳定的两种构型分别为N···C四配位键（tetrel bond）构型，以及CO2位于咪唑环π云上方的构型。与未取代的咪唑-CO2参考体系相比，甲基取代显著调控了这些能量极小点的相对稳定性。多甲基化不仅提升了结合能，还能够通过增强π–π相互作用克服优势结合位点产生的空间位阻效应，这一机制得到了对称适配微扰理论（SAPT-DFT）的能量分解结果支持。此外，研究提供了微波波段的光谱参数信息，并应用机器学习方法对体系的势能面进行了探索与建模。

该研究发表于《Journal of Molecular Spectroscopy》，针对咪唑及其甲基衍生物与CO₂的弱相互作用机制展开深入探讨。研究背景聚焦于当前CO₂等温室气体捕集技术的迫切需求，以及咪唑类杂环分子在生物医药、天体化学及材料科学中的重要地位。现有研究表明，咪唑可通过四配位键（tetrel bond）、氢键及π–π堆积等多种非共价作用结合小分子，但在气相中甲基取代如何调控其与CO₂的结合模式、尤其是π–π相互作用强度的变化规律尚不明确，缺乏系统的理论解析与光谱预测数据。为此，研究人员选取咪唑、单甲基咪唑、1,4-二甲基咪唑及1,2,4,5-四甲基咪唑为模型体系，结合高精度量子化学计算与机器学习方法，阐明甲基化对结合位点偏好、相互作用本质及微波光谱特征的影响，为设计新型气体吸附材料提供理论依据。

关键技术方法方面，研究人员采用分层优化策略：首先利用MP2/6–311++G(d,p)水平优化单体与复合物几何结构并验证驻点性质，随后在CCSD与CCSD(T)水平进行单点能与结构精修；通过对称适配微扰理论（SAPT-DFT）在PBE0/aug-cc-pVTZ水平分解相互作用能，量化静电、交换、诱导与色散贡献；基于量子理论中的原子（QTAIM）分析电子密度拓扑性质，结合约化密度梯度（NCI）指数可视化弱相互作用区域；最后构建包含3348个构象的数据集，利用高斯过程回归（GPR）等监督机器学习算法建立势能面预测模型，并通过五折交叉验证评估模型性能。

研究结果部分，在3.1节对咪唑-CO₂复合物的基准研究中，研究人员确认了四个能量极小点：四配位键（TB）、π–π堆积、NH···O氢键及CH···O弱氢键构型。CCSD(T)/6–311++G(d,p)计算显示，TB构型结合能为?20.1 kJ·mol^?1，略高于π–π构型的?16.8 kJ·mol^?1，两者差距为3.3 kJ·mol^?1；QTAIM分析表明所有相互作用均为闭壳层弱作用，键临界点电子密度介于0.006–0.012 au之间。

在3.2节单甲基咪唑复合物的研究中，甲基化使结合能整体提升，且TB与π–π构型的能量差显著缩小至0.7–2.7 kJ·mol^?1，其中1-甲基咪唑的差距仅为0.7 kJ·mol^?1，表明π–π作用已具备与TB竞争的能力。此外，4-甲基与5-甲基咪唑为互变异构体，前者孤立态能量低2.2 kJ·mol^?1，形成TB复合物后差距扩大至3.0 kJ·mol^?1。研究人员还计算了各复合物的转动常数（A、B、C）与Ray不对称参数（k），发现CO₂结合导致转动常数显著降低，且甲基取代位置对光谱参数具有灵敏调控作用，可为实验光谱指认提供依据。

在3.3节多甲基化效应的研究中，1,4-二甲基咪唑的π–π构型结合能首次超过TB构型（?23.0 vs ?22.1 kJ·mol^?1）；而在1,2,4,5-四甲基咪唑体系中，π–π构型结合能达到?28.3 kJ·mol^?1，显著高于TB构型的?24.3 kJ·mol^?1。SAPT能量分解显示，多甲基化大幅增强了色散作用贡献，四甲基体系的π–π构型色散能达?23.6 kJ·mol^?1，占总吸引能的60%以上；NCI分析进一步证实，甲基化扩大了π体系的电子云分布，强化了堆叠区域的相互作用密度。

在机器学习建模部分，研究人员对比了随机森林、支持向量回归及高斯过程回归（GPR）等七种算法，最终确定Rational Quadratic核的GPR模型表现最优，在全能量区间测试集的均方根偏差（RMSD）为2.37 kJ·mol^?1，当限定能量低于0 kJ·mol^?1时RMSD降至0.87 kJ·mol^?1，实现了对复合物势能面的高精度预测。

讨论与结论部分指出，甲基化通过增强色散作用与扩大接触界面，成功逆转了咪唑-CO₂体系中TB与π–π作用的稳定性顺序，突破了传统认知中氮位点主导结合的局限。这一发现不仅解释了多甲基化材料在CO₂吸附中的优异性能，也为通过取代基调控非共价作用提供了新思路。研究同时证实，MP2与CCSD优化的几何结构与CCSD(T)结果高度一致，可采用CCSD(T)//MP2复合方案在保证精度的前提下降低计算成本。此外，预测的转动光谱参数可指导后续的微波实验探测，为天体化学中咪唑类分子的星际介质搜寻提供理论支撑。