来源论文: https://arxiv.org/abs/2603.04498v1 生成时间: Mar 06, 2026 21:46

0. 执行摘要

本研究提出了一种在超冷原子系统中实现手性与配对超流体的创新平台，该平台基于态相关的 Kronig-Penney 晶格，通过巧妙设计的三足式（tripod-type）原子-光耦合方案构建。这种晶格设计独特之处在于它能同时引入强烈的几何挫折和非标准的二体相互作用，特别是可控的配对跃迁（pair hopping）和密度诱导跃迁（density-induced tunneling），而这些在传统光学晶格中通常可以忽略。通过大规模密度矩阵重正化群（DMRG）计算，我们首次绘制了相应扩展 Bose-Hubbard 模型的量子相图。结果表明，配对跃迁能够稳定鲁棒的配对超流体（PSF），其特征是配对关联函数呈现准长程衰减，并有效抑制单粒子关联。此外，由竞争性的最近邻（NN）和次最近邻（NNN）跃迁导致的几何挫折，诱导了时间反演对称性自发破缺，从而产生了具有长程电流关联的手性超流体（CSF）。在高势垒极限下，该系统可以精确映射到反铁磁 XXZ 旋量链，这不仅使得相变点能够通过解析方法确定，也为数值模拟结果提供了重要的基准，极大地验证了理论模型的准确性与DMRG计算的可靠性。

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

核心科学问题

传统光学晶格中的超冷原子系统是量子模拟的理想平台，但其在实现强关联效应、特定几何结构（如三角梯形）中的几何挫折以及可控的非标准多体相互作用（如配对跃迁和密度诱导跃迁）方面存在局限性。这些特殊的物理机制是生成非常规量子相的关键驱动力。本研究的核心科学问题是：如何设计并实现一个能够在超冷原子系统中同时具备显著几何挫折和可调控的关联跃迁（特别是配对跃迁和密度诱导跃迁）的平台，并探究这些独特物理机制如何共同作用，诱导产生新颖的量子物态？

具体而言，研究人员希望回答以下问题：

  1. 如何通过原子-光耦合技术构建一个具有特定几何（如三角梯形）和态相关势垒的晶格，从而在系统中引入几何挫折？
  2. 如何利用态相关的原子-原子相互作用，自然地产生可调控的配对跃迁和密度诱导跃迁，并量化其强度？
  3. 在同时存在几何挫折和关联跃迁的复杂多体系统中，会涌现出哪些新颖的量子相？这些相的特征是什么？
  4. 这些新奇量子相（如配对超流体和手性超流体）是如何稳定形成的？它们的相变边界如何界定？
  5. 在某些极限情况下，系统能否映射到已知的可解模型（如XXZ自旋模型），从而提供解析洞察力并验证数值结果？

理论基础

本研究的核心理论基础在于态相关光学晶格的设计与有效哈密顿量的推导。

态相关 Kronig-Penney 晶格

传统的Kronig-Penney模型描述了周期性势垒对粒子运动的影响。这里，研究人员引入了“态相关”的概念，意味着势垒对不同内部态的原子散射方式不同。具体实现方式是通过一个三足式原子-光耦合方案（tripod atom-light coupling scheme），该方案能够产生两个独立的“暗态”（dark states）。暗态是原子-光耦合哈密顿量的零本征态，意味着处于这些态的原子不受光场影响，从而减少了自发辐射和加热，是超冷原子实验中常用的技术。

论文中，原子内部态被选为三个基态 $|1 angle, |2 angle, |3 angle$ 中的两个暗态 $|D_1 angle$ 和 $|D_2 angle$。这些基态通过拉比频率 $\Omega_l(x)$ 耦合到一个共同的激发态 $|e angle$。拉比频率的特定空间依赖性是实现态相关势垒的关键：

$$ \Omega_1(x) = \Omega_p \ \Omega_2(x) = \Omega_p \cos(2\pi x/a) \ \Omega_3(x) = \Omega_c \sin(2\pi x/a) $$

其中，$a$ 是光学波长，$k$ 是波数。通过设定 $\Omega_{p,2}(x)$ 远小于 $\Omega_3(x)$（即 $\epsilon = \Omega_p/\Omega_c \ll 1$），可以形成宽度为 $\epsilon a$ 的亚波长势垒。这些势垒在 $\Omega_3(x)=0$ 的位置（即 $x=na/2$）处形成，它们对暗态 $|D_2 angle$ 中的原子产生影响，导致原子在奇偶 $n$ 对应的势垒处分别经历对称和反对称的原子基态叠加态 $|+ angle = (|1 angle+|2 angle)/\sqrt{2}$ 和 $|- angle = (|1 angle-|2 angle)/\sqrt{2}$ 的散射。这使得 Wannier 函数 $w_n(x)$ 呈现出交错的性质，即 $n \in 2Z+1$ 对应 $|+ angle$ 态，而 $n \in 2Z$ 对应 $|- angle$ 态。

有效哈密顿量：扩展 Bose-Hubbard 模型

在最低 Bloch 能带的 Wannier 函数基下，整个原子系统（包括原子-原子相互作用）可以被一个扩展的 Bose-Hubbard 型哈密顿量描述：

$$ H = H_{ ext{tun}} + H_{ ext{int}} $$

1. 隧穿项 ($H_{ ext{tun}}$): 传统光学晶格中，通常只有最近邻隧穿。而本方案通过独特的态相关势垒，引入了最近邻（NN）和次最近邻（NNN）隧穿，且 NNN 隧穿强度 $|J_2|$ 大于 NN 隧穿强度 $|J_1|$，并且 $J_2$ 为负值。这种 $J_2 < 0$ 和 $|J_2| > |J_1|$ 的特性是产生几何挫折的关键。具体的，NN 隧穿是由于不同内部态 $|+ angle$ 和 $|- angle$ 之间 Wannier 函数的微小混合重叠产生，而 NNN 隧穿则是原子跨越势垒而不改变内部态的隧穿。这种特殊的隧穿结构将系统映射到一个带有 $\pi$ 磁通的三角梯形结构上，如图3(a)所示。

$$ H_{ ext{tun}} = -J_1\sum_j (\hat{a}_j^\dagger \hat{a}_{j+1} + ext{h.c.}) - J_2\sum_j (\hat{a}_j^\dagger \hat{a}_{j+2} + ext{h.c.}) $$

2. 相互作用项 ($H_{ ext{int}}$): 原子-原子相互作用 $H_{ ext{int}}$ 的独特之处在于其态相关性。通过散射长度 $g_{pq}$ 和 Wannier 函数重叠积分 $G_{ijk}$ 的参数化，可以自然地导出以下几种重要的相互作用：

  • 在位相互作用 ($U$)：原子在同一格点上的相互作用，由 $G_{000}$ 决定。
  • 最近邻相互作用 ($V$)：原子在相邻格点上的相互作用，由 $G_{011}$ 决定。
  • 密度诱导跃迁 ($D$)：一种关联隧穿，其中原子在跃迁时会受到相邻格点密度的影响。这在论文中被设为零 ($D=0$)，对应于原子在相同内部态之间的散射长度相等的情况。
  • 配对跃迁 ($P$)：这是本研究的亮点之一，指的是一对原子同时从一个格点跃迁到另一个格点，或者从相邻格点分别跃迁到同一格点。其强度由 $P = g_x G_{011}/2$ 决定。这是由原子-原子相互作用的态相关性自然产生的。

相互作用哈密顿量表示为：

$$ H_{ ext{int}} = rac{U}{2}\sum_j \hat{n}_j (\hat{n}_j-1) + V\sum_j \hat{n}_j \hat{n}_{j+1} - D\sum_j \hat{a}_j^\dagger (\hat{a}_{j-1} + \hat{a}_{j+1}) + ext{h.c.} - P\sum_j \hat{a}_{j+1}^\dagger \hat{a}_{j+1}^\dagger \hat{a}_j \hat{a}_j + ext{h.c.} $$

其中，参数 $U, V, D, P$ 均通过散射长度 $g_0, g_x$ 和 Wannier 函数的重叠积分 $G_{000}, G_{001}, G_{011}$ 进行控制。值得注意的是，为了稳定配对超流体并防止亮孤子塌缩，研究人员额外施加了一个三体约束，即每个格点最多只能容纳两个玻色子。这可以通过三体损耗或三体相互作用在实验中实现。

技术难点

  1. 精确实现态相关 Kronig-Penney 晶格：

    • 拉比频率的精密控制： 需要精确控制三个拉比频率 $\Omega_1(x), \Omega_2(x), \Omega_3(x)$ 的空间依赖性，以确保在正确位置形成亚波长势垒。特别是 $\Omega_p/\Omega_c \ll 1$ 的条件需要精细调节。
    • 暗态流形的限制： 要求原子的质心能量（频率单位）远小于总拉比频率 $\Omega(x)$，以确保绝热近似的有效性，将系统动力学限制在暗态流形内。
    • Wannier 函数的交错性质： 确保原子在偶数和奇数势垒处分别散射到对称和反对称叠加态中，从而产生所需的 Wannier 函数交错性质，是实现三角梯形结构和挫折的关键。

  2. 生成和控制非标准相互作用：

    • 态相关散射长度： 原子-原子相互作用的强度 $g_{pq}$ 依赖于内部态，这通常需要通过 Feshbach 共振等技术进行精确调控。
    • Wannier 函数重叠积分： 相互作用参数 $U, V, D, P$ 依赖于 Wannier 函数的重叠积分 $G_{ijk}$。这些积分对势垒形状和宽度高度敏感，需要精确的晶格参数控制。
    • 配对跃迁的强度： 实现足够强的配对跃迁 $P$ 是稳定配对超流体的关键。这要求 $g_x$ 和 $G_{011}$ 具有适当的值。论文指出，$p$ 波段的 Wannier 函数可能更有利于实现更强的配对跃迁，但实验上更难实现。

  3. 几何挫折与关联跃迁的共存与竞争：

    • NNN 隧穿的主导作用： $J_2 < 0$ 且 $|J_2| > |J_1|$ 是引入几何挫折的先决条件。在实验中精确控制 $J_1$ 和 $J_2$ 的相对强度和符号是一项挑战，尤其是在 $J_1$ 和 $J_2$ 都相对较小的情况下。
    • 多体效应的复杂性： 几何挫折和关联跃迁（尤其是配对跃迁）的同时存在，使得系统的多体行为变得异常复杂，难以直观理解，需要强大的数值模拟工具。

  4. 相图的精确确定：

    • DMRG 计算的精度和效率： 对于强关联系统，DMRG 需要足够大的键维度（bond dimension）才能捕获量子纠缠，这会带来高昂的计算成本。在具有开放边界条件和单位填充的较大系统上，需要仔细控制截断误差，并使用中心区域数据来最小化边缘效应。
    • 相变点的识别： 区分不同量子相（如 MI, DW, PSF, CSF）需要计算多种序参量和关联函数，并精确识别其行为变化。特别是 BKT 相变（如 DW-PSF 转变）的识别，其特征是关联函数的代数衰减而非突然消失，这比一级相变或二级相变更具挑战性。

  5. 实验观测：

    • 手性序的测量： 手性超流体中的时间反演对称性破缺通过电流-电流关联函数 $K_2(\Delta j)$ 来表征。在实验中，这可能需要通过动量分布测量或干涉测量来推断。
    • 配对关联的测量： 配对超流体中的长程配对关联 $C_2(\Delta j)$ 可以通过光缔合测量等技术进行探测。

方法细节

本研究主要采用密度矩阵重正化群（DMRG）方法，并结合高势垒极限下的自旋映射来分析系统的基态性质。

1. 密度矩阵重正化群（DMRG）模拟

DMRG 是一种强大的数值方法，特别适用于研究一维（1D）强关联量子系统的基态性质。它可以有效地处理系统中大量的量子纠缠，并能精确计算各种序参量和关联函数。

  • 模型： 扩展 Bose-Hubbard 模型，如前所述，包括 NN 和 NNN 隧穿项，以及在位、NN 和配对跃迁相互作用项。密度诱导跃迁 $D$ 在本次研究中被设为零，以简化问题并聚焦于几何挫折和配对跃迁。
  • 系统设置：
    • 格点数 ($L$)： 200个格点。
    • 原子数： 200个原子，对应于单位填充（unit filling），即每个格点平均一个原子。
    • 边界条件： 开放边界条件（open boundary conditions）。
    • 玻色子维度： 每个格点最多容纳2个玻色子。这是一个关键的三体约束，用于稳定配对超流体并防止亮孤子态的形成。论文指出，如果没有这个约束，在吸引性在位相互作用下可能会形成亮孤子。
    • 键维度 ($\chi$)： 最大键维度设为 500。键维度是 DMRG 算法中保留的态的数量，它决定了计算的精度和计算成本。更大的键维度可以更好地捕获系统的纠缠，但计算量更大。
    • 截断误差： 达到最大键维度时，截断误差（truncation error）低于 $10^{-6}$。这个误差量化了在每次重正化步骤中丢弃的小本征值（即信息损失）的总和，是 DMRG 结果精度的重要指标。
    • 期望值计算： 所有的期望值（包括关联函数和序参量）都在链的中心 160 个格点上进行评估，以最小化开放边界条件带来的边缘效应。
  • 软件： 计算使用基于 Julia 编程语言的 ITensors.jl 库 [37, 38, 39] 完成。这是一个专为张量网络计算而设计的开源库，支持 DMRG 等多种算法。

2. 量子相的表征

为了识别和区分不同的量子相，研究人员计算了以下几种序参量和关联函数：

  • 单粒子格林函数 ($C_1(\Delta j)$)： $C_1(\Delta j) = \langle \hat{a}_j^\dagger \hat{a}_{j+\Delta j} angle$。在绝缘相（Mott 绝缘体 MI 和密度波 DW）中，它会指数衰减到零。在超流体相（SF, PSF, CSF）中，由于准凝聚现象，它会以代数形式衰减（准长程关联）。
  • 配对关联函数 ($C_2(\Delta j)$)： $C_2(\Delta j) = \langle \hat{a}_j^\dagger \hat{a}_j^\dagger \hat{a}_{j+\Delta j} \hat{a}_{j+\Delta j} angle$。配对超流体（PSF）的特征是 $C_2(\Delta j)$ 呈代数衰减，而 $C_1(\Delta j)$ 呈指数衰减。SF 和 CSF 也可能表现出 $C_2(\Delta j)$ 的代数衰减，但伴随着 $C_1(\Delta j)$ 的代数衰减。
  • 偶宇称序参量 ($O_{ ext{even}}^2(\Delta j)$)： 用于表征 Mott 绝缘体（MI）相，它测量粒子数密度围绕平均值的波动模式。
  • 结构因子 ($S(k=\pi)$)： $S(k) = rac{1}{L^2}\sum_{j,l} \langle \hat{n}_j \hat{n}_l angle e^{ik(j-l)}$。在 $k=\pi$ 处达到最大值通常预示着周期为2的密度波（DW）相的出现，对应于离散平移对称性自发破缺。
  • 电流-电流关联函数 ($K_2(\Delta j)$)： $K_2(\Delta j) = \langle \hat{k}_j \hat{k}_{j+\Delta j} angle$，其中 $\hat{k}_j = i(\hat{a}_j^\dagger \hat{a}_{j+1} - \hat{a}_{j+1}^\dagger \hat{a}_j)$ 是格点 $j$ 和 $j+1$ 之间的电流算符。非零的 $K_2(\Delta j)$（长程关联）表明时间反演对称性自发破缺，这是手性超流体（CSF）的特征。

3. 高势垒极限下的自旋映射

为了提供解析洞察力并验证 DMRG 结果，研究人员在高势垒极限下进行了分析。在这种极限下，单粒子隧穿 ($J_1, J_2$) 和密度诱导跃迁 ($D$) 可以被忽略 ($H_{ ext{tun}} = 0$ 或 $J_1=J_2=0, D=0$)。系统的低能物理可以由每个格点上零个或两个玻色子占据的态（即 $|0 angle_j$ 和 $|2 angle_j$）精确描述。通过将这些玻色子态映射到赝自旋-1/2系统：$|0 angle_j ightarrow |\downarrow angle_j$ 和 $|2 angle_j ightarrow |\uparrow angle_j$，可以将原玻色子哈密顿量映射到一个反铁磁 XXZ 自旋-1/2 链哈密顿量：

$$ H_{ ext{hb}} = -J_{ ext{hb}}\sum_j (\hat{S}_j^x \hat{S}_{j+1}^x + \hat{S}_j^y \hat{S}_{j+1}^y + \Delta_{ ext{hb}} \hat{S}_j^z \hat{S}_{j+1}^z) $$

其中 $J_{ ext{hb}} = 4P = 2g_x G_{011}$ 且 $\Delta_{ ext{hb}} = -V/P = -4g_0/g_x$。由于 $g_0>0$ 和 $g_x<0$，因此 $J_{ ext{hb}}<0$，这对应于反铁磁 XXZ 模型。

  • 相对应： 在这个映射下，玻色子系统中的密度波（DW）相对应于自旋系统中的 Néel 相（$\Delta_{ ext{hb}} > 1$），而配对超流体（PSF）相对应于 XY 相（$|\Delta_{ ext{hb}}| \le 1$）。
  • 解析确定相变点： 利用 XXZ 模型的 Bethe ansatz 解，可以解析计算 Néel 相的能量，并与 MI 相的能量进行比较，从而精确确定 MI-DW 相变点。这个解析结果与 DMRG 模拟结果吻合得非常好，为数值方法提供了强有力的验证。

这些方法共同使得研究人员能够全面地探索该复杂量子系统的基态性质，从数值模拟到解析验证，确保了结果的可靠性。

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

关键 Benchmark 体系

本研究主要关注一个独特的三角梯形（triangular ladder）结构中的扩展 Bose-Hubbard 模型。这个模型通过态相关的 Kronig-Penney 晶格在超冷原子中实现。以下是主要的 benchmark 体系及其特点：

  1. 态相关 Kronig-Penney 晶格：

    • 原子-光耦合方案： 基于三足式（tripod-type）方案，将三个原子基态 $|1 angle, |2 angle, |3 angle$ 耦合到一个激发态 $|e angle$（图1a）。
    • 拉比频率： 具有特定空间依赖性 $\Omega_1(x) = \Omega_p, \Omega_2(x) = \Omega_p \cos(2\pi x/a), \Omega_3(x) = \Omega_c \sin(2\pi x/a)$（图1b）。
    • 势垒形成： 在 $\Omega_3(x)=0$ 处形成亚波长势垒，对原子对称态 $|+ angle$ 和反对称态 $|- angle$ 产生交替作用（图2a）。
    • Wannier 函数： 最低 Bloch 能带的 Wannier 函数 $w_n(x)$ 具有交错性质，中心在 $x=na/2$，其中 $n \in 2Z+1$ 对应 $|+ angle$ 态， $n \in 2Z$ 对应 $|- angle$ 态（图2b）。

  2. 扩展 Bose-Hubbard 模型（有效哈密顿量）：

    • 隧穿项： 包含最近邻（NN）隧穿 $J_1$ 和次最近邻（NNN）隧穿 $J_2$。关键特征是 $J_2 < 0$ 且 $|J_2| > |J_1|$，这引入了几何挫折并使得每个三角格点具有一个有效 $\pi$-磁通（图3a）。在数值模拟中，$J_2/J_1$ 固定为 $\approx -3.565$。
    • 相互作用项： 包括在位相互作用 $U$，最近邻相互作用 $V$，密度诱导隧穿 $D$ (本研究设为 $D=0$)，以及配对跃迁 $P$。这些参数由原子散射长度 $g_0, g_x$ 和 Wannier 函数重叠积分 $G_{000}, G_{001}, G_{011}$ 决定（公式11）。
    • 三体约束： 为防止亮孤子形成和稳定配对超流体，每个格点最多允许有2个玻色子。这意味着模型中考虑的最大玻色子占据数为 $n_{max}=2$。

  3. 高势垒极限下的 XXZ 自旋-1/2 链：

    • 近似条件： 单粒子隧穿 $J_1, J_2$ 和密度诱导隧穿 $D$ 均可忽略（即 $H_{ ext{tun}}=0$ 或 $J_1=J_2=0$）。
    • 自旋映射： 将玻色子占据态 $|0 angle_j$ 和 $|2 angle_j$ 映射到赝自旋-1/2 态 $|\downarrow angle_j$ 和 $|\uparrow angle_j$。这个近似在 $J_1=0$ 时最准确。
    • 映射模型： 得到一个反铁磁 XXZ 自旋-1/2 链哈密顿量，其参数 $J_{ ext{hb}}$ 和 $\Delta_{ ext{hb}}$ 与 $P, V, g_0, g_x$ 相关。具体地，$J_{ ext{hb}} = 2g_x G_{011}$ 和 $\Delta_{ ext{hb}} = -4g_0/g_x$。
    • 目的： 作为 DMRG 结果的解析验证，特别是在 MI-DW 和 DW-PSF 转变点。

计算所得数据

1. 量子相图 (图4a)

相图展示了基态量子相如何随参数比率 $g_x/g_0$ 和 $J_1/g_0G_{000}$ 变化。主要的量子相包括：

  • Mott 绝缘体 (MI)： 在 $J_1/g_0G_{000}$ 和 $g_x/g_0$ 都较小（强在位排斥）的区域。特征是粒子数固定，单粒子和配对关联函数呈指数衰减。
  • 密度波 (DW)： 在 $g_x/g_0$ 减小（最近邻相互作用 $V$ 减弱）时，MI 区域会转变为 DW。特征是 $S(k=\pi)$ 有最大值，表示周期为2的密度调制。
  • 配对超流体 (PSF)： 当 $g_x/g_0$ 进一步减小，导致配对跃迁 $P$ 增强时，DW 会转变为 PSF。特征是配对关联函数 $C_2(\Delta j)$ 呈代数衰减，而单粒子关联函数 $C_1(\Delta j)$ 呈指数衰减。
  • 手性超流体 (CSF)： 在 $J_1/g_0G_{000}$ 增加（隧穿增强）时，MI 和 DW 区域会转变为 CSF。这是由几何挫折引起的，特征是电流-电流关联函数 $K_2(\Delta j)$ 呈代数衰减，表明时间反演对称性自发破缺。

相图中的转变线：

  • 实线： 锐利相变。
  • 虚线： BKT 相变（如 DW-PSF 转变）。
  • 黄 × 标记： $J_1=0$ 时由自旋映射得到的精确相变点。

2. 序参量和关联函数行为 (图4b-e)

图4的底部面板展示了在固定 $J_1/g_0G_{000} = 0.075$ 时，各序参量和关联函数随 $g_x/g_0$ 的变化：

  • 图4b (偶宇称序参量 $O_{ ext{even}}$)： 在 MI 相 ($g_x/g_0 \approx 0$ 附近) 具有显著非零值，在其他相中接近零。
  • 图4c (手性-手性关联 $\kappa_2$)： 在 CSF 相 ($g_x/g_0 < -2$) 中显示出非零的准长程关联，而在其他相中为零。
  • 图4d (结构因子 $S(k=\pi)$)： 在 DW 相（例如 $g_x/g_0 \approx -1.5$ 到 $-0.5$ 之间）中表现出尖锐的峰值，表明密度波序的存在。
  • 图4e (配对关联函数 $C_2$)： 在 PSF 相（例如 $g_x/g_0 < -2$）中显示出显著的非零值，与 CSF 区域的 $C_2$ 值接近。论文指出，$C_1(\Delta j)$ 在 CSF 区域也非零但未展示。

3. 相变点和自旋映射验证

  • MI-DW 相变： 在 $J_1/g_0G_{000}=0$ 的极限下（即高势垒 regime），通过将 MI 和 XXZ 模型的 Néel 相能量进行比较，解析计算出 MI-DW 相变点为 $g_x/g_0 \approx -1.26829$（图5）。这个解析结果与 DMRG 模拟的黄×标记完美吻合。
  • DW-PSF 相变： 在高势垒极限下，DW-PSF 转变对应于 XXZ 模型中的 Néel-XY 相变，这是一种 BKT 相变。这个转变是由 $g_x/g_0$ 值控制的。

4. 子空间近似的准确性 (图6)

  • 子空间概率 $P_Q$： 定量评估了将系统限制在每个格点最多2个玻色子的子空间中的准确性。$P_Q = \langle Q angle$，其中 $Q$ 是投影算符，将系统投影到占据数为0或2的格点空间。
  • 结果： 在 DW 和 PSF 区域，$P_Q \approx 1$，表明子空间近似非常准确。然而，随着隧穿强度 $J_1/g_0G_{000}$ 的增加，近似准确性会下降（$P_Q < 1$），意味着单占据格点的数量变得不可忽略。

性能数据

DMRG 模拟的性能数据如下：

  • 系统规模： $L=200$ 个格点，200个原子（单位填充）。
  • 局部玻色子维度： 每个格点最大占据数为2。
  • 最大键维度 ($\chi$)： 500。
  • 最大截断误差： $10^{-6}$。这个低误差表明了计算的高精度，也意味着需要足够的计算资源来实现。
  • 计算区域： 期望值在中心 160 个格点上评估，以减少边缘效应。这表明模拟系统足够大，以确保本体性质的准确性。

这些计算数据和结果共同描绘了该特殊三角梯形晶格中丰富的量子相图，并提供了关于几何挫折与配对跃迁如何协同作用以产生新奇物态的深刻理解。高势垒极限下的解析结果进一步增强了数值模拟的可靠性。

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

代码实现细节

本研究的数值模拟主要通过**密度矩阵重正化群（DMRG）**算法实现，并依赖于专门的量子计算库。虽然论文没有直接提供完整的源代码，但其详细描述了实现的关键要素和DMRG的通用工作流程。

  1. 哈密顿量构建：

    • 首先，需要构建如论文公式(6)-(8)所示的扩展 Bose-Hubbard 哈密顿量 $H = H_{ ext{tun}} + H_{ ext{int}}$。这包括最近邻和次最近邻隧穿项，以及在位、最近邻和配对跃迁相互作用项。
    • 参数 $J_1, J_2, U, V, D, P$ 需根据实验设定的散射长度 $g_0, g_x$ 和 Wannier 函数重叠积分 $G_{000}, G_{001}, G_{011}$ 进行计算。特别地，$D$ 在本次研究中设为0。
    • 在哈密顿量中，玻色子算符 $\hat{a}_j, \hat{a}_j^\dagger$ 需要在每个格点上满足粒子数最大为2的约束（即 $\hat{n}_j \le 2$）。这通常通过在构造局部格点希尔伯特空间时限制其维度来实现。

  2. DMRG 算法核心：

    • 系统初始化： 通常将系统初始化为某种平均填充的基态，或随机初始化并进行塌缩（sweeping）迭代以找到基态。论文指出是单位填充（200个原子在200个格点上）。
    • 张量网络表示： DMRG 将量子态表示为矩阵乘积态（MPS）。哈密顿量通常表示为矩阵乘积算符（MPO）。
    • 迭代过程（Sweeping）： DMRG 算法通过迭代地优化 MPS 的键维度来找到基态。这个过程涉及在系统内部“扫描”，通过左右块的收缩，不断优化中间键上的密度矩阵，并保留最大的本征值以最小化截断误差。
    • 键维度控制： 键维度 $\chi$ 是DMRG的关键参数，它直接影响计算精度和计算资源。本研究使用了最大键维度 $\chi=500$。DMRG通常从较小的 $\chi$ 开始，逐步增加，直到达到最大值或截断误差满足要求。
    • 截断误差： 每次迭代中，在密度矩阵对角化后，会丢弃一部分小本征值。截断误差是这些被丢弃本征值的总和，应保持在足够低的水平（本研究为 $10^{-6}$），以确保结果的精度。

  3. 序参量和关联函数计算：

    • 在 DMRG 找到基态 MPS 后，可以计算各种物理量。这些量通过在 MPS 上收缩相应的算符 MPO 来获得。
    • 单粒子关联函数 $C_1(\Delta j)$： $\langle \hat{a}_j^\dagger \hat{a}_{j+\Delta j} angle$。
    • 配对关联函数 $C_2(\Delta j)$： $\langle \hat{a}_j^\dagger \hat{a}_j^\dagger \hat{a}_{j+\Delta j} \hat{a}_{j+\Delta j} angle$。
    • 结构因子 $S(k)$： $ rac{1}{L^2}\sum_{j,l} \langle \hat{n}_j \hat{n}_l angle e^{ik(j-l)}$。
    • 电流-电流关联函数 $K_2(\Delta j)$： $\langle \hat{k}_j \hat{k}_{j+\Delta j} angle$。
    • 偶宇称序参量 $O_{ ext{even}}^2(\Delta j)$： $\left\langle \prod_{j_0 \le j \le j_0+\Delta j} e^{i\pi \delta \hat{n}_j} ight angle$。
    • 数据提取： 为了避免开放边界条件造成的边缘效应，所有期望值均在系统中间 160 个格点上进行评估。

  4. 相图绘制：

    • 通过系统地改变哈密顿量参数（例如 $g_x/g_0$ 和 $J_1/g_0G_{000}$），重复上述DMRG计算。
    • 根据计算得到的序参量和关联函数的行为（如是否非零、衰减形式、峰值位置），识别不同的量子相。
    • 确定不同相之间的相变边界，并根据关联函数的行为来区分锐利相变和 BKT 相变。

复现指南

要复现本研究的DMRG结果，需要遵循以下步骤。虽然没有公开的代码仓库，但通用DMRG库提供了实现所需的功能。

  1. 安装 Julia 和 ITensors.jl：

    • 首先，确保安装了 Julia 编程语言（建议使用最新稳定版本）。
    • 通过 Julia 包管理器安装 ITensors.jl：using Pkg; Pkg.add("ITensors")。

  2. 定义格点和物理空间：

    • 创建一个 SiteSet，定义每个格点的局部希尔伯特空间。由于每个格点最多允许2个玻色子，可以定义一个 Boson 类型的 SiteSet，并设置 dim=3（对应占据数 0, 1, 2）。或者更精确地，定义一个自定义的 SiteSet 来严格强制 $n_{max}=2$。
    • 设置系统大小 L = 200。

  3. 构建哈密顿量 MPO：

    • 将玻色子哈密顿量的每一项（隧穿、相互作用）转换为 ITensors.jl 的 OpSum 对象。
    • 单粒子隧穿： add!.(ampo, "Ad", j, "A", j+1) 和 add!.(ampo, "Ad", j, "A", j+2) 等，并乘以相应的系数 $-J_1, -J_2$。
    • 相互作用： add!.(ampo, "N", j, "N", j, "N", j, "N", j+1) 等，并乘以相应的系数 $U, V, P$。注意在位项 $ rac{U}{2}\hat{n}(\hat{n}-1)$ 的表示。
    • 将 OpSum 转换为 MPO。

  4. DMRG 参数设置：

    • sweeps 对象： 配置 DMRG 的迭代次数、键维度和截断误差。
    • maxdim： 键维度的序列，例如 [20, 50, 100, 200, 500]，表示逐渐增加键维度。
    • cutoff： 截断误差的序列，例如 [1E-8, 1E-8, 1E-7, 1E-7, 1E-6]，通常在迭代后期降低。
    • noise： 在早期迭代中添加少量噪声有助于避免陷入局部极小值。

  5. 运行 DMRG：

    • 调用 ITensors.jl 的 dmrg(H_mpo, init_mps, sweeps) 函数来计算基态 MPS。
    • init_mps 可以是一个随机初始化的 MPS，或者一个由平均填充得到的 MPS。

  6. 计算期望值和关联函数：

    • 使用 expect(psi, "N") 计算粒子数密度 $\langle \hat{n}_j angle$。
    • 使用 correlation_matrix(psi, "Ad", "A") 计算单粒子关联 $C_1(\Delta j)$。
    • 配对关联 $C_2(\Delta j)$ 和电流-电流关联 $K_2(\Delta j)$ 需要自定义算符并计算四点或更复杂的关联函数。
    • 结构因子 $S(k)$ 可以通过对粒子数密度关联函数的傅里叶变换来计算。
    • 确保只在中心区域（例如格点 21 到 180）提取数据。

  7. 遍历参数空间：

    • 编写一个循环，系统地改变哈密顿量中的 $g_x/g_0$ 和 $J_1/g_0G_{000}$ 参数。
    • 对每个参数组合执行 DMRG 模拟并计算所有相关的序参量和关联函数。
    • 将结果保存到文件中，以便后续进行数据分析和相图绘制。

  8. 数据分析和可视化：

    • 使用 Julia 或 Python（例如 matplotlib）等工具对保存的数据进行分析和可视化，绘制相图。
    • 根据序参量的行为和关联函数的衰减特性来识别不同的量子相和相变点。

  • 主要软件包：

  • 开源代码库：

    • 论文引用了 ITensors.jl 库作为其数值模拟的基础工具，但没有提供专门针对本研究的开源代码仓库链接。这意味着研究人员使用了该库的功能来执行他们的模拟，但并未公开发布他们具体的哈密顿量构建脚本、DMRG运行脚本或数据分析脚本。因此，复现需要基于对论文方法的理解和 ITensors.jl 库的通用功能自行实现。要找到该库的最新开发版本和文档，请访问其 GitHub 仓库。

注意事项：

  • 虽然 DMRG 是一个非常通用的算法，但精确复现特定物理系统的结果通常需要对模型参数、初始化、收敛标准以及如何处理边界条件有深入的理解。
  • 高精度 DMRG 模拟（如本研究中的键维度和截断误差）可能需要大量的计算资源，建议在高性能计算（HPC）集群上运行。

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

关键引用文献

本研究构建于量子模拟和超冷原子物理领域的大量先行工作之上，以下是一些被引用且对理解该工作至关重要的文献类别：

  1. 超冷原子与光学晶格基础 [1-4]： 奠定了超冷原子作为量子模拟平台的基石，介绍了 Bose-Hubbard 模型和光学晶格的基本概念。例如 Jaksch 等人的开创性工作 [1]。
  2. 态相关光学晶格与人造规范势 [5-14]： 这是本研究实现其独特晶格设计的核心技术。这些文献探讨了如何利用原子内部态与光场的耦合，创建出对不同内部态原子产生不同势垒的晶格，以及如何产生人造规范势。特别是关于三足式原子-光耦合方案 [9, 29, 30] 和亚波长势垒 [12, 13] 的工作，为本研究提供了直接的技术基础。
  3. 几何挫折与手性序 [15-17, 27, 51-54]： 阐述了几何挫折在自旋系统和玻色子系统中的重要性，以及它如何导致手性序和时间反演对称性破缺。McCulloch 等人 [15] 和 Moessner & Ramirez [16] 是该领域的经典文献。Dasgupta 等人 [27] 的近期工作也探讨了三角梯形结构中的手性相。
  4. 关联跃迁与配对超流体 [18-21, 25, 26]： 探讨了配对跃迁和密度诱导跃迁等非标准关联项的物理机制及其在形成配对超流体中的作用。Tovmasyan 等人 [18] 和 Eckholt & García-Ripoll [19] 的工作是这方面的代表。Cuzzuol 等人 [26] 提出了区分配对超流体的方法。
  5. DMRG 方法与 Bose-Hubbard 模型 [36, 46, 49, 70-72]： 引用了 White [36] 提出的密度矩阵重正化群算法，以及它在研究一维 Bose-Hubbard 模型及其扩展版本中的应用。Giamarchi [50] 的专著提供了关于一维量子系统更广泛的理论背景。ITensors.jl 库 [37-39] 是具体实现该方法的工具。
  6. XXZ 自旋模型与 Bethe Ansatz [55-57, 75]： 在高势垒极限下，本研究将系统映射到 XXZ 自旋链，并利用其已知的相图和 Bethe ansatz 解进行解析验证。Yang & Yang 的系列工作 [55-57] 是 XXZ 模型解析解的基石，而 Takahashi [75] 的专著提供了 Bethe ansatz 能量的详细推导。
  7. 实验技术 [60-63]： 提及了 Feshbach 共振用于调控散射长度 [60]，以及通过时间飞行（time-of-flight）测量结构因子 [61-63] 和光缔合（photoassociation）测量配对关联 [19, 20] 的技术。
  8. 三体约束与亮孤子 [40-42]： 解释了为何需要三体约束来稳定配对超流体并避免亮孤子的形成。Daley 等人 [40] 和 Bonnes & Wessel [42] 的工作探讨了三体损失和三体约束对玻色子系统的影响。

对这项工作局限性的评论

尽管本研究在理论上取得了显著进展，并在数值模拟中揭示了丰富的量子相图，但仍存在一些潜在的局限性值得讨论：

  1. 三体约束的实验挑战性：

    • 限制： 为了稳定配对超流体并防止亮孤子形成，研究强制每个格点最多只能容纳两个玻色子。这在实验中需要通过三体损失或三体相互作用（通常是强排斥）来实现。虽然论文提到了可以通过三体损失实现 [40-42]，但如何在保持系统相干性的同时精确控制三体损失，以确保其作为有效约束而非简单损耗机制，仍然是一个实验上的挑战。
    • 影响： 如果这个约束未能完美实现，那么相图可能会发生变化，特别是在吸引性在位相互作用区域，可能会形成亮孤子，从而改变 PSF 的性质甚至使其消失。

  2. 密度诱导隧穿 $D=0$ 的假设：

    • 限制： 为了简化模型，本研究将密度诱导隧穿强度 $D$ 设为零，这对应于原子在相同内部态之间的散射长度相等的情况（$g_{11}=g_{22}$ 或 $g_x=0$）。
    • 影响： 实际系统中，$D$ 通常是非零的。论文在未来方向中也提及了允许有限 $D$ 的情况。如果 $D$ 足够强，它会引入额外的关联隧穿机制，可能与配对跃迁和几何挫折发生复杂竞争，从而进一步丰富相图，甚至可能导致共存的超流体相 [19] 或改变现有相的稳定性。

  3. NNN 隧穿与 NN 隧穿比率 $J_2/J_1$ 的固定值：

    • 限制： 论文将 $J_2/J_1$ 比率固定为 $\approx -3.565$，这个值对应于一个实验上可行的参数区域 [12]。
    • 影响： $J_2/J_1$ 的精确值对几何挫折的强度和有效磁通的性质至关重要。探索不同 $J_2/J_1$ 比率下相图的变化，特别是当 $J_2/J_1$ 接近 $-1$ 或其他关键值时，可能会揭示出更复杂的相变行为或新的量子相。

  4. 限制在最低 Bloch 能带（s 波段）：

    • 限制： 论文假设原子仅占据最低 Bloch 能带的 Wannier 函数（s 波段）。
    • 影响： 高能带（如 p 波段）的 Wannier 函数具有不同的空间结构，可能导致更强的配对跃迁（如论文附录A中提及的，$p$ 波段更有利于实现大 $G_{011}$）或其他类型的相互作用。虽然实验上实现 p 波段物理更具挑战性 [43-45]，但理论上探索这些高能带的影响可能会揭示出新的物理现象，例如不同的手性序或更强的配对凝聚。

  5. 绝热近似的有效性：

    • 限制： 系统动力学被限制在暗态流形内，要求原子的质心能量远小于总拉比频率 $\Omega(x)$。
    • 影响： 如果这个条件不严格满足，非绝热效应可能会导致原子泄漏到亮态，增加加热，并破坏系统的相干性，从而影响实验的可行性和理论模型的准确性。

  6. 自旋映射的适用范围：

    • 限制： 将玻色子模型映射到 XXZ 自旋链的近似只在高势垒极限下（即单粒子隧穿可以忽略时）才成立。论文也明确指出，随着隧穿振幅比率 $J_1/g_0G_{000}$ 的增加，近似准确性会下降（图6）。
    • 影响： 这意味着在更广泛的参数区域内（特别是隧穿项显著时），不能简单地依赖自旋映射的解析结果。DMRG 在这些区域仍然是必要的，但缺乏解析验证可能使得对某些相变机制的理解变得更加困难。

  7. 实验实现的复杂性：

    • 限制： 尽管论文认为所需工具“在当前超冷原子技术可及范围内”，但同时实现所有条件仍然复杂。包括精确控制亚波长势垒、精细调节态相关散射长度、以及同时测量多种复杂的关联函数（如手性关联和配对关联）。
    • 影响： 多个实验挑战的叠加可能会使得该理论预测的相难以在短期内完全实现和探测。需要更进一步的技术突破来同时满足所有实验要求。

总而言之，这项工作为在超冷原子中模拟复杂量子现象开辟了新途径，但其理论模型的某些简化和实验实现的固有挑战，也为未来的研究指明了方向。

5. 其他你认为必要的补充

实验可行性与现有技术

本研究提出的方案虽然复杂，但其各个组成部分所依赖的实验技术大多已在超冷原子物理领域实现或正在积极开发中，这使得理论预测的实验实现具备坚实基础。

  1. 亚波长势垒与态相关晶格的创建：

    • 技术基础： 论文引用的 [58, 59] 表明，简化的 $\Lambda$ 型构型暗态原子亚波长势垒已经实现。本研究的三足式耦合方案 [9, 29, 30] 是 $\Lambda$ 型构型的推广，原理相似，因此创建具有亚波长空间结构的态相关光学晶格在技术上是可行的。关键在于精确控制激光的频率、相位和强度，以形成所需空间依赖性的拉比频率 $\Omega_l(x)$。
    • 挑战： 维持 $\Omega_p/\Omega_c \ll 1$ 的条件，以及确保绝热近似，需要高度稳定的激光系统和精细的原子参数控制。

  2. 态相关散射长度的调控：

    • 技术基础： 通过 Feshbach 共振 [60] 技术可以精确调控原子之间的散射长度，包括自旋和内部态相关的散射长度。这意味着可以独立调整 $g_0$ 和 $g_x$（或 $g_{11}, g_{12}, g_{22}$），从而控制相互作用参数 $U, V, D, P$ 的相对强度。
    • 挑战： 寻找同时满足所有散射长度要求且具有合适 Feshbach 共振的应用原子种类。

  3. Wannier 函数的重叠控制：

    • 技术基础： 势垒高度（通过 $\Omega_p/\Omega_c$ 比值控制 [12, 13, 34, 35]）影响 Wannier 函数的形状和重叠，进而影响隧穿振幅 $J_1, J_2$ 和相互作用参数 $U, V, D, P$。这为调控这些参数提供了灵活的手段。
    • 挑战： 精确控制亚波长势垒的宽度和深度，以实现所需的 Wannier 函数重叠积分，尤其是在 $J_2 < 0$ 且 $|J_2| > |J_1|$ 这种特殊比例下。

  4. 可观测量的测量：

    • 密度波序： 可通过测量时间飞行（time-of-flight, TOF）图像的结构因子 [61-63] 来检测。TOF 图像反映了系统的动量分布，其中的峰值可以指示密度波序的存在。
    • 配对关联： 可通过光缔合（photoassociation）测量技术 [19, 20] 进行探测。光缔合可以形成分子，其形成速率与原子对的密度相关。
    • 手性序： 可通过测量动量分布来推断 [14, 27]。手性超流体中的净电流会导致动量分布的特定不对称性或偏移。
    • 挑战： 同时测量多个复杂的关联函数，并且需要高信噪比的探测方法来捕获准长程关联的细微特征。

  5. 三体约束的实现：

    • 技术基础： 可以通过引入三体损失或调节三体相互作用强度来实现每个格点最大占据数为2的约束 [40-42]。
    • 挑战： 确保三体损失在提供约束的同时，不显著影响系统的相干性和寿命。精确控制三体相互作用的强度和符号也是必要的。

未来方向

本研究为后续工作奠定了坚实基础，并指明了多个令人兴奋的未来研究方向：

  1. 考虑有限的密度诱导隧穿 ($D e 0$)：

    • 复杂性： 允许非零的密度诱导隧穿 $D$ 将引入额外的关联跃迁机制，使几何挫折、配对跃迁和密度诱导隧穿之间发生更复杂的竞争。这可能导致新的量子相或改变现有相的稳定性。
    • 潜在物态： 理论上可能出现密度波和配对超流体的共存相 [19]，或者其他具有新颖关联特征的相。

  2. 放松三体约束：

    • 物态： 如果移除每个格点最多2个玻色子的约束，在吸引性在位相互作用下，系统可能会形成亮孤子态或聚集态（clusters）。这与 Bose-Einstein 凝聚体的稳定性密切相关。
    • 动力学： 探索在没有约束的情况下，系统如何演化，以及这些聚集态的动力学行为。

  3. p 波段 Wannier 态的探索：

    • 增强配对跃迁： 论文附录提及，p 波段的 Wannier 函数可能导致更强的配对跃迁强度 $P$，因为 $G_{011}$ 在 p 波段中可能更大。
    • 新奇手性序： p 波段的轨道自由度也可能导致不同类型的手性序或轨道手性超流体 [51, 53]。
    • 挑战： 实验上在光学晶格中激发和控制 p 波段原子通常比 s 波段更具挑战性 [43-45]。

  4. Floquet 工程与多能带模型：

    • 动态控制： 通过 Floquet 工程（周期性驱动） [64]，可以实现对隧穿参数和相互作用参数的动态调控，甚至耦合不同的能带 [65-67]。这可能为实现更多样化的有效哈密顿量和相图提供途径。
    • 多能带物理： 研究多能带玻色子模型中的关联效应，特别是当能带之间存在耦合时，可能会出现新奇的相，如多轨道超流体。

  5. Dublon 和空穴的动力学：

    • 关联输运： 研究双占据态（dublons）和空穴（holes）在晶格中的输运动力学，这在其他超冷原子系统中已有探索 [68, 69]。
    • 激发态物理： 探索系统的激发态性质，而不仅仅是基态。

  6. 引入长程相互作用：

    • 偶极原子： 如果使用具有长程偶极相互作用的原子，可以引入除了最近邻之外的更长程相互作用，这会极大地改变多体物理，可能导致新的超流体相或空间关联模式。
    • 几何挫折与长程相互作用： 长程相互作用与几何挫折的结合是一个非常活跃的研究领域，有望发现更多新奇的量子物态 [14]。

broader Impact/Significance

这项工作具有深远的科学意义和广泛的影响：

  1. 推进量子模拟前沿： 提供了一个在超冷原子中实现复杂哈密顿量的新颖平台，能够同时模拟几何挫折和关联跃迁。这扩展了超冷原子在量子模拟领域的能力，使其能够探索超出标准 Bose-Hubbard 模型范围的物理现象。
  2. 揭示新奇量子物态： 通过理论预测并数值确认了配对超流体和手性超流体的存在，丰富了人们对强关联玻色子系统基态相图的理解。这些物态可能在凝聚态物理中具有对应物，例如高温超导体中的配对机制或手性磁体中的拓扑序。
  3. 连接不同物理领域： 成功地将超冷原子物理、凝聚态物理（挫折系统、关联效应）、量子光学（原子-光耦合）和计算物理（DMRG）等多个领域联系起来，促进了跨学科研究。
  4. 为实验提供指导： 详细的理论预测和参数依赖性为未来的实验设计提供了明确的指导。特别是高势垒极限下的自旋映射，为实验验证提供了关键的解析基准，有助于快速定位实验参数。
  5. 基础科学理解： 通过深入分析几何挫折和关联跃迁的竞争，加深了我们对这些基本物理机制如何塑造物质形态的理解，这对于理解更复杂的材料系统至关重要。
  6. 张量网络方法的发展与应用： 利用 ITensors.jl 库进行大规模 DMRG 模拟，也展示了张量网络方法在处理复杂量子多体系统中的强大能力。

总而言之，这项研究不仅提出了一个创新性的实验方案，更通过严谨的理论分析和数值计算，揭示了强关联量子世界中丰富的物理现象，为未来的量子技术和基础科学探索奠定了重要基础。