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

U(1)ゲージ場を含む量子階層Fokker-Planck方程式：Aharonov-Bohm環への応用

基本情報

論文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電流のシミュレーションを研究している。

研究背景と動機

問題背景

巨視的量子現象の研究需要: ボーズ・アインシュタイン凝縮と巨視的トンネル効果などの巨視的量子現象は量子物理学の中心的課題であり、デコヒーレンスと緩和が量子特性の寿命を決定している。
SQUID系の理論的課題: 超伝導量子干渉装置(SQUID)の研究は量子超伝導技術とナノデバイスに重要な意義を持つが、既存の理論研究のほとんどは平衡条件に限定されており、動力学的分析が不足している。
Caldeira-Leggettモデルの限界:
- CLモデルは回転対称性を欠いており、熱浴座標がθとθ+2πで一致しない
- 減衰自由回転子の回転スペクトル帯を再現できない
- 回転不変系-浴モデル(RISB)の導入が必要

研究動機

量子ノイズは不確定性原理によって制約され、ノイズ相関時間とエネルギー変動幅の間に関係が存在するため、系-浴相互作用は非摂動的かつ非マルコフ的枠組みの下で処理され、系と浴間の量子エンタングルメントを考慮する必要がある。

核心的貢献

U(1)-HEOM理論枠組みの確立: 階層運動方程式(HEOM)をU(1)ゲージ場を含む場合に拡張し、ゲージ対称性を保持した。
U(1)-QHFPE方程式の導出: 離散化Wigner変換を通じてU(1)-HEOMを相空間の量子階層Fokker-Planck方程式に変換した。
重要な対称性の保持: 提案された方程式はU(1)ゲージ不変性と回転対称性の両方を同時に保持している。
AB環系の完全な分析: 非マルコフ条件下でのAB環の持続電流を初めて研究し、低温非マルコフ条件下での持続電流の存在を証明した。
RISBモデルの優位性の検証: RISBとCLモデルの比較を通じて、RISBモデルがナノ環系の記述において優れていることを証明した。

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)

技術的革新点

ゲージ不変性の保持: ゲージ変換A'(x;t) = A(x;t) + ∇χ(x,t)の下で、U(1)-HEOMは不変性を保持する。
回転対称性: RISBモデルを通じてシステムの回転対称性を確保し、CLモデルの欠陥を克服した。
非マルコフ処理: Padé周波数分解を使用してDrudeスペクトル密度関数を処理し、マルコフ近似の高温限界を回避した。
離散化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法