A new paradigm for spectral localisation of operator pencils and analytic operator-valued functions

Lead Research Organisation: CARDIFF UNIVERSITY
Department Name: Sch of Mathematics

Abstract

The natural frequencies of vibration of an elastic structure (bridge; airframe; bicycle frame) determine its mechanical behaviour. Frequencies which are likely to be excited by external loads should be designed out, in order to avoid mechanical failures. The yellow light emitted by sodium street-lamps (now mostly replaced by LCDs) was also an example of a natural frequency, determined by the law E = h f, where h is Planck's constant and E is determined by finding the allowed energy levels of electrons in a sodium atom: effectively, by solving something called an eigenvalue problem, for the Schroedinger equation.

The importance of eigenvalue problems in so many different applications has resulted in large amounts of public and industrial money being spent on the development of appropriate software, such as ARPACK and LAPACK. These deal with problems in which the eigen-parameter appears linearly. Similarly, huge efforts have been made to develop theoretical tools to understand, at the level of pure mathematics, where eigenvalues lie. Are they real or complex? Do they lie inside the unit disc, or outside? Do they have positive or negative real part? The numerical range is one of the most widely used tools for making such estimates.

However, eigen-parameters often appear polynomially or rationally. When this happens, a common approach has been to apply clever transformations to linearise the problem. In this proposal, we intend to show that this is generally not the best approach, by proving theorems and constructing estimates which treat the problem in its original form and get better estimates. In fact, we even propose that in most cases one should do the opposite of the usual approach: transform the original linear problem into a family of rational problems and, by considering estimates for the whole family together, obtain better estimates than can be obtained by direct treatment of the original. In fact, we shall show that this approach can, in theory, yield *all* information about the eigenvalues and, more generally, the spectrum, of the original problem.

Planned Impact

The pathways from this project lead to two categories of impact: Knowledge and People, and Economy and Society.

Knowledge and People

We aim to train a postdoc with skills at the interface between analysis, applied mathematics and computational mathematics, opening for them a wide panorama
of career options by equipping them with skills which are valued both in industry and academia. Our project will link areas in which the UK is already an international leader, such as computational mathematics and numerical linear algebra, with areas in which we are internationally competitive and in which IRM 2004 and IRM 2010 recommended investment, including mathematical analysis. EPSRC has made what are arguably its most strategically significant investments in
doctoral education at this interface, funding CCA in Cambridge and the PDE CTD in Oxford. Funding our proposal would target the next level up, postdoctoral training. Many of the skills which would be acquired by the PDRA would also be applicable in other rapidly growing areas such as inverse problems and imaging, particularly in view of the fact that the PI has held continuous EPSRC funding on inverse-problem projects since January 2013.

Economy and Society

Ideas from spectral analysis of large matrices are now ubiquitous in daily life and data analysis. In Cardiff, our colleague Bertrand Gauthier of the
Statistics Group is working on a new technique for low-dimensional approximation of spectra of operators in reproducing kernel Hilbert spaces, which has shown great promise in the development of techniques in machine learning. Gauthier is actively discussing this work with the PI and other members of the PI's group. The
impact of these techniques will be felt by society and by our economy, even if those who experience the impact will be unaware of its origin. Cardiff School of
Mathematics also provides an outstanding environment in which to conduct this work as we have invested in new lectureships, both in mathematics and computer
science, to underpin these sorts of collaborations with a critical mass of resident expertise. As Director of Research for the School, the PI has had a key role in setting these strategic directions and choosing new appointees.

An important part of the fourth work package in this proposal will be to seek out the new applications such as these and provide the tools which allow techniques
initially developed by practitioners to be widened in their applicability. One enticing possible outcome here is a proper understanding of just why certain techniques in machine learning function well in certain situations and not in others. If achieved this would be comparable to the work of Hansen et al. in Cambridge
on compressed sensing, which provided the first properly rigorous analysis explaining the hidden limitations behind methods which, a priori, seemed too good to be
true.

Publications

10 25 50
publication icon
Ben-Artzi J (2023) On the complexity of the inverse Sturm-Liouville problem in Pure and Applied Analysis

publication icon
Ben-Artzi J (2021) Computing the Sound of the Sea in a Seashell in Foundations of Computational Mathematics

publication icon
Ben-Artzi J (2022) Universal algorithms for computing spectra of periodic operators in Numerische Mathematik

publication icon
Ben-Artzi J (2022) Computing scattering resonances in Journal of the European Mathematical Society

publication icon
Bögli S (2023) Spectral analysis and domain truncation for Maxwell's equations in Journal de Mathématiques Pures et Appliquées

publication icon
Ferraresso F (2022) On the spectral instability for weak intermediate triharmonic problems in Mathematical Methods in the Applied Sciences

publication icon
Ferraresso F (2023) On the eigenvalues of the biharmonic operator with Neumann boundary conditions on a thin set in Bulletin of the London Mathematical Society

 
Description Maxwell's equations govern the behaviour of electromagnetic radiation described by vector-valued functions called `electric field' and `magnetic field'. Just as sound engineers analyse sound waves by decomposing them into the parts which oscillate at different frequencies (Fourier analysis),
spectral theorists analyse Maxwell equations by finding special resonant electromagnetic frequencies which are characteristic of a given material. The most interesting materials are heterogeneous and anisotropic, composed of conductors, semi-conductors and insulators, all characterised by different values of the permeability, permittivity and conductivity in the Maxwell equations. As a result of our work on `essential spectra' of Maxwell systems we now understand much bettter the presence of short-duration spatially trapped waves near material interfaces, and we have mathematically rigorous proofs that the presence of such phenomena does not contribute to the way that waves radiate to infinity. This allows us to be certain, when performing numerical experiments, that the results will be free of spectral pollution in a much wider region of spectral space than could be guaranteed by previous theorems. Indeed, we ourselves were surprised by the results: three years ago, we would have conjectured that they should be false. The results also hold for materials in which the electrical properties depend upon the frequencies of the waves to which they are exposed (described by Drude-Lorentz models).

At the same time that we did this work, we also worked on questions of `computability' for waves in Helmholtz resonators and in periodic media: what can realistically be simulated reliably on a computer, and what would require a superposition of infinite processes in order to generate a result?
The answers to these questions depend partly upon understanding sources of `spectral pollution', but relied in essential ways on the use of operator-valued analytic functions and their Schatten-class properties, as we had conjectured in the grant application.

We are now in a place where we should be able to study computability problems for Maxwell equations and indeed for other vector-valued systems.

Our work on computing the sound of seashells was mentioned in the popular science press:

https://au.news.yahoo.com/curious-kids-why-seashells-sound-003934597.html?guccounter=1
Exploitation Route Elliptic systems with vector-valued solutions are now the subject of significant effort, particularly in the area of microlocal analysis, in several groups in the EU, as well as in UCL and Heriot-Watt. Our work should be of interest to them. Our descriptions of interface modes and other mechanisms for generation of essential spectrum should also give ways of understanding EM cloaking in the presence of material interfaces.
Sectors Education