埋藏界面分子态的亚单层分辨追踪 Sub-monolayer tracking of molecular states at buried interfaces
ML精修DFT光学模型让RSoXS突破衍射极限，定量解析埋藏界面分子态至亚单层级
前置知识
NEXAFS 与跃迁偶极矩（TDM）
近边X射线吸收精细结构谱（NEXAFS）探测芯层电子向未占据分子轨道的共振跃迁，每个共振对应一个跃迁偶极矩（TDM）方向和振子强度。由于吸收对X射线电场与TDM的相对取向敏感，变角（AR-）NEXAFS 可用来测定分子平均取向；本文用 DFT 计算单分子所有跃迁的能量、TDM 与振子强度，再经机器学习聚类精修后与实验谱定量对齐。
本文光学模型的全部化学指纹和方向信息都来自 TDM 簇，理解 NEXAFS/TDM 才能明白共振能量为何能同时指认化学身份与分子倾角。
RSoXS 反射率与复折射率张量
共振软X射线散射（RSoXS）把光子能量调到碳 K 边附近的电子共振上，测量镜面反射率随散射矢量 q 的振荡（Kiessig 条纹）。材料的复折射率 $n(E)=1-\delta(E)+i\beta(E)$ 决定临界角 $\theta_C=\sqrt{2\delta}$ 与 Brewster 角等特征；对取向分子，$\delta,\beta$ 是各向异性张量，其实部虚部由 Kramers-Kronig 关系耦合。反射率对深度方向的光学剖面敏感，是本文解析埋藏结构的数据来源。
整篇论文的测量对象与正演模型都是偏振、能量分辨的 RSoXS 反射率，不懂张量折射率就无法理解'锁定能量与偏振信息'的含义。
传输矩阵法与层状薄膜拟合
把薄膜近似为若干均匀平板（slab）叠在衬底上，在每个界面用 Fresnel 系数匹配边界条件，通过 4×4 传输矩阵算出总的 s/p 偏振反射率；界面粗糙度以 Nevot-Croce 指数衰减计入。拟合时调节每层厚度、粗糙度和光学常数使计算曲线逼近实验。本文用 refloxide 计算正演、refnx 管理参数与目标函数，差分进化+L-BFGS-B 求全局最优。
'自由模型 vs ML-DFT 模型 vs 梯度模型'的对比都发生在这一正演框架内，参数从 350 个减到 101 个的来龙去脉要靠它理解。
模型选择准则（AIC/BIC/约化χ²）
比较不同复杂度模型时不能只看拟合残差，还要惩罚参数数量：约化卡方 $\chi^2_\nu$ 度量残差水平，Akaike 信息准则（AIC）和贝叶斯信息准则（BIC）在似然基础上各加一个与参数数 $N$ 相关的惩罚项，MSC（模型选择准则）是另一常用指标。值越优说明模型在'拟合好'与'不过拟合'之间取得更好平衡。
本文用这些统计量严格证明 ML-DFT 模型不仅物理上可解释，而且在统计意义上全面优于传统自由张量模型。
Ginzburg-Landau 序参量与梯度模型
二级相变的 Ginzburg-Landau 框架把系统状态描述为序参量场（本文取密度偏差 $\Delta\rho(z)$ 与取向偏差 $\Delta\gamma(z)$），自由能泛函包含局域势能项、梯度惩罚项 $\frac{1}{2}c(\Delta\rho')^2$ 和界面锚定势。对泛函变分得到耦合的二阶微分方程，其解是从界面钉扎态向体相弛豫的指数衰减，特征长度 $\tau$ 即界面层厚度。
论文的梯度结构模型直接由这套物理导出，它替代离散三层板，使亚单层界面剖面的拟合既更优又参数更少。
Fisher 信息与 Cramér-Rao 下界
Fisher 信息度量数据对参数的敏感程度，其逆给出参数估计方差的下界（Cramér-Rao 界，$\sigma_{CR}$）。把 $|\Delta\theta|/\sigma_{CR}$ 作为'可探测度'，可以预先判断某个界面参数在给定测量（能量、偏振、q 点）下能否与体相区分开来。
作者用它定量证明亚单层界面钉扎态以 37σ（取向）和 9σ（密度）置信度可测，这是'突破衍射极限'主张的统计学根基。
研究动机
有机电子器件的核心功能源自器件内部无法直接触及的埋藏异质界面上的分子态：在有机场效应晶体管（OFET）中，紧邻电介质的第一个分子单层直接决定场效应电荷传输；有机太阳能电池与 OLED 中界面取向控制电荷产生与复合；有机自旋阀中单层构象决定自旋杂化；EUV 光刻胶潜像边缘的化学梯度则限定微处理器的最小工艺尺寸。然而现有表征手段均不理想：原子探针层析和截面电子显微术破坏样品且无法分辨键合态与取向；在沉积覆盖率检查点测量的间接方法存在争议——暴露出来的分子层未必代表最终埋藏态。RSoXS 虽无损且对局域化学态敏感，但传统'自由模型'分析需为每个探测能量、每层拟合双折射、二向色性和各向同性 δ、β 共 4 个张量分量，本文的自由模型膨胀到 252 个光学分量、总参数 266 个，拟合不稳定；分析者被迫退化为单一诊断共振，无法同时区分化学种类与取向，更糟的是通常必须假设埋藏界面存在一层未知污染或氧化层才能拟合数据，而无法判断这层是真实存在还是软X射线衍射极限（约数纳米）造成的假象。
本文的目标是本文的具体目标是建立一套非破坏性、定量、可直接解释的分析框架，使共振软X射线散射（RSoXS）能以亚单层分辨率——即突破衍射极限——直接量化埋藏界面处分子的化学身份、取向与堆积密度，并用它回答一个具体的材料科学问题：ZnPc/SiOx 有机晶体管沟道-电介质界面处的分子到底处于什么状态？可量化的子目标包括：（1）把传统自由模型中每能量每层 84 个自由张量参数压缩为每层仅 2 个物理参数（分子倾角 $\gamma$ 与密度 $\rho$），使总参数从 350 个降到 101 个、参数不确定度缩小 10 倍以上；（2）通过 Kramers-Kronig 一致性和化学共振指纹把不同能量、不同偏振的散射特征锁定为一个自洽物理模型，恢复被衍射极限（实空间约 23 Å）掩盖的约 4 Å 尺度界面结构；（3）用所得深度分辨结构裁决长期争论：界面层究竟是污染/氧化，还是本征分子构型；薄膜堆积是被基底图案化外延决定，还是由表面吸附与分子间相互作用的竞争决定。
与已有工作不同的是，本文的独特切入角度是把'光学常数从哪里来'这一问题从拟合器手中彻底接管。此前的 RSoXS 反射率分析要么以经验（Edisonian）方式逐能量自由拟合光学张量，再用二次分析猜测其化学与取向含义——参数自由度过大、不满足 Kramers-Kronig 一致性、且依赖未知的界面污染层假设；要么从 DFT 直接构建光学模型——但受限于基础计算的精度不足。作者团队在前期工作（Murcia et al., Phys. Rev. Lett. 137, 158001, 2026）中证明机器学习聚类精修的 DFT 光学张量能定量复现多晶薄膜的变角 NEXAFS，并编码每个共振的键身份与 TDM 方向；本文首次把这一 ML-DFT 单分子光学张量作为硬约束嵌入 RSoXS 反射率的 4×4 传输矩阵正演中，使不同能量、不同偏振的散射特征通过同一个分子的物理模型相互锁定，从而找回被衍射极限'丢失'的信息。另一个新颖角度在结构模型本身：不用离散三层板，而从 Ginzburg-Landau 二级相变物理导出连续的取向-密度梯度剖面，直接对应生长前沿的物理图像，这在X射线反射率分析中前所未有。
核心方法
方法的直觉是：反射率曲线中每一条 Kiessig 条纹的振幅与相位同时编码了薄膜的几何纳米结构和材料的光学相互作用；如果用第一性原理加机器学习把'光学相互作用'这部分用极少量物理参数（分子倾角、密度）完全固定，那么反射率数据里剩下的自由度就全部归于结构，埋藏界面的亚单层细节就能被'挤出'数据。技术路线分四步：（1）用 StoBe 库对单个 ZnPc 分子做转态势能 DFT 计算，得到 4 个不等价碳激发中心的跃迁能量、振子强度与 TDM 方向；（2）用层次凝聚聚类把跃迁归并为 TDM 簇，再依次精修振子强度、跃迁能量与 TDM 取向，使模型定量复现实测变角 NEXAFS（平均误差 2.7%、最大 5.3%），据此投影出分子坐标系下的复折射率张量 $n(E)=1-\delta(E)+i\beta(E)$，实部 $\delta(E)$ 由 Kramers-Kronig 变换得到；（3）将张量按倾角 $\gamma$ 旋转并做方位平均后，嵌入三层（表面/体相/界面）ZnPc + SiO2/Si 衬底的传输矩阵反射率模型，在 21 个光子能量、s 与 p 两种偏振下对全部数据做单一全局目标联合拟合；（4）进一步用 Ginzburg-Landau 梯度模型替换离散三层板，解出连续的取向-密度深度剖面，并以 Fisher 灵敏度分析和第二片独立薄膜交叉验证其真实性。
核心创新是把全部化学与取向信息'焊死'在光学模型里，而不是留给拟合器自由支配。传统自由模型为每个能量、每层拟合 4 个独立张量分量，21 个能量×3 层共 252 个光学参数：参数爆炸、无物理约束（甚至出现 287 eV 处的非物理二向色性）、无法保证 Kramers-Kronig 一致性，不同能量还会给出互相矛盾的取向结论。ML-DFT 模型则用每层仅有的两个参数，通过旋转-方位平均公式 $n^{xx}_f(E;\gamma,\rho)=\frac{\rho}{2}\left(n^{xx}_m(1+\cos^2\gamma)+n^{zz}_m\sin^2\gamma\right)$ 与 $n^{zz}_f=\rho\left(n^{xx}_m\sin^2\gamma+n^{zz}_m\cos^2\gamma\right)$，从单分子张量确定性生成所有能量、所有偏振的光学响应——一个共振强度的变化必然经 K-K 关系牵动其余能量的 $\delta$，能量间自洽自动满足。配合'双诊断'测量策略：在前边（低于第一个吸收共振）高密度采能定位临界角消光 $\theta_C=\sqrt{2\delta}$，以低辐射剂量读取膜平均取向；再在精细结构共振处采集，由 DFT 指认的 4 个不等价碳环境提供化学身份。前边报告'取向'、共振报告'化学'，两者通过同一分子模型交叉锁定，最终使约 4 Å 的亚单层特征虽远小于衍射极限（$2\pi/q_{max}\approx 23$ Å）仍能通过条纹振幅与相位的微小改变被恢复。
方法步骤详情
第一步：样品制备。ZnPc（C32H16N8Zn，纯度≥99.5%）在 PVD 腔内逐级加热纯化，以 0.2 Å/s 蒸镀约 20 nm 到 UV-臭氧清洗的 Si/SiOx 片上（室温 24 °C），氮气惰性保存送同步辐射。第二步：预表征。VASE 椭偏（65°/70°/75° 三入射角）测得厚度 21.4(5) nm，并确认鱼骨堆叠的劈裂 Q 带（1.98/1.80 eV）；GIXRD 给出面外 (200) d=12.6 Å、面内 (002) 11.8 Å、π-π 堆叠 3.3 Å；AFM 测得 RMS 粗糙度 0.8 nm。第三步：AR-NEXAFS（澳大利亚同步辐射 SXR 线站，总电子产额模式，20°–90° 五个角度、270–350 eV），定出膜平均倾角 74°、约化 χ²=0.213。第四步：DFT+ML 光学张量。StoBe 转态势能计算（HCH 方法给跃迁流形、ΔKS 方法定绝对能量），仅取 4 个碳激发中心；层次凝聚聚类超参数经网格搜索定为振子强度阈值 OS=0.25、重叠百分数 OVP=52.5；TDM 簇精修至复现 AR-NEXAFS（平均误差 2.7%）。第五步：RSoXS 采集（ALS 11.0.1.2 线站）：从碳 K 边下方 250 eV 起步，21 个能量（含 283.7 eV 诊断能量、δ-交叉与前边/后边系列）×s/p 偏振，分段曝光与衰减片拼接出多个量级动态范围的完整反射率曲线，全部数据 12 小时内完成。第六步：全局拟合与验证。refloxide 计算 4×4 传输矩阵反射率（Nevot-Croce 粗糙度），refnx 提供目标函数与先验，差分进化+L-BFGS-B 求最大似然，Hessian 逆矩阵给出不确定度；先做三层板对比（ML-DFT 17 参数 vs 自由张量 266 参数），再换成 Ginzburg-Landau 梯度模型，最后做 Fisher 灵敏度分析与第二片 15 nm 独立薄膜的交叉验证。
技术新颖性
技术新颖性体现在四个层面。其一，光学模型层面：前期 ML-DFT 张量（Murcia et al., PRL 2026）只用于解释谱学测量，本文首次将其作为反射率正演的硬约束嵌入，使每层光学自由度从 84 个骤降到 2 个，三 ZnPc 层描述从 252 个光学分量降到 7 个参数，并用 MSC/AIC/BIC 统计指标严格证明优越性（共振能量平均偏差 1.7σ vs 2.2σ，参数不确定度小 10 倍以上）。其二，信息论层面：把不同能量的共振特征锁定到同一分子的同一组 TDM 上后，原本互相独立的能量通道变成互相约束的证据链——Fisher 信息/Cramér-Rao 分析定量证明埋藏界面钉扎态相对体相的可探测度达 37σ（取向）与 9σ（密度），这在传统逐能量拟合中不可想象。其三，结构模型层面：以 Ginzburg-Landau 自由能泛函（含密度-取向双线性耦合项 $b\Delta\rho\Delta\gamma$）导出从界面钉扎态出发的双指数梯度剖面，替代离散三层板——物理上对应生长前沿的连续再取向，参数更少、拟合更优（消除条纹相位漂移、正确复现 Brewster 角附近结构）。其四，科学结论层面：用统计严格的终态分析推翻了领域内长期默认的'基底图案化外延'与'界面污染/氧化层'两种解释，代之以 2D 表面吸附与 3D 分子间 π-π 凝聚竞争的生长图像，并指出覆盖率检查点式研究测到的只是生长前沿的瞬态而非埋藏终态。
实验结果
核心发现有五点。（1）光学模型验证：ML-DFT 张量复现 AR-NEXAFS 的平均误差 2.7%、最大 5.3%，膜平均倾角 74°（edge-on），约化 χ²=0.213，四个不等价碳环境的精细结构共振全部指认。（2）三层板模型对比：在共振能量系列上 ML-DFT 平均偏差 1.7σ，显著优于自由模型的 2.2σ；描述三 ZnPc 层仅需 7 个光学参数对 252 个，总参数 17 vs 266；参数不确定度小 10 倍以上；MSC/AIC/BIC 全面胜出。自由模型暴露出 287 eV 非物理二向色性、跨能量取向结论矛盾等缺陷；而 ML-DFT 模型在'只允许 ZnPc 存在'的约束下仍拟合成功，直接否定了膜内存在可测量污染或氧化——传统分析归咎于污染的'界面层'实为 ZnPc 自身的独特构型。（3）梯度模型参数（Table I）：总厚 192.7(4) Å、表面粗糙度 6.8(4) Å；体相 $\gamma$=69.6(2)°、$\rho$=1.64(1) g/cm³，与鱼骨晶型文献一致；表面层 $\gamma$=58.7(6)°、$\rho$=2.0(1) g/cm³、$\tau$=13(1) Å（更致密、接近魔角 54.7°、可能为结晶前的无序态）；埋藏界面层 $\gamma$=38(1)°、$\rho$=0.87(7) g/cm³、$\tau$=14.9(1) Å——显著更平躺（face-on）且疏松。（4）亚单层分辨：全膜仅约 15 个分子层，两个界面区各相当于单层厚度，界面钉扎分子高约 4 Å，远小于衍射极限（约 23 Å）；Fisher 分析给出界面态相对体相 37σ（取向，$\sigma_{\gamma}^{Si}$=1.28°）与 9σ（密度，$\sigma_{\rho}^{Si}$=0.084 g/cm³）的可探测度；梯度模型还消除了三层板在 283.7 eV 处的条纹相位漂移并正确复现 Brewster 角结构。（5）生长机制：界面剖面否定基底图案化模型（其预测界面为 edge-on 鱼骨堆积、载流子密度与横向跳跃频率将显著不同），支持表面吸附（favor face-on）与分子间 π-π 内聚（favor edge-on）的竞争；第二片独立薄膜（15 nm）剖面一致，说明界面参数由真空与二氧化硅界面的自由能本征决定。
查看结构化数据
| 任务 | 指标 | 本文 | 基线 | 提升 |
|---|---|---|---|---|
| 光学张量复现变角 NEXAFS（模型验证） | 平均误差 |ε|（最大 εmax） | 2.7%（εmax=5.3%），全部实验特征复现 | 经验/Edisonian 光学模型（需辅助二阶分析，有效性不确定） | 首次实现定量级复现并直接给出化学归属与 TDM 指认 |
| 共振能量系列反射率拟合（三 ZnPc 层结构） | 平均拟合偏差 | 1.7σ（ML-DFT 模型） | 2.2σ（自由张量模型） | 残差更小且总参数从 266 降到 17（>15 倍简约） |
| 模型复杂度与统计优越性 | 总参数数 N 及 AIC/BIC/MSC | 101（含 84 个实验修正项，梯度模型） | 350（自由张量模型） | 参数减少 17 倍，参数不确定度缩小 >10 倍，AIC/BIC 全面胜出 |
| 埋藏界面层分子态分辨 | 可分辨结构尺寸 / 参数置信度 | 约 4 Å 亚单层（界面 γ=38(1)°、ρ=0.87(7) g/cm³，37σ/9σ 可探测度） | 传统衍射极限 2π/qmax≈23 Å | 分辨能力超越衍射极限约 6 倍，进入亚单层区间 |
| 界面层化学归属判定 | 模型约束下的物种判定 | 界面层为 ZnPc 本征构型（单一物种模型拟合成功），否定污染/氧化 | 传统分析需假设未知污染/氧化界面层 | 用化学共振指纹首次无损裁决界面层本质 |
| 体相结构自洽性（与独立技术交叉验证） | 与文献晶体学/MD 值的偏差 | γ=69.6(2)°、ρ=1.64(1) g/cm³，GIXRD d-spacing 吻合约 2% 以内 | 文献鱼骨 ZnPc：72.9(4)°、1.64(3) g/cm³；MD 松弛 72(1)°、1.657 g/cm³ | 体相参数与独立表征一致，支撑界面参数的可信度 |
局限与改进
作者承认的局限：梯度模型本质上是'经验近似'——指数衰减形式虽然连续、参数少、拟合最优，但界面锚定势究竟服从 Dirichlet 还是 Robin 边界条件未做严格判定，这决定了指数幅值与界面势之间的解析关系仍不确定；模型假设单分子 DFT 张量足以代表薄膜全部电子态，忽略了多分子聚集可能产生的新态（本体系中经与实验光学模型比对确认不主导，但对强电子耦合体系未必成立）；MD 可视化中的鱼骨张开角、位置与取向无序仅为视觉连续性设置，非经验或推导因子。我自己的观察：第一，整套流程依赖同步辐射设施——ALS 11.0.1.2 的偏振切换、12 小时专班光束时间、精细能量校准（聚苯乙烯标样），普通实验室无法复制，推广受线站资源限制；第二，方法演示于高度成熟的模型分子 ZnPc（晶体结构已知、DFT 与实验高度吻合、近似单轴），对取向自由度更多、构象可变、多种化学环境重叠的共轭聚合物或混合体系，TDM 聚类的可分性与'每层仅 2 参数'的充分性都待检验；第三，梯度模型假设沿深度单调的双指数过渡且每个界面由单一模式主导，若真实界面存在振荡分层或面内各向异性排列（单轴假设失效），模型会引入系统性偏置；第四，界面层低密度 0.87 g/cm³ 被解释为部分覆盖的平躺分子，但覆盖率并无独立直接测量，该图像与 MD 重构属于自洽性论证而非独立证据。
独立分析的弱点
弱点一：参数-结构简并风险未完全消除。反射率对界面参数的敏感度集中于 p 偏振 Brewster 角附近与 283.7 eV 诊断能量两个通道（Fisher 分析明确显示 s 偏振对界面几乎不敏感），若某体系在这两个通道上灵敏度天然较弱，界面取向与密度参数可能退化或与表面粗糙度交叉耦合。改进方向：对每个新体系先行 Fisher 信息/最优实验设计计算，自适应挑选诊断能量、偏振与 q 采样点组合。弱点二：结论是'终态反推动力学'而非直接观测。本文测量的是生长完成后的冻结剖面，'表面吸附 vs π-π 内聚竞争'是从终态痕迹反推的动力学叙事，缺少沉积过程的实时观测；改进方向是与原位蒸镀 RSoXS 或实时 GIWAXS 联用，直接追踪生长前沿的取向演化。弱点三：体系可推广性存疑。ZnPc 是高对称、刚性、近似单轴的分子，对低对称分子或面内各向异性排列（摩擦对齐膜、单晶），方位平均假设失效，需要展开方位分布函数，参数简约的优势会缩水；改进方向是把单轴旋转-平均推广为完整取向分布（如球谐展开）并与谱学取向重建方法结合。弱点四：低密度界面层的物理解释缺乏独立验证。0.87 g/cm³ 的'部分覆盖平躺分子'图像在光学上与'空隙+空气'难以严格区分，作者用共振各向异性排除了污染，但覆盖率本身未被独立测量；改进方向是用 AFM/STM 对数埃级超薄覆盖薄膜直接成像，或用元素特异的 NEXAFS/同位素标记交叉验证覆盖率。弱点五：误差模型中存在人为设定的经验包络（如 q 依赖不确定度 $\sigma(E,q)=aR_{exp}+bq$，a=0.08、b=0.5×10⁻⁶），不确定度估计的鲁棒性依赖这些选择；改进方向是用重复测量与自助法（bootstrap）经验校准误差模型。
未来方向
作者明确提出的方向：把分析推广到具有更多取向与化学自由度的分子体系——构象变化、轨道杂化影响未占据轨道及其 TDM 的情形；推广到反射率之外的散射几何，实现对不可及埋藏纳米结构的局域分析；应用于自组装单层（作者点名钙钛矿太阳能电池中的分子单层，其取向与覆盖率是效率与热稳定性的决定因素）、有机太阳能电池/OLED 的电荷产生与复合界面、有机自旋阀的自旋杂化界面、电化学器件的离子输运界面。基于本文成果可自然延伸的方向：（1）与原位/operando 蒸镀或溶液加工腔联用，把'终态梯度'变成'动力学电影'，直接检验 2D/3D 竞争生长模型并测量锚定势与特征长度 τ 的时间演化；（2）把 ML-DFT 张量构建流程扩展为可检索的分子光学常数数据库，让新体系无需重做昂贵 DFT 即可套用，形成'分子结构→光学张量→RSoXS 反演'的标准管线；（3）从差分进化+Hessian 的最大似然估计升级为完整贝叶斯后验采样，量化多解风险与参数相关性；（4）把拟合出的钉扎态与 τ 与分子动力学/相场生长模拟闭环，反演界面锚定能与分子间相互作用势，使表征数据直接输入理论模型；（5）拓展到引言中点名的 EUV 光刻胶潜像化学梯度和电池电极-电解质埋藏界面等应用场景；（6）工程化：聚类超参数网格搜索（OS、OVP）与全局拟合已随代码开源，可进一步封装为面向同步辐射普通用户的一站式自动化流水线。
复现评估
复现条件评估总体良好但门槛偏高。代码高度开源：反射率结构分析（wsu-carbon-lab/refl-analysis）、RSoXS 数据归约（als-rsoxs/pxr-reduce）、反射率计算与拟合（als-rsoxs/refloxide）、DFT 聚类与光学模型生成（wsu-carbon-lab/dft-learn）均在 GitHub 公开；NEXAFS 数据已存入 X-ray Atlas（accession 193bc0c5）及 Zenodo（DOI 10.5281/zenodo.14758696、10.5281/zenodo.22742990），GIWAXS/GIXRD 处理用开源 NIKA/Gwyddion，MD 用开源 LAMMPS（GAFF2+UFF 力场、B3LYP/def2-SVP 电荷）。减分项：RSoXS 反射率曲线、非共振 XRR、VASE、AFM 数据及拟合参数仅'合理请求可从通讯作者获得'，未完全公开。硬件需求：关键测量需要同步辐射——ALS 11.0.1.2 线站的偏振切换 RSoXS（一个 12 小时班次可完成全部数据）与澳大利亚同步辐射 SXR 线站的 AR-NEXAFS；DFT 的 StoBe 转态势能计算在 NERSC 超算完成；样品制备为常规 PVD 蒸镀（0.2 Å/s、室温）加标准溶剂清洗+UV-臭氧，难度低。总体难度中等偏高：对拥有同步辐射线站资源的团队，三层板对比与梯度拟合可直接复现——SI 对对准流程、分段曝光拼接、误差建模（Fano 因子、q 依赖包络）、修正参数边界（归一化 0.8–1.2、入射角 ±0.1°）有详尽描述；但复现 Fisher 灵敏度分析需要理解 Cramér-Rao 界与三明治协方差，对使用者的统计功底要求较高。
论文图表
