来源论文: https://arxiv.org/abs/2606.01116v1 生成时间: Jun 07, 2026 01:00

里德堡原子三角棱柱阵列中的基态相图：DMRG模拟与对称性破缺的多重演生相

0. 执行摘要

强关联量子多体系统中的基态相图和量子临界现象一直是凝聚态物理和量子信息科学的核心研究方向。近年来，基于光镊阵列的里德堡原子系统凭借其高度可编程的空间几何结构、超强的长程范德华相互作用以及里德堡能级固有的“里德堡阻碍（Rydberg blockade）”机制，已成为研究新型量子相、拓扑有序以及非平庸量子熔融过程的重要平台。

目前，学术界对一维里德堡原子链以及二维正方形、三角形晶格中的量子相变已有了较为深入的理解。然而，介于一维和二维之间的准一维多腿梯子（Multi-leg ladders）和棱柱（Prism）几何结构，由于其特殊的空间限制和离散对称性，表现出更为丰富且非平庸的物理图像。本文将深度解析最新发表的学术成果——《Ground-state phase diagram of Rydberg atoms in a triangular-prism array》。该工作系统性地采用**密度矩阵重整化群（DMRG）**算法，研究了排布在三角棱柱光镊阵列中的里德堡原子的基态相图。该系统在空间上具有高度非平凡的 $D_3$ 点群腿交换对称性（Leg-exchange symmetry），这与传统的两腿里德堡梯子（具有 $\mathbb{Z}_2$ 腿交换对称性）有着本质区别。

本研究的核心发现包括：

1. 多样化的密度波相：通过精细调节激光解调-拉比频率比（$\Delta/\Omega$）和里德堡阻碍半径（$R_b/a$），系统展现出由空间平移对称性和 $D_3$ 腿交换对称性自发破缺所主导的多种密度波（DW）相。
2. $\mathbb{Z}_6$ 钟模型物理的涌现：在较小的阻碍半径和高解调区，发现了一种交替双-单里德堡占位的共占密度波相（$\mathbb{Z}_2 \times \mathbb{Z}_3$）。该相通过两次连续的 Berezinskii-Kosterlitz-Thouless (BKT) 相变熔融进入无序相，在低能极限下由 $\mathbb{Z}_6$ 钟模型（Clock model）完美描述。
3. 第一阶 $D_3$ 对称性破缺相变：在较大的阻碍半径下，每个三角形单元仅允许单个里德堡激发，系统通过一阶相变进入完全破缺 $D_3$ 对称性的晶格相。
4. Rung-Trimerized（横向三聚化）密度波相：当强阻碍机制抑制了相邻三角形单元上的双重激发时，随着解调度的增加，系统会自发形成横向三聚化的 $\mathbb{Z}_p^+$ 密度波相。这些相在熔融时的临界行为与一维里德堡链极为相似，表现出包含 Ising 临界线、Potts 临界点、Ashkin-Teller (AT) 临界点以及手征临界线（Chiral critical lines）和浮动相（Floating phases）的复杂临界结构。

本博客将站在量子化学与强关联理论物理研究者的视角，对该工作的理论基础、数值计算细节、关键 benchmark 数据、Julia ITensor 代码复现方案以及该工作的物理局限性进行全方位的深度剖析。

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

1.1 核心科学问题

本工作的核心科学问题在于：非阿贝尔 $D_3$ 腿交换对称性如何与沿棱柱轴向的一维平移对称性破缺发生竞争与耦合，从而在准一维里德堡原子系统中演生出不同于传统一维链和两腿梯子的新奇量子相和相变机制？

在传统的两腿里德堡梯子中，腿交换对称性为 $\mathbb{Z}_2$。而在三角棱柱结构中，横截面为等边三角形（如 图 1(d) 所示），其点群对称性为 $D_3$（包含 3 阶旋转对称性 $\mathbb{Z}_3$ 和 3 个反射对称性 $\mathbb{Z}_2$）。这种更高的离散对称性极大地扩展了对称性自发破缺的模式空间，使得系统能够容纳非共度的浮动相、手征熔融相变以及由 $\mathbb{Z}_6$ 对称性主导的拓扑相变。

1.2 理论基础与哈密顿量

1.2.1 原始里德堡哈密顿量

在实验可实现的等离激元光镊阵列中，三角棱柱上的里德堡原子动力学由以下紧束缚类型的哈密顿量（以 $\hbar=1$ 为单位）描述：

其中：

• $L$ 为三角棱柱的轴向层数（每一层包含 3 个原子，排布在等边三角形的顶点上，共 $3L$ 个格点）。
• $i \in [1, L]$ 表示轴向层数索引，$s \in [1, 3]$ 表示同一层等边三角形上的三个顶点（腿索引）。
• $|g_{i,s}\rangle$ 和 $|r_{i,s}\rangle$ 分别代表格点 $(i,s)$ 处的基态和里德堡激发态。$\hat{n}_{i,s} = |r_{i,s}\rangle\langle r_{i,s}|$ 为里德堡数算符。
• $\Omega$ 是拉比频率（Rabi frequency），$\Delta$ 是激光解调度（Laser detuning）。
• $V_{r,r'} = C_6 / |\mathbf{r} - \mathbf{r'}|^6$ 代表各向同性的长程范德华（van der Waals）排斥相互作用，其中 $C_6$ 为色散系数。定义里德堡阻碍半径为 $R_b = (C_6/\Omega)^{1/6}$。

1.2.2 强阻碍极限下的有效哈密顿量

当阻碍半径 $R_b/a$ 足够大（$a$ 为相邻格点间距），使得单个等边三角形单元内绝对无法容纳多于一个里德堡激发（即约束条件 $\sum_{s=1}^3 \hat{n}_{i,s} \le 1$ 严格成立）时，可以通过投影算符将局部希尔伯特空间维度从 $2^3 = 8$ 压缩为 $1 + 3 = 4$。在这种投影子空间内，系统的有效低能哈密顿量可以写为：

在此有效公式中：

• $\hat{N}_i = \sum_{s=1}^3 \hat{n}_{i,s}$ 为第 $i$ 层三角形上的总里德堡激发数算符。
• $\hat{A}_{\text{eff},i}$ 为有效湮灭算符，其作用满足 $\hat{A}_{\text{eff}}|100\rangle = |000\rangle$（对其他分量类似）。
• $V_l$ 控制层与层之间的总电荷密度调制相互作用。
• $U_l$ 代表层间非对角相互作用，表现为有效反铁磁耦合。其中引入了风味-符号算符（Flavor-sign operator） $\hat{D}_{s,i}$，定义为： $$\hat{D}_{s,i} \hat{n}_{i,s} = -\hat{n}_{i,s}, \quad \hat{D}_{s,i} \hat{n}_{i,s'} = \hat{n}_{i,s'} \quad (s \neq s')$$ 该算符对当前腿上的激发赋予 $-1$ 权重，对另外两腿的激发赋予 $+1$ 权重。这种耦合倾向于使相邻两层占用的里德堡原子处于不同的腿上，从而降低范德华排斥能。

1.3 技术难点与应对策略

1. 长程范德华作用的截断与精确模拟：里德堡相互作用随距离呈 $r^{-6}$ 衰减。在 DMRG 计算中，直接处理真正长程的相互作用会显著增加张量网络的键维度。作者通过将长程交互作用截断在 $l_{\text{max}} = 20$ 个轴向单位，既保证了物理上的精确性（误差 $< 0.1\%$），又维持了矩阵乘积态（MPS）表示的可控性。
2. BKT 相变的准确定位：BKT 相变常伴随极强的系统尺寸修正和对数修正，常规的能隙外推或磁化强度分析法极难准确捕捉相变点。本工作通过引入 ESCO（特征态交叉重叠，Eigenstate Crosscap Overlap） 方法，在周期性边界条件（PBC）下提取 Luttinger 参数 $K$。通过寻找 $K = 1/8$ 和 $K = 1/18$ 的精确相交点，成功锁定了 $\mathbb{Z}_2 \times \mathbb{Z}_3$ 相的熔融边界。
3. 多腿系统中的纠缠熵局域涨落：由于等边三角形结构的内秉几何限制，开边界条件（OBC）下 DMRG 计算的纠缠熵 $S_{\text{vN}}(l)$ 会在层内和层间产生强烈的奇偶震荡。本工作利用一维保角场论（CFT）的经典标度公式并辅以对数修正项，对纠缠熵数据进行了精细拟合，从而准确提取了中心电荷 $c$。

1.4 计算方法细节与诊断算符

1.4.1 保角场论与纠缠熵标度

在 $(1+1)$ 维共形临界点，系统在 OBC 下的冯·诺伊曼纠缠熵服从如下标度行为：

其中 $c$ 为 CFT 中心电荷。对于高度涨落的临界相或手征相，拟合公式需引入对数修正项 $\propto 1/\ln L$。

1.4.2 三聚化密度波序参量

为了探测保留 $D_3$ 对称性但破缺平移对称性的 $\mathbb{Z}_p^+$ 晶格相，定义局域三聚算符：

相应的 $\mathbb{Z}_p^+$ 序参量为：

通过 Binder 累积量（Binder Cumulant）分析其临界交叉行为：

1.4.3 $D_3$ 腿交换对称性破缺算符

为了量化三条腿上密度的不均匀分布，定义复数局域场 $\eta_i$ 以及傅里叶模式算符 $\Psi_q$：

对于完全破缺 $D_3$ 的相，当 ordering 矢分别为 $q = 2\pi/3$ 和 $q = \pi/3$ 时，$\Psi_q$ 的突变行为分别指示了 $D_3$ 相和 $\mathbb{Z}_2 \times D_3$ 相的一阶转变。

2. 关键 Benchmark 体系、计算所得数据及性能数据

本工作通过严谨的数值计算，构建了非常详实的基态相图（如 图 1 所示）。以下将分类详述关键相和临界点的定量计算结果。

2.1 $\mathbb{Z}_2 \times \mathbb{Z}_3$ 钟模型相与双 BKT 熔融

• 体系参数：$R_b/a = 1.14$，采用哈密顿量 (1)，截断长程作用至第 4 邻单元，系统尺寸 $L = 24$ 至 $100$（使用 PBC 以进行 ESCO 测量），以及 $L = 109$ 至 $349$（使用 OBC 进行能隙计算）。
• 物理图像：该相展示出周期为 2 的轴向电荷调制（一个三角形内占有 2 个里德堡原子，相邻三角形内占有 1 个），且同时破缺 $\mathbb{Z}_3$ 旋转对称性，但保留了反射对称性。这产生了六重简并的基态。
• 计算所得性能数据：
  • 能隙交叉外推（图 7(a)）：利用有限尺寸标度形式 $(\Delta/\Omega)_L = A/[\ln(L) + C]^2 + (\Delta/\Omega)_c$，拟合得出：
    • 左边界（无序相至临界相）：$(\Delta/\Omega)_{c1} = 3.1542$
    • 右边界（临界相至有序相）：$(\Delta/\Omega)_{c2} = 3.7449$
  • ESCO 导出的 Luttinger 参数 $K$（图 6, 图 7(b)）：
    • 理论预测，左边界（Rabi 项驱动）对应 $K_c = 1/8 = 0.125$；右边界（$\mathbb{Z}_6$ 各向异性驱动）对应 $K_c = 1/18 \approx 0.0556$。
    • 通过有限尺寸 $1/L$ 外推拟合交叉点，得到极佳的自洽结果：
      • 左边界：$(\Delta/\Omega)_{c1} = 3.1435$（与能隙外推相差仅 $0.3\%$）
      • 右边界：$(\Delta/\Omega)_{c2} = 3.7661$（与能隙外推相差仅 $0.5\%$）
    • 在临界相内部（例如选取 $\Delta/\Omega = 3.4$, $R_b/a = 1.14$），纠缠熵拟合直接给出中心电荷 $c \approx 1.023$（如 图 9 紫色圆圈数据），与 $\mathbb{Z}_6$ 高斯临界理论预期的 $c = 1$ 高度契合。

2.2 $\mathbb{Z}_2^+$ 晶格相与 Ising 临界行为

• 体系参数：选取 $R_b/a = 1.9$ 和 $2.0$ 横切线，系统尺寸 $L = 529, 589, 649$。
• 计算所得性能数据（图 2）：
  • 对于原始哈密顿量 (1)，在 $R_b/a = 1.9$ 处，Binder 累积量 $U_4$ 数据崩塌（Data Collapse）给出：
    • 临界点：$(\Delta/\Omega)_c = 2.5605$
    • 相关长度临界指数：$\nu = 0.981$
  • 对于有效哈密顿量 (2)，在 $R_b/a = 1.9$ 处，崩塌给出：
    • 临界点：$(\Delta/\Omega)_c = 2.5805$
    • 相关长度临界指数：$\nu = 0.980$
  • 在 $R_b/a = 2.0$ 处，原始哈密顿量和有效哈密顿量的 $\nu$ 分别为 $0.968$ 和 $0.964$。
  • 提取的临界纠缠熵中心电荷为 $c \approx 0.505$（如 图 9 红色折线），与标准二维 Ising 普适类预期值（$\nu = 1, c = 0.5$）完全一致。

2.3 $\mathbb{Z}_3^+$ 和 $\mathbb{Z}_4^+$ 晶格相的 CFT 多体临界点

• $\mathbb{Z}_3^+$ 三聚化相（图 3）：
  • 在 $R_b/a = 2.6545$、$\Delta/\Omega = 3.6372$ 处发现三态 Potts 共度临界点。
  • 数据崩塌得出：$\nu = 0.842$（与 3-state Potts 理论值 $\nu = 5/6 \approx 0.833$ 极近），动力学临界指数 $z = 0.993 \approx 1$。
  • 结合对数修正拟合，提取的中心电荷 $c \approx 0.793$，与 Potts 预言值 $c = 0.8$ 误差小于 $1\%$。
• $\mathbb{Z}_4^+$ 共度相（图 4）：
  • 在 $R_b/a = 3.6754$、$\Delta/\Omega = 4.2886$ 处定位出 Ashkin-Teller 临界点。
  • 数据标度给出：$\nu = 0.815$，中心电荷 $c \approx 1.008$（理论预期 $c = 1$），动力学指数 $z = 1.008 \approx 1$。

2.4 $D_3$ 对称性自发破缺的一阶相变证据

在 $R_b/a = 1.6$ 处（图 8(a)）和 $R_b/a = 2.8$ 处（图 8(b)），系统分别经历向 $D_3$ 相和 $\mathbb{Z}_2 \times D_3$ 相的转变：

• 纠缠熵 $S_{\text{vN}}$ 与相应的腿交换序参量 $\Psi_{2\pi/3}$、$\Psi_{\pi/3}$ 表现出极其陡峭且不随系统尺寸 $L$ 减小的跃迁间断行为（Jumps）。
• 能隙（Energy Gap）在转变点处表现出极窄的尖峰，且当尺寸 $L$ 增加时，其极小值很快在 $10^{-3}$ 量级饱和，不满足连续相变的 $L^{-z}$ 标度，明确指示了其一阶相变（Level Crossing）本质。

3. 代码实现细节、复现指南及开源 Repo 链接

本工作使用 Julia 语言开发，基于著名的张量网络开源库 ITensors.jl 进行数值模拟。为了让量子计算与凝聚态物理同行能够快速复现相图中的关键基态数据，下面给出针对有效投影哈密顿量 (2) 的 DMRG 模拟复现指南及核心 Julia 代码实现。

3.1 几何构建与哈密顿量 MPO 编写逻辑

要在 ITensors 中模拟三角棱柱几何，核心在于将一个准一维格点（三腿，共 $3L$ 个物理格点）投影为局部维数为 4 的层算符。定义每一层 $i$ 的局部基组为： |0⟩ (无激发), |1⟩ (第1腿激发), |2⟩ (第2腿激发), |3⟩ (第3腿激发)。

核心复现 Julia 代码

using ITensors
using LinearAlgebra

# 1. 定义自定义层物理 site (4-level system: Empty, Leg1, Leg2, Leg3)
function ITensors.space(::SiteType"RungTriangle")
  return 4
end

# 定义常用的局部物理算符
function ITensors.op(::OpName"N", ::SiteType"RungTriangle")
  # 总激发数算符 N = |1><1| + |2><2| + |3><3|
  return [0.0 0.0 0.0 0.0;
          0.0 1.0 0.0 0.0;
          0.0 0.0 1.0 0.0;
          0.0 0.0 0.0 1.0]
end

function ITensors.op(::OpName"A_eff", ::SiteType"RungTriangle")
  # 湮灭算符 A_eff = |0><1| + |0><2| + |0><3|
  return [0.0 1.0 1.0 1.0;
          0.0 0.0 0.0 0.0;
          0.0 0.0 0.0 0.0;
          0.0 0.0 0.0 0.0]
end

function ITensors.op(::OpName"Adag_eff", ::SiteType"RungTriangle")
  return transpose(op(OpName"A_eff", SiteType"RungTriangle"))
end

# 定义风味符号算符 D_s (s=1, 2, 3)
function ITensors.op(::OpName"D1", ::SiteType"RungTriangle")
  return [0.0 0.0  0.0  0.0;
          0.0 -1.0 0.0  0.0;
          0.0 0.0  1.0  0.0;
          0.0 0.0  0.0  1.0]
end

function ITensors.op(::OpName"D2", ::SiteType"RungTriangle")
  return [0.0 0.0  0.0  0.0;
          0.0 1.0  0.0  0.0;
          0.0 0.0 -1.0  0.0;
          0.0 0.0  0.0  1.0]
end

function ITensors.op(::OpName"D3", ::SiteType"RungTriangle")
  return [0.0 0.0  0.0  0.0;
          0.0 1.0  0.0  0.0;
          0.0 0.0  1.0  0.0;
          0.0 0.0  0.0 -1.0]
end

# 2. 主体计算函数：构建有效哈密顿量并执行 DMRG
function run_rydberg_prism_dmrg(L::Int, Rb_over_a::Float64, detuning::Float64; l_max::Int=20)
  sites = siteinds("RungTriangle", L)
  
  # 相互作用强度参数 V0 (在等边三角形边长 a=1 下，设 Rabi 频率 Omega = 1)
  # Rb/a = V0^(1/6) => V0 = (Rb/a)^6
  V0 = Rb_over_a^6
  Omega = 1.0

  # 初始化 Hamiltonian 的 OpSum
  ampo = OpSum()

  # 单体项: (Omega/2)*(A_eff + A_eff_dag) - Delta * N
  for i in 1:L
    ampo += (Omega / 2.0), "A_eff", i
    ampo += (Omega / 2.0), "Adag_eff", i
    ampo += -detuning, "N", i
  end

  # 双体项: 考虑轴向长程相互作用并截断至 l_max
  for l in 1:l_max
    # 计算有效层间排斥常数 V_l 和非对角反铁磁常数 U_l
    V_l = V0 * ( 1.0 / (4.0 * l^6) + 3.0 / (4.0 * (l^2 + 1.0)^3) )
    U_l = V0 * ( 1.0 / (4.0 * l^6) - 1.0 / (4.0 * (l^2 + 1.0)^3) )

    for i in 1:(L - l)
      # 密度相互作用项
      ampo += V_l, "N", i, "N", i + l
      
      # 风味-反铁磁非对角耦合项 
      for s in 1:3
        D_name = "D$s"
        ampo += U_l, D_name, i, D_name, i + l
      end
    end
  end

  H = MPO(ampo, sites)

  # 3. 设定 DMRG 扫描参数
  sweeps = Sweeps(15)
  maxdim!(sweeps, 50, 100, 200, 400, 600)
  cutoff!(sweeps, 1e-12)

  # 设定初始状态：简单的空态/无里德堡激发
  state = ["0" for n in 1:L]
  psi0 = MPS(sites, state)

  # 运行 DMRG
  energy, psi = dmrg(H, psi0, sweeps)
  
  println("DMRG 收敛基态能量: ", energy)
  return psi, sites
end

# 4. 纠缠熵测量辅助函数
function measure_entanglement_entropy(psi::MPS, b::Int)
  # 对键 b 进行奇异值分解(SVD)以提取奇异值谱
  orthogonalize!(psi, b)
  U, S, V = svd(psi[b], (siteind(psi, b), linkind(psi, b-1)))
  SvN = 0.0
  for n in 1:dim(S, 1)
    λ = S[n, n]^2
    if λ > 1e-16
      SvN -= λ * log(λ)
    end
  end
  return SvN
end

3.2 复现实验设计推荐

想要成功复现文中经典的 Figure 2(a) 数据崩溃和 Figure 9 的中心电荷拟合，建议执行以下步骤：

1. 克隆主流 Julia 包：使用 julia -e 'using Pkg; Pkg.add("ITensors")' 安装最新版稳定库。
2. Ising 临界线验证：固定 $R_b/a = 1.9$，在 $\Delta/\Omega \in [2.50, 2.62]$ 区间内以 $0.005$ 为步长，针对尺寸 $L = 100, 150, 200$ 计算纠缠熵，寻找电荷流最大值点并套用公式 (3) 以验证 $c = 0.5$ 的超精细临界值。
3. ITensors 开源生态链接：
   • ITensors 官方 GitHub 仓库
   • ITensors 官方网站文档与示例

4. 关键引用文献及对该项工作局限性的客观学术评论

4.1 关键参考文献

1. [12] P. Fendley, K. Sengupta, and S. Sachdev, Phys. Rev. B 69, 075106 (2004).
   • 科学贡献：奠定了里德堡一维链中约束哈密顿量的基本临界理论框架，首次详细刻画了 $\mathbb{Z}_p$ 空间对称性破缺中的 Potts 临界和手征过渡机制，是本工作的核心理论支撑。
2. [2] H. Bernien et al., Nature 551, 579 (2017).
   • 科学贡献：实验展示了 51 个可编程光学镊子中里德堡原子的晶格相变，为准一维梯子和棱柱的研究奠定了强有力的实验基础。
3. [72] B.-Y. Tan, Y. Zhang, H.-H. Tu, and Y.-H. Wu, Phys. Rev. Lett. 134, 076501 (2020).
   • 科学贡献：系统开发了利用 MPS 构造的特征态交叉重叠（ESCO）方案来在有限尺寸下提取 Luttinger 参数 $K$，解决了高精度定位高斯临界 BKT 转变点的数值难题。
4. [29] J. Zhang, S. H. Cantu, F. Liu, et al., Nature Communications 16, 712 (2025).
   • 科学贡献：细致刻画了双腿里德堡梯子（Ladders）中的相图、横向 dimerized 状态以及手征熔融线，是本项工作向三腿“棱柱”推广的最直接前驱研究。

4.2 本工作局限性与局外人批判评论

作为一个理性的量子化学和张量网络方向的科学研究者，我们在高度赞赏本工作对 $D_3$ 对称性相变完美刻画的同时，必须客观看待其存在的学术局限性：

1. 一维张量网络的维度限制：
   • 该工作本质上仍然是将三角棱柱“拉直”为一维的长条物理轨道来使用一维 DMRG 求解。这导致其纠缠熵在面对更大的棱柱系统（例如横截面扩展为正六边形棱柱、乃至真二维晶格）时，将无可避免地面临**纠缠熵面积律（Area Law）**惩罚。一维 MPS 的键维度 $D$ 必须随腿数呈指数增长才能维持同等精度，阻碍了该方案直接推广到真正二维具有高点群对称性的系统中。
2. 有效哈密顿量投影的误差敏感性：
   • 公式 (2) 依赖于严格的投影约束 $\sum n_{i,s} \le 1$。但在真实的实验温度（通常为微开尔文级，$\mu\text{K}$）及激光起伏下，高能态（单层容纳双激发的态）依然会有非零占位。这意味着该投影理论相图在与实验直接比对时，高解调区的一阶相变临界行为可能会由于“热软化（thermal softening）”或“激发泄漏”而发生显著偏移。
3. 周期性边界条件（PBC）的尺寸惩罚：
   • 在第 III 部分，为了使用 ESCO 测量 Luttinger 参数 $K$，作者不得不选用 PBC 进行 DMRG 计算。正如 DMRG 领域的常识，由于哈密顿量矩阵带宽增加了一倍，PBC 下的 DMRG 键维度需求呈二次方上涨，系统最大可模拟尺寸被严重限制在 $L = 100$ 左右。这极大限制了 BKT 临界点标度外推在热力学极限（$L \to \infty$）下的准确度上限，产生了约为 $0.5\%$ 的微观不确定性。
4. 色散力各项异性的缺失：
   • 论文假设原子间的范德华相互作用是完全各向同性的 $C_6/r^6$。但是在某些具体实验设计中，由于里德堡原子的强电偶极激发方向受偏振激光的限制，不同朝向（例如沿着棱柱轴向和沿着横截面三角形边方向）的 $C_6$ 实际上具有微小的各向异性（Anisotropy）。该各向异性如何修正当前的 $D_3$ 相图，论文并未予以深入讨论。

5. 补充内容：从 $\mathbb{Z}_6$ 钟模型相变到 Rung-Trimerization 物理图像

为了给致力于里德堡量子计算和量子相变研究的同行提供更完整的学术背景，本节对论文中最为精彩的两处物理内涵进行拓展补充。

5.1 $\mathbb{Z}_6$ 钟模型的双 BKT 相变微观起源

在三角棱柱物理中，$\mathbb{Z}_2 \times \mathbb{Z}_3$ 相的熔融展示出极具吸引力的双连续 BKT 转变。为什么一个纯粹的强关联格点模型会演生出连续的 $\mathbb{Z}_6$ 钟模型物理？

1. 对称性匹配：该有序相具有两部分特征：首先是轴向电荷呈周期为 2 的 $\mathbb{Z}_2$ 平移破缺；其次是横截面上里德堡激发指向特定顶点的 $\mathbb{Z}_3$ 旋转对称性自发破缺。两者相结合，形成了基态的 $2 \times 3 = 6$ 重简并。
2. 低能有效相位标量场：在场论描述中，这 6 个简并基态可以映射为一个相位角 $\theta \in \{ 0, \frac{\pi}{3}, \frac{2\pi}{3}, \pi, \frac{4\pi}{3}, \frac{5\pi}{3} \}$。其低能动力学由如下二维正弦-高斯（Sine-Gordon）类有效场论控制： $$S_{\text{eff}} = \int d\tau dx \left[ \frac{1}{2K} (\partial_\mu \theta)^2 + g_6 \cos(6\theta) + g_R \cos(\tilde{\theta}) \right]$$ 其中 $g_6$ 代表系统的六折晶格各向异性项，$g_R$ 代表拉比（Rabi）相干激发的扰动项。根据重整化群（RG）理论：
   • 当系统从低温（大 $\Delta$）出发时，六折各向异性算符 $\cos(6\theta)$ 的标度维度为 $9K$。若 $K < 1/18$，该算符为相关算符（Relevant），锁定系统进入有能隙的 $\mathbb{Z}_2 \times \mathbb{Z}_3$ 有序相。
   • 随着 $\Delta$ 减小，涨落增强，当 $K > 1/18$ 时，该各向异性算符变为不相关（Irrelevant），系统进入无能隙的卢廷格液体临界相（Gaussian Critical Phase），这也对应着第一个 BKT 转变点。
   • 当 $\Delta$ 进一步减小，拉比相干项 $\cos(\tilde{\theta})$（其标度维度为 $1/(4K)$）在 $K > 1/8$ 时转变为相关算符，从而将临界高斯相彻底摧毁，系统通过第二个 BKT 相变点熔融进入平庸的无序相。这就是 图 5 能隙交叉和 图 6 $K$ 参数分界线的深层物理本质。

5.2 Rung-Trimerization (横向三聚化) 与 Rung-Dimerization 的对比

在量子化学和关联物性中，我们熟悉“二聚化”（Dimerization，如 SSH 模型）。那么“横向三聚化”（Rung-Trimerization）在里德堡棱柱中是如何工作的？

• 物理构造对比（见 图 1(d)）：
  • 两腿梯子的 Dimerization：由于每一横梁（Rung）仅有两个格点，在强相互作用下形成形如 $\frac{1}{\sqrt{2}}(|rg\rangle + |gr\rangle)$ 的二聚纠缠共振态，该态完美保持了横梁的 $\mathbb{Z}_2$ 交换反射对称性。
  • 三角棱柱的 Trimerization：每一个横梁（等边三角形）包含 3 个格点。为了保持 $D_3$ 腿交换对称性，里德堡原子选择在三个格点之间发生高度相干流失，自发形成一个对称化的三聚态： $$|\Psi_{\text{trimer}}\rangle = \frac{1}{\sqrt{3}} (|rgg\rangle + |grg\rangle + |ggr\rangle)$$ 每个原子的平均里德堡占位数为完美的 $1/3$。这种特殊的态不仅没有破缺任何横截面 D3 对称性，而且通过在横梁上的空间离域降低了系统的动力学动能。这一机制是系统在 $R_b/a \ge 2.5$ 区域形成广泛稳定的 $\mathbb{Z}_2^+$、$\mathbb{Z}_3^+$ 和 $\mathbb{Z}_4^+$ 相的物理基石，充分展示了人工量子模拟平台构筑高度可控新型凝聚态物质的独特魅力。

通过对本论文的透彻解析，可以看出利用精确的张量网络方法与先进的冷原子实验平台紧密结合，正不断刷新人类对低维量子关联物理、自发对称性破缺及复杂拓扑熔融相变行为的认知疆界。对于量子化学和强关联理论工作者而言，该三角棱柱中的多体关联机制也为理解复杂大分子内的激子相干传输和多通道电荷转移提供了非常宝贵的低维模型启示。