Numerical instability of the Akhmediev breather and a finite-gap model of it

In this paper we study the numerical instabilities of the NLS Akhmediev breather, the simplest space periodic, one-mode perturbation of the unstable background, limiting our considerations to the simplest case of one unstable mode. In agreement with recent theoretical findings of the authors, in the situation in which the round-off errors are negligible with respect to the perturbations due to the discrete scheme used in the numerical experiments, the split-step Fourier method (SSFM), the numerical output is well-described by a suitable genus 2 finite-gap solution of NLS. This solution can be written in terms of different elementary functions in different time regions and, ultimately, it shows an exact recurrence of rogue waves described, at each appearance, by the Akhmediev breather. We discover a remarkable empirical formula connecting the recurrence time with the number of time steps used in the SSFM and, via our recent theoretical findings, we establish that the SSFM opens up a vertical unstable gap whose length can be computed with high accuracy, and is proportional to the inverse of the square of the number of time steps used in the SSFM. This neat picture essentially changes when the round-off error is sufficiently large. Indeed experiments in standard double precision show serious instabilities in both the periods and phases of the recurrence. In contrast with it, as predicted by the theory, replacing the exact Akhmediev Cauchy datum by its first harmonic approximation, we only slightly modify the numerical output. Let us also remark, that the first rogue wave appearance is completely stable in all experiments and is in perfect agreement with the Akhmediev formula and with the theoretical prediction in terms of the Cauchy data.Comment: 27 pages, 8 figures, Formula (30) at page 11 was corrected, arXiv admin note: text overlap with arXiv:1707.0565

On scattering of solitons for the Klein-Gordon equation coupled to a particle

We establish the long time soliton asymptotics for the translation invariant nonlinear system consisting of the Klein-Gordon equation coupled to a charged relativistic particle. The coupled system has a six dimensional invariant manifold of the soliton solutions. We show that in the large time approximation any finite energy solution, with the initial state close to the solitary manifold, is a sum of a soliton and a dispersive wave which is a solution of the free Klein-Gordon equation. It is assumed that the charge density satisfies the Wiener condition which is a version of the ``Fermi Golden Rule''. The proof is based on an extension of the general strategy introduced by Soffer and Weinstein, Buslaev and Perelman, and others: symplectic projection in Hilbert space onto the solitary manifold, modulation equations for the parameters of the projection, and decay of the transversal component.Comment: 47 pages, 2 figure

The exact rogue wave recurrence in the NLS periodic setting via matched asymptotic expansions, for 1 and 2 unstable modes

The focusing Nonlinear Schr\"odinger (NLS) equation is the simplest universal model describing the modulation instability (MI) of quasi monochromatic waves in weakly nonlinear media, the main physical mechanism for the generation of rogue (anomalous) waves (RWs) in Nature. In this paper we investigate the $x$-periodic Cauchy problem for NLS for a generic periodic initial perturbation of the unstable constant background solution, in the case of $N=1,2$ unstable modes. We use matched asymptotic expansion techniques to show that the solution of this problem describes an exact deterministic alternate recurrence of linear and nonlinear stages of MI, and that the nonlinear RW stages are described by the N-breather solution of Akhmediev type, whose parameters, different at each RW appearence, are always given in terms of the initial data through elementary functions. This paper is motivated by a preceeding work of the authors in which a different approach, the finite gap method, was used to investigate periodic Cauchy problems giving rise to RW recurrence.Comment: 20 pages. arXiv admin note: text overlap with arXiv:1708.00762 and substantial text overlap with arXiv:1707.0565

XMDS2: Fast, scalable simulation of coupled stochastic partial differential equations

XMDS2 is a cross-platform, GPL-licensed, open source package for numerically integrating initial value problems that range from a single ordinary differential equation up to systems of coupled stochastic partial differential equations. The equations are described in a high-level XML-based script, and the package generates low-level optionally parallelised C++ code for the efficient solution of those equations. It combines the advantages of high-level simulations, namely fast and low-error development, with the speed, portability and scalability of hand-written code. XMDS2 is a complete redesign of the XMDS package, and features support for a much wider problem space while also producing faster code.Comment: 9 pages, 5 figure

Boundary effects on the dynamics of chains of coupled oscillators

We study the dynamics of a chain of coupled particles subjected to a restoring force (Klein-Gordon lattice) in the cases of either periodic or Dirichlet boundary conditions. Precisely, we prove that, when the initial data are of small amplitude and have long wavelength, the main part of the solution is interpolated by a solution of the nonlinear Schr\"odinger equation, which in turn has the property that its Fourier coefficients decay exponentially. The first order correction to the solution has Fourier coefficients that decay exponentially in the periodic case, but only as a power in the Dirichlet case. In particular our result allows one to explain the numerical computations of the paper \cite{BMP07}