Researcher profile

Marcus Webb

Marcus Webb contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 17 - UnverifiedVerification L1Unclaimed author
4works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

4 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.