来源论文: https://arxiv.org/abs/2606.30176v1 生成时间: Jul 04, 2026 07:03

扩展哈伯德模型中三聚体链的磁性与相分离：Lieb定理的稳定性及DMRG深度解析

0. 执行摘要

在强关联凝聚态物理领域中，寻找并稳定具有非平庸磁序的基态是一项长期而艰巨的任务。Lieb定理（Lieb’s Theorem）作为量子磁学理论的奠基石之一，指出在半满（Half-filling）状态下的双分格点（Bipartite Lattice）上，若仅存在格点内部的排斥力 $U$，且子格点数目不相等（$N_B \neq N_A$），系统将展现出总自旋为 $S = |N_B - N_A|/2$ 的亚铁磁（Ferrimagnetic）基态。然而，真实的物理系统往往存在非局域的库仑相互作用。最邻近库仑排斥力 $V$ 的引入是否会迅速破坏这一精妙的磁性基态，以及它如何诱导电荷密度波或相分离，一直是理论物理学家关注的核心课题。

本研究针对一维双分耦合三聚体链（Coupled Trimer Chain）上的扩展哈伯德模型（Extended Hubbard Model）展开了系统性的研究。利用高精度的密度矩阵重正化群（Density Matrix Renormalization Group, DMRG）方法，本文系统地绘制了该模型在半满状态下的基态相图。研究结果表明：

Lieb亚铁磁相的鲁棒性：尽管最邻近库仑相互作用 $V$ 会显著增加系统内的双占据（Doublon）密度，但在达到相分离（Phase Separation, PS）线之前，Lieb亚铁磁绝缘相表现出惊人的鲁棒性。
相分离临界线：相分离发生在 $V \gtrsim U/4$ 的临界区域，这一临界值与原子极限（Atomic Limit）下的简并能量分析高度吻合。
多相共存行为：在相分离区以上，系统分裂为两个空间共存的相。其中一相是由 $B$ 子格点上的双占据主导的电荷绝缘相；而另一共存相在强关联区（$U \gtrsim t$）表现为金属化的、未饱和的铁磁相（Unsaturated Ferromagnetic Phase），在弱关联区（$U \lesssim t$）则表现为单态（Singlet）绝缘相。

本博客将从核心物理问题出发，详尽推导其理论边界，剖析DMRG的数值实现，并结合前沿文献对该工作的局限性与未来应用前景进行深度评论。

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

1.1 核心科学问题：非局域相互作用对Lieb定理的冲击

经典的哈伯德模型仅考虑格点上的原位库仑排斥作用 $U$：

$$H_{Hubbard} = -\sum_{\langle i,j \rangle, \sigma} t_{ij} (c^{\dagger}_{i\sigma} c_{j\sigma} + H.c.) + U \sum_{i} n_{i\uparrow} n_{i\downarrow}$$

当系统处于半满（每个格点平均有一个电子）时，若晶格具有双分拓扑结构（即格点可分为 $A$ 和 $B$ 两套子格点，且跳跃仅发生在 $A$ 和 $B$ 之间），Lieb定理保证了其基态具有宏观总自旋。对于本文研究的三聚体链，每个原胞（Trimer Unit）包含三个格点，标记为 $B_1 - A - B_2$。如图1所示，格点 $A$ 作为中心，仅与本单元及相邻单元的 $B_1$、$B_2$ 格点发生跳跃耦合。因此，$A$ 构成了子格点 $A$（每个原胞1个格点，$N_A = L$），而 $B_1, B_2$ 构成了子格点 $B$（每个原胞2个格点，$N_B = 2L$）。整个系统在半满时的总格点数为 $N = 3L$，电子数为 $N_e = 3L$。根据Lieb定理，当 $V = 0, U > 0$ 时，基态总自旋严格为：

$$S = \frac{|N_B - N_A|}{2} = \frac{2L - L}{2} = \frac{L}{2}$$

即每个原胞贡献 $1/2$ 的自旋，系统呈现亚铁磁长程序。

然而，当引入最邻近库仑相互作用 $V \sum_{\langle i,j \rangle} n_i n_j$ 时，相互作用能不仅取决于局域的双占据，还取决于相邻格点之间的电荷排布。这一非局域项不仅打破了原有单粒子算符在子格点变换下的对称性，还倾向于使电荷空间分布不均匀，从而导致电荷密度波（CDW）或相分离（PS）。核心科学问题在于：Lieb亚铁磁相在面对 $V$ 的扰动时能坚持多久？其瓦解的微观机制是连续的电荷转移还是急剧的一阶相分离？

1.2 理论基础：原子极限下的严格边界推导

为了直观理解这一竞争，我们首先在跳跃振幅 $t = t' = 0$ 的原子极限（Atomic Limit）下，分析几种具有代表性的静态电荷排布结构，从而为相互作用驱动的相变建立精确的解析基准。

考虑三种典型的静态构型（如图2所示）：

构型 (i)：标准Lieb亚铁磁电荷态 每个格点上恰好占据一个电子（无双占据，即双占据数 $\langle d \rangle = 0$）。此时，原位相互作用 $U$ 的贡献为 0。每个原胞内，中心格点 $A$ 与四个最邻近的 $B$ 格点（同单元内的 $B_1, B_2$，以及相邻单元耦合的两个 $B$ 格点）相连。因此，每个原胞产生的最邻近排斥能为 $4V$。每个原胞的总能量为：
$$E_{(i)}/L = 4V$$
构型 (ii)：部分电荷转移态 原胞内的 $A$ 格点被彻底排空（$\langle n_A \rangle = 0$），而其多余的电荷转移到 $B$ 子格点上，导致其中一个 $B$ 格点双占据（例如 $B_1$ 占据2个电子），另一个 $B$ 格点单占据（$B_2$ 占据1个电子）。在此构型下，由于 $A$ 格点为空，所有连接 $A$ 与 $B$ 的最邻近键上的电荷密度乘积 $n_A n_B = 0$。因此，最邻近相互作用能降为 0。但代价是产生了一个双占据，贡献了原位库仑能 $U$。每个原胞的总能量为：
$$E_{(ii)}/L = U$$
构型 (iii)：最大双占据态 $A$ 格点依然为空（$\langle n_A \rangle = 0$），而 $B$ 子格点上的电荷达到最大双占据。这意味着在 $2L$ 个 $B$ 格点上均匀分布了 $3L$ 个电子，从而产生 $3L/2$ 个双占据，平均每个原胞包含 $3/2$ 个双占据。最邻近排斥能同样为 0。每个原胞的总能量为：
$$E_{(iii)}/L = \frac{3}{2}U$$

通过比较构型 (i) 和构型 (ii) 的能量，我们可以立刻发现：当 $E_{(i)} > E_{(ii)}$，即 $4V > U$ 时，系统将倾向于将电荷从 $A$ 格点剥离并注入到 $B$ 格点，以消除巨大的最邻近排斥能。这给出了相不稳定性发生的临界值：

$$V_c^{atomic} = \frac{U}{4}$$

这一极为简洁的原子极限推导，不仅预言了不稳定性发生的精确物理位置，而且指明了这一相变伴随着电荷由 $A$ 向 $B$ 的转移。然而，在有跳跃项（$t \neq 0, t' \neq 0$）的真实体系中，量子涨落和动能项会显著修正这一边界，必须依赖精确的数值计算手段。

1.3 技术难点与DMRG计算策略

在一维强关联系统中，DMRG是公认最精确的无符号问题（Minus-sign-problem-free）数值方法。然而，在本模型中应用DMRG面临以下三大技术瓶颈：

局域极小值与亚稳态陷阱：由于相分离相对应于实空间中两种截然不同的电荷密度区域的共存（一侧极富 doublon，另一侧极贫），在重正化扫频（Sweeping）过程中，若初始波函数选择不当，DMRG极易陷入亚稳态（Metastable States）。特别是在 $V \approx U/4$ 的临界点附近，不同实空间电荷排布的态能量极其接近。
- 对策：作者采用了两种重正化初始化策略。第一种是标准的随机自旋排布（具有 $S^z = 0$ 约束）；第二种是人工构造的非均匀电荷排布（例如将左侧 $3L/4$ 的格点设为双占据，右侧排空）。通过对比两种方法在数十次扫频后的基态能量，确保获取真正的全局基态。
可压缩性的实空间表征：相分离的经典特征是电荷可压缩性 $\kappa \to \infty$。在有限尺寸晶格中，利用能量的二阶差分定义局部可压缩性会带来巨大的边界效应。定义如下：
$$\frac{1}{n^2\kappa} = \frac{N}{4} [E(N_e + 2) + E(N_e - 2) - 2E(N_e)]$$
由于相分离区域存在畴壁（Domain Wall），边界和畴壁的相互作用使得可压缩性的数值表现出高度的波动性。
- 对策：引入实空间局部电荷偏差最大值 $\Delta_{max}$。通过监测偏离均值（3个电子/原胞）的最大局域电荷涨落，精确标定相分离的实空间边界，并进行系统的有限尺寸外推（Finite-size Scaling）。
计算对称性破缺：本模型中的两个跳跃项 $t$ 与 $t'$ 若完全相等，会引入高度的局域对称性。这在DMRG扫频时可能导致局域纠缠熵异常增加，进而减缓收敛速度。为了破坏这种非物理的局域对称性，作者特意选取了不对称的跳跃振幅组合：$t = 1.0, t' = 0.8$。

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

本节深度剖析研究中三个关键 Benchmark 体系的 DMRG 计算结果，这些数据是勾勒完整相图的基石。

2.1 Benchmark 1: $U = 0$ 的非相互作用/弱关联极限 ($t-t'-V$ 模型)

当原位排斥力 $U = 0$ 时，系统仅由最邻近排斥力 $V$ 驱动。直觉上，任何微小的 $V$ 都会促使电子避免占据最邻近格点。

电荷与双占据分布（见图3）：
- 在 $V = 0$ 时，各格点上的平均电子数严格为 $\langle n_A \rangle = \langle n_{B1} \rangle = \langle n_{B2} \rangle = 1.0$。此时双占据 $\langle d \rangle$ 受动能驱动，保持在一个常数水平。
- 随着 $V$ 增加至 $V = 4$，$\langle n_A \rangle$ 迅速下降至 $\approx 0.25$，而 $B$ 子格点的总电荷 $\langle n_B \rangle = \langle n_{B1} \rangle + \langle n_{B2} \rangle$ 增加至 $\approx 2.75$。与之相伴的是 $B$ 上的双占据密度 $\langle d_B \rangle$ 陡增至 $\approx 1.25$。
- 结论：在 $U=0$ 极限下，系统没有展现出一阶相变，而是展现出一个由 $V$ 驱动的连续电荷转移过程。然而，该状态下的实空间分布（图4a-d）表明，当 $V > 1.0$ 时，系统在空间上已经分裂成富 doublon 的绝缘相区域（平均每原胞4个电子）与富空穴的金属相区域（平均每原胞2个电子）。
可压缩性外推（见图4e）：
- 通过对 $L = 20, 40, 60$ 进行有限尺寸外推，发现当 $V > 1.0$ 时，反比可压缩性 $1/\kappa$ 迅速降为 0，这确证了 $U=0$ 系统在热力学极限下对于 $V > 0$ 也是电荷不稳定的，一经扰动即发生相分离。

2.2 Benchmark 2: $U \neq 0$ 的有限关联体系——Lieb相的崩溃与电荷不均匀性

当引入有限的 $U$ 后，原位排斥会强烈抑制双占据的产生，从而与 $V$ 的电荷转移效应产生激烈竞争。作者选取了 $U = 2.0, 4.0, 8.0, 16.0, 32.0$ 等典型强度进行详尽扫频。

关联强度 $U$	相分离临界点 $V_c$ (DMRG)	原子极限预测值 $U/4$	Lieb相区间 (鲁棒区)	磁性行为 ($V > V_c$)
2.0	$\approx 0.52$	0.50	$V < 0.52$	单态 (Singlet, $S/L \approx 0$)
4.0	$\approx 1.08$	1.00	$V < 1.08$	亚铁磁与铁磁共存 ($S/L \approx 0.3$)
8.0	$\approx 2.15$	2.00	$V < 2.15$	未饱和铁磁与单态共存 ($S/L \gtrsim 0.3$)
16.0	$\approx 4.20$	4.00	$V < 4.20$	强铁磁有序

数据深度解读：
- 精确相变点的偏离：DMRG 计算得到的 $V_c$ 始终略高于原子极限预测的 $U/4$。例如当 $U=8.0$ 时，$V_c \approx 2.15$。这一微小的正向偏离完全来自于系统的动能项（$t, t'$ 引起的量子跳跃）。量子动能倾向于使电荷局域化均匀分布，从而对相分离起到了稳定和延缓作用。
- 相分离区的实空间波函数特征（图7）：在 $U=8.0, V=4.0$（即 $4V/U = 2.0$）的强相分离区中，DMRG精确还原了两个共存相的实空间特征：
  1. 黄色区域 (富电荷相)：$\langle n_A \rangle \approx 0.05$，$\langle n_B \rangle \approx 2.0$，且双占据 $\langle d_B \rangle \approx 1.0$。这完美对应了静态构型 (iii)，即电荷在 $B$ 子格点上完全两两配对，形成电荷绝缘体，其局域总自旋为 0。
  2. 蓝色区域 (贫电荷相)：$\langle n_A \rangle < 0.2$，$\langle n_B \rangle \approx 1.33$，双占据 $\langle d_B \rangle \approx 0.36$。由于电荷不满足整数填充，此区域具有极高的电荷迁移率，表现出金属特征，且其局域自旋由于双分结构的破坏未被完全淬灭，展现出未饱和的铁磁有序。

2.3 Benchmark 3: 自旋长程序与磁相图的定量标定

为了判定磁性转变，作者通过计算自旋-自旋关联函数定量提取了总自旋平方 $S(S+1)$（式6）：

$$S(S+1) = \sum_{i,j} \langle \mathbf{S}_i \cdot \mathbf{S}_j \rangle$$

自旋密度的空间演化（图8b）：
- 在整个均匀相区（$V < V_c$），无论 $V$ 如何增加，每个原胞的平均自旋 $S/L$ 始终坚挺地保持在严格的 $1/2$。这表明即使最邻近排斥力 $V$ 已经在不断把 $A$ 格点的电荷往 $B$ 格点上“驱赶”，只要系统没有整体崩溃发生相分离，量子纠缠仍能维持经典的 Lieb 亚铁磁长程序。这一鲁棒性超出了以往学术界的普遍预期。
- 一旦跨越相分离线（$V > V_c$），对于弱关联 $U \lesssim 1.0$，总自旋急剧跌落至 $S/L \approx 0$（表现为非磁性单态共存）；而对于强关联 $U \gtrsim 1.0$，由于相分离产生的金属区内存在强未饱和铁磁矩，系统总自旋保持在 $S/L \gtrsim 0.3$ 的高位。

3. 代码实现细节，复现指南，所用的软件包及开源 repo link

本研究的数值计算完全基于著名的开源张量网络库 ITensor。ITensor 提供了高效的矩阵乘积态（MPS）和矩阵乘积算符（MPO）数据结构，非常适合一维和准一维量子多体系统的精确求解。以下提供基于 ITensor Julia 版本的完整复现指南。

3.1 晶格映射与哈密顿量构建

在一维 DMRG 中，必须将三聚体链的 2D 拓扑结构映射为一维链。对于包含 $L$ 个原胞的链，我们将其映射为包含 $3L$ 个格点的一维 MPS。原胞 $l$ ($l = 1, \dots, L$) 的内部格点映射关系为：

$B_{1, l} \to 3l - 2$
$A_l \to 3l - 1$
$B_{2, l} \to 3l$

基于此映射，哈密顿量中的跳跃项和相互作用项可以精确写为一维链上格点索引的形式。

3.2 完整复现代码（Julia + ITensor）

请在 Julia 环境下安装 ITensor（using Pkg; Pkg.add("ITensors")）后运行以下脚本：

using ITensors

function run_trimer_dmrg(L::Int, U::Float64, V::Float64; t::Float64=1.0, t_prime::Float64=0.8)
    N = 3 * L  # 总格点数
    sites = siteinds("Electron", N; conserve_qns=true) # 守恒粒子数和总 Sz

    # 1. 构建哈密顿量算符和 (OpSum)
    os = OpSum()

    for l in 1:L
        # 确定一维格点索引
        b1 = 3*l - 2
        a  = 3*l - 1
        b2 = 3*l

        #--- 原位相互作用 U ---
        os += U, "Nupdn", b1
        os += U, "Nupdn", a
        os += U, "Nupdn", b2

        #--- 内跳跃项 t (A_l 与 B1_l, B2_l) ---
        for sigma in ["Up", "Dn"]
            os += -t, "Cdag$sigma", a, "C$sigma", b1
            os += -t, "Cdag$sigma", b1, "C$sigma", a
            os += -t, "Cdag$sigma", a, "C$sigma", b2
            os += -t, "Cdag$sigma", b2, "C$sigma", a
        end

        #--- 内最邻近排斥 V (A_l 连 B1_l, B2_l) ---
        os += V, "Ntot", a, "Ntot", b1
        os += V, "Ntot", a, "Ntot", b2

        #--- 跨原胞项 (仅在非边界处) ---
        if l < L
            next_b1 = 3*(l+1) - 2
            next_a  = 3*(l+1) - 1

            # 跨原胞跳跃 t'
            for sigma in ["Up", "Dn"]
                os += -t_prime, "Cdag$sigma", a, "C$sigma", next_b1
                os += -t_prime, "Cdag$sigma", next_b1, "C$sigma", a
                os += -t_prime, "Cdag$sigma", b2, "C$sigma", next_a
                os += -t_prime, "Cdag$sigma", next_a, "C$sigma", b2
            end

            # 跨原胞最邻近排斥 V
            os += V, "Ntot", a, "Ntot", next_b1
            os += V, "Ntot", b2, "Ntot", next_a
        end
    end

    H = MPO(os, sites)

    # 2. 设置半满初始化状态 (3L 个电子，总 Sz = L/2 对应 Lieb 相，或 Sz = 0 对应测试相分离)
    # 为了探寻真实的基态，我们初始化一个具有 Sz = 0 且在空间上分布均匀的态
    state = String[]
    for i in 1:N
        if mod(i, 2) == 1
            push!(state, "Up")
        else
            push!(state, "Dn")
        end
    end
    psi0 = randomMPS(sites, state, 10)

    # 3. 设置DMRG扫频参数
    sweeps = Sweeps(15)
    setmaxdim!(sweeps, 50, 100, 200, 400, 800, 1200, 2000)
    setcutoff!(sweeps, 1E-8)
    # 引入噪声项 (Noise term) 极其关键，能有效防止 DMRG 锁死在局域相分离的亚稳态中
    setnoise!(sweeps, 1E-5, 1E-6, 1E-7, 1E-8, 0.0)

    # 4. 运行DMRG
    energy, psi = dmrg(H, psi0, sweeps)

    # 5. 测量物理量：实空间电荷密度与双占据数
    n_profile = Float64[]
    d_profile = Float64[]
    for i in 1:N
        # 测量电荷数 <n_i>
        n_val = expect(psi, "Ntot"; sites=i)
        push!(n_profile, n_val)
        # 测量双占据 <d_i> = <n_i,up * n_i,dn>
        d_val = expect(psi, "Nupdn"; sites=i)
        push!(d_profile, d_val)
    end

    println("Calculated ground state energy per unit: ", energy / L)
    return n_profile, d_profile, psi
end

# 示例运行：U=8.0, V=1.0 (处于鲁棒的 Lieb 亚铁磁相)
n_prof, d_prof, ground_psi = run_trimer_dmrg(40, 8.0, 1.0)

3.3 数值收敛监控要点

在复现相分离区域（如 $U = 8.0, V = 3.0$）时，需特别监控以下指标：

最大截断误差（Truncation Error）：在最后几次扫频中，截断误差需严格小于 $10^{-6}$。若截断误差显著偏高，需将最大保留维数（maxdim）提升至 3000 或更高。
电荷畴壁的漂移：在相分离状态下，由于体系具有平移对称性，电荷富集区可能位于晶格的左侧、右侧或中间。若发现多次独立计算得到的电荷密度分布图形状不一，但总能量极度接近，此乃正常现象，表明相分离共存态的畴壁能极低。

开源库链接：
- ITensor 官方 Github 仓库
- ITensor 官方文档与 C++/Julia 指南

4. 关键引用文献，以及你对这项工作局限性的评论

4.1 关键历史文献背书

本项研究工作立足于数十年强关联多体理论的演进，其学术谱系可追溯至以下里程碑式文献：

Lieb定理的奠基：
- Lieb, E. H. Two theorems on the Hubbard model. Phys. Rev. Lett. 62, 1201 (1989).
- 贡献：严格证明了半满下双分晶格哈伯德模型的独特亚铁磁基态，为量子磁性提供了非微扰的解析基础。本文的核心立论均基于此定理。
DMRG 方法的创立：
- White, S. R. Density matrix formulation for quantum renormalization groups. Phys. Rev. Lett. 69, 2863 (1992).
- 贡献：克服了传统量子蒙特卡洛（QMC）的负符号问题，使一维强关联体系的精确求解成为可能。
$AB_2$链模型的早期探索：
- Macédo, A. M. S., dos Santos, M. C., Coutinho-Filho, M. D. Magnetism and Phase Separation in Polymeric Hubbard Chains. Phys. Rev. Lett. 74, 1851 (1995).
- 贡献：首次探讨了双分高分子链哈伯德模型的相分离倾向，是本文研究三聚体扩展哈伯德模型的直接理论前身。

4.2 局限性深度剖析与学术批判

尽管本文通过精湛的DMRG数值模拟，全面刻画了三聚体链的相分离和磁性不稳定性，但从前沿强关联材料物理的角度审视，该工作仍存在以下不容忽视的局限性：

1. 维度灾难与真实材料的跨越（Dimensionality Gap）

本文所探讨的物理模型是一维耦合三聚体链。然而，自然界中实际存在的、具有三聚体结构的量子磁性材料，如：

顺磁体 $\text{Cu}_3(\text{P}_2\text{O}_6\text{OD})_2$
著名的蓝铜矿化合物蓝铜矿 $\text{Cu}_3(\text{CO}_3)_2(\text{OH})_2$
磷酸盐系统 $A_3\text{Cu}_3(\text{PO}_4)_4$ ($A = \text{Ca, Sr, Pb}$)

这些材料在极低温度下，链间的弱三维耦合（Inter-chain Coupling）会发生不可忽视的作用。在一维体系中，根据 Mermin-Wagner 定理，在有限温度下不存在真正的长程序，相分离和磁有序仅能在零温基态下稳定存在。三维耦合的存在可能会极大地拓宽或压缩相分离区，甚至在有限温度下诱导出截然不同的磁结构（例如非共线螺旋磁序）。因此，将一维DMRG结论直接套用于三维真实材料时需持高度审慎态度。

2. 最邻近相互作用近似的过分简化（Truncation of Coulomb Tail）

扩展哈伯德模型将库仑相互作用截断在最邻近项 $V$。但在真实的过渡金属氧化物或低维有机超导体中，库仑相互作用的屏蔽效应往往不完全，表现为具有缓慢衰减尾部的长程库仑势（如欧姆或 Yukawa 型势，包含次邻近 $V_2$、三邻近 $V_3$ 等）：

$$V_{lr} \sim \sum_{i \neq j} \frac{V}{|i-j|}$$

已有研究表明，长程库仑相互作用会强烈抑制一阶相分离（一阶相分离伴随着电荷的宏观聚集，这会带来巨大的长程静电能惩罚），转而稳定各种微观的条纹相（Stripes）或电荷晶体态（Wigner Crystal）。因此，本文观察到的宏观相分离在更真实的电荷相互作用下可能会“崩塌”为精细的纳米级电荷调制结构。

3. 动力学谱函数证据的缺失（Absence of Dynamical Spectra）

本工作完全聚焦于基态静态关联（自旋密度、双占据、静态可压缩性）。然而，实验探测手段如角分辨光电子能谱（ARPES）和非弹性中子散射（INS）直接测量的是动力学关联函数，即单粒子谱函数 $A(k, \omega)$ 和自旋动力学结构因子 $S(k, \omega)$。

缺失的科学信息：在进入相分离区前，Lieb 相中由于 $V$ 诱导的电荷转移是否会导致能隙（Charge Gap）的变窄？在相分离区中，共存的金属相和绝缘相是如何在自旋动力学谱上叠加并表现出独特的“双能隙”特征的？缺乏动力学谱计算，使得该工作与现代前沿谱学实验的对话能力大打折扣。

5. 其他必要补充：物理机制的深度洞察与未来展望

为了给强关联领域的科研工作者提供最透彻的物理图景，本节对三聚体链中的电荷转移机制及亚铁磁鲁棒性的微观起源进行深度物理剖析，并展望其在量子器件领域的应用。

5.1 亚铁磁鲁棒性的微观本质：动能驱动的有效超交换

为什么在 $V$ 已经开始将电子从 $A$ 挤向 $B$（导致 $\langle n_A \rangle$ 显著减小）时，Lieb 亚铁磁序仍能完好地保持，甚至没有出现自旋大小的衰减？

这可以通过**动能驱动的有效超交换（Kinetic-energy-driven Effective Superexchange）**来理解。在哈伯德模型中，半满下的磁序本质上是电子为了最大化其虚跳跃空间（从而降低动能）而产生的。在三聚体单元中，空着的 $A$ 格点或高度占据的 $B$ 格点看似削弱了局域的单粒子磁矩，但在量子力学层面上：

即使 $\langle n_A \rangle < 1.0$，由于 $t$ 和 $t'$ 项的存在，电子在 $A$ 和 $B$ 子格点之间处于强烈的相干相干相干相干相干跃迁状态。
$A$ 格点作为超交换媒介的角色没有变。只要系统保持空间均匀性，虚跃迁过程 $B_1 \to A \to B_2$ 所建立的有效反铁磁超交换相互作用 $J_{eff} \sim t^2/U$ 依然主导着整个晶格。
$B_1$ 和 $B_2$ 上的自旋方向相同，与 $A$ 上的自旋方向相反。当电子发生部分转移，这些自旋磁矩在空间上进行了量子重新分配，但其自旋配对的位相相干性（Phase Coherence）完好无损。因此，直到一阶相分离彻底撕裂晶格的均匀性之前，整体的 Lieb 亚铁磁拓扑序被量子相干牢牢地“锚定”住。

5.2 实验观测的指路明灯：如何验证相分离？

对于致力于合成低维晶格的实验学家，本研究提供了一套清晰的物理量演化特征，可用于鉴别晶格是否发生了相分离：

红外光谱中的光学导电率 (Optical Conductivity, $\sigma(\omega)$)：
- 在均匀 Lieb 绝缘相中，$\sigma(\omega)$ 表现为单一的绝缘体能隙吸收边。
- 进入相分离区后，由于金属相与绝缘相的共存，光谱将同时展现出代表金属区自由电子的德鲁德峰 (Drude Peak) 以及代表富电荷绝缘区的中红外激发峰 (Mid-infrared Peak)。这种独特的双峰共存是空间不均匀相分离的直接证据。
核磁共振 (NMR) 谱线分裂：
- 在空间均匀的亚铁磁相中，由于每个原胞等价，NMR 探测到的局部磁场（Knigth Shift）应当是离散且对称的。
- 一旦系统跨入相分离区，实空间被分割成电荷密度截然不同的两个宏观区域，NMR 谱线将发生显著的分裂，直接映射出两个不同局部环境下的电荷/自旋密度丰度。

5.3 未来展望：拓扑边缘态与量子器件设计

三聚体链由于其天然的亚晶格非对称性，往往伴随着非平庸的拓扑边界态（Topological Edge States），在超-Su-Schrieffer-Heeger (super-SSH) 晶格中有着广泛的应用。本项研究揭示的“相互作用调控电荷排布”机制，为设计拓扑自旋开关器件开辟了全新道路：

自旋开关设计概念：通过施加外加局部栅极电压（等效于调控最邻近相互作用 $V$ 或局部势能），我们可以将系统精确地控制在 Lieb 均匀相与不均匀相分离相的临界点边缘。微小的电场扰动即可诱导系统在“宏观亚铁磁绝缘态”与“局部非磁性单态 + 金属导电态”之间发生巨幅跃变。这为开发新型非易失性拓扑自旋电子学器件（Spintronic Devices）提供了理论上可行的崭新方案。

未来，随着硅基人工冷原子晶格（Cold Atoms in Optical Lattices）和纳米自旋链组装技术（SP-STM）的蓬勃发展，在高度可控的实验平台上复现并操纵本文揭示的 Lieb 亚铁磁相及其相分离不稳定性，必将成为量子模拟领域的一颗璀璨明珠。