⚛️ quantum physics

Work distribution for strongly coupled many-body open quantum systems

This paper extends the time-dependent numerical renormalization group (TDNRG) method to calculate the quantum work distribution function for strongly coupled many-body open systems, revealing non-perturbative power-law threshold behavior and universal scaling collapse driven by the Anderson orthogonality catastrophe.

Original authors: H. T. M. Nghiem, T. A. Costi, S. Campbell, A. K. Mitchell

Published 2026-08-12

📖 4 min read🧠 Deep dive

CC0 1.0

Original authors: H. T. M. Nghiem, T. A. Costi, S. Campbell, A. K. Mitchell

Original paper dedicated to the public domain under CC0 1.0 (http://creativecommons.org/publicdomain/zero/1.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

Imagine the universe as a giant, bustling dance floor. In this dance, tiny particles like electrons are the dancers, and they are constantly moving to the rhythm of their own energy. Sometimes, we want to change the music suddenly—maybe we speed up the beat or change the key. In the world of quantum physics, this sudden change is called a "quench." But here's the tricky part: these dancers don't just move alone; they are often holding hands with a massive crowd of other particles (the "bath" or environment). When you suddenly change the music for one dancer, the whole crowd has to react. Because they are all so tightly connected, the reaction isn't simple; it's a chaotic, collective shuffle that creates a ripple effect through the entire system.

Scientists have long been trying to figure out exactly how much "work" (energy) is required to make these sudden changes and how the system reacts. Usually, they could only calculate the average amount of work, like guessing the average height of everyone in the crowd. But they wanted to know the whole story: the full range of possibilities, from the smallest shuffle to the wildest jump. This is called the "work distribution." It's like knowing not just the average height, but the exact probability of finding a dancer who is 4 feet tall versus one who is 7 feet tall. Understanding this is crucial because it helps us build better quantum computers and tiny energy machines, but it's incredibly hard to solve when the dancers are holding hands tightly (strong coupling).

This paper is a breakthrough in how we solve that puzzle. The authors, a team of physicists from Vietnam, Germany, and Ireland, have developed a new, super-powerful mathematical tool to watch these quantum dances in extreme detail. They focused on two famous models of quantum systems: one where an electron interacts with a sea of other electrons (the Anderson impurity model) and another where a tiny quantum bit interacts with a field of vibrations (the spin-boson model). Using a method called Time-Dependent Numerical Renormalization Group (TDNRG), they simulated what happens when these systems are suddenly "quenched."

Their main finding is that when you look at the energy distribution right at the very bottom—the minimum amount of work needed—the system doesn't behave randomly. Instead, it follows a very specific, predictable pattern called a "power-law." Think of it like a waterfall: as the water gets closer to the bottom, it doesn't just stop; it speeds up in a very precise way. The authors discovered that this "waterfall" shape is caused by something called the "Anderson orthogonality catastrophe." In simple terms, when the system changes, the environment is so sensitive that the "before" and "after" states become completely different, like two songs that sound nothing alike even if they share a few notes. This difference forces the energy distribution to follow a strict mathematical rule.

The paper also shows that this behavior is "universal," meaning it doesn't matter exactly how strong the connections are between the particles; as long as you look at the energy low enough, the pattern is always the same. They found that this pattern is controlled by a new, hidden energy scale that emerges from the interactions themselves. For the electron model, this scale is related to the "Kondo temperature" (a measure of how strongly the electrons are screening a magnetic impurity), and for the vibration model, it's a "renormalized tunneling amplitude" (how easily the particle can jump between states after being slowed down by the environment).

To make sure their math was right, the team also checked their results against a famous rule called the "Crooks fluctuation theorem," which relates the work done in a forward process to the work done if you ran the movie backward. Their simulations matched this rule perfectly, even at different temperatures, proving their method is incredibly accurate. They didn't just guess; they simulated the entire system with such fine resolution that they could see energy scales that are exponentially small, something previous methods struggled to do.

In short, this paper doesn't just give us a new way to calculate energy; it reveals that even in the chaotic, messy world of strongly connected quantum particles, there is a hidden, elegant order waiting to be found right at the edge of the minimum energy. The authors suggest that this new tool could be used to explore even more complex quantum systems in the future, helping us understand the fundamental rules of how energy moves in the quantum world.

Technical Summary: Work Distribution for Strongly Coupled Many-Body Open Quantum Systems

Problem Statement
The paper addresses the challenge of defining and calculating thermodynamic work in strongly coupled many-body open quantum systems. While operational definitions of work exist via two-time projective energy measurements, obtaining exact results for the full quantum work distribution function (WDF) in the strong-coupling regime is notoriously difficult. In this regime, system-bath correlations are non-perturbative, and interactions generate emergent low-energy scales. Standard perturbative approaches or Markovian approximations fail to capture the collective, infrared-singular rearrangements of the environment (the Anderson orthogonality catastrophe) that occur following a sudden quench. The authors focus on "quantum impurity models"—specifically the Anderson Impurity Model (AIM) and the Spin-Boson Model (SBM)—where the system and bath are treated as a single composite entity.

Methodology
To solve the underlying non-equilibrium quantum impurity problem non-perturbatively, the authors extend the Time-Dependent Numerical Renormalization Group (TDNRG) method. Key technical aspects of their approach include:

Non-Perturbative Framework: The method treats the system and bath on an equal footing in the thermodynamic limit, avoiding weak-coupling approximations.
Complete-Basis Formulation: The authors utilize a complete-basis formulation of TDNRG. They construct a many-body basis from the discarded states of all NRG iterations, supplemented by environment states. This allows for the calculation of the characteristic function $G(u) = \langle e^{iuH_f} e^{-iuH_i} \rangle$ , from which the WDF $P(W)$ is obtained via Fourier transform.
Resolution and Scaling: By leveraging the logarithmic discretization of the Wilson chain, the method provides exponentially fine energy resolution at low energies. This is crucial for resolving the universal scaling behavior near the work threshold.
Model Mapping: For the Ohmic Spin-Boson Model, the authors employ a mapping to the equivalent spinless Interacting Resonant Level Model (IRLM) to facilitate more efficient TDNRG treatment.

Key Results
The study analyzes the WDF following a sudden quench in two paradigmatic models:

Anderson Impurity Model (AIM): The authors investigate level quenches where the dot potential $\epsilon_d$ $ϵ_{d}$ is changed.
- Threshold Behavior: At zero temperature, the WDF exhibits a power-law threshold behavior as the work $W$ approaches its minimum value $W_{min}$ : $P(W) \sim |W - W_{min}|^{-\beta}$ .
- Orthogonality Exponent: The exponent $\beta$ is determined by the Anderson orthogonality catastrophe (OC). Specifically, $\beta = 1 - (n_f - n_i)^2/2$ , where $n_i$ and $n_f$ are the initial and final local level occupations.
- Universal Scaling: Below an emergent low-energy scale (the Kondo temperature $T_K$ ), the distribution exhibits universal scaling collapse. The scaled distribution $P(W) \times (W - W_{min})^\beta$ collapses onto a single curve for varying interaction strengths $U$ or occupation changes.
Spin-Boson Model (SBM): The authors examine bias quenches from a fully polarized state to an unbiased state.
- Dissipation Dependence: The threshold exponent is found to be $\beta = 1 - \alpha/2$ , where $\alpha$ is the dimensionless dissipation strength. This result is derived via the mapping to the IRLM, where the orthogonality exponent relates to the change in scattering phase shifts.
- Scaling: Similar to the AIM, universal data collapse is observed below the renormalized tunneling amplitude $\Delta_r$ , which acts as the emergent low-energy scale.
Finite Temperature Validation: The numerical approach is validated at finite temperatures by verifying the Crooks fluctuation theorem and checking that the first and second moments of the WDF match analytical expectations derived from the initial state properties.

Significance and Claims
The paper claims to provide a numerically exact solution for the full WDF of strongly coupled open quantum systems, a task previously intractable for true many-body problems. The primary significance lies in:

Direct Probe of RG Flow: The WDF is shown to act as a direct thermodynamic probe of the impurity renormalization group (RG) flow, resolving both the infrared fixed-point structure and the interaction-renormalized crossover scales ( $T_K$ and $\Delta_r$ ).
Resolution of Infrared Singularities: The TDNRG method's logarithmic energy separation allows for the direct access of exponentially small many-body scales and the resolution of the universal many-body continuum of work processes near the threshold, which is computationally demanding for other methods like Monte Carlo or real-time process tensors.
Universality: The work demonstrates that the threshold behavior of the WDF is universally governed by the Anderson orthogonality catastrophe, with exponents fixed by the change in low-energy scattering phase shifts, regardless of the specific details of the quench amplitude.

The authors note that this framework opens the door to studying work statistics in more complex quantum impurity problems, including multichannel systems and quenches involving interaction parameters, and suggests that finite-time driving could be explored via sequences of quenches within this formalism.

Technical Summary: Work Distribution for Strongly Coupled Many-Body Open Quantum Systems

More like this