Lyonell Boulton

Lyonell Boulton, Aatef Hobiny

The convergence of the so-called quadratic method for computing eigenvalue enclosures of general self-adjoint operators is examined. Explicit asymptotic bounds for convergence to isolated eigenvalues are found. These bounds turn out to improve significantly upon those determined in previous investigations. The theory is illustrated by means of several numerical experiments performed on particularly simple benchmark models of one-dimensional Schrodinger operators.

Lyonell Boulton, Gabriel Lord

For $q>12/11$ the eigenfunctions of the non-linear eigenvalue problem associated to the one-dimensional $q$-Laplacian are known to form a Riesz basis of $L^2(0,1)$. We examine in this paper the approximation properties of this family of functions and its dual, in order to establish non-orthogonal spectral methods for the $p$-Poisson boundary value problem and its corresponding parabolic time evolution initial value problem. The principal objective of our analysis is the determination of optimal values of $q$ for which the best approximation is achieved for a given $p$ problem.

Lyonell Boulton, Nabile Boussaid

We discuss a novel strategy for computing the eigenvalues and eigenfunctions of the relativistic Dirac operator with a radially symmetric potential. The virtues of this strategy lie on the fact that it avoids completely the phenomenon of spectral pollution and it always provides two-side estimates for the eigenvalues with explicit error bounds on both eigenvalues and eigenfunctions. We also discuss convergence rates of the method as well as illustrate our results with various numerical experiments.

Lyonell Boulton, Michael Strauss

Let A be a self-adjoint operator acting on a Hilbert space. The notion of second order spectrum of A relative to a given finite-dimensional subspace L has been studied recently in connection with the phenomenon of spectral pollution in the Galerkin method. We establish in this paper a general framework allowing us to determine how the second order spectrum encodes precise information about the multiplicity of the isolated eigenvalues of A. Our theoretical findings are supported by various numerical experiments on the computation of inclusions for eigenvalues of benchmark differential operators via finite element bases.

Lyonell Boulton, Michael Strauss

We discuss how to compute certified enclosures for the eigenvalues of benchmark linear magnetohydrodynamics operators in the plane slab and cylindrical pinch configurations. For the plane slab, our method relies upon the formulation of an eigenvalue problem associated to the Schur complement, leading to highly accurate upper bounds for the eigenvalue. For the cylindrical configuration, a direct application of this formulation is possible, however, it cannot be rigourously justified. Therefore in this case we rely on a specialized technique based on a method proposed by Zimmermann and Mertins. In turns this technique is also applicable for finding accurate complementary bounds in the case of the plane slab. We establish convergence rates for both approaches.

Lehel Banjai, Lyonell Boulton

We establish a strategy for finding sharp upper and lower numerical bounds of the Poincaré constant on a class of planar domains with piecewise self-similar boundary. The approach consists of four main components: W1) tight inner-outer shape interpolation, W2) conformal mapping of the approximate polygonal regions, W3) grad-div system formulation of the spectral problem and W4) computation of the eigenvalue bounds. After describing the method, justifying its validity and determining general convergence estimates, we show concrete evidence of its effectiveness by computing lower and upper bound estimates for the constant on the Koch snowflake.

Lyonell Boulton, Monika Winklmeier

A certified strategy for determining sharp intervals of enclosure for the eigenvalues of matrix differential operators with singular coefficients is examined. The strategy relies on computing the second order spectrum relative to subspaces of continuous piecewise linear functions. For smooth perturbations of the angular Kerr-Newman Dirac operator, explicit rates of convergence due to regularity of the eigenfunctions are established. Existing benchmarks are validated and sharpened by several orders of magnitude in the unperturbed setting.

Lyonell Boulton, Michael Levitin

This paper addresses the problem of computing the eigenvalues lying in the gaps of the essential spectrum of a periodic Schrodinger operator perturbed by a fast decreasing potential. We use a recently developed technique, the so called quadratic projection method, in order to achieve convergence free from spectral pollution. We describe the theoretical foundations of the method in detail, and illustrate its effectiveness by several examples.

Lyonell Boulton, Peter Lancaster

The theme of this paper was motivated by the question: How effective are path-following procedures for tracing the pseudospectral boundary? The present study of the mathematical properties of the boundary of the pseudospectrum is the result. This boundary is generally made up piecewise smooth curves. We shown how the Schur triangular form of the matrix can be used to analyse dynamical properties of the singular points on these curves.

Lyonell Boulton, Peter Lancaster, Panayiotis Psarrakos

In the first part of this paper, the main concern is with smoothness properties of the boundary of the pseudospectrum of a matrix polynomial. In the second part, results are obtained concerning the number of connected components of pseudospectra, as well as results concerning matrix polynomials with multiple eigenvalues, or the proximity to such polynomials.

Lyonell Boulton

We establish sufficiency conditions in order to achieve approximation to discrete eigenvalues of self-adjoint operators in the second-order projection method suggested recently by Levitin and Shargorodsky, [math.SP/0212087]. We find explicit estimates for the eigenvalue error and study in detail two concrete model examples. Our results show that, unlike the majority of the standard methods, second-order projection strategies combine non-pollution and approximation at a very high level of generality.