来源论文: https://arxiv.org/abs/2606.25411v1 生成时间: Jun 25, 2026 16:24

动态晶格上的自旋不平衡费米子：基于 DMRG 的一维关联磁性与双模式自旋密度波（SDW）深度解析

0. 执行摘要

在一维强关联电子系统中，电荷与自旋自由度的相互耦合和竞争是产生新奇量子相（如非平庸磁有序、超导配对以及拓扑相）的核心机制。传统的固体物理通常将晶格视为静态背景，或者将其振动简化为声子。然而，近年来“动态晶格”（Dynamical Lattice）的概念在冷原子物理和量子模拟中脱颖而出：晶格本身的自由度（如化学键上的局域自旋或规范场）与流动费米子的运动发生直接的动力学反馈，构建出一种高度非线性的、各向异性的动态耦合体系。这种耦合最典型的代表之一便是 $\mathbb{Z}_2$ 规范场耦合的费米子系统及量子链接模型（Quantum Link Models）。

本文针对山西大学量子光学与光电器件国家重点实验室（State Key Laboratory of Quantum Optics Technologies and Devices, Shanxi University）的刘杰、周晓凡和贾锁堂等人的最新研究成果进行深度学术解析。该研究关注一个一维自旋不平衡（Spin-imbalanced）的自旋-1/2费米子动态晶格模型。在这个模型中，流动的费米子通过一种类似于 Ising 耦合的自旋依赖性跃迁（Spin-dependent hopping），直接与键中心（Bond-centered）的局域自旋发生耦合。流动费米子的跃迁振幅受到其穿过键上局域自旋构型的直接调制，反过来，费米子的自旋密度分布也作为有效磁场，反作用于局域自旋。

利用先进的**密度矩阵重正化群（DMRG, Density Matrix Renormalization Group）数值模拟方法，研究团队系统地绘制了不同参数空间下的基态相图。最为引人注目的发现是，在部分自旋极化费米子相（Partially spin-polarized fermion phase）**中，局域自旋的结构因子 $S^z(k)$ 展现出了两种性质截然不同的自旋密度波（SDW）有序波矢。它们分别完美地对应于自旋分解费米面（Spin-resolved Fermi surfaces）的嵌套波矢（Nesting vectors） $2k_{F\uparrow}$ 和 $2k_{F\downarrow}$。这两种自旋密度波相（命名为 SDW-I 和 SDW-II）在引入排斥性的 Hubbard-$U$ 相互作用后表现出极强的鲁棒性，甚至 Hubbard-$U$ 还通过抑制双占位有效增强了费米子的自旋极化倾向。这一研究揭示了一种在低维关联体系中人工调控磁性调制的新途径，深化了学术界对动态晶格磁性及一维皮尔斯（Peierls）类似转变的微观理解。

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

1.1 核心科学问题

在一维多体系统中，量子涨落极为剧烈，导致了传统的费米液体理论失效。取而代之的是 Tomonaga-Luttinger 液体（TLL） 理论，其特征是自旋-电荷分离（Spin-charge separation）和各种具有幂律关联的电荷、自旋密度波波动。对于二组分自旋平衡（即自旋向上和向下的粒子数相等， $N_\uparrow = N_\downarrow$）的一维费米子系统，其自旋密度波（SDW）的主导有序波矢通常遵循经典的普适关系：

$$k_{max} = 2k_F$$

其中，当总填充密度为 $\rho = N/L$ 时，费米波矢定义为 $k_F = \pi \rho / 2$。因此，主导波矢与总密度具有简单的线性关系 $k_{max} = \pi \rho$。

然而，一旦引入自旋不平衡（Spin imbalance, $N_\uparrow \neq N_\downarrow$），系统的物理图像将发生深刻改变。这不仅为体系注入了额外的新奇自由度（如可以稳定类似 FFLO 的非零动量配对超导态或部分极化的金属态），更对自旋空间调制提出了挑战：

在引入动态晶格耦合（Bond-centered 局域自旋调制跃迁）后，费米子自旋不平衡如何反过来调控局域自旋的 SDW 有序波矢？
这种由动态晶格介导的相互作用是否会分裂出不同的 SDW 模式？
系统的临界行为和全局相图在参数耦合下如何演化？

该研究试图解决的核心科学问题正是：在自旋不平衡流动费米子与化学键上局域自旋相互耦合的动态晶格模型中，其基态自旋调制波矢的微观机制是什么？它如何与自旋分解的费米面发生定量关联？

1.2 理论基础与模型构建

研究团队考虑的一维自旋-1/2 费米子动态晶格模型由以下哈密顿量定义（其中设置普朗克常数 $\hbar = 1$）：

$$\hat{H} = -t \sum_{i,\sigma} (\hat{c}_{i,\sigma}^\dagger \hat{c}_{i+1,\sigma} + H.c.) - \alpha \sum_{i,\sigma} (\hat{c}_{i,\sigma}^\dagger \hat{S}_i^z \hat{c}_{i+1,\sigma} + H.c.) + h \sum_i (\hat{n}_{i,\uparrow} - \hat{n}_{i,\downarrow}) + \frac{\Delta}{2} \sum_i \hat{S}_i^z + \beta \sum_i \hat{S}_i^x + U \sum_i \hat{n}_{i,\uparrow} \hat{n}_{i,\downarrow}$$

其中各物理项与算符含义如下：

流动费米子项： $\hat{c}_{i,\sigma}^\dagger$ ( $\hat{c}_{i,\sigma}$ ) 表示在晶格格点 $i$ 上创建（湮灭）一个自旋为 $\sigma \in \{\uparrow, \downarrow\}$ 的费米子。 $\hat{n}_{i,\sigma} = \hat{c}_{i,\sigma}^\dagger \hat{c}_{i,\sigma}$ 是相应的粒子数算符。 $t$ 代表常规的最近邻跃迁速率。
局域自旋度规： $\hat{S}_i^z$ 和 $\hat{S}_i^x$ 是居住在格点 $i$ 与 $i+1$ 之间**化学键中心（Bond-center）**的自旋-1/2 算符。它们代表动态晶格的局部量子化自由度。
自旋依赖的动态跃迁耦合（关键项）：第二项是耦合常数为 $\alpha$ 的项，它将流动费米子的跃迁直接与键中心局域自旋 $\hat{S}_i^z$ 联系起来。这种耦合导致沿该键的有效跃迁振幅变为：
$$t_{eff} = -(t \pm \alpha)$$
具体取决于局域自旋是处于向上（$+1/2$）还是向下（$-1/2$）的状态。因此，费米子的空间运动受到了局域自旋构型的动力学调制，带有典型的动态晶格特征。
流动费米子塞曼场：第三项 $h$ 为施加在流动费米子上的外部塞曼场，用于调节费米子自旋向上和向下组分的布居数不平衡（即改变极化度 $P = (N_\uparrow - N_\downarrow)/(N_\uparrow + N_\downarrow)$）。
局域自旋磁场：第四项 $\Delta$ 和第五项 $\beta$ 分别为施加在键中心局域自旋上的纵向磁场（偏向自旋 $z$ 方向排列）和横向量子涨落场（诱导自旋在 $x-z$ 面发生翻转，驱使相变发生）。
Hubbard 项：最后一项 $U$ 是原位排斥相互作用，控制两组分费米子之间的电荷关联。

这一模型在数学上与 $\mathbb{Z}_2$ 规范场耦合费米子系统及量子链接模型高度关联，其中局域自旋扮演了介导费米子间非线性相互作用的动态规范场角色。

1.3 技术难点

处理该模型面临三大显著的技术与计算难点：

非局部量子纠缠与局域希尔伯特空间膨胀：由于晶格点和键中心自旋交替排列（物理构型为 ... - Fermion_i - Spin_i - Fermion_{i+1} - Spin_{i+1} - ...），在计算中，必须将每个格点和邻近的键联合定义为一个复合的局域希尔伯特空间。对于原位有自旋-1/2 费米子（4维空间：空位、单占向上、单占向下、双占）和自旋-1/2 局域自旋（2维空间：向上、向下），两者的联合张量积空间维数为 $4 \times 2 = 8$。如果引入 Hubbard 相互作用并保留所有电荷守恒，在构建矩阵乘积态（MPS）时需要极为精确地处理各种对称性分区。
多体项 MPO 复杂性：哈密顿量中的自旋依赖跃迁项 $\hat{c}_{i,\sigma}^\dagger \hat{S}_i^z \hat{c}_{i+1,\sigma}$ 涉及三个格点算符（格点 $i$ 的费米子算符、键 $i$ 的自旋算符、格点 $i+1$ 的费米子算符）。在标准的 DMRG 矩阵乘积算符（MPO）表示中，这类三体交互作用项会显著增加 MPO 的键维数（Bond dimension），对扫描收敛性提出了更高的要求。
大参数空间的精细相变判定：系统包含 $t, \alpha, h, \Delta, \beta, U, \rho$ 共 7 个自变量，相互之间存在极强的竞争。要在海量的参数组合中精确定位从顺磁（PM）到自旋密度波（SDW）的二阶量子相变点，并排除有限尺寸效应，计算开销极其庞大。

1.4 方法细节：极高精度 DMRG 模拟

为了攻克上述难点，作者采用了基于矩阵乘积态（MPS）架构的最先进的 密度矩阵重正化群（DMRG） 方法。模拟技术细节如下：

边界条件：采用开边界条件（OBC, Open Boundary Conditions），这虽然会在边界处引入局域反射，但在 1D 强关联计算中，开边界条件可以使 MPO 的形式更加紧凑，同时避免了周期性边界条件（PBC）下键维数呈指数增长的劣势。
尺度选取：选择了一系列典型的晶格尺寸： $L = 24, 36, 48, 60$。系统总格点数在 DMRG 内部由于自旋与费米子双链结构被转化为交替排列，因此实际优化的 MPS 长度最长可达 $2L-1 = 119$ 个位置。
精度控制与截断：
- 保持每个 DMRG 子块（Block）的保留态数（Kept States / Bond Dimension）高达 $M = 600$。
- 在每次扫描中进行了至少 60 次完整的扫频（Sweeps），并辅以细微的噪声扰动（Noise term）以避免 DMRG 优化落入局域极小值。
- 截断误差（Truncation Error）上限严格控制在极其精确的 $10^{-9}$ 以下。
物理可观测量的提取：
- 局域自旋结构因子 $S^z(k)$：
  $$S^z(k) = \frac{1}{L} \sum_{i,j} e^{ik(i-j)} \langle (\hat{S}_i^z - \bar{S}^z)(\hat{S}_j^z - \bar{S}^z) \rangle$$
  其中 $\bar{S}^z = \sum_i \langle \hat{S}_i^z \rangle / L$ 是局域自旋的平均期望值。结构因子最大峰值的位置 $k^{max}$ 和高度 $S^{max}$ 是判定 SDW 有序的关键依据。
- 实空间自旋-自旋关联函数 $S^z(r)$：
  $$S^z(r) = \frac{1}{L} \sum_l \langle (\hat{S}_l^z - \bar{S}^z)(\hat{S}_{l+r}^z - \bar{S}^z) \rangle$$
  用于区分长程有序的幂律振荡（SDW 相）与无序的指数衰减（PM 相）。
- 费米子动量分布函数 $n_\sigma(k)$：
  $$n_\sigma(k) = \frac{1}{L} \sum_{l,j} \langle \hat{c}_{l,\sigma}^\dagger \hat{c}_{j,\sigma} \rangle e^{i(l-j)k}$$
  用于直接探测费米面（Fermi Surface）位置及其对应的嵌套波矢。

2. 关键 Benchmark 体系、计算数据与物理相图深度剖析

2.1 长程有序判定与有限尺寸外推 (Benchmark of Long-Range Order)

研究团队首先分析了无 Hubbard 相互作用（ $U=0$）情况下系统基态的磁性特征。选取总填充 $\rho = 1$ （半填充），费米子塞曼场 $h = 1.0$，自旋依赖跃迁强度 $\alpha = 0.5$：

强纵向场极限（例如 $\Delta = 1.0$）：
- 如图 2(a1) 所示，随着系统尺寸 $L$ 从 24 增大到 60，缩标后的局域自旋结构因子 $S^z(k)/L$ 在有限波矢 $k \approx 0.67\pi$ 处形成一个极其陡峭、尖锐的 Bragg 峰。而在 $k=0$ 和 $k=\pi$ 处值极小。
- 图 2(a2) 展示了在 $1/L \to 0$ 极限下的有限尺寸外推（Finite-size scaling）。明显可见，在外推至热力学极限（ $L \to \infty$）后，主导峰值强度 $S^z(k^{max})/L$ 收敛到一个确定的大于零的有限值（ $\approx 0.1$），表明系统存在非零的自旋长程自组织调制。
- 与此对应，图 2(a3) 描绘了实空间关联函数 $S^z(r)$。结果呈现强烈的空间余弦型持续振荡，没有发生任何空间衰减。这无懈可击地证明了系统处于一个长期稳定的空间调制自旋密度波（SDW）相。
弱纵向场极限（例如 $\Delta = 0.1$）：
- 当 $\alpha > \Delta$ 导致自旋跃迁机制完全主导局域自旋分布时，物理图像突变。如图 2(b1) 和 2(b2) 所示，在 $k=0.67\pi$ 处的峰值随 $1/L$ 迅速萎缩，在热力学极限下完全消失为零。
- 图 2(b3) 显示其实空间关联函数 $S^z(r)$ 表现出近乎完美的指数衰减，几乎在距离 $r > 3$ 后便衰竭至零。这表明此时系统的长程磁有序被完全破坏，系统进入了磁无序的顺磁（PM）相。

2.2 磁相变与塞曼场调控 (Phase Transitions & Zeeman Field Effect)

为了更定量地描绘这一磁性量子相变，研究人员固定参数 $\rho = 1/2$（ 1/4 填充）， $\alpha = 0.5$, $\Delta = 1.0$，以流动费米子塞曼场 $h$ 作为控制变量进行了精细扫频：

主导有序波矢 $k^{max}$ 的连续漂移：如图 3(a) 所示，在零温无塞曼场（ $h=0$）时， $k^{max} = \pi / 2$。随着 $h$ 从 0 开始连续增加，流动费米子开始发生自旋极化（自旋向上和向下的密度不再对称），主导波矢 $k^{max}$ 从 $\pi/2$ 开始单调且平滑地减小。这直接说明塞曼场引入的自旋极化连续拓宽了自旋调制的空间周期（即 $2\pi / k^{max}$ 变大）。
临界点 $h_c$ 的确定：图 3(b) 绘制了外推的热力学有序强度 $S^{max} = \lim_{L \to \infty} S^z(k^{max})/L$。在小 $h$ 区域内， $S^{max}$ 保持在 $\approx 0.08$ 的稳健平台。随着 $h$ 接近临界场 $h_c \approx 0.7$ ， $S^{max}$ 呈幂律形式连续下降并最终归零。在 $h > h_c$ 后， $S^{max}$ 严格为零，标志着系统从共存的、空间调制的 SDW 态，经历了一个连续的量子相变，过渡到了无序的 PM 金属态。这一跃迁可以被精确归类为二阶量子相变。

2.3 双模式自旋密度波的发现与嵌套机制（SDW-I & SDW-II）

本研究最核心的物理发现展示在图 4 和图 5 中。通过大范围改变费米子的局域自旋密度 $\rho_\uparrow$ 和 $\rho_\downarrow$，研究人员在部分极化区域内观测到了两种迥异的 SDW 响应规律：

SDW-I 模式（自旋向上密度 $\rho_\uparrow$ 线性主导）：
- 当自旋依赖跃迁耦合较弱（如 $\alpha = 0.4$, 图 4(a1, a2)）时，SDW 相仅出现在中间极化区。随着 $\rho_\uparrow$ 的增加，主导波矢 $k^{max}$ 表现出与 $\rho_\uparrow$ 强烈的、具有正斜率的线性依赖关系。这一依赖关系精确地遵循以下定量公式：
  $$\frac{k^{max}}{\pi} \approx 2\rho_\uparrow = 2\rho - 2\rho_\downarrow$$
- 微观起源分析（图 5(a)）：对该区域费米子动量分布 $n_\sigma(k)$ 的计算表明，由于自旋向上分量（多数派）占主导，多数派费米子在 $k_{F\uparrow} = \pi \rho_\uparrow$ 处形成一个清晰的准费米面。通过键中心自旋介导的倒格矢过程，局域自旋感受到了多数派费米子费米面的嵌套不稳定性（Fermi Surface Nesting），从而激发了波矢为 $2k_{F\uparrow}$ 的自旋密度波。这正是 SDW-I 模式。
SDW-II 模式（自旋向下密度 $\rho_\downarrow$ 线性主导）：
- 然而，当增大自旋依赖耦合（如 $\alpha = 0.5$, 图 4(b1, b2)）时，令人惊异的现象发生了：SDW 相被顺磁相从中段截断，并在高极化（即 $\rho_\uparrow$ 趋向 0.5 且 $\rho_\downarrow$ 极小）和低极化两个两极区域稳定下来。在这两个极性区中， $k^{max}$ 展现出了与 $\rho_\uparrow$ 呈负斜率的线性关系，等价于它与少数派自旋密度 $\rho_\downarrow$ 呈正斜率关系：
  $$\frac{k^{max}}{\pi} \approx 2\rho_\downarrow = 2\rho - 2\rho_\uparrow$$
- 微观起源分析（图 5(b)）：在此强耦合区域，尽管少数派自旋向下费米子密度极低，但其动量分布 $n_\downarrow(k)$ 在更窄的边界 $k_{F\downarrow} = \pi \rho_\downarrow$ 处形成高度局域化的费米突变。这种少数派费米子的嵌套机制超越了多数派，成功接管了局域自旋调制，产生具有 $2k_{F\downarrow}$ 空间周期的磁有序。此即为 SDW-II 模式。

这种在一个简洁哈密顿量中由耦合常数 $\alpha$ 调控的两种嵌套波矢之间的磁性竞争（从多数派嵌套向少数派嵌套的转换），是一维磁性理论上的一个重大突破。这可以被视为在动态局域自旋扇区中发生的、由费米面嵌套不稳定性驱动的皮尔斯不稳定性（Peierls-like instability）的磁学对应物。

2.4 全局相图及其参数依赖性 (Global Phase Diagrams)

研究人员在图 6 和图 7 中通过大尺度 DMRG 模拟完整绘制了该动态晶格系统的二维全局量子相图：

$\alpha-\Delta$ 平面相图（图 6(a, b)）：
- 在部分极化（PP, $h=0.5$）和完全极化（FP, $h=2.0$）下，相图均由边界清晰的顺磁相（PP-PM, FP-PM）和自旋密度波相（PP-SDW, FP-SDW）组成。这揭示了纵向磁场 $\Delta$ 与量子耦合 $\alpha$ 的直接对抗： $\Delta$ 强制局域自旋朝向单一取向对齐（导致均匀相，即 PM）；而 $\alpha$ 则通过诱导有效的非均匀空间耦合来降低能量，从而促进非均匀调制的 SDW 有序。
$h-\Delta$ 与 $h-\alpha$ 平面相图（图 7）：
- 在填充因子分别处于 $\rho = 1/2$、 $2/3$、 $1$ （从下往上趋于半填充）时，相图展示了丰富的物理现象。特别是在 $\rho < 1$ 的分数填充区，相图中稳定地共存着 SDW-I 与 SDW-II 两个独立区域。而当系统达到半填充（ $\rho = 1$）时，由于电荷自由度受到强烈限制，自旋向下费米子的物理空间被极大压缩，SDW-II 区域随之萎缩乃至完全消亡，整个相图被 SDW-I 完全主导。

2.5 Hubbard-$U$ 排斥相互作用的物理影响 (Repulsive Hubbard-$U$ Effect)

为了评估真实凝聚态实验中必不可少的费米子原位库仑排斥对该效应的影响，作者研究了 $U > 0$ 时系统的演化（图 8）：

相区的稳定性：图 8(a) 和 8(b) 描绘了在 $h-U$ 平面中，对于 $\alpha = 0.4$ （对应 SDW-I）和 $\alpha = 0.5$ （对应 SDW-II）的相图。可以清晰地看出，即使施加极强的排斥相互作用（如 $U/t = 3.0$），SDW-I 和 SDW-II 的基态有序区依然稳固存在。这证明双模式 SDW 有序并不依赖于弱相互作用近似，而是对强关联具有极高免疫力的拓扑学鲁棒结构。
自旋极化的增强效应：随着 $U$ 的增大，向完全极化（FP）顺磁相转变的边界向下移动（即临界塞曼场 $h_c$ 随着 $U$ 的增加而显著降低）。这是因为强排斥 Hubbard-$U$ 严重惩罚了任何双占位状态（ $\uparrow\downarrow$ 构型），迫使费米子自旋自发地朝相同方向（单占位）极化以最小化电荷库仑能量。这种效应等效于给系统施加了一个额外的、由电荷关联诱导的强交换磁场，从而大大拓宽了系统的自旋极化相区。

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

为了使量子化学和凝聚态计算的科研工作者能够快速、准确地复现该论文中的多体物理结果，本节提供基于著名张量网络计算框架 ITensor 的详细代码实现细节和算法复现指南。

3.1 软件栈与计算环境准备

我们选用由 Julia 语言驱动的 ITensors.jl 软件包。Julia 的即时编译（JIT）特性与高性能 C++ 后端张量库相结合，是目前运行一维 DMRG 的最佳选择。

安装 Julia (推荐 v1.9 及以上版本)。
安装 ITensors：在 Julia REPL 中输入 using Pkg; Pkg.add("ITensors")。

3.2 动态交替格点物理对象（SiteType）的设计与映射

在 ITensor 中，无法直接在一条一维链的同一个位置放置两种不同类型的物理自由度（这里是 Fermion 和 Spin-1/2）。因此，最优雅、最高效的复现方案是构建一条交替格点链（Bipartite Chain），链的总长度设为 $2L-1$：

奇数格点（ $2i-1$ ）定义为流动费米子格点（物理类型为 "Electron" 或 "Fermion"，推荐使用带自旋的 "Electron" 以利用电荷量子数 $N_e$ 和自旋对称性 $S^z$ 守恒）。
偶数格点（ $2i$ ）定义为化学键中心的局域自旋（物理类型为 "SpinHalf"，利用自旋分量 $S^z$ 守恒）。

这种巧妙的拓扑映射不仅完美契合了一维 MPS 的一维拓扑邻近性，而且大幅降低了单个位置的局域物理维度（物理空间从 8 降至交替的 4 和 2），从而能够节省数倍的 SVD 分解计算时间。

3.3 Julia 复现核心代码：哈密顿量构建与 DMRG 求解

以下给出一个高度优化且功能完备的 Julia 脚本，用于复现论文中的 Hamiltonian 构建和 DMRG 模拟过程：

using ITensors

function run_dynamical_lattice_dmrg(L::Int, t::Float64, α::Float64, h::Float64, Δ::Float64, β::Float64, U::Float64; maxdim_max=600)
    # 1. 构建交替链格点系统
    # 奇数点为 Electron (费米子)，偶数点为 SpinHalf (局域自旋)
    N_total = 2 * L - 1
    sites = Vector{Index}(undef, N_total)
    for i in 1:N_total
        if isodd(i)
            sites[i] = siteind("Electron", i; conserve_qns=true)
        else
            sites[i] = siteind("SpinHalf", i; conserve_qns=true)
        end
    end

    # 2. 利用 AutoMPO 自动构建复杂的量子多体算符 MPO
    # 注意格点物理映射关系：
    # 真实格点 i -> 转换为 MPS 链上的奇数点 2i-1
    # 真实键 i   -> 转换为 MPS 链上的偶数点 2i
    # 真实格点 i+1 -> 转换为 MPS 链上的奇数点 2i+1
    ampo = OpSum()

    for i in 1:(L-1)
        f1 = 2*i - 1
        s_bond = 2*i
        f2 = 2*i + 1

        # --- (a) 自由费米子跃迁项: -t * (c_i^† c_{i+1} + H.c.) ---
        for σ in ["Up", "Dn"]
            ampo += -t, "Cdag$(σ)", f1, "C$(σ)", f2
            ampo += -t, "Cdag$(σ)", f2, "C$(σ)", f1
        end

        # --- (b) 局域自旋依赖跃迁项: -α * (c_i,σ^† S_i^z c_{i+1,σ} + H.c.) ---
        # 这是三体项，在 ITensor 中直接连写表示算符乘积
        for σ in ["Up", "Dn"]
            ampo += -α, "Cdag$(σ)", f1, "Sz", s_bond, "C$(σ)", f2
            ampo += -α, "Cdag$(σ)", f2, "Sz", s_bond, "C$(σ)", f1
        end

        # --- (c) 局域自旋纵向场与横向场 ---
        ampo += Δ / 2.0, "Sz", s_bond
        if β != 0.0
            # 横向场涉及自旋翻转 Sx = 0.5*(S+ + S-)
            # 注：若 conserve_qns=true，由于自旋守恒，直接添加 Sx 会破坏量子数守恒。
            # 这里为演示方便，若 β!=0，需将 SpinHalf 初始化为 conserve_sz=false
            ampo += β, "Sx", s_bond
        end
    end

    # --- (d) 流动费米子的塞曼场 h 与原位排斥 U ---
    # 作用在奇数点 (费米子格点)
    for i in 1:L
        f_idx = 2*i - 1
        # 塞曼场: h * (n_up - n_dn)
        ampo += h, "Nup", f_idx
        ampo += -h, "Ndn", f_idx
        # Hubbard 相互作用: U * n_up * n_dn
        ampo += U, "Nupdn", f_idx
    end

    H = MPO(ampo, sites)

    # 3. 设置守恒的初态并构建随机 MPS
    # 设定粒子数 N_up = N_dn = L/4 (例如无极化的 1/4 填充)
    state = Vector{String}(undef, N_total)
    n_up_count = 0
    n_dn_count = 0
    target_up = div(L, 4)
    target_dn = div(L, 4)

    for i in 1:N_total
        if isodd(i)  # 费米子点
            if n_up_count < target_up && n_dn_count < target_dn
                state[i] = "UpDn"
                n_up_count += 1
                n_dn_count += 1
            elseif n_up_count < target_up
                state[i] = "Up"
                n_up_count += 1
            elseif n_dn_count < target_dn
                state[i] = "Dn"
                n_dn_count += 1
            else
                state[i] = "Emp"
            end
        else  # 键自旋点，初始随机设置为 Up/Dn
            state[i] = isodd(i÷2) ? "Up" : "Dn"
        end
    end

    psi0 = randomMPS(sites, state, 10)

    # 4. 设定高级 DMRG 扫频规划 (Sweeps)
    sweeps = Sweeps(60)
    maxdim!(sweeps, 10, 20, 100, 200, 400, maxdim_max)
    cutoff!(sweeps, 1e-10)
    noise!(sweeps, 1e-6, 1e-7, 1e-8, 1e-10, 0.0)

    # 5. 执行 DMRG 计算
    energy, psi = dmrg(H, psi0, sweeps)

    println("计算完成! 基态能量 E = ", energy)
    return energy, psi, sites
end

3.4 测量结构因子与相关物理量的后处理

在 DMRG 求解得到高质量基态 MPS psi 后，我们需要提取局域自旋在动量空间的结构因子。下面的后处理代码展示了如何利用 ITensor 极其便捷地计算实空间关联并执行离散傅里叶变换：

function calculate_spin_structure_factor(psi, sites, L::Int)
    # 保存实空间关联矩阵
    corr_matrix = zeros(Float64, L-1, L-1)
    
    # 首先提取每个键中心自旋的 Sz 期望值 <Sz_i>
    sz_vals = zeros(Float64, L-1)
    for i in 1:(L-1)
        s_bond = 2*i
        sz_vals[i] = expectation_value(psi, "Sz", s_bond)
    end
    sz_avg = sum(sz_vals) / (L-1)

    # 计算两点关联 <Sz_i * Sz_j>
    for i in 1:(L-1)
        for j in 1:(L-1)
            s_i = 2*i
            s_j = 2*j
            if i == j
                corr_matrix[i, j] = 0.25 - sz_avg^2 # 自旋-1/2 的 Sz^2 = 1/4
            else
                # 计算关联函数
                op_i = op("Sz", sites[s_i])
                op_j = op("Sz", sites[s_j])
                # ITensor 内置的 correlation_matrix 函数适合快速批处理，
                # 此处展示单点精确提取：
                corr_matrix[i, j] = correlation_matrix(psi, "Sz", "Sz")[s_i, s_j] - sz_vals[i]*sz_vals[j]
            end
        end
    end

    # 执行空间傅里叶变换，求出 1D Brillouin 区内的 S^z(k)
    k_points = range(-π, π, length=200)
    S_k = zeros(Float64, length(k_points))

    for (idx, k) in enumerate(k_points)
        sum_val = 0.0
        for i in 1:(L-1)
            for j in 1:(L-1)
                sum_val += cos(k * (i - j)) * corr_matrix[i, j]
            end
        end
        S_k[idx] = sum_val / L
    end

    return k_points, S_k
end

通过运行上述后处理函数，绘制出不同 α 值和自旋密度下的 S_k 图像，科研人员即可完美地复现论文图 2 和图 4 中令人振奋的 SDW-I 和 SDW-II 的峰值连续演化和分裂特征。

4. 关键引用文献及该研究的局限性客观评论

4.1 关键引用文献

本研究的工作扎根于一维多体物理和拓扑晶格调控的经典理论基础之上。为了全面理解其学术脉络，以下 6 篇关键文献是必读的基石：

一维量子多体关联经典：Giamarchi, T. Quantum Physics in One Dimension (Oxford University Press, 2003). 【奠定了自旋不平衡一维费米系统 Tomonaga-Luttinger 液体理论及 SDW 关联幂律衰减分析的理论基石】
一维 Hubbard 关联理论：Schulz, H. J. Correlation exponents and the metal-insulator transition in the one-dimensional Hubbard model, Phys. Rev. Lett. 64, 2831 (1990). 【阐明了费米面嵌套（Nesting）在一维费米子系统自旋/电荷不稳定性中的本质物理机制】
DMRG 算法奠基之作：White, S. R. Density matrix formulation for quantum renormalization groups, Phys. Rev. Lett. 69, 2863 (1992). 【提供了本研究赖以实现极高计算精度的算法核心工具】
动态晶格概念起源：Cuadra, D. G., Grzybowski, P. R., Dauphin, A. & Lewenstein, M. Strongly Correlated Bosons on a Dynamical Lattice, Phys. Rev. Lett. 121, 090402 (2018). 【首次提出了流动玻色子与键中心经典/量子局域自由度动态耦合的概念】
一维自旋不平衡实验：Liao, Y.-A. et al. Spin-imbalance in a one-dimensional Fermi gas, Nature 467, 567 (2010). 【为本模型中考虑的费米子不平衡极化提供了直接的冷原子实验参考背景】
$\mathbb{Z}_2$ 规范场耦合费米子研究：Borla, U., Verresen, R., Grusdt, F. & Moroz, S. Confined Phases of One-Dimensional Spinless Fermions Coupled to $\mathbb{Z}_2$ Gauge Theory, Phys. Rev. Lett. 124, 120503 (2020). 【建立了键中心局域自旋与流动费米子之间的动力学耦合在规范场理论下的等价诠释】

4.2 研究局限性评论

虽然该项工作物理图像清晰，数值计算极其扎实，但作为一个面向物理前沿的探索，它在以下几个重要维度上仍存在明显的局限性，值得未来深入研究：

一维特殊性与向高维推广的断层：该研究所有的物理效应（特别是 SDW-I 和 SDW-II 的分裂）都高度依赖于一维系统完美的费米面“点状”特征。在一维中，费米面只是动量空间中的两个独立点，因而任何微小的扰动和嵌套都会导致极强的发散不稳定性。然而，一旦体系拓宽到二维（2D）或三维（3D）晶格，费米面将变成一维闭合曲线或二维闭合曲面。此时，局域的费米面嵌套不稳定性会被严重稀释。动态晶格介导的双模式 SDW 调制是否能在二维层状材料（如铜氧化物、铁基超导体或转角石墨烯）中幸存？这是一个巨大的问号，也是理论计算亟待突破的瓶颈。
准静态极限（Quasiadiabatic Limit）的约束：作者在文中明确指出，模拟主要集中在 $\beta \ll t$ 的准静态物理极限下（设定横向涨落场 $\beta = 0.02$）。在这种情况下，键中心局域自旋的翻转被严重抑制，自旋几乎等效于静态的非均匀磁背景。然而，最本质的动态晶格物理应当发生在完全动态区（Fully Dynamical Regime, $\beta \sim t$）。在强量子涨落场下，局域自旋会发生高频的相干翻转，其动力学行为是否会彻底抹杀 SDW-I 和 SDW-II 的长程有序？抑或会催生出高度新奇的、具有时空分数化激发的规范自旋液体？当前的静态相图未能提供关于这个强动态区域的解答。
非平衡动力学（Nonequilibrium Dynamics）的缺失：作为动态晶格，其最迷人、最适合超快泵浦探测（Pump-probe）实验的特性应当在时间演化和非平衡激发中展现。例如，通过含时超快激光脉冲猝灭（Quench）费米子极化度，局域自旋的 SDW-I 和 SDW-II 模式如何实现超快的时间域相干转换？在时间演化中是否存在非热力学非平庸稳态？目前的工作仅局限于静态基态分析，缺乏时域动力学的刻画。
实验可行性挑战：尽管冷原子和里德堡（Rydberg）原子晶格提供了强大的模拟潜力，但要完美实现哈密顿量中第二项涉及流动费米子自旋、格点位置和中间化学键自旋的三体动力学耦合项 $\hat{c}_{i,\sigma}^\dagger \hat{S}_i^z \hat{c}_{i+1,\sigma}$，在实验设计上具有极高的门槛。这需要极其精密的辅助拉曼光路或复杂的里德堡能级魔术调控（Magic wavelength tuning），目前尚未有成熟、普适的成熟方案能无损地实现这种自旋依赖性跃迁调制。

5. 补充探讨：动态晶格与规范场论、量子链接模型的深层融合

为了帮助读者将本研究置于更宏大的近代凝聚态与高能物理交叉学科视角中，我们有必要深入探讨其背后的**规范场理论（Gauge Field Theories）与量子链接模型（Quantum Link Models, QLM）**的深层对应关系。

5.1 动态晶格中的“局部涌现规范对称性”

在粒子物理和格点规范场论（Lattice Gauge Theory）中，电荷的移动必须伴随着规范玻色子（如光子或胶子）在化学键上的产生或湮灭，以严格遵循电荷守恒（即高斯定律）。

在本论文考虑的哈密顿量中，如果我们考虑无自旋费米子在 $\alpha = t$ 极限下的行为，跃迁算符变为：

$$\hat{c}_i^\dagger (\frac{1}{2} - \hat{S}_i^z) \hat{c}_{i+1}$$

这表明费米子只有在键上局域自旋处于向下（$-1/2$）的状态时才能进行跃迁，并在跃迁后，局域自旋的状态与费米子的密度分布发生纠缠。这在数学上等价于将局域自旋 $\hat{S}_i^z$ 视为定义在格点键上的 $\mathbb{Z}_2$ 量子规范场 $A_{i,i+1}$。费米子的跃迁项本质上构成了一个一维的 $\mathbb{U}(1)$ 或者是 $\mathbb{Z}_2$ 量子链接模型。高斯定律在每一个格点上涌现，将电荷密度与相邻键上的通量差强力锁定：

$$\hat{G}_i = \hat{S}_i^z - \hat{S}_{i-1}^z - \hat{n}_i = \text{const}$$

因此，该工作所发现的双模式 SDW 现象，其物理本质不仅是简简单单的低维磁性，更是在自旋不平衡费米子调制下，格点规范场中的量子链接流发生的一种新奇的‘局域规范场凝聚不稳定性’（Gauge condensates instability）。这一视角将一维冷原子多体模拟的成果，直接嫁接到了高能格点规范场论的非微扰现象研究中，展现了多学科大交叉的非凡魅力。

5.2 实验方案展望：里德堡 dressed 原子的量子模拟

为了克服第 4 节中提到的实验挑战，这里提供一个极具前景的、基于里德堡 dressed 费米原子的实验模拟架构建议：

载体选择：选择超冷费米子碱土金属原子（如 $^{173}\text{Yb}$ 或 $^{87}\text{Sr}$）。利用其极其长寿命的亚稳激发态（“钟”态）作为费米子的自旋向上和向下状态。
动态键构建：将原子捕获在精心设计的、间距可调的交替光晶格中。通过引入周期性振荡晶格势（Lattice shaking, 弗洛凯工程），将常规跃迁 $t$ 压制到零。
三体跃迁的实现：利用一束紧聚焦的紫外激光，将化学键中心的虚拟能级弱耦合到具有极强范德华相互作用的里德堡态（Rydberg-dressed states）。此时，费米子穿过中间化学键的能垒高度将直接取决于邻近原子是否处于特定的里德堡激发态。这就巧妙地将空间跃迁振幅与中间化学键的局域量子状态锁定，在实验室中真实还原该动态晶格哈密顿量中的自旋依赖跃迁耦合项。

5.3 结语

山西大学研究团队的这项工作通过极高精度的数值手段，将一维自旋不平衡多体物理推进到了前沿的动态晶格范畴。其发现的基于 Fermi 嵌套的分裂式双模式自旋密度波（SDW-I 和 SDW-II），为人们调控低维材料的非平庸磁调制周期提供了崭新的维度。随着多体数值方法（如张量网络、变分蒙特卡洛）以及冷原子模拟技术的不断跃升，我们有理由相信，动态晶格磁性研究将成为连接强关联材料磁性与非微扰规范场论的一座宏伟桥梁。