来源论文: https://arxiv.org/abs/2606.01254v1 生成时间: Jun 02, 2026 17:28

非厄米自旋弛豫下 GZ 量子多体伤疤态的动力学稳定化：从理论、iTEBD 模拟到哈伯德模型映射的深度解析

0. 执行摘要

自热力学与统计力学建立以来，孤立量子多体系统的热力学化过程一直是凝聚态物理和量子信息科学的核心议题。本征态热力学化假说（Eigenstate Thermalization Hypothesis, ETH）为这一过程奠定了现代量子力学基础。然而，近年来物理学家们发现，在某些强非积分、混沌的多体系统内部，存在极少数不符合 ETH 预测的特殊高能激发本征态，它们被称为量子多体伤疤态（Quantum Many-Body Scars, QMBS）。这些状态在单一时间演化下表现出弱遍历性破缺，能长期保持局域量子信息，因而具有重大的量子计算与量子精密测量应用价值。

在自旋为 $S \ge 1$ 的各向异性 XYZ 量子自旋链中，存在一类精确的自旋螺旋形积态本征态——Granovskii-Zhedanov (GZ) 伤疤态。尽管 GZ 态在理想的幺正动力学下能完美保持局域自旋相干性，但在实际物理平台（如冷原子光晶格、超导量子比特或固态自旋链）中，环境噪声、晶格无序及不可避免的哈密顿量微扰会迅速破坏其脆弱的相干性，导致其无可挽回地热力学化。为了克服这一瓶颈，由 Dhiman Bhowmick 撰写的最新研究论文《Stabilization of Granovskii–Zhedanov Scars of the XYZ Quantum Spin Chain via Non-Hermitian Spin Relaxation》提出了一种创新的物理方案：通过引入非厄米自旋弛豫过程配合外加螺旋磁场，可以极大地增强 GZ 伤疤态的生存寿命，并将其稳定在具有高保真度的非平衡稳态（Non-Equilibrium Steady State, NESS）中。

本文将面向系统研究量子多体物理、量子化学模拟及开放量子系统的科研人员，对该工作进行全方位、深层次的技术解构。我们将系统推导 GZ 态的代数构造，详细分析非厄米稳定化算符的微观物理起源，阐明基于张量网络（MPS/iTEBD）的非厄米演化模拟技术，展示哈伯德模型向非厄米 XYZ 自旋模型的 Schrieffer-Wolff 微扰映射，并对比 Lindblad 耗散主方程与非厄米哈密顿量描述的等效性。最后，我们将剖析该理论方案在实验实现上的潜在局限，并指明未来的研究方向。

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

1.1 核心科学问题：扰动下量子多体伤疤态的退相干与稳定化

对于一个混沌且满足 ETH 的孤立多体系统，几乎所有高能本征态的局域物理量观测值都与微正则系综下的热力学平衡态一致。量子多体伤疤（QMBS）作为一种弱遍历性破缺现象，其本征能量在光谱中并非孤立，而是隐匿在密集的“热本征态海洋”之中，但其纠缠熵显著低于周围的热态，呈现对数或常数标度特征。这种特殊的代数结构使它们在特定初始状态下能展现出长期的量子相干振荡。

然而，伤疤态的稳定性极易受到系统微扰的破坏。在现实实验中，理想的 $XYZ$ 相互作用往往会叠加多体跳跃微扰（如费米子流动）、局域磁场无序或退相干通道。这些微扰会在伤疤子空间与热力学本征态子空间之间建立起非零的基元跃迁矩阵元，从而将局域信息快速扩散到整个多体希尔伯特空间中。如何建立一种实验可行的、主动的动力学稳定化协议（Active Stabilization Protocol），在混沌多体背景下无限期或显著延长伤疤态的相干寿命，是该领域亟待解决的核心科学问题。

1.2 理论基础：GZ 伤疤态的严格代数构造

在 $S \ge 1$ 的各向异性一维 $XYZ$ 自旋链中，系统哈密顿量为：

其中，$\hat{S}_n^{\alpha}$ 为位置 $n$ 处的自旋-$S$ 算符。利用雅可比椭圆函数（Jacobi Elliptic Functions）$cn(u, \kappa)$、$sn(u, \kappa)$ 和 $dn(u, \kappa)$，我们可以精确地构造出一类螺旋自旋积态（Product States），即 Granovskii-Zhedanov (GZ) 态。首先，参数 $q$ 和椭圆模数 $\kappa$ 由哈密顿量的各向异性耦合系数决定：

GZ 态的积态形式定义为：

其中 $|\uparrow_n\rangle$ 是算符 $\hat{S}_n^z$ 的最大特征值本征态（即 $|S, S\rangle$ 状态）。局域旋转角 $\theta_n$ 和 $\phi_n$ 的引入，使得每一个格点处的自旋取向都沿着由经典三维自旋向量 $\langle \hat{\mathbf{S}}_n \rangle$ 决定的方向：

而自旋期望值的空间分布则由如下的雅可比椭圆轨线刻画：

这里 $\alpha = \sqrt{1 - \gamma^2}$，$\beta = \sqrt{1 - \gamma^2 + \kappa^2 \gamma^2}$，$\gamma \in [-1, 1]$ 是调控三维螺旋结构取向的连续变参。通过调节自由参数 $\gamma$ 和 $\phi$，我们实际上可以定义一个极其庞大的 GZ 伤疤子空间，其维度通常具有多项式标度。对于有限周期边界条件的自旋链，波矢 $q$ 被限制为离散值 $q = 4pK(\kappa)/N$（其中 $p \in \mathbb{Z}$，$K(\kappa)$ 为第一类完全椭圆积分）。

在理想情况下，由于 GZ 态中的每一个局域自旋都被转动到了由哈密顿量决定的局部各向异性有效场方向上，多体算符作用在该积态上仅产生局域守恒流（详见附录 A 的严格证明），使其成为 $\hat{H}_0$ 的严格本征态：

1.3 核心技术方案：非厄米弛豫与螺旋磁场的协同控制

为了对抗有害微扰 $\hat{H}_P$，作者提出了两个层面的控制手段：

1. 幺正控制（螺旋磁场 $\hat{H}_h$）： 施加一个完全贴合 GZ 态局部自旋纹理方向的局域螺旋磁场：

   $$\hat{H}_h = -h \sum_{n=1}^{N} \mathbf{B}_n \cdot \hat{\mathbf{S}}_n$$

   其中，$\mathbf{B}_n = \left( \alpha cn(nq, \kappa), \beta sn(nq, \kappa), \gamma dn(nq, \kappa) \right)$。显然，由于 GZ 态在格点 $n$ 处的局部自旋方向与 $\mathbf{B}_n$ 完全对齐，该螺旋磁场作用于 GZ 态仅产生恒定的能量标移，即 $|\Psi_{\text{GZ}}\rangle$ 仍是 $\hat{H}_0 + \hat{H}_h$ 的本征态。在物理上，该磁场极大地加宽了伤疤态与相邻热态之间的有效能隙，抑制了局域多体微扰引发的共振跃迁。然而，数值模拟表明，纯幺正磁场仅仅是延缓了热力学化，并不能在 $t \to \infty$ 极限下阻止信息的衰减。

2. 非厄米控制（耗散弛豫 $\hat{H}_\lambda$）： 引入非厄米项：

   $$\hat{H}_\lambda = i\lambda \sum_{n=1}^{N} \hat{S}_n^{z'}$$

   其中，$\hat{S}_n^{z'}$ 是格点 $n$ 处的局域转动基底自旋算符，其 $z'$ 轴指向局域自旋期望值 $\langle \hat{\mathbf{S}}_n \rangle$ 的方向（参见 Fig. 1(a)）。若采用 Schwinger 玻色子（或自旋子 Spinon）表符化定义，我们可以将该局域自旋算符写为：

   $$\hat{S}_n^{z'} = \frac{1}{2} (\hat{c}_{n\uparrow}^{\dagger} \hat{c}_{n\uparrow} - \hat{c}_{n\downarrow}^{\dagger} \hat{c}_{n\downarrow})$$

   这里，$\hat{c}_{n\uparrow}^{\dagger}$ 和 $\hat{c}_{n\downarrow}^{\dagger}$ 分别创造上自旋子和下自旋子（分别对应 GZ 局域对齐基态和激发态）。带入后，非厄米项在自旋子表象下表示为：

   $$\hat{H}_\lambda = \frac{i\lambda}{2} \sum_{n=1}^{N} (\hat{c}_{n\uparrow}^{\dagger} \hat{c}_{n\uparrow} - \hat{c}_{n\downarrow}^{\dagger} \hat{c}_{n\downarrow})$$

   物理图景（Fig. 1(b)）：当 $\lambda > 0$ 时，虚数能量项对下自旋子（激发态）引入了负虚部（即衰减率），而对上自旋子（基态）引入了正虚部（即放大率）。这意味着系统的耗散动力学在不断地将偏离 GZ 轨线的热激发粒子“泵浦”回局部自旋对齐的基态（即 GZ 态的格点局域真空）。因此，非厄米项扮演了“自旋冷凝器”的角色，不断提取系统由于微扰产生的多体熵，从而达到长时间稳定伤疤态的目的。

将该局部转动算符逆变换回实验室参考系，非厄米稳定项可以用常规自旋算符的内积简洁表示：

1.4 技术难点与方法细节

1. 非厄米时间演化的数值挑战： 在厄米哈密顿量的时间演化中，态振幅的模长始终守恒。但在包含非厄米项 $\hat{H}_\lambda$ 之后，演化算符 $\hat{U}(t) = \exp(-i \hat{H} t)$ 不再是幺正的。在张量网络方法（如 iTEBD）中，非幺正算符的作用会导致矩阵乘积态（MPS）的奇异值谱发生剧烈畸变，若直接应用标准幺正 TEBD 代码，数值误差将随时间呈指数放大。为此，必须在每一次 Trotter 演化步后进行显式的波函数归一化，并动态调整奇异值截断阈值以维持 MPS 的物理代表性。论文中采用的有效保真度定义为转移矩阵的最大特征值，这巧妙地规避了热力学极限下传统波函数重叠积恒等于零的问题。

2. 物理机制的可实现性证明： 纯数学上的非厄米算符（例如 $i\lambda \hat{S}_n^z$）很容易在理论上写出，但其实验可行性一向饱受质疑。如何从一个包含实物粒子（如费米子、声子）的厄米系统，严格推导出该非厄米项，是凝聚态理论上的一个难点。作者通过构建双轨道费米子哈伯德模型，借助 Schrieffer-Wolff 变换，成功完成了这一严格的理论映射，证明了非厄米稳定化算符可以自然地源于声子诱导的自旋弛豫与 Purcell 效应。

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

为了系统检验非厄米自旋弛豫对 GZ 伤疤态的稳定效果，作者设计了针对无限大自旋链（Thermodynamic Limit）与有限尺寸自旋链（Finite Size）的 Benchmark 模拟。以下是详细的数值实验结果与性能分析：

2.1 Benchmark 系统参数配置

• 基底各向异性 XYZ 自旋链：设定 $\kappa = 0.5$, $\gamma = 0.8$, 椭圆波矢 $q = 4K(\kappa)/5$, 自旋量子数 $S=1$。哈密顿量系数均以 $J_y = 1$ 为基准单位进行标度。
• 初始状态：一般初始化为严格的 GZ 积态 $|\Psi_{\text{GZ}}\rangle$。
• 微扰类型：
  • 局域单体横场微扰 (Onsite Perturbation)：$\hat{H}_P^{(1)} = \mu \sum_n \hat{S}_n^x$（设定标称值 $\mu = 0.3$）。
  • 自旋跳跃微扰 (Hopping Perturbation)：$\hat{H}_P^{(2)} = t \sum_n (\hat{S}_n^+ \hat{S}_{n+1}^- + \text{H.c.})$（设定标称值 $t = 0.3$）。

2.2 核心动力学数据与性能表现 (Fig. 2 深度解构)

2.2.1 局域微扰下的动力学稳定表现

在存在局域微扰 $\mu = 0.3$ 且无任何稳定机制时（$\lambda = 0, h = 0$，见 Fig. 2(a) 中的黑色实线），GZ 态保真度 $\mathcal{F}_{GZ}$ 迅速发生指数衰减，在极短时间内 ($t < 5$) 跌落至 0.25 以下，表现出混沌多体系统典型的热力学化行为。而在引入非厄米自旋弛豫项（$\lambda > 0$）后，系统的动力学特征发生了根本性的转变：

• $\lambda = 0.2$（蓝色虚线）：保真度在经历短暂小幅下降后，于 $t \approx 10$ 时迅速拉平并稳定在 $\mathcal{F}_{GZ} \approx 0.65$ 的稳态平台。
• $\lambda = 0.4$（蓝色实线）：稳定化性能进一步提升，稳态保真度在全时段内被牢牢锁定在 $\mathcal{F}_{GZ} \ge 0.90$。这无可辩驳地证明了非厄米自旋弛豫机制在混沌多体背景下对局域量子相干性的高保真度保护能力。
• $\lambda < 0$ 的非物理对立情况：当引入负的弛豫率（如 $\lambda = -0.2, -0.4$，见 Fig. 2(a) 的红色曲线）时，热力学化过程不仅没有被抑制，反而被剧烈加速，保真度在 $t < 3$ 内便直奔零点，这从反面印证了耗散的方向性（即指向局域自旋对齐方向的泵浦）对于伤疤态稳定性的决定性作用。

2.2.2 动能跳跃微扰下的稳定表现 (Fig. 2(b))

相较于局域微扰，非局域的自旋跳跃微扰（$t = 0.3$）通常会引入多体流动，其解相干效应更为剧烈。在无稳定化机制时，保真度瞬时归零。然而，得益于局域非厄米算符在能谱上对多体激发态的强烈压制：

• 当设定 $\lambda = 0.4$ 时，即使在强跳跃微扰下，稳态保真度仍能维持在 $\mathcal{F}_{GZ} \approx 0.85$ 左右（Fig. 2(b) 蓝色实线）。
• 同样，反向流 $\lambda = -0.4$（红色虚线）会导致态保真度瞬间坍塌。

2.2.3 幺正螺旋磁场的局限性比较 (Fig. 2(c))

为了凸显非厄米控制方案的优越性，作者单独模拟了仅外加幺正局域螺旋磁场 $h$（无非厄米耗散，即 $\lambda = 0$）的情况：

• 当不加磁场 $h = 0$ 时，系统热力学化（黑色曲线）。
• 当施加正磁场 $h = 0.2, 0.4$ 时，保真度的初阶段指数衰减被抑制，且保真度曲线上移（蓝色曲线）。然而，在较长时间极限下（$t \approx 15$），保真度仍呈现出持续单调下滑的趋势。这表明幺正能隙控制只能延缓多体纠缠向环境的扩散，而无法建立动力学自适应机制来消除熵的持续产生。只有非厄米耗散项的存在，才能将系统引向一个确定的、有限保真度的非平衡稳态。

2.2.4 GZ 态的动力学泵浦产生 (Fig. 2(d) & (e))

除了维持既有的伤疤态，非厄米自旋弛豫机制还具备动力学态制备（Dynamical State Preparation）的能力。作者模拟了从完全平凡的磁性积态出发的演化过程：

• 铁磁全对齐态 $|1111\dots\rangle$（因为其与 GZ 伤疤子空间存在较大的初始重叠 $\mathcal{F}_{GZ}(0) = 0.5$）：在正向自旋弛豫 $\lambda = 0.4$ 的驱动下，即使在强微扰下，态保真度也能在 $t \approx 5$ 时自适应地攀升并稳定在 1.0 附近（实线蓝/红线）。这代表着一个完美的自动泵浦制备过程。
• 反铁磁对齐态 $|1\bar{1}1\bar{1}\dots\rangle$ 与结构复杂的积态 $|1\bar{1}\bar{1}1\bar{1}\bar{1}\dots\rangle$：由于它们本身处于高能热本征态空间，与 GZ 伤疤子空间的重叠度微乎其微（$\mathcal{F} \approx 0$），即便施加 $\lambda = 0.4$ 稳定项，也无法实现向 GZ 态的动态汇聚（保真度始终保持在 0 附近）。这深刻揭示了该泵浦机制的选择性只对“具有一定相干前驱体”的局域状态生效。

2.3 稳定性相图分析 (Fig. 5 深度解构)

为了给出完整的参数生存空间，作者在附录中利用大量 iTEBD 数据绘制了稳态保真度 $\mathcal{F}_{GZ}(t=15)$ 在 $(\mu, \lambda)$ 及 $(t, \lambda)$ 参数空间中的相图（Fig. 5(c)-(f)）。他们将 $\mathcal{F}_{GZ} \ge 0.9$ 定义为“完美稳定区”，并拟合出了分界线（相边界）：

我们对关键性能参数进行定量提取，制成下表：

这些 Benchmark 数据为物理学家在实验上评估微扰容忍度提供了极其精确的参考。特别是，大自旋系统虽然对抗局域退相干的能力更强，但对空间自旋流动的微扰则显著更加敏感，这一物理机制在设计光晶格自旋链实验时具有极强的指导价值。

3. 代码实现细节、复现指南与开源工具链

为便于同行复现论文中的关键动力学结果，本节提供一份完整的技术路线图，详细说明如何使用主流张量网络软件包 TenPy 和量子动力学工具箱 QuTiP 来编写数值模拟代码。

3.1 理论基础：非厄米 iTEBD 算法细节

要在无限大矩阵乘积态（MPS）上应用 iTEBD 模拟非厄米演化，必须对标准的二阶 Suzuki-Trotter 分解进行修改。由于 $\hat{H}$ 包含非厄米项，时间步长 $\delta \tau$ 对应的演化算符：

是不保模长的。复现该算法的核心步骤如下：

1. 非等效相互作用键算符构造：将哈密顿量分解为奇数键 $A$ 和偶数键 $B$ 的局域两格点算符 $\hat{h}_{i, i+1}$。对于非厄米部分，需将单体项（如 $i\lambda \mathbf{B}_n \cdot \hat{\mathbf{S}}_n$）平均分配到相邻的键算符中。
2. Trotter 算符指数化：计算非幺正矩阵 $U_{bond} = \exp(-i \hat{h}_{bond} \delta\tau)$，这会产生奇异值不为 1 的局域演化算符。
3. SVD 压缩与规范化（Crucial Step）：作用 $U_{bond}$ 于 MPS 后，通过奇异值分解（SVD）更新兰姆达矩阵（Singular Values $\Lambda$）。在非厄米演化中，$\Lambda$ 的平方和不再等于 1。必须手动计算 $\mathcal{N} = \sqrt{\sum_k \lambda_k^2}$ 并执行 $\Lambda \leftarrow \Lambda / \mathcal{N}$ 进行显式物理归一化。
4. 有效保真度计算：对于无限大系统，无法直接使用 $\langle \Psi(0) | \Psi(t) \rangle$，而必须通过计算左、右固定 MPS 之间的转移矩阵（Transfer Matrix）的最大本征值来获取单位长度的有效保真度 $\mathcal{F}_{GZ}$。

3.2 TenPy 代码复现框架：构建非厄米自旋模型

以下是基于 Python 开源张量网络库 TenPy (Tensor Network Python) 复现该论文核心动力学演化（Fig. 2(a)）的实验级代码框架：

import numpy as np
from tenpy.models.model import CouplingModel, NearestNeighborModel
from tenpy.networks.site import SpinSite
from tenpy.models.lattice import Chain
from tenpy.algorithms import tebd
from tenpy.networks.mps import MPS

class NonHermitianXYZChain(CouplingModel, NearestNeighborModel):
    def __init__(self, model_params):
        # 1. 提取物理参数
        L = model_params.get('L', np.inf) # 设定格点数，inf 代表热力学极限
        S = model_params.get('S', 1.0)
        Jx = model_params.get('Jx', 0.86)   # 各向异性相互作用强度
        Jy = model_params.get('Jy', 1.0)
        Jz = model_params.get('Jz', 0.5)
        h = model_params.get('h', 0.0)      # 幺正磁场强度
        lmbda = model_params.get('lmbda', 0.4) # 非厄米弛豫强度
        mu = model_params.get('mu', 0.3)    # 局域单体微扰强度
        
        # 构造自旋基底
        site = SpinSite(S=S, conserve=None) # 非厄米系统通常不守恒自旋对称性
        lattice = Chain(L, site, bc='infinite' if np.isinf(L) else 'periodic')
        CouplingModel.__init__(self, lattice)
        
        # 2. 构造局域相互作用 (H0)
        for u in range(lattice.N_sites):
            self.add_coupling(Jx, u, 'Sx', u, 'Sx', dx=1)
            self.add_coupling(Jy, u, 'Sy', u, 'Sy', dx=1)
            self.add_coupling(Jz, u, 'Sz', u, 'Sz', dx=1)
            
        # 3. 添加微扰项 (HP)
        if mu != 0:
            self.add_onsite(mu, 0, 'Sx')
            
        # 4. 重点：添加非厄米项 H_lambda 与 幺正螺旋磁场 H_h
        # 为简化，在实验室框架下直接加入复杂耦合常数的局域磁场项
        # 以 n=0 格点为例（实验中需根据雅可比椭圆公式计算每个格点的 Bn 方向）
        # 对应的总场算符为 (-h + i * lmbda) * (Bx * Sx + By * Sy + Bz * Sz)
        # 注意：此处必须使用 complex 复数耦合参数以实现非厄米性
        for u in range(lattice.N_sites):
            # 此处数值仅作为椭圆局部磁场的示意，具体椭圆坐标根据公式 (3) 生成
            Bx, By, Bz = 0.57, 0.0, 0.82 
            heff = -h + 1j * lmbda
            
            self.add_onsite(heff * Bx, u, 'Sx')
            self.add_onsite(heff * By, u, 'Sy')
            self.add_onsite(heff * Bz, u, 'Sz')
            
        NearestNeighborModel.__init__(self, lattice)

# 5. 时间演化执行与保真度监控
def run_simulation():
    model_params = {
        'L': np.inf, 'S': 1.0,
        'Jx': 0.86, 'Jy': 1.0, 'Jz': 0.5,
        'h': 0.0, 'lmbda': 0.4, 'mu': 0.3
    }
    model = NonHermitianXYZChain(model_params)
    
    # 初始化 MPS 态为非厄米的 GZ 积态
    # 通过旋转算符将 |1111...⟩ 变换到本地螺旋轴向
    # 此处省略具体局域酉变换旋转生成，直接使用模型基态或定义积态
    psi = MPS.from_product_state(model.lat.mps_sites(), ["up"] * model.lat.N_sites, bc=model.lat.bc)
    
    # 设定演化参数
    tebd_params = {
        'order': 2, 
        'dt': 0.001,      # Trotter 步长，必须足够小
        'N_steps': 10,     # 每次输出的步数
        'trunc_params': {'chi_max': 120, 'svd_min': 1.e-10}, # MPS 截断维度
        'verbose': 1
    }
    
    eng = tebd.Engine(psi, model, tebd_params)
    
    times = []
    fidelities = []
    for t in range(150): # 演化到 t = 15
        eng.run() # 内部自动处理非厄米演化的波函数重归一化
        current_t = eng.evolved_time
        
        # 计算与初始 GZ 态的有效重叠（保真度）
        # 对于无限大自旋链，通过转移矩阵求解最大特征值
        fid = eng.psi.overlap(psi, charge_sector=None) # 此处在有限截断下做近似评估
        
        times.append(current_t)
        fidelities.append(np.abs(fid))
        print(f"Time: {current_t:.3f}, Fidelity: {np.abs(fid):.6f}")

if __name__ == "__main__":
    run_simulation()

3.3 开源工具链及 Repo 推荐

1. TenPy 官方仓库：
   • 链接：https://github.com/tenpy/tenpy
   • 说明：基于矩阵乘积态（MPS）算法最强大的 Python 开源张量网络库。其内部的 tebd.Engine 在处理非幺正实时间演化（Non-Unitary Time Evolution）时能自动重归一化奇异值，是实现本研究非厄米 iTEBD 算法的最佳选择。
2. QuTiP 官方仓库：
   • 链接：https://github.com/qutip/qutip
   • 说明：用于模拟开放量子系统动力学的开源 Python 库，包含强大的 ODE 求解器。可用于极小尺寸（如 $N=9$）的精准对角化（Exact Diagonalization）演化以及 Lindblad 主方程（Eq. 27）的数值积分。

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

4.1 关键引用文献及其在科学树中的位置

1. [36, 37] Y. I. Granovskii and A. S. Zhedanov (1985, 1986):
   • 位置：理论基石。首次在各向异性一维 $XYZ$ 量子自旋链中发现并解析构造出了特殊的螺旋积态。本工作正是建立在这一精确代数解的基础之上。
2. [12, 13] C. J. Turner, A. A. Michailidis, D. A. Abanin, M. Serbyn, and Z. Papić (2018):
   • 位置：QMBS 现代复兴的起点。在里德堡原子非正规链（PXP 模型）中首次通过实验和数值证实了量子多体伤疤的存在，点燃了多体遍历性破缺的研究热潮。
3. [53-56] V. Popkov, C. Presilla, J. Schmidt et al. (2016, 2017):
   • 位置：方法学先驱。率先探索了利用非厄米边界耗散在 $XXZ$ 自旋链中稳定自旋螺旋态（Spin-Helix States）的机制。本工作将这一思想成功地推广到了更为普适且混沌的各向异性 $XYZ$ 自旋-$S$（$S\ge 1$）体系及体积耗散通道中。
4. [58] J. Hauschild and F. Pollmann (Scipost Phys. 2018):
   • 位置：计算工具。发布了极其完备的开源张量网络库 TenPy。本研究所使用的 iTEBD 核心算法和热力学极限保真度提取技术，均深度依赖该开源库。

4.2 本工作局限性之批判性评论

尽管本论文在物理学理论上极具美感，且数值模拟展示了惊人的稳定化性能，但站在客观、冷峻的凝聚态物理实验与工程实现视角，该方案仍存在以下不可忽视的物理局限性：

1. 对“后选择”（Post-selection）或非物理监测的高度依赖

非厄米哈密顿量 $\hat{H}_\lambda$ 的本质是对开放量子系统在无量子跃迁（No-Jump）子空间中动力学的有效近似描述。在标准的开放量子系统（如 Lindblad 描述）中，除了非厄米有效演化外，还必然伴随着随机发生的“量子跃迁”（Quantum Jump）过程（即 $L_j \rho L_j^{\dagger}$ 项，它会将部分概率流无规地散射到其他多体态中）。 要实现严格的非厄米哈密顿量演化，通常需要对系统实施强连续弱测量并进行“后选择”，即丢弃所有发生过量子跃迁的实验轨迹。随着系统格点数 $N \to \infty$ 以及时间 $t \to \infty$，完全不发生跃迁的轨迹概率呈指数衰减（$P_{\text{no-jump}} \sim e^{-N \lambda t}$）。这意味着在大系统和长演化尺度下，直接后选择的效率将降为零，极大限制了该协议在多比特信息保持中的直接工业应用。

2. 空间不均匀螺旋磁场的实验设计难度

稳定化方案要求施加一个空间连续旋转的螺旋外磁场（$\hat{H}_h$ 和 $\hat{H}_\lambda$ 均正比于局域向量 $\mathbf{B}_n$），其周期性由雅可比椭圆波矢 $q$ 极其精确地定义。在固体自旋链器件或超导比特阵列中，要在纳米级尺度下任意调控、实现空间上每个格点磁场方向及强度的个性化定制，工艺难度呈几何级数增长。 虽然作者在论文中提出了利用“环状光晶格+倾斜均匀场”（Fig. 3(c)）来在局部参考系中模拟螺旋场的简化方案，但该方案仅适用于非常特定且单一的自旋波矢 $q = 2\pi/N$ 的特殊情形，无法灵活稳定更具一般代数性质的 GZ 伤疤子空间。

3. 声子诱导弛豫的不可控环境背噪

在利用磁致伸缩效应（Section IV.A）推导自旋弛豫微观机制时，作者假设晶格产生的声子能迅速且无阻碍地逸散到外部无限大衬底介质中，不发生二次重吸收激发。但在实际的三维多晶或低温超导衬底中，声子局域反射与过热声子瓶颈（Phonon Bottleneck）广泛存在。这些无法排除的声子反馈会产生等效的热噪声，相当于在原本定向的非厄米泵浦上叠加了无序的“热泵浦”，进而导致 NESS 稳态保真度发生显著退化。

5. 关键补充：哈伯德微观映射与 Lindblad 主方程等同性证明

为了展示该工作的完整物理图像，本节深度拆解两个关键理论支撑点：哈伯德模型在强关联极限下的微扰映射，以及非厄米哈密顿量与 Lindblad 主方程的等价性区间。

5.1 深度拆解：双轨道费米子哈伯德模型向非厄米 XYZ 自旋链的映射

在实验实现上，单纯构造一个相互作用各项异性且自带非厄米局域弛豫的自旋-1 哈密顿量是极其困难的。作者指出，这一系统可以通过一个更底层的、完全符合实物粒子特征的双轨道费米子哈伯德模型（Two-orbital Fermionic Hubbard Model）在半满及强关联极限下，通过第二阶微扰理论（Schrieffer-Wolff 变换）严格推导出来。

微观费米子哈密顿量 $\hat{H}$ 如等式（20）所示，包含两套轨道 $m=1,2$。在相互作用极强（$U, U' \gg \tau$）且满足半满（每格点填充两个电子）时，极强的排斥力将电荷自由度完全冻结，只留下自旋和轨道自由度。更精妙的是，格点内的超交换相互作用（Hund’s coupling $J_H$）会将低能能谱强行分裂成三态（对称的自旋转角自旋-1 三重态，Eq. (21)）和高能单重态。三重态的形式为：

这三个状态构成了一个完美的等效自旋-1 自由度！当我们在这套系统上施加一个微观的虚数自旋项（源于电极或光学微腔耦合导致的非等性损耗）：

通过二阶 Schrieffer-Wolff 变换投影到自旋-1 三重态低能流形（Low-Energy Manifold）中，其等效的自旋交换相互作用耦合常数（Eq. (26)）变为：

其中，超交换耦合系数直接与微观费米子跳跃项相关：

这完美地证明了：微观跳跃振幅的自旋各向异性（$\tau_{\uparrow} \neq \tau_{\downarrow}$）直接导致了等效低能自旋链中各向异性 XYZ 相互作用的产生，而微观上对费米子施加的不等号虚数损耗，在投影后精确地转变为了局域指向 GZ 本地自旋对齐方向的非厄米自旋弛豫算符 $\hat{H}_\lambda$！ 这一令人叹为观止的映射，为在冷原子光晶格（如双轨道碱土金属原子 $^{173}\text{Yb}$ 或 $^{87}\text{Sr}$）中人工合成非厄米量子伤疤态扫清了理论障碍。

5.2 Lindblad 主方程与非厄米演化的等同性与区别

为了消除学术界对非厄米有效哈密顿量“非幺正”物理基础的疑虑，作者在论文第 V 节将数值结果与完整的 Lindblad 耗散主方程描述进行了详细比对。系统的标准开放主方程具有如下形式：

通过对其进行代数变换，主方程可以等价改写为：

其中，有效哈密顿量为 $\hat{H}_{\text{eff}} = \hat{H}_0 - \frac{i}{2} \sum_j \Gamma_j L_j^{\dagger} L_j$。这展示了两种框架的内在关联：

• 等同性：若系统只进行短时间演化，或者对测量设备进行严苛的后选择（监测到系统未发生量子跃迁，即把 $\sum_j \Gamma_j L_j \rho L_j^{\dagger}$ 这一“跳跃项”完全滤除），系统的密度矩阵演化便完全由非厄米有效哈密顿量 $\hat{H}_{\text{eff}}$ 决定。

• 区别与等效稳定机制：当无法进行后选择，必须考虑完整物理系统的长期 NESS 稳态动力学时，跳跃项无法被忽略。在 GZ 态的特殊结构下，由于 GZ 本身就是局域上自旋子的真空态，局域自旋阶梯算符满足极高阶消灭性质。作者将非厄米自旋弛豫项进行特制谱分解，构造出如下形式的 Lindblad 跳跃算符（Eq. (30)）：

  当 $\lambda > 0$ 时，选择 $L_n = \hat{S}_n^{+\prime}$ 作为主方程的耗散通道。

  在这一耗散作用下，即使包含跳跃项（其扮演的是将系统从其他无关激发态退激发回到 GZ 本质子空间的角色），长期演化的 NESS 仍能表现出与非厄米演化定性高度契合的稳态保真度保护能力（Fig. 4 数值证实）。这也强有力地向实验界宣告：即便后选择不能完美执行，直接依靠量子通道耗散（如 Purcell Spontaneous Emission）也能够无损地实现对 GZ 伤疤态强有力的物理稳定！

6. 结语

Dhiman Bhowmick 的这项工作不仅具有极高的学术价值，更为量子多体物理和量子信息控制提供了一条崭新的思路：“化敌为友，借力打力”——我们不应一味追求将多体系统与有害环境绝对隔离（这在热力学上是不可能的），相反，应当主动、精确地设计、引入环境耗散，在非厄米和开放系统框架下，利用自组织动力学将混沌系统中的高相干状态（如多体伤疤态）锚定为系统长期的动力学不动点。 本研究所揭示的通过微观哈伯德模型合成非厄米自旋自适应弛豫的方案，以及对 GZ 态保真度的量化评估，必将极大地推动多体相干调控及新一代量子拓扑材料的研究进程。