超导量子比特阵列中的量子孤子及其量子行走

arXiv: 2606.19339v1

论文信息

标题: 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

论文背景和研究动机

量子模拟已成为多体物理研究的重要支柱,它借助高度可控的实验平台,探索经典计算难以处理的复杂量子体系。在众多候选架构中,由 transmon 量子比特构成的超导电路阵列凭借其灵活的设计与操控能力,成为模拟格点模型的突出平台。传统观念中,弱非谐 transmon 的高能级是实现高保真量子门时需抑制的 “泄漏” 通道,但若善加利用,这些高能级可将 transmon 扩展为 qutrit 或 qudit,从而大大丰富局域希尔伯特空间,为更复杂的量子模拟提供资源。

在这一背景下,探索具有吸引相互作用的玻色-哈伯德模型的物理意义尤为深远。与排斥相互作用体系不同,吸引玻色子的低能动力学以粒子在特定格点 “聚集” 为特征,导致局域希尔伯特空间急剧膨胀,并伴随着强关联效应,给经典数值方法带来巨大挑战。这种 “聚集” 行为在低能本征态的局域密度中并不明显,但在基态的密度-密度关联函数中表现为指数衰减,这暗示基态可以理解为空间局域化量子孤子的叠加态。此前已有工作借助 “钉扎” 位点从基态中挑选孤子并研究其释放后的短时演化,并发现该演化展现出量子行走的特征。量子行走不仅是波粒二象性的生动体现,也在量子搜索和通用量子计算中扮演关键角色,尤其是多粒子量子行走有望实现量子增强测量。因此,在超导电路中制备并深入分析量子孤子及其量子行走,不仅具有基础物理意义,也为量子技术的应用开辟了新途径。

核心方法和技术细节

本文以包含位点依赖化学势和跃迁系数的玻色-哈伯德模型为出发点:

H^BH=i=1M1Ji,i+1(b^ib^i+1+b^i+1b^i)+U2i=1Mn^i(n^i1)+i=1Mμi(t)n^i,\frac{\hat{\mathcal{H}}_{\rm BH}}{\hbar}= \sum_{i=1}^{M-1}J_{i,i+1}(\hat{b}_{i}^{\dagger}\hat{b}_{i+1}+\hat{b}_{i+1}^{\dagger}\hat{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<0U<0 为吸引相互作用强度。对于均匀系统,相互作用足够强时,谱中会出现一个包含 MM 个本征态的孤子能带,其能量显著低于其他能态。在这些态中,所有玻色子倾向于占据同一位点,构成有效复合粒子,其有效跃迁幅度为 J~JN/(N1)!(J/U)N1\tilde{J}\simeq J\,N/(N-1)!\, (J/|U|)^{N-1},并伴随由边界效应导致的等效化学势偏移。

为了明确定义空间局域化的量子孤子,作者提出并比较了三种构造方案:将 NN 个玻色子全部置于位点 ii 的福克态(称为玻色子堆栈)向孤子能带投影得到投影玻色子堆栈;要求构造的叠加态仅与指定位点的玻色子堆栈有重叠,从而定义单堆栈孤子;以及在孤子能带内通过优化最小化空间宽度得到的最窄孤子。数值计算表明,对于强相互作用参数(例如 U=10JU=-10J),这三种态之间的保真度极高(差异小于 10310^{-3}),且具有几乎相同的衰减长度和关联长度,因此后续研究可统一采用投影玻色子堆栈作为局域量子孤子的代表。

制备孤子的协议针对超导 transmon 阵列设计。Transmon 的等离子体激发可映射为玻色子,其跃迁频率对应化学势,电荷能提供吸引相互作用,电容耦合则产生相邻位点的跃迁。通过外加磁通可实现化学势的快速调节。传统上,可以通过制备最强钉扎下的基态再绝热缓释钉扎来逼近孤子态,但分析表明,由于能谱中存在能级间距极小的区域,绝热过程所需时长往往超过系统的相干时间,实验上不可行。

为此,作者提出了一种非绝热的短时演化协议:首先在目标位点上以强钉扎 Δμ\Delta\mu 制备一个玻色子堆栈,保持该钉扎演化一个精心选择的短时间 t0t_0。在仅考虑堆栈态与相邻位点的单粒子跳跃态构成的三维子空间中,可解析求解系统演化。通过要求演化后相邻位点的粒子数占比与目标孤子态中的占比匹配,得到最优钉扎强度 μoptU(N1)\mu^{\text{opt}}\simeq U(N-1) 和相应的最优保持时间 t0optπ/(2U(N1))t_0^{\text{opt}}\simeq \pi/(2|U|(N-1))。在该时刻迅速将钉扎撤去(淬火),即可获得高保真度的局域孤子态。进一步考虑现实实验中的有限淬火速率,通过微调保持时间,仍能恢复与理想瞬时淬火相当的高保真度,展示了该协议在实验条件下的鲁棒性。

创新点和贡献

本文的创新之处体现在多个层面。首先,在理论构建上,系统给出了玻色-哈伯德模型中局域量子孤子的三种精确定义,并证明了它们在强相互作用极限下的等价性,这为后续研究提供了统一且清晰的量子孤子概念。

其次,所提出的非绝热制备协议充分发掘了 transmon 阵列的动态调控能力,巧妙地绕过了绝热方案遇到的能级拥挤瓶颈,将制备时间缩短到远小于典型弛豫时间的尺度,使得在现有实验条件下实现高保真度(最高可达 99.999% 以上)的量子孤子成为可能。该协议对参数的选择给出了明晰的解析指引,并考虑了有限时间淬火的补偿方法,具有极强的实验可操作性。

最后,在量子行走分析中,作者首次定量比较了量子孤子、玻色子堆栈和单粒子的演化行为。通过引入 Hellinger 距离以及监测阵列端部 transmon 的各激发态占据概率,鲜明地区分了三者:堆栈中会混入快速传播的单粒子成分,而量子孤子的演化则表现出清晰的复合粒子干涉图案,且端部出现激发的时间显著滞后。这些判据为实验验证量子孤子的多体相干性提供了切实可行的观测方案,也深化了对量子多粒子行走中关联效应的理解。

实验结果分析

文章的数值实验涵盖了 N=2N=2N=3N=3 的玻色子数,相互作用强度覆盖从强相互作用 (U=10JU=-10J) 到略高于孤子能带形成的临界值 (U=2JU=-2J),阵列长度 M=11M=111515,充分验证了方法的普适性。

在态制备部分,最优钉扎下 N=2N=2 时保真度达 99.8%,N=5N=5 时超过 99.999%,且最佳演化时间随 NN 增大而缩短。当引入有限淬火时间 tq=0.01/Jt_q=0.01/J 后,仅需将保持时间从 0.1377/J0.1377/J 微调至约 0.134/J0.134/JN=2N=2 情况),便可将保真度恢复至与瞬时淬火相同的水平,表明该协议对实际硬件上的频率调制速率具有很好的适应性。

量子行走的对比模拟直观地展示了三者在动力学上的差异。以 N=2N=2U=10JU=-10J 为例,若将时间按有效跃迁 J~0\tilde{J}_0 缩放,孤子的期望粒子数演化图案与无相互作用单粒子的图案高度相似,而玻色子堆栈的图案则夹杂了前沿快速传播的微弱信号。Hellinger 距离的量化分析显示,在短时尺度内孤子与单粒子分布的距离明显小于堆栈与单粒子的距离。更为显著的是,观察端部 transmon 的第一激发态概率,堆栈在大约 5/J5/J 的时间内便出现明显峰,而孤子的相应峰值则延迟到约 17/J17/J,并且其间无早期杂散峰,这为实验区分两者提供了清晰的 “延迟到达” 信号。在较弱相互作用 U=2JU=-2JN=3N=3 的情况下,虽然干涉图案总体模糊,但孤子的优势依然保持,端部第二激发态和第一激发态的演化一致性更彰显了其复合粒子特性。

实践应用建议和未来发展方向

从实验实现角度,本文建议在具有可调频 transmon 的阵列上直接验证所提协议。实验中可用微波脉冲将某位点激发至 NN 光子福克态,再利用单独磁通偏置线将该位点频率按最优钉扎偏移并保持 t0t_0 时间,最后快速扫回共振频率。读出则可借助单发色散读取技术区分 transmon 的 0\vert 0\rangle1\vert 1\rangle2\vert 2\rangle 态,从而重构出各格点的粒子数分布乃至各激发态占据概率,直接观察孤子量子行走的干涉图样以及端部到达时间。这些测量并不需要复杂的量子层析,使得验证相对简便。

在应用层面,量子孤子作为具有确定粒子数且空间局域的纠缠态,在量子计量学中可充当类似 NOON 态的探针,有望提升对相位、微弱力或惯性效应的测量灵敏度。同时,其稳健的复合粒子行走特性也为构建基于量子行走的搜索算法或多粒子干涉仪提供了新的资源。进一步,探索更大规模的阵列、更多粒子数的孤子,以及引入无序或规范场后的孤子输运行为,将有望在超导电路平台上模拟更广泛的非平衡量子多体现象。

理论和方法上,可进一步研究环境中退相干与耗散对孤子制备保真度和量子行走的影响,发展具有主动纠错或动态解耦的增强方案。此外,将本文的制备思想推广到二维阵列或具有拓扑性质的晶格中,可能催生对拓扑量子孤子或分数统计粒子模拟的新研究。

总结与展望

该工作系统性地在吸引玻色-哈伯德模型中构建了局域量子孤子理论,并给出了与当前超导 transmon 技术完美兼容的制备方案,通过量子行走的详细分析揭示了这类复合粒子的独特动力学特征。论文不仅清晰定义了量子孤子的多种等价构造,更提供了具有解析指导的非绝热制备流程,弥合了抽象模型与实验实现之间的鸿沟。数值结果有力地证明,即便在有限淬火速率下,所提协议仍能以极高保真度产生目标态,并利用简单可观的端部信号甄别量子孤子与平庸的福克态堆栈。

展望未来,随着超导量子处理器规模的扩展与相干性能的持续提升,在实验上实现并操控量子孤子将水到渠成。这一平台将为研究量子多体行走的干涉、束缚态输运以及量子增强测量提供独特窗口,并有望推动基于量子模拟的精密测量技术与新型量子算法的发展。