URVAvis：基于统一反应谷方法（Unified Reaction Valley Approach, URVA）化学反应可视化机理分析的图形用户界面（Graphical User Interface, GUI）

《The Journal of Physical Chemistry A》：URVAvis: A Graphical User Interface for the Visual Mechanistic Analysis of Chemical Reactions Based on the Unified Reaction Valley Approach

　　

统一反应谷方法（Unified Reaction Valley Approach, URVA）通过对反应复合物（Reaction Complex, RC）在势能面（Potential Energy Surface, PES）上追踪的反应路径（Reaction Pathway, RP）进行详细考察，实现对化学反应的精确机理分析。在该方法中，RP的曲率被分解为与给定反应中RC内部坐标变化（如键长、键角或二面角变化）相关的曲率分量。随后对RP曲率分解进行分析，以识别反应的关键机理事件。然而，这种曲率分解分析涉及碎片化且通常耗时的计算工作流程，需要多种可视化软件。为解决这一分析局限性，本文介绍了URVAvis图形用户界面（Graphical User Interface, GUI）。URVAvis使用户能够在单一集成平台上可视化曲率分解分析所需的数据，从而简化该分析工作流程，并增强URVA在机理研究中的实际可用性。URVAvis的能力通过乙腈（CH3CN）与一氯化碳离子（CCl+）之间具有天体化学相关性的气相反应基元步骤得以展示。通过URVAvis促进的该基元步骤RP曲率分解分析，阐明了反应的基本机理事件。

该研究背景源于量子化学自20世纪80年代以来被广泛用于化学反应机理的理论研究，机理洞察来源于对势能面（Potential Energy Surface, PES）的分析，化学反应沿反应复合物（Reaction Complex, RC）的PES上的反应路径（Reaction Pathway, RP）发生转化，常包含过渡态（Transition State, TS）和中间体等驻点。传统方法仅关注RP上的驻点，导致对驻点间RC演化认识不完整，限制了机理分析深度。统一反应谷方法（Unified Reaction Valley Approach, URVA）通过关注RP的曲率，将其分解为内部坐标曲率分量，识别由局部曲率极大值和相邻极小值构成的反应相（Reaction Phase），揭示非驻点处的重要机理变化，但URVA的RP曲率分解分析需多款软件协同，流程碎片化且耗时。为此研究人员开展本研究，开发URVAvis图形用户界面（Graphical User Interface, GUI），集成RP曲率、RC几何结构与电子能量的同步交互可视化，简化分析流程，并以乙腈（CH₃CN）与一氯化碳离子（CCl⁺）的气相反应基元步骤为案例验证其能力，得出URVAvis可有效阐明反应关键机理事件的结论，对提升URVA在机理研究中的实用性具有重要意义，论文发表在《The Journal of Physical Chemistry A》。

作者为开展研究用到的主要关键技术方法包括：采用Fukui的内禀反应坐标（Intrinsic Reaction Coordinate, IRC） formalism确定RP，利用URVA2025程序计算RP各点的标量曲率、识别反应相并进行曲率分解至内部坐标分量，通过质量加权梯度向量g?(s)、质量加权Hessian矩阵f?^x(s)与RP方向向量η(s)按公式κ(s)=?1/∥g?(s)∥ (f?^x(s)η(s)?[(η(s))^?f?^x(s)η(s)]η(s))计算标量曲率，曲率分量κ_n(s)=u_n(s)^?κ(s)（u_n(s)为内部坐标向量），输入文件含XYZ格式的RC几何文件、CSV格式的RC能量文件与DAT格式的RP曲率文件；案例反应为乙腈（CH₃CN）与一氯化碳离子（CCl⁺）气相反应的基元步骤，其 stationary point单点能在DLPNO–CCSD(T)/aug-cc-pVTZ//MP2/cc-pVTZ理论水平评估并校正零点振动能，IRC计算采用Hratchian和Kraka的反应路径跟踪流程。

研究结果部分保留小标题并说明如下：

Introduction：介绍量子化学机理研究背景与PES、RP、驻点分析的局限，阐述URVA聚焦RP曲率、定义反应相的原理及曲率与曲率分量计算公式，指出URVA RP曲率分解分析的三步流程（识别显著内部坐标、同步可视化RC几何、能量与曲率、结合变化判定机理事件），以及现有分析需多软件的瓶颈，引出URVAvis的开发目的。

Methods：说明URVAvis的技术细节见于补充信息，案例为乙腈（CH₃CN）与一氯化碳离子（CCl⁺）反应；介绍该反应的天体化学背景，CH₃CN在星际分子云等被探测到，CCl⁺未探测到但其前体HCl与同类CF⁺在PKS 1830–211分子吸收体中被探测到，两者可能反应生成的C₂H₂是复杂有机分子（Complex Organic Molecules, COMs）前体，关联地外生命起源假说；回顾该反应已有机理研究，其势能面（PES）显示CH₃CN与CCl⁺经无势垒反应形成中间体3，再经异构化、氢迁移等生成多种产物，重点关注中间体5到6的基元步骤，仅比较驻点会误判仅C4N6断键与H3-C1-C4的3c-2e键形成，无法解释TS附近能量软凸起，需URVA分析；该基元步骤URVA结果显示8个反应相，通过曲率分解分析判定各相机理，如相1-2为C4N6断键形成范德华复合物，相3为H2与C1、C4形成3c-2e键生成B–C₂H₃⁺及碎片重排，相4-7为B–C₂H₃⁺与Y–C₂H₃⁺互变解释能量凸起，相8为碎片进一步取向形成中间体6，URVA揭示的瞬态事件弥补驻点分析不足。

Results and Discussion：含三个子部分，

URVAvis Input Files：说明URVAvis需三类输入文件：RC几何文件（XYZ格式，含RP各点RC几何，由IRC输出后处理脚本提取）、RC能量文件（CSV格式，URVA2025自动输出）、RP曲率文件（DAT格式，URVA2025自动输出），IRC计算采用Hratchian和Kraka流程以跟踪进出通道。

URVAvis Functionalities：介绍GUI三大功能：2D RC视觉（基于ASE库由RC几何文件生成各RP点2D动画帧，自动选最优动画平面并存储复用）；RC视觉控制面板（含滑块、指定s跳转、正反向动画、步进、首尾与过渡态（TS）定位、下载XYZ文件、调用ASE GUI查看3D结构并测量内部坐标）；曲率与能量图（交互式绘制RP标量曲率与RC能量随反应坐标s的变化，含红色示踪垂直线与RC视觉联动，动画时禁用示踪防崩溃可切换）。

URVA RP Curvature Decomposition Analysis with URVAvis：以反应相3为例演示分析，第一步从曲率分解（图4b）识别相3最显著的5个内部坐标：H2C4、C4N6、∠H2C1C4、∠C1C4N6、∠H3C1C4；第二步用URVAvis同步可视化相3起止点（图8）与动画，通过ASE GUI量化变化：H2C4从1.6 ?降至1.4 ?，∠H2C1C4从86.2°降至67.2°，C4N6从2.3 ?升至2.6 ?，∠H3C1C4从152.4°升至173.1°，∠C1C4N6从121.2°微升至121.7°，能量近乎平台仅微升数kcal/mol；第三步判定机理事件：一是C₂H₃⁺碎片的H2转向C4形成3c-2e键生成B–C₂H₃⁺，对应H2C4、∠H2C1C4减小与∠H3C1C4趋近180°；二是C₂H₃⁺与NCCl碎片的初始相互远离（C4N6增大）与重取向（∠C1C4N6微增），需克服离子-偶极与偶极-偶极吸引致能量微升，两事件为后续相4-7异构与终态构型奠基；URVAvis将多软件碎片化分析整合为单一平台，简化流程并阐明相3近能量平台成因与关键事件。

总结讨论部分：URVAvis是基于URVA的化学反应可视化机理分析GUI，旨在简化URVA2025输出数据的RP曲率分解分析，消除多软件需求，通过交互同步可视化RP曲率、RC能量与几何随反应坐标s的变化，简化反应相研究与机理事件判定，促进URVA机理研究；未来版本将增加非能量标量RC属性（内部坐标、自然键轨道（Natural Bond Orbital, NBO）电荷、偶极矩等）随s变化图，以及曲率图的反应相 demarcation线与曲率分量曲线，使其可生成可直接交流的能垒与曲率图用于报告与出版，扩展实用价值。

翻译研究结论部分：URVAvis，一种基于URVA用于化学反应可视化机理分析的GUI，在本文中被呈现。该GUI设计用于简化URVA2025程序生成的输出数据的URVA RP曲率分解分析，消除了在分析每个URVA反应相曲率分解时使用多种分子与数据可视化软件的需求。该程序对RP曲率、RC能量与RC几何结构随反应坐标的交互式同步可视化，为研究反应相提供了一种简化的、用户友好的方法。这种简化反过来使用户能够更容易地判定每个相的机理事件，从而促进基于URVA的反应机理研究。在URVAvis的未来版本中，将包含两个新的GUI功能集。其中第一个将是呈现各种非能量标量RC属性随反应坐标变化的图（例如内部坐标值、源自NBO分析的原子电荷、偶极矩等）。将这些属性的图与当前URVAvis的能量和曲率图一起包含在内，将使研究人员能够从URVA RP曲率分解分析揭示的机理事件中获得更详细的见解。将添加到URVAvis的第二组功能涉及GUI的 energy 和曲率图。与图4中的对应图不同，这些图目前不包括划分反应相的线，并且在曲率图的情况下也不包括曲率分量曲线。然而，这些功能将包含在未来的GUI更新中。一旦两个图的反应相划分线和曲率图的曲率分解曲线被包含，这些图将类似于图4的图。因此，这些添加将把URVAvis的实用性扩展到促进URVA RP曲率分解分析之外，并允许用户也容易地生成可以容易地传播到研究社区的能量和曲率图（例如通过演示或出版物）。