富氧石榴石Li7La3M2O13(M = Ta, Nb) 的微观结构解析与锂离子传导抑制机制

《JOURNAL OF POWER SOURCES》：Garnets Li7La3M2O13 (M = Ta, Nb): Stable phases and superior ionic conductors?

【编辑推荐】为解决锂镧钽/铌氧“O13”型石榴石固态电解质离子电导率远低于预期、其结构与稳定性和导电机理不明的难题，研究人员通过结合扩散模型MatterGen、SOSS结构搜索、密度泛函理论（DFT）和机器学习势（CHGNet）分子动力学模拟，揭示了此类材料为亚稳态，其额外氧原子引入的局部扭曲改变了锂亚晶格连通性，形成阻碍长程扩散的瓶颈，从而解释了其低电导率。此项研究不仅阐明了“O13”化合物的导电短板，其方法学也为高性能锂填充固态电解质的理性设计提供了数据驱动的指导。

　　

研究背景与亟待破解的“石榴石”之谜

自1991年诞生以来，锂离子电池经历了飞速发展，但传统液态电解质带来的安全隐患和能量密度瓶颈，促使科学家们将目光投向固态电解质。其中，石榴石型陶瓷电解质，特别是Li₇La₃Zr₂O₁₂(LLZO)，以其高离子电导率和良好的电化学稳定性，被视为极具潜力的候选者。然而，在石榴石家族中，还存在一类鲜为人知的成员——化学式为Li₇La₃M₂O₁₃(M = Nb, Ta) 的富氧“O₁₃”型石榴石。理论上，它们拥有与明星材料LLZO相似的立方框架和更高的锂含量，理应表现出优异的离子传导性能。但令人困惑的是，已报道的实验数据显示，无论是Li₇La₃Nb₂O₁₃(LLNO₁₃) 还是Li₇La₃Ta₂O₁₃(LLTO₁₃)，其室温离子电导率都比LLZO低了一个数量级以上。这就像一个“石榴石悖论”：结构与锂含量相似，性能却天差地别。此外，这些材料中“额外”氧原子的确切占位、材料本身的热力学和电化学稳定性、以及导致其低电导率的微观机制，都笼罩在迷雾之中。为了揭开“O₁₃”石榴石的神秘面纱，并从中汲取对设计下一代固态电解质至关重要的普适性见解，一项由Henry A. Cortés等人领导的综合性计算研究在《JOURNAL OF POWER SOURCES》上展开。

研究方法概览：人工智能与计算模拟的深度联姻

为系统探究LLTO₁₃和LLNO₁₃，研究团队构建了一套先进的计算工作流。首先，他们利用生成式人工智能扩散模型MatterGen，在Li-La-Ta-O化学空间中广泛生成候选结构，并结合“最优稳定结构搜索”（SOSS）算法，高效探索了LLTO₁₂中锂空位的最佳排列。随后，通过向LLTO₁₂晶格中引入Li₂O分子，构建了“O₁₃”结构模型。所有候选结构的能量和优化通过密度泛函理论（DFT）计算完成，并用于评估材料的热力学稳定性（构建相图与凸包）和电化学稳定窗口。为了深入分析锂离子动力学和由过量氧引入的微观结构变化，研究采用了基于机器学习势（MLP）CHGNet的分子动力学（MD）模拟，并在500-1300 K的温度范围内进行了长时间模拟。最后，利用“晶体学位点发现器”（CrySF）这一无监督数据分析方法，对MD轨迹进行自动化的位点识别、表征和锂离子扩散网络分析，精准定位了离子跳跃路径和潜在的捕获位点。

研究结果：层层剥茧，揭示低电导率的根源

3.1. 晶格膨胀

DFT计算预测，LLTO₁₃的晶胞参数相对于LLTO₁₂平均膨胀约4.5%，这与实验观察到的Li₂O掺入导致的轻微膨胀趋势一致，但计算值偏大，可能源于模型中Li₂O分布尚未完全优化或实验样品中存在其他物相。

3.2. 热力学稳定性

通过MatterGen生成相图并结合DFT能量计算，研究发现LLTO₁₃和LLNO₁₃均位于凸包之上，但能量分别仅高出0.058 eV/atom和0.0958 eV/atom，处于通常认为可合成的亚稳态窗口内。这表明它们虽然在热力学上不是最稳定的相，但通过实验手段合成是可能的。

3.3. 电化学稳定窗口

计算得到的电化学稳定窗口显示，两种“O₁₃”化合物的稳定窗口均比其对应的“O₁₂”母相更窄。具体而言，LLTO₁₃的窗口为0.35-2.75 V，而LLTO₁₂为0.45-3.30 V；LLNO₁₃为1.15-2.15 V，而LLNO₁₂为1.15-2.90 V。这表明“O₁₃”材料在充放电过程中更易发生分解，尤其是对脱锂（高电压侧）更敏感。

3.4. LLTO₁₃的结构

通过CrySF对MD模拟轨迹的分析，揭示了过量氧原子（O_m）在晶格中的最终归宿：它们并未停留在初始假设的锂八面体位点中心，而是通过两种方式重组。一部分直接迁移到间隙位（O_int）；另一部分则“挤占”邻近的96h氧位点，导致原有氧原子被踢出形成新的间隙氧。这种重组伴随着局部TaO₆八面体的扭曲。更重要的是，锂离子扩散网络发生了显著改变。每个O_int位点附近都有一个锂离子与其紧密结合（键长约1.9 ?），形成了稳定的“Li-O_int”单元，这些锂离子的位点占据率高达~0.99，表明它们几乎被“钉扎”在原地，成为不参与长程迁移的“死锂”。此外，部分原本在LLTO₁₂中连通的锂位点消失或发生位移，导致扩散网络出现瓶颈和阻断。X射线衍射（XRD）模拟表明，尽管存在这些微观结构变化，LLTO₁₃的衍射谱与LLTO₁₂的谱图极为相似，仅因晶胞膨胀有微小峰位偏移，这从侧面说明仅凭XRD难以明确区分和证实“O₁₃”相的形成，需要结合中子衍射、热重分析等其他表征手段。

3.5. LLTO₁₃的离子电导率

基于CHGNet势的MD模拟计算了离子电导率。结果清晰地显示，在整个500-1300 K的温度范围内，LLTO₁₃的离子电导率均低于LLTO₁₂。将结果外推至室温，LLTO₁₃的电导率估计值仅为LLTO₁₂的80%。对LLNO体系的计算也得出了相同结论，且Nb基材料的电导率普遍高于同构的Ta基材料。对扩散网络的深入分析（图8）直观地揭示了原因：在LLTO₁₂中，锂离子可以在整个三维网络中进行相对均匀的扩散，包括偶尔发生的八面体-八面体（Oh?Oh）直接跳跃；而在LLTO₁₃中，高频的离子跳跃被限制在一些孤立的局部区域内，区域之间仅通过低频的跳跃连接，扩散的“长程连通性”被严重破坏。被O_int捕获的锂离子以及因结构重组而堵塞的迁移通道，共同抑制了锂离子的宏观输运。

结论与意义：从“失败案例”中汲取的设计智慧

本研究通过多层次、多尺度的计算模拟，首次对富氧“O₁₃”型石榴石LLTO₁₃和LLNO₁₃的结构、稳定性和离子传输性质提供了系统而深入的理解。核心结论在于，这些材料并非如早期猜测那样是潜在的高性能导体，反而是离子传导的“反面教材”。其根本原因在于，为容纳过量氧（Li₂O）而进行的锂“填充”（Li-stuffing），在微观上引发了氧亚晶格的重排和局部结构扭曲。这不仅没有创造出更多畅通的锂离子迁移通道，反而通过形成强局域化的“Li-O_int”束缚对和改变锂位点网络拓扑结构，引入了阻碍离子长程迁移的瓶颈和陷阱，最终导致了较低的体相离子电导率。同时，研究确认了它们是亚稳态材料，且电化学稳定窗口较窄，在电池应用中面临分解风险。

尽管LLTO₁₃和LLNO₁₃本身作为固态电解质前景黯淡，但这项工作的意义却十分深远。首先，它破解了长期存在的“石榴石性能悖论”，为理解化学计量与结构微调对离子传输的复杂影响提供了关键机理认识。其次，研究揭示了仅凭XRD难以唯一确定此类复杂掺杂相，强调了多技术联用表征的必要性。最后，也是最重要的，研究所开发和集成的先进计算方法学（融合了生成式AI、全局结构搜索、机器学习势和无监督轨迹分析），为未来高性能固态电解质的“理性设计”树立了范本。它指明，单纯增加锂含量并非提高电导率的万能钥匙，必须精细调控晶体结构以避免产生离子迁移的微观瓶颈。这些从“失败”化合物中获得的深刻见解，正指引着科学家们避开设计陷阱，朝着开发真正高效、稳定的下一代全固态电池电解质迈出更坚实的步伐。