2025-11-26T03:01:17.371794

Interplay of localization and topology in disordered dimerized array of Rydberg atoms

Prodius, Aramthottil, Zakrzewski

Rydberg tweezer arrays provide a platform for realizing spin-1/2 Hamiltonians with long-range tunneling that decays as a power law with distance. We numerically investigate the effects of positional disorder and dimerization on the properties of excited states in such a one-dimensional system. Our model allows for continuous tuning of both the dimerization pattern and the disorder strength. Within the parameter space constrained by our geometry, we show that both mechanisms lead to a localized phase that does not resemble standard MBL. This phase can be understood as an ensemble of distinct Hilbert space fragmented realizations induced by small inter-spin separations. As dimerization is commonly associated with Symmetry Protected Topological (SPT) physics, we also examine the SPT states across the entire energy spectrum. Despite a partial spin-glass order, we argue that the system hosts an extensive fraction of SPT states.

academic

Interplay of localization and topology in disordered dimerized array of Rydberg atoms

基本信息

论文ID: 2505.07720
标题: Interplay of localization and topology in disordered dimerized array of Rydberg atoms
作者: Maksym Prodius, Adith Sai Aramthottil, Jakub Zakrzewski
单位: Jagiellonian University, Kraków, Poland
分类: cond-mat.dis-nn, cond-mat.quant-gas, quant-ph
发表时间: November 26, 2025
论文链接: https://arxiv.org/abs/2505.07720v2

摘要

本文研究了Rydberg原子光镊阵列中位置无序和二聚化对激发态性质的影响。研究系统为具有长程幂律衰减隧穿的一维自旋-1/2哈密顿量。通过数值模拟发现，在参数空间约束下，无序和二聚化机制导致的局域化相不同于标准的多体局域化(MBL)，而是由小自旋间距诱导的希尔伯特空间碎片化(Hilbert space fragmentation)集合。尽管存在部分自旋玻璃序，系统在整个能谱中仍承载大量对称性保护拓扑(SPT)态。

研究背景与动机

研究问题

核心问题: 在具有长程相互作用的Rydberg原子阵列中，位置无序和二聚化如何共同影响系统的局域化性质和拓扑特性？
关键挑战: 理解非标准多体局域化机制与对称性保护拓扑序的共存关系

研究重要性

实验可行性: Rydberg原子光镊阵列提供了可控的实验平台，可精确调控原子位置和相互作用
理论意义: 探索超越标准MBL范式的新型非遍历性机制
拓扑物理: 研究激发态中SPT序的存在性，这对理解拓扑保护机制至关重要

现有方法局限性

标准MBL理论: 基于局域积分运动(LIOMs)的标准MBL理论不适用于键无序和长程相互作用系统
拓扑态检测: 传统的纠缠谱简并度在存在cat态时无法正确识别拓扑特性
实验局限: 大多数理论研究集中在基态，而激发态的拓扑性质研究不足

研究动机

通过结合实空间重整化群(RSRG-X)方法和数值精确对角化，系统研究位置无序二聚化链的完整能谱性质，揭示希尔伯特空间碎片化、自旋玻璃序和拓扑序的相互作用。

核心贡献

揭示非标准局域化机制: 证明位置无序和二聚化导致的局域化相不是标准MBL，而是由自旋配对引起的希尔伯特空间碎片化集合
发现自旋玻璃序: 在没有ZZ相互作用项的纯XY模型中观察到Edwards-Anderson自旋玻璃序参数的有限值
建立RSRG-X理论框架: 提出适用于二聚化长程相互作用系统的实空间重整化群方案，成功预测本征态结构
提出改进的拓扑探测方法: 引入新的断开熵(disconnected entropy) $S_2^D$ 定义，更有效地识别激发态中的拓扑非平凡态
预测实验可观测量: 提出通过特定初态的时间演化来验证激发态拓扑性质的实验方案，振荡周期分别按 $L^3$ 和 $L^6$ 标度

方法详解

任务定义

研究一维长程XY模型在位置无序和二聚化下的性质：

输入: 无序强度 $W$ 和二聚化参数 $\delta$
输出: 本征态的局域化性质、拓扑特性、自旋玻璃序
约束: 开边界条件，保持 $U(1) \times Z_2^T$ 对称性

模型架构

哈密顿量

$H = J\sum_{i>j} \frac{1}{r_{ij}^3}(s_i^x s_j^x + s_i^y s_j^y)$

其中自旋位置定义为： $r_i = \left(i + \delta(-1)^{i+1} - \frac{1}{2}\right)a + iD_{rb} + W_i$

$\delta \in [-1, 1]$ : 二聚化参数
$W_i \in [-W, W]$ : 均匀分布的位置无序
$D_{rb} = 0.2$ : Rydberg阻塞直径
$a = 1$ : 平均格点间距

参数空间约束

几何约束条件： $W < a \cdot \min\left(\frac{1-\delta}{2}, \frac{\delta+1}{2}\right)$

形成三角形参数空间，如图1所示。

技术创新点

1. RSRG-X方案适配

针对二聚化系统的改进：

强键识别: 优先消除最强耦合 $J_{ij} = J/r_{ij}^3$ 对应的两格点算符
基态选择策略:
- 对 $|\pm\rangle_{ij} = \frac{1}{\sqrt{2}}(|\uparrow\downarrow\rangle \pm |\downarrow\uparrow\rangle)$ 态：生成有效隧穿
- 对 $|Z^\pm\rangle_{ij} = |\uparrow\uparrow\rangle/|\downarrow\downarrow\rangle$ 态：生成化学势项或配对跃迁项
路径依赖: 不同消除路径 $\mathcal{P}$ 对应不同的近似本征态

2. 改进的断开熵定义

提出 $S_2^D$ 定义（图9b）： $S_2^D = S_A^{ent} + S_B^{ent} - S_{A\cup B}^{ent} - S_{A\cap B}^{ent}$

其中：

$A = \{1, ..., L/2 + (L/2 \mod 2)\}$
$B = \{2, ..., L/2 + (L/2 \mod 2), L\}$

优势:

对 $\delta < 0$ : $S_2^D = 2$ 表示拓扑非平凡态， $S_2^D = 1$ 表示体内畴壁态
对 $\delta > 0$ : $S_2^D = 0$ 表示拓扑平凡态
避免了传统 $S_1^D$ 定义的假阳性问题

3. 量子化态分数统计

定义量子化态比例： $P_Q(S^D = k) = \langle N_k/N \rangle$

其中 $N_k$ 是 $S^D \in [k-\epsilon, k+\epsilon]$ （ $\epsilon = 0.05$ ）的态数， $N$ 是希尔伯特空间维度。

实验设置

数值方法

精确对角化(ED): $L \leq 16$
POLFED算法: $L > 16$ ，多项式滤波精确对角化
对称性扇区: 零磁化、偶宇称扇区

系统参数

链长度: $L = 10, 12, 14, 16, 18, 20$
无序实现数: 100-5000（取决于系统尺寸）
能谱采样: 谱中心最多1000个本征值

评价指标

1. 平均能隙比

$\langle r_n \rangle = \left\langle \frac{\min\{\delta_k, \delta_{k+1}\}}{\max\{\delta_k, \delta_{k+1}\}} \right\rangle$

GOE值: $r_{GOE} \approx 0.53$ （遍历相）
Poisson值: $r_{Poi} = 0.389$ （标准局域化）

2. 半链纠缠熵

$S_{ent}(L/2) = -\text{Tr}(\rho \log_2 \rho)$

归一化为Page值 $S_P$ 。

3. Edwards-Anderson自旋玻璃序参数

$m_{EA} = \frac{1}{L^2}\sum_n \sum_{i \neq j} \langle n|s_i^z s_j^z|n\rangle^2$

辅助量： $\chi_{EA} = L \cdot m_{EA}$

对比分析

水平切片: 固定 $\delta = 0$ ，变化 $W$
垂直切片: 固定 $W = 0.1$ ，变化 $\delta$
边界扫描: 沿三角参数空间边界

实验结果

主要结果

1. 相图结构（图2）

在 $(\delta, W)$ 参数空间发现：

遍历区域: 中心低 $W$ 、 $\delta > 0$ 区域， $\langle r_n \rangle \approx 0.53$ ， $S_{ent}/S_P \approx 1$
局域化区域: 大 $W$ $W$ 或 $|\delta|$ $∣ δ ∣$ 区域，特征：
- 上部区域（大 $W$ 小 $|\delta|$ ）: $\langle r_n \rangle \approx 0.4$ （略高于Poisson）
- 两侧区域（大 $|\delta|$ ）: $\langle r_n \rangle < 0.389$ （强烈亚Poisson）

2. 无序驱动的局域化转变（图3a-b）

固定 $\delta = 0$ ：

交叉点: $W_c \approx 0.33$
遍历侧( $W < W_c$ $W < W_{c}$ ):
- $\langle r_n \rangle \to 0.53$ （ $L \geq 14$ ）
- $S_{ent}/S_P$ 超体积标度
局域化侧( $W > W_c$ $W > W_{c}$ ):
- $\langle r_n \rangle \approx 0.4$ （未达Poisson值）
- $S_{ent}/S_P$ 亚体积标度

3. 二聚化驱动的转变（图3c-d）

固定 $W = 0.1$ ：

两个交叉点: $\delta \approx -0.35$ 和 $\delta \approx 0.4$
异常特征: 局域化侧纠缠熵曲线不按系统尺寸单调排序（因切割位置在强/弱键上）
亚Poisson行为: 大 $|\delta|$ 时 $\langle r_n \rangle$ 显著低于0.389

4. 能谱碎片化（图6）

对 $W = 0.5, \delta = 0.0$ ，1000个无序实现的能量密度：

整体分布: 类Lorentzian而非Gaussian
个别实现:
- 某些呈Gaussian（小位移或集体平移）
- 某些显示2-3个"岛"（自旋配对导致的碎片化）
物理机制: 不同无序实现产生不同的希尔伯特空间碎片化模式

消融实验

1. 纠缠熵和能隙比分布（图5）

三个典型点的详细分布：

点A（ $W=0.1, \delta=-0.8$ ，左下角）:

能隙比: 强烈的能级吸引， $r_n \approx 0.45$ 处共振峰
纠缠熵: 峰值在 $S_{ent} = 1, 2$ （整数值）
解释: cat态和边缘-体交换态

点B（ $W=0.5, \delta=0.0$ ，上部）:

能隙比: 小间隙处有凹陷（残余能级排斥）
纠缠熵: 多峰结构， $S_{ent} \gtrsim 1, 2$
解释: 冻结自旋畴和最大关联子系统

点C（ $W=0.1, \delta=0.8$ ，右下角）:

能隙比: 类似点A的能级吸引和共振峰
纠缠熵: 主峰在 $S_{ent} = 1, 2$
解释: 平凡拓扑态的二聚体交换

2. 自旋玻璃序分析（图7）

垂直切片（ $\delta=0$ $δ = 0$ ）:
- $m_{EA}$ 随 $W$ 增长
- $\chi_{EA}$ 在 $W \approx 0.37$ 交叉
- 外推表明 $m_{EA}(L\to\infty) \neq 0$
水平切片（ $W=0.1$ $W = 0.1$ ）:
- 两个 $\chi_{EA}$ 交叉点对应二聚化转变
- 大 $|\delta|$ 区域 $m_{EA}$ 有限
标度行为: 局域化区域 $\chi_{EA} \propto L$ ，遍历区域 $\chi_{EA} \sim 1/L$

3. 拓扑态分数（图11, 13）

$\delta = -0.8$ （非平凡相）:

$P_Q(S_2^D = 2)$ : 从 $L=10$ 的~~15%降至 $L=18$ 的~~10%
$P_Q(S_2^D = 1)$ : 从~~20%降至~~15%
外推（图12）: 有限尺寸拟合显示非零渐近值

$\delta = 0.8$ （平凡相）:

$P_Q(S_2^D = 0)$ : 主导，约70-80%
$P_Q(S_2^D = 1,2)$ : 可忽略

边界扫描（图13）:

左边界（负 $\delta$ ）: 非平凡态分数在强二聚化时显著
右边界（正 $\delta$ ）: 几乎全为平凡态

案例分析

时间动力学验证（图14-15）

案例1: 边缘振荡（图14）

初态: $|\Psi\rangle = |\uparrow\rangle_1 |+\rangle_{2,3} |-\rangle_{4,5} \cdots |\downarrow\rangle_L$
物理图像: 两个近似本征态 $|\pm\rangle_{1,L}$ 的叠加
观测量: 边缘磁化 $s_1^z, s_L^z$
结果:
- 周期 $T = 2\pi/(E_1 - E_2)$
- RSRG-X预测与精确演化一致
- 无序平均: $T \propto L^3$ （拟合证实）

案例2: 边缘-体交换（图15）

初态: $|\uparrow\rangle_1 \otimes |-\rangle_{2,3} \otimes \cdots \otimes |Z^-\rangle_{L/2,L/2+1} \otimes \cdots \otimes |\uparrow\rangle_L$
物理图像: 边缘 $|Z^+\rangle_{1,L}$ 与体内 $|Z^-\rangle_{L/2,L/2+1}$ 交换
结果:
- 振荡周期 $T \propto L^6$ （对数拟合: $\ln T = 6.17\ln L - 0.35$ ）
- 能量尺度: $\Delta E \sim J_{1,L/2+1}J_{L/2+1,L}/J_{L/2,L/2+1}$

实验发现

非标准局域化: 局域化相的能隙比行为（亚Poisson或略超Poisson）与标准MBL显著不同
希尔伯特空间碎片化: 能谱的"岛"结构表明系统是碎片化哈密顿量的集合
拓扑-玻璃共存: 尽管存在自旋玻璃序，仍有大量SPT态（ $\sim 25-30\%$ 在 $\delta=-0.8$ ）
整数纠缠熵: 局域化区域本征态纠缠熵集中在整数值，反映cat态和冻结畴结构
RSRG-X有效性: 重整化群预测的态结构和动力学与精确计算高度吻合

结论与讨论

主要结论

局域化机制: 位置无序和二聚化在Rydberg原子阵列中产生的局域化不是标准MBL，而是希尔伯特空间碎片化的集合体现
相图结构: 在 $(\delta, W)$ 参数空间中，中心遍历区域被局域化区域包围，形成三角形相图
自旋玻璃序: 纯XY模型（无ZZ项）在局域化区域展现有限的Edwards-Anderson序参数
拓扑态存在性: 尽管部分能谱呈玻璃态，在 $\delta < 0$ 区域仍有约25-30%的激发态保持拓扑非平凡特性
RSRG-X适用性: 实空间重整化群方法成功预测本征态结构和时间动力学，证明其在长程相互作用系统的有效性

局限性

1. 系统尺寸限制

最大链长: $L = 20$ （ED）， $L = 18$ （POLFED全谱）
外推不确定性: 拓扑态分数的热力学极限行为需要更大尺寸验证
有限尺寸效应: 图3c中纠缠熵曲线的非单调性反映切割位置依赖

2. 参数空间约束

几何约束: 式(3)限制了可探索的 $W$ - $\delta$ 区域
固定参数: $a=1$ ， $D_{rb}=0.2$ ，不同参数可能改变相图
一维限制: 高维系统可能有不同物理

3. 数值方法局限

POLFED精度: 对于大系统只能获取部分谱
无序平均: 即使5000个实现，统计误差仍存在
对称性扇区: 仅研究零磁化偶宇称扇区

4. 理论理解

能隙比共振峰: $r_n \approx 0.45$ 处的峰值缺乏理论解释
玻璃-拓扑界面: 两种序如何在同一能谱共存的微观机制不清楚
非局域积分运动: 具体形式未明确构造

未来方向

1. 实验验证

Rydberg平台: 利用光镊阵列实现可控位置无序
动力学观测: 测量图14-15预测的边缘振荡（ $T \propto L^3, L^6$ ）
拓扑探测: 通过多体干涉测量断开熵

2. 理论扩展

二维系统: 研究平面阵列的拓扑相变
其他相互作用: 考虑van der Waals项（附录B显示 $\Delta \sim 0.1$ 为临界点）
有限温度: 研究热激发对拓扑序的影响

3. 数值改进

更大系统: 开发tensor network方法处理 $L > 20$
全谱分析: 改进算法获取完整能谱
其他对称性扇区: 验证四重简并性

4. 物理机制

积分运动构造: 显式构造非局域积分运动
玻璃态微观结构: 理解自旋玻璃序的起源
共振峰解释: 发展理论解释能隙比的共振特征

深度评价

优点

1. 方法创新性

RSRG-X适配: 首次将实空间重整化群系统应用于二聚化长程系统，基态选择策略（ $|\pm\rangle$ vs $|Z^\pm\rangle$ ）体现深刻物理洞察
改进拓扑探测: $S_2^D$ 定义巧妙解决cat态干扰问题，分区设计针对性强
多方法结合: ED、POLFED、RSRG-X、统计分析形成完整工具链

2. 实验充分性

参数空间全覆盖: 图2的三角形扫描、图3的正交切片、图4/13的边界分析提供全面图景
多尺度验证: $L = 10-20$ 的系统尺寸标度分析严谨
统计可靠性: 100-5000个无序实现保证统计显著性
多角度刻画: 能隙比、纠缠熵、自旋玻璃序参数、断开熵、时间动力学相互印证

3. 结果说服力

定量预测: 时间动力学 $T \propto L^3, L^6$ 的精确验证（图14-15）
理论自洽: RSRG-X预测与数值结果高度一致
新物理发现:
- 亚Poisson能隙比（图3d，5a,c）
- 整数纠缠熵分布（图5d-f）
- 能谱碎片化（图6）
- 无ZZ项的自旋玻璃序（图7）

4. 写作清晰度

逻辑结构: 从相图→统计性质→玻璃序→RG理论→拓扑分析→动力学验证，层层递进
图表质量: 22个主图+附录图，信息密度高且易读
物理图像: RSRG-X方案的图解（图8）、分区定义（图9）清晰直观
文献综述: 引言全面回顾MBL、碎片化、SPT等领域

不足

1. 方法局限性

RSRG-X近似: 路径依赖性未充分讨论，不同路径的误差估计缺失
外推不确定性: 图12的线性拟合在 $L=10-18$ 范围，外推到 $L\to\infty$ 可能有系统偏差
对称性扇区: 仅分析零磁化偶宇称，四重简并验证成本高（第6节提及但未系统执行）

2. 实验设置缺陷

参数选择: $D_{rb}=0.2$ 的物理动机不足，敏感性分析缺失
能谱采样: "谱中心最多1000个本征值"可能遗漏边缘态信息
无序分布: 仅考虑均匀分布 $W_i \in [-W, W]$ ，其他分布（如Gaussian）未探索

3. 分析不足

共振峰: 图5a,c中 $r_n \approx 0.45$ 的共振峰"目前无解释"（第3节末），削弱理论完整性
玻璃-拓扑界限: 哪些态是玻璃态、哪些是SPT态？两者如何区分？（第6节讨论模糊）
ZZ项影响: 附录B显示 $\Delta=0.1$ 附近质变，但主文未深入讨论物理机制
有限尺寸趋势: 图11c插图显示 $P_Q(S_2^D=2)$ 随 $L$ 递减，但未明确外推结论

4. 表述问题

"extensive fraction": 摘要和结论中声称"广延分数的SPT态"，但图12外推并不完全支持（ $L\to\infty$ 时可能趋零）
"does not resemble standard MBL": 与标准MBL的定量对比不足，仅定性描述
术语混用: "localized phase"有时指MBL，有时指碎片化，需更明确区分

影响力

1. 对领域的贡献

理论: 建立键无序+长程相互作用系统的RSRG-X框架，可推广到其他平台（如偶极分子、离子阱）
实验: 为Rydberg原子实验提供明确预测（动力学振荡周期），可直接验证
概念: 挑战"MBL需要位置无序"的传统观念，拓展非遍历性机制的认知

2. 实用价值

量子模拟: 指导Rydberg平台设计拓扑量子态
拓扑探测: $S_2^D$ 方法可应用于其他开链SPT系统
误差诊断: 理解位置误差对量子计算的影响

3. 可复现性

代码: 未提供（常见问题），但方法描述详细
参数: 所有超参数明确（ $a=1$ , $D_{rb}=0.2$ , $\epsilon=0.05$ 等）
数据: 图表数据点清晰，误差棒标注（如图7插图）

4. 局限性

计算成本: POLFED算法对一般研究者门槛高
实验难度: $L=18$ 的Rydberg阵列技术要求高
理论gap: 缺乏解析结果，依赖数值

适用场景

1. 直接适用

Rydberg原子阵列: 论文设计的主要平台， $1/r^3$ 偶极相互作用天然实现
极性分子: 类似长程相互作用，可调控位置和二聚化
离子阱: 可编程长程Ising/XY模型

2. 需修改

短程模型: RSRG-X方案需调整（原始Fisher方法64更适用）
高维系统: 碎片化模式更复杂，拓扑不变量需重新定义
有限温度: 当前零温结果，热激发影响需另行研究

3. 不适用

标准MBL系统: 本文方法针对键无序，位置无序需不同处理
体能隙闭合系统: RSRG-X要求能隙分布宽，金属相不适用
强相互作用极限: $U \gg J$ 时需超交换微扰论

潜在扩展

机器学习: 用神经网络学习无序实现→碎片化模式映射
量子控制: 利用拓扑态的鲁棒性设计量子存储
多体疤痕: 研究碎片化与量子疤痕的关联（参考文献18-22）
开放系统: 加入耗散，研究拓扑序的稳定性

参考文献（精选关键文献）

MBL基础: Nandkishore & Huse, Annu. Rev. Condens. Matter Phys. 6, 15 (2015) - 综述
希尔伯特空间碎片化: Rakovszky et al., Phys. Rev. B 101, 125126 (2020) - 统计局域化
SPT与MBL: Vasseur et al., Phys. Rev. B 93, 134207 (2016) - 粒子-空穴对称性
Rydberg实验: de Léséleuc et al., Science 365, 775 (2019) - SPT相观测
RSRG-X: Pekker et al., Phys. Rev. X 4, 011052 (2014) - 原始方法
键无序: Aramthottil et al., Phys. Rev. Lett. 133, 196302 (2024) - 同组前期工作
断开熵: Zeng & Zhou, Europhys. Lett. 113, 56001 (2016) - 拓扑不变量
自旋玻璃: Parisi, Phys. Rev. Lett. 50, 1946 (1983) - Edwards-Anderson参数

总评: 这是一篇高质量的理论物理论文，在方法创新（RSRG-X适配、改进断开熵）、数值严谨性（多尺度多指标）和物理洞察（碎片化机制、玻璃-拓扑共存）方面表现突出。主要贡献是揭示了Rydberg原子阵列中非标准局域化机制，并提出可实验验证的预测。不足之处在于有限尺寸外推的不确定性和部分现象（如能隙比共振峰）缺乏理论解释。论文对凝聚态理论（非遍历性、拓扑物态）和冷原子实验领域均有重要影响，为理解长程相互作用无序系统提供了新范式。建议未来工作聚焦于更大系统尺寸模拟、实验验证和理论gap的填补。