Circular Rosenzweig-Porter random matrix ensemble

本文提出并数值验证了一个由 Dyson 布朗运动定义的 Rosenzweig-Porter 随机矩阵系综的酉（圆）模拟，用于刻画周期驱动（Floquet）系统中多体局域化的能级统计与特征态分形性。

要点

想象你有一个巨大的、混乱的舞池，里面挤满了成千上万的舞者。在量子物理的世界里，这些舞者代表了一个复杂系统（例如一组相互作用的原子）的可能状态。

本文旨在理解当音乐以两种不同方式变化时，这些舞者如何移动和相互作用：

1. 静态系统：音乐是稳定、不变的嗡嗡声（就像一间正常、静止的房间）。
2. 弗洛凯（周期性驱动）系统：音乐是有节奏、重复的节拍，每隔几秒就不断改变规则（就像频闪灯或脉冲激光）。

长期以来，物理学家为第一种场景（静态房间）掌握了一本出色的“规则手册”（称为罗森茨韦格 - 波特系综）。这本手册帮助他们预测舞者们是会自由混合（混沌），还是会困在自己小小的角落里（局域化）。

然而，对于第二种场景（脉冲、有节奏的房间），没有人拥有一本好的规则手册。由于当系统受到节奏驱动时，量子力学的规则会发生变化，旧的数学公式并不完全适用。

新想法：圆形舞池

本文的作者问道：“我们能否构建一个适用于有节奏、脉冲系统的旧版规则手册的变体？”

他们创建了一个他们称之为圆形罗森茨韦格 - 波特系综的新模型。

以下是他们构建该模型的过程，使用了一个简单的类比：

• 旧方法（布朗运动）：想象舞者在一条平坦、笔直的线上随机移动。如果你随时间随机地推动他们，他们会以可预测的方式散开。这就是旧模型的运作方式。
• 新方法（圆周运动）：对于有节奏的系统，作者意识到舞者并非在直线上移动，而是在圆上移动。想象舞者在圆形跑道上奔跑。他们的位置由角度（如钟面）而非直线距离来衡量。

他们将新模型定义为专门发生在这个圆上的“随机游走”的结果。他们并非凭空猜测数学公式，而是在计算机上模拟了这一过程，以观察会发生什么。

他们的发现

作者进行了大规模的计算机模拟（涉及多达 1,000 名舞者），以观察他们新的“圆形”模型是否表现得像旧的“直线”模型。他们主要检查了两点：

1. 舞者的间距（能级）
他们观察了舞者之间的间隙。

• 在“混沌”区域：舞者均匀散开，他们之间的间隙遵循一种特定的、复杂的模式（就像拥挤的派对中每个人都在推挤）。
• 在“局域化”区域：舞者聚集在一起，或者以非常可预测、简单的方式彼此远离（就像人们排成一队）。
• 结果：他们的新圆形模型显示出与旧模型完全相同的从“混沌”到“聚集”的转变。行为发生变化的“临界点”出现在相同的位置。

2. 舞者的形状（本征态）
他们观察了单个舞者的影响力是如何“扩散”的。

• 扩散：一个舞者的能量被许多其他人共享。
• 分形：一个舞者处于一种奇怪的中间状态——既扩散了，又未完全扩散。它就像一团具有模糊、自相似形状的云。
• 局域化：一个舞者被困在一个点上。
• 结果：圆形模型重现了这些完全相同的形状。无论舞者是完全混合、部分混合（分形）还是被困住，新模型都与旧模型完美匹配。

结论

本文声称，他们成功构建了著名的罗森茨韦格 - 波特模型的幺正（圆形）版本。

通过将系统视为圆形而非直线，他们创造了一种能够准确描述周期性驱动（脉冲）量子系统行为的工具。就像旧模型是静态系统的“唯象”（描述性）工具一样，这个新的圆形模型充当了受节奏性摇晃或驱动系统的描述性工具。

他们通过证明其新模型的统计“指纹”（能级如何分布以及状态如何成形）与原始、易于理解的模型的指纹无法区分，从而证实了这一点。这为物理学家提供了一种新的、可靠的方法来研究复杂的、有节奏的量子系统，而无需从头开始求解极其困难的方程。

技术摘要

技术摘要：圆 Rosenzweig-Porter 随机矩阵系综

问题陈述
Rosenzweig-Porter (RP) 随机矩阵系综是一个成熟的唯象模型，用于描述静态无序相互作用量子系统中多体局域化 (MBL) 跃迁过程中的能级统计和特征态分形性。在静态系统中，动力学由厄米哈密顿量支配。然而，周期驱动（Floquet）系统由酉 Floquet 算符描述，其特征值位于复平面的单位圆上。尽管 MBL 已推广至 Floquet 系统，但此前尚未构建出能够捕捉周期驱动系统相同唯象学的直接酉（圆）类比 RP 系综。本文探讨了此类酉类比是否可被定义，以及它是否具备原始 RP 系综的关键统计特性。

方法论
作者基于标准 RP 系综可被视为有限时间 Dyson 布朗运动 (DBM) 过程的结果这一观察，提出了圆 Rosenzweig-Porter (CRP) 系综的构建方案。

1. 构建：作者将 CRP 系综定义为含时对称酉矩阵 $S(t)$（维度为 $N$）随机演化的结果。该演化遵循适用于酉矩阵的 Dyson 布朗运动过程：
   $$dS(t) = i\sqrt{dt} U^T(t) X U(t)$$
   其中 $X$ 是在每个无穷小时间步采样的高斯正交系综 (GOE) 矩阵，$U(t)$ 源自分解 $S(t) = U^T(t)U(t)$。
2. 初始条件：该过程从对角矩阵 $S(0) = \text{diag}(e^{i\theta_1(0)}, \dots, e^{i\theta_N(0)})$ 开始，其中初始相位 $\theta_n(0)$ 独立采样自 $[-\pi, \pi)$ 上的均匀分布。这确保了态密度均匀，与通用 Floquet 算符一致。
3. 系综定义：CRP 系综定义为矩阵 $S(t)$ 在特定时间 $t = \epsilon^2 / N^\gamma$ 的状态，其中 $\epsilon$ 为微扰强度，$\gamma$ 控制对角项与非对角项的相对强度。
4. 数值实现：作者使用一种通过 $S(t+dt) = S(t)e^{i\sqrt{dt}X}$ 更新矩阵的算法来数值评估该随机过程。该方法在有限时间步长下保持酉性，并避免了显式的矩阵分解。模拟在矩阵维度 $N \leq 1000$ 下进行，对至少 $10^6$ 个特征值或特征态进行了平均。

主要贡献

• CRP 的定义：本文引入了圆 Rosenzweig-Porter 系综作为标准 RP 系综的酉类比，其通过圆 Dyson 布朗运动过程明确定义。
• 理论映射：作者提供了一个启发式论证，即只要态密度经过适当缩放（具体而言，将 RP 参数 $\epsilon$ 映射为 CRP 语境下的 $\epsilon/\sqrt{2\pi}$），CRP 系综在谱中间的谱统计应与 RP 系综相匹配。
• 数值验证：该研究提供了全面的数值证据，表征了 CRP 系综特征值和特征态的统计特性。

结果
数值调查证实，CRP 系综在不同参数 $\gamma$ 区间内表现出与 RP 系综类似的统计行为：

1. 能级间距统计：连续能级间距的平均比值 $r$ 作为量子混沌的诊断指标。
   • 当 $\gamma < 2$ 时，系统表现出 Wigner-Dyson 统计 ($r \approx 0.530$)，这是混沌系统的特征。
   • 当 $\gamma \geq 2$ 时，系统过渡到泊松统计 ($r \approx 0.386$)，这是可积系统的特征。
   • 相变发生在临界点 $\gamma_c = 2$，临界指数为 $\nu = 1$。数据显示 $r$ 随 $(\gamma - \gamma_c) \ln(N)^{1/\nu}$ 的标度坍缩，与 RP 系综一致。
2. 特征态分形性：在 $S(0)$ 为对角的基底下分析了特征态的逆参与比 (IPR)。
   • 当 $\gamma \leq 1$ 时，特征态是扩展的 ($IPR_2 \sim N^{-1}$)。
   • 当 $\gamma \geq 2$ 时，特征态是局域化的 ($IPR_2 \sim O(1)$)。
   • 在中间区间 $1 < \gamma < 2$ 内，特征态是分形的，标度为 $IPR_2 \sim N^{-(2-\gamma)}$。
3. 特征态形状：特征态振幅 $|\langle m | \psi_n \rangle|^2$ 的分布被发现遵循 Breit-Wigner 统计。数据与理论预测一致：
   $$|\langle m | \psi_n \rangle|^2 \sim \frac{1}{(\theta_n - \theta_m^{(0)})^2 + \Gamma^2}$$
   其中 $\Gamma$ 为展宽宽度。虽然有限尺寸效应和连续近似导致了轻微的定量不匹配，但 $\gamma = 0.5$ 和 $\gamma = 1.5$ 区间之间明显的定性一致性和结构差异已被清晰观测到。

意义与主张
作者主张，圆 Rosenzweig-Porter 系综可作为描述周期驱动（Floquet）系统中多体局域化跃迁过程中观测到的能级统计和特征态分形性的有效唯象模型。通过确立 CRP 系综与 RP 系综共享关键统计特性，该工作为研究驱动量子系统中的非遍历相提供了一种可行的理论工具。

本文谦逊地将此项工作表述为构建更广泛推广（例如包含多重分形性的模型或其他具有多重分形特征态的 Floquet 模型）的“第一步”。作者建议，未来的推广可以通过考虑具有相关增量的随机过程来实现，这类似于针对标准 RP 系综的最新提议。此外，他们指出 CRP 系综在随机量子电路研究中作为非最大随机构建块的潜在用途。