Power-law distributions in nonequilibrium open quantum systems

1. Problem Statement

Power-law probability distributions are ubiquitous in complex classical systems (e.g., earthquakes, stock market crashes), often associated with "black swan" events—low-frequency, high-impact occurrences. While extensively studied in classical stochastic systems, the emergence of such heavy-tailed statistics in open quantum systems has remained largely unexplored.

The central problem addressed is whether a single dissipative quantum system, driven out of equilibrium, can naturally exhibit power-law tails in its steady-state energy distribution without fine-tuning. Specifically, the paper investigates the physical mechanism behind this phenomenon and whether it is a purely classical effect (multiplicative noise) or an intrinsic quantum feature.

2. Methodology

The author employs a combination of analytical derivation and numerical simulation using two primary models:

• The M-boson Model (Analytical Approach):

  • A minimal toy model consisting of a 1D quantum harmonic oscillator interacting with $M$ independent thermal baths.
  • The system exchanges $m$ excitation quanta with the $m$-th bath at dissipation rate $\gamma_m$ and pumping rate $\kappa_m$.
  • The dynamics are governed by a Lindblad master equation. The author derives the steady-state population distribution $\rho_{n,n}$ in the excitation number basis.
  • The analysis utilizes the quantum Itô equation in the Heisenberg picture to explicitly identify the nature of quantum noise (additive vs. multiplicative) arising from the dissipation terms.

• Quantum Liénard Systems (Numerical Approach):

  • To test the generality of the findings beyond the specific M-boson toy model, the author quantizes a general class of classical nonlinear dynamical systems known as Liénard systems (which include van der Pol, Rayleigh, and Duffing oscillators).
  • Using a quantization method that ensures the resulting master equation generates completely positive and trace-preserving (CPTP) maps, the author simulates the steady-state density matrix.
  • Numerical simulations are performed using the QuTiP package with Hilbert space truncation to observe the scaling of both populations (diagonal elements) and coherences (off-diagonal elements).

3. Key Contributions

A. Identification of the Mechanism: Multiplicative Quantum Noise

The paper establishes a general mechanism for power-law emergence in open quantum systems:

1. Nonlinear Dissipation: When a quantum system undergoes nonlinear dissipation (e.g., $m$-photon loss/gain where $m > 1$), the associated quantum noise is multiplicative.
2. Self-Amplification: Unlike additive noise, multiplicative noise is coupled to the system's state (e.g., proportional to $\hat{a}^{\dagger m-1}$). This causes fluctuations to self-amplify as the excitation energy increases.
3. Competition: The power-law tail arises from the competition between this noise amplification and the nonlinear dissipation. At high excitation numbers, the damping forces become insufficient to suppress the amplified fluctuations, leading to a heavy tail.
4. Quantum Constraint: Crucially, the form of this noise is not arbitrary; it is enforced by the constraints of quantum mechanics (specifically, the preservation of bosonic commutation relations $[\hat{a}, \hat{a}^\dagger] = 1$).

B. Analytical Proof in the M-boson Model

The author analytically proves that for the M-boson model, the steady-state energy distribution follows a truncated power law:
$$\rho_{n,n} \propto \frac{e^{-n\delta/M}}{n^\nu}$$
where $\delta = 1 - \kappa_M/\gamma_M$ represents the deviation from a critical point (infinite temperature bath).

• Critical Point ($\delta = 0$): The distribution becomes a pure power law $\rho_{n,n} \sim n^{-\nu}$.
• Exponent $\nu$: The exponent is determined by the rates of the $(M-1)$-boson processes: $\nu = \frac{M-1}{M^2} \frac{\gamma_{M-1} - \kappa_{M-1}}{\gamma_M}$.
• Extreme Events: In the regime where $\nu \le 2$, the mean excitation number diverges, leading to "extreme events" where the system occasionally occupies states with orders of magnitude more energy than the typical value.

C. Discovery of Power Laws in Coherences

A significant finding is that power-law tails appear not only in the populations (diagonal elements $\rho_{n,n}$) but also in the coherences (off-diagonal elements $\rho_{n,m}$) of the density matrix in the quantum Liénard systems. This indicates that the heavy-tailed statistics are an intrinsic feature of the quantum state, not just a classical probability distribution artifact.

4. Key Results

• Robustness: Unlike classical systems where power laws often require fine-tuning, the quantum Liénard systems exhibit power-law tails for typical parameters in the classically stable regime (limit cycle regime).
• Infinite Mean Excitation: The simulations show regimes where the steady state is normalizable but possesses an infinite mean excitation number, a hallmark of heavy-tailed distributions.
• Physical Signatures:
  • Photon Superbunching: The second-order correlation function $g^{(2)}(0)$ can become arbitrarily large, indicating extreme photon bunching far exceeding thermal light sources.
  • Phase Space: The Wigner quasiprobability distribution exhibits a heavy tail $W(q,p) \sim (q^2+p^2)^{-\nu}$, implying power-law decay in position and momentum distributions.
• Comparison with Classical: Even when the corresponding classical system is stable and does not exhibit power laws, the quantum system does, highlighting the essential role of quantum noise.

5. Significance and Implications

• Fundamental Physics: The work bridges the gap between complex systems theory and quantum optics, demonstrating that extreme statistical events are a natural consequence of nonlinear quantum dissipation. It challenges the standard view that quantum fluctuations become negligible at large amplitudes in dissipative systems.
• Quantum Engineering: The phenomenon offers a new route to generate extreme photon sources. Systems exhibiting these power laws can produce single pulses containing $\sim 10^6$ photons (as observed in related experiments), which is useful for:
  • Sensing: Applications like ghost imaging that rely on high-order correlations.
  • Light-Matter Interaction: Exploring extreme regimes of interaction where standard perturbative methods fail.
• Quantum Error Correction: The ability to engineer nonlinear dissipation (reservoir engineering) to stabilize specific states (like Schrödinger cat states) is directly linked to the mechanisms described, suggesting new pathways for fault-tolerant quantum computing.

In summary, the paper provides a rigorous theoretical framework showing that nonlinear dissipation in open quantum systems inevitably generates multiplicative quantum noise, which drives the system toward a nonequilibrium steady state characterized by power-law energy distributions and extreme statistical events.