Exact diagonalization study of energy level statistics in harmonically confined interacting bosons

本文采用精确对角化方法证明，虽然在适度相互作用机制下，受谐振势束缚的相互作用玻色子表现出规则（泊松）或弱混沌统计特性，但在强相互作用机制下会向以高斯正交系综（GOE）分布为特征的强混沌行为转变，且旋转会进一步增强这种混沌性。

要点

想象一个在一个巨大的、隐形的碗状容器内，挤满了舞者的舞池。在这个舞池上，有一群完全相同的舞者（玻色子），他们都在试图跟随同一个节奏起舞。这个舞池的形状像一个扁平的煎饼（“准二维平面”），而音乐是持续的、有节奏的嗡鸣声（“谐振势阱”）。

舞者的动作主要受两件事影响：

1. 碗的形状： 碗壁将他们推向中心。这就是“陷阱能量”。
2. 舞者的个人空间： 舞者不喜欢互相碰撞。他们有一种排斥力，就像隐形的泡泡膜一样，当彼此靠得太近时会将对方推开。这就是“相互作用能”。

科学家们想要知道：当这些舞者移动时，他们的模式是有序可预测的，还是混乱随机的？

为了弄清这一点，他们不仅观察了舞蹈本身，还观察了“能级”（即舞者可以采取的特定步法或音符）。他们使用了一种特殊的数学工具箱，来观察这些步阶之间的间隙是随机的，还是遵循严格的规则。

两种主要的场景

研究人员测试了两种不同的舞池“情绪”设置：

1. “冷静”之舞（中等相互作用）

• 设置： 舞者们很有礼貌。与其他力量相比，将他们推开的力量相对较弱。
• 结果（不旋转）： 当碗没有旋转时，舞者的移动方式非常有序、可预测。他们的步阶遵循“泊松分布”（Poisson distribution）。
  • 类比： 想象一群人在排队等公交车。他们以随机的时间间隔站立，但并不在意彼此。有时两人站得很近，有时又离得很远。这里不存在“能级排斥”（即他们不会主动避开彼此）。这是一个规则的、非混沌的系统。
• 结果（旋转）： 如果你开始缓慢旋转这个碗（产生一个单个涡旋），舞者们会变得有点坐立不安。他们开始表现出“弱混沌”的迹象。他们还没有完全陷入随机状态，但也并非完全有序。

2. “狂野”之舞（强相互作用）

• 设置： 舞者们变得非常有攻击性。现在，将他们推开的力量与碗壁的力量一样强大。
• 结果（不旋转）： 突然之间，舞池变得混乱了。步阶不再看起来是随机的，而是看起来像一个复杂的、混沌的系统。
  • 类比： 现在，舞者们在积极地避开彼此。如果一个人迈出一步，其他人会立即调整位置以避免碰撞。这被称为“能级排斥”。步阶的模式现在符合“高斯正交系综（GOE）分布”，这是混沌的数学特征。
• 结果（旋转）： 当你在舞者表现得非常有攻击性时旋转这个碗，混沌会进入过载状态。系统变得强混沌。

转折点：有多少舞者？

研究人员还改变了舞者的数量（12、16 或 20 个）。

• 在**“冷静”*场景中，增加更多舞者实际上使系统变得更加有序*（更像随机的公交车队伍）。
• 在**“狂野”**场景中，增加更多舞者会导致混沌发生波动。有时它变得更混沌，有时又会稍微平复下来，但它总体上保持在混沌区域。

“旋转”因素

论文发现，旋转是终极的混沌放大器。

• 即使在舞者只是表现出中等程度的相互作用时，旋转这个碗也会让他们表现得更加混沌。
• 当舞者已经表现得非常有攻击性时，旋转这个碗会让混沌变得更强。
• 他们甚至测试了快速旋转这个碗（产生 2 或 3 个涡旋）。在这种情况下，无论舞池上有多少舞者，系统都是纯粹混沌的。

他们使用的工具（简化版）

为了测量这一点，科学家们使用了四种不同的“尺子”：

1. 最近邻间距 (NNSD)： 测量一个步阶与紧接着的下一个步阶之间的距离。
2. 间距比例： 比较步阶 A 与 B 的距离，以及 B 与 C 的距离。（这是一个巧妙的技巧，可以避免数学误差）。
3. 长程尺子（Dyson-Mehta 与 能级涨落）： 这些工具观察一段较长步阶范围内的模式，以观察整个舞池是僵硬的还是灵活的。

核心结论

论文得出结论，这些被束缚的原子行为是陷阱（希望有序）与相互作用（创造复杂性）之间的一场拉锯战。

• 弱相互作用 + 不旋转 = 有序（规则）。
• 强相互作用 或 旋转 = 混沌。
• 强相互作用 + 旋转 = 最大程度的混沌。

本质上，这项研究表明，通过简单地改变粒子相互推挤的力度或系统旋转的速度，你可以将一个量子系统从一台可预测的精密钟表切换到一场狂野的混沌风暴。这有助于科学家理解“量子混沌”是如何在现实世界中涌现的，特别是对于像铷原子这样构成的超冷气体。

技术摘要

技术摘要：谐振势约束相互作用玻色子能级统计的精确对角化研究

问题陈述
本文研究了量子多体系统的能量谱统计特性，特别关注准二维平面内受谐振势约束的相互作用玻色子。虽然随机矩阵理论（RMT）提供了一个框架，通过 Bohigas-Giannoni-Schmit (BGS) 和 Berry-Tabor 猜想来区分正则（可积）系统与混沌（不可积）系统，但关于改变二维相互作用强度、粒子数 ($N$) 以及旋转对受陷相互作用玻色子的具体影响尚未得到充分探索。以往的研究通常采用近似方法或侧重于特定极限。本研究旨在通过进行精确对角化研究，来理解相互作用能、陷阱能以及旋转之间的相互作用如何驱动从正则到混沌的谱行为转变，从而填补这一空白。

方法论
作者采用精确对角化方法求解了由 $N = 12, 16,$ 和 $20$个$^{87}\text{Rb}$ 原子组成的许多体哈密顿量。

• 系统模型： 玻色子被约束在一个准二维谐振陷阱中，各向异性参数为 $\lambda_z = \omega_z/\omega_\perp = 4$。相互作用由范围 $\sigma = 0.1 a_\perp$ 的排斥高斯势建模。研究在非旋转（$L_z = 0$）和旋转框架（单涡旋态 $L_z = N$，以及更高角动量 $L_z = 2N, 3N$）下进行分析。
• 相互作用机制： 研究考察了两种不同的机制：
  1. 中等相互作用： 相互作用能相对于陷阱能较小（$g_2 = 0.3669$，对应 $a_{sc} = 1000 a_0$）。
  2. 强相互作用： 相互作用能与陷阱能相当（$g_2 = 3.669$，对应 $a_{sc} = 10000 a_0$）。
• 统计工具： 研究通过短程和长程相关性度量分析了最低 100 个能级：
  • 短程： 最近邻间距分布 (NNSD) $P(s)$ 和连续能级间距比值的分布 $P(r)$。作者强调了 $P(r)$ 的重要性，因为它避免了 $P(s)$ 所需的展开（unfolding）过程，从而最大限度地减少了系统误差。使用 Brody 分布来拟合数据，其参数 $b$ 在 Poisson 分布（$b=0$）与高斯正交系综（GOE, $b=1$）之间进行插值。
  • 长程： Dyson-Mehta $\Delta_3(L)$ 统计量和能级数方差 $\Sigma^2(L)$。

主要贡献与结果

1. 非旋转情况 ($L_z = 0$)：

   • 中等相互作用： 系统表现出正则行为。NNSD 和 $P(r)$ 分布符合 Poisson 分布，表明缺乏能级排斥。Brody 参数 $b$ 随 $N$ 的增加而系统性下降（从 $N=12$ 时的 $0.33$降至$N=20$时的$0.04$），表明在该机制下，随着玻色子数量的增加，谱相关性逐渐减弱。长程统计（$\Delta_3$和$\Sigma^2$）在小能量区间也遵循 Poisson 行为。
   • 强相互作用： 系统向混沌行为转变。分布与 GOE 对齐，标志着强烈的能级排斥。然而，混沌程度受 $N$ 的调制。对于 $N=16$，系统表现出最强的混沌特征（$b=0$，$\langle r \rangle \approx 0.538$）；而对于 $N=12$ 和 $N=20$，系统则回归到弱混沌机制，带有一定的正则性。在强相互作用机制下，$\Delta_3(L)$ 统计量在饱和时的数值显著低于中等相互作用机制，证实了向不可积系统的转变。

2. 旋转情况（单涡旋态 $L_z = N$）：

   • 中等相互作用： 旋转诱导了从正则到弱混沌行为的转变。系统表现出带有正则性的混沌特征（Brody 参数 $b \approx 0.52 - 0.64$）。平均间距比值 $\langle r \rangle$ 接近 $0.5$，处于 Poisson 与 GOE 之间。
   • 强相互作用： 系统表现出与 GOE 统计一致的强混沌行为（$b \approx 1.0, \langle r \rangle \approx 0.55$）。研究发现，旋转相比于相同相互作用机制下的非旋转情况，增强了系统的混沌本质。

3. 更高角动量 ($L_z = 2N, 3N$)：

   • 在强相互作用机制下，具有 $L_z = 2N$ 和 $L_z = 3N$ 的系统表现出强烈的量子混沌特征，其短程（$P(r)$）和长程（$\Delta_3$）统计均与 GOE 分布一致。

意义与主张
作者声称，其研究为如何通过相互作用强度、粒子数和旋转来控制受陷玻色子的谱相关性提供了全面的理解。

• 转变机制： 相互作用能与陷阱能之间的相互作用决定了从正则（Poisson）到混沌（GOE）的行为转变。
• 旋转的作用： 旋转被确定为在中等和强相互作用机制下都能增强混沌行为的因素。
• 方法论的稳健性： 本研究验证了使用连续能级间距比值 $P(r)$ 作为短程相关性分析可靠工具的有效性，因为它避免了与 NNSD 相关的潜在展开过程误差。
• 普适性： 研究结果支持随机矩阵理论在描述量子多体系统（特别是超冷玻色系统）谱相关性方面的普适适用性。

论文结论较为谦逊，指出这些发现推进了对超冷玻色子系统中量子混沌的理解，并指出未来的工作可以扩展到谱形式因子以及具有长程相互作用的扩展 Bose-Hubbard 模型。