本文是对 arXiv:2609.19105《Characterization of the CSST Survey Camera CCDs: I. basic electro-optical performance》的精读记录,同时整理了整套探测器基础电光性能表征的方法与判据口径,供课题组内的红外/光学焦面探测器验收与标定参考。
原文信息:罗遵、詹虎等(中国科学院国家天文台 空间天文与技术重点实验室;中国科学院大学;北京大学科维理天文与天体物理研究所),arXiv:2609.19105v1 [astro-ph.IM],2026-09-16,22 页 17 图,投 RAA。
一、这项工作在做什么
CSST(中国空间站巡天望远镜)是 2 米口径宽视场望远镜,与中国空间站同轨,计划执行 17,500 deg² 的成像与无缝光谱巡天,覆盖 255–1000 nm,配置九个测光带(NUV、u、g、r、i、z、WV、WI、v)与三个光谱带(GU、GV、GI)。
本文是该系列的第一篇,对 CSST 巡天相机(CSST-SC)的 31 片飞行型号 CCD 做地面实验室表征,只报告基础电光性能。所谓"基础",论文界定得很明确:只用偏置帧、平场、暗帧三类最简单数据就能得到的那部分。
它对应的是探测器的光→电转换与电荷处理性能,与"光学性能"(波前误差、PSF、畸变、杂散光,用干涉仪或星点投影测)分属两条线。
哪些不属于基础电光性能(需要专门成像装置,属后续论文):串扰、像元内响应非均匀性、电荷扩散,另加辐照退化。测试台上有一段可拆挡管,拆掉后电动平移台可把成像测试装置推入光路——这条"分支"就是标题里那个"I."的物理含义。
摘要给出七项关键结果:
- 满阱容量 $\ge 80\ \mathrm{ke^-}$(受测试电子学限制)
- $7\text{–}70\ \mathrm{ke^-}$ 范围内非线性 $< 1.4\%$
- 读出噪声中位数约 $3.6\ \mathrm{e^-}$
- 工作温度 $-85\ ^\circ\mathrm{C}$ 下暗电流 $\le 0.02\ \mathrm{e^- \, pix^{-1} \, s^{-1}}$
- 在 255–1000 nm 内量子效率高
- 典型像元响应非均匀性(PRNU)$< 1\%$
- $10\ \mathrm{ke^-}$ 平场水平下电荷转移效率(CTE)$> 99.9994\%$
这七项要分三类读:结果值(读出噪声、暗电流、PRNU、CTE)可直接引用;判据值(满阱、非线性)必须先知道定义才有意义;条件性表述(满阱"受测试电子学限制"、QE 只说"高")都带限定语。
二、焦面布局与器件命名
焦面由 31 片 9k × 9k CCD 拼成,编号 CCD00–CCD30,来自两个厂商、五个子型。论文为每个波段选了在该波段量子效率最高的子型——这是因为没有一种器件能覆盖全波段。
图注里有一处值得留意:三个无缝光谱波段的 CCD 与光栅,在焦面左右两侧彼此相对旋转了 2°。正文未提、也未给出理由。从无缝光谱的几何约束看,最可能是要让色散方向与探测器的读出通道边界相容(每条光谱若斜跨通道边界会引入增益不连续与伪影);图注说的是"左右两侧"而非"全部",也支持这是局部排布而非全局色散方向的考虑。引用时须分开讲:这是设计事实,原因属推断。
五种子型的命名
| 代号 | 厂商 | 硅型 | 涂层 | 装机 | 强项波段 |
|---|---|---|---|---|---|
| 44NUV | 中电 44 所 | 标准硅 | NUV 涂层 | 是 | 250–400 nm(QE 74–76%) |
| 44VIS | 中电 44 所 | 标准硅 | VIS 涂层 | 是 | 400–700 nm(QE 86–94%) |
| e2vNUV | Teledyne e2v | 标准硅 | NUV 涂层 | 否(只测不装) | 250–450 nm |
| e2vVIS | Teledyne e2v | 标准硅 | astro-multi 2 | 是 | 400–700 nm(QE 89–93%) |
| e2vDD | Teledyne e2v | 深耗尽硅 | astro-multi 2 | 是 | 700–1000 nm(850 nm 达 77–80%) |
读法:44 与 e2v 是并列的厂商代号,后缀 NUV/VIS/DD 指涂层或硅型。
标准硅 vs 深耗尽硅是理解全篇的钥匙。 标准硅感光层只有几到十几微米厚,而 850 nm 的光子在硅中要穿透 15–20 µm,未及收集就损失掉了,所以长波 QE 断崖下跌。深耗尽用高电阻率衬底加反向偏压,把耗尽区延伸到几十微米深,长波光子直接产生在耗尽区内、靠电场漂移收集,850 nm 因此能到 77–80%。代价是紫外:光子在表面产生、要穿过更厚的死层,e2vDD 在 250–400 nm 只有 31–42%,远低于 44NUV 的 74–76%。
深耗尽用紫外端的牺牲换近红外的翻身。 这一条还顺带解释了两处正文:为什么深耗尽的亮-胖效应更严重(衬底厚、电子收集路径长,与已积累电荷的横向相互作用时间更长),以及为什么相互作用系数 $a_{ij}$ 随波长下降(长波吸收更深、相互作用时间更短)。
单器件结构:16 个读出通道
单片 9216 × 9216 像元、像元 10 µm 见方、16 个读出通道、每通道含 27 个 prescan 像元。CCD290-99 每通道 4616 × 1152;GM1201A 的通道 1、8、9、16 为 4620 × 1154,其余为 4620 × 1152。两者差异是 CCD290-99 有 dummy 输出、GM1201A 没有。
这件事决定了全文的组织方式。 论文明确要求:转换增益必须为每片 CCD 的每个通道单独标定。16 个通道各有独立的读出电路,增益各不相同,所以电子单位的换算必须在通道一级做。这就是为什么后面所有结果都先按通道给出、再按子型汇总成区间。
三、测试台
光路是:光源(激光驱动)→ 单色仪 → 快门 → 光纤 → 积分球 → 挡板 → 杜瓦窗口 → CCD。另有两条控制/监测通道。
通道一:时间基准。 积分球上挂监视光电二极管,配皮安计以 70 Hz 采样。指标是:单点噪声 $< 0.2\%$,曝光长于 0.1 s 时平均电流噪声 $< 0.1\%$。
这是全篇的精度天花板。 论文要报零点几个百分点量级的非线性,参考通道就必须比它更稳。反过来也成立:要测多小的非线性,参考通道就必须有多稳——这条原则可以直接搬进我们自己相机的测试规格。
通道二:温控。 CCD 温度由加热其安装板来控制,稳定在 $-85\ ^\circ\mathrm{C}$,约 10 小时工作期内波动 $\pm 0.1\ ^\circ\mathrm{C}$;暗电流另测 $-95\ ^\circ\mathrm{C}$ 与 $-75\ ^\circ\mathrm{C}$ 两个温度点。
其余关键条件:积分球面内均匀性 $< 2\%$(照明不均与器件 PRNU 在空间频率上重叠且无法完全解耦,装置越均匀解耦越干净);背景 $< 10^{-3}\ \mathrm{e^- \, pix^{-1} \, s^{-1}}$(比暗电流要求低两个数量级,才不污染暗电流测量);读出参数用在轨默认设置(像元速率 150 kHz、行转移时间 200 µs、全帧读出约 39 s)。
顺带算清一个容易误读的地方:150 kHz 是每通道串行寄存器的像元读出速率,读出一行 1152 个像元需 $1152/150000 \approx 7.7\ \mathrm{ms}$,4616 行共约 35.5 s,加上并行转移与行开销即约 39 s。而 200 µs 是行转移时间,与串行像元读出交替进行,不是那 39 s 的组成部分——两个时钟不能相乘混算。
论文有一句关键限定:非线性与满阱都受读出电路影响,最终结果必须在焦面接上真实飞行电子学之后重测。所以满阱"受测试电子学限制"、预期换飞行件后 $> 90\ \mathrm{ke^-}$,这个外推只对读出侧的非线性合理,对像元侧的本征非线性不成立。
四、七项测量
4.1 PTC 与转换增益
PTC 是像元值方差对均值的关系。假设泊松统计,斜率为
$$V = \frac{\alpha}{G} F \quad \Longrightarrow \quad \text{斜率} = \frac{V}{F} = \frac{\alpha}{G} $$
其中 $V$ 是仅由光子散粒噪声引起的方差,$F$ 是平均信号,$\alpha$ 是量子产额(每个光子产生的平均光电子数),$G$ 是转换增益(单位 $\mathrm{e^-/ADU}$)。
为什么斜率不等于 $1/G$。 因为要得到 $G$ 必须已知 $\alpha$;而独立测 $\alpha$ 需要绝对电子数,那又需要 $G$——这是个闭环。所以 $\alpha = 1$ 只能是物理前提,不能是测量结果。论文的做法是:在足够长的波长处标定增益(标准硅用 650 nm,e2vDD 用 850 nm),因为光子能量降低时 $\alpha$ 渐近趋于 1,之后量子产额由 PTC 斜率与增益的乘积给出。
$\alpha > 1$ 的物理是:一个被吸收的光子激发出多个电子-空穴对。硅的带隙约 1.1 eV,若光子能量超过约 $2E_g$,激发出的"热"电子仍有足够动能通过碰撞电离再激发第二对。250 nm 的光子有 4.96 eV(约 4.5 倍 $E_g$),余量很大;400 nm 是 3.10 eV(约 2.8 倍),余量已经不多,再往长波走 $\alpha$ 迅速收敛到 1。
论文还强调第二层依赖:量子产额受硅内电场强度影响,而电场强度取决于背照减薄与钝化工艺,因此它在不同子型间差异显著。强电场中光电子可在两次碰撞之间被电场加速、额外获取动能,相当于降低碰撞电离的有效阈值;而短波光子恰好被吸收在紧邻背面偏压区的最强场里,长波光子则进入场更弱的深处。两个机制同向叠加,所以 $\alpha$ 单调下降并趋于 1。图 6 中 44NUV 在 250 nm 处 $\alpha$ 最高达约 1.28。
必须用成对平场的差分图像。 单帧方差由散粒噪声、读出噪声与 PRNU 共同组成,而 PRNU(约 1%)在 $45\ \mathrm{ke^-}$ 平场下比散粒噪声($1/\sqrt{45000} \approx 0.47\%$)还大。两帧相减,跨帧不变的 PRNU、坏点、照明不均全部抵消;又因
$$\mathrm{Var}(A-B) = \mathrm{Var}(A) + \mathrm{Var}(B) = 2\sigma^2 $$
单帧的散粒噪声方差就是差分方差的一半——这正是论文里相关矩阵元素"取中位数的一半"的来源。此外还须用第 0 档零曝光帧扣掉读出噪声的方差贡献。
亮-胖效应:公式一行,改正好几段。 像元里积累的电荷产生的横向电场会把后续光电子推向邻像元,使高流量端 PTC 向下弯。论文按 Astier 等(2019)的处理,用相关矩阵元素之和代替方差,把被推到邻像元的电荷协方差加回来。经验是标准硅取三邻域、深耗尽 e2vDD 取五邻域(后者衬底厚、相互作用更强)。
4.2 非线性与满阱
非线性和满阱由同一批平场数据(即光子响应曲线 PRC)算出,不额外采数据。
非线性定义中,每一步都要先修正"有效曝光时间"这个自变量,三项改正按影响从小到大:
- 光源强度漂移:$T_i = (U_i/\bar{U})\, t_i$,用监视二极管把每次曝光的实际光强归一到全序列均值;
- 快门时延:机械快门开合约 10 ms,短曝光时积分球总输出越来越非线性。做法是取 10 帧 30 ms(快门允许的最短曝光)的平场,把其平均信号 $F_0$ 与有效曝光时间 $T_0$ 从其他曝光中减掉;
- 采样率限制:皮安计最小采样间隔约 14 ms,所以 PRC 的低端无法用短曝光准确恢复(作者明确申报)。
然后逐点算相对偏差,区间内的非线性取 $|NL_i|$ 的最大值:
$$NL_i = \frac{F_i - F_0}{f\,(T_i - T_0)} - 1 $$
取最大值而不是 RMS 是关键:非线性的物理形态是单调往一个方向偏,RMS 会把端点处的系统性偏差平均掉、给出虚假的乐观值;而在误差预算里非线性是系统项,不会因多次测量下降,必须用上界记账。
满阱容量的定义是本文最值得照抄的一条:
$$\text{满阱容量} = NL\ \text{首次低于} -2\%\ \text{处的电子数} $$
不用饱和点的理由很直白:饱和之前响应已经相当非线性,那部分动态范围在科学上不可用。这样定义的好处是把"满阱"从一个器件物理参数变成了数据质量参数——满阱以内线性度是有保证的。
结果是:$14\text{–}56\ \mathrm{ke^-}$ 内非线性 $< 1\%$(0.03–0.96%),$7\text{–}70\ \mathrm{ke^-}$ 内 $< 1.4\%$(要求 $\le 2\%$);满阱 80–93 $\mathrm{ke^-}$。注意别混引两种口径:Table 2 里"要求 $\ge 70\ \mathrm{ke^-}$“是要求,摘要里”$\ge 80\ \mathrm{ke^-}$"是结果。
作者的两条自我限定:① 最终结果必须用飞行电子学复测;② 低信号端非线性更显著,达到亚百分比精度需要更精心设计的测试,另文报道。
4.3 暗电流
暗电流由中性体材料、耗尽材料、以及 Si–SiO₂ 界面处的少数载流子热产生引起——第三条(界面态)是背照 CCD 的重要来源,也是钝化工艺能压暗电流的原因。
做法是暗帧积分 1000 s 以累积足够信号;$-85\ ^\circ\mathrm{C}$ 取 5 帧暗加 60 帧偏置,$-95$ 与 $-75\ ^\circ\mathrm{C}$ 各取 3 帧暗加 10 帧偏置;偏置帧在暗帧前后各取一次以监视偏置电平漂移;逐像元跨帧取中位数以去除宇宙线。
结果:$-85\ ^\circ\mathrm{C}$ 下 $0.37\text{–}18 \times 10^{-3}\ \mathrm{e^- \, pix^{-1} \, s^{-1}}$,全部满足 $\le 0.02$,但器件间跨度超过一个数量级。$-95$ 与 $-85\ ^\circ\mathrm{C}$ 的差别在均值与分布上都不大,而 $-75\ ^\circ\mathrm{C}$ 显著增大,最差子型已贴到要求线——由此论证 $-85\ ^\circ\mathrm{C}$ 是性能与工程成本的平衡点。
自己验算一下。 暗电流近似 $I \propto T^3 \exp(-E_g/kT)$。用 e2vVIS 在 188.15 K(取区间中点 $11.3\times10^{-3}$)与 198.15 K($5.75\times10^{-2}$)的数据,比值为 5.1,解得 $E_g \approx 0.85\ \mathrm{eV}$。这个数字接近硅的带隙,说明测到的确实是本征热产生而非缺陷泄漏——若是缺陷主导,活化能会明显偏小。用同一关系把 188 K 外推到 178 K,暗电流只降到约 0.27 倍,与实测的 0.37–1.3 与 1.0–1.6(均 $10^{-3}$)同量级。这就定量解释了为什么从 $-85\ ^\circ\mathrm{C}$ 再往下压、收益很小却要多付制冷代价。
论文另发现两个器件效应,需要与真实暗电流区分:暗电流抬升斑(部分 e2v 器件读出寄存器附近,中心像素可超过 $0.1\ \mathrm{e^- \, pix^{-1} \, s^{-1}}$,成因未定)与放大器辉光(峰值出现在通道边缘或极近处)。
4.4 读出噪声
做法有三个细节都要照搬:数据来源是 $-85\ ^\circ\mathrm{C}$ 的偏置帧序列(不是暗帧);逐像元算标准差,计算前去掉该像元时序中最高 4 个与最低 4 个值;通道读出噪声取该通道所有像元的中位数。
这样设计的原因分别是:去掉极值抑制宇宙线与偏置波动;取中位数避免少数异常像元拉偏;用偏置帧则测到的是读出噪声本身,不含暗电流散粒噪声。
结果全部 $\le 5\ \mathrm{e^-}$(要求),中位数约 $3.6\ \mathrm{e^-}$,各子型 3.0–4.9 $\mathrm{e^-}$。作者再次强调需用飞行电子学复测。
4.5 量子效率
QE 测量的难点不在仪器读数,而在**"绝对通量"这一环**。论文用传递标准的思路:监视二极管与标定二极管处在同一积分球照明下,两者电流的比值只取决于响应率与位置,一次测准即可反复使用。入射光子数为
$$N_{\mathrm{ph}} = \frac{I_{\mathrm{mon}} R_\lambda \lambda}{hc} \cdot \frac{I_{\mathrm{cal}}}{I_{\mathrm{mon,cal}}} \cdot \frac{A_{\mathrm{pix}}}{A_{\mathrm{phot}}} \cdot T_\lambda \, t_{\mathrm{exp}} $$
逐项看:$I_{\mathrm{mon}}$ 是成像时监视二极管的平均电流;$R_\lambda$(单位 $\mathrm{W\,A^{-1}}$)是标定二极管的响应率,作用是把电流换算成绝对光功率;$\lambda/(hc)$ 把功率换成每秒光子数;$I_{\mathrm{cal}}/I_{\mathrm{mon,cal}}$ 完成从监视电流到绝对通量的传递;$A_{\mathrm{pix}}/A_{\mathrm{phot}}$ 把 $1\ \mathrm{cm^2}$ 上的通量缩到单个像元面积($100\ \mu\mathrm{m^2}$);$T_\lambda$ 是杜瓦窗口透过率;$t_{\mathrm{exp}}$ 是曝光时间。
两个容易忽略的点。 一是几何:标定二极管面积 $1\ \mathrm{cm^2}$,而 $N_e$ 取中心 $1000 \times 1000$ 像元 $= 1\ \mathrm{cm^2}$——两端取样面积相同,面积比项才成立,这正是必须用激光位移传感器先测出焦面中心位置、移走杜瓦后把标定二极管放到同一位置的原因。二是窗口透过率 $T_\lambda$:标定二极管测的是窗口之前的光,CCD 收到的是穿过窗口之后的光,不乘 $T_\lambda$ 会把窗口损失算进 QE、导致系统性低估;而它随波长变化(石英或蓝宝石在紫外端明显下降),所以必须写成 $T_\lambda$ 而不是常数。
结果(Table 2 平均 QE):250–400 nm 44NUV 达 74.1–75.7%(最高),e2vVIS 仅 34.8–39.3%;550–700 nm e2vDD 达 94.3–96.1%;800–900 nm e2vDD 77.4–79.5%;900–1000 nm e2vDD 仍有 38.6–40.9%,而 e2vVIS 只剩 19.1–20.4%。焦面为每个波段选择了 QE 最高的子型。
顺带纠正一个常见误解:“$\alpha \approx 1$"不等于"QE 不受波长影响”。$\alpha$ 是内禀的电荷产生效率,QE 还要乘涂层透过率与电荷收集效率。所以 250 nm 处 $\alpha$ 可以高达 1.28,而 QE 仍只有 74%——那里涂层与收集的损失很大。论文原话:“QE 由涂层透过率、量子产额与电荷收集效率共同决定。”
4.6 PRNU:三类图案与误差申报
做法是每个 CCD 在两个波长各取 25 帧 $45\ \mathrm{ke^-}$ 平场;因为积分球照明在 CCD 上只有缓慢变化,把归一化平场除以它自身的四阶多项式拟合以去掉照明不均;再用 $5\sigma$ 裁剪,取标准差作为 PRNU。
为什么用 $45\ \mathrm{ke^-}$(论文未明说理由,逻辑如下):光子噪声 $1/\sqrt{45000} \approx 0.47\%$,读出噪声贡献约 $0.009\%$,让被测量的 PRNU(约 1%)成为主导项;同时它落在"非线性 $\le 1\%$"的 $14\text{–}56\ \mathrm{ke^-}$ 区间内,远离饱和与亮-胖影响的端点。往下走光子噪声会反超 PRNU,往上走非线性增大,两侧都有代价。
最值得学的态度是误差申报:论文明确写出,四阶多项式在扣照明不均的同时也把 CCD 自身大尺度 PRNU 部分扣掉了,导致 PRNU 被低估 0.05–0.3%。
三类图案及其成因:
- 砖墙图案(CCD290-99,300 与 470 nm 显著,$>400$ nm 迅速消失):激光退火用脉冲激光逐块扫描背照表面,使表面薄层瞬时熔化再结晶以激活掺杂、修复离子注入损伤并钝化——它会留下周期性的非均匀印记,而 NUV 光子恰好只在这最表层几纳米内被吸收。长波光子穿透该层,图案随之消失。GM1201A 采用原子层沉积(ALD)工艺,逐层自限制生长,厚度精确、覆盖均匀,不产生砖墙图案——这是国产器件在该项上的实际优势。
- 同心环(e2vDD 850 nm):既不是 fringing(出现在更长波长)也不是更细的"年轮",论文判断为长波光子在前表面反射所致。
- 边缘窄特征:杜瓦内部杂散光;所有面板共有的中心偏亮:积分球照明不均。
典型 PRNU $< 1\%$。对测光的含义(§4 明确写出):常规平场无法完全改正这些波长相关 PRNU,会带来不可忽略的光度误差;作者引用自家 Luo 等(2024),并计划在飞行件上取更多波长的平场以引入半物理模型重建砖墙图案。更长波长下的 fringing 将在未来论文表征。
4.7 CTE
CTE 用扩展像元边缘响应(EPER)法测量。核心物理假设是:陷阱若持续遇到大小相近的密集电荷包,它释放一个电子后会立刻再俘获一个(净效果是一直饱和);一旦电荷包进入 overscan 区(那里本征光电子为零),释放的电子不再有替代来源,俘获概率随之下降。所以只要 overscan 足够长,成像区被俘获的电子最终全部释放到 overscan:
$$\mathrm{CTE} = 1 - \frac{N_{\mathrm{trap}}}{N_{\mathrm{flat}}\, n_{\mathrm{trans}}} = 1 - \frac{N_{\mathrm{trail}}}{N_{\mathrm{flat}}\, n_{\mathrm{trans}}} $$
其中 $n_{\mathrm{trans}}$ 并行 4616、串行 1152。操作上:每条 overscan 行(列)先离群裁剪再平均,成为拖尾剖面的一个像素;剖面末 50 像素的均值作为偏置;前 20 像素含 >90% 的俘获电子,取它们求和即得 $N_{\mathrm{trail}}$。
结论必须绑定信号水平。 陷阱在像元内均匀分布,因此陷阱密度恒定,$N_{\mathrm{trap}}$ 正比于电荷包占据的体积,而该体积随 $N_{\mathrm{flat}}^{\,\beta}$ 标度($\beta \approx 0.5$),故
$$1 - \mathrm{CTE} \propto N_{\mathrm{flat}}^{\,\beta - 1} = N_{\mathrm{flat}}^{-0.5} $$
信号越低,被俘获比例越大,CTE 越差。 结果印证:$10\ \mathrm{ke^-}$ 下串行与并行均 $> 99.9994\%$(并行最好达 99.99999%);$1\ \mathrm{ke^-}$ 下 e2vVIS 串行最低降到 99.99257%。串行拖尾比并行严重(LSSTCam 也一样)。
为什么报两个信号水平:一是覆盖实际使用的信号范围,两点还能验证 $\beta \approx 0.5$ 的幂律;二是 $1\ \mathrm{ke^-}$ 是辐照退化的灵敏探针——CTE 是空间 CCD 辐照退化最敏感的指标,而陷阱对低信号的影响远大于高信号(信号低十倍,$1-\mathrm{CTE}$ 约增大 $\sqrt{10} \approx 3$ 倍),用 $10\ \mathrm{ke^-}$ 作对照可把"陷阱增加"与"工作点变化"分开。三是与要求口径对应:Table 2 中 CTE 要求写的是"$10\ \mathrm{ke^-}$ 串行/并行 $\ge 99.999\%$",而 $1\ \mathrm{ke^-}$ 那一行标注为 “Information only”(仅供参考)。所以不要把 $1\ \mathrm{ke^-}$ 的结果当成"未达标"。
温度通过指数支配一切。 陷阱俘获电子后靠热激发释放,特征时间常数为
$$\tau_e \propto \frac{1}{\sigma_n v_{\mathrm{th}} N_c} \exp\!\left(\frac{\Delta E}{kT}\right) $$
能级深度 $\Delta E$ 决定量级(在指数上,主导),温度通过同一个指数放大它。取 $\Delta E = 0.3\ \mathrm{eV}$(典型辐照诱生陷阱),从 300 K 降到 188 K,$\exp(\Delta E/kT)$ 从约 $1.1\times10^5$ 增到约 $1.1\times10^8$——增大近三个数量级。而读出时标(约 7.7 ms 一行)是固定的,于是判据是:$\tau_e \ll$ 读出时标则无害,$\tau_e \gtrsim$ 读出时标则产生拖尾与电荷损失。这就是低温下 CTE 变差的原因,也解释了论文为什么把辐照退化留到后续(“To be continued”)。
五、总结与讨论:Table 2 与后续计划
5.1 亮-胖效应
本文用相关矩阵元素之和代替方差来消除亮-胖引起的 PTC 非线性。LSSTCam 的做法更进一步:用平场相关矩阵拟合亮-胖模型,同时求出转换增益与相互作用系数 $a_{ij}$(Astier 等 2019)。作者引用 Broughton 等(2024)指出,该模型的改正可恢复"与流量无关的 PSF"——对弱透镜分析是必需项,对测光精度也有二阶影响。但精确测量相关矩阵的各个元素需要多个流量水平下的数百张平场,并配以细致的误差控制;本文只给出三片 CCD 在若干波长下的初步结果。
三条结论:e2vDD 的 $a_{ij}$ 与 LSSTCam 深耗尽器件同量级;标准硅器件的 $a_{ij}$ 更低,因为其硅衬底更薄,新生成的光电子在被势阱收集前与像元内已有电荷发生相互作用的时间更短;$a_{ij}$ 随波长略微下降,因为长波光子在硅中吸收更深、留给相互作用的时间更短。
注意 Fig.17 的横轴是 $\sqrt{i^2+j^2}$(离中心像元的距离),而 $a_{ij}$ 随该距离增大而下降,不是正比关系。
5.2 Table 2:可直接借用的验收表范式
| 指标 | 条件 | 要求 | 结果(子型区间) |
|---|---|---|---|
| 非线性 | 14–56 / 7–70 ke⁻ | ≤1% / ≤2% | 0.03–0.96% / 0.08–1.37% |
| 满阱 | — | ≥70 ke⁻ | 80–93 ke⁻(受测试电子学限制) |
| 暗电流 | −85 ℃ | ≤0.02 e⁻/pix/s | 0.37–18 ×10⁻³(跨一个数量级) |
| 读出噪声 | — | ≤5 e⁻ | 中位数 ≈3.6 e⁻;各子型 3.0–4.9 |
| PRNU | 300/470/625/850 nm | ≤3% | 典型 <1%;e2vNUV 300 nm 达 2.94% |
| CTE | 10 ke⁻,串行/并行 | ≥99.999% | >99.9994%(并行最好 99.99999%) |
三条设计要点:每一项都写明条件与单位;要求只给必要的那一档、结果给区间;用脚注标明"仅供参考"的项、未装机的子型、以及每个子型只在自己的适用波段内需满足要求。
5.3 其他效应
串扰、像元内响应非均匀性、电荷扩散、辐照退化均推到后续论文。PRNU 方面,VIS/NIR 的年轮图案与深耗尽器件的 fringing 已有其他研究建立模型,本文计划在飞行件上取更多波长平场以纳入。
六、对 80 cm 红外望远镜 J/K 相机的移植要点
可直接照搬
- PTC 增益标定:成对平场差分方差 + 分块相关矩阵 + 扣读出噪声方差 + $< 0.8$ 倍满阱拟合;均值取中位数;裁掉放大器附近行列。注意红外阵列的邻域选择需按 IPC 与电荷扩散的实际范围重新实验确定,不能照抄三或五。
- 非线性判据与满阱口径:取区间内 $|NL|$ 最大值;满阱 = $NL$ 首次低于 $-2\%$ 处。对 K 波段尤其重要,因为 HgCdTe 的满阱与非线性是亮源测光饱和的根源。
- 曝光时间三项改正:光强漂移归一化、快门时延相减、采样率对低端不可恢复的申报。我们单帧只有 20 s,快门时延的相对影响比论文大 10–50 倍,这一项必须做。
- 暗电流三点温度法:用于论证工作温度与制冷余量。
- 读出噪声统计口径:去 4 高 4 低 + 通道中位数。
- PRNU 口径与误差申报:多项式扣照明不均 + $5\sigma$ 裁剪 + 标准差,并申报多项式导致的低估量。
- Table 2 验收表范式,以及用区间而非标准差别呈现器件间散布。
不可照搬:CTE/EPER 与串并行拖尾(混合式红外阵列没有逐行逐列电荷转移,也就没有 CTE 概念);亮-胖相关矩阵(红外主导耦合是 IPC 与电荷扩散);砖墙图案与 NUV 量子产额结论;$-85\ ^\circ\mathrm{C}$、满阱 $80\ \mathrm{ke^-}$、读出噪声 $3.6\ \mathrm{e^-}$ 等绝对数值。
必须自补:fringing、persistence、IPC,以及暗电流与噪声的偏压/温度二维扫描——前四项论文都推到了未来论文,而 fringing 与 persistence 正是我们 K 波段最痛的项。
我们额外要做的一项(论文没有的):J 与 K 必须各自独立标定增益波长、各自重做参考探测器的绝对定标。论文的参考二极管响应率与窗口透过率都是可见波段数据,我们的蓝宝石或石英窗口在 $1.25\ \mu\mathrm{m}$ 与 $2.2\ \mu\mathrm{m}$ 的透过率必须实测、不能由可见波段外推。
两个值得警惕的反向效应
- 降温压暗电流,但会让余像变长。 persistence 与 CTE 是同一套物理——缺陷能级俘获电荷后缓慢释放,时标由 $\tau_e \propto \exp(\Delta E/kT)$ 控制。我们为压暗电流要降温,余像反而更长。这两项的取舍必须实测,不能只引硅 CCD 的经验。
- 背景限 ≠ 暗电流可忽略。 J 相机暗电流约 $14\ \mathrm{e^- \, s^{-1} \, pix^{-1}}$、读出噪声约 $11\ \mathrm{e^-}$,20 s 单帧下暗电流累积约 $280\ \mathrm{e^-}$,其散粒噪声约 $16.7\ \mathrm{e^-}$,已超过读出噪声约 1.5 倍。所以必须分别测量读出噪声与暗电流,才能判断降噪的着力点在哪里。
七、几处待核实与论文未给出的地方
- 视场口径不一致:若按 $0.15''$/像元理解,31 片 9k×9k 覆盖约 $4.57\ \mathrm{deg^2}$,是论文所述 $1.1^\circ \times 1.2^\circ$($1.32\ \mathrm{deg^2}$)的 3.5 倍。更可能是"可用成像视场"与"焦面外接尺寸"两种定义混用。可查 CSST Collaboration: Gong 等(2026)。同时可确认:2 m 口径在 0.6 µm 的衍射极限为 $0.0755''$,所以 $0.15''$ 是采样/像质口径而非衍射极限。
- 串行/并行与像元速率口径:150 kHz 与全帧 39 s 自洽(1152 像元/150 kHz ≈ 7.7 ms 一行 × 4616 行 ≈ 35.5 s),但行转移时间 200 µs 单独就是 0.92 s,两者不是同一件事,引用时须分别标注。
- 工作温度的来源:是"先定温度再推暗电流要求"还是"先有科学要求再反推温度",论文未给出。从逻辑与 “operating temperature” 的措辞看更像后者(系统级交换)。
- 原文未给出:坏点与缺陷统计、串扰数值、像元内响应、电荷扩散、persistence、fringing、Fe55 测试、辐照后数据、低信号端非线性的定量结果。
八、小结
这篇论文的价值不在于产出天体物理结果,而在于把"探测器表征"写成了一份完整、可复现、可照抄的方法模板:只用偏置帧、平场、暗帧三类最普通的数据,就得到转换增益、量子产额、非线性、满阱、暗电流、读出噪声、QE、PRNU、CTE 九项结果,每一项都写明装置、判据、改正项与系统误差。
它的缺口同样明确:硅 CCD、红外混合式阵列的物理不可迁移,而 fringing、persistence、IPC 被作者明确推到未来论文。但正因为缺口被写清楚,它作为模板反而更好用——我们知道哪里能抄、哪里必须自己补。