4. 噪声源分析与建模

量子传感器的读出链路上，真正决定“能不能测到”的往往不是信号有多大，而是噪声有多小。一个 SQUID 磁强计、一个原子自旋进动探测器或一个超导纳米线单光子探测器，其输出端总是叠加着若干来源不同、频谱形状各异的噪声。工程上要做的第一件事，是把总噪声拆成可辨识的分量，再判断哪一项在目标频段内占主导，从而决定是改电路、降温度、换放大器，还是调整调制频率。

一、五类噪声的物理来源与谱特征

热噪声（Johnson–Nyquist 噪声）来自导体中载流子的随机热运动，任何有电阻的元件都会产生。其电压谱密度为 S_V = 4k_BT R，电流谱密度为 S_I = 4k_BT/R，在远低于 k_BT/h 的频率范围内是白噪声。室温下 1 kΩ 电阻的电压噪声密度约为 4 nV/√Hz，这是很多读出电路的本底参照。降低热噪声只有两条路：降低温度或降低等效电阻，但后者往往与带宽、阻抗匹配冲突。

散粒噪声（shot noise）源于电荷或光子到达事件的离散性，谱密度 S_I = 2qI（q 为电子电荷，I 为平均电流），同样是白噪声。在光电探测、隧道结和单电子器件中，散粒噪声常是主导项。若探测器量子效率为 η，光电流散粒噪声对应的等效光功率噪声为 NEP_shot = √(2hνP/η)，可见提高量子效率和增大光功率都能压低该项。

1/f 噪声（闪烁噪声）的谱密度近似为 S(f) = A/f^α，α 通常在 0.8～1.2 之间。它来自材料缺陷、界面陷阱、载流子俘获-释放等慢过程，在低频段迅速抬升。对直流或准直流工作的量子传感器（如磁强计、加速度计），1/f 噪声决定了低频灵敏度的“地板”。工程上常用斩波、调制-解调或相关双采样把工作点搬到 1/f 拐点之上。

放大器噪声包括输入级电压噪声 e_n、电流噪声 i_n 和等效输入噪声。它与源阻抗共同作用：源阻抗 Z_s 上的总输入噪声为 √(e_n² + (i_nZ_s)²)。存在一个最佳源阻抗使总噪声最小，这就是噪声匹配。低温量子读出常用 SiGe HBT、HEMT 或 SQUID 前置放大器，其 e_n 可低至 0.1 nV/√Hz 量级。

量子反作用噪声是量子极限读出特有的：测量本身会扰动被测系统。以光机或微波腔探测为例，位移测量精度受标准量子极限（SQL）约束，散粒噪声（测量噪声）与辐射压力噪声（反作用噪声）之和在 SQL 处取极小值。当探测功率增大时，测量噪声下降但反作用噪声上升，二者交叉点即为最优探测功率。这一项无法通过常规降噪消除，只能通过压缩态注入、背作用规避（back-action evasion）等量子技术突破。

二、噪声等效功率（NEP）与灵敏度指标

NEP 是衡量探测器噪声性能的核心指标，定义为使输出信噪比等于 1 时所需的输入信号功率，单位 W/√Hz。对光探测器，NEP = i_n / R，其中 i_n 为总电流噪声密度（A/√Hz），R 为响应度（A/W）。若同时考虑散粒噪声和热噪声，则 NEP_total = √(NEP_shot² + NEP_thermal² + NEP_1/f² + …)，各项按平方和叠加。

灵敏度则更贴近系统层面，通常指可分辨的最小信号变化，如磁强计的 pT/√Hz、加速度计的 μg/√Hz。灵敏度与 NEP 的关系取决于换能系数：若磁强计的输出电压噪声为 S_V（V/√Hz），磁场响应为 ∂V/∂B（V/T），则磁场灵敏度 S_B = S_V / (∂V/∂B)，单位 T/√Hz。工程上常把灵敏度换算为等效输入噪声，便于与信号幅度直接比较。

三、建模与计算示例

假设一个原子磁强计读出链路：光电流 I = 100 μA，量子效率 η = 0.8，波长 780 nm，探测器负载电阻 R = 50 kΩ，前置放大器 e_n = 2 nV/√Hz，i_n = 10 fA/√Hz，1/f 噪声在 1 Hz 处为 50 nV/√Hz，拐点频率 10 Hz。目标工作频段 100 Hz～1 kHz。

散粒噪声电流密度：
  i_shot = sqrt(2 * q * I) = sqrt(2 * 1.6e-19 * 100e-6)
         ≈ 5.66e-12 A/√Hz = 5.66 pA/√Hz

热噪声电流密度（R = 50 kΩ, T = 300 K）：
  i_thermal = sqrt(4 * k_B * T / R)
            = sqrt(4 * 1.38e-23 * 300 / 50e3)
            ≈ 1.82e-14 A/√Hz = 0.0182 pA/√Hz

放大器电流噪声贡献（源阻抗 50 kΩ）：
  i_amp = e_n / R = 2e-9 / 50e3 = 4e-14 A/√Hz = 0.04 pA/√Hz

总电流噪声（100 Hz，1/f 已可忽略）：
  i_total = sqrt(i_shot² + i_thermal² + i_amp²)
          ≈ sqrt((5.66e-12)² + (1.82e-14)² + (4e-14)²)
          ≈ 5.66 pA/√Hz

NEP = i_total / R_resp，若响应度 R_resp = 0.5 A/W：
  NEP ≈ 5.66e-12 / 0.5 ≈ 1.13e-11 W/√Hz = 11.3 pW/√Hz

计算表明，在 100 Hz 以上，散粒噪声比热噪声和放大器电流噪声高约两个数量级，是绝对主导项。此时降低 NEP 最有效的措施是提高光电流或量子效率，而不是更换更低噪声的放大器。若把工作频率降到 1 Hz，1/f 噪声折算到输入端约为 50 nV/√Hz，对应电流噪声 1 pA/√Hz，仍小于散粒噪声，但已接近同一量级，此时调制解调的必要性就凸显出来。

四、量子反作用噪声的建模要点

量子反作用噪声不能简单当作加性噪声处理，它与测量强度耦合。以腔光机位移探测为例，测量噪声谱密度 S_meas ∝ 1/P_c，反作用噪声 S_ba ∝ P_c，总噪声 S_total = S_meas + S_ba 在 S_meas = S_ba 时取最小值，即 SQL。建模时需要同时标定探测功率与噪声谱，找到总噪声最低点对应的功率，而不是一味增大激光功率。若系统已接近 SQL，进一步降低噪声只能借助压缩真空态注入，将某一正交分量的噪声压到 SQL 以下。

实际建模流程建议：先测量暗噪声本底，再逐项屏蔽或改变条件（如改变温度、改变光功率、改变调制频率）观察噪声谱变化，从而分离各分量；然后用平方和模型拟合总噪声曲线，提取 1/f 拐点、散粒噪声系数和放大器噪声参数；最后代入 NEP 和灵敏度公式，得到目标频段内的预期指标，并与实测信噪比对比验证模型是否自洽。