The origin of KPZ-scaling in arrays of polariton condensates

This paper establishes that the Kardar-Parisi-Zhang (KPZ) scaling observed in the phase dynamics of polariton condensate arrays originates from fluctuations in Goldstone modes resulting from spontaneous $U(1)$ symmetry breaking, thereby linking microscopic system parameters to the coherent properties of emitted light.

The Gist

Imagine a vast, crowded dance floor where thousands of tiny particles (called polaritons) are dancing. These aren't just any particles; they are a hybrid of light and matter, trapped inside a semiconductor chip. When you shine a laser on them, they start to synchronize, moving in perfect unison like a single, giant wave. This is called a condensate.

For a long time, scientists noticed something strange about how these dancers move. Even though they are trying to dance in perfect sync, there are tiny, random wobbles and jitters in their rhythm. When they measured how "in step" the dancers were over time and distance, the pattern of these wobbles followed a very specific, universal rule known as KPZ scaling.

Think of KPZ scaling like the way a pile of sand grows. If you pour sand onto a pile, the surface gets rougher and rougher in a predictable way, no matter if the sand is fine or coarse. Scientists found this same "roughness" pattern in the quantum dance of polaritons, but they didn't know why it happened. They knew the math worked, but they didn't understand the microscopic cause.

The "Ghost Dancers" (Goldstone Modes)

This new paper solves the mystery by introducing a cast of "ghost dancers."

In physics, when a system breaks symmetry (like when the dancers suddenly decide to all face North instead of spinning randomly), it creates massless ripples called Nambu-Goldstone modes. You can think of these as the "ghosts" of the dance floor. They are low-energy excitations that don't have a fixed position but float around, causing tiny fluctuations in the phase (the timing) of the main condensate.

The authors argue that the KPZ scaling isn't caused by the main dancers themselves, but by the fluctuations of these ghost dancers.

The Analogy: The Orchestra and the Background Noise

Imagine a massive orchestra playing a single, perfect note (the condensate).

• The Main Note: This is the strong, coherent light the system emits.
• The Ghosts (Goldstone Modes): These are the tiny, random whispers and hiccups in the air around the orchestra.

In the past, scientists thought the "KPZ roughness" was a fundamental property of the orchestra itself. This paper says: No, it's the background noise.

The authors created a mathematical model showing that when the "ghosts" (the excited states) are numerous and active, their random jitters interfere with the main note. This interference creates a specific type of "roughness" in the timing of the music.

• In a 1D line (a single row of dancers): The wobbles follow a specific rhythm (exponent 1/3).
• In a 2D grid (a full dance floor): The wobbles follow a slightly different rhythm (exponent ~0.24).

These rhythms match the KPZ predictions perfectly.

The Volume Knob (Pump Power)

The paper also explains why this effect is only seen at certain times. Imagine the laser pumping energy into the system is a volume knob.

• Low Volume (Low Power): The main orchestra is playing softly. The "ghosts" are loud and active relative to the main note. Their random jitters dominate, and you see the KPZ pattern clearly.
• High Volume (High Power): You crank the volume up. The main orchestra becomes so loud and powerful that the "ghosts" are drowned out. The system becomes too stable, and the KPZ pattern disappears, replaced by a different, more orderly type of behavior.

Why Does This Matter?

This discovery is like finding the "source code" for a complex video game. Instead of just observing that the graphics look a certain way, the authors found the specific line of code (the fluctuation of Goldstone modes) that generates the effect.

The Takeaway:
The strange, universal "roughness" seen in these quantum fluids isn't a mystery of the universe's deep laws, but a direct result of how the "ghosts" of the system interact with the main wave. This gives scientists a new tool: by controlling the population of these ghost modes (by adjusting the laser power or the material properties), they can tune the quality of the light emitted by these condensates. This could lead to better quantum light sources for future technologies, allowing us to design light with very specific, tailored properties.

Technical Summary

1. Problem Statement

The Kardar–Parisi–Zhang (KPZ) universality class, originally formulated to describe the stochastic growth of interfaces, has recently been observed in the phase dynamics of non-equilibrium quantum fluids, specifically exciton-polariton condensates. While previous studies established a connection between phase dynamics and KPZ scaling using phenomenological driven-dissipative Gross-Pitaevskii equations (GPE) with added white noise, a fundamental microscopic explanation remained lacking.

Specifically, the paper addresses the need to derive the characteristic power-law decay of the first-order coherence function, $g^{(1)}$, from first principles. The goal is to link the emergent KPZ scaling and its critical exponents directly to the microscopic parameters of polariton systems (interaction strength, energy dispersion, effective temperature) rather than treating noise as an external, phenomenological input.

2. Methodology

The authors propose a simple analytical model based on the stochastic dynamics of Nambu-Goldstone modes arising from spontaneous $U(1)$ symmetry breaking.

• Theoretical Framework:

  • The system is modeled as an infinite periodic lattice of exciton-polariton condensates described by a scalar space-dependent order parameter $\psi(r,t)$.
  • The order parameter is expanded using a Bogolyubov ansatz, separating the macroscopically occupied ground state (condensate) from the fluctuating excited states (Nambu-Goldstone modes).
  • The population of excited states is governed by the Bose-Einstein distribution at an effective temperature $T$.
  • The Nambu-Goldstone modes obey a conventional Bogolyubov dispersion relation: $\varepsilon(k) = \sqrt{E(k)(E(k) + 2\mu)}$, where $\mu$ is the effective chemical potential derived from polariton-polariton interactions.

• Calculation of Correlation:

  • The first-order correlation function $g^{(1)}(\Delta r, \Delta t)$ is calculated by performing ensemble averaging over the random phases of the harmonics in the Bogolyubov expansion.
  • Due to the finite size of experimental arrays ($L$), the $k$-space is discretized, and the summation is performed over the first Brillouin zone.

• Numerical Simulations:

  • The model was applied to both one-dimensional (1D) linear chains and two-dimensional (2D) triangular lattices of polariton condensates.
  • Specific energy band structures were assumed for the lattices (e.g., cosine-based dispersion for 1D chains and geometric structure factors for 2D triangular arrays).
  • Simulations varied key parameters, including the interaction energy $\mu$ (proxy for pump power) and system dimensions, to observe the transition between scaling regimes.

3. Key Contributions

• Microscopic Origin of KPZ Scaling: The paper identifies the fluctuation of the population of Nambu-Goldstone modes as the primary mechanism driving KPZ scaling. It demonstrates that the critical exponents are a direct consequence of the dynamics of these massless excitations, rather than an artifact of phenomenological noise terms.
• Derivation from First Principles: The authors provide a derivation linking microscopic parameters (interaction strength, dispersion, effective temperature) to the emergent macroscopic scaling laws, bridging the gap between microscopic quantum fluid dynamics and universal statistical physics.
• Pump Power Dependence: The study elucidates the specific regime where KPZ scaling occurs. It establishes that KPZ dynamics dominate only when the integrated population of Nambu-Goldstone modes is comparable to the condensate population (typically near the bosonic condensation threshold).

4. Key Results

• Critical Exponents:
  • 1D Case: The simulations reproduced the temporal and spatial critical exponents $\beta = 1/3$ and $\chi = 1/2$, respectively, matching the 1D KPZ universality class.
  • 2D Case: The model yielded $\beta \approx 0.241$ and $\chi \approx 0.387$, consistent with the 2D KPZ universality class and recent experimental observations (e.g., Widmann et al., 2025).
• Role of Pump Power:
  • As pump power increases, the condensate density rises, and the effective chemical potential $\mu$ increases.
  • This leads to a steeper slope in the Nambu-Goldstone dispersion and a reduction in the population of these excited modes.
  • Consequently, the time interval over which KPZ scaling prevails shrinks. At high pump powers, the system transitions away from KPZ scaling toward equilibrium Berezinskii-Kosterlitz-Thouless (BKT) dynamics.
• Finite Size Effects: The model explains the oscillations and limited range of KPZ regions observed in simulations as artifacts of Brillouin zone discretization due to finite system size, noting that these are often unresolved in experiments.

5. Significance

• Theoretical Validation: The work resolves a long-standing question regarding the microscopic origin of KPZ universality in driven-dissipative quantum fluids, moving beyond phenomenological descriptions to a mechanism rooted in symmetry breaking and Goldstone modes.
• Experimental Guidance: It provides a clear criterion for experimentalists: to observe KPZ scaling, one must operate near the condensation threshold where excited state populations are significant. High-power regimes will suppress these fluctuations, masking the KPZ behavior.
• Technological Application: By establishing a direct link between microscopic parameters and coherent light properties, the paper offers a pathway for controlling scaling behaviors in bosonic condensate arrays. This is crucial for developing tailored quantum light sources with specific photon correlation properties.
• Generalizability: The mechanism proposed (fluctuating Goldstone modes) suggests that KPZ scaling may be a universal feature of any system containing bosonic condensates where symmetry is spontaneously broken, provided the excited state populations are sufficiently populated.