芳香亲核加成中迈森heimer中间体形成的计算热力学与动力学分析

《Molecules》：Computational Thermodynamic and Kinetic Analysis of Meisenheimer Intermediate Formation in Aromatic Nucleophilic Addition

　　

迈森海默（Meisenheimer）中间体为亲核芳香取代（SNAr）与脱芳香官能化提供了切入点，但在强活化芳香体系之外，仍难以确定可成键的底物-亲核试剂组合。研究人员采用密度泛函理论（DFT）对95个取代苯与吡啶底物同4种碳亲核试剂（共380种组合）的迈森海默中间体形成过程进行映射。计算所得反应自由能跨度为?90.1至+33.9 kcal/mol。亲电性与亲核性指数虽能捕捉宽泛的热力学趋势，但均不足以精确预测单个反应能量。以实验参照为中心的热力学窗口筛选出84种组合进行过渡态搜寻；定位到44个加成过渡态，计算所得加成能垒为0.4–11.4 kcal/mol。反应自由能不能对这些能垒排序，表明显式动力学分析改变了仅由热力学筛选得到的优先级。依次施加热力学窗口与“能垒不高于实验参照（8.0 kcal/mol）”的判据，确定32种组合供实验评估。结果阐明了芳香骨架、取代拓扑与亲核试剂强度如何共同控制迈森海默中间体的可及性，并为实验测试提供了一组聚焦的底物-亲核试剂配对。

研究背景方面，芳香亲核加成是亲核芳香取代（S_NAr）与脱芳香官能化的基元步骤；经典S_NAr中亲核试剂进攻缺电子芳环π体系生成阴离子σ-络合物即迈森海默（Meisenheimer）中间体，再经离去基团消除恢复芳香性。若重芳构化被抑制或延迟，该中间体可继续成键积累，其积累取决于σ-络合物形成自由能与亲核加成能垒这两个未必平行变化的量。传统迈森海默化学依赖硝基等强吸电子基稳定阴离子σ-络合物，将化学拓展至弱活化芳底物可拓宽脱芳香转化范围。Li等人近期分离出工作台稳定、酯稳定的迈森海默络合物，证明取代模式与亲核试剂匹配时弱于硝基的取代基亦可支撑持久脱芳中间体，为该研究提供实验锚点。多取代芳底物中单取代规则无法捕捉骨架、位置与亲核试剂身份的联合效应，已有计算关联σ-络合物稳定性与S_NAr位点选择性，但缺乏完整集上的热力学-动力学联立评估。该研究即在M06-2X/def2-SVP（PCM(DME)，298.15 K）水平上对380组合做反应自由能计算，再以实验参照反应（ΔG_rxn=?30.7 kcal/mol）为中心开?40~?20 kcal/mol热力学窗口筛84组合做过渡态搜寻，借亲电性ω、亲核性N及含邻位诱导项的Σσ_H做描述符分析，从而回答“哪些底物-亲核试剂配对可动力学可及地形成迈森海默中间体”。

主要关键技术方法上，研究人员构建四、五、六取代苯与三、四、五取代吡啶共95底物，配4种碳亲核试剂（1,3-二羰基稳定碳负离子Nu1、苯磺酰溴甲基碳负离子Nu2、α-氰基苄基碳负离子Nu3、烯丙基碳亲核试剂Nu4）得380组合；在M06-2X/def2-SVP几何优化与频率计算、M06-2X/def2-TZVPP单点能、PCM(DME)溶剂模型下算溶液相Gibbs自由能，双分子能加?1.89 kcal/mol校正至1 mol/L标准态；用概念DFT算底物亲电性ω与亲核性N，用文献Hammett常数并扩为Σσ_H=Σσ_m,p+n_o·σ_I处理邻位诱导；以实验酯稳定迈森海默参照反应定ΔG_rxn=?30.7 kcal/mol并设±10 kcal/mol热力学窗口选84组合做过渡态扫描，对定位鞍点验虚频与IRC。

研究结果部分，2.1节“迈森海默中间体形成中的热力学趋势”通过对380组合算ΔG_rxn（跨度?90.1~+33.9 kcal/mol，中位?6.8 kcal/mol）发现，固定亲核试剂时ΔG_rxn随底物亲电性ω增大而降低（R²=0.69–0.70，MAE 6.9–8.8 kcal/mol）；固定底物时ΔG_rxn随亲核性N增大而更负；匹配Hammett系列中引入邻位诱导项Σσ_H使20个含邻位系列R²提升0.019–0.144。结论：描述符抓主趋势但留大残差，个体判别需显式算反应能。

2.2节“迈森海默中间体形成的动力学可及性”在?40~?20 kcal/mol窗口内84组合中搜得44个一阶鞍点（成功率52%），其中36例含Nu4、29例未定位者含Nu1–Nu3；ΔG^?范围0.4–11.4 kcal/mol（中位5.1 kcal/mol）。窗口内ΔG_rxn与ΔG^?无有用线性关系。结论：热力学窗口定可居集，能垒定动力学可达性，二者互补且反应自由能不能给能垒排序。

2.3节“迈森海默中间体形成中的底物-亲核试剂兼容性”顺序施加热力学窗口与“ΔG^?≤参照8.0 kcal/mol”得32组合（24例Nu4、6例Nu2、各1例Nu1/Nu3）；例举四取代苯CF₃位+Nu4（ΔG_rxn=?30.6，ΔG^?=0.4）、四取代苯NO₂位+Nu2（?33.8，0.6）、四取代吡啶NO₂位+Nu1（?23.3，0.7）、五取代吡啶CN位+Nu2（?27.7，1.3）。结论：相似ΔG_rxn可伴不同能垒，环氮改变入窗配对，显式过渡态计算改写纯热力学筛选优先级。

讨论与结论部分，研究人员指出380组合中骨架、取代拓扑、反应位点与亲核试剂身份联合决定迈森海默中间体计算可及性；ω与N复现热力学主响应方向但个体判别需显式反应能，匹配Hammett系列证邻位短程诱导在多取代底物有可测贡献。参照中心窗口84组合中定位44过渡态、能垒0.4–11.4 kcal/mol，该子集内反应自由能不能排能垒序，故纯热力学筛选不能识别动力学优选组合；序贯用热力学窗口与参照相对能垒判据挑出32个苯/吡啶组合作实验评估集。结论译文：在380个底物-亲核试剂组合中，芳香骨架、取代拓扑、反应位点与亲核试剂身份共同决定迈森海默中间体的计算可及性；亲电性与亲核性复现主热力学响应方向，而显式反应能计算是区分个体组合之所需；匹配Hammett系列进一步表明短程邻位诱导在多取代底物中有可测贡献。参照中心热力学窗口84组合中44例定位到过渡态，其计算加成能垒跨0.4–11.4 kcal/mol；该定位子集中反应自由能不能对加成能垒排序，故仅热力学筛选无法识别动力学优选组合。施加热力学窗口与参照相对能垒判据确定32个苯基与吡啶基组合供实验评估。所选集提供具体实验靶标，并表明可及迈森海默中间体源于骨架活化、取代模式与亲核试剂强度的匹配组合。论文发表于《Molecules》。