锂硫电池里吉布斯自由能怎么算？

引言

锂硫电池里说吉布斯自由能不是单独的吸附能数值，最常见的误差不是公式不会写，而是把吸附能、反应能、放电电位和 Li₂S 分解势垒混成同一个指标。Li₂S₈、Li₂S₆、Li₂S₄、Li₂S₂ 和 Li₂S 对应不同硫化程度，VASP 给出的总能只是起点；真正要比较催化剂对转化的作用，必须把每一步反应的 ΔG 放在同一套参照态中。

Li-S 自由能到底算哪类量？

反应步和吸附态怎样区分？

Li-S 体系里至少有两层能量。第一层是 Li₂Sₓ 在催化剂表面的吸附稳定性，回答某个多硫化锂是否更容易被锚定；第二层是相邻硫物种之间的转化自由能，回答放电或充电过程中哪一步更难推进。首张配图同时给出有机框架模型、Li₂Sₓ 构型、结合能和态密度信息，说明一个完整判断往往要把结构、吸附和电子性质放到同一个模型中。

吸附态可以写作 Cat* + Li₂Sₓ → Cat-Li₂Sₓ*，反应态则要写成 Li₂S₈* → Li₂S₆* → Li₂S₄* → Li₂S₂* → Li₂S* 或相应的配平反应。吸附能与反应自由能的分工必须先讲清；前者的参照是清洁催化剂和孤立 Li₂Sₓ，后者的参照是相邻反应中间体、Li⁺/e⁻ 或金属 Li 等效库。两个口径可以共用同一批优化结构，但后处理表不能共用同一列公式。

图1 TM-rTCNQ 单层、多硫化锂构型、结合能和态密度示意。该图说明 Li-S 自由能分析需要先明确催化剂模型、Li₂Sₓ 对象和电子结构参照。DOI: 10.3390/molecules28052389。

为什么不能只看吸附能？

只看 E_ads 容易把“锚定强”误读成“转化快”。若 Li₂S₈* 在某表面吸附很强，但 Li₂S₂* → Li₂S* 的 ΔG 很高，放电末端仍可能出现转化迟缓；若 Li₂S* 结合太强，充电时 Li₂S 分解势垒又可能升高。吸附能是入口条件，不是完整反应路径。

同一催化剂上，长链锚定与短链转化需要分开看。长链 Li₂S₆/Li₂S₈ 需要足够吸附以降低溶解损失，短链 Li₂S₂/Li₂S 则还要能继续转化或分解。换到另一类表面时，即使所有 Li₂Sₓ 的 E_ads 更负，只要 ΔG 台阶出现更高上坡，该表面也未必更适合 Li-S 转化。Li-S 催化筛选至少要同时比较吸附强度、反应 ΔG 和关键势垒。

ΔG 由哪些项组成？

总能差从哪里取？

VASP 里先得到每个端点的电子总能 E_DFT。对表面反应，端点通常是 Cat-Li₂S₈*、Cat-Li₂S₆*、Cat-Li₂S₄*、Cat-Li₂S₂*、Cat-Li₂S*，以及必要的 Cat*、S₈、Li 金属或电化学 Li 参照。每个端点必须来自同样的泛函、赝势、ENCUT、k 点、色散校正、真空层和自旋设置。

ΔE_rxn = ΣE_DFT(products) – ΣE_DFT(reactants)　　式(1)

式(1)是反应能的骨架。若采用单个吸附态路径比较，可以直接把相邻端点相减；若采用严格配平方程，就要把所有产物和反应物的化学计量写完整。对 Li₂S₈ → Li₂S₆ → Li₂S₄ → Li₂S₂ → Li₂S 这类简化路径，同一组吸附态端点必须使用同一个表面位点、覆盖度和超胞大小，否则相邻 ΔE 会混入模型差异。若中途吸附物从金属位迁到缺陷位，应把迁移前后的端点拆成两个可说明的状态，而不是直接把能量差写进同一步转化。

修正项怎样加入？

从 ΔE 到 ΔG，需要加入零点能、熵和环境修正。固体表面和吸附态的平动/转动熵通常被抑制，常用振动频率计算得到 ZPE 与振动熵；气相 S₈ 或溶剂环境中的 Li₂Sₓ 参照更复杂，若没有显式溶剂或实验热化学数据，必须说明采用的是近似口径。

ΔG = ΔE_rxn + ΔZPE – TΔS + ΔG_solv + ΔG_U　　式(2)

式(2)里的 ΔG_solv 可来自隐式溶剂、显式溶剂团簇或文献热力学修正；ΔG_U 代表电位项，放电和充电方向符号相反。MoS₂、N 掺杂碳和 MoS₂/NCO 的配图把 Gibbs 自由能台阶与 Li₂S 分解势垒并列放置，说明低 ΔG 上坡和低分解势垒需要分开读取，二者共同决定转化和再氧化表现。单独给出 ΔG 上坡只能讨论放电转化，单独给出分解势垒只能讨论充电端动力学。若加入显式溶剂团簇，还应保持分子数、初始取向和是否参与配位的判定一致。

图2 MoS₂、N 掺杂碳和 MoS₂/NCO 模型中的多硫化锂转化 Gibbs 自由能与 Li₂S 分解势垒对照。DOI: 10.1002/advs.202201579。

电位项怎样进入放电台阶？

Li 参照怎样定义？

Li-S 放电可看作 Li⁺ + e⁻ 不断加入硫物种。若用计算氢电极类似思想处理，可把 Li⁺ + e⁻ 合并到相对于 Li/Li⁺ 的 Li 化学势中；若用金属 Li 总能作为参照，则电位 U 会改变每转移一个电子的自由能。两种写法要在同一篇计算里保持一致。

μ_Li(U) = μ_Li,metal – eU　　式(3)

式(3)说明，电位项来自 Li 电子化学势。对 n 个电子转移步骤，可写成 ΔG(U) = ΔG(0) + neU，具体符号取决于反应方向定义。实际作图时，先在 U = 0 V vs Li/Li⁺ 得到各步 ΔG，再平移到目标电位，得到对应的上坡步和限制电位。放电电位越高，加入 Li⁺/e⁻ 对反应自由能的贡献越不同。

平台和限速步怎样读取？

自由能台阶图不等同于实验充放电曲线，但可以解释平台附近哪一步存在热力学上坡。金属硫化物算例给出 Li₂S 氧化电压与相对能量曲线，旁边结构图显示 Li₂S 与表面的结合方式。该能量剖面适合说明充电端的 Li₂S 分解不能只看最终态能量，还要看沿反应坐标的能量变化。

在放电方向，若 Li₂S₄* → Li₂S₂* 或 Li₂S₂* → Li₂S* 上坡最大，这一步通常被视为热力学限制步；在充电方向，台阶图和 NEB 的分工要分清。Li₂S* 分解成 LiS* 或更高硫化物时，还要引入 NEB 端点、插值图像和过渡态频率确认，两类结果不能互相替代。

图3 不同金属硫化物表面 Li₂S 氧化电压、相对能量曲线和吸附构型。该图用于区分充电端 Li₂S 分解电位、能量剖面和局域结合方式。DOI: 10.1073/pnas.1615837114。

不同催化剂怎样保持同一计算口径？

同一反应步怎样建端点？

以 Li₂S₂* → Li₂S* 为例，端点必须保持同一表面、同一吸附区域和相近覆盖度。初态是优化后的 Li₂S₂* 吸附构型，终态是优化后的 Li₂S* 吸附构型；若反应涉及断裂或重排，需要另建 NEB 路径。若初态用 Li 端吸附，终态换成远离该位点的构型，得到的 ΔG 会同时包含反应和扩散贡献。

官能团化 MXene 算例比较了不同表面终止下的 Li₂S 分解势垒，插图显示吸附构型。端点定义必须一致。同一个 Li₂S 分解问题，在 F 终止和 O 终止表面上可能给出不同能垒；若端点本身也改变，材料差异会被建模差异放大。

图4 不同官能团化 MXene 上 Li₂S 分解能垒和吸附构型。该图说明反应端点和表面终止环境会共同影响 Li₂S 分解动力学。DOI: 10.1039/d5ra01387a。

不同材料怎样统一端点？

比较多种催化剂时，端点数据矩阵应覆盖 Cat*、Cat-Li₂S₈*、Cat-Li₂S₆*、Cat-Li₂S₄*、Cat-Li₂S₂*、Cat-Li₂S*、Li₂S 分解初末态，以及必要的 S₈ 和 Li 参照。每行对应一个材料或位点，每列对应一个物种或反应步。矩阵用于区分差异来自吸附、转化自由能，还是 Li₂S 分解势垒。

单原子催化算例用 Li-S 键 ICOHP、LiS* 吸附能和势垒建立相关关系，说明电子结构描述符可以解释为什么某些位点降低 Li₂S 氧化势垒。只有当吸附能、反应自由能、势垒和电子结构指标指向同一类位点作用时，催化解释才具备可比性。描述符不是独立结论，它应回到 ΔG 和 E_b 中验证。

图5 单原子催化 Li₂S 氧化中的 Li-S 键 ICOHP、反应能、LiS* 吸附能和势垒相关关系。该图说明描述符需要回到自由能和能垒中验证。DOI: 10.1002/advs.202102809。

VASP后处理可以怎样组织？

计算步骤先后怎样排？

第一步优化清洁催化剂表面和孤立参照物，二维或 slab 模型应检查真空层、偶极修正和底层固定策略；第二步枚举 Li₂Sₓ 初构型，包括 Li 端、S 端、桥位和边缘位；第三步对每个低能构型做静态总能，并对代表性吸附态做频率计算，得到 ZPE 和振动熵近似。

第四步把每个端点写入同一张端点数据矩阵，列出 E_DFT、ZPE、TS、溶剂近似、电位 U、ΔG(0) 和 ΔG(U)；第五步绘制放电 SRR 台阶图和充电 Li₂S 分解能垒图。频率任务只需对吸附物和邻近活性原子充分放松的构型做代表性校正，远离反应中心的底层原子可保持固定，但这个选择应与结构优化边界一致。参数可从 PBE-D3 或 optB86b-vdW、PAW 赝势、按体系收敛的 ENCUT/KPOINTS、自旋极化和 EDIFF = 10⁻⁵ 到 10⁻⁶ eV 量级开始测试，最终阈值以材料收敛结果为准。

哪些结论不能直接推出？

ΔG 更低不必然代表循环寿命更长，因为穿梭抑制、电解液溶剂化、孔结构、导电网络和真实硫载量都会影响实验表现。VASP 结果能提供的是原子尺度热力学与局域动力学解释：某个表面是否稳定 Li₂Sₓ，某个转化步是否上坡，某个 Li₂S 分解坐标是否存在较高势垒。

Li-S 吉布斯自由能计算的核心，是让每个 ΔG 台阶都能追溯到明确端点、参照态和校正项。比较结果时，应先指出同一模型口径下哪一步 ΔG 最大，再说明该步骤是否对应 Li₂S₂/Li₂S 末端转化或 Li₂S 分解，并用 Bader 电荷、差分电荷密度、PDOS 或 ICOHP 解释键合变化。若后续要和实验放电平台对应，还应说明温度、硫载量、电解液和传质因素没有完全包含在静态 DFT 模型中。

本章要点总结

要点：Li-S 中的吸附能、反应自由能、放电电位和 Li₂S 分解势垒回答不同问题，不能混作一个指标。

要点：ΔG 由 VASP 总能差、ZPE、熵、溶剂近似和电位项共同组成，所有端点应使用同一计算口径。

要点：Li 化学势决定电位项，台阶图应先在统一参照下建立，再换算到目标电位。

要点：不同催化剂比较时，要统一 Li₂Sₓ 物种、吸附构型筛选范围、覆盖度、端点定义和 NEB 反应坐标。