Researcher profile

P. G. Kevrekidis

P. G. Kevrekidis contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
38works
0followers
16topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

38 published item(s)

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.