Discrete and Semi-Discrete Multidimensional Solitons and Vortices: Established Results and Novel Findings

This article presents a concise survey of basic discrete and semi-discrete nonlinear models, which produce two- and three-dimensional (2D and 3D) solitons, and a summary of the main theoretical and experimental results obtained for such solitons. The models are based on the discrete nonlinear Schrödinger (DNLS) equations and their generalizations, such as a system of discrete Gross-Pitaevskii (GP) equations with the Lee-Huang-Yang corrections, the 2D Salerno model (SM), DNLS equations with long-range dipole-dipole and quadrupole-quadrupole interactions, a system of coupled discrete equations for the second-harmonic generation with the quadratic (χ(2)) nonlinearity, a 2D DNLS equation with a superlattice modulation opening mini-gaps, a discretized NLS equation with rotation, a DNLS coupler and its PT-symmetric version, a system of DNLS equations for the spin-orbit-coupled (SOC) binary Bose-Einstein condensate, and others. The article presents a review of the basic species of multidimensional discrete modes, including fundamental (zero-vorticity) and vortex solitons, their bound states, gap solitons populating mini-gaps, symmetric and asymmetric solitons in the conservative and PT-symmetric couplers, cuspons in the 2D SM, discrete SOC solitons of the semi-vortex and mixed-mode types, 3D discrete skyrmions, and some others.

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.

Rogue waves and instability arising from long-wave-short-wave resonance beyond the integrable regime

We consider instability and localized patterns arising from the long-wave-short-wave resonance in the nonintegrable regime numerically. We study the stability and instability of elliptic-function periodic waves with respect to subharmonic perturbations, whose period is a multiple of the period of the elliptic waves. We thus find the modulational instability (MI) of the corresponding dnoidal waves. Upon varying parameters of dnoidal waves, spectrally unstable ones can be transformed into stable states via the Hamiltonian Hopf bifurcation. For snoidal waves, we find a transition of the dominant instability scenario between the MI and the instability with a bubblelike spectrum. For cnoidal waves, we produce three variants of the MI. Evolution of the unstable states is also considered, leading to formation of rogue waves on top of the elliptic-wave and continuous-wave backgrounds.

Bimetal Thin Film, Semiconductors, and 2D Nanomaterials in SPR Biosensors: An Approach to Enhanced Urine Glucose Sensing

This work introduces a systematic approach for the development of Kretschmann configuration-based biosensors designed for non-invasive urine glucose detection. The methodology encompasses the utilization of various semiconductors, including Silicon (Si), Germanium (Ge), Gallium Nitride (GaN), Aluminum Nitride (AlN), and Indium Nitride (InN), in combination with a bimetallic layer (comprising Au and Ag films of equal thickness) to enhance the biosensor sensitivity. Additionally, 2D nanomaterials, such as Black Phosphorus and Graphene, are integrated into the semiconductor layers to enhance performance further. These configurations are meticulously optimized through the application of the transfer matrix method (TMM), and the sensing parameters are assessed using the angular modulation method. Among the semiconductors, AlN and GaN exhibit superior results. On these substrates, Graphene and Black phosphorous (BP) layers are applied, resulting in four final structures (thicknesses in nm): BK7/Au(26)/Ag(26)/Si(6)/BP(0.53)/Biosample, BK7/Au(26)/Ag(26)/AlN(14)/BP(0.53)/Biosample, BK7/Au(26)/Ag(26)/GaN(12)/BP(0.53)/Biosample, and BK7/Au(26)/Ag(26)/GaN(12)/Graphene(0.34)/Biosample. These biosensors achieve Sensitivity(° /RIU) and Figure of Merit (FoM) (1/RIU) of 380, 360, 440, 400, and 58.5, 90, 90.65, and 82.4, respectively. Subsequently, these high-performing sensors undergo testing with actual urine glucose samples. Among them, two biosensors, BK7/Au(26)/Ag(26)/AlN(14)/BP (0.53)/Biosample and BK7/Au(26)/Ag(26)/GaN(14)/Graphene(0.34)/Biosample, exhibit outstanding performance, with sensitivities (° /RIU) and FoM (1/RIU) of 394.44 & 294.44, and 112.6 & 92.01 respectively. A comparison is also made with relevant previously published work, revealing improved performance in glucose detection.

Analysis of high-contrast all-optical dual-wavelength switching in asymmetric dual-core fibers

We systematically present experimental and theoretical results for the dual-wavelength switching of 1560 nm, 75 fs signal pulses (SPs) driven by 1030 nm, and 270 fs control pulses (CPs) in a dual-core fiber (DCF). We demonstrate a switching contrast of 31.9 dB, corresponding to a propagation distance of 14 mm, achieved by launching temporally synchronized SP-CP pairs into the fast core of the DCF with moderate inter-core asymmetry. Our analysis employs a system of three coupled propagation equations to identify the compensation of the asymmetry by nonlinearity as the physical mechanism behind the efficient switching performance.

Semidiscrete optical vortex droplets in quasi-phase-matched photonic crystals

What we believe is a new scheme for producing semidiscrete self-trapped vortices ("swirling photon droplets") in photonic crystals with competing quadratic (χ(2)) and self-defocusing cubic (χ(3)) nonlinearities is proposed. The photonic crystal is designed with a striped structure, in the form of spatially periodic modulation of the χ(2) susceptibility, which is imposed by the quasi-phase-matching technique. Unlike previous realizations of semidiscrete optical modes in composite media, built as combinations of continuous and arrayed discrete waveguides, the semidiscrete vortex "droplets" are produced here in the fully continuous medium. This work reveals that the system supports two types of semidiscrete vortex droplets, viz., onsite- and intersite-centered ones, which feature, respectively, odd and even numbers of stripes, N. Stability areas for the states with different values of N are identified in the system's parameter space. Some stability areas overlap with each other, giving rise to the multistability of states with different N. The coexisting states are mutually degenerate, featuring equal values of the Hamiltonian and propagation constant. An experimental scheme to realize the droplets is outlined, suggesting new possibilities for the long-distance transmission of nontrivial vortex beams in nonlinear media.

Three-dimensional solitons in Rydberg-dressed cold atomic gases with spin-orbit coupling

We present numerical results for three-dimensional (3D) solitons with symmetries of the semi-vortex (SV) and mixed-mode (MM) types, which can be created in spinor Bose-Einstein condensates of Rydberg atoms under the action of the spin-orbit coupling (SOC). By means of systematic numerical computations, we demonstrate that the interplay of SOC and long-range spherically symmetric Rydberg interactions stabilize the 3D solitons, improving their resistance to collapse. We find how the stability range depends on the strengths of the SOC and Rydberg interactions and the soft-core atomic radius.

Ring-shaped quantum droplets with hidden vorticity in a radially periodic potential

We study the stability and characteristics of two-dimensional circular quantum droplets (QDs) with embedded hidden vorticity (HV), i.e., opposite angular momenta in two components, formed by binary Bose-Einstein condensates (BECs) trapped in a radially periodic potential. The system is modeled by the Gross-Pitaevskii equations with the Lee-Huang-Yang terms, which represent the higher-order self-repulsion induced by quantum fluctuations around the mean-field state, and a potential which is a periodic function of the radial coordinate. Ring-shaped QDs with high winding numbers (WNs) of the HV type, which are trapped in particular circular troughs of the radial potential, are produced by means of the imaginary-time-integration method. Effects of the depth and period of the potential on these QD states are studied. The trapping capacity of individual circular troughs is identified. Stable compound states in the form of nested multiring patterns are constructed too, including ones with WNs of opposite signs. The stably coexisting ring-shaped QDs with different WNs can be used for the design of BEC-based data-storage schemes.

Stable higher-charge vortex solitons in the cubic-quintic medium with a ring potential

We put forward a model for trapping stable optical vortex solitons (VSs) with high topological charges m. The cubic-quintic nonlinear medium with an imprinted ring-shaped modulation of the refractive index is shown to support two branches of VSs, which are controlled by the radius, width, and depth of the modulation profile. While the lower-branch VSs are unstable in their nearly whole existence domain, the upper branch is completely stable. Vortex solitons with m ≤  12 obey the anti-Vakhitov-Kolokolov stability criterion. The results suggest possibilities for the creation of stable narrow optical VSs with a low power, carrying higher vorticities.

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.

Formation of Rogue Waves and Modulational Instability with Zero-Wavenumber Gain in Multicomponent Systems with Coherent Coupling

It is known that rogue waves (RWs) are generated by the modulational instability (MI) of the baseband type. Starting with the Bers-Kaup-Reiman system for three-wave resonant interactions, we identify a specific RW-building mechanism based on MI which includes zero wavenumber in the gain band. An essential finding is that this mechanism works solely under a linear relation between the MI gain and a vanishingly small wavenumber of the modulational perturbation. The same mechanism leads to the creation of RWs by MI in other multicomponent systems-in particular, in the massive Thirring model.

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.

One-dimensional Townes solitons in dual-core systems with localized coupling

The recent creation of Townes solitons (TSs) in binary Bose-Einstein condensates and experimental demonstration of spontaneous symmetry breaking (SSB) in solitons propagating in dual-core optical fibers has drawn renewed interest in the TS and SSB phenomenology in these and other settings. In particular, stabilization of TSs, which are always unstable in free space, is a relevant problem with various ramifications. We introduce a system which admits exact solutions for both TSs and SSB of solitons. It is based on a dual-core waveguide with quintic self-focusing and fused (localized) coupling between the cores. The respective system of coupled nonlinear Schrödinger equations gives rise to exact solutions for full families of symmetric and asymmetric solitons, which are produced by the supercritical SSB bifurcation (i.e., the symmetry-breaking phase transition of the second kind). Stability boundaries of asymmetric solitons are identified by dint of numerical methods. Unstable solitons spontaneously transform into robust moderately asymmetric breathers or strongly asymmetric states with small intrinsic oscillations. The setup can be used in the design of photonic devices operating in coupling and switching regimes.

Topological solitonic macromolecules

Being ubiquitous, solitons have particle-like properties, exhibiting behaviour often associated with atoms. Bound solitons emulate dynamics of molecules, though solitonic analogues of polymeric materials have not been considered yet. Here we experimentally create and model soliton polymers, which we call "polyskyrmionomers", built of atom-like individual solitons characterized by the topological invariant representing the skyrmion number. With the help of nonlinear optical imaging and numerical modelling based on minimizing the free energy, we reveal how topological point defects bind the solitonic quasi-atoms into polyskyrmionomers, featuring linear, branched, and other macromolecule-resembling architectures, as well as allowing for encoding data by spatial distributions of the skyrmion number. Application of oscillating electric fields activates diverse modes of locomotion and internal vibrations of these self-assembled soliton structures, which depend on symmetry of the solitonic macromolecules. Our findings suggest new designs of soliton meta matter, with a potential for the use in fundamental research and technology.

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.

Dichromatic "Breather Molecules" in a Mode-Locked Fiber Laser

Bound states of solitons ("molecules") occur in various settings, playing an important role in the operation of fiber lasers, optical emulation, encoding, and communications. Soliton interactions are generally related to breathing dynamics in nonlinear dissipative systems, and maintain potential applications in spectroscopy. In the present work, dichromatic breather molecules (DBMs) are created in a synchronized mode-locked fiber laser. Real-time delay-shifting interference spectra are measured to display the temporal evolution of the DBMs, that cannot be observed by means of the usual real-time spectroscopy. As a result, robust out-of-phase vibrations are found as a typical intrinsic mode of DBMs. The same bound states are produced numerically in the framework of a model combining equations for the population inversion in the mode-locked laser and cross-phase-modulation-coupled complex Ginzburg-Landau equations for amplitudes of the optical fields in the fiber segments of the laser cavity. The results demonstrate that the Q-switching instability induces the onset of breathing oscillations. The findings offer new possibilities for the design of various regimes of the operation of ultrafast lasers.

Vortex Solitons in Quasi-Phase-Matched Photonic Crystals

We report solutions for stable compound solitons in a three-dimensional quasi-phase-matched photonic crystal with the quadratic (χ^{(2)}) nonlinearity. The photonic crystal is introduced with a checkerboard structure, which can be realized by means of the available technology. The solitons are built as four-peak vortex modes of two types, rhombuses and squares (intersite- and onsite-centered self-trapped states, respectively). Their stability areas are identified in the system's parametric space (rhombuses occupy an essentially broader stability domain), while all bright vortex solitons are subject to strong azimuthal instability in uniform χ^{(2)} media. Possibilities for experimental realization of the solitons are outlined.

Formations and dynamics of two-dimensional spinning asymmetric quantum droplets controlled by a PT-symmetric potential

In this paper, vortex solitons are produced for a variety of 2D spinning quantum droplets (QDs) in a PT-symmetric potential, modeled by the amended Gross-Pitaevskii equation with Lee-Huang-Yang corrections. In particular, exact QD states are obtained under certain parameter constraints, providing a guide to finding the respective generic family. In a parameter region of the unbroken PT symmetry, different families of QDs originating from the linear modes are obtained in the form of multipolar and vortex droplets at low and high values of the norm, respectively, and their stability is investigated. In the spinning regime, QDs become asymmetric above a critical rotation frequency, most of them being stable. The effect of the PT-symmetric potential on the spinning and nonspinning QDs is explored by varying the strength of the gain-loss distribution. Generally, spinning QDs trapped in the PT-symmetric potential exhibit asymmetry due to the energy flow affected by the interplay of the gain-loss distribution and rotation. Finally, interactions between spinning or nonspinning QDs are explored, exhibiting elastic collisions under certain conditions.

Experimental realisations of the fractional Schrödinger equation in the temporal domain

The fractional Schrödinger equation (FSE)-a natural extension of the standard Schrödinger equation-is the basis of fractional quantum mechanics. It can be obtained by replacing the kinetic-energy operator with a fractional derivative. Here, we report the experimental realisation of an optical FSE for femtosecond laser pulses in the temporal domain. Programmable holograms and the single-shot measurement technique are respectively used to emulate a Lévy waveguide and to reconstruct the amplitude and phase of the pulses. Varying the Lévy index of the FSE and the initial pulse, the temporal dynamics is observed in diverse forms, including solitary, splitting and merging pulses, double Airy modes, and "rain-like" multi-pulse patterns. Furthermore, the transmission of input pulses carrying a fractional phase exhibits a "fractional-phase protection" effect through a regular (non-fractional) material. The experimentally generated fractional time-domain pulses offer the potential for designing optical signal-processing schemes.

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.

Interactions of solitary waves in the Adlam-Allen model

We study the interactions of two or more solitary waves in the Adlam-Allen model describing the evolution of a (cold) plasma of positive and negative charges, in the presence of electric and transverse magnetic fields. In order to show that the interactions feature an exponentially repulsive nature, we elaborate two distinct approaches: (a) using energetic considerations and the Hamiltonian structure of the model, and (b) using the so-called Manton method. We compare these findings with results of direct simulations, and we identify adjustments necessary to achieve a quantitative match between them. Additional connections are made, such as with solitons of the Korteweg-de Vries equation. New challenges are identified in connection to this model and its solitary waves.

Stability limits for modes held in alternating trapping-expulsive potentials

We elaborate a scheme of trapping-expulsion management (TEM), in the form of the quadratic potential periodically switching between confinement and expulsion, as a means of stabilization of two-dimensional dynamical states against the backdrop of the critical collapse driven by the cubic self-attraction with strength g. The TEM scheme may be implemented, as spatially or temporally periodic modulations, in optics or BEC, respectively. The consideration is carried out by dint of numerical simulations and variational approximation (VA). In terms of the VA, the dynamics amounts to a nonlinear Ermakov equation, which, in turn, is tantamount to a linear Mathieu equation. Stability boundaries are found as functions of g and parameters of the periodic modulation of the trapping potential. Below the usual collapse threshold, which is known, in the numerical form, as g<g_{c}^{(num)}≈5.85 (in the standard notation), the stability is limited by the onset of the parametric resonance. This stability limit, including the setup with the self-repulsive sign of the cubic term (g<0), is accurately predicted by the VA. At g>g_{c}^{(num)}, the collapse threshold is found with the help of full numerical simulations. The relative increase of g_{c} above g_{c}^{(num)} is ≈1.5%. It is a meaningful result, even if its size is small, because the collapse threshold is a universal constant which is difficult to change.

Trapping wave fields in an expulsive potential by means of linear coupling

We demonstrate the existence of confined states in one- and two-dimensional (1D and 2D) systems of two linearly coupled components, with the confining harmonic-oscillator (HO) potential acting upon one component and an expulsive anti-HO potential acting upon the other. The systems can be implemented in optical and BEC dual-core waveguides. In the 1D linear system, codimension-one solutions are found in an exact form for the ground state (GS) and dipole mode (the first excited state). Generic solutions are produced by means of the variational approximation and are found in a numerical form. Exact codimension-one solutions and generic numerical ones are also obtained for the GS and vortex states in the 2D system (the exact solutions are found for all values of the vorticity). Both the trapped and antitrapped components of the bound states may be dominant ones, in terms of the norm. The localized modes may be categorized as bound states in continuum, as they coexist with delocalized ones. The 1D states, as well as the GS in 2D, are weakly affected and remain stable if the self-attractive or repulsive nonlinearity is added to the system. The self-attraction makes the vortex states unstable against splitting, while they remain stable under the action of the self-repulsion.

Fermi-Pasta-Ulam phenomena and persistent breathers in the harmonic trap

We consider the long-term weakly nonlinear evolution governed by the two-dimensional nonlinear Schrödinger (NLS) equation with an isotropic harmonic oscillator potential. The dynamics in this regime is dominated by resonant interactions between quartets of linear normal modes, accurately captured by the corresponding resonant approximation. Within this approximation, we identify Fermi-Pasta-Ulam-like recurrence phenomena, whereby the normal-mode spectrum passes in close proximity of the initial configuration, and two-mode states with time-independent mode amplitude spectra that translate into long-lived breathers of the original NLS equation. We comment on possible implications of these findings for nonlinear optics and matter-wave dynamics in Bose-Einstein condensates.