GITT计算锂离子扩散系数：从电化学测试到数据解析的关键步骤

GITT（恒电流间歇滴定技术）是锂离子电池领域测量材料体相扩散系数的经典方法。原理上，它通过给一个短时间的恒电流脉冲，然后测开路电压的弛豫过程，利用半无限扩散模型的解来反推扩散系数。但”原理简单”和”结果可靠”之间隔着好几道坎。

测试条件的设定：脉冲时间和电流密度

GITT测试的输入参数有三个需要事先确定的量：脉冲电流密度（I）、脉冲时间（τ）、弛豫时间（休息时间）。脉冲电流密度的选取取决于电池的C-rate能力——如果电流太大，电压响应里会混入浓差极化和欧姆极化，解析出来的扩散系数会偏大；如果电流太小，电压变化信号太弱，信噪比不够。

团队在做LiFePO₄正极的GITT测试时，脉冲电流密度设在了C/20（0.05C）。这个值对LiFePO₄来说是一个保守选择——在这么低的倍率下，浓差极化的贡献可以控制在总电压变化的10%以内。对应的脉冲时间τ设为30分钟，弛豫时间设为4小时（至少是脉冲时间的8倍，确保电压弛豫到稳态）。

脉冲时间τ的选择也有讲究。太短的脉冲（比如5分钟）会让扩散前沿还没进入材料体相，测到的扩散系数反映的是界面过程而非体相扩散；太长的脉冲（比如2小时）会违背半无限扩散模型的前提——扩散前沿到达了颗粒中心。对LiFePO₄来说，颗粒半径一般在100-500 nm，30分钟的脉冲时间大致让锂离子的扩散深度在50-200 nm之间，对大部分颗粒来说还在半无限扩散的范围内。

数据处理：公式里的简化假设

GITT计算锂离子扩散系数的标准公式是：

D = (4/πτ) × (mV_M/MS)² × (ΔE_s/ΔE_t)²

其中ΔE_s是稳态电压变化（弛豫后的电压差），ΔE_t是脉冲期间的瞬态电压变化（扣除IR降之后的电压差）。这个公式的推导基于几个关键假设：颗粒是球形、扩散各向同性、脉冲期间扩散前沿没有到达颗粒中心。

LiFePO₄有一个特殊问题：它有两相共存区（LiFePO₄和FePO₄），在两相区内的GITT数据会出现一个”平台+尖峰”的特征——脉冲瞬间电压急剧上升（相界移动），然后弛豫到稳态电压。在两相共存区内用标准GITT公式算出来的扩散系数可能偏高一到两个数量级，因为公式假设的是单相固溶体扩散，而两相区的控制步骤实际上是相界迁移。

团队在数据处理时，把充电曲线分成了两个区域：单相区（SOC<20%和>80%）用标准GITT公式处理，两相区（20-80%）改用相界迁移模型（修改过的Weppner-Huggins公式）。单相区的扩散系数在10^-14到10^-16 cm²/s之间，与文献报道的LiFePO₄体相扩散系数一致。两相区的等效扩散系数在10^-12到10^-13 cm²/s之间，反映了相界迁移的动力学特征——这两个数值不能直接对比，但各自在物理上是有意义的。

IR降的扣除：一个常被忽视的步骤

脉冲瞬间的电压跳升包含了两个成分：欧姆电阻引起的IR降和浓差极化。GITT公式里用的是扣除IR降之后的ΔE_t。IR降的判定方法是在脉冲开始后的前几秒到几十秒内做线性外推，外推到t=0时的电压就是扣除IR降后的起始电压。

这一步的操作精度对结果影响很大。如果外推的时间窗口太短（比如前5秒），可能还没完全排除电流建立瞬间的非稳态效应；如果太长（比如前2分钟），浓差极化的贡献已经开始混入。团队对LiFePO₄的数据用前30秒做了线性外推，跟前10秒和前60秒的结果做了交叉验证——三种外推窗口给出的ΔE_t差异在3%以内，说明结果对这一步的操作不敏感，结果是可靠的。

与DFT计算的交叉验证

GITT测出的扩散系数是宏观有效值，包含了晶界、缺陷和各向异性的综合贡献。如果需要区分不同晶格方向的本征扩散势垒，需要借助DFT计算——用NEB（微弹性带）方法计算Li⁺沿不同迁移路径的能垒，再用过渡态理论估算扩散系数。

团队在这个项目里配合做了DFT计算：LiFePO₄沿b轴方向（[010]通道）的Li⁺迁移能垒是0.27 eV，沿a轴方向是0.55 eV——扩散的各向异性非常显著。用Nudged Elastic Band方法计算的扩散系数在300 K下约为2×10^-14 cm²/s，与GITT单相区测到的量级一致，数值上偏低约一个数量级——这个差异可以用多晶样品中晶界和缺陷的增强效应来解释。

GITT计算锂离子扩散系数是一项需要实验技巧和物理理解并重的测试工作。在VASP/第一性原理的框架下，DFT可以提供本征扩散势垒的理论预测，与GITT的实测值形成互补。更多电池材料的计算与测试案例，可参阅科研学术网首页。