A bound state attractor in optical turbulence

Abstract: We study numerically the nonintegrable dynamics of coherent, solitonic, nonlinear waves, in a spatially nonlocal nonlinear Schrodinger equation relevant to realistic modelling of optical systems: the Schrodinger-Helmholtz equation. We observe a single oscillating, coherent solitary wave emerging from a variety of initial conditions. Using the direct scattering transform of the (integrable) cubic nonlinear Schrodinger equation, we find that this structure is a bound state, comprising of a primary and secondary soliton whose amplitudes oscillate out of phase. We interpret this as the solitons periodically exchanging mass. We observe this bound state self-organising from a state of incoherent turbulence, and from solitonic structures launched into the system. When a single (primary) solitonic structure is launched, a resonance process between it and waves in the system generates the secondary soliton, resulting in the bound state. Further, when two solitons are initially launched, we show that they can merge if their phases are synchronised when they collide. When the system is launched from a turbulent state comprised of many initial solitons, we propose that the bound state formation is preceded a sequence of binary collisions, in which the mass is transferred on average from the weak soliton to the strong one, with occasional soliton mergers. Both processes lead to increasingly stronger and fewer dominant solitons. The final state -- a solitary double-solitonic bound state surrounded by weakly nonlinear waves -- is robust and ubiquitous. We propose that for nonlocal media, it is a more typical statistical attractor than a single-soliton attractor suggested in previous literature.