EPSRC Reference: 
EP/I00761X/1 
Title: 
Robust Eigenvalue Computation 
Principal Investigator: 
Boulton, Professor L 
Other Investigators: 

Researcher CoInvestigators: 

Project Partners: 

Department: 
S of Mathematical and Computer Sciences 
Organisation: 
HeriotWatt University 
Scheme: 
First Grant  Revised 2009 
Starts: 
01 February 2011 
Ends: 
16 June 2013 
Value (£): 
101,427

EPSRC Research Topic Classifications: 
Mathematical Analysis 
Numerical Analysis 

EPSRC Industrial Sector Classifications: 
No relevance to Underpinning Sectors 


Related Grants: 

Panel History: 
Panel Date  Panel Name  Outcome 
09 Sep 2010

Mathematics Prioritisation Panel

Announced


Summary on Grant Application Form 
The eigenvalues of linear operators manifest themselves almost everywhere in our everyday life. The colour of light we see is related to eigenvalues of models describing atoms and molecules. The tones and overtones of musical instruments we hear are determined by the eigenvalues of mathematical models describing strings and drums. The resonances produced by cars running over bridges are predicted by the eigenvalue analysis of beams and suspension cables. This proposal consists of four parts all connected by the central topic of computation of eigenvalues of linear operators. The eigenvalues of linear operators can be found analytically only for a few simple models, hence the importance of procedures for approximating them in practical applications. The projection method is by far the most successful robust tool for computing onesided bounds of eigenvalues of socalled selfadjoint operators. It is based on the RayleighRitz variational principle whose roots can be traced back over a century. This principle is robust in the sense that only minimal a priori information about the mathematical problem is required in order to obtain trustworthy information about the eigenvalues. In 1928 Kellner and Hylleraas used this approach to compute bounds on the ionisation energy of the helium atom. In the early days of quantum theory, the agreement of about 0.01% (now regarded as crude ) between the model calculations and the experimentally measured ionisation energy, represented important evidence supporting the validity of quantum mechanics. Remarkably, all these calculations were made long before the invention of the first digital computers. A realisation of the projection method via the finite element method underlies most of the modern computational packages, commercial or otherwise.The problem of obtaining robust bounds on eigenvalues complementary to those provided by the RayleighRitz principle, hence producing an enclosure, was begun with the work of Temple, Lehmann and Kato. This problem has a long and complicated history, and it belongs to a very active area of mathematical research. The problem is closely connected with a remarkable numerical phenomenon known to specialists as spectral pollution. The later has received a substantial amount of attention in the last 15 years. The main goal of the first part of this proposal is to investigate robust procedures for calculating intervals of enclosure for eigenvalues of selfadjoint operators. We will mainly focus in two methods which have recently been identified as successful for avoiding spectral pollution: the second order method and the Fnmethod. A most successful strategy for eigenvalue computation is to combine different numerical tools. The second part of the proposal aims at applying a combined strategy for calculating the critical magnetic field strength required for spontaneous electronpositron pair creation in a model of hydrogenic atoms. These are atoms with only one electron. It has been conjectured that this field strength is huge and can only be encountered in extreme objects such as magnetars (neutron stars with an intense magnetic field).When an operator is nonselfadjoint, the RayleighRitz principle cannot be applied and the numerical estimation of its eigenvalues is usually a highly nontrivial task. The third and fourth part of the project address the problem of computation of eigenvalues in this regime. In the third parts, we will consider extensions of a wellknown theorem by H. Weyl to the projection method. In the fourth part we will investigate a geometrical approach to eigenvalue computation which can be regarded as an extension of the Fnmethod.

Key Findings 
This information can now be found on Gateway to Research (GtR) http://gtr.rcuk.ac.uk

Potential use in nonacademic contexts 
This information can now be found on Gateway to Research (GtR) http://gtr.rcuk.ac.uk

Impacts 
Description 
This information can now be found on Gateway to Research (GtR) http://gtr.rcuk.ac.uk 
Summary 

Date Materialised 


Sectors submitted by the Researcher 
This information can now be found on Gateway to Research (GtR) http://gtr.rcuk.ac.uk

Project URL: 

Further Information: 

Organisation Website: 
http://www.hw.ac.uk 