交叉电容矩阵

交叉电容矩阵是一张"谁影响谁"的局部地图：对角元是栅极对目标量子点的控制效率，非对角元量化它对其他量子点的串扰。

物理图像

门控定义的半导体量子点中，每根金属栅极的电场在空间上没有硬边界。一根柱塞栅（plunger gate）除了决定”自己的”量子点的电化学势，还会通过寄生电容（parasitic capacitance）拉动邻近量子点的电化学势，势垒栅则会同时抬高/压低左右两侧点间的势垒。在电荷稳定图上，这种栅极串扰（capacitive crosstalk）表现为倾斜的电荷隧穿线——理想情况下某点的隧穿线应只沿自身柱塞栅方向延伸，串扰使其对邻近栅压也产生依赖。

单点、双点时代，实验者凭经验手动”回拧”其他栅压即可补偿；但当量子点阵列扩展到几十根栅极时，串扰使控制空间变成一个高维、强耦合的电压空间，人工试探式调谐很快失效。交叉电容矩阵（cross-capacitance matrix，又称串扰矩阵）就是用来量化这种耦合的线性响应矩阵：在某个局部工作点附近，物理栅压的微小变化如何映射到各量子点的电化学势变化。它的逆矩阵直接给出虚拟电极的线性组合系数，矩阵条件数过大则预警控制轴接近线性相关。该概念最早由代尔夫特理工大学 Vandersypen 组在九点 Fermi–Hubbard 阵列等工作中系统化，2019 年普林斯顿大学 Petta 组又将其纳入自动校准流程，已成为自动调控与阵列扩展的标准接口。

理论模型

局部线性化与线性响应

对一个工作在弱阻塞区（$(N_1,N_2,\dots ,N_n)$ 固定）的阵列，把各量子点电化学势 $\mathbf{\mu }=({\mu }_1,\dots ,{\mu }_n{)}^{\mathrm{T}}$ 与各物理栅压 $\mathbf{V}=(V_1,\dots ,V_m{)}^{\mathrm{T}}$ 的局部依赖泰勒展开到一阶：

$$\Delta \mathbf{\mu }={\left(\frac{\partial \mathbf{\mu }}{\partial \mathbf{V}}\right)}_{{\mathbf{V}}_0}\Delta \mathbf{V}+\mathcal{O}(\Vert \Delta \mathbf{V}{\Vert }^2).$$

Jacobian 矩阵 $\partial \mathbf{\mu }/\partial \mathbf{V}$ 即交叉电容矩阵 $\mathbf{T}$ 的核心，其元素 $t_{ij}=\partial {\mu }_i/\partial V_j$ 表征”栅极 $j$ 改变 $1$ 单位时量子点 $i$ 电化学势的一阶响应”。在常相互作用模型框架下，栅压与电化学势本就是线性关系，因此 $\mathbf{T}$ 在整个工作区内都是常数（不限于一阶近似）；势垒栅对隧穿耦合的影响则呈指数律，不在 $\mathbf{T}$ 的管辖范围。这一点决定了 $\mathbf{T}$ 只能独立消除电化学势方向的串扰，势垒方向必须另立”虚拟势垒栅极”。

从 CI 模型到串扰矩阵

把 $n$ 量子点的静电能写成多维电路的形式

$$\mathbf{Q}=\mathbf{CV},U=\frac{1}{2}{\mathbf{V}}^{\mathrm{T}}\mathbf{CV}$$

其中电容矩阵 $\mathbf{C}$ 按”量子点”与”电极”分块。计入门控栅与点间互电容后，对一个双量子点

$$\mathbf{C}=\left(\begin{array}{cc}{C}_1& -C_{12}\\ -C_{21}& C_2\end{array}\right)$$

其对角元 $C_i=C_{L(R)}+C_{g,i}+C_m$ 是点 $i$ 的等效总电容，非对角元 $-C_{ij}$ 则是点 $i$ 与点 $j$ 的互电容。这是 CI 模型”电容网络”视角下的精确结果。要把它转成”栅压 → 电化学势”的映射，只需对 ${\mu }_i=\partial U/\partial N_i$ 求导并代入节点电荷与栅压的线性关系，即可得到

$$\Delta \mathbf{\mu }=\alpha \mathbf{T}\Delta \mathbf{V}$$

其中 $\alpha$ 是与器件几何相关的整体杠杆臂因子，$\mathbf{T}$ 是无量纲的归一化矩阵。这就是把经典电路图转译为量子点控制语言的标准通道。对双点再写一次分量形式，可得

$$\Delta {\mu }_1=\alpha (t_{11}\Delta V_1+t_{12}\Delta V_2),\Delta {\mu }_2=\alpha (t_{21}\Delta V_1+t_{22}\Delta V_2)$$

矩阵元 $t_{ij}=C_{g,j}{C}_{ij}^{-1}/{\sum }_k{C}_{g,k}$ 等无量纲比例，其中 $C_{g,k}$ 是栅极 $k$ 与各点的互电容展开，$C^{-1}$ 是 $\mathbf{C}$ 矩阵的逆。这与在 Si/SiGe 2×2 阵列中采用的写法 $\mathbf{vG}=\mathbf{MG}$ 等价——只是坐标系约定略有差别：把”对角归一”留给变换矩阵、把”杠杆臂”放进变换矩阵或显式提出都不会影响本质。

对角归一化约定

实验上一般把对角元固定为 $1$，让矩阵的非对角元直接等于串扰强度。具体做法是取 $\Delta {\mu }_i$ 对 $V_i$ 的响应为单位响应，定义

$$t_{ii}=1,t_{ij}=\frac{\partial {\mu }_i/\partial V_j}{\partial {\mu }_i/\partial V_i}(j\ne i).$$

归一化后矩阵 $\mathbf{T}$ 的逆 ${\mathbf{T}}^{-1}$ 给出一组虚拟电极变换，把”只移动一个 ${\mu }_i$ 而不扰动其他 ${\mu }_j$“变成若干物理栅压的同时变化。这也是 2×2 阵列中 11×11 矩阵所采用的约定——7 行势垒栅 + 4 行能级栅，对应 11 根物理栅极串扰到 4 个电化学势的耦合：前 7 行记录势垒栅对势垒的串扰（点间隧穿耦合不参与时除对角线外为零），后 4 行记录所有栅极对 4 个 ${\mu }_i$ 的耦合。

从隧穿线斜率提取矩阵元

非对角元不需要逐点拟合电容，可直接从电荷稳定图读出。关键观察是：沿量子点 $i$ 的电荷隧穿线移动时，该点的电化学势不变，即 $\Delta {\mu }_i=0$。代入线性映射并在只扫描两根栅极（横轴 $p$、纵轴 $q$）的二维相图中化简：

$$0=\Delta {\mu }_i=\alpha (t_{i,q}\Delta V_q+t_{i,p}\Delta V_p).$$

约定对角归一后立即得到

$$t_{i,p}=-\frac{\Delta V_q}{\Delta V_p}\equiv -r_i$$

其中 $r_i$ 即量子点 $i$ 的隧穿线斜率。因此整个矩阵可以通过一次双点相图获得：在 (1,1) 电荷区附近，对每根”非目标”栅极 $j$ 与目标点柱塞栅 $i$ 一起扫描稳定图，从隧穿线斜率读出 $t_{i,j}$。在双量子点器件上给出完整例子：测得 $r_1=-8$、$r_2=-0.325$，于是

$$\mathbf{T}=\left(\begin{array}{cc}1& 1/8\\ -0.325& 1\end{array}\right)=\left(\begin{array}{cc}1& 0.125\\ 0.325& 1\end{array}\right)$$

反演得

$${\mathbf{T}}^{-1}=\left(\begin{array}{cc}1.042& -0.13\\ -0.339& 1.042\end{array}\right)$$

扫描虚拟电极而非物理栅极重新测量，同一对电荷隧穿线即变正交，“独立控制”立即获得实验验证。在 Si/SiGe 2×2 阵列中则把 (1,1,1,1) 电荷区作为参考状态，逐栅扫描 $V_{B_i/P_i}$ 与目标柱塞栅 $V_{P_i}$ 测电荷稳定图，对每根隧穿线拟合峰值位置得到 ${\alpha }_{B_1-P_1}=|k_{B_1-P_1}|$ 等系数，依次填入矩阵；整个串扰矩阵 $\mathbf{M}$ 的维度为 $11\times 11$，前 7 行描述势垒栅对势垒的串扰，后 4 行描述所有栅对电化学势的串扰。

矩阵反演、条件数与正则化

虚拟电极的存在完全依赖于 $\mathbf{T}$ 可逆。实验上需要关注的几个量：

• 条件数 $\kappa (\mathbf{T})={\sigma }_{\max }/{\sigma }_{\min }$（最大与最小奇异值之比）直接衡量控制轴的线性无关程度。$\kappa$ 越大，矩阵越接近奇异，逆变换会把测量噪声放大为控制动作的剧烈起伏；
• 正则化伪逆 ${\mathbf{T}}^{+}=({\mathbf{T}}^{\mathrm{T}}\mathbf{T}+\lambda \mathbf{I}{)}^{-1}{\mathbf{T}}^{\mathrm{T}}$ 是处理接近奇异矩阵的标准手段：引入 Tikhonov 正则化参数 $\lambda$ 后，最小二乘意义下的小特征值被压缩，反演不再发散，但虚拟电极对相应方向的控制会被牺牲；
• 重新选择栅极组合是更彻底的方案：去掉高度相关的栅、改用更稀疏的电极版图，或显式引入屏蔽栅；
• 迭代更新：当工作点漂移、电容随栅压变化时，旧矩阵需要重新测量。附录给出 ${\mathbf{M}}_3={\mathbf{M}}_2{\mathbf{M}}_1$ 的迭代更新关系，新矩阵等于新测量矩阵与旧矩阵的乘积，这一乘法更新结构在多点阵列中只引入线性测量开销。

多点扩展：迭代更新的物理根据

把虚拟电极扩展到 $n$ 个量子点时不需要一次性测量完整 $n\times m$ 矩阵。相邻双点依次建立局部变换矩阵后，按乘法复合即可：

$${\mathbf{M}}_n={\mathbf{M}}_{\mathrm{update}}{\mathbf{M}}_{n-1}$$

其物理基础是：在对新点更新之前，旧虚拟电极对新点的影响与对应物理栅极完全相同。指出：“在虚拟电极没有对 P3 进行更新之前，虚拟电极 VP2 和物理栅极 P2 对 QD3 的影响是相同的，这意味着在新加入量子点后……可以直接使用物理栅极下测量的电荷稳定图计算变换矩阵并使用式 4.12 的关系进行更新”，从而避免对已有栅极全部重新测量。这一递推结构是阵列可扩展性的关键：在四量子点器件上，自动程序给出的实际变换矩阵示例是

$$\mathbf{T}=\left(\begin{array}{cccc}1& 0.125& 0& 0\\ 0.325& 1& 0.167& 0\\ 0& 0.091& 1& 0.257\\ 0& 0& 0.27& 1\end{array}\right).$$

为简化操作，程序只针对四条柱塞栅测量对应电荷稳定图，势垒栅串扰因调控范围小而被忽略——这是工程上对”完整矩阵 vs. 测量开销”的常见折中。

势垒方向：为什么串扰矩阵救不了虚拟势垒

电化学势方向的串扰是电容性的、近似线性，所以 $\mathbf{T}$ 可以干净地消除。但势垒栅对隧穿耦合的调控不是线性的，而是强指数依赖：

$$t_{ij}=t_0+t_1\exp \left({\beta }_{ij}{\mathrm{vB}}_{ij}\right)$$

其中 ${\beta }_{ij}$ 表征该势垒栅的调控能力。在 2×2 阵列中测得四个最近邻势垒的 ${\beta }_{ij}$ 从 $1.42\times 10^{-2}$ 到 $6.85\times 10^{-2}{\mathrm{mV}}^{-1}$ 不等，最弱者仅为最强者的 $21\%$。此时交叉调控系数 ${\beta }^{\prime }$ 不能由 $\mathbf{T}$ 给出，必须从数据另行提取。例：${\mathrm{vB}}_{41}$ 增大会把 QD4 推离 QD3 而显著压低 $t_{34}$，对应的 ${\beta }^{\prime }=-1.03\pm 0.11\times 10^{-2}{\mathrm{mV}}^{-1}$；据此在扫 ${\mathrm{vB}}_{41}$ 时按此系数补偿 ${\mathrm{vB}}_{34}$，即可维持 $t_{34}$ 不变——这是指数型”虚拟势垒栅极”的手动实现，$\mathbf{T}$ 完全不参与。

参数与量级

量	典型值	来源
矩阵维度	$11\times 11$（Si/SiGe 2×2：7 势垒 + 4 柱塞）
隧穿线斜率实例 $r_1,r_2$	$-8$、$-0.325$（GaAs 双点）
非对角元量级	$0.091$–$0.325$（相邻点串扰约 $9\%$–$33\%$）	（式 4.14）
反演矩阵元	$1.042$、$-0.13$、$-0.339$	（式 4.10）
串扰系数 ${\alpha }_{B_1-P_1}$	由隧穿线斜率绝对值给定
平均杠杆臂 $\alpha$	$\approx 0.12eV/V$
势垒调控系数 ${\beta }_{ij}$	$(1.42$–$6.85)\times 10^{-2}{\mathrm{mV}}^{-1}$
耦合串扰补偿系数 ${\beta }^{\prime }$	$-1.03\pm 0.11\times 10^{-2}{\mathrm{mV}}^{-1}$（vB41→t34）

注意非对角元 $0.091$–$0.325$ 表示相邻栅极在另一量子点上引起约 $9\%$–$33\%$ 的串扰；这种串扰在三维 Si 量子点阵列中会进一步放大，是矩阵条件数偏大的主要原因。

实验特征与测量

标准测量流程

测量流程以阵列已稳定在 (1,1,1,1) 等参考电荷区为前提：在 Wang Ning 的 2×2 阵列中，其余量子点被置于深阻塞区以避免干扰；对 QD1 的串扰系数，依次扫描 $V_{B_i/P_i}$ 与 $V_{P_1}$ 测稳定图，从隧穿线斜率拟合 ${\alpha }_{B_i/P_i-P_1}=|k_{B_i/P_i-P_1}|$；其余三个量子点用相同方式处理，最终拼成完整矩阵。锁定放大器参数典型为积分时间 $300\mathrm{ms}$、激励频率 $73\mathrm{Hz}$、幅值 $0.5\mathrm{mV}$（分压后），两台锁定放大器分别配给两个 SET 并做信号同步以提高探测效率。SET 直流偏置 $0.2\mathrm{mV}$，交流激励 $200\mu \mathrm{V}$，可在保证信噪比的同时抑制功率展宽。

自动化：霍夫线变换

人工逐条拟合斜率在大阵列中不现实。霍夫线变换（Hough line transform）是标准自动方案：把笛卡尔坐标系下的每个像素 $(x_i,y_i)$ 映射到极坐标

$$\rho =x_i\cos \theta +y_i\sin \theta$$

同一直线上的点在 $(\rho ,\theta )$ 平面上相交于一点，统计相交曲线最多的 $(\rho ,\theta )$ 即给出原图中直线位置与斜率。该算法对残缺或被遮挡的直线同样鲁棒，特别适合从含噪声、电荷隧穿线段被传感器信号截断的实测相图中提取斜率。霍夫线变换把虚拟电极建立完全无人工化，与电荷态识别的卷积神经网络串联，即构成自动调控流水线的两个核心模块。

矩阵反演的实验判据

虚拟电极是否生效不靠理论计算，要看电荷稳定图：隧穿线由倾斜变正交即直接证明控制轴已正交化。同时物理栅压电压范围、积分时间和最大电压限制都必须重新评估——虚拟电极下的电压绝对值没有物理意义，只有变化量 $\Delta {\mathrm{vP}}_i$ 是真实的； 2×2 阵列测得所有四个量子点的电荷隧穿线在虚拟电极下均与坐标轴平行，等价于 ${\mathbf{M}}^{-1}$ 的实验确认。

失效与维护

虚拟电极是局部线性化的产物，其有效性有明确边界：

• 工作点漂移：栅压响应随电子数、势垒状态变化，线性矩阵只在建立时的工作点附近成立。增大隧穿耦合会使量子点位置明显移动，原本正交的电荷隧穿线重新倾斜—— 5.A 附录明确指出：“因为耦合强度的增加使得量子点的位置发生较大变化，原来的虚拟栅极变得不那么有效”；
• 周期性重标定：必须重新测量串扰系数并按 ${\mathbf{M}}_n={\mathbf{M}}_{n-1}{\mathbf{M}}_{\mathrm{update}}$ 复合更新；阵列运行时间越长、耦合调得越强，更新越频繁；
• 斜率法与几何法偏差：库仑菱形三种杠杆臂提取方法不一致时，往往提示栅压在扫描范围内改变了点的形状（电容随栅压漂移）或存在未被计入的交叉电容；
• 方向选择：电化学势方向用 $\mathbf{T}$ 即可；势垒方向因指数律无法靠同一矩阵描述，必须配合逐点提取的 $\beta$ 系数手动补偿。

与其他概念的关系

• 常相互作用模型给出”电化学势对栅压为线性”的微观根据；$\mathbf{T}$ 就是这条线性依赖的矩阵化表达。Hubbard 模型参数 $U_i$、$U_{ij}$ 也用 CI 模型从蜂窝图标定，串扰矩阵则是同一网络在”控制侧”的对应物；
• 虚拟电极就是 ${\mathbf{T}}^{-1}$ 的行向量：给定矩阵，虚拟电极立即写出来；矩阵条件数过大时控制轴近线性相关，反演会放大噪声；
• 电荷稳定图上的隧穿线斜率直接给出 $\mathbf{T}$ 的非对角元；虚拟电极生效的判据也画在同一张图上（隧穿线变正交）；
• 电化学势方向的线性来自 CI 模型；势垒方向由隧穿耦合的指数律支配，是 $\mathbf{T}$ 唯一管不到的地方；
• 量子点阵列的调控难度按点数线性增长，串扰矩阵是该线性增长的”记账本”：它的维度随阵列规模线性上升，条件数随结构非对称性放大；
• 自动调控靠卷积神经网络读图、霍夫线变换提斜率、用 ${\mathbf{T}}^{-1}$ 写虚拟电极三步串联完成，因此 $\mathbf{T}$ 是自动调控从”看图”到”拧电压”之间必须有的接口；
• 读出串扰发生在测量侧（频分复用色散读出中圆斑整体平移），交叉电容矩阵发生在控制侧（栅压同时牵动多个 ${\mu }_i$），二者都是多比特扩展必须面对但完全独立的两类串扰；
• Si/SiGe 异质结器件中二维阵列要靠四层铝栅（屏蔽 / 能级 / 两层势垒）把柱塞、势垒、屏蔽三类电极扇出，正是为减少 $\mathbf{T}$ 的非对角元；COMSOL 静电模拟显示屏蔽板的非对称性是 2×2 阵列四点充电能 $U_1/U_2/U_3/U_4=2.83/5.02/3.05/4.63\mathrm{meV}$ 不均匀的直接来源——也即同一矩阵在不同点有不同的对角比例。