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

Basic Information

• Paper ID: 2510.11981
• Title: Open Quantum Dynamics Theory for Coulomb Potentials: Hierarchical Equations of Motion for Atomic Orbitals (AO-HEOM)
• Authors: Yankai Zhang, Yoshitaka Tanimura (Department of Chemistry, Kyoto University)
• Classification: quant-ph, cond-mat.stat-mech, physics.chem-ph
• Publication Date: Updated October 15, 2025
• Paper Link: https://arxiv.org/abs/2510.11981

Abstract

This paper investigates the quantum dynamics of Coulomb potential systems in thermal baths. Within the framework of open quantum dynamics theory, emphasis is placed on maintaining rotational symmetry of the entire system including the thermal bath. Employing the three-dimensional rotationally invariant system-bath (3D-RISB) model, numerically "exact" solutions of the hierarchical equations of motion for atomic orbitals (AO-HEOM) are derived, enabling non-perturbative and non-Markovian treatment of system-bath interactions at finite temperatures. The formalism is systematically evaluated by calculating linear absorption spectra of atomic systems in isotropic thermal environments, with systematic variation of system-bath coupling strength and temperature.

Research Background and Motivation

Core Problems

1. Limitations of Traditional Methods: Standard open quantum dynamics approaches (such as Lindblad master equations and Markovian-Redfield equations) exhibit fundamental defects when treating Coulomb potential systems with rotational symmetry, particularly exhibiting non-physical anomalies under low-temperature and strong-coupling conditions.
2. Rotational Symmetry Breaking: Although the classical Caldeira-Leggett (CL) model is classically equivalent to Langevin dynamics, it fails to reproduce discrete rotational energy bands in quantum systems, leading to destruction of "bath entanglement."
3. Non-Markovian Nature of Quantum Thermal Noise: Due to the uncertainty principle, quantum thermal noise is inherently non-Markovian, with correlation time and amplitude necessarily satisfying uncertainty relations, requiring non-perturbative treatment.

Research Significance

• Fundamental Physics: Coulomb potentials encompass a broad range of physical systems from hydrogen atoms to color center problems in ionic crystals
• Practical Applications: In condensed matter physics, thermal excitations play crucial roles, requiring non-Markovian fluctuation treatment
• Cavity Quantum Electrodynamics: Provides a non-perturbative treatment framework for strongly-coupled cavity QED systems

Core Contributions

1. Proposed 3D-RISB Model: Developed a three-dimensional rotationally invariant system-bath model that maintains rotational symmetry of the entire system including the thermal bath
2. Established AO-HEOM Formalism: Derived hierarchical equations of motion for atomic orbitals, enabling numerically "exact" non-perturbative and non-Markovian treatment
3. GPU-Optimized Implementation: Developed fully GPU-accelerated AO-HEOM computational code, significantly enhancing computational efficiency
4. Systematic Verification: Validated the theoretical framework across different temperatures and coupling strengths through linear absorption spectrum calculations
5. Theoretical Insights: Revealed the critical influence of "bath entanglement" phenomena on quantum system dynamics

Methodology Details

Task Definition

Investigation of quantum dynamics of Coulomb potential systems in thermal bath environments, with inputs being the system Hamiltonian, bath parameters, and coupling strength, and outputs being the time evolution and spectral response of the system.

Model Architecture

1. 3D-RISB Hamiltonian

The total Hamiltonian is expressed as:

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

where the system part is the Coulomb potential:

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

and the thermal baths in three directions are:

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

2. AO-HEOM Equations

Using Kα-1/Kα Padé approximation, the AO-HEOM equations are:

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. Linear Absorption Spectrum Calculation

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

Technical Innovations

1. Rotational Symmetry Preservation: The 3D-RISB model ensures rotational invariance of the entire system including the thermal bath
2. Non-Perturbative Treatment: The HEOM method can handle strong system-bath coupling, avoiding limitations of perturbative approximations
3. Bath Entanglement Preservation: The hierarchical structure describes quantum correlations arising from multiple system-bath interactions
4. Spherical Coordinate Representation: System-bath interactions are represented as (V̂x, V̂y, V̂z) ≡ (r sin θ cos φ, r sin θ sin φ, r cos φ)

Experimental Setup

Computational Parameters

• Number of Basis States: 55 energy eigenstates (principal quantum numbers n=1 to n=5)
• Hierarchy Depth: Nα=2 for each thermal bath, using 0/1 form Padé approximation
• Time Integration: Fourth-order Runge-Kutta method with 3000 time steps
• Isotropic Bath: γx=γy=γz=1, ηx=ηy=ηz=η

Computational Environment

• Hardware: Intel Core i9-13900KF CPU + NVIDIA GeForce RTX 4090 GPU
• Software: Python 3.13.5, CuPy 13.4.1, NumPy 2.2.5, SciPy 1.15.3
• Resource Consumption: Each simulation requires 2.4GB GPU memory and 20,000 seconds computational time

Evaluation Metrics

• Linear absorption spectrum intensity and peak positions
• Comparison with Fermi's golden rule
• Spectral characteristics under different temperatures and coupling strengths

Experimental Results

Main Results

1. Strong Coupling Case (η = 0.001-0.01)

• High Temperature (β=1.0): Broadened single peak appearing mainly in the Lyman region
• Intermediate Temperature (β=2.0): Balmer series peaks begin to emerge
• Low Temperature (β=5.0): Clear Lyman-α (ω≈0.4) and Balmer-α/β (ω≈0.08) peaks

2. Weak Coupling Case (η = 0.0001-0.0005)

• Series peaks are sharper and more distinct
• Brackett series peaks observed even at low temperatures
• Peak positions show good agreement with Fermi's golden rule predictions

Important Findings

1. High-n' Transition Suppression: Under strong thermal bath coupling and high temperature, large-n to n' transitions are suppressed, exhibiting semi-classical behavior
2. Spectral Line Broadening Mechanism: Bath entanglement leads to effective continuum-like atomic energy levels, producing broadened peaks
3. Temperature Dependence: At low temperatures, thermal excited state populations decrease, with low-frequency transitions becoming more prominent

Ablation Studies

Verification of theoretical predictions through variation of maximum principal quantum number n': transitions satisfying βℏωnn'≪1 exhibit semi-classical characteristics under strong coupling and high temperature.

Related Work

Limitations of Traditional Methods

1. Lindblad Master Equation: Produces non-physical equal population anomalies under Ohmic spectral density
2. Quantum Fokker-Planck Equation: Reliable at high temperatures but violates positive-definiteness at low temperatures
3. Time-Convolutionless Redfield Equation: Although capturing non-Markovian noise, it destroys bath entanglement

Related Developments

• QHFPE (Quantum Hierarchical Fokker-Planck Equation)
• MB-HEOM (Multi-Bath HEOM)
• U(1)-HEOM (HEOM with U(1) gauge field)

Conclusions and Discussion

Main Conclusions

1. AO-HEOM successfully handles quantum dynamics of Coulomb potential systems in thermal baths
2. The 3D-RISB model maintains rotational symmetry, avoiding fundamental defects of traditional methods
3. Quantum-to-classical transitions observed under strong coupling validate theoretical predictions

Limitations

1. Computational Cost: The 3D-RISB model is computationally intensive, requiring GPU acceleration
2. Temperature Range: Non-physical high-temperature conditions were employed to highlight thermal effects
3. Single-Electron Restriction: Currently applicable only to single-electron systems

Future Directions

1. Extension to multi-electron systems and Fock space
2. Calculation of higher-order nonlinear spectra, such as 2D spectroscopy
3. Inclusion of gauge field effects for cavity QED research
4. Application to practical systems such as color centers and ionic liquids

In-Depth Evaluation

Strengths

1. Theoretical Rigor: Derived from first principles while maintaining physical symmetries
2. Methodological Innovation: Extension of HEOM to rotationally symmetric systems with high technical difficulty
3. Numerical Implementation: GPU-optimized code to be made publicly available, promoting field development
4. Systematic Verification: Multi-parameter systematic calculations validate theoretical validity

Weaknesses

1. Practical Limitations: Currently demonstrates only simple systems such as hydrogen atoms
2. Parameter Selection: Parameters chosen to highlight effects may deviate from actual conditions
3. Insufficient Comparison: Lacks direct comparison with experimental data

Impact

1. Academic Value: Provides important extension to open quantum system theory
2. Application Prospects: Offers new tools for cavity QED, quantum transport, and related fields
3. Methodological Contribution: Demonstrates the importance of maintaining physical symmetries in quantum dynamics

Applicable Scenarios

• Radiation damping studies in atomic physics
• Color center problems in condensed matter physics
• Strongly-coupled systems in cavity quantum electrodynamics
• Coulomb interactions in ionic liquids

References

The paper cites 75 relevant references covering key works in open quantum system theory, HEOM method development, cavity QED, and related fields, providing a solid theoretical foundation for this research.