Statistical mechanics of beam self-cleaning in GRIN multimode optical fibers

Random matrix model of Kolmogorov-Zakharov turbulence

We introduce and study a random matrix model of Kolmogorov-Zakharov turbulence in a nonlinear purely dynamical finite-size system with many degrees of freedom. For the case of a direct cascade, the energy and norm pumping takes place at low energy scales with absorption at high energies. For a pumping strength above a certain chaos border, a global chaotic attractor appears with a stationary energy flow through a Hamiltonian inertial energy interval. In this regime, the steady-state norm distribution is described by an algebraic decay with an exponent in agreement with the Kolmogorov-Zakharov theory. Below the chaos border, the system is located in the quasi-integrable regime similar to the Kolmogorov-Arnold-Moser theory and the turbulence is suppressed. For the inverse cascade, the system rapidly enters a strongly nonlinear regime where the weak turbulence description is invalid. We argue that such a dynamical turbulence is generic, showing that it is present in other lattice models with disorder and Anderson localization. We point out that such dynamical models can be realized in multimode optical fibers.

Nature of Optical Thermodynamic Pressure Exerted in Highly Multimoded Nonlinear Systems

The theory of optical thermodynamics provides a comprehensive framework that enables a self-consistent description of the intricate dynamics of nonlinear multimoded photonic systems. This theory, among others, predicts a pressurelike intensive quantity (p[over ^]) that is conjugate to the system's total number of modes (M)-its corresponding extensive variable. Yet at this point, the nature of this intensive quantity is still nebulous. In this Letter, we elucidate the physical origin of the optical thermodynamic pressure and demonstrate its dual essence. In this context, we rigorously derive an expression that splits p[over ^] into two distinct components, a term that is explicitly tied to the electrodynamic radiation pressure and a second entropic part that is responsible for the entropy change. We utilize this result to establish a formalism that simplifies the quantification of radiation pressure under nonlinear equilibrium conditions, thus eliminating the need for a tedious evaluation of the Maxwell stress tensor. Our theoretical analysis is corroborated by numerical simulations carried out in highly multimoded nonlinear optical structures. These results may provide a novel way in predicting and controlling radiation pressure processes in a variety of nonlinear electromagnetic settings.

High-temperature wave thermalization spoils beam self-cleaning in nonlinear multimode GRIN fibers

In our experiments, we reveal a so-far unnoticed power limitation of beam self-cleaning in graded-index nonlinear multimode optical fibers. As the optical pulse power is progressively increased, we observed that the initial Kerr-induced improvement of the spatial beam quality is eventually lost. Based on a holographic mode decomposition of the output field, we show that beam spoiling is associated with high-temperature wave thermalization, which depletes the fundamental mode in favor of a highly multimode power distribution.

Obstruction to ergodicity in nonlinear Schrödinger equations with resonant potentials

We identify a class of trapping potentials in cubic nonlinear Schrödinger equations (NLSEs) that make them nonintegrable, but prevent the emergence of power spectra associated with ergodicity. The potentials are characterized by equidistant energy spectra (e.g., the harmonic-oscillator trap), which give rise to a large number of resonances enhancing the nonlinearity. In a broad range of dynamical solutions, spanning the regimes in which the nonlinearity may be either weak or strong in comparison with the linear part of the NLSE, the power spectra are shaped as narrow (quasidiscrete), evenly spaced spikes, unlike generic truly continuous (ergodic) spectra. We develop an analytical explanation for the emergence of these spectral features in the case of weak nonlinearity. In the strongly nonlinear regime, the presence of such structures is tracked numerically by performing simulations with random initial conditions. Some potentials that prevent ergodicity in this manner are of direct relevance to Bose-Einstein condensates: they naturally appear in 1D, 2D, and 3D Gross-Pitaevskii equations (GPEs), the quintic version of these equations, and a two-component GPE system.

Nonlinear Perturbation of Random Matrix Theory

We consider a system of linear oscillators, or quantum states, described by random matrix theory and analyze how its time evolution is affected by a nonlinear perturbation. Our numerical results show that above a certain chaos border a weak or moderate nonlinearity leads to a dynamical thermalization of a finite number of degrees of freedom with energy equipartition over linear eigenmodes as expected from the laws of classical statistical mechanics. The system temperature is shown to change in a broad range from positive to negative values, and the dependence of system characteristics on the initial injected energy is determined. Below the chaos border the dynamics is described by the Kolmogorov-Arnold-Moser integrability. Owing to universal features of random matrix theory we argue that the obtained results describe the generic properties of its nonlinear perturbation.

Modal phase-locking in multimode nonlinear optical fibers

Spatial beam self-cleaning, a manifestation of the Kerr effect in graded-index multimode fibers, involves a nonlinear transfer of power among modes, which leads to robust bell-shaped output beams. The resulting mode power distribution can be described by statistical mechanics arguments. Although the spatial coherence of the output beam was experimentally demonstrated, there is no direct study of modal phase evolutions. Based on a holographic mode decomposition method, we reveal that nonlinear spatial phase-locking occurs between the fundamental and its neighboring low-order modes, in agreement with theoretical predictions. As such, our results dispel the current belief that the spatial beam self-cleaning effect is the mere result of a wave thermalization process.

Thermalization of the Ablowitz-Ladik lattice in the presence of non-integrable perturbations

We investigate the statistical mechanics of the photonic Ablowitz-Ladik lattice, the integrable version of the discrete nonlinear Schrödinger equation. In this regard, we demonstrate that in the presence of perturbations, the complex response of this system can be accurately captured within the framework of optical thermodynamics. Along these lines, we shed light on the true relevance of chaos in the thermalization of the Ablowitz-Ladik system. Our results indicate that when linear and nonlinear perturbations are incorporated, this weakly nonlinear lattice will thermalize into a proper Rayleigh-Jeans distribution with a well-defined temperature and chemical potential, in spite of the fact that the underlying nonlinearity is non-local and hence does not have a multi-wave mixing representation. This result illustrates that in the supermode basis, a non-local and non-Hermitian nonlinearity can in fact properly thermalize this periodic array in the presence of two quasi-conserved quantities.

Analysis of an image noise sensitivity mechanism for matrix-operation-mode-decomposition and a strong anti-noise method

Mode decomposition (MD) based on the matrix operation (MDMO) is one of the fastest mode decomposition methods in fiber laser which has great potential for optical communications, nonlinear optics and spatial characterization applications. However, we found that the image noise sensitivity is the main limit to the accuracy of the original MDMO method, but improving the decomposition accuracy by using conventional image filtering methods is almost ineffective. By using the norm theory of matrices, the analysis result shows that both the image noise and the coefficient matrix condition number determine the total upper-bound error of the original MDMO method. Besides, the greater the condition number, the more sensitive of MDMO method is to noise. In addition, it is found that the local error of each mode information solution in the original MDMO method is different, which depends on the L2-norm of each row vector of the inverse coefficient matrix. Moreover, a more noise-insensitive MD method is achieved by screening out the information corresponding to large L2-norm. In particular, selecting the higher accuracy among the original MDMO method and such noise-insensitive method as the result in a single MD process, a strong anti-noise MD method was proposed in this paper, which displays high MD accuracy in strong noise for both near-filed and far-filed MD cases.

Coherence properties of light in highly multimoded nonlinear parabolic fibers under optical equilibrium conditions

We study the coherence characteristics of light propagating in nonlinear graded-index (GRIN) multimode fibers after attaining optical thermal equilibrium conditions. The role of optical temperature on the spatial mutual coherence function and the associated correlation area is systematically investigated. In this respect, we show that the coherence properties of the field at the output of a multimode nonlinear fiber can be controlled through its optical thermodynamic properties.

Observation of Light Thermalization to Negative-Temperature Rayleigh-Jeans Equilibrium States in Multimode Optical Fibers

Although the temperature of a thermodynamic system is usually believed to be a positive quantity, under particular conditions, negative-temperature equilibrium states are also possible. Negative-temperature equilibriums have been observed with spin systems, cold atoms in optical lattices, and two-dimensional quantum superfluids. Here we report the observation of Rayleigh-Jeans thermalization of light waves to negative-temperature equilibrium states. The optical wave relaxes to the equilibrium state through its propagation in a multimode optical fiber-i.e., in a conservative Hamiltonian system. The bounded energy spectrum of the optical fiber enables negative-temperature equilibriums with high energy levels (high-order fiber modes) more populated than low energy levels (low-order modes). Our experiments show that negative-temperature speckle beams are featured, in average, by a nonmonotonic radial intensity profile. The experimental results are in quantitative agreement with the Rayleigh-Jeans theory without free parameters. Bringing negative temperatures to the field of optics opens the door to the investigation of fundamental issues of negative-temperature states in a flexible experimental environment.

Universality of light thermalization in multimoded nonlinear optical systems

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.

Nonlinear dynamics of beam self-cleaning on LP11 mode in multimode fibers

We investigate the modal energy flow of the femtosecond-pulsed beam self-cleaning on LP₁₁ mode with the influence of different factors such as the initial fraction of LP₁₁ mode, initial peak power, distribution of high-order modes and the numerical aperture of the fiber. It is interesting that there is a critical value of the initial peak power, Pcr, which is the watershed, not only in the quantitatively dominant transverse mode converting from LP₁₁ mode to LP₀₁ mode, but also in the behavior of HOMs of the transition from Attractor to chaos. Our simulation results may provide a novel perspective to understanding the beam self-cleaning on LP₁₁ mode.

Interplay of Thermalization and Strong Disorder: Wave Turbulence Theory, Numerical Simulations, and Experiments in Multimode Optical Fibers

We address the problem of thermalization in the presence of a time-dependent disorder in the framework of the nonlinear Schrödinger (or Gross-Pitaevskii) equation with a random potential. The thermalization to the Rayleigh-Jeans distribution is driven by the nonlinearity. On the other hand, the structural disorder is responsible for a relaxation toward the homogeneous equilibrium distribution (particle equipartition), which thus inhibits thermalization (energy equipartition). On the basis of the wave turbulence theory, we derive a kinetic equation that accounts for the presence of strong disorder. The theory unveils the interplay of disorder and nonlinearity. It unexpectedly reveals that a nonequilibrium process of condensation and thermalization can take place in the regime where disorder effects dominate over nonlinear effects. We validate the theory by numerical simulations of the nonlinear Schrödinger equation and the derived kinetic equation, which are found in quantitative agreement without using any adjustable parameter. Experiments realized in multimode optical fibers with an applied external stress evidence the process of thermalization in the presence of strong disorder.

Thermal control of the topological edge flow in nonlinear photonic lattices

The chaotic evolution resulting from the interplay between topology and nonlinearity in photonic systems generally forbids the sustainability of optical currents. Here, we systematically explore the nonlinear evolution dynamics in topological photonic lattices within the framework of optical thermodynamics. By considering an archetypical two-dimensional Haldane photonic lattice, we discover several prethermal states beyond the topological phase transition point and a stable global equilibrium response, associated with a specific optical temperature and chemical potential. Along these lines, we provide a consistent thermodynamic methodology for both controlling and maximizing the unidirectional power flow in the topological edge states. This can be achieved by either employing cross-phase interactions between two subsystems or by exploiting self-heating effects in disordered or Floquet topological lattices. Our results indicate that photonic topological systems can in fact support robust photon transport processes even under the extreme complexity introduced by nonlinearity, an important feature for contemporary topological applications in photonics.

The nonlinear evolution dynamics in topological photonic lattices is systematically investigated within the framework of optical thermodynamics. This approach allows for the precise prediction of topological currents even under the extreme complexity introduced by nonlinearity.

Thermalization of Orbital Angular Momentum Beams in Multimode Optical Fibers

We report on the thermalization of light carrying orbital angular momentum in multimode optical fibers, induced by nonlinear intermodal interactions. A generalized Rayleigh-Jeans distribution of asymptotic mode composition is obtained, based on the conservation of the angular momentum. We confirm our predictions by numerical simulations and experiments based on holographic mode decomposition of multimode beams. Our work establishes new constraints for the achievement of spatial beam self-cleaning, giving previously unforeseen insights into the underlying physical mechanisms.