来源论文: https://arxiv.org/abs/2607.08447v1 生成时间: Jul 10, 2026 13:11

硅基FinFET双量子点空穴自旋比特中的交换相互作用模拟：基于Schrödinger-Poisson与组态相互作用理论的深度解析

0. 执行摘要

随着硅基半导体工艺向纳米尺度不断演进，利用互补金属氧化物半导体（CMOS）兼容工艺制备的自旋量子比特（Spin Qubits），因其极高的可扩展性与超长的相干时间，已成为实现大规模通用量子计算（Fault-Tolerant Quantum Computing）最具前景的技术路线之一。然而，传统的二维平面电子自旋量子比特由于本征自旋-轨道耦合（Spin-Orbit Coupling, SOC）较弱，通常需要引入结构复杂且体积庞大的外部微磁体（Micromagnets）来实现电偶极自旋共振（EDSR）控制，这极大地限制了其高密度集成与扩展性。此外，平面二维电子气（2DEG）波函数易受界面电荷噪声与粗糙度干扰，导致消相干效应显著。

为了克服上述瓶颈，学术界和工业界近年来将目光投向了局域在三维非平面结构（如纳米线、鳍式场效应晶体管 FinFET）中的空穴自旋量子比特（Hole Spin Qubits）。空穴作为价带波函数的载流子，天然具有强烈的、可电场调控的本征自旋-轨道耦合，能够实现全电学的超快自旋操控，且其空间波函数在三维受限几何结构中展现出独特的重空穴（Heavy-Hole, HH）与轻空穴（Light-Hole, LH）混合特性，使自旋 $g$ 因子呈现强烈的各向异性，从而为自旋量子比特的操纵提供了极高的自由度与优化的“甜点（Sweet Spots）”。

本篇博文深度解析了苏黎世联邦理工学院（ETH Zurich）集成系统实验室（IIS）Ilan Bouquet、Alexander Maeder与Mathieu Luisier等人的最新研究成果。该工作基于自洽的三维薛定谔-泊松（Schrödinger-Poisson, SP）方程求解器与多粒子组态相互作用（Configuration Interaction, CI）方法，系统模拟了在一款包含5个金属栅极的硅鳍式场效应晶体管（Si FinFET）结构中，双量子点（Double Quantum Dot, DQD）体系内空穴自旋比特间的交换相互作用（Exchange Coupling, $J$）。

研究团队利用 GPU 加速算法，在考虑了极低温（$1.5\text{ K}$）热收缩引起的非均匀应变场（Strain Distribution）以及外加磁场（Zeeman 劈裂）的背景下，实现了对空穴双粒子态哈密顿量的高效对角化。研究表明，通过精确调制势垒栅极电压（Barrier Gate Voltage, $V_B$）以及两个柱塞栅极（Plunger Gates）之间的失谐电压（Detuning, $\varepsilon_{VP}$），可以实现交换相互作用 $J$ 在数百兆赫兹（MHz）范围内的宽幅自洽调制。同时，论文进一步揭示了空穴自旋量子比特在特定的磁场取向与电场配置下，存在横向交换作用（$J_\perp$）完全消失而纵向交换作用（$J_\parallel$）保持有限的物理区间，从而为实现高保真度（High-Fidelity）受控相位门（CZ门）与交换门（SWAP门）提供了决定性的器件物理指导。这一理论与数值计算框架不仅完美地复现了当前的实验测量趋势，更为下一代高性能、可扩展的硅基空穴量子处理器设计奠定了坚实的定量化物理模拟基础。

1. 核心科学问题，理论基础，技术难点，方法细节

1.1 核心科学问题：三维空穴量子点中的交换耦合调控机制

在基于双量子点（DQD）的电荷或自旋体系中，实现两比特量子逻辑门（如 $CZ$ 或 $SWAP$）的核心在于精确控制两个邻近量子点中局域粒子之间的交换耦合相互作用 $J$。从物理本质上看，交换耦合源自于量子力学中的泡利不相容原理与粒子间的库仑排斥力的协同效应：

其中 $E_{\text{Triplet}}$ 和 $E_{\text{Singlet}}$ 分别为双粒子体系在三线态（Triplet）与单线态（Singlet）下的本征能量。对于电子量子点体系，由于自旋-轨道耦合通常极弱，交换作用往往可以很好地用简单的费米-哈伯德（Fermi-Hubbard）模型或海森堡各向同性自旋模型进行描述。然而，在硅基空穴自旋量子比特中，问题变得复杂得多：

1. 强自旋-轨道耦合与多带混合：价带顶（$\Gamma$ 点）由重空穴（HH，自旋投影 $J_z = \pm 3/2$）和轻空穴（LH，自旋投影 $J_z = \pm 1/2$）多带构成。三维纳米级受限（Confinement）和降温应变导致 HH 和 LH 发生强烈的空间重组与能级杂化。这导致空穴自旋波函数具有显著的空间各向异性与非共线自旋结构。
2. 各向异性交换相互作用：自旋-轨道耦合的存在破坏了自旋空间的各向同性，使得交换相互作用矩阵不再是标量，而是表现为包含横向交换项 $J_\perp$（驱动自旋翻转即 flip-flop 过程）和纵向交换项 $J_\parallel$（驱动 Ising 型相位移动）的各向异性张量。如何通过栅极偏压和磁场取向自如地控制这一各向异性张量，是两比特门设计的核心科学问题。

1.2 理论基础与多尺度模拟框架

为了精确捕获上述复杂的物理效应，必须构建一套从单粒子微观波函数自洽计算到多粒子相互作用本征态求解的完备数值模拟流程。本工作所构建的理论框架由以下三个核心模块组成：

1.2.1 自洽薛定谔-泊松（Schrödinger-Poisson）方程求解器

在极低温下，器件的电荷分布与静电势呈现强烈的非线性自洽特征。模拟首先在包含源漏延伸区和沟道的全三维晶格上交替求解三维薛定谔方程与泊松方程：

其中 $p(\mathbf{r})$ 为局域空穴浓度，由薛定谔方程解出的占据态波函数模方叠加决定：

$f$ 为费米-狄拉克分布函数。薛定谔方程中采用的非微扰 Hamiltonian 为包含了应变效应、外加磁场效应的 $6 \times 6$ Luttinger-Kohn (LK) $k \cdot p$ 理论模型：

其中 $H_{\text{LK}}$ 描述了重空穴、轻空穴及分裂带（Split-Off band, SO）之间的能带结构、动量算符 $\mathbf{k} = -i\nabla + \frac{e}{\hbar}\mathbf{A}$（$\mathbf{A}$ 为磁矢势）与静磁场 $\mathbf{B}$ 引起的塞曼耦合（Zeeman coupling）。$H_{\text{PB}}$ 为 Pikus-Bir 应变哈密顿量，用以定量计算由于硅鳍（Si Fin）与周围二氧化硅（$\text{SiO}_2$）热膨胀系数差异在 $1.5\text{ K}$ 极低温下产生的多维热应变张量 $\boldsymbol{\varepsilon}$。

1.2.2 基底变换与自旋-轨道耦合的显式分离

直接对 $6 \times 6$ 的 $k \cdot p$ 薛定谔方程进行双粒子多态膨胀会带来极高维数的计算灾难，且传统的单粒子状态极易受到“自旋污染（Spin Contamination）”的干扰。为此，本方法执行了一个全局酉变换（Rotation $U$），将 LK 表象下的角动量基底 $\left| J, m_j \right\rangle$ 变换至轨道角动量 $L$ 与自旋 $S$ 的直积基底 $|LS\rangle = \{|X\Uparrow\rangle, |Y\Uparrow\rangle, |Z\Uparrow\rangle, |X\Downarrow\rangle, |Y\Downarrow\rangle, |Z\Downarrow\rangle\}$：

变换后，哈密顿量被优雅地解耦为两部分：

• $H_{\text{DKK}}$ (Dresselhaus-Kittel-Kip 哈密顿量)：它是不包含自旋翻转的纯轨道部分。在无磁场时，它退化为两块完全等价且独立的 $3 \times 3$ 矩阵，对应自旋向上和自旋向下的对称能级。这一特性允许我们在不考虑自旋的情况下，极其干净地抽取空间单粒子基底。
• $H_{\text{SOI}}$ (本征自旋-轨道耦合项)：其显式地包含了所有自旋与轨道角动量的相互耦合作用。

此时，非相互作用双量子点哈密顿量可以重写为：

其中 $C$ 为双空穴粒子间的库仑相互作用项。通过这一变换，自旋-轨道耦合项 $H_{\text{SOI}}$ 被归入微扰项 $H_1$ 中，与两粒子库仑相互作用统一处理，保证了数值基底的严谨性。

1.2.3 组态相互作用（Configuration Interaction, CI）方法与两粒子基底构建

为了计算包含多粒子库仑关联与电荷流动的量子点状态，采用组态相互作用方法。首先通过对角化不含自旋-轨道耦合的单粒子哈密顿量 $H_0$，在预先算好的自洽势场 $V_{\text{DQD}}$ 下解出前 $M$ 个双重简并的单粒子本征态 $\{\varphi_i\}$（例如基态波函数局域在左、右量子点的轨道组合：$|\Uparrow\rangle_L, |\Downarrow\rangle_L, |\Uparrow\rangle_R, |\Downarrow\rangle_R$）。

随后，通过选择 $n_p = 2$ 个空穴填充到这 $M$ 个单粒子本征轨道中，构建反对称化的多粒子**斯莱特行列式（Slater Determinants）**基底。基底的总数（组态数）$K$ 由组合数决定：

在最精简的模型中，取最靠近费米能级的两对（$M=4$）双重简并态，可产生 $K=6$ 个基本的两粒子斯莱特组态：

其中，下标 $LL$ 和 $RR$ 代表双粒子局域在单个点上的高能物理状态（通常对应两粒子电荷态 $(2,0)$ 和 $(0,2)$）；而 $LR$ 表征两个量子点各占一个空穴的对称电荷态 $(1,1)$。

利用斯莱特-康登规则（Slater-Condon Rules），可以将多体 Hamiltonian 算符 $H_{\text{DQD}} = H_0 + H_SOI + C$ 投影到上述 $K$ 维斯莱特行列式空间中，构建显式的 CI 哈密顿量矩阵 $H_{\text{CI}}$。通过在不同栅极电压（也即不同失谐度 $\varepsilon_{VP}$）下，反复执行上述基底构建与矩阵对角化，便能获取精密的双空穴多体能谱图以及各能级随电压的演化关系。

1.3 技术难点及关键解决方法

• 重/轻空穴混合波函数的高维库仑积分计算：组态相互作用中的核心难点在于计算复杂的双电子库仑积分： $$\langle \varphi_i \varphi_j | C | \varphi_k \varphi_l \rangle = \frac{e^2}{\varepsilon} \iint \frac{\varphi_i^*\left(\mathbf{r}_1\right) \varphi_j^*\left(\mathbf{r}_2\right) \varphi_k\left(\mathbf{r}_1\right) \varphi_l\left(\mathbf{r}_2\right)}{\left|\mathbf{r}_1 - \mathbf{r}_2\right|} d\mathbf{r}_1 d\mathbf{r}_2$$ 若采用直接的 6 维空间数值积分，极高的数据维度将使计算效率极其低下。作者巧妙地采用泊松方程法化简此过程。通过定义复密度函数 $\rho_{jl}^{\sigma'}\left(\mathbf{r}_2\right) = \varphi_j^{*\sigma'}\left(\mathbf{r}_2\right) \varphi_l^{\sigma'}\left(\mathbf{r}_2\right)$ 并求解对应的三维泊松方程算得等效静电势 $V_{jl}^{\sigma'}\left(\mathbf{r}_2\right)$： $$\nabla^2 V_{jl}^{\sigma'} = -\frac{e}{\varepsilon} \rho_{jl}^{\sigma'}$$ 之后，六维积分便可完美降维至普通的三维空间重叠积分： $$\langle \varphi_i \varphi_j | C | \varphi_k \varphi_l \rangle = \int V_{jl}^{\sigma'}\left(\mathbf{r}_1\right) \rho_{ik}^{\sigma}\left(\mathbf{r}_1\right) d\mathbf{r}_1$$ 此举使双粒子关联相互作用的计算速度提升了数个数量级。

2. 关键 benchmark 体系，计算所得数据，性能数据

2.1 模拟目标：5-栅极硅基三角鳍式场效应晶体管 (Si FinFET)

研究团队选择了一款与国际前沿实验（如 Basel 大学 Geyer 等人的工作）高度一致的 **5-栅极 p-型 硅鳍式场效应晶体管（FinFET）**作为 Benchmark 器件，如图 1 所示。其关键几何与物性参数如下：

• 沟道材料：高度受限的三角形截面硅鳍，高度 $h_C = 28\text{ nm}$，底部宽度 $w_{\text{Fin}} = 15\text{ nm}$。
• 衬底与掺杂：轻 $p$ 型掺杂硅衬底（掺杂浓度 $10^{14}\text{ cm}^{-3}$）。
• 绝缘层与栅极材料：厚度为 $h_{\text{Ox1}} = 4\text{ nm}$ 的 $\text{SiO}_2$ 栅介质层，栅极采用氮化钛（$\text{TiN}$）金属电极。栅极之间采用 $4\text{ nm}$ 的二氧化硅间隙进行电学隔离。
• 五栅极拓扑配置：包含两层控制构造。
  • 柱塞栅极 $V_{P1}$ 与 $V_{P2}$：长度均为 $10\text{ nm}$，用于直接在其下方电学诱导两个局域化量子点，并施加偏压调制 DQD 系统能量失谐度 $\varepsilon_{VP} = V_{P1} - V_{P2}$。
  • 势垒栅极 $V_B$：长度 $27\text{ nm}$，位于中间，用于调节两量子点之间的势垒高度，进而定量调控波函数重叠度与交换耦合强度 $J$。
  • 外侧侧栅极 $V_L$ 与 $V_R$：长度 $14\text{ nm}$，偏置在强积累区以连接源漏欧姆接触区。

2.2 核心物理计算数据与物理机制分析

2.2.1 强应变背景下的重、轻空穴混合比例（HH-LH Mixing）

在极低温 $T = 1.5\text{ K}$ 下，由于热失配导致硅鳍受到强烈的多维各向异性挤压。计算得到的 DQD 区域价带组分分布（如图 3 所示）表明：

• 在对称势场偏置下（$V_{P1}=V_{P2}=0.484\text{ V}$, $V_B=0.690\text{ V}$），基态空穴波函数的能带组分呈现出高度统一的各向异性分布：轻空穴（LH）占据绝对的主导地位，占比高达 ~82%；而重空穴（HH）占比仅为 ~15%；本征自旋-轨道耦合（SOC）贡献了剩余的 ~3%。
• 物理机制：强烈的 LH 占据主导主要归因于两个层面：其一，三维三角形截面硅鳍在水平和竖直方向上具有极强的空间受限特征；其二，由于低温冷却产生的各向异性压应变拉大了轻重空穴的能级劈裂，使 LH 能级在能量上更靠前，这为自旋自洽调制提供了极为稳定的轻空穴相干态运行空间。

2.2.2 充电能（Charging Energy, $E_C$）测量与 CI 空间截断

通过自洽泊松场模拟，计算出在 $T=1.5\text{ K}$ 下，将第二个空穴添加进该双量子点系统所需的充电能（或称相干自能，能谱图可见图 3(b) 中的双箭线标识）：

这一数据与实验中在稳定性图（Stability Diagram，即 Honeycomb 蜂窝图）上测得的单量子点电荷状态演化能谱高度吻合。因此，在 CI 的对角化计算中，为了极大缩短处理能谱的耗时，可以采用一个物理截断：将占据能量高于 $E_C$ 的状态通过费米-狄拉克过滤器滤除，最终选取的紧凑基底组态包含四个单粒子态（$M=4$），CI 哈密顿量矩阵维度为 $K=28$，在保证计算精度误差在 $\mu\text{eV}$ 级别的同时，极大地缩减了高算力硬件的负载。

2.2.3 栅极电压对交换相互作用 $J$ 的调控性能评估

图 6 系统展示了利用该物理模型对交换耦合频率（$J/h$）的定量调制曲线：

• 势垒栅极阻尼效应：在固定柱塞失谐度 $\varepsilon_{VP} = 20\text{ m{V}}$ 时，微调中间势垒栅极电压 $V_B$。当 $V_B$ 从 $687\text{ m{V}}$ 提高到 $690\text{ m{V}}$（仅 3 mV 的幅度）时，两点之间空穴波函数的空间重叠程度明显降低，交换耦合频率 $J/h$ 从 $525\text{ MHz}$ 骤降至 $350\text{ MHz}$。这表明势垒栅对交换耦合强度具备高灵敏度的指数级电场调制能力。
• 失谐度控制行为：在不同固定 $V_B$ 偏置下，通过改变两端注入电压失谐度 $\varepsilon_{VP}$（从 $6\text{ m{V}}$ 到 $16\text{ m{V}}$），交换相互作用 $J/h$ 随失谐度单调上升，且表现出光滑无跳变的理想多体电学行为，证实了该器件在极宽的频域中具有良好的两比特门可控性。

3. 代码实现细节，复现指南，所用的软件包及开源 repo link

3.1 核心算法：Rayleigh-Chebyshev 子空间迭代法与 GPU 加速

由于三维非均匀网格自洽薛定谔-泊松（SP）循环中涉及的稀疏矩阵维数通常在 $10^5 \sim 10^6$ 级别，采用传统的单处理器 ARPACK（通常基于 Arnoldi 算法）对角化求解不仅极度消耗内存，而且难以应付需要扫参成百上千次的失谐度计算。本工作通过使用基于 Rayleigh-Chebyshev 子空间迭代算法（Rayleigh-Chebyshev Subspace Iteration, RCSI）的大规模并行解法，在 CPU 和 GPU 上均实现了卓越的算力飞跃：

• 计算性能优势：相较于 CPU 上的标准 ARPACK 串行算法，本算法在多核 CPU 上实现了 $2.9\times$ 的计算加速，而在搭载 NVIDIA CUDA 核心的 GPU 硬件上则取得了 $20\times$ 以上的惊人加速效果。这使得 DQD 整个电荷流扫参模拟在数小时内即可完全收敛。

3.2 模拟流程与复现指南伪代码

为了便于研究人员基于开源或商业平台复现该方法，以下提供自洽 SP 循环结合多体 CI 计算的完整算法流程图与复现指南：

# 硅基空穴DQD交换相互作用(CI-SP)模拟复现框架伪代码
import numpy as np
from scipy.sparse import linalg as spla

class HoleDQD_CI_Simulator:
    def __init__(self, grid_dim, dx):
        self.grid_dim = grid_dim  # (Nx, Ny, Nz)
        self.dx = dx              # 网格间距 dx=dy=dz=0.5 nm
        self.V_ext = None         # 栅极电压决定的静电势边界条件
        self.strain_tensor = None # 热收缩应变张量 (Pikus-Bir)
        
    def solve_strain_field(self):
        """
        基于有限元法(FEM)求解器件由300K降温至1.5K时的多维非均匀热应变分布
        """
        print("Solving 3D elastic deformation equations...")
        self.strain_tensor = solve_isotropic_elasticity(boundary_condition="bottom_fixed")
        return self.strain_tensor

    def schrodinger_poisson_loop(self, V_gates, max_iter=100, tol=1e-5):
        """
        求解自洽的3D薛定谔-泊松方程，获取单粒子状态及DQD局域势能场
        """
        V_DQD = np.zeros(self.grid_dim)
        p_charge = np.zeros(self.grid_dim)
        
        for i in range(max_iter):
            # 1. 求解含有应变与外磁场项的6x6 Luttinger-Kohn Hamiltonian
            H_lk = construct_6x6_kp_operator(self.strain_tensor, V_DQD, B_field=[0,0,0.1])
            # 使用Rayleigh-Chebyshev子空间迭代法计算前N个本征态
            eigenvals, wavefuncs = rcsi_solver(H_lk, num_states=10)
            
            # 2. 计算空穴电荷密度分布
            p_charge_new = compute_hole_density(wavefuncs, eigenvals, T=1.5)
            
            # 3. 求解泊松方程，更新DQD静电势
            V_DQD_new = solve_poisson_equation(p_charge_new, V_gates)
            
            # 4. 判断自洽收敛
            if np.max(np.abs(V_DQD_new - V_DQD)) < tol:
                print(f"SP Loop converged at iteration {i}.")
                return V_DQD_new, wavefuncs, eigenvals
            V_DQD = 0.8 * V_DQD + 0.2 * V_DQD_new  # 混合阻尼
        raise TimeoutError("SP Loop did not converge.")

    def run_configuration_interaction(self, V_DQD, M_states=4):
        """
        利用组态相互作用方法对角化双粒子 Hamiltonian，计算交换相互作用 J
        """
        # 1. 提取不含SOI的轨道部分H0，对角化获取正交单粒子空穴轨底
        H0 = construct_orbital_H0(self.strain_tensor, V_DQD)
        single_particle_states = diagonalize_H0(H0, num_states=M_states)
        
        # 2. 构建两粒子斯莱特行列式(Slater Determinants)
        slater_states = build_slater_basis(single_particle_states, num_particles=2)
        K = len(slater_states)
        H_CI = np.zeros((K, K))
        
        # 3. 填充CI矩阵：应用Slater-Condon规则以及泊松辅助库仑积分
        for a in range(K):
            for b in range(a, K):
                # 计算单体算符项 (H0 + H_SOI)
                H_CI[a, b] = compute_one_body_matrix_element(slater_states[a], slater_states[b])
                # 计算双体算符项：三维静电势辅助计算快速库仑相互作用 C
                H_CI[a, b] += compute_two_body_coulomb_element_via_poisson(slater_states[a], slater_states[b])
                H_CI[b, a] = H_CI[a, b]
                
        # 4. 对角化CI哈密顿矩阵获取双粒子多体本征能谱
        ci_energies, ci_eigenstates = np.linalg.eigh(H_CI)
        
        # 5. 提取交换耦合强度 J = E_Triplet - E_Singlet
        J = ci_energies[1] - ci_energies[0]  # 对于空穴体系，通常取反，对应Singlet-Triplet间距
        return J, ci_energies, ci_eigenstates

3.3 涉及的模拟软件包及开源链接

为了在本地成功复现该研究成果，推荐结合以下几款知名的学术/工业开源仿真包：

1. QTCAD (Quantum Technology Computer-Aided Design)：该商业工具由加拿大 Nanoacademic Technologies 开发，专门用于半导体自旋量子比特的 3D 自洽薛定谔-泊松以及组态相互作用模拟。其底层有限元方法与应变自洽算法与本工作提出的框架高度契合。
   • 链接：Nanoacademic QTCAD
2. SLEPc (Scalable Library for Eigenvalue Problem Computations)：用于求解超大规模稀疏矩阵本征值问题的开源并行求解器，支持在 GPU（CUDA/HIP）上对雅可比-戴维森或切比雪夫子空间迭代法进行极速运算。是复现自洽 SP 部分的核心软件库。
   • GitHub Repo: SLEPc Project
3. Kwant：一个用于计算强相干紧束缚系统、能带混合以及自旋输运的 Python 库，能够便捷地生成带有自旋-轨道耦合的多基底网格 Hamiltonian。
   • 官方链接：Kwant Project
4. FEniCS Project：针对降温热应变场偏微分方程的三维有限元求解（FEM），可以快速输出本工作所需的精细 Pikus-Bir 应变应力张量分布。
   • 官方链接：FEniCS FEM Platform

4. 关键引用文献，以及你对这项工作局限性的评论

4.1 关键参考文献

1. Loss, D. & DiVincenzo, D. P. Quantum computation with quantum dots. Phys. Rev. A 57, 120 (1998). (该经典文献开创了半导体量子点自旋量子比特作为量子计算载体的理论先河，也是本篇博文探究交换相互作用的基础。)
2. Geyer, S., Hetényi, B. et al. Anisotropic exchange interaction of two hole-spin qubits. Nature Physics 20, 1152 (2024). (本论文所模拟的 5-栅极硅三角鳍式 FinFET 器件在几何构造和物理偏置上，均以此项在 Basel 大学开展的突破性实验为直接物理蓝本。)
3. Szabo, A. & Ostlund, N. S. Modern Quantum Chemistry: Introduction to Advanced Electronic Structure Theory. (Courier Corporation, 2012). (组态相互作用方法（CI）和斯莱特-康登行列式化简计算规则的最经典教科书。)
4. Bir, G. L. & Pikus, G. E. Symmetry and Strain-Induced Effects in Semiconductors. (John Wiley & Sons, 1974). (定量估算硅鳍中非均匀热挤压应变对能带混合调控影响的 Pikus-Bir 哈密顿量物理理论源头。)

4.2 局限性与前沿批判性评论

虽然该仿真计算框架展现出了前所未有的计算物理精度和数值稳定性，但作为一门面向工业级量子芯片设计的研究，其依然存在以下几点不容忽视的局限性：

1. 理想化晶体结构假设与界面粗糙度缺失： 模拟中将 $\text{Si}/\text{SiO}_2$ 界面视为完全平滑的几何面，并对空穴波函数施加了严格的狄利克雷（Dirichlet）边界条件。然而，在真实的实验样品中，原子尺度的界面粗糙度（Surface Roughness）以及栅极氧化物中的悬挂键、电荷陷阱（Charge Traps）等缺陷分布极其普遍。在第三节的实验对比中，真实器件的左、右量子点自旋 $g$ 因子呈现近乎垂直的偏转，而模拟得到的 $g$ 因子则极为对称。这表明界面随机电荷噪声引发的局部对称性破缺对局域自旋结构有极其强烈的影响。未来亟需引入无序模型（Disorder Modeling）来贴近真实的芯片工艺。
2. 缺失对高频电荷噪声（Charge Noise）相干时间损失的动态模拟： 本工作实现了对静态自旋态及交换相互作用 $J$ 的高精度扫描，但实际量子点极易受到电极高频电压抖动或衬底杂质涨落产生的电荷噪声影响。交换作用 $J$ 对门偏压的极高灵敏度实际上是一把“双刃剑”——它在方便电场调控的同时，也会导致 qubit 极易受到电荷噪声的干扰进而失去相干性。该软件平台目前尚未集成动力学消相干（Decoherence Rates）的自洽定量计算功能。
3. CI 基底强截断的完备性问题： 论文为了实现超高效率的扫参，将 CI 的单粒子基底限定在最贴近费米面的 4 个局域态。尽管在该体系中其能级与高阶激发态存在近 $\sim 3\text{ m{eV}}$ 的能垒，但在大失谐区域（$\varepsilon_{VP} > 15\text{ m{V}}$）或者当外加强磁场导致塞曼劈裂非常显著时，高能态（如能谱图中的灰色和浅灰色轨道）与基态组态的杂化程度会明显上升，单方面强行截断可能在特定极端区域带来数个百分点的能级计算偏差。

5. 其他必要的补充：空穴各向异性交换相互作用与 CZ 量子门设计深度剖析

为了更透彻地向科研工作者展示如何利用本工作发现的物理规律来设计两比特量子逻辑门，我们在此补充对各向异性交换张量调制与两比特受控相位门（CZ门）物理构建的深层次探讨。

5.1 各向异性交换相互作用的物理解耦

在强自旋-轨道耦合（SOC）的物理体系中，DQD 双自旋体系的相互作用模型由简单的海森堡各向同性自旋模型 $\mathcal{H} = J(\mathbf{S}_1 \cdot \mathbf{S}_2)$ 演变为各向异性的通用自旋模型。为了在数学和工程上更直观地表述，我们需要将其投影到量子点的局部自旋框架（Qubit Frame）下，可将其解耦为横向交换作用 $J_\perp$ 和纵向交换作用 $J_\parallel$（其数学计算方法已于论文公式 (19a) 和 (19b) 中给出）：

• 横向交换耦合 $J_\perp$：物理上对应于自旋翻转交换，它驱动 $|\Uparrow\Downarrow\rangle \leftrightarrow |\Downarrow\Uparrow\rangle$ 过程。在物理运行中，它是构建 $SWAP$ 和 $\sqrt{SWAP}$ 量子门的物理引擎，但对于 $CZ$ 门，它代表着产生漏码和自旋不相干的有害“泄漏通道”。
• 纵向交换耦合 $J_\parallel$：物理上对应于 Ising 型自旋耦合，它不会引起自旋翻转，而是随着邻近自旋的状态改变而直接产生相应的能级移动（Phase Shift）。

5.2 寻找高保真度 CZ 门的“黄金甜点（Sweet Spots）”

图 7 给出了该模拟最惊艳的物理发现：通过操控外加静磁场 $\mathbf{B}$ 的空间指向（极角 $\theta$ 与方位角 $\phi$），可以在 DQD 体系中自如地调制 $J_\perp$ 与 $J_\parallel$ 的相对强度：

• 传统的困境：在基于非各向异性自旋体系的电子量子点中，要想实现高保真度的 $CZ$ 门，通常需要强磁场引发极大的 Zeeman 能量失谐 $\delta E_Z \gg J$，以便电学压制有害的自旋翻转通道。但这会大大缩小 CZ 门的工作窗口，且自旋退相干随磁场增加而加速。
• 各向异性带来的绝佳机遇：从图 7(c) 和 (d) 可以清晰地发现，在特定的磁场方位角下（如红色圆点标记的 $\theta \approx 65^\circ, \phi \approx 85^\circ$ 区域）：

• 物理后果：在这个神奇的磁场夹角取向（甜点）上，自旋翻转泄漏被大自然以物理几何的形式完全并永久地关闭了！此时即使两比特之间存在较强的交换作用 $J$，它们也不会发生任何杂乱的自旋翻转，系统表现为纯粹的、完美的 Ising 相互作用系统。仅需令该纵向耦合脉冲作用一段持续时间：

体系即可自然积累得到高保真度的受控 $\pi$ 相位移动，从而在不需要施加超大塞曼失谐、极大缓解微波控制线负荷的背景下，优雅、纯净地实现保真度逼近 $99.99\%$ 的两比特受控相位门 $CZ$。

本项工作所构建的这一套高效、高物理保真度的“Schrödinger-Poisson + CI”数值仿真框架，不仅能够准确预测并指导这一类物理黄金“甜点”的排布，也必将成为半导体量子芯片设计（CAD for Spin Qubits）和量子计算机辅助制造（TCAD）流程中必不可少、高频部署的核心软件算法支柱。