Localized states in discrete nonlinear Schrödinger equations

Spatial structure of the non-integrable discrete defocusing Hirota equation

In this paper, we investigate the spatial property of the non-integrable discrete defocusing Hirota equation utilizing a planar nonlinear discrete dynamical map method. We construct the periodic orbit solutions of the stationary discrete defocusing Hirota equation. The behavior of the orbits in the vicinity of the special periodic solution is analyzed by taking advantage of the named residue. We characterize the effects of the parameters on the aperiodic orbits with the aid of numerical simulations. A comparison with the non-integrable discrete defocusing nonlinear Schrödinger equation case reveals that the non-integrable discrete defocusing Hirota equation has more abundant spatial properties. Rather an interesting and novel thing is that for any initial value, there exists triperiodic solutions for a reduced map.

Propagating intrinsic localized mode in a cyclic, dissipative, self-dual one-dimensional nonlinear transmission line

A well-known feature of a propagating localized excitation in a discrete lattice is the generation of a backwave in the extended normal mode spectrum. To quantify the parameter-dependent amplitude of such a backwave, the properties of a running intrinsic localized mode (ILM) in electric, cyclic, dissipative, nonlinear 1D transmission lines, containing balanced nonlinear capacitive and inductive terms, are studied via simulations. Both balanced and unbalanced damping and driving conditions are treated. The introduction of a unit cell duplex driver, with a voltage source driving the nonlinear capacitor and a synchronized current source, the nonlinear inductor, provides an opportunity to design a cyclic, dissipative self-dual nonlinear transmission line. When the self-dual conditions are satisfied, the dynamical voltage and current equations of motion within a cell become the same, the strength of the fundamental, resonant coupling between the ILM and the lattice modes collapses, and the associated fundamental backwave is no longer observed.

How close are integrable and nonintegrable models: A parametric case study based on the Salerno model

In the present work we revisit the Salerno model as a prototypical system that interpolates between a well-known integrable system (the Ablowitz-Ladik lattice) and an experimentally tractable, nonintegrable one (the discrete nonlinear Schrödinger model). The question we ask is, for "generic" initial data, how close are the integrable to the nonintegrable models? Our more precise formulation of this question is, How well is the constancy of formerly conserved quantities preserved in the nonintegrable case? Upon examining this, we find that even slight deviations from integrability can be sensitively felt by measuring these formerly conserved quantities in the case of the Salerno model. However, given that the knowledge of these quantities requires a deep physical and mathematical analysis of the system, we seek a more "generic" diagnostic towards a manifestation of integrability breaking. We argue, based on our Salerno model computations, that the full spectrum of Lyapunov exponents could be a sensitive diagnostic to that effect.

Thermalization in the one-dimensional Salerno model lattice

The Salerno model constitutes an intriguing interpolation between the integrable Ablowitz-Ladik (AL) model and the more standard (nonintegrable) discrete nonlinear Schrödinger (DNLS) one. The competition of local on-site nonlinearity and nonlinear dispersion governs the thermalization of this model. Here, we investigate the statistical mechanics of the Salerno one-dimensional lattice model in the nonintegrable case and illustrate the thermalization in the Gibbs regime. As the parameter interpolating between the two limits (from DNLS toward AL) is varied, the region in the space of initial energy and norm densities leading to thermalization expands. The thermalization in the non-Gibbs regime heavily depends on the finite system size; we explore this feature via direct numerical computations for different parametric regimes.

Nonautonomous discrete bright soliton solutions and interaction management for the Ablowitz-Ladik equation

We present the nonautonomous discrete bright soliton solutions and their interactions in the discrete Ablowitz-Ladik (DAL) equation with variable coefficients, which possesses complicated wave propagation in time and differs from the usual bright soliton waves. The differential-difference similarity transformation allows us to relate the discrete bright soliton solutions of the inhomogeneous DAL equation to the solutions of the homogeneous DAL equation. Propagation and interaction behaviors of the nonautonomous discrete solitons are analyzed through the one- and two-soliton solutions. We study the discrete snaking behaviors, parabolic behaviors, and interaction behaviors of the discrete solitons. In addition, the interaction management with free functions and dynamic behaviors of these solutions is investigated analytically, which have certain applications in electrical and optical systems.

Nonintegrable semidiscrete Hirota equation: gauge-equivalent structures and dynamical properties

In this paper, we investigate nonintegrable semidiscrete Hirota equations, including the nonintegrable semidiscrete Hirota(-) equation and the nonintegrable semidiscrete Hirota(+) equation. We focus on the topics on gauge-equivalent structures and dynamical behaviors for the two nonintegrable semidiscrete equations. By using the concept of the prescribed discrete curvature, we show that, under the discrete gauge transformations, the nonintegrable semidiscrete Hirota(-) equation and the nonintegrable semidiscrete Hirota(+) equation are, respectively, gauge equivalent to the nonintegrable generalized semidiscrete modified Heisenberg ferromagnet equation and the nonintegrable generalized semidiscrete Heisenberg ferromagnet equation. We prove that the two discrete gauge transformations are reversible. We study the dynamical properties for the two nonintegrable semidiscrete Hirota equations. The exact spatial period solutions of the two nonintegrable semidiscrete Hirota equations are obtained through the constructions of period orbits of the stationary discrete Hirota equations. We discuss the topic regarding whether the spatial period property of the solution to the nonintegrable semidiscrete Hirota equation is preserved to that of the corresponding gauge-equivalent nonintegrable semidiscrete equations under the action of discrete gauge transformation. By using the gauge equivalent, we obtain the exact solutions to the nonintegrable generalized semidiscrete modified Heisenberg ferromagnet equation and the nonintegrable generalized semidiscrete Heisenberg ferromagnet equation. We also give the numerical simulations for the stationary discrete Hirota equations. We find that their dynamics are much richer than the ones of stationary discrete nonlinear Schrödinger equations.

Mechanical energy fluctuations in granular chains: the possibility of rogue fluctuations or waves

The existence of rogue or freak waves in the ocean has been known for some time. They have been reported in the context of optical lattices and the financial market. We ask whether such waves are generic to late time behavior in nonlinear systems. In that vein, we examine the dynamics of an alignment of spherical elastic beads held within fixed, rigid walls at zero precompression when they are subjected to sufficiently rich initial conditions. Here we define such waves generically as unusually large energy fluctuations that sustain for short periods of time. Our simulations suggest that such unusually large fluctuations ("hot spots") and occasional series of such fluctuations through space and time ("rogue fluctuations") are likely to exist in the late time dynamics of the granular chain system at zero dissipation. We show that while hot spots are common in late time evolution, rogue fluctuations are seen in purely nonlinear systems (i.e., no precompression) at late enough times. We next show that the number of such fluctuations grows exponentially with increasing nonlinearity whereas rogue fluctuations decrease superexponentially with increasing precompression. Dissipation-free granular alignment systems may be possible to realize as integrated circuits and hence our observations may potentially be testable in the laboratory.

Extreme events in discrete nonlinear lattices

We perform statistical analysis on discrete nonlinear waves generated through modulational instability in the context of the Salerno model that interpolates between the integrable Ablowitz-Ladik (AL) equation and the nonintegrable discrete nonlinear Schrödinger equation. We focus on extreme events in the form of discrete rogue or freak waves that may arise as a result of rapid coalescence of discrete breathers or other nonlinear interaction processes. We find power law dependence in the wave amplitude distribution accompanied by an enhanced probability for freak events close to the integrable limit of the equation. A characteristic peak in the extreme event probability appears that is attributed to the onset of interaction of the discrete solitons of the AL equation and the accompanied transition from the local to the global stochasticity monitored through the positive Lyapunov exponent of a nonlinear map.

Nonlinear lattice model for spatially guided solitons in nonlinear photonic crystals

Numerical simulations have shown the existence of transversely localized guided modes in nonlinear two-dimensional photonic crystals. These soliton-like Bloch waves induce their own waveguide in a photonic crystal without the presence of a linear defect. By applying a Green's function method which is limited to within a strip perpendicular to the propagation direction, we are able to describe these Bloch modes by a nonlinear lattice model that includes the long-range site-to-site interaction between the scattered fields and the non-local nonlinear response of the photonic crystal. The advantages of this semi-analytical approach are discussed and a comparison with a rigorous numerical analysis is given in different configurations. Both monoatomic and diatomic nonlinear photonic crystals are considered.

Nonintegrable Schrodinger discrete breathers

In an extensive numerical investigation of nonintegrable translational motion of discrete breathers in nonlinear Schrödinger lattices, we have used a regularized Newton algorithm to continue these solutions from the limit of the integrable Ablowitz-Ladik lattice. These solutions are shown to be a superposition of a localized moving core and an excited extended state (background) to which the localized moving pulse is spatially asymptotic. The background is a linear combination of small amplitude nonlinear resonant plane waves and it plays an essential role in the energy balance governing the translational motion of the localized core. Perturbative collective variable theory predictions are critically analyzed in the light of the numerical results.

Spatiotemporal discrete multicolor solitons

We have found various families of two-dimensional spatiotemporal solitons in quadratically nonlinear waveguide arrays. The families of unstaggered odd, even, and twisted stationary solutions are thoroughly characterized and their stability against perturbations is investigated. We show that the twisted and even solitons display instability, while most of the odd solitons show remarkable stability upon evolution.

Perturbation-induced radiation by the Ablowitz-Ladik soliton

An efficient formalism is elaborated to analytically describe dynamics of the Ablowitz-Ladik soliton in the presence of perturbations. This formalism is based on using the Riemann-Hilbert problem and provides the means of calculating evolution of the discrete soliton parameters, as well as shape distortion and perturbation-induced radiation effects. As an example, soliton characteristics are calculated for linear damping and quintic perturbations.

Using the finite domain remnant of the continuous spectrum to examine integrability: effect of boundary conditions

The aim of this work is to propose a method for testing the integrability of a model partial differential (PDE) and/or differential difference equation (DDE), by examining it in a finite but large domain. For monoparametric families of PDE/DDE's, that are known to possess isolated integrable points, we find that very special features occur in the finite domain remnant of the continuous ("phonon") spectrum at these "singular" points. We identify these features in the case example of a PDE and a DDE (that sustain front and pulselike solutions, respectively) for different types of boundary conditions. The key finding of the work is that such spectral features are generic near the singular, integrable points and hence we propose to explore a given PDE/DDE in a finite but large domain for such traits, as a means of assessing its potential integrability.

Study of intrinsic localized vibrational modes in micromechanical oscillator arrays

Intrinsic localized modes (ILMs) have been observed in micromechanical cantilever arrays, and their creation, locking, interaction, and relaxation dynamics in the presence of a driver have been studied. The micromechanical array is fabricated in a 300 nm thick silicon-nitride film on a silicon substrate, and consists of up to 248 cantilevers of two alternating lengths. To observe the ILMs in this experimental system a line-shaped laser beam is focused on the 1D cantilever array, and the reflected beam is captured with a fast charge coupled device camera. The array is driven near its highest frequency mode with a piezoelectric transducer. Numerical simulations of the nonlinear Klein-Gordon lattice have been carried out to assist with the detailed interpretation of the experimental results. These include pinning and locking of the ILMs when the driver is on, collisions between ILMs, low frequency excitation modes of the locked ILMs and their relaxation behavior after the driver is turned off.

Continuum approach to discreteness

We study analytically and numerically continuum models derived on the basis of Padé approximations and their effectiveness in modeling spatially discrete systems. We not only analyze features of the temporal dynamics that can be captured through these continuum approaches (e.g., shape oscillations, radiation effects, and trapping) but also point out ones that cannot be captured (such as Peierls-Nabarro barriers and Bloch oscillations). We analyze the role of such methods in providing an effective "homogenization" of spatially discrete, as well as of heterogeneous continuum equations. Finally, we develop numerical methods for solving such equations and use them to establish the range of validity of these continuum approximations, as well as to compare them with other semicontinuum approximations.

Methods for discrete solitons in nonlinear lattices

A method to find discrete solitons in nonlinear lattices is introduced. Using nonlinear optical waveguide arrays as a prototype application, both stationary and traveling-wave solitons are investigated. In the limit of small wave velocity, a fully discrete perturbative analysis yields formulas for the mode shapes and velocity.

Blocking and routing discrete solitons in two-dimensional networks of nonlinear waveguide arrays

It is shown that discrete solitons can be navigated in two-dimensional networks of nonlinear waveguide arrays. This can be accomplished via vector interactions between two classes of discrete solitons: signals and blockers. Discrete solitons in such two-dimensional networks can exhibit a rich variety of functional operations, e.g., blocking, routing, logic functions, and time gating.

Discrete solitons and breathers with dilute Bose-Einstein condensates

We study the dynamical phase diagram of a dilute Bose-Einstein condensate (BEC) trapped in a periodic potential. The dynamics is governed by a discrete nonlinear Schrödinger equation: intrinsically localized excitations, including discrete solitons and breathers, can be created even if the BEC's interatomic potential is repulsive. Furthermore, we analyze the Anderson-Kasevich experiment [Science 282, 1686 (1998)], pointing out that mean field effects lead to a coherent destruction of the interwell Bloch oscillations.

Twisted localized modes

In a number of recent papers the so-called twisted localized mode of the discrete nonlinear Schrödinger equation has been proposed. Herein, we study the existence and stability properties of such modes. We analyze the persistence of quasiperiodic modes and study the domains of existence and numerical stability of the exact form of such solutions. We identify the bifurcations through which they lose their stability and follow the behavior of the intrinsic localized modes and their eigenmodes even in the unstable regime.

Electron-vibron-breather interaction

We study the interaction of breathers in the context of a coupled electron-vibron lattice system. Starting with single-site excitations, it is demonstrated that constellations exist for which the coexistence of electronic and vibronic breathers is assured. The energy exchange between the vibrational and electronic subsystems and its impact on the breather formation are discussed in detail. The coupled electron-vibron dynamics shows a tendency toward energy redistribution into the vibronic degrees of freedom at the expense of the electronic energy content. Attention is paid to the relaxation dynamics in the energy exchange and we discuss the attainment of a steady regime for the coupled electron-vibron dynamics starting from a nonequilibrium state. It is demonstrated that the presence of breathers has a strong impact on the relaxation dynamics. Breathers can assist the relaxation process. With the help of a linear stability analysis, we show why the electronic subsystem acts as an energy donor while the vibron system serves as the energy acceptor. To this end we investigate the existence and stability of localized breathing eigenmodes capable of energy trapping. A frequency analysis reveals that strong exchange also occurs due to a temporal transition from single-frequency breathers to those oscillating with two frequencies and their temporal resonance interaction. Finally, the self-stabilized electron-vibron system relaxes to a combined electron-vibron breather. On increasing the electron-vibron coupling strength, only a vibronic phonobreather of large amplitude survives, whereas the electronic subsystem tends to energy equipartition.

Perturbative study of classical Ablowitz-Ladik type soliton dynamics in relation to energy transport in alpha-helical proteins

Classical Ablowitz-Ladik type soliton dynamics from three closely related classical nonlinear equations is studied using a perturbative method. Model nonintegrable equations are derived by assuming nearest neighbor hopping of an exciton(vibron) in the presence of a full exciton(vibron)-phonon interaction in soft molecular chains in general and spines of alpha-helices in particular. In all cases, both trapped and moving solitons are found implying activation energy barrier for propagating solitons. Analysis further shows that staggered and nearly staggered trapped solitons will have a negative effective mass. In some models the exciton(vibron)-phonon coupling affects the hopping. For these models, when the conservation of probability is taken into account, only propagating solitons with a broad profile are found to be acceptable solutions. Of course, for the soliton to be a physically meaningful entity, total nonlinear coupling strength should exceed a critical value. On the basis of the result, a plausible modification in the mechanism for biological energy transport involving conformational change in alpha-helix is proposed. Future directions of the work are also mentioned.

Solitonic energy transfer in a coupled exciton-vibron system

We consider the exciton transfer along a one-dimensional molecular chain. The exciton motion is influenced by longitudinal vibrations evolving in a Toda lattice potential. It is shown how the soliton solutions of the vibron system coupled to the exciton system induce solitonic exciton transfer. To this aim the existence of a regime of suppressed energy exchange between the coupled excitonic and vibrational degrees of freedom is established in the case of which a nonlinear Schrodinger equation for the exciton variable is derived. The nonlinear Schrodinger equation possesses soliton solutions corresponding to coherent transfer of the localized exciton.