来源论文: https://arxiv.org/abs/2607.08278v1 生成时间: Jul 10, 2026 23:47

关联诱导的奇迹：从非局域库仑相互作用中涌现的 1D 拓扑相——D-TRILEX 方法深度解析

0. 执行摘要

拓扑物态通常是在非相互作用或弱相互作用框架下进行分类和研究的，此时拓扑不变量由单粒子 Bloch 哈密顿量的本征态直接决定。然而，在强关联多体系统中，电子与电子之间的强库仑排斥会产生丰富的集体激发和自发对称性破缺。一个极其重要且具有挑战性的科学问题是：电子关联本身能否从一个完全拓扑平庸的非相互作用能带出发，自发地产生出非平庸的拓扑物态？

最近发表于 arXiv:2607.08278v1 的工作《Emergent Topology from Nonlocal Electronic Correlations in One Dimension》（作者：Félix Fossati, Erik Linner, 并在强关联多体物理领域深耕的 Evgeny A. Stepanov）给出了令人振奋的肯定回答。该研究以经典的一维扩展哈伯德模型（1D Extended Hubbard Model）为研究对象，采用了一种最近开发的先进团簇图符多体方法——双三重不可约局部展开（D-TRILEX，Dual Triply Irreducible Local Expansion）团簇扩展方法，系统地研究了非局域相互作用 $V$ 驱动的键序波（BOW）相和电荷密度波（CDW）相。研究表明，这些由相互作用驱动和自发对称性破缺产生的物相，在低能极限下可以通过一个有效的单粒子描述（拓扑哈密顿量，Topological Hamiltonian）来精确刻画。具体而言，自发对称性破缺的 BOW 相可以精确地映射到经典的 Su-Schrieffer-Heeger (SSH) 模型，而伴随有次领头键序波关联的 CDW 相（即新发现的 CDW+BOW 混合相）则可以映射到经典的 Rice-Mele (RM) 模型。两者的低能有效单粒子能谱在开放边界条件下都显现出强局域化的边缘态，从而无可辩驳地证实了非局域电子关联能够直接自发诱导涌现拓扑相。这一发现为寻找关联拓扑量子物态开辟了全新的理论视角。

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

1.1 核心科学问题

传统的拓扑物态分类依赖于单粒子 Bloch 态的几何相位（如 Berry 相位、Zak 相位）。在引入微弱相互作用时，只要能带间隙不闭合，拓扑特征通常是稳健的。然而，在强关联极限下，尤其是一维（1D）系统中，剧烈的量子涨落和电荷/自旋自由度的强耦合会导致单粒子能带图像完全失效。以往关于强关联系统中拓扑物态的研究，大多集中在相互作用如何调制、破坏或稳定那些在非相互作用极限下就已经存在的拓扑相（例如，研究 Hubbard 相互作用 $U$ 对 SSH 模型或 Kane-Mele 模型中拓扑边缘态的抑制或重整化）。

相反，本研究聚焦于一个本质上更具颠覆性的机制：从一个完全平庸（即不具备任何空间交替跃迁或势垒调制、单粒子哈密顿量平庸）的一维紧束缚格点模型出发，仅通过调节电子间的非局域排斥力，系统是否能自发产生拓扑相？ 这一问题的难点在于，要同时处理自发对称性破缺（产生非均匀的电荷分布或跃迁重整化）和强多体关联效应，而传统的平均场理论在低维系统中往往因为低估涨落而得出定性错误的结论。

1.2 理论基础：模型哈密顿量与拓扑映射

研究的基石是一维扩展哈伯德模型，其哈密顿量为：

$$H = -t \sum_{\langle i, j \rangle, \sigma} c^{\dagger}_{i,\sigma} c_{j,\sigma} + U \sum_{j} n_{j,\uparrow} n_{j,\downarrow} + \frac{V}{2} \sum_{\langle i, j \rangle, \sigma, \sigma'} n_{i,\sigma} n_{j,\sigma'}$$

其中 $t=1$ 为最近邻格点间的跃迁系数，$U$ 为局域库仑排斥，$V$ 为最近邻非局域库仑排斥，$\sigma$ 代表电子自旋。在半满（half-filling）条件下，该模型随着 $U$ 和 $V$ 的竞争，展现出丰富的相图。当 $V$ 占主导时，系统倾向于克服 $U$ 的局域单态排斥，自发形成不同的有序相：

  1. 电荷密度波相（CDW）：电荷在奇偶格点上交替排列，自发破缺了平移对称性。
  2. 键序波相（BOW）：跃迁振幅发生自发二聚化（dimerization），形成交替的强键和弱键，自发破缺了空间反演对称性。

拓扑不变量的定义和边缘态的涌现可以通过与以下两个经典的单粒子拓扑模型进行映射来理解：

  • SSH 模型：描述了一维二聚化格点，其单粒子跃迁在 $t_1$ 和 $t_2$ 之间交替。当 $t_1 < t_2$ 时，系统处于拓扑非平庸相，具有 Zak 相位 $\pi$，在开放边界下存在零能边缘态。
  • Rice-Mele (RM) 模型：在 SSH 模型的基础上引入了交替的格点势 $\pm \Delta_d$。RM 模型由于破缺了反演对称性，其边缘态不再被锁定在绝对零能，但在能隙中依然存在强局域化的边缘态，并伴随有有限的子格点极化。

1.3 技术难点与现有方法的局限性

在一维多体系统的研究中，常用的数值方法各有优缺点，但在处理有限温度下的非局域相互作用时遇到了巨大阻碍：

  • **密度矩阵重整化群（DMRG）和基态量子蒙特卡洛（QMC）**极其擅长零温（$T=0$）下的计算，但在高精度处理有限温动力学、自能结构以及大尺度谱函数时计算成本极高。
  • **动力学平均场理论（DMFT）**通过将格点模型映射到单杂质 Anderson 模型，能够精确捕获局域电子关联，但单格点 DMFT 的自能 $\Sigma(\omega)$ 是空间均匀且无动量依赖的。由于 BOW 相的序参量 $O_{BOW} = \frac{1}{2} \text{Re}[\Sigma_{intra}(\nu_0) - \Sigma_{inter}(\nu_0)]$ 本质上是非局域的，单格点 DMFT 完全无法描述此类物相。
  • 传统的团簇 DMFT（如 CDMFT）虽然引入了短程空间关联，但在有限温度下，如何在图符微扰论中一致地将非局域库仑相互作用 $V$ 引入到自能的重整化中，一直是个悬而未决的难题，通常会遭遇非因果性（noncausality）问题，即自能虚部在某些频率下变为正值（$\text{Im}\Sigma(i\nu_0) > 0$）。

1.4 方法细节：D-TRILEX 团簇扩展与局部场自洽（Local-Field Surrogate）

为了克服这些技术瓶颈，作者采用并扩展了最近开发的 D-TRILEX 团簇框架（参见论文引文 [75]）。其基本思想是以二格点（dimer）DMFT 为参考系统，将二聚体内部的局域和最近邻关联进行精确的非微扰处理，而将二聚体之外的更长程关联通过基于双玻色子/双费米子图符微扰论的 D-TRILEX 展开进行处理。具体步骤如下：

  1. 参考系统的构建与求解： 参考系统选择为在 $V=0$ 处求解的二格点 DMFT 杂质模型。这保证了二聚体内最强烈的局域 $U$ 关联被精确求解，获得了非微扰的格林函数 $g_{\nu}$ 和杂质杂化函数 $\Delta_{\nu}$。

  2. 外自洽（Outer Self-Consistency）与静态对称性破缺场 $h$： 由于 BOW 和 CDW 都是自发对称性破缺相，在参考系统中必须引入一个静态的对称性破缺场 $h$ 来诱导有序态。在二聚体空间中，$h$ 表示为一个 2x2 的厄米矩阵：

    $$h = h_x \sigma^x + h_z \sigma^z$$

    其中，对角成分 $h_z$ 充当与 CDW 共轭的交替格点势，而对角外成分 $h_x$ 充当与 BOW 共轭的键二聚化场。外自洽条件要求参考系统的密度矩阵与晶格系统的密度矩阵在二聚体内完全一致：

    $$\langle n^{ref}_{l_1 l_2}[h] \rangle \stackrel{!}{=} \langle n^{latt}_{l_1 l_2}[h] \rangle$$

  3. 线性响应近似与图符展开： 由于在每个外自洽循环中重新求解带非零 $h$ 的强关联二格点杂质模型极其昂贵，作者引入了线性响应近似（Linear Response Approximation）。自能对外加静态场 $h$ 的偏导数可以通过参考系统中的电荷通道三点顶点函数 $\lambda^{ch}_{\nu,\omega=0}$ 来计算：

    $$\partial \Sigma_{\nu, l_1 l_2} = \sum_{l_3, l_4} \bar{\lambda}^{ch}_{\nu,\omega=0, l_1 l_2, l_3 l_4} h_{l_3 l_4}$$

    其中，经过扣除自由场贡献后的顶点函数 $\bar{\lambda}^{ch}$ 完整捕获了多体相互作用对响应的重整化。经过场修饰后的二聚体格林函数 $g^h_{\nu}$ 为：

    $$g^h_{\nu} = \left[ g_{\nu}^{-1} - \lambda^{ch}_{\nu,\omega=0} h \right]^{-1}$$

    这一局部场替代方案极大地降低了计算复杂度，同时通过线性响应理论严格保证了因果性。

  4. 低能有效哈密顿量（Topological Hamiltonian Approach）： 为了在强关联下定义拓扑不变量，作者利用了拓扑哈密顿量理论。在低能极限下，由于动力学涨落被有序相的能隙所抑制，自能的虚部在零频极限下趋于零（$\text{Im}\Sigma(k, \nu \to 0) \to 0$）。因此，可以定义一个等效的非相互作用哈密顿量：

    $$H_{eff}(k) = \epsilon_k + \text{Re} \Sigma(k, \nu_0)$$

    其中 $\nu_0 = \pi / \beta$ 为最低费米子马苏巴拉（Matsubara）频率。$H_{eff}(k)$ 与全格林函数 $G(k, 0)$ 具有完全相同的拓扑结构（如零点和极点分布），可以直接用于计算拓扑不变量（如 Zak 相位）以及在开放边界条件下对边缘态进行精确对角化。

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

2.1 物理参数设置

为了展示电子关联如何诱导拓扑物相，作者对半满一维扩展哈伯德模型进行了详尽的数值模拟，设定温度为 $T=0.1$（该温度足够低，能够清晰展现强关联相变的特征，同时避免了绝对零温下可能遇到的数值发散）。主要探索的参数空间为：

  • 局域排斥强度：$U \in [0, 3.5]$
  • 最近邻排斥强度：$V \in [0, 3.0]$
  • 跃迁强度：$t = 1$ 作为能量单位。

2.2 相图分析（Figure 1）与关键数据点

通过执行 D-TRILEX 团簇自洽计算，作者绘制了极其精确的 $U-V$ 相图（见主文 Figure 1）。相图主要由三个相组成：

  1. 无序相（Unordered Phase）：位于低 $V$ 区域。随着 $U$ 的增大，无序相的范围向高 $V$ 方向延伸。这是由于增强的局域单排斥抑制了电荷涨落，使得系统需要更强的最近邻相互作用 $V$ 才能稳定电荷或键的有序排列。
  2. 键序波拓扑相（BOW Phase）：在弱到中等 $U$ 强度（例如 $U=1.5$）下，随着 $V$ 增加，系统首先跨越临界值 $V^{BOW}_c = 0.99$ 进入 BOW 相。在这个相内，由于自发的键二聚化，格点间的等效跃迁系数交替变化，单粒子自能显现出显著的非局域分量，使得等效哈密顿量映射为经典的拓扑 SSH 模型。
  3. 电荷密度波拓扑相（CDW+BOW Phase）：当 $V$ 进一步增大跨越 $V^{CDW}_c$（例如在 $U=1.5$ 时，$V^{CDW}_c = 1.07$；在 $U=3.0$ 时，$V^{CDW}_c = 1.95$）后，系统进入 CDW 相。本工作的一个重大发现是：该 CDW 相绝非纯粹的电荷有序，而是伴随着强烈的领头或次领头键序波（BOW）关联。因此，它实际上是一个 CDW+BOW 混合拓扑相。在这个相中，格点势能的交替（由 $h_z$ 驱动）与跃迁系数的交替（由 $h_x$ 驱动）同时存在，使得等效哈密顿量在低能下完美映射为 Rice-Mele 模型。

2.3 序参量的演化与相变特征（Figure 2）

Figure 2 展示了在不同 $U$ 条件下，CDW 序参量 $O_{CDW}$ 和 BOW 序参量 $O_{BOW}$ 随 $V$ 的演化规律：

  • 对于 $U = 1.5$（中等耦合极限）：

    • 在 $V < 0.99$ 处，系统处于无序相，但由于二聚体参考系统的几何截断效应，存在一个微小的准线性背景 $O_{BOW}$。作者通过定义序参量曲率 $\kappa(V) = \partial^2 O_{BOW} / \partial V^2$ 的极大值，确定了真实的物理相变点 $V^{BOW}_c = 0.99$。
    • 在 $V^{BOW}_c < V < V^{CDW}_c = 1.07$ 范围内，$O_{CDW} = 0$，而 $O_{BOW}$ 显现出显著的非线性自发增长，这是纯粹由关联驱动的拓扑 SSH 相。
    • 一旦 $V > 1.07$，$O_{CDW}$ 陡然跳跃到一个有限值（展现出一阶相变特征），与此同时，$O_{BOW}$ 不仅没有受到抑制，反而展现出一个断裂式的显著增强（不连续的跳跃），直接证明了在 CDW 区域内，空间反演对称性和格点平移对称性同时被强烈破缺，从而证实了 CDW+BOW 混合相的存在。

  • 对于 $U = 3.0$（强耦合极限）：

    • 由于 $U$ 强烈抑制了局域电荷涨落，纯 BOW 相被完全挤压消失。系统直接在 $V^{CDW}_c = 1.95$ 处从无序相直接转变到 CDW+BOW 拓扑相，表现为 $O_{CDW}$ 伴随着 $O_{BOW}$ 的突变同时涌现。

2.4 单核与团簇 D-TRILEX 结果对比（方法性能展示）

在 Figure 1 中，作者对比了单格点 D-TRILEX（红虚线）和二格点团簇 D-TRILEX（红实线）的 CDW 相界：

  • 在弱耦合区（$U \le 2.0$），单格点与二格点方法的 CDW 相界几乎重合。这表明在该区域内，非局域关联是弱微扰的，可以通过图符展开得到良好描述。
  • 在强耦合区（$U \ge 3.0$），单格点计算显著高估了 CDW 的稳定性，其相界维持准线性增长。而二格点团簇相界则显著向上弯曲，即需要极大的最近邻排斥 $V$ 才能诱导 CDW。这表明此时二聚体内部的非微扰短程自旋单态（singlet）涨落开始主导，这种强烈的磁不稳定性与电荷有序产生竞争，抑制了 CDW。单格点方法由于缺乏二聚体内非微扰关联的处理，完全忽略了这一竞争效应，从而凸显了团簇 D-TRILEX 方案的优越性与必要性。

2.5 拓扑边缘态谱分析（Figure 3）

为了提供决定性的拓扑证据，作者利用自洽求解出的有效哈密顿量 $H_{eff}(k)$，在具有开放边界条件（OBC）的 $L=64$ 格点长链上进行了精确对角化（Figure 3）。

  • 在 $V < V^{BOW}_c$ 区域：所有能态均属于体态（bulk states），呈现为均匀分布在整个格点上的波函数，没有在边缘局域化的态。边缘局域化权重 $|\phi^{(n)}_{edge}|^2 \approx 0$，子格点极化 $|P^{(n)}| \approx 0$。
  • 在 BOW 相区（$V^{BOW}_c < V < V^{CDW}_c$）：有两个能态迅速从体能带中脱离，向绝对零能移动。在 BOW 相内，这两个态的能量被精确锁定在 $E=0$ 处。这两个零能模在物理边界上表现出极强的指数衰减局域化特征（Figure 3 中的红色/黑色标记，局部化权重极高，但子格点极化 $|P^{(n)}| \approx 0$），这正是经典 SSH 拓扑非平庸相中受对称性保护的零能边缘态的完美再现。
  • 在 CDW+BOW 相区（$V > V^{CDW}_c$）：随着格点对称性破缺，格点势能交替开启，这两个边缘态被推离零能，向体能隙的两侧漂移。然而，它们依然维持在体能隙之中，并且依然保持着极强的边缘局域化特征（表现为深色/黑色三角形）。与此同时，由于交替格点势的作用，这两个边缘态自发获得了强烈的子格点极化（$|P^{(n)}| \to 1$，三角形方向指向左或右，代表电荷强烈局域在 A 或 B 子格点上）。这与 Rice-Mele 模型中的拓扑边缘电荷极化特征完全吻合，从而在数值上无可辩驳地证实了该相的拓扑本质。

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

要复现本论文的核心结果，研究人员需要一个集成有高性能量子杂质求解器（如连续时间杂化膨胀算法 CT-HYB）的团簇 DMFT 框架，并在此基础上编写 D-TRILEX 图符校正和局部场外自洽循环程序。

3.1 核心算法复现流程图

 ┌────────────────────────────────────────────────────────┐
 │ 1. 求解 V=0 处的 2-site (Dimer) DMFT 杂质模型         │
 │    - 得到格林函数 g_v, 杂化函数 Delta_v                │
 │    - 计算非微扰的三点顶点函数 \lambda^ch_{v, \omega=0}   │
 └───────────────────────────┬────────────────────────────┘
                             │
                             ▼
 ┌────────────────────────────────────────────────────────┐
 │ 2. 初始化对称性破缺场 h = h_x * \sigma^x + h_z * \sigma^z│
 └───────────────────────────┬────────────────────────────┘
                             │
                             ▼
                       ┌───────────┐
                ┌─────>│ 迭代开始  │
                │      └─────┬─────┘
                │            │
                │            ▼
                │┌────────────────────────────────────────┐
                ││ 3. 根据场修饰格林函数:                   │
                ││    g^h_v = [g_v^-1 - \lambda^ch * h]^-1 │
                │└───────────────────┬────────────────────┘
                │                    │
                │                    ▼
                │┌────────────────────────────────────────┐
                ││ 4. 运行 D-TRILEX 晶格图符求解器:        │
                ││    - 计算动量依赖的自能 \Sigma(k, iw_n) │
                ││    - 获得格林函数 G(k, iw_n)            │
                ││    - 计算晶格电荷/键密度分布 n_latt       │
                │└───────────────────┬────────────────────┘
                │                    │
                │                    ▼
                │┌────────────────────────────────────────┐
                ││ 5. 求解局部自洽方程以更新场 h:         │
                ││    <n_ref[h]> = <n_latt[h]>            │
                │└───────────────────┬────────────────────┘
                │                    │
                │          [未收敛]  │
                └────────────────────┴────────────────────
                                     │ [收敛]
                                     ▼
 ┌────────────────────────────────────────────────────────┐
 │ 6. 提取低能有效哈密顿量:                              │
 │    H_eff(k) = \epsilon_k + Re \Sigma(k, iv_0)            │
 └───────────────────────────┬────────────────────────────┘
                             │
                             ▼
 ┌────────────────────────────────────────────────────────┐
 │ 7. 在 L=64 的 OBC 长链上对角化 H_eff(k)               │
 │    - 绘制能谱演化, 计算局域化权重和子格点极化          │
 └────────────────────────────────────────────────────────┘

3.2 关键数学物理方程的代码化实现指南

A. 场修饰格林函数的矩阵构建

在马苏巴拉频率轴上，对每个 $\nu_n$，二聚体格林函数 $g_{\nu_n}$ 和顶点函数 $\lambda^{ch}_{\nu_n, 0}$ 均为 2x2 矩阵。场修饰方程（方程 S3）需要采用精确的矩阵逆运算。在 Python 中可利用 numpy.linalg.inv 实现：

import numpy as np

def get_field_dressed_green(g_v, lambda_ch, h_x, h_z):
    """
    g_v: numpy array of shape (N_freq, 2, 2) - Dimer green function
    lambda_ch: numpy array of shape (N_freq, 2, 2, 2, 2) - Three-point vertex
    h_x, h_z: float - static fields
    """
    N_freq = g_v.shape[0]
    g_h_v = np.zeros_like(g_v, dtype=complex)
    
    # 构建场矩阵 h_matrix = h_x * sigma_x + h_z * sigma_z
    sigma_x = np.array([[0, 1], [1, 0]], dtype=complex)
    sigma_z = np.array([[1, 0], [0, -1]], dtype=complex)
    h_matrix = h_x * sigma_x + h_z * sigma_z
    
    for n in range(N_freq):
        g_inv = np.linalg.inv(g_v[n])
        # 收缩顶点函数 lambda_ch[n] 与 h_matrix
        # \lambda_corr_{l1, l2} = \sum_{l3, l4} \lambda_{l1, l2, l3, l4} * h_{l3, l4}
        lambda_corr = np.einsum('ijkl,kl->ij', lambda_ch[n], h_matrix)
        
        g_h_inv = g_inv - lambda_corr
        g_h_v[n] = np.linalg.inv(g_h_v_inv)
        
    return g_h_v

B. 求解自洽方程（方程 S4）

自洽条件寻找合适的 $h_x, h_z$ 使得二聚体参考密度的对角和外对角成分与通过 D-TRILEX 算出的晶格密度一致。可以采用 scipy.optimize.root 中的多维非线性方程组求解器（如 Powell 或 Levenberg-Marquardt 算法）：

from scipy.optimize import root

def self_consistency_objective(h_vector, dimer_solver, trilex_solver):
    h_x, h_z = h_vector
    
    # 1. 计算当前场 h 下的参考系统密度
    n_ref = dimer_solver.compute_density(h_x, h_z)
    
    # 2. 运行 D-TRILEX 获得当前场下的晶格密度
    n_latt = trilex_solver.solve(h_x, h_z)
    
    # 3. 目标：差异趋于 0
    diff_x = n_ref.off_diagonal - n_latt.off_diagonal # 用于确定 h_x (BOW)
    diff_z = n_ref.staggered - n_latt.staggered       # 用于确定 h_z (CDW)
    
    return [diff_x, diff_z]

# 执行自洽求解
sol = root(self_consistency_objective, x0=[0.1, 0.0], method='hybr')
h_x_opt, h_z_opt = sol.x

3.3 推荐开源软件包及 Repo 链接

  1. TRIQS (Toolbox for Research on Interacting Quantum Systems): 强关联多体物理计算的黄金标准框架，包含高精度的 C++ 库和 Python 接口，内置了基础的 DMFT 和 CT-HYB 求解器。
  2. D-TRILEX / Dual-Fermion 相关开源实现: 由于 D-TRILEX 是近期开发的高级微扰论方法，完整的代码库通常由各研究组维护。研究人员可以参考典型的双费米子（Dual Fermion）或图符自洽计算代码进行二次开发：

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

4.1 关键历史引用文献

  1. Su-Schrieffer-Heeger (SSH) 模型原著:
    • Su, W. P., Schrieffer, J. R., & Heeger, A. J. (1979). Solitons in polyacetylene. Physical Review Letters, 42(16), 1698. (拓扑一维模型的奠基之作，定义了非平庸 Zak 相和零能边缘态)。
  2. Rice-Mele (RM) 模型原著:
    • Rice, M. J., & Mele, E. J. (1982). Elementary excitations of a linearly conjugated diatomic polymer. Physical Review Letters, 49(19), 1455. (引入交替势的二聚化长链模型，描述了拓扑电荷极化泵浦)。
  3. 拓扑哈密顿量方法奠基:
    • Wang, Z., & Zhang, S. C. (2012). Simplified topological invariants for interacting insulators. Physical Review X, 2(3), 031008. (证明了低能极限下格林函数的拓扑性质可以用零频有效自能构建的 $H_{eff}$ 精确描述)。
  4. D-TRILEX 核心方法文献:
    • Fossati, F., & Stepanov, E. A. (2026). Dual-space cluster-diagrammatic approach to nonlocal electronic correlations. Physical Review B, 113, 075149. (系统阐述了二聚体 D-TRILEX 微扰展开的数学框架)。

4.2 本工作局限性剖析与客观评论

虽然该项工作在理论上取得了突破性进展，深刻揭示了关联涌现拓扑的物理图像，但从高精度量子化学和凝聚态物理的严苛视角来看，仍存在以下局限性：

  1. 团簇尺度的截断效应（Dimer Approximation）： 该方法的微扰参考系统仅为一个二格点（dimer）系统。虽然它相比单格点 DMFT 迈出了关键一步（能够处理最近邻关联），但在严格的一维系统中，长程量子涨落和电荷/自旋关联（例如 $1/r$ 衰减的幂律关联，即 Luttinger 液体行为）具有全局性。仅用 2-site 团簇作为起点，可能在极低温度和临界点附近低估更长程涨落对拓扑相的破坏作用。未来有必要将参考系统扩展到 4-site 或 8-site 团簇，以验证相界的系统收敛性。

  2. 有限温度下一维连续与离散对称性破缺的物理限制： 根据 Mermin-Wagner 定理，在有限温度下一维系统无法自发破缺任何连续对称性。本研究涉及的 CDW（电荷交替）和 BOW（键二聚化）属于离散的 $Z_2$ 对称性破缺。虽然离散对称性在一维零温下可以完全破缺，但在有限温度下，严格的长程有序依然会受到热涨落激发的孤子（soliton）和反孤子的强烈热激发干扰，从而退化为准长程有序或有限关联长度。论文中观察到的清晰拓扑边缘态，很大程度上得益于计算是在 $T=0.1$ 的有限尺寸格点链上进行的，且采用了有效平均场性质的静态 $h$ 场自洽。若在热力学极限下完全考虑动态涨落，边缘态的寿命（相干时间）和退相干效应需要更精细的动力学分析。

  3. 马苏巴拉零频近似的系统误差： 将 $\text{Re}\Sigma(k, \omega=0)$ 近似为最低马苏巴拉频率下的值 $\text{Re}\Sigma(k, i\nu_0)$ 是一种极其便利的数值权宜之计，避免了病态的解析延拓。然而，在强关联区，自能在实频低能区可能存在复杂的非平庸动力学结构（例如，由于低能自旋激发导致的非费米液体行为）。用最低马苏巴拉频率代替零频，可能会掩盖真正的低能激发动力学，从而影响对有效哈密顿量 $H_{eff}$ 精度边界的严谨评估。

5. 补充探讨：物理机理深挖与前沿展望

为了帮助量子化学和固体物理背景的读者更透彻地理解该工作的物理内核，以下对两个深层物理机制进行补充探讨。

5.1 从非局域库仑排斥 $V$ 到跃迁二聚化的物理机制

在非相互作用紧束缚模型中，要产生二聚化，必须人工设定不均匀的键长（例如，令奇数键跃迁为 $t - \delta$，偶数键跃迁为 $t + \delta$）。那么，在均匀的长链上，最近邻库仑排斥 $V$ 是如何自发催生出交替的跃迁重整化的？

我们可以从一阶微扰论的 Hartree-Fock 交换项（Fock 关联）来直观理解这一物理过程。在关联系统里，非局域相互作用项可以做如下解耦：

$$V n_i n_j \approx V \langle n_i \rangle n_j + V n_i \langle n_j \rangle - V \chi_{ij} \left( c^{\dagger}_i c_j + c^{\dagger}_j c_i \right)$$

其中，$\chi_{ij} = \langle c^{\dagger}_i c_j \rangle$ 为格点 $i$ 和 $j$ 之间的键关联函数（即键电荷密度）。

  • Hartree 项 ($V \langle n_i \rangle n_j$) 驱动了 CDW 有序：如果奇数格点电荷密度高，偶数格点电荷密度低，Hartree 项就会产生一个交替的势能，使这一电荷不均匀性自洽地维持下去。

  • Fock 项 ($- V \chi_{ij} c^{\dagger}_i c_j$) 直接重整化了等效跃迁系数：

    $$t_{eff, ij} = t + V \chi_{ij}$$

    这意味着，如果由于量子涨落，某一个键上的电子重叠积分 $\chi_{ij}$ 稍大，非局域排斥 $V$ 就会通过 Fock 交换项使该键的等效跃迁 $t_{eff}$ 变得更大；相反，邻近的键等效跃迁会变小。这种正反馈机制在 $V$ 足够大时，会自发导致 $\chi_{ij}$ 在奇数和偶数键上交替，从而在无需任何晶格畸变（Peierls 不稳定性）的情况下，纯粹通过非局域电子关联自发产生跃迁的二聚化（BOW 相）。本研究正是利用 D-TRILEX 方法，在超越 Hartree-Fock 的非微扰水平上，精确证明了这一涌现拓扑相的存在。

5.2 实验实现前景与冷原子平台

在真实的固体材料（如过渡金属氧化物或一维有机导体 Bechgaard 盐）中，电子-声子耦合往往与电子-电子相互作用交织在一起，很难说拓扑相是由电子关联主导还是由传统的晶格 Peierls 畸变主导。

然而，**光晶格中的冷原子气体（Cold Atoms in Optical Lattices）**为复现和验证本工作提供了一个近乎完美的物理平台。利用超冷极性分子（如 $^{23}\text{Na}^{87}\text{Rb}$）或具有大磁偶极矩的强磁性原子（如 $^{162}\text{Dy}$ 或 $^{166}\text{Er}$），实验物理学家可以在光晶格中实现具有长程偶极排斥作用的一维扩展哈伯德模型。通过调节外加磁场（利用 Feshbach 共振改变局域 $U$）和调整偶极子的取向（调节最近邻排斥 $V$），可以精确控制相互作用参数。由于光晶格中完全不存在固体中的晶格缺陷和声子干扰，超冷原子实验可以纯粹在电子（原子）关联水平上复现本工作预言的自发 BOW 拓扑相和 CDW+BOW 混合拓扑相，通过飞行时间原位成像（in-situ imaging）和布拉格光谱技术直接观测到强局域化的边界原子态。

5.3 展望：向更高维度的推广

本工作建立的“关联自发诱导拓扑”的理论框架（D-TRILEX + 局部自洽场 + 拓扑哈密顿量）具有极强的普适性。一个极具前景的方向是将该方法推广至二维格点（如蜂窝状 Honeycomb 晶格或 Kagome 晶格）。在二维情况下，非局域相互作用 $V$ 有望自发诱导产生时空反演对称性破缺的自发环流（Loop Currents），从而在零磁场下自发产生非零的 Chern 数，诱导涌现出拓扑莫特绝缘体（Topological Mott Insulators）或分数量子反常霍尔相。这不仅能深化对强关联多体物理的理解，也为设计新型的无耗散拓扑电子学器件提供了全新的材料设计逻辑。