2. 量子测量中的噪声与漂移源分类
量子测量的精度极限并不只由量子力学本身决定。在真实实验平台上,被测量子态往往要在室温或低温环境中经历一段有限的相干时间,而测量仪器、样品架、线缆、屏蔽体乃至建筑本身都在持续向系统注入扰动。这些扰动可以粗略地分成两类:噪声表现为围绕均值的快速随机起伏,通常具有零均值、可用功率谱密度(PSD)描述;漂移则是慢变的、非零均值的系统性偏移,它不会因为多次平均而消失,反而会在长时间测量中主导误差预算。理解各类噪声与漂移的物理机制和频谱特征,是后续设计补偿策略的前提。
一、热噪声:白噪声的物理下限
热噪声(Johnson–Nyquist 噪声)来自导体中载流子的热运动。任何处于温度 T 的电阻 R,其两端都会出现电压涨落,单边功率谱密度为 SV = 4kBTR,其中 kB 为玻尔兹曼常数。在量子测量中,热噪声主要出现在三个位置:探测器的前端放大器、连接量子器件与室温电子学的同轴线、以及低温恒温器内部的衰减器与滤波器。
热噪声的关键特征是频谱平坦(白噪声),在远低于 kBT/h 的频率范围内与频率无关。以 300 K 为例,kBT/h ≈ 6.25 THz,因此在微波频段(几 GHz)热噪声几乎是纯白的。这意味着它无法通过降低测量带宽以外的办法消除,只能通过冷却、降低电阻或使用量子极限放大器来压制。例如,将前端从 300 K 冷却到 4 K,热噪声功率下降约 18.8 dB;若进一步到 20 mK,则下降约 42 dB。这也是超导量子比特普遍工作在稀释制冷机中的原因之一。
二、1/f 噪声:低频漂移的普遍来源
1/f 噪声(闪烁噪声)的功率谱密度近似为 S(f) ∝ 1/fα,其中 α 通常在 0.8~1.2 之间。它广泛存在于半导体器件、超导约瑟夫森结、量子点接触以及材料缺陷中。其物理机制并不唯一:在半导体中常归因于载流子被缺陷态随机俘获与释放,产生洛伦兹型涨落的叠加;在超导量子比特中,则与衬底、界面处的两能级系统(TLS)有关。
1/f 噪声对量子测量的影响体现在两个方面。第一,它直接导致量子比特频率的慢漂移,使 Ramsey 干涉条纹在长时间平均后失去对比度。第二,它使测量信号的基线缓慢移动,若不加处理,长时间积分会得到错误均值。判断 1/f 噪声是否主导,可以观察 Allan 方差曲线:当平均时间 τ 超过某一拐点后,Allan 偏差不再随 τ-1/2 下降,而是趋于平坦甚至上升,这正是 1/f 漂移的典型标志。
三、电磁干扰:外部耦合的窄带与宽带成分
电磁干扰(EMI)来自外部电磁场通过传导或辐射途径耦合进测量链路。常见来源包括:开关电源的开关频率及其谐波(几十 kHz 到几 MHz)、无线通信信号(MHz~GHz)、实验室中其他脉冲设备的瞬态辐射。EMI 的特征是窄带尖峰或周期性脉冲,在频谱上表现为离散谱线,而非平滑连续谱。
与热噪声和 1/f 噪声不同,EMI 往往具有明确的频率标签,因此可以通过屏蔽、滤波、接地和同步触发来抑制。例如,若发现频谱在 100 kHz 处有稳定尖峰,可先检查附近是否有开关电源;若尖峰频率随设备位置变化,则更可能是辐射耦合而非传导耦合。在量子测量中,EMI 的另一个危险是它可能被非线性器件整流,转化为低频漂移,从而伪装成 1/f 噪声。
四、振动与温度波动:机械与热学的慢变扰动
振动通过两种方式影响量子测量:一是直接改变器件与天线、腔体或磁体的相对位置,从而改变耦合强度;二是通过压电效应或摩擦生电在电缆中产生电荷。振动的频谱通常集中在几十 Hz 到几百 Hz,与建筑机械、泵组、人员走动有关。在超导量子比特中,振动还会引起磁通噪声,因为约瑟夫森结环路面积对机械形变敏感。
温度波动则通过热膨胀、超导临界电流变化、以及热噪声功率变化影响测量。温度漂移的时间尺度通常从秒到小时,属于典型的慢漂移。例如,稀释制冷机混合室温度若在 10 mK 附近有 ±0.5 mK 的缓慢起伏,就足以使量子比特频率发生可观测的偏移。温度波动与 1/f 噪声在低频段常常叠加,难以从单一频谱中区分,需要结合温度传感器同步记录来分离。
五、磁场漂移:对相位敏感的量子系统尤为致命
磁场漂移来自地磁场变化、附近磁性物体移动、以及磁屏蔽层本身的磁化弛豫。对于超导量子比特、原子磁力计、金刚石 NV 色心等对磁通敏感的体系,磁场漂移会直接改变能级分裂,表现为频率的慢漂移。其时间尺度可从秒到小时,频谱上同样呈现 1/f 特征,但物理来源与材料缺陷无关,而是外部环境。
区分磁场漂移与材料 1/f 噪声的一个实用方法是:在磁屏蔽筒内放置一个参考磁力计,同步记录磁场。若量子比特频率漂移与参考磁力计读数高度相关,则磁场漂移主导;若相关性弱,则更可能是材料或电荷噪声。磁场漂移的抑制通常依靠多层磁屏蔽(如坡莫合金加超导屏蔽)和主动反馈线圈。
六、从频谱到误差预算:如何判断谁在主导
面对一台具体的量子测量装置,第一步不是立刻上补偿算法,而是先做噪声表征。推荐的操作顺序如下:
- 采集时域数据:在测量系统处于稳定工作点时,连续记录探测器的输出电压或量子比特的读出信号,采样率至少覆盖关心的最高频率,记录时长应远大于最慢漂移的周期(例如 1 小时以上)。
- 计算功率谱密度:用 Welch 方法分段平均,得到 S(f)。观察低频段是否呈 1/f 斜率,中频段是否平坦,高频段是否出现离散尖峰。
- 计算 Allan 方差:对同一组时域数据做 Allan 偏差分析。若曲线在 τ 大于某值后偏离 τ-1/2,说明存在 1/f 或更慢的漂移。
- 同步记录环境量:同时记录温度、磁场、振动加速度。计算它们与测量信号的相关系数。若相关系数绝对值大于 0.5,则该环境量很可能是主导漂移源。
- 做受控实验:例如暂时关闭附近开关电源,观察 EMI 尖峰是否消失;或改变磁屏蔽状态,观察低频漂移是否变化。通过“关掉一个源、看谱线变化”的方式逐一排除。
七、各类扰动的特征对照
| 扰动类型 | 物理机制 | 频谱特征 | 典型时间尺度 | 主要抑制手段 |
|---|---|---|---|---|
| 热噪声 | 载流子热运动 | 白噪声,平坦 | ns~μs | 冷却、降低 R、量子极限放大 |
| 1/f 噪声 | 缺陷俘获/释放、TLS | 1/fα,α≈1 | μs~s | 动态解耦、反馈、材料优化 |
| 电磁干扰 | 外部场传导/辐射耦合 | 离散尖峰或脉冲 | ns~ms | 屏蔽、滤波、接地、同步 |
| 振动 | 机械位移改变耦合 | 几十~几百 Hz 窄带 | ms~s | 隔振平台、刚性固定 |
| 温度波动 | 热膨胀、临界电流变化 | 低频 1/f 型 | s~h | 恒温控制、热锚定 |
| 磁场漂移 | 地磁、磁性物体移动 | 低频 1/f 型 | s~h | 磁屏蔽、主动反馈线圈 |