Linear response and exact hydrodynamic projections in Lindblad equations with decoupled Bogoliubov hierarchies

This paper investigates a class of spinless-fermion Lindblad equations with decoupled BBGKY hierarchies, utilizing a mapping to non-Hermitian Schrödinger equations to derive exact hydrodynamic projections and linear response functions for systems exhibiting diffusive dynamics and infinite temperature steady states.

The Gist

Imagine a crowded dance floor where everyone is trying to move, but there's a twist: the music is glitching, and occasionally, people are randomly bumped by invisible ghosts. This is a bit like what happens in the quantum world when particles interact with their environment. This interaction is called dissipation or noise, and it usually makes predicting how the system behaves incredibly difficult.

This paper by Patrik Penc and Fabian Essler is like finding a secret cheat code that allows us to predict exactly how this noisy dance floor behaves, even when things get chaotic.

Here is a breakdown of their discovery using everyday analogies:

1. The Problem: The "Impossible" Puzzle

Usually, when you have a system of many particles (like electrons) that are losing energy to their surroundings, the math becomes a tangled mess. It's like trying to predict the path of every single drop of water in a crashing wave while a hurricane is blowing. The equations get so complex that scientists usually have to rely on approximations or supercomputers that can only guess the answer.

2. The Solution: The "Decoupled" Ladder

The authors studied a specific class of these noisy systems. They discovered that for these particular models, the complex web of interactions untangles itself.

The Analogy: Imagine a giant, multi-story building where every floor depends on the one below it. Usually, if you want to know what's happening on the 10th floor, you have to solve the equations for floors 1 through 9 first. It's a nightmare.
In these specific models, the authors found that the building is actually made of separate, independent elevators. The 10th floor doesn't care about the 9th. This is what they call a "decoupled hierarchy." Because the floors don't talk to each other, you can solve the math for each "floor" (or group of particles) independently.

3. The Magic Trick: Turning Time Backwards

Once they separated the problem, they found a way to turn the messy "noise" equations into something much simpler: a Schrödinger equation (the basic rulebook of quantum mechanics) but with a twist.

The Analogy: Normally, quantum mechanics is like a movie playing forward. But because of the noise, the authors realized they could rewrite the rules so the movie plays in reverse (imaginary time) and the "actors" (particles) have negative health points (non-Hermitian Hamiltonians).
This sounds scary, but it's actually a superpower. It turns a complex, messy problem into a standard math problem that physicists know how to solve perfectly.

4. The Big Discoveries

A. The "Hydrodynamic Projection" (The Slow Drift)

In a noisy system, some things move fast and die out quickly, while others move very slowly and stick around. The authors figured out exactly how to spot the "slow movers."

The Analogy: Imagine throwing a handful of confetti into a windy room. Most of it flies away instantly. But a few pieces get caught in a slow, swirling draft near the ceiling.
The authors calculated exactly how much of the original confetti ends up in that slow draft. They found that these slow pieces don't just fade away; they spread out in a specific, predictable pattern (like a drop of ink spreading in water). This is called diffusion. They gave a precise formula for how this spreading happens, which is a huge deal because usually, we can only guess this behavior.

B. The "Pump and Probe" (Linear Response)

Scientists often study materials by hitting them with a pulse of energy (a "pump") and watching how they react (the "probe").
The Analogy: Think of tapping a bell. In a perfect world, the bell rings with a clear, sharp tone. But in a noisy, dampened world, the sound is muffled and blurry.
The authors showed how to calculate exactly how the "bell" sounds when it's being tapped while it's being shaken by the wind (the noise). They found that the noise smooths out the sharp edges of the sound, turning a clear ring into a dull thud, but they can predict exactly what that thud will look like.

C. Integrable vs. Non-Integrable (The Orderly vs. The Chaotic)

Some of the models they studied are "integrable," meaning they have hidden symmetries that keep them orderly (like a perfectly choreographed dance). Others are "non-integrable" and chaotic.
The Surprise: Even though the chaotic models usually behave very differently from the orderly ones, in this specific noisy environment, they end up behaving almost the same way in the long run. They both settle into the same "infinite temperature" state (a state of maximum chaos where everything is mixed up). This suggests that noise can sometimes wash away the differences between orderly and chaotic systems.

Why Does This Matter?

This paper is a toolkit. By showing that these specific noisy systems can be solved exactly, the authors have given scientists a way to:

1. Predict how quantum computers might behave when they get noisy (which is a huge problem for building real quantum computers).
2. Understand how heat and electricity move through materials that are constantly interacting with their environment.
3. Design better experiments to measure these effects without needing to guess.

In short, they took a jumbled, impossible-to-solve puzzle, found a way to separate the pieces, and showed us exactly how the picture forms over time, even in the middle of a storm.

Technical Summary

1. Problem Statement

The dynamics of open quantum many-body systems, described by Lindblad equations (LEs), are generally intractable due to the exponential growth of the Hilbert space and the complexity of non-Markovian or non-Gaussian interactions. While perturbative and numerical methods exist, there is a scarcity of exact analytical results for non-equilibrium dynamics, particularly regarding:

• Hydrodynamic Projections: Determining how microscopic lattice operators project onto slow, diffusive hydrodynamic modes at late times. This is crucial for understanding power-law tails in correlation functions.
• Linear Response in Non-Equilibrium: Calculating response functions (e.g., pump-probe scenarios) on top of dissipative time evolution, which is generally difficult.
• Integrability vs. Dissipation: Understanding how Yang-Baxter integrability manifests in dissipative systems, especially since many integrable LEs relax to infinite-temperature steady states (lacking traditional conservation laws).

The authors focus on a specific class of spinless-fermion Lindblad equations where the Bogoliubov–Born–Green–Kirkwood–Yvon (BBGKY) hierarchy decouples. Unlike standard free-fermion systems, these models are interacting (no Wick's theorem applies to the density matrix), yet the hierarchy of equations of motion for $n$-point functions closes at a finite order $n$.

2. Methodology

The authors employ a vectorization technique to map the Heisenberg-picture time evolution of operators onto a few-particle imaginary-time Schrödinger equation with a non-Hermitian Hamiltonian.

• Vectorization: Operators of the form $\hat{\Psi}_{n,m} = \sum \Psi(\vec{j}, \vec{\ell}) c^\dagger_{j_1} \dots c^\dagger_{j_n} c_{\ell_1} \dots c_{\ell_m}$ are mapped to states in a doubled Hilbert space (spin-up and spin-down fermions).
• Non-Hermitian Mapping: The dual Lindbladian $\mathcal{L}^*$ acting on operators becomes a Hamiltonian $H$ acting on these vectorized states:
  $$\frac{d}{dt} |\Psi(t)\rangle = H |\Psi(t)\rangle$$
  The Hamiltonian $H$ describes interacting spin-1/2 fermions with purely imaginary hopping terms and short-range interactions.
• Spectral Analysis: The time evolution is solved by expanding the initial state in the eigenbasis of $H$:
  $$|\Psi(t)\rangle = \sum_\alpha \langle \Phi_{L,\alpha} | \Psi(0) \rangle e^{\lambda_\alpha t} |\Phi_{R,\alpha}\rangle$$
  The authors analyze the eigenvalue spectrum $\lambda_\alpha$ of these non-Hermitian Hamiltonians for various models (I–V), distinguishing between integrable (Model I) and non-integrable cases.
• Hydrodynamic Projection: By identifying eigenmodes with eigenvalues $\lambda(p) \approx -Dp^2$ (diffusive modes) near zero momentum, the authors calculate the overlap between local operators and these slow modes to derive exact late-time asymptotics.

3. Key Contributions and Results

A. Exact Hydrodynamic Projections

The paper derives closed-form expressions for the hydrodynamic projections of operators quadratic in fermions (e.g., $c^\dagger_x c_y$).

• Diffusive Eigenoperators: The authors explicitly construct the diffusive eigenmodes of the Lindbladian. These correspond to "particle-hole bound states" in the vectorized picture, where the wavefunction decays exponentially with distance.
• Late-Time Tails: For spatially local operators, the time-evolved operator takes the form:
  $$(c^\dagger_{x_0} c_{y_0})(t) \simeq \sum_{x,y} \psi_{x_0,y_0}(x,y;t) c^\dagger_x c_y$$
  The amplitude $\psi$ exhibits power-law tails at late times ($t \gg \gamma^{-1}$):
  $$\psi \sim (Dt)^{-1 + X/2} \quad (\text{for even } Y)$$
  $$\psi \sim (Dt)^{-2 + X/2} \quad (\text{for odd } Y)$$
  where $X = |x-y| + |x_0-y_0|$ and $Y = (x-x_0) + (y-y_0)$. The diffusion constant $D$ is model-dependent (e.g., $D = J^2/\gamma$ for Model I).
• Momentum Dependence: Operators carrying definite non-zero momentum decay exponentially and do not exhibit hydrodynamic power-law tails, a result of the specific kinematics of the bound states.

B. Linear Response in Non-Equilibrium

The authors apply their formalism to calculate linear response functions for systems driven by a quantum quench (starting from a dimerized tight-binding ground state) and subsequently subjected to dissipation.

• Dissipative Smearing: They show that dissipation washes out sharp features (like square-root singularities) present in the unitary response functions.
• Emergence of Diffusive Modes: At low frequencies and large dissipation rates, signatures of the diffusive hydrodynamic mode emerge in the response function, distinct from the ballistic features of the unitary limit.

C. Integrability and Spectral Structure

The study compares Yang-Baxter integrable models (Model I: on-site dephasing) with non-integrable models (Models II–IV: bond dephasing).

• Integrable Case (Model I): The spectrum of the non-Hermitian Hamiltonian can be solved exactly using the Bethe Ansatz. The diffusive modes correspond to specific "k-$\Lambda$ string" solutions (analogous to bound states in the Hubbard model).
• Non-Integrable Cases: While exact Bethe ansatz solutions are unavailable, the decoupled hierarchy allows for exact numerical diagonalization of the few-body problem. The authors find that the qualitative structure of the diffusive modes (bound states with real eigenvalues vanishing as $p^2$) persists even in non-integrable models, though the exact wavefunctions differ.
• Higher-Order Operators:
  • Cubic operators: Decay exponentially (no hydrodynamic tails) due to odd fermion parity.
  • Quartic operators: In the integrable case, they are expected to show hydrodynamic tails, but the wavefunctions are no longer simple Bethe ansatz states (requiring nested Bethe ansatz), making analytical treatment significantly harder.

D. Models Studied

The authors analyze five specific models of spinless fermions with different jump operators:

1. Model I: On-site dephasing (Integrable, maps to imaginary Hubbard model).
2. Model II: Two-channel bond dephasing (Non-integrable).
3. Model III & IV: Single-channel bond dephasing variations (Non-integrable).
4. Model V: Jump operators breaking U(1) symmetry (Particle number not conserved; no hydrodynamic tails, purely exponential decay).

4. Significance

• Exactness in Non-Equilibrium: This work provides one of the few instances where exact, closed-form results for hydrodynamic projections and linear response are derived for interacting, dissipative many-body systems.
• Beyond Gaussianity: Unlike previous exact solutions for open systems (which often rely on quadratic Lindbladians preserving Gaussianity), this work handles models where the density matrix remains non-Gaussian, yet the equations of motion for operators decouple.
• Universality of Hydrodynamics: The results suggest that the existence of diffusive hydrodynamic modes and their associated power-law tails is a robust feature of systems with conserved quantities and decoupled hierarchies, regardless of whether the underlying model is Yang-Baxter integrable.
• New Perspective on Integrability: The paper clarifies that in dissipative systems, integrability simplifies the form of the eigenfunctions (Bethe ansatz) but does not necessarily imply the existence of non-trivial conservation laws (as steady states are often infinite temperature).
• Experimental Relevance: The derived linear response functions provide theoretical benchmarks for pump-probe experiments in open quantum systems, predicting how dissipation modifies spectral lineshapes and transport properties.

In summary, the paper establishes a powerful framework linking the decoupling of the BBGKY hierarchy to exact solvability in the Heisenberg picture, enabling the precise characterization of long-time hydrodynamic behavior and non-equilibrium response in a class of interacting open quantum systems.