GITT计算锂离子扩散系数：脉冲-弛豫数据的Warburg拟合与边界效应修正

GITT计算锂离子扩散系数在锂电材料研发里是评估倍率性能的标配手段。恒电流间歇滴定技术——名字听起来有点复杂，操作上就是”给一小段恒电流脉冲→让电压弛豫回平衡→再给下一段脉冲”。从电压响应曲线里提取扩散系数DLi，物理上走的是Warburg扩散控制的暂态响应公式。但公式推导有三个假设条件，实际测试几乎没有完全满足的，这就导致了D值的数量级偏差。

Weppner-Huggins公式的隐含条件

GITT扩散系数由Weppner-Huggins公式给出：D = (4/πτ)(m_B V_M / M_B S)² (ΔE_s / ΔE_t)²。前提条件有三个：①脉冲足够短，扩散层仅限于颗粒表面；②弛豫足够长，电压回到平衡态；③颗粒内部的锂浓度均匀，脉冲前后的浓度分布可以用Fick第二定律的半无限边界解。

团队在NCM811正极材料的GITT测试中吃过亏——初始脉冲时间设为30 min，每脉冲弛豫4 h。看起来够充分了，但出来的DLi在脱锂初期高达10⁻⁹ cm²/s，和文献报道的10⁻¹⁰-10⁻¹¹ cm²/s差了一个量级。排查后发现问题是脉冲时间太长：30 min的恒电流在二次颗粒尺度（~10 μm）上已经穿透了颗粒表面，扩散层深度超过了D/τ^{-1/2}对应半无限扩散的假设区间。把脉冲时间缩短到10 min、对应电流增大3倍以保持相同SOC步长，D值回归到~3.2×10⁻¹⁰ cm²/s的正常范围。

弛豫区间的选取与边界效应

弛豫区间的终态电压稳定性是第二个关键。GITT计算锂离子扩散系数时，ΔEs是脉冲前后电压的稳态值之差——如果弛豫没回到稳态，ΔEs偏小，D值被高估。NCM811的电压弛豫尾端非常慢：4 h后电压变化率约0.05 mV/min，对高镍材料来说离真正的稳态还有距离。团队把弛豫时间延长到8 h，四个SOC点下的ΔEs相比4 h增大了3-7%，D值对应降低5-12%。

但弛豫时间不能无限拉长——超过12 h后，低电流自放电开始混杂进弛豫响应，U(t)曲线上出现斜率反转。处理方法是在弛豫响应的末端做线性外推：取最后2 h的线性dV/dt趋势，向更长时间做外推，估算真正的稳态电压。这个修正方案虽然是近似，但在NCM811上让GITT得到的D值和EIS低频Warburg阻抗推导的D值一致性从偏差50%缩小到15%以内。

τ^{-1/2}线性区与非线性的分界

公式里还有一个容易被跳过的条件：脉冲期间的电压降ΔEt必须和√τ成正比。GITT公式推导的前提是Fick第二定律的半无限扩散解，在τ½坐标系下，dV/d(t½)是常数。但实际数据绘制在E vs √τ图上，前段和后半段的斜率往往不一致。

NCM811在50% SOC处的脉冲数据揭示：前30 s内dV/d(√τ)稳定在0.82 mV/s½，30 s之后斜率逐渐偏离线性——后半段进入球面扩散域（不再是半无限边界）。物理原因是二次颗粒（~10 μm）由一次颗粒（~1 μm）聚集体构成，锂离子在二次颗粒间隙的扩散对应快响应（前段），在一次颗粒晶格内的扩散对应慢响应（后段）。取哪一段做拟合决定了D值对应的是哪一级扩散通道。团队取前30 s的线性段计算二次颗粒尺度的D，取60-300 s段计算一次颗粒尺度的D，二者相差约5倍。

与EIS的一致校验

单一方法得到的D值必须用另一种独立方法做校验。EIS低频区的Warburg阻抗给出了化学扩散系数，公式D = ½[(V_M/FSA)(dE/dx)]²(1/σ_w)²，其中σ_w是Warburg系数（Zre vs ω^{-1/2}斜率）。NCM811在50% SOC下EIS得到的D = 2.8×10⁻¹⁰ cm²/s，GITT经过上述三项修正后得到D = 3.2×10⁻¹⁰ cm²/s，偏差约14%——在电化学扩散系数测量的可接受范围内（不同方法之间差2-3倍并不罕见）。

GITT计算锂离子扩散系数不是一个套公式的任务，是对恒电流暂态条件下扩散物理边界的判断——脉冲时间、弛豫时长、√τ线性区间，三个边界定义决定了D值的物理含义。锂电材料研发中对不同倍率下D值变化的趋势分析，更多讨论在站内相关案例中有整理。

更多内容请访问 https://www.keyanxueshu.com/