Symmetric and asymmetric solitons in twin-core nonlinear optical fibers

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.

A solvable model for symmetry-breaking phase transitions

Analytically solvable models are benchmarks in studies of phase transitions and pattern-forming bifurcations. Such models are known for phase transitions of the second kind in uniform media, but not for localized states (solitons), as integrable equations which produce solitons do not admit intrinsic transitions in them. We introduce a solvable model for symmetry-breaking phase transitions of both the first and second kinds (alias sub- and supercritical bifurcations) for solitons pinned to a combined linear-nonlinear double-well potential, represented by a symmetric pair of delta-functions. Both self-focusing and defocusing signs of the nonlinearity are considered. In the former case, exact solutions are produced for symmetric and asymmetric solitons. The solutions explicitly demonstrate a switch between the symmetry-breaking transitions of the first and second kinds (i.e., sub- and supercritical bifurcations, respectively). In the self-defocusing model, the solution demonstrates the transition of the second kind which breaks antisymmetry of the first excited state.

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.

Zigzag Solitons and Spontaneous Symmetry Breaking in Discrete Rabi Lattices with Long-Range Hopping

A new type of discrete soliton, which we call zigzag solitons, is founded in two-component discrete Rabi lattices with long-range hopping. The spontaneous symmetry breaking (SSB) of zigzag solitons is also studied. Through numerical simulation, we found that by enhancing the intensity of the long-range linearly-coupled effect or increasing the total input power, the SSB process from the symmetric soliton to the asymmetric soliton will switch from the supercritical to subcritical type. These results can help us better understand both the discrete solitons and the Rabi coupled effect.

Moving Bragg grating solitons in a semilinear dual-core system with dispersive reflectivity

The existence, stability and collision dynamics of moving Bragg grating solitons in a semilinear dual-core system where one core has the Kerr nonlinearity and is equipped with a Bragg grating with dispersive reflectivity, and the other core is linear are investigated. It is found that moving soliton solutions exist as a continuous family of solutions in the upper and lower gaps of the system's linear spectrum. The stability of the moving solitons are investigated by means of systematic numerical stability analysis, and the effect and interplay of various parameters on soliton stability are analyzed. We have also systematically investigated the characteristics of collisions of counter-propagating solitons. In-phase collisions can lead to a variety of outcomes such as passage of solitons through each other with increased, reduced or unchanged velocities, asymmetric separation of solitons, merger of solitons into a quiescent one, formation of three solitons (one quiescent and two moving ones) and destruction of both solitons. The outcome regions of in-phase collisions are identified in the plane of dispersive reflectivity versus frequency. The effects of coupling coefficient, relative group velocity in the linear core, soliton velocity and dispersive reflectivity and the initial phase difference on the outcomes of collisions are studied.

PT-symmetric couplers with competing cubic-quintic nonlinearities

We introduce a one-dimensional model of the parity-time ( PT)-symmetric coupler, with mutually balanced linear gain and loss acting in the two cores, and nonlinearity represented by the combination of self-focusing cubic and defocusing quintic terms in each core. The system may be realized in optical waveguides, in the spatial and temporal domains alike. Stationary solutions for PT-symmetric solitons in the systems are tantamount to their counterparts in the ordinary coupler with the cubic-quintic nonlinearity, where the spontaneous symmetry breaking of solitons is accounted for by bifurcation loops. A novel problem is stability of the PT-symmetric solitons, which is affected by the competition of the PT symmetry, linear coupling, cubic self-focusing, and quintic defocusing. As a result, the solitons become unstable against symmetry breaking with the increase of the energy (alias integral power, in terms of the spatial-domain realization), and they retrieve the stability at still larger energies. Above a certain value of the strength of the quintic self-defocusing, the PT symmetry of the solitons becomes unbreakable. In the same system, PT-antisymmetric solitons are entirely unstable. We identify basic scenarios of the evolution of unstable solitons, which may lead to generation of additional ones, while stronger instability creates expanding quasi-turbulent patterns with limited amplitudes. Collisions between stable solitons are demonstrated to be quasi-elastic.

Two-dimensional solitons in conservative and parity-time-symmetric triple-core waveguides with cubic-quintic nonlinearity

We analyze a system of three two-dimensional nonlinear Schrödinger equations coupled by linear terms and with the cubic-quintic (focusing-defocusing) nonlinearity. We consider two versions of the model: conservative and parity-time (PT) symmetric. These models describe triple-core nonlinear optical waveguides, with balanced gain and losses in the PT-symmetric case. We obtain families of soliton solutions and discuss their stability. The latter study is performed using a linear stability analysis and checked with direct numerical simulations of the evolutional system of equations. Stable solitons are found in the conservative and PT-symmetric cases. Interactions and collisions between the conservative and PT-symmetric solitons are briefly investigated, as well.

Coexisting rogue waves within the (2+1)-component long-wave-short-wave resonance

The coexistence of two different types of fundamental rogue waves is unveiled, based on the coupled equations describing the (2+1)-component long-wave-short-wave resonance. For a wide range of asymptotic background fields, each family of three rogue wave components can be triggered by using a slight deterministic alteration to the otherwise identical background field. The ability to trigger markedly different rogue wave profiles from similar initial conditions is confirmed by numerical simulations. This remarkable feature, which is absent in the scalar nonlinear Schrödinger equation, is attributed to the specific three-wave interaction process and may be universal for a variety of multicomponent wave dynamics spanning from oceanography to nonlinear optics.

Multisoliton Newton's cradles and supersolitons in regular and parity-time-symmetric nonlinear couplers

We demonstrate the existence of stable collective excitation in the form of "supersolitons" propagating through chains of solitons with alternating signs (i.e., Newton's cradles built of solitons) in nonlinear optical couplers, including the parity-time-symmetric (PT-symmetric) version thereof. In the regular coupler, stable supersolitons are created in the cradles composed of both symmetric solitons and asymmetric ones with alternating polarities. Collisions between moving supersolitons are investigated too, by means of direct simulations in both the regular and PT-symmetric couplers.

Stability boundary and collisions of two-dimensional solitons in PT-symmetric couplers with the cubic-quintic nonlinearity

We introduce one- and two-dimensional (1D and 2D) models of parity-time (PT)-symmetric couplers with the mutually balanced linear gain and loss applied to the two cores and cubic-quintic (CQ) nonlinearity acting in each one. The 2D and 1D models may be realized in dual-core optical wave guides in the spatiotemporal and spatial domains, respectively. Stationary solutions for PT-symmetric solitons in these systems reduce to their counterparts in the usual coupler. The most essential problem is the stability of the solitons, which become unstable against symmetry breaking with the increase of the energy (norm) and retrieve the stability at still larger energies. The boundary value of the intercore-coupling constant, above which the solitons are completely stable, is found by means of an analytical approximation, based on the cw (zero-dimensional) counterpart of the system. The approximation demonstrates good agreement with numerical findings for the 1D and 2D solitons. Numerical results for the stability limits of the 2D solitons are obtained by means of the computation of eigenvalues for small perturbations, and verified in direct simulations. Although large parts of the soliton families are unstable, the instability is quite weak. Collisions between 2D solitons in the PT-symmetric coupler are studied by means of simulations. Outcomes of the collisions are inelastic but not destructive, as they do not break the PT symmetry.

Interface solitons in locally linked two-dimensional lattices

Existence, stability, and dynamics of soliton complexes, centered at the site of a single transverse link connecting two parallel two-dimensional (2D) lattices, are investigated. The system with the onsite cubic self-focusing nonlinearity is modeled by the pair of discrete nonlinear Schrödinger equations linearly coupled at the single site. Symmetric, antisymmetric, and asymmetric complexes are constructed by means of the variational approximation (VA) and numerical methods. The VA demonstrates that the antisymmetric soliton complexes exist in the entire parameter space, while the symmetric and asymmetric modes can be found below a critical value of the coupling parameter. Numerical results confirm these predictions. The symmetric complexes are destabilized via a supercritical symmetry-breaking pitchfork bifurcation, which gives rise to stable asymmetric modes. The antisymmetric complexes are subject to oscillatory and exponentially instabilities in narrow parametric regions. In bistability areas, stable antisymmetric solitons coexist with either symmetric or asymmetric ones.

Symmetry breaking of solitons in two-component Gross-Pitaevskii equations

We revisit the problem of the spontaneous symmetry breaking (SSB) of solitons in two-component linearly coupled nonlinear systems, adding the nonlinear interaction between the components. With this feature, the system may be realized in new physical settings, in terms of optics and the Bose-Einstein condensate (BEC). SSB bifurcation points are found analytically, for both symmetric and antisymmetric solitons (the symmetry between the two components is meant here). Asymmetric solitons, generated by the bifurcations, are described by means of the variational approximation (VA) and numerical methods, demonstrating good accuracy of the variational results. In the space of the self-phase-modulation (SPM) parameter and soliton's norm, a border separating stable symmetric and asymmetric solitons is identified. The nonlinear coupling may change the character of the SSB bifurcation, from subcritical to supercritical. Collisions between moving asymmetric and symmetric solitons are investigated too. Antisymmetric solitons are destabilized by a supercritical bifurcation, which gives rise to self-confined modes featuring Josephson oscillations, instead of stationary states with broken antisymmetry. An additional instability against delocalized perturbations is also found for the antisymmetric solitons.

Spontaneous symmetry breaking in coupled parametrically driven waveguides

We introduce a system of linearly coupled parametrically driven damped nonlinear Schrödinger equations, which models a laser based on a nonlinear dual-core waveguide with parametric amplification symmetrically applied to both cores. The model may also be realized in terms of parallel ferromagnetic films, in which the parametric gain is provided by an external field. We analyze spontaneous symmetry breaking (SSB) of fundamental and multiple solitons in this system, which was not studied systematically before in linearly coupled dissipative systems with intrinsic nonlinearity. For fundamental solitons, the analysis reveals three distinct SSB scenarios. Unlike the standard dual-core-fiber model, the present system gives rise to a vast bistability region, which may be relevant to applications. Other noteworthy findings are restabilization of the symmetric soliton after it was destabilized by the SSB bifurcation, and the existence of a generic situation with all solitons unstable in the single-component (decoupled) model, while both symmetric and asymmetric solitons may be stable in the coupled system. The stability of the asymmetric solitons is identified via direct simulations, while for symmetric and antisymmetric ones the stability is verified too through the computation of stability eigenvalues, families of antisymmetric solitons being entirely unstable. In this way, full stability maps for the symmetric solitons are produced. We also investigate the SSB bifurcation of two-soliton bound states (it breaks the symmetry between the two components, while the two peaks in the shape of the soliton remain mutually symmetric). The family of the asymmetric double-peak states may decouple from its symmetric counterpart, being no longer connected to it by the bifurcation, with a large portion of the asymmetric family remaining stable.

Symmetry breaking in linearly coupled dynamical lattices

We examine one- and two-dimensional models of linearly coupled lattices of the discrete-nonlinear-Schrödinger type. Analyzing ground states of the system with equal powers (norms) in the two components, we find a symmetry-breaking phenomenon beyond a critical value of the total power. Asymmetric states, with unequal powers in their components, emerge through a subcritical pitchfork bifurcation, which, for very weakly coupled lattices, changes into a supercritical one. We identify the stability of various solution branches. Dynamical manifestations of the symmetry breaking are studied by simulating the evolution of the unstable branches. The results present the first example of spontaneous symmetry breaking in two-dimensional lattice solitons. This feature has no counterpart in the continuum limit because of the collapse instability in the latter case.

Solitary pulses in linearly coupled Ginzburg-Landau equations

This article presents a brief review of dynamical models based on systems of linearly coupled complex Ginzburg-Landau (CGL) equations. In the simplest case, the system features linear gain, cubic nonlinearity (possibly combined with cubic loss), and group-velocity dispersion (GVD) in one equation, while the other equation is linear, featuring only intrinsic linear loss. The system models a dual-core fiber laser, with a parallel-coupled active core and an additional stabilizing passive (lossy) one. The model gives rise to exact analytical solutions for stationary solitary pulses (SPs). The article presents basic results concerning stability of the SPs; interactions between pulses are also considered, as are dark solitons (holes). In the case of the anomalous GVD, an unstable stationary SP may transform itself, via the Hopf bifurcation, into a stable localized breather. Various generalizations of the basic system are briefly reviewed too, including a model with quadratic (second-harmonic-generating) nonlinearity and a recently introduced model of a different but related type, based on linearly coupled CGL equations with cubic-quintic nonlinearity. The latter system features spontaneous symmetry breaking of stationary SPs, and also the formation of stable breathers.

Optical solitons as quantum objects

The intensity of classical bright solitons propagating in linearly coupled identical fibers can be distributed either in a stable symmetric state at strong coupling or in a stable asymmetric state if the coupling is small enough. In the first case, if the initial state is not the equilibrium state, the intensity may switch periodically from fiber to fiber, while in the second case the asymmetrical state remains forever, with most of its energy in either fiber. The latter situation makes a state of propagation with two exactly reciprocal realizations. In the quantum case, such a situation does not exist as an eigenstate because of the quantum tunneling between the two fibers. Such a tunneling is a purely quantum phenomenon without counterpart in the classical theory. We estimate the rate of tunneling by quantizing a simplified dynamics derived from the original Lagrangian equations with test functions. This tunneling could be within reach of the experiments, particularly if the quantum coherence of the soliton can be maintained over a sufficient amount of time.

Quasisymmetric and asymmetric gap solitons in linearly coupled Bragg gratings with a phase shift

We introduce a model including two linearly coupled Bragg gratings, with a mismatch (phase shift theta) between them. The model may be realized as parallel-coupled fiber Bragg gratings (FBGs), or, in the spatial domain, as two parallel planar waveguides carrying diffraction gratings. The phase shift induced by a shear stress may be used to design a different type of FBG sensor. In the absence of the mismatch, the symmetry-breaking bifurcation of gap solitons (GSs) in this model was investigated before. Our objective is to study how mismatch theta affects families of symmetric and asymmetric GSs, and the bifurcation between them. We find that the system's band gap is always filled with solitons (for theta not equal to 0, the gap's width does not depend on coupling constant lambda if it exceeds some minimum value). The largest velocity of the moving soliton, cmax, is found as a function of theta and lambda (cmax grows with theta). The mismatch transforms symmetric GSs into quasisymmetric (QS) ones, in which the two components are not identical, but their peak powers and energies are equal. The mismatch also breaks the spatial symmetry of the GSs. The QS solitons are stable against symmetry-breaking perturbations as long as asymmetric (AS) solutions do not exist. If theta is small, AS solitons emerge from their QS counterparts through a supercritical bifurcation. However, the bifurcation may become subcritical at larger theta. The condition for the stability against oscillatory perturbations (unrelated to the symmetry breaking) is essentially the same as in the ordinary FBG model: both QS and AS solitons are stable if their intrinsic frequency is positive (i.e., in a half of the band gap).

Vector solitons with an embedded domain wall

We present a class of soliton solutions to a system of two coupled nonlinear Schrödinger equations, with an intrinsic domain wall (DW) which separates regions occupied by two different fields. The model describes a binary mixture of two Bose-Einstein condensates (BECs) with interspecies repulsion. For the attractive or repulsive interactions inside each species, we find solutions which are bright or dark solitons in each component, while for the opposite signs of the intraspecies interaction, a bright-dark soliton pair is found (each time, with the intrinsic DW). These solutions can arise in the context of discrete lattices, and most of them can be supported in continuum settings by an external parabolic trap. The stability of the solitons with intrinsic DWs is examined, and the evolution of unstable ones is analyzed. We also briefly discuss the possibility of generating such families of solutions in the presence of linear coupling between the components, and an application of the model to bimodal light propagation in nonlinear optics.

Self-similar propagation and compression of chirped self-similar waves in asymmetric twin-core fibers with nonlinear gain

Ultrashort-pulse propagation in asymmetric twin-core fiber amplifiers is studied with the aid of self-similarity analysis of the nonlinear Schrödinger-type equation interacting with a source, variable dispersion, variable Kerr nonlinearity, variable gain or loss, and nonlinear gain. Exact chirped pulses that can propagate self-similarly subject to simple scaling rules of this model have been found. It is reported that the pulse position of these chirped pulses can be precisely piloted by appropriately tailoring the dispersion profile. This fact is profitably exploited to achieve optimal pulse compression of these chirped self-similar solutions.

Nonlinear compression of solitary waves in asymmetric twin-core fibers

We demonstrate a different pulse compression technique based on exact solutions to the nonlinear Schrödinger-type equation interacting with a source, variable dispersion, variable Kerr nonlinearity, and variable gain or loss. We show that this model is appropriate for the pulse propagation in asymmetric twin-core fibers. The chirped pulses are compressed due to the nonlinearity as well as dispersion management as also due to the space dependence of the gain coefficient. We also obtain singular solitary wave solutions, pertaining to extreme increase of the amplitude due to self-focusing.

Symmetric and asymmetric solitons in linearly coupled Bragg gratings

We demonstrate that a symmetric system of two linearly coupled waveguides, with Kerr nonlinearity and resonant grating in both of them, gives rise to a family of symmetric and antisymmetric solitons in an exact analytical form, a part of which exists outside of the bandgap in the system's spectrum, i.e., they may be regarded as embedded solitons (ES's, i.e., the ones partly overlapping with the continuous spectrum). Parameters of the family are the soliton's amplitude and velocity. Asymmetric ES's, unlike the regular (nonembedded) gap solitons (GS's), do not exist in the system. Moreover, ES's exist even in the case when the system's spectrum contains no bandgap. The main issue is the stability of the solitons. We demonstrate that some symmetric ES's are stable, while all the antisymmetric solitons are unstable; an explanation is given to the latter property, based on the consideration of the system's Hamiltonian. We produce a full stability diagram, which comprises both embedded and regular solitons, quiescent and moving. A stability region for ES's is found around the point where the constant of the linear coupling between the two cores is equal to the Bragg-reflectivity coefficient accounting for the linear conversion between the right- and left-traveling waves in each core, i.e., the ES's are the "most endemic" solitary solitons in this system. The stability region quickly shrinks with the increase of the soliton's velocity c, and completely disappears when c exceeds half the maximum velocity. Collisions between stable moving solitons of various types are also considered, with a conclusion that the collisions are always quasielastic.

Localization in physical systems described by discrete nonlinear Schrodinger-type equations

Following a short introduction on localized modes in a model system, namely the discrete nonlinear Schrodinger equation, we present explicit results pertaining to three different physical systems described by similar equations. The applications range from the Raman scattering spectra of a complex electronic material through intrinsic localized vibrational modes, to the manifestation of an abrupt and irreversible delocalizing transition of Bose-Einstein condensates trapped in two-dimensional optical lattices, and to the instabilities of localized modes in coupled arrays of optical waveguides.

Discrete vector solitons in two-dimensional nonlinear waveguide arrays: solutions, stability, and dynamics

We identify and investigate bimodal (vector) solitons in models of square-lattice arrays of nonlinear optical waveguides. These vector self-localized states are, in fact, self-induced channels in a nonlinear photonic-crystal matrix. Such two-dimensional discrete vector solitons are possible in waveguide arrays in which each element carries two light beams that are either orthogonally polarized or have different carrier wavelengths. Estimates of the physical parameters necessary to support such soliton solutions in waveguide arrays are given. Using Newton relaxation methods, we obtain stationary vector-soliton solutions, and examine their stability through the computation of linearized eigenvalues for small perturbations. Our results may also be applicable to other systems such as two-component Bose-Einstein condensates trapped in a two-dimensional optical lattice.

Bragg-grating solitons in a semilinear dual-core system

We investigate the existence and stability of gap solitons in a double-core optical fiber, where one core has the Kerr nonlinearity and the other one is linear, with the Bragg grating (BG) written on the nonlinear core, while the linear one may or may not have a BG. The model considerably extends the previously studied families of BG solitons. For zero-velocity solitons, we find exact solutions in a limiting case when the group-velocity terms are absent in the equation for the linear core. In the general case, solitons are found numerically. Stability borders for the solitons are found in terms of an internal parameter of the soliton family. Depending on the frequency omega, the solitons may remain stable for large values of the group velocity in the linear core. Stable moving solitons are also found. They are produced by interaction of initially separated solitons, which shows a considerable spontaneous symmetry breaking in the case when the solitons attract each other.

Soliton interaction with an external traveling wave

The dynamics of soliton pulses in the nonlinear Schrodinger equation (NLSE) driven by an external traveling wave is studied analytically and numerically. The Hamiltonian structure of the system is used to show that, in the adiabatic approximation for a single soliton, the problem is integrable despite the large number of degrees of freedom. Fixed points of the system are found, and their linear stability is investigated. The fixed points correspond to a Doppler shifted resonance between the external wave and the soliton. The structure and topological changes of the phase space of the soliton parameters as functions of the strength of coupling are investigated. A physical derivation of the driven NLSE is given in the context of optical pulse propagation in asymmetric, twin-core optical fibers. The results can be applied to soliton stabilization and amplification.