GITT法计算锂离子扩散系数的实验设计与数据处理要点

GITT计算锂离子扩散系数，原理上就三行公式：施加恒定电流脉冲→记录电压随时间的变化→把稳态和瞬态的电压差代入Fick第二定律的简化解，得到扩散系数D。但这个”三行公式”背后，每次拿到新材料的原始数据，真正花时间的都不是公式本身，而是判断哪些数据是可信的、哪些段要扔掉。

去年做过一个NCM811正极材料的扩散系数测量。客户的要求很直接：在不同的锂含量（x=0.2到0.9之间）给出一组D值，误差控制在±20%以内。这个精度放在电池动力学领域不算苛刻，但实际数据跑出来的分散程度远超预期。

脉冲参数设置：电流大小和时间的”组合拳”

GITT的核心假设之一是：在脉冲期间，锂离子在电极颗粒中的扩散满足半无限扩散条件。这个条件成立需要两个前提：脉冲时间足够短（扩散层厚度远小于颗粒半径），同时电流不能太大（避免副反应和欧姆过电位淹扩散过电位）。

对于NCM811这种二次颗粒直径大约5-10 μm的体系，如果按教科书推荐的C/20电流、脉冲30分钟来跑，扩散层厚度大约是1-2 μm——和颗粒半径在同一量级，半无限扩散条件事实上不成立，算出来的D值会严重低估。

团队对比了两种脉冲方案：方案A（C/20，30分钟）和方案B（C/10，10分钟）。方案B扩散层厚度控制在~0.5 μm，满足半无限假设，测出来的D值在10⁻¹¹到10⁻¹⁰ cm²/s之间。方案A的D值系统偏低0.3-0.5个数量级——不是扩散系数真的差那么多，而是边界条件失效了。

脉冲时间还有一个被忽略的影响：短脉冲对弛豫判据更敏感。每次脉冲后要让电池弛豫到新平衡态才能读取稳态电压E_s，弛豫时间设得太短，读到的E_s其实不是真正的平衡态，D值自然不准。

弛豫判据：dV/dt < 什么时候才算”稳定”

弛豫的终点判定通常写成”电压变化率小于某个阈值”，常见的是dV/dt < 0.5 mV/h或1 mV/h。这个阈值的合理性取决于材料的本征动力学特征。

NCM811在锂含量较高（x>0.6）时Li层间距较大，扩散较快，弛豫到dV/dt < 0.5 mV/h只需大约30-40分钟。但在低锂含量（x<0.3）时——此时Li层已经收缩，扩散变慢——达到同样的电压稳定性需要2小时以上。按固定弛豫时间来跑的话，低锂区读到的E_s会偏低，算出来的D值随锂含量变化的趋势曲线就可能出现人为的拐点。

团队后来调整了弛豫策略：不用固定时间，而是用自适应判据——连续3次采样间隔（每次10分钟）的电压变化均小于0.2 mV才算弛豫完成。这样在低锂区自动延长弛豫时间，在高锂区不浪费时间。代价是单组GITT测试的总时长从12小时拉到了18小时，但数据的一致性好了一个档次。

D值提取：两种经典公式的选择

有了瞬态电位E_τ和稳态电位E_s之后，锂离子扩散系数的计算通常有两个公式：

简化公式（适用于脉冲时间远小于弛豫时间常数τ ~ L²/D的情况）：

D = (4/π) × (I·V_m / FA)² × (dE/dδ / dE/d√t)²

完整公式（对脉冲时间和弛豫行为没有简化假设）：

需要考虑抛物线型电压响应的完整时间序列。

实际操作中，大部分文献用的是简化公式——因为它只需要读两个电压值和线性拟合dE/d√t的斜率。但如果脉冲时间和弛豫时间常数不满足τ >> t_pulse的条件，简化公式会引入系统偏差。对于NCM这类扩散速率中等的材料，τ大概在10³到10⁴秒，30分钟脉冲（1800秒）的简化公式误差在5-15%之间。

团队最终采用了完整公式处理——多花了一下午做数据拟合，但把双电层充电的初始段（首秒）和非扩散的电压贡献都剔掉了。最终的D值在10⁻¹¹到2×10⁻¹⁰ cm²/s之间，随Li含量变化呈现典型的”M形”趋势。这组数据后来成了团队内部GITT数据处理的标准流程模板，在科研学术网首页上能找到相应的数据处理笔记本和参考文献清单。