Recently, we analysed spontaneous symmetry breaking (SSB) of solitons in linearly coupled dual-core waveguides with fractional diffraction and cubic nonlinearity. In a practical context, the system can serve as a model for optical waveguides with the fractional diffraction or Bose–Einstein condensate of particles with Lévy index $\alpha <2$. In an earlier study, the SSB in the fractional coupler was identified as the bifurcation of subcritical type, becoming extremely subcritical in the limit of $\alpha \rightarrow 1$. There, the moving solitons and collisions between them at low speeds were also explored. In the present paper, we present new numerical results for fast solitons demonstrating restoration of symmetry in post-collision dynamics.

A nonlinear evolution equation correct to fourth order is developed for gravity-capillary waves on linear shear currents in finite water depth. Therefore, this equation covers both effects of depth uniform currents and uniform vorticity. Starting from this equation, an instability analysis is then made for narrow banded uniform Stokes waves. The notable feature is that our investigation due to fourth order shows a remarkable improvement compared with the third-order one, and produces an excellent result compatible with the exact result of Longuet-Higgins. We observe that linear shear currents considerably change the modulational instability properties of capillary-gravity waves, such as the growth rate and bandwidth of instability.

We derive a nonlinear Schrödinger equation for the propagation of the three-dimensional broader bandwidth gravity-capillary waves including the effect of depth-uniform current. In this derivation, the restriction of narrow bandwidth constraint is extended, so that this equation will be more appropriate for application to a realistic sea wave spectrum. From this equation, an instability condition is obtained and then instability regions in the perturbed wavenumber space for a uniform wave train are drawn, which are in good agreement with the exact numerical results. As it turns out, the corrections to the stability properties that occur at the fourth-order term arise from an interaction between the mean flow and the frequency-dispersion term. Since the frequency-dispersion term, in the absence of depth-uniform current, for pure capillary waves is of opposite sign for pure gravity waves, so too are the corrections to the instability properties.

We consider the Cauchy problem for the nonlinear Schrödinger equation on the whole space. After introducing a weaker concept of finite speed of propagation, we show that the concatenation of initial data gives rise to solutions whose time of existence increases as one translates one of the initial data. Moreover, we show that, given global decaying solutions with initial data u0, v0, if |y| is large, then the concatenated initial data u0 + v0(· − y) gives rise to globally decaying solutions.

We are concerned with the existence of positive weak solutions, as well as the existence of bound states (i.e. solutions in W1, p (ℝN)), for quasilinear scalar field equations of the form

$$ - \Delta _pu + V(x) \vert u \vert ^{p - 2}u = K(x) \vert u \vert ^{q - 2}u + \vert u \vert ^{p^ * - 2}u,\qquad x \in {\open R}^N,$$

In this paper, we introduce the Hamiltonian boundary value method (HBVM) to solve nonlinear Hamiltonian PDEs. We use the idea of Fourier pseudospectral method in spatial direction, which leads to the finite-dimensional Hamiltonian system. The HBVM, which can preserve the Hamiltonian effectively, is applied in time direction. Then the nonlinear Schrödinger (NLS) equation and the Korteweg-de Vries (KdV) equation are taken as examples to show the validity of the proposed method. Numerical results confirm that the proposed method can simulate the propagation and collision of different solitons well. Meanwhile the corresponding errors in Hamiltonian and other intrinsic invariants are presented to show the good preservation property of the proposed method during long-time numerical calculation.

In this paper, we develop a novel energy-preserving wavelet collocation method for solving general multi-symplectic formulations of Hamiltonian PDEs. Based on the autocorrelation functions of Daubechies compactly supported scaling functions, the wavelet collocation method is conducted for spatial discretization. The obtained semi-discrete system is shown to be a finite-dimensional Hamiltonian system, which has an energy conservation law. Then, the average vector field method is used for time integration, which leads to an energy-preserving method for multi-symplectic Hamiltonian PDEs. The proposed method is illustrated by the nonlinear Schrödinger equation and the Camassa-Holm equation. Since differentiation matrix obtained by the wavelet collocation method is a cyclic matrix, we can apply Fast Fourier transform to solve equations in numerical calculation. Numerical experiments show the high accuracy, effectiveness and conservation properties of the proposed method.

We study infinite soliton trains solutions of nonlinear Schrödinger equations, i.e. solutions behaving as the sum of infinitely many solitary waves at large time. Assuming the composing solitons have sufficiently large relative speeds, we prove the existence and uniqueness of such a soliton train. We also give a new construction of multi-solitons (i.e. finite trains) and prove uniqueness in an exponentially small neighbourhood, and we consider the case of solutions composed of several solitons and kinks (i.e. solutions with a non-zero background at infinity).

This paper studies a local discontinuous Galerkin method combined with fourth order exponential time differencing Runge-Kutta time discretization and a fourth order conservative method for solving the nonlinear Schrödinger equations. Based on different choices of numerical fluxes, we propose both energy-conserving and energy-dissipative local discontinuous Galerkin methods, and have proven the error estimates for the semi-discrete methods applied to linear Schrödinger equation. The numerical methods are proven to be highly efficient and stable for long-range soliton computations. Extensive numerical examples are provided to illustrate the accuracy, efficiency and reliability of the proposed methods.

This paper is concerned with the Cauchy problem for a nonlinear Schrödinger equation with a harmonic potential and exponential growth nonlinearity in two space dimensions. In the defocusing case, global well-posedness is obtained. In the focusing case, existence of nonglobal solutions is discussed via potential-well arguments.

In this paper, we are concerned with the existence of multi-bump solutions for anonlinear Schrödinger equations with electromagnetic fields. We prove under some suitableconditions that for any positive integer m, there existsε(m) > 0 such that, for0 < ε < ε(m),the problem has an m-bump complex-valued solution. As a result, whenε → 0, the equation has more and more multi-bumpcomplex-valued solutions.

We study numerically the semiclassical limit for the nonlinearSchrödinger equation thanks to a modification of the Madelungtransform due to Grenier. This approach allows for the presence ofvacuum. Even if the mesh size and the time step do not depend on the Planck constant, we recover the position and current densities in thesemiclassical limit, with a numerical rate of convergence inaccordance with the theoreticalresults, before shocks appear in the limiting Eulerequation. By using simple projections, the mass and the momentum ofthe solution are well preserved by the numerical scheme,while the variation of the energy is not negligiblenumerically. Experiments suggest that beyond the critical time for theEuler equation, Grenier's approach yields smooth but highlyoscillatory terms.

We consider the existence and multiplicity of standing-wave solutions

of nonlinear Schrödinger equations with electromagnetic fields and critical nonlinearity

Under suitable assumptions, we prove that it has at least one solution and that, for any m ∈ ℕ, it has at least m pairs of solutions.

The plane wave stability properties of the conservative schemes of Besse [SIAM J. Numer. Anal.42 (2004) 934–952] and Fei et al. [Appl. Math. Comput.71 (1995) 165–177] for the cubic Schrödinger equation are analysed. Although the two methods possess many of the same conservation properties, we show that their stability behaviour is very different.An energy preserving generalisation of the Fei method with improved stability is presented.

We prove global internal controllability in large time for the nonlinear Schrödinger equation on a bounded interval with periodic, Dirichlet or Neumann conditions. Our strategy combines stabilization and local controllability near 0. We use Bourgain spaces to prove this result on L2 . We also get a regularity result about the control if the data are assumed smoother.

We consider the focusing nonlinear Schrödinger equation on the quarter plane. Initial data vanish at infinity while boundary data are time-periodic. The main goal of this paper is to introduce a Riemann–Hilbert problem whose solution gives the solution of our initial–boundary-value problem. This is a preliminary step to obtain uniform long-time asymptotics for the solution of this equation using both the stationary phase and the Deift–Zhou methods.

We discretize the nonlinear Schrödinger equation,with Dirichlet boundary conditions, by a linearlyimplicit two-step finite element method which conserves the L 2 norm. We prove optimal order a priori error estimates in the L 2 and H 1 norms, under mild mesh conditions for two and three space dimensions.

The aim of this work is to establish, from amathematical point of view, the limit α → +∞ in the system $i \partial_t E+\nabla (\nabla . E)-\alpha^2 \nabla \times\nabla \times E =-|E|^{2\sigma}E,$ where $E:{\ensuremath{{\Bbb R}}}^3\rightarrow{\mathbb C}^3$ . This corresponds to an approximationwhich is made in the context of Langmuir turbulence in plasmaPhysics. The L 2-subcritical σ (that is σ ≤ 2/3)and the H 1-subcritical σ (that is σ ≤ 2) arestudied. In the physical case σ = 1, the limit is then studied for the $H^1({\ensuremath{{\Bbb R}}}^3)$ norm.