3.3 保持热化学校正与标准态一致

这些输入和起始模型为未执行的教学示例。原创配图为示意图，不是计算结果。使用前请验证版本相关语法、授权数据、数值收敛及科学模型。

3.3.1 模型、单位与来源记录

结构长度使用 Å，电子能使用 hartree，振动波数使用 cm⁻¹。核对各输出字段，区分热修正与标准态项。

3.3.2 未执行输入与明确修改片段

按照下方说明确认每个片段的母计算与替换位置；片段不自动构成独立输入。保留空行与文件来源要求。

3.3.2.1 输入块 1

%oldchk=case10_water_min.chk
%chk=case11_water_298.chk
%mem=4GB
%nprocshared=2
#p B3LYP/6-31G(d) Freq=ReadFC Geom=AllCheck Guess=Read
   Temperature=298.15 Pressure=1.0

3.3.2.2 输入块 2

%oldchk=case10_water_min.chk
%chk=case11_water_350.chk
%mem=4GB
%nprocshared=2
#p B3LYP/6-31G(d) Freq=ReadFC Geom=AllCheck Guess=Read
   Temperature=350.0 Pressure=1.0

3.3.3 实作研究

3.3.3.1 问题与直观理解

同一个优化分子为什么会有电子能、零点校正能、焓和 Gibbs 自由能等不同数值？电子能只是账目的一部分，核的零点运动以及与温度有关的转动、平移、振动还会带来其他贡献，自由能又包含熵项。可以把这些量看成账本里标签明确的不同合计：每个合计已包含一些项目，将两个合计再次相加可能造成重复。本例沿用案例 10 验证后的水分子，先不引入化学反应，专注理解记账关系。

前提是 case10_water_min.chk 中确实包含同一计算水平下已完成的频率计算，不能只有 SCF 波函数，也不能只有优化过程中的近似 Hessian。重用已计算的 Hessian 是改变统计热力学分析条件，而不是重新求解电子势能面。修改温度关键字不会自动执行分子动力学、生成液体环境或优化有限温度自由能面。

3.3.3.2 检查点续算输入

下面两个输入分别运行，各自将已验证源检查点复制到新的输出检查点。使用 Geom=AllCheck 时，标题、电荷与多重度、坐标区全部省略，仅以空行结束 route 区。源文件因此保留供后续使用。

参见上方输入块 1。

参见上方输入块 2。

3.3.3.3 操作步骤

确认源极小点没有内部虚频，读取原输出中的质量、同位素和转动对称数；即使电子能不变，这些信息仍会影响热化学结果。
分别运行两次分析，核对输出中的温度单位 K 与压力单位 atm；文件名含有 350K 不代表关键字一定按预期生效。
抄录电子能、零点校正、内能/焓/Gibbs 自由能热校正及对应合计，作差前保留完整精度。多数校正和合计行使用 hartree；热化学表中的熵通常使用 cal mol−1 K−1，应看表头而非凭印象判断。
用 H = Eelec + Hcorr、G = Eelec + Gcorr 重建结果，统一单位后检查 G = H − TS。所列热校正已包含零点能，不要重复再加一次零点校正。
比较两种温度；几何和 Hessian 未变时，变化来自热力学处理。每个 G 数值都应注明温度、压力、同位素、对称性约定和缩放约定。

3.3.3.4 标准态核对

Gaussian 常规频率热化学采用理想气体的刚性转子—谐振子近似，默认条件为 298.15 K、1 atm。项目里另有溶液计算，并不会自动把平移配分函数替换成液体模型。如果所选热力学循环需要将理想气体 1 atm 标准化学势转换到 1 mol L−1 标准化学势，则每个物种的校正为 ΔG° = RT ln(C°RT/p°)，对数中的组合必须无量纲。在 298.15 K 时约为 +7.93 kJ mol−1。此数值由常数推导，是标准态转换量，不是 Gaussian 对某个分子的计算结果。

对于配平反应，总标准态转换量是 Δν RT ln(C°RT/p°)，其中 Δν 为需要转换的物种在产物侧与反应物侧化学计量系数之和的差。一分子变为一分子的异构化会抵消，缔合反应通常不会。若某溶剂化自由能已按目标标准态定义，不能机械地再加一次转换校正。应先把热力学循环各端的气相与溶液参考态写清楚，实验采用 1 bar 时也要与 1 atm 区分。

3.3.3.5 检查、错误与练习

本例的验收要求是能重建输出合计，并给出明确标准态说明。常见错误包括把 350 K 的焓与 298 K 的熵混用、把卡与焦耳混用、把热校正当成总能量，以及把分子 Gibbs 能误认为从元素标准态生成的标准生成自由能。低频扭转可能使 RRHO 熵不可靠；水分子避开了大部分此类复杂性，但不意味着近似绝对准确。

中文练习：从 μ = μ° + RT ln a 推导 1 atm→1 M 校正的正负号，再判断 A+B→AB 的总校正符号。随后说明该校正在案例 14–16 的氨翻转中为什么抵消。回答这些问题无需编造任何分子能量。