非Hermitian皮肤效应如何点亮莫尔超导？深度解析三角Hubbard系统中的配对关联增强机制

来源论文: https://arxiv.org/abs/2606.20425v1 生成时间: Jun 20, 2026 00:45

0. 执行摘要

非Hermitian物理与强关联电子学是当代凝聚态物理中最活跃的两大前沿领域。非Hermitian皮肤效应（Non-Hermitian Skin Effect, NHSE）作为非Hermitian系统特有的拓扑物理现象，表现为在开边界条件（OBC）下，系统所有的体态波函数均会发生指数级的退化与边界坍缩。与此同时，莫尔超晶格（Moiré Superlattices，如魔角双层石墨烯MATBG、过渡金属硫族化合物tTMDs）因其高度可调的平带结构，为研究强关联超导、Mott绝缘体等非常规超导物理提供了完美的实验平台。然而，如何利用非Hermitian不平衡驱动（如非互易跃迁）作为一种新型的量子调控手段去积极调控甚至增强强关联莫尔系统中的超导配对关联，一直是一个极具挑战性且悬而未决的课题。

本博文深度解析了最新研究成果《Non-Hermitian Skin Effect Enhances Pairing Correlations in Moiré Hubbard Systems》。该工作巧妙地将非Hermitian Hatano-Nelson非互易跃迁引入到三角莫尔Hubbard模型中，利用精确对角化（ED）与非Hermitian密度矩阵重整化群（NH-DMRG）方法，首次揭示了NHSE能够通过**通道选择性机制（Channel-Selective Mechanism）**大幅度增强系统局域配对关联。研究表明，在适度的非互易性区间 $\gamma \in [0.5, 1.2]t$（被称为“黄金窗口”）内，边界电荷的指数级堆积使得局域态密度（LDOS）获得极大的放大，进而强力催化了局域的s波配对；与此同时，NHSE对竞争的反铁磁（AFM）关联起到了显著的压制作用。这两者的协同效应使得系统的总配对易受性 $\chi_{SC}$ 实现了高达 +98% 的惊人增长。这一发现不仅为非平衡态超导调控指明了全新的物化路径，也为未来的莫尔纳米器件设计提供了理论支撑。

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

1.1 核心科学问题

在平衡态强关联物理中，调控超导配对关联（Pairing Correlations）主要依赖于调整载流子浓度、施加外加压力、或通过调谐扭角来改变带宽与库仑能之比 $U/t$。这些手段均局限于Hermitian框架。而在非平衡开放量子系统中，不平衡的能量或粒子注入通常被视为超导的敌人，因为耗散与非相干去相位容易破坏脆弱的库珀对。

本研究探讨的核心科学问题是：非Hermitian皮肤效应（NHSE）引起的电荷空间重分布，究竟是仅仅简单地重标度了局域电荷密度（平庸的电荷堆积），还是能够作为一种积极的量子相干机制，从根本上重组多体相互作用通道，从而增强真实的许多体配对关联？ 进一步地，在强库仑排斥（Hubbard $U$）与非互易跃迁产生的电荷坍缩之间，存在怎样的竞争与协同？

1.2 理论基础：非互易性三角莫尔Hubbard模型

本工作的理论大厦建立在Hatano-Nelson非互易跃迁与莫尔Hubbard模型的结合之上。考虑二维三角莫尔超晶格上的电子系统，其非Hermitian多体哈密顿量由三部分构成：

1. 动力学项（Kinetic Term, $H_{kin}$）：描述相邻莫尔超晶格点之间的各向同性跃迁以及化学势调节：

   $$H_{kin} = -t \sum_{\langle ij \rangle, \sigma} \left( c^\dagger_{i\sigma} c_{j\sigma} + H.c. \right) - \mu \sum_{i, \sigma} n_{i\sigma}$$

   其中 $t$ 代表莫尔格点间的有效跃迁（在莫尔超晶格中，每个格点对应一个莫尔超胞，其有效跃迁通常在 $1\sim 10 \text{ meV}$ 数量级），$\mu$ 为化学势，用于调节系统处于半满（Half-filling）状态。

2. 强关联相互作用项（Hubbard Term, $H_U$）：表征莫尔平带带来的强电子-电子排斥能：

   $$H_U = U \sum_{i} n_{i\uparrow} n_{i\downarrow}$$

   这里 $U/t \in [0, 12]$，完整地覆盖了从弱耦合金属到强关联Mott绝缘体的物理区间。

3. 非Hermitian非互易跃迁项（Non-Hermitian Term, $H_{NH}$）：沿三角晶格的 $\hat{x}$ 方向引入Hatano-Nelson型非互易跳跃：

   $$H_{NH} = \gamma \sum_{\langle ij \rangle_x, \sigma} c^\dagger_{i\sigma} c_{j\sigma}$$

   此处 $\gamma > 0$。这导致向右跃迁的有效振幅变为 $t + \gamma$，而向左跃迁的有效振幅衰减为 $t - \gamma$。这种跃迁的非对称性直接打破了哈密顿量的Hermiticity（$H^\dagger \neq H$），在开边界条件（OBC）下会诱导全部单粒子及多体本征态向右侧边界指数级坍缩。

1.3 技术难点与非Hermitian多体物理的挑战

将强关联Hubbard模型与非Hermitian物理相结合，在理论和计算上面临极大的技术瓶颈：

• 双正交基底的失效与右本征态的选择：在非Hermitian量子力学中，哈密顿量的左本征态 $\langle \psi_L |$ 与右本征态 $|\psi_R \rangle$ 不再互为共轭。传统的厄米变分原理和热力学期望值定义面临修正。对于多体本征值问题： $$H |\psi_{R, n}\rangle = E_n |\psi_{R, n}\rangle, \quad H^\dagger |\psi_{L, n}\rangle = E_n^* |\psi_{L, n}\rangle$$ 研究必须在“双正交表象”（Biorthogonal Formalism, $\langle \psi_L | \hat{O} | \psi_R \rangle$）与“右本征态表象”（Right-eigenstate Formalism, $\langle \psi_R | \hat{O} | \psi_R \rangle / \langle \psi_R | \dots \rangle$）之间做出严格的论证。当系统靠近非Hermitian退化点（Exceptional Point, EP）或发生局域能级交叉时，双正交重叠 $\langle \psi_L | \psi_R \rangle \to 0$，会导致计算数值发散。如何稳健地提取多体基态物理量是一大难点。
• 电荷过度局域化（Over-localization）的压制：当非互易跃迁强度 $\gamma$ 过大时，所有电子会被极端地压缩至最右侧的一两个格点上。这种局域的电荷饱和（即格点被双重占满，无法进行相干跃迁）会导致泡利不相容原理主导的物理“死锁”，反而会剧烈压制相干的配对关联。因此，必须寻找一个平衡电荷堆积与相干跃迁的“黄金窗口”。
• 指数级的希尔伯特空间与非对称Arnoldi算法的收敛性：即使在固定的粒子数扇区（如半满且 $S_z = 0$ 守恒），三角晶格的多体基底数也随格点数指数增长。对于非对易稀疏矩阵，常用的Hermitian Lanczos方法无法保证收敛到具有最小实部的“多体基态”。必须采用基于Krylov子空间的Arnoldi算法（如ARPACK中的 eigs），但复数本征值的排序和追踪极易受到近简并能级的干扰。

1.4 方法细节

为了确保物理结论的坚实可靠，作者采用了多重先进的数值方案相互印证：

精确对角化（Exact Diagonalization, ED）

作者在 $L_x \times L_y$ 的三角格点系统（点数 $N_s = 4 \sim 12$）上实现了完全粒子数守恒的ED计算。通过固定向上自旋和向下自旋的电子数（$N_\uparrow, N_\downarrow$），构建出完全正交的多体职业数基底。哈密顿量被拆解为：

利用位运算（XOR、popcount）极速构建稀疏哈密顿量。通过ARPACK库，在复平面内精确搜寻具有最小能量实部的右本征态 $|\psi_R\rangle$ 作为非Hermitian多体基态。同时求解 $H^\dagger$ 以获取对应的左本征态 $|\psi_L\rangle$ 进行双正交基分析。

非Hermitian密度矩阵重整化群（NH-DMRG）

为了克服ED的尺寸限制，作者引入了基于矩阵乘积态（MPS）的非Hermitian DMRG算法。该算法核心在于采用双正交MPS表象（Biorthogonal MPS），即分别用两套不同的张量矩阵 $\{A^{[i]}_L\}$ 和 $\{A^{[i]}_R\}$ 来表示左、右多体波函数，从而能够高精度、无偏地计算中等尺寸（可达 $N_s = 24$ 的三角格点柱状体）的边界密度分布和双占有率。

耗散开放系统：Lindblad主方程验证

为了排除单纯哈密顿量非Hermiticity的数学人工痕迹，作者还采用了一种更具实验基础的微观描述——Lindblad主方程方法，在 $2 \times 2$ 晶格上求解非平衡稳态（NESS）：

通过引入具有空间梯度的粒子注入与抽失跳跃算符 $L_i^{\pm}$，在纯正的开放系统框架下，成功复现了NHSE驱动的电荷边界堆积与配对关联演化行为。

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

2.1 Benchmark体系配置

本工作的核心测试体系是一个带有开边界条件（OBC）的 $3 \times 3$ 三角超晶格模型（格点数 $N_s = 9$），在半满扇区（$N_\uparrow = 5, N_\downarrow = 4$ 或自旋对称区 $N_\uparrow = 4, N_\downarrow = 4$），其对应的多体希尔伯特空间维度为 $15,876$。该体系虽然尺寸有限，但由于三角晶格的高度配位数（每个内格点有6个邻居），能够完美保留超导配对与反铁磁磁性关联之间的强烈竞争。

2.2 核心物理量定义

为了量化分析，定义了以下两个关键阶参量：

1. 双占有率（Double Occupancy, $D$）：作为局域s波配对易受性的直接代用指标： $$D = \frac{1}{N_s} \sum_i \langle n_{i\uparrow} n_{i\downarrow} \rangle_R$$
2. 皮肤效应阶参量（Skin Order, $S$）：用于度量电荷在非互易跃迁下的非对称重分布： $$S = \frac{1}{N_s} \sum_i (x_i - \bar{x}) \langle n_i \rangle_R$$ 其中 $x_i$ 是格点的横坐标，$\bar{x}$ 是系统几何中心。

2.3 关键计算数据分析

黄金窗口中的非单调双占有率增强（$D(\gamma)$）

在传统的Hermitian系统中，增加电子跃迁的不对称性通常会削弱关联效应。但在该系统中，作者观察到了惊人的非单调（Dome-shaped）行为（参见论文 Fig. 3(a)）：

物理机制解析：在 $\gamma < \gamma^*$ 区间，非互易跃迁驱动电子向右侧边界堆积，导致右边界的局域态密度（LDOS）发生指数级爆发（如 $L=9, \gamma=0.3t$ 时，边界LDOS放大系数高达 $10^2$）。这种极高的局域电荷密度在库仑作用 $U$ 的协助下，强力促进了局域电子配对，双占率随之攀升。然而，当 $\gamma > \gamma^*$ 时，电子被过度压缩（Over-localization）在极少数格点上，导致格点完全被双重占满（电荷死锁），泡利不相容原理阻止了电子的相干流动，导致配对关联迅速崩塌。因此，系统存在一个最优调控的**“黄金窗口”**。

配对易受性 $\chi_{SC}$ 的翻倍增长与通道选择性（Channel-selective）机制

为了证明双占率的提升并非平庸的静电电荷聚集，作者计算了s波自旋单态配对涨落的关联函数：

以及总配对易受性（Pairing Susceptibility）：

将 $\chi_{SC}$ 分解为原位（On-site）配对通道 $\chi_{on} = N_s D$ 与离位（Off-site）相干通道 $\chi_{off} = \sum_{i \neq j} P(i,j)$。计算所得核心性能数据如下（在 $U=4t$ 体系下，参见论文 Fig. 4）：

• 总配对易受性 $\chi_{SC}$：从 $\gamma = 0$ 时的 0.44 一路飙升至 $\gamma^* = 1.05t$ 时的 0.87，相对增强幅度高达惊人的 +98%！
• 原位通道 $\chi_{on}$：在 $\gamma^*$ 处增长了 +21%（贡献了正向的配对配给）。
• 离位通道 $\chi_{off}$：在三角晶格中，由于反铁磁涨落的竞争，离位配对关联通常为负值（表示排斥）。令人振奋的是，在非互易性作用下，该负值的绝对值收缩了 22%（从负得较多变为负得较少），这意味着竞争的反铁磁通道受到了极大的压制。

这一结果完美揭示了NHSE的通道选择性物理机制：非非对易的皮肤效应空间局域化不仅在空间上堆积了电荷（增强了原位配对），而且极大地打乱并压制了在三角晶格中占据主导地位的竞争性磁性关联（反铁磁涨落）。这两者的建设性干涉，导致总超导配对易受性实现了近乎翻倍的飞跃。相比之下，在 $U=0$ 的无相互作用极限下，原位通道仅增长 +10.9%，且无大 $\gamma$ 下的陡峭崩塌。这强力证明了黄金窗口是一个真正的多体强关联相干效应，无法被平庸的单粒子密度红利所解释。

NH-DMRG 大尺寸 Benchmark 数据

为排除有限尺寸效应，作者利用 NH-DMRG 对 $3 \times L_y$ 和 $4 \times 4$ 的三角柱状格进行计算（参见 Fig. 7 & Fig. 8）：

• 在 $3 \times 12$ ($N_s=36$) 体系上，双占率依然在 $\gamma^* \approx 0.9t$ 处展示出明显的非单调“峰”结构，增强幅度维持在 +15% ~ +18%。
• 在 $4 \times 4$ ($N_s=16$) 体系下（Fig. 8），多 $U$ 扫描展现出高精度的峰值结构：$U=2t$ 时配对增强为 +13.1%，$U=4t$ 时为 +15.5%，而强耦合极限 $U=6t$ 时受到强烈磁性锁止，增强被压制至 +1.7%。
• 利用不同几何尺寸（$N_s = 4 \sim 12$ 的7种几何形状）进行 $1/N_s$ 有限尺寸标度外推（Fig. 6），在热力学极限下（$1/N_s \to 0$），该黄金窗口带来的配对增强仍稳定外推至约 +13%，确立了其热力学稳定性。

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

3.1 代码复现整体架构与核心逻辑

要复现该非Hermitian三角格点Hubbard模型的ED计算，核心任务在于建立支持非互易跃迁的稀疏矩阵，并利用非对称特征值求解器。以下基于 Python 的主流科学计算库 SciPy 给出核心实现步骤和伪代码：

步骤 1：构建 Fock 空间基底与位运算编码

对于 $N_s = 9$ 且处于半满自旋守恒扇区的系统（如 $N_\uparrow=4, N_\downarrow=4$），我们可以使用一个 int64 整数来编码电子状态。前 9 位存储向上自旋（$\uparrow$）的职业数，接下来的 9 位存储向下自旋（$\downarrow$）的职业数。

import numpy as np
from scipy.sparse import csr_matrix
import scipy.sparse.linalg as jla

def get_basis(n_sites, n_particles):
    """生成特定粒子数扇区的所有二进制状态"""
    basis = []
    for state in range(1 << n_sites):
        if bin(state).count('1') == n_particles:
            basis.append(state)
    return np.array(basis, dtype=np.int64)

步骤 2：定义格点几何与非互易邻接矩阵

在三角格点上，确定横向（$\hat{x}$）非互易键和其余厄米键。定义在横向键上向右跳跃系数为 $t+\gamma$，向左跳跃为 $t-\gamma$，纵向或斜向跳跃系数保持 $t$：

def get_hopping_amplitude(i, j, gamma, t=1.0):
    # 假设格点位置 (x_i, y_i), (x_j, y_j) 已知
    # 检测是否为横向的最近邻键
    if is_x_neighbor(i, j):
        if x[j] > x[i]: # 向右跃迁
            return -(t + gamma)
        else:           # 向左跃迁
            return -(t - gamma)
    elif is_other_neighbor(i, j):
        return -t
    return 0.0

步骤 3：稀疏哈密顿量矩阵组装（极其关键）

利用多体基底之间的单体跃迁算符作用，计算哈密顿量矩阵元。注意非对易性导致哈密顿量不再是对称矩阵：

def build_hamiltonian(basis_up, basis_down, U, gamma, n_sites):
    dim_up = len(basis_up)
    dim_down = len(basis_down)
    dim = dim_up * dim_down
    
    # 用于构建CSR稀疏矩阵的临时数组
    data, rows, cols = [], [], []
    
    # 1. 构建对角相互作用项 HU
    for i_up, s_up in enumerate(basis_up):
        for i_down, s_down in enumerate(basis_down):
            row_idx = i_up * dim_down + i_down
            # 计算双占格点数
            double_occ = bin(s_up & s_down).count('1')
            u_energy = U * double_occ
            if u_energy != 0:
                data.append(u_energy)
                rows.append(row_idx)
                cols.append(row_idx)
                
    # 2. 构建非对易动力学与非互易项 (H_kin + H_NH)
    # 这里分别对自旋向上和自旋向下的基底作用跃迁，利用Kronecker积的思想加速
    # 伪代码：遍历所有多体态跃迁...
    # 跃迁作用细节需要考虑费米子交换反对称化（负号确定）
    # ... (略去繁琐的费米负号追踪过程)
    
    H = csr_matrix((data, (rows, cols)), shape=(dim, dim), dtype=np.complex128)
    return H

步骤 4：非对称特征值求解与本征态追踪

在非Hermitian系统中，我们寻找具有**最小实部（Smallest Real Part）**的本征态：

def solve_ground_state(H):
    # 使用 ARPACK 求解器，which='SR' 代表寻找实部最小的本征值
    eigenvalues, eigenvectors = jla.eigs(H, k=6, which='SR')
    
    # 排序
    idx = np.argsort(np.real(eigenvalues))
    gs_energy = eigenvalues[idx[0]]
    psi_R = eigenvectors[:, idx[0]]
    
    # 归一化右本征态
    psi_R /= np.linalg.norm(psi_R)
    
    # 求解对应的左本征态
    H_adj = H.conjugate().transpose()
    evals_L, evecs_L = jla.eigs(H_adj, k=6, which='SR')
    # 在简并或能级交叉时需要小心匹配左本征态，使其满足 <psi_L|psi_R> != 0
    # ...
    return gs_energy, psi_R

步骤 5：观测计算

def compute_double_occupancy(psi_R, basis_up, basis_down, n_sites):
    # 计算 <psi_R | n_i_up * n_i_down | psi_R>
    dim_down = len(basis_down)
    n_occ = np.zeros(n_sites)
    for i_up, s_up in enumerate(basis_up):
        for i_down, s_down in enumerate(basis_down):
            idx = i_up * dim_down + i_down
            coeff = np.abs(psi_R[idx])**2
            for site in range(n_sites):
                if ((s_up >> site) & 1) and ((s_down >> site) & 1):
                    n_occ[site] += coeff
    return np.mean(n_occ)

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

对于更大型或更高效的复现，建议直接使用社区内已经高度成熟的非Hermitian多体计算框架：

1. ITensors (Julia/C++)：功能强大的张量网络库。其支持高度自定义的算符，非常适合用于实现本工作中的非Hermitian DMRG（NH-DMRG）。利用其自带的非对称Lanczos/Krylov方法，可高效更新双正交矩阵积态（Biorthogonal MPS）。
   • Link: ITensors.jl GitHub
2. QuSpin (Python)：专门针对一维和二维量子多体系统（包括Hubbard、Spin、Bose-Hubbard等）的精确对角化和动力学演化的成熟Python包。QuSpin原生支持非Hermitian算符的构建与对角化，其底层的C++加速机制能大幅节约内存与计算时间。
   • Link: QuSpin GitHub
3. Biorthogonal MPS Solver：论文作者在NH-DMRG中参考了基于矩阵乘积态的双正交方法，该方法在张量网络学界已有很多开源的实现基础。可以参考凝聚态物理领域的非Hermitian强关联算法库：
   • Link: TeNPy (Tensor Network Python) Library

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

4.1 关键引用文献

本工作立足于以下几项奠基性研究之上，构成其核心立论链条：

• [1] N. Hatano and D. R. Nelson, Phys. Rev. Lett. 77, 570 (1996)：非互易跳跃（Hatano-Nelson模型）的开山之作。首次论证了不对称跳跃如何导致局域波函数的指数级定位，这是NHSE的理论源头。
• [2] S. Yao and Z. Wang, Phys. Rev. Lett. 121, 086803 (2018)：非Bloch能带理论的建立，彻底奠定了在开边界条件下如何处理非Hermitian皮肤效应（NHSE）的量子力学数学框架。
• [6, 7] Y. Cao et al., Nature 556, 43 & 80 (2018)：魔角双层石墨烯（MATBG）强关联与非常规超导的实验发现，为本工作提供了真实的莫尔Hubbard物理载体。
• [16] X.-J. Yu et al., Phys. Rev. Lett. 132, 116503 (2024)：研究了周期性边界条件（PBC）下非Hermitian蜂窝格点Hubbard模型的强关联物理。该研究发现非Hermiticity能够通过“费米速度还原”促进反铁磁性。本工作与之形成鲜明对比——本工作强调在**开边界条件（OBC）**下，NHSE带来的空间重新分布反而会压制反铁磁性，促进超导配对，充分展示了边界条件在非Hermitian强关联系统中的决定性作用。
• [20] P. Zhong et al., arXiv:2401.15000：提出了双正交矩阵积态（Biorthogonal MPS）的DMRG新算法，为本工作大尺寸NH-DMRG计算提供了直接的方法论支撑。

4.2 局限性深度点评（Critical Commentary）

作为一项发表于前沿的理论探索，该工作机制新颖，计算严谨，但在走向真实器件应用与完整理论建构的过程中，仍存在以下值得深入探讨和批判的局限性：

1. 物理效应在“热力学极限”下的稀释问题

非Hermitian皮肤效应（NHSE）本质上是一个**边界驱动（Boundary-driven）**的现象。电荷堆积主要发生在距离最右侧边界 $\lambda_{skin}$ 的范围内。在作者计算的小尺寸（$3 \times 3$ 到 $4 \times 4$）体系中，边界格点占比极大（如 $3 \times 3$ 中有 8 个格点位于边界或紧邻边界），因此整个系统都深度感受到了局域态密度的放大。然而，随着二维系统尺寸 $L \to \infty$，边界占体相比例 $\sim 1/L$ 会趋于零。这意味着，在宏观的热力学极限下，局域超导增强效应会被“ bulk ”严重稀释，这并非真正意义上的体相超导电性（Bulk Superconductivity）。这限制了其在传统宏观超导材料中的应用，使其更偏向于一种介观纳米器件（Mesoscopic Moiré Devices）的局域配对增强。

2. Floquet 相干驱动引发的“热死”（Floquet Heating）效应

论文提出，可以通过相干激光驱动（高频 Floquet 调制）来动态工程化出 Hatano-Nelson 型非互易跃迁 $\gamma$。然而，在具有强库仑相互作用 $U$ 的强关联莫尔Hubbard系统中，高频周期驱动会引发严重的 Floquet 吸收与电热效应（Floquet Heating）。强烈的热激发和多体态加热会极快地破坏极为脆弱的超导配对（莫尔系统的超导能隙通常仅有几毫开尔文 $\sim \text{mK}$）。如何平衡非对易非平衡驱动的相干性与系统升温，是该方案在实验上落地时最难跨越的鸿沟。

3. 右本征态与双正交观测量的物理诠释分歧

非Hermitian多体系统中，物理可观测量的最佳数学定义在物理学界仍有争议。作者在论文中主要采用右本征态期望值 $\langle \hat{O} \rangle_R$，并辅以双正交期望值 $\langle \hat{O} \rangle_{bi}$ 进行印证。如 Table I 所示，在 $\gamma/t = 0.5$ 和 $1.0$ 附近，由于系统的多体近简并导致双正交重叠矩阵极其接近奇异值，双正交范数 $\langle \psi_L | \dots \rangle \to 0$，发生了严重的数值发散（Divergence）。尽管作者声称这不影响物理 dome 区间的存在，但这种数学上的非解析性也暗示了在接近 Exceptional Point 附近时，传统的量子测量和物理解释面临逻辑上的不完备性。

4. 真实实验中 $\gamma$ 的高精度调谐可行性

在二维材料器件中，实现单向非互易跃迁需要精准打破空间反演对称性与时间反演对称性。利用外加电极提供栅极偏压和外加不对称耗散层虽在理论上可行，但在纳米级精度的莫尔格点上维持均一且稳定的非互易跳跃能 $\gamma \in [0.5, 1.2]t$，其工艺要求达到了当前凝聚态实验的最前沿，短期内面临极大的工艺瓶颈。

5. 其他必要的补充：非Hermitian BCS平均场与拓扑不变量

5.1 非Hermitian BCS平均场理论推导

为了从物理图像上直观理解为何双占率和配对易受性会呈现 dome 结构，可以使用非对易BCS平均场理论进行定性推导。

在非互易跃迁作用下，单粒子能带结构被重构为复数能谱。考虑有效色散关系 $E(\mathbf{k})$ 引入虚部（由非互易性跃迁产生）：

在 BCS 框架下，自旋单态配对的自洽能隙方程为：

引入非对易Bogoliubov-de Gennes（BdG）哈密顿量：

求解该 BdG 哈密顿量，得到的准粒子激发谱为：

从这个解析表达式可以看出，非互易项 $\gamma \sin(k_x)$ 作为一个纯虚数偏置，直接对应于光谱在虚轴上的拉伸。在开边界条件下，非Bloch band collapse 使得单粒子光谱被极大地压缩并坍缩到实轴，这在热力学和态密度层面导致了一个有效的边界态密度放大：

在 $\gamma < \gamma^*$ 时，随着 $\gamma$ 增加，实部能带扁平化，$\rho_{eff}(\gamma)$ 陡然上升，直接导致超导转变温度 $T_c(\gamma) \sim T_{c0} e^{-1/[V \rho_{eff}(\gamma)]}$ 呈非线性指数上升。而当 $\gamma > \gamma^*$ 时，系统能谱实部发生坍缩并引发过局域化，使得有效费米面处的实际可用态密度骤降，从而使得能隙 $\Delta$ 和转变温度 $T_c$ 迅速下降。这在理论上优雅地解释了 dome 形配对响应的物理起源（参见论文 Fig. 3(b)）。

5.2 单粒子拓扑不变量：Generalized Brillouin Zone 上的非零 Winding Number

非Hermitian皮肤效应（NHSE）并非平庸的耗散，其本质受到单粒子能带拓扑的严格保护。该不变量不能在常规实数动量空间（Brillouin Zone, BZ）中定义，而必须在**广义布里渊区（Generalized Brillouin Zone, GBZ）**上定义。

在周期性边界条件（PBC）下，一维非对易跃迁的 Bloch 哈密顿量为：

当 $k$ 从 $0$ 扫到 $2\pi$ 时，复能量 $E(k)$ 在复平面内扫出一个椭圆轨迹。只要 $\gamma \neq 0$，该轨迹就一定会绕过复平面的原点（或参考能量点 $E_0$），其单粒子拓扑非平凡性由 Winding Number（绕数 $\nu$） 表征：

在该莫尔三角系统中，横向（$\hat{x}$）键上的非互易性赋予了系统非零的绕数 $\nu = 1$（如 Fig. 10(a) 所示）。

根据非Hermitian体-边对应关系（Bulk-Boundary Correspondence），只要 PBC 下复能谱的绕数 $\nu \neq 0$，那么在切断边界（改为 OBC）时，复能谱轨道就必然会全部塌缩至实轴内，并伴随产生向边界堆积的皮肤模。

因此，三角莫尔格点最右侧边界的电荷海量堆积和局域态密度爆发，是受到广义布里渊区非平凡拓扑绕数 $\nu = 1$ 严格物理保护的。这使得该增强机制具有对抗弱无序（Disorder）、格点缺陷及外界非相干噪声的强大拓扑鲁棒性，也是该工作的核心学术亮点。

5.3 对未来量子化学/计算材料研究的启示

本工作对于广大量子化学与多体计算材料领域的科研人员具有重要的借鉴意义。目前，大多数从头算（Ab-initio）量子化学和固体物性计算，依然基于 Hermiticity 假设，重点开发密度泛函理论（DFT）或后Hartree-Fock关联算法。然而：

1. 非互易分子激子传输：在光合作用有机捕光复合物及非共轭高分子体系中，分子间激子能量传输也存在天然的非互易非对称行为。引入非Hermitian多体对角化方法，能够帮助我们设计出具有超高效、定向能量传输能力的量子分子器件。
2. 分子晶体与开放量子体系计算：对于外接电极的分子导线、分子开关，由于电荷的动态注入和抽出，体系是一个天生的不平衡开放多体系统。使用本工作所展示的非Hermitian ED、NH-DMRG及分子级 Lindblad 模拟方法，能够极大提高电导率、自旋选择性传输以及局域激子配对关联的计算精度。