Hadamard regularization of open quantum systems coupled to unstructured environments in the Schwinger-Keldysh formalism

This paper proposes a Hadamard regularization-based separation-of-scales ansatz within the Schwinger-Keldysh formalism to develop a computationally efficient time-stepping algorithm for the Kadanoff-Baym equations, enabling the simulation of damped quantum harmonic oscillators in unstructured environments while capturing non-Markovian and renormalization effects without prohibitive cubic scaling.

The Gist

The Big Picture: The "Noisy Room" Problem

Imagine you are trying to listen to a single, delicate violin playing a beautiful melody in a massive, crowded concert hall. The violin is your quantum system (like a tiny computer chip or an atom). The crowd is the environment (the air, the walls, other people).

In the real world, the violin never plays in total silence. The crowd makes noise, bumps into the violin, and absorbs its energy. This causes the music to fade (damping) and the notes to lose their perfect pitch (decoherence).

Physicists want to predict exactly how the violin sounds in this noisy room. To do this, they use a powerful mathematical tool called the Schwinger-Keldysh formalism. Think of this tool as a super-advanced recording device that can capture every tiny interaction between the violin and the crowd.

The Problem: The "Too-Fast" Crowd

The problem is that the "crowd" in these quantum models is often treated as featureless and infinitely fast. It's like the crowd is made of millions of tiny, invisible particles buzzing around at speeds we can't even imagine.

When the physicists try to simulate this using a computer, they hit a wall:

1. The Speed Trap: Because the crowd moves so fast, the computer has to take tiny, tiny steps in time to keep up. It's like trying to film a hummingbird's wings with a camera that can only take one photo per second; you miss everything.
2. The Explosion: To get an accurate picture, the computer needs to take steps so small that the calculation time explodes. If the crowd is 1,000 times faster than the violin, the computer takes $1,000^3$ (one billion) times longer to run the simulation. This makes it impossible to study complex scenarios.
3. The Singularity: Mathematically, when you try to ignore the crowd's speed (make it infinite), the equations break down and give you "infinity" as an answer. It's like dividing by zero.

The Solution: The "Smart Blur" (Hadamard Regularization)

The author, Jakob Dolgner, proposes a clever new way to solve this. Instead of trying to track every single buzzing particle in the crowd, he suggests separating the fast from the slow.

Here is the analogy:
Imagine you are watching a car race.

• The Slow Part: You care about where the cars are on the track (the system).
• The Fast Part: You don't care about the individual pixels of the asphalt or the vibration of the engine bolts (the environment).

Dolgner's method is like using a smart blur.

1. The Blur: He acknowledges that the "fast" part of the environment is too fast to see. Instead of trying to resolve it, he mathematically "blurs" it out.
2. The Fix (Hadamard Regularization): When you blur something, you usually lose detail. But in math, this "blur" creates a specific kind of error (a singularity). Dolgner uses a technique called Hadamard Regularization. Think of this as a special "mathematical eraser" that removes the infinite noise but keeps the shape of the effect. It's like saying, "We know the crowd is infinitely fast, so let's just calculate the average push they give the violin, ignoring the impossible details."

The New Algorithm: The "Two-Speed Camera"

The paper introduces a new computer algorithm that acts like a two-speed camera:

• Fast Mode: It instantly calculates the "average push" (the renormalization) that the fast environment gives the system. This fixes the pitch of the violin immediately.
• Slow Mode: It then simulates the violin's movement using normal, manageable time steps. It doesn't waste time watching the invisible, fast particles.

This allows the computer to simulate the system as if the environment were infinitely fast, without actually having to calculate every single fast step.

Why This Matters

1. Solving the "Infinity" Problem: The paper proves that even though the math looks like it blows up (goes to infinity), the "blur" technique makes it finite and solvable. It turns a broken equation into a working one.
2. Capturing the "Ghost" Effects: Even though the environment is fast, it leaves a "ghost" on the system. The system remembers the past interactions (non-Markovianity). This new method captures these memories perfectly, even when the system is very cold (where quantum effects are strongest).
3. Real-World Use: This is a game-changer for designing quantum computers and new materials. It allows scientists to simulate how quantum devices behave in the real, noisy world without needing a supercomputer the size of a city.

Summary in One Sentence

The author invented a mathematical "smart blur" that lets computers simulate how fast, invisible noise affects delicate quantum systems, solving a problem that previously made the calculations take forever or break entirely.

Technical Summary

1. Problem Statement

The paper addresses a critical computational bottleneck in the numerical simulation of open quantum systems, specifically those coupled to "featureless" (unstructured) environments like Ohmic baths.

• Context: The Schwinger-Keldysh (SK) formalism is a rigorous framework for non-equilibrium quantum systems, capable of capturing non-Markovian dynamics and exact correlation functions. However, its numerical solution via the Kadanoff-Baym equations (KBE) typically scales cubically with the number of time steps ($O(N^3)$).
• The Stiffness Issue: In standard treatments of unstructured environments (e.g., the Wide-Band Limit or Ohmic limit), the environment is modeled with a high-frequency cutoff $\omega_c$ to regularize divergences. To accurately resolve the physics, the time step $\Delta t$ must be small enough to resolve the fastest scale, $1/\omega_c$.
• The Conflict: In many physical scenarios, the environment scale $\omega_c$ is orders of magnitude larger than the system's characteristic frequency $\omega_0$ (i.e., $\omega_c \gg \omega_0$). Resolving $\omega_c$ forces an extremely small time step, making simulations prohibitively expensive ($O((\omega_c/\omega_0)^3)$) for long-time dynamics.
• Divergences: Taking the limit $\omega_c \to \infty$ directly leads to mathematical singularities, specifically a logarithmic divergence in the momentum variance $\langle \pi^2 \rangle$ and a singularity in the second time derivative of the symmetric correlator on the time-diagonal.

2. Methodology

The author proposes a novel regularization and time-stepping algorithm based on Hadamard regularization and scale separation.

A. Hadamard Finite Part (HFP) Interpretation

The paper reinterprets the embedding self-energies ($\Sigma_{em}$) arising from the environment as distributions.

• Symmetric Self-Energy ($\Sigma_S$): In the Ohmic limit, the symmetric component of the self-energy (related to fluctuations) behaves like a distribution with a singularity of order $1/t^2$. The author identifies this as a nascent Hadamard Finite Part (HFP) distribution.
• Regularization: Instead of introducing an artificial high-frequency cutoff $\omega_c$ that requires fine time steps, the theory treats the singular integrals using the HFP definition. This allows for an analytical treatment of the $\omega_c \to \infty$ limit for off-diagonal times.

B. Scale Separation Ansatz

The core algorithmic innovation is the decomposition of the convolution integral $(\Sigma \cdot G)(t_1, t_2)$ into "fast" and "slow" components:

1. Taylor Expansion: The correlator $G(t', t_2)$ is expanded around the critical point $t' = t_1$ (where the singularity occurs):
   $$G(t', t_2) \approx G(t_1, t_2) + \partial_{t_1}G(t_1, t_2)(t' - t_1) + R_2(t', t_1, t_2)$$
2. Term Isolation:
   • Zeroth and First Order Terms: These terms ($P$ and $Q$) contain the explicit dependence on the cutoff $\omega_c$ and the singular behavior. They are integrated analytically.
   • Remainder Term ($R_2$): This term scales as $(t' - t_1)^2$, which cancels the singularity of the kernel, rendering the integral regular and slowly varying.
3. Result: The fast, stiff, $\omega_c$-dependent part is isolated into local terms ($P$ and $Q$), while the history-dependent part becomes smooth and can be resolved on the coarse system time scale ($\Delta t \sim 1/\omega_0$).

C. Numerical Implementation

• Time-Stepping: The KBEs are solved using a time-stepping algorithm (similar to Verlet or Exponential Time Differencing) adapted for the coarse grid.
• Filon-Type Quadrature: To handle the integration of the stiff local terms ($P$ and $Q$) against the oscillating Green's functions, the author employs Filon-type quadrature. This method is specifically designed for integrals with oscillatory kernels and allows for high accuracy without requiring the grid to resolve the high-frequency cutoff.
• Renormalization: The method effectively renormalizes the momentum variance. The logarithmic divergence is absorbed into the definition of the initial conditions or the bare parameters, allowing the simulation to proceed with $\omega_c \to \infty$ (or effectively infinite) without numerical instability.

3. Key Contributions

1. Resolution of the Stiffness Problem: The paper provides a rigorous method to simulate open quantum systems with unstructured environments using time steps determined solely by the system's dynamics ($\omega_0$), rather than the environment's cutoff ($\omega_c$). This reduces computational cost from cubic scaling with the scale separation to a manageable level.
2. Hadamard Regularization Framework: It establishes that the singularities in the KBEs for Ohmic baths are naturally handled by Hadamard Finite Part distributions, offering a mathematically consistent way to take the wide-band limit ($\omega_c \to \infty$).
3. Non-Markovian Dynamics at Low Temperatures: The method successfully captures strong non-Markovian effects and quantum fluctuations in the ultra-cold regime ($T \ll \gamma$), where standard Markovian approximations (like Lindblad equations) fail.
4. Analytical Derivations: The paper derives exact analytical forms for the time-domain self-energies and the coefficients ($P$ and $Q$) required for the numerical scheme, including their limits for large $\omega_c$.

4. Results

• Benchmarking: The algorithm was tested on a damped quantum harmonic oscillator (DQHO).
  • Thermalization: The system correctly thermalizes to the bath temperature $T$ across a wide range of parameters (damping $\gamma$, temperature $T$, and cutoff $\omega_c$).
  • Momentum Variance: The numerical results for the momentum variance $\langle \pi^2 \rangle$ match the exact analytical predictions, confirming that the logarithmic divergence is correctly regularized.
  • Correlation Decay: In the ultra-cold regime ($T < \gamma$), the symmetric correlator $G_S(t, t)$ exhibits a power-law decay ($t^{-2}$) at intermediate times, a feature absent in Markovian (Lindblad) simulations which predict immediate exponential decay.
• Efficiency: The method allows for simulations where $\omega_c \sim 10^5 \omega_0$ without increasing the computational cost relative to the isolated system, effectively bypassing the stiffness constraint.

5. Significance

• Bridging Scales: This work enables the study of high-bandwidth, low-temperature environments (common in condensed matter and quantum optics) using first-principles field theory without resorting to uncontrolled approximations like the Born-Markov approximation.
• General Applicability: While demonstrated on bosonic systems, the authors note that the regularization procedure is directly applicable to fermionic systems (wide-band limit) and can be generalized to non-Ohmic baths (sub-Ohmic and super-Ohmic), provided the singularity order is handled correctly.
• Theoretical Insight: It clarifies the nature of divergences in open quantum systems, showing that they are artifacts of the continuum limit that can be systematically removed via coarse-graining and renormalization, rather than fundamental physical obstructions.

In summary, Dolgner presents a robust numerical framework that leverages the mathematical properties of Hadamard finite parts to decouple the simulation time scale from the environmental cutoff, enabling efficient and exact simulations of non-Markovian open quantum systems.