Analytical classification of Majorana zero-mode spatial profiles in extended Kitaev chains: probability maxima can shift inward

本文提出了一个扩展的 Kitaev 链的分析框架，揭示了马约拉纳零模可以呈现多样的空间分布特征，包括内部概率极大值和不同的衰减行为，这些特征完全由导出的递推关系的特征根所决定。

要点

想象一条长长的、一维的原子链，就像一串珠子。在量子物理世界中，这被称为Kitaev 链。科学家们对这条链非常感兴趣，因为它可以容纳一种特殊的“幽灵粒子”，称为马约拉纳零能模（MZMs）。

可以将这些 MZM 想象为不可见的、零能量的幽灵，它们喜欢躲藏在链的最两端。由于它们位于两端，彼此相距甚远，这使得它们非常稳定，并且对于构建未来需要免受错误干扰的量子计算机非常有用。

通常，物理学家使用“拓扑不变量”（一个复杂的数学数值）来统计这些幽灵的数量。如果该数值为 1，则链的两端各有一个幽灵；如果为 2，则有两个。但这里有个问题：这个数值只能告诉你有多少个幽灵，却无法告诉你它们确切躲在哪里，或者它们长什么样。

本文就像一部侦探故事，深入观察这些幽灵的实际“形状”和“位置”，揭示了一些令人惊讶的秘密。

主要发现：幽灵并不总是守在门口

在最简单的模型中，科学家曾假设这些幽灵总是完美地粘在链的第一个或最后一个珠子上。他们认为“概率”（即发现幽灵的几率）在边缘处最高，并随着向内部深入而平滑地衰减。

本文证明这并不总是正确的。

通过一个巧妙的数学技巧（将问题转化为一组重复模式，或称“递推关系”），作者发现，根据链的“设置”不同，这些幽灵可以表现出三种截然不同的行为：

1. 单调幽灵：它表现得符合预期。它在边缘处最强，并随着向链深处深入而平滑衰减。
2. 振荡幽灵：它在衰减过程中会波动。想象一波海浪，随着远离海岸而变得越来越小。幽灵的存在感在穿透链的过程中上下起伏。
3. 完美局域化幽灵：在某些特殊情况下，幽灵根本不会逐渐衰减。它严格局限于第一个或前两个珠子上，就像一束聚光灯只照亮楼梯的最第一级台阶，而照不到其他地方。

大惊喜：“偏移”的幽灵

最令人兴奋的发现是，幽灵不必在边缘处最强。

想象你在一条长长的走廊里寻找一只走失的猫。你本以为它就在前门。但在这篇论文中，作者表明，这只猫（马约拉纳模）实际上最有可能出现在走廊深处两三个房间的位置，尽管它仍然属于前门。

• 隐喻：将幽灵想象成从隧道尽头扬声器发出的声波。通常，声音在扬声器处最响亮。但在这些特定的量子链中，声波可以相互干涉，从而在隧道内部几米处产生一个“响亮点”（概率最大值），即使声源位于墙壁处。
• 包络：尽管“响亮点”在内部，但随着你向更深处走去，以及向墙壁方向退去，声音仍然会衰减。它仍然是一个“边界”幽灵，但其峰值已向内偏移。

这对真实实验为何重要

在现实世界中，我们无法建造无限长的链；我们的链是有限的（短的）。当链较短时，来自左端和右端的幽灵会相互“碰撞”，混合它们的身份，使它们变得略微不完美。

作者提供了一个数学“尺子”（基于他们方程的根），告诉科学家：

• 链需要多长，才能在不被两端干扰的情况下观察到幽灵的“真实”形状。
• 去哪里寻找。如果你是一位试图寻找这些幽灵的实验人员，你不应该只盯着第一个原子看。你可能需要向链内部深入几个原子去寻找，因为幽灵的“峰值”可能正躲藏在那里。

一句话总结

• 问题：我们知道这些链中存在多少个量子幽灵，但我们不知道它们确切长什么样，或者它们的“心脏”在哪里。
• 解决方案：作者通过求解数学方程，描述了这些幽灵的确切形状。
• 转折：这些幽灵并不总是粘在边缘。它们可以波动，可以完美地粘在第一个珠子上，或者在即使没有缺陷的完全均匀系统中，其“最强点”也会向链深处偏移。
• 启示：如果你在搜寻这些粒子，不要只盯着边缘看。向深处多看一点，因为粒子的“峰值”可能正躲藏在那里，等待被发现。

技术摘要

技术摘要：扩展 Kitaev 链中马约拉纳零模空间分布的解析分类

问题陈述
尽管拓扑不变量（如缠绕数）能够成功预测一维超导系统中马约拉纳零模（MZMs）的数量，但它们无法确定这些模式的空间结构或局域化性质。在本质上有限的实验实现中，区分 MZM 的内在属性与有限尺寸效应（如边界诱导的杂化和能级分裂）仍然是一个挑战。现有研究通常依赖于有限系统的数值对角化或特定极限情况，缺乏一个通用的解析框架来刻画具有次近邻相互作用的扩展 Kitaev 链中 MZM 的内在空间分布。具体而言，尚不清楚系统参数如何影响 MZM 是呈单调衰减、振荡衰减，还是在远离链边缘处出现概率峰值的偏移。

方法论
作者分析了具有最近邻（$\lambda_1$）和次近邻（$\lambda_2$）耦合的扩展 Kitaev 链，重点关注跃迁振幅与配对振幅相等（$t_i = \Delta_i = \lambda_i$）的特定参数区域。

1. 马约拉纳基变换：将哈密顿量重写为马约拉纳算符基（$a_i, b_i$）。该变换通过将马约拉纳振幅的运动方程解耦，简化了问题。
2. 递推关系推导：针对零能解（$E=0$），作者推导出了马约拉纳模式振幅 $A_j$ 的线性递推关系。通解表示为涉及特征二次方程 $\lambda_2 q^2 + \lambda_1 q - \mu = 0$ 的根（$q_\pm$）的项的线性组合。
3. 半无限极限分析：为了将内在属性与有限尺寸杂化隔离开来，分析在半无限极限下进行。MZM 的存在及其数量由归一化条件 $|q| < 1$ 决定。
4. 根的归类：根据特征根的性质（实数与复数、互异与简并）及其模长对空间分布进行分类。作者推导出了各种参数区域下振幅 $A_j$ 的闭式表达式。
5. 有限尺寸对比：将半无限极限的解析预测与有限链的数值结果进行比较，以确定有限系统重现内在空间结构所需的长度尺度。

主要贡献与结果

• 空间分布的解析表征：本文确立了 MZM 的空间结构完全由递推关系的特征根决定。这建立了哈密顿量参数（$\mu, \lambda_1, \lambda_2$）与模式衰减行为之间的直接联系。
• 多样的衰减行为：在单一拓扑相（固定缠绕数）内，MZM 可表现出截然不同的空间行为：
  • 单调衰减：概率从边缘开始指数衰减。
  • 振荡衰减：概率在衰减过程中发生振荡。
  • 完美局域化：模式被限制在有限数量的格点上（例如第一个或前两个格点），没有指数拖尾。
• 向内偏移的概率峰值：一个核心发现是，源于边界的 MZM 其最大概率不一定位于链边缘（$j=1$）。取决于递推关系两个独立解之间的干涉（由 $A_1$ 和 $A_2$ 的相对权重决定），概率峰值可以偏移至内部格点。这些模式在峰值两侧仍保持指数衰减包络，证实了其源于边界的本质，尽管发生了偏移。
• 相图分类：作者根据可归一化根的数量（$w=0, 1, 2$）将参数空间映射为区域（$G_1$ 至 $G_4$）。他们确定了根从实数转变为复数或其模长跨越单位圆的临界线，这些标志着拓扑相变。
• 有限尺寸交叉：该研究确定了由 $|q|$ 导出的、依赖于参数的长度尺度，这些尺度决定了有限链何时足够大以展现出内在的半无限分布。它阐明了小链中的“大质量边缘模式”是有限尺寸伪影，随着系统尺寸超过穿透深度，它们会演变为理想的 MZM。

意义与主张
本文声称提供了首个通用的解析框架，能够直接表征扩展 Kitaev 链中 MZM 的空间结构，且独立于系统尺寸。

• 与拓扑不变量的区别：该工作表明，虽然拓扑不变量固定了 MZM 的数量，但空间分布是系统参数的连续函数。这种变异性包括峰值偏移的可能性，这对仅通过第一个格点处的峰值来识别边缘模式的实验启发式方法提出了挑战。
• 内在效应与有限尺寸效应：通过求解半无限极限，作者区分了内在哈密顿量特征（如由干涉引起的向内偏移峰值）与有限尺寸效应（如杂化诱导的能级分裂）。这种区分对于解释系统尺寸受限的实验数据至关重要。
• 鲁棒性：作者指出，只要体带隙保持开启且保护对称性得以保留，定性特征（边界模式的存在以及峰值偏移的可能性）预计会在弱无序下持续存在，尽管具体的空间分布可能会发生改变。
• 方法论优势：与以往依赖数值扫描或特定约束的工作不同，这种解析方法允许对整个相图中的空间分布进行完整分类，揭示了如向内偏移峰值等此前在有限链的解析处理中未被观察到的现象。

本文得出结论：MZM 的空间结构由多个衰减分量的组合与干涉决定，即使在不改变拓扑相的情况下，也能产生一系列截然不同的分布。该框架建立了哈密顿量参数与实验可观测的空间分布之间的直接联系。