钙钛矿太阳能电池反向工程 — 系统性研究笔记
作者: 夔嵬
整理时间: 2026-06-04
项目时间: 2026-03-09 ~ 2026-04-09
核心工具: SIMsalabim v2.1(漂移-扩散模拟器)、Python
一、项目总览
本项目的核心目标是建立一套 从 J-V 曲线反推物理参数和退化机制的系统方法。以 FAPI+ALD 钙钛矿太阳能电池为研究对象,通过 SIMsalabim 漂移-扩散模拟器,对数百条 J-V 曲线进行了从宏观(PCE/Voc/Jsc/FF)→ 微观(μn/μp/Nt/Rsh)的深度诊断。
研究流程可分为四层递进:
方法开发(exp1 → exp5) → 找到最优反向工程策略
↓
三级分析体系(原型) → 建立测量→电路→物理的分析框架
↓
老化报告 v1(328 曲线) → 完整老化退化诊断
↓
老化对比 v2(328 vs 172) → 两种老化模式对比
关键数字
| 指标 | 数值 |
|---|---|
| 研究周期 | 2026-03-09 ~ 04-09(30 天) |
| 总 J-V 曲线分析数 | ~550 条 |
| 实验批次 | FAPI+ALD 54 器件 + JV_beginner_guide 基础实验 |
| SIMsalabim eval 次数 | ~5000-8000 次 |
| 最终报告 | 2 份完整 Word 报告 + 25 张分析图表 |
二、研究方法链:从随机搜索到物理约束坐标下降
2.1 问题定义
给定: 某老化阶段的实验 J-V 曲线
目标: 找到一组物理参数 (μn, μp, Nt_bulk, Nt_int, R_shunt),使 SIMsalabim 模拟的 J-V 曲线与实验的 fitError 最小
5 个自由参数及其物理含义:
| 参数 | 符号 | 单位 | 物理含义 |
|---|---|---|---|
| 电子迁移率 | μn | m²/V·s | 电子在钙钛矿层中的迁移能力 |
| 空穴迁移率 | μp | m²/V·s | 空穴在钙钛矿层中的迁移能力 |
| 体相陷阱密度 | Nt_bulk | m⁻³ | 钙钛矿体内 SRH 复合中心密度 |
| 界面陷阱密度 | Nt_int | m⁻² | 电荷传输层界面的缺陷态密度 |
| 并联电阻 | R_shunt | Ω·cm² | 器件的分流/漏电路径 |
文献参考(Heester et al., APL Energy 3, 036106, 2025):
- μn: 1e-6 ~ 1e-3 m²/V·s
- μp: 1e-6 ~ 1e-3 m²/V·s
- Nt_bulk: 1e19 ~ 1e23 m⁻³
- Nt_int: 1e11 ~ 1e15 m⁻²
2.2 五轮实验对比
| 实验 | 方法 | avg fitError | vs 基线 | 物理合理? | 关键发现 |
|---|---|---|---|---|---|
| exp1 | 固定网格搜索(10⁵ 点) | 0.2416 | 基线 | ✅ | 搜索保守,改善有限 |
| exp2 | 多种方法混合 | — | — | ❌ | 收敛不稳定,结果不可靠 |
| exp3 | 网格加密(4 版本迭代) | 0.2376 | -1.7% | ✅ | 改善微小(<2%),性价比低 |
| ⭐exp4 | 坐标下降(无约束) | 0.1425 | -41.2% | ❌ | Late/Final 严重过拟合 |
| ⭐exp5 | 坐标下降(物理约束) | 0.1721 | -29.1% | ✅ | 推荐方案 |
2.3 关键突破:坐标下降 + 物理硬约束
核心思路
坐标下降(Coordinate Descent)的思想是:每次只优化一个参数,保持其他参数不变,一轮一轮迭代。
轮次 1: μn → μp → Nt_bulk → Nt_int → R_shunt
轮次 2: μn → μp → Nt_bulk → Nt_int → R_shunt
轮次 3: ...(直到收敛或达到最大轮次)
每轮在参数空间中用 15 个对数间距采样点搜索最优值。
物理约束(防止过拟合)
exp4 暴露出严重的过拟合问题:
| 阶段 | μn 超范围 | Nt_bulk 超范围 |
|---|---|---|
| Initial | ↑ 47× (4.75e-2) | ↑ 3× |
| Mid | ↑ 34× (3.40e-2) | ↑ 6× |
| Final | ↑ 130× (1.30e-1) | ↑ 32,400× (3.24e27) |
问题: 没有约束时,算法会通过”提高 μn + 疯狂提高 Nt_bulk”来获得低 fitError——参数之间相互补偿,但毫无物理意义。
exp5 的解决方案——硬约束:
- μn ≤ 5e-3 m²/V·s
- μp ≤ 1e-3 m²/V·s
- Nt_bulk ≤ 1e24 m⁻³
- Nt_int ≤ 1e16 m⁻²
- R_shunt ≥ 1e-3 Ω·cm²
最佳实践总结
| 要点 | 说明 |
|---|---|
| 搜索起点 | 必须从已知收敛点出发(随机起点收敛率 0%) |
| 采样密度 | 每轮每参数 15 个对数间距点 |
| 轮次数 | 2-3 轮通常足够 |
| 并行化 | 各老化阶段独立并行,节省 3-4 倍时间 |
| 单次 eval | ~10-12 秒 |
| 总计算量 | ~5000-8000 次 eval |
2.4 V2 映射法的失败(策略对比)
在「混合策略研究」中,尝试了一种替代方案——V2 物理映射法:
V2 映射法: 实验 J-V → 二极管拟合(n, Rs, Rsh, J₀) → 物理公式映射 → (μn, μp, Nt_bulk, Nt_int, R_shunt)
但结果证明:间接映射(电路→物理)在所有 5 个老化阶段均无法超越 Heester 的直接网格搜索:
- V2 avg fitError = 0.3345
- Heester 基线 = 0.2416
这验证了一个重要原则:在 SIMsalabim 反向工程中,直接优化优于间接映射。
三、三级反向工程分析体系
在「反向工程副本」中,建立了从宏观到微观的三级分析框架(原型),后续在两个老化报告中得到完善。
三级结构
Level 1 (测量参数) PCE / Jsc / Voc / FF
↓
Level 2 (电路参数) n / Rs / Rsh / J₀(单二极管拟合)
↓
Level 3 (物理机制) 复合类型 / 损失分解 / 暗电流评估 / SIMsalabim 物理参数
原型结果(102 条曲线)
Level 1 — 测量参数
| 参数 | 最大 | 最小 | 平均 |
|---|---|---|---|
| PCE (%) | 13.54 | 0.00 | 5.64 |
| Jsc (mA/cm²) | 18.03 | 0.00 | 8.81 |
| Voc (V) | 1.045 | 0.004 | 0.588 |
| FF | 0.774 | 0.128 | 0.549 |
Level 2 — 电路参数
| 参数 | 范围 | 平均 |
|---|---|---|
| 理想因子 n | 1.30 ~ 3.00 | 1.60 |
| 串联电阻 Rs | 1.1 ~ 8.9 Ω·cm² | 2.3 |
| 并联电阻 Rsh | 2857 ~ 20000 Ω·cm² | 13539 |
Level 3 — 复合机制分布
| 类型 | 比例 |
|---|---|
| SRH 体相复合 | 90.2% |
| 高陷阱辅助复合 | 8.2% |
| SRH + 界面复合 | 1.6% |
关键发现:SRH 体相复合主导(90.2%),表明体相陷阱是钙钛矿器件的主要复合路径,与后续 FAPI+ALD 分析中 94.3% SRH/界面复合的结论一致。
四、老化退化研究
4.1 v1:完整老化报告(328 条曲线)
核心发现:激进老化模式
| 指标 | 值 |
|---|---|
| PCE 总损失 | 88.1%(~13.5% → ~1.6%) |
| Voc 损失 | -33% |
| Jsc 损失 | -72% |
| Nt_int 变化 | +50 倍 |
| μn 下降 | 96% |
退化特征:
- 复合损耗主导(60-70%):J₀ 增加 → Voc 下降
- 界面崩溃:Nt_int 增加 50 倍 → 载流子传输严重受阻
- 指数衰减动力学:Initial → Mid 损失 ~50%,后期趋于饱和
- 并联电阻 Rsh 下降 96%,分流是主要损耗路径
报告结构(5 章 + 4 附录,11 张图表):
- 实验方法(器件制备、测量条件)
- 结果(11 张图,系统性分析)
- 讨论(退化机制推导 + 局限性)
- 改进建议(工艺/表征/模拟/分析)
- 结论
- 附录(基础方法/计算方法/术语/修改对照)
审核状态: ✅ R4 终审接受,114 项参数、5 张数据表、11 张图全部逐项验证通过。
4.2 v2:对比分析(328 vs 172 曲线)
v2 的核心价值在于对比两种老化模式,从”描述退化”升级为”诊断退化机制”。
两种老化模式对比
| 特征 | 模式 A(激进老化) | 模式 B(温和老化) |
|---|---|---|
| 曲线数 | 328 | 172 |
| PCE 损失 | 88.1% | 33.2% |
| 主导机制 | 复合损耗(60-70%) | 传输+复合均衡 |
| Nt_int | +50× | +2.2× |
| μn 下降 | 96% | 57% |
| n 终值 | 2.5 | 3.0 |
| 退化模式 | 指数衰减→饱和 | 近线性缓降 |
| fitError | 较高(非理想 J-V 形态) | 较低(接近理想行为) |
共性发现(两种模式均验证)
- Nt_bulk 持续递增 → 体相缺陷生成是钙钛矿老化的普遍机制
- μ 持续下降 → 与陷阱散射增强一致
- n 持续增大 → 反映扩散复合→SRH 陷阱辅助复合的转变
- J₀ 增加与 Voc 下降的对数关系 → 验证复合理论框架
差异性发现
- 退化主因:激进 = 复合;温和 = 二者均衡
- 界面崩溃阈值:Nt_int 增加 50× vs 2.2×,存在临界转变点
- 工程意义:
- 温和老化阶段 → 优化传输层(提高 μ)
- 激进老化阶段 → 优先降低陷阱密度、改善界面钝化
4.3 L2/L3 分层拟合(v2 新增)
v2 引入了 L2(电路参数)与 L3(物理参数)的对应关系分析:
Level 2 电路参数拟合:n(t), Rs(t), Rsh(t), J₀(t)
↕ 关联分析
Level 3 物理参数优化:μn(t), μp(t), Nt_bulk(t), Nt_int(t), R_shunt(t)
新增了 6 张补充图:
- S1: 残差阴影图(定位拟合偏差区域)
- S2: 相关矩阵(参数间耦合关系)
- S3: 散点图(L2-L3 对应关系)
- S4: 敏感性雷达图(各参数对 fitError 的贡献)
- S5: 参数摘要表
- S6: 退化贡献分解图
五、技能与工具链
本项目相关的 3 个核心技能(位于 skills/):
| 技能 | 功能 |
|---|---|
| JV Curve Generator | 基于单二极管方程的 J-V 曲线生成器,随机参数+牛顿-拉夫森求解 |
| Solar Cell Physics | 太阳能电池物理知识库,含效率损失机制可视化 |
| Solar Cell Physics Enhanced | 增强版,整合文献理论+J-V 分析工具包 |
核心脚本(位于 scripts/)
| 脚本 | 用途 |
|---|---|
fapi_ald_reverse_engineering_tool_v2.py (23.5 KB) | FAPI+ALD 批量反向工程主工具 |
fapi_ald_generate_physics_report.py (9.6 KB) | 物理报告自动生成 |
fapi_ald_shunt_analysis.py (8.4 KB) | 并联电阻专项分析 |
fapi_ald_rsh_verification.py (8.1 KB) | Rsh=1000 误判验证 |
reverse_engineer_experiment_two_v4.py (9.0 KB) | 基础实验反向工程 v4 |
deeply_reverse_engineer_wpd.py (14.1 KB) | WebPlotDigitizer 数据深度反向工程 |
solar_practitioner.py / solar_team_launcher.py | 太阳能学习协作框架 |
simss_context.md (20.1 KB) | SIMsalabim 使用上下文说明 |
六、文件组织
C:\Users\23012\Desktop\solar\
│
├── 方法开发 ──────────────────────────────
│ ├── method_testing_exp1-5/ exp1→exp5 五轮系统对比
│ ├── hybrid_strategy_simsalabim/ 混合策略 + V2 映射法验证
│ └── reverse_engineering_3levels/ 三级分析体系原型(102 曲线)
│
├── 老化研究 ──────────────────────────────
│ ├── aging_analysis_r1-r4/ v1 完整老化报告(328 曲线)
│ └── aging_analysis_v2_comparison/ v2 对比分析(328 vs 172)
│
├── 技能与工具 ────────────────────────────
│ ├── skills/ JV曲线生成 / 太阳能电池物理 (×2)
│ ├── scripts/ FAPI+ALD 工具链 / 反向工程脚本
│ ├── solar_cell_knowledge/ 知识库
│ └── notes-website/ HTML 笔记(实验报告/文献/可视化)
│
└── 笔记 ────────────────────────────────
└── 钙钛矿太阳能电池反向工程研究笔记.md ← 你在这里
七、关键结论
方法论的
- 坐标下降 + 物理硬约束 = 最优反向工程策略(avg fitError 0.1721,比基线改善 29.1%)
- 间接映射(电路→物理)不可靠,直接优化优于间接推导
- 无约束优化必然过拟合——参数相互补偿制造虚假低误差
- 并行化各老化阶段可节省 3-4 倍时间
物理的
- SRH 体相复合是钙钛矿器件的主要复合路径(≥90%)
- Nt_bulk 持续递增 + μ 持续下降是老化普适特征
- 存在”界面崩溃”临界点:Nt_int 增幅从 2.2× 跳跃到 50×
- 不同老化阶段需要不同的优化策略:温和→传输,激进→复合
研究产出的
- 两份完整老化报告(Word + 25 张图)
- 一套可复用的 SIMsalabim 反向工程工具链
- 五轮实验系统对比 → 方法论证明
整理完成于:2026-06-04 10:20 GMT+8
文件数量:321 文件 / 9 目录
正在载入讨论…