【纯计算】npj Comput. Mater.：通过硫取代缓解Li1.2Ni0.6Mn0.2O2阴极材料中晶格氧演化

研究背景

为了满足对大规模储能设备日益增长的需求，提高商业大规模储能设备的能量密度将是必要的。因此，设计和开发具有高比容量和氧化还原电位的创新阴极材料至关重要。在此，台湾科技大学Jyh-Chiang Jiang等利用密度泛函理论（DFT）计算和从头算分子动力学（AIMD）模拟，证明了通过S-阴离子取代可以成功地抑制高容量阴极Li_1.2Ni_0.6Mn_0.2O₂（LNMO）表面的晶格氧释放。这为提高LNMO（003）表面的结构稳定性提供了见解，为下一代锂开发富锂LNMO阴极材料铺平了道路。

计算方法

本文通过VASP量子计算软件包实现催化性质的计算，电子-离子相互作用通过投影增强平面波赝势（PAW）来描述。通过GGA-PBE泛函进行电子交互关联的分析，采用optB88-vdW函数中包含的长程色散校正来描述范德瓦尔斯（vdW）的相互作用。将总能量的收敛准则设为10⁻⁴eV。所有的原子位置和晶格结构都被完全放松，最大力的阈值为0.01eVÅ⁻¹，以保证计算的准确性和效率，通过之前研究，选定LNMO（003）进行计算，同时考虑O端，TM端及Li端表面模型。

结果与讨论

为了进一步研究LNMO（003）表面在脱氧过程中的氧演化机制，我们利用AIMD模拟方法比较了O²⁻阴离子和TM金属的各脱氧结构和氧化态。在析出过程中，考虑了原始Li_1.2−xNi_0.2Mn_0.6O₂（003）表面中从x = 0到0.8的不同析出状态。在Li_1.2−xNi_0.2Mn_0.6O₂（003）表面上，一次去除4个锂原子。在解离前模拟的原始模型如图1a所示，电解质组分的优化几何形状如图1b所示。在从LNMO（003）表面去除Li离子后，系统平衡5 ps，在此期间它表现出结构变化，达到新的平衡。

图1. AIMD计算和平均插层电压分布的仿真模型

作者展示了阴极/电解质界面从低到高的解离水平x = 0.13, 0.26, 0.4, 0.53, 0.67和0.80的构型，如图2所示。从图1c中可以看出，LNMO正极材料在x = 0.07的析出水平下表现出~4 V的高电压。电压分布表明，随着析出程度的增加，电压逐渐增大。

此外，还分析了LNMO(003)在脱除过程中的表面构型，其构型如图2所示。LiPF₆和溶剂分子(EC和DEC)未发生分解，表明它们在LNMO(003)表面具有良好的稳定性。LNMO(003)表面在x = 0.26时表现出较高的结构稳定性；在x = 0.4时，LNMO(003)表面与电解质界面形成氧分子。

当解离水平达到x = 0.67时，Li_1.2−xNi_0.2Mn_0.6O₂(003)表面的第二层发现了第二次氧分子的形成，这表明部分O原子在高解离程度下是欠配位的(少于3个阳离子)。这种配位不足的O原子可能会发生析氧并进一步降解在LNMO(003)表面，从而导致不可逆反应和容量衰减。

图2. 针对原始LNMO系统的优化结构

晶格氧的参与与O^{2 –}阴离子的氧化活性有关。Bader电荷分析用于表征和可视化LNMO(003)表面原子电荷的变化过程。计算得到的O阴离子在不同脱出水平时的Bader荷如图3a所示。在x = 0的初始状态下，O阴离子的平均电荷约为- 1.3|e|。随着脱除水平的提高，O负离子的平均原子电荷负电荷变小，表明晶格氧负离子显著参与了脱除过程中的电荷补偿过程。当LNMO(003)由初始状态变为x = 0.33时，未检测到O阴离子的过度氧化。然而，在x = 0.4的解离水平上，原子电荷从- 1.0到- 0.2|e|突然增加，表明氧化状态从O²⁻到O₂ⁿ⁻发生了显著变化。

计算得出的巴德电荷表明，氧是通过原子电荷突然增加(大于- 0.4|e|)而形成的，这可以作为氧气生成的基线。此外，观察到第二次氧生成在沉降水平x = 0.67，原子电荷的值也超过了氧生成的基线，从- 0.9到- 0.3|e|(图3a中的绿线)。

为了进一步了解LNMO(003)表面氧的详细析氧机制，作者计算了LNMO(003)表面O阴离子和TM离子(Ni和Mn)在不同析氧水平下的原子电荷变化(ΔQ)，结果如图3b所示。ΔQ的正值表明原子失去电子，氧化态增加。在整个脱除过程中，Mn原子的ΔQ值略有增加，之后几乎保持不变，说明Mn原子没有参与电荷补偿过程。与Mn原子相比，Ni原子的ΔQ值在析出过程中有较明显的增加，说明Ni原子主要参与氧化过程，这与我们之前的研究和其他研究一致。

图3. LNMO电子性质分析

为了进一步探究在原始和S取代的LNMO表面提取Li原子后Ni电子性质的变化，在x = 0和0.4的解离水平下，投影了第一层和第二层TM层Ni和Mn原子三维轨道的态密度(PDOS)，如图4所示。S取代的LNMO表面PDOS(图4c)显示，与原始体系相比，在x = 0处Ni价带明显上移，接近费米能级(图4a)，表明在S取代体系中Ni氧化比原始体系更有利。如图4d所示，与原始体系(图4b)相比，S取代体系中Ni原子的三维态在解离后(x = 0.4)从价带向导带转移更多。

此外，在x = 0和x = 0.4的解离水平下，S取代的LNMO(003)表面Ni原子的三维态在费米能级以上的积分面积分别为16.2和17.4，表明与原始体系相比，S取代可以提高Ni原子的失去电子的能力。总的来说，增加Ni原子的氧化程度可以显著提高容量保持率。

图4. LNMO的PDOS图像

在富锂的LNMO正极材料中，晶格氧的形成通常会产生不需要的O空位，而O空位过多会导致LNMO正极的结构退化，从而导致lib中容量/电压的快速衰减。为了研究s取代对LNMO表面O空位形成的影响，对原始Li_1.2−xNi_0.2Mn_0.6O₂(003)表面进行了逐步的脱除过程。

然后，作者考虑了原始和S取代的LNMO表面第一层和第二层O空位位置的所有可能性，选择最有利的O空位位置进行下面的计算。图5a显示了x = 0和0.4时原始和s取代的LNMO(003)表面的优化结构。黑色虚线圈表示单位单元中O空位的位置。

然后评估了单位细胞中O空位的形成，Li_1.2Ni_0.2Mn_0.6O₂的脱除状态为Li_0.4Ni_0.2Mn_0.6O₂, O空位与脱除状态的形成能如图5b所示。如图5b所示，从原始LNMO(003)表面去除Li原子后，O空位的形成能急剧降低，表明O空位的形成在热力学上变得更加有利。当脱除态达到x = 0.4时，O空位的形成能为负，这意味着O空位可以在原始表面自发形成。然而，S阴离子的取代抑制了S取代的LNMO(003)表面O空位的形成。

图5. O空位的优化构型和形成能

循环过程中氧空位和TM迁移引起的不可逆的层状结构向新相的相变影响了LNMO阴极的结构稳定性，导致电压衰减和容量衰退。为了理解这一点，作者探索了原始和s取代的LNMO(003)表面在脱除过程中每个TM离子的详细构型。图7为原始Li_1.2−xNi_0.2Mn_0.6O₂(003)表面和s取代Li_1.2−xNi_0.2Mn_0.6O₂(003)表面的结构变形。观察到原始LNMO表面的结构变形发生在x = 0.4的精化水平，这是由于Mn和Ni原子的数量从八面体增加到四面体。对于原始体系，MnO₆和NiO₆的八面体转变为MnO₄和NiO₄的四面体，Mn-O环境在x = 0.4时形成了表面氧。

此外，当沉降状态达到x = 0.67时，可以观察到S取代表面的结构变形，这种结构变化的程度小于较高沉降程度的原始表面。这一结果表明，S取代的LNMO(003)在低程度的析出状态下结构变化不明显。此外，S取代的LNMO(003)在x = 0.4时保持层状结构(图6)，表明S取代可以有效防止装药过程中的结构变形。由此可见，在Li_1.2−xNi_0.2Mn_0.6O₂正极材料中，S-取代在电子缓冲层中发挥了重要作用，参与电荷补偿，从而进一步避免了晶格氧的形成和释放，有效地减小了锂离子分解过程中的结构变形。

图6. 结构变形的说明

图7. 研究中进行AIMD的设计流程

结论与展望

综上所述，作者利用DFT和AIMD模拟方法，系统地研究和比较了电解质混合物存在下原始和S取代LNMO（003）表面的O空位形成和晶格氧演化机理。本文还分析了原始LNMO（003）表面的结构和电子性质，在高水平的解化过程中观察到一些过氧化的O阴离子，这些O阴离子发生氧演化，导致O空位的形成。

该结果有助于理解LNMO正极材料的析氧机制，特别强调抑制晶格氧形成的策略，以开发具有高可逆性和高结构稳定性的高性能锂离子正极材料。

原创文章，作者：计算搬砖工程师，如若转载，请注明来源华算科技，注明出处：https://www.v-suan.com/index.php/2024/03/29/0f14c59636/

【纯计算】npj Comput. Mater.：通过硫取代缓解Li1.2Ni0.6Mn0.2O2阴极材料中晶格氧演化

相关推荐