Source author record

Panayotis G. Kevrekidis

Panayotis G. Kevrekidis appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

37works
15topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

37 published item(s)

preprint2026arXiv

Rogue waves in extended Gross-Pitaevskii Models with a Lee-Huang-Yang correction

We explore the existence and dynamical generation of rogue waves (RWs) within a one dimensional quantum droplet bearing environment. RWs are computed by deploying a spacetime fixed point scheme to the relevant extended Gross Pitaevskii equation (eGPE). Parametric regions where the ensuing RWs are different from their counterparts in the nonlinear Schroedinger equation are identified. To corroborate the controllable generation, relevant to ultracold atom experiments, of these rogue patterns we exploit two different protocols. The first is based on interfering dam break flows emanating from Riemann initial conditions and the second refers to the gradient catastrophe of a spatially localized waveform. A multitude of possible RWs are found in this system, spanning waveforms reminiscent of the Peregrine soliton, its spatially periodic variants, namely, the Akhmediev breathers, and other higher order RW solutions of the nonlinear Schroedinger equation. Key elements of the shape of the corresponding eGPE RWs traced back to nonintegrability and the presence of competing interactions are discussed. Our results set the stage for probing a multitude of unexplored rogue like waveforms in such mixtures with competing interactions and should be accessible to current ultracold atom experiments.

preprint2022arXiv

A Modified PINN Approach for Identifiable Compartmental Models in Epidemiology with Applications to COVID-19

A variety of approaches using compartmental models have been used to study the COVID-19 pandemic and the usage of machine learning methods with these models has had particularly notable success. We present here an approach toward analyzing accessible data on Covid-19's U.S. development using a variation of the "Physics Informed Neural Networks" (PINN) which is capable of using the knowledge of the model to aid learning. We illustrate the challenges of using the standard PINN approach, then how with appropriate and novel modifications to the loss function the network can perform well even in our case of incomplete information. Aspects of identifiability of the model parameters are also assessed, as well as methods of denoising available data using a wavelet transform. Finally, we discuss the capability of the neural network methodology to work with models of varying parameter values, as well as a concrete application in estimating how effectively cases are being tested for in a population, providing a ranking of U.S. states by means of their respective testing.

preprint2021arXiv

Discrete embedded solitary waves and breathers in one-dimensional nonlinear lattices

For a one-dimensional linear lattice, earlier work has shown how to systematically construct a slowly-decaying linear potential bearing a localized eigenmode embedded in the continuous spectrum. Here, we extend this idea in two directions: The first one is in the realm of the discrete nonlinear Schrodinger equation, where the linear operator of the Schrodinger type is considered in the presence of a Kerr focusing or defocusing nonlinearity and the embedded linear mode is continued into the nonlinear regime as a discrete solitary wave. The second case is the Klein-Gordon setting, where the presence of a cubic nonlinearity leads to the emergence of embedded-in-the-continuum discrete breathers. In both settings, it is seen that the stability of the modes near the linear limit turns into instability as nonlinearity is increased past a critical value, leading to a dynamical delocalization of the solitary wave (or breathing) state. Finally, we suggest a concrete experiment to observe these embedded modes using a bi-inductive electrical lattice.

preprint2021arXiv

Instabilities of vortex-ring-bright soliton in trapped binary 3D Bose-Einstein condensates

Instabilities of vortex-ring-bright coherent structures in harmonically trapped two-component three-dimensional Bose-Einstein condensates are studied numerically within the coupled Gross-Pitaevskii equations and interpreted analytically. Interestingly, the filled vortex core with a sufficiently large amount of the bright component is observed to reduce the parametric interval of stability of the vortex ring. We have identified the mechanisms of several linear instabilities and one nonlinear parametric instability in this connection. Two of the linear instabilities are qualitatively different from ones reported earlier, to our knowledge, and are associated with azimuthal modes of $m=0$ and $m=1$, i.e., deviations of the vortex from the stationary ring shape. Our nonlinear parametric resonance instability occurs between the $m=0$ and $m=2$ modes and signals the exchange of energy between them.

preprint2021arXiv

Neural Networks Enforcing Physical Symmetries in Nonlinear Dynamical Lattices: The Case Example of the Ablowitz-Ladik Model

In this work we introduce symmetry-preserving, physics-informed neural networks (S-PINNs) motivated by symmetries that are ubiquitous to solutions of nonlinear dynamical lattices. Although the use of PINNs have recently attracted much attention in data-driven discovery of solutions chiefly to partial differential equations, we demonstrate that they fail at enforcing important physical laws including symmetries of solutions and conservation laws. Through the correlation of parity symmetries in both space and time of solutions to differential equations with their group equivariant representation, we construct group-equivariant NNs which respect spatio-temporal parity symmetry. Moreover, we adapt the proposed architecture to enforce different types of periodicity (or localization) of solutions to nonlinear dynamical lattices. We do so by applying S-PINNs to the completely integrable Ablowitz-Ladik model, and performing numerical experiments with a special focus on waveforms that are related to rogue structures. These include the Kuznetsov-Ma soliton, and Akhmediev breather as well as the Peregrine soliton. Our numerical results demonstrate the superiority and robustness of the proposed architecture over standard PINNs.

preprint2021arXiv

Rogue and solitary waves in coupled phononic crystals

In this work we present an analytical and numerical study of rogue and solitary waves in a coupled one-dimensional nonlinear lattice that involves both axial and rotational degrees of freedom. Using a multiple-scale analysis we derive a system of coupled nonlinear Schrödinger type equations in order to approximate solitary waves and rogue waves of the coupled lattice model. Numerical simulations are found to agree with the analytical approximations. We also consider generic initialization data in the form of a Gaussian profile and observe that they can result in the spontaneous formation of rogue wave-like patterns in the lattice. The solitary and rogue waves in the lattice demonstrate both energy isolation and exchange between the axial and rotational degrees of freedom of the system. This suggests that the studied coupled lattice has the potential to be an efficient energy isolation, transfer, and focusing medium.

preprint2020arXiv

Collision of $ϕ^4$ kinks free of the Peierls-Nabarro barrier in the regime of strong discreteness

The two major effects observed in collisions of the continuum $ϕ^4$ kinks are (i) the existence of critical collision velocity above which the kinks always emerge from the collision and (ii) the existence of the escape windows for multi-bounce collisions with the velocity below the critical one, associated with the energy exchange between the kink's internal and translational modes. The potential merger (for sufficiently low collision speeds) of the kink and antikink produces a bion with oscillation frequency $ω_B$, which constantly radiates energy, since its higher harmonics are always within the phonon spectrum. Similar effects have been observed in the discrete $ϕ^4$ kink-antikink collisions for relatively weak discreteness. Here we analyze kinks colliding with their mirror image antikinks in the regime of strong discreteness considering an exceptional discretization of the $ϕ^4$ field equation where the static Peierls-Nabarro potential is precisely zero and the not-too-fast kinks can propagate practically radiating no energy. Several new effects are observed in this case, originating from the fact that the phonon band width is small for strongly discrete lattices and for even higher discreteness an inversion of the phonon spectrum takes place with the short waves becoming low-frequency waves. When the phonon band is narrow, not a bion but a discrete breather with frequency $ω_{DB}$ and all higher harmonics outside the phonon band is formed. When the phonon spectrum is inverted, the kink and antikink become mutually repulsive solitary waves with oscillatory tails, and their collision is possible only for velocities above a threshold value sufficient to overcome their repulsion.

preprint2020arXiv

Decay of two-dimensional quantum turbulence in binary Bose-Einstein condensates

We study two-dimensional quantum turbulence in miscible binary Bose-Einstein condensates in either a harmonic trap or a steep-wall trap through the numerical simulations of the Gross-Pitaevskii equations. The turbulence is generated through a Gaussian stirring potential. When the condensates have unequal intra-component coupling strengths or asymmetric trap frequencies, the turbulent condensates undergo a dramatic decay dynamics to an interlaced array of vortex-antidark structures, a quasi-equilibrium state, of like-signed vortices with an extended size of the vortex core. The time of formation of this state is shortened when the parameter asymmetry of the intra-component couplings or the trap frequencies are enhanced. The corresponding spectrum of the incompressible kinetic energy exhibits two noteworthy features: (i) a $k^{-3}$ power-law around the range of the wave number determined by the spin healing length (the size of the extended vortex-core) and (ii) a flat region around the range of the wave number determined by the density healing length. The latter is associated with the small scale phase fluctuation relegated outside the Thomas-Fermi radius and is more prominent as the strength of intercomponent interaction approaches the strength of intra-component interaction. We also study the impact of the inter-component interaction to the cluster formation of like-signed vortices in an elliptical steep-wall trap, finding that the inter-component coupling gives rise to the decay of the clustered configuration.

preprint2020arXiv

Exploring Critical Points of Energy Landscapes: From Low-Dimensional Examples to Phase Field Crystal PDEs

In the present work we explore the application of a few root-finding methods to a series of prototypical examples. The methods we consider include: (a) the so-called continuous-time Nesterov (CTN) flow method; (b) a variant thereof referred to as the squared-operator method (SOM); and (c) the the joint action of each of the above two methods with the so-called deflation method. More traditional methods such as Newton's method (and its variant with deflation) are also brought to bear. Our toy examples start with a naive one degree-of-freedom (dof) system to provide the lay of the land. Subsequently, we turn to a 2-dof system that is motivated by the reduction of an infinite-dimensional, phase field crystal (PFC) model of soft matter crystallisation. Once the landscape of the 2-dof system has been elucidated, we turn to the full PDE model and illustrate how the insights of the low-dimensional examples lead to novel solutions at the PDE level that are of relevance and interest to the full framework of soft matter crystallization.

preprint2020arXiv

Solitary and Periodic Waves in Collisionless Plasmas: The Adlam-Allen Model Revisited

We consider the Adlam-Allen (AA) system of partial differential equations which, arguably, is the first model that was introduced to describe solitary waves in the context of propagation of hydrodynamic disturbances in collisionless plasmas. Here, we identify the solitary waves of the model by implementing a dynamical systems approach. The latter suggests that the model also possesses periodic wave solutions --which reduce to the solitary wave in the limiting case of infinite period-- as well as rational solutions which are obtained herein. In addition, employing a long-wave approximation via a relevant multiscale expansion method, we establish the asymptotic reduction of the AA system to the Korteweg-de Vries equation. Such a reduction, is not only another justification for the above solitary wave dynamics, but also may offer additional insights for the emergence of other possible plasma waves. Direct numerical simulations are performed for the study of multiple solitary waves and their pairwise interactions. The stability of solitary waves is discussed in terms of potentially relevant criteria, while the robustness of spatially periodic wave solutions is touched upon by our numerical experiments.

preprint2020arXiv

Stability of topological edge states under strong nonlinear effects

We examine the role of strong nonlinearity on the topologically-robust edge state in a one-dimensional system. We consider a chain inspired from the Su-Schrieffer-Heeger model, but with a finite-frequency edge state and the dynamics governed by second-order differential equations. We introduce a cubic onsite-nonlinearity and study this nonlinear effect on the edge state's frequency and linear stability. Nonlinear continuation reveals that the edge state loses its typical shape enforced by the chiral symmetry and becomes generally unstable due to various types of instabilities that we analyze using a combination of spectral stability and Krein signature analysis. This results in an initially-excited nonlinear-edge state shedding its energy into the bulk over a long time. However, the stability trends differ both qualitatively and quantitatively when softening and stiffening types of nonlinearity are considered. In the latter, we find a frequency regime where nonlinear edge states can be linearly stable. This enables high-amplitude edge states to remain spatially localized without shedding their energy, a feature that we have confirmed via long-time dynamical simulations. Finally, we examine the robustness of frequency and stability of nonlinear edge states against disorder, and find that those are more robust under a chiral disorder compared to a non-chiral disorder. Moreover, the frequency-regime where high-amplitude edge states were found to be linearly stable remains intact in the presence of small amount of disorder of both types.

preprint2016arXiv

Energy criterion for the spectral stability of discrete breathers

Discrete breathers are ubiquitous structures in nonlinear anharmonic models ranging from the prototypical example of the Fermi-Pasta-Ulam model to Klein-Gordon nonlinear lattices, among many others. We propose a general criterion for the emergence of instabilities of discrete breathers analogous to the well-established Vakhitov-Kolokolov criterion for solitary waves. The criterion involves the change of monotonicity of the discrete breather's energy as a function of the breather frequency. Our analysis suggests and numerical results corroborate that breathers with increasing (decreasing) energy-frequency dependence are generically unstable in soft (hard) nonlinear potentials.

preprint2016arXiv

Experimental Study of Nonlinear Resonances and Anti-resonances in a Forced, Ordered Granular Chain

We experimentally study a one-dimensional uncompressed granular chain composed of a finite number of identical spherical beads with Hertzian interactions. The chain is harmonically excited by an amplitude- and frequency-dependent boundary drive at its left end and has a fixed boundary at its right end. Such ordered granular media represent an interesting new class of nonlinear acoustic metamaterials, since they exhibit essentially nonlinear acoustics and have been designated as 'sonic vacua' due to the fact that their corresponding speed of sound (as defined in classical acoustics) is zero.This paves the way for essentially nonlinear and energy-dependent acoustics with no counterparts in linear theory. We experimentally detect time-periodic, strongly nonlinear resonances whereby the particles (beads) of the granular chain respond at integer multiples of the excitation period, and which correspond to local peaks of the maximum transmitted force at the chain's right, fixed end.In between these resonances we detect a local minimum of the maximum transmitted forces corresponding to an anti-resonance in the stationary-state dynamics. The experimental results of this work confirm previous theoretical predictions, and verify the existence of strongly nonlinear resonance responses in a system with a complete absence of any linear spectrum; as such, the experimentally detected nonlinear resonance spectrum is passively tunable with energy and sensitive to dissipative effects such as internal structural damping in the beads, and friction or plasticity effects. The experimental results are verified by direct numerical simulations and by numerical stability analysis.

preprint2016arXiv

Performing Hong-Ou-Mandel-type Numerical Experiments with Repulsive Condensates: The case of Dark and Dark-bright Solitons

The Hong-Ou-Mandel experiment leads indistinguishable photons simultaneously reaching a 50:50 beam splitter to emerge on the same port through two-photon interference. Motivated by this phenomenon, we consider numerical experiments of the same flavor for classical, {\it wave} objects in the setting of repulsive condensates. We examine dark solitons interacting with a repulsive barrier, a case in which we find no significant asymmetries in the emerging waves after the collision, presumably due to their topological nature. We also consider case examples of two-component systems, where the dark solitons trap a bright structure in the second-component (dark-bright solitary waves). For these, pronounced asymmetries upon collision {\it are} possible for the non-topological bright component. We also show an example of a similar phenomenology for ring dark-bright structures in two dimensions.

preprint2015arXiv

Dynamics of vortex dipoles in anisotropic Bose-Einstein condensates

We study the motion of a vortex dipole in a Bose-Einstein condensate confined to an anisotropic trap. We focus on a system of ordinary differential equations describing the vortices' motion, which is in turn a reduced model of the Gross-Pitaevskii equation describing the condensate's motion. Using a sequence of canonical changes of variables, we reduce the dimension and simplify the equations of motion. We uncover two interesting regimes. Near a family of periodic orbits known as guiding centers, we find that the dynamics is essentially that of a pendulum coupled to a linear oscillator, leading to stochastic reversals in the overall direction of rotation of the dipole. Near the separatrix orbit in the isotropic system, we find other families of periodic, quasi-periodic, and chaotic trajectories. In a neighborhood of the guiding center orbits, we derive an explicit iterated map that simplifies the problem further. Numerical calculations are used to illustrate the phenomena discovered through the analysis. Using the results from the reduced system we are able to construct complex periodic orbits in the original, partial differential equation, mean-field model for Bose-Einstein condensates, which corroborates the phenomenology observed in the reduced dynamical equations.

preprint2015arXiv

Existence and stability of PT-symmetric vortices in nonlinear two-dimensional square lattices

Vortices symmetric with respect to simultaneous parity and time reversing transformations are considered on the square lattice in the framework of the discrete nonlinear Schrödinger equation. The existence and stability of vortex configurations is analyzed in the limit of weak coupling between the lattice sites, when predictions on the elementary cell of a square lattice (i.e., a single square) can be extended to a large (yet finite) array of lattice cells. Our analytical predictions are found to be in good agreement with numerical computations.

preprint2015arXiv

Formation of rarefaction waves in origami-based metamaterials

We investigate the nonlinear wave dynamics of origami-based metamaterials composed of Tachi-Miura polyhedron (TMP) unit cells. These cells exhibit strain softening behavior under compression, which can be tuned by modifying their geometrical configurations or initial folded conditions. We assemble these TMP cells into a cluster of origami-based metamaterials, and we theoretically model and numerically analyze their wave transmission mechanism under external impact. Numerical simulations show that origami-based metamaterials can provide a prototypical platform for the formation of nonlinear coherent structures in the form of rarefaction waves, which feature a tensile wavefront upon the application of compression to the system. We also demonstrate the existence of numerically exact traveling rarefaction waves. Origami-based metamaterials can be highly useful for mitigating shock waves, potentially enabling a wide variety of engineering applications.

preprint2015arXiv

High energy density in multi-soliton collisions

Solitons are very effective in transporting energy over great distances and collisions between them can produce high energy density spots of relevance to phase transformations, energy localization and defect formation among others. It is then important to study how energy density accumulation scales in multi-soliton collisions. In this study, we demonstrate that the maximal energy density that can be achieved in collision of $N$ slowly moving kinks and antikinks in the integrable sine-Gordon field, remarkably, is proportional to $N^2$, while the total energy of the system is proportional to $N$. This maximal energy density can be achieved only if the difference between the number of colliding kinks and antikinks is minimal, i.e., is equal to 0 for even $N$ and 1 for odd $N$ and if the pattern involves an alternating array of kinks and anti-kinks. Interestingly, for odd (even) $N$ the maximal energy density appears in the form of potential (kinetic) energy, while kinetic (potential) energy is equal to zero. The results of the present study rely on the analysis of the exact multi-soliton solutions for $N=1,2,$ and 3 and on the numerical simulation results for $N=4,5,6,$ and 7. Based on these results one can speculate that the soliton collisions in the sine-Gordon field can, in principle, controllably produce very high energy density. This can have important consequences for many physical phenomena described by the Klein-Gordon equations.

preprint2015arXiv

Nonlinear instabilities of multi-site breathers in Klein-Gordon lattices

In the present work, we explore the possibility of excited breather states in a nonlinear Klein--Gordon lattice to become nonlinearly unstable, even if they are found to be spectrally stable. The mechanism for this fundamentally nonlinear instability is through the resonance with the wave continuum of a multiple of an internal mode eigenfrequency in the linearization of excited breather states. For the nonlinear instability, the internal mode must have its Krein signature opposite to that of the wave continuum. This mechanism is not only theoretically proposed, but also numerically corroborated through two concrete examples of the Klein--Gordon lattice with a soft (Morse) and a hard ($ϕ^4$) potential. Compared to the case of the nonlinear Schr{ö}dinger lattice, the Krein signature of the internal mode relative to that of the wave continuum may change depending on the period of the excited breather state. For the periods for which the Krein signatures of the internal mode and the wave continuum coincide, excited breather states are observed to be nonlinearly stable.

preprint2015arXiv

Solitary waves of a PT-symmetric Nonlinear Dirac equation

In the present work, we consider a prototypical example of a PT-symmetric Dirac model. We discuss the underlying linear limit of the model and identify the threshold of the PT-phase transition in an analytical form. We then focus on the examination of the nonlinear model. We consider the continuation in the PT-symmetric model of the solutions of the corresponding Hamiltonian model and find that the solutions can be continued robustly as stable ones all the way up to the PT-transition threshold. In the latter, they degenerate into linear waves. We also examine the dynamics of the model. Given the stability of the waveforms in the PT-exact phase we consider them as initial conditions for parameters outside of that phase. We find that both oscillatory dynamics and exponential growth may arise, depending on the size of the corresponding "quench". The former can be characterized by an interesting form of bi-frequency solutions that have been predicted on the basis of the SU(1,1) symmetry. Finally, we explore some special, analytically tractable, but not PT-symmetric solutions in the massless limit of the model.

preprint2015arXiv

Staggered $\mathcal{PT}$-symmetric ladders with cubic nonlinearity

We introduce a ladder-shaped chain with each rung carrying a $\mathcal{PT}$ -symmetric gain-loss dipole. The polarity of the dipoles is staggered along the chain, meaning that a rung bearing gain-loss is followed by one bearing loss-gain. This renders the system $\mathcal{PT}$-symmetric in both horizontal and vertical directions. The system is governed by a pair of linearly coupled discrete nonlinear Schrödinger (DNLS) equations with self-focusing or defocusing cubic onsite nonlinearity. Starting from the analytically tractable anti-continuum limit of uncoupled rungs and using the Newton's method for identifying solutions and parametric continuation in the inter-rung coupling for following the associated branches, we construct families of $\mathcal{PT}$-symmetric discrete solitons and identify their stability regions. Waveforms stemming from a single excited rung, as well as ones from multiple rungs are identified. Dynamics of unstable solitons is presented too.

preprint2015arXiv

Traveling waves for the Mass in Mass model of granular chains

In the present work, we consider the mass in mass (or mass with mass) system of granular chains, namely a granular chain involving additionally an internal resonator. For these chains, we rigorously establish that under suitable "anti-resonance" conditions connecting the mass of the resonator and the speed of the wave, bell-shaped traveling wave solutions continue to exist in the system, in a way reminiscent of the results proven for the standard granular chain of elastic Hertzian contacts. We also numerically touch upon settings where the conditions do not hold, illustrating, in line also with recent experimental work, that non-monotonic waves bearing non-vanishing tails may exist in the latter case.

preprint2015arXiv

Weakly Nonlinear Analysis of Vortex Formation in a Dissipative Variant of the Gross-Pitaevskii Equation

For a dissipative variant of the two-dimensional Gross-Pitaevskii equation with a parabolic trap under rotation, we study a symmetry breaking process that leads to the formation of vortices. The first symmetry breaking leads to the formation of many small vortices distributed uniformly near the Thomas-Fermi radius. The instability occurs as a result of a linear instability of a vortex-free steady state as the rotation is increased above a critical threshold. We focus on the second subsequent symmetry breaking, which occurs in the weakly nonlinear regime. At slightly above threshold, we derive a one-dimensional amplitude equation that describes the slow evolution of the envelope of the initial instability. We show that the mechanism responsible for initiating vortex formation is a modulational instability of the amplitude equation. We also illustrate the role of dissipation in the symmetry breaking process. All analyses are confirmed by detailed numerical computations.

preprint2014arXiv

$\mathcal{PT}$-symmetric ladders with a scattering core

We consider a $\mathcal{PT}$-symmetric chain (ladder-shaped) system governed by the discrete nonlinear Schrödinger equation where the cubic nonlinearity is carried solely by two central "rungs" of the ladder. Two branches of scattering solutions for incident plane waves are found. We systematically construct these solutions, analyze their stability, and discuss non-reciprocity of the transmission associated with them. To relate these results to finite-size wavepacket dynamics, we also perform direct simulations of the evolution of the wavepackets, which confirm that the transmission is indeed asymmetric in this nonlinear system with the mutually balanced gain and loss.

preprint2014arXiv

A Mean-Field Analogue of the Hong-Ou-Mandel Experiment With Bright Solitons

In the present work, we theoretically propose and numerically illustrate a mean-field analogue of the Hong-Ou-Mandel experiment with bright solitons. More specifically, we scatter two solitons off of each other (in our setup, the bright solitons play the role of a classical analogue to the quantum photons of the original experiment), while the role of the beam splitter is played by a repulsive Gaussian barrier. In our classical scenario, distinguishability of the particles yields, as expected, a $0.5$ split mass on either side. Nevertheless, for very slight deviations from the completely symmetric scenario a near-perfect transmission i.e., a $|2,0>e$ or a $|0,2 >$ state can be constructed instead, very similarly to the quantum mechanical output. We demonstrate this as a generic feature under slight variations of the relative soliton speed, or of the relative amplitude in a wide parametric regime. We also explore how variations of the properties of the "beam splitter" (i.e., the Gaussian barrier) affect this phenomenology.

preprint2014arXiv

Asymmetric Wave Propagation Through Saturable Nonlinear Oligomers

In the present paper we consider nonlinear dimers and trimers (more generally, oligomers) embedded within a linear Schr{ö}dinger lattice where the nonlinear sites are of saturable type. We examine the stationary states of such chains in the form of plane waves, and analytically compute their reflection and transmission coefficients through the nonlinear oligomer, as well as the corresponding rectification factors which clearly illustrate the asymmetry between left and right propagation in such systems. We examine not only the existence but also the dynamical stability of the plane wave states. Lastly, we generalize our numerical considerations to the more physically relevant case of Gaussian initial wavepackets and confirm that the asymmetry in the transmission properties also persists in the case of such wavepackets.

preprint2014arXiv

Dynamics of Three Non-co-rotating Vortices in Bose-Einstein Condensates

In this work we use standard Hamiltonian-system techniques in order to study the dynamics of three vortices with alternating charges in a confined Bose-Einstein condensate. In addition to being motivated by recent experiments, this system offers a natural vehicle for the exploration of the transition of the vortex dynamics from ordered to progressively chaotic behavior. In particular, it possesses two integrals of motion, the {\it energy} (which is expressed through the Hamiltonian $H$) and the {\it angular momentum} $L$ of the system. By using the integral of the angular momentum, we reduce the system to a two degree-of-freedom one with $L$ as a parameter and reveal the topology of the phase space through the method of Poincaré surfaces of section. We categorize the various motions that appear in the different regions of the sections and we study the major bifurcations that occur to the families of periodic motions of the system. Finally, we correspond the orbits on the surfaces of section to the real space motion of the vortices in the plane.

preprint2014arXiv

Self trapping transition for a nonlinear impurity within a linear chain

In the present work we revisit the issue of the self-trapping dynamical transition at a nonlinear impurity embedded in an otherwise linear lattice. For our Schrödinger chain example, we present rigorous arguments that establish necessary conditions and corresponding parametric bounds for the transition between linear decay and nonlinear persistence of a defect mode. The proofs combine a contraction mapping approach applied in the fully dynamical problem in the case of a 3D-lattice, together with variational arguments for the derivation of parametric bounds for the creation of stationary states associated with the expected fate of the self-trapping dynamical transition. The results are relevant for both power law nonlinearities and saturable ones. The analytical results are corroborated by numerical computations.

preprint2014arXiv

Spectral stability analysis for standing waves of a perturbed Klein-Gordon equation

In the present work, we introduce a new $\mathcal{PT}$-symmetric variant of the Klein-Gordon field theoretic problem. We identify the standing wave solutions of the proposed class of equations and analyze their stability. In particular, we obtain an explicit frequency condition, somewhat reminiscent of the classical Vakhitov-Kolokolov criterion, which sharply separates the regimes of spectral stability and instability. Our numerical computations corroborate the relevant theoretical result.

preprint2013arXiv

Nonlinear dynamics in PT-symmetric lattices

We consider nonlinear dynamics in a finite parity-time-symmetric chain of the discrete nonlinear Schr{ö}dinger (dNLS) type. We work in the range of the gain and loss coefficient when the zero equilibrium state is neutrally stable. We prove that the solutions of the dNLS equation do not blow up in a finite time and the trajectories starting with small initial data remain bounded for all times. Nevertheless, for arbitrary values of the gain and loss parameter, there exist trajectories starting with large initial data that grow exponentially fast for larger times with a rate that is rigorously identified. Numerical computations illustrate these analytical results for dimers and quadrimers.

preprint2013arXiv

Nonlinear stationary states in PT-symmetric lattices

In the present work we examine both the linear and nonlinear properties of two related PT-symmetric systems of the discrete nonlinear Schrodinger (dNLS) type. First, we examine the parameter range for which the finite PT-dNLS chains have real eigenvalues and PT-symmetric linear eigenstates. We develop a systematic way of analyzing the nonlinear stationary states with the implicit function theorem at an analogue of the anti-continuum limit for the dNLS equation. Secondly, we consider the case when a finite PT-dNLS chain is embedded as a defect in the infinite dNLS lattice. We show that the stability intervals of the infinite PT-dNLS lattice are wider than in the case of a finite PT-dNLS chain. We also prove existence of localized stationary states (discrete solitons) in the analogue of the anti-continuum limit for the dNLS equation. Numerical computations illustrate the existence of nonlinear stationary states, as well as the stability and saddle-center bifurcations of discrete solitons.

preprint2013arXiv

Solitary Matter Waves in Combined Linear and Nonlinear Potentials: Detection, Stability, and Dynamics

We study statically homogeneous Bose-Einstein condensates with spatially inhomogeneous interactions and outline an experimental realization of compensating linear and nonlinear potentials that can yield constant-density solutions. We illustrate how the presence of a step in the nonlinearity coefficient can only be revealed dynamically and consider, in particular, how to reveal it by exploiting the inhomogeneity of the sound speed with a defect-dragging experiment. We conduct computational experiments and observe the spontaneous emergence of dark solitary waves. We use effective-potential theory to perform a detailed analytical investigation of the existence and stability of solitary waves in this setting, and we corroborate these results computationally using a Bogoliubov-de Gennes linear stability analysis. We find that dark solitary waves are unstable for all step widths, whereas bright solitary waves can become stable through a symmetry-breaking bifurcation as one varies the step width. Using phase-plane analysis, we illustrate the scenarios that permit this bifurcation and explore the dynamical outcomes of the interaction between the solitary wave and the step.

preprint2012arXiv

Dark lattice solitons in one-dimensional waveguide arrays with defocusing saturable nonlinearities and alternating couplings

In the present work, we examine "binary" waveguide arrays, where the coupling between adjacent sites alternates between two distinct values $C_1$ and $C_2$ and a saturable nonlinearity is present on each site. Motivated by experimental investigations of this type of system in fabricated LiNbO$_3$ arrays, we proceed to analyze the nonlinear wave excitations arising in the self-defocusing nonlinear regime, examining, in particular, dark solitons and bubbles. We find that such solutions may, in fact, possess a reasonably wide, experimentally relevant parametric interval of stability, while they may also feature both prototypical types of instabilities, namely exponential and oscillatory ones, for the same configuration. The spectral properties of the linearization around the solutions and the presence of two bands and associated mini-gaps are discussed and the nontrivial differences from the standard cubic discrete nonlinear Schr{ö}dinger model are highlighted.

preprint2011arXiv

Dark-bright gap solitons in coupled-mode one-dimensional saturable waveguide arrays

In the present work, we consider the dynamics of dark solitons as one mode of a defocusing photorefractive lattice coupled with bright solitons as a second mode of the lattice. Our investigation is motivated by an experiment which illustrates that such coupled states can exist with both components in the first gap of the linear band spectrum. This finding is further extended by the examination of different possibilities from a theoretical perspective, such as symbiotic ones where the bright component is supported by states of the dark component in the first or second gap, or non-symbiotic ones where the bright soliton is also a first-gap state coupled to a first or second gap state of the dark component. While the obtained states are generally unstable, these instabilities typically bear fairly small growth rates which enable their observation for experimentally relevant propagation distances.

preprint2011arXiv

Dynamics of bright solitons and soliton arrays in the nonlinear Schrodinger equation with a combination of random and harmonic potentials

We report results of systematic simulations of the dynamics of solitons in the framework of the one-dimensional nonlinear Schrödinger equation (NLSE), which includes the harmonic-oscillator (HO) potential and a random potential. The equation models experimentally relevant spatially disordered settings in Bose-Einstein condensates (BECs) and nonlinear optics. First, the generation of soliton arrays from a broad initial quasi-uniform state by the modulational instability (MI) is considered, following a sudden switch of the nonlinearity from repulsive to attractive. Then, we study oscillations of a single soliton in this setting, which models a recently conducted experiment in BEC. Basic characteristics of the MI-generated array, such as the number of solitons and their mobility, are reported as functions of the strength and correlation length of the disorder, and of the total norm. For the single oscillating soliton, its survival rate is found. Main features of these dependences are explained qualitatively.

preprint2009arXiv

On the stability of multibreathers in Klein-Gordon chains

In the present paper, a theorem, which determines the linear stability of multibreathers in Klein-Gordon chains, is proven. Specifically, it is shown that for soft nonlinearities, and positive inter-site coupling, only structures with adjacent sites excited out-of-phase may be stable, while only in-phase ones may be stable for negative coupling. The situation is reversed for hard nonlinearities. This theorem can be applied in $n$-site breathers, where $n$ is any finite number and provides an $\cal{O}(\sqrtε)$ estimation of the characteristic exponents of the solution. To complement the analysis, we perform numerical simulations and establish that the results are in excellent agreement with the theoretical predictions, at least for small values of the coupling constant $ε$.

preprint2004arXiv

Motion of discrete solitons assisted by nonlinearity management

We demonstrate that periodic modulation of the nonlinearity coefficient in the discrete nonlinear Schrödinger (DNLS) equation can strongly facilitate creation of traveling solitons in the lattice. We predict this possibility in an analytical form, and test it in direct simulations. Systematic simulations reveal several generic dynamical regimes, depending on the amplitude and frequency of the time modulation, and on initial thrust which sets the soliton in motion. These regimes include irregular motion, regular motion of a decaying soliton, and regular motion of a stable one. The motion may occur in both the straight and reverse directions, relative to the initial thrust. In the case of stable motion, extremely long simulations in a lattice with periodic boundary conditions demonstrate that the soliton keeps moving as long as we can monitor without any visible loss. Velocities of moving stable solitons are in good agreement with the analytical prediction, which is based on requiring a resonance between the ac drive and motion of the soliton through the periodic potential. All the generic dynamical regimes are mapped in the model's parameter space. Collisions between moving stable solitons are briefly investigated too, with a conclusion that two different outcomes are possible: elastic bounce, or bounce with mass transfer from one soliton to the other. The model can be realized experimentally in a Bose-Einstein condensate trapped in a deep optical lattice.