Li-S电池里多硫化锂吸附能怎样算？

引言

Li-S 电池计算里，多硫化锂吸附能经常被用来评价催化剂是否能抑制穿梭、促进转化。真正难的不是写出 E_ads，而是保证不同材料、不同位点和不同 Li₂Sₓ 物种之间有相同口径。Li₂S₄、Li₂S₆ 和 Li₂S₈ 的尺寸、极性和可弯曲程度都不同；同一个分子用 Li 端、S 端或桥位靠近表面，优化后也可能得到完全不同的接触方式。多硫化锂吸附能只有在物种、构型集合、表面面积和参照态一致时，才有横向比较意义。

要让这个数值真正可比，计算时应先固定 Li₂Sₓ 对象，再枚举合理构型，随后用同一套 slab 和同一套能量参照求吸附能。这样得到的结果才适合进入催化剂筛选、差分电荷密度、Bader 电荷或自由能台阶分析。

图1 多硫化物-正极结合能景观。图中展示不同表面对 Li-S 物种结合强度的差异，说明吸附能需要放在统一物种和统一表面口径下比较。DOI: 10.1002/advs.202206057。

先确定比较的是哪一个 Li₂Sₓ

Li-S 转化链条中常见对象包括 S₈、Li₂S₈、Li₂S₆、Li₂S₄、Li₂S₂ 和 Li₂S。它们不是同一个吸附物的不同名字，而是电化学转化过程中电子数、Li/S 比和几何尺寸逐步变化的物种。把 Li₂S₄ 在材料 A 上的吸附能，与 Li₂S₆ 在材料 B 上的吸附能直接比较，不能回答哪个表面更能锚定多硫化锂。

同一篇比较中，建议先选定一到两个核心物种。若研究穿梭抑制，Li₂S₆ 或 Li₂S₄ 常更能代表可溶性多硫化物；若研究沉积和分解，Li₂S₂ 或 Li₂S 更接近固相终态。吸附能矩阵应按“同一 Li₂Sₓ 物种 × 不同表面/位点”组织，而不是把整条转化链混成一个排序。

气相或孤立 Li₂Sₓ 参照也要一致。孤立分子需要放在足够大的盒子里优化，保持相同泛函、赝势、截断能、展宽和自旋设置；若分子存在多个低能构型，应固定最低能或明确选择规则。否则表面吸附能会把分子参照误差一起带进来。

E_ads(Li₂Sₓ) = E_{slab+Li₂Sₓ} – E_slab – E_Li₂Sₓ 式(1)

式(1)中的三个能量必须来自同一计算口径。E_{slab+Li₂Sₓ} 是吸附后优化结构总能，E_slab 是同一表面裸 slab 总能，E_Li₂Sₓ 是孤立多硫化锂分子总能。若 E_slab 来自不同超胞或不同固定层，吸附能就不再只反映吸附强弱。

图2 Li₂S 生长与多硫化物结合能之间的关系。图中强调吸附强弱会影响后续成核和沉积过程，不能只看单个末态能量。DOI: 10.1002/advs.202206057。

构型枚举怎样才不偏向某个材料

多硫化锂不是小分子 CO 或 H*，它有多个 Li 原子和 S-S 链段，能以 Li 端、S 端、平躺、侧立、桥连等方式接触表面。极性氧化物、硫化物、氮化物、碳基单原子位点给出的低能构型往往不同。若材料 A 枚举了 12 个构型，材料 B 只放了一个 Li 端吸附构型，最后排序会明显偏向枚举更充分的一侧。

比较不同催化剂时，应给每个材料相同层级的初始构型集合。例如对 Li₂S₄，可以至少考虑 Li-Li 双端接触、末端 S 接触、Li-S 桥连、链段平躺和靠近金属位点的弯折构型。每个构型都应充分优化，再取最低能构型进入横向比较。公平比较不是固定一个初始姿势，而是固定同样的构型搜索规则。

优化后还要检查结构身份。若 Li₂S₆ 在某表面断裂成 Li₂S₄ + S₂，或者 Li 原子转移到表面缺陷位，此时它不再是完整分子吸附能，而接近解离吸附或表面转化反应能。这样的结果不能和完整 Li₂S₆ 吸附能放在同一列，需要单独命名和讨论。

图3 高指数晶面纳米晶在 Li-S 电池中的双功能催化行为。图中把吸附、催化转化和电化学表现联系起来，说明单一吸附能需要与转化证据共同判断。DOI: 10.1007/s40820-021-00769-2。

表面面积和覆盖度为什么要统一

Li₂Sₓ 分子体积较大，小超胞容易产生周期性镜像相互作用。若一个材料用 3×3 表面超胞，另一个材料用 5×5 超胞，单个 Li₂S₄ 的覆盖度已经不同。吸附能差异可能来自分子-分子镜像排斥，而不是材料本身的锚定能力。

更可靠的做法是让每个模型满足两个条件：吸附物之间的周期性距离足够大，且比较组中的表面面积尽量接近。对二维材料或碳基载体，还要保证真空层、层数和固定策略一致；对有缺陷或单原子位点的表面，需要说明一个超胞中有几个活性中心。多硫化锂吸附能比较要同时标注表面超胞、活性位数量和等效覆盖度。

如果研究高覆盖或多分子协同，能量口径要从“单分子吸附能”切换到增量吸附能。例如先放一个 Li₂S₄，再放第二个 Li₂S₄，第二个分子的吸附能应以已经吸附一个分子的表面为参照。这可以区分孤立锚定能力和拥挤条件下的多分子容纳能力。

E_ads,2 = E_{slab+2Li₂Sₓ} – E_{slab+Li₂Sₓ} – E_Li₂Sₓ 式(2)

式(2)是第二个 Li₂Sₓ 的增量吸附能，适合讨论覆盖度升高后的吸附变化。它不能直接和式(1)混在一起排序，因为参照表面已经不同。若需要比较材料筛选趋势，应优先使用同一覆盖度下的同类能量。

图4 金属氧化物上多硫化物吸附与电催化转化示意。图中体现了表面化学组成、吸附构型和转化动力学之间的联系。DOI: 10.1002/advs.202204930。

吸附过强和过弱怎样判断

Li-S 催化不是吸附越强越好。吸附太弱，可溶性多硫化物容易离开正极区域；吸附太强，Li₂Sₓ 或 Li₂S 可能难以继续转化，导致沉积和分解过程受阻。因此吸附能应与 Li-S 键长、S-S 键变化、差分电荷密度、Bader 电荷和后续反应能垒一起判断。

一个实用判断是把吸附能分成三类：弱吸附说明锚定不足；中等吸附说明表面可稳定多硫化物并保留转化空间；过强吸附则需要继续看 Li₂S 分解、S-S 断裂或 Li+ 迁移是否困难。吸附能给出锚定强弱，不等于直接给出 Li-S 催化活性排序。

对不同材料的筛选，最好同时报告最低能构型和关键几何量，例如 Li-O、Li-N、Li-S、M-S 距离，以及 S-S 链是否被拉长。若吸附能更负且 S-S 键适度活化，才更适合作为转化促进证据；若只是 Li 原子被表面强烈固定，后续转化可能反而受限。

图5 金属硫化物表面 Li₂S 氧化的计算模型与能量关系。吸附后的终态还需要继续连接 Li₂S 分解或氧化步骤，才能判断转化动力学。DOI: 10.1073/pnas.1615837114。

VASP任务怎样安排

建议把任务分成四层。第一层优化裸 slab 和孤立 Li₂Sₓ 分子，得到统一参照。第二层对每个表面生成相同构型集合，优化后筛出最低能构型。第三层做静态单点、差分电荷密度、Bader 或 PDOS，用来解释吸附来源。第四层选择关键体系继续计算 Li₂S 成核、Li₂S 分解或多硫化物转化能垒。

INCAR 口径要保持一致。吸附优化通常需要 ISPIN = 2，过渡金属体系要检查磁矩初值；范德华修正对多硫化锂和碳基表面尤其重要；大分子吸附要控制 EDIFFG、真空层和偶极修正。若不同材料的收敛设置不同，应说明这是收敛测试后的必要差异，而不是随意更改。

SYSTEM = Li2S4_adsorption
ISPIN = 2
ENCUT = 520
EDIFF = 1E-5
EDIFFG = -0.03
IBRION = 2
NSW = 250
IVDW = 12

# 参考示例：参数需按材料、赝势和超胞大小收敛测试；多硫化锂构型应成组枚举。

最终结果建议按“材料-位点-Li₂Sₓ-最低能构型-E_ads-关键键长-电荷转移-后续转化计算”组织。这样既能比较吸附强弱，也能避免把不同物种、不同覆盖度和不同结构身份混成一个数字。有比较意义的吸附能结果，应能让别人复现同一组参照和同一组构型筛选规则。

若要进一步连接电池性能，可以把吸附能和 Li₂S 成核能垒、Li₂S 分解能垒、硫还原自由能台阶放在同一逻辑里。吸附能负责回答“能不能抓住”，能垒负责回答“能不能继续转化”，电荷和结构证据负责回答“为什么这个表面有效”。三者合在一起，才适合作为 Li-S 催化剂筛选依据。

还有一个容易被忽略的细节是溶剂和电解液环境。真空 slab 中得到的 Li₂Sₓ 吸附能适合做材料内部分组，但若要和实验穿梭抑制直接对应，应说明没有显式溶剂、盐离子和液相构型熵。真空吸附能是表面锚定能力的第一层指标，不是完整电池环境下的自由能。

因此，正式比较时可以先用真空模型筛掉明显弱吸附表面，再对少数候选加入隐式溶剂、显式 DOL/DME 片段或 Li+ 配位环境做对照。若加入溶剂后排序发生变化，结论应转向“界面环境改变吸附稳定性”，而不是继续沿用真空模型排序。多硫化锂的可比性要随着模型层级一起说明。

本章要点总结

要点：多硫化锂吸附能比较前必须固定 Li₂Sₓ 物种、孤立分子参照、slab 口径和构型搜索规则。

要点：Li 端、S 端、桥连和平躺构型都可能成为低能吸附态，不能只用单一初始构型代表材料能力。

要点：表面面积、活性位数量和覆盖度会改变周期性相互作用，必须随吸附能一起说明。

要点：吸附过强可能影响后续转化，吸附能需要和键长、电荷、自由能或能垒共同判断。

要点：若吸附过程中发生断键或分解，应改用解离吸附能或反应能口径，不能继续写成完整 Li₂Sₓ 吸附能。