钙钛矿太阳能电池反向工程 — 系统性研究笔记

作者: 夔嵬
整理时间: 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 个自由参数及其物理含义:

参数符号单位物理含义
电子迁移率μnm²/V·s电子在钙钛矿层中的迁移能力
空穴迁移率μpm²/V·s空穴在钙钛矿层中的迁移能力
体相陷阱密度Nt_bulkm⁻³钙钛矿体内 SRH 复合中心密度
界面陷阱密度Nt_intm⁻²电荷传输层界面的缺陷态密度
并联电阻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 fitErrorvs 基线物理合理?关键发现
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.540.005.64
Jsc (mA/cm²)18.030.008.81
Voc (V)1.0450.0040.588
FF0.7740.1280.549

Level 2 — 电路参数

参数范围平均
理想因子 n1.30 ~ 3.001.60
串联电阻 Rs1.1 ~ 8.9 Ω·cm²2.3
并联电阻 Rsh2857 ~ 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 张图表):

  1. 实验方法(器件制备、测量条件)
  2. 结果(11 张图,系统性分析)
  3. 讨论(退化机制推导 + 局限性)
  4. 改进建议(工艺/表征/模拟/分析)
  5. 结论
  6. 附录(基础方法/计算方法/术语/修改对照)

审核状态: ✅ R4 终审接受,114 项参数、5 张数据表、11 张图全部逐项验证通过。

4.2 v2:对比分析(328 vs 172 曲线)

v2 的核心价值在于对比两种老化模式,从”描述退化”升级为”诊断退化机制”。

两种老化模式对比

特征模式 A(激进老化)模式 B(温和老化)
曲线数328172
PCE 损失88.1%33.2%
主导机制复合损耗(60-70%)传输+复合均衡
Nt_int+50×+2.2×
μn 下降96%57%
n 终值2.53.0
退化模式指数衰减→饱和近线性缓降
fitError较高(非理想 J-V 形态)较低(接近理想行为)

共性发现(两种模式均验证)

  1. Nt_bulk 持续递增 → 体相缺陷生成是钙钛矿老化的普遍机制
  2. μ 持续下降 → 与陷阱散射增强一致
  3. n 持续增大 → 反映扩散复合→SRH 陷阱辅助复合的转变
  4. J₀ 增加与 Voc 下降的对数关系 → 验证复合理论框架

差异性发现

  1. 退化主因:激进 = 复合;温和 = 二者均衡
  2. 界面崩溃阈值:Nt_int 增加 50× vs 2.2×,存在临界转变点
  3. 工程意义
    • 温和老化阶段 → 优化传输层(提高 μ)
    • 激进老化阶段 → 优先降低陷阱密度、改善界面钝化

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    ← 你在这里

七、关键结论

方法论的

  1. 坐标下降 + 物理硬约束 = 最优反向工程策略(avg fitError 0.1721,比基线改善 29.1%)
  2. 间接映射(电路→物理)不可靠,直接优化优于间接推导
  3. 无约束优化必然过拟合——参数相互补偿制造虚假低误差
  4. 并行化各老化阶段可节省 3-4 倍时间

物理的

  1. SRH 体相复合是钙钛矿器件的主要复合路径(≥90%)
  2. Nt_bulk 持续递增 + μ 持续下降是老化普适特征
  3. 存在”界面崩溃”临界点:Nt_int 增幅从 2.2× 跳跃到 50×
  4. 不同老化阶段需要不同的优化策略:温和→传输,激进→复合

研究产出的

  1. 两份完整老化报告(Word + 25 张图)
  2. 一套可复用的 SIMsalabim 反向工程工具链
  3. 五轮实验系统对比 → 方法论证明

整理完成于:2026-06-04 10:20 GMT+8
文件数量:321 文件 / 9 目录