来源论文: https://arxiv.org/abs/2603.10238v1 生成时间: Mar 12, 2026 03:15

执行摘要

在强关联电子系统研究中，Mott 绝缘体与无序导致的局域化（如 Anderson 局域化）是两个核心范式。传统观点往往将两者分开讨论，但当两者并存时，系统展现出的物理特性极为复杂。本文深度解析了一篇关于“Mott-Hubbard 模型中无序诱导局域化”的科研工作。该研究利用“关联层次结构（Hierarchy of Correlations）”理论框架，系统探讨了在二维点阵（六角、正方形、蜂窝状）中，电荷无序（随机原位势）和自旋无序（随机静态自旋背景）对准粒子激发（空穴子 holon 和双电子 doublon）演化的影响。

研究的核心发现包括：

  1. 电荷无序会导致能谱分裂，并在空间和能量上产生局域态与离域态的清晰分离，表现为明显的双峰参与率（Participation Ratio）分布。
  2. 自旋无序则在整个能带内产生局域态，且不引起明显的能带分裂。
  3. 1/Z 关联层次理论在捕捉准粒子细微结构方面优于传统的一阶强耦合微扰理论。

1. 核心科学问题，理论基础与方法细节

1.1 核心科学问题：Mott 与 Anderson 的博弈

Mott 绝缘体是由电子间的强库仑排斥作用（$U$）引起的，而 Anderson 局域化则是由单体势的无序分布（$V$）引起的。在强关联极限下（$T \ll U$），系统的基本激发不再是自由电子，而是所谓的“准粒子”：holon（空穴子，带正电荷，无自旋）和 doublon（双电子，带负电荷，无自旋）。研究的问题在于：这些由强关联产生的准粒子，在面对环境无序时，其定域化行为遵循什么规律？多体局域化（MBL）在二维系统中是否稳定？

1.2 理论基础：Fermi-Hubbard 模型与无序项

研究的出发点是二维 Fermi-Hubbard 哈密顿量（公式 1）：

$$\hat{H} = -\frac{1}{Z} \sum_{\mu,\nu,s} T_{\mu\nu}\hat{c}^\dagger_{\mu,s}\hat{c}_{\nu,s} + U \sum_{\mu} \hat{n}_{\mu,\uparrow}\hat{n}_{\mu,\downarrow}$$

其中 $Z$ 是配位数。无序通过两种方式引入：

  • 电荷无序（Charge Disorder）：引入原位势 $\hat{H}_c = \sum_{\mu} V_\mu (\hat{n}_{\mu,\uparrow} + \hat{n}_{\mu,\downarrow})$，其中 $V_\mu$ 以概率 $r$ 取 $V$，否则为 0。这种二值分布与传统的连续 Anderson 无序不同。
  • 自旋无序（Spin Disorder）：引入自旋分裂项 $\hat{H}_s = \sum_{\mu} v_\mu (\hat{n}_{\mu,\uparrow} - \hat{n}_{\mu,\downarrow})$，其中 $v_\mu \in \{-v, v\}$。在时间尺度 $\sim 1/T$ 内，自旋背景被假定为静态，这在实验中对应于自旋动力学被冻结的状态。

1.3 技术难点：强关联下的动力学求解

直接处理强关联无序系统的多体动力学在计算上是指数级爆炸的。传统的动力学平均场理论（DMFT）在处理空间非均匀性时成本极高。因此，本文采用了 1/Z 关联层次理论（Hierarchy of Correlations）。这种方法通过 $1/Z$ 的幂次展开来截断相关的密度矩阵方程。在 $Z \gg 1$ 的极限下，系统倾向于平均场行为，而 $1/Z$ 项捕获了邻近格点间的关联效应。

1.4 方法细节：准粒子算符与方程闭合

为了描述准粒子，定义了“着装（Dressed）”费米子算符（公式 6）：

  • $I=1$ (Doublon): $\hat{C}^1_{\mu,s} = \hat{n}_{\mu,\bar{s}} \hat{c}_{\mu,s}$
  • $I=0$ (Holon): $\hat{C}^0_{\mu,s} = (1 - \hat{n}_{\mu,\bar{s}}) \hat{c}_{\mu,s}$

通过对二点关联函数进行因式分解（Factorisation ansatz），研究者推导出了准粒子概率幅 $p^I_{\mu,s}$ 的演化方程（公式 8）：

$$i\partial_t p^I_{\mu,s} = U^I p^I_{\mu,s} - \frac{1}{Z} \sum_{\nu,X} T_{\nu\mu} \langle \hat{N}^I_{\mu,\bar{s}} \rangle^0 p^X_{\nu,s}$$

该方程组的优势在于它将复杂的多体问题转化为了具有有效势能项和有效跳迁项的线性方程组，从而允许在包含数千个格点的系统上进行数值求解。

2. 关键 Benchmark 体系与数据分析

2.1 模拟体系设置

研究对比了三种具有不同对称性和配位数的二维格点：

  1. 六角格点（Hexagonal, Z=6）：提供最强的平均场近似效果。
  2. 正方形格点（Square, Z=4）：最常见的凝聚态模拟平台。
  3. 蜂窝状格点（Honeycomb, Z=3）：量子波动效应最强，接近理论下限。 计算设定 $U = 0.35$ eV，$T/Z = 0.05$ eV，格点规模通常为 $25 \times 25 = 625$ 个格点。

2.2 关键数据：态密度（DOS）分裂

在电荷无序情况下（图 1），系统表现出显著的子带分裂。当无序势 $V = -0.20$ eV 时，原始的 holon 带和 doublon 带各自裂解为两部分。这两部分的能量间距恰好由 $V$ 决定。这意味着准粒子能够“感知”原位势的阶梯，从而在能量空间形成隔离区。

2.3 参与率（Participation Ratio, PR）分析

这是区分局域化与离域态的关键性能数据。研究者计算了逆参与率（IPR），并提取了 PR（即 $1/IPR$）。

  • 电荷无序（图 4 & 6）：PR 分布呈现明显的“双峰结构”。第一个峰值靠近 0，对应于局域在少数格点上的状态（由掺杂引起的团簇）；第二个峰值对应于在整个格点上延伸的准平面波状态。这种二元性表明系统在空间上发生了局域与离域的共存。
  • 自旋无序（图 7 & 8）：与电荷无序不同，自旋无序的 PR 分布是连续且单峰的，所有状态都表现出一定程度的局域化倾向。这意味着自旋随机性造成的相干性破坏在全能带范围内是均匀的。

2.4 格点几何效应

对比图 6 可以发现，随着 $Z$ 的降低（从 6 到 3），局域态的比例显著增加。在蜂窝状格点中，由于连通性差，准粒子更容易被无序势捕获，这符合渗流理论（Percolation Theory）的预测：配位数越低，形成无限扩展团簇的临界掺杂率 $r$ 越高。

3. 代码实现细节与复现指南

3.1 核心算法：线性代数与对角化

复现该工作的核心在于求解方程 (10) 和 (14) 的特征值问题。这是一个典型的非厄米（在某些近似下可转化为厄米）矩阵对角化问题。

计算步骤：

  1. 格点生成：构建 $L \times L$ 的拓扑邻接矩阵 $T_{\mu\nu}$。需要支持周期性边界条件（PBC）。
  2. 无序生成：按概率 $r$ 生成二值随机向量 $V_\mu$ 或 $v_\mu$。
  3. 哈密顿矩阵构造：
    • 矩阵维度为 $2N^2 \times 2N^2$（考虑 holon 和 doublon 的耦合）。
    • 对于电荷无序，哈密顿矩阵元素由方程 (10) 给定：对角项为 $U^I + V_\mu$，非对角项为 $-T_{\nu\mu}/2Z$。
  4. 求解器选择：对于 $N=25$ 的体系，矩阵大小约为 $1250 \times 1250$，标准 LAPACK（如 numpy.linalg.eigh）即可处理。若要扩展到 $N > 100$，需使用 ARPACK 或 Krylov 子空间法。
  5. 统计平均：由于无序的随机性，需要对 100-500 个样本进行蒙特卡洛平均，以获得平滑的 DOS 曲线。

3.2 软件包建议

  • 编程语言：推荐使用 Julia（其线性代数性能与 C 相当，且语法接近 Python）或 Python (NumPy + SciPy)。
  • 库支持：
    • LinearAlgebra.jl (Julia)
    • scipy.sparse.linalg (Python)
    • QuantPylib (假设存在此类量子模拟库)
  • 可视化：Matplotlib 或 Plotly 用于绘制 DOS 和格点局域化状态图（如图 3）。

3.3 关键参数复现

  • $U = 0.35$
  • $T = 0.05 \times Z$
  • $r = 0.26$（这是一个关键点，接近某些格点的渗流阈值）
  • 时间尺度限制：$t < U/T^2$（确保自旋静态假设成立）。

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

4.1 关键引用

  1. Anderson (1958) [1]：局域化理论的奠基石。
  2. Hubbard (1963) [15]：定义了强关联电子系统的基本框架。
  3. Basko et al. (2006) [11]：关于多体局域化（MBL）的现代阐述。
  4. Navez & Schützhold (2010/2014) [16, 17]：本文所用 $1/Z$ 关联层次理论的数学来源。

4.2 工作评论与局限性

优势：

  • 该工作成功避开了昂贵的 DMFT 计算，利用 $1/Z$ 展开提供了一个直观的准粒子动力学图像。
  • 明确区分了电荷无序和自旋无序对局域化的不同物理机制（能带分裂 vs 全带局域化）。

局限性（技术角度）：

  1. 静态背景假设：研究假设自旋背景是静态的。在真实的量子系统中，即使 $T \ll U$，超交换作用（Superexchange, $J \sim T^2/U$）依然会导致自旋翻转，这会破坏本研究描述的相干局域化。因此，该结论仅在极短的时间尺度内有效。
  2. 1/Z 截断误差：仅保留到一阶。在低配位数系统（如蜂窝格点 $Z=3$）中，高阶项（如环状跃迁）可能非常重要，论文对此的讨论略显不足。
  3. 有限尺寸效应：$N=25$ 在二维系统中虽然不算小，但对于接近相变点的局域化研究，长程关联可能受到边界条件的限制。图 4 中不同尺寸的 PR 分布虽有重合，但趋势外推至无限大体系时仍需谨慎。

5. 补充：深度原理解析与物理直觉

5.1 为什么电荷无序会导致子带分裂？

从物理直觉上看，电荷无序引入的是一个“二值”势。这意味着一部分格点能量为 0，另一部分格点能量为 $V$。对于准粒子来说，这相当于在格点阵列中挖出了大量的“坑”或“障。如果准粒子的动能（由 $T/Z$ 决定）小于势垒 $V$，它们就会被困在特定的掺杂团簇内。这种能量上的显著不连续性直接导致了 DOS 的断裂。

5.2 自旋无序的特殊性

自旋无序并不改变格点的原位势能总量，而是改变了自旋向上和向下电子的相对能量。因为准粒子（holon/doublon）在移动时必须保持或交换背景自旋，随机的自旋排列就像是一个充满随机相位因子的介质。这种随机性对所有动量状态的准粒子都是“平等的”，因此不会产生能隙，而是导致整个能带的弥散化和局域化。

5.3 $1/Z$ 展开与平均场理论的关系

读者可以将 $1/Z$ 关联层次结构理解为“超越平均场”的第一步。在 $Z \to \infty$ 时，每个格点感受到的是周围格点的平均势场（Weiss Field）。当 $Z$ 为有限值时，格点 $\mu$ 与格点 $\nu$ 之间的瞬时关联不可忽略。本文的方法通过保留两点关联函数 $\rho^{corr}_{\mu\nu}$，实际上捕获了准粒子在跳迁过程中的受激激发现象。这解释了为什么在图 9 和图 10 的对比中，关联层次理论能看到微扰论看不见的“细微结构”。

5.4 对量子模拟实验的指导意义

对于冷原子实验（Optical Lattice），该研究建议，通过精确控制人工规范势或随机光晶格，可以分别模拟这两类无序。如果实验观测到的 DOS 出现了明显的 sub-band 分裂，则可以判定为电荷主导的无序；若能谱保持连续但输运特性消失，则更倾向于自旋无序引起的局域化。

总结：这项工作为我们理解强关联物质中的局域化提供了宝贵的理论工具，特别是证明了在处理非均匀无序体系时，基于密度矩阵的关联层次展开比传统的费曼图微扰论更具表现力。对于计算化学家而言，这种将复杂关联因式分解的思路，也为开发更高效的多体动力学算法提供了灵感。