2025-11-15T05:28:11.845417

Quantum hierarchical Fokker-Planck equations with U(1) gauge fields: Application to the Aharonov-Bohm ring

Yang, Koyanagi, Tanimura
We investigate a three-dimensional subsystem under a time-dependent U(1) gauge field coupled to rotationally invariant environments. To capture the dynamic behavior of the subsystem under thermal excitations and dissipations, it is imperative to treat the bath in a non-Markovian and nonperturbative manner. This is because quantum noise is constrained by the uncertainty principle, which dictates the relationship between the noise correlation time and the amplitude of the energy fluctuation. To this end, we derive the hierarchical equations of motion (HEOM) incorporating the gauge field, enabling a rigorous investigation of the dynamics of the reduced subsystem. Transforming the HEOM into the Wigner representation yields quantum hierarchical Fokker-Planck equations [U(1)-QHFPE] with U(1) gauge fields. These equations incorporate vector fields into the damping operators while preserving both gauge invariance and rotational symmetry. To demonstrate the practical use of the formalism, the effects of a heat bath in the Aharonov-Bohm (AB) ring. Our investigation includes simulations of the equilibrium distribution, linear absorption spectra, and AB currents under thermal conditions. Within a rotationally invariant system-bath (RISB) model, we predict the emergence of a persistent current even in dissipative environments, provided the bath is non-Markovian and the temperature is sufficiently low. We also assessed the validity of the Caldeira-Leggett model in this context.
academic

Quantum hierarchical Fokker-Planck equations with U(1) gauge fields: Application to the Aharonov-Bohm ring

基本信息

  • 论文ID: 2509.11462
  • 标题: Quantum hierarchical Fokker-Planck equations with U(1) gauge fields: Application to the Aharonov-Bohm ring
  • 作者: Hyeonseok YANG, Shoki KOYANAGI, Yoshitaka TANIMURA (京都大学化学系)
  • 分类: quant-ph cond-mat.mes-hall cond-mat.stat-mech
  • 发表时间: 2025年9月16日更新
  • 论文链接: https://arxiv.org/abs/2509.11462

摘要

本文研究了三维子系统在时间依赖U(1)规范场作用下与旋转不变环境的耦合问题。为了捕获子系统在热激发和耗散作用下的动力学行为，必须以非马尔可夫和非微扰的方式处理热浴。作者推导了包含规范场的层次运动方程(HEOM)，通过Wigner表示变换得到了带有U(1)规范场的量子层次Fokker-Planck方程U(1)-QHFPE。这些方程在阻尼算子中引入矢量场，同时保持规范不变性和旋转对称性。作为实际应用，文章研究了Aharonov-Bohm(AB)环中热浴的效应，包括平衡分布、线性吸收谱和热条件下的AB电流的模拟。

研究背景与动机

问题背景

  1. 宏观量子现象的研究需求: 宏观量子现象如玻色-爱因斯坦凝聚和宏观隧穿效应是量子物理的核心问题，而退相干和弛豫决定了量子性质的寿命。
  2. SQUID系统的理论挑战: 超导量子干涉器件(SQUID)的研究对量子超导技术和纳米器件具有重要意义，但现有理论研究大多局限于平衡条件，缺乏动力学分析。
  3. Caldeira-Leggett模型的局限性:
    • CL模型缺乏旋转对称性，导致热浴坐标在θ和θ+2π处不一致
    • 无法重现阻尼自由转子的旋转光谱带
    • 需要引入旋转不变系统-浴模型(RISB)

研究动机

量子噪声受不确定性原理约束，噪声相关时间与能量涨落幅度之间存在关系，因此系统-浴相互作用必须在非微扰和非马尔可夫框架下处理，考虑系统与浴之间的量子纠缠。

核心贡献

  1. 建立了U(1)-HEOM理论框架: 将层次运动方程(HEOM)扩展到包含U(1)规范场的情况，保持规范对称性。
  2. 推导了U(1)-QHFPE方程: 通过离散化Wigner变换，将U(1)-HEOM转换为相空间中的量子层次Fokker-Planck方程。
  3. 保持关键对称性: 所提出的方程同时保持了U(1)规范不变性和旋转对称性。
  4. AB环系统的完整分析: 首次在非马尔可夫条件下研究了AB环的持续电流，证明了在低温非马尔可夫条件下持续电流的存在。
  5. 验证了RISB模型的优势: 通过对比RISB和CL模型，证明了RISB模型在描述纳米环系统方面的优越性。

方法详解

任务定义

研究带电粒子在U(1)规范场中与三维旋转不变热浴耦合的开放量子系统动力学，特别关注AB环系统的平衡分布、线性响应谱和持续电流。

模型架构

1. 总哈密顿量

H^RISB_tot(t) = H_S(t) + H^RISB_I+B

其中系统哈密顿量为：

H_S(t) = 1/(2m_S)[p̂ - qA(x̂;t)]² + U_φ(x̂;t)

RISB模型的浴哈密顿量：

H^RISB_I+B = Σ_α Σ_k [(p̂^α_k)²/(2m^α_k) + m^α_k(ω^α_k)²/2 (q̂^α_k - c^α_k V̂_α/(m^α_k(ω^α_k)²))²]

2. U(1)-HEOM方程

∂ρ̂_{n_α}(t)/∂t = -[i/ℏ Ĥ^×_S(t) + Σ_α Σ_j n^α_j ν^α_j] ρ̂_{n_α}(t)
                  + Σ_α Σ_j Φ̂_α ρ̂_{n_α+e^α_j}(t) + Σ_α Σ_j n^α_j Θ̂^α_j ρ̂_{n_α-e^α_j}(t)

3. U(1)-QHFPE方程

通过离散化Wigner变换(DWT-PBC)：

W_{n_α}(p_n, r) = Π_α [1/(2πℏ) ∫^L_α_-L_α dξ_α e^(-ip^α_n ξ_α/ℏ)] ρ_{n_α}(r + ξ/2, r - ξ/2)

得到U(1)-QHFPE：

∂W_{n_α}(p,r)/∂t = -[L̂_qm(t) + Σ_α Σ_j n^α_j ν^α_j] W_{n_α}(p,r)
                   + Σ_α Σ_j Φ̃_α W_{n_α+e^α_j}(p,r) + Σ_α Σ_j n^α_j Θ̃^α_j(t) W_{n_α-e^α_j}(p,r)

技术创新点

  1. 规范不变性保持: 在规范变换A'(x;t) = A(x;t) + ∇χ(x,t)下，U(1)-HEOM保持不变性。
  2. 旋转对称性: 通过RISB模型确保了系统的旋转对称性，克服了CL模型的缺陷。
  3. 非马尔可夫处理: 使用Padé频率分解处理Drude谱密度函数，避免了马尔可夫近似的高温限制。
  4. 离散化Wigner变换: 引入周期边界条件的离散化Wigner变换，适用于纳米尺度系统。

实验设置

AB环系统参数

  • 系统参数：ℏ = 1, q = -1, m_S = 0.5, r_0 = 1
  • 特征频率：ω_0 = ℏ/(2I_S) = 1.0
  • 磁通量：以通量量子Φ_0 = h/q为单位

数值计算设置

  • 网格大小:
    • RISB模型：|p_n|/ℏ ≤ 31/2, N_θ = 64
    • CL模型：N_p = 128, dp = 0.25, N_θ = 64
  • 层次截断: N = 2(弱耦合), N = 8(强耦合)
  • Padé频率数: K = 1(高温), K = 4(低温)
  • 时间积分: Runge-Kutta-Fehlberg方法

对比方法

  1. RISB模型: 保持旋转对称性的系统-浴模型
  2. CL模型: 传统Caldeira-Leggett模型
  3. 马尔可夫vs非马尔可夫: 对比不同近似下的结果

实验结果

平衡分布

在马尔可夫极限下，两种模型都显示以p = qr_0A为中心的高斯分布：

  • RISB模型：由于DWT-PBC显示离散化能级
  • CL模型：连续分布
  • 由于径向-角向耦合，RISB结果的动量分布更窄

线性吸收谱

无磁场时(Φ̄ = 0)：

  • RISB模型：显示反映离散角动量态的旋转带
  • CL模型：由于缺乏旋转对称性仅显示宽峰

有磁场时：

  • 每个旋转带分裂为两个峰
  • 分裂距离随Φ̄增加先增后减
  • 在Φ̄ = 0.5时分裂消失

转移能量公式：ΔE_{±|n|→±(|n|+1)} = ℏω_0(2|n| + 1 ∓ 2Φ̄)

持续电流(非马尔可夫)

使用Byers-Yang关系：I_n(Φ) = -∂E_n(Φ)/∂Φ = (2ℏω_0/Φ_0)(n - Φ̄)

主要发现：

  1. 低温条件下(β = 2.5)：观察到明显的持续电流
  2. 高温条件下(β = 0.2)：持续电流完全被抑制
  3. 耦合强度影响：强耦合导致净电流减少
  4. 周期性：电流具有Φ_0的周期性

实验发现

  1. 马尔可夫近似的局限性: 高温近似无法描述持续电流等重要量子效应
  2. RISB模型的优势: 在纳米尺度系统中更准确地描述量子相干效应
  3. 温度的关键作用: 低温是观察持续电流的必要条件

相关工作

主要研究方向

  1. 开放量子系统理论: HEOM、QUAPI、伪模方法
  2. 介观物理: AB效应、持续电流、量子相干
  3. 规范场理论: 腔QED、超强耦合系统

本文贡献

相比现有工作，本文首次：

  • 将HEOM扩展到包含U(1)规范场
  • 在非马尔可夫条件下研究AB环持续电流
  • 系统比较RISB和CL模型在规范场系统中的表现

结论与讨论

主要结论

  1. 成功建立了保持规范不变性和旋转对称性的U(1)-HEOM/QHFPE理论框架
  2. 证明了在低温非马尔可夫条件下AB环中持续电流的存在
  3. 验证了RISB模型在描述纳米尺度量子系统方面的优越性

局限性

  1. 计算复杂度: 非马尔可夫计算需要大量计算资源
  2. 模型适用范围: RISB模型适用于亚相干尺寸的纳米环，CL模型适用于大尺寸环
  3. 温度限制: 需要足够低的温度才能观察到量子效应

未来方向

  1. 扩展应用: 量子霍尔效应、Berry相位、拓扑相系统
  2. 任意时变外场: 研究时间依赖外场下的AB效应
  3. 量子计算实现: 利用量子计算机处理更复杂的现实场景

深度评价

优点

  1. 理论严谨性: 完整推导了保持对称性的方程体系
  2. 方法创新性: 首次将HEOM扩展到规范场系统
  3. 实验充分性: 从平衡态到动力学的全面分析
  4. 物理洞察: 揭示了温度和耦合强度对量子相干的影响

不足

  1. 计算成本: 非马尔可夫计算的高复杂度限制了应用范围
  2. 模型假设: Drude谱密度的特定假设可能限制普适性
  3. 实验验证: 缺乏与实际实验数据的直接比较

影响力

  1. 理论贡献: 为研究规范场中的开放量子系统提供了新工具
  2. 实用价值: 对量子器件设计和纳米电子学具有指导意义
  3. 可复现性: 作者承诺提供计算代码，增强了可复现性

适用场景

  1. 纳米电子学: 量子点、纳米环、超导器件
  2. 量子材料: 拓扑绝缘体、量子霍尔系统
  3. 量子信息: 量子相干和退相干研究

参考文献

论文引用了106篇相关文献，涵盖了开放量子系统理论、介观物理、AB效应等多个领域的重要工作，为本研究提供了坚实的理论基础。