多波长-变温 PL 钙钛矿膜厚测量方案：可行性分析

时间：2026-06-09 初稿 — 2026-06-11 文献验证更新 (BR19/BR21/Wong24/DG26) 目标：金字塔绒面硅基底上钙钛矿膜厚的非破坏性光学测量 方案层级：可行性验证 → 理论建模 → 实验设计 文献支撑：12篇WOS搜索+精读论文 + 6篇新精读论文（详见附录B及本文引用）

---

1. 方案架构

核心思想

多波长解决”载流子在哪里产生” → 变温提供独立 α 调制通道 → 比值法消除系统误差

┌─────────────────────────────────────────────────┐
│  激发波长 λ₁ (浅穿透, ~50nm)  ───→  I_PL(λ₁,T₁) │
│  激发波长 λ₂ (深穿透, ~300nm) ───→  I_PL(λ₂,T₁) │  } T₁ = 室温
│  激发波长 λ₁ (浅穿透, ~50nm)  ───→  I_PL(λ₁,T₂) │
│  激发波长 λ₂ (深穿透, ~300nm) ───→  I_PL(λ₂,T₂) │  } T₂ = 升温
└─────────────────────────────────────────────────┘
         ↓ 4个独立测量值 ↓
    → 过确定方程组 → 膜厚 d + 误差评估

测量量与被测量

符号	含义	已知/未知
I_PL(λ,T)	测量的 PL 强度	测量值
α_exc(λ)	激发光吸收系数	✅ 已知 (文献)
α_PL(T)	PL 波长处的吸收系数	✅ 已知 (文献+变温)
η_rad(T)	辐射复合效率	✅ 可通过 Arrhenius 校准
d	钙钛矿膜厚	❓ 待求
G	金字塔几何因子	⚠️ 需射线光学建模

---

2. 物理模型

2.1 单个体素内的 PL 信号

位置 z 处产生的 PL 光子在逃逸到表面探测器途中按 Beer-Lambert 衰减：

dI_PL(z) = Φ_exc × α_exc(λ) × exp[-α_exc(λ)·z] × η_rad(T) × exp[-α_PL(T)·z·G]

其中 G 是金字塔几何因子（平面时 G=1）。

2.2 膜厚内的积分信号

I_PL(λ_ex, T) = Φ₀ × η_rad(T) × α_exc × ∫₀ᵈ exp[-(α_exc + α_PL·G)·z] dz
                ──────────────────  ─────────────────────────────
                   前置因子                  深度积分

解析解：

           Φ₀ × η_rad(T) × α_exc
I = ──────────────────────────────── × { 1 - exp[-(α_exc + α_PL·G)·d] }
         α_exc + α_PL(T)·G

2.3 双波长比值（消除 Φ₀ 和 η_rad 的绝对值依赖）

    I(λ₁, T)   α₁    α₂ + α_PL·G   1 - exp[-(α₁ + α_PL·G)·d]
R = ──────── = ── × ──────────── × ──────────────────────────
    I(λ₂, T)   α₂    α₁ + α_PL·G   1 - exp[-(α₂ + α_PL·G)·d]

Φ₀ 被完全消去。但 η_rad(T) 仅在其空间均匀时被完全消去。 如果 η_rad 存在深度依赖 η_rad(z)（Ahmad et al. 2024 缺陷形成能深度衰减模型），浅激发（405nm）和深激发（640nm）采样的深度加权平均 η 不同，导致残余偏倚。Monte Carlo 模拟（2026-06-11）给出：d_true=500nm、L_surf=30nm 时，忽略 η_rad(z) 的比值法系统偏倚 = +176 nm (+35%)。详见 §5.4 新增的深度非均匀性分析。

2.4 变温引入第二通道

在 T₁ 和 T₂ 下各测量 R₁ 和 R₂：

R₁ = f(d, G, α_PL(T₁))
R₂ = f(d, G, α_PL(T₂))    ← α_PL 变了（因为 Eg 红移）

两个方程解两个未知数 (d, G) → 理论上封闭可解。

变温导致的 α_PL 变化量级：

• ΔEg/ΔT ≈ -0.3 meV/K（钙钛矿典型值）
• ΔT = 60K → ΔEg ≈ -18 meV
• 导致 Δα_PL(780nm) ≈ 5-10% → 可测量

---

2A. 文献支撑与关键假设验证

以下内容基于2026-06-10 WOS系统文献检索（8组搜索，31篇RIS记录）与12篇全文精读。详见附录B。

2A.1 方案核心假设的文献验证状态

关键假设	验证状态	支撑文献	置信度
钙钛矿可在金字塔绒面硅上保形生长	✅ 已验证	Sahli 2018 (Nat. Mater.) 首次展示；Kan 2025 (Nat. Photon.) CuSCN共沉积31.46%；Luo 2023 (Adv. Mater.) >3μm金字塔上28.6%	🔴🔴🔴 极高
金字塔高度范围2-5μm	⚠️ 需修正为3-10μm	Chen 2020 (Joule): 标准工业金字塔3-10 μm	🔴🔴🔴 极高
PL重吸收是绒面PL光谱变化主因	✅ 已由实验证实	Callies 2025 (Solar RRL): −47 meV/μm峰位偏移，排除应力假说	🔴🔴🔴 极高
钙钛矿n,k数据可从文献获取	✅ 已有可靠数据	van Eerden 2017 (Adv. Opt. Mater.): 椭偏+PDS交叉验证MAPI的n,k(300-900nm)	🔴🔴🔴 极高
共聚焦PL可在绒面样品上做3D扫描	✅ 已有先例	Callies 2025: 405nm共聚焦(NA=0.9, 100×), z步进500nm	🔴🔴🔴 极高
PL强度 vs 膜厚呈单调关系	✅ 物理可预期	Beer-Lambert自吸收，Callies Fig.4c&4d PL强度与金字塔高度相关。BR19验证d≥200nm时相干/非相干误差<2%，支持该假设	🔴🔴 高
变温下PL可稳定测量	✅ 已有先例	Fiorentino 2022 (ACS AMI): 85°C in situ PL实时监测2D/3D界面	🔴🔴 高
粗糙表面可用射线追踪建模	✅ 已有方法+工具	BR19 (MC+TMM框架) → RayFlare (Pearce 2021, LGPL v3) → Wong24 (金字塔绒面工业验证)	🔴🔴🔴 极高
工业在线PL检测	❌ 未见报道	He 2023 (Nanomaterials): in situ综述确认仅实验室研究，产线inline空白	🔴 低（opportunity）

2A.2 最直接的文献对标：Callies 2025 精读

Callies et al., Solar RRL 2025, DOI: 10.1002/solr.202500048

• 团队: Fraunhofer ISE + AMOLF (欧洲顶级光伏研究所联队)
• 核心发现: 绒面硅上钙钛矿PL峰位随金字塔高度系统红移（−47 meV/μm at 50% height），主因为光学重吸收（光子循环导致PL光谱中短波长成分被再吸收，等效峰位红移），排除应力假说
• 方法: 405nm共聚焦3D PL + 532nm宏观PL + Raytrace3D模拟
• 金字塔: 6个梯度高度（<1 ~ >6 μm）
• 未覆盖: 仅分析PL峰位，未做定量膜厚反演、未用多波长、未做变温

与本方案的关系（详见 Callies_2025_精读对比分析.md）：

维度	Callies 2025	本方案	差异化
PL分析量	峰位 (meV)	强度 (counts) → 定量d	将物理机制转化为工程测量
激发	单波长 (405nm共聚焦; 532nm宏观)	双波长 (405/635nm)	增加穿透深度维度
温度	室温	变温 (25/85°C)	消除系统误差+α调制
目标	解释重吸收物理机制	定量提取膜厚d	工程应用
模型	复杂Raytrace3D	解析Beer-Lambert + 射线追踪验证	降低计算成本

✅ Callies证明了你方案的前提（重吸收显著），但没有涉及你的创新（比率法、定量膜厚、变温）。本方案可以定位为”Callies工作的工程化延伸”。

2A.3 光学常数数据来源

来源	材料	方法	波长范围	可靠性
Dasgupta/Stranks/Snaith 2026 ⭐	FA₁₋ₓCsₓPb(I₁₋ᵧBrᵧ)₃ 全系列 (Br 0-100%, 5档)	椭偏仪 + 插值	210-2500 nm	⭐⭐⭐⭐⭐ 首选
van Eerden 2017	CH₃NH₃PbI₃	椭偏仪 + 分光光度法 + PDS验证	300-900 nm	⭐⭐⭐⭐⭐
Callies 2025 Supp. Fig.S8	CsFAPb(IBr)₃ (MA-free)	椭偏 + PL反推Urbach尾	full range	⭐⭐⭐⭐

推荐使用 Dasgupta 2026 nk_repo（已克隆至 C:\Users\23012\Desktop\产业项目\nk_repo\）。该数据库覆盖 FA₁₋ₓCsₓPb(I₁₋ᵧBrᵧ)₃ 全系列（Br 0%, 25%, 50%, 75%, 100%），数据格式为 CSV（Stitched_Not_Interpolated/ 原始测量 + Interpolated/ 插值后），直接可用。van Eerden 2017 的 MAPI 数据（经 PDS 验证）可作为交叉校验。

2A.4 粗糙表面光学建模方法（三级方案，经文献更新）

层级	方法	精度	计算成本	来源	适用场景
Level 1	解析 Beer-Lambert + 有效光程因子 G	中	低	本方案核心公式	快速初筛、d≥200nm 验证有效 [BR19]
Level 2	RayFlare RT_TMM（解析 RT + Lambertian 近似）	较高	中	Pearce 2021 (JOSS) + Wong 2024 (IEEE PVSC)	金字塔 G 因子计算、LUT 生成
Level 3	RayFlare 全 MC 光迹追踪 + TMM 验证	高	高	Bonnin-Ripoll 2019 (Sol. Energy Mater.) → RayFlare 实现	关键点精细化校准

Level 1 — Beer-Lambert 的物理有效性（来自 BR19 验证）： Bonnin-Ripoll et al. (2019) [BR19] 通过 MC+TMM 对比实验验证了：当钙钛矿膜厚 d ≥ 200 nm 时，TMM（相干薄膜光学，等价于 Beer-Lambert + 干涉修正）与全 MC（非相干射线追踪）的误差 < 2%。这意味着本方案的 Level 1 解析公式在绝大多数实际膜厚（300-1000 nm）范围内是物理有效的。只有在 d < 200 nm 的极薄膜区才需要升级到 MC 方法。

BR19 的完整 G(z) 公式为 [BR19, Eq.(6)]：

G(z) = α(λ) · ∫ N₀(0,θ) · exp(−αz / cosθ) dθ

其中 N₀(0,θ) 是 MC 光线的表面角分布。当膜厚较薄且金字塔光俘获显著时，G(z) 偏离简单指数衰减（光俘获均质化效应）——这会在方案的多波长比值法中引入额外的深度加权偏倚。

Level 2 — RayFlare 作为实操工具（替代 Kumar 2021 sPSF）：

• RayFlare (Pearce et al. 2021, JOSS, LGPL v3) 已开源并克隆至 C:\Users\23012\Desktop\产业项目\rayflare_repo\
• 提供 RT_TMM 混合方法：专门处理绒面基底 + 薄膜层的光学耦合
• 内置 regular_pyramids() 函数 → 原生金字塔绒面支持
• 已有 perovskite_Si_pyramids_tandem.py 示例脚本 → 可作为模板修改
• Wong et al. (2024) [Wong24] 已用 RayFlare 成功拟合 1000 片金字塔 TOPCon 电池的 EQE（27 LED 波长, 350-1200 nm），验证了其在该几何上的可靠性
• 用解析 RT + Lambertian 近似可大幅加速（Wong 2024 报道 “substantial speed increase”）→ 适用于快速 G 因子 LUT 生成

Level 3 — 全 MC 精化： 当 Level 2 的 G 因子与 FIB-SEM 交叉验证偏差 > 20% 时，升级到 RayFlare 全 MC 模式（n_rays=20000），精化角分布 N₀(θ)。

HTM 层光学效应的处理（来自 BR21）： Bonnin-Ripoll et al. (2021) [BR21] 通过对比 4 种 HTM（Spiro-OMeTAD, P3HT, CuSCN, CuPc）证实：宽禁带 HTM（CuSCN, E_g=3.8eV; Spiro-OMeTAD, E_g≈3.0eV）在 500-850 nm PL 波段近似透明（IQE ≈ 1）[BR21, Fig.3 & Fig.5]。因此，若你的器件使用 CuSCN 或 Spiro-OMeTAD 作为 HTM，在光学建模中可忽略 HTM 层的吸收。若使用窄禁带 HTM（如 P3HT 或 CuPc），则需在 G(z) 计算中加入 HTM 的吸收修正。

---

3. 实验方案设计

3.1 设备需求

设备	用途	替代方案
共聚焦 PL 显微镜	空间分辨 PL 成像	标准倒置荧光显微镜 + 光纤光谱仪
多激光线 (405/532/640nm)	多波长激发	LED 光源 + 带通滤光片
温控样品台 (RT~100°C)	变温测量	Peltier 热台 (TEC)
光谱仪 (500-850nm)	PL 光谱采集	EMCCD 光谱仪

大学实验室可行性：⭐⭐⭐⭐（共聚焦 PL + 热台是半导体物理实验室标准配置）

3.2 测量流程

Step 1: 校准 (平面钙钛矿参考样品)
  ├── 已知膜厚 d_ref (台阶仪/截面SEM确认)
  ├── 测 I_PL(λ₁, T₁), I_PL(λ₂, T₁)
  ├── 测 I_PL(λ₁, T₂), I_PL(λ₂, T₂)
  └── → 提取 G=1 时的校正曲线 R(d)∣G=1

Step 2: 绒面样品测量
  ├── 空间扫描 (40×40 μm² map, 每像素 1 μm)
  ├── 每点测4个 I_PL(λ, T)
  ├── 每点计算 R₁, R₂
  └── → 解出每个像素的 (d, G)

Step 3: 验证
  └── 截面 SEM/FIB 随机 3-5 个点 → 交叉验证

3.3 采集矩阵

┌─────────┬─────────────┬─────────────┬──────────────┐
│ 测量 #  │  λ_ex (nm)  │  T (°C)     │  物理量       │
├─────────┼─────────────┼─────────────┼──────────────┤
│   M1    │   405       │   25 (T₁)   │  I₁₁         │
│   M2    │   635       │   25 (T₁)   │  I₂₁ → R₁    │
│   M3    │   405       │   85 (T₂)   │  I₁₂         │
│   M4    │   635       │   85 (T₂)   │  I₂₂ → R₂    │
└─────────┴─────────────┴─────────────┴──────────────┘
         ↓ 每个空间像素 × 4 次测量 ↓
    → d(x,y) 二维膜厚分布图

---

4. 关键参数估计

4.1 钙钛矿关键参数 (以 CsFAPbIBr ~1.6eV 为例)

参数	值	来源
α_exc(405nm)	~8×10⁴ cm⁻¹	van Eerden 2017 MAPI n,k → 转换
α_exc(635nm)	~1.2×10⁴ cm⁻¹	van Eerden 2017 MAPI n,k → 转换
α_PL(775nm) @ 25°C	~1.5×10⁴ cm⁻¹	van Eerden 2017 + Eg修正
α_PL(790nm) @ 85°C	~1.1×10⁴ cm⁻¹	推算 (Eg红移~18meV, ΔEg/ΔT≈−0.3 meV/K)
η_rad @ RT	0.1-0.9	取决于膜质量
d 范围	0.2-2.0 μm	预期范围（Sahli 2018; Chen 2020 典型值）
金字塔高度	3-10 μm（标准工业）	Chen 2020 (Joule); <1 μm（亚微米工艺可溶液加工）

n,k数据获取方法：椭偏仪 + 分光光度法，用PDS（光热偏转光谱）交叉验证α值（van Eerden 2017）。对非MAPI组分，用Manzoor 2018的波长平移法修正。Callies 2025 Supp. Fig.S8也给出了CsFAPb(IBr)₃的n,k和Urbach尾k值。

4.2 信号灵敏度估算

对于 d = 500 nm 的钙钛矿层：

测量量	值	对膜厚的灵敏度
I(405)/I(635) @ RT	~0.6	R 变化 ~15% per 100nm Δd
I(405)/I(635) @ 85°C	~0.55	R 变化 ~12% per 100nm Δd
ΔR = R₁ - R₂	~0.05	唯一对应一组 (d, G)

→ 在共聚焦 PL 的信噪比下（>100:1），100 nm 的膜厚变化是可检测的。

---

5. 误差源与缓解策略

5.1 系统误差

误差源	影响	缓解
α_exc, α_PL 文献值不确定度	系统偏移	平面样品校准补偿；优先使用van Eerden 2017的PDS验证数据
金字塔几何模型简化	G 估计偏差	三级光学建模验证（解析→sPSF→MC+TMM; 见§2A.4）
光子循环（重吸收→再发射）	等效光程变长，PL峰位红移，有效α_PL降低	Callies 2025已定量：−47 meV/μm峰移，100%归因于光学重吸收（非应力）。本方案通过变温调制可进一步分离出循环的等效G修正。Raytrace3D模拟辅助标定。
热台温度不均匀	局部 T 偏差	IR 相机温度校验

5.2 随机误差

误差源	缓解
激光功率波动	实时功率监测归一化
探测器噪声	长时间积分/多次平均
钙钛矿光致降解	低功率激发 + 光斑快速扫描

5.3 最大不确定度来源：η_rad(z) 的深度非均匀性 ⚠️ 2026-06-11 更新

本节的原始内容（非辐射复合空间分布）已整合到下方 §5.4 的完整分析中。

理论模型最大的脆弱点。η_rad 在缺陷密布的区域（晶界、金字塔峰应变区）会显著降低。更关键的是，2026-06-11 的 Monte Carlo 模拟揭示了一个此前忽略的效应：即使面内 η_rad 均匀，深度方向的 η_rad(z) 非均匀性也会污染比值法。

变温仍然是重要的诊断工具：

变温 Arrhenius 分析 → 提取每个像素的 E_a → 绘制缺陷分布图
→ 将"非辐射热区"从膜厚反演中剔除或加权

5.4 η_rad(z) 深度非均匀性：Monte Carlo 定量分析与修正方案（2026-06-11 新增）

5.4.1 物理根源

Ahmad et al. (2024, Phys. Rev. Materials 8, 125402) 的第一性原理计算证明：

缺陷形成能 DFE(x) = a + b · exp(−x / l)

其中 l 是衰减长度（MAPI 碘空位: ~0.7 Å，影响 ~3-4 nm；CsPbI₃ 缺陷键长偏移影响 ~10 层原子，~20-30 nm）。

DFE 降低 → n_def ∝ exp(−DFE/kT) 增加 → τ_nr ∝ 1/n_def 缩短 → η_rad = τ_r/(τ_r+τ_nr) 指数衰减：

η_rad(z) = η_bulk × [1 − Δη · exp(−z / L_surf)]

其中 L_surf 是有效表面影响层厚度（~20-50 nm），Δη 是表面 η_rad 的相对降低分数（~0.5-0.8）。

5.4.2 为什么比值法会被污染——反直觉发现

传统直觉认为比值法消去了 η_rad 的绝对值依赖，因此对 η_rad 分布不敏感。但这是错的：

激发波长	穿透深度 (~1/α_exc)	采样区域	有效平均 η_rad
405 nm	~67 nm	主要表面层，η_rad 被淬灭	低
640 nm	~400 nm	均匀采样表面对体相	接近 η_bulk

→ 两束光的有效平均 η 不同 → 比值 R = I₄₀₅/I₆₄₀ 被压低 → 均质求解器将压低比值误读为”膜更厚”→ d 系统性高估。

5.4.3 Monte Carlo 定量结果

模拟参数：d_true=500nm, L_surf=30nm, Δη=0.70 (η_surf=0.234), α_PL(300K)=2000cm⁻¹, α_PL(350K)=3200cm⁻¹, f_PR=0.92, G=1.15。10,000 样本，含参数不确定度和测量噪声。

求解方法	噪声-free 偏倚	MC 均值 ± σ	偏倚
比值法 + 均质模型	+176.2 nm (35.2%)	677.4 ± 33.3 nm	+177.4 nm
比值法 + 匹配模型	0.0 nm	512.9 ± 65.7 nm	+12.9 nm
绝对法 + 均质模型	0.0 nm	499.4 ± 16.5 nm	−0.6 nm
绝对法 + 匹配模型	0.0 nm	511.8 ± 16.5 nm	+11.8 nm

模拟脚本：mc_eta_rad_v2.py | 结果图表：mc_eta_rad_full_v2.png | 数据文件：mc_v2_results.npz

偏倚对 L_surf 和膜厚的依赖：

偏倚场景	d (nm)	L_surf (nm)	比值法偏倚	相对偏倚
弱表面淬灭	500	10	+55 nm	+11%
基准	500	30	+176 nm	+35%
强表面淬灭	500	50	+230 nm	+46%
极端	500	100	+310 nm	+62%
薄膜效应	200	30	+105 nm	+53%
厚膜效应	800	30	+215 nm	+27%

5.4.4 修正策略（按推荐优先级）

策略 A：将 (L_surf, Δη) 作为额外求解参数（推荐首选）

修改比值法求解器：4 方程（R₁, R₂, 两个绝对 I_PL 或额外温度点）→ 求解 4 未知数 (d, G, L_surf, Δη)。形式上封闭，但 L_surf 和 Δη 对 R 的敏感度低 → 需要先验约束（Ahmad 文献值 ±50% 作为 Bayesian 先验）。

策略 B：TRPL 先验表征 η_rad(z)（最稳健）

平面钙钛矿 TRPL → τ_eff(T) → Arrhenius → τ_r 分离
  → η_rad(T) 绝对基准
绒面峰区/谷区 TRPL mapping → η_rad 面内分布
Ahmad 2024 → L_surf 先验范围 (10-50 nm)
→ 将 η_bulk, L_surf, Δη 作为已知量输入求解器
→ d 偏倚从 35% 降至 < 5%

策略 C：FIB-SEM 标定偏倚校正曲线（最务实）

即使不物理修正 η_rad(z)，偏倚也是系统性的。通过 FIB-SEM 或双光子 z-scan 做 5-10 点标定 → 建立 d_PL(比值法) vs d_SEM(真值) 的校正曲线 → 软件偏倚补偿 → 精度恢复至 ±10%。

策略 D：切换到单波长绝对法（最简化但需额外归一化）

放弃比值法，每个温度点仅用 405nm 绝对 I_PL 反演 d。405nm 浅激发本身限制了采样深度（~70nm），使有效平均 η 接近 η_bulk，偏倚仅 −0.6 nm。代价是失去了 Φ₀ 的自动消去，需要 (a) 标准荧光参考样归一化 或 (b) 实时功率监测。

---

6. 可行性等级

评估维度	评级	说明
物理基础	⭐⭐⭐⭐⭐	Beer-Lambert + 直接带隙 PL 物理稳固
理论自洽性	⭐⭐⭐⭐	4 方程 2 未知数，过确定。需加入 η_rad(z) 修正（Ahmad 2024 框架）
设备可及性	⭐⭐⭐⭐	共聚焦 PL + 热台，大学实验室标配
信噪比	⭐⭐⭐⭐	100nm 膜厚变化可检测
校准难度	⭐⭐⭐	需平面参考样品 + 金字塔几何建模 + η_rad(z) 表征或标定
最大风险	⚠️	η_rad(z) 深度非均匀性已定量（MC: +35%偏倚）——可修正但需额外实验
综合	⭐⭐⭐⭐	可行，建议推进——但需将 η_rad(z) 表征纳入实验计划

---

7. 推荐实验路线

第一轮：平面校准（2-5天）

• 制备 3-4 个平面钙钛矿样品，膜厚 200/500/800/1500 nm
• 截面 SEM 验证真实膜厚
• 多波长+变温 PL 测量 → 建立校正曲线
• 验证 Beer-Lambert 模型的线性区 → 扁平硅上成立即推广到绒面

第二轮：绒面验证（1-2周）

• 金字塔绒面硅上保形沉积钙钛矿
• FIB-SEM 3-5 点位取截面 → 提供”真实值”
• PL 测量同一点位 → 提取 d_PL
• 对比 d_SEM vs d_PL → 评估几何模型准确度

第三轮：全场扫描（1周）

• 40×40 μm² 2D 扫描 → d(x,y) 厚度分布图
• 统计金字塔峰/谷/斜面的膜厚差异
• 与 Callies 2025 的 PL 光谱结果做交叉验证

---

8. 风险与后备

风险	概率	后备方案
η_rad(z) 深度非均匀性导致比值法偏倚 >20%	高（MC 已定量：35%）	① 求解器中加入 η_rad(z) 模型（Ahmad 2024 框架）；② TRPL 先验表征；③ FIB-SEM 标定偏倚校正曲线；④ 切换单波长绝对法
金字塔光学建模偏差 >20%	中	Kumar 2021 sPSF方法 + Bonnin-Ripoll MC+TMM做G因子验证；引入CL做交叉校准
变温 α_PL 变化量 <5%（测量阈值以下）	低	扩大 ΔT 到 100K 或用更多波长（如增加532nm线）
光子循环使等效光程不可预测	低 ⬇	Callies 2025已证明该效应可定量模拟(Raytrace3D)；变温Arrhenius分离直接PL与循环PL；Monte Carlo光子输运模拟
热致相变/降解	低	限温 <100°C（Fiorentino 2022 已验证85°C可行），氮气环境
金字塔峰区膜厚<100nm导致信号太弱	中	采用更高NA物镜(1.4 oil)增加收集效率；增加积分时间；用谷区数据倒推峰区

---

9. 替代激发方式扩展：双光子激发方案

核心问题：方案 B 最大的脆弱点是”载流子在哪里产生”——多波长穿透深度差被载流子扩散抹平，变温依赖带来热退化风险。双光子激发从物理机制上直接绕过这些问题。

9.1 为什么双光子能解决根本问题

单光子 vs 双光子的本质差异

单光子激发 (方案B):
  405nm 光 → 在材料中按 exp(-αz) 衰减
  → e-h 对产生在整个光路所有深度
  → 无法确定"这个 PL 光子来自哪个深度"
  → 需要多波长 + 变温做间接区分

双光子激发:
  双光子吸收只在焦斑 (~1 μm³) 内发生
  → 焦斑外的材料透明（近红外激发，钙钛矿不吸收）
  → e-h 对只在焦点处产生
  → PL 信号精确地来自已知的深度坐标
  → z-scan 直接给出膜厚

双光子对方案B五个弊端的规避

方案B的弊端	双光子能否规避	机制
① 载流子扩散均匀化	✅ 消除	焦斑 ~1 μm³ 是真正的点激发。扩散仍然发生，但起点精确已知 — z 坐标就是焦平面位置。扩散只影响横向分辨率，不影响 z 定位
② 光子循环	❌ 不可规避	材料本征性质。但可用 z-scan 数据量化循环贡献：焦点在上表面时 I_PL(z) vs 焦点接近界面时 I_PL(z) 的差异可分离直接PL和循环PL
③ 变温热降解	✅ 不再需要变温	z-scan 直接提供深度信息，不需要 α 的温度调制
④ α_PL 变化不足	✅ 不再需要变温	同上
⑤ 非辐射复合分布	⚠️ 缓解	焦斑小于金字塔特征尺寸 → 可在峰区/谷区选择性采集，规避缺陷热区

规避了 4/5 个弊端，唯一不可规避的光子循环也可被量化。

9.2 物理原理

双光子吸收条件

钙钛矿带隙 Eg ≈ 1.60 eV。双光子激发需要 2hν ≥ Eg：

• 激发波长 λ_ex ≤ 2 × 1240/1.60 ≈ 1550 nm
• 常用飞秒激光：1030 nm (Yb:KGW)，~800 nm (Ti:Sapphire)
• 1030 nm 正好在双光子吸收窗口内（2hν ≈ 2.41 eV > 1.60 eV）

关键：1030 nm 光子在钙钛矿中的单光子吸收可以忽略（hν = 1.20 eV < Eg），所以只有焦斑处才有激发。

双光子 PL 的 z 轴分辨率

双光子激发效率 ∝ I² → 信号严格局限在焦斑
焦斑尺寸: Δx × Δy × Δz ≈ 0.3 × 0.3 × 1.0 μm³ (NA=0.9, λ=1030nm)
z 轴分辨率: Δz ≈ 1.0 μm → 对 200-2000 nm 的膜厚，需要亚焦斑定位

膜厚显然小于 z 轴分辨率，但这不是问题——信号峰值位置的提取精度远优于 Δz：

I_PL(z) 峰值拟合精度 ≈ Δz / √(SNR) ≈ 1000 nm / 10 = 100 nm

加上上下界面两个信号变化点（上表面：信号从0→有；下界面：信号从有→0更快衰减），通过双边界定位可进一步提高到 50 nm 量级。

测量逻辑：z-scan 直接膜厚

z-scan 从空气进入钙钛矿:

空气: I_PL = 0
      ↓ 焦点触及上表面
z₁:   I_PL 开始上升        ← 上表面位置
      ↓ 焦点深入钙钛矿
z:    I_PL 持续（探测到体相PL）
      ↓ 焦点接近硅界面
z₂:   I_PL 迅速下降（Si 是间接带隙，PL 效率极低）
      ↓ 焦点进入硅
硅:   I_PL ≈ 0

→ 膜厚 d = z₂ - z₁  （直接，不需要任何反演公式）

在金字塔绒面上的应用

金字塔峰区:
  上表面 = 金字塔顶端（尖）
  界面   = Si 金字塔顶端
  d     = 峰区膜厚（最薄处，~200-300 nm）

金字塔谷区:
  上表面 = 谷底平坦区
  界面   = Si 谷底
  d     = 谷区膜厚（最厚处，可能 >1000 nm）

扫描策略:
  对选定金字塔做离散 z-scan:
    (1) 峰顶: 上下边界 → d_peak
    (2) 谷底: 上下边界 → d_valley
    (3) 斜面: 上下边界 (需几何修正) → d_slope(x)
  → 得到 d(位置) 的离散采样

9.3 实验方案

设备需求

设备	规格	用途
双光子显微镜	倒置或正置，含飞秒激光器	激发源 + 扫描系统
飞秒激光	1030 nm (Yb:KGW) 或 ~800 nm (Ti:Sapphire)，~100 fs，~80 MHz	双光子激发
物镜	高 NA (≥0.9)，长工作距离	聚焦 + 收集
压电 z 扫描台	步进 ≤ 50 nm	精确定位焦平面
GaAsP PMT 或光谱仪	500-850 nm 响应	PL 检测
二向色镜 + 滤光片组	分离激发光 (1030 nm) 和 PL (700-800 nm)	光学滤波

大学实验室可行性：⭐⭐⭐（双光子显微镜不如共聚焦普及，但多数985/211大学有共享平台）

低成本替代方案：如果没有专用双光子显微镜，可用改装方案：

飞秒激光器（共享平台）
  + 自制光路（扩束 → 振镜扫描 → 物镜 → 样品）
  + 收集光路（二向色镜 → 滤光片 → PMT/光谱仪）
  → 总改装成本 ≈ 5-15万元（飞秒激光共享使用费除外）

测量流程

Step 1: 系统校准
  ├── 用已知厚度的平面钙钛矿验证 z-scan 精度
  ├── 确定 z 扫描范围和步进参数
  └── 标定 z 轴漂移（重复扫描同一位置）

Step 2: 绒面 z-scan
  ├── 选择 5-10 个代表性金字塔
  ├── 每个金字塔: 峰顶 × 5 个 z-scan + 谷底 × 5 个 z-scan
  ├── 每条 z-scan 提取上表面和界面位置
  └── 计算局部膜厚

Step 3: 全场扫描（可选，2D z-scan map）
  ├── x-y 网格扫描，每点做完整 z-scan
  ├── 生成 d(x,y) 全场膜厚分布
  └── 统计峰/谷/斜面膜厚差异

Step 4: 交叉验证
  └── FIB-SEM 相同位置取截面 → d_SEM vs d_2P-PL

采集参数

x-y 像素数: 256 × 256 (全场) 或 10-15 个定点 (快速)
z 范围:    5 μm（覆盖金字塔高度 + 膜厚）
z 步进:    50-100 nm
单点积分:  0.1-1 s（取决于 PL 效率和激光功率）
总时长:    定点模式 ~30 分钟，全场扫描 ~6-8 小时

9.4 关键技术问题与解决方案

问题 1：金字塔斜面上的焦斑变形

NA=0.9 的光束在金字塔斜面（~54.7°）上:
  → 折射后焦斑在 z 方向拉长
  → z 定位精度下降

解决：

• 峰顶和谷底区域相对平坦（接近平面），优先做定量测量
• 斜面区域用峰顶/谷底插值，或做几何修正因子
• 如果必需斜面定量，可改用更低 NA 物镜（NA=0.5）增加焦深，减少像差

问题 2：上下表面信号的精确判定

上表面: 信号从 0 → 非0 的跃迁可能有 100-200 nm 模糊带
        （焦斑边缘触及表面时已有微弱信号）

下界面: 硅虽然 PL 极弱，但可能有界面反射回的光子
        被误读为 PL → 界面位置误判

解决：

• 不只依赖强度阈值，同时监控 PL 光谱形状：钙钛矿 PL (~775 nm) vs 硅拉曼/荧光，光谱不同
• 用导数 dI/dz 的最大斜率点作为边界位置（比阈值法精确）
• Si 界面的反射光可通过时间门控（TCSPC）排除（PL 有 ns 级寿命，反射是即时的）

问题 3：激光功率选择

太低:  信号弱，z-scan 信噪比不足
太高:  可能损伤钙钛矿（热效应 + 多光子电离）

参考值：

• 钙钛矿双光子 PL 典型测量功率：0.5-5 mW（样品面）
• 对应焦斑功率密度：~10⁸-10⁹ W/cm²（远低于损伤阈值 ~10¹¹ W/cm²）
• 测试流程：从低功率开始，检查 I_PL ∝ P²（确保是双光子过程，排除热效应）

问题 4：光致卤化物分相

双光子近红外激发（1030 nm）→ 光子能量 1.2 eV < Eg
→ 钙钛矿在激发波长处不吸收
→ 光化学损伤风险远低于单光子 UV/蓝光激发
→ 几乎不存在卤化物分相问题 ✅

这是一个显著优于方案 B 多波长（405 nm 激发）的安全特性。

9.5 误差与精度分析

膜厚测量精度估算

误差源              贡献 (nm)    说明
─────────────────────────────────────────
z 台定位精度        ±10        压电台标准精度
焦斑 z 尺寸         ±50-100    通过峰值拟合减少
上表面判定          ±30-50     导数法定位
下界面判定          ±30-50     光谱辅助去歧义
温度漂移             ±20        恒温控制可忽略
─────────────────────────────────────────
总 RSS 不确定度     ±70-130 nm

→ 对 300-1000 nm 的膜厚，相对不确定度 7-25%。

与方案 B 的精度对比

误差源	方案 B（多波长+变温）	双光子	双光子优势
深度定位	间接（Beer-Lambert反演）	直接（z台位置）	消除反演误差
扩散影响	中等（假设不成立时偏很大）	低（扩散在z方向影响小）	降维
温度相关	有（α(T)不确定度）	无	消除一维误差
几何因子	需建模（5-20%不确定性）	仅斜面需要	峰区/谷区自由
综合不确定度	估计 ±100-200 nm	估计 ±70-130 nm	精度提升 20-60%

9.6 双光子已有文献支撑

虽然直接用于”金字塔绒面膜厚定量”未见报道，但相关技术已有成熟先例：

已有工作	与本法案的关系
de Quilettes et al. (Science 2015): 双光子PL用于钙钛矿晶粒成像	验证了双光子PL对钙钛矿的可行性
双光子z-scan测薄膜厚度 (多种材料体系)	验证了z-scan测厚的基本方法
Walters et al. (JPCL 2015): 钙钛矿双光子吸收截面测量	提供了关键光学参数
Callies 2025: 共聚焦单光子PL的3D扫描方法	提供了绒面PL分析的方法论基础
Fiorentino 2022 (ACS AMI): 85°C in situ PL实时监测钙钛矿界面	证明PL在升温条件下可行
He 2023 (Nanomaterials): in situ/operando综述 (30页)	综述了现有原位PL技术栈
Sahli 2018 (Nat. Mater.): 首次展示钙钛矿在5μm金字塔上保形生长	验证了绒面硅体系的存在和可研究性
Chen 2020 (Joule): 标准金字塔3-10μm, 亚微米金字塔刮涂	提供了绒面工艺基线参数
Kumar 2021 (MRS Adv.): pencil beam ray tracing处理粗糙薄膜	可辅助解释z-scan在粗糙表面的信号混合

→ 双光子方案有技术积累支撑，非空中楼阁。而且S5（WOS搜索：双光子PL钙钛矿）返回365条结果，领域活跃。

9.7 可行性评级

评估维度	评级	说明
物理基础	⭐⭐⭐⭐⭐	双光子吸收定域激发 + z-scan 直接定位
理论自洽性	⭐⭐⭐⭐⭐	z-scan 不需要间接反演，物理上最直接
设备可及性	⭐⭐⭐	双光子显微镜不如共聚焦普及，但共享平台可预约
精度	⭐⭐⭐⭐	预期 ±70-130 nm，优于方案B
规避方案B弊端	⭐⭐⭐⭐⭐	消除扩散/变温/深度反演4个问题
对金字塔适应性	⭐⭐⭐⭐	峰区谷区直接测，斜面需几何修正
光损伤风险	⭐⭐⭐⭐⭐	近红外激发极低损伤
发表新颖性	⭐⭐⭐⭐⭐	双光子测金字塔绒面膜厚，未见报道
综合	⭐⭐⭐⭐½	非常可行，强烈推荐——如果设备可达

9.8 三步方案的关系（更新）

方案 B 原版 (多波长+变温):
  优势: 设备门槛最低
  劣势: 5个弊端，间接反演

方案 A (Callies 扩展):
  优势: 有发表工作验证
  劣势: 依赖模拟工具

双光子方案 (本扩展):
  优势: 物理上最直接，规避4/5个弊端
  劣势: 设备门槛最高

推荐优先级:
  有双光子 → 双光子做"金标准"定量 + 方案A做全场统计
  无公共平台 → 方案A主力 + 方案B增强 + 双光子(外部合作验证)
  仅标准荧光显微镜 → 方案B降配版

---

结论：双光子激发方案是物理上最直接的解决路径——用真正的三维空间分辨取代多波长+变温的间接反演，从根本上规避了方案B最脆弱的4个弊端。如果设备可达，应作为首选的”金标准”方案。次选方案A（Callies扩展），方案B作为设备受限时的后备。

---

10. 从实验室到产线：双光子数据的产业化转化

核心问题：双光子设备昂贵（¥200万+）、操作复杂、速度慢（秒级），不能直接部署到产线。但它的高精度数据可以用来训练一个廉价、快速的代理测量系统——这是一次性的知识转移。

10.1 转化逻辑：金标准 → 代理模型 → 产线部署

┌─────────────────────────────────────────────────────┐
│                   阶段 1：实验室                      │
│                                                     │
│  同一批样品，两种方式同时测量：                         │
│                                                     │
│  ┌──────────────────┐    ┌──────────────────────┐   │
│  │ 双光子 z-scan     │    │ 代理测量 (待选型)      │   │
│  │ → d(x,y) 真值     │    │ → 快速光学信号         │   │
│  │ (金标准,±70 nm)   │    │ (产线级别的速度)        │   │
│  └────────┬─────────┘    └──────────┬───────────┘   │
│           │                         │               │
│           └──── 配对数据集 ─────────┘               │
│                    ↓                                │
│           训练: 代理信号 → d                          │
└─────────────────────────────────────────────────────┘
                         ↓
┌─────────────────────────────────────────────────────┐
│                   阶段 2：产线部署                     │
│                                                     │
│  只部署代理测量设备 + 已训练的校准模型                   │
│  → 实时输出 d(x,y)                                   │
│  → 不需要双光子、不需要飞秒激光                         │
│  → 每台设备成本 <¥5万                                  │
└─────────────────────────────────────────────────────┘

本质：双光子做一次（训练校准模型）→ 模型固化到产线设备里 → 产线永远不需要双光子。与半导体行业的椭偏仪（实验室）→ 反射率膜厚计（产线）是同一逻辑。

10.2 候选代理测量方案

方案 ①：单波长 PL 强度 — 最直接 ⭐ 推荐优先

要素	说明
测量量	单一激发波长（405 nm）的 PL 积分强度
设备	405nm 激光二极管 + 聚焦透镜 + 硅光电探测器 + 长通滤光片
成本	<¥3万
速度	单点 ~0.1 ms
物理原理	PL 强度 ∝ e^(-α_PL·G·d)（Beer-Lambert 自吸收），校准后 I → d

为什么这个方案最合理：

双光子:  像素 #1256 → d = 520 ± 70 nm
同点位 PL强度: I_405 = 1.2×10⁵ counts
双光子:  像素 #1410 → d = 680 ± 65 nm
同点位 PL强度: I_405 = 8.7×10⁴ counts
...
→ 1000+ 个配对点 → 校准曲线 I_405 → d
→ 校准后新样品: 测 I_405 → 模型输出 d ± 80-120 nm

产线硬件清单：

组件	规格	价格 (¥)
405nm 激光二极管	50 mW, 温控	~2,000
聚焦透镜组	光斑 ~5 μm	~3,000
二向色分束器	405R/500-850T	~2,000
长通滤光片	>450 nm	~1,000
Si 光电探测器	放大模块	~3,000
微控制器 + ADC	Arduino/STM32	~500
机械外壳 + XY台	定制	~10,000
总计		~¥21,500

方案 ②：多光谱反射成像 — 信息密度高

要素	说明
测量量	白光照明 → 反射光谱 R(λ)
设备	线光源 + 线扫高光谱相机（400-900 nm）
成本	¥15-30万（工业高光谱相机）
速度	线扫 10-100 mm/s
物理原理	干涉色 → 膜厚，类似 SiO₂ 膜厚测量的成熟技术

优势：

• 每像素有完整光谱 → 信息冗余度高 → 神经网络可学习复杂映射
• 线扫覆盖整片 wafer → 适合大面积在线检测
• 半导体行业已有先例（SiO₂/SiN 膜厚）

与方案①的关系：

• 方案①更低成本，适合定点/多点采样
• 方案②适合全片扫描，检测面积更大
• 双光子数据可同时训练两者的校准模型

方案 ③：双波长反射率比 — 极简化

要素	说明
测量量	双波长 (405/635 nm) 反射光强比
设备	双 LED + 两个光电二极管
成本	<¥5万
速度	单点 ~1 ms
原理	膜厚干涉调制反射率比
局限	金字塔随机性导致单点噪音大，需要多点平均

10.3 推荐产线转化路线

实验室阶段 (1-2月):
  ├── 双光子 z-scan → d(x,y) 金标准 (5-10个金字塔)
  ├── 同一区域采集方案① PL强度 + 方案② 多光谱反射
  ├── 配对: 每个双光子像素的 d vs PL强度/反射光谱
  ├── 训练回归模型 (多项式/神经网络)
  └── → 校准曲线 + 不确定度评估

验证阶段 (2-4周):
  ├── 独立验证样品 (不同批次/组分)
  ├── 代理系统测量 + FIB-SEM 5点交叉验证
  ├── 重复性评估 (同一位置10次测量)
  └── → 生成校准报告和操作规范

产线集成 (1月):
  ├── 选择最终代理方案 (方案① 或 ①+② 组合)
  ├── 工业级封装 (防尘/防振/温控)
  ├── 集成到 MES 系统 (输出 d 值到产线数据库)
  └── → 操作培训 + SOP

10.4 产线与实验室的指标对比

指标	双光子（实验室金标准）	PL代理（产线方案①）	产线需求
单点测量时间	1 s	<10 ms	<100 ms ✅
设备成本	¥200万+	<¥3万	<¥30万 ✅
单点精度	±70 nm	±80-120 nm	±50-100 nm ⚠️ 边界
操作人员	博士生/工程师	产线工人	普通员工 ✅
工作环境	光学平台/恒温	车间	车间 ✅
维护周期	每月	每季度	每半年 ⚠️
可否集成 MES	否	是	是 ✅
采样模式	定点	多点/线扫	全检或抽检 ✅
是否需要液氮	否	否	否 ✅

精度从 ±70 nm 降到 ±80-120 nm，这是校准传递的代价。但产线通常接受 ±10-15% 膜厚的精度——对 500 nm 膜厚即 ±50-75 nm。方案①的 ±80-120 nm 略超此界。改进路径：

(1) 增加采样点数: 峰区+谷区 各5点平均 → 统计降噪 → ±60-90 nm
(2) 引入方案②多光谱补充: 双通道融合 → ±50-80 nm
(3) 同批次多点参考: 用SEM验证的同批次样品做内标 → 系统误差补偿

10.5 校准模型的维护策略

产线运行后，钙钛矿组分或工艺微调可能导致校准偏移。维护方案：

触发条件	维护动作	周期
同批次定期验证	取1片做FIB-SEM测3-5点，对比PL输出	每批次/每周
组分变更	新组分 → 重新双光子校准（一次性）	仅变更时
激光衰减	测标准参考样品，自动温补	每月
工艺窗口漂移 >10%	重新双光子校准	按需（预计2-4次/年）

→ 双光子仅在初始校准和重大变更时使用。产线日常运行完全独立。

10.6 产线可行性总结

评估维度	评级	说明
代理方案精度	⭐⭐⭐	±80-120 nm，可通过多点平均增强到 ±60-90 nm
产线兼容性	⭐⭐⭐⭐⭐	体积小、速度快、免维护、可集成MES
成本	⭐⭐⭐⭐⭐	<¥3万/台，可复制到多条产线
校准维护	⭐⭐⭐⭐	定期FIB-SEM验证，组分变更时重新双光子
产业化风险	🟢 低	半导体行业已有成熟先例（椭偏→反射率膜厚计）
综合	⭐⭐⭐⭐	双光子一次投入，产线长期受益——可行

---

结论：双光子是产线工具的”老师”而非产线工具本身。通过一次性的配对校准（双光子真值 + 廉价代理信号），将实验室精度转化为产线能力。核心代理方案（405nm单色PL强度）成本不足3万元，可在多条产线上复制部署，仅在校准和维护时需要双光子支持。这是典型的”实验室金标准赋能产线”路径。

---

附录A：WOS文献检索概览

搜索日期：2026-06-10 | 数据库：Web of Science All Databases | WOS账号：黄鼎铭/四川大学CARSI

搜索	主题	结果数	核心发现
S1	钙钛矿光学常数 (n,k, absorption coefficient, ellipsometry)	9	van Eerden 2017 MAPI n,k (PDS验证), Manzoor 2018 波长平移法
S2	光子循环/自吸收 (photon recycling, self-absorption, reabsorption)	244	Callies 2025 重吸收定量 (−47 meV/μm)
S3	变温光学性质 (Urbach, temperature-dependent absorption)	171	Urbach能量提取方法, α(T)行为
S4	共聚焦PL膜厚 + 钙钛矿 (confocal, thickness measurement)	12	Callies 2025 — 最直接对标
S5	双光子PL钙钛矿 (two-photon, nonlinear)	365	CsPbBr₃/FAPbBr₃激发光谱，领域活跃
S6	绒面硅钙钛矿涂层形貌	大量	Sahli 2018, Chen 2020, Bonnin-Ripoll 2019等
S7	粗糙面光学建模 (rough surface, scattering, ray tracing, Beer-Lambert)	多篇	van Eerden 2017, Kumar 2021, Bonnin-Ripoll 2019
S8	工业在线PL检测 (inline, in-situ, in-line, monitoring, thickness)	少量	He 2023综述 — 产线inline属明显空白

附录B：全文精读文献清单（12篇）

P0 — 与方案直接相关 ⭐⭐⭐⭐⭐

#	论文	DOI	精读要点	文件
1	Callies 2025, Solar RRL	10.1002/solr.202500048	绒面硅PL重吸收定量, 3D共聚焦, Raytrace3D	✅ 已下载精读
2	van Eerden 2017, Adv. Opt. Mater.	10.1002/adom.201700151	MAPI n,k(300-900nm, PDS验证), TMM+ray-tracing混合	✅ 已下载精读
3	Bonnin-Ripoll 2019, Sol. Energy Mater. Sol. Cells	10.1016/j.solmat.2019.109984	MC ray tracing + TMM, 50-1000nm膜厚扫描	✅ 已下载精读

P1 — 重要支撑 ⭐⭐⭐⭐

#	论文	DOI	精读要点	文件
4	Kumar 2021, MRS Advances	10.1557/s43580-021-00114-0	Pencil beam ray tracing, sPSF, 粗糙薄膜采样体积校正	✅ 已下载
5	Fiorentino 2022, ACS Appl. Mater. Interfaces	10.1021/acsami.2c07122	85°C in situ PL实时监测2D/3D界面演化	✅ 已下载
6	Chen 2020, Joule	10.1016/j.joule.2020.01.008	金字塔3-10μm标准值, 亚微米金字塔刮涂, 保形vs平面化	✅ 已下载
7	Sahli 2018, Nat. Mater.	10.1038/s41563-018-0115-4	首次绒面硅保形钙钛矿, 共蒸发+旋涂, 25.2%	✅ 已下载
8	Jost 2018, Energy Environ. Sci.	10.1039/C8EE02469C	OPTOS形式化, LM foil, 绒面vs平面Jsc差	✅ 已下载
9	He 2023, Nanomaterials	10.3390/nano13131983	In situ/operando综述(30页), 确认产线inline空白	✅ 已下载

P2 — 背景参考 ⭐⭐⭐

#	论文	DOI	精读要点	文件
10	Luo 2023, Adv. Mater.	10.1002/adma.202207883	MACl+MASCN添加剂, >3μm金字塔保形, 28.6%	✅ 已下载
11	Kan 2025, Nat. Photon.	10.1038/s41566-024-01561-5	CuSCN共沉积替代HTL, 31.46%	✅ 已下载
12	Aydin 2023, Nature	10.1038/s41586-023-06631-y	超薄非晶IZO(5nm)解决SAM绒面不均, 32.5%	✅ 已下载

待获取（非必需，已有替代）

论文	DOI	缺的原因	替代方案
Manzoor 2018 (Opt. Express)	10.1364/OE.26.027441	未下载	用van Eerden 2017 MAPI n,k + 波长平移近似
Tucher 2019 (Opt. Express)	10.1364/OE.27.0A1419	未下载	引用Jost 2018的OPTOS方法即可

---

方案文档创建：2026-06-09 初稿 — 2026-06-10 文献验证更新 — 2026-06-11 η_rad(z) MC定量分析更新 文献搜索报告：wos_search_report_2026-06-10.md Callies精读对比：Callies_2025_精读对比分析.md 精读汇总报告：论文精读汇总报告.md MC 模拟：mc_eta_rad_v2.py | 结果图表：mc_eta_rad_full_v2.png | 摘要：mc_v2_summary.txt ArXiv 论文分析：arxiv_analysis_report.md