Source author record

Arieh Iserles

Arieh Iserles appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

7works
3topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

7 published item(s)

preprint2022arXiv

Approximation of wave packets on the real line

In this paper we compare three different orthogonal systems in $\mathrm{L}_2(\mathbb{R})$ which can be used in the construction of a spectral method for solving the semi-classically scaled time dependent Schrödinger equation on the real line, specifically, stretched Fourier functions, Hermite functions and Malmquist--Takenaka functions. All three have banded skew-Hermitian differentiation matrices, which greatly simplifies their implementation in a spectral method, while ensuring that the numerical solution is unitary -- this is essential in order to respect the Born interpretation in quantum mechanics and, as a byproduct, ensures numerical stability with respect to the $\mathrm{L}_2(\mathbb{R})$ norm. We derive asymptotic approximations of the coefficients for a wave packet in each of these bases, which are extremely accurate in the high frequency regime. We show that the Malmquist--Takenaka basis is superior, in a practical sense, to the more commonly used Hermite functions and stretched Fourier expansions for approximating wave packets

preprint2022arXiv

Positivity-preserving methods for population models

Many important applications are modelled by differential equations with positive solutions. However, it remains an outstanding open problem to develop numerical methods that are both (i) of a high order of accuracy and (ii) capable of preserving positivity. It is known that the two main families of numerical methods, Runge-Kutta methods and multistep methods, face an order barrier: if they preserve positivity, then they are constrained to low accuracy: they cannot be better than first order. We propose novel methods that overcome this barrier: our methods are of second order, and they are guaranteed to preserve positivity. Our methods apply to a large class of differential equations that have a special graph Laplacian structure, which we elucidate. The equations need be neither linear nor autonomous and the graph Laplacian need not be symmetric. This algebraic structure arises naturally in many important applications where positivity is required. We showcase our new methods on applications where standard high order methods fail to preserve positivity, including infectious diseases, Markov processes, master equations and chemical reactions.

preprint2022arXiv

Sobolev-Orthogonal Systems with Tridiagonal Skew-Hermitian Differentiation Matrices

We introduce and develop a theory of orthogonality with respect to Sobolev inner products on the real line for sequences of functions with a tridiagonal, skew-Hermitian differentiation matrix. While a theory of such L2-orthogonal systems is well established, Sobolev orthogonality requires new concepts and their analysis. We characterise such systems completely as appropriately weighed Fourier transforms of orthogonal polynomials and present a number of illustrative examples, inclusive of a Sobolev-orthogonal system whose leading N coefficients can be computed in $\mathcal{O}(N \log N)$ operations.

preprint2022arXiv

Solving the linear semiclassical Schrödinger equation on the real line

The numerical solution of a linear Schrödinger equation in the semiclassical regime is very well understood in a torus $\mathbb{T}^d$. A raft of modern computational methods are precise and affordable, while conserving energy and resolving high oscillations very well. This, however, is far from the case with regard to its solution in $\mathbb{R}^d$, a setting more suitable for many applications. In this paper we extend the theory of splitting methods to this end. The main idea is to derive the solution using a spectral method from a combination of solutions of the free Schrödinger equation and of linear scalar ordinary differential equations, in a symmetric Zassenhaus splitting method. This necessitates detailed analysis of certain orthonormal spectral bases on the real line and their evolution under the free Schrödinger operator.

preprint2016arXiv

Efficient methods for time-dependence in semiclassical Schrödinger equations

We build an efficient and unitary (hence stable) method for the solution of the semi-classical Schrödinger equation subject with explicitly time-dependent potentials. The method is based on a combination of the Zassenhaus decomposition (Bader, Iserles, Kropielnicka & Singh 2014) with the Magnus expansion of the time-dependent Hamiltonian. We conclude with numerical experiments.

preprint1994arXiv

Biorthogonal polynomials and zero-mapping transformations

The authors have presented in \cite{IN2} a technique to generate transformations $\cal T$ of the set ${\Bbb P}_n$ of $n$th degree polynomials to itself such that if $p\in{\Bbb P}_n$ has all its zeros in $(c,d)$ then ${\cal T}\{p\}$ has all its zeros in $(a,b)$, where $(a,b)$ and $(c,d)$ are given real intervals. The technique rests upon the derivation of an explicit form of biorthogonal polynomials whose Borel measure is strictly sign consistent and such that the ratio of consecutive generalized moments is a rational $[1/1]$ function of the parameter. Specific instances of strictly sign consistent measures that have been debated in \cite{IN2} include $x^μ\Dψ(x)$, $μ^x\Dψ(x)$ and $x^{\log_qμ}\Dψ(x)$, $q\in(0,1)$. In this paper we identify all measures $ψ$ such that their consecutive generalized moments have a rational $[1/1]$ quotient, thereby characterizing all possible zero-mapping transformations of this kind.

preprint1993arXiv

From Schrödinger spectra to orthogonal polynomials, via a functional equation

The main difference between certain spectral problems for linear Schrödinger operators, e.g. the almost Mathieu equation, and three-term recurrence relations for orthogonal polynomials is that in the former the index ranges across $\ZZ$ and in the latter only across $\Zp$. We present a technique that, by a mixture of Dirichlet and Taylor expansions, translates the almost Mathieu equation and its generalizations to three term recurrence relations. This opens up the possibility of exploiting the full power of the theory of orthogonal polynomials in the analysis of Schrödinger spectra. Aforementioned three-term recurrence relations share the property that their coefficients are almost periodic. We generalize a method of proof, due originally to Jeff Geronimo and Walter van Assche, to investigate essential support of the Borel measure of associated orthogonal polynomials, thereby deriving information on the underlying absolutely continuous spectra of Schrödinger operators.