Metastable soliton necklaces supported by fractional diffraction and competing nonlinearities

Basic fractional nonlinear-wave models and solitons

This review article provides a concise summary of one- and two-dimensional models for the propagation of linear and nonlinear waves in fractional media. The basic models, which originate from Laskin's fractional quantum mechanics and more experimentally relevant setups emulating fractional diffraction in optics, are based on the Riesz definition of fractional derivatives, which are characterized by the respective Lévy indices. Basic species of one-dimensional solitons, produced by the fractional models which include cubic or quadratic nonlinear terms, are outlined too. In particular, it is demonstrated that the variational approximation is relevant in many cases. A summary of the recently demonstrated experimental realization of the fractional group-velocity dispersion in fiber lasers is also presented.

One-dimensional Lévy quasicrystal

Space-fractional quantum mechanics (SFQM) is a generalization of the standard quantum mechanics when the Brownian trajectories in Feynman path integrals are replaced by Lévy flights. We introduce Lévy quasicrystal by discretizing the space-fractional Schrödinger equation using the Grünwald-Letnikov derivatives and adding on-site quasiperiodic potential. The discretized version of the usual Schrödinger equation maps to the Aubry-André (AA) Hamiltonian, which supports localization-delocalization transition even in one dimension. We find the similarities between Lévy quasicrystal and the AA model with power-law hopping, and show that the Lévy quasicrystal supports a delocalization-localization transition as one tunes the quasiperiodic potential strength and shows the coexistence of localized and delocalized states separated by mobility edge. Hence, a possible realization of SFQM in optical experiments should be a new experimental platform to test the predictions of AA models in the presence of power-law hopping.

Symmetry-breaking transitions in quiescent and moving solitons in fractional couplers

We consider phase transitions, in the form of spontaneous symmetry breaking (SSB) bifurcations of solitons, in dual-core couplers with fractional diffraction and cubic self-focusing acting in each core, characterized by Lévy index α. The system represents linearly coupled optical waveguides with the fractional paraxial diffraction or group-velocity dispersion (the latter system was used in a recent experiment [Nat. Commun. 14, 222 (2023)10.1038/s41467-023-35892-8], which demonstrated the first observation of the wave propagation in an effectively fractional setup). By dint of numerical computations and variational approximation, we identify the SSB in the fractional coupler as the bifurcation of the subcritical type (i.e., the symmetry-breaking phase transition of the first kind), whose subcriticality becomes stronger with the increase of fractionality 2-α, in comparison with very weak subcriticality in the case of the nonfractional diffraction, α=2. In the Cauchy limit of α→1, it carries over into the extreme subcritical bifurcation, manifesting backward-going branches of asymmetric solitons which never turn forward. The analysis of the SSB bifurcation is extended for moving (tilted) solitons, which is a nontrivial problem because the fractional diffraction does not admit Galilean invariance. Collisions between moving solitons are studied too, featuring a two-soliton symmetry-breaking effect and merger of the solitons.

Domain walls in fractional media

Currently, much interest is drawn to the analysis of optical and matter-wave modes supported by the fractional diffraction in nonlinear media. We predict a new type of such states in the form of domain walls (DWs) in the two-component system of immiscible fields. Numerical study of the underlying system of fractional nonlinear Schrödinger equations demonstrates the existence and stability of DWs at all values of the respective Lévy index (α<2), which determines the fractional diffraction, and at all values of the XPM/SPM ratio β in the two-component system above the immiscibility threshold. The same conclusion is obtained for DWs in the system which includes the linear coupling, alongside the XPM interaction between the immiscible components. Analytical results are produced for the scaling of the DW's width. The DW solutions are essentially simplified in the special case of β=3, as well as close to the immiscibility threshold. In addition to symmetric DWs, asymmetric ones are constructed too, in the system with unequal diffraction coefficients and/or different Lévy indices of the two components.

Fractional Integrable Nonlinear Soliton Equations

Nonlinear integrable equations serve as a foundation for nonlinear dynamics, and fractional equations are well known in anomalous diffusion. We connect these two fields by presenting the discovery of a new class of integrable fractional nonlinear evolution equations describing dispersive transport in fractional media. These equations can be constructed from nonlinear integrable equations using a widely generalizable mathematical process utilizing completeness relations, dispersion relations, and inverse scattering transform techniques. As examples, this general method is used to characterize fractional extensions to two physically relevant, pervasive integrable nonlinear equations: the Korteweg-deVries and nonlinear Schrödinger equations. These equations are shown to predict superdispersive transport of nondissipative solitons in fractional media.

Generation of random soliton-like beams in a nonlinear fractional Schrödinger equation

We investigate the generation of random soliton-like beams based on the Kuznetsov-Ma solitons in a nonlinear fractional Schrödinger equation (NLFSE). For Lévy index α = 1, the Kuznetsov-Ma solitons split into two nondiffracting beams during propagation in linear regime. According to the different input positions of the Kuznetsov-Ma solitons, the diffraction-free beams can be divided into three different types: bright-dark, dark-bright and bright-bright beams. In the nonlinear regime, the Kuznetsov-Ma solitons can be evolved into random soliton-like beams due to the collapse. The number of soliton-like beams is related to the nonlinear coefficient and the Lévy index. The bigger the nonlinear coefficient, the more beams generated. Moreover, the peak intensity of soliton-like beams presents a Gaussian distribution under the large nonlinear effect. In practice, the evolution of KM soliton can be realized by a plane wave with a Gaussian perturbation, which can be confirmed that they have the similar dynamics of propagation. In two dimensions, the plane wave with a Gaussian perturbation can be evolved into a bright-dark axisymmetric ring beam in the linear regime. Under the nonlinear modulation, the energy accumulates to the center and finally breaks apart into random beam filaments.

Optical Solitons and Vortices in Fractional Media: A Mini-Review of Recent Results

The article produces a brief review of some recent results which predict stable propagation of solitons and solitary vortices in models based on the nonlinear Schrödinger equation (NLSE) including fractional one-dimensional or two-dimensional diffraction and cubic or cubic-quintic nonlinear terms, as well as linear potentials. The fractional diffraction is represented by fractional-order spatial derivatives of the Riesz type, defined in terms of the direct and inverse Fourier transform. In this form, it can be realized by spatial-domain light propagation in optical setups with a specially devised combination of mirrors, lenses, and phase masks. The results presented in the article were chiefly obtained in a numerical form. Some analytical findings are included too, in particular, for fast moving solitons and the results produced by the variational approximation. Moreover, dissipative solitons are briefly considered, which are governed by the fractional complex Ginzburg–Landau equation.

Existence, symmetry breaking bifurcation and stability of two-dimensional optical solitons supported by fractional diffraction

We study existence, bifurcation and stability of two-dimensional optical solitons in the framework of fractional nonlinear Schrödinger equation, characterized by its Lévy index, with self-focusing and self-defocusing saturable nonlinearities. We demonstrate that the fractional diffraction system with different Lévy indexes, combined with saturable nonlinearity, supports two-dimensional symmetric, antisymmetric and asymmetric solitons, where the asymmetric solitons emerge by way of symmetry breaking bifurcation. Different scenarios of bifurcations emerge with the change of stability: the branches of asymmetric solitons split off the branches of unstable symmetric solitons with the increase of soliton power and form a supercritical type bifurcation for self-focusing saturable nonlinearity; the branches of asymmetric solitons bifurcates from the branches of unstable antisymmetric solitons for self-defocusing saturable nonlinearity, featuring a convex shape of the bifurcation loops: an antisymmetric soliton loses its stability via a supercritical bifurcation, which is followed by a reverse bifurcation that restores the stability of the symmetric soliton. Furthermore, we found a scheme of restoration or destruction the symmetry of the antisymmetric solitons by controlling the fractional diffraction in the case of self-defocusing saturable nonlinearity.