方案B不成立时的后备路线

生成日期：2026-06-10
问题：如果载流子扩散跨膜厚（L_d > 700nm）或 η_rad 空间严重不均匀，怎么办？
关联：方案不确定性分析.md §二

---

路线总览

               实验1结果 → 实验2结果 → 最终路线
                    │
        ┌───────────┴───────────┐
        │                       │
   L_d << 700nm            L_d >> 700nm
   (扩散在浅层停止)         (扩散跨过全膜厚)
        │                       │
   ┌────┴────┐             ┌────┴────┐
   │         │             │         │
η_rad均匀  η_rad不均   η_rad均匀  η_rad不均
   │         │             │         │
   ▼         ▼             ▼         ▼
 路线A     路线B         路线C     路线D
(方案B)  (方案B+修正)  (方案B改)  (方案A或双光子)

---

场景1：L_d >> 700nm（扩散跨膜厚）——方案B的物理基础被破坏

为什么致命

方案B的核心假设:
  405nm → 载流子在 z≈0-125nm 产生 (浅)
  635nm → 载流子在 z≈0-833nm 产生 (深)
  → 两个群体来自不同深度 → PL强度差异含深度信息

但如果 L_d > 1μm:
  405nm产生的载流子扩散到全膜厚
  635nm产生的载流子也扩散到全膜厚
  → 两个群体在空间上完全混匀
  → 比值R只反映"哪个波长激发了更多总载流子"
  → 与膜厚d失去关联 ❌

后备方案

备选C1：更多波长 + 更极端的穿透深度差

思路：如果405nm和635nm的穿透深度差不够（700nm），用更极端的波长组合。

激发波长	α (cm⁻¹)	穿透深度 1/α	说明
355 nm (UV)	~2×10⁵	50 nm	极浅穿透，几乎仅在表面
405 nm	~8×10⁴	125 nm	当前λ₁
635 nm	~1.2×10⁴	833 nm	当前λ₂
785 nm (近红外)	~5×10³	2000 nm	带边激发，穿透最深

新的穿透深度差：2000 - 50 = 1950 nm。即使 L_d = 1μm，也还有 ~950nm 的有效深度差。

代价：

• 785nm靠近带边，α对组分/T极度敏感 → 不确定度大
• 355nm可能引发光降解（高能光子）
• 需要更多激光线 → 设备成本/复杂度增加

评级：可行，但降低了方案的优雅性

---

备选C2：时间分辨替代空间分辨（扩散时间 = 深度信息）

思路：载流子在不同深度产生 → 扩散到表面需要不同时间 → PL上升沿包含深度信息。

405nm脉冲激发 (ps/fs):
  → 浅层载流子: 几乎瞬间到达表面 → PL上升沿快 (τ_rise ≈ 0)
  → 深层载流子: 扩散 t_diffuse = L_d²/D  → PL上升沿慢
  
  测量 PL(t) 的上升沿 → 提取 t_rise → 反演产生深度

设备需求：TCSPC（时间相关单光子计数）+ 皮秒脉冲激光。比CW PL复杂，但比双光子便宜。

文献先例：

• 钙钛矿的TRPL已经有大量工作（载流子寿命测量）
• PL上升沿分析在GaAs等半导体中已有成熟方法论
• 尚未在绒面钙钛矿膜厚测量中应用 → 仍然是新的

优点：不依赖α的空间选择性，直接利用载流子输运物理

缺点：需要ps脉冲激光 + TCSPC；信号分析比CW PL复杂一个量级

评级：⭐⭐⭐ 物理上更优雅但实验上更复杂

---

备选C3：PL光谱形状替代PL强度（Callies路线）

思路：放弃强度比，转向PL峰位或光谱形状。Callies 2025已证明PL峰位对膜厚/金字塔高度敏感（−47 meV/μm）。

测量PL光谱（不是积分强度）:
  → 提取峰位 E_peak
  → 提取半高宽 FWHM  
  → 提取短波长侧斜率 (Urbach尾)
  
  E_peak(d) = f(d)   ← 用Raytrace3D或实验校准确定f(d)

优点：

• Callies已为这个路线提供了方法论基础和实验验证
• 峰位不受 η_rad 绝对值影响（只受相对光谱形状影响）
• 扩散不影响PL光谱形状（载流子热化后在带边复合，光谱形状由联合态密度和Urbach尾决定）

缺点：

• 峰的提取精度通常不如积分强度
• 依赖Raytrace3D或等效模拟工具
• 峰位变化量级较小（~10 meV for Δd=200nm），需高光谱分辨率

评级：⭐⭐⭐⭐ 最自然的后备——Callies已铺好了路

---

场景2：η_rad 空间严重不均匀（面内 + 深度两个维度）—— 2026-06-11 更新

2026-06-11 MC 定量更新：10,000 样本 Monte Carlo 模拟揭示了一个反直觉发现：比值法对 η_rad(z) 深度方向的非均匀性极其敏感（而非免疫）。即使面内 η_rad 完全均匀，仅深度方向的 η_rad(z) 衰减就导致 +176 nm (+35%) 的系统偏倚（d_true=500nm, L_surf=30nm）。详见 mc_eta_rad_v2.py 和 mc_v2_summary.txt。

为什么致命（深度维度）

原分析（仅面内）:
  405nm偏向峰区激发 → I₁ ∝ η_rad(峰区)
  635nm偏向谷区激发 → I₂ ∝ η_rad(谷区)
  → R = η_rad(峰)/η_rad(谷) × (深度积分比)
  → η_rad 没有被消去

新增分析（深度维度——Ahmad 2024 物理根源）:
  η_rad(z) 在表面被淬灭（缺陷形成能低），在体相恢复:
    η_rad(z) = η_bulk × [1 − Δη · exp(−z / L_surf)]
  
  405nm 浅穿透 (~67nm) → 主要采样表面低 η 区
  640nm 深穿透 (~400nm) → 均匀采样表面对体相
  
  → 两束光的有效平均 η 不同
  → 比值被压低
  → 均质求解器将"压低比值"误读为"膜更厚，深穿透信号更多"
  → d 系统性高估 10-60%（取决于 L_surf 和膜厚）
  
  **即使面内 η_rad 完全均匀，深度方向差异也会污染比值法。**

Monte Carlo 定量结果

偏倚场景	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%

对比：单波长绝对法（405nm 浅激发）对 η_rad(z) 几乎免疫——噪声-free 偏倚 0 nm，MC 偏倚 −0.6 nm。

后备方案（按优先级）

备选B1：TRPL mapping 作为先验输入（修正，不改方案）

思路：不试图让 η_rad 约分，而是先表征它，再把它作为已知量输入模型。

Step 1: 全场 TRPL mapping
  → 每个像素的载流子寿命 τ(x,y)
  → η_rad(x,y) = τ(x,y) / τ_radiative  (τ_radiative 从文献或平面样品获取)
  → 得到 η_rad(x,y) 二维分布图

Step 2: PL强度测量（与方案B相同）

Step 3: 修正模型
  I(λ₁) = Φ₀ × η_rad(x,y) × α₁ × ∫₀ᵈ exp[-(α₁ + α_PL·G)·z] dz
                            ↑
                     不再是常数，而是每个像素的已知值

Step 4: 反演 d
  η_rad 已知 → 不再需要约分 → d 可直接反演

优点：不需要改变方案B的基本测量流程 缺点：增加TRPL测量步（设备门槛↑，时间↑） 评级：⭐⭐⭐⭐ 如果TRPL设备可达，这是最直接的后备

---

备选B2：变温Arrhenius分离辐射与非辐射通道

思路：利用变温区分 η_rad(T) 的热激活行为。

η_rad(T) = k_rad / [k_rad + k_nonrad(T)]

其中 k_nonrad(T) = k₀ × exp(-E_a/kT)

测量 I_PL(T) 在两个波长下:
  → 拟合 Arrhenius 曲线 → 提取 E_a 和 k₀
  → 推算出每个像素在RT下的 k_nonrad / k_rad 比值
  → 得到 η_rad(x,y) 的相对分布
  → 用于修正 R

优点：利用方案B已有的变温数据，不增加新设备 缺点：

• 需要至少3-4个温度点（不是原来的2个）
• Arrhenius拟合在低温区（RT附近）灵敏度有限
• 如果非辐射通道也是多机制的，拟合可能不稳定

评级：⭐⭐⭐ 加分项，但不宜作为唯一后备

---

备选B3：选择性测量——只在”安全区”取值

思路：既然η_rad在谷区较均匀（应变小、缺陷少），只相信谷区的数据。

空间筛选举例:
  → TRPL mapping 识别 η_rad 高且均匀的区域
  → PL测厚只在这些区域内进行
  → 峰区膜厚通过 SEM 标定的峰/谷膜厚比推算:
      d_peak = d_valley × r_peak/valley
  
  r_peak/valley 来自:
    - FIB-SEM 3-5个金字塔的截面统计
    - 或双光子z-scan一次标定（如果可用）

优点：极其简单，不增加测量复杂度 缺点：

• 只能给出谷区d的统计值 + 峰区的推算值，不是全场的d(x,y)
• 需要SEM或双光子做峰/谷比的单次标定

评级：⭐⭐⭐ 最务实的后备——牺牲空间分辨率换可靠性

---

备选B4：η_rad(z) 深度模型直接嵌入求解器（NEW — Ahmad 2024 + MC 驱动）

思路：在比值法求解器中加入 η_rad(z) 的深度衰减模型作为待拟合参数。

传统比值法（2 未知数）:
  R₁(d, G) = I(405, 300K) / I(640, 300K)
  R₂(d, G) = I(405, 350K) / I(640, 350K)
  → 2 方程解 (d, G)  ✅ 封闭

修正比值法（4 未知数）:
  R₁(d, G, L_surf, Δη) = ...
  R₂(d, G, L_surf, Δη) = ...
  → 2 方程解 4 未知数 ❌ 不封闭

解决方案:
  增加两个新通道:
  (a) 引入第3激发波长（如 532nm）→ R₃ → 3 方程
  (b) 引入绝对 PL 强度（I₄₀₅）作为第 4 方程
  → 4 方程解 4 未知数 ✅ 封闭
  
  或使用 Bayesian 框架:
  L_surf 和 Δη 从 Ahmad 2024 获得先验分布
  → 2 方程 + 2 先验约束 → 可解

所需新增数据：

• Ahmad 2024：L_surf 先验范围（MAPI ~3-30nm 原子级 + CsPbI₃ ~10-50nm 键长偏移级）
• 平面钙钛矿 TRPL：η_bulk 基准值
• 可选：532nm 激光线（¥5,000-10,000 设备增量）

优点：

• 直接在方案 B 的数学框架内解决，不改变测量流程
• 利用已有文献（Ahmad 2024）提供物理先验
• MC 显示修正后偏倚从 35% 降至 <5%

缺点：

• L_surf 和 Δη 对 R 的敏感度较低 → 需要较宽的 L_surf 才能稳定收敛
• Bayesian 先验约束是新的代码工作

评级：⭐⭐⭐⭐ 最完整的理论修正——代价最小（仅增加计算复杂度）

---

备选B5：单波长绝对法降级（简化降维）

思路：放弃比值法，仅用 405nm 浅激发的绝对 PL 强度反演 d。

MC 模拟的意外发现：405nm 浅激发本身将采样限制在表面 ~70nm 内 → 有效平均 η_rad 恰好在 η_bulk 附近 → 偏倚仅 −0.6 nm。

测量: I_PL(405nm, 300K)，I_PL(405nm, 350K)
求解: d 从 I_PL(405nm) 绝对强度反演

代价:
  × 失去 Φ₀ 的自动消去 → 需要:
    (a) 标准荧光参考样（罗丹明薄膜）每次测量前归一化
    (b) 激光功率实时监测 + 归一化
  × 失去双波长深度分辨力
  × 仅表面层对 d 敏感 (d > 200nm 时 I_PL 对 d 趋于饱和)
  
优势:
  √ 对 η_rad(z) 几乎免疫
  √ 设备最简单（仅需 405nm 激光）
  √ 变温维度的贡献保留

精度预估：对 d=200-600nm，绝对法定标后精度 ±10-15%。

评级：⭐⭐⭐ 极简化——但牺牲膜厚测量范围和精度

---

场景3：双杀（L_d 长 + η_rad 不均）——方案B在两个维度上同时失效

唯一出路：不走间接反演，走直接测量

方法	测量量	原理	设备门槛
双光子z-scan	I_PL(z)	焦斑三维定位 → 直接读d	⭐⭐⭐⭐
共聚焦反射z-scan	R(z)	折射率突变定位上下界面	⭐⭐
FIB-SEM	d直接	截面成像	⭐ (破坏性)
台阶仪	Δh	机械探针划过钙钛矿台阶	⭐ (需制备台阶)
椭偏仪	Δ, Ψ vs λ	薄膜干涉反演n,k,d	⭐⭐ (要求平面)

推荐降级路径：

双光子可达 → 双光子z-scan作为金标准
     │
     ├── 只需5-10个金字塔的峰区/谷区离散测量
     ├── 耗时 ~30分钟
     └── 辅助: PL mapping全场统计 (方案A——Callies扩展)

双光子不可达 → 方案A（Callies PL峰位路线）主力
     │
     ├── 测量PL峰位 → 校准曲线 → d
     ├── 用Raytrace3D做正向模拟验证
     └── FIB-SEM 3-5点交叉验证

双光子不可达 + 无Raytrace3D能力 →
     │
     ├── 选择谷区（η_rad较均匀处）
     ├── TRPL mapping 做 η_rad 先验修正
     ├── PL光谱峰位 + 强度双通道
     └── 底线：FIB-SEM做5点，PL做全场统计趋势

---

总后备矩阵（2026-06-11 更新版）

场景	L_d	η_rad 深度	η_rad 面内	首选后备	精度损失	设备增量
理想	✓ 短	✓ 均匀	✓ 均匀	方案B原版	0	无
场景A1	✓ 短	✗ 不均	✓ 均匀	方案B + η_rad(z) 模型修正 (备选B4)	+5-10%	无 (仅计算)
场景A2	✓ 短	✓ 均匀	✗ 不均	方案B + TRPL先验修正 (备选B1)	+15-30%	TRPL
场景B	✗ 长	✓ 均匀	✓ 均匀	PL光谱峰位 (Callies路线)	+20-40%	Raytrace3D
场景C1	✓ 短	✗ 不均	✗ 不均	方案B + η_rad(z)模型 + TRPL (B4+B1)	+15-25%	TRPL
场景C2	✗ 长	✗ 不均	✗ 不均	双光子z-scan 或 Callies路线+TRPL	+30-50%	双光子或TRPL+Raytrace3D
场景D	✗ 长	✗ 不均	✗ 不均	底线：FIB-SEM 5点 + PL全场统计趋势	+50-100%	FIB-SEM
简并降级	✓ 短	✗ 不均	—	单波长绝对法 (备选B5)	+10-15% (d<600nm)	无

新增维度：η_rad 的深度方向非均匀性（Ahmad 2024 指数衰减模型）和面内非均匀性是独立的不确定性源，需要不同的修正策略。MC 模拟证明深度方向的影响可达 35% 量级，不可忽略。

---

最终结论

即使方案B完全不成立（L_d长 + η_rad不均），方案仍然是可行的——只是换了实现路径。

• Callies路线（PL峰位）已经由FhG ISE+AMOLF验证了物理基础，只是精度和定量性不如强度比
• 双光子路线是物理上最直接的，精度最高，只是设备门槛最高
• TRPL修正是成本适中的中间路线

三者不互斥。推荐的策略是：

      双光子 (一次, 共享平台)
           │
           ├→ 金标准 d 数据 (5-10个金字塔)
           │
     ┌─────┴─────┐
     │           │
  PL峰位      PL强度
(Callies)    (方案B)
     │           │
     └─────┬─────┘
           │
    产线代理模型 (405nm单色PL)

双光子做一次性的”真相来源”，两条PL路线（峰位+强度）做主力测量，产线代理做日常。

这比押注单一方案稳健得多——你的方案不会因为实验1或实验2的结果而”失败”，只会在路线树的不同分支上前进。

---

生成时间：2026-06-10 14:30 GMT+8 — 2026-06-11 10:30 GMT+8 更新：η_rad(z) 深度非均匀性 MC 定量分析 关联：方案不确定性分析.md MC 模拟：mc_eta_rad_v2.py | 图表：mc_eta_rad_full_v2.png | 摘要：mc_v2_summary.txt