Search

We investigate time-independent solutions of a discrete optical cavity model featuring saturable Kerr nonlinearity, a discrete version of the Lugiato-Lefever equation. This model supports continuous wave (uniform) and localized (discrete soliton) solutions. Stationary bright solitons arise through the interaction of dark and bright uniform states, forming a homoclinic snaking bifurcation diagram within the Pomeau pinning region. As the system approaches the anti-continuum limit (weak coupling), this snaking bifurcation widens and transitions into $\subset$-shaped isolas. We propose a one-active-site approximation that effectively captures the system's behavior in this regime. The approximation also provides insight into the stability properties of soliton states. Numerical continuation and spectral analysis confirm the accuracy of this semianalytical method, showing excellent agreement with the full model., Comment: published more...

This study investigates the existence and stability of localized states in the discrete nonlinear Schr\"odinger (DNLS) equation with quadratic and cubic nonlinearities, describing the so-called quantum droplets and bubbles. Those states exist within an interval known as the pinning region, as we vary a control parameter. Within the interval, multistable states are connected through multiple hysteresis, called homoclinic snaking. In particular, we explore its mechanism and consider two limiting cases of coupling strength: weak (anti-continuum) and strong (continuum) limits. We employ an asymptotic and a variational method for the weak and strong coupling limits, respectively, to capture the pinning region's width. The width exhibits an algebraic and an exponentially small dependence on the coupling constant for the weak and strong coupling, respectively. This finding is supported by both analytical and numerical results, which show excellent agreement. We also consider the modulational instability of spatially uniform solutions. Our work sheds light on the intricate interplay between multistability and homoclinic snaking in discrete quantum systems, paving the way for further exploration of complex nonlinear phenomena in this context., Comment: Chaos, Solitons and Fractals more...

In this work, we study a model of a one-dimensional magnetic metamaterial formed by a discrete array of nonlinear resonators. We focus on periodic and localized traveling waves of the model, in the presence of loss and an external drive. Employing a Melnikov analysis we study the existence and persistence of such traveling waves, and study their linear stability. We show that, under certain conditions, the presence of dissipation and/or driving may stabilize or destabilize the solutions. Our analytical results are found to be in good agreement with direct numerical computations., Comment: 12 pages more...

We consider the discrete Allen-Cahn equation with cubic and quintic nonlinearity on the Lieb lattice. We study localized nonlinear solutions of the system that have linear multistability and hysteresis in their bifurcation diagram. In this work, we investigate the system's homoclinic snaking, i.e. snaking-like structure of the bifurcation diagram, particularly the effect of the lattice type. Numerical continuation using a pseudoarclength method is used to obtain localized solutions along the bifurcation diagram. We then develop an active-cell approximation to classify the type of solution at the turning points, which gives good agreement with the numerical results when the sites are weakly coupled. Time dynamics of localized solutions inside and outside the pinning region is also discussed., Comment: published more...

We study continuations of topological edge states in the Su-Schrieffer-Heeger model with on-site cubic (Kerr) nonlinearity, which is a 1D nonlinear photonic topological insulator (TI). Based on the topology of the underlying spatial dynamical system, we establish the existence of nonlinear edge states (edge solitons) for all positive energies in the topological band gap. We discover that these edge solitons are stable at any energy when the ratio between the weak and strong couplings is below a critical value. Above the critical coupling ratio, there are energy intervals where the edge solitons experience an oscillatory instability. Though our paper focuses on a photonic system, we also discuss the broader relevance of our methods and results to 1D nonlinear mechanical TIs., Comment: 7 pages, 6 figures; accepted version in Physical Review E after addressing referee comments more...

We consider a damped, parametrically driven discrete nonlinear Klein-Gordon equation, that models coupled pendula and micromechanical arrays, among others. To study the equation, one usually uses a small-amplitude wave ansatz, that reduces the equation into a discrete nonlinear Schr\"odinger equation with damping and parametric drive. Here, we justify the approximation by looking for the error bound with the method of energy estimates. Furthermore, we prove the local and global existence of {solutions to the discrete nonlinear} Schr\"odinger equation. To illustrate the main results, we consider numerical simulations showing the dynamics of errors made by the discrete nonlinear equation. We consider two types of initial conditions, with one of them being a discrete soliton of the nonlinear Schr\"odinger equation, that is expectedly approximate discrete breathers of the nonlinear Klein-Gordon equation., Comment: arXiv admin note: text overlap with arXiv:1911.01631 more...

We present a study of time-independent solutions of the two-dimensional discrete Allen-Cahn equation with cubic and quintic nonlinearity. Three different types of lattices are considered, i.e., square, honeycomb, and triangular lattices. The equation admits uniform and localised states. We can obtain localised solutions by combining two different states of uniform solutions, which can develop a snaking structure in the bifurcation diagrams. We find that the complexity and width of the snaking diagrams depend on the number of "patch interfaces" admitted by the lattice systems. We introduce an active-cell approximation to analyse the saddle-node bifurcation and stabilities of the corresponding solutions along the snaking curves. Numerical simulations show that the active-cell approximation gives good agreement for all of the lattice types when the coupling is weak. We also consider planar fronts that support our hypothesis on the relation between the complexity of a bifurcation diagram and the number of interface of its corresponding solutions., Comment: To appear in Nonlinearity more...

We consider the interactions of traveling localized wave solutions with a vertex in a star graph domain that describes multiple Josephson junctions with a common/branch point (i.e., tricrystal junctions). The system is modeled by the sine-Gordon equation. The vertex is represented by boundary conditions that are determined by the continuity of the magnetic field and vanishing total fluxes. When one considers small-amplitude breather solutions, the system can be reduced into the nonlinear Schr\"odinger equation posed on a star graph. Using the equation, we show that a high-velocity incoming soliton is split into a transmitted component and a reflected one. The transmission is shown to be in good agreement with the transmission rate of plane waves in the linear Schr\"odinger equation on the same graph (i.e., a quantum graph). In the context of the sine-Gordon equation, small-amplitude breathers show similar qualitative behaviors, while large-amplitude ones produce complex dynamics., Comment: 10 pages, 4 figures more...

We study static nonlinear waves in networks described by a nonlinear Schrodinger equation with point-like nonlinearities on metric graphs. Explicit solutions fulfilling vertex boundary conditions are obtained. Spontaneous symmetry breaking caused by bifurcations is found. more...

In this work, we consider a semi-infinite discrete nonlinear Schr\"odinger equation with saturable nonlinearity driven at one edge by a driving force. The equation models the dynamics of coupled photorefractive waveguide arrays. It has been reported that when the frequency of the driving force is in the forbidden band, energy can be transmitted along the lattices provided that the driving amplitude is above a critical value. This nonlinear tunneling is called supratransmission. Here, we explain the source of supratransmission using geometric illustrations. Approximations to the critical amplitude for supratransmission are presented as well., Comment: 7 pages, 2 figures, presented at the 6th International Symposium on Modern Optics and Its Applications (ISMOA) 2007 more...

We consider an array of dual-core waveguides, which represent an optical realization of a chain of dimers, with an active (gain-loss) coupling between the cores, opposite signs of discrete diffraction in the parallel arrays, and a phase-velocity mismatch between them (which is necessary for the stability of the system). The array provides an optical emulation of the charge-parity ($\mathcal{CP}$) symmetry. The addition of the intracore cubic nonlinearity gives rise to several species of fundamental discrete solitons, which exist in continuous families, although the system is non-Hermitian. The existence and stability of the soliton families are explored by means of analytical and numerical methods. An asymptotic analysis is presented for the case of weak intersite coupling (i.e., near the anticontinuum limit), as well as weak coupling between cores in each dimer. Several families of fundamental discrete solitons are found in the semi-infinite gap of the system's spectrum, that have no counterparts in the continuum limit, as well as a branch which belongs to the finite band gap and carries over into a family of stable gap solitons in that limit. One branch develops an oscillatory instability above a critical strength of the intersite coupling, others being stable in their entire existence regions. Unlike solitons in conservative lattices, which are controlled solely by the strength of the intersite coupling, here fundamental-soliton families have several control parameters, one of which, viz., the coefficient of the intercore coupling in the active host medium, may be readily adjusted in the experiment by varying the gain applied to the medium., Comment: 11 pages, 8 figures more...

We consider linearly coupled discrete nonlinear Schr\"odinger equations with gain and loss terms and with a cubic-quintic nonlinearity. The system models a parity-time ($\cal{PT}$)-symmetric coupler composed by a chain of dimers. Particularly we study site-centered and bond-centered spatially-localized solutions and present that each solution has a symmetric and antisymmetric configuration between the arms. When a parameter is varied, the resulting bifurcation diagrams for the existence of standing localized solutions have a snaking behaviour. The critical gain/loss coefficient above which the $\cal{PT}-$symmetry is broken corresponds to the condition when bifurcation diagrams of symmetric and antisymmetric states merge. Past the symmetry breaking, the system no longer has time-independent states. Nevertheless, equilibrium solutions can be analytically continued by defining a dual equation that leads to so-called ghost states associated with growth or decay, that are also identified and examined here. We show that ghost localized states also exhibit snaking bifurcation diagrams. We analyse the width of the snaking region and provide asymptotic approximations in the limit of strong and weak coupling where good agreement is obtained., Comment: To appear in PRE more...

In this work we investigate a one-dimensional parity-time (PT)-symmetric magnetic metamaterial consisting of split-ring dimers having gain or loss. Employing a Melnikov analysis we study the existence of localized travelling waves, i.e. homoclinic or heteroclinic solutions. We find conditions under which the homoclinic or heteroclinic orbits persist. Our analytical results are found to be in good agreement with direct numerical computations. For the particular nonlinearity admitting travelling kinks, numerically we observe homoclinic snaking in the bifurcation diagram. The Melnikov analysis yields a good approximation to one of the boundaries of the snaking profile., Comment: 15 pages, 4 figures, Wave Motion 76 (2018) 9-18. arXiv admin note: text overlap with arXiv:1304.0556 by other authors more...

An array of a finite number waveguides, driven laterally by injecting light at the outer waveguides, is considered. The array is modeled by a discrete nonlinear Schr\"{o}dinger equation. It has been shown [Phys. Rev. Lett. 94, 243902 (2005)] that, when the injected light is in the proximity of a bifurcation point, such a system can be sensitive to small disturbances, making it possible to act as a light detector. Here, the optimum intensity of the injected light is discussed, and an analytical approximation is presented., Comment: 3 pages, 2 figures more...

We investigate the flow of a one-dimensional nonlinear Schrodinger model with periodic boundary conditions past an obstacle, motivated by recent experiments with Bose--Einstein condensates in ring traps. Above certain rotation velocities, localized solutions with a nontrivial phase profile appear. In striking difference from the infinite domain, in this case there are many critical velocities. At each critical velocity, the steady flow solutions disappear in a saddle-center bifurcation. These interconnected branches of the bifurcation diagram lead to additions of circulation quanta to the phase of the associated solution. This, in turn, relates to the manifestation of persistent current in numerous recent experimental and theoretical works, the connections to which we touch upon. The complex dynamics of the identified waveforms and the instability of unstable solution branches are demonstrated., Comment: Resubmitted (2016) more...

A Parity-Time (PT)-symmetric system with periodically varying-in-time gain and loss modeled by two coupled Schrodinger equations (dimer) is studied. It is shown that the problem can be reduced to a perturbed pendulum-like equation. This is done by finding two constants of motion. Firstly, a generalized problem using Melnikov type analysis and topological degree arguments is studied for showing the existence of periodic (libration), shift periodic (rotation), and chaotic solutions. Then these general results are applied to the PT-symmetric dimer. It is interestingly shown that if a sufficient condition is satisfied, then rotation modes, which do not exist in the dimer with constant gain-loss, will persist. An approximate threshold for PT-broken phase corresponding to the disappearance of bounded solutions is also presented. Numerical study is presented accompanying the analytical results., Comment: submitted (2014) more...

We study collisions of moving nonlinear-Schr\"{o}dinger solitons with a $\mathcal{PT}$-symmetric dipole embedded into the one-dimensional self-focusing or defocusing medium. Accurate analytical results are produced for bright solitons, and, in a more qualitative form, for dark ones. In the former case, an essential aspect of the approximation is that it must take into regard the intrinsic chirp of the soliton, thus going beyond the framework of the simplest quasi-particle description of the soliton's dynamics. Critical velocities separating reflection and transmission of the incident bright solitons are found by means of numerical simulations, and in the approximate semi-analytical form. An exact solution for the dark soliton pinned by the complex $\mathcal{PT}$-symmetric dipole is produced too., Comment: 12 pages, 5 figures more...

Two effectively one-dimensional parallel coupled Bose-Einstein condensates in the presence of external potentials are studied. The system is modelled by linearly coupled Gross-Pitaevskii equations. In particular, the interactions of grey-soliton-like solutions representing analogues of superconducting Josephson fluxons as well as coupled dark solitons are discussed. A theoretical approximation based on variational formulations to calculate the oscillation frequency of the grey-soliton-like solution is derived and a qualitatively good agreement is obtained., Comment: 26 pages, 17 figures. Appeared in "Spontaneous Symmetry Breaking, Self-Trapping, and Josephson Oscillations", B.A. Malomed, ed. (Springer, Berlin, 2013) more...

We consider propagating, spatially localised waves in a class of equations that contain variational and non-variational terms. The dynamics of the waves is analysed through a collective coordinate approach. Motivated by the variational approximation, we show that there is a natural choice of projection onto collective variables for reducing the governing (nonlinear) partial differential equation (PDE) to coupled ordinary differential equations (ODEs). This projection produces ODEs whose solutions are exactly the stationary states of the effective Lagrangian that would be considered from applying the variational approximation method. We illustrate our approach by applying it to a modified Fisher equation for a travelling front, containing a non-constant-coefficient nonlinear term. We present numerical results that show that our proposed projection captures both the equilibria and the dynamics of the PDE much more closely than previously proposed projections., Comment: To appear in PRE (2013) more...

We study numerically a parametrically driven discrete nonlinear Schr\"odinger equation modelling periodically curved waveguide arrays. We show that discrete surface solitons persist, but their threshold power is altered by the drive. There are critical drives at which the threshold values vanish. We also show that parametric drives can create resonance with a phonon making a new barrier for discrete solitons. By calculating the corresponding Floquet multipliers, we find that the stability of symmetric and antisymmetric off-side discrete surface solitons switches approximately at the critical drives for thresholdless solitons. more...

We consider a parametrically driven damped discrete nonlinear Schr\"odinger (PDDNLS) equation. Analytical and numerical calculations are performed to determine the existence and stability of fundamental discrete bright solitons. We show that there are two types of onsite discrete soliton, namely onsite type I and II. We also show that there are four types of intersite discrete soliton, called intersite type I, II, III, and IV, where the last two types are essentially the same, due to symmetry. Onsite and intersite type I solitons, which can be unstable in the case of no dissipation, are found to be stabilized by the damping, whereas the other types are always unstable. Our further analysis demonstrates that saddle-node and pitchfork (symmetry-breaking) bifurcations can occur. More interestingly, the onsite type I, intersite type I, and intersite type III-IV admit Hopf bifurcations from which emerge periodic solitons (limit cycles). The continuation of the limit cycles as well as the stability of the periodic solitons are computed through the numerical continuation software Matcont. We observe subcritical Hopf bifurcations along the existence curve of the onsite type I and intersite type III-IV. Along the existence curve of the intersite type I we observe both supercritical and subcritical Hopf bifurcations., Comment: to appear in "Spontaneous Symmetry Breaking, Self-Trapping, and Josephson Oscillations in Nonlinear Systems", B.A. Malomed, ed. (Springer, Berlin, 2012) more...

A nonlinear Schr\"odinger equation with repulsive (defocusing) nonlinearity is considered. As an example, a system with a spatially varying coefficient of the nonlinear term is studied. The nonlinearity is chosen to be repelling except on a finite interval. Localized standing wave solutions on a non-zero background, e.g., dark solitons trapped by the inhomogeneity, are identified and studied. A novel instability criterion for such states is established through a topological argument. This allows instability to be determined quickly in many cases by considering simple geometric properties of the standing waves as viewed in the composite phase plane. Numerical calculations accompany the analytical results., Comment: 20 pages, 11 figures more...

In the present work, we generalize earlier considerations for intrinsic localized modes consisting of a few excited sites, as developed in the one-component discrete nonlinear Schrodinger equation model, to the case of two-component systems. We consider all the different combinations of "up" (zero phase) and "down" ({\pi} phase) site excitations and are able to compute not only the corresponding existence curves, but also the eigenvalue dependences of the small eigenvalues potentially responsible for instabilities, as a function of the nonlinear parameters of the model representing the self/cross phase modulation in optics and the scattering length ratios in the case of matter waves in optical lattices. We corroborate these analytical predictions by means of direct numerical computations. We infer that all the modes which bear two adjacent nodes with the same phase are unstable in the two component case and the only solutions that may be linear stable are ones where each set of adjacent nodes, in each component is out of phase. more...

A periodically inhomogeneous Schrodinger equation is considered. The inhomogeneity is reflected through a non-uniform coefficient of the linear and non-linear term in the equation. Due to the periodic inhomogeneity of the linear term, the system may admit spectral bands. When the oscillation frequency of a localized solution resides in one of the finite band gaps, the solution is a gap soliton, characterized by the presence of infinitely many zeros in the spatial profile of the soliton. Recently, how to construct such gap solitons through a composite phase portrait is shown. By exploiting the phase-space method and combining it with the application of a topological argument, it is shown that the instability of a gap soliton can be described by the phase portrait of the solution. Surface gap solitons at the interface between a periodic inhomogeneous and a homogeneous medium are also discussed. Numerical calculations are presented accompanying the analytical results., Comment: Comments and suggestions are welcome more...

Travelling solitary waves in the one-dimensional discrete nonlinear Schr\"{o}dinger equation (DNLSE) with saturable onsite nonlinearity are studied. A variational approximation (VA) for the solitary waves is derived in an analytical form. The stability is also studied by means of the VA, demonstrating that the solitons are stable, which is consistent with previously published results. Then, the VA is applied to predict parameters of travelling solitons with non-oscillatory tails (\textit{embedded solitons}, ESs). Two-soliton bound states are considered too. The separation distance between the solitons forming the bound state is derived by means of the VA. A numerical scheme based on the discretization of the equation in the moving coordinate frame is derived and implemented using the Newton--Raphson method. In general, a good agreement between the analytical and numerical results is obtained. In particular, we demonstrate the relevance of the analytical prediction of characteristics of the embedded solitons., Comment: This paper is to appear in Journal of Physics A more...

Localised structures appear in a wide variety of systems, arising from a pinning mechanism due to the presence of a small-scale pattern or an imposed grid. When there is a separation of lengthscales, the width of the pinning region is exponentially small and beyond the reach of standard asymptotic methods. We show how this behaviour can be obtained using a variational method, for two systems. In the case of the quadratic-cubic Swift-Hohenberg equation, this gives results that are in agreement with recent work using exponential asymptotics. Secondly, the method is applied to a discrete system with cubic-quintic nonlinearity, giving results that agree well with numerical simulations., Comment: submitted. Comments are welcome more...

We introduce a system of two linearly coupled discrete nonlinear Schr\"{o}dinger equations (DNLSEs), with the coupling constant subject to a rapid temporal modulation. The model can be realized in bimodal Bose-Einstein condensates (BEC). Using an averaging procedure based on the multiscale method, we derive a system of averaged (autonomous) equations, which take the form of coupled DNLSEs with additional nonlinear coupling terms of the four-wave-mixing type. We identify stability regions for fundamental onsite discrete symmetric solitons (single-site modes with equal norms in both components), as well as for two-site in-phase and twisted modes, the in-phase ones being completely unstable. The symmetry-breaking bifurcation, which destabilizes the fundamental symmetric solitons and gives rise to their asymmetric counterparts, is investigated too. It is demonstrated that the averaged equations provide a good approximation in all the cases. In particular, the symmetry-breaking bifurcation, which is of the pitchfork type in the framework of the averaged equations, corresponds to a Hopf bifurcation in terms of the original system., Comment: 6 pages, 3 figures more...

We investigate the snaking of localised patterns, seen in numerous physical applications, using a variational approximation. This method naturally introduces the exponentially small terms responsible for the snaking structure, that are not accessible via standard multiple-scales asymptotic techniques. We obtain the symmetric snaking solutions and the asymmetric 'ladder' states, and also predict the stability of the localised states. The resulting approximate formulas for the width of the snaking region show good agreement with numerical results., Comment: 4 pages, 3 figures, submitted more...

A long Josephson junction containing regions with a phase shift of pi is considered. By exploiting the defect modes due to the discontinuities present in the system, it is shown that Josephson junctions with phase-shift can be an ideal setting for studying localized mode interactions. A phase-shift configuration acting as a double-well potential is considered and shown to admit mode tunnelings between the wells. When the phase-shift configuration is periodic, it is shown that localized excitations forming bright and dark solitons can be created. Multi-mode approximations are derived confirming the numerical results., Comment: 4 pages, to appear in Phys. Rev. B more...

We study the dynamics of matter waves in an effectively one-dimensional Bose-Einstein condensate in a double well potential. We consider in particular the case when one of the double wells confines excited states. Similarly to the known ground state oscillations, the states can tunnel between the wells experiencing the physics known for electrons in a Josephson junction, or be self-trapped. As the existence of dark solitons in a harmonic trap are continuations of such non-ground state excitations, one can view the Josephson-like oscillations as tunnelings of dark solitons. Numerical existence and stability analysis based on the full equation is performed, where it is shown that such tunneling can be stable. Through a numerical path following method, unstable tunneling is also obtained in different parameter regions. A coupled-mode system is derived and compared to the numerical observations. Regions of (in)stability of Josephson tunneling are discussed and highlighted. Finally, we outline an experimental scheme designed to explore such dark soliton dynamics in the laboratory., Comment: submitted more...

A nonlinear Schrodinger equation arising from light propagation down an inhomogeneous medium is considered. The inhomogeneity is reflected through a non-uniform coefficient of the non-linear term in the equation. In particular, a combination of self-focusing and self-defocusing nonlinearity, with the self-defocusing region localized in a finite interval, is investigated. Using numerical computations, the extension of linear eigenmodes of the corresponding linearized system into nonlinear states is established, particularly nonlinear continuations of the fundamental state and the first excited state. The (in)stability of the states is also numerically calculated, from which it is obtained that symmetric nonlinear solutions become unstable beyond a critical threshold norm. Instability of the symmetric states is then investigated analytically through the application of a topological argument. Determination of instability of positive symmetric states is reduced to simple geometric properties of the composite phase plane orbit of the standing wave. Further the topological argument is applied to higher excited states and instability is again reduced to straightforward geometric calculations. For a relatively high norm, it is observed that asymmetric states bifurcate from the symmetric ones. The stability and instability of asymmetric states is also considered., Comment: Suggestions and comments are welcome more...

We consider a parametrically driven Klein--Gordon system describing micro- and nano-devices, with integrated electrical and mechanical functionality. Using a multiscale expansion method we reduce the system to a discrete nonlinear Schrodinger equation. Analytical and numerical calculations are performed to determine the existence and stability of fundamental bright and dark discrete solitons admitted by the Klein--Gordon system through the discrete Schrodinger equation. We show that a parametric driving can not only destabilize onsite bright solitons, but also stabilize intersite bright discrete solitons and onsite and intersite dark solitons. Most importantly, we show that there is a range of values of the driving coefficient for which dark solitons are stable, for any value of the coupling constant, i.e. oscillatory instabilities are totally suppressed. Stability windows of all the fundamental solitons are presented and approximations to the onset of instability are derived using perturbation theory, with accompanying numerical results. Numerical integrations of the Klein--Gordon equation are performed, confirming the relevance of our analysis. more...

In the present work, we propose a new set of coherent structures that arise in nonlinear dynamical lattices with more than one components, namely interlaced solitons. These are waveforms in which in the relevant anti-continuum limit, i.e. when the sites are uncoupled, one component has support where the other component does not. We illustrate systematically how one can combine dynamically stable unary patterns to create ones such for the binary case of two-components. In the one-dimensional setting, we provide also a detailed theoretical analysis of the existence and stability of these waveforms, while in higher dimensions, where such analytical computations are far more involved, we resort to corresponding numerical computations. Lastly, we perform direct numerical simulations to showcase how these structures break up, when exponentially or oscillatorily unstable, to structures with a smaller number of participating sites., Comment: 11 pages, 10 figures more...

In this report, we consider a semi-infinite discrete nonlinear Schr\"odinger equation driven at one edge by a driving force. The equation models the dynamics of coupled waveguide arrays. When the frequency of the forcing is in the allowed-band of the system, there will be a linear transmission of energy through the lattice. Yet, if the frequency is in the upper forbidden band, then there is a critical driving amplitude for a nonlinear tunneling, which is called supratransmission, of energy to occur. Here, we discuss mathematically the mechanism and the source of the supratransmission. By analyzing the existence and the stability of the rapidly decaying static discrete solitons of the system, we show rigorously that two of the static solitons emerge and disappear in a saddle-node bifurcation at a critical driving amplitude. One of the emerging solitons is always stable in its existence region and the other is always unstable. We argue that the critical amplitude for supratransmission is then the same as the critical driving amplitude of the saddle-node bifurcation. We consider as well the case of the forcing frequency in the lower forbidden band. It is discussed briefly that there is no supratransmission because in this case there is only one rapidly decaying static soliton that exists and is stable for any driving amplitude., Comment: Some figures are removed due to file size limit more...

Solitons and necklaces in the first band-gap of a two-dimensional optically induced honeycomb defocusing photonic lattice are theoretically considered. It is shown that dipoles, soliton necklaces, and vortex necklaces exist and may possess regions of stable propagation through a photorefractive crystal. Most of the configurations disappear in bifurcations close to the upper edge of the first band. Solutions associated with such bifurcations are also numerically examined, and it is found that they are often asymmetric and more exotic. The dynamics of the relevant unstable structures are also examined through direct numerical simulations revealing either breathing oscillations or, in some cases, destruction of the original waveform. more...

Dipole and quadrupole solitons in a two-dimensional optically induced defocus- ing photonic lattice are theoretically predicted and experimentally observed. It is shown that in-phase nearest-neighbor dipole and out-of-phase next-nearest-neighbor dipoles exist and can be stable in the intermediate intensity regime. The other types of dipoles are always unstable. In-phase nearest-neighbor quadrupoles are also numerically obtained, and may also be linearly stable. Out-of-phase, nearest-neighbor quadrupoles are found to be typically unstable. These numerical results are found to be aligned with the main predictions obtained analytically in the discrete nonlinear Schroedinger model. Finally, experimental results are presented for both dipole and quadrupole structures, indicating that self-trapping of such structures in the defocusing lattice can be realized for the length of the nonlinear crystal (10 mm). more...

We consider a spatially non-autonomous discrete sine-Gordon equation with constant forcing and its continuum limit(s) to model a 0-$\pi$ Josephson junction with an applied bias current. The continuum limits correspond to the strong coupling limit of the discrete system. The non-autonomous character is due to the presence of a discontinuity point, namely a jump of $\pi$ in the sine-Gordon phase. The continuum models admits static solitary waves which are called $\pi$-kinks and are attached to the discontinuity point. For small forcing, there are three types of $\pi$-kinks. We show that one of the kinks is stable and the others are unstable. There is a critical value of the forcing beyond all static $\pi$-kinks fail to exist. Up to this value, the (in)stability of the $\pi$-kinks can be established analytically in the strong coupling limits. Applying a forcing above the critical value causes the nucleation of $2\pi$-kinks and -antikinks. Besides a $\pi$-kink, the unforced system also admits a static $3\pi$-kink. This state is unstable in the continuum models. By combining analytical and numerical methods in the discrete model, it is shown that the stable $\pi$-kink remains stable, and that the unstable $\pi$-kinks cannot be stabilized by decreasing the coupling. The $3\pi$-kink does become stable in the discrete model when the coupling is sufficiently weak., Comment: figures are not included due to file size limit more...

While fundamental-mode discrete solitons have been demonstrated with both self-focusing and defocusing nonlinearity, high-order-mode localized states in waveguide lattices have been studied thus far only for the self-focusing case. In this paper, the existence and stability regimes of dipole, quadrupole and vortex soliton structures in two-dimensional lattices induced with a defocusing nonlinearity are examined by the theoretical and numerical analysis of a generic envelope nonlinear lattice model. In particular, we find that the stability of such high-order-mode solitons is quite different from that with self-focusing nonlinearity. As a simple example, a dipole (``twisted'') mode soliton which may be stable in the focusing case becomes unstable in the defocusing regime. Our results may be relevant to other two-dimensional defocusing periodic nonlinear systems such as Bose-Einstein condensates with a positive scattering length trapped in optical lattices., Comment: 14 pages, 10 figures more...

In the present work, we study dark solitons in dynamical lattices with the saturable nonlinearity and compare them with those in lattices with the cubic nonlinearity. This comparison has become especially relevant in light of recent experimental developments in the former context. The stability properties of the fundamental waves, for both on-site and inter-site modes, are examined analytically and corroborated by numerical results. Furthermore, their dynamical evolution when they are found to be unstable is obtained through appropriately crafted numerical experiments., Comment: 15 pages, 5 figures more...

We investigate fundamental localized modes in 2D lattices with an edge (surface). Interaction with the edge expands the stability area for ordinary solitons, and induces a difference between perpendicular and parallel dipoles; on the contrary, lattice vortices cannot exist too close to the border. Furthermore, we show analytically and numerically that the edge stabilizes a novel wave species, which is entirely unstable in the uniform lattice, namely, a "horseshoe" soliton, consisting of 3 sites. Unstable horseshoes transform themselves into a pair of ordinary solitons., Comment: 6 pages, 4 composite figures more...

We study the mobility of solitons in second-harmonic-generating lattices. Contrary to what is known for their cubic counterparts, discrete quadratic solitons are mobile not only in the one-dimensional (1D) setting, but also in two dimensions (2D). We identify parametric regions where an initial kick applied to a soliton leads to three possible outcomes, namely, staying put, persistent motion, or destruction. For the 2D lattice, it is found that, for the solitary waves, the direction along which they can sustain the largest kick and can attain the largest speed is the diagonal. Basic dynamical properties of the discrete solitons are also discussed in the context of an analytical approximation, in terms of an effective Peierls-Nabarro potential in the lattice setting., Comment: 4 pages more...

In this brief report, we consider parametrically driven bright solitons in the vicinity of the anti-continuum limit. We illustrate the mechanism through which these solitons become unstable due to the collision of the phase mode with the continuous spectrum, or eigenvelues bifurcating thereof. We show how this mechanism typically leads to complete destruction of the bright solitary wave., Comment: 4 pages, 4 figures more...

Dynamics of a fluxon in a stack of coupled long Josephson junctions is studied numerically. Based on the numerical simulations, we show that the dependence of the propagation velocity $c$ on the external bias current $\gamma$ is determined by the ratio of the critical currents of the two junctions $J$., Comment: Tex-file with 8 pages and 13 figures more...