The exponentially convergent trapezoidal rule

It is well known that the trapezoidal rule converges geometrically when applied to analytic functions on periodic intervals or the real line. The mathematics and history of this phenomenon are reviewed and it is shown that far from being a curiosity, it is linked with computational methods all across scientific computing, including algorithms related to inverse Laplace transforms, special functions, complex analysis, rational approximation, integral equations, and the computation of functions and eigenvalues of matrices and operators

Convergence of a Boundary Integral Method for Water Waves

We prove nonlinear stability and convergence of certain boundary integral methods for time-dependent water waves in a two-dimensional, inviscid, irrotational, incompressible fluid, with or without surface tension. The methods are convergent as long as the underlying solution remains fairly regular (and a sign condition holds in the case without surface tension). Thus, numerical instabilities are ruled out even in a fully nonlinear regime. The analysis is based on delicate energy estimates, following a framework previously developed in the continuous case [Beale, Hou, and Lowengrub, Comm. Pure Appl. Math., 46 (1993), pp. 1269–1301]. No analyticity assumption is made for the physical solution. Our study indicates that the numerical methods must satisfy certain compatibility conditions in order to be stable. Violation of these conditions will lead to numerical instabilities. A breaking wave is calculated as an illustration

Representations of Coherent and Squeezed States in a ff-deformed Fock Space

We establish some of the properties of the states interpolating between number and coherent states denoted by $| n >_{\lambda}$; among them are the reproducing of these states by the action of an operator-valued function on $| n>$ (the standard Fock space) and the fact that they can be regarded as $f$-deformed coherent bound states. In this paper we use them, as the basis of our new Fock space which in this case are not orthogonal but normalized. Then by some special superposition of them we obtain new representations for coherent and squeezed states in the new basis. Finally the statistical properties of these states are studied in detail.Comment: 13 pages, 4 Figure

Fredholm determinants for the stability of travelling waves

This thesis investigates both theoretically and numerically the stability of travelling wave solutions using Fredholm determinants, on the real line. We identify a class of travelling wave problems for which the corresponding integral operators are of trace class. Based on the geometrical interpretation of the Evans function, we give an alternative proof connecting it to (modified) Fredholm determinants. We then extend that connection to the case of front waves by constructing an appropriate integral operator. In the context of numerical evaluation of Fredholm determinants, we prove the uniform convergence associated with the modified/regularised Fredholm determinants which generalises Bornemann's result on this topic. Unlike in Bornemann's result, we do not assume continuity but only integrability with respect to the second argument of the kernel functions. In support to our theory, we present some numerical results. We show how to compute higher order determinants numerically, in particular for integral operators belonging to classes I3 and I4 of the Schatten-von Neumann set. Finally, we numerically compute Fredholm determinants for some travelling wave problems e.g. the `good' Boussinesq equation and the fth-order KdV equation.UK EPSRC (Engineering and Physical Sciences Research Council) grant EP/G03613

Numerical verification of a gap condition for a linearized nonlinear Schrödinger equation

We make a detailed numerical study of the spectrum of two Schrödinger operators L± arising from the linearization of the supercritical nonlinear Schrödinger equation (NLS) about the standing wave, in three dimensions. This study was motivated by a recent result of the second author on the conditional asymptotic stability of solitary waves in the case of a cubic nonlinearity. Underlying the validity of this result is a spectral condition on the operators L±, namely that they have no eigenvalues nor resonances in the gap (a region of the positive real axis between zero and the continuous spectrum), which we call the gap property. The present numerical study verifies this spectral condition and shows further that the gap property holds for NLS exponents of the form 2 β + 1, as long as β* < β ≤ 1, where $$\begin{equation*}\beta_{\ast} = 0.913\,958\,905 \pm 1e-8.\end{equation*}$$ Our strategy consists of rewriting the original eigenvalue problem via the Birman–Schwinger method. From a numerical analysis viewpoint, our main contribution is an efficient quadrature rule for the kernel 1/|x - y| in {\mathbb R}^3 , i.e. proved spectrally accurate. As a result, we are able to give similar accuracy estimates for all our eigenvalue computations. We also propose an improvement in Petviashvili's iteration for the computation of standing wave profiles which automatically chooses the radial solution. All our numerical experiments are reproducible. The Matlab code can be downloaded from http://www.acm.caltech.edu/~demanet/NLS/