来源论文: https://arxiv.org/abs/2607.04876v1 生成时间: Jul 07, 2026 11:21

0. 执行摘要

半导体硅基自旋量子比特（Silicon Spin Qubits）由于其与现代互补金属氧化物半导体（CMOS）工艺的高度兼容性，被视为实现大规模集成量子芯片（Up to Millions of Qubits）最具前景的物理路线之一。然而，在纳米级 CMOS 芯片中，紧凑且高度不对称的栅极堆叠（Gate Stacks）结构给器件建模带来了极大的挑战。设计技术协同优化（Design-Technology Co-Optimization, DTCO）的核心诉求在于：如何建立一个透明、高效且可审计的物理模型，将晶圆厂的几何版图（Layout）与底层的量子多体自旋填充和量子相干特性直接联系起来？

传统的器件建模存在严重的“两阶段断裂”：

  1. 静电学阶段：依赖于重型、黑盒化的技术计算机辅助设计（TCAD）有限体积法（FV）求解器，计算资源消耗极大，且难以直观解析几何参数对势场的物理制约关系。
  2. 量子力学阶段：直接抛弃真实的静电势场，代之以人为设定的简谐势（Harmonic Potential）或简单的拟合势，从而割裂了“工艺制造-多体自旋”的因果链条。

本工作开发了一种统一的静电到自旋（Electrostatic-to-Spin）降低阶解析计算框架，填补了这一空白。其核心静电引擎——泊松核耦合界面格林函数（Poisson-kernel Coupled-interface Green-function, PK-GF）模型，在不对任何数值解进行参数拟合的前提下，实现了与独立有限体积法（FV）解在毫伏（millivolt）尺度上的高度一致性。在此基础上，该框架将真实的栅极约束势直接输入自旋-谷多体计算模块，对含有 $N = 2-17$ 个电子的延展型“软糖豆”（Jellybean）量子点进行了模拟。计算表明，不设限哈特里-福克（UHF）方法能精确捕捉到由于强库仑排斥引起的类 Wigner 分子电荷局域化，但会系统性高估自旋极化；而完全活性空间组态相互作用（CASCI）计算则支持了与实验观测相符的低自旋分支。这一端到端的工程链条为工业级 CMOS qubit 芯片的设计与 DTCO 提供了强大的理论工具。

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

1.1 核心科学问题

在硅自旋量子点芯片中，微型栅极不仅用于界定电荷的物理边界，还决定了隧道势垒高度、轨道能级差、电荷占有率以及自旋壳层填充规律。在非对称多栅结构中，栅极宽度不等、分布于不同的介质界面上，且通过复杂的连续介电层（如 $\text{SiO}_2$、$\text{Al}_2\text{O}_3$、空气等）相互作用。这导致量子点内部的约束势（Confinement Potential）表现出极强的不对称性和空间不均匀性。如何建立一个能够自洽处理多导体电荷再分配、多介质界面极化反射，同时在数学上解析可解的静电模型，是当前该领域的核心科学问题。

1.2 理论基础与公式推导

PK-GF 模型从第一原理出发，通过多介质层格林算符和边界积分方法，构建了非对称栅极系统的静电场解析解。

1.2.1 基础传播子与泊松核（Poisson Kernel, PK）

对于厚度为 $H$、底部接地的硅衬底（介电常数为 $\epsilon_{\text{Si}}$），在 $y < 0$ 的半无边界域中，由单个有限宽度栅极 $j$（区间为 $[L_j, R_j]$）贡献的二维半平面响应可表示为泊松核形式。对于 $y < 0$，该有限区间响应函数 $A_j(x,y)$ 定义为：

$$A_j(x, y) = \frac{1}{\pi} \left[ \arctan\left(\frac{R_j - x}{-y}\right) - \arctan\left(\frac{L_j - x}{-y}\right) \right]$$

在傅里叶空间中，有限厚度衬底的接地边界条件通过衬底传播子（Substrate Propagator）$G_H(k,y)$ 予以修正：

$$G_H(k, y) = \frac{\sinh[|k|(y + H)]}{\sinh(|k|H)}, \quad G_H(0, y) = \frac{y+H}{H}$$

基础 PK 模型通过对 4 个孤立栅极在平面对称多层介质中的多重反射进行算符传播，但忽略了未覆金属区域的电极再分配效应。其一阶近似电势为：

$$\phi_{\text{PK}}(x, y) = \mathcal{F}^{-1} \left\{ G_H(k, y) \sum_{j=1}^4 V_j \hat{\chi}_j(k) T_j(k) \right\}$$

其中，$\hat{\chi}_j(k)$ 是栅极空间窗口函数的傅里叶变换，$T_j(k)$ 是通过传输矩阵方法算得的多层介质传输因子。

1.2.2 耦合界面格林函数（PK-GF）拓展

基础 PK 模型最大的物理局限性在于其**“零连续假设”（Zero-continuation Assumption）**：它假设在栅极所在的有效平面上，除了金属栅极之外，其余未覆盖的介质区域电势均为零。这显然违背了连续介质极化的物理规律。为了修正这一偏差，PK-GF 模型引入了耦合界面格林算符。

设定系统的 4 个关键电学界面为：硅/二氧化硅界面、下层氧化铝界面、栅极平面及上层介质界面。记其界面势向量为 $\mathbf{u}(k) = [u_0, u_6, u_{10}, u_{14}]^{\top}$（数字对应 $y$ 轴坐标）。通过齐次层消元法（Homogeneous-layer Elimination），可将多层介质中的一维泊松方程转化为四界面狄利克雷-诺伊曼（Dirichlet-to-Neumann）矩阵 $\mathbf{D}(k)$：

$$\mathbf{q}(k) = \mathbf{D}(k)\mathbf{u}(k) \implies \mathbf{u}(k) = \mathbf{D}^{-1}(k)\mathbf{q}(k) = \mathbf{G}(k)\mathbf{q}(k)$$

这里，$\mathbf{G}(k)$ 即为多界面格林算符。$\mathbf{D}(k)$ 的矩阵元素由各介质层的电学参数和空间波速确定：

$$A_{\ell}(k) = \epsilon_{\ell} |k| \coth(|k| d_{\ell}), \quad B_{\ell}(k) = \epsilon_{\ell} |k| \text{csch}(|k| d_{\ell})$$

其中 $\epsilon_{\ell}$ 和 $d_{\ell}$ 分别为第 $\ell$ 层的介电常数和厚度。

1.2.3 边界奇异性与切比雪夫电荷展开

导体边缘存在众所周知的电荷累积和静电场奇异性（Edge Singularity），其电荷密度在边界处呈 $1/\sqrt{1-\xi^2}$ 散逸。为了确保高精度求解，PK-GF 在金属电极上使用了**边缘感知的切比雪夫基底函数（Edge-aware Chebyshev Basis Functions）**来展开未知电荷密度 $\sigma_j(x)$：

$$\sigma_j(x) \approx \sigma_j^{\text{eff}, (M_q)}(x) = \frac{\chi_j(x)}{\sqrt{1 - \xi_j^2(x)}} \sum_{m=0}^{M_q-1} a_{jm} T_m[\xi_j(x)]$$

其中，$\xi_j(x) = \frac{2(x - x_{c, j})}{W_j}$ 为映射至 $[-1, 1]$ 区间内的无量纲局部坐标，$W_j$ 为栅极宽度，$x_{c,j}$ 为其中心坐标。$T_m(\xi)$ 为第一类切比雪夫多项式。

通过对格林算符施加这些基底函数，并在空间中配置（Collocate）给定的栅极偏置电压，可获得紧凑的线性代数方程组：

$$\mathbf{K}\mathbf{a} = \mathbf{V}$$

其中 $\mathbf{a} = \{a_{jm}\}$ 为物理电荷展开系数，$\mathbf{V}$ 为给定的栅极电压向量。求解该矩阵即可得到全局势场的解析表示，其形式高度闭合，不再依赖繁琐的数值网络求解。

+--------------------------------------------------------------+
|               PK-GF 降低阶静电学建模工作流                      |
+--------------------------------------------------------------+
                               |
      [ 输入: 栅极位置 L_j, R_j, 介质常数 \epsilon_l, 偏置 V_j ]
                               |
                               v
              +---------------------------------+
              |  计算傅里叶空间 Dirichlet-to-   |
              |  Neumann 耦合矩阵 D(k) 和 G(k)  |
              +---------------------------------+
                               |
                               v
              +---------------------------------+
              | 引入边缘奇异性切比雪夫电荷展开   |
              | \sigma_j(x) ∝ T_m(\xi)/√(1-\xi^2)|
              +---------------------------------+
                               |
                               v
              +---------------------------------+
              | 求解非奇异配置矩阵方程 K a = V   |
              |      (SVD 最小二乘法求解)        |
              +---------------------------------+
                               |
                               v
              +---------------------------------+
              | 执行快速傅里叶逆变换 (IFFT)     |
              | 得到 y = -2 nm 处的 1D 约束电势   |
              +---------------------------------+
                               |
                               v
+--------------------------------------------------------------+
|                    对接多体自旋填充计算                       | 
+--------------------------------------------------------------+

1.3 技术难点与物理瓶颈

  1. 多介质界面反射的非指数衰减：在传统的等效氧化层厚度（EOT）近似中，多层介质通常被简化为单一呈指数衰减的衰减常数。然而，本工作发现，在高空间频率 $|k|$ 区域，多层介质带来的界面反射和高频滤波效应无法被单一特征长度捕获，导致常规方法在通道深度处（$y = -2\text{ nm}$）产生巨大误差。
  2. 导体电荷的协同再分配：当给某一栅极施加偏置时，周围未偏置的金属电极以及裸露的介质界面上电荷会发生重分布，形成全局非定域极化。PK-GF 模型必须自洽地在算符层面解耦各电极之间的互感电荷，这对非对称三维堆叠器件而言在数学上极具挑战。
  3. 多体理论中的自旋极化过高问题（Spin Over-polarization）：在量子化学层次，采用平均场近似（如不设限哈特里-福克 UHF）处理延展量子点时，由于电子数较多（$N \ge 12$），其断裂对称性的行列式波函数会严重混入高自旋态（Spin Contamination），从而给出过于极化的磁性基态，这与实验观测不符。

2. 关键 Benchmark 体系、计算数据与性能展示

为了全面验证该方法的精度和计算性能，研究团队定义了一个极具代表性的非对称四栅极硅基量子器件（基于真实的 Jellybean 耦合器几何设计），并开展了精细的数值基准测试。

2.1 测试体系配置

  • 几何结构：
    • 硅衬底：厚度 $H = 2000\text{ nm}$，相对介电常数 $\epsilon_{\text{Si}} = 11.7$。
    • 连续氧化层：$6\text{ nm}$ $\text{SiO}_2$（$\epsilon = 3.9$）。
    • 下层介质层：$4\text{ nm}$ $\text{Al}_2\text{O}_3$（$\epsilon = 9.0$）。
    • 柱塞栅极（Plunger Gates）：G2 和 G4 置于 $y = 10\text{ nm}$ 平面。
    • 栅间绝缘层：$4\text{ nm}$ $\text{Al}_2\text{O}_3$。
    • 交换栅极（Exchange Gates）：G1 和 G3 置于 $y = 14\text{ nm}$ 平面。
    • 栅极覆盖区间：$G_1 = [0, 80)\text{ nm}$, $G_2 = [83, 233)\text{ nm}$, $G_3 = [236, 276)\text{ nm}$, $G_4 = [279, 319)\text{ nm}$。这一分布呈现典型的不对称和间隙分布特征。
  • 验证电压偏置：$V_1 = 0.80\text{ V}$, $V_2 = 0.35\text{ V}$, $V_3 = 1.05\text{ V}$, $V_4 = 0.55\text{ V}$（非对称偏置，最大化考量静电耦合的复杂性）。

2.2 性能对比：PK vs. PK-prior vs. PK-GF vs. 数值 FV 法

作为参照的数值模型是利用成熟半导体工业 TCAD 软件构建的经典五点有限体积法（Finite Volume, FV）模型。其计算网格划分极其致密：在整个验证区内利用 $dx = dy = 1\text{ nm}$ 的超细网格进行离散，共有大约 649 万个节点。

测试在 $x \in [-100, 420]\text{ nm}$，$y \in [-80, -3]\text{ nm}$ 这一宽广的二维区域（覆盖了量子点形成和输运的核心物理通道）进行，其势场均方根误差（RMSE）统计结果如下：

静电模型方案势场均方根误差 (RMSE)最大绝对误差 (Max Error)物理机制特征
基础泊松核模型 (PK)$239.453\text{ mV}$$> 400\text{ mV}$采用零连续性假设，严重低估未覆盖区域势场
先验知识模型 (PK-prior)$46.017\text{ mV}$$\sim 120\text{ mV}$引入几何背景权重，但未求解导体电荷再分配
本工作模型 (PK-GF)$1.134\text{ mV}$$2.699\text{ mV}$完美考虑多界面耦合与切比雪夫电荷累积，无需拟合

在最为严苛的量子点通道深度（$y = -2\text{ nm}$ 截线）处，不同模型对真实有限体积物理场的复现能力对比如下：

  • PK 模型 RMSE：$264.797\text{ mV}$
  • PK-prior 模型 RMSE：$62.284\text{ mV}$
  • PK-GF 模型 RMSE：$1.283\text{ mV}$

电场强度 RMSE 验证：

  • PK-GF 的电场均方根误差仅为 $5.68 \times 10^{-6}\text{ V/nm}$，而传统 PK 模型的误差则高达 $1.46 \times 10^{-3}\text{ V/nm}$。这表明在力场精度上，PK-GF 相比前代方案实现了近三个数量级的跨越。

2.3 模式数收敛与稳定性分析

为了定量探究在切比雪夫多项式中展开项数 $M_q$（每个栅极的模式数）对最终势场精度的边际贡献，论文系统扫描了不同的基底维度：

  • $M_q = 10 \to 12$：势场平均变化仅为 $0.081\text{ mV}$ RMS；
  • $M_q = 12 \to 14$：变化量进一步收敛到 $0.034\text{ mV}$；
  • 条件数（Condition Number）随基底扩大的变化：
    • 在 $M_q = 12$ 时，矩阵 $\mathbf{K}$ 的条件数为 $5.35 \times 10^3$；
    • 在 $M_q = 14$ 时，条件数上升到 $1.24 \times 10^4$。

这组数据清晰地揭示了量子比特设计在计算效率与数值条件稳定性之间的“工程拐点”（Engineering Knee）：每个栅极设置 12 个切比雪夫模式（整个系统共 48 个求解未知数），即可在毫伏精度下实现系统的绝对收敛，且配置矩阵保持了极好的稳定性。 相比求解数百万个节点的传统 FV 有限体积矩阵，PK-GF 成功将计算自由度降低了 5 个数量级！

3. 代码实现细节与复现指南

本章节旨在为量子化学计算及器件物理方向的科研人员提供一套可以直接部署和编程复现的实用指南。

3.1 静电核心 PK-GF 的数学算法流程

要复现 PK-GF，首先需完成以下五个核心计算步骤：

  1. 傅里叶域狄利克雷-诺伊曼矩阵 $\mathbf{D}(k)$ 的构建： 针对任意给定层厚 $d_\ell$ 和介电常数 $\epsilon_\ell$，预先在空间频率 $k$ 轴（网格范围取 $k \in [0, k_{\text{max}}]$）上构建 $4 \times 4$ 传输阻抗矩阵。

  2. 切比雪夫空间到光谱域的变换： 对于任意栅极 $j$ 的第 $m$ 阶模式，其空间基底 $\theta_{jm}(x) = T_m[\xi_j(x)] / \sqrt{1 - \xi_j^2(x)}$ 可通过解析公式进行一维快速傅里叶变换（FFT），得到其在 $k$ 空间内的表示 $\hat{\theta}_{jm}(k)$。由于切比雪夫多项式与贝塞尔函数（Bessel Functions）的内生对应，其解析变换公式为：

    $$\hat{\theta}_{jm}(k) = \pi (-i)^m J_m\left(\frac{W_j k}{2}\right) e^{-i k x_{c, j}}$$

    其中 $J_m$ 是 $m$ 阶第一类贝塞尔函数。这一解析性质极大提升了代码执行速度，规避了数值积分误差。

  3. 配置矩阵 $\mathbf{K}$ 的组装： 通过在 4 个栅极的金属表面均匀布设 $M_{\text{coll}} = M_q$ 个实空间探测点，将格林函数卷积积分转化为谱空间求和。配置矩阵中的第 $(p, q)$ 个元素表示第 $q$ 个电荷模式在第 $p$ 个空间探测点所产生的电势：

    $$K_{p, q} = \mathcal{F}^{-1} \left\{ G_H(k, y_{\text{probe}}) \mathbf{e}_{\text{probe}}^{\top} \mathbf{G}(k) \mathbf{e}_{\text{source}} \hat{\theta}_{jm}(k) \right\}_{x = x_{\text{probe}}}$$

  4. 代数系统奇异值分解（SVD）求解： 由于多重电极间的强静电屏蔽，配置矩阵 $\mathbf{K}$ 在模式数较高时可能出现轻微的病态（Ill-conditioned）。复现代码中应统一采用 SVD 最小二乘算法（如 Python 中的 scipy.linalg.lstsq 或 Julia 的 SVD 求解器），以保证数值求解的鲁棒性。

  5. 计算一维势场切片： 提取计算通道线 $y = -2\text{ nm}$ 处的势场分布 $\phi_{\text{PK-GF}}(x, y_0)$，并乘以电荷转换系数 $\alpha_{\text{lever}} = 0.052\text{ eV/V}$，得到单电子势能：

    $$U_x(x) = -\alpha_{\text{lever}} \phi_{\text{PK-GF}}(x, y_0)$$

3.2 量子多体计算流（UHF 与 CASCI）

量子物理计算部分的核心在于将前面得到的真实物理电势 $U_x(x)$ 代入自旋-谷多体哈密顿量中求解：

3.2.1 空间谷-轨道自洽求解

硅具有两个沿输运轴投影的低能传导带谷（Valley $v_{\pm}$）。在势能 $U_x(x)$ 和横向限域势（由 $z$ 方向的约束孔径定义，一般取 $z_{\text{aperture}} = \pm 20\text{ nm}$）共同作用下，对二维单电子有效质量方程进行离散求解：

$$\left[ -\frac{\hbar^2}{2 m_x^*} \frac{\partial^2}{\partial x^2} - \frac{\hbar^2}{2 m_z^*} \frac{\partial^2}{\partial z^2} + U_{xz}(x, z) \pm \frac{\Delta_{\text{valley}}}{2} \right] \psi_i(x, z) = \epsilon_i \psi_i(x, z)$$

其中 $m_x^* = 0.916 m_0$，$m_z^* = 0.190 m_0$，谷分裂能级设定为经典的 $\Delta_{\text{valley}} = 0.12\text{ m{e}V}$。空间哈密顿量在实空间格点上通过切比雪夫滤波子空间迭代法（Chebyshev-filtered Subspace Iteration）快速对角化，提取出前 45 个最低的谷-轨道空间态作为后续多体计算的活性轨道（Active Orbitals）。

3.2.2 不设限哈特里-福克（UHF）自洽场与 DIIS

在单粒子空间谷-轨道基底（45 个空间函数 $\times$ 2 个 spin 自由度 $=$ 90 个自旋轨道）上，组装两体库仑矩阵元素 $V_{ijkl}$ 和交换项 $K_{ijkl}$。多体计算的单粒子自洽 Fock 算符可紧凑写为：

$$\hat{H}_s [\mathbf{n}] = \hat{H}_{\text{kin}} + U_{\text{PK-GF}}(x,z) + \hat{H}_{\text{valley}} + V_H[\mathbf{n}] + V_x^{(s)}[\mathbf{n}] - \frac{s}{2} g \mu_B B$$

其中 $s = \pm 1$。使用**直接反转子空间迭代（DIIS）**方法加速 UHF 自洽场计算。收敛判定标准极其严苛：单步能量波动必须小于 $10^{-4}\text{ m{e}V}$，且电荷密度矩阵和Fock算符的残差模长均须在 $10^{-4}$ 以下。

3.2.3 CASCI 多体组态相互作用

为纠正 UHF 所引入的自旋对称性破缺（Spin Contamination），本框架开发了基于**非冷冻芯（No-frozen-core）**的完全活性空间组态相互作用（CASCI）对角化算法。

  • 对特定的电子数 $N$（如 $N=2, 4, 6, 8, 12$），选取费米面附近的 $n_{\text{act}} = 12$ 个单粒子空间能级（包含 valley 轨道）作为活性空间。
  • 构造在总自旋 $S^2$ 及磁分量 $S_z$ 下具有严格对称性的量子化学 Slater 行列式态。
  • 在整个活性空间中构建精确的多体哈密顿矩阵，并通过 Lanczos 稀疏矩阵求解器直接计算低自旋（Singlet, $S=0$）与高自旋（Triplet 或更高）状态的激发谱，从而确认物理真实的基态自旋属性。

3.3 开源推荐与可复现接口

虽然该论文中的 PK-GF 模型为研究团队的自研代码，但量子化学计算模块（UHF 和 CASCI）可与现有的开源高性能科学计算工具包完美融合。推荐读者采用以下开源接口复现该框架：

  1. 静电学模块 (PK-GF)： 可基于 Julia 的核心线性代数库 LinearAlgebra.jl 以及快速傅里叶变换包 FFTW.jl 编写。用不超过 300 行 Julia 代码即可实现极速的四界面 $\mathbf{K}$ 配置矩阵组装与 SVD 求解。
  2. 电子结构与多体量子化学计算 (UHF/CASCI)： 推荐使用 PySCF (Python Chemistry Framework) 软件。PySCF 提供了高度优化的 UHF 求解器、DIIS 收敛器以及完整的 CASCI 接口。具体复现流程为：
    • 通过 PK-GF 得到一维静电势，在实空间数值格点上计算前 45 个单粒子波函数；
    • 计算这些波函数之间的两电子积分（ERI, Electron Repulsion Integrals），将其打包写入 NumPy .npy 文件中；
    • 将该两电子积分矩阵及单粒子核心 Hamiltonian 矩阵通过 pyscf.gto.mole 模块导入 PySCF；
    • 初始化 pyscf.scf.UHF 求解自洽场，随后使用 pyscf.mcscf.CASCI 定义活性空间 $n_{\text{act}} = 12$。这样能无缝复现论文中如图 4h 和图 5 所示的所有物理曲线。

4. 关键引用文献与局限性评论

4.1 关键引用文献

本项工作在理论上根植于以下经典文献：

  1. Davies et al., J. Appl. Phys. 77, 4504 (1995): 首次建立了经典 Patterned 二维电子气电极的静电泊松核计算模型。本工作在数学上是对其一维电极边界条件的重大突破，成功将其推广到了多界面极化反射的全新体系。
  2. Loss & DiVincenzo, Phys. Rev. A 57, 120 (1998): 提出了利用量子点中自旋态进行量子计算的奠基性设想。本工作开发的静电-自旋闭环正是实现这一设想的工程利器。
  3. Zwerver et al., Nat. Electron. 5, 184 (2022): 报道了采用先进工业半导体产线制造的自旋量子比特芯片。该文献指明了半导体 qubit 进行大规模 DTCO 的必要性，也是本工作静电叠层参数的物理原型基础。
  4. Wang et al., Adv. Mater. 35, 2208557 (2023): 首次探索了硅基“软糖豆”（Jellybean）量子点在长程自旋耦合及车载模拟中的应用。本论文多体计算部分直接继承并深化了该工作的耦合物理学。

4.2 工作局限性深度评论

尽管 PK-GF 展现出了极高的物理保真度和令人印象深刻的毫伏级精度，作为一名资深的技术作者与审稿人，我必须指出该工作在向实用化量子芯片 CAD 软件演进过程中存在的几项物理局限性：

1. 二维静电场（2D Electrostatics）局限

目前的 PK-GF 模型本质上是一个二维垂直截面模型。它默认栅极在垂直于截面的第三轴方向（$z$ 轴）具有平移不变性。然而，现代前沿 CMOS 比特（如 FinFET、GAA 纳米线比特等）具有极其强烈的三维空间特征。三维情况下的边缘电荷积累不仅发生在栅极的侧边缘，还发生在栅端角部。这种三维边缘奇异性会严重改变三维量子点的真实限域高度。2D 简化使得该模型目前无法处理三维角落效应（Corner Effects）以及复杂的电极间串扰问题。

2. 线性静电局限，缺失自洽载流子屏蔽效应

PK-GF 是一个纯线性经典泊松方程求解器（$\nabla \cdot [\epsilon(x,y)\nabla \phi] = 0$）。它在求解静电势时，假设硅衬底内完全没有移动电荷（Mobile Charge）。这在器件偏置极低、电子数极少（如 $N < 2$）时是合理的。然而，随着注入的电子数增多（本工作模拟到了 $N=17$），通道底部的反型层会积累高密度的电子气。这些移动电荷会对栅极势场产生强烈的量子屏蔽（Quantum Screening）。如果不引入自洽的 Poisson-Schrödinger 循环求解，纯电势场在电子数较多时必然会高估约束势深度，从而影响多体积分的准确度。

3. 简化且理想化的物理界面

本模型的数学基石建立在边界条件完全规则的平面异质结层上。但在实际工艺中：

  • 金属电极具有非零的物理厚度，而非本模型中设定的“零厚度薄片”；
  • 极复杂的二氧化硅/硅界面处存在高密度的电荷陷阱、无序缺陷（Disorder）以及表面粗糙度（Surface Roughness）。 缺陷产生的随机电势涨落通常可以达到 $10-50\text{ m{e}V}$ 量级，这远大于 PK-GF 模型努力消除的几毫伏静电误差。在未来的模型版本中，必须引入带有噪声和无序扰动算符的随机介电层项，才有可能对真实的成品器件进行精准预测。

4. 自旋-谷相互作用的唯象学近似

在量子多体计算中，磁场仅通过 Zeeman 劈裂项引入，忽略了高磁场下（$B=5\text{ T}$）轨道受到的洛伦兹力偏转效应。此外，谷轨道相互作用被简化为了常数。真实的硅自旋比特中，由于界面尖锐度不足，谷能级极易与空间轨道能级发生剧烈杂化。这会导致计算给出的自旋-谷壳层排布规律（Shell Filling）在特定磁场下产生系统性偏差。

5. 补充技术洞察：Jellybean 耦合器中的物理学与 DTCO 应用指南

为了让科研人员在实际 DTCO 流程中最大化利用这套框架，本节将深入探讨“软糖豆”量子点独特的物理学特征，并提供一份系统级优化指南。

5.1 类 Wigner 分子局域化的物理实质

在常规的极小、紧凑量子点中，电子空间分布呈圆形，且壳层填充严格遵循类似于原子核外电子排布的“费米液体”图像（Fermi Liquid Picture）。但是在长达 $150-300\text{ nm}$ 的延展型 Jellybean 量子点中，由于纵向一维约束极弱，电子间的库仑斥力（静电排斥）开始压倒单粒子轨道能级差。

本工作的多体模拟清晰地展示了这一竞争过程：

  • 相互作用能与能级差之比 $R$：

    $$R = \frac{E_C}{\Delta}$$

    在测试中，这一比值始终处于 $2.5$ 到 $4.4$ 之间。这表明系统进入了典型的中度到强关联区域。

  • 多峰局域化现象： 如图 5 所示，随着电子数 $N$ 从 2 逐步增加到 17，一维电子电荷密度曲线不再是一个平滑的、展宽的包络面，而是分裂成了若干个分立的空间峰（Lobes）。至 $N=17$ 时，系统已经自发形成了 8 个高度局域化的电荷簇。这种电子在静电斥力驱动下发生晶格化排布的现象，正是经典的**“类 Wigner 分子”（Wigner-molecule-like）**状态。

  • 自旋极化的物理根源： 在 UHF 平均场层次下，电荷局域化通常伴随着自旋极化的剧烈抬升（图 4f 所示的 $S_z$ 呈现极高不规则极化状态，最高达 8.5）。这一高自旋行为是“Stoner 转变”的一种非局域物理呈现，因为局域化的轨道重叠度降低，体系自发倾向于自旋铁磁对齐来规避两体排斥。然而，CASCI 相关计算表明，这些局域化电子之间由于超交换（Superexchange）作用，在没有外界轨道各向异性破缺时，真实的基态自旋会高度倾向于“自旋安静”（Spin-quiet）的低自旋分支（图 5r，$S=0$ 保持在每一活性空间）。这就要求器件设计师在设计比特输运通道时，必须充分考虑这一多体量子相关纠缠效应，不能直接采用 UHF 进行简单预测。

+-------------------------------------------------------------+
|                        DTCO 优化闭环                        |
+-------------------------------------------------------------+
                               |
          [ 工艺选择 / 几何参数版图 (Layout, Stack) ]
                               |
                               v
              +---------------------------------+
              |      通过 PK-GF 模拟势场限域    |
              |    (快速扫描栅极间距、介质层厚度) |
              +---------------------------------+
                               |
                               v
              +---------------------------------+
              | 评估 Wigner 局域化指数 R = Ec/Δ  |
              |   (确保通道在注入电子时维持相干) |
              +---------------------------------+
                               |
                               v
              +---------------------------------+
              | 评估低自旋分支激发能级 (CASCI)  |
              | (使 singlet-triplet 能量差最大化)|
              +---------------------------------+
                               |
                               v
          [ 优化输出: 提高量子比特读取纯度与调控保真度 ]

5.2 面向芯片设计者的 DTCO 实践指南

半导体量子比特制造的核心瓶颈在于极低的成品率和器件间性能的剧烈漂移。借助 PK-GF，版图设计师和工艺工程师可以无需经过繁重的工艺流片和漫长的 TCAD 计算，直接实施以下优化设计：

  1. 设计规则审查（Design Rule Manual, DRM）制定： 通过快速扫描柱塞栅与交换栅的间距（例如从 $1\text{ nm}$ 变化到 $10\text{ nm}$），可以立刻评估其对电子势垒高度的影响。这有助于确定工艺中允许的最小光刻套准偏差（Overlay Budget）。
  2. 叠层介质工程（Dielectric Stack Engineering）： PK-GF 模型可以极其方便地修改各氧化层和阻挡层的厚度和介电常数（如扫描 $\text{SiO}_2$ 与 $\text{Al}_2\text{O}_3$ 的物理配比）。设计目标在于：在不改变光刻版图的情况下，通过调整叠层成分组合，使量子点处的横向静电约束电场（Electric Field Slope）达到最陡峭状态，从而最大限度地压制谷-轨道杂化，提高量子调控速度。
  3. 利用电荷局域化作为设计信号（Design Signal）： 在设计用于连接两个远距离处理器单元的 Jellybean 长程耦合器（Quantum Coherent Link）时，核心诉求是实现“自旋相干的长程无损传输”。如果 Wigner 局域化过强（$R > 4.5$），电子将彻底钉扎在固定的势能格子上，转变为类似 Hubbard 链的多体强局域状态，此时电荷输运效率将呈现指数衰减；若 $R$ 过小，则无法实现电子的单发控制。利用本框架计算 $R$ 指数对偏置电压和电学结构的响应，可以协助建立具有最佳自旋中介（Spin-quiet Mediator）传输能力的自洽偏置工作窗口。

总而言之，本工作不仅是一项优雅的解析静电学物理突破，更搭建了首条连接 “半导体物理制造” 与 “量子物理 observables” 的数字化桥梁。对于任何致力于半导体硅自旋量子芯片商业化进程的团队来说，本工作所蕴含的理论框架都是极具借鉴意义的重要工程指引。