RC-HEOM Hybrid Method for Non-Perturbative Open System Dynamics

1. Problem Statement

Open quantum system dynamics often involve strong system-bath coupling and non-Markovian memory effects, which render standard perturbative approaches (like the Lindblad master equation) invalid. Two primary non-perturbative methods exist, but both have significant limitations:

• Hierarchical Equations of Motion (HEOM): Provides numerically exact solutions for systems coupled to Gaussian baths by encoding bath memory into Auxiliary Density Operators (ADOs). However, it treats the bath as a "black box," making it difficult to access specific details about the bath's collective modes or system-bath correlations.
• Reaction Coordinate (RC) Mapping: Transforms the system-bath Hamiltonian to explicitly map a dominant collective bath mode (the RC) onto the system. This allows direct access to the system-RC joint state. However, standard implementations treat the residual bath (the remainder of the bath after RC extraction) using perturbative Markovian approximations (Lindblad master equations), which fails when the residual bath coupling is strong or highly structured.

The Core Challenge: There is a lack of a method that simultaneously offers numerical exactness (handling strong coupling/non-Markovian effects of the residual bath) and explicit access to the dominant bath mode (RC) for analyzing system-environment correlations.

2. Methodology: RC–HEOM

The authors propose RC–HEOM, a hybrid framework that unifies the strengths of both approaches.

• Step 1: Reaction Coordinate Mapping:
  The system is linearly coupled to a bosonic or fermionic bath. A Bogoliubov transformation is applied to the bath operators to define a new set of modes. This maps the original Hamiltonian into:
  $$H_{tot} = H_{sys+RC} + H_{residual}$$
  Where $H_{sys+RC}$ describes the system coupled to a single collective mode (the RC), and $H_{residual}$ describes the RC coupled to the remaining "residual" bath modes.

• Step 2: Non-Perturbative Treatment of Residual Bath:
  Instead of applying a Markovian master equation to the $sys+RC$ system interacting with the residual bath, the authors apply the HEOM formalism.

  • The residual bath's influence is captured via a two-time correlation function, which is decomposed into a sum of exponentials.
  • A hierarchy of Auxiliary Density Operators (ADOs) is constructed to track the system-RC correlations with the residual bath exactly.
  • The equations of motion for the ADOs are solved to obtain the reduced density matrix of the $sys+RC$ joint state ($\rho_{sys+RC}$).

• Key Advantage: This approach retains the exact non-Markovian memory of the residual bath (via HEOM) while providing explicit access to the RC mode, enabling the calculation of system-RC correlations.

3. Key Contributions

1. Development of RC–HEOM: A novel hybrid algorithm that combines RC mapping with HEOM, overcoming the limitations of both the standard RC-ME (perturbative) and standard HEOM (lack of bath mode resolution).
2. Validation Framework: The method is benchmarked against:
   • HEOM:* Standard HEOM without RC mapping (serves as the ground truth for exact dynamics).
   • RC–ME: RC mapping with a Lindblad master equation (standard approximate approach).
3. Application to Anderson Impurity Models: The method is applied to two distinct physical regimes to demonstrate its utility in accessing previously unobservable quantities.

4. Results

Example 1: Single-Impurity Anderson Model (SIAM) & Kondo Physics

• Context: Studying the formation of the Kondo singlet at low temperatures.
• Findings:
  • Density of States: RC–HEOM reproduces the correct Kondo resonance and Hubbard resonance, matching the HEOM* benchmark. In contrast, RC–ME fails to produce the correct resonance due to the breakdown of the Born-Markov approximation in the specific bath spectral density regime.
  • Singlet Fraction: The authors calculated the singlet fraction ($F$), a measure of how close the $sys+RC$ state is to a spin singlet. They observed that $F$ increases as temperature decreases, correctly tracking the emergence of the Kondo singlet.
  • Significance: RC–HEOM successfully captures strong-coupling Kondo physics where RC–ME fails.

Example 2: Two-Impurity Anderson Model (TIAM) & Coherence Revival

• Context: Investigating bath-mediated coherence revival between two impurities.
• Findings:
  • Coherence Dynamics: Both HEOM* and RC–HEOM correctly predict a sharp revival in the $l_1$ norm of coherence at a specific time ($Wt' \approx 97.9$), whereas RC–ME fails to capture this.
  • Mechanism Unveiling: By accessing the $sys+RC$ joint state, the authors decomposed the coherence into "one-fermion" (OF) coherence terms. They identified that the revival is driven by the interference of four distinct paths mediated by the RC (vacuum, spin-up, spin-down, and doubly occupied RC states).
  • Interference Factor: An interference factor ($I$) was defined to quantify the constructive/destructive interference of these paths. The results showed that the revival occurs when the interference shifts from destructive to constructive.
  • Significance: This provides a microscopic understanding of how the bath mediates coherence, a level of detail impossible with standard HEOM (which lacks RC resolution) or RC–ME (which is inaccurate).

Bosonic Systems (Supplemental Material)

• The method was also validated on a spin-boson model (Jaynes-Cummings). RC–HEOM correctly recovered short-time dynamics that RC–ME failed to capture, confirming its applicability to bosonic environments.

5. Significance and Impact

• Bridging the Gap: RC–HEOM fills a critical gap in open quantum system theory by allowing researchers to study regimes where the system-bath coupling is strong and non-Markovian, while still retaining physical insight into the bath's collective modes.
• New Observables: It enables the calculation of quantities like the "singlet fraction" and "system-RC coherence," which are essential for understanding quantum thermodynamics, entanglement, and transport in strongly coupled regimes.
• Future Applications: The authors suggest this method is particularly valuable for:
  • Quantum thermodynamic devices and heat engines in strong coupling.
  • Quantum transport and ultrastrong light-matter interactions.
  • Non-Markovian heat and work statistics.

In conclusion, RC–HEOM represents a significant advancement in computational quantum dynamics, offering a versatile and exact tool for characterizing complex open quantum systems where conventional master-equation approaches are insufficient.