文章目录
科学原理导论
反应能把分子电子结构与化学变化联系起来。在最基础的量子化学层面，电子反应能是产物与反应物电子能量按照化学计量数加权后的差。公式虽然很短，可信计算却必须同时满足多项条件：反应式元素守恒，每个分子的电荷和多重度正确，坐标来源明确，所有物种采用一致的电子结构方法，自洽场确实收敛，单位转换透明，并且报告准确说明该数值代表什么、不代表什么。
本案例使用 Psi4 1.10.2 的限制性 Hartree–Fock 方法研究固定几何构型的 H₂ + F₂ → 2 HF，并比较 STO-3G 与 6-31G。小分子让真实计算能够在 CPU 上快速完成，而含氟体系又能清楚展示基组敏感性。目标不是提出生产级热化学方案，而是建立一个可审计流程：每个分子能量、电子数、基函数数、化学计量求和、单位换算、原生日志、表格和图片都能被独立检查。
固定几何 SCF 能量差只属于电子能贡献。它没有包含结构弛豫、零点振动能、有限温度下的振动和转动项、平动贡献、压力体积功、熵、溶剂效应或标准态修正。因此结果应称为 ΔE_elec，不能自动写成焓变 ΔH、吉布斯自由能 ΔG、活化能或实验反应热。这个边界是验证合同的一部分，而不是文章末尾可以忽略的附注。
反应能公式与守恒
对于化学计量系数为 νᵢ 的反应，电子反应能写成 ΣᵢνᵢEᵢ，产物系数为正，反应物系数为负。本例为 ΔE = 2E(HF) − E(H₂) − E(F₂)。所有能量项必须使用相同方法、基组、几何约定、电荷约定和数值设置。即使每项都收敛，混用理论级别也会使差值难以解释。
计算前先检查原子守恒。反应式两侧均含两个氢原子和两个氟原子。电荷与自旋同样重要：H₂、F₂ 和 HF 都按中性闭壳层单重态表示，所以 RHF 与输入状态一致。波函数元数据提供第二层审计，每个物种的 α 与 β 电子数之和必须等于总电子数。
Psi4 以 Hartree 输出能量。流程先完成化学计量求和，再使用 1 Hartree = 627.509474 kcal/mol 转换，并同时保留两种单位。Hartree 数值便于与原生日志核对，千卡每摩尔则便于化学解释。提前四舍五入单个总能量会放大差分误差，因此计算使用全精度，展示时才格式化。
测试进度
| 验证门槛 | 结果 | 证据 |
|---|---|---|
| 技能安装 | 通过 | Skill Hub 复制完整 reaction-energy 技能 |
| Psi4 环境 | 通过 | Psi4 1.10.2 真实执行六次 RHF |
| 原子与电子审计 | 通过 | H/F 守恒且 α/β 数一致 |
| 数值验证 | 通过 | 两个确定性反应能范围匹配 |
| 实时聊天端到端 | 通过 | 智能体读数据、写代码、运行并检查产物 |
| 产物验证 | 通过 | JSON、两个 CSV、日志、PNG 和 SVG 完整 |
| 聚焦界面证据 | 必须 | 科学结果在聊天结果区渲染 |
原来的真实聊天尝试已经产生有效文件，但它早于当前语义积分字段。我们没有因为旧报告写着 passed 就直接发布，而是新增独立验证器重新读取保留的 JSON 与 CSV，核对精确数值容差、六条计算记录、电子算术、基函数数量、原生 Psi4 日志和图片。只有这个检查成功后，文章阶段才解锁。
演示用户请求
使用已安装的 reaction-energy 技能和
hydrogen-fluorination.json，通过 Psi4 1.10.2 计算 H₂ + F₂ → 2 HF 的固定几何电子反应能。验证原子守恒、电荷、多重度、电子数与基函数数；分别运行 RHF/STO-3G 和 RHF/6-31G；保存原生日志；报告 Hartree 与 kcal/mol；比较基组敏感性；输出 JSON、CSV、PNG 和 SVG。不要把结果描述为焓、吉布斯自由能、活化能或实验反应热。
这是一条自然语言科研请求，而不是写死答案的脚本。端到端测试要求智能体加载已安装技能和数据，针对任务编写程序，通过受管 Psi4 环境执行，再打开和解释结果。提示词规定科学合同，但工具选择、程序实现和产物审阅由智能体完成。
演示数据
输入文件 hydrogen-fluorination.json 保存笛卡尔坐标、分子名称、电荷、多重度、反应系数、方法、基组和换算设置。显式坐标让计算可以重现，也清楚暴露固定几何边界。这些坐标用于功能验证，不宣称是在任一基组下优化得到的最低点。
| 物种 | 化学计量角色 | 电荷 | 多重度 | 电子数 |
|---|---|---|---|---|
| H₂ | 反应物，系数 −1 | 0 | 1 | 2 |
| F₂ | 反应物，系数 −1 | 0 | 1 | 18 |
| HF | 产物，系数 +2 | 0 | 1 | 10 |
生产研究还应记录几何来源、是否优化、频率分析是否确认极小值，以及优化与最终单点是否采用相同理论级别。本测试故意固定坐标，使核心能量计算与记账路径能在无 GPU 的机器上快速完成。
安装与可复现执行
技能使用 Python 3.11 受管环境和 Psi4 1.10.2。环境边界防止系统中偶然存在的其他版本冒充测试依赖。检测会真正调用受管解释器与 Psi4，而不是只检查目录。开发者可用以下命令复现边界：
python scripts/install.py
python scripts/tool.py probe
python scripts/tool.py run -- reaction_energy_analysis.py
交互用户不需要手写这些命令或分析程序。智能体根据聊天请求编写并运行任务代码。命令块用于审计、故障诊断和独立复现，而不是把编程负担转交给用户。
物种电子能量
六次 RHF 计算全部收敛。从 STO-3G 增大到 6-31G 后，每个分子的基函数数量上升，总电子能量下降。不同分子的总能量不能直接作为稳定性排名；只有按照平衡反应式组合后才形成目标差值。
| 基组 | 物种 | 电子能量（Hartree） | 基函数数 | α / β 电子 |
|---|---|---|---|---|
| STO-3G | H₂ | −1.1166843871 | 2 | 1 / 1 |
| STO-3G | F₂ | −195.9674606982 | 10 | 9 / 9 |
| STO-3G | HF | −98.5707575917 | 6 | 5 / 5 |
| 6-31G | H₂ | −1.1267339671 | 4 | 1 / 1 |
| 6-31G | F₂ | −198.6460967419 | 18 | 9 / 9 |
| 6-31G | HF | −99.9834071596 | 11 | 5 / 5 |
独立验证器确认每一行 α 加 β 等于预期总电子数，并要求基函数数为正，而不是相信复制到表格中的文字。这些关系检查可以发现物种标签交换、错误多重度和不完整日志解析，即使错误结果看上去仍像正常数字。
结果与产物
| 基组 | ΔE（Hartree） | ΔE（kcal/mol） | 解释 |
|---|---|---|---|
| STO-3G | −0.0573700980 | −36.0003 | 固定几何下电子能差为放热 |
| 6-31G | −0.1939836102 | −121.7266 | 固定几何下电子能差为放热 |
两个结果都为负，但它们之间的巨大差异才是主要信息。最小基组 STO-3G 与分裂价层基组 6-31G 对电子分布，尤其是氟原子的孤对电子与成键区域，提供不同的表示能力。这个比较证明体系对基组敏感，却不能证明 6-31G 已经收敛或与实验一致。生产方案应继续考察极化函数、弥散函数、电子相关、几何优化、振动修正和参考数据。
截图只保留科学结果、交付物与范围警告，隐藏文件浏览器、原始 JSON 编辑器和无关面板，使读者可以直接看到两个反应能及其含义。
主图来自保留的真实计算结果，直观显示电子反应能随基组明显变化。JSON 和 CSV 是权威数据源，图片只用于解释。
辅助图根据机器可读产物生成，并在来源清单中保存源文件与 SHA-256。它不是为了凑数量而绘制的虚构科研图。
如何解释基组敏感性
基组限制了分子轨道可使用的一电子函数。STO-3G 是紧凑的最小表示，6-31G 为价层提供分裂灵活性。反应能是多个很大总能量的差，反应物与产物之间不均匀的基组误差可能显著改变较小差值。氟的成键和孤对电子密度让这个演示尤其敏感。
不能仅因为基组更大就把更负的 6-31G 结果称为“更正确”。Hartree–Fock 忽略电子相关，结构没有优化，两套基组都不是基组完备极限。系统改进需要预先设计的方法与基组层级，并在合适情况下采用相关方法。对于其他反应类别，还可能需要审查基组重叠误差、尺寸一致性和多参考特征。
基组比较还说明为什么只报告最终数字不够。研究者应同时查看每个物种的收敛状态、基函数数、电子状态与原生日志。若某个物种使用了错误多重度，或者一个计算没有真正收敛，化学计量求和仍会给出格式正确的数字。关系型验证比简单文件存在检查更能保护科学结论。
可复现性与产物审计
完整复现需要精确坐标与化学计量、Psi4 版本、方法、基组、SCF 参考与算法、能量和密度收敛阈值、内存和线程设置、单位换算常数、生成程序及原生输出。JSON 保存这些设置，species-energies.csv 和 reaction-energy-summary.csv 提供便于人工审阅的表格，原生日志证明 Psi4 确实执行。
语义验证器检查科学关系，而不只检查文件名。它核对元素守恒、数值容差、记录数量、电子算术、日志体积和图片有效性。未来如果回归版本写出空文件、交换基组标签，或者用无关的序列统计冒充量子化学，即使进程退出码为零也会失败。
可复现还不等于科学有效。完全相同的输入和软件可以稳定重现一个错误多重度或不恰当理论级别。因此审计应分成两层：先确认计算与记录可重现，再确认模型假设适合研究问题。本文提供第一层的强证据，同时明确第二层仍需要专业判断。
独立审阅时应直接从物种能量表重新计算化学计量和，而不是只相信汇总列。审阅者还应核对单位换算常数、查看每次 SCF 的最终收敛信息，并确认所有物种确实采用目标基组。这些检查成本很低，却能发现严重问题，因为反应能来自多个大数相减；一次行选择错误也可能生成看似合理的化学数值。JSON 适合精确关系检查，CSV 适合快速人工阅读，原生日志则在软件升级后数值变化时提供诊断证据。
还应检查符号约定。本文把产物系数设为正、反应物设为负，因此负 ΔE 表示在指定固定几何和电子理论下产物电子能更低。如果其他软件或论文使用相反定义，数值绝对值可以一致而符号相反。报告应同时给出反应式、公式和单位，避免读者只看到一个负数却不知道计算方向。
对于基组比较，审阅重点不只是两个最终柱形高度。应逐项检查 H₂、F₂ 与 HF 的基函数数是否按预期增加，确认电子数未随基组改变，并检查两个基组分别完成完整的三物种循环。若只重算了部分物种，混合级别差分可能仍能输出，却失去清晰物理含义。
数据治理同样影响可复现性。输入坐标、生成程序、环境清单、结果文件和图片应放在同一测试档案中，并记录摘要校验值。文章中的表格应来自机器可读文件，而不是人工再次录入。这样能够追踪每个展示数字的来源，也能让自动回归在未来版本中判断差异来自数值容差、软件行为还是文件损坏。
真实研究还需要考虑构象。小型双原子分子没有复杂构象空间，但一般反应物和产物可能存在多个低能构象。如果两侧构象选择不一致，单一固定几何差值可能主要反映构象偏差。生产流程通常需要构象搜索、相同温度下的统计处理以及明确的构象选择规则。
误差讨论也不能只写“更大基组更好”。Hartree–Fock 的相关误差与基组不完备误差可能部分抵消，也可能同向累积。合理做法是设计方法与基组矩阵，查看趋势是否稳定，并在资源允许时与高等级参考或可靠实验量比较。每一步都应保留原始数据，而不是只保留最终推荐值。
最后，文章的负反应能不能证明反应会快速发生。热力学方向与动力学速率是不同问题；即使产物电子能更低，反应仍可能具有很高势垒。要讨论速率，需要寻找过渡态并建立自由能势垒，而不能从 ΔE_elec 直接推断。
审阅者还应保存失败与警告信息。SCF 虽然可能返回最后一次迭代能量，但未收敛数值不能与正常收敛结果等价使用。自动流程应把收敛状态提升为显式字段，并在任一物种失败时阻止反应能发布。对于接近简并、强相关或键断裂体系，单一闭壳层参考可能从物理上不合适，即便数值迭代看似稳定；这类情况需要轨道稳定性、自旋状态和多参考特征诊断。
局限性
本案例没有优化结构、进行频率计算、加入零点或热修正、估计熵、模拟溶剂、定义标准态、寻找过渡态，也没有与实验比较。它不是 Psi4 精度基准，只是在小型 CPU 体系上测试平衡固定几何电子反应能的关键机制。
更完整的热化学研究应在合理理论级别优化每个物种，用频率确认极小值，计算目标温度下的零点与热贡献，应用正确标准态修正，考察相关方法和基组收敛，并量化相对参考的不确定度。对于开壳层物种，还需要选择 UHF 或 ROHF 等参考并审查自旋污染。
如果研究目标是反应动力学，则必须定位过渡态、确认唯一虚频、连接内禀反应坐标，并计算反应物与过渡态之间的自由能差。电子反应能本身不能替代活化势垒。若目标是溶液反应，还必须说明溶剂模型、浓度标准态与构象采样。
参考文献
无需编程即可尝试此工作流
MindPlot 已内置支持这种科学工作流。上传分子几何并用自然语言描述反应，智能体会编写和执行受管计算，保留原生与发布产物，验证化学记账并解释科学边界。可在 mindplot.ai 在线体验，也可下载桌面版以获得更强的本地数据隐私和环境保留能力。