Transmon 阵列中的量子孤子及其量子行走

论文信息

标题: Quantum solitons and their quantum walks in transmon arrays

作者: Ben Blain, Giampiero Marchegiani, Luigi Amico, et al.

发布日期: 2026-06-17

arXiv ID: 2606.19339v1

PDF 链接: 下载 PDF

3 分钟速览

• 研究问题：这篇论文研究如何在超导量子比特（transmon）阵列中构造和制备空间局域化的量子孤子态，并分析其量子漫步行为，以区分真正的量子孤子与简单的 “玻色子堆叠” 态。

• 核心方法：利用 attractive Bose-Hubbard 模型描述电容耦合的 transmon 阵列，通过孤子带本征态构造三种空间局域化孤子；重点提出了一种非绝热制备协议，在极短时间内（$t_0^{\text{opt}}\ll 1/J$）将玻色子堆叠态演化成高保真度的量子孤子态。

• 关键结果：通过优化常值钉扎强度和演化时间，可制备出与目标量子孤子保真度高达 99.999% 的态（$N=5$ 时），且该方法能有效补偿有限频率调谐速率带来的保真度损失。

• 主要局限：实验上需要实现对单个 transmon 频率的快速动态调控；研究限于少体系统（$N=2$ 到 $N=5$）和小型阵列（$M=11$）；量子孤子的演化时间受限于量子比特的退相干时间 $T_1$。

• 适合读者：从事超导量子计算、量子模拟、Bose-Hubbard 模型实验、以及量子孤子物理和量子漫步理论研究的科研人员和研究生。

论文背景和研究动机

量子模拟已成为现代多体物理的核心工具，可控平台使得经典算法难以处理的问题可以被实验探索。在众多候选平台中，超导电路——特别是 transmon qubit 阵列——展现了模拟复杂格点模型的多功能性。传统上，弱非谐 transmon 的高能级被视为引起量子门泄漏的 “祸害”，但论文指出，这些高能级可以被视作 qubit 或 qudit，从而扩展局部 Hilbert 空间，成为宝贵的量子模拟资源。

本文关注的焦点是 attractive Bose-Hubbard 模型在 transmon 阵列中的实现。与 repulsive 相互作用系统不同，这类系统中粒子倾向于 “堆积” 在特定格点上，形成强关联。这种局域化的粒子堆积并非体现在简单的局域密度观测量上，而是隐含在基态的密度-密度关联中——关联函数的指数衰减信号表明，基态可被理解为多个孤子的叠加。

论文承接了前期工作的重要发现：通过在 Bose-Hubbard 基态中引入钉扎位点，可以选择出单个孤子；释放钉扎后的短时演化展现出量子漫步的特征。本文进一步将量子孤子定义得更加精确，提出适用于当前超导硬件的制备方案，并对孤子的量子漫步进行定量分析，系统比较它与单粒子漫步和玻色子堆叠漫步的区别。

研究者对量子漫步的兴趣来自两个层面。基础层面，量子干涉模式生动体现了波粒二象性；对于不可区分多粒子漫步，它能揭示相互作用引起的关联。应用层面，量子漫步是量子搜索算法的自然框架（可达平方加速），甚至能实现通用量子计算。对于量子孤子而言，多粒子漫步还可能实现量子增强测量。

核心方法和技术细节

1. Bose-Hubbard 模型与孤子带

论文从一维吸引性 Bose-Hubbard 哈密顿量出发：

$$\frac{\hat{\mathcal{H}}_{\rm BH}}{\hbar}= \sum_{i=1}^{M-1}J_{i,i+1}(b_i^{\dagger}b_{i+1}+b_{i+1}^{\dagger}b_i)+\frac{U}{2}\sum_{i=1}^{M}\hat{n}_i(\hat{n}_i-1) +\sum_{i=1}^{M}\mu_i(t)\hat{n}_i$$

其中 $U<0$ 为吸引相互作用，$\mu_i$ 为位置相关的化学势，可随时间动态调节（这对 transmon 平台的协议至关重要）。当相互作用强度超过临界值 $U_C\simeq J[2N/(N-1)!]^{1/(N-1)}$（见论文 Sec. II），$M$ 个所有玻色子占据同一位点的本征态形成一个分离的 “孤子带”。在强吸引区（$-U\gg J$），该孤子带可等效地用复合粒子的有效跃迁 $J̃\simeq J\frac{N}{(N-1)!}(J/|U|)^{N-1}$ 和边界诱导的有效化学势描述。

2. 空间局域化量子孤子的三种构造方法

论文提出了三种仅用孤子带本征态构造空间局域化孤子的协议（见 Sec. II.1）：

• 投影玻色子堆叠（projected boson stack）：将 Fock 态 $|N_i\rangle$ 投影到孤子带上并归一化： $|\psi_P^{(i)}\rangle=\frac{\mathcal{\hat{P}}_B|N_i\rangle}{\langle N_i|\mathcal{\hat{P}}_B|N_i\rangle}$

• 单堆叠孤子（single-stack soliton）：构造孤子带本征态的线性叠加，只与选定的位点 $i$ 上的玻色子堆叠有非零交叠。

• 最窄孤子（narrowest soliton）：通过 Trust-Region 优化算法最小化态的宽度 $R[|\psi\rangle]=\sqrt{\frac{1}{N}\sum_{l=1}^{M}\langle\psi|\hat{n}_l|\psi\rangle(l-i)^2}$。

三种构造给出的态在长度尺度（关联长度 $\xi_C$、密度衰减长度 $\xi_N$、宽度 $R$）和彼此间的 infidelity 上差异极小（$<10^{-3}$ 对 $N=2$，$\sim 10^{-6}$ 对 $N=3$，见表 1），证实了它们物理等效。

密度衰减长度可通过变分法解析近似（见论文公式(8)及附录 A）：

$$\frac{1}{\xi_N}\simeq 2\ln\left[\frac{|U|(N-1)}{J}\right]$$

这表明该局域孤子的衰减长度小于晶格常数，在连续极限中没有对应。

3. Transmon 阵列中的非绝热制备协议

Transmon 芯片中，电容耦合的阵列可用上述 BH 模型描述，其中化学势 $\mu_i$ 对应于 $|0\rangle\to|1\rangle$ 跃迁频率，$U\simeq -E_C/\hbar$ 为充电能。关键在于 transmon 频率可通过磁通快速调谐。

论文指出，传统绝热方案（从极强钉扎绝热演化到最优钉扎）面临严重挑战：需要穿越能级间距极小的区域（$\sim\pi J̃/(M+3-2\lceil N/2\rceil)$），这要求绝热时间 $\gg M/J̃$，可能与退相干时间 $T_1$ 相当。

因此论文提出非绝热方案：利用玻色子堆叠 $|N_i\rangle$ 与相邻态 $|N-1_i 1_{i\pm 1}\rangle$ 间通过跃迁的直接耦合。在三态子空间中解析求解 Schrödinger 演化（见附录 B），发现通过选择最优钉扎强度 $\mu^{\rm opt}\simeq U(N-1)$ 和演化时间 $t_0^{\rm opt}\simeq\pi/2|U|(N-1)$，可使态演化至与目标局域孤子最高保真度（见图 4、图 5）。

该协议通过在时间 $t_0$ 处将钉扎突然切换至 0（quench），得到均匀系统中的高保真度孤子态。即使考虑到实验上有限频率调谐速率，通过适当调整准备时间 $t_0$，保真度损失可被完全补偿（见附录 C）。

创新点和贡献

1. 明确构造空间局域态：论文首次系统地用三种等效方法从孤子带本征态构造真正的量子孤子态，并与简单玻色子堆叠、钉扎基态进行了定量区分（见表 1 中的 infidelity、关联长度等）。

2. 提出了硬件适配的制备协议：非绝热协议是本文的核心实用贡献。它避免了传统方案中极缓慢的绝热演化需求，在 $\ll 1/J$ 的时间尺度内实现 $>99.8\%$ 的保真度（$N=5$ 时 $>99.999\%$，见图 5(a)），且对有限 quench 速率具有鲁棒性。

3. 定量量子漫步分析：利用 Hellinger 距离（公式(16)）和端点位点粒子数演化，定量比较了量子孤子、玻色子堆叠和单粒子的量子漫步。发现量子孤子在整个演化过程中保持 “复合粒子” 特性——尤其在监测端点位点 transmon 的 $|1\rangle$ 态概率时，孤子表现为单峰，而玻色子堆叠则有更早的单粒子成分到达（见图 8(c)、图 9(b-c)）。

4. 可实验验证的判别方法：提出通过监测量子漫步的端点位点占据数来区分孤子与堆叠态。这不仅基于期望值，还可扩展至单个 transmon 的能级分辨读出（$\mathcal{P}_i(\mathcal{N})$），实验中可通过多能级单发读出实现。

实验结果分析

论文的 “实验” 主要是数值精确对角化和 Trotterized 动力学模拟，但所有参数均对标实际 transmon 参数选取（$E_J/E_C\simeq 50-100$，$J\simeq 1-40$ MHz）。

1. 最优钉扎与非绝热制备

图 2 展示了投影玻色子堆叠与钉扎基态间的 infidelity 随钉扎强度 $\Delta\mu_6$ 的变化。infidelity 呈非单调行为，在最优钉扎处达到极小。最优钉扎强度的绝对值随 $N$ 增大而减小（$N=2$ 时约为 $-3.66J$，$N=3$ 时约 $-1.13J$，见论文描述）。

图 4(a)展示了 $N=2$ 玻色子堆叠在 $\Delta\mu_6=-11.6J$ 常值钉扎下的动力学：中心位点粒子数仅微弱减少，相邻位点在极短时间内（$t_0^{\rm opt}\ll 1/J$）被布居。图 4(b-c)的 infidelity 曲线显示，与目标孤子的 infidelity 在 $t_0$ 处达到极低极小值，且与孤子带投影的 infidelity 几乎一致，证实了主要布居在孤子带内。

图 5(a)展示了在时间范围 $0\ldots 1/J$ 内演化后与目标孤子的最小 infidelity。当钉扎接近 $U(N-1)$ 时，infidelity 达到极小。随 $N$ 增大，infidelity 显著降低——$N=2$ 时为 $0.2\%$，$N=5$ 时达到 $\sim 10^{-5}$ 量级，展示了方案的优越扩展性。

2. 量子漫步对比分析

图 6 和图 7 分别对 $N=2, U=-10J$ 和 $N=3, U=-2J$ 两种情况，展示了单粒子、玻色子堆叠和量子孤子的量子漫步。通过将时间用有效跃迁 $J̃_0$ 缩放，三种态的干涉图案定性相似。但关键区别在于：

• 玻色子堆叠的演化（图 6(b)）上叠加了快速传播（速度 $\sim J$）的单粒子组分；
• 量子孤子（图 6(c)）完全消除了这些 “宽带” 组分，干涉图案更清晰。

Hellinger 距离（图 6(d)）定量证实了这一点：在短时范围内，孤子与单粒子的 Hellinger 距离明显小于堆叠与单粒子的距离。

端点位点占据数（图 6(e)、图 7(e)）提供了最清晰的实验判别途径。堆叠态的单粒子组分在时间 $\sim (M-1)/2\sqrt{2}J$ 到达端点，而孤子态则需更长时间 $\sim (M-1)/2\sqrt{2}J̃_0$。纸中 “A” 曲线（橙线）展示了按附录 C 优化协议制备的实际态——其端点密度两个数量级低于堆叠态，与理想孤子几乎一致。

3. Transmon 本征态分辨测量

图 8 和图 9 进一步分析了单发读出的能级分辨结果。量子孤子的 $\mathcal{P}_i(1)$ 和 $\mathcal{P}_i(2)$（或 $\mathcal{P}_i(2)$ 和 $\mathcal{P}_i(3)$，对 $N=3$）干涉图案几乎相同，体现了复合粒子行为。而玻色子堆叠的 $\mathcal{P}_i(1)$ 包含明显的快速移动组分。端点 transmon 的 $\mathcal{P}_1(1)$ 监测清晰捕捉到这一差异。

实践建议

对于希望在超导量子芯片上实现这一工作的实验团队，可以按以下步骤推进：

1. 器件参数选择：选择 $E_J/E_C$ 比值在 50-100 范围的 transmon，确保非谐性 $-E_C$（即 $U$）远大于跃迁强度 $J$，但又不致于过大而使 $U_C$ 难以达到。建议 $U/J$ 设置在-10 到-2 这个范围（论文中用的两个典型值）。

2. 钉扎标定：利用磁通脉冲对单个 transmon 的频率进行快速调控。根据论文公式 $\mu_{\rm opt}\simeq U(N-1)$ 设置钉扎强度初值，然后通过测量演化后与目标孤子的 infidelity（可用态层析或密度分布测量）精细调优。注意 $N=2$ 时实验上需要 $\Delta\mu \approx -11.6J$（论文数值结果），而非简单的 $-10J$。

3. 时序控制：非绝热方案要求演化时间 $t_0^{\rm opt}\simeq\pi/2|U|(N-1)$。对于典型的 $E_C/h=200$ MHz，$N=2$ 时 $t_0^{\rm opt}\approx 1.25$ ns。制备完成后立即将频率 quench 回均匀值。如果频率调谐速率有限（典型值 1 GHz/10-20 ns），需按附录 C 所示减短准备时间以补偿。

4. 数据采集与判别：不只是测量平均粒子数 $\langle\hat{n}_i\rangle$，更要利用高保真单发读出获得 $\mathcal{P}_i(\mathcal{N})$ 分布。比较端点位点第 1 激发态概率的时间演化——孤子态为单峰慢速到达，堆叠态则先出现快速前兆。这比整体密度分布的 Hellinger 距离更易实现。

5. 进一步方向：论文未讨论但是可行的是，利用这一平台探索量子孤子对旋转、磁场的增强测量灵敏度（论文 Sec. V 中提及），以及多孤子态的结合制备——这些都可能受益于本文建立的坚实基础。