怎么计算反应的吉布斯自由能变？

说明：本文华算科技主要介绍反应吉布斯自由能变在第一性原理计算中的定义、单个物种自由能的组成、表面与电化学反应的参照态写法，以及模型条件改变 ΔG 时的判读范围。

反应ΔG的定义和计算来源

在 DFT文章里，反应吉布斯自由能变通常写成 ΔG，它回答的是同一温度、压力和参照态下，反应物变成产物后自由能升高还是降低。计算表达式并不神秘：ΔG = Σν_iG_i,产物 – Σν_jG_j,反应物。配平系数 ν、物种 G 和参照态共同决定 ΔG 是否可比。

DFT 给出的最低能量常是电子总能 E_DFT，它接近 0 K 下某个优化结构的势能面最低点。反应在 298 K、气相分压、溶液活度或电极电位下讨论时，E_DFT要接上零点能、热焓、熵、压力或活度修正。对一个气体分子，这些修正来自平动、转动、振动和标准态；对表面吸附物，平动和转动自由度被表面束缚，振动项更常用频率计算获得。

自由能对应一个明确的热力学状态，结构截图本身无法限定温度、压力和参照态。温度从 0 K 改到 298 K，气体从 1 atm 改到 0.1 atm，水溶液里 H⁺ 活度改变，电极电位从 0 V 推到 -0.7 V，同一个反应式都会得到不同的 ΔG。计算反应 ΔG 的第一步，是把反应式、物种状态和参照态固定下来，再去谈总能差和修正项。

单个物种的G由哪些能量项组成

常见写法是 G(T,p) = E_DFT + E_ZPE + ΔH_thermal – TS + pV。固体和吸附物里的 pV 项通常很小，气体分子在标准热化学表里会把热焓和熵项列得更完整。E_ZPE 来自振动频率的 1/2 hν，轻元素参与的反应会放大这一项，H*、OH*、OOH*、NH_x* 这类中间体尤其敏感。

气相 CO₂、H₂、H₂O 的 G 可以使用热化学表或频率计算给出的标准态修正；表面 *COOH、*OCHO、*OH 的 ZPE 与振动熵要来自该吸附构型的频率，并且结构要处在同一 slab、同一覆盖度和同一计算设置下。若反应式里有电子、质子、溶剂分子或离子，它们不能直接用真空单分子能量代替，参照态要另行定义。

在 Cu-ZIF 这类多金属位点材料中，CO₂ 可以吸附在不同 Cu 位点，局域配位、孔道限域和自旋极化都会改变 *OCHO 或 *COOH 的稳定程度。吸附位点一旦改变，进入 ΔG 的 G_ads 也随之改变，后续路径比较要使用同一套结构模型。

表面反应和电化学步骤怎样选参照态

表面反应常把干净表面写成 *，吸附物写成 *A。吸附自由能可以写成 ΔG_ads = G_*A – G_* – G_A。干净表面的 slab 能量在前后相减时抵消，剩下的是吸附物和气相或溶液物种之间的自由能差。若 slab 尺寸、真空层、k 点、磁性初态或 Hubbard U 改变，抵消关系会被破坏，ΔG_ads 会混入模型差异。

计算氢电极把电催化步骤中的 H⁺ + e^– 化学势接到 1/2 H₂，在 0 V vs. RHE、pH = 0 的标准写法下，质子电子对不用单独算真空质子和电子。外加电位进入能量时常写成 ΔG(U) = ΔG(0) – n e U，每转移一个电子，电位改变 1 V 会移动 1 eV。

HER 里常用的 ΔG_H* 就是一个参照态很清楚的例子：H⁺ + e^– + * → H* 通过 CHE 转成 1/2 H₂ + * → H*。ΔG_H* 接近 0 eV 时，H* 不会过分稳定，也不会过分难吸附；但它只描述热力学吸附强度，水解离能垒、双电层结构和局部离子分布仍会改变真实速率。

CO₂ 还原到 CH₃COOH 时，*CO、*COH、*OCCOH、*OCCH₂、*CH₃COOH 等中间体依次进入反应式。每一步 ΔG 都要按照同一电位、同一质子电子参照和同一表面模型计算，随后才能比较不同路径分支。路径上某一步上坡很大，只能说明该步热力学代价高；真正的速率由过渡态能垒控制。

哪些模型条件会改变 ΔG

低覆盖度下互相隔开的 *CO 可能稳定，高覆盖度下偶极排斥、氢键网络或 C-C 偶联几何会让同一中间体移动几百 meV。金属氧化物、单原子位点和含缺陷表面的价态、磁矩和局域电荷同样会进入数值；如果初始磁性收敛到不同电子态，ΔG 的差别就来自电子结构状态。

隐式溶剂会改变极性中间体的稳定度，显式水分子可以给 *COOH、*CHO、*CH₂OH 提供氢键，Li⁺、Na⁺、K⁺、Cs⁺ 这类阳离子会改变双电层电场和质子转移环境。强偶极、带电中间体或质子转移步骤通常要把真空 slab 的 ΔG 当作基线，再比较溶剂和离子模型带来的位移。

pH、气体分压和溶质活度通过化学势进入 ΔG。理想气体常写成 μ(T,p) = μ°(T) + k_BT ln(p/p°)，溶液物种常写成 μ = μ° + k_BT ln a。CO₂ 分压降低、CO 覆盖度升高、H₂ 活度改变，都会把对应物种的自由能推高或拉低。写结果时，温度、压力、pH、电位标尺和活度假设要和 ΔG 数值放在同一段里。

反应吉布斯自由能变的可靠度，取决于每个能量项是否指向同一个热力学状态。反应式要配平，物种 G 要来自同一计算设置，气体、溶液、表面和电极电位要使用同一参照态。若 ΔG = -0.35 eV，它表示在这些条件下产物侧自由能低 0.35 eV；换成另一种覆盖度、另一种 pH 或另一组显式水结构，这个数就要重新计算。

判读 ΔG 时，可以把结果拆成四行：反应式、物种 G 来源、电位与化学势修正、模型限定条件。四行齐全，数值才适合比较位点、路径或材料；缺其中一行，差异可能来自参照态或模型设置。