Source author record

Marcus Webb

Marcus Webb 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
5topics
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)

preprint2026arXiv

Simulation of Radiation Chemistry by a One-Shot Hybrid Continuum / Monte Carlo Method

Understanding the spatio-temporal evolution of radiolytic species created by high-energy electrons in water underpins key applications from radiotherapy and nuclear safety to environmental processing and electron microscopy. Here, using the Manchester Inhomogeneous Radiation Chemistry by Linear Expansions (MIRaCLE) toolkit, we introduce and benchmark a novel approach to simulating these processes. Although the initial conditions are determined stochastically, the subsequent time evolution is calculated deterministically using a continuum representation, derived from those initial conditions. This hybrid approach essentially averages over many chemistry ``trajectories'' simultaneously, often converging to the 1% level in one shot, not requiring multiple runs. We demonstrate this new approach through the calculation of time-dependent G-values for e_{aq}^-$, \dot{\mathrm{OH}} and other radiolytic products, including at unprecedented dose rates where calculations which would take years with a conventional Monte Carlo approach can be performed in mere hours on a commercial laptop. We demonstrate that the main artifact of continuum modelling can be mitigated by a correction term. These results establish MIRaCLE as a flexible and efficient platform for modelling long-timescale radiolysis, providing a bridge between Monte Carlo approaches and macroscopic reaction--diffusion schemes, with broad implications for radiation chemistry in medicine, energy, and materials science.

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

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

Fast polynomial transforms based on Toeplitz and Hankel matrices

Many standard conversion matrices between coefficients in classical orthogonal polynomial expansions can be decomposed using diagonally-scaled Hadamard products involving Toeplitz and Hankel matrices. This allows us to derive $\smash{\mathcal{O}(N(\log N)^2)}$ algorithms, based on the fast Fourier transform, for converting coefficients of a degree $N$ polynomial in one polynomial basis to coefficients in another. Numerical results show that this approach is competitive with state-of-the-art techniques, requires no precomputational cost, can be implemented in a handful of lines of code, and is easily adapted to extended precision arithmetic.

preprint2015arXiv

The Acyclicity of the Frobenius Functor for Modules of Finite Flat Dimension

Let $R$ be a commutative Noetherian local ring of prime characteristic $p$ and $f:R\to R$ the Frobenius ring homomorphism. For $e\ge 1$ let $R^{(e)}$ denote the ring $R$ viewed as an $R$-module via $f^e$. Results of Peskine, Szpiro, and Herzog state that for finitely generated modules $M$, $M$ has finite projective dimension if and only if $\operatorname{Tor}_i^R(R^{(e)},M)=0$ for all $i>0$ and all (equivalently, infinitely many) $e\ge 1$. We prove this statement holds for arbitrary modules using the theory of flat covers and minimal flat resolutions.

preprint2015arXiv

Volume Preservation by Runge-Kutta Methods

It is a classical theorem of Liouville that Hamiltonian systems preserve volume in phase space. Any symplectic Runge-Kutta method will respect this property for such systems, but it has been shown that no B-Series method can be volume preserving for all volume preserving vector fields (BIT 47 (2007) 351-378 and IMA J. Numer. Anal. 27 (2007) 381-405). In this paper we show that despite this result, symplectic Runge-Kutta methods can be volume preserving for a much larger class of vector fields than Hamiltonian systems, and discuss how some Runge-Kutta methods can preserve a modified measure exactly.