一维量子基态中 U(1) 对称性能否自发破缺？

一维基态为何被认为不会破缺连续对称性

Hohenberg-Mermin-Wagner 定理排除了在一维和二维、相互作用足够短程的经典平衡系统中，非零温度下连续对称性的自发破缺。通过量子-经典对应，即把 d 维量子系统映射为 d+1 维经典系统，这一结果让人预期连续对称性在一维量子基态中同样不会破缺。Tomonaga-Luttinger 液体是标准例证：对数增长的相位涨落取代真正的长程有序，只留下代数衰减的准长程有序。

论文强调，这一对应并不自动继承经典定理所依赖的前提。在一维量子晶格系统中，对连续对称性破缺的严格限制只在附加假设下才能导出。各向同性铁磁 Heisenberg 链是熟悉的例外：即使在一维，完全极化的基态也破缺 SU(2) 自旋旋转对称性；但在那里，序参量与哈密顿量对易，基态流形包含极化乘积态。较新的无阻挫构造可以在序参量不守恒的情况下破缺连续对称性，因此一般机制此前仍属未决。

经典理论缺失的 Berry 相位要素

论文的场论出发点是自旋相干态路径积分。相邻虚时上相干态的重叠会产生一个 Berry 相位项，它由自旋取向扫过的立体角几何地决定。这一项在欧几里得作用量中纯为虚数，在经典统计力学的 Boltzmann 权重中没有直接对应物。

在铁磁体中，Berry 相位使方位角相位与纵向磁化强度成为正则共轭量，从而改变长距离涨落，使其不同于相对论性理论。基于 Nahum 的提议，作者考虑把 SU(2) 自旋旋转对称性约化为 U(1)，同时保留磁化反转对称性和空间对称性。这一约化允许纵向磁化强度出现局域四次相互作用，该相互作用在一维空间中是相关的，会破坏高斯不动点的稳定性。在固定正四次相互作用下调节二次密度项，便可接近一个相互作用不动点，在那里相位涨落及其共轭磁化强度即使在长距离上仍保持非平凡耦合。

子系统 Binder 标度定位相变点

自旋-1 模型守恒总纵向磁化强度，因此本征态可取确定的 M_z，非简并有限链基态必然位于 M_z = 0 扇区。在该扇区中，磁化强度密度的二阶矩和四阶矩都为零，常规 Binder 比因此无定义。作者改用半链的磁化强度，它在总磁化强度固定时仍可涨落，并把 R_half 定义为其四阶矩与二阶矩平方之比。

在 M_z = 0 扇区用 DMRG 得到若干系统尺寸下 R_half 随各向异性参数 g_1 的变化。曲线在 g_1 约 -0.826 附近的狭窄区域相交，尺寸 L 与 L+2 之间的交点 g_x(L, L+2) 随 L 增大而收敛。用 R_half = F_p(x) + L^{-omega}(b_0 + b_1 x)、x = (g_1 - g_1c) L^{1/nu}、多项式次数 p = 3、4、5，拟合窗口 |g_1 + 0.8258| <= 0.01，得到 g_1c = -0.8258 ± 0.0002 和 nu = 0.60 ± 0.02。所给不确定度是反映拟合选择敏感性的保守系统范围，而非单次拟合的统计误差。

nu = 0.60 指向相互作用而非高斯临界性

估计值 nu = 0.60 ± 0.02 不同于高斯值 nu_G = 1/2。等价地说，相关微扰的 RG 本征值 1/nu 不同于其高斯值 2，作者将其作为相互作用临界行为的数值证据。该值也落在微扰 RG 结果 nu 约 0.5833（epsilon 阶）与 nu 约 0.6258（epsilon 平方阶）之间。

在自旋-1/2 链中，调节 g_1 会同时改变相干态拉格朗日量中二次项和梯度项的系数，但其行为在定性上相同。其半链 Binder 比在狭窄耦合区间相交，用同一标度形式、L_max = 68，得到 g_1c = 0.3826 ± 0.0003 和 nu = 0.61 ± 0.02。该 nu 与自旋-1 结果一致，同样大于高斯值，说明相互作用临界性并不局限于单一自旋长度。

热力学极限下横向磁化强度为正

由于基态保持 U(1) 对称性，M_x 和 M_y 的期望值各自为零。作者因此追踪横向磁化强度平方密度 m_perp^2(L) = <M_x^2 + M_y^2>_0 / L^2，在 M_z = 0 扇区的最低能量态中计算。正的极限 m_perp^2 = lim_{L to infinity} m_perp^2(L) > 0 足以判定 U(1) 对称性自发破缺，因为在取完热力学极限后再撤去无穷小横向场时，非零横向磁化强度仍然存在。

在临界耦合 g_1 = -0.8258 处，对周期自旋-1 链做 DMRG，并用二次拟合 m_perp^2(L) = m_perp^2 + a_1/L + a_2/L^2、20 <= L <= 72，得到 m_perp^2 约 0.9473。固定 L_max = 72、把 L_min 从 20 变到 40，估计值保持在 0.9469 到 0.9473 之间；在 L_min = 20 下比较线性、二次和三次拟合，得到 0.94687 到 0.94889 之间的值。采用的估计是 m_perp^2 = 0.947 ± 0.003，明显为正。作为对照，铁磁 Heisenberg 链在 S_tot = L、M_z = 0 的基态中有精确有限尺寸结果 m_perp^2(L) = 1 + 1/L。

有限纠缠标度在无限长链上确认有序

为直接研究无限长链，作者使用带 U(1) 对称矩阵乘积态拟设的 VUMPS。这里的截断不是系统尺寸，而是键维 chi，它限制 MPS 能表示的纠缠，因此可观测量必须通过有限纠缠标度外推。由于 S_i^+ 携带 U(1) 荷 q = 1，横向关联探测的是 q = 1 转移矩阵扇区，该扇区中最长的关联长度 xi_1(chi) 便作为标度变量。

作者先在固定块长度下消除有限纠缠效应，再把块长度外推到无穷，即在探测长程有序之前先取 chi 趋于无穷的极限。所得估计为 m_perp^2 约 0.949，与有限尺寸外推一致。对自旋-1/2 链，在 g_1 = 0.3826 处对 m_perp^2(L) 做有限尺寸标度，得到热力学极限估计 0.212，保守区间 0.209 <= m_perp^2 <= 0.215，对照的 Heisenberg 精确结果为 m_perp^2(L) = 1/4 + 1/(2L)。

z 小于 2 排除无阻挫实现

作者分析最低中性激发能隙的有限尺寸标度，得到动力学临界指数 z = 1.986 ± 0.005，它刻画低能色散 omega 正比于 |k|^z。这相对高斯值 z_G = 2 是微小但显著的偏离。论文还用 VUMPS 检验纵向自旋关联：比较不同试探标度维数时，由 z 推出的 Delta_psi = 0.507 使重标度关联曲线相对其均值的相对均方根散布比高斯值 Delta_psi = 0.5 减小约 29.5%。作者提醒，这种塌缩改善并不构成对 Delta_psi 的精确独立测定，因为最小化散布得到的有效维数会随 chi_min 从 512 增至 1536（固定 chi_max = 3072）而从约 0.520 漂移到 0.531。

作者得出的关键推论是：z < 2 与无阻挫模型对该临界理论的微观实现不相容，因为无阻挫哈密顿量要求 z >= 2。因此一维连续对称性破缺并不止于无阻挫构造，暗示存在一种不同于无阻挫性的机制。论文还指出，同一经典场论在非平衡 Ising 晶格气体中有微观实现，其 Monte Carlo 模拟给出 nu = 0.62 ± 0.03 和 nu = 0.61 ± 0.03，与 DMRG 值一致；而对应的有效动力学指数 1.94 ± 0.02 和 1.965 ± 0.005 与本文报告的 z 有偏离，作者将此差异标记为值得进一步研究。

判断破缺结论前应比较什么

本文最具决策价值的比较，不是某条一维链是否破缺连续对称性，而是哪个诊断量确立了破缺、哪个指数区分了机制。破缺结论依赖外推的横向磁化强度平方密度：自旋-1 链为 m_perp^2 = 0.947 ± 0.003，自旋-1/2 链为 0.212，区间 0.209 到 0.215；这些都在保持对称性的态中得到，而非通过加入对称破缺场。相互作用性质依赖 nu = 0.60 ± 0.02 和 nu = 0.61 ± 0.02，二者都高于高斯值 1/2，并被微扰 RG 值夹在中间。机制论断依赖 z = 1.986 ± 0.005 低于 2，即无阻挫哈密顿量所要求的阈值。

读者在权衡这一结果时应留意三条边界。这些指数是保守的系统范围，与拟合窗口和多项式次数选择有关，而非统计置信区间。自旋-1 链的临界耦合 g_1c = -0.8258 ± 0.0002 与早先的 g_c = -0.826 ± 0.005 一致，但报告的 z 与早先的 z = 1.50 ± 0.04 差别很大；把早先的 beta = 0.30 ± 0.01 与其 z 通过标度关系结合，得到 nu = 0.40 ± 0.02，与本文的值不相容。Delta_psi 塌缩的改善只是提示性的，并非精确测定。因此论文自身的重点落在 z < 2 这一阈值上：它是把该临界理论与无阻挫实现区分开来的唯一数字，也支撑了一维连续对称性破缺存在不止一条路径的论断。

本文引用的文献