一级衰减与 Arrhenius 从半衰期到加速老化的边界
钙钛矿太阳能电池一文提醒我们,短时效率保持率与多年户外寿命不是同一个指标。若想进一步理解“加速老化”,先要能回答:下降的是哪一个量,速率怎样定义,温度改变后沿用了哪些假设?
本文从最简单的一级衰减出发。所有数值曲线都是明确参数下的教学合成结果,不是材料或器件实测,也没有用几条高温曲线预测真实产品的服役年限。
区分每小时减固定量与减固定比例。
使用开尔文,明确速率常数的时间单位。
同一机制和可比较测试条件需要额外证据。
一级衰减为什么不是匀速下降
设 c(t) 是某种反应物的浓度,温度等条件固定,且满足 dc/dt=−kc,k>0。分离变量并积分:
积分后的关系:ln[c(t)/c(0)] = −kt,因此归一化保持率 **r(t) = exp(−kt)**。
浓度越低,同一时刻消耗的绝对量越少;相同时间段留下相同比例,而非减去相同数量。每小时固定减去 0.01 的线性模型会在足够长时间后变成负数,与这个方程不同。
令 r=1/2,得到半衰期 t₁/₂=ln2/k,与初始浓度无关。IUPAC 的半衰期条目给出了相关术语。这里的积分推导同时说明:只有在一级、恒定 k 的假设下,固定半衰期才自然成立。
将 r 换为阈值 q(0<q<1),得到 t_q=−ln(q)/k。到 80% 的时间是 −ln(0.8)/k,不是半衰期的 80%,也不是效率从 20% 下降到 80%。q 是相对于同一基准的保持比例。
如果把归一化器件输出拟合为 exp(−kt),这是额外的经验建模假设:输出不一定与某一种反应物浓度成正比。不能仅凭看起来像指数下降,就宣布已识别唯一降解反应。
flowchart TD
A[定义被观测量与初始基准] --> B{一级且速率恒定可用吗}
B -- 作为教学假设 --> C[写出 dc除以dt等于负kc]
C --> D[积分得到指数保持率]
D --> E[给定阈值求到达时间]
B -- 实测尚未支持 --> F[比较机制 噪声 与备选模型]
E --> G[模型计算不等于器件寿命认证]
温度怎样进入速率常数
IUPAC 的 Arrhenius 方程定义使用绝对温度 T:k=A exp[−Eₐ/(RT)],在原始形式中将 A 与 Eₐ 视为温度无关参数。R 取 8.314462618 J·mol⁻¹·K⁻¹,Eₐ 用 J/mol。
对同一机制、同一参数适用区间,用参考温度 T_ref 消去 A,计算更方便:
速率比:ln[k(T)/k_ref] = (Eₐ/R) × (1/T_ref − 1/T)。
两温度反算:Eₐ = R × [ln k₂ − ln k₁] / (1/T₁ − 1/T₂)。
当 Eₐ>0、T₂>T₁ 时,分母为正,较高温度的速率更大。摄氏度必须先加 273.15 转成 K;对摄氏温度直接取倒数会改变关系。两个 k 也必须使用相同时间单位。
两点总能确定一条直线,不能据此检验 ln k 对 1/T 是否线性。实际拟合还需要多个温度、重复样品、误差模型和机制证据;温度太接近时,分母很小,测量误差会被放大。
两个温度怎样手算
人为设定 T_ref=300 K,k_ref=ln2/100 h⁻¹,并设 Eₐ/R=3300ln2 K。则从 300 K 升至 330 K:
代入后,ln(k₃₃₀/k₃₀₀)=ln2,所以 330 K 的速率是300 K的两倍。反算 Eₐ=3300Rln2,约 19.02 kJ/mol。
| 温度 | k / h⁻¹ | 50 h 保持率 | 100 h 保持率 | 半衰期 / h | 到 80% 的时间 / h |
|---|---|---|---|---|---|
| 300 K | ln2/100 ≈ 0.00693147 | 2⁻¹ᐟ² ≈ 0.707107 | 0.5 | 100 | 32.1928 |
| 330 K | ln2/50 ≈ 0.01386294 | 0.5 | 0.25 | 50 | 16.0964 |
两温度速率比恰好为 2 是合成参数的设计结果,不是“每升高 30 K 所有材料都快两倍”的经验规律。
温度变化时先累计暴露
如果仍假设只有当前温度影响 k,没有恢复或路径依赖,那么分段恒温的累计暴露量 H=Σk(T_i)Δt_i,r=exp(−H)。每段按自己的速率推进;H 无量纲。
先 300 K 放 50 h,再 330 K 放 50 h,H=0.5ln2+ln2=1.5ln2,因此 r=2^(-1.5)≈0.353553。调换两段顺序在这个模型中得到相同结果。真实材料若出现不可逆相变、裂纹或恢复过程,这个可交换性可能不成立。
把两段简化为“平均 315 K、共 100 h”,一般会得到不同结果,因为 k(T) 非线性。这里逐项比较:分段计算得到0.353553,平均温度近似得到0.369143。后者高估了这组合成条件下的保持率,不能当作精确积分。
flowchart LR A[分段时间与开尔文温度] --> B[逐段计算 k] B --> C[逐段求 k乘时间] C --> D[求和得到无量纲 H] D --> E[保持率等于 exp负H] A --> F[只计算平均温度] F --> G[会丢失非线性温度响应]
把实验条件放回老化曲线
只比较“多少小时后还剩多少”会漏掉重要背景。读图时,至少同时核对下面几项:
- 测的是什么:浓度、光电转换效率,还是相对初值的输出保持率。
- 怎样施加应力:温度和湿度是否恒定,有无光照、昼夜循环或电偏压。
- 样品如何工作:封装状态、负载与最大功率点跟踪条件是否相同。
- 起点与恢复怎么处理:有没有初始快速下降、停光恢复或多阶段变化。
- 观察覆盖多长时间:有多少重复样品,外推部分占结论的多大比例。
哪些结果仍需另外证明
上述手算与曲线已经过独立高精度计算核对。它们回答的是“这些参数与方程会算出什么”,没有验证某种真实材料满足一级模型,也没有据此拟合噪声、给出置信区间或认证户外寿命。
阅读真实光伏稳定性结果时,还应记录光照、湿度、气氛、负载、封装、归一化起点和恢复过程。2020 年 ISOS 共识原文的机构存档强调了测试和报告条件的可比性,并讨论多种应力及可恢复变化;它并不把实验室协议当作户外寿命认证。本文只用该文说明报告边界,没有把它的器件数据拿来拟合上述合成模型。
自检:两条高温曲线拟合很好,就能写成户外可用二十年吗?
不能。还需要确认外推区间的失效机制、湿度和光照等应力、封装与负载条件、恢复和多阶段变化,并处理样品差异及预测不确定性。公式算得准确,只说明在给定假设下计算正确。
回看钙钛矿效率与稳定性,把保持率与绝对效率分开;写读书笔记时,可以用引用溯源与主张边界区分原文报告、自己的计算和外推判断。