Dynamical thermalization, Rayleigh-Jeans condensate, vortexes and wave collapse in quantum chaos fibers and fluid of light

This paper investigates the time evolution of nonlinear fields in chaotic D-shaped billiards, revealing that strong nonlinearity drives dynamical thermalization into a Rayleigh-Jeans condensate, while also characterizing phenomena such as wave collapse, vortex dynamics, and superfluidity in both focusing and defocusing regimes relevant to optical fibers and fluid light.

The Gist

Imagine a long, twisted glass fiber, shaped like a "D" (a circle with one flat side cut off). Inside this fiber, we send a beam of light. Usually, light travels in straight lines, but inside this specific fiber, the walls are shaped so that the light bounces around in a wildly unpredictable, chaotic way, much like a pinball in a machine with no flippers.

This paper explores what happens when we turn up the "volume" on the light's own personality. In physics terms, this is called nonlinearity. When the light is intense enough, it starts to interact with itself, changing how it moves. The researchers wanted to see: Does this chaotic, self-interacting light eventually settle down into a predictable pattern, or does it stay wild forever?

Here is the story of their findings, broken down into simple concepts:

1. The Two Worlds: Order vs. Chaos

Think of the light inside the fiber as a crowd of people in a room.

• The Quiet Room (Low Nonlinearity): If the light is weak, the "people" (light waves) don't really talk to each other. They just bounce around the walls. If the room is shaped just right, they might get stuck in specific patterns, never mixing with everyone else. This is called being "quasi-integrable." It's like a dance where everyone stays in their own lane.
• The Wild Party (High Nonlinearity): If the light is strong, the waves start bumping into each other, pushing and shoving. This creates chaos. The researchers found that once the chaos gets strong enough (crossing a specific "chaos border"), the system stops being a dance and becomes a mosh pit. But here is the surprise: this mosh pit eventually settles down into a very specific, organized state.

2. The Great Migration: The "Rayleigh-Jeans" Condensate

When the chaos settles, something magical happens. Almost all the energy (about 80% to 90% of it) decides to move to the very lowest, calmest spot in the fiber—the "ground state."

Imagine a crowded stadium where everyone is running around wildly. Suddenly, without any external force telling them to, 90% of the crowd spontaneously sits down in the front row, leaving the rest of the stadium almost empty. The researchers call this a Rayleigh-Jeans Condensate.

• Why is this special? In the quantum world (like cold atoms), you expect things to spread out or behave differently. But here, because the light acts like a "classical fluid" (a wave of water rather than tiny particles), it follows different rules. It piles up in the lowest energy state, creating a super-dense, calm core of light.

3. The "Fröhlich" Confusion

The paper makes a clear distinction between this new discovery and an old idea called the "Fröhlich condensate."

• The Old Idea (Fröhlich): Imagine a machine that keeps pumping energy into a system while also draining it away (like a leaky bucket being filled). In this scenario, energy can pile up at high temperatures.
• The New Discovery (Rayleigh-Jeans): The fiber in this experiment is a closed system. No energy is being pumped in or drained out. It's a self-contained universe. The light piles up in the lowest state only when the system is "cool" (low energy relative to the number of modes). It's a spontaneous gathering, not a forced one.

4. The "Collapse" and the "Vortex"

The researchers also looked at what happens if the light tries to focus too hard or if it spins.

• The Collapse: If the light tries to focus too intensely (like a magnifying glass focusing sunlight), it can theoretically "collapse" into a single point of infinite density. In an open field, this is a known danger. But inside this chaotic "D-shaped" fiber, the chaos actually fights the collapse, creating a weird, unstable dance between the two forces.
• The Vortex: When the light is defocused (spread out), it can form swirling patterns, like water going down a drain. The researchers found that even in this chaotic fiber, these swirls (vortices) can survive for a long time, acting like tiny, stable tornadoes of light.

5. The Entropy Puzzle (The "Messiness" Meter)

In physics, "entropy" is a measure of messiness. Usually, when things settle down, they get messier (entropy goes up).

• The Twist: The researchers tracked a specific type of "messiness" called quantum entropy. They found it went up (the system got messy as waves mixed), reached a peak, and then went down as the system settled into the condensate.
• The Analogy: Imagine a messy room where you throw everything in the air (entropy goes up). Then, instead of leaving it messy, everyone suddenly agrees to put everything back in perfect, neat piles (entropy goes down). The system found a new kind of order that is very different from the initial chaos.

Summary

The paper proves that in a chaotic, D-shaped optical fiber, strong light waves don't just stay chaotic. They undergo a process of dynamical thermalization. They shake themselves out, and then, surprisingly, they all migrate to the lowest energy state, forming a massive, stable "condensate" of light.

This isn't just a theory; the math and computer simulations show that this happens naturally in these fibers. It suggests that we can use these "quantum chaos fibers" to study how complex systems organize themselves, potentially leading to new ways to control light in telecommunications or understanding how fluids behave at a microscopic level.

In short: Chaos leads to a specific kind of order where almost all the light gathers in one spot, creating a stable, calm core in the middle of a wild, swirling storm.

Technical Summary

Technical Summary: Dynamical Thermalization, Rayleigh-Jeans Condensate, Vortexes, and Wave Collapse in Quantum Chaos Fibers and Fluid of Light

Problem Statement
The paper investigates the time evolution of a nonlinear field governed by the Nonlinear Schrödinger Equation (NSE) within a chaotic D-shaped billiard geometry. This model serves a dual physical purpose: describing longitudinal light propagation in multimode optical fibers with a chaotic cross-section and modeling the dynamics of a fluid of light in a Kerr nonlinear medium (atomic vapor). The central problem addresses the mechanism of dynamical thermalization in Hamiltonian systems without external thermal baths. Specifically, the authors aim to distinguish between regimes of quasi-integrable dynamics (governed by Kolmogorov-Arnold-Moser (KAM) theory at weak nonlinearity) and chaotic thermalization (occurring above a specific "chaos border"). A key objective is to determine whether the resulting thermal distribution follows the quantum Bose-Einstein (BE) statistics or the classical Rayleigh-Jeans (RJ) statistics, and to analyze the formation of a condensate at the ground state under these conditions.

Methodology
The study employs a combination of analytical derivations and extensive numerical simulations:

• Model: The system is defined by the NSE with Dirichlet boundary conditions in a D-shaped billiard (a circle with a straight cut). The nonlinearity parameter $\beta$ controls the interaction strength, with $\beta > 0$ representing defocusing and $\beta < 0$ representing focusing nonlinearity.
• Numerical Integration: Two methods are utilized. For short-time dynamics ($t \le 40$), a Trotter decomposition method is used, conserving norm and energy with high precision. For long-time evolution ($t$ up to 1000), the authors employ the FREEFEM package using finite element methods with up to $1.36 \times 10^5$ mesh nodes to ensure stability and accuracy over extended periods.
• Initial Conditions: Simulations typically start with superpositions of linear eigenmodes (e.g., $(\psi_m + \psi_{m+1})/\sqrt{2}$) to break parity symmetry and induce chaotic mixing.
• Analytical Framework: The authors derive analytical estimates for the Rayleigh-Jeans condensate fraction in a simplified model with equidistant energy levels (RJS model) using microcanonical ensemble statistics and saddle-point approximations, comparing these with Monte Carlo simulations.

Key Contributions and Results

1. Dynamical Thermalization and the Chaos Border:
   The authors demonstrate that for nonlinearity strengths exceeding a specific chaos border, the system undergoes dynamical thermalization. Below this border, the dynamics remain quasi-integrable (KAM regime). Above the border, the system evolves toward a steady state described by the classical Rayleigh-Jeans thermal distribution:
   $$\rho_m = \frac{T}{E_m - \mu}$$
   where $T$ is temperature and $\mu$ is the chemical potential. This confirms that the thermalization is a result of purely Hamiltonian chaos, not external dissipation.

2. Rayleigh-Jeans (RJ) Condensate:
   A primary finding is the emergence of an RJ condensate, where approximately 80–90% of the total probability (or particle number) accumulates in the ground state ($m=1$).

   • Mechanism: This condensation occurs at low effective temperatures ($T \ll E_m$) in systems with a large number of modes ($N$). The fraction of non-condensed probability scales as $W_{ex} \propto (\ln N)/N$, vanishing as $N \to \infty$.
   • Distinction from Fröhlich Condensate: The paper explicitly distinguishes this RJ condensate from the Fröhlich condensate. The Fröhlich model requires high temperatures ($T \gg E_m$) and external pumping/dissipation. In contrast, the RJ condensate arises in a conservative Hamiltonian system at low temperatures without external pumping.
   • Distinction from Bose-Einstein (BE) Condensate: The numerical results show a drastic difference between BE and RJ distributions. While BE statistics would predict a ground state population of $\approx 0.38$ for the studied parameters, the RJ statistics predict $\approx 0.98$. The simulations align with the RJ prediction, confirming the classical nature of the field thermalization.

3. Entropy Evolution:
   The study tracks both classical Boltzmann entropy ($S_B$) and quantum von Neumann entropy ($S_q$).

   • $S_B(t)$ grows monotonically, consistent with the Boltzmann H-theorem.
   • $S_q(t)$ exhibits a non-monotonic behavior: it initially increases as energy spreads to higher modes, reaches a maximum, and then decreases significantly as the system settles into the RJ condensate state. This decrease is attributed to the high concentration of probability in the ground state, which reduces the effective number of populated modes.

4. Wave Collapse (Focusing Nonlinearity):
   For the focusing case ($\beta < 0$), the authors investigate wave collapse. Unlike the infinite-space case where collapse occurs at finite times, the finite size of the D-shaped billiard and the presence of chaos lead to an interplay between collapse and thermalization. The study shows that wave collapse can occur even at sufficiently high positive energies, a behavior distinct from open-space scenarios.

5. Vortex Dynamics (Defocusing Nonlinearity):
   In the defocusing regime ($\beta > 0$) with strong nonlinearity, the system exhibits a superfluid regime characterized by long-lived vortexes. This establishes the D-shaped fiber as a platform for studying vortex dynamics in a chaotic "fluid of light."

6. Experimental Relevance:
   The paper discusses parameters for realizing these effects in multimode optical fibers and atomic vapor systems. It notes that while previous experiments on "self-cleaning" in fibers may have occurred in degenerate mode spectra (where KAM theory fails), the D-shape fiber offers a generic case where the KAM regime is valid, allowing for the observation of the chaos border and subsequent RJ thermalization.

Significance
The paper claims to provide a fundamental understanding of dynamical thermalization in generic nonlinear oscillator systems. By utilizing the D-shape billiard, it isolates the effects of chaos from spectral degeneracy, demonstrating that thermalization leads to a classical Rayleigh-Jeans distribution rather than a quantum Bose-Einstein distribution. The identification of the RJ condensate as a low-temperature, conservative phenomenon distinct from the Fröhlich condensate clarifies the statistical mechanics of classical fields. Furthermore, the work bridges the gap between quantum chaos theory, nonlinear optics, and fluid dynamics, offering a theoretical basis for experimental investigations of thermalization, condensates, and vortex dynamics in multimode fibers and atomic vapors.