Local fermion density in inhomogeneous free-fermion chains: a discrete WKB approach

本文提出了一种基于离散 WKB 近似的解析方法，通过直接求解单粒子本征函数的递推关系，推导出了适用于任意填充、跃迁强度和磁场分布的均匀自由费米子链局域费米子密度闭式表达式，从而为理解非均匀系统中的纠缠熵抑制现象提供了超越传统场论技术的理论框架。

要点

这篇论文讲述了一个关于**“如何在拥挤且地形复杂的房间里安排客人”**的数学故事。

想象一下，你有一排很长的座位（这就叫“一维量子链”），上面坐满了客人（这就是“费米子”，一种微观粒子）。在普通的房间里，座位间距和灯光都是一样的，客人分布得很均匀。但在这篇论文研究的模型里，座位的间距忽宽忽窄（跳跃强度 $J_n$ 变化），而且头顶的灯光忽明忽暗（磁场 $B_n$ 变化）。

科学家们想知道：在这种复杂的环境下，如果你只让一定数量的客人坐下（比如只让前 30% 的座位有人），那么每个座位上坐人的概率（局部密度）到底是多少？

1. 核心难题：为什么以前很难算？

以前，科学家想算这个问题，就像试图用**“天气预报”的模型去预测“微观粒子”**的行为。

• 传统方法（场论）：就像用气象卫星看大范围的云层。这种方法在“天气平稳”（均匀系统）或“临界点”（半满状态）时很准。
• 问题所在：一旦系统变得不均匀（座位忽宽忽窄），或者客人填得很少/很满，传统的“气象卫星”就失灵了。之前的研究只能猜出“低密度”时的情况，却算不出“高密度”时会出现什么奇怪现象（比如某些区域突然完全没人，或者突然挤满人）。

2. 新武器：离散版的"WKB 近似”

作者发明了一种新工具，可以把它想象成**“微观地形测绘仪”**。

• WKB 是什么？ 在物理学里，WKB 是一种处理“快速波动”的方法。想象你在听一段忽高忽低的音乐，WKB 不关心每一个具体的音符，而是关心旋律的整体起伏趋势。
• 作者的创新：他们直接把这种“听旋律”的方法用在了离散的座位上（而不是连续的波浪）。他们不依赖复杂的“气象卫星”（场论），而是直接分析每个座位上的“客人波函数”（即客人出现在那里的数学概率波）。

3. 核心发现：神奇的“密度公式”

通过这种新方法，作者推导出了一个万能公式。这个公式就像一把**“地形透视镜”**，只要告诉你：

1. 这个位置的座位间距（$J$）；
2. 这个位置的灯光亮度（$B$）；
3. 你总共请了多少客人（费米能级 $\epsilon_F$）；

它就能立刻告诉你：这个座位上坐满人的概率是多少？

公式给出了三种状态：

• 完全空置（Depletion/耗尽）：如果灯光太暗或座位太挤，客人根本进不去，概率为 0。就像在沙漠里，水根本存不住。
• 完全坐满（Saturation/饱和）：如果灯光太亮或座位太松，客人挤得满满当当，概率为 1。就像在早高峰的地铁里，连缝隙都塞满了。
• 中间状态：概率在 0 到 1 之间平滑过渡，由一个反余弦函数决定。

4. 有趣的比喻：彩虹链与克劳特库链

作者用几个具体的例子来验证这个公式：

• 彩虹链（Rainbow Chain）：想象一个房间，中间座位很宽，两边座位很窄（像彩虹拱门）。
  • 现象：当客人不多时，他们只敢待在中间宽敞的地方，两边完全没人（耗尽）。
  • 公式预测：完美吻合！
• 克劳特库链（Krawtchouk Chain）：想象座位分布像一个椭圆。
  • 现象：随着客人增加，他们会先从左边坐满，然后中间开始有人，最后右边坐满。
  • 公式预测：不仅算出了哪里坐满，还算出了哪里是空的，甚至算出了“半满”时的对称性。

5. 为什么这很重要？（纠缠熵的“消音器”）

这篇论文最酷的地方在于，它解释了为什么在某些不均匀的系统中，“纠缠”（量子纠缠）会突然消失。

• 什么是纠缠？ 想象两个客人，无论隔多远，他们的动作都是同步的。
• 发生了什么？ 当某些区域完全没人（耗尽）或完全坐满（饱和）时，这些区域就像**“隔音墙”**。
  • 如果左边区域完全没人（真空），它就像一张白纸，和右边没有任何联系。
  • 如果左边区域完全坐满，它也像一堵实墙，切断了联系。
• 结论：作者的新公式能精确画出这些“隔音墙”的位置。这意味着，我们终于能从数学上理解，为什么在不均匀的环境中，量子系统的“混乱程度”（纠缠熵）会被抑制。

总结

这就好比以前我们只能用模糊的望远镜看星星，知道大概有几颗，但看不清细节。
这篇论文给科学家发了一副**“高精度眼镜”**。戴上它，我们不仅能看清在不均匀的地形下，粒子（客人）到底分布在哪里，还能预测哪里会形成“真空区”，哪里会形成“拥堵区”。

这不仅解决了物理学的一个老难题，还为未来设计量子模拟器和量子计算机提供了重要的理论地图——告诉我们如何控制粒子的分布，从而控制量子信息的传递。

技术摘要

这是一份关于论文《Local fermion density in inhomogeneous free-fermion chains: a discrete WKB approach》（非均匀自由费米链中的局域费米子密度：离散 WKB 方法）的详细技术总结。

1. 研究背景与问题 (Problem)

• 研究对象：一维非均匀自由费米子链（即非均匀 XX 自旋链），其跳跃振幅 $J_n$ 和外部磁场 $B_n$ 随格点位置变化。
• 核心挑战：
  • 在热力学极限下，如何解析地计算非均匀链中的局域费米子密度（local fermion density）$\langle c^\dagger_n c_n \rangle$。
  • 现有的场论方法（如共形场论）通常仅适用于均匀系统、半填充（half-filling）或零磁场情况。对于任意填充率（filling fraction）和非零磁场的非均匀系统，缺乏通用的解析处理手段。
  • 数值模拟观察到了“耗尽”（depletion，密度趋近于 0）和“饱和”（saturation，密度趋近于 1）现象，这些现象与纠缠熵的抑制密切相关，但缺乏普适的解析公式来解释和预测这些区域的分布。
• 现有局限：之前的解析尝试（如 Mula 等人）仅适用于低填充率和零磁场，无法处理高填充率下的饱和效应或任意磁场情况。

2. 方法论 (Methodology)

作者提出了一种直接应用于离散模型本征函数递推关系的离散 WKB（Wentzel-Kramers-Brillouin）近似方法，避免了中间场论近似。

• 基本设定：

  • 考虑哈密顿量 $H = \sum J_n(c^\dagger_n c_{n+1} + h.c.) + \sum B_n c^\dagger_n c_n$。
  • 引入晶格常数 $a$，在热力学极限下（$N \to \infty, a \to 0, Na = \ell$ 固定），将离散格点 $n$ 映射为连续坐标 $x = na$。
  • 假设 $J_n$ 和 $B_n$ 平滑收敛为连续函数 $J(x)$ 和 $B(x)$。

• WKB 近似推导：

  1. 递推关系：单粒子本征函数满足三项递推关系。将其重写为对称形式，并取连续极限。
  2. WKB 试探解：假设波函数形式为 $\psi(x) = e^{i S(x)/a}$，其中 $S(x)$ 是展开为 $a$ 的幂级数的相位函数。
  3. 相位求解：通过逐阶分析，得到相位 $S_0(x)$ 和 $S_1(x)$ 的表达式。
     • 定义无量纲参数 $\xi(x, \epsilon) = \frac{\epsilon - B(x)}{2J(x)}$。
     • 当 $|\xi| < 1$ 时，波函数呈振荡形式（允许区）；当 $|\xi| > 1$ 时，波函数呈指数衰减或增长（禁戒区）。
  4. 波函数包络：推导出单粒子本征函数的 WKB 近似形式：
     $$\phi(x, \epsilon) \approx \frac{A(\epsilon) \sin \varphi^*(x, \epsilon)}{\sqrt{J(x) |1 - \xi(x, \epsilon)^2|^{1/2}}}$$
     其中 $\varphi^*$ 是累积相位，$A(\epsilon)$ 是归一化常数。在禁戒区，波函数近似为零。

• 密度计算：

  • 利用态密度 $D(\epsilon)$ 和本征函数的平方，将求和转化为积分。
  • 通过对快速振荡项 $\sin^2 \varphi^*$ 进行平均，导出局域费米子密度的解析表达式。

3. 关键贡献与核心公式 (Key Contributions & Results)

论文导出了一个通用的闭式解析公式，用于计算任意填充率、任意跳跃振幅分布和任意磁场分布下的局域费米子密度。

• 核心公式 (Eq. 5 & Eq. 50)：
  局域费米子密度（以格点占据数形式表示）近似为：
  $$\langle c^\dagger_n c_n \rangle \simeq \begin{cases} 0, & \epsilon_F \le B_n - 2J_n \quad (\text{耗尽区}) \\ \frac{1}{\pi} \arccos\left( \frac{B_n - \epsilon_F}{2J_n} \right), & B_n - 2J_n \le \epsilon_F \le B_n + 2J_n \\ 1, & \epsilon_F \ge B_n + 2J_n \quad (\text{饱和区}) \end{cases}$$
  其中 $\epsilon_F$ 是费米能级。

• 主要发现：

  1. 耗尽与饱和的解析判据：
     • 耗尽区 (Depletion)：当费米能级低于局部能带底 ($B(x) - 2J(x)$) 时，密度为 0。
     • 饱和区 (Saturation)：当费米能级高于局部能带顶 ($B(x) + 2J(x)$) 时，密度为 1。
     • 这一发现解释了之前数值模拟中观察到的纠缠熵抑制现象：在耗尽或饱和区域，子系统的态是纯态（真空或全满），因此纠缠熵为零。
  2. 任意填充率的适用性：该公式不仅适用于低填充（与之前文献一致），还成功描述了高填充下的“镜像”耗尽（即饱和）现象，这是之前启发式公式无法做到的。
  3. 态密度与填充率：推导了态密度 $D(\epsilon)$ 和填充率 $\nu(\epsilon_F)$ 的解析表达式，建立了费米能级与总粒子数的关系。

4. 验证与示例 (Examples & Validation)

作者在多种典型的非均匀链模型中验证了该公式的准确性，并与数值模拟结果进行了对比：

1. 均匀链 (Homogeneous chain)：公式退化为常数密度，与精确解一致，并解释了边界处的 Friedel 振荡。
2. Krawtchouk 链：展示了在不同填充率下，链两端出现的耗尽区或饱和区，以及中间区域的平滑过渡。
3. 彩虹链 (Rainbow chain)：针对指数衰减的跳跃振幅，公式准确预测了低能级下的双端耗尽区和高能级下的无耗尽区。
4. 余弦链 (Cosine chain)：展示了跳跃振幅周期性变化导致的复杂密度分布，包括多个耗尽/饱和区间。
5. 非对称余弦链：在同时存在非对称跳跃和抛物线磁场的情况下，公式依然表现出极高的精度，能够处理多个不连续的势阱（wells）。

在所有案例中，解析公式（红线）与数值计算结果（蓝点）在整个能量和填充率范围内均表现出极好的一致性。

5. 意义与展望 (Significance & Outlook)

• 理论突破：提供了一种不依赖共形场论（CFT）的解析工具，能够处理非均匀、非临界（away from criticality）且存在磁场的自由费米子系统。
• 纠缠熵研究：由于自由费米子系统的纠缠熵完全由关联矩阵决定，而关联矩阵又由局域密度和本征函数决定，该 WKB 方法为解析计算非均匀系统的纠缠熵提供了基础。这对于理解非均匀环境下的纠缠面积律（Area Law）及其修正至关重要。
• 普适性：该方法不仅适用于一维链，其离散 WKB 的思想有望推广到更复杂的场景，如包含非局域效应的交换关联能计算等。
• 物理洞察：清晰地揭示了局域磁场和跳跃振幅的空间变化如何直接调控费米子的局域分布，从而控制量子纠缠的分布。

总结：这篇论文通过引入一种创新的离散 WKB 方法，成功解决了非均匀自由费米链中局域费米子密度的解析计算难题，给出了一个普适、精确且物理图像清晰的公式，为深入研究非均匀量子多体系统的纠缠性质奠定了坚实的理论基础。