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

Basic Information

• Paper ID: 2509.11462
• Title: Quantum hierarchical Fokker-Planck equations with U(1) gauge fields: Application to the Aharonov-Bohm ring
• Authors: Hyeonseok YANG, Shoki KOYANAGI, Yoshitaka TANIMURA (Department of Chemistry, Kyoto University)
• Classification: quant-ph cond-mat.mes-hall cond-mat.stat-mech
• Publication Date: Updated September 16, 2025
• Paper Link: https://arxiv.org/abs/2509.11462

Abstract

This paper investigates the coupling of a three-dimensional subsystem under time-dependent U(1) gauge fields with a rotationally invariant environment. To capture the subsystem dynamics under thermal excitation and dissipation, the thermal bath must be treated in a non-Markovian and non-perturbative manner. The authors derive the hierarchical equations of motion (HEOM) incorporating gauge fields and obtain quantum hierarchical Fokker-Planck equations with U(1) gauge fields U(1)-QHFPE through Wigner representation transformation. These equations introduce vector fields in the damping operator while maintaining gauge invariance and rotational symmetry. As a practical application, the paper studies the effects of thermal baths in the Aharonov-Bohm (AB) ring, including equilibrium distributions, linear absorption spectra, and simulations of AB currents under thermal conditions.

Research Background and Motivation

Problem Background

1. Research needs for macroscopic quantum phenomena: Macroscopic quantum phenomena such as Bose-Einstein condensation and macroscopic tunneling effects are central to quantum physics, while decoherence and relaxation determine the lifetime of quantum properties.
2. Theoretical challenges in SQUID systems: The study of superconducting quantum interference devices (SQUIDs) is important for quantum superconducting technology and nanodevices, but existing theoretical research is mostly limited to equilibrium conditions, lacking dynamical analysis.
3. Limitations of the Caldeira-Leggett model:
   • The CL model lacks rotational symmetry, leading to inconsistencies in thermal bath coordinates at θ and θ+2π
   • Cannot reproduce the rotational spectral bands of damped free rotors
   • Requires the introduction of the rotationally invariant system-bath (RISB) model

Research Motivation

Quantum noise is constrained by the uncertainty principle, with a relationship between the noise correlation time and the magnitude of energy fluctuations. Therefore, system-bath interactions must be treated within non-perturbative and non-Markovian frameworks, accounting for quantum entanglement between the system and bath.

Core Contributions

1. Established U(1)-HEOM theoretical framework: Extended the hierarchical equations of motion (HEOM) to include U(1) gauge fields while maintaining gauge symmetry.
2. Derived U(1)-QHFPE equations: Converted U(1)-HEOM to quantum hierarchical Fokker-Planck equations in phase space through discretized Wigner transformation.
3. Maintained key symmetries: The proposed equations simultaneously preserve U(1) gauge invariance and rotational symmetry.
4. Complete analysis of AB ring systems: First study of persistent currents in AB rings under non-Markovian conditions, demonstrating the existence of persistent currents at low temperatures under non-Markovian conditions.
5. Verified advantages of RISB model: Through comparison of RISB and CL models, demonstrated the superiority of the RISB model in describing nanoscale ring systems.

Methodology Details

Problem Definition

Study the dynamics of open quantum systems with charged particles in U(1) gauge fields coupled to three-dimensional rotationally invariant thermal baths, with particular focus on equilibrium distributions, linear response spectra, and persistent currents in AB ring systems.

Model Architecture

1. Total Hamiltonian

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

where the system Hamiltonian is:

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

The RISB model bath Hamiltonian:

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 Equations

∂ρ̂_{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 Equations

Through discretized Wigner transformation (DWT-PBC):

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

The U(1)-QHFPE is obtained:

∂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)

Technical Innovations

1. Gauge invariance preservation: U(1)-HEOM remains invariant under gauge transformation A'(x;t) = A(x;t) + ∇χ(x,t).
2. Rotational symmetry: The RISB model ensures rotational symmetry of the system, overcoming deficiencies of the CL model.
3. Non-Markovian treatment: Padé frequency decomposition is used to handle the Drude spectral density function, avoiding high-temperature limitations of Markovian approximations.
4. Discretized Wigner transformation: Introduction of discretized Wigner transformation with periodic boundary conditions, applicable to nanoscale systems.

Experimental Setup

AB Ring System Parameters

• System parameters: ℏ = 1, q = -1, m_S = 0.5, r_0 = 1
• Characteristic frequency: ω_0 = ℏ/(2I_S) = 1.0
• Magnetic flux: Expressed in units of flux quantum Φ_0 = h/q

Numerical Computation Setup

• Grid size:
  • RISB model: |p_n|/ℏ ≤ 31/2, N_θ = 64
  • CL model: N_p = 128, dp = 0.25, N_θ = 64
• Hierarchy truncation: N = 2 (weak coupling), N = 8 (strong coupling)
• Padé frequency terms: K = 1 (high temperature), K = 4 (low temperature)
• Time integration: Runge-Kutta-Fehlberg method

Comparison Methods

1. RISB model: System-bath model preserving rotational symmetry
2. CL model: Traditional Caldeira-Leggett model
3. Markovian vs. non-Markovian: Comparison of results under different approximations

Experimental Results

Equilibrium Distribution

In the Markovian limit, both models display Gaussian distributions centered at p = qr_0A:

• RISB model: Shows discretized energy levels due to DWT-PBC
• CL model: Continuous distribution
• Due to radial-angular coupling, the momentum distribution from RISB results is narrower

Linear Absorption Spectra

Without magnetic field (Φ̄ = 0):

• RISB model: Shows rotational bands reflecting discrete angular momentum states
• CL model: Only shows broad peaks due to lack of rotational symmetry

With magnetic field:

• Each rotational band splits into two peaks
• Split distance initially increases then decreases with Φ̄
• Split vanishes at Φ̄ = 0.5

Transition energy formula: ΔE_{±|n|→±(|n|+1)} = ℏω_0(2|n| + 1 ∓ 2Φ̄)

Persistent Current (Non-Markovian)

Using the Byers-Yang relation: I_n(Φ) = -∂E_n(Φ)/∂Φ = (2ℏω_0/Φ_0)(n - Φ̄)

Key findings:

1. Low temperature conditions (β = 2.5): Significant persistent currents observed
2. High temperature conditions (β = 0.2): Persistent currents completely suppressed
3. Coupling strength effects: Strong coupling leads to reduced net current
4. Periodicity: Current exhibits Φ_0 periodicity

Experimental Findings

1. Limitations of Markovian approximation: High-temperature approximation cannot describe important quantum effects such as persistent currents
2. Advantages of RISB model: More accurately describes quantum coherence effects in nanoscale systems
3. Critical role of temperature: Low temperature is necessary for observing persistent currents

Related Work

Main Research Directions

1. Open quantum system theory: HEOM, QUAPI, pseudomode methods
2. Mesoscopic physics: AB effect, persistent currents, quantum coherence
3. Gauge field theory: Cavity QED, ultrastrong coupling systems

Contributions of This Work

Compared to existing work, this paper is the first to:

• Extend HEOM to include U(1) gauge fields
• Study persistent currents in AB rings under non-Markovian conditions
• Systematically compare RISB and CL model performance in gauge field systems

Conclusions and Discussion

Main Conclusions

1. Successfully established a U(1)-HEOM/QHFPE theoretical framework maintaining gauge invariance and rotational symmetry
2. Demonstrated the existence of persistent currents in AB rings at low temperatures under non-Markovian conditions
3. Verified the superiority of the RISB model in describing nanoscale quantum systems

Limitations

1. Computational complexity: Non-Markovian calculations require substantial computational resources
2. Model applicability: RISB model applies to sub-coherence-length nanorings, CL model to larger rings
3. Temperature constraints: Sufficiently low temperatures are required to observe quantum effects

Future Directions

1. Extended applications: Quantum Hall effect, Berry phase, topological phase systems
2. Arbitrary time-varying external fields: Study AB effects under time-dependent external fields
3. Quantum computing implementation: Utilize quantum computers to handle more complex realistic scenarios

In-Depth Evaluation

Strengths

1. Theoretical rigor: Complete derivation of equation systems maintaining symmetries
2. Methodological innovation: First extension of HEOM to gauge field systems
3. Experimental comprehensiveness: Complete analysis from equilibrium to dynamics
4. Physical insights: Reveals effects of temperature and coupling strength on quantum coherence

Limitations

1. Computational cost: High complexity of non-Markovian calculations limits applicability
2. Model assumptions: Specific Drude spectral density assumptions may limit universality
3. Experimental validation: Lacks direct comparison with actual experimental data

Impact

1. Theoretical contribution: Provides new tools for studying open quantum systems in gauge fields
2. Practical value: Offers guidance for quantum device design and nanoelectronics
3. Reproducibility: Authors commit to providing computational code, enhancing reproducibility

Applicable Scenarios

1. Nanoelectronics: Quantum dots, nanorings, superconducting devices
2. Quantum materials: Topological insulators, quantum Hall systems
3. Quantum information: Quantum coherence and decoherence research

References

The paper cites 106 relevant references covering multiple fields including open quantum system theory, mesoscopic physics, and AB effects, providing a solid theoretical foundation for this research.