Revealing Short- and Long-range Li-ion diffusion in Li$_2$MnO$_3$ from finite-temperature dynamical mean field theory

通过将有限温度动力学平均场理论与 DFT+$U$及爬行弹性带计算相结合，本研究揭示动力学关联显著降低了顺磁 Li$_2$MnO$_3$中的锂离子迁移势垒，从而在不依赖外源无序或簇状空位的情况下，为短程和长程输运的活化能提供了统一的解释。

要点

想象一下，将锂离子电池比作一座繁忙的城市，其中微小的锂离子是通勤者，试图从城市的一端前往另一端，以为你的手机或汽车提供动力。它们行驶的“道路”位于一种名为Li₂MnO₃的材料内部。

长期以来，科学家们一直困惑于这些通勤者的移动速度。一些实验（观察极短距离）称道路超级平滑且快速。另一些实验（观察长距离）则称道路拥堵不堪且非常缓慢。这就像有人说：“你可以在 10 秒内跑完冲刺！”但同时又说：“你无法跑完马拉松，因为赛道是坏的。”

本文通过一种超先进的计算机模拟，以新的方式观察“交通”，从而解开了这一谜团。

旧地图与新地图

此前，科学家们使用一种标准的计算机模型（称为 DFT+U）来绘制道路图。该模型就像一个基础 GPS：它看到了锂离子试图翻越墙壁，但计算出的墙壁非常高（约 0.6 至 0.9 eV）。这表明离子移动非常缓慢，这与短距离实验中观察到的快速“冲刺”数据不符。

作者意识到，旧模型缺失了一个关键要素：热量与混乱。在现实世界中，电池内的原子并非静止不动；由于热量（温度），它们正在颤动和振动。材料中的锰原子还具有微小的自旋磁矩，它们正在随机翻转。旧模型将这些自旋视为完美排列的冻结状态，但这并不符合实际工作的电池。

“动态”模拟

为了解决这个问题，作者使用了一种更强大的工具，称为DFT+DMFT。这可以想象成从静态的二维地图升级为考虑了热量和磁自旋随机翻转的 3D 实时模拟。

他们模拟了锂城市中的一个“空座位”（空位）。锂离子需要跳入这个空座位才能向前移动。

两种移动速度

当他们运行新的“热且混乱”的模拟时，发现了一些惊人的现象。能量势垒（离子必须翻越的墙壁）显著降低，但仅针对特定类型的跳跃。

1. 短途跳跃（冲刺）：
   对于两个相邻位置之间的极短跳跃，新模拟显示墙壁高度仅为0.18 eV。

   • 结果： 这与短距离实验中观察到的“快速冲刺”数据完美吻合。
   • 类比： 想象一名通勤者跨过一个小路缘。这既容易又快速。旧模型认为路缘是 10 英尺高的栅栏；而新模型意识到那仅仅是一小步。

2. 长途跋涉（马拉松）：
   然而，要穿越整个城市进行长途旅行，通勤者不能永远只走轻松的步伐。他们最终必须迈出一小步更难的步伐。模拟发现了一个稍高的第二道墙壁，高度为0.50 eV。

   • 结果： 这与长距离实验中观察到的“缓慢马拉松”数据相符。
   • 类比： 要穿越全城，你必须迈出许多轻松的步伐，但偶尔会遇到一座小山丘。即使大多数步伐都很轻松，你的整体速度仍受那座山丘的限制。

为何这很重要

重大发现在于，你无需发明复杂的解释来解决速度问题。你无需假设电池充满了空座位的“团块”，也无需假设材料是破损的。

该论文表明，Li₂MnO₃实际上是一种非常好的材料（近乎完美，或称“化学计量”）。我们在不同实验中观察到不同速度的原因仅仅是：

• 短程实验只看到了容易的低矮山丘（0.18 eV）。
• 长程实验看到了整个旅程，而旅程偶尔会被较高的山丘（0.50 eV）所减缓。

结语

通过考虑原子的热量和磁“抖动”，作者构建了一个统一的单一故事。他们证明了锂离子在局部尺度上可以轻松穿梭，但其整体旅程受少数稍难的步骤控制。这解释了为什么电池的表现取决于测量方式，而无需归咎于材料中的缺陷或杂质。

简而言之：电池并没有坏；我们只是需要一张更好的地图，该地图考虑了原子的热量和磁舞，以理解锂离子真正的移动方式。

技术摘要

技术摘要：基于有限温度动力学平均场理论揭示 Li2MnO3 中锂离子的短程与长程扩散

问题陈述
Li2MnO3 是富锂层状正极材料的关键组分，但其在限制锂离子传输中的作用仍存在争议。实验上对 Li2MnO3 中锂离子迁移的估算因探测尺度不同而表现出显著差异。在近乎化学计量的粉末样品上进行的微观μ+SR 测量显示，在顺磁态下短程扩散的活化能较低（Ea ≈ 0.156 eV）。相比之下，宏观交流阻抗测量则得出长程传输的表观势垒显著更高（Ea ≈ 0.46 eV）。

此前基于 Hubbard U 修正的密度泛函理论（DFT+U）的第一性原理研究难以调和这些数值。针对单个锂空位的标准 DFT+U 计算通常预测 0.5–0.8 eV 的势垒，这过高而无法解释短程实验数据。为了匹配较低的实验势垒，先前的模型往往引入聚集的空位构型（双空位或三空位），或假设特定的磁有序。然而，在近乎化学计量的样品中，此类空位团簇并不常见，且 Li2MnO3 在电池工作条件下呈顺磁性，而许多 DFT+U 研究假设了铁磁态或非磁态，从而抑制了局域锰磁矩。

方法论
本研究采用多尺度计算方法，在室温（300 K）下研究含单个锂空位的顺磁 Li2MnO3 中的 Li+ 迁移。该方法整合了三个组成部分：

1. DFT+U：用于通过nudged-elastic-band（NEB）方法优化迁移路径和鞍点几何结构。采用 PBEsol 泛函和 1×2×2 超胞进行自旋极化 DFT+U 计算。
2. DFT+DMFT：应用有限温度动力学平均场理论（DMFT），结合连续时间量子蒙特卡洛（CTQMC）杂质求解器，在固定的 NEB 几何结构上评估总能量。该框架捕捉了顺磁相中的动力学电子关联。
3. 电子结构处理：针对 Li2MnO3 的低晶体对称性（C2/m），该对称性导致标准 Wannier 投影中出现较大的非对角矩阵元，作者对 Mn-d 块进行了对角化处理。这构建了一个有效的低能“仅 d 轨道”模型，其中 Wannier 基组隐式编码了 Mn-d–O-p 杂化，从而能够准确复现绝缘能隙。

该研究通过计算初始、鞍点和终态图像的 DMFT 总能量，评估了六种对称性不等价的迁移路径（四条层内路径和两条层间路径），并在每个点上允许关联电子态弛豫。

主要结果
有限温度动力学关联的应用显著重正化了特定迁移路径的活化能，在无需引入聚集空位的情况下，解决了微观与宏观实验数据之间的差异：

• 短程扩散：对于势垒最低的层内跳跃（路径 ii：$4h \to T_{2b}^{Li} \to 4h$），DFT+DMFT 计算的活化能降低至0.18 eV。该数值定量复现了μ+SR 测量中观察到的短程活化能（0.156 eV）。这种降低主要由 DMFT 总能量中的本征值修正项驱动，该修正项稳定了鞍点构型。
• 长程扩散：次低势垒路径（路径 iv：$4h \to T_{2b}^{Li} \to 2c$）在 DMFT 下的活化能为0.50 eV。该势垒控制着锂离子的长程渗流。它与从交流阻抗测量中提取的宏观活化能（≈ 0.46 eV）一致。
• 路径依赖性：势垒降低具有高度的路径依赖性。虽然路径 (ii) 从 DFT+U 值（约 0.68 eV）急剧下降至 0.18 eV，但其他路径（例如路径 iii）变化甚微。作者将此与迁移锂离子相对于 MnO6 网络的局域排列相关联；路径 (ii) 穿过更开阔的区域，具有更长的 Li-Mn 距离，因此能更有效地被动力学关联所稳定。
• 电子结构：与 DFT+U 在产生空位后显示出微弱的位点依赖性 Mn-d 占据变化不同，DMFT 揭示了谱权重的强位点依赖性重分布，表现为远离空位的 Mn 位点上占据态 Mn-d 权重的显著耗竭。

意义与主张
该论文声称提供了对 Li2MnO3 中局部和宏观锂离子传输的统一微观解释。通过利用单空位、顺磁性的 DMFT 描述，作者证明活化能的层级结构（短程较低，长程较高）自然源于特定迁移路径与有限温度动力学关联之间的相互作用。

研究结论如下：

1. 实验观察到的短程活化能可以通过顺磁宿主中的单空位迁移来合理化，无需引入外源无序或聚集空位构型来解释数据。
2. 长程传输由跳跃网络中较高势垒的步骤（0.50 eV）控制，这与宏观阻抗数据一致。
3. 有限温度动力学关联是预测关联氧化物电极中离子迁移能量的关键要素，因为静态 DFT+U 方法无法捕捉势垒高度所需的必要重正化。

作者指出，这些发现强调了在关联氧化物中进行离子扩散研究时超越 DFT 处理的必要性，并建议未来的工作应将此框架扩展至更高空位浓度以及混合阳离子的富锂层状正极，特别是在氧氧化还原变得活跃的体系中。