本研究由美国阿贡国家实验室纳米尺度材料中心的Maria K. Y. Chan和Jeffrey P. Greeley，以及西北大学材料科学系的C. Wolverton共同完成。研究成果以题为《First Principles Simulations of the Electrochemical Lithiation and Delithiation of Faceted Crystalline Silicon》的论文形式，于2012年7月21日发表在《Journal of the American Chemical Society》上。
该研究属于计算材料科学与电化学储能领域，聚焦于锂离子电池下一代高性能硅负极材料的原子尺度机理探索。硅因其极高的理论容量（约为石墨的10倍）而被视为极具前景的负极材料。然而，在锂化过程中，硅会经历复杂的结构演变，包括室温下的电化学固态非晶化、在完全锂化时形成Li15Si4晶相、首次充放电循环中的电压滞后现象，以及在晶体样品中观察到的显著锂化各向异性。尽管这些现象已被实验观测到，但其在原子层面上的详细过程和根本物理原因尚不清楚。为此，本研究旨在开发和应用一种基于第一性原理的、具有历史依赖性的锂离子嵌入/脱出算法，来模拟晶体硅的锂化和随后的脱锂过程，以期从最基础的理论层面揭示上述现象的原子机理，并为未来设计性能更优的硅基负极提供坚实的理论依据。
研究的核心工作流程围绕自主研发的锂化/脱锂模拟算法展开，并辅以详细的原子结构分析和能量计算。
首先，研究团队构建了不同晶面的晶体硅表面模型。他们选取了硅晶体最常见的（100）、（110）和（111）三个晶面作为研究对象。对于每个晶面，均考虑了不同重构方式（如（100）面的2×1和2×2对称/弯曲二聚体重构）以及不同尺寸的周期性原胞（每层包含2个或4个硅原子）。这些表面模型底部通过氢原子钝化，以模拟块体基底，并在垂直方向留有至少17 Å的真空层。所有计算均基于密度泛函理论（Density Functional Theory, DFT），使用VASP软件包，采用PBE泛函描述交换关联作用。
其次，关键的研究工具——历史依赖的锂化/脱锂算法被开发并应用。该算法的核心是逐步、顺序地插入或移除锂原子，以模拟非平衡的室温电化学过程。其具体步骤如下： 1. 起始点确定：首先针对每个表面模型，确定表面锂吸附位点，并找到表面锂化电压最高（最易发生）的构型，以此作为后续体相锂化的起点。 2. 模拟区域定义：确定一个厚度为5层（约7-10 Å）的“可锂化”硅薄层区域。该区域下方的硅原子保持其晶格位置固定，以模拟未受扰动的晶体基底对上方区域的约束。 3. 迭代锂化过程：对于每一个待插入的锂原子，执行以下子步骤： a. 候选位点生成：在可锂化区域内，基于一个规则的空间网格，生成所有距离最近原子至少d_min（1.9 Å）的间隙位点。 b. 能量计算与筛选：将锂原子依次放入每个候选位点，使用DFT“快速”参数对包含该锂原子的整个结构进行离子弛豫，计算其总能。 c. 最优构型选择：选择所有尝试中总能量最低的弛豫后构型，并使用DFT“精确”参数再次弛豫，得到该锂含量下的最稳定（或亚稳态）构型。该构型将作为插入下一个锂原子的起点。 4. 脱锂模拟：脱锂过程被模拟为上述锂化过程的逆过程。从锂化程度最高的稳定构型开始，逐一尝试移除每一个锂原子，计算移除后的结构弛豫能量，并选择能量最低的构型进行下一步脱锂。 为了区分表面效应和晶体取向本身的内在影响，研究团队还进行了另一组模拟：使用不含真空层的周期性原胞（沿<100>, <110>, <111>方向），仅固定横向尺寸，允许纵向膨胀，以模拟特定晶向的几何约束下的锂化，而无显式界面。
再次，基于模拟生成的原子构型和总能量，进行了一系列深入的分析： 1. 结构分析：通过可视化锂化过程中的原子构型演变，直观观察非晶化进程。计算了硅-硅、锂-硅、锂-锂的径向分布函数（Radial Distribution Function, RDF），并与实验数据以及已知晶态LixSi相的RDF进行对比。 2. 能量与电压分析：根据每个锂含量构型的形成能（相对于金属锂和未锂化硅表面），构建了能量-成分关系的凸包线，并基于一系列两相反应的假设，推导出相应的理论电压曲线。这是判断不同晶面锂化热力学倾向性的关键。 3. 扩散能垒计算：使用攀爬图像微动弹性带法（Climbing Image Nudged Elastic Band, CI-NEB）计算了锂在不同表面覆盖层、亚表层以及锂化-未锂化界面附近的扩散能垒，以评估动力学各向异性的可能贡献。 4. 应变与键合能分解：将锂化构型的形成能分解为锂原子应变能、硅骨架应变能和锂-硅键合能，以量化不同晶面在锂化过程中能量变化的来源。
本研究获得了丰富且相互印证的原子尺度结果。
在原子结构演变方面，模拟成功再现了晶体硅在室温锂化下的非晶化过程。结果显示，锂化是一个明显的两相反应过程，随着锂的插入，清晰的“非晶-晶体”界面向内推进。更重要的是，不同晶面的非晶化路径存在显著差异：（110）面和（100）面倾向于形成“之字形”硅链或二聚体，其中（110）面形成的链更为有序；而（111）面则先形成相互连接的硅原子层，随后在高锂含量下才进一步破碎。径向分布函数分析表明，模拟得到的结构在关键特征（如高锂含量下依然存在的Si-Si最近邻键合）上与实验测量结果高度吻合，验证了算法的可靠性。此外，在模拟得到的最高锂含量（Li~4.1Si）构型中，发现了孤立硅原子被约12个锂原子配位的局域结构，这与Li15Si4晶相中的硅配位环境相似，为实验观察到的“非晶向Li15Si4结晶”转变提供了原子层面的线索。
在热力学与电压特性方面，模拟得出的电压曲线揭示了关键的各向异性现象。对于具有真空层的表面模型，（110）面的锂化电压平台明显高于（100）和（111）面。这一差异在无真空层、仅考虑晶向几何约束的模拟中依然存在，说明电压差异源于硅晶体自身在不同方向上的结构对称性和键合特性，而非表面本身。通过能量分解发现，（110）面虽然硅骨架应变能较大，但其锂-硅键合能也更强（更负），从而在净形成能上占据优势，导致更高的锂化电压。这一热力学倾向性具有重要推论：在由不同晶面构成的硅微米/纳米结构中，锂将优先从热力学更有利的（110）面嵌入。为了证实这一点，研究者求解了具有各向异性表面锂溶解度的扩散方程，结果显示即使假定扩散系数各向同性，热力学上的差异也足以导致实验观测到的严重各向异性膨胀形貌（例如，微柱优先沿某些方向膨胀），从而为实验现象提供了一个更合理的、基于热力学的解释，替代了之前可能过度依赖扩散动力学各向异性的解释。
在动力学方面，锂扩散能垒的计算表明，尽管在单层锂覆盖的表面，锂从表面向亚表层扩散的能垒存在差异（（110）面较低），但在进入更深层或通过锂化-未锂化界面时，不同晶向的扩散能垒变得相似。这表明扩散各向异性在微米尺度结构中的作用有限，进一步支持了热力学主导机制。
在脱锂与滞后现象方面，模拟成功捕捉到了首次充放电循环中的电压滞后。脱锂过程的电压平台整体高于首次锂化。这是因为脱锂后的硅无法恢复到原始的晶体状态，而是形成了能量更高的非晶硅层。滞后回线的大小与非晶硅与晶体硅之间的能量差相关。模拟显示，（100）面的滞后最大，暗示其脱锂后非晶化程度更高或晶体恢复程度更低。对脱锂后结构进行的二次锂化模拟显示，第二次锂化的电压曲线更接近第一次脱锂曲线，滞后减小，这与实验观察到的首次循环后库仑效率改善的趋势一致。
本研究的结论是，通过自主开发的、物理动机明确的DFT模拟方法，首次在原子尺度上详细描绘了晶体硅电化学锂化/脱锂的复杂图景。研究阐明了不同晶面非晶化路径的结构差异，并从热力学角度首次明确解释了硅负极锂化各向异性的根本原因——即（110）晶面具有更高的锂化电压，这源于其硅骨架在特定应变模式下与锂原子更强的键合相互作用。这一热力学优势足以驱动微米/纳米结构中的选择性锂化与形变，无需引入显著的扩散各向异性。同时，研究揭示了电压滞后的起源与非晶-晶体能量差直接相关，并为Li15Si4相的形成提供了原子层面的见解。这些发现将实验观测到的宏观现象与原子尺度的电子结构和化学键合联系起来，为理解硅负极的工作机理建立了坚实的理论基础。
本研究的亮点突出体现在以下几个方面：第一，方法创新：开发了首个基于第一性原理的、历史依赖的锂化/脱锂模拟算法，能够有效地模拟室温非平衡电化学过程，填补了传统基态搜索方法在此类问题上的空白。第二，机理突破：首次从热力学角度而非动力学角度，合理解释了硅锂化各向异性这一关键实验现象，提出了更普适和深刻的物理图像，解决了原有解释的困境。第三，预测与验证：模拟得到的原子结构（如RDF）与实验数据高度吻合，不仅验证了方法的可靠性，也提供了实验难以直接获取的锂原子局域环境信息。第四，系统性与深度：研究涵盖了从锂化、脱锂到二次锂化的完整循环，结合了结构、能量、动力学和电压的全面分析，构成了对硅负极电化学行为的系统性原子尺度阐释。
此外，研究还隐含了重要的应用指导价值：理解各向异性有助于设计特定的电极微结构（如取向生长的纳米线或图案化薄膜），以优化锂离子传输路径和缓解应力；而对滞后和非晶化过程的认识，则为通过表面工程、界面调控或复合结构设计来改善硅负极的循环稳定性提供了理论线索。