2025-11-19T16:13:13.264208

Open Quantum Dynamics Theory for Coulomb Potentials: Hierarchical Equations of Motion for Atomic Orbitals (AO-HEOM)

Zhang, Tanimura

We investigate the quantum dynamics of Coulomb potential systems in thermal baths. We study these systems within the framework of open quantum dynamics theory, focusing on preserving the rotational symmetry of the entire system, including the baths. Thus, we employ a three-dimensional rotationally invariant system-bath (3D-RISB) model to derive numerically ``exact'' hierarchical equations of motion for atomic orbitals (AO-HEOM) that enable a non-perturbative and non-Markovian treatment of system-bath interactions at finite temperatures. To assess the formalism, we calculated the linear absorption spectrum of an atomic system under isotropic thermal environment, with systematic variation of system-bath coupling strength and temperature.

academic

クーロンポテンシャルのための開放量子動力学理論：原子軌道の階層運動方程式(AO-HEOM)

基本情報

論文ID: 2510.11981
タイトル: Open Quantum Dynamics Theory for Coulomb Potentials: Hierarchical Equations of Motion for Atomic Orbitals (AO-HEOM)
著者: Yankai Zhang、Yoshitaka Tanimura(京都大学化学科)
分類: quant-ph、cond-mat.stat-mech、physics.chem-ph
発表日時: 2025年10月15日更新
論文リンク: https://arxiv.org/abs/2510.11981

従来手法の限界: 標準的な開放量子動力学手法(Lindblad主方程式、マルコフ-Redfield方程式など)は、回転対称性を有するクーロンポテンシャル系を処理する際に根本的な欠陥を有しており、特に低温および強結合条件下で非物理的な異常挙動を示す。
回転対称性の破れ: 従来のCaldeira-Leggett(CL)モデルは古典的にはランジュバン動力学と等価であるが、量子系では離散的な回転エネルギーバンドを再現できず、「浴エンタングルメント」の破壊につながる。
量子熱ノイズの非マルコフ性: 不確定性原理により、量子熱ノイズは本質的に非マルコフ的であり、その相関時間と振幅は不確定性関係を満たす必要があり、非摂動的処理が必要である。

研究の重要性

基礎物理的意義: クーロンポテンシャルは水素原子からイオン結晶中の色心問題に至るまで、広範な物理系を包含している
実用的応用: 凝縮系物理において、熱励起は重要な役割を果たし、非マルコフ的な揺らぎ処理が必要である
キャビティ量子電動力学: 強結合キャビティQED系に対する非摂動的処理枠組みを提供する

核心的貢献

3D-RISBモデルの提案: 三次元回転不変系-浴モデルを開発し、熱浴を含む全系の回転対称性を保持した
AO-HEOM形式化手法の確立: 原子軌道の階層運動方程式を導出し、数値「厳密」な非摂動的・非マルコフ的処理を実現した
GPU最適化実装: GPU を完全に活用した AO-HEOM 計算コードを開発し、計算効率を大幅に向上させた
系統的検証: 線形吸収スペクトル計算を通じて、異なる温度と結合強度下での理論枠組みの有効性を検証した
理論的洞察: 「浴エンタングルメント」現象が量子系動力学に及ぼす重要な影響を明らかにした

Ĥtot = ĤS + ∑α=x,y,z Ĥα I+B

系の部分はクーロンポテンシャルである：

ĤS = ∑α=x,y,z p̂²α/2me - Zpe²/4πε₀r

三方向の熱浴は以下の通りである：

Ĥα I+B = ∑j [(p̂αj)²/2mαj + ½mαj(ωαj)²(q̂αj - cαjV̂α/mαj(ωαj)²)²]

2. AO-HEOM方程式

Kα-1/Kα Padé近似を採用し、AO-HEOM方程式は以下の通りである：

d/dt ρ̂{nα} = -[i/ℏ ĤS× + ∑α=x,y,z ∑k=0^Kα (nαkναk)] ρ̂{nα}
             - i/ℏ ∑α=x,y,z ∑k=0^Kα nαk Θ̂αk ρ̂{nα-ekα}
             - i/ℏ ∑α=x,y,z ∑k=0^Kα V̂α× ρ̂{nα+ekα}

3. 線形吸収スペクトル計算

Iα'α(ω) = Im(i/ℏ)∫₀^∞ dt e^(iωt) Tr{μ̂α' Ĝ(t)μ̂α× ρ̂eq}

技術的革新点

回転対称性の保持: 3D-RISBモデルを通じて全系(熱浴を含む)の回転不変性を確保する
非摂動的処理: HEOM手法は強い系-浴結合を処理でき、摂動近似の限界を回避する
浴エンタングルメント保持: 階層構造は複数回の系-浴相互作用から生じる量子相関を記述する
球座標表現: 系-浴相互作用は(V̂x, V̂y, V̂z) ≡ (r sin θ cos φ, r sin θ sin φ, r cos φ)として表現される