skip to main content


Title: Universality of light thermalization in multimoded nonlinear optical systems
Abstract

Recent experimental studies in heavily multimoded nonlinear optical systems have demonstrated that the optical power evolves towards a Rayleigh–Jeans (RJ) equilibrium state. To interpret these results, the notion of wave turbulence founded on four-wave mixing models has been invoked. Quite recently, a different paradigm for dealing with this class of problems has emerged based on thermodynamic principles. In this formalism, the RJ distribution arises solely because of ergodicity. This suggests that the RJ distribution has a more general origin than was earlier thought. Here, we verify this universality hypothesis by investigating various nonlinear light-matter coupling effects in physically accessible multimode platforms. In all cases, we find that the system evolves towards a RJ equilibrium—even when the wave-mixing paradigm completely fails. These observations, not only support a thermodynamic/probabilistic interpretation of these results, but also provide the foundations to expand this thermodynamic formalism along other major disciplines in physics.

 
more » « less
Award ID(s):
1711230
NSF-PAR ID:
10392446
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Nature Publishing Group
Date Published:
Journal Name:
Nature Communications
Volume:
14
Issue:
1
ISSN:
2041-1723
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. We show that in highly multimoded nonlinear photonic systems, the optical thermodynamic pressures emerging from different species of the optical field obey Dalton’s law of partial pressures. In multimode settings, the optical thermodynamic pressure is defined as the conjugate to the extensive variable associated with the system’s total number of modes and is directly related to the actual electrodynamic radiation forces exerted at the physical boundaries of the system. Here, we extend this notion to photonic configuration supporting different species of the optical field. Under thermal equilibrium conditions, we formally derive an equation that establishes a direct link between the partial thermodynamic pressures and the electrodynamic radiation pressures exerted by each polarization species. Our theoretical framework provides a straightforward approach for quantifying the total radiation pressures through the system’s thermodynamic variables without invoking the Maxwell stress tensor formalism. In essence, we show that the total electrodynamic pressure in such arrangements can be obtained in an effortless manner from initial excitation conditions, thus avoiding time-consuming simulations of the utterly complex multimode dynamics. To illustrate the validity of our results, we carry out numerical simulations in multimoded nonlinear optical structures supporting two polarization species and demonstrate excellent agreement with the Maxwell stress tensor method.

     
    more » « less
  2. Abstract

    Cosmic rays produced by young stellar objects can potentially alter the ionization structure, heating budget, chemical composition, and accretion activity in circumstellar disks. The inner edges of these disks are truncated by strong magnetic fields, which can reconnect and produce flaring activity that accelerates cosmic radiation. The resulting cosmic rays can provide a source of ionization and produce spallation reactions that alter the composition of planetesimals. These reconnection and particle acceleration processes are analogous to the physical processes that produce flaring in and the heating of stellar coronae. Flaring events on the surface of the Sun exhibit a power-law distribution of energy, reminiscent of those measured for earthquakes and avalanches. Numerical lattice reconnection models are capable of reproducing the observed power-law behavior of solar flares under the paradigm of self-organized criticality. One interpretation of these experiments is that the solar corona maintains a nonlinear attractor—or “critical”—state by balancing energy input via braided magnetic fields and output via reconnection events. Motivated by these results, we generalize the lattice reconnection formalism for applications in the truncation region of magnetized disks. Our numerical experiments demonstrate that these nonlinear dynamical systems are capable of both attaining and maintaining criticality in the presence of Keplerian shear and other complications. The resulting power-law spectrum of flare energies in the equilibrium attractor state is found to be nearly universal in magnetized disks. This finding indicates that magnetic reconnection and flaring in the inner regions of circumstellar disks occur in a manner similar to the activity on stellar surfaces. These results, in turn, have ramifications for the spallation-driven injection of radionuclides in planetesimals, disk ionization, and the subsequent planetary formation process.

     
    more » « less
  3. Abstract: There is no theoretical underpinning that successfully explains how turbulent mixing is fed by wave breaking associated with nonlinear wave-wave interactions in the background oceanic internal wavefield. We address this conundrum using one-dimensional ray tracing simulations to investigate interactions between high frequency internal waves and inertial oscillations in the extreme scale separated limit known as “Induced Diffusion”. Here, estimates of phase locking are used to define a resonant process (a resonant well) and a non-resonant process that results in stochastic jumps. The small amplitude limit consists of jumps that are small compared to the scale of the resonant well. The ray tracing simulations are used to estimate the first and second moments of a wave packet’s vertical wavenumber as it evolves from an initial condition. These moments are compared with predictions obtained from the diffusive approximation to a self-consistent kinetic equation derived in the ‘Direct Interaction Approximation’. Results indicate that the first and second moments of the two systems evolve in a nearly identical manner when the inertial field has amplitudes an order of magnitude smaller than oceanic values. At realistic (oceanic) amplitudes, though, the second moment estimated from the ray tracing simulations is inhibited. The transition is explained by the stochastic jumps obtaining the characteristic size of the resonant well. We interpret this transition as an adiabatic ‘saturation’ process which changes the nominal background wavefield from supporting no mixing to the point where that background wavefield defines the normalization for oceanic mixing models. 
    more » « less
  4. Abstract

    Large-amplitude internal solitary waves were recently observed in a coastal plain estuary and were hypothesized to evolve from an internal lee wave generated at the channel–shoal interface. To test this mechanism, a 3D nonhydrostatic model with nested domains and adaptive grids was used to investigate the generation of the internal solitary waves and their subsequent nonlinear evolution. A complex sequence of wave propagation and transformation was documented and interpreted using the nonlinear wave theory based on the Korteweg–de Vries equation. During the ebb tide a mode-2 internal lee wave is generated by the interaction between lateral flows and channel–shoal topography. This mode-2 lee wave subsequently propagates onto the shallow shoal and transforms into a mode-1 wave of elevation as strong mixing on the flood tide erases stratification in the bottom boundary layer and the lower branch of the mode-2 wave. The mode-1 wave of elevation evolves into an internal solitary wave due to nonlinear steepening and spatial changes in the wave phase speed. As the solitary wave of elevation continues to propagate over the shoaling bottom, the leading edge moves ahead as a rarefaction wave while the trailing edge steepens and disintegrates into a train of rank-ordered internal solitary waves, due to the combined effects of shoaling and dispersion. Strong turbulence in the bottom boundary layer dissipates wave energy and causes the eventual destruction of the solitary waves. In the meantime, the internal solitary waves can generate elevated shear and dissipation rate in local regions.

    Significance Statement

    In the coastal ocean nonlinear internal solitary waves are widely recognized to play an important role in generating turbulent mixing, modulating short-term variability of nearshore ecosystem, and transporting sediment and biochemical materials. However, their effects on shallow and stratified estuaries are poorly known and have been rarely studied. The nonhydrostatic model simulations presented in this paper shed new light into the generation, propagation, and transformation of the internal solitary waves in a coastal plain estuary.

     
    more » « less
  5. By utilizing notions from statistical mechanics, we develop a general and self-consistent theoretical framework capable of describing any weakly nonlinear optical multimode system involving conserved quantities. We derive the fundamental relations that govern the grand canonical ensemble through maximization of the Gibbs entropy at equilibrium. In this classical picture of statistical photo-mechanics, we obtain analytical expressions for the probability distribution, the grand partition function, and the relevant thermodynamic potentials. Our results universally apply to any other weakly nonlinear multimode bosonic system.

     
    more » « less