《Chemical Science》：Recovering the causal kinetic–thermodynamic relationship in asynchronous reactionsOpen Access
编辑推荐：
动力学–热力学关系（Kinetic–Thermodynamic Relationship, KTR）是反应活性解释与设计的基础，然而强异步反应常表现出微弱、无规律甚至反向的能量响应，此类行为常被解读为异常甚至“反热力学”现象，看似违背反应能量与能垒之间预期的单
动力学–热力学关系（Kinetic–Thermodynamic Relationship, KTR）是反应活性解释与设计的基础，然而强异步反应常表现出微弱、无规律甚至反向的能量响应，此类行为常被解读为异常甚至“反热力学”现象，看似违背反应能量与能垒之间预期的单调关系。本研究表明，在异步体系中直接将活化能垒与整体反应能量相关联的标准做法，会将能量上彼此独立的阶段混并为单一表观步骤；由此产生的负梯度或任意大的梯度是表征方式造成的伪迹，而非奇特反应活性的证据。通过解耦反应阶段，可在协同异步反应中恢复其内在的能量–能垒关系，使其与同步反应直接可比。本研究通过展示初级与次级阶段能量响应的差异如何扭曲整体能量–能垒关系，对上述伪梯度进行了合理化解释。在包括1,2-氢负迁移（1,2-hydride shifts）、炔醇锂–硝酮（3 + 2）环加成以及金(I)（gold(I)）介导的6-endo-dig环化在内的代表性体系中，通过受约束内坐标（constrained internal coordinates）隔离反应阶段，恢复了具有物理意义的动力学–热力学关系。该框架表明，一旦各组成过程被恰当解析，即使是高度异步的反应也遵循明确的动力学–热力学关系。由此得到的诊断方法为实验上评估同步性以及提取各组成反应阶段的相对能量学提供了实用手段。
一、研究背景与科学问题
动力学–热力学关系（Kinetic–Thermodynamic Relationship, KTR）是物理有机化学中解释反应活性与设计新反应的核心概念，其本质在于活化能垒与反应能量之间存在的单调关联，这一关联的理论基础可追溯至Bell–Evans–Polanyi原理、Marcus理论以及Hammond假设。在这些经典框架中，过渡态（Transition State, TS）被视为反应物与产物极小值之间能量加权插值的结果，因此能垒对反应能量的响应梯度应处于0到1之间。然而，大量实验与理论研究表明，在协同异步反应中，直接以整体反应能量对活化能作图往往得到弱相关、噪声大甚至负斜率或大于1的异常梯度，这类现象曾被归因于“反热力学”行为或量子隧穿等非经典效应。Bordwell等人报道的硝基烷烃异常（nitroalkane anomaly）即为典型例证，其Br?nsted斜率（α）出现小于0或大于1的异常值。尽管已有研究从价键理论等角度对此类现象进行解释，但异步反应中动力学–热力学关系系统性失真的根源仍未得到统一阐明。
本文的核心科学问题在于：协同异步反应中，过渡态在几何与电子结构上可能与某一极小值不再保持插值关系，导致Hammond假设失效，从而使直接的能量–能垒关联沦为纯相关性而非因果性描述。因此，如何从异步反应中提取具有物理意义的动力学–热力学关系，成为亟待解决的方法学难题。
二、研究内容与主要结论
研究人员提出了一种通过反应阶段解耦恢复异步反应中动力学–热力学关系的方法框架。其核心思路是：沿内禀反应坐标（Intrinsic Reaction Coordinate, IRC）识别化学上独立的反应阶段，通过约束内坐标（constrained internal coordinates）阻止次级无垒过程的干扰，从而恢复主过程过渡态与有效极小值之间的插值关系，使Hammond假设重新成立。
研究选取了三类具有代表性的反应体系进行验证。第一类是质子化二苯乙烯的1,2-氢负迁移反应，研究发现当碳正离子高度不稳定时，邻位苯环会通过ipso亲电进攻发生无垒螺环化，形成与主过程能量响应方向相反的次级过程，导致表观动力学–热力学图出现负梯度或大于1的局部梯度。通过将Ar–CH?–阳离子角度约束为108.5°以阻止螺环化，研究人员恢复了R2=0.991的典型弯曲动力学–热力学关系，局部梯度恢复至0到1之间。第二类是炔醇锂–硝酮的（3 + 2）环加成反应，其表观动力学–热力学斜率为?0.5，呈现明显的“反相”响应。通过将烯酮CCO角冻结为180°以获取隐藏中间体（hidden intermediate）几何，研究人员恢复了初始亲核进攻阶段的单调动力学–热力学关系，局部梯度大于0.5，表明过渡态更接近烯酮而非炔醇锂，由此可指导取代基设计以降低能垒。第三类是金(I)介导的6-endo-dig环化反应，其中8个最不利的底物表现为协同异步机制，直接形成α-环丙基卡宾配合物而绕过乙烯基金(I)中间体。通过IRC分析确定各底物CCC角的最大值并用于约束优化，研究人员恢复了缺失的苄基阳离子结构，使解耦能量与同步反应数据落在同一动力学–热力学曲线上（R2=0.90）。此外，对第二步环丙烷环形成反应的动力学–热力学分析表明，其渐近线约为?2.5 kcal mol?1，解释了步级反应在特定能量范围内塌缩为协同异步过程的临界行为。
三、关键技术与方法
本研究主要采用密度泛函理论（Density Functional Theory, DFT）计算，在ORCA 6.0程序中以ωB97X-D3(BJ)泛函配合def2-TZVPPD基组进行几何优化与单点能计算，部分体系采用SMD隐式溶剂模型模拟四氢呋喃或二氯甲烷溶剂环境。反应阶段分离通过内禀反应坐标（IRC）计算实现，采用局部二次近似算法。隐藏中间体的识别依据最小能量路径（Minimum Energy Pathway, MEP）上的拐点与曲率变化。反应阶段解耦通过约束特定内坐标（如CCC角或CCO角）进行受限优化实现，以获取非平稳几何的能量参考值。
四、研究结果
2.1 异步性与表观动力学–热力学响应中的化学意义
对质子化二苯乙烯1,2-氢负迁移的研究表明，直接以整体反应能量对活化能作图得到?4.7至4.7 kcal mol?1范围内的浅曲线，放热端出现极大负梯度，吸热端局部梯度超过1，均超出动力学–热力学关系的理论允许范围。IRC分析揭示该反应为协同异步过程：氢迁移发生在过渡态附近的局域区段，而螺环化仅在更大的反应坐标位移处无垒发生。通过约束CCC角为108.5°阻止螺环化后，恢复了R2=0.991的典型弯曲动力学–热力学关系，证明初始的非单调图与大于1的梯度完全源于异步性伪迹。
2.2 协同异步与步级过程过渡区的动力学–热力学关系
对金(I)介导6-endo-dig环化的25个反应的计算表明，8个最不利底物呈现协同异步机制，直接形成α-环丙基卡宾配合物，其隐藏中间体对应乙烯基金(I)物种。从3-三氟甲基苯基到3,4,5-三氟苯基的取代变化仅使能垒改变不足0.1 kcal mol?1，即构成步级过程塌缩为异步过程的临界点。通过IRC分析确定各底物CCC角最大值并用于约束优化，恢复了全部8个缺失的苄基阳离子结构，其能量与同步反应数据拟合于同一动力学–热力学曲线（R2=0.90）。对第二步环丙烷环形成的动力学–热力学分析显示渐近线约为?2.5 kcal mol?1，解释了步级反应在特定能量范围内塌缩为协同异步过程的临界行为。
五、讨论与结论
本研究表明，协同异步反应中表观动力学–热力学关系的失真源于将能量上独立的反应阶段混并为单一表观步骤。当次级过程与主过程的能量响应方向相反时，表观梯度可超出0到1的理论范围，甚至出现负斜率，这些现象并非“反热力学”行为的证据，而是表征方式的伪迹。通过IRC分析识别反应阶段并采用约束内坐标解耦各过程，可恢复具有物理意义的动力学–热力学关系，使异步反应与同步反应直接可比。该框架同时提供了实验上评估反应同步性的诊断工具：当实验测得的局部梯度超出0到1范围时，可提示反应可能经历协同异步过程或存在未识别的额外极小值。此外，解耦后的动力学–热力学关系可揭示过渡态的电子结构特征（如更接近烯酮而非炔醇锂），为理性设计取代基以调控能垒提供直接指导。该方法的普适性不限于本文所涉及的反应类型，可推广至任意协同异步化学反应。
