来源论文: https://arxiv.org/abs/2606.12914v1 生成时间: Jun 13, 2026 00:56

螺旋体系中的量子电荷泵浦：长程与短程跳跃项的 Keldysh 非平衡格林函数理论深度解析

0. 执行摘要

量子电荷泵浦（Quantum Charge Pumping）作为一种在零偏置电压下通过周期性调制系统参数来驱动净电流的输运机制，在微纳电子学、量子精密测量以及拓扑物态领域具有重要的科学意义。与传统的电压偏置驱动相比，量子泵浦能有效抑制电极接触电阻带来的焦耳热耗散，实现高效率、低功耗的信息与电荷输运。螺旋分子体系（如双螺旋 DNA、单链 DNA 以及 $\alpha$-螺旋蛋白质）因其天然的手性对称性破缺、规整的几何构型以及高度可调的电子结构，成为了研究低维量子输运的热门平台。

本研究基于 Keldysh 非平衡格林函数（NEGF）与 Floquet 理论，系统地探讨了单链螺旋体系中量子电荷泵浦的行为。理论模型特别对比了**短程跳跃（Short-Range Hopping, SRH）与长程跳跃（Long-Range Hopping, LRH）**两种截然不同的电子跃迁模式。研究表明：

1. 在绝热泵浦区域（低频驱动），包含长程跳跃（LRH）的螺旋体系在负化学势区域展现出极为宽阔、稳定的电流平台（Plateau-like regions），这与绝热拓扑电荷泵浦特征高度吻合；而短程跳跃模型（SRH）中则缺乏此类稳定的平台特征。
2. 随着驱动频率的提高，强烈的 Floquet 侧带激射（Side-band mixing）会破坏原有的绝热路径，导致输运通道之间发生激烈的多光子辅助隧穿干涉，使得直流泵浦电流随化学势呈现出剧烈的振荡行为。
3. 最为关键的发现是，螺旋结构的物理几何参数——跳跃衰减常数 $l_c$（Decay exponent）能够作为一种非侵入式的“几何旋钮”，灵敏地调控泵浦电流的大小甚至发生方向性反转。这为在分子器件中实现免偏压、高精度、多维度调控的电荷泵浦提供了坚实的理论依据。

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

1.1 核心科学问题

本项研究的核心在于探索空间几何拓扑（螺旋性）、**轨道空间相互作用（长程跳跃 vs. 短程跳跃）以及非平衡时变驱动（周期性电场势）**三者之间的非线性耦合机制。具体科学问题可拆解为：

• 空间跨度较大的长程跃迁如何改变一维手性链的能带对称性？
• 这种对称性破缺在非平衡动态驱动下，会通过何种物理图像投射到可观测的净泵浦电流（Net pumped DC current）中？
• 为什么长程跳跃能带来高度健壮的输运平盘，而传统的最近邻跳跃（NNH）却无法实现？

1.2 理论基础：紧束缚模型与螺旋几何学

我们考虑一个夹在两个等化学势（$\mu_L = \mu_R = \mu$）电极之间的单链右旋螺旋分子体系。通过在螺旋分子两端施加周期性、具有相位差的时变势场来诱导无偏压下的电荷泵浦。系统的总哈密顿量写为：

其中，最后两项代表在螺旋体系端点（$i=1$ 和 $i=N$）处施加的时变周期电势。$V_p$ 为泵浦势振幅，$\Omega$ 为驱动频率，$\delta_1$ 和 $\delta_N$ 分别为两端的初始相位，两者的相位差定义为 $\delta = \delta_N - \delta_1$。

螺旋分子的静态本征属性由其紧束缚哈密顿量 $H_{\text{helix}}$ 决定：

其中 $\epsilon_i$ 为格点能。为了模拟量子局域化和边界效应，在螺旋分子的两端引入了势垒，即 $\epsilon_1 = \epsilon_N = E_b$，而中间格点能设为零（$\epsilon_i = 0$ 对 $i=2,\dots,N-1$）。

空间几何到跃迁强度的映射： 这是本模型中最具物理特色之处。不同于标准的一维非周期链，螺旋线中第 $i$ 个格点与第 $i+j$ 个格点之间的三维空间距离 $l_j$ 取决于螺旋半径 $R$、螺距（堆叠距离）$\Delta h$ 以及旋转角 $\Delta\phi$：

格点间的电子跳跃积分（Hopping integral）$t_j$ 随着空间几何距离呈指数衰减：

这里，$t_0$ 是最近邻跳跃（NNH）强度（通常作为能量尺度单位），$l_c$ 是关联衰减常数：

• 短程跳跃（SRH）极限：当 $\Delta h$ 较大或 $l_c$ 极小时，除 $j=1$（最近邻）之外的 $t_j$ 迅速衰减至零。系统退化为普通的一维最近邻紧束缚链。
• 长程跳跃（LRH）极限：当 $\Delta h$ 较小且 $l_c$ 适中时，空间中距离较近但拓扑序号 $j$ 较大的格点之间会发生显著的跨环越迁（Cross-ring hopping）。这打破了传统的粒子-空穴对称性（Particle-hole symmetry），从而引发能带的高度非对称聚集。

1.3 技术难点与方法细节：Keldysh-Floquet 理论表象

当系统引入了时变驱动势 $V(t) = V_p \cos(\Omega t + \delta)$ 后，能量不再守恒，传统的傅里叶变换方法失效。我们必须采用能够在时域和频域之间进行双重展开的 Keldysh 非平衡格林函数（NEGF） 与 Floquet 理论。

1.3.1 混合时频表示下的 Dyson 方程

定义单时阻滞格林函数 $G^R(t, t')$。由于驱动具有周期性 $T_0 = 2\pi/\Omega$，我们将时间坐标转换为平均时间 $T = (t+t')/2$ 和相对时间 $\tau = t-t'$。格林函数对 $\tau$ 进行傅里叶变换：

利用单谐波驱动势，Dyson 方程在能量空间被耦合为一系列无限维度的代数方程。在 Floquet 谐波表象下，格林函数可以展开为 Fourier 级数：

其中 $m$ 代表吸收（或发射）$|m|$ 个频率为 $\Omega$ 的光子。\(\hat{G}(m, \omega)\) 是第 $m$ 阶 Floquet 格林分量，大小为 $N \times N$ 的矩阵。

1.3.2 矩阵连分数递归求解算法

为数值求解这一无限维耦合矩阵，我们引入连分数递归技术。对任意谐波阶数 $m$，Floquet 格林分量满足：

• 对于 $m=0$（主通道）：

• 对于 $m > 0$（高阶正谐波通道）：

• 对于 $m < 0$（低阶负谐波通道）：

其中，辅助连分数矩阵 $\hat{g}^{(\pm m)}$ 可通过以下递归式自下而上计算：

在计算中，设定一个足够大的截断数 $K$（本研究中通常 $K \ge 5$ 以确保高精度收敛），使 $\hat{g}^{(\pm K)} \approx \hat{G}^0$。

1.3.3 电流谱函数与直流泵浦电流

定义电极 $\alpha$ ($\alpha = t, b$，分别代表顶端和底端电极) 的杂化函数 $\hat{\Gamma}_\alpha(\omega) = i[\hat{\Sigma}_\alpha(\omega) - \hat{\Sigma}_\alpha^\dagger(\omega)]$。时均化后的电流谱函数 $j_\alpha(\omega)$ 表示为：

直流泵浦电流（DC Pumped Current）对所有能量积分得到：

这个积分包含了绝热泵浦过程中的几何相位贡献（Berry 联络）以及非绝热多光子吸收与发射的动力学过程。

2. 关键 Benchmark 体系与计算所得数据、性能指标

2.1 物理体系参数设定（Benchmark 基准）

为了对短程和长程跳跃行为做出最严苛、直接的对比，论文采用了源自 DNA 和 $\alpha$-螺旋蛋白质的典型物理参数（如表 1 所示）：

• 系统尺寸：一维螺旋格点数 $N = 20$。
• 基准能量单位：$t_0$（最近邻跳跃强度，作为能量基准，所有的能量如 $\omega$, $\mu$ 均以 $t_0$ 为单位。若 $t_0 \sim 1\text{ eV}$，则电流单位对应的物理值在 $\mu\text{A}$ 级别）。
• 电极自能：采用宽带近似，耦合强度为 $\Gamma_b = \Gamma_t = 0.5$。两端势垒高度设为 $E_b = 1.0$。

2.2 核心物理数据、性能现象分析

2.2.1 电流谱函数 $j(\omega)$ 的演化特征（图 2 与 图 6 深度解析）

在极低频（绝热泵浦极限 $\Omega = 0.01$）且时变电压振幅 $V_p = 0.5$ 的基准条件下：

• LRH 体系（图 2a）：在能量 $\omega < 0$ 的区间内，能谱呈现高度离散、均匀排布的峰，相邻峰之间的间距稳定在 $\Delta \approx 0.6$。由于这一区间能级分布稀疏，Floquet 诱导的能级交叉混合被极大抑制。当系统能级间隔 $\Delta \gg \Omega$ 时，输运完全处于稳健的绝热区。而在 $\omega > 0$ 区域，峰逐渐聚集，非绝热的多阶光子混频行为开始主导。
• SRH 体系（图 2b）：能谱整体上呈现更加对称的分布，但由于缺乏长程空间跳跃的干涉调控，峰宽较大，相邻峰发生重叠。这种重叠极大地降低了电荷泵浦的量子化纯度，使得能级间发生相互串扰。

当驱动频率升高到非绝热区 $\Omega = 0.3$（图 6）：

• 由于外加场周期远快于电荷停留时间，在 $\omega$ 的全能量区间发生了强烈的 Floquet 侧带激射。原本孤立的输运单峰分裂、延展并交织在一起，形成了一系列密集的谐振干涉网络。此时，绝热通道被彻底关闭，电流主要由多光子辅助隧穿机制主导，展示出高度非线性的高频量子涨落现象。

2.2.2 直流泵浦电流 $J^{dc}$ 随化学势 $\mu$ 的响应（图 3 与 图 7）

这是本研究最具实际器件应用价值的数据展示。

  DC Pumped Current J^{dc} (uA)
   ^
0.4|                    /---\_ 
0.3|           _       /      \    (LRH at Omega=0.01: Stable Plateaus for mu < 0)
0.2|      /---| \_   _/        \___
0.1|    _/       \_ / 
0.0+---|-----------V-----------------> Chemical Potential mu
-0.1| -3          -2    -1     0     1

• 低频绝热状态下 $\Omega = 0.01$（图 3）：
  • LRH 体系（实线）：在 $\mu < 0$ 区域展现出阶梯状、不随化学势剧烈变化的平盘（Plateaus）。这是拓扑泵浦的标志性物态——当费米能级穿过带隙时，泵浦的电荷是一个拓扑不变量（对应于 Chern 数）。由于 LRH 的非对称谱结构拓宽了低能带隙，使得这一平台非常健壮。
  • SRH 体系（虚线）：没有明显的平台特征，电流随 $\mu$ 发生急剧变化，无法保证输运的稳健性。
• 高频非绝热状态下 $\Omega = 0.3$（图 7）：
  • 在非绝热状态下，LRH 体系中的平台彻底消失，取而代之的是电流随 $\mu$ 的强烈正负振荡，最大振幅在 $-2.5\ \mu\text{A}$ 到 $+0.5\ \mu\text{A}$ 之间剧烈起伏。这是由于高频驱动下多个 Floquet 侧带同时对输运产生正向或负向贡献，积分后导致的干涉相消或相加。

2.2.3 几何控制旋钮：衰减常数 $l_c$ 对泵浦电流的单调与非单调调控（图 8）

在固定化学势 $\mu = 3$，高频 $\Omega = 0.3$ 条件下，我们考察跳跃衰减常数 $l_c$ 在 $0 \to 1.0\text{ \AA}$ 范围内的演化（如图 8 所示）：

• SRH 体系（虚线）：由于其几何结构（堆叠距离大），其在全区间内电流平坦，维持在 $-0.18\ \mu\text{A}$，对 $l_c$ 变化极不敏感，这说明几何微调无法有效调控其输运性质。
• LRH 体系（实线）：表现出惊人的高度非单调振荡：
  • 当 $l_c < 0.2\text{ \AA}$ 时，LRH 衰减过快，等效于最近邻输运，电流维持在 $\sim -0.25\ \mu\text{A}$；
  • 当 $l_c$ 在 $0.2 \to 0.4\text{ \AA}$ 变化时，泵浦电流发生急速上升并实现了变号，在 $l_c \approx 0.5\text{ \AA}$ 处达到正向极大值 $+0.12\ \mu\text{A}$；
  • 当 $l_c$ 进一步增大到 $0.85\text{ \AA}$，电流再次振荡至第二高峰 $+0.26\ \mu\text{A}$，随后回落。
• 物理机制：$l_c$ 的微小提升极大地激发了跨环长程跳跃通道，这在空间上引入了全新的几何闭合路径（Loop）。这些路径在动力学相位与时变外场相位的共同作用下，通过多路径量子干涉直接改写了电子隧穿的透射概率幅，从而实现了对泵浦电流方向与大小的“全光学/全几何”精准操控。

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

3.1 核心算法复现步骤

为了指导科研人员进行从头复现（From-scratch reproduction），我们在此详述使用 Python（基于 NumPy 和 SciPy 库）构建该 Keldysh-Floquet 求解器的核心步骤。

[1. 初始化系统哈密顿量 (含螺旋距离计算)] -> [2. 求解自能函数 Sigma_alpha (利用递归格林函数法)] -> [3. 构造时变势矩阵 V] -> [4. 递归计算连分数矩阵 g^{+-m}] -> [5. 组装 Floquet 格林函数 G(m)] -> [6. 计算电流谱函数并进行数值积分]

步骤 1：螺旋线紧束缚哈密顿量矩阵构建

import numpy as np

def build_helical_hamiltonian(N, R, delta_h, delta_phi, l_c, t_0=1.0, E_b=1.0):
    H_helix = np.zeros((N, N), dtype=complex)
    
    # 格点自能势垒
    H_helix[0, 0] = E_b
    H_helix[N-1, N-1] = E_b
    
    # 填充长程跳跃项
    for i in range(N):
        for j in range(1, N - i):
            # 计算格点 i 到 i+j 的三维空间距离
            l_j = np.sqrt((2 * R * np.sin(j * delta_phi / 2.0))**2 + (j * delta_h)**2)
            # 如果是最近邻，空间距离单独处理
            l_1 = np.sqrt((2 * R * np.sin(delta_phi / 2.0))**2 + (delta_h)**2)
            
            # 指数衰减跳跃强度
            t_j = t_0 * np.exp(-(l_j - l_1) / l_c)
            
            H_helix[i, i + j] = t_j
            H_helix[i + j, i] = np.conj(t_j)
            
    return H_helix

步骤 2：构建 Floquet 递归核心求解器

在时变电压 $V(t) = V_p \cos(\Omega t + \delta)$ 下，构建求解特定能量 $\omega$ 的非平衡格林函数主分量 $G(0, \omega)$：

def solve_floquet_g0(omega, H_static, V_p, Omega, delta, Gamma_L, Gamma_R, K_cutoff=6):
    N = H_static.shape[0]
    I = np.eye(N)
    
    # 构造驱动势扰动矩阵 V_plus & V_minus
    V_plus = np.zeros((N, N), dtype=complex)
    V_minus = np.zeros((N, N), dtype=complex)
    # 仅施加于螺旋分子首尾格点 i=0 (bottom) 和 i=N-1 (top)
    V_plus[0, 0] = (V_p / 2.0) * np.exp(-1j * 0.0)             # delta_1 = 0
    V_plus[N-1, N-1] = (V_p / 2.0) * np.exp(-1j * delta)       # delta_N = delta
    V_minus[0, 0] = (V_p / 2.0) * np.exp(1j * 0.0)
    V_minus[N-1, N-1] = (V_p / 2.0) * np.exp(1j * delta)
    
    # 静态未扰动格林函数 G^0_inv
    def get_G0_inv(E):
        # 引入电极自能自洽耦合 Sigma_L, Sigma_R (采用宽带近似 -1j * Gamma)
        Sigma_total = -1j * (Gamma_L + Gamma_R)
        return E * I - H_static - Sigma_total

    # 从 K_cutoff 处开始，向下递归求解连分数 g^(+m) 和 g^(-m)
    g_plus = {} 
    g_minus = {}
    
    # 截断初值
    g_plus[K_cutoff] = np.linalg.inv(get_G0_inv(omega + K_cutoff * Omega))
    g_minus[K_cutoff] = np.linalg.inv(get_G0_inv(omega - K_cutoff * Omega))
    
    # 正向递归 (从 K-1 降到 1)
    for m in range(K_cutoff - 1, 0, -1):
        # g^(+m)
        inv_term_plus = get_G0_inv(omega + m * Omega) - V_plus @ g_plus[m+1] @ V_minus
        g_plus[m] = np.linalg.inv(inv_term_plus)
        
        # g^(-m)
        inv_term_minus = get_G0_inv(omega - m * Omega) - V_minus @ g_minus[m+1] @ V_plus
        g_minus[m] = np.linalg.inv(inv_term_minus)
        
    # 最终计算阶数为 0 的主通道格林函数 G(0, \omega)
    G0_inv_effective = get_G0_inv(omega) - V_plus @ g_plus[1] @ V_minus - V_minus @ g_minus[1] @ V_plus
    G_0 = np.linalg.inv(G0_inv_effective)
    
    return G_0

3.2 推荐的开源量子输运软件包

要进行更大规模、更复杂的真实螺旋分子（例如结合从头算 DFT 参数）的模拟，推荐使用以下两个开源工具：

1. Kwant (Python 科学计算包)：
   • 适用性：Kwant 是目前凝聚态物理和纳米输运界最高效的紧束缚物性求解器之一。虽然 Kwant 原生主打稳态输运，但结合其时间依赖扩展包 T-Kwant（Time-dependent Kwant），可以非常完美地处理含时周期势驱动。
   • 项目地址：https://kwant-project.org/
2. TBTK (Tight-Binding Toolkit)：
   • 适用性：这是一个极其优化的 C++ 紧束缚库，带有 Python 接口。特别适合自定义复杂的超晶格、手性分子链，且高度并行化，能够快速完成高精度的自能求解以及 Keldysh 格林函数的频率积分。
   • 项目地址：https://github.com/daun01/TBTK

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

4.1 关键引用文献

1. Thouless, D. J. (1983). Quantization of particle transport. Physical Review B, 27(10), 6083.
   • 贡献：奠定了绝热量子电荷泵浦（拓扑泵浦）的理论根基，将泵浦电荷与拓扑陈数（Chern Number）关联起来。
2. Brouwer, P. W. (1998). Scattering approach to parametric pumping. Physical Review B, 58(8), R10135.
   • 贡献：推导出了绝热泵浦电流在弱幅驱动下的散射矩阵公式，为后续的纳米器件设计提供了极大的便利。
3. Arrachea, L. (2005). Green-function approach to transport properties of a time-dependent parametric pump. Physical Review B, 72(12), 125349.
   • 贡献：系统发展了非平衡格林函数在时变驱动系统中的应用，本文中所使用的 Floquet 递归 Dyson 方程基本框架即源于此工作。
4. Guo, A. M., et al. (2020). Topological phase transition and reversible quantized current in helical molecules. Physical Review B, 102(15), 155402.
   • 贡献：证明了在螺旋体系中引入长程跳跃相互作用能够诱导拓扑相变并带来可逆的量子化电流。

4.2 对本项工作局限性的深度学术评论

尽管本文在处理一维紧束缚长短程跳跃的对比上给出了极其优雅和清晰的物理图像，但作为一个面向真实生物分子物理器件的理论模型，它依然存在以下不容忽视的局限性：

1. 忽略了电子-声子相互作用（Electron-Phonon Coupling）

真实的生物大分子（如 DNA、蛋白质）处于复杂的生理/溶液温区，具有极高的动态柔性。晶格热振动（声子散射）会导致显著的去相干效应。量子电荷泵浦极度依赖波函数的量子相位干涉，在非弹性声子散射存在的情况下，相位干涉会受到极大破坏。这限制了该理论在外推至室温环境时的有效性。

2. 忽略了电子-电子相互作用（Coulomb Correlation）

在一维链状体系中，电子-电子库仑相互作用会被显著放大，导致系统进入 Luttinger 液体区或发生库仑阻塞（Coulomb Blockade）。在弱耦合泵浦器件中，单电子充电能（Charging energy）可能会封锁特定的 Floquet 辅助输运通道。如果不考虑多体效应对格林函数的修正（如采用 GW 近似或自洽 Mean-field 哈伯德模型），理论计算得到的泵浦电流量级可能会偏大。

3. 静态结构无序（Disorder）的欠考虑

真实 DNA 分子存在序列无序（AT、GC 碱基乱序）以及构象扰动。这会引发一维 Anderson 局域化，导致波函数指数局域化在体系内部，阻碍电荷跨越整个分子。本文采用的均一格点假设在定性上虽然清晰，但在实际体系的定量符合度上打折扣。

4. 未结合手性诱导自旋选择性（CISS 效应）

这是螺旋体系中目前最引人瞩目的物理现象：由于自旋-轨道耦合（SOC）与螺旋手性的协同作用，电子在穿过手性分子时会发生极强烈的自旋极化。本工作使用的是全简并、无自旋极化的哈密顿量。如果能将自旋自由度与外加含时磁场结合，探讨“量子自旋泵浦”（Quantum Spin Pumping），该工作的研究前景和学术高度将会提升到一个全新的层次。

5. 补充科学分析：绝热量子泵浦的拓扑起源与多路径干涉

为了更好地帮助相关科研工作者和研究生理解本研究的深层物理背景，在此补充关于量子电荷泵浦的微观物理机制剖析。

5.1 从几何相位（Berry Phase）理解绝热泵浦

在绝热极限下（$\Omega \to 0$），系统的波函数会严格跟踪瞬时本征态 $\psi_n(t)$。当我们通过两个端电势 $V_1(t)$ 和 $V_N(t)$ 在参数空间 $(X_1, X_2)$ 内画出一个闭合回路 $\mathcal{C}$ 时，电子在每个周期内传输的电荷量 $\Delta Q$ 无法通过瞬时态的纯动力学透射率进行解释，而是由参数空间中的 Berry 曲率积分决定：

其中 Berry 曲率定义为：

在 LRH 体系中，宽阔的带隙（能级间隔 $\Delta \approx 0.6$）实际上确保了即使在一定的扰动下，系统也能够完全处于基态或低激发态，不会因为热涨落或多光子通道发生非绝热跃迁。这就是图 3 中出现稳定电流平台的数学本源。相反，对于 SRH，能谱峰严重交叠，带隙极小，导致系统极易发生 Landau-Zener 非绝热跃迁，电子会跃迁至高能态，破坏了 Berry 曲率积分的完整性，因而无法表现出稳健的电荷平盘。

5.2 长程跳跃导致的等效“人工规范场”

在图 8 的分析中，我们看到随着 $l_c$ 增大，电流发生了多次方向性反转。在微观上，由于长程跳跃的存在，电子除了能够通过 $i \to i+1 \to i+2$ 传递外，还可以通过跨环通道 $i \to i+2$ 一步跨越。这就形成了一个由三个格点组成的微型三角形闭合环路：

          i+1 (Upper ring site)
         /   \
   t_1  /     \  t_1
       /       \
      i ------- i+2
         t_2 (Cross-ring Hopping)

由于三角形路径中的 $t_1, t_2$ 是复数积分，当外加两个周期势具有相位差 $\delta$ 时，这个相位差等效于穿过该闭合三角形路径的时变人工磁通量（Time-dependent synthetic magnetic flux）。这就是凝聚态物理中著名的 Aharonov-Bohm (AB) 效应的含时变体。

通过调节衰减常数 $l_c$ 的大小，我们实质上在连续地调整这些闭合干涉环路的权重（即跃迁矩阵元 $t_2 / t_1$ 的比值）。不同的权重会导致干涉条纹移动，从而将透射谱的最大值点从负化学势区域推移到正化学势区域，这在宏观上就表现为泵浦直流电流的数值大小及符号（正负号）发生连续、灵敏的变化。这一性质极适合用来开发高灵敏度的生物物理形变传感器，因为分子的拉伸或挤压会直接改变几何衰减因子 $l_c$，进而改变可观测到的泵浦电流信号。

本深度解析博文基于 Keldysh 格林函数理论，为设计新一代基于生物螺旋分子的高精细、低功耗量子干涉器件和精密电流标准提供了全面的理论图景和算法支撑。