钙钛矿纳米晶超晶格薄膜的偶极取向反推

Kumar 等（2022）：弱腔薄膜的 p 偏振角分辨光致发光与反推得到的水平偶极占比

作者Luke Cole

Anisotropic nanocrystal superlattices overcoming intrinsic light outcoupling efficiency limit in perovskite quantum dot light-emitting diodes

作者： Sudhir Kumar、Tommaso Marcato、Frank Krumeich、Yen-Ting Li、Yu-Cheng Chiu 等（通讯作者 Chih-Jen Shih）

期刊： Nature Communications 13, 2106（2022） · 对照目标： Figure 3a

发光层里的偶极是躺平还是竖起，直接决定同一套器件叠层能取出多少光。实验上量它的通行做法，是把薄膜贴到半球玻璃透镜上逐角度记录 p 偏振光致发光，再用薄膜光学模型把这条曲线反推回取向占比。

本案例建的是一个两层弱腔模型：18 nm 空穴传输层加 30 nm 发光层，夹在玻璃与空气之间。跑一次角度扫描得到论文 Figure 3a 的曲线，再用图中的 18 个实验点反推水平偶极占比。

Kumar Figure 3a——薄膜 p 偏振光致发光强度随视角的分布，圆圈为实测点，四条计算曲线对应水平偶极占比 0.67、0.72、0.91 和 1.00。Kumar 等，Nature Communications 13, 2106（2022），Figure 3a CC BY

实测点在 0° 归一化到 1，随视角单调下降到 40° 处的极小 0.498，越过玻璃与空气的临界角后跳到 50° 的 1.034，再单调衰减到 85° 的 0.042。水平偶极占比为 1 的那条曲线在 41.9° 处严格归零；占比越低，41° 附近的极小越浅、41° 以上进入基底模式的份额越高。

背景

半球透镜与折射率匹配液消除了基底背面的界面，转台读数因此等于玻璃内的出射角，基底在光学上按半无限介质处理。这样 41° 以上的基底模式区也能被探测到，而常规测量里这部分光被全反射困在基底内。

p 偏振曲线的形状对跃迁偶极矩的空间取向敏感：平行于膜面的偶极把能量集中在小角度，垂直于膜面的偶极把能量送向大角度。41° 附近的尖锐特征来自玻璃与空气界面的全反射，玻璃 $n=1.52$ 对应临界角 41.14°；超过该角度后空气一侧转为倏逝场，水平偶极的 p 偏振发射几乎消失，垂直偶极则在同一角度强耦合进基底。这正是论文用来分辨取向的物理机制。

论文报告的

\Theta_{\mathrm H}

是水平分量占比，软件的 Vertical Fraction 是垂直分量占比，两者相加为 1。论文的 0.72 对应填入 0.28。取用任何文献取向值之前都要先确认它报的是水平占比、垂直占比还是序参量。

结构

叠层只有空气、发光层与空穴传输层，再加上被视为半无限介质的玻璃。没有金属反射镜，属于弱腔体系：结果对绝对效率和发光深度都不敏感，只依赖两层的折射率、厚度和取向占比。

位置	材料	厚度	光学输入
`Incidence Medium`（探测侧）	玻璃	半无限边界	$n=1.52$ ， $k=0$
层 1	X-F6-TAPC 空穴传输层	`18 nm`	逐波长 $n(\lambda)$ ，`528 nm` 处 1.5648； $k=0$
层 2	纳米晶超晶格发光层	`30 nm`	$n$ 见下； $k=0.03548$
`Transmission Medium`	空气	半无限边界	$n=1$ ， $k=0$

把玻璃放在 Incidence Medium 是这个案例的关键一步：角分布的角度定义在入射侧介质内，玻璃在入射侧时视角可以取到 90°，覆盖论文用半球透镜测到的整个基底模式区。把玻璃改放到 Transmission Medium 一侧得到的是空气侧分布，41° 以上没有数据，临界角特征随之消失。

发光层的消光系数由论文给出的 528 nm 光密度 0.011 与 30 nm 厚度反推得 0.0355。发光层折射率则有两个取值，论文在两处分别给出：

取值（`528 nm`）	论文中的出处	有效性
`1.760`	正文与补充材料明写的 $n=1.76\pm0.05$ （`530 nm`），来自椭偏实测	直接给出，波长与本案例的 `528 nm` 相邻
`1.680`	补充材料色散参数在 `528 nm` 求值	该拟合只用了非吸收区 `550–850 nm` 的数据，`528 nm` 属于外推

本案例主结果采用 1.760——论文自己在补充材料中的效率计算用的也是这个值。两个取值都印在论文里，选择其一不引入论文之外的数字；两种取值下的反推结果都会列出。

结构表里只有两行有限厚度薄膜，两侧边界分别是玻璃与空气，发光层折射率为常数 1.76：

论文正文折射率模型的结构页，玻璃入射介质、18 nm X-F6-TAPC、30 nm 发光层与空气透射介质

用补充材料折射率的两个模型只有发光层这一行不同，改为由色散参数求值得到的逐波长数据：

补充材料折射率模型的结构页，发光层改用逐波长折射率

光学与发光设置

五个模型只有两处不同：发光层折射率与 Vertical Fraction。

模型	发光层 $n$ （`528 nm`）	`Vertical Fraction`	$\Theta_{\mathrm H}$
论文取向，正文折射率	1.760	`0.28`	0.72
反推取向，正文折射率	1.760	`0.24655`	0.753
纯水平偶极，正文折射率	1.760	`0`	1.00
论文取向，补充材料折射率	1.68005	`0.28`	0.72
反推取向，补充材料折射率	1.68005	`0.2192`	0.781

其余设置全部共用：

设置	值
发射谱	峰 `527.4 nm`、半高全宽 `23.2 nm` 的高斯谱，`490–570 nm`
偶极位置与分布	相对位置 `0.5`，`Delta`
探测器	Emission 的 `Intensity`（另需勾选一个 Propagation 探测器，本案例用 `R`）
波长模式	`Single`，`528 nm`
角度模式	`Sweep`，`0°–89.5°`，步长 `0.5°`
对照量	`Normalized Angular Distribution` 的 `TM` 曲线

发光体面板应显示 Custom 取向、对应的 Vertical Fraction、相对位置 0.5 和 Delta 分布：

论文取向模型的发光体设置，Custom 取向与 0.28 的垂直占比

角度范围要取到 89.5°，才能覆盖整个基底模式区：

528 nm 单波长、0 到 89.5 度角度扫描的 Intensity 探测器设置

归一化前的角分布对取向严格线性。需要比较一批取向值时，先把 Vertical Fraction 设成 0 与 1 各跑一次，再按占比线性合成即可，不必逐个重跑。

仿真结果与 Figure 3a 对照

论文的原图：

论文四条计算曲线对应水平偶极占比 0.67、0.72、0.91 和 1.00。把这四个取向逐个填进 Vertical Fraction 跑出来，与论文曲线画在同一坐标里：

论文 Figure 3a 的四条计算曲线与实测点，与 Dreapex TMM 角度扫描结果同轴对照 — 虚线空心方块为论文 Figure 3a 的独立读图值，实线为 Dreapex TMM 的真实运行结果，绿色空心圆为论文实测点。论文数据来自 Kumar 等 Figure 3a；Dreapex TMM 数据来自本案例真实计算；图形独立绘制 CC BY 4.0（独立绘制）

四条曲线的走向全部复现：取向占比越低，41° 附近的极小越浅、52° 附近的基底模式峰越高。

水平偶极占比 $\Theta_{\mathrm H}$	与论文同一条曲线的均方根偏差
`0.67`	`0.038`
`0.72`	`0.087`
`0.91`	`0.080`
`1.00`	`0.072`

下面几张是软件里的原始结果。

直接代入论文取向

把 0.72 换算成 Vertical Fraction 0.28 填进去跑一次，曲线在 39.5° 取极小 0.601、在 52.5° 取基底模式峰 1.118：

论文取向、正文折射率下的归一化角分布

放开取向、用 18 个实验点最小二乘反推，最佳值为 0.753（1σ 区间 0.746–0.761），极小移到 41.5° 的 0.547、峰移到 53.0° 的 1.047：

反推取向、正文折射率下的归一化角分布

反推值与论文的 0.72 相差 +0.033，全角度均方根偏差从 5.25 % 降到 3.57 %。

零自由参数的纯水平偶极曲线

把 Vertical Fraction 设为 0 得到纯水平偶极曲线，没有任何自由参数，对应论文 Figure 3a 中 $\Theta_{\mathrm H}=1$ 的那条：

纯水平偶极、正文折射率下的归一化角分布

模型的零点落在 41.5°，论文标注在 41.9°，两者都比 41.14° 的理论临界角偏大，偏移来自偶极到薄膜上表面 15 nm 的传播相位。峰值 0.691 位于 55.5°，在 28°–84° 区间上与论文同一条曲线的均方根偏差为 4.68 %。

换成补充材料折射率

把发光层折射率换成外推得到的 1.68005，其余输入不动。论文取向下的极小抬到 39.0° 的 0.636、峰升到 51.5° 的 1.201：

论文取向、补充材料折射率下的归一化角分布

放开取向反推得到 0.781，此时曲线本身与实验点吻合得很好：

反推取向、补充材料折射率下的归一化角分布

五个模型的特征值汇总。论文的对照值为：实验点极小 0.498（40°）、基底模式峰 1.034（50°）；论文自己的 $\Theta_{\mathrm H}=0.72$ 计算曲线极小 0.523（38.7°）、峰 1.061（51.9°）。

模型	`41°` 附近的极小	基底模式峰	相对 18 个实验点的 RMS
论文取向，正文折射率	`0.601` @ `39.5°`	`1.118` @ `52.5°`	`5.25 %`
反推取向，正文折射率	`0.547` @ `41.5°`	`1.047` @ `53.0°`	`3.57 %`
纯水平偶极，正文折射率	`0.000` @ `41.5°`	`0.691` @ `55.5°`	不适用
论文取向，补充材料折射率	`0.636` @ `39.0°`	`1.201` @ `51.5°`	`8.32 %`
反推取向，补充材料折射率	`0.555` @ `41.5°`	`1.055` @ `52.5°`	`3.10 %`

五个模型的方向都与论文一致：取向占比越低，41° 附近的极小越浅、基底模式峰越高。

发光层 $n$ （`528 nm`）	反推 $\Theta_{\mathrm H}$	1σ 区间	与论文 `0.72` 的偏差
1.760（正文椭偏实测值）	0.753	0.746–0.761	`+0.033`
1.680（补充材料外推）	0.781	0.774–0.787	`+0.061`

反推值几乎完全由发光层折射率决定：折射率每增加 0.02，反推的取向占比降低约 0.007。其余输入都是弱决定量——X-F6-TAPC 厚度 5–35 nm、发光层厚度 24–40 nm、偶极层内位置 0.3–0.7 与谱积分替代单波长四项加在一起只让反推值动 0.02 上下。读图精度也不是瓶颈：单点归一化强度的定位误差 ±0.0036 传播到反推值的 3σ 上限只有 ±0.002。

偏差说明

论文对发光层折射率给出两个不同取值，这是残余偏差的主要来源。对论文色散曲线的像素读数给出 528 nm 处约 1.69–1.72，介于两者之间，而论文未提供该曲线的拟合参数表。论文自己给出的实测带宽是 $1.76\pm0.05$ 。

在 52°–68° 的基底模式区内，论文与本模型的水平通道强度比为 0.918、垂直通道为 0.887，垂直通道多弱 2.5 个百分点。这个不对称使本模型只需较少的垂直分量就能复现同样高度的基底模式峰，反推的取向占比因此系统性偏高。

角分辨光致发光样品的空穴传输层厚度论文未单独标注，本案例取器件值 18 nm；发光层消光系数取由光密度反推的常数，未按波长色散，而 528 nm 紧贴该材料的吸收边。

返回案例复现