电荷噪声

电荷噪声是量子点实验里最"会变形"的噪声：它既能移动电荷跃迁线，也能通过电荷混合、交换作用或自旋轨道耦合转化为比特相位噪声。

物理图像与定义

电荷噪声（charge noise）指器件中随机涨落的电势 $\delta V(t)$ 作用在量子点上，使点内电化学势、双点失谐（detuning）以及由此派生的一切比特参数发生随时间的漂移。它不是某一种具体的微观机制，而是”所有把随机电场送到量子点位置上的过程”的统称。

工程上更方便的做法是把它折算成能量涨落 $\Delta \epsilon$：栅极电压涨落 $\Delta V$ 经杠杆臂（lever arm）$\alpha$ 转换为点内能量涨落

$$\Delta \epsilon =\alpha \Delta V.$$

之所以要折算，是因为不同器件的隧穿率、$\alpha$ 各不相同，只有能量涨落才是可以横向比较的”器件噪声性能”指标。掺杂 GaAs 门控量子点上实测的典型值是数 $\mu \mathrm{eV}$，工艺优化后可压到亚 $\mu \mathrm{eV}$。

电荷噪声的关键特征是低频主导：其涨落的特征时间尺度在毫秒量级，功率谱密度在实验可及的频段近似正比于 $1/f$。这决定了它的两个实验后果——单次测量之间的稳定图会跳变、工作点会漂移；以及自由感应衰减（Ramsey）比自旋回波衰减得快得多。

微观起源

学界至今没有完全确定 $1/f$ 噪声的起源，但普遍认为其主要来源是随机电荷涨落：某个局域能级俘获或释放一个载流子，产生一次电势阶跃；大量具有不同特征时间的这类过程叠加，谱形自然趋近 $1/f$。在门控量子点中，被反复确认的具体通道有：

• 肖特基栅极到二维电子气的漏电流。掺杂 GaAs 器件中，栅极电压涨落驱动电子从肖特基电极经掺杂层漏向二维电子气，被认为是纳米器件噪声的主要来源之一。相应的抑制思路也很直接：移除栅极下方的二维电子气（浅刻蚀），或干脆移除掺杂层（非掺杂异质结 + 顶栅积累）。
• 掺杂层中电离施主的电荷涨落。这是GaAs/AlGaAs 调制掺杂结构的固有代价，也是转向非掺杂结构的主要动机。
• 半导体–介质界面的界面缺陷。氧化层中的悬挂键在量子尺度上可视为一个二能级系统（two-level system, TLS）。这些 TLS 既向耦合系统贡献电荷噪声、限制比特相干，又能与微波谐振腔交换能量、压低品质因子。同一批缺陷还通过”表面隧穿–俘获”过程屏蔽栅压，造成阈值电压漂移与工作电流衰减，使工作点校准变得困难。
• 栅氧厚度与工艺。过厚的氧化铝会显著增加电荷噪声，因此栅氧厚度需要在绝缘可靠性与噪声之间折中——Si/SiGe 重叠栅器件中通常控制在 $7\mathrm{nm}$ 左右。

材料平台的差异主要体现在界面质量上：Si-MOS 的电子被限制在 Si–SiO${}_2$ 界面，界面质量低于外延结构，均匀性有限、电荷噪声更大；Si/SiGe 与应变锗的外延界面则把载流子推离表面缺陷，噪声环境相对干净。

噪声谱：1/f 模型

把栅极 $G_i$ 上的电压噪声写成自相关函数的傅里叶变换

$$S_{G_i}(f)={\int }_{-\infty }^{+\infty }dt{e}^{-2\pi ift}⟨\delta V_{G_i}(t)\delta V_{G_i}(0)⟩$$

实验上常用的两种参数化形式为

$$S_q(\omega )=\frac{2\pi A_q}{|\omega |}S(f)=S_{\mathrm{hf}}{\left(\frac{f_0}{f}\right)}^{\alpha }$$

前者中 $A_q$ 是失谐参数 $q$ 的噪声在 $1\mathrm{Hz}$ 处的功率，后者把谱指数 $\alpha$ 留作可拟合量、$f_0$ 为参考频率。理想 $1/f$ 对应 $\alpha =1$；实测值通常略小于 1（例如平面锗空穴器件上拟合得到 $\alpha \approx 0.907$），仍属于典型的 $1/f$ 型电荷噪声。

$\alpha$ 不能直接测量，但可以从回波衰减的形状反推。Hahn 回波幅度随等待时间的衰减写成

$$\exp \left[-{\left(\frac{{\tau }_{\mathrm{idle}}}{T_2^{\mathrm{Echo}}}\right)}^{\beta }\right]$$

拟合出的指数 $\beta$ 与谱指数满足 $\alpha =\beta -1$。这条关系把”波形长什么样”和”噪声谱多陡”直接联系起来，是仅凭比特实验就能读出噪声谱形的实用技巧。

噪声如何进入比特：三类耦合通道

电荷噪声本身只涨落电势。它能否伤害比特，取决于比特频率对电势的导数有多大。

电荷比特：一阶失谐耦合

电荷量子比特的最小模型是

$$H=\frac{\epsilon }{2}{\sigma }_z+t_c{\sigma }_x\hslash {\omega }_{01}=\sqrt{{\epsilon }^2+4t_c^2}$$

其中 $\epsilon$ 是双点失谐、$t_c$ 是隧穿耦合。比特频率对失谐的敏感度

$$\frac{\partial {\omega }_{01}}{\partial \epsilon }=\frac{1}{\hslash }\frac{\epsilon }{\sqrt{{\epsilon }^2+4t_c^2}}$$

在 $\epsilon =0$ 处为零、在 $|\epsilon |\gg t_c$ 处饱和到 $1/\hslash$。这就是电荷比特”在反交叉点最抗噪、远离对称点退相干最快”的定量来源，也解释了为什么电荷比特虽有最大的电偶极矩、最容易做到电荷–光子强耦合，相干时间却最短。

交换类比特：二阶展开与协方差

对共振交换比特这类由两个失谐参数 $(\epsilon ,\Delta )$ 共同控制的比特，需要把哈密顿量展开到二阶。设电荷噪声使两个失谐各自偏移 $\delta q$（$q=\epsilon ,\Delta$），均值为零、服从高斯分布，则

$$\delta {\omega }_x=\frac{\partial {\omega }_{\mathrm{RX}}}{\partial \epsilon }\delta \epsilon +\frac{\partial {\omega }_{\mathrm{RX}}}{\partial \Delta }\delta \Delta +\frac{1}{2}\frac{{\partial }^2{\omega }_{\mathrm{RX}}}{\partial {\epsilon }^2}\delta {\epsilon }^2+\frac{1}{2}\frac{{\partial }^2{\omega }_{\mathrm{RX}}}{\partial {\Delta }^2}\delta {\Delta }^2+\frac{{\partial }^2{\omega }_{\mathrm{RX}}}{\partial \epsilon \partial \Delta }\delta \epsilon \delta \Delta +\cdots$$

假设所有电荷涨落具有 $1/f$ 谱，只保留导致退相位的纵向项，退相位率的主项为

$${\gamma }_{\varphi }=\frac{1}{2\pi \hslash }{\left[\sum _{q=\epsilon ,\Delta }\left(\frac{1}{4}{\left(\frac{\partial {\omega }_{\mathrm{RX}}}{\partial q}\right)}^2\mathrm{var}(q)+\frac{1}{16}{\left(\frac{{\partial }^2{\omega }_{\mathrm{RX}}}{\partial q^2}\right)}^2\mathrm{var}(q{)}^2\right)+\cdots \right]}^{1/2}$$

完整表达式还含 $\epsilon$ 与 $\Delta$ 的交叉项，其权重由相关系数 $c=\mathrm{cov}(\epsilon ,\Delta )/({\sigma }_{\epsilon }{\sigma }_{\Delta })$ 决定。这个形式的实用价值在于：当一阶导数被工作点设计压到零时，二阶项和交叉项就成为剩余退相干的主导，展宽随失谐的依赖曲线因此不是简单的 V 形，而是可以用上式整体拟合——拟合与实验数据吻合，即可判定电荷噪声是主导噪声源。

作为对照，同一体系中电子–声子相互作用引起的弛豫可用形变势哈密顿量与费米黄金定则估算：硅中声子贡献相对电荷噪声可忽略，而 GaAs 中电声相互作用占比相当大，不能忽略。

自旋比特：纵向自旋–电场敏感度

自旋本身不带电，电荷噪声必须”借道”才能影响它。借道方式有二：

1. 自旋轨道耦合 / $g$ 张量的电压依赖。把空穴比特写成有效自旋 $1/2$ 系统

$$H_s=\mathbf{S}\cdot {\mathbf{\omega }}_L({\mathbf{V}}_G){\mathbf{\omega }}_L({\mathbf{V}}_G)=\frac{{\mu }_B}{\hslash }\mathsf{g}({\mathbf{V}}_G)\cdot \mathbf{B}$$

$g$ 张量显含栅压，于是拉莫尔频率 $f_L=\frac{{\mu }_B}{h}\left|\mathsf{g}({\mathbf{V}}_G^0)\cdot \mathbf{B}\right|$ 也随栅压变化。定义纵向自旋–电场敏感度（longitudinal spin–electric susceptibility, LSES）

$$D_{G_i}=\frac{\partial f_L}{\partial V_{G_i}^0}$$

它的大小直接反映比特对电荷噪声的敏感程度，实验上就是”比特频率随柱塞栅压的斜率”。

2. 微磁体纵向梯度场。微磁体的杂散场与电驱动位移结合给出横向梯度 $b_{tr}=({\vec{e}}_{MW}\cdot \nabla )B_{MM}^{\perp }$，驱动EDSR；而其纵向分量与涨落电场结合给出 $b_{long}=({\vec{e}}_{noise}\cdot \nabla )B_{MM}^{//}$，直接把电荷噪声转成频率噪声。二者的比值定义了品质因子

$$Q=\frac{\mathrm{d}{B}_{tr}^{tot}/\mathrm{d}y}{\sqrt{{\left(\mathrm{d}{B}_{long}/\mathrm{d}x\right)}^2+{\left(\mathrm{d}{B}_{long}/\mathrm{d}y\right)}^2}}$$

最大化 $Q$ 的面内磁场角并不沿磁体轴向，而是偏离一定角度。

在平面锗中，有效 $g$ 因子的各向异性源于重空穴–轻空穴（HH–LH）混合，而混合程度又由量子点的电势分布决定，因此 $g$ 与 LSES 都强烈依赖磁场取向。面内旋转磁场时

$$g(\theta )=\sqrt{g_a^2{\cos }^2(\theta -{\theta }_0)+g_b^2{\sin }^2(\theta -{\theta }_0)}$$

对应的自旋–电场敏感度解析近似为

$$\frac{\partial g}{\partial V}=\frac{1}{g}\left(g_x\frac{\partial g_x}{\partial V}{\cos }^2\theta +g_y\frac{\partial g_y}{\partial V}{\sin }^2\theta \right).$$

退相干的滤波函数理论

把上述”敏感度 × 噪声谱”整合成相干时间，需要滤波函数（filter function）形式。栅压噪声在自由演化中累积随机相位

$$\delta \varphi (t)=2\pi {\int }_{-\infty }^{+\infty }d{t}^{\prime }\sum _i{D}_{G_i}\delta V_{G_i}(t^{\prime }){\eta }_t(t^{\prime })$$

其中 ${\eta }_t(t^{\prime })$ 描述脉冲序列的符号翻转。对含 $N_{\pi }$ 个 $\pi$ 脉冲的 CPMG 序列

$${\eta }_t^{\mathrm{CPMG}}(t^{\prime })=\sum _{k=0}^{N_{\pi }}(-1{)}^k\Theta (t_{k+1}-t^{\prime })\Theta (t^{\prime }-t_k)$$

$\Theta$ 为 Heaviside 阶跃函数，Hahn 回波对应 $N_{\pi }=1$。旋转坐标系中密度矩阵非对角元的衰减为

$$⟨{\tilde{\rho }}_{01}(t)⟩={\tilde{\rho }}_{01}(0)e^{-\frac{1}{2}⟨\delta {\varphi }^2(t)⟩}⟨\delta {\varphi }^2(t)⟩=4{\pi }^2{\int }_{-\infty }^{+\infty }df\sum _i{D}_{G_i}^2{S}_{G_i}(f){\left|{\tilde{\eta }}_t(f)\right|}^2$$

即退相干速率 = 噪声谱与滤波函数的重叠积分。Ramsey 序列的滤波函数

$${\left|{\tilde{\eta }}_t^R(f)\right|}^2={\left(\frac{\sin (\pi ft)}{\pi f}\right)}^2$$

在 $f\to 0$ 处不衰减，因此对低频分量全盘接收；对 $1/f$ 谱做积分会在低频端发散，须用实验时长与带宽设定截止 $f_l$、$f_h$，得到

$$\frac{1}{T_2^{\ast }}\approx 2\pi \sqrt{\ln \left(\frac{f_h}{f_l}\right)f_0\sum _i{\left(\frac{\partial f_Q}{\partial V_{P_i}}\right)}^2{S}_{P_i}^{\mathrm{lf}}}.$$

$\pi$ 脉冲把滤波函数的低频响应压掉，回波时间因而由较高频段的噪声幅度决定；取 $\alpha =0.907$ 时

$$\frac{1}{T_2^{\mathrm{Echo}}}\approx 7.6f_0^{1/\alpha }{\left(\sum _i{\left(\frac{\partial f_Q}{\partial V_{P_i}}\right)}^2{S}_{P_i}^{\mathrm{hf}}\right)}^{0.695}.$$

两式的对比解释了实验中反复出现的现象：$T_2^{\ast }$ 常常对工作点变化不敏感（被幅度很大的低频分量钉住），而 $T_2^{\mathrm{Echo}}$ 会随敏感度 $D_{G_i}$ 明显起伏。反过来，用不同 $N_{\pi }$ 的 CPMG 测一组 $T_2^{\mathrm{CPMG}}$，再配合长时间 Ramsey 频率监测，就能反演出跨越多个数量级的噪声谱——这是动力学解耦序列作为”噪声谱仪”的用法。

最优工作点

既然伤害正比于导数，压制电荷噪声的通用策略就是把工作点放在导数为零的位置，即最优工作点（sweet spot）。三个典型例子：

• 电荷比特的对称点 $\epsilon =0$：$\partial {\omega }_{01}/\partial \epsilon =0$，一阶失谐噪声被完全抑制，代价是失去了用失谐调频的手段。
• 交换门的对称操作点。交换耦合与失谐 $ϵ=\alpha (V_{P2}-V_{P1})$ 的关系为

$$J=\frac{2tU}{U^2-{ϵ}^2}$$

$t$ 为点间隧穿耦合、$U$ 为充电能。$J$ 在 $ϵ=0$ 处取极小值、对 $ϵ$ 一阶不敏感，因此实验上不再用失谐脉冲、而是用快速电压脉冲直接抬降点间势垒来开关 $J$，可有效抑制电荷噪声对 $J$ 的干扰。需要注意的是，正反脉冲会使两点电化学势非对称偏移，实测对称点常偏离失谐零点，需用 Hahn 回波型序列单独标定（引入虚拟栅极可缓解）。

• 空穴自旋比特的磁场取向甜点。由于 LSES 随磁场角呈抛物线型变化，选择使 $|D_{G_i}|$ 最小的磁场方向即可显著抑制退相干。理论进一步指出，LSES 为零的方向在描述磁场取向的单位球面上不是孤立点，而是连成连续轨迹（“最优操控线”，sweet lines），且轨迹位置可由栅压调控——这为在多比特阵列中寻找共同最优磁场配置提供了可能。

代价是普遍存在的：“操控速度–相干时间”权衡。强自旋轨道耦合是一把双刃剑：它带来超快全电学操控，也放大自旋–电场敏感度。实验上确实看到 $f_{\mathrm{Rabi}}$ 与 $T_2^{\mathrm{Echo}}$ 随磁场角呈相反趋势，使品质因子 $Q_E=2T_2^{\mathrm{Echo}}{f}_{\mathrm{Rabi}}$ 不再有单调规律。

参数与量级

量	典型值	来源
能量涨落 $\Delta \epsilon$（传统掺杂 GaAs 单点）	$3.1\pm 1.1\mu \mathrm{eV}$
能量涨落 $\Delta \epsilon$（浅刻蚀 GaAs 单点）	$0.50\pm 0.23\mu \mathrm{eV}$
能量涨落 $\Delta \epsilon$（非掺杂 GaAs 单点）	$0.55$–$1.4\mu \mathrm{eV}$
杠杆臂 $\alpha$（用于 $\Delta \epsilon =\alpha \Delta V$）	$0.05$（浅刻蚀）/ $0.03meV/mV$（传统）
噪声谱指数 $\alpha$	$0.907$（由 $\beta =1.907\pm 0.23$ 反推）
电荷涨落特征时间尺度	毫秒量级（低频主导）
比特频率抖动（55 h Ramsey 监测）	近似高斯分布，$\sigma =0.794\mathrm{MHz}$
内禀失谐噪声幅度 $\delta f_{\mathrm{Larmor}}$	$\approx 1.66\mathrm{MHz}$（由 $T_2^{\ast }$ 估算），拟合取 $1.80\mathrm{MHz}$
平面锗空穴比特 $T_2^{\ast }$	$169.9\mathrm{ns}$，随面内磁场角近似不变
平面锗空穴比特 $T_2^{\mathrm{Echo}}$	Q1 $500\to 2200\mathrm{ns}$；Q2 $400\to 800\mathrm{ns}$
甜点处单比特门保真度增益	$+2\%$，最高 $99.82\%$
RX 比特退相干率 ${\gamma }_{\mathrm{RX}}/2\pi$	$16.9\mathrm{MHz}$
Si-MOS 微磁体器件 $T_2^{\ast }$ / $T_2^{\mathrm{echo}}$	$624\mathrm{ns}$ / $60.2\mu \mathrm{s}$	楚宁 2025
电流噪声积分频段	$1$–$9\mathrm{Hz}$（浅刻蚀）；$5$–$45\mathrm{Hz}$（非掺杂）
温度依赖	$250\mathrm{mK}\to 1\mathrm{K}$，峰顶电流涨落上升约 $20\%$

实验测量方法

输运法：库仑峰上的电流噪声

这是把电荷噪声折算成能量的经典方案，四个步骤：

1. 定标。先测库仑菱形提取杠杆臂 $\alpha$。
2. 测谱。把量子点调到库仑阻塞区（避开能量量子化的影响），源极加直流偏压（如 $200\mu \mathrm{V}$，用直流而非交流以排除激励信号干扰），漏极电流经前置放大器转成电压后送入频谱仪，得到电流噪声谱 $S_I(f)$。在库仑峰的峰底（阻塞区）、斜率最大处、峰顶三个位置分别取谱：峰底零电流处测到的实际上是隧穿几率涨落 $\Delta \Gamma$ 引起的噪声，可作为本底。
3. 积分。对低频段积分得到电流涨落

$$\Delta I=\sqrt{{\int }_{f_1}^{f_2}\left[S_I^2(f)-S_{BG}^2(f)\right]df}$$

$S_{BG}$ 为系统噪声本底。 4. 换算。$\Delta V=\Delta I/\left|dI/dV\right|$，再乘 $\alpha$ 得 $\Delta \epsilon$。

一个重要的自洽性判据：$\Delta I$ 沿库仑峰扫描应出现两个峰，与微分电导 $|dI/dV|$ 的两个峰一一对应，即电流噪声正比于微分电导——这说明噪声确实来自栅极电压涨落，而不是别的机制。峰顶附近 $|dI/dV|\to 0$ 会让 $\Delta V$ 发散，这几个点必须舍弃。

比特谱学法

• $T_2^{\ast }$ 与多脉冲解耦时间之比。$T_2^{\mathrm{Echo}}/T_2^{\ast }$ 越大，说明噪声越集中在低频；配合 CPMG 阶数扫描可重构噪声谱。
• 回波衰减指数。拟合 $\beta$ 得 $\alpha =\beta -1$，直接给出谱形。
• 长时间频率监测。固定等待时间做 Ramsey，连续数十小时跟踪共振频率漂移，得到极低频端（$10^{-5}\mathrm{Hz}$ 量级）的谱点。
• LSES 扫描。三段式（Empty–Load–Readout）脉冲中改变 Load 阶段的柱塞栅压，测比特频率随 $\delta V_P$ 的线性斜率，即 $\partial f_Q/\partial V_P$；不同栅极给出的 LSES 可正可负，符号反映波函数在该栅极作用下是被拉伸还是被压缩。
• 腔谱线宽。在cQED 构型中，用双色调制谱提取比特展宽随失谐的依赖关系，再用上文的二阶电荷噪声模型整体拟合。

需要强调：不同方法测到的噪声投影不完全相同。输运法测的是量子点自身能量涨落，比特谱学法测的是”噪声谱 × 该比特敏感度 × 该序列滤波函数”的加权积分。此外，比特模型往往只包含能施加快速脉冲的柱塞栅极，势垒栅极相关的噪声通道容易被遗漏——这正是数值模拟能复现 $T_2^{\ast }$ 却低估 $T_2^{\mathrm{Echo}}$ 各向异性的原因之一。

抑制手段

层次	手段	效果
材料	移除掺杂层（非掺杂异质结 + 顶栅积累）	GaAs 单点能量涨落降至 $0.55$–$1.4\mu \mathrm{eV}$
工艺	浅刻蚀移除栅极下方二维电子气，切断漏电流通道	相对传统结构降低约一个数量级
工艺	界面钝化处理、控制栅氧厚度	降低界面态密度与 TLS 参与度
工作点	电荷比特对称点、交换门对称操作点、磁场取向甜点	消去一阶敏感度
波形	势垒脉冲代替失谐脉冲控制 $J$	抑制电荷噪声对 $J$ 的干扰
序列	Hahn 回波 / CPMG 动力学解耦	滤除低频分量，$T_2$ 可提升数倍至数十倍
器件设计	优化微磁体几何、优化腔电场分布远离缺陷区	减小纵向梯度、降低缺陷耦合

与其他概念的关系

• 电荷噪声直接移动电化学势，因此电荷稳定图中的跃迁线会随时间漂移、库仑阻塞峰会跳变，这也是自动调控必须处理的稳定性问题。
• 它是电荷比特相干时间的主要限制因素；对单自旋比特和空穴自旋比特，则经 $g$ 因子、自旋轨道耦合或微磁体梯度间接作用。
• 通过交换作用 $J(ϵ)$，电荷噪声成为单重态–三重态比特、共振交换比特和杂化比特两比特门误差的来源。
• 界面缺陷既是电荷噪声的微观载体，也是高阻抗谐振腔内损耗的来源，二者需一起优化才能满足强耦合判据。
• 测量层面，QPC 电荷传感与射频反射测量既是探测电荷噪声的工具，其自身工作点也受电荷噪声漂移影响。
• 在翻转模式比特与自旋–光子耦合中，比特被有意置于电荷–自旋混合区以增大偶极矩，代价正是更强的电荷噪声敏感性。

参考文献

• 电荷噪声对门保真度的影响与对策：Burkard et al., RMP 95, 025003 (2023)、Hendrickx et al., Nat. Mater. 23, 920 (2024)。

完整文献库（含各篇站内全文页）见 参考文献库。