Source author record

P. G. Kevrekidis

P. 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

248works

26topics

4close collaborators

Actions

Connect this record

Open graph Browse works

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.

preprint2022arXiv

A Spectral Analysis of the Nonlinear Schroedinger Equation in the Co-Exploding Frame

The nonlinear Schroedinger model is a prototypical dispersive wave equation that features finite time blowup, either for supercritical exponents (for fixed dimension) or for supercritical dimensions (for fixed nonlinearity exponent). Upon identifying the self-similar solutions in the so-called "co-exploding frame", a dynamical systems analysis of their stability is natural, yet is complicated by the mixed Hamiltonian-dissipative character of the relevant frame. In the present work, we study the spectral picture of the relevant linearized problem. We examine the point spectrum of 3 eigenvalue pairs associated with translation, $U(1)$ and conformal invariances, as well as the continuous spectrum. We find that two eigenvalues become positive, yet are attributed to symmetries and are thus not associated with instabilities. In addition to a vanishing eigenvalue, 3 more are found to be negative and real, while the continuous spectrum is nearly vertical and on the left-half (spectral) plane. Finally, the subtle effects of the boundaries are also assessed and their role in the observed weak eigenvalue oscillations is clarified.

preprint2022arXiv

An Ising machine based on networks of subharmonic electrical resonators

We explore a case example of networks of classical electronic oscillators evolving towards the solution of complex optimization problems. We show that when driven into subharmonic response, a network of such nonlinear electrical resonators can minimize the Ising Hamiltonian on non-trivial graphs such as antiferromagnetically coupled rewired-M{ö}bius ladders. In this context, the spin-up and spin-down states of the Ising machine are represented by the oscillators' response at the even or odd driving cycles. Our experimental setting of driven nonlinear oscillators coupled via a programmable switch matrix leads to a unique energy minimizer when one such exists, and probes frustration where appropriate. Theoretical modeling of the electronic oscillators and their couplings allows us to accurately reproduce the qualitative features of the experimental results. This suggests the promise of this setup as a prototypical one for exploring the capabilities and limitations of such an unconventional computing platform.

preprint2022arXiv

Backcasting COVID-19: A Physics-Informed Estimate for Early Case Incidence

It is widely accepted that the number of reported cases during the first stages of the COVID-19 pandemic severely underestimates the number of actual cases. We leverage delay embedding theorems of Whitney and Takens and use Gaussian Process regression to estimate the number of cases during the first 2020 wave based on the second wave of the epidemic in several European countries, South Korea, and Brazil. We assume that the second wave was more accurately monitored and hence that it can be trusted. We then construct a manifold diffeomorphic to that of the implied original dynamical system, using fatalities or hospitalizations only. Finally, we restrict the diffeomorphism to the reported cases coordinate of the dynamical system. Our main finding is that in the European countries studied, the actual cases are under-reported by as much as 50\%. On the other hand, in South Korea -- which had an exemplary and proactive mitigation approach -- a far smaller discrepancy between the actual and reported cases is predicted, with an approximately 17\% predicted under-estimation. We believe that our backcasting framework is applicable to other epidemic outbreaks where (due to limited or poor quality data) there is uncertainty around the actual cases.

preprint2022arXiv

Dark solitons in a trapped gas of long-range interacting bosons

We consider the interplay of repulsive short-range and same-sign long-range interactions in the dynamics of dark solitons, as prototypical coherent nonlinear excitations in a trapped 1D Bose gas. First, the form of the ground state is examined, and then both the existence of the solitary waves and their stability properties are explored, and corroborated by direct numerical simulations. We find that single- and multiple-dark-soliton states can exist and are generically robust in the presence of long-range interactions. We analyze the modes of vibration of such excitations and find that their respective frequencies are significantly upshifted as the strength of the long-range interactions is increased. Indeed, we find that a prefactor of the long-range interactions considered comparable to the trap strength may upshift the dark soliton oscillation frequency by {\it an order of magnitude}, in comparison to the well established one of $Ω/\sqrt{2}$ in a trap of frequency $Ω$.

preprint2022arXiv

Dispersive Shock Waves in Lattices: A Dimension Reduction Approach

Dispersive shock waves (DSWs), which connect states of different amplitude via a modulated wave train, form generically in nonlinear dispersive media subjected to abrupt changes in state. The primary tool for the analytical study of DSWs is Whitham's modulation theory. While this framework has been successfully employed in many space-continuous settings to describe DSWs, the Whitham modulation equations are virtually intractable in most spatially discrete systems. In this article, we illustrate the relevance of the reduction of the DSW dynamics to a planar ODE in a broad class of lattice examples. Solutions of this low-dimensional ODE accurately describe the orbits of the DSW in self-similar coordinates and the local averages in a manner consistent with the modulation equations. We use data-driven and quasi-continuum approaches within the context of a discrete system of conservation laws to demonstrate how the underlying low dimensional structure of DSWs can be identified and analyzed. The connection of these results to Whitham modulation theory is also discussed.

preprint2022arXiv

Floquet solitons in square lattices: Existence, Stability and Dynamics

In the present work, we revisit a recently proposed and experimentally realized topological 2D lattice with periodically time-dependent interactions. We identify the fundamental solitons, previously observed in experiments and direct numerical simulations, as exact, exponentially localized, periodic in time solutions. This is done for a variety of phase-shift angles of the central nodes upon a period oscillation of the coupling strength. Subsequently, we perform a systematic Floquet stability analysis of the relevant structures. We analyze both their point and their continuous spectrum and find that the solutions are generically stable, aside from the possible emergence of complex quartets due to the collision of bands of continuous spectrum. The relevant instabilities become weaker as the lattice size gets larger. Finally, we also consider multi-soliton analogues of these Floquet states, inspired by the corresponding discrete nonlinear Schrödinger (DNLS) lattice. When exciting initially multiple sites in phase, we find that the solutions reflect the instability of their DNLS multi-soliton counterparts, while for configurations with multiple excited sites in alternating phases, the Floquet states are spectrally stable, again in analogy to their DNLS counterparts.

preprint2022arXiv

Forced symmetry breaking as a mechanism for rogue bursts in a dissipative nonlinear dynamical lattice

We propose an alternative to the standard mechanisms for the formation of rogue waves in a non-conservative, nonlinear lattice dynamical system. We consider an ODE system that features regular periodic bursting arising from forced symmetry breaking. We then connect such potentially exploding units via a diffusive lattice coupling and investigate the resulting spatio-temporal dynamics for different types of initial conditions (localized or extended). We find that in both cases, particular oscillators undergo extremely fast and large amplitude excursions, resembling a rogue wave burst. Furthermore, the probability distribution of different amplitudes exhibits bimodality, with peaks at both vanishing and very large amplitude. While this phenomenology arises over a range of coupling strengths, for large values thereof the system eventually synchronizes and the above phenomenology is suppressed. We use both distributed (such as a synchronization order parameter) and individual oscillator diagnostics to monitor the dynamics and identify potential precursors to large amplitude excursions. We also examine similar behavior with amplitude-dependent diffusive coupling.

preprint2022arXiv

Mixed dispersion nonlinear Schrödinger equation in higher dimensions: theoretical analysis and numerical computations

In the present work we provide a characterization of the ground states of a higher-dimensional quadratic-quartic model of the nonlinear Schr{ö}dinger class with a combination of a focusing biharmonic operator with either an isotropic or an anisotropic defocusing Laplacian operator (at the linear level) and power-law nonlinearity. Examining principally the prototypical example of dimension $d=2$, we find that instability arises beyond a certain threshold coefficient of the Laplacian between the cubic and quintic cases, while all solutions are stable for powers below the cubic. Above the quintic, and up to a critical nonlinearity exponent $p$, there exists a progressively narrowing range of stable frequencies. Finally, above the critical $p$ all solutions are unstable. The picture is rather similar in the anisotropic case, with the difference that even before the cubic case, the numerical computations suggest an interval of unstable frequencies. Our analysis generalizes the relevant observations for arbitrary combinations of Laplacian prefactor $b$ and nonlinearity power $p$.

preprint2022arXiv

Nonlinear Anti-(Parity-Time) symmetric dimer

In the present work we propose a nonlinear anti-$\mathcal{PT}$-symmetric dimer, that at the linear level has been experimentally created in the realm of electric circuit resonators. We find four families of solutions, the so-called upper and lower branches, both in a symmetric and in an asymmetric (symmetry-broken) form. We unveil analytically and confirm numerically the critical thresholds for the existence of such branches and explore the bifurcations (such as saddle-node ones) that delimit their existence, as well as transcritical ones that lead to their potential exchange of stability. We find that out of the four relevant branches, only one, the upper symmetric branch, corresponds to a spectrally and dynamically robust solution. We subsequently leverage detailed direct numerical computations in order to explore the dynamics of the different states, corroborating our spectral analysis results.

preprint2022arXiv

Stationary multi-kinks in the discrete sine-Gordon equation

We consider the existence and spectral stability of static multi-kink structures in the discrete sine-Gordon equation, as a representative example of the family of discrete Klein-Gordon models. The multi-kinks are constructed using Lin's method from an alternating sequence of well-separated kink and antikink solutions. We then locate the point spectrum associated with these multi-kink solutions by reducing the spectral problem to a matrix equation. For an $m$-structure multi-kink, there will be $m$ eigenvalues in the point spectrum near each eigenvalue of the primary kink, and, as long as the spectrum of the primary kink is imaginary, the spectrum of the multi-kink will be as well. We obtain analytic expressions for the eigenvalues of a multi-kink in terms of the eigenvalues and corresponding eigenfunctions of the primary kink, and these are in very good agreement with numerical results. We also perform numerical time-stepping experiments on perturbations of multi-kinks, and the outcomes of these simulations are interpreted using the spectral results.

preprint2022arXiv

Superfluid vortex multipoles and soliton stripes on a torus

We study the existence, stability, and dynamics of vortex dipole and quadrupole configurations in the nonlinear Schrödinger (NLS) equation on the surface of a torus. For this purpose we use, in addition to the full two-dimensional NLS on the torus, a recently derived [Phys. Rev. A 101, 053606 (2021)] reduced point-vortex particle model which is shown to be in excellent agreement with the full NLS evolution. Horizontal, vertical, and diagonal stationary vortex dipoles are identified and continued along the torus aspect ratio and the chemical potential of the solution. Windows of stability for these solutions are identified. We also investigate stationary vortex quadrupole configurations. After eliminating similar solutions induced by invariances and symmetries, we find a total of 16 distinct configurations ranging from horizontal and vertical aligned quadrupoles, to rectangular and rhomboidal quadrupoles, to trapezoidal and irregular quadrupoles. The stability for the least unstable and, potentially, stable quadrapole solutions is monitored both at the NLS and the reduced model levels. Two quadrupole configurations are found to be stable on small windows of the torus aspect ratio and a handful of quadrupoles are found to be very weakly unstable for relatively large parameter windows. Finally, we briefly study the dark soliton stripes and their connection, through a series of bifurcation cascades, with steady state vortex configurations.

preprint2022arXiv

The role of mobility in the dynamics of the COVID-19 epidemic in Andalusia

Metapopulation models have been a popular tool for the study of epidemic spread over a network of highly populated nodes (cities, provinces, countries) and have been extensively used in the context of the ongoing COVID-19 pandemic. In the present work, we revisit such a model, bearing a particular case example in mind, namely that of the region of Andalusia in Spain during the period of the summer-fall of 2020 (i.e., between the first and second pandemic waves). Our aim is to consider the possibility of incorporation of mobility across the province nodes focusing on mobile-phone time dependent data, but also discussing the comparison for our case example with a gravity model, as well as with the dynamics in the absence of mobility. Our main finding is that mobility is key towards a quantitative understanding of the emergence of the second wave of the pandemic and that the most accurate way to capture it involves dynamic (rather than static) inclusion of time-dependent mobility matrices based on cell-phone data. Alternatives bearing no mobility are unable to capture the trends revealed by the data in the context of the metapopulation model considered herein.

preprint2021arXiv

Dark solitons under higher order dispersion

We show theoretically that dark solitons can exist in the presence of pure quartic dispersion, and also in the presence of both quadratic and quartic dispersive effects, displaying a much greater variety of possible solutions and dynamics than for pure quadratic dispersion. The interplay of the two dispersion orders may lead to oscillatory non-vanishing tails, which enables the possibility of bound, potentially stable, multi-soliton states. Dark soliton-like states which connect to low amplitude oscillations are also shown to be possible. Dynamical evolution results corroborate the stability picture obtained, and possible avenues for dark soliton generation are explored.

preprint2021arXiv

Dynamics of two languages competing on a network: a case study

A language dynamics model on a square lattice, which is an extension of the one popularized by Abrams and Strogatz [1], is analyzed using ODE bifurcation theory. For this model we are interested in the existence and spectral stability of structures such as stripes, which are realized through pulses and/or the concatenation of fronts, and spots, which are a contiguous collection of sites in which one language is dominant. Because the coupling between sites is nonlinear, the boundary between sites containing speaking two different languages is "sharp"; in particular, in a PDE approximation it allows for the existence of compactly supported pulses (compactons). The dynamics are considered as a function of the prestige of a language. In particular, it is seen that as the prestige varies, it allows for a language to spread through the network, or conversely for its demise.

preprint2021arXiv

Existence, Stability and Dynamics of Monopole and Alice Ring Solutions in Anti-Ferromagnetic Spinor Condensates

In this work we study the existence, stability, and dynamics of select topological point and line defects in anti-ferromagnetic, polar phase, $F=1$ $^{23}$Na spinor condensates. Specifically, we leverage fixed-point and numerical continuation techniques in three spatial dimensions to identify solution families of monopole and Alice rings as the chemical potential (number of atoms) and trapping strengths are varied within intervals of realizable experimental parameters. We are able to follow the monopole from the linear limit of small atom number all the way to the Thomas-Fermi regime of large atom number. Additionally, and importantly, our studies reveal the existence of {\em two} Alice ring solution branches, corresponding to, relatively, smaller and larger ring radii, that bifurcate from each other in a saddle-center bifurcation as the chemical potential is varied. We find that the monopole solution is always dynamically unstable in the regimes considered. In contrast, we find that the larger Alice ring is indeed stable close to the bifurcation point until it destabilizes from an oscillatory instability bubble for a larger value of the chemical potential. We also report on the possibility of dramatically reducing, yet not completely eliminating, the instability rates for the smaller Alice ring by varying the trapping strengths. The dynamical evolution of the different unstable waveforms is also probed via direct numerical simulations.

preprint2021arXiv

Kink-antikink stripe interactions in the two-dimensional sine-Gordon equation

The main focus of the present work is to study quasi-one-dimensional kink-antikink stripes embedded in the two-dimensional sine-Gordon equation. Using variational techniques, we reduce the interaction dynamics between a kink and an antikink stripe on their respective time and space dependent widths and locations. The resulting reduced system of coupled equations is found to accurately describe the width and undulation dynamics of a single kink stripe as well as that of interacting ones. As an aside, we also discuss two related topics: the computational identification of the kink center and its numerical implications and alternative perturbative and multiple scales approaches to the transverse direction induced dynamics for a single kink stripe in the two-dimensional realm.

preprint2021arXiv

Measurement and Memory in the Periodically Driven Complex Ginzburg-Landau equation

In the present work we illustrate that classical but nonlinear systems may possess features reminiscent of quantum ones, such as memory, upon suitable external perturbation. As our prototypical example, we use the two-dimensional complex Ginzburg-Landau equation in its vortex glass regime. We impose an external drive as a perturbation mimicking a quantum measurement protocol, with a given "measurement rate" (the rate of repetition of the drive) and "mixing rate" (characterized by the intensity of the drive). Using a variety of measures, we find that the system may or may not retain its coherence, statistically retrieving its original glass state, depending on the strength and periodicity of the perturbing field. The corresponding parametric regimes and the associated energy cascade mechanisms involving the dynamics of vortex waveforms and domain boundaries are discussed.

preprint2021arXiv

Moving discrete breathers in a $β$-FPU lattice revisited

In the present work we revisit the existence, stability and dynamical properties of moving discrete breathers in $β$-FPU lattices. On the existence side, we propose a numerical procedure, based on a continuation along a sequence of velocities, that allows to systematically construct breathers traveling more than one lattice site per period. On the stability side, we explore the stability spectrum of the obtained waveforms via Floquet analysis and connect it to the energy-frequency bifurcation diagrams. We illustrate in this context examples of the energy being a multivalued function of the frequency, showcasing the coexistence of different moving breathers at the same frequency. Finally, we probe the moving breather dynamics and observe how the associated instabilities change their speed, typically slowing them down over longtime simulations.

preprint2021arXiv

Phase Diagram, Stability and Magnetic Properties of Nonlinear Excitations in Spinor Bose-Einstein Condensates

We present the phase diagram, the underlying stability and magnetic properties as well as the dynamics of nonlinear solitary wave excitations arising in the distinct phases of a harmonically confined spinor $F=1$ Bose-Einstein condensate. Particularly, it is found that nonlinear excitations in the form of dark-dark-bright solitons exist in the antiferromagnetic and in the easy-axis phase of a spinor gas, being generally unstable in the former while possessing stability intervals in the latter phase. Dark-bright-bright solitons can be realized in the polar and the easy-plane phases as unstable and stable configurations respectively; the latter phase can also feature stable dark-dark-dark solitons. Importantly, the persistence of these types of states upon transitioning, by means of tuning the quadratic Zeeman coefficient from one phase to the other is unravelled. Additionally, the spin-mixing dynamics of stable and unstable matter waves is analyzed, revealing among others the coherent evolution of magnetic dark-bright, nematic dark-bright-bright and dark-dark-dark solitons. Moreover, for the unstable cases unmagnetized or magnetic droplet-like configurations and spin-waves consisting of regular and magnetic solitons are seen to dynamically emerge remaining thereafter robust while propagating for extremely large evolution times. Interestingly, exposing spinorial solitons to finite temperatures, their anti-damping in trap oscillation is showcased. It is found that the latter is suppressed for stronger bright soliton component "fillings". Our investigations pave the wave for a systematic production and analysis involving spin transfer processes of such waveforms which have been recently realized in ultracold experiments.

preprint2020arXiv

A quantitative framework for exploring exit strategies from the COVID-19 lockdown

Following the highly restrictive measures adopted by many countries for combating the current pandemic, the number of individuals infected by SARS-CoV-2 and the associated number of deaths is steadily decreasing. This fact, together with the impossibility of maintaining the lockdown indefinitely, raises the crucial question of whether it is possible to design an exit strategy based on quantitative analysis. Guided by rigorous mathematical results, we show that this is indeed possible: we present a robust numerical algorithm which can compute the cumulative number of deaths that will occur as a result of increasing the number of contacts by a given multiple, using as input only the most reliable of all data available during the lockdown, namely the cumulative number of deaths.

preprint2020arXiv

Controlled Generation of Dark-Bright Soliton Complexes in Two-Component and Spinor Bose-Einstein Condensates

We report on the controlled creation of multiple soliton complexes of the dark-bright type in one-dimensional two-component, three-component and spinor Bose-Einstein condensates. The formation of these states is based on the so-called matter wave interference of separated condensate fragments. In all three cases a systematic numerical study is carried out upon considering different variations of each systems' parameters both in the absence and in the presence of a harmonic trap. It is found that the judicious selection of the initial separation or the chemical potential of the participating components can be utilized to tailor the number of nucleated states. The latter increases as the former parameters are increased. Similarities and differences of the distinct models considered herein are showcased while the robustness of the emerging states is illustrated via the numerical experiments demonstrating their long time propagation. Importantly, for the spinorial system, we unravel the existence of beating dark soliton arrays that are formed due to the spin-mixing dynamics. These states persist in the presence of a parabolic trap, often relevant for associated experimental realizations.

preprint2020arXiv

Kink-Antikink Interaction Forces and Bound States in a Biharmonic $ϕ^4$ Model

We consider the interaction of solitons in a biharmonic, beam model analogue of the well-studied $ϕ^4$ Klein-Gordon theory. Specifically, we calculate the force between a well separated kink and antikink. Knowing their accelerations as a function of separation, we can determine their motion using a simple ODE. There is good agreement between this asymptotic analysis and numerical computation. Importantly, we find the force has an exponentially-decaying oscillatory behaviour (unlike the monotonically attractive interaction in the Klein-Gordon case). Corresponding to the zeros of the force, we predict the existence of an infinite set of field theory equilibria, i.e., kink-antikink bound states. We confirm the first few of these at the PDE level, and verify their anticipated stability or instability. We also explore the implications of this interaction force in the collision between a kink and an oppositely moving antikink.

preprint2020arXiv

Kuznetsov-Ma breather-like solutions in the Salerno model

The Salerno model is a discrete variant of the celebrated nonlinear Schrödinger (NLS) equation interpolating between the discrete NLS (DNLS) equation and completely integrable Ablowitz-Ladik (AL) model by appropriately tuning the relevant homotopy parameter. Although the AL model possesses an explicit time-periodic solution known as the Kuznetsov-Ma (KM) breather, the existence of time-periodic solutions away from the integrable limit has not been studied as of yet. It is thus the purpose of this work to shed light on the existence and stability of time-periodic solutions of the Salerno model. In particular, we vary the homotopy parameter of the model by employing a pseudo-arclength continuation algorithm where time-periodic solutions are identified via fixed-point iterations. We show that the solutions transform into time-periodic patterns featuring small, yet non-decaying far-field oscillations. Remarkably, our numerical results support the existence of previously unknown time-periodic solutions {\it even} at the integrable case whose stability is explored by using Floquet theory. A continuation of these patterns towards the DNLS limit is also discussed.

preprint2020arXiv

Many-body effects on second-order phase transitions in spinor Bose-Einstein condensates and breathing dynamics

We unravel the correlation effects of the second-order quantum phase transitions emerging on the ground state of a harmonically trapped spin-1 Bose gas, upon varying the involved Zeeman terms, as well as its breathing dynamics triggered by quenching the trapping frequency. It is found that the boundaries of the associated magnetic phases are altered in the presence of interparticle correlations for both ferromagnetic and anti-ferromagnetic spin-spin interactions, an effect which becomes more prominent in the few-body scenario. Most importantly, we unveil a correlation-induced shrinking of the anti-ferromagnetic and broken-axisymmetry phases implying that ground states with bosons polarized in a single spin-component are favored. Turning to the dynamical response of the spinor gas it is shown that its breathing frequency is independent of the system parameters while correlations lead to the formation of filamentary patterns in the one-body density of the participating components. The number of filaments is larger for increasing spin-independent interaction strengths or for smaller particle numbers. Each filament maintains its coherence and exhibits an anti-correlated behavior while distinct filaments show significant losses of coherence and are two-body correlated. Interestingly, we demonstrate that for an initial broken-axisymmetry phase an enhanced spin-flip dynamics takes place which can be tuned either via the linear Zeeman term or the quench amplitude.

preprint2020arXiv

Nonlinear edge modes in a honeycomb electrical lattice near the Dirac points

We examine - both experimentally and numerically - a two-dimensional nonlinear driven electrical lattice with honeycomb structure. Drives are considered over a range of frequencies both outside (below and above) and inside the band of linear modes. We identify a number of discrete breathers both existing in the bulk and also (predominantly) ones arising at the domain boundaries, localized either along the arm-chair or along the zig-zag edges. The types of edge-localized breathers observed and computed emerge in distinct frequency bands near the Dirac-point frequency of the dispersion surface while driving the lattice subharmonically (in a spatially homogeneous manner). These observations/computations can represent a starting point towards the exploration of the interplay of nonlinearity and topology in an experimentally tractable system such as the honeycomb electrical lattice.

preprint2020arXiv

Nonlinearity and Topology

The interplay of nonlinearity and topology results in many novel and emergent properties across a number of physical systems such as chiral magnets, nematic liquid crystals, Bose-Einstein condensates, photonics, high energy physics, etc. It also results in a wide variety of topological defects such as solitons, vortices, skyrmions, merons, hopfions, monopoles to name just a few. Interaction among and collision of these nontrivial defects itself is a topic of great interest. Curvature and underlying geometry also affect the shape, interaction and behavior of these defects. Such properties can be studied using techniques such as, e.g. the Bogomolnyi decomposition. Some applications of this interplay, e.g. in nonreciprocal photonics as well as topological materials such as Dirac and Weyl semimetals, are also elucidated.

preprint2020arXiv

Observation and Analysis of Multiple Dark-Antidark Solitons in Two-Component Bose-Einstein Condensates

We report on the static and dynamical properties of multiple dark-antidark solitons (DADs) in two-component, repulsively interacting Bose-Einstein condensates. Motivated by experimental observations involving multiple DADs, we present a theoretical study which showcases that bound states consisting of dark (antidark) solitons in the first (second) component of the mixture exist for different values of interspecies interactions. It is found that ensembles of few DADs may exist as stable configurations, while for larger DAD arrays, the relevant windows of stability with respect to the interspecies interaction strength become progressively narrower. Moreover, the dynamical formation of states consisting of alternating DADs in the two components of the mixture is monitored. A complex dynamical evolution of these states is observed, leading either to sorted DADs or to beating dark-dark solitons depending on the strength of the interspecies coupling. This study demonstrates clear avenues for future investigations of DAD configurations.

preprint2020arXiv

On-demand generation of dark soliton trains in Bose-Einstein condensates

Matter-wave interference mechanisms in one-dimensional Bose-Einstein condensates that allow for the controlled generation of dark soliton trains upon choosing suitable box-type initial configurations are described. First, the direct scattering problem for the defocusing nonlinear Schrödinger equation with nonzero boundary conditions and general box-type initial configurations is discussed, and expressions for the discrete spectrum corresponding to the dark soliton excitations generated by the dynamics are obtained. It is found that the size of the initial box directly affects the number, size and velocity of the solitons, while the initial phase determines the parity of the solutions. The analytical results are compared to those of numerical simulations of the Gross-Pitaevskii equation, both in the absence and in the presence of a harmonic trap. The numerical results bear out the analytical results with excellent agreement.

preprint2020arXiv

Parametrically excited star-shaped patterns at the interface of binary Bose-Einstein condensates

A Faraday-wave-like parametric instability is investigated via mean-field and Floquet analysis in immiscible binary Bose-Einstein condensates. The condensates form a so-called \textit{ball-shell} structure in a two-dimensional harmonic trap. To trigger the dynamics, the scattering length of the core condensate is periodically modulated in time. We reveal that in the dynamics the interface becomes unstable towards the formation of oscillating patterns. The interface oscillates sub-harmonically exhibiting an $m$-fold rotational symmetry that can be controlled by maneuvering the amplitude and the frequency of the modulation. Using Floquet analysis we are able to predict the generated interfacial tension of the mixture and derive a dispersion relation for the natural frequencies of the emergent patterns. A heteronuclear system composed of $^{87}$Rb-$^{85}$Rb atoms can be used for the experimental realization of the phenomenon, yet our results are independent of the specifics of the employed atomic species {and of the parameter at which the driving is applied.

preprint2020arXiv

Quench induced vortex-bright-soliton formation in binary Bose-Einstein condensates

We unravel the spontaneous generation of vortex-bright-soliton structures in binary Bose-Einstein condensates with a small mass imbalance between the species confined in a two-dimensional harmonic trap where one of the two species has been segmented into two parts by a potential barrier. To trigger the dynamics the potential barrier is suddenly removed and subsequently the segments perform a counterflow dynamics. We consider a relative phase difference of $π$ between the segments, while a singly quantized vortex may be imprinted at the center of the other species. The number of vortex structures developed within the segmented species following the merging of its segments is found to depend on the presence of an initial vortex on the other species. In particular, a $π$ phase difference in the segmented species and a vortex in the other species result in a single vortex-bright-soliton structure. However, when the non-segmented species does not contain a vortex the counterflow dynamics of the segmented species gives rise to a vortex dipole in it accompanied by two bright solitary waves arising in the non-segmented species. Turning to strongly mass imbalanced mixtures, with a heavier segmented species, we find that the same overall dynamics takes place, while the quench-induced nonlinear excitations become more robust. Inspecting the dynamics of the angular momentum we show that it can be transferred from one species to the other, and its transfer rate can be tuned by the strength of the interspecies interactions and the mass of the atomic species.

preprint2020arXiv

Solitary waves in the resonant nonlinear Schrödinger equation: stability and dynamical properties

The stability and dynamical properties of the so-called resonant nonlinear Schrödinger (RNLS) equation, are considered. The RNLS is a variant of the nonlinear Schrödinger (NLS) equation with the addition of a perturbation used to describe wave propagation in cold collisionless plasmas. We first examine the modulational stability of plane waves in the RNLS model, identifying the modifications of the associated conditions from the NLS case. We then move to the study of solitary waves with vanishing and nonzero boundary conditions. Interestingly the RNLS, much like the usual NLS, exhibits both dark and bright soliton solutions depending on the relative signs of dispersion and nonlinearity. The corresponding existence, stability and dynamics of these solutions are studied systematically in this work.

preprint2020arXiv

Stability of traveling waves in a driven Frenkel-Kontorova model

In this work we revisit a classical problem of traveling waves in a damped Frenkel-Kontorova lattice driven by a constant external force. We compute these solutions as fixed points of a nonlinear map and obtain the corresponding kinetic relation between the driving force and the velocity of the wave for different values of the damping coefficient. We show that the kinetic curve can become non-monotone at small velocities, due to resonances with linear modes, and also at large velocities where the kinetic relation becomes multivalued. Exploring the spectral stability of the obtained waveforms, we identify, at the level of numerical accuracy of our computations, a precise criterion for instability of the traveling wave solutions: monotonically decreasing portions of the kinetic curve always bear an unstable eigendirection. We discuss why the validity of this criterion in the {\it dissipative} setting is a rather remarkable feature offering connections to the Hamiltonian variant of the model and of lattice traveling waves more generally. Our stability results are corroborated by direct numerical simulations which also reveal the possible outcomes of dynamical instabilities.

preprint2020arXiv

Vortex Pairs in the Discrete Nonlinear Schrödinger Equation

The aim of the present work is to examine the role of discreteness in the interaction of both co-winding and counter-winding vortices in the context of the nonlinear Schr{ö}dinger equation. Contrary to the well-known rotation of same charge vortices, and translation of opposite charge vortices, we find that strong discreteness is able to halt both types of pairs into stationary, potentially stable configurations up to a critical inter-site coupling strength. Past the relevant critical point the behavior is also somewhat counter-intuitive as, for instance, counterwinding vortices start moving but also approach each other. This lateral motion becomes weaker as the continuum limit is approached and we conjecture that genuine traveling appears only at the continuum limit. Analogous features arise in the cowinding where the discrete coherent structure pair spirals outward, with rigid rotation being restored only in the continuum limit.

preprint2019arXiv

Bifurcation analysis of stationary solutions of two-dimensional coupled Gross-Pitaevskii equations using deflated continuation

Recently, a novel bifurcation technique known as the deflated continuation method (DCM) was applied to the single-component nonlinear Schrödinger (NLS) equation with a parabolic trap in two spatial dimensions. The bifurcation analysis carried out by a subset of the present authors shed light on the configuration space of solutions of this fundamental problem in the physics of ultracold atoms. In the present work, we take this a step further by applying the DCM to two coupled NLS equations in order to elucidate the considerably more complex landscape of solutions of this system. Upon identifying branches of solutions, we construct the relevant bifurcation diagrams and perform spectral stability analysis to identify parametric regimes of stability and instability and to understand the mechanisms by which these branches emerge. The method reveals a remarkable wealth of solutions: these do not only include some of the well-known ones including, e.g., from the Cartesian or polar small amplitude limits of the underlying linear problem but also a significant number of branches that arise through (typically pitchfork) bifurcations. In addition to presenting a ``cartography'' of the landscape of solutions, we comment on the challenging task of identifying {\it all} solutions of such a high-dimensional, nonlinear problem.

preprint2019arXiv

Discrete breathers in $ϕ^4$ and related models

We touch upon the wide topic of discrete breather formation with a special emphasis on the the $ϕ^4$ model. We start by introducing the model and discussing some of the application areas/motivational aspects of exploring time periodic, spatially localized structures, such as the discrete breathers. Our main emphasis is on the existence, and especially on the stability features of such solutions. We explore their spectral stability numerically, as well as in special limits (such as the vicinity of the so-called anti-continuum limit of vanishing coupling) analytically. We also provide and explore a simple, yet powerful stability criterion involving the sign of the derivative of the energy vs. frequency dependence of such solutions. We then turn our attention to nonlinear stability, bringing forth the importance of a topological notion, namely the Krein signature. Furthermore, we briefly touch upon linearly and nonlinearly unstable dynamics of such states. Some special aspects/extensions of such structures are only touched upon, including moving breathers and dissipative variations of the model and some possibilities for future work are highlighted.

preprint2019arXiv

Existence and spectral stability of multi-pulses in discrete Hamiltonian lattice systems

In the present work, we consider the existence and spectral stability of multi-pulse solutions in Hamiltonian lattice systems. We provide a general framework for the study of such wave patterns based on a discrete analogue of Lin's method, previously used in the continuum realm. We develop explicit conditions for the existence of $m$-pulse structures and subsequently develop a reduced matrix allowing us to address their spectral stability. As a prototypical example for the manifestation of the details of the formulation, we consider the discrete nonlinear Schrödinger equation. Different families of $2$- and $3$-pulse solitary waves are discussed, and analytical expressions for the corresponding stability eigenvalues are obtained which are in very good agreement with numerical results.

preprint2019arXiv

Rogue Waves and Periodic Solutions of a Nonlocal Nonlinear Schrödinger Model

In the present work, a nonlocal nonlinear Schrödinger (NLS) model is studied by means of a recent technique that identifies solutions of partial differential equations, by considering them as fixed points in {\it space-time}. This methodology allows to perform a continuation of well-known solutions of the local NLS model to the nonlocal case. Four different examples of this type are presented, namely (a) the rogue wave in the form of the Peregrine soliton, (b) the generalization thereof in the form of the Kuznetsov-Ma breather, as well as two spatio-temporally periodic solutions in the form of elliptic functions. Importantly, all four waveforms can be continued in intervals of the parameter controlling the nonlocality of the model. The first two can be continued in a narrower interval, while the periodic ones can be extended to arbitrary nonlocalities and, in fact, present an intriguing bifurcation whereby they merge with (only) spatially periodic structures. The results suggest the generic relevance of rogue waves and related structures, as well as periodic solutions, in nonlocal NLS models.

preprint2019arXiv

Two-dimensional rogue waves on zero background of the Davey-Stewartson II equation

A prototypical example of a rogue wave structure in a two-dimensional model is presented in the context of the Davey-Stewartson~II (DS~II) equation arising in water waves. The analytical methodology involves a Taylor expansion of an eigenfunctionof the model's Lax pair which is used to form a hierarchy of infinitely many new eigenfunctions. These are used for the construction of two-dimensional (2D) rogue waves (RWs) of the DS~II equation by the even-fold Darboux transformation (DT). The obtained 2D RWs, which are localized in both space and time, can be viewed as a 2D analogue of the Peregrine soliton and are thus natural candidates to describe oceanic RW phenomena,as well as ones in 2D fluid systems and water tanks.

preprint2017arXiv

Existence, Stability & Dynamics of Nonlinear Modes in a 2d Partially $\mathcal{PT}$ Symmetric Potential

It is known that multidimensional complex potentials obeying $\mathcal{PT}$-symmetry may possess all real spectra and continuous families of solitons. Recently it was shown that for multi-dimensional systems these features can persist when the parity symmetry condition is relaxed so that the potential is invariant under reflection in only a single spatial direction. We examine the existence, stability and dynamical properties of localized modes within the cubic nonlinear Schrödinger equation in such a scenario of partially $\mathcal{PT}$-symmetric potential.

preprint2016arXiv

$SO(2)$-induced breathing patterns in multi-component Bose-Einstein condensates

In this work, we employ the $SO(2)$-rotations of a two-component, one-, two- and three-dimensional nonlinear Schrödinger system at and near the Manakov limit, to construct vector solitons and vortex structures. This way, stable stationary dark-bright solitons and their higher-dimensional siblings are transformed into robust oscillatory dark-dark solitons (and generalizations thereof), with and without a harmonic confinement. By analogy to the one-dimensional case, vector higher-dimensional structures take the form of vortex-vortex states in two dimensions and, e.g., vortex ring-vortex ring ones in three dimensions. We consider the effects of unequal (self- and cross-) interaction strengths, where the $SO(2)$ symmetry is only approximately satisfied, showing the dark-dark soliton oscillation is generally robust. Similar features are found in higher dimensions too, although our case examples suggest that phenomena such as phase separation may contribute to the associated dynamics. These results, in connection with the experimental realization of one-dimensional variants of such states in optics and Bose-Einstein condensates (BECs), suggest the potential observation of the higher-dimensional bound states proposed herein.

preprint2016arXiv

A PT-symmetric dual-core system with the sine-Gordon nonlinearity and derivative coupling

We propose a system of sine-Gordon equations, with the $\mathcal{PT}$ symmetry represented by balanced gain and loss in them. The equations are coupled by sine-field terms and first-order derivatives. The sinusoidal coupling stems from local interaction between adjacent particles in the coupled Frenkel-Kontorova (FK) chains and related sine-lattices, while the cross-derivative coupling, which was not considered before, is induced by \emph{three-particle} interactions, provided that the particles in the parallel FK\ chains move in different directions. Nonlinear wave structures are then studied in this model. In particular, kink-kink (KK) and kink-antikink (KA) complexes are explored by means of analytical and numerical methods. It is predicted analytically and confirmed numerically that the complexes are unstable for one sign of the sinusoidal coupling, and stable for another. Stability regions are delineated in the underlying parameter space. Unstable complexes split into free kinks/antikinks that may propagate or become stationary, depending on whether they are subject to gain or loss, respectively.

preprint2016arXiv

Dark spherical shell solitons in three-dimensional Bose-Einstein condensates: Existence, stability and dynamics

In this work we study spherical shell dark soliton states in three-dimensional atomic Bose-Einstein condensates. Their symmetry is exploited in order to analyze their existence, as well as that of topologically charged variants of the structures, and, importantly, to identify their linear stability Bogolyubov-de Gennes spectrum. We compare our effective 1D spherical and 2D cylindrical computations with the full 3D numerics. An important conclusion is that such spherical shell solitons can be stable sufficiently close to the linear limit of the isotropic condensates considered herein. We have also identified their instabilities leading to the emergence of vortex line and vortex ring cages. In addition, we generalize effective particle pictures of lower dimensional dark solitons and ring dark solitons to the spherical shell solitons concerning their equilibrium radius and effective dynamics around it. In this case too, we favorably compare the resulting predictions such as the shell equilibrium radius, qualitatively and quantitatively, with full numerical solutions in 3D.

preprint2016arXiv

On the spontaneous time-reversal symmetry breaking in synchronously-pumped passive Kerr resonators

We study the spontaneous temporal symmetry breaking instability in a coherently-driven passive optical Kerr resonator observed experimentally by Xu and Coen in Opt.~Lett.~{\bf 39}, 3492 (2014). We perform a detailed stability analysis of the Lugiato-Lefever model for the optical Kerr resonators and analyze the temporal bifurcation structure of stationary symmetric and the emerging asymmetric states as a function of the pump power. For intermediate pump powers a pitchfork loop is responsible for the destabilization of symmetric states towards stationary asymmetric ones while at large pump powers we find the emergence of periodic asymmetric solutions via a Hopf bifurcation. From a theoretical perspective, we use local bifurcation theory in order to analyze the most unstable eigenmode of the system. We also explore a non-conservative variational approximation capturing, among others, the evolution of the solution's amplitude, width and center of mass. Both methods provide insight towards the pitchfork bifurcations associated with the symmetry breaking.

preprint2016arXiv

Rogue Waves in Ultracold Bosonic Seas

In this work, we numerically consider the initial value problem for nonlinear Schrödinger (NLS) type models arising in the physics of ultracold boson gases, with generic Gaussian wavepacket initial data. The corresponding Gaussian's width and, wherever relevant also its amplitude, serve as control parameters. First we explore the one-dimensional, standard NLS equation with general power law nonlinearity, in which large amplitude excitations reminiscent of Peregrine solitons or regular solitons appear to form, as the width of the relevant Gaussian is varied. Furthermore, the variation of the nonlinearity exponent aims at a first glimpse of the interplay between rogue or soliton formation and collapse features. The robustness of the main features to noise in the initial data is also confirmed. To better connect our study with the physics of atomic condensates, and explore the role of dimensionality effects, we also consider the nonpolynomial Schrödinger equation (NPSE), as well as the full three-dimensional NLS equation, and examine the degree to which relevant considerations generalize.

preprint2016arXiv

Scattering of Waves by Impurities in Precompressed Granular Chains

We study scattering of waves by impurities in strongly precompressed granular chains. We explore the linear scattering of plane waves and identify a closed-form expression for the reflection and transmission coefficients for the scattering of the waves from both a single impurity and a double impurity. For single-impurity chains, we show that, within the transmission band of the host granular chain, high-frequency waves are strongly attenuated (such that the transmission coefficient vanishes as the wavenumber $k\rightarrow \pmπ$), whereas low-frequency waves are well-transmitted through the impurity. For double-impurity chains, we identify a resonance --- enabling full transmission at a particular frequency --- in a manner that is analogous to the Ramsauer--Townsend (RT) resonance from quantum physics. We also demonstrate that one can tune the frequency of the RT resonance to any value in the pass band of the host chain. We corroborate our theoretical predictions both numerically and experimentally, and we directly observe complete transmission for frequencies close to the RT resonance frequency. Finally, we show how this RT resonance can lead to the existence of reflectionless modes even in granular chains (including disordered ones) with multiple double impurities.

preprint2016arXiv

Stability of solitary waves and vortices in a 2D nonlinear Dirac model

We explore a prototypical two-dimensional model of the nonlinear Dirac type and examine its solitary wave and vortex solutions. In addition to identifying the stationary states, we provide a systematic spectral stability analysis, illustrating the potential of spinor solutions consisting of a soliton in one component and a vortex in the other to be neutrally stable in a wide parametric interval of frequencies. Solutions of higher vorticity are generically unstable and split into lower charge vortices in a way that preserves the total vorticity. These results pave the way for a systematic stability and dynamics analysis of higher dimensional waveforms in a broad class of nonlinear Dirac models and a comparison revealing nontrivial differences with respect to their better understood non-relativistic analogue, the nonlinear Schrödinger equation.

preprint2016arXiv

The dynamical playground of a higher-order cubic Ginzburg-Landau equation: from orbital connections and limit cycles to invariant tori and the onset of chaos

The dynamical behavior of a higher-order cubic Ginzburg-Landau equation is found to include a wide range of scenarios due to the interplay of higher-order physically relevant terms. We find that the competition between the third-order dispersion and stimulated Raman scattering effects, gives rise to rich dynamics: this extends from Poincaré-Bendixson--type scenarios, in the sense that bounded solutions may converge either to distinct equilibria via orbital connections, or space-time periodic solutions, to the emergence of almost periodic and chaotic behavior. One of our main results is that the third-order dispersion has a dominant role in the development of such complex dynamics, since it can be chiefly responsible (i.e., even in the absence of the other higher-order effects) for the existence of the periodic, quasi-periodic and chaotic spatiotemporal structures. Suitable low-dimensional phase space diagnostics are devised and used to illustrate the different possibilities and identify their respective parametric intervals over multiple parameters of the model.

preprint2016arXiv

The emergence and analysis of Kuramoto-Sakaguchi-like models as an effective description for the dynamics of coupled Wien-bridge oscillators

We derive the Kuramoto-Sakaguchi model from the basic circuit equations governing two coupled Wien-bridge oscillators. A Wien-bridge oscillator is a particular realization of a tunable autonomous oscillator that makes use of frequency filtering (via a RC band-pass filter) and positive feedback (via an Op-Amp). In the last few years, such oscillators have started to be utilized in synchronization studies. We first show that the Wien-bridge circuit equations can be cast in the form of a coupled pair of Duffing - Van der Pol equations. Subsequently, by applying the method of multiple time scales, we derive the differential equations that govern the slow evolution of the oscillator phases and amplitudes. These equations are directly reminiscent of the Kuramoto-Sakaguchi type models for the study of synchronization. We analyze the resulting system in terms of existence and stability of various coupled oscillator solutions and explain on that basis how their synchronization emerges. The phase-amplitude equations are also compared numerically to the original circuit equations, and good agreement is found. Finally, we report on experimental measurements on two coupled Wien-bridge oscillators and relate the results back to the theoretical predictions.

preprint2016arXiv

Vector Dark-Antidark Solitary Waves in Multi-Component Bose-Einstein condensates

Multi-component Bose-Einstein condensates exhibit an intriguing variety of nonlinear structures. In recent theoretical work, the notion of magnetic solitons has been introduced. Here we generalize this concept to vector dark-antidark solitary waves in multi-component Bose-Einstein condensates. We first provide concrete experimental evidence for such states in an atomic BEC and subsequently illustrate the broader concept of these states, which are based on the interplay between miscibility and inter-component repulsion. Armed with this more general conceptual framework, we expand the notion of such states to higher dimensions presenting the possibility of both vortex-antidark states and ring-antidark-ring (dark soliton) states. We perform numerical continuation studies, investigate the existence of these states and examine their stability using the method of Bogolyubov-de Gennes analysis. Dark-antidark and vortex-antidark states are found to be stable for broad parametric regimes. In the case of ring dark solitons, where the single-component ring state is known to be unstable, the vector entity appears to bear a progressively more and more stabilizing role as the inter-component coupling is increased.

preprint2016arXiv

Vortex-soliton complexes in coupled nonlinear Schrödinger equations with unequal dispersion coefficients

We consider a two-component, two-dimensional nonlinear Schrödinger system with unequal dispersion coefficients and self-defocusing nonlinearities, chiefly with equal strengths of the self- and cross-interactions. In this setting, a natural waveform with a nonvanishing background in one component is a vortex, which induces an effective potential well in the second component, via the nonlinear coupling of the two components. We show that the potential well may support not only the fundamental bound state, but also multi-ring excited radial state complexes for suitable ranges of values of the dispersion coefficient of the second component. We systematically explore the existence, stability, and nonlinear dynamics of these states. The complexes involving the excited radial states are weakly unstable, with a growth rate depending on the dispersion of the second component. Their evolution leads to transformation of the multi-ring complexes into stable VB solitons ones with the fundamental state in the second component. The excited states may be stabilized by a harmonic-oscillator trapping potential, as well as by unequal strengths of the self- and cross-repulsive nonlinearities.

preprint2015arXiv

Bose-Einstein condensation: Twenty years after

The aim of this introductory article is two-fold. First, we aim to offer a general introduction to the theme of Bose-Einstein condensates, and briefly discuss the evolution of a number of relevant research directions during the last two decades. Second, we introduce and present the articles that appear in this Special Volume of Romanian Reports in Physics celebrating the conclusion of the second decade since the experimental creation of Bose-Einstein condensation in ultracold gases of alkali-metal atoms.

preprint2015arXiv

Bright Discrete Solitons in Spatially Modulated DNLS Systems

In the present work, we revisit the highly active research area of inhomogeneously nonlinear defocusing media and consider the existence, spectral stability and nonlinear dynamics of bright solitary waves in them. We use the anti-continuum limit of vanishing coupling as the starting point of our analysis, enabling in this way a systematic characterization of the branches of solutions. Our stability findings and bifurcation characteristics reveal the enhanced robustness and wider existence intervals of solutions with a broader support, culminating in the "extended" solution in which all sites are excited. Our eigenvalue predictions are corroborated by numerical linear stability analysis. Finally, the dynamics also reveal a tendency of the solution profiles to broaden, in line with the above findings. These results pave the way for further explorations of such states in discrete systems, including in higher dimensional settings.

preprint2015arXiv

Collapse for the higher-order nonlinear Schrödinger equation

We examine conditions for finite-time collapse of the solutions of the higher-order nonlinear Schrödinger (NLS) equation incorporating third-order dispersion, self-steepening, linear and nonlinear gain and loss, and Raman scattering; this is a system that appears in many physical contexts as a more realistic generalization of the integrable NLS. By using energy arguments, it is found that the collapse dynamics is chiefly controlled by the linear/nonlinear gain/loss strengths. We identify a critical value of the linear gain, separating the possible decay of solutions to the trivial zero-state, from collapse. The numerical simulations, performed for a wide class of initial data, are found to be in very good agreement with the analytical results, and reveal long-time stability properties of localized solutions. The role of the higher-order effects to the transient dynamics is also revealed in these simulations.

preprint2015arXiv

Conical Wave Propagation and Diffraction in 2D Hexagonally Packed Granular Lattices

Linear and nonlinear mechanisms for conical wave propagation in two-dimensional lattices are explored in the realm of phononic crystals. As a prototypical example, a statically compressed granular lattice of spherical particles arranged in a hexagonal packing configuration is analyzed. Upon identifying the dispersion relation of the underlying linear problem, the resulting diffraction properties are considered. Analysis both via a heuristic argument for the linear propagation of a wavepacket, as well as via asymptotic analysis leading to the derivation of a Dirac system suggests the occurrence of conical diffraction. This analysis is valid for strong precompression i.e., near the linear regime. For weak precompression, conical wave propagation is still possible, but the resulting expanding circular wave front is of a non-oscillatory nature, resulting from the complex interplay between the discreteness, nonlinearity and geometry of the packing. The transition between these two types of propagation is explored.

preprint2015arXiv

Conservation laws, exact travelling waves and modulation instability for an extended nonlinear Schrödinger equation

We study various properties of solutions of an extended nonlinear Schrödinger (ENLS) equation, which arises in the context of geometric evolution problems -- including vortex filament dynamics -- and governs propagation of short pulses in optical fibers and nonlinear metamaterials. For the periodic initial-boundary value problem, we derive conservation laws satisfied by local in time, weak $H^2$ (distributional) solutions, and establish global existence of such weak solutions. The derivation is obtained by a regularization scheme under a balance condition on the coefficients of the linear and nonlinear terms -- namely, the Hirota limit of the considered ENLS model. Next, we investigate conditions for the existence of traveling wave solutions, focusing on the case of bright and dark solitons. The balance condition on the coefficients is found to be essential for the existence of exact analytical soliton solutions; furthermore, we obtain conditions which define parameter regimes for the existence of traveling solitons for various linear dispersion strengths. Finally, we study the modulational instability of plane waves of the ENLS equation, and identify important differences between the ENLS case and the corresponding NLS counterpart. The analytical results are corroborated by numerical simulations, which reveal notable differences between the bright and the dark soliton propagation dynamics, and are in excellent agreement with the analytical predictions of the modulation instability analysis.

preprint2015arXiv

Dark solitons near potential and nonlinearity steps

We study dark solitons near potential and nonlinearity steps and combinations thereof, forming rectangular barriers. This setting is relevant to the contexts of atomic Bose-Einstein condensates (where such steps can be realized by using proper external fields) and nonlinear optics (for beam propagation near interfaces separating optical media of different refractive indices). We use perturbation theory to develop an equivalent particle theory, describing the matter-wave or optical soliton dynamics as the motion of a particle in an effective potential. This Newtonian dynamical problem provides information for the soliton statics and dynamics, including scenarios of reflection, transmission, or quasi-trapping at such steps. The case of multiple such steps and its connection to barrier potentials is also touched upon. Our analytical predictions are found to be in very good agreement with the corresponding numerical results.

preprint2015arXiv

Dark-bright solitons in coupled nonlinear Schrödinger equations with unequal dispersion coefficients

We study a two-component nonlinear Schr{ö}dinger system with equal, repulsive cubic interactions and different dispersion coefficients in the two components. We consider states that have a dark solitary wave in one component. Treating it as a frozen one, we explore the possibility of the formation of bright-solitonic structures in the other component. We identify bifurcation points at which such states emerge in the bright component in the linear limit and explore their continuation into the nonlinear regime. An additional analytically tractable limit is found to be that of vanishing dispersion of the bright component. We numerically identify regimes of potential stability, not only of the single-peak ground state (the dark-bright soliton), but also of excited states with one or more zero crossings in the bright component. When the states are identified as unstable, direct numerical simulations are used to investigate the outcome of the instability development. Although our principal focus is on the homogeneous setting, we also briefly touch upon the counterintuitive impact of the potential presence of a parabolic trap on the states of interest.

preprint2015arXiv

Discrete Solitons and Vortices in Anisotropic Hexagonal and Honeycomb Lattices

In the present work, we consider the self-focusing discrete nonlinear Schrodinger equation on hexagonal and honeycomb lattice geometries. Our emphasis is on the study of the effects of anisotropy, motivated by the tunability afforded in recent optical and atomic physics experiments. We find that important classes of solutions, such as the so-called discrete vortices, undergo destabilizing bifurcations as the relevant anisotropy control parameter is varied. We quantify these bifurcations by means of explicit analytical calculations of the solutions, as well as of their spectral linearization eigenvalues. Finally, we corroborate the relevant stability picture through direct numerical computations. In the latter, we observe the prototypical manifestation of these instabilities to be the spontaneous rearrangement of the solution, for larger values of the coupling, into localized waveforms typically centered over fewer sites than the original unstable structure. For weak coupling, the instability appears to result in a robust breathing of the relevant waveforms.

preprint2015arXiv

Effects of interactions on the generalized Hong-Ou-Mandel effect

We numerically investigate the influence of interactions on the generalized Hong-Ou-Mandel (HOM) effect for bosonic particles and show results for the cases of $N=2$, $N=3$ and $N=4$ bosons interacting with a beam splitter, whose role is played by a $δ$-barrier. In particular, we focus on the effect of attractive interactions and compare the results with the repulsive case, as well as with the analytically available results for the non-interacting case (that we use as a benchmark). We observe a fermionization effect both for growing repulsive and attractive interactions, i.e., the dip in the HOM coincidence count is progressively smeared out, for increasing interaction strengths. The role of input asymmetries is also explored.

preprint2015arXiv

Exciting and Harvesting Vibrational States in Harmonically Driven Granular Chains

This article explores the excitation of different vibrational states in a spatially extended dynamical system through theory and experiment. As a prototypical example, we consider a one-dimensional packing of spherical particles (a so-called granular chain) that is subject to harmonic boundary excitation. The combination of the multi-modal nature of the system and the strong coupling between the particles due to the nonlinear Hertzian contact force leads to broad regions in frequency where different vibrational states are possible. In certain parametric regions, we demonstrate that the Nonlinear Schrodinger (NLS) equation predicts the corresponding modes fairly well. We propose that nonlinear multi-modal systems can be useful in vibration energy harvest- ing and discuss a prototypical framework for its realization. The electromechanical model we derive predicts accurately the conversion from mechanical to electrical energy observed in the experiments.

preprint2015arXiv

Existence, Stability and Dynamics of Discrete Solitary Waves in a Binary Waveguide Array

Recent work has explored binary waveguide arrays in the long-wavelength, near-continuum limit, here we examine the opposite limit, namely the vicinity of the so-called anti-continuum limit. We provide a systematic discussion of states involving one, two and three excited waveguides, and provide comparisons that illustrate how the stability of these states differ from the monoatomic limit of a single type of waveguide. We do so by developing a general theory which systematically tracks down the key eigenvalues of the linearized system. When we find the states to be unstable, we explore their dynamical evolution through direct numerical simulations. The latter typically illustrate, for the parameter values considered herein, the persistence of localized dynamics and the emergence for the duration of our simulations of robust quasi-periodic states for two excited sites. As the number of excited nodes increase, the unstable dynamics feature less regular oscillations of the solution's amplitude.

preprint2015arXiv

Generating and Manipulating Quantized Vortices On-Demand in a Bose-Einstein Condensate: a Numerical Study

We numerically investigate an experimentally viable method, that we will refer to as the "chopsticks method", for generating and manipulating on-demand several vortices in a highly oblate atomic Bose-Einstein condensate (BEC) in order to initialize complex vortex distributions for studies of vortex dynamics. The method utilizes moving laser beams (the "chopsticks") to generate, capture and transport vortices inside and outside the BEC. We examine in detail this methodology and show a wide parameter range of applicability for the prototypical two-vortex case, and show case examples of producing and manipulating several vortices for which there is no net circulation, equal numbers of positive and negative circulation vortices, and for which there is one net quantum of circulation. We find that the presence of dissipation can help stabilize the pinning of the vortices on their respective laser beam pinning sites. Finally, we illustrate how to utilize laser beams as repositories that hold large numbers of vortices and how to deposit individual vortices in a sequential fashion in the repositories in order to construct superfluid flows about the repository beams with several quanta of circulation.

preprint2015arXiv

Multidimensional Discrete Compactons in Nonlinear Schrödinger Lattices with Strong Nonlinearity Management

The existence of multidimensional lattice compactons in the discrete nonlinear Schrödinger equation in the presence of fast periodic time modulations of the nonlinearity is demonstrated. By averaging over the period of the fast modulations, a new effective averaged dynamical equation arises with coupling constants involving Bessel functions of the first and zeroth kind. These terms allow one to solve, at this averaged level, for exact discrete compacton solution configurations in the corresponding stationary equation. We focus on seven types of compacton solutions: single site and vortex solutions are found to be always stable in the parametric regimes we examined. Other solutions such as double site in- and out-of-phase, four site symmetric and anti-symmetric, and a five site compacton solution are found to have regions of stability and instability in two-dimensional parametric planes, involving variations of the strength of the coupling and of the nonlinearity. We also explore the time evolution of the solutions and compare the dynamics according to the averaged with those of the original dynamical equations without the averaging. Possible observation of compactons in the BEC loaded in a deep two-dimensional optical lattice with interactions modulated periodically in time is discussed.

preprint2015arXiv

Multifrequency and edge breathers in the discrete sine-Gordon system via subharmonic driving: theory, computation and experiment

We consider a chain of torsionally-coupled, planar pendula shaken horizontally by an external sinusoidal driver. It has been known that in such a system, theoretically modeled by the discrete sine-Gordon equation, intrinsic localized modes, also known as discrete breathers, can exist. Recently, the existence of multifrequency breathers via subharmonic driving has been theoretically proposed and numerically illustrated by Xu {\em et al.} in Phys. Rev. E {\bf 90}, 042921 (2014). In this paper, we verify this prediction experimentally. Comparison of the experimental results to numerical simulations with realistic system parameters (including a Floquet stability analysis), and wherever possible to analytical results (e.g. for the subharmonic response of the single driven-damped pendulum), yields good agreement. Finally, we report on period-1 and multifrequency edge breathers which are localized at the open boundaries of the chain, for which we have again found good agreement between experiments and numerical computations

preprint2015arXiv

Non-Conservative Variational Approximation for Nonlinear Schroedinger Equations

Recently, Galley [Phys. Rev. Lett. {\bf 110}, 174301 (2013)] proposed an initial value problem formulation of Hamilton's principle applied to non-conservative systems. Here, we explore this formulation for complex partial differential equations of the nonlinear Schrodinger (NLS) type, examining the dynamics of the coherent solitary wave structures of such models by means of a non-conservative variational approximation (NCVA). We compare the formalism of the NCVA to two other variational techniques used in dissipative systems; namely, the perturbed variational approximation and a generalization of the so-called Kantorovich method. All three variational techniques produce equivalent equations of motion for the perturbed NLS models studied herein. We showcase the relevance of the NCVA method by exploring test case examples within the NLS setting including combinations of linear and density dependent loss and gain. We also present an example applied to exciton polariton condensates that intrinsically feature loss and a spatially dependent gain term.

preprint2015arXiv

Nonlinear Resonances and Antiresonances of a Forced Sonic Vacuum

We consider a harmonically driven acoustic medium in the form of a (finite length) highly nonlinear granular crystal with an amplitude and frequency dependent boundary drive. Remarkably, despite the absence of a linear spectrum in the system, we identify resonant periodic propagation whereby the crystal responds at integer multiples of the drive period, and observe that they can lead to local maxima of transmitted force at its fixed boundary. In addition, we identify and discuss minima of the transmitted force ("antiresonances") in between these resonances. Representative one-parameter complex bifurcation diagrams involve period doublings, Neimark-Sacker bifurcations as well as multiple isolas (e.g. of period-3, -4 or -5 solutions entrained by the forcing). We combine them in a more detailed, two-parameter bifurcation diagram describing the stability of such responses to both frequency and amplitude variations of the drive. This picture supports an unprecedented example of a notion of a (purely) "nonlinear spectrum" in a system which allows no sound wave propagation (due to zero sound speed: the so-called sonic vacuum). We rationalize this behavior in terms of purely nonlinear building blocks: apparent traveling and standing nonlinear waves.

preprint2015arXiv

On the Bifurcation and Stability of Single and Multiple Vortex Rings in Three-Dimensional Bose-Einstein Condensates

In the present work, we investigate how single- and multi-vortex-ring states can emerge from a planar dark soliton in three-dimensional (3D) Bose-Einstein condensates (confined in isotropic or anisotropic traps) through bifurcations. We characterize such bifurcations quantitatively using a Galerkin-type approach, and find good qualitative and quantitative agreement with our Bogoliubov-de Gennes (BdG) analysis. We also systematically characterize the BdG spectrum of the dark solitons, using perturbation theory, and obtain a quantitative match with our 3D BdG numerical calculations. We then turn our attention to the emergence of single- and multi-vortex-ring states. We systematically capture these as stationary states of the system and quantify their BdG spectra numerically. We find that although the vortex ring may be unstable when bifurcating, its instabilities weaken and may even eventually disappear, for sufficiently large chemical potentials and suitable trap settings. For instance, we demonstrate the stability of the vortex ring for an isotropic trap in the large chemical potential regime.

preprint2015arXiv

Positive and negative mass solitons in spin-orbit coupled Bose-Einstein condensates

We present a unified description of different types of matter-wave solitons that can emerge in quasi one-dimensional spin-orbit coupled (SOC) Bose-Einstein condensates (BECs). This description relies on the reduction of the original two-component Gross-Pitaevskii SOC-BEC model to a single nonlinear Schrödinger equation, via a multiscale expansion method. This way, we find approximate bright and dark soliton solutions, for attractive and repulsive interatomic interactions respectively, for different regimes of the SOC interactions. Beyond this, our approach also reveals "negative mass" regimes, where corresponding "negative mass" bright or dark solitons can exist for repulsive or attractive interactions, respectively. Such a unique opportunity stems from the structure of the excitation spectrum of the SOC-BEC. Numerical results are found to be in excellent agreement with our analytical predictions.

preprint2015arXiv

PT Meets Supersymmetry and Nonlinearity: An Analytically Tractable Case Example

In the present work, we combine the notion of $\mathcal{PT}$-symmetry with that of super-symmetry (SUSY) for a prototypical case example with a complex potential that is related by SUSY to the so-called P{ö}schl-Teller potential which is real. Not only are we able to identify and numerically confirm the eigenvalues of the relevant problem, but we also show that the corresponding nonlinear problem, in the presence of an arbitrary power law nonlinearity, has an exact bright soliton solution that can be analytically identified and has intriguing stability properties, such as an oscillatory instability, which the corresponding solution of the regular nonlinear Schr{ö}dinger equation with arbitrary power law nonlinearity does not possess. The spectral properties and dynamical implications of this instability are examined. We believe that these findings may pave the way towards initiating a fruitful interplay between the notions of $\mathcal{PT}$-symmetry, super-symmetric partner potentials and nonlinear interactions.

preprint2015arXiv

Robust Vortex Lines, Vortex Rings and Hopfions in 3D Bose-Einstein Condensates

Performing a systematic Bogoliubov-de Gennes spectral analysis, we illustrate that stationary vortex lines, vortex rings and more exotic states, such as hopfions, are robust in three-dimensional atomic Bose-Einstein condensates, for large parameter intervals. Importantly, we find that the hopfion can be stabilized in a simple parabolic trap, without the need for trap rotation or inhomogeneous interactions. We supplement our spectral analysis by studying the dynamics of such stationary states; we find them to be robust against significant perturbations of the initial state. In the unstable regimes, we not only identify the unstable mode, such as a quadrupolar or hexapolar mode, but we also observe the corresponding instability dynamics. Furthermore, deep in the Thomas-Fermi regime, we investigate the particle-like behavior of vortex rings and hopfions.

preprint2015arXiv

Scattering of matter-waves in spatially inhomogeneous environments

We study scattering of quasi one-dimensional matter-waves at an interface of two spatial domains, one with repulsive and one with attractive interatomic interactions. It is shown that the incidence of a Gaussian wavepacket from the repulsive to the attractive region gives rise to generation of a soliton train. More specifically, the number of emergent solitons can be controlled e.g. by the variation of the amplitude or the width of the incoming wavepacket. Furthermore, we study the reflectivity of a soliton incident from the attractive region to the repulsive one. We find the reflection coefficient numerically and employ analytical methods, that treat the soliton as a particle (for moderate and large amplitudes) or a quasi-linear wavepacket (for small amplitudes), to determine the critical soliton momentum - as function of the soliton amplitude - for which total reflection is observed.

preprint2015arXiv

Solitons in Multi-Component Nonlinear Schrodinger Models: A Survey of Recent Developments

In this review we try to capture some of the recent excitement induced by experimental developments, but also by a large volume of theoretical and computational studies addressing multi-component nonlinear Schrodinger models and the localized structures that they support. We focus on some prototypical structures, namely the dark-bright and dark-dark solitons. Although our focus will be on one-dimensional, two-component Hamiltonian models, we also discuss variants, including three (or more)-component models, higher-dimensional states, as well as dissipative settings. We also offer an outlook on interesting possibilities for future work on this theme.

preprint2015arXiv

Some case example exact solutions for quadratically nonlinear optical media with $\mathcal{PT}$-symmetric potentials

In the present paper we consider an optical system with a $χ^{(2)}$-type nonlinearity and unspecified $\mathcal{PT}$-symmetric potential functions. Considering this as an inverse problem and positing a family of exact solutions in terms of cnoidal functions, we solve for the resulting potential functions in a way that ensures the potentials obey the requirements of $\mathcal{PT}$-symmetry. We then focus on case examples of soliton and periodic solutions for which we present a stability analysis as a function of their amplitude parameters. Finally, we numerically explore the nonlinear dynamics of the associated waveforms to identify the outcome of the relevant dynamical instabilities of localized and extended states.

preprint2015arXiv

Stabilization of ring dark solitons in Bose-Einstein condensates

Earlier work has shown that ring dark solitons in two-dimensional Bose-Einstein condensate are generically unstable. In this work, we propose a way of stabilizing the ring dark soliton via a radial Gaussian external potential. We investigate the existence and stability of the ring dark soliton upon variations of the chemical potential and also of the strength of the radial potential. Numerical results show that the ring dark soliton can be stabilized in a suitable interval of external potential strengths and chemical potentials. We also explore different proposed particle pictures considering the ring as a moving particle and find, where appropriate, results in very good qualitative and also reasonable quantitative agreement with the numerical findings.

preprint2015arXiv

Superdiffusive Transport and Energy Localization in Disordered Granular Crystals

We study the spreading of initially localized excitations in 1D disordered granular crystals. We thereby investigate localization phenomena in strongly nonlinear systems, which we demonstrate to be fundamentally different from localization in linear and weakly nonlinear systems. We compare wave dynamics in chains with 3 different types of disorder: an uncorrelated (Anderson-like) disorder and 2 types of correlated disorders (random dimer arrangements), and for 2 types of initial conditions: displacement excitations and velocity excitations. For strongly precompressed chains, the dynamics depend strongly on the initial condition. For displacement excitations, the long-time behavior of the second moment $\tilde{m}_2$ has oscillations that depend on the type of disorder, with a complex trend that differs markedly from a power law and which is particularly evident for an Anderson disorder. For velocity excitations, we find a scaling $\tilde{m}_2\sim t^γ$ (for a constant $γ$) for all 3 types of disorder. For weakly precompressed (strongly nonlinear) chains, $\tilde{m}_2$ and the inverse participation ratio $P^{-1}$ satisfy $\tilde{m}_2\sim t^γ$ and $P^{-1}\sim t^{-η}$, and the dynamics is superdiffusive for all examined cases. When precompression is strong, the IPR decreases slowly for all 3 types of disorder, and we observe a partial localization around the core and the leading edge of the wave. For an Anderson disorder, displacement perturbations lead to localization of energy primarily in the core, and velocity perturbations cause the energy to be divided between the core and the leading edge. This localization phenomenon does not occur in the sonic-vacuum regime, which yields the surprising result that the energy is no longer contained in strongly nonlinear waves but instead is spread across many sites. In this regime, the exponents are very similar in all cases.

preprint2014arXiv

A Tale of Two Distributions: From Few To Many Vortices In Quasi-Two-Dimensional Bose-Einstein Condensates

Motivated by the recent successes of particle models in capturing the precession and interactions of vortex structures in quasi-two-dimensional Bose-Einstein condensates, we revisit the relevant systems of ordinary differential equations. We consider the number of vortices $N$ as a parameter and explore the prototypical configurations ("ground states") that arise in the case of few or many vortices. In the case of few vortices, we modify the classical result of Havelock [Phil. Mag. ${\bf 11}$, 617 (1931)] illustrating that vortex polygons in the form of a ring are unstable for $N \geq7$. Additionally, we reconcile this modification with the recent identification of symmetry breaking bifurcations for the cases of $N=2,\dots,5$. We also briefly discuss the case of a ring of vortices surrounding a central vortex (so-called $N+1$ configuration). We finally examine the opposite limit of large $N$ and illustrate how a coarse-graining, continuum approach enables the accurate identification of the radial distribution of vortices in that limit.

preprint2014arXiv

Beating dark-dark solitons and Zitterbewegung in spin-orbit coupled Bose-Einstein condensates

We present families of beating dark-dark solitons in spin-orbit (SO) coupled Bose-Einstein condensates. These families consist of solitons residing simultaneously in the two bands of the energy spectrum. The soliton components are characterized by two different spatial and temporal scales, which are identified by a multiscale expansion method. The solitons are "beating" ones, as they perform density oscillations with a characteristic frequency, relevant to Zitterbewegung (ZB). We find that spin oscillations may occur, depending on the parity of each soliton branch, which consequently lead to ZB oscillations of the beating dark solitons. Analytical results are corroborated by numerical simulations, illustrating the robustness of the solitons.

preprint2014arXiv

Characterizing Traveling Wave Collisions in Granular Chains Starting from Integrable Limits: the case of the KdV and the Toda Lattice

Our aim in the present work is to develop approximations for the collisional dynamics of traveling waves in the context of granular chains in the presence of precompression. To that effect, we aim to quantify approximations of the relevant Hertzian FPU-type lattice through both the Korteweg-de Vries (KdV) equation and the Toda lattice. Using the availability in such settings of both 1-soliton and 2-soliton solutions in explicit analytical form, we initialize such coherent structures in the granular chain and observe the proximity of the resulting evolution to the underlying integrable (KdV or Toda) model. While the KdV offers the possibility to accurately capture collisions of solitary waves propagating in the same direction, the Toda lattice enables capturing both co-propagating and counter-propagating soliton collisions. The error in the approximation is quantified numerically and connections to bounds established in the mathematical literature are also given.

preprint2014arXiv

Dark-bright Solitons and their Lattices in Atomic Bose-Einstein Condensates

In the present contribution, we explore a host of different stationary states, namely dark-bright solitons and their lattices, that arise in the context of multi-component atomic Bose-Einstein condensates. The latter, are modeled by systems of coupled Gross-Pitaevskii equations with general interaction (nonlinearity) coefficients $g_{ij}$. It is found that in some particular parameter ranges such solutions can be obtained in analytical form, however, numerically they are computed as existing in a far wider parametric range. Many features of the solutions under study, such as their analytical form without the trap or the stability/dynamical properties of one dark-bright soliton even in the presence of the trap are obtained analytically and corroborated numerically. Additional features, such as the stability of soliton lattice homogeneous states or their existence/stability in the presence of the trap, are examined numerically.

preprint2014arXiv

Dynamic and Energetic Stabilization of Persistent Currents in Bose-Einstein Condensates

We study conditions under which vortices in a highly oblate harmonically trapped Bose-Einstein condensate (BEC) can be stabilized due to pinning by a blue-detuned Gaussian laser beam, with particular emphasis on the potentially destabilizing effects of laser beam positioning within the BEC. Our approach involves theoretical and numerical exploration of dynamically and energetically stable pinning of vortices with winding number up to $S=6$, in correspondence with experimental observations. Stable pinning is quantified theoretically via Bogoliubov-de Gennes excitation spectrum computations and confirmed via direct numerical simulations for a range of conditions similar to those of experimental observations. The theoretical and numerical results indicate that the pinned winding number, or equivalently the winding number of the superfluid current about the laser beam, decays as a laser beam of fixed intensity moves away from the BEC center. Our theoretical analysis helps explain previous experimental observations, and helps define limits of stable vortex pinning for future experiments involving vortex manipulation by laser beams.

preprint2014arXiv

Exploring Vortex Dynamics in the Presence of Dissipation: Analytical and Numerical Results

In this paper, we systematically examine the stability and dynamics of vortices under the effect of a phenomenological dissipation used as a simplified model for the inclusion of the effect of finite temperatures in atomic Bose-Einstein condensates. An advantage of this simplified model is that it enables an analytical prediction that can be compared directly (and favorably) to numerical results. We then extend considerations to a case of considerable recent experimental interest, namely that of a vortex dipole and observe good agreement between theory and numerical computations in both the stability properties (eigenvalues of the vortex dipole stationary states) and the dynamical evolution of such configurations.

preprint2014arXiv

Generalized Dimers and their Stokes-variable Dynamics

In the present work, we generalize the setting of dimers with potential gain and loss which have been extensively considered recently in $\mathcal{P T}$-symmetric contexts. We consider a pair of waveguides which are evanescently coupled but may also be actively coupled and may possess onsite gain and loss, as well as (possibly non-uniform) nonlinearity. We identify (and where appropriate review from earlier work) a plethora of interesting dynamical scenaria ranging from the existence of stable and unstable fixed points and integrable dynamics, to the emergence of pitchfork or Hopf bifurcations and the generation of additional fixed points and limit cycles, respectively, as well as the potential deviation of trajectories to infinity. Thus, a catalogue of a large number of possible cases is given and their respective settings physically justified (where appropriate).

preprint2014arXiv

Generating Functions, Polynomials and Vortices with Alternating Signs in Bose-Einstein Condensates

In this work, we construct suitable generating functions for vortices of alternating signs in the realm of Bose-Einstein condensates. In addition to the vortex-vortex interaction included in earlier fluid dynamics constructions of such functions, the vortices here precess around the center of the trap. This results in the generating functions of the vortices of positive charge and of negative charge satisfying a modified, so-called, Tkachenko differential equation. From that equation, we reconstruct collinear few-vortex equilibria obtained in earlier work, as well as extend them to larger numbers of vortices. Moreover, particular moment conditions can be derived e.g. about the sum of the squared locations of the vortices for arbitrary vortex numbers. Furthermore, the relevant differential equation can be generalized appropriately in the two-dimensional complex plane and allows the construction e.g. of polygonal vortex ring and multi-ring configurations, as well as ones with rings surrounding a vortex at the center that are again connected to earlier bibliography.

preprint2014arXiv

Highly Nonlinear Wave Propagation in Elastic Woodpile Periodic Structures

In the present work, we experimentally implement, numerically compute with and theoretically analyze a configuration in the form of a single column woodpile periodic structure. Our main finding is that a Hertzian, locally-resonant, woodpile lattice offers a test bed for the formation of genuinely traveling waves composed of a strongly-localized solitary wave on top of a small amplitude oscillatory tail. This type of wave, called a nanopteron, is not only motivated theoretically and numerically, but are also visualized experimentally by means of a laser Doppler vibrometer. This system can also be useful for manipulating stress waves at will, for example, to achieve strong attenuation and modulation of high-amplitude impacts without relying on damping in the system.

preprint2014arXiv

Instability dynamics and breather formation in a horizontally shaken pendulum chain

Inspired by the experimental results of Cuevas et al. (Physical Review Letters 102, 224101 (2009)), we consider theoretically the behavior of a chain of planar rigid pendulums suspended in a uniform gravitational field and subjected to a horizontal periodic driving force applied to the pendulum pivots. We characterize the motion of a single pendulum, finding bistability near the fundamental resonance, and near the period-3 subharmonic resonance. We examine the development of modulational instability in a driven pendulum chain and find both a critical chain length and critical frequency for the appearance of the instability. We study the breather solutions and show their connection to the single pendulum dynamics, and extend our analysis to consider multi-frequency breathers connected to the period-3 periodic solution, showing also the possibility of stability in these breather states. Finally we examine the problem of breather generation and demonstrate a robust scheme for generation of on-site and off-site breathers.

preprint2014arXiv

Interaction of sine-Gordon kinks and breathers with a $\mathcal{PT}$-symmetric defect

The scattering of kinks and low-frequency breathers of the nonlinear sine-Gordon (SG) equation on a spatially localized $\mathcal{PT}$-symmetric perturbation (defect) with a balanced gain and loss is investigated numerically. It is demonstrated that if a kink passes the defect, it always restores its initial momentum and energy and the only effect of the interaction with the defect is a phase shift of the kink. A kink approaching the defect from the gain side always passes, while in the opposite case it must have sufficiently large initial momentum to pass through the defect instead of being trapped in the loss region. The kink phase shift and critical velocity are calculated with the use of the collective variable method. Kink-kink (kink-antikink) collisions at the defect are also briefly considered, showing how their pairwise repulsive (respectively, attractive) interaction can modify the collisional outcome of a single kink within the pair with the defect. For the breather, the result of its interaction with the defect strongly depends on the breather parameters (velocity, frequency and initial phase) and on the defect parameters. The breather can gain some energy from the defect and as a result potentially even split into a kink-antikink pair or it can lose a part of its energy. Interestingly, the breather translational mode is very weakly affected by the dissipative perturbation, so that a breather penetrates more easily through the defect when it comes from the lossy side, than a kink. In all studied soliton-defect interactions the energy loss to radiation of small-amplitude extended waves is negligible.

preprint2014arXiv

Lattice Three-Dimensional Skyrmions Revisited

In the continuum a skyrmion is a topological nontrivial map between Riemannian manifolds, and a stationary point of a particular energy functional. This paper describes lattice analogues of the aforementioned skyrmions, namely a natural way of using the topological properties of the three-dimensional continuum Skyrme model to achieve topological stability on the lattice. In particular, using fixed point iterations, numerically exact lattice skyrmions are constructed; and their stability under small perturbations is explored by means of linear stability analysis. While stable branches of such solutions are identified, it is also shown that they possess a particularly delicate bifurcation structure, especially so in the vicinity of the continuum limit. The corresponding bifurcation diagram is elucidated and a prescription for selecting the branch asymptoting to the well-known continuum limit is given. Finally, the robustness of the spectrally stable solutions is corroborated by virtue of direct numerical simulations .

preprint2014arXiv

Nonlinear dressed states at the miscibility-immiscibility threshold

The dynamical evolution of spatial patterns in a complex system can reveal the underlying structure and stability of stationary states. As a model system we employ a two-component rubidium Bose-Einstein condensate at the transition from miscible to immiscible with the additional control of linear interconversion. Excellent agreement is found between the detailed experimental time evolution and the corresponding numerical mean-field computations. Analyzing the dynamics of the system, we find clear indications of stationary states that we term nonlinear dressed states. A steady state bifurcation analysis reveals a smooth connection of these states with dark-bright soliton solutions of the integrable two-component Manakov model.

preprint2014arXiv

On the spectral stability of kinks in some PT-symmetric variants of the classical Klein-Gordon Field Theories

In the present work we consider the introduction of PT-symmetric terms in the context of classical Klein-Gordon field theories. We explore the implication of such terms on the spectral stability of coherent structures, namely kinks. We find that the conclusion critically depends on the location of the kink center relative to the center of the PT-symmetric term. The main result is that if these two points coincide, the kink's spectrum remains on the imaginary axis and the wave is spectrally stable. If the kink is centered on the "lossy side" of the medium, then it becomes stabilized. On the other hand, if it becomes centered on the "gain side" of the medium, then it is destabilized. The consequences of these two possibilities on the linearization (point and essential) spectrum are discussed in some detail.

preprint2014arXiv

PT-Symmetric dimer in a generalized model of coupled nonlinear oscillators

In the present work, we explore the case of a general PT-symmetric dimer in the context of two both linearly and nonlinearly coupled cubic oscillators. To obtain an analytical handle on the system, we first explore the rotating wave approximation converting it into a discrete nonlinear Schrödinger type dimer. In the latter context, the stationary solutions and their stability are identified numerically but also wherever possible analytically. Solutions stemming from both symmetric and anti-symmetric special limits are identified. A number of special cases are explored regarding the ratio of coefficients of nonlinearity between oscillators over the intrinsic one of each oscillator. Finally, the considerations are extended to the original oscillator model, where periodic orbits and their stability are obtained. When the solutions are found to be unstable their dynamics is monitored by means of direct numerical simulations.

preprint2014arXiv

PT-symmetric sine-Gordon breathers

In this work, we explore a prototypical example of a genuine continuum breather (i.e., not a standing wave) and the conditions under which it can persist in a $\mathcal{P T}$-symmetric medium. As our model of interest, we will explore the sine-Gordon equation in the presence of a $\mathcal{P T}$- symmetric perturbation. Our main finding is that the breather of the sine-Gordon model will only persist at the interface between gain and loss that $\mathcal{P T}$-symmetry imposes but will not be preserved if centered at the lossy or at the gain side. The latter dynamics is found to be interesting in its own right giving rise to kink-antikink pairs on the gain side and complete decay of the breather on the lossy side. Lastly, the stability of the breathers centered at the interface is studied. As may be anticipated on the basis of their "delicate" existence properties such breathers are found to be destabilized through a Hopf bifurcation in the corresponding Floquet analysis.

preprint2014arXiv

Quasi-energies, parametric resonances, and stability limits in ac-driven $\mathcal{PT}$-symmetric systems

We introduce a simple model for implementing the concepts of quasi-energy and parametric resonances (PRs) in systems with the $\mathcal{PT}$ symmetry, i.e., a pair of coupled and mutually balanced gain and loss elements. The parametric (ac) forcing is applied through periodic modulation of the coefficient accounting for the coupling of the two degrees of freedom. The system may be realized in optics as a dual-core waveguide with the gain and loss applied to different cores, and the thickness of the gap between them subject to a periodic modulation. The onset and development of the parametric instability for a small forcing amplitude ($V_{1}$) is studied in an analytical form. The full dynamical chart of the system is generated by systematic simulations. At sufficiently large values of the forcing frequency, $ω$, tongues of the parametric instability originate, with the increase of $V_{1}$, as predicted by the analysis. However, the tongues following further increase of $V_{1}$ feature a pattern drastically different from that in usual (non-$\mathcal{PT}$)\ parametrically driven systems: instead of bending down to larger values of the dc coupling constant, $V_{0}$, they maintain a direction parallel to the $V_{1}$ axis. The system of the parallel tongues gets dense with the decrease of $ω$, merging into a complex small-scale structure of alternating regions of stability and instability. The cases of $ω\rightarrow 0$ and $ω\rightarrow \infty $ are studied analytically by means of the adiabatic and averaging approximation, respectively. The cubic nonlinearity, if added to the system, alters the picture, destabilizing many originally robust dynamical regimes, and stabilizing some which were unstable.

preprint2014arXiv

Rogue Waves In Nonlinear Schrödinger Models with Variable Coefficients: Application to Bose-Einstein Condensates

We explore the form of rogue wave solutions in a select set of case examples of nonlinear Schrödinger equations with variable coefficients. We focus on systems with constant dispersion, and present three different models that describe atomic Bose-Einstein condensates in different experimentally relevant settings. For these models, we identify exact rogue wave solutions. Our analytical findings are corroborated by direct numerical integration of the original equations, performed by two different schemes. Very good agreement between numerical results and analytical predictions for the emergence of the rogue waves is identified. Additionally, the nontrivial fate of small numerically induced perturbations to the exact rogue wave solutions is also discussed.

preprint2014arXiv

Scattering and leapfrogging of vortex rings in a superfluid

The dynamics of vortex ring pairs in the homogeneous nonlinear Schrödinger equation is studied. The generation of numerically-exact solutions of traveling vortex rings is described and their translational velocity compared to revised analytic approximations. The scattering behavior of co-axial vortex rings with opposite charge undergoing collision is numerically investigated for different scattering angles yielding a surprisingly simple result for its dependence as a function of the initial vortex ring parameters. We also study the leapfrogging behavior of co-axial rings with equal charge and compare it with the dynamics stemming from a modified version of the reduced equations of motion from a classical fluid model derived using the Biot-Savart law.

preprint2014arXiv

Solitary Waves in a Discrete Nonlinear Dirac equation

In the present work, we introduce a discrete formulation of the nonlinear Dirac equation in the form of a discretization of the Gross-Neveu model. The motivation for this discrete model proposal is both computational (near the continuum limit) and theoretical (using the understanding of the anti-continuum limit of vanishing coupling). Numerous unexpected features are identified including a staggered solitary pattern emerging from a single site excitation, as well as two- and three-site excitations playing a role analogous to one- and two-site, respectively, excitations of the discrete nonlinear Schrödinger analogue of the model. Stability exchanges between the two- and three-site states are identified, as well as instabilities that appear to be persistent over the coupling strength $ε$, for a subcritical value of the propagation constant $Λ$. Variations of the propagation constant, coupling parameter and nonlinearity exponent are all examined in terms of their existence and stability implications and long dynamical simulations are used to unravel the evolutionary phenomenology of the system (when unstable).

preprint2014arXiv

Solitons and Vortices in Two-dimensional Discrete Nonlinear Schrodinger Systems with Spatially Modulated Nonlinearity

We consider a two-dimensional (2D) generalization of a recently proposed model [Phys. Rev. E 88, 032905 (2013)], which gives rise to bright discrete solitons supported by the defocusing nonlinearity whose local strength grows from the center to the periphery. We explore the 2D model starting from the anti-continuum (AC) limit of vanishing coupling. In this limit, we can construct a wide variety of solutions including not only single-site excitations, but also dipole and quadrupole ones. Additionally, two separate families of solutions are explored: the usual "extended" unstaggered bright solitons, in which all sites are excited in the AC limit, with the same sign across the lattice (they represent the most robust states supported by the lattice, their 1D counterparts being what was considered as 1D bright solitons in the above-mentioned work), and the vortex cross, which is specific to the 2D setting. For all the existing states, we explore their stability (analytically, whenever possible). Typical scenarios of instability development are exhibited through direct simulations.

preprint2014arXiv

Stability and Tunneling Dynamics of a Dark-Bright Soliton Pair in a Harmonic Trap

We consider a binary repulsive Bose-Einstein condensate in a harmonic trap in one spatial dimension and investigate particular solutions consisting of two dark-bright (DB) solitons. There are two different stationary solutions characterized by the phase difference in the bright component, in-phase and out-of-phase states. We show that above a critical particle number in the bright component, a symmetry breaking bifurcation of the pitchfork type occurs that leads to a new asymmetric solution whereas the parental branch, i.e., the out-of-phase state becomes unstable. These three different states support different small amplitude oscillations, characterized by an almost stationary density of the dark component and a tunneling of the bright component between the two dark solitons. Within a suitable effective double-well picture, these can be understood as the characteristic features of a Bosonic Josephson Junction (BJJ), and we show within a two-mode approach that all characteristic features of the BJJ phase space are recovered. For larger deviations from the stationary states, the simplifying double-well description breaks down due to the feedback of the bright component onto the dark one, causing the solitons to move. In this regime we observe intricate anharmonic and aperiodic dynamics, exhibiting remnants of the BJJ phase space.

preprint2014arXiv

Time-Periodic Solutions of Driven-Damped Trimer Granular Crystals

In this work, we consider time-periodic structures of trimer granular crystals consisting of alternate chrome steel and tungsten carbide spherical particles yielding a spatial periodicity of three. The configuration at the left boundary is driven by a harmonic in-time actuation with given amplitude and frequency while the right one is a fixed wall. Similar to the case of a dimer chain, the combination of dissipation, driving of the boundary, and intrinsic nonlinearity leads to complex dynamics. For fixed driving frequencies in each of the spectral gaps, we find that the nonlinear surface modes and the states dictated by the linear drive collide in a saddle-node bifurcation as the driving amplitude is increased, beyond which the dynamics of the system become chaotic. While the bifurcation structure is similar for solutions within the first and second gap, those in the first gap appear to be less robust. We also conduct a continuation in driving frequency, where it is apparent that the nonlinearity of the system results in a complex bifurcation diagram, involving an intricate set of loops of branches, especially within the spectral gap. The theoretical findings are qualitatively corroborated by the experimental full-field visualization of the time-periodic structures.

preprint2014arXiv

Transitions From Order to Disorder in Multi-Dark and Multi-Dark-Bright Soliton Atomic Clouds

We have performed a systematic study quantifying the variation of solitary wave behavior from that of an ordered cloud resembling a "crystalline" configuration to that of a disordered state that can be characterized as a soliton "gas". As our illustrative examples, we use both one-component, as well as two-component, one dimensional atomic gases very close to zero temperature, where in the presence of repulsive inter-atomic interactions and of a parabolic trap, a cloud, respectively of dark (dark-bright) solitons can form in the one- (two-) component system. We corroborate our findings through three distinct types of approaches, namely a Gross-Pitaevskii type of partial differential equation, particle-based ordinary differential equations describing the soliton dynamical system and Monte-Carlo simulations for the particle system. We define an "empirical" order parameter to characterize the order of the soliton lattices and study how this changes as a function of the strength of the "thermally" (i.e., kinetically) induced perturbations. As may be anticipated by the one-dimensional nature of our system, the transition from order to disorder is gradual without, apparently, a genuine phase transition ensuing in the intermediate regime.

preprint2014arXiv

Traveling Waves and their Tails in Locally Resonant Granular Systems

In the present study, we revisit the theme of wave propagation in locally resonant granular crystal systems, also referred to as Mass-in-Mass systems. We use 3 distinct approaches to identify relevant traveling waves. The first consists of a direct solution of the traveling wave problem. The second one consists of the solution of the Fourier tranformed variant of the problem. or, more precisely, of its convolution reformulation (upon an inverse Fourier transform) of the problem in real space. Finally, our third approach will restrict considerations to a finite domain, utilizing the notion of Fourier series for important technical reasons, namely the avoidance of resonances, that will be discussed in detail. All three approaches can be utilized in either the displacement or the strain formulation. Typical resulting computations in finite domains result in the solitary waves bearing symmetric non-vanishing tails at both ends of the computational domain. Importantly, however, a countably infinite set of resonance conditions is identified for which solutions with genuinely monotonic decaying tails arise.

preprint2014arXiv

Traveling waves of the regularized short pulse equation

In the present work, we revisit the so-called regularized short pulse equation (RSPE) and, in particular, explore the traveling wave solutions of this model. We theoretically analyze and numerically evolve two sets of such solutions. First, using a fixed point iteration scheme, we numerically integrate the equation to find solitary waves. It is found that these solutions are well approximated by a truncated series of hyperbolic secants. The dependence of the soliton's parameters (height, width, etc) to the parameters of the equation is also investigated. Second, by developing a multiple scale reduction of the RSPE to the nonlinear Schrödinger equation, we are able to construct (both standing and traveling) envelope wave breather type solutions of the former, based on the solitary wave structures of the latter. Both the regular and the breathing traveling wave solutions identified are found to be robust and should thus be amenable to observations in the form of few optical cycle pulses.

preprint2014arXiv

Vector rogue waves and dark-bright boomeronic solitons in autonomous and non-autonomous settings

In this work, we consider the dynamics of vector rogue waves and dark-bright solitons in two-component nonlinear Schrödinger equations with various physically motivated time-dependent nonlinearity coefficients, as well as spatio-temporally dependent potentials. A similarity transformation is utilized to convert the system into the integrable Manakov system and subsequently the vector rogue and dark-bright boomeron-like soliton solutions of the latter are converted back into ones of the original non-autonomous model. Using direct numerical simulations we find that, in most cases, the rogue wave formation is rapidly followed by a modulational instability that leads to the emergence of an expanding soliton train. Scenarios different than this generic phenomenology are also reported.

preprint2014arXiv

Vortex Nucleation in a Dissipative Variant of the Nonlinear Schrödinger Equation under Rotation

In the present work, we motivate and explore the dynamics of a dissipative variant of the nonlinear Schr{ö}dinger equation under the impact of external rotation. As in the well established Hamiltonian case, the rotation gives rise to the formation of vortices. We show, however, that the most unstable mode leading to this instability scales with an appropriate power of the chemical potential $μ$ of the system, increasing proportionally to $μ^{2/3}$. The precise form of the relevant formula, obtained through our asymptotic analysis, provides the most unstable mode as a function of the atomic density and the trap strength. We show how these unstable modes typically nucleate a large number of vortices in the periphery of the atomic cloud. However, through a pattern selection mechanism, prompted by symmetry-breaking, only few isolated vortices are pulled in sequentially from the periphery towards the bulk of the cloud resulting in highly symmetric stable vortex configurations with far fewer vortices than the original unstable mode. These results may be of relevance to the experimentally tractable realm of finite temperature atomic condensates.

preprint2014arXiv

When Does Linear Stability Not Exclude Nonlinear Instability ?

We describe a mechanism that results in the nonlinear instability of stationary states even in the case where the stationary states are linearly stable. This instability is due to the nonlinearity-induced coupling of the linearization's internal modes of negative energy with the wave continuum. In a broad class of nonlinear Schr{ö}dinger (NLS) equations considered, the presence of such internal modes guarantees the nonlinear instability of the stationary states in the evolution dynamics. To corroborate this idea, we explore three prototypical case examples: (a) an anti-symmetric soliton in a double-well potential, (b) a twisted localized mode in a one-dimensional lattice with cubic nonlinearity, and (c) a discrete vortex in a two-dimensional saturable lattice. In all cases, we observe a weak nonlinear instability, despite the linear stability of the respective states.

preprint2013arXiv

${\cal PT}$-symmetric coupler with $χ^{(2)}$ nonlinearity

We introduce the notion of a ${\cal PT}$-symmetric dimer with a $χ^{(2)}$ nonlinearity. Similarly to the Kerr case, we argue that such a nonlinearity should be accessible in a pair of optical waveguides with quadratic nonlinearity and gain and loss, respectively. An interesting feature of the problem is that because of the two harmonics, there exist in general two distinct gain/loss parameters, different values of which are considered herein. We find a number of traits that appear to be absent in the more standard cubic case. For instance, bifurcations of nonlinear modes from the linear solutions occur in two different ways depending on whether the first or the second harmonic amplitude is vanishing in the underlying linear eigenvector. Moreover, a host of interesting bifurcation phenomena appear to occur including saddle-center and pitchfork bifurcations which our parametric variations elucidate. The existence and stability analysis of the stationary solutions is corroborated by numerical time-evolution simulations exploring the evolution of the different configurations, when unstable.

preprint2013arXiv

Coupled backward- and forward-propagating solitons in a composite right/left-handed transmission line

We study the coupling between backward- and forward-propagating wave modes, with the same group velocity, in a composite right/left-handed nonlinear transmission line. Using an asymptotic multiscale expansion technique, we derive a system of two coupled nonlinear Schr{ö}dinger equations governing the evolution of the envelopes of these modes. We show that this system supports a variety of backward- and forward propagating vector solitons, of the bright-bright, bright-dark and dark-bright type. Performing systematic numerical simulations in the framework of the original lattice that models the transmission line, we study the propagation properties of the derived vector soliton solutions. We show that all types of the predicted solitons exist, but differ on their robustness: only bright-bright solitons propagate undistorted for long times, while the other types are less robust, featuring shorter lifetimes. In all cases, our analytical predictions are in a very good agreement with the results of the simulations, at least up to times of the order of the solitons' lifetimes.

preprint2013arXiv

Damped-Driven Granular Crystals: An Ideal Playground for Dark Breathers and Multibreathers

By applying an out-of-phase actuation at the boundaries of a uniform chain of granular particles, we demonstrate experimentally that time-periodic and spatially localized structures with a nonzero background (so-called dark breathers) emerge for a wide range of parameter values and initial conditions. Importantly, the number of ensuing breathers within the multibreather pattern produced can be "dialed in" by varying the frequency or amplitude of the actuation. The values of the frequency (resp. amplitude) where the transition between different multibreather states occurs are predicted accurately by the proposed theoretical model, which is numerically shown to support exact dark breather solutions. The existence, linear stability, and bifurcation structure of the theoretical dark breathers are also studied in detail. Moreover, the distributed sensing technologies developed herein enable a detailed space-time probing of the system and a systematic favorable comparison between theory, computation and experiments.

preprint2013arXiv

Directed Ratchet Transport in Granular Crystals

Directed-ratchet transport (DRT) in a one-dimensional lattice of spherical beads, which serves as a prototype for granular crystals, is investigated. We consider a system where the trajectory of the central bead is prescribed by a biharmonic forcing function with broken time-reversal symmetry. By comparing the mean integrated force of beads equidistant from the forcing bead, two distinct types of directed transport can be observed ---spatial and temporal DRT. Based on the value of the frequency of the forcing function relative to the cutoff frequency, the system can be categorized by the presence and magnitude of each type of DRT. Furthermore, we investigate and quantify how varying additional parameters such as the biharmonic weight affects DRT velocity and magnitude. Finally, friction is introduced into the system and is found to significantly inhibit spatial DRT. In fact, for sufficiently low forcing frequencies, the friction may even induce a switching of the DRT direction.

preprint2013arXiv

Dynamics of Few Co-rotating Vortices in Bose-Einstein Condensates

We study the dynamics of small vortex clusters with few (2--4) co-rotating vortices in Bose-Einstein condensates by means of experiments, numerical computations, and theoretical analysis. All of these approaches corroborate the counter-intuitive presence of a dynamical instability of symmetric vortex configurations. The instability arises as a pitchfork bifurcation at sufficiently large values of the angular momentum that induces the emergence and stabilization of asymmetric rotating vortex configurations. The latter are quantified in the theoretical model and observed in the experiments. The dynamics is explored both for the integrable two-vortex system, where a reduction of the phase space of the system provides valuable insight, as well as for the non-integrable three- (or more) vortex case, which additionally admits the possibility of chaotic trajectories.

preprint2013arXiv

Eigenstates and instabilities of chains with embedded defects

We consider the eigenvalue problem for one-dimensional linear Schrödinger lattices (tight-binding) with an embedded few-sites linear or nonlinear, Hamiltonian or non-conservative defect (an oligomer). Such a problem arises when considering scattering states in the presence of (generally complex) impurities as well as in the stability analysis of nonlinear waves. We describe a general approach based on a matching of solutions of the linear portions of the lattice at the location of the oligomer defect. As specific examples we discuss both linear and nonlinear, Hamiltonian and $\cP \cT$-symmetric dimers and trimers. In the linear case, this approach provides us a handle for semi-analytically computing the spectrum [this amounts to the solution of a polynomial equation]. In the nonlinear case, it enables the computation of the linearization spectrum around the stationary solutions. The calculations showcase the oscillatory instabilities that strongly nonlinear states typically manifest.

preprint2013arXiv

Few-Particle Vortex Cluster Equilibria in Bose-Einstein Condensates: Existence and Stability

Motivated by recent experimental and theoretical studies of fewparticle vortex clusters in Bose-Einstein condensates, we consider the ordinary differential equations of motion and systematically examine settings for up to N = 6 vortices. We analyze the existence of corresponding stationary state configurations and also consider their spectral stability properties. We compare our particle model results with the predictions of the full partial differential equation system. Whenever possible, we propose generalizations of these results in the context of clusters of N vortices. Some of these, we can theoretically establish, especially so for the N-vortex polygons, while others we state as conjectures, e.g. for the N-vortex line equilibrium.

preprint2013arXiv

Going beyond the double well: complex mode dynamics of effective coupled oscillators in infinite dimensional systems

In this work we explore how nonlinear modes described by a dispersive wave equation (in our example, the nonlinear Schrodinger equation) and localized in a few wells of a periodic potential can act analogously to a chain of coupled mechanical oscillators. We identify the small-amplitude oscillation modes of these `coupled wave oscillators' and find that they can be extended into the large amplitude regime, where some can `ring' for long times. We also identify prototypical case examples of more complex dynamical behaviour that can arise in such systems beyond the double well paradigm, including the breakdown of Josephson-like oscillations and of internal modes more generally, the transfer of energy out of/destabilization of fundamental oscillation modes and the emergence of chaotic oscillations for large amplitude excitations. We provide details of the phase perturbations required for experimental observations of such dynamics, and show that the oscillator formalism can be extended to predict large amplitude excitations in genuinely two-dimensional configurations.

preprint2013arXiv

Inelastic Collisions of Solitary Waves in Anisotropic Bose-Einstein Condensates: Sling-Shot Events and Expanding Collision Bubbles

We study experimentally and theoretically the dynamics of apparent dark soliton stripes in an elongated Bose-Einstein condensate. We show that for the trapping strengths corresponding to our experimental setup, the transverse confinement along one of the tight directions is not strong enough to arrest the formation of solitonic vortices or vortex rings. These solitonic vortices and vortex rings, when integrated along the transverse direction, appear as dark soliton stripes along the longitudinal direction thereby hiding their true character. The latter significantly modifies the interaction dynamics during collision events and can lead to apparent examples of inelasticity and what may appear experimentally even as a merger of two dark soliton stripes. We explain this feature by means of the interaction of two solitonic vortices leading to a sling shot event with one of the solitonic vortices being ejected at a relatively large speed. Furthermore we observe expanding collision bubbles which consist of repeated inelastic collisions of a dark soliton stripe pair with an {\it increasing} time interval between collisions.

preprint2013arXiv

Interactions of solitons with a Gaussian barrier: Splitting and recombination in quasi-1D and 3D

The interaction of matter-wave solitons with a potential barrier is a fundamentally important problem, and the splitting and subsequent recombination of the soliton by the barrier is the essence of soliton matter-wave interferometry. We demonstrate the three-dimensional (3D) character of the interactions in the case relevant to ongoing experiments, where the number of atoms in the soliton is sufficiently large that the soliton is close to collapse. The mean-field description is quite accurate, but the proximity to the collapse threshold makes the use of the 1D Gross-Pitaevskii equation (GPE) irrelevant. We examine the soliton dynamics in the framework of the effectively 1D nonpolynomial Schr{ö}dinger equation (NPSE), which admits the collapse in a modified form, and in parallel we use the full 3D GPE. Both approaches produce results which are very different from those produced in recent work which used the 1D cubic GPE. The features observed in both the NPSE and the 3D GPE include (a) an increase in the first reflection coefficient for increasing barrier height and decreasing atom number; (b) strong modulation in the secondary reflection/recombination probability vs. barrier height; (c) a pronounced asymmetry in the oscillation amplitudes of the transmitted and reflected fragments; (d) an enhancement of the transverse excitations as the number of atoms is increased. We also explore effects produced by variations of the barrier width and outcomes of the secondary collision upon phase imprinting on the fragment in one arm of the interferometer.

preprint2013arXiv

Justification of Leading Order Quasicontinuum Approximations of Strongly Nonlinear Lattices

We consider the leading order quasicontinuum limits of a one-dimensional granular medium governed by the Hertz contact law under precompression. The approximate model which is derived in this limit is justified by establishing asymptotic bounds for the error with the help of energy estimates. The continuum model predicts the development of shock waves, which are also studied in the full system with the aid of numerical simulations. We also show that existing results concerning the Nonlinear Schrodinger (NLS) and Korteweg de-Vries (KdV) approximation of FPU models apply directly to a precompressed granular medium in the weakly nonlinear regime.

preprint2013arXiv

Matter-wave dark solitons and their excitation spectra in spin-orbit coupled Bose-Einstein condensates

We present three types of dark solitons in quasi-one-dimensional spin-orbit coupled repulsive Bose-Einstein condensates. Among these families, two are always stable, while the third one is only stable sufficiently close to the linear regime. The solitons' excitation spectra reveal the potential existence of a second anomalous mode. While the first such mode describes the soliton oscillatory motion in a parabolic trap, the second, when present, reflects the double well structure of the underlying single-particle spectrum. This novel mode results in moving density stripes in the vicinity of the soliton core, or in an out-of-phase oscillation of the constituent components, with little effect on the nearly stationary striped total density of the composite soliton.

preprint2013arXiv

Matter-wave solitons in the counterflow of two immiscible superfluids

We study formation of solitons induced by counterflows of immiscible superfluids. Our setting is based on a quasi-one-dimensional binary Bose-Einstein condensate (BEC), composed of two immiscible components with large and small numbers of atoms in them. Assuming that the "small" component moves with constant velocity, either by itself, or being dragged by a moving trap, and intrudes into the "large" counterpart, the following results are obtained. Depending on the velocity, and on whether the small component moves in the absence or in the presence of the trap, two-component dark-bright solitons, scalar dark solitons, or multiple dark solitons may emerge, the latter outcome taking place due to breakdown of the superfluidity. We present two sets of analytical results to describe this phenomenology. In an intermediate velocity regime, where dark-bright solitons form, a reduction of the two-component Gross-Pitaevskii system to an integrable Mel'nikov system is developed, demonstrating that solitary waves of the former are very accurately described by analytically available solitons of the latter. In the high-velocity regime, where the breakdown of the superfluidity induces the formation of dark solitons and multi-soliton trains, an effective single-component description, in which a strongly localized wave packet of the "small" component acts as an effective potential for the "large" one, allows us to estimate the critical velocity beyond which the coherent structures emerge in good agreement with the numerical results.

preprint2013arXiv

Nonlinear localized modes in two-dimensional electrical lattices

We report the observation of spontaneous localization of energy in two spatial dimensions in the context of nonlinear electrical lattices. Both stationary and traveling self-localized modes were generated experimentally and theoretically in a family of two-dimensional square, as well as hon- eycomb lattices composed of 6x6 elements. Specifically, we find regions in driver voltage and frequency where stationary discrete breathers, also known as intrinsic localized modes (ILM), exist and are stable due to the interplay of damping and spatially homogeneous driving. By introduc- ing additional capacitors into the unit cell, these lattices can controllably induce traveling discrete breathers. When more than one such ILMs are experimentally generated in the lattice, the interplay of nonlinearity, discreteness and wave interactions generate a complex dynamics wherein the ILMs attempt to maintain a minimum distance between one another. Numerical simulations show good agreement with experimental results, and confirm that these phenomena qualitatively carry over to larger lattice sizes.

preprint2013arXiv

Nonlinear modes and symmetries in linearly-coupled pairs of PT-invariant dimers

The subject of the work are pairs of linearly coupled PT-symmetric dimers. Two different settings are introduced, namely, straight-coupled dimers, where each gain site is linearly coupled to one gain and one loss site, and cross-coupled dimers, with each gain site coupled to two lossy ones. The latter pair with equal coupling coefficients represents a "PT-hypersymmetric" quadrimer. We find symmetric and antisymmetric solutions in these systems, chiefly in an analytical form, and explore the existence, stability and dynamical behavior of such solutions by means of numerical methods. We thus identify bifurcations occurring in the systems, including spontaneous symmetry breaking and saddle-center bifurcations. Simulations demonstrate that evolution of unstable branches typically leads to blowup. However, in some cases unstable modes rearrange into stable ones.

preprint2013arXiv

Nonlinear PT-symmetric models bearing exact solutions

We study the nonlinear Schr$\ddot{o}$dinger equation with a PT-symmetric potential. Using a hydrodynamic formulation and connecting the phase gradient to the field amplitude, allows for a reduction of the model to a Duffing or a generalized Duffing equation. This way, we can obtain exact soliton solutions existing in the presence of suitable PT-symmetric potentials, and study their stability and dynamics. We report interesting new features, including oscillatory instabilities of solitons and (nonlinear) PT-symmetry breaking transitions, for focusing and defocusing nonlinearities.

preprint2013arXiv

PT-Symmetric Dimer of Coupled Nonlinear Oscillators

We provide a systematic analysis of a prototypical nonlinear oscillator system respecting PT-symmetry i.e., one of them has gain and the other an equal and opposite amount of loss. Starting from the linear limit of the system, we extend considerations to the nonlinear case for both soft and hard cubic nonlinearities identifying symmetric and anti-symmetric breather solutions, as well as symmetry breaking variants thereof. We propose a reduction of the system to a Schrödinger type PT-symmetric dimer, whose detailed earlier understanding can explain many of the phenomena observed herein, including the PT phase transition. Nevertheless, there are also significant parametric as well as phenomenological potential differences between the two models and we discuss where these arise and where they are most pronounced. Finally, we also provide examples of the evolution dynamics of the different states in their regimes of instability.

preprint2013arXiv

PT-symmetry Management in Oligomer Systems

We study the effects of management of the PT-symmetric part of the potential within the setting of Schrödinger dimer and trimer oligomer systems. This is done by rapidly modulating in time the gain/loss profile. This gives rise to a number of interesting properties of the system, which are explored at the level of an averaged equation approach. Remarkably, this rapid modulation provides for a controllable expansion of the region of exact PT-symmetry, depending on the strength and frequency of the imposed modulation. The resulting averaged models are analyzed theoretically and their exact stationary solutions are translated into time-periodic solutions through the averaging reduction. These are, in turn, compared with the exact periodic solutions of the full non-autonomous PT-symmetry managed problem and very good agreement is found between the two.

preprint2013arXiv

Revisiting the PT-symmetric Trimer: Bifurcations, Ghost States and Associated Dynamics

In this paper, we revisit one of the prototypical PT-symmetric oligomers, namely the trimer. We find all the relevant branches of "regular" solutions and analyze the bifurcations and instabilities thereof. Our work generalizes the formulation that was proposed recently in the case of dimers for the so-called "ghost states" of trimers, which we also identify and connect to symmetry-breaking bifurcations from the regular states. We also examine the dynamics of unstable trimers, as well as those of the ghost states in the parametric regime where the latter are found to exist. Finally, we present the current state of the art for optical experiments in PT-symmetric trimers, as well as experimental results in a gain-loss-gain three channel waveguide structure.

preprint2013arXiv

The Drastic Role of Beyond Nearest-Neighbor Interactions on Two-Dimensional Dynamical Lattices: A Case Example

In the present work, we highlight the significant effect that the simplest beyond nearest neighbor interactions can have on two-dimensional dynamical lattices. To do so, we select as our case example the closest further neighbor, namely the diagonal one, and a prototypical nonlinear lattice, the discrete nonlinear Schrodinger equation. Varying solely the strength of this extra neighbor interaction, we see examples of (a) destabilization of states that start out as stable in the nearest neighbor limit; (b) stabilization of states that start out as unstable in that limit; (c) bifurcation of novel states that do not exist in the nearest neighbor case. These dramatic changes are first theoretically highlighted through an analysis of a reduction of the problem to a few excited sites and the associated set of conditions that govern their existence and their dynamical stability. Then, they are corroborated numerically through fixed point computations, spectral analysis and nonlinear dynamical evolution simulations.

preprint2013arXiv

Traveling Waves in 2D Hexagonal Granular Crystal Lattices

We describe the dynamic response of a two-dimensional hexagonal packing of uncompressed stainless steel spheres excited by localized impulsive loadings. After the initial impact strikes the system, a characteristic wave structure emerges and continuously decays as it propagates through the lattice. Using an extension of the binary collision approximation (BCA) for one-dimensional chains, we predict its decay rate, which compares well with numerical simulations and experimental data. While the hexagonal lattice does not support constant speed traveling waves, we provide scaling relations that characterize the power law decay of the wave velocity. Lastly, we discuss the effects of weak disorder on the directional amplitude decay rates.

P. G. Kevrekidis

What is connected

Connect this record

See the researcher in context

Building this map preview

248 published item(s)

A Spectral Analysis of the Nonlinear Schroedinger Equation in the Co-Exploding Frame

An Ising machine based on networks of subharmonic electrical resonators

Backcasting COVID-19: A Physics-Informed Estimate for Early Case Incidence

Dark solitons in a trapped gas of long-range interacting bosons

Dispersive Shock Waves in Lattices: A Dimension Reduction Approach

Floquet solitons in square lattices: Existence, Stability and Dynamics

Forced symmetry breaking as a mechanism for rogue bursts in a dissipative nonlinear dynamical lattice

Mixed dispersion nonlinear Schrödinger equation in higher dimensions: theoretical analysis and numerical computations

Nonlinear Anti-(Parity-Time) symmetric dimer

Stationary multi-kinks in the discrete sine-Gordon equation

Superfluid vortex multipoles and soliton stripes on a torus

The role of mobility in the dynamics of the COVID-19 epidemic in Andalusia

Dark solitons under higher order dispersion

Dynamics of two languages competing on a network: a case study

Existence, Stability and Dynamics of Monopole and Alice Ring Solutions in Anti-Ferromagnetic Spinor Condensates

Kink-antikink stripe interactions in the two-dimensional sine-Gordon equation

Measurement and Memory in the Periodically Driven Complex Ginzburg-Landau equation

Moving discrete breathers in a $β$-FPU lattice revisited

Phase Diagram, Stability and Magnetic Properties of Nonlinear Excitations in Spinor Bose-Einstein Condensates

A quantitative framework for exploring exit strategies from the COVID-19 lockdown

Controlled Generation of Dark-Bright Soliton Complexes in Two-Component and Spinor Bose-Einstein Condensates

Kink-Antikink Interaction Forces and Bound States in a Biharmonic $ϕ^4$ Model

Kuznetsov-Ma breather-like solutions in the Salerno model

Many-body effects on second-order phase transitions in spinor Bose-Einstein condensates and breathing dynamics

Nonlinear edge modes in a honeycomb electrical lattice near the Dirac points

Nonlinearity and Topology

Observation and Analysis of Multiple Dark-Antidark Solitons in Two-Component Bose-Einstein Condensates

On-demand generation of dark soliton trains in Bose-Einstein condensates

Parametrically excited star-shaped patterns at the interface of binary Bose-Einstein condensates

Quench induced vortex-bright-soliton formation in binary Bose-Einstein condensates

Solitary waves in the resonant nonlinear Schrödinger equation: stability and dynamical properties

Stability of traveling waves in a driven Frenkel-Kontorova model

Vortex Pairs in the Discrete Nonlinear Schrödinger Equation

Bifurcation analysis of stationary solutions of two-dimensional coupled Gross-Pitaevskii equations using deflated continuation

Discrete breathers in $ϕ^4$ and related models

Existence and spectral stability of multi-pulses in discrete Hamiltonian lattice systems

Rogue Waves and Periodic Solutions of a Nonlocal Nonlinear Schrödinger Model

Two-dimensional rogue waves on zero background of the Davey-Stewartson II equation

Existence, Stability & Dynamics of Nonlinear Modes in a 2d Partially $\mathcal{PT}$ Symmetric Potential

$SO(2)$-induced breathing patterns in multi-component Bose-Einstein condensates

A PT-symmetric dual-core system with the sine-Gordon nonlinearity and derivative coupling

Dark spherical shell solitons in three-dimensional Bose-Einstein condensates: Existence, stability and dynamics

On the spontaneous time-reversal symmetry breaking in synchronously-pumped passive Kerr resonators

Rogue Waves in Ultracold Bosonic Seas

Scattering of Waves by Impurities in Precompressed Granular Chains

Stability of solitary waves and vortices in a 2D nonlinear Dirac model

The dynamical playground of a higher-order cubic Ginzburg-Landau equation: from orbital connections and limit cycles to invariant tori and the onset of chaos

The emergence and analysis of Kuramoto-Sakaguchi-like models as an effective description for the dynamics of coupled Wien-bridge oscillators

Vector Dark-Antidark Solitary Waves in Multi-Component Bose-Einstein condensates

Vortex-soliton complexes in coupled nonlinear Schrödinger equations with unequal dispersion coefficients

Bose-Einstein condensation: Twenty years after

Bright Discrete Solitons in Spatially Modulated DNLS Systems

Collapse for the higher-order nonlinear Schrödinger equation

Conical Wave Propagation and Diffraction in 2D Hexagonally Packed Granular Lattices

Conservation laws, exact travelling waves and modulation instability for an extended nonlinear Schrödinger equation

Dark solitons near potential and nonlinearity steps

Dark-bright solitons in coupled nonlinear Schrödinger equations with unequal dispersion coefficients

Discrete Solitons and Vortices in Anisotropic Hexagonal and Honeycomb Lattices

Effects of interactions on the generalized Hong-Ou-Mandel effect

Exciting and Harvesting Vibrational States in Harmonically Driven Granular Chains

Existence, Stability and Dynamics of Discrete Solitary Waves in a Binary Waveguide Array

Generating and Manipulating Quantized Vortices On-Demand in a Bose-Einstein Condensate: a Numerical Study

Multidimensional Discrete Compactons in Nonlinear Schrödinger Lattices with Strong Nonlinearity Management

Multifrequency and edge breathers in the discrete sine-Gordon system via subharmonic driving: theory, computation and experiment

Non-Conservative Variational Approximation for Nonlinear Schroedinger Equations

Nonlinear Resonances and Antiresonances of a Forced Sonic Vacuum

On the Bifurcation and Stability of Single and Multiple Vortex Rings in Three-Dimensional Bose-Einstein Condensates

Positive and negative mass solitons in spin-orbit coupled Bose-Einstein condensates

PT Meets Supersymmetry and Nonlinearity: An Analytically Tractable Case Example

Robust Vortex Lines, Vortex Rings and Hopfions in 3D Bose-Einstein Condensates

Scattering of matter-waves in spatially inhomogeneous environments

Solitons in Multi-Component Nonlinear Schrodinger Models: A Survey of Recent Developments

Some case example exact solutions for quadratically nonlinear optical media with $\mathcal{PT}$-symmetric potentials

Stabilization of ring dark solitons in Bose-Einstein condensates