🔬 mesoscale physics

Quasiparticle-resolved variational theory of Andreev spin qubits

本文开发了一种针对安德烈耶夫自旋量子比特（Andreev spin qubits）的连续变分理论，该理论揭示了相干双准粒子构型如何主导自旋依赖的约瑟夫森耦合，并提供了关键参数的闭式表达式，从而在强关联量子点模型与有效电路哈密顿量之间搭建了桥梁。

原作者： Teodor Iličin, Rok Žitko

发布于 2026-10-07

📖 1 分钟阅读☕ 轻松阅读

CC BY 4.0

原作者： Teodor Iličin, Rok Žitko

原始论文采用 CC BY 4.0 许可（http://creativecommons.org/licenses/by/4.0/）。 ✨ 这是对下方论文的AI生成解释。它不是由作者撰写或认可的。如需技术准确性，请参阅原始论文。阅读完整免责声明

在量子计算的微观世界中，信息存储在最脆弱的状态之中：单个电子的自旋。这种自旋就像一根微小的指南针，通过指向“上”或“下”来代表零或一。为了保护这些脆弱的信息免受嘈杂环境的影响，科学家们通常将电子捕获在超导电路中，这是一种电阻为零的导电材料。当电子被限制在这种电路中时，它会形成一种特殊的束缚态，被称为安德烈耶夫束缚态（Andreev bound state）。这些状态是一种新型量子比特——即安德烈耶夫自旋比特（Andreev spin qubit）的核心。研究人员面临的挑战是，准确理解这个被捕获的电子是如何与其周围的超导电子海进行相互作用的。虽然电子承载着信息，但它并非孤立存在；它在不断地与超导体交换能量和虚粒子，这种复杂的量子力学之舞决定了量子比特的行为方式以及如何对其进行控制。

来自斯洛文尼亚约热夫·斯特凡研究所（Jožef Stefan Institute）的一个物理学家团队开发了一种能够高精度映射这种相互作用的新方法。他们专注于一个特定的装置，其中一个量子点（充当电子的微小岛屿）连接着两个超导引线。在这个系统中，岛屿上的电子可以通过两条不同的路径与超导引线发生相互作用：一条是通过岛屿本身的直接路径，另一条是绕过岛屿的背景路径。研究人员希望了解这些相互作用如何产生量子比特自旋向上和自旋向下状态之间的能量差，这一特性对于读取和写入量子信息至关重要。为此，他们建立了一个详细的数学模型，用于追踪电子及其周围“准粒子”（在超导体中表现得像粒子的激发态）的行为。

其工作的核心涉及一种称为变分理论（variational theory）的方法，这本质上是一种通过猜测量子波函数最可能的形状，然后不断优化该猜测，使其尽可能符合物理定律的方法。在这种情况下，研究人员构建的模型包含了系统在任何给定时刻可能拥有零个、一个或两个超导引线中准粒子的可能性。虽然发现两个准粒子的概率非常小，但团队发现这些罕见的构型却出人意料地重要。尽管它们出现的频率很低，但它们以一种显著改变量子比特能量的方式干扰了更常见的零准粒子状态。具体而言，即使在电子间排斥力较弱的情况下，这两个准粒子构型也贡献了导致自旋向上和自旋向下状态产生能量差异的主要效应的一半。随着电子间排斥力的增强，这些两准粒子构型成为了决定量子比特行为的主导因素。

研究人员使用被称为数值重整化群（numerical renormalization group）计算的高性能计算机模拟对他们的模型进行了测试，这被认为是解决此类问题的金标准。他们的结果与模拟结果高度吻合，证实了他们简化的模型捕捉到了本质的物理特性。他们发现，量子比特与超导体之间的相互作用强度很大程度上取决于库仑排斥力（即推开电子的力量）。当这种排斥力很强时，驱动超电流的传统能量比自旋相关能量下降得更快。这意味着在具有强电子排斥力的系统中，自旋相关效应变得相对更加重要，尽管整体与超导体的耦合度减弱了。这一发现解释了为什么某些实验装置在具有强电子相互作用的情况下，仍能显示出强烈的自旋相关信号。

除了计算能量之外，团队还绘制出了这些虚粒子在空间中的分布图。他们发现，这些粒子并不局限于微小的量子点，而是分布在与超导体相干长度（衡量超导序延伸距离的度量）相当的距离内。此外，他们还发现电子的自旋并非完全局限在量子点上；一部分自旋与超导引线共享。这种自旋共享受到一种被称为自旋轨道耦合（将电子自旋与其运动联系起来）的性质的影响。当施加磁场时，这种自旋共享会产生一种独特的效果：即使在接头两端没有电压差的情况下，也会产生电流。这种电流对于自旋向上和自旋向下的状态都是相同的，这意味着它不依赖于量子比特的信息，而是取决于自旋本身的物理重新分布。

这项研究在复杂的、相互作用电子的微观现实与用于设计量子电路的更简单的有效模型之间架起了一座清晰的桥梁。通过展示波函数结构如何决定耦合关系，研究人员为工程师提供了一种预测这些量子比特在实际器件中如何表现的工具。他们证明了，如果忽略掉那些罕见的两准粒子事件，将会导致对量子比特性能预测出现显著误差。这项工作证实，要真正理解并控制这些量子比特，必须考虑被捕获的电子与超导体中虚激发之间微妙且相干的相互作用，这种关系比此前假设的要复杂得多，也深远得多。

技术摘要：安德烈耶夫自旋比特的准粒子分辨变分理论

问题陈述
安德烈耶夫自旋比特（ASQs）将量子信息编码在受限于超导点接触或混合异质结构中安德烈耶夫束缚态（Andreev bound state）内的单个准粒子（QP）自旋中。虽然唯象哈密顿量在许多机制下都能成功描述 ASQs，但对于自旋依赖约瑟夫森耦合（ $E_S$ ）以及虚准粒子（virtual quasiparticle, QP）构型的作用，微观理解仍不完整。具体而言，涉及库仑相互作用、自旋轨道耦合（SOC）与超导连续态之间的相互作用，需要一种超越简单有效模型或低阶微扰论的描述。开放性问题包括：携带自旋的束缚态的空间延展性、部分近藤屏蔽（Kondo screening）的起始点，以及在存在连续态的情况下，自旋分裂的精确微观起源。

方法论
作者为由两个超导引线耦合的安德森杂质模型（AIM）所构成的量子点约瑟夫森结中的奇宇称双重态（odd-parity doublet）开发了一种连续变分描述。该系统包含：

哈密顿量： 一个相互作用的杂质轨道（库仑排斥 $U$ ，失谐 $\delta$ ），通过自旋守恒的跳跃（ $V$ ）以及包含自旋依赖振幅和 SOC 的直接引线间隧穿项（ $W_\sigma$ ）与 BCS 引线耦合。
变分波函数（Variational Ansatz）： 波函数构建为包含最多两个 Bogoliubov 准粒子的构型叠加。这包括真空态（0-QP）、单准粒子态（1-QP）以及进一步分解为单态和三线态浴对（singlet and triplet bath pairs）的双准粒子态（2-QP）。
对称性与通道： 作者利用幺正对称性（ $\eta$ ）并引入对称/反对称准粒子通道（ $\gamma_{s,k,\sigma}$ ）以简化变分方程。
近似方案： 通过三个层级的核截断（Kernel truncation）来求解导出的耦合积分方程（第 1 层：仅 1-QP；第 2 层：无反馈的 2-QP；第 3 层：带有混合核反馈的 2-QP）。
验证： 变分结果通过使用开源 qdjj solver 与数值重整化群（NRG）和密度矩阵重整化群（DMRG）计算进行了基准测试。

核心贡献与结果

约瑟夫森耦合的解析表达式：
作者推导了在有限超导能隙 $\Delta$ 下，常规约瑟夫森能量（ $E_J$ ）和自旋依赖约瑟夫森能量（ $E_S$ ）的闭合形式解析表达式。
- 常规耦合（ $E_J$ ）： 在杂化强度 $\Gamma^2$ 的二阶项处恢复，这与多能级量子点的四阶微扰论结果一致。
- 自旋依赖耦合（ $E_S$ ）： 在直接隧穿振幅（ $W$ ）和杂化强度（ $\Gamma$ ）的领先阶推导。表达式表明， $E_S$ 相对于原子极限被 $\Delta/(\Delta + U/2)$ 因子所抑制。
双准粒子构型的作用：
一个核心发现是，尽管 2-QP 构型的概率权重很小，但它们对于自旋分裂的大小至关重要。
- 截面相干性（Intersector Coherence）： 自旋依赖耦合源于 0-QP（真空）与 2-QP 截面以及 1-QP 截面之间的相干性。
- 定量贡献： 在零库仑排斥（ $U=0$ ）时，0–2-QP 干涉项贡献了恰好一半的领先自旋依赖耦合。随着 $U$ 的增加，该干涉项主导了 1-QP 的贡献，提供了耦合的领先 $1/U$ 行为。
准确性与收敛性：
- 基准测试： 第 3 层变分近似（包含 2-QP 反馈）在能量和波函数分量上与 NRG 和 DMRG 表现出极佳的一致性，其差异主要源于忽略了 3-QP 及更高阶截面。
- 截断误差： 在 2-QP 空间内由核截断引入的误差被发现很小（在 $\phi=0$ 时能量误差小于 0.255%）。在强杂化区域与 NRG 产生的较大偏差归因于排除了高阶 QP 截面，而非核近似本身。
自旋重新分布与磁响应：
- SOC 诱导的转移： SOC 的引入导致了量子点与引线之间自旋的转移。这种重新分布通过局部矩分数 $q_d$ 的变化来体现。
- 共同约瑟夫森电流： 当量子点上的塞曼能与引线中的塞曼能不同（ $b_{rel} \neq 0$ ）时，SOC 诱导的自旋转移会产生一个对两种比特态都通用的约瑟夫森电流贡献。即使在零相位偏置下，该电流也是有限的，这代表了由 SOC 与磁场梯度相互作用驱动的反常约瑟夫森效应（ $\phi_0$ -结行为）。
空间结构：
作者绘制了虚准粒子（virtual QPs）的定位以及引线中物理自旋密度的分布图。
- 定位： 虚准粒子延伸至 BCS 相干长度 $\xi_0$ 。2-QP 分量由每个通道中的单个自然轨道主导，其占据概率超过 95%。
- 自旋密度： 引线中的物理自旋密度为正，但点-引线自旋相关函数随距离发生符号改变。0–2-QP 干涉项主导了自旋的空间重新分布，但并未改变积分自旋值。

意义
本文在相互作用量子点模型与用于设计 ASQ 电路的有效 $E_J$ – $E_S$ 哈密顿量之间建立了严密的微观联系。它阐明了波函数结构——特别是不同 QP 截面之间的相干性——是如何决定耦合强度的。研究结果为电路设计提供了定量约束，表明库仑排斥对常规耦合和自旋依赖耦合的抑制程度不同（在大 $U$ 极限下， $E_J \propto U^{-2}$ 而 $E_S \propto U^{-1}$ ）。这意味着虽然相对自旋依赖性（ $E_S/E_J$ ）可能会随相互作用增强，但自旋依赖耦合的绝对量级会下降。该工作还识别了负责自旋分裂的具体机制（0–2-QP 相干性），提供了一个比以往唯象或低阶微扰方法更准确的预测框架。

类似论文