Three-Dimensional Kardar--Parisi--Zhang Scaling in Polariton Condensates

Technical Summary: Three-Dimensional Kardar–Parisi–Zhang Scaling in Polariton Condensates

Problem Statement
The Kardar–Parisi–Zhang (KPZ) universality class describes macroscopic scaling arising from microscopic violations of detailed balance in driven-dissipative systems. While KPZ scaling has been extensively explored in one- and two-dimensional driven condensates and growing interfaces, demonstrating this scaling in three spatial dimensions ($d=3$) remains a significant challenge. In $d=3$, the KPZ equation exhibits a weak-coupling smooth phase and a strong-coupling rough phase; experimental realization requires driving the system into the latter regime. Furthermore, the growth exponent ($\beta \approx 0.18$) is small, necessitating a broad scaling window to observe coherence decay, while finite-size effects and multiple relaxation channels complicate the observation of asymptotic behavior. Existing platforms, such as nanowires or planar microcavities, are inherently one- or two-dimensional, failing to support the necessary three-dimensional Bloch dispersion.

Methodology
The authors propose a three-dimensional exciton-polariton crystal formed in an inverse-opal photonic crystal doped with excitonic emitters (e.g., semiconductor quantum dots) as a platform for observing 3D KPZ universality.

1. Theoretical Derivation:

   • The study begins with a stochastic driven-dissipative Gross–Pitaevskii equation (GPE) describing a polariton condensate in a 3D photonic crystal lower-polariton band. The model includes anisotropic effective masses ($m_i$), a reservoir density ($n_R$), and noise terms representing pumping and decay.
   • By eliminating the massive condensate-density and reservoir-density fluctuations in a local adiabatic approximation, the authors derive an effective equation for the condensate phase ($\theta$).
   • This reduction yields an effective anisotropic KPZ equation: $\partial_t \theta = \sum_i \nu_i \partial_i^2 \theta + \frac{1}{2} \sum_i \lambda_i (\partial_i \theta)^2 + \eta_{\text{eff}}$. Under conditions of an approximately isotropic valley, this reduces to the standard 3+1-dimensional KPZ equation.
   • The coefficients ($\nu, \lambda$) are shown to be tunable via pump distance from threshold, effective mass, noise strength, and decay rates.

2. Numerical Simulations:

   • Control Calculation: The authors first simulate the standard 3+1 KPZ equation (Eq. 10) to establish benchmark scaling exponents ($\beta \approx 0.1845$, $\chi \approx 0.3135$) and dynamic exponent $z \approx 1.699$.
   • Full Model Simulation: They then simulate the full stochastic open-dissipative polariton model (Eqs. 2 and 3). To mimic stochastic growth and avoid the suppression of phase fluctuations seen in uniform pumping, they introduce a weakly disordered extended pump $P(r,t)$.
   • Observables: The study focuses on the first-order correlation function $g^{(1)}(\Delta r, \Delta t)$. In the phase-dominated regime, $-\ln |g^{(1)}|$ is proportional to the mean-square phase fluctuation. The authors analyze temporal coherence ($-\ln |g^{(1)}(0, \Delta t)| \propto |\Delta t|^{2\beta}$) and spatial coherence ($-\ln |g^{(1)}(\Delta r, 0)| \propto |\Delta r|^{2\chi}$).

Key Results

• Benchmark Validation: Numerical simulations of the pure KPZ equation confirm the expected 3+1 scaling. The temporal exponent derived from compact coherence is $\approx 0.375$ (corresponding to $2\beta \approx 0.369$), and spatial correlations show movement toward the benchmark $2\chi \approx 0.627$. Spatiotemporal correlators exhibit data collapse consistent with the Family-Vicsek scaling form.
• Polariton Model Performance: Simulations of the full driven-dissipative polariton model with a disordered pump reveal an intermediate-asymptotic regime where the system exhibits 3D KPZ scaling.
  • The temporal coherence slope is found to be $0.367$ (close to the benchmark $0.369$).
  • The spatial coherence slope is $0.618$ (close to the benchmark $0.627$).
  • The spatiotemporal correlation function $C(\Delta r, \Delta t)$ collapses onto a universal curve when plotted in KPZ variables ($C/\Delta t^{2\beta}$ vs. $\Delta r/\Delta t^{1/z}$), supporting the dynamic scaling hypothesis.
• Spectral Characteristics: The analysis of the 3D Bogoliubov dispersion indicates that, unlike in 2D, the 3D Goldstone mode does not produce a low-energy white spectral plateau, with the population spectrum scaling linearly with energy ($S_{\text{pop}}(E) \propto E$) in the classical low-energy limit.

Significance and Claims
The paper claims to identify three-dimensional polariton crystals as a controllable quantum fluid route to higher-dimensional nonequilibrium universality. The proposed platform addresses the primary obstacles to observing 3D KPZ scaling:

1. It provides a genuine three-dimensional Bloch dispersion, avoiding the reduction of phase dynamics to lower dimensions.
2. The band structure allows for the engineering of effective mass tensors and the tuning of KPZ coupling parameters without altering dimensionality.
3. The driven-dissipative nature of the system is intrinsic to the generation of KPZ phase dynamics, rather than being a parasitic perturbation.

The authors assert that their results provide a concrete route to observing nonequilibrium universality beyond two dimensions. They propose that the predicted scaling can be measured via momentum-resolved photoluminescence and interferometric measurements of the emitted field to reconstruct the full three-dimensional correlation function. The work suggests that 3D polariton crystals serve as a unique platform where photonic band engineering and stochastic many-body scaling can be studied within a single experimental setup.