来源论文: https://arxiv.org/abs/2606.06605v1 生成时间: Jun 08, 2026 12:32

一维自旋-声子耦合系统中的反常受限与解禁闭量子临界性：张量网络与双正弦-Gordon场论深度解析

0. 执行摘要

传统的兰道-金兹堡-威尔逊（Landau-Ginzburg-Wilson, LGW）对称破缺范式认为，两个破缺不同群对称性的无序-有序相之间的转变必须是相继的一阶相变，或者需要精细调节。然而，解禁闭量子临界点（Deconfined Quantum Critical Points, DQCP） 的提出打破了这一教条，允许在两个截然不同的对称破缺相之间存在连续的量子相变。尽管二维 DQCP 的研究已经历了二十余年，但真实材料中不可避免的晶格畸变（即声子自由度）如何影响临界点附近的物理，仍然是凝聚态物理界的核心前沿课题。

本研究针对一个具有代表性的一维自旋-声子耦合模型，通过高精度无限系统密度矩阵重整化群（iDMRG）以及时变变分原理（TDVP）等张量网络数值模拟，系统性地解答了“晶格涨落如何重整化并破坏解禁闭临界性”的关键科学问题。研究表明：

1. 当动态 Einstein 声子频率 $\omega_0$ 大于临界频率 $\omega_c \approx 0.45J_1$ 时，一维 DQCP 依然能够保持稳定，系统展现出完整的 $U(1)$ 涌现对称性。
2. 当声子频率低于 $\omega_c$ 时，自旋-声子耦合通过降低 Luttinger 参数 $K$，使得原本不可重整的相关项转化为重整化群流中的相关算符。这导致临界线终止于一个属于 四态 Potts（4-state Potts）共性类 的关键终点，并在此之下转化为强一阶相变。
3. 通过分析双正弦-Gordon（double sine-Gordon, DSG）有效场论，本工作澄清了这一一阶相变起源于双重余弦项产生的势垒，这与经典的 Ashkin-Teller 模型具有深刻的对偶性，但其势垒预因子正负号相反。
4. 工作定量计算了声子谱函数 $A(q, \omega)$，在连续临界点处观测到显著的 Kohn 异常（Kohn anomaly），并给出了其普适的幂律标度行为 $A(q=\pi, \omega) \sim \omega^{-2K}$，为中子散射和共振非弹性X射线散射（RIXS）等实验手段探测自旋-声子解禁闭临界性提供了直接的定量指引。

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

1.1 核心科学问题

本研究的核心问题是：动态晶格涨落（声子）是否会破坏一维解禁闭量子临界性？如果破坏，其演化的场论机制和普适性类是什么？

在一维自旋-1/2链中，无能隙的 Luttinger 液体相在晶格畸变下极易发生不稳定性，这被称为 自旋-Peierls 不稳定性（Spin-Peierls Instability）。晶格通过对称破缺（二聚化）降低系统能量，形成有能隙的价键固体（VBS）相。类似地，一维解禁闭量子临界点（在反铁磁 Néel 相与 VBS 相之间）本身是一个无能隙点，因而理论上极易受到晶格自由度的干扰。此前的静态晶格理论预测，静态晶格畸变将直接使 DQCP 转化为一阶相变，但真实的固体材料中，声子是具有动态能量尺度的量子涨落。本工作致力于从定量和定性两个维度，精确给出动态声子对一维 DQCP 稳定性的完整物理图像。

1.2 理论基础与模型构建

1.2.1 微观哈密顿量

研究考虑了一维耦合自旋-晶格系统，其总哈密顿量定义为：

其中，自旋部分 $H_s$ 是典型的具有强易轴各向异性（$\Delta > 1$）和次近邻相互作用的 $J_1-J_2$ 自旋-1/2 链：

在无声子耦合的极限下（$J_2=0$），该系统表现出 Néel 反铁磁序。随着 frustration 调控参数 $J_2/J_1$ 的增大，基态过渡到二聚化的价键固体（VBS）相，其中自旋单态交替排布。对于 $\Delta = 4$，相变点位于 $J_{2,c}/J_1 \approx 0.3923$，表现为无二聚化扰动的连续 DQCP。

声子部分 $H_p$ 采用无耗散、无耦合的 Einstein 声子模型描述：

其中 $u_j$ 是第 $j$ 个格点的位移。自旋-声子耦合项 $H_{sp}$ 起源于近邻自旋交换作用对键长的依赖性（即二阶展开的线性项）：

引入产生消灭算符 $a_j, a_j^\dagger$ 后，耦合项写为：

其中 $g = J_1 \tilde{g} / \sqrt{2m}$。为了简便计算，自旋相互作用被各向同性地推广为 $\mathbf{S}_j \cdot \mathbf{S}_{j+1}$（各向同性自旋-声子耦合），这不改变物理本质，但提升了数值收敛性。

1.2.2 玻色化与连续场论

利用 Jordan-Wigner 变换和玻色化（Bosonization）技术，自旋链的低能物理可描述为 compact 玻色场 $\phi(x) \in [0, 2\pi]$ 的正弦-Gordon（sine-Gordon）作用量：

这里：

• $K$ 为 Luttinger 参数。对于纯自旋链，若无二聚化，相变发生在 $\mu = 0$ 处，此时拥有涌现的 $U(1)$ 对称性（表现为 $\phi \to \phi + \text{const.}$ 的低能对称性）。
• Néel 序和 VBS 序参数分别对应于 $M_j^z \sim \cos(2\phi)$ 和 $\Psi_j \sim \sin(2\phi)$。根据自旋对称性，$\cos(4\phi) \sim \Psi^2 - M_z^2$。当 $\mu > 0$ 时稳定 Néel 相，$\mu < 0$ 时稳定 VBS 相。其缩放维度为 $[\mu] = 4K$。当 $K > 1/2$ 时，该算符是不相关的（irrelevant）；而在临界点 $\mu=0$ 且 $K < 1/2$ 时，它是唯一的对称性允许的相关（relevant）微扰。

1.2.3 动态声子的重整化效应

将动态声子场积掉（integrate out）之后，低能有效作用量中生成了一个重整化的 Luttinger 参数 $K'(\omega_0)$。通过重整化群计算得到：

这表明：有限频率的动态声子总是会减小有效 Luttinger 参数（$K' < K$）。随着声子频率 $\omega_0$ 的降低，$\gamma$ 迅速减小，导致 $K'$ 不断被压低。

当 $K'$ 降至 $1/8$ 以下时，更高阶的对称性允许项 $-\rho \cos(8\phi)$ 成为相关算符（$[\rho] = 16K'$，在 $K' < 1/8$ 时 $[\rho] < 2$）。此时，系统必须用 双正弦-Gordon（double sine-Gordon, DSG）作用量 来描述：

这一项正是破坏连续临界性并将其转化为强一阶相变的罪魁祸首。

1.3 技术难点与应对策略

1. 玻色自由度截断与计算复杂性： 声子具有无限维的局部希尔伯特空间。为了在张量网络中处理声子，必须对声子占用数进行合理截断。本工作在局部空间中引入截断 $N_{\text{max}} = 6$。通过系统性检验不同截断下基态物理的偏差，确认在 $\omega_0 \ge 0.4$ 时其系统相对偏差低于 $10^{-3}$，验证了截断的精确性。
2. 相变点附近的极强迟滞效应： 在一阶相变和临界终点附近，常规的自洽 DMRG 极易陷入亚稳态（Metastable state）。为此，本工作设计了双向扫频协议（Sweep protocol），分别从 Néel 相和 VBS 相出发，并在模拟中引入衰减密度混合器（Density mixer），其幅值从 $10^{-5}$ 衰减至 $10^{-8}$，以确保生成足够的纠缠并避免系统过早陷入亚稳态，从而精确提取两个相分支的交叉点。
3. 动态格林函数的实时演化： 利用时变变分原理（TDVP），在包含 $L=256$ 个声子/自旋格点并外接无限边界环境的杂化大系统上进行长时演化（最大时间 $t_m \approx 80/J_1$）。之后通过线性预测（Linear prediction）外推至原时间的两倍，并应用高斯窗口消除吉布斯振荡（Gibbs ringing），高精度还原出声子谱函数的动量-频率细节。

2. 关键 Benchmark 体系与物理性能数据

2.1 零声子耦合极限（无晶格畸变 benchmark）

对于纯自旋链（$\Delta=4$），在调控参数 $J_2/J_1$ 变化时：

• 连续量子临界点确定为：$J_{2,c}/J_1 = 0.3923$。
• 在临界点处，Néel 序与 VBS 序的残余值随着 MPS 截断键维度 $\chi$ 的增加，呈现出普适的标度消失规律：$O_{J_{2c}} \sim \xi(\chi)^{-\beta/\nu_\mu} \sim \xi(\chi)^{-K}$（见图2(a)插图），其中 Luttinger 参数精确提取为 $K = 0.372$。这完美证实了临界点的二阶连续特征。
• 序参数关联函数在长程表现出代数衰减：$\langle M^z(0) M^z(r) \rangle \sim r^{-2K}$，其幂律指数与 $K = 0.372$ 完美吻合（见图2(b)）。

2.2 静态声子不稳定性极限（Adiabatic limit）

在声子频率 $\omega_0 \to 0$ 且晶格刚度 $\kappa = 3.5$ 固定的情况下：

• 系统的能量面 $E_0(\delta)$（作为晶格交替畸变幅值 $\delta$ 的函数）在 $J_2/J_1$ 较小时最小点位于 $\delta=0$（对称 Néel 相）。
• 随着 $J_2/J_1$ 的增加，能量在有限二聚化位移 $\delta_{\text{opt}} \neq 0$ 处发展出新的全局极小值，导致最优畸变 $\delta_{\text{opt}}$ 发生突变（见图3(c)），同时伴随着 Néel 序和 VBS 序的断崖式跳跃。
• 关联长度 $\xi$ 在相变点处不再发散，而是保持在有限值，直接证实了静态声子场必然将一维 DQCP 转化为强一阶相变（自旋-Peierls 不稳定性）。

2.3 动态声子体系（核心数据）

对于 $\Delta_1=\Delta_2=4$ 且耦合强度 $g=0.25J_1$ 的动态体系：

2.3.1 弱耦合连续相变区（$\omega_0 = 0.6J_1 > \omega_c$）

• 相变点左移至：$J_{2,c} \approx 0.2573J_1$，这是由于声子涨落倾向于稳定 VBS 基态。
• 序参数连续消失（见图4(a)），残余标度分析给出临界点 Luttinger 参数为 $K' \approx 0.221$ 到 $0.237$。这显著小于无声子时的 $K=0.372$，定量证实了动态声子对 Luttinger 参数的压低效应。
• 纠缠熵 $S$ 与关联长度 $\xi$ 表现出共形场论（CFT）行为：$S \sim \frac{c}{6} \log \xi$，拟合得到中心电荷 $c \approx 0.98$ 到 $1.12$（普适值 $c=1$），对应于无质量自由玻色子。
• 临界指数确定为：$\beta = 0.189$，$\nu = 0.78$。这也与标度关系 $K' = \beta/\nu \approx 0.242$ 一致。

2.3.2 强一阶相变区（$\omega_0 = 0.4J_1 < \omega_c$）

• 序参数在相变点发生高达 $0.4$ 数量级的非连续性跳跃（见图5(a)）。
• 关联长度 $\xi$ 迅速饱和至有限值，能量一阶导数 $\partial E/ \partial J_2$ 出现不连续台阶（见图5(c)）。
• VBS 序的 Binder 累积量 $U_B$ 标度分析（见图5(d)）：随着系统尺寸 $L$（$L = 64, 128, 256, 512$）的增大，在相变点处出现急剧发散的负尖峰。这是由于相相干（Phase coexistence）导致的双峰概率分布，是一阶相变极为显著和无可争议的数值判据。

2.3.3 临界终点（$\omega_0 = \omega_c \approx 0.45J_1$）

• 在此临界终点上，重整化 Luttinger 参数刚好触及临界值 $K' = 1/8 = 0.125$（见图1(c)）。

• 有限尺寸标度分析（Finite-size scaling） 拟合（公式 (16)）：

  $$\hat{O} \chi^{\beta\kappa/\nu} = F(\delta J_2 \chi^{\kappa/\nu})$$

  在磁化强度和 VBS 序参数下展示出近乎完美的标度崩塌（Collapse，见图6(a,b)），提取到的临界指数为：

  $$\beta_M \approx 0.08, \quad \beta_\Psi \approx 0.087, \quad \nu \approx 0.635$$

  这与经典的 二维四态 Potts（4-state Potts）共性类 的解析值（$\beta = 1/12 \approx 0.083$，$\nu = 2/3 \approx 0.667$）在误差范围内高度一致，成功将此一维自旋-声子系统终点划归为此高对称性普适类。

2.3.4 声子谱函数响应数据（Kohn 异常）

• 有序相内（Néel 或 VBS 相）： 声子谱表现出能隙，且在布里渊区边缘 $q=\pi$ 处由于自旋交替起伏发生显著变软（Softening），呈现出 Kohn 异常特征。

• 连续临界点（DQCP）： 声子模式在 $q=\pi$ 处完全无能隙，在低频展现出连续体（Continuum）特征（见图8中间子图）。

• 谱函数幂律标度： 在临界点处，声子谱函数表现出非平凡的幂律衰减（见图9(a)）：

  $$A(q=\pi, \omega) \sim \omega^{-2K'}$$

  对于 $\omega_0 = 0.6J_1$，其谱权重在低频极限下的扫频行为与此解析幂律高度重合，证实了声子与临界二聚化起伏的强杂化。

3. 代码实现细节与复现指南

3.1 软件栈与核心软件包

数值计算全部基于 Python 开源张量网络库 TeNPy (Tensor Network Python)（由 Johannes Hauschild 和 Frank Pollmann 维护）。该库提供了高效的基于矩阵乘积态（MPS）和矩阵乘积算符（MPO）的数值算法，并原生支持任意一维局部对称性的守恒（本工作严格守恒总 $S^z$ 算符）。

• 开源代码与原始数据仓库： 本工作的复现数据和仿真脚本均已在 Zenodo 开源发布，读者可以通过以下永久链接直接下载： https://doi.org/10.5281/zenodo.20037182

3.2 自旋-声子晶格的 MPS 表示

由于包含自旋和声子两种不同的格点，哈密顿量构建中，晶格采用 自旋与玻色子交替格点链 表示：

• 奇数格点 $2j$ 放置 Spin-1/2 局部空间，维数为 2。
• 偶数格点 $2j+1$ 放置声子局部空间，采用被截断的 Boson 空间，维数为 $d_{\text{boson}} = N_{\text{max}} + 1 = 7$。对应的局部湮灭算符 $a$ 表达为 $7 \times 7$ 的实数矩阵。

3.3 iDMRG 扫频复现步骤

要重构本工作中发现的一阶/二阶转变和确定临界点：

1. 初始化态： 为追踪迟滞区，分别用 tenpy.networks.mps.MPS.from_product_state 建立两种初始态：
   • Néel 反铁磁态：['up', '0', 'down', '0'] * (L/2)
   • VBS 二聚二聚态：['up', '0', 'up', '0'] 的自旋二聚对（在局部键上投影成 Singlet）。
2. 密度混合参数调优（Density Mixer）：

   mixer_params = {
       'amplitude': 1e-5,
       'decay_rate': 0.5,
       'decay_steps': 20,
       'disable_after': 30
   }
   

3. 扫频与交叉点提取： 对于给定的声子频率 $\omega_0$，以步长 $\Delta (J_2/J_1) = 0.001$ 逐步改变哈密顿量。通过计算当前步收敛态在下一步哈密顿量下的期望能（作为初始态），从而描绘出两条亚稳能流曲线 $E_{\text{AFM}}(J_2)$ 与 $E_{\text{VBS}}(J_2)$。相变点 $J_{2,c}$ 由这两条能流曲线的交叉点决定：$\Delta E = E_{\text{VBS}} - E_{\text{AFM}} = 0$。

3.4 时变 TDVP 算子演化

1. 格林函数测量： 演化在基态 $\lvert \psi_0 \rangle$ 上施加局部算符 $A = (a + a^\dagger)$。
2. 二格点 TDVP (2-site TDVP)： 使用 tenpy.algorithms.tdvp.TwoSiteTDVPEngine 执行实时演化，步长设为 $dt = 0.2/J_1$。允许最大辅助键维 $\chi$ 增长至 $600$。演化持续进行，直到最大截断时间 $t_{\text{max}} = 80/J_1$。
3. 线性预测外推： 对得到的离散关联函数 $G(r, t)$，提取时间截取区数据并使用 Python 中的自回归线性外推算法，将其外推一倍至 $t = 160/J_1$，而后利用 Gaussian 窗函数进行离散傅里叶变换，求得高精度动量谱 $A(q, \omega)$。

4. 关键引用文献与局限性评论

4.1 关键引用文献

1. Senthil et al., Science 303, 1490 (2004) [文献1]： 首次提出二维解禁闭量子临界点（DQCP）理论体系，奠定了超越兰道对称破缺范式的物理基础。
2. Roberts et al., PRB 99, 165143 (2019) [文献21]： 提出了易轴各向异性自旋链中的一维 DQCP 模型，并论证了其在无声子极限下的连续相变特征，本工作的自旋哈密顿量设计直接继承于此。
3. Hofmeier, Willsher, Seifert, Knolle, PRB 110, 085121 (2024) [文献32]： 基于重整化群场论分析，预测了晶格动态畸变对一维 DQCP 的压低和潜在的一阶化机制，为本工作提供了直接的场论出发点。
4. Giamarchi, Quantum physics in one dimension (2003) [文献22]： 一维玻色化与 Luttinger 液体理论的权威教科书，本工作附录 A 中自旋链向 sine-Gordon 场论映射的全部技术细节均基于此书。
5. Haegeman et al., PRL 107, 070601 (2011) [文献62]： 时变变分原理（TDVP）在矩阵乘积态中的核心算法实现，本工作计算声子谱函数的关键数值技术。

4.2 本工作局限性评论

虽然该项工作在理论和数值上极为详实，但从真实的量子化学和实验物理角度看，仍存在以下重要局限性：

1. Einstein 声子模型的局限： 本工作将晶格振动简化为互不耦合的 Einstein 声子（即各格点声子频率独立且为常数 $\omega_0$），在低频近似下这相当于一条平带。而在实际的三维晶体中，晶格振动具有高度色散的声学支（Acoustic branch）与光学支（Optical branch）。声学声子在 $q \to 0$ 时能量趋于零，会与自旋体系发生更长程、更复杂的非局域杂化。如何将真实的声子色散关系纳入张量网络模拟，仍是一个巨大的计算瓶颈。
2. 玻色截断 $N_{\text{max}}=6$ 的潜在失效区： 虽然作者通过能量偏差论证了在 $\omega_0 \ge 0.4$ 时截断足够精准，但当声子频率进一步降至极低频区域（例如更接近静态极限的 $\omega_0 \sim 0.1$）或者在极强的自旋-声子耦合 $g$ 下，声子实空间占用数将显著增大，高能态的玻色涨落无法忽略。此时强制截断为 6 会导致物理量的非物理饱和，从而扭曲一阶相变的真实演化边界规律。
3. 维度的鸿沟（1D to 2D）： 本文针对的是一维自旋链，虽然这允许高精度的张量网络模拟，但一维体系中是不存在真正的长程 Néel 序的（仅在各向异性 $\Delta > 1$ 的有限能隙相中存在有效序）。而在真实的强关联电子材料（如 Shastry-Sutherland 晶格材料 $\text{SrCu}_2(\text{BO}_3)_2$）中，DQCP 发生于二维体系。将本工作中的动态声子积掉技术推广至二维，由于纠缠面积律的限制，目前的二维张量网络（PEPS）或大尺寸 DMRG 还很难实现包含动态自旋-声子耦合的精确实时演化，这限制了该结论在二维材料中的直接定量应用。
4. 场论分析中双正弦-Gordon 势预因子 $\rho$ 的推导缺乏微观自洽性： 工作在第 III B 节中虽然指出了相变一阶化依赖于 $-\rho \cos(8\phi)$ 项中 $\rho > 0$ 的势垒性质，但这一结论是通过数值结果“倒推”得出的。如何通过全非线性玻色化技术，从自旋-声子耦合哈密顿量的微观参数（如 $J_1, J_2, g, \omega_0$）自洽、解析地导出来 $\rho$ 的正负号和精确大小，依然是一个悬而未决的理论难题。

5. 补充解析：双正弦-Gordon 势能面的微观拓扑结构与实验探测图景

5.1 双正弦-Gordon（DSG）势能极小值微观物理图景

为了深化对该物理图像的理解，我们有必要详细剖析双正弦-Gordon 势能函数（等式 (17)）的拓扑极小值结构如何决定不同的基态物理：

利用简单的三角恒等式，我们将极小值解（等式 (18)）在相空间中展开，并绘制其物理对应：

• 情况 1：$\mu > 0$ 且 $\rho > -\mu/4$ 极小值位于 $\phi_0 = 0, \pi/2$，对应于 $\cos(2\phi)$ 取极大值。这正是 Néel 反铁磁序的特征。此时自旋取向交替排列。
• 情况 2：$\mu < 0$ 且 $\rho > -\mu/4$ 极小值转移至 $\phi_0 = \pi/4, 3\pi/4$，对应于 $\sin(2\phi)$ 取极大值。这对应于二聚化的 VBS 序，空间平移对称性破缺。
• 情况 3：$\mu = 0$ 且 $\rho > 0$（一阶临界点） 此时势能极小值同时位于 $\phi_0 = 0, \pi/4, \pi/2, 3\pi/4$，意味着 Néel 序与 VBS 序的四个简并态完全共存（见图7中间的四个势阱图）。然而，由于 $\rho > 0$，在相邻势阱之间（例如 $\phi = \pi/8$ 处）存在一个势垒。正是这个经典的势垒阻止了简并基态之间的量子隧穿和连续过度，直接迫使相变变为一阶突变。
• 情况 4：$\rho < 0$（经典的 Ashkin-Teller 行为） 在经典的 Ashkin-Teller 模型中，由于其对偶性，预因子 $\rho$ 总是负的。在这种情况下，势垒不复存在，原有的单一临界点分裂为两个相继的 Ising 相变点，其中极小值在相变点之间可以连续滑移（即中间存在共存中介相）。

本工作最重要的贡献之一就是在数值上无可辩驳地证实了自旋-动态声子相互作用产生的是 $\rho > 0$ 的非 Ashkin-Teller 类型有效场，这在凝聚态理论界具有极高的系统重要性。

5.2 实验探测：如何在实验室中寻找 1D DQCP 的声子印记？

对于实验凝聚态物理和量子化学工作者来说，该项研究指明了利用现有的谱学技术探测解禁闭量子临界性的重要方向：

5.1.1 实验候选材料建议

实验研究人员可重点关注以下具有准一维特征且具有强自旋-声子耦合作用的低维反铁磁各向异性材料：

1. $\text{BaCo}_2\text{V}_2\text{O}_8$ 及 $\text{SrCo}_2\text{V}_2\text{O}_8$： 这是目前国际上公认研究一维 XXZ 链及 Ising 型临界性的标准候选材料。其强易轴各向异性 $\Delta \approx 2$ 且自旋-晶格耦合极其显著，可通过施加外加磁场调控 $J_2$。在这类材料中，通过**非弹性中子散射（INS）**测量声子能量在倒空间边缘的变软，可直接映射出本文所述的 Kohn 异常。
2. $\text{CuGeO}_3$（经典的自旋-Peierls 材料）： 虽然该材料在极低温下会发生二聚化自旋-Peierls 相变，但通过高压或者杂质掺杂可以压制二聚化温度。在压制临界点附近，声子谱在低能限是否呈现非相干连续体以及谱指数是否偏离常数，是检验本文 DQCP-声子杂化理论（$A \sim \omega^{-2K'}$）的极佳实验切入点。