Search

We investigate the propagation characteristics of the chirped self-similar solitary waves in non-Kerr nonlinear media within the framework of the generalized nonlinear Schrödinger equation with distributed dispersion, two-power-law nonlinearities, and gain or loss. This model contains many special types of nonlinear equations that appear in various branches of contemporary physics. We extend the self-similar analysis presented for searching chirped self-similar structures of the cubic model to a more general problem involving two nonlinear terms of arbitrary power. A variety of exact linearly chirped localized solutions with interesting properties are derived in the presence of all physical effects. The solutions comprise bright, kink and antikink, and algebraic solitary wave solutions, illustrating the potentially rich set of self-similar pulses of the model. It is shown that these optical pulses possess a linear chirp that leads to efficient compression or amplification, and thus are particularly useful in the design of optical fiber amplifiers, optical pulse compressors, and solitary wave based communication links. Finally, the stability of the self-similar solutions is discussed numerically under finite initial perturbations.

We consider the reflectionless transport of solitons in networks. The system is modeled in terms of the nonlinear Schrödinger equation on metric graphs, for which transparent boundary conditions at the branching points are imposed. This approach allows to derive simple constraints, which link the equivalent usual Kirchhoff-type vertex conditions to the transparent ones. Our method is applied to a metric star graph. An extension to more complicated graph topologies is straightforward.

For waves described by the focusing nonlinear Schrödinger equation (FNLS), we present an effective dispersion relation (EDR) that arises dynamically from the interplay between the linear dispersion and the nonlinearity. The form of this EDR is parabolic for a robust family of "generic" FNLS waves and equals the linear dispersion relation less twice the total wave action of the wave in question multiplied by the square of the nonlinearity parameter. We derive an approximate form of this EDR explicitly in the limit of small nonlinearity and confirm it using the wave-number-frequency spectral (WFS) analysis, a Fourier-transform based method used for determining dispersion relations of observed waves. We also show that it extends to the FNLS the universal EDR formula for the defocusing Majda-McLaughlin-Tabak (MMT) model of weak turbulence. In addition, unexpectedly, even for some spatially periodic versions of multisolitonlike waves, the EDR is still a downward shifted linear-dispersion parabola, but the shift does not have a clear relation to the total wave action. Using WFS analysis and heuristic derivations, we present examples of parabolic and nonparabolic EDRs for FNLS waves and also waves for which no EDR exists.

We formulate, using the discrete nonlinear Schroedinger equation (DNLS), a general approach to encode and process information based on reservoir computing. Reservoir computing is a promising avenue for realizing neuromorphic computing devices. In such computing systems, training is performed only at the output level, by adjusting the output from the reservoir with respect to a target signal. In our formulation, the reservoir can be an arbitrary physical system, driven out of thermal equilibrium by an external driving. The DNLS is a general oscillator model with broad application in physics and we argue that our approach is completely general and does not depend on the physical realisation of the reservoir. The driving, which encodes the object to be recognised, acts as a thermodynamical force, one for each node in the reservoir. Currents associated to these thermodynamical forces in turn encode the output signal from the reservoir. As an example, we consider numerically the problem of supervised learning for pattern recognition, using as reservoir a network of nonlinear oscillators., Comment: 5 pages, 4 figures

An integrable nonlocal nonlinear Schr\"odinger (NLS) equation with clear physical motivations is proposed. This equation is obtained from a special reduction of the Manakov system, and it describes Manakov solutions whose two components are related by a parity symmetry. Since the Manakov system governs wave propagation in a wide variety of physical systems, our nonlocal equation has clear physical meanings. Solitons and multisolitons in this nonlocal equation are also investigated in the framework of Riemann-Hilbert formulations. Surprisingly, symmetry relations of discrete scattering data for this equation are found to be very complicated, where constraints between eigenvectors in the scattering data depend on the number and locations of the underlying discrete eigenvalues in a very complex manner. As a consequence, general $N$-solitons are difficult to obtain in the Riemann-Hilbert framework. However, one- and two-solitons are derived, and their dynamics investigated. It is found that two-solitons are generally not a nonlinear superposition of one-solitons, and they exhibit interesting dynamics such as meandering and sudden position shifts. As a generalization, other integrable and physically meaningful nonlocal equations are also proposed, which include NLS equations of reverse-time and reverse-space-time types as well as nonlocal Manakov equations of reverse-space, reverse-time, and reverse-space-time types.

Multiscale turbulence naturally develops and plays an important role in many fluid, gas, and plasma phenomena. Statistical models of multiscale turbulence usually employ Kolmogorov hypotheses of spectral locality of interactions (meaning that interactions primarily occur between pulsations of comparable scales) and scale-invariance of turbulent pulsations. However, optical turbulence described by the nonlinear Schrodinger equation exhibits breaking of both the Kolmogorov locality and scale-invariance. A weaker form of spectral locality that holds for multi-scale optical turbulence enables a derivation of simplified evolution equations that reduce the problem to a single scale modeling. We present the derivation of these equations for Kerr media with random inhomogeneities. Then, we find the analytical solution that exhibits a transition between inverse and direct energy cascades in optical turbulence.

We present one- and two-breather solutions of the fourth-order nonlinear Schr\"odinger equation. With several parameters to play with, the solution may take a variety of forms. We consider most of these cases including the general form and limiting cases when the modulation frequencies are 0 or coincide. The zero-frequency limit produces a combination of breather-soliton structures on a constant background. The case of equal modulation frequencies produces a degenerate solution that requires a special technique for deriving. A zero-frequency limit of this degenerate solution produces a rational second-order rogue wave solution with a stretching factor involved. Taking, in addition, the zero limit of the stretching factor transforms the second-order rogue waves into a soliton. Adding a differential shift in the degenerate solution results in structural changes in the wave profile. Moreover, the zero-frequency limit of the degenerate solution with differential shift results in a rogue wave triplet. The zero limit of the stretching factor in this solution, in turn, transforms the triplet into a singlet plus a low-amplitude soliton on the background. A large value of the differential shift parameter converts the triplet into a pure singlet.

A lossless nonlinear $LC$ transmission network is considered. With the use of the reductive perturbation method in the semidiscrete limit, we show that the dynamics of matter-wave solitons in the network can be modeled by a one-dimensional Gross-Pitaevskii (GP) equation with a time-dependent linear potential in the presence of a chemical potential. An explicit expression for the growth rate of a purely growing modulational instability (MI) is presented and analyzed. We find that the potential parameter of the GP equation of the system does not affect the different regions of the MI. Neglecting the chemical potential in the GP equation, we derive exact analytical solutions which describe the propagation of both bright and dark solitary waves on continuous-wave (cw) backgrounds. Using the found exact analytical solutions of the GP equation, we investigate numerically the transmission of both bright and dark solitary voltage signals in the network. Our numerical studies show that the amplitude of a bright solitary voltage signal and the depth of a dark solitary voltage signal as well as their width, their motion, and their behavior depend on (i) the propagation frequencies, (ii) the potential parameter, and (iii) the amplitude of the cw background. The GP equation derived in this paper with a time-dependent linear potential opens up different ideas that may be of considerable theoretical interest for the management of matter-wave solitons in nonlinear $LC$ transmission networks.

We present rogue wave solutions of the integrable nonlinear Schrödinger equation hierarchy with an infinite number of higher-order terms. The latter include higher-order dispersion and higher-order nonlinear terms. In particular, we derive the fundamental rogue wave solutions for all orders of the hierarchy, with exact expressions for velocities, phase, and "stretching factors" in the solutions. We also present several examples of exact solutions of second-order rogue waves, including rogue wave triplets.

The freak oscillation in one-dimensional dusty plasma is studied numerically by particle-in-cell method. Using a perturbation method, the basic set of fluid equations is reduced to a nonlinear Schrödinger equation (NLSE). The rational solution of the NLSE is presented, which is proposed as an effective tool for studying the rogue waves in dusty plasma. Additionally, the application scope of the analytical solution of the rogue wave described by the NLSE is given.

We present a numerical study of the cubic-quintic nonlinear Schr\"odinger equation in two transverse dimensions, relevant for the propagation of light in certain exotic media. A well known feature of the model is the existence of flat-top bright solitons of fixed intensity, whose dynamics resembles the physics of a liquid. They support traveling wave solutions, consisting of rarefaction pulses and vortex-antivortex pairs. In this work, we demonstrate how the vortex-antivortex pairs can be generated in bright soliton collisions displaying destructive interference followed by a snake instability. We then discuss the collisional dynamics of the dark excitations for different initial conditions. We describe a number of distinct phenomena including vortex exchange modes, quasielastic flyby scattering, soliton-like crossing, fully inelastic collisions and rarefaction pulse merging.

The effect of derivative nonlinearity and parity-time- (PT-) symmetric potentials on the wave propagation dynamics is investigated in the derivative nonlinear Schrodinger equation, where the physically interesting Scarff-II and hamonic-Hermite-Gaussian potentials are chosen. We study numerically the regions of unbroken/broken linear PT-symmetric phases and find some stable bright solitons of this model in a wide range of potential parameters even though the corresponding linear PT-symmetric phases are broken. The semi-elastic interactions between exact bright solitons and exotic incident waves are illustrated such that we find that exact nonlinear modes almost keep their shapes after interactions even if the exotic incident waves have evidently been changed. Moreover, we exert the adiabatic switching on PT-symmetric potential parameters such that a stable nonlinear mode with the unbroken linear PT-symmetric phase can be excited to another stable nonlinear mode belonging to the broken linear PT-symmetric phase., 11 pages, 11 figures

Explicit symplectic integrators have been important tools for accurate and efficient approximations of mechanical systems with separable Hamiltonians. For the first time, the article proposes for arbitrary Hamiltonians similar integrators, which are explicit, of any even order, symplectic in an extended phase space, and with pleasant long time properties. They are based on a mechanical restraint that binds two copies of phase space together. Using backward error analysis, KAM theory, and additional multiscale analysis, an error bound of $\mathcal{O}(T\delta^l \omega)$ is established for integrable systems, where $T$, $\delta$, $l$ and $\omega$ are respectively the (long) simulation time, step size, integrator order, and some binding constant. For non-integrable systems with positive Lyapunov exponents, such an error bound is generally impossible, but satisfactory statistical behaviors were observed in a numerical experiment with a nonlinear Schr\"{o}dinger equation., Comment: Accepted by Phys. Rev. E

Analysis of short-pulse propagation in positive dispersion media, e.g., in optical fibers and in shallow water, requires assorted high-order derivative terms. We present an infinite-order "dark" hierarchy of equations, starting from the basic defocusing nonlinear Schrödinger equation. We present generalized soliton solutions, plane-wave solutions, and periodic solutions of all orders. We find that "even"-order equations in the set affect phase and "stretching factors" in the solutions, while "odd"-order equations affect the velocities. Hence odd-order equation solutions can be real functions, while even-order equation solutions are complex. There are various applications in optics and water waves.

We study a variable-coefficient nonlinear Schr\"odinger (vc-NLS) equation with higher-order effects. We show that the breather solution can be converted into four types of nonlinear waves on constant backgrounds including the multi-peak solitons, antidark soliton, periodic wave and W-shaped soliton. The transition condition requiring the group velocity dispersion (GVD) and third-order dispersion (TOD) to scale linearly is obtained analytically. We display several kinds of elastic interactions between the transformed nonlinear waves. We discuss the dispersion management of multi-peak soliton, which indicates that the GVD coefficient controls the number of peaks of the wave while the TOD coefficient has compression effect. The gain or loss has influence on the amplitudes of the multi-peak soliton. We further derive the breather multiple births by using multiple compression points of Akhmediev breathers in optical fiber systems with periodic dispersion. The number of ABs depends on the amplitude of the modulation but not on its wavelength, which affects their separation distance. In the limiting case, the breather multiple births reduce to the Peregrine combs. We discuss the effects of TOD coefficient on the spatiotemporal characteristics of Peregrine combs. When the amplitude of the modulation is equal to 1, the Peregrine comb is converted into a Peregrine wall that can be seen as intermediate state between rogue wave and W-shaped soliton. We finally find that the modulational stability regions with zero growth rate coincide with the transition condition using rogue wave eigenvalues. Our results could be useful for the experimental control and manipulation of the formation of generalized Peregrine rogue waves in diverse physical systems modeled by vc-NLS equation with higher-order effects.

We consider accurate numerical solutions of the one-dimensional time-dependent Schrödinger equation when the potential is piecewise continuous. Spatial step sizes are defined for each of the regions between the discontinuities and a matching condition at the boundaries of the regions is employed. The Numerov method for spatial integration is particularly appropriate to this approach. By employing Padé approximants for the time-evolution operator, we obtain solutions with significantly improved precision without increased CPU time. This approach is also appropriate for adaptive changes in spatial step size even when there is no discontinuity of the potential.

PubMed ID: 27415263 In this paper we analyze the rogue wave spectra of the Kundu-Eckhaus equation (KEE). We compare our findings with their nonlinear Schrodinger equation (NLSE) analogs and show that the spectra of the individual rogue waves significantly differ from their NLSE analogs. A remarkable difference is the one-sided development of the triangular spectrum before the rogue wave becomes evident in time. Also we show that increasing the skewness of the rogue wave results in increased asymmetry in the triangular Fourier spectra. Additionally, the triangular spectra of the rogue waves of the KEE begin to develop at earlier stages of their development compared to their NLSE analogs, especially for larger skew angles. This feature may be used to enhance the early warning times of the rogue waves. However, we show that in a chaotic wave field with many spectral components the triangular spectra remain as the main attribute as a universal feature of the typical wave fields produced through modulation instability and characteristic features of the KEE's analytical rogue wave spectra may be suppressed in a realistic chaotic wave field. Publisher's Version

This article explores the excitation of different vibrational states in a spatially extended dynamical system through theory and experiment. As a prototypical example, we consider a one-dimensional packing of spherical particles (a so-called granular chain) that is subject to harmonic boundary excitation. The combination of the multimodal nature of the system and the strong coupling between the particles due to the nonlinear Hertzian contact force leads to broad regions in frequency where different vibrational states are possible. In certain parametric regions, we demonstrate that the nonlinear Schrödinger equation predicts the corresponding modes fairly well. The electromechanical model we apply predicts accurately the conversion from the obtained mechanical energy to the electrical energy observed in experiments.

We study the dynamics of self-trapping in Bose-Einstein condensates (BECs) loaded in deep optical lattices with Gaussian initial conditions, when the dynamics is well described by the discrete nonlinear Schrödinger equation (DNLSE). In the literature an approximate dynamical phase diagram based on a variational approach was introduced to distinguish different dynamical regimes: diffusion, self-trapping, and moving breathers. However, we find that the actual DNLSE dynamics shows a completely different diagram than the variational prediction. We calculate numerically a detailed dynamical phase diagram accurately describing the different dynamical regimes. It exhibits a complex structure that can readily be tested in current experiments in BECs in optical lattices and in optical waveguide arrays. Moreover, we derive an explicit theoretical estimate for the transition to self-trapping in excellent agreement with our numerical findings, which may be a valuable guide as well for future studies on a quantum dynamical phase diagram based on the Bose-Hubbard Hamiltonian.

We inquire how structure emerges during the process of protein folding. For this we scrutinise col- lective many-atom motions during all-atom molecular dynamics simulations. We introduce, develop and employ various topological techniques, in combination with analytic tools that we deduce from the concept of integrable models and structure of discrete nonlinear Schroedinger equation. The example we consider is an alpha-helical subunit of the HIV envelope glycoprotein gp41. The helical structure is stable when the subunit is part of the biological oligomer. But in isolation the helix becomes unstable, and the monomer starts deforming. We follow the process computationally. We interpret the evolving structure both in terms of a backbone based Heisenberg spin chain and in terms of a side chain based XY spin chain. We find that in both cases the formation of protein super-secondary structure is akin the formation of a topological Bloch domain wall along a spin chain. During the process we identify three individual Bloch walls and we show that each of them can be modelled with a very high precision in terms of a soliton solution to a discrete nonlinear Schroedinger equation., 20 pages 29 figures

In this paper we present a general framework for solving the stationary nonlinear Schr\"odinger equation (NLSE) on a network of one-dimensional wires modelled by a metric graph with suitable matching conditions at the vertices. A formal solution is given that expresses the wave function and its derivative at one end of an edge (wire) nonlinearly in terms of the values at the other end. For the cubic NLSE this nonlinear transfer operation can be expressed explicitly in terms of Jacobi elliptic functions. Its application reduces the problem of solving the corresponding set of coupled ordinary nonlinear differential equations to a finite set of nonlinear algebraic equations. For sufficiently small amplitudes we use canonical perturbation theory which makes it possible to extract the leading nonlinear corrections over large distances., Comment: 26 pages

The reduced dynamics for dark and bright soliton chains in the one-dimensional nonlinear Schrödinger equation is used to study the behavior of collective compression waves corresponding to Toda lattice solitons. We coin the term hypersoliton to describe such solitary waves riding on a chain of solitons. It is observed that in the case of dark soliton chains, the formulated reduction dynamics provides an accurate an robust evolution of traveling hypersolitons. As an application to Bose-Einstein condensates trapped in a standard harmonic potential, we study the case of a finite dark soliton chain confined at the center of the trap. When the central chain is hit by a dark soliton, the energy is transferred through the chain as a hypersoliton that, in turn, ejects a dark soliton on the other end of the chain that, as it returns from its excursion up the trap, hits the central chain repeating the process. This periodic evolution is an analog of the classical Newton's cradle.

The dynamics of solitons is considered in the framework of an extended nonlinear Schrödinger equation (NLSE), which is derived from a Zakharov-type model for wind-driven high-frequency surface waves in the ocean, coupled to damped low-frequency internal waves. The drive gives rise to a convective (but not absolute) instability in the system. The resulting NLSE includes a pseudo-stimulated-Raman-scattering (pseudo-SRS) term, which is a spatial-domain counterpart of the SRS term, a well-known ingredient of the temporal-domain NLSE in optics. Analysis of the field-momentum balance and direct simulations demonstrates that wave-number downshift by the pseudo-SRS may be compensated by the upshift induced by the wind traction, thus maintaining robust bright solitons in both stationary and oscillatory forms; in particular, they are not destroyed by the underlying convective instability. Analytical soliton solutions are found in an approximate form and are verified by numerical simulations. Solutions for soliton pairs are obtained in the numerical form.

In the present work, we combine the notion of parity-time (PT ) symmetry with that of supersymmetry (SUSY) for a prototypical case example with a complex potential that is related by SUSY to the so-called P¨oschl-Teller potential which is real. Not only are we able to identify and numerically confirm the eigenvalues of the relevant problem, but we also show that the corresponding nonlinear problem, in the presence of an arbitrary power-law nonlinearity, has an exact bright soliton solution that can be analytically identified and has intriguing stability properties, such as an oscillatory instability,which is absent for the corresponding solution of the regular nonlinear Schr¨odinger equation with arbitrary power-law nonlinearity. The spectral properties and dynamical implications of this instability are examined. We believe that these findings may pave the way toward initiating a fruitful interplay between the notions of PT symmetry, supersymmetric partner potentials, and nonlinear interactions. This work was supported in part by the US Department of Energy. A.K. wishes to thank the Indian National Science Academy (INSA) for financial support at Pune University. P.G.K. gratefully acknowledges support from NSF-DMS- 1312856, the Binational Science Foundation under Grant No. 2010239, and the ERC under FP7, Marie Curie Actions, People, International Research Staff Exchange Scheme (IRSES-605096)

We study integrable coupled nonlinear Schrödinger equations with pair particle transition between components. Based on exact solutions of the coupled model with attractive or repulsive interaction, we predict that some new dynamics of nonlinear excitations can exist, such as the striking transition dynamics of breathers, new excitation patterns for rogue waves, topological kink excitations, and other new stable excitation structures. In particular, we find that nonlinear wave solutions of this coupled system can be written as a linear superposition of solutions for the simplest scalar nonlinear Schrödinger equation. Possibilities to observe them are discussed in a cigar-shaped Bose-Einstein condensate with two hyperfine states. The results would enrich our knowledge on nonlinear excitations in many coupled nonlinear systems with transition coupling effects, such as multimode nonlinear fibers, coupled waveguides, and a multicomponent Bose-Einstein condensate system.

In this paper, a simple and constructive method is presented to find the generalized perturbation (n,M)-fold Darboux transformations (DTs) of the modified nonlinear Schrodinger (MNLS) equation in terms of fractional forms of determinants. In particular, we apply the generalized perturbation (1,N-1)-fold DTs to find its explicit multi-rogue-wave solutions. The wave structures of these rogue-wave solutions of the MNLS equation are discussed in detail for different parameters, which display abundant interesting wave structures, including the triangle and pentagon, etc., and may be useful to study the physical mechanism of multirogue waves in optics. The dynamical behaviors of these multi-rogue-wave solutions are illustrated using numerical simulations. The same Darboux matrix can also be used to investigate the Gerjikov-Ivanov equation such that its multi-rogue-wave solutions and their wave structures are also found. The method can also be extended to find multi-rogue-wave solutions of other nonlinear integrable equations.