Exact SU($N$) symmetry resolution of bosonic Hubbard models: general boson-fermion mappings and flat-band singlet-cluster states

技术摘要：玻色型哈伯德模型中 SU(N) 对称性的精确解析

问题陈述
本文解决了求解 SU(N) 对称玻色型哈伯德模型时的数值挑战，特别是在希尔伯特空间维度随 $N^{N_b}$（其中 $N$ 为味数，$N_b$ 为粒子数）增长的机制下。虽然 SU(N) 费米子模型已通过适配于不可约表示（irrep，如半标准杨表）的基组得到了成功处理，但玻色系统在很大程度上仍受限于阿贝尔 U(1)$^N$ 子群或小 $N$ 近似。这种局限性非常显著，因为近期在超冷极性分子和类碱土金属原子领域取得的实验进展，使得具有吸引或排斥相互作用的大 $N$ SU(N) 玻色系统在实验上变得具有相关性。理解这些玻色系统的基态物理学仍存在关键的理论空白，尤其是在标准 Eisenberg-Lieb 定理（预测铁磁性）可能不再适用的受挫晶格中。

方法论
作者利用作用于轨道自由度的酉群 $U(L)$ 与作用于颜色/味自由度的 $U(N)$ 之间的 Howe 对偶性来进行精确对角化。

• 颜色因子化： $N_b$ 个玻色子的希尔伯特空间被分解为不可约表示张量积的直和：$\mathcal{H}_b = \bigoplus_{|\alpha|=N_b} V^{(L)}_\alpha \otimes V^{(N)}_\alpha$。至关重要的是，对于玻色子，轨道不可约表示 $\alpha$ 与颜色不可约表示 $\alpha$ 相同（而对于费米子，两者互为转置）。
• 基组构建： 该方法利用填充有位点标签的形状为 $\alpha$ 的半标准杨表（SSYT）来构成轨道空间 $V^{(L)}_\alpha$ 的正交归一基。哈密顿量生成元 $E_{i,j}$ 在该基组上的作用仅取决于轨道形状 $\alpha$，这使得整个 SU(N) 问题可以简化为一系列较小的、不含 SU(N) 简并性的轨道问题 $H_\alpha$。
• 约束条件： 该框架能够自然地容纳软核和硬核约束（例如，每个位点最多有 $n$ 个玻色子），方法是简单地丢弃位点标签出现次数超过 $n$ 次的杨表。

主要贡献与结果

1. 一般结构性质：
   本文推导了任何由 $U(L)$ 生成元构建的 SU(N) 不变量哈密顿量的三个基本性质：

   • N-无关性： 在特定轨道扇区 $\alpha$ 内，能谱与 $N$ 无关（只要 $N \ge \ell(\alpha)$），且与粒子统计性质无关。只有颜色简并因子取决于 $N$。
   • 玻色-费米对应关系： 玻色子与费米子可容许扇区之间存在严格的包含关系。具体而言，如果 $L \le N$，则每个费米子扇区都是玻色子扇区；如果 $N_b \le N$，则每个玻色子扇区都是费米子扇区。
   • 配分函数恒等式： 玻色子（$Z_b$）与费米子（$Z_f$）的正则配分函数满足精确恒等式 $Z_b(-N) = (-1)^{N_b} Z_f(N)$。这意味着 SU(N) 不变算符的热力学平均值遵循 $\langle O \rangle_b(-N) = \langle O \rangle_f(N)$。

2. 少位点环的应用（大 N 极限）：
   将 $N \to -N$ 对应关系应用于固定填充度 $\nu = N_b/N$ 的环状吸引玻色子系统：

   • 零温： 基态是一个 SU(N) 铁磁体（扇区 $(N_b)$），其能量与无自旋玻色子相同。
   • 有限温度： 本文推导了从无序相到自陷相（self-trapped phase）转变的相图。对于 $L=2$，该转变是连续的，临界耦合为 $\beta_c U N_b = -4/(2+\nu)$。对于 $L \ge 3$，该转变变为一阶转变。通过设置 $\nu = -1$，可以精确恢复费米子系统的结果（文献 [41]）。

3. 受挫跳跃与平带物理：
   本文研究了 Eisenberg-Lieb 定理不适用的系统（受挫跳跃），特别是锯齿链（sawtooth chain）和卡格梅（kagome）晶格。

   • 平带定理： 对于 $U > 0$，$N_b \le D_{flat}$（平带简并度）且 $N_b \le N$，证明其基态是一个颜色反对称的 SU(N) 单态（扇区 $(1^{N_b})$）。这是费米子平带铁磁性的玻色对应物。
   • 单态簇态（Singlet-Cluster States）： 在有限密度下，作者构造了作为“单态簇”乘积的精确基态。这些簇是通过将紧凑局部态（CLS）填充进平带，并使各簇之间不共享位点而形成的。
     • 锯齿链： 在临界密度 $\rho_c(N) = N/[2(N+1)]$ 以下，存在精确的单态基态。超过此密度后，基态倾向于使用反对称不可约表示以抑制双占据。
     • 卡格梅晶格： 该问题映射为在三角形晶格上填充最多 $N$ 个六边形的连通簇的堆积问题。本文提供了 $N=1$ 到 $4$的临界密度下界（例如，$\rho_{c,2} \ge 2/15$）。
   • 能量竞争： 在整数填充和强耦合下，系统回归为 SU(N) 铁磁体（Mott 限）。在中间耦合和非整数填充下，基态位于中间不可约表示中，代表了动能与相互作用能之间的竞争。

意义与主张
本文声称提供了第一个用于解析玻色型哈伯德模型全 SU(N) 对称性的精确数值工具。其主要意义在于：

• 统一玻色子与费米子： 通过 $N \to -N$ 对应关系建立了玻色子与费米子能谱之间的严谨“字典”，使得一类统计的结果可以为另一类提供参考。
• 计算效率： 证明了与标准的 U(1)$^N$ 方法相比，解析 SU(N) 对称性可以大幅降低颜色单态扇区的希尔伯特空间维度（如表 I 所示，降幅达数个数量级）。
• 受挫系统中的新物理： 确定了“平带单态簇态”是玻色子平带铁磁性的对应物。本文断言，虽然受挫平带会导致无自旋玻色子的定域晶体和费米子的铁磁性，但 SU(N) 玻色子展现出一种独特的颜色反对称单态簇相。

作者指出，目前的数值结果受限于精确对角化方法以及玻色子希尔伯特空间固有的规模。本文谦虚地建议，将经过 SU(N) 解析的密度矩阵重整化群（DMRG）算法扩展到玻色子领域，是获取更大一维系统的必要步骤，并且对于更大的受挫系统，可能需要使用二维张量网络。