Quantum Statistics and Non-Equilibrium Dynamics

• Quantum Statistics and Non-Equilibrium Dynamics is a framework examining quantum ensembles and emergent behaviors from external driving and environmental coupling.
• The field employs advanced statistical ensembles, fluctuation relations like the Jarzynski equality, and experimental techniques such as programmable photonic simulations to analyze work distributions.
• Recent findings reveal that bosonic interference, adiabatic–nonadiabatic crossovers, and non-Markovian effects critically shape entropy production and system relaxation.

Quantum statistics and non-equilibrium dynamics constitute a fundamental research axis at the intersection of condensed matter, AMO physics, and statistical mechanics. This topic encompasses both the microscopic statistical structure of quantum systems and the emergent behavior induced by external driving, interactions, and environmental coupling, with a particular emphasis on the breakdown of equilibrium, entropy production, and the appearance of novel dynamical steady states or critical phenomena. Experimental advances in quantum simulation, programmable hardware, and fast measurement now allow direct access to quantum work statistics, fluctuation relations, and dynamical indicators across a variety of platforms.

1. Foundations: Quantum Statistical Ensembles and Non-Equilibrium Concepts

Quantum statistics describes the probabilistic structure of quantum ensembles—including Gibbs, generalized Gibbs (GGE), and their extension to quantum jump trajectories—and underlies any definition of thermodynamic observables such as work, heat, or entropy. In equilibrium, the maximum-entropy principle specifies the density matrix: $\rho_{\rm eq} = \frac{e^{-\beta H}}{Z}\,,$ where $\beta=1/k_B T$ and $Z$ is the partition function. In integrable systems with extensive local conservation laws $I_\alpha$, equilibration is to a GGE: $$\rho_{\rm GGE} = \frac{1}{Z_{\rm GGE}} \exp\left(-\sum_\alpha \lambda_\alpha I_\alpha\right) \,,$$ with $\lambda_\alpha$ fixed by initial conditions (Altman, 2015, Garrahan, 2017). In open systems, quantum statistics extends to density operators evolving under master equations or Keldysh field theory.

Non-equilibrium dynamics arise in protocols where the system is driven out of equilibrium by parameter quenches, external fields, or coupling to engineered baths. These settings demand statistical frameworks capturing non-adiabatic transitions, entropy production, fluctuation theorems, and full distribution functions of thermodynamic observables.

2. Fluctuation Relations and Quantum Work Statistics

Out-of-equilibrium quantum thermodynamics is fundamentally structured by quantum fluctuation relations, such as the Jarzynski equality: $\langle e^{-\beta W}\rangle = e^{-\beta \Delta F}$ where $W$ is the work, and $\Delta F$ the free energy difference. The experimental measurement of $P(W)$ employs two-point projective schemes: the work distribution is

$\beta=1/k_B T$0

with $\beta=1/k_B T$1 the transition probability, determined by the quantum dynamics (Krishna et al., 11 Mar 2026).

Recent experiments simulate a two-boson piston protocol on programmable photonic hardware, leveraging quasi-unitary embedding to reconstruct $\beta=1/k_B T$2 for both expansion and compression scenarios, directly observing the crossover from adiabatic to non-adiabatic regimes, and confirming the Jarzynski equality at the quantum level (Krishna et al., 11 Mar 2026). Bosonic statistics—encoded in matrix permanents for indistinguishable particles—directly shape the statistics of work and the protocol-to-protocol variation in dissipation and irreversibility.

3. Non-Equilibrium Dynamics: Adiabatic–Nonadiabatic Crossover, Dissipation, and Entropy Production

Non-equilibrium protocols such as finite-speed boundary deformations induce complex quantum dynamics, including non-adiabatic transitions and entropy production. In the quantum piston experiment, the wall speed $\beta=1/k_B T$3 tunes the crossover:

• Quasi-adiabatic regime ($\beta=1/k_B T$4): Populations remain close to their initial (thermal) weights; $\beta=1/k_B T$5 is sharply peaked and dissipation negligible.
• Strongly non-adiabatic regime ($\beta=1/k_B T$6): Significant off-diagonal Fock state occupation, broadening, and skewing of $\beta=1/k_B T$7, with positive-work tails in expansions and high-work tails in compressions.
• Bosonic interference modifies transition probabilities, giving rise to unique many-body statistical signatures inaccessible to distinguishable particle dynamics (Krishna et al., 11 Mar 2026).

Dissipated work $\beta=1/k_B T$8 quantifies irreversibility, with entropy production $\beta=1/k_B T$9. Irreversibility can also be assessed via geometric measures such as the Bhattacharyya coefficient between initial and final distributions, with more rapid protocols leading to greater dissipation and reduced state overlap.

4. Approaches: Full Counting Statistics, Large Deviations, and Quantum Trajectory Methods

Theoretical characterization of quantum non-equilibrium fluctuations relies on several complementary frameworks:

• Full Counting Statistics (FCS): FCS encodes all moments (cumulants) of time-integrated observables such as work, heat, or transferred charge. In open systems, the characteristic function is constructed via two-point measurement protocols or by dressing Hamiltonians with counting fields, leading to expressions such as

$Z$0

or via path-ordered exponentials in Keldysh space (Cerrillo et al., 2016, Guarnieri et al., 2017). In this way, cumulant generating functions, large deviation rate functions, and trajectory-level fluctuation relations can be established.

• Large Deviation Theory: The exponential suppression of rare events is captured by scaled cumulant generating functions and their Legendre transforms, revealing dynamical phase transitions in trajectory space and first-order kink singularities in the dynamical free energy (Garrahan, 2017, Guarnieri et al., 2017).
• Quantum Trajectories: For open quantum systems (Lindblad evolution), unraveling the master equation yields stochastic quantum trajectories. Statistical mechanics of these trajectories—counting quantum jumps—directly analogizes classical non-equilibrium processes, with tilted generators providing access to the full trajectory statistics and dynamical phase transitions (Garrahan, 2017).

5. Strong Coupling, Non-Markovianity, and Memory Effects

Non-equilibrium quantum dynamics in the presence of strong system-reservoir coupling or structured (non-Markovian) environments exhibit behaviors inaccessible within simple Markovian or weak-coupling approximations.

• Non-Markovian Master Equations: The exact reduced dynamics can be written in terms of time-nonlocal memory kernels, with coefficients dependent on system-bath spectral properties. Strong memory effects support the emergence of non-thermal steady states ("qumemory" and oscillating qumemory) retaining initial state information, in contrast to conventional relaxation to equilibrium (Xiong et al., 2013).
• Transport and Fluctuations: In transport problems (e.g., energetic conductance between baths), non-Markovian and strong-coupling regimes cause substantial deviations from fluctuation-dissipation relations and lead to rich transient behavior. Hierarchical equations of motion (HEOM) enable accurate computation of time-dependent cumulants and transient fluctuation theorems beyond the steady-state (Cerrillo et al., 2016).

6. Spectral Statistics and Dynamical Signatures in Many-Body Systems

Spectral statistics—specifically the statistics of energy level spacings—play a critical role in the temporal behavior of many-body systems. Quantum-chaotic (Wigner-Dyson) level statistics imply rapid dephasing, characteristic Gaussian or Bessel-like decay of survival probabilities, and the presence of "correlation holes" in return amplitude, reflecting spectral rigidity. Integrable or many-body localized systems display Poissonian level statistics, power-law relaxation, and anomalously slow entanglement growth, challenging the reach of equilibrium quantum statistical mechanics (Santos et al., 2018, Zvyagin, 2017, Altman, 2015).

The identification of dynamical quantum phase transitions (DQPTs), signaled by nonanalyticities in Loschmidt rate functions, connects non-equilibrium dynamics with equilibrium criticality and topological features. DQPTs can be analyzed via real-time analogs of Fisher zero crossings, with scaling exponents and universality classes determined by the underlying quantum statistical properties (Zvyagin, 2017).

7. Experimental Access and Programmable Quantum Platforms

Recent advances enable direct experimental interrogation of quantum non-equilibrium statistics:

• Photonic Quantum Processors: Precise engineering of multi-mode photonic networks encodes truncated many-body propagators, directly reconstructing work distributions $Z$1, verifying fluctuation relations, and realizing genuine multi-particle bosonic interference in a programmable and scalable architecture (Krishna et al., 11 Mar 2026).
• Quantum Simulators and Hybrid Algorithms: Quantum-classical hybrid schemes implement nonequilibrium DMFT protocols for correlated electron systems, simulating Keldysh Green's function dynamics, and enabling the characterization of quantum statistics deep into the thermodynamically large, strongly-interacting regime (Kreula et al., 2015).

These platforms facilitate systematic study of how quantum indistinguishability, interference, protocol geometry, and environmental structure interlock to determine work, entropy production, thermalization pathways, and the boundaries of applicability for equilibrium statistical mechanics.

---

References:

• "Experimental simulation of non-equilibrium quantum piston on a programmable photonic quantum computer" (Krishna et al., 11 Mar 2026)
• "Aspects of non-equilibrium in classical and quantum systems: slow relaxation and glasses, dynamical large deviations, quantum non-ergodicity, and open quantum dynamics" (Garrahan, 2017)
• "Non-Markovian Dynamics Impact on the Foundations of Statistical Mechanics" (Xiong et al., 2013)
• "Non-Equilibrium Dynamics of a Noisy Quantum Ising Chain: statistics of the work and prethermalization after a sudden quench of the transverse field" (Marino et al., 2013)
• "Non-equilibrium quantum transport coefficients and the transient dynamics of full counting statistics in the strong coupling and non-Markovian regimes" (Cerrillo et al., 2016)
• "Nonequilibrium many-body quantum dynamics: from full random matrices to real systems" (Santos et al., 2018)
• "Full counting statistics approach to the quantum non-equilibrium Landauer bound" (Guarnieri et al., 2017)
• "Dynamical quantum phase transitions" (Zvyagin, 2017)
• "Non equilibrium quantum dynamics in ultra-cold quantum gases" (Altman, 2015)
• "Non-linear quantum-classical scheme to simulate non-equilibrium strongly correlated fermionic many-body dynamics" (Kreula et al., 2015)