-
Dynamic inverse problems, such as those arising in dynamic MRI, often involve reconstructing time-evolving processes from highly under-sampled measurements. Classical approaches often rely on infimal convolution regularisation which combines multiple priors, such as low-rank and sparsity, to separate static and dynamic components. While effective, these methods can be limited in their flexibility and fail to fully exploit complex data structures. In this talk, we explore how Plug-and-Play (PnP) techniques—where learned denoisers act as implicit priors—can be integrated into the infimal convolution framework. By using multiple denoisers, each tailored to different aspects of the reconstruction, we aim to enhance both adaptability and reconstruction quality in dynamic imaging settings.
-
Ergodic theory concerns the long-term behaviour of dynamical systems, particularly through the study of existence and uniqueness of invariant measures. Things get a little tricky in the case of systems with memory where the Markov property fails e.g. fractional Brownian motion and in the case of infinite dimensions e.g. SPDEs, but there are existing frameworks.
In this talk we will present ongoing work on combining these two cases to study the ergodicity of non-Markovian systems in infinite dimensions. A key motivation is the 2-D Stochastic Navier-Stokes equation with degenerate fractional noise.
P.S. I will start with an overview of measure theory!
-
Interacting particle systems are often exact solutions to singular SPDEs. Despite their singularity, numerical approximations can be used to regularise them. The resulting equations can then be used to estimate fluctuations of the underlying particle system, often with a much lower computational cost (in specific scaling regimes). In this talk, we will discuss these aspects for the Dean-Kawasaki equation and present preliminary work on discontinuous Galerkin methods for this particular SPDE.
-
The spirochetal bacterium Borrelia burgdorferi is a tick-borne zoonosis that circulates in various wildlife populations in temperate rural regions of Europe, North America and Asia. Humans are not usually competent for transmission, but spillover infections can lead to Lyme disease (LD). The infection is passed to human hosts via the bite of an infected tick. Ticks have multiple life stages and complex phenology (seasonal pattern). Over the last decade, there has been a sustained increase in Borrelia prevalence in wildlife, leading to an increase in spillover events, often via residential areas that back onto woodland. Understanding tick ecology is essential for predicting the spread of Lyme disease, informing control strategies, and assessing impacts of environmental change.
In this talk, I will discuss the development of a tick population model in a fragmented peridomestic environment. We will consider a metapopulation framework of residential patches, where humans might encounter ticks. Our principal goal is to understand the impact of deer movement on the tick ecological dynamics. Deer are the primary host of adult ticks and a necessary component of tick reproduction. They visit the residential patches in very small numbers (1 or 2 per hectare/ patch) and can disperse over large distances, transporting any feeding ticks with them. Consequently, the deer population and any events that concern them are inherently stochastic. I will discuss how we can incorporate this feature into the model whilst maintaining computational efficiency using a hybrid modelling approach. Time permitting, I will discuss the use of acaricide treatment in deer feeders as a control measure against ticks.
-
Conditional extreme value models have proven useful for analysing the joint tail behaviour of random vectors. Conditional extreme value models describe the distribution of components of a random vector conditional on at least one exceeding a suitably high threshold, and they can flexibly capture a variety of structures in the distribution tails. One drawback of these methods is that model estimates tend to be highly uncertain due to the natural scarcity of extreme data. This motivates the development of clustering methods for this class of models; pooling similar within-cluster data drastically reduces parameter estimation uncertainty.
While an extensive amount of work to estimate conditional extremes models exists in multivariate and spatial applications, the prospect of clustering for models of this type has not yet been explored. As a motivating example, we explore tail dependence of meteorological variables across multiple spatial locations and seek to identify sites which exhibit similar multivariate tail behaviour. To this end, we introduce a clustering framework for conditional extremes models which provides a novel and principled, parametric methodology for summarising multivariate extremal dependence.
In a first step, we define a dissimilarity measure for conditional extremes models based on the Jensen-Shannon divergence and common working assumptions made when fitting these models. One key advantage of our measure is that it can be applied in arbitrary dimension and, as opposed to existing methods for clustering extremal dependence, is not restricted to a bivariate setting. Clustering is then performed by applying the k-medoids algorithm to our novel dissimilarity matrix, which collects the dissimilarity between all pairs of spatial sites.
A detailed simulation study shows our technique to be superior to the leading competitor in the bivariate case across a range of possible dependence structures and uniquely provides a tool for clustering in the multivariate extremal dependence setting. We also outline a methodology for selecting the number of clusters to use in a given application. Finally, we apply our clustering framework to Irish meteorological data.
-
In this talk, we consider the 3D Euler equations in vorticity formulation for an incompressible, homogeneous fluid. Our focus is on constructing solutions that exhibit a high concentration of vorticity along a time-evolving curve, known as vortex filaments. We firstly apply the so-called Inner-Outer Gluing Scheme to construct a compactly supported helical filament, associated with a rotating-translating helix evolving by its binormal curvature flow. We then use this construction to explore the intriguing Leapfrogging Phenomenon of Vortex Helices, which describes the dynamic interaction of multiple vortex helices sharing a common axis of symmetry (e.g interaction of tornadoes). Numerical simulations available online that showcase this phenomenon will also be presented.
-
We consider a rooted tree where each vertex is labelled by an independent and identically distributed (i.i.d.) uniform(0,1) random variable, plus a parameter theta times its distance from the root. We study paths from the root to infinity along which the vertex labels are increasing. The existence of such increasing paths depends on both the structure of the tree and the value of theta. The goal is to determine the critical value of theta such that, above this value, increasing paths occur with positive probability, while below it, no such paths exist. Additionally, we extend this problem to consider the case of the integer lattice.
-
The (linear) Boltzmann transport equation plays a fundamental role in mathematical physics, with applications ranging from cancer treatment to nuclear reactor design. A key challenge in solving this equation numerically arises from the integral term over an independent angular variable, which introduces global coupling and leads to large linear systems. On top of this complexity, classical iterative solvers can fail to converge in certain regimes due to specific relationships between the equation's coefficients. In this talk, we'll explore one such challenging regime and discuss how a diffusive approximation can be used as a preconditioner to improve convergence toward a numerical solution. Expect a mix of (possibly butchered) functional analysis, hand-waving numerics, and the occasional "you'll just have to take my word for it".
-
Hypocoercivity is a property of certain dissipative evolution equations that implies exponential convergence of solutions to the steady state. We apply this theory to a kinetic toy model for heat transfer in a gas which approximates the linearised Boltzmann equation. We prove the model is hypocoercive in L2 after splitting the operator appropriately which enables us to use a recipe framework provided by Dolbeaut, Mouhot and Schmeiser.
-
Inverse problems arise in various applications involving particle transport, such as photoacoustic imaging (PAI) and proton therapy. These problems often involve inferring unknown coefficients of a PDE governing the transport process. One example is the Boltzmann equation for photon transport, where the goal is to recover optical tissue properties from observed data.
We formulate this inverse problem in a Bayesian framework and seek to infer the unknown PDE coefficients. A common choice is to use Monte Carlo simulations to obtain the PDE solution. Three key challenges need to be tackled: First, the high dimensionality of the model problem, second, the substantial computational cost, and third, stochastic nature of forward solves. We discuss several methods that can be used to address this problem.
To reduce computational cost of solving the inverse problem while maintaining accuracy, we propose a framework in which a computationally inexpensive, approximate model is corrected using an accurate forward solver. The approximate model may take various forms, such as a neural network surrogate, a diffusion approximation, or another reduced-complexity representation of particle transport. To achieve this we embed this methodology in the framework of Delayed Acceptance MCMC.
-
As in many areas of industry, there is significant interest in how mathematical modelling techniques may be applied at various stages of chocolate production, to enhance yields and improve efficiency. In my brief academic career to date, I've been fortunate to work on not one but two separate research projects in this (somewhat niche) field. My talk will begin with a discussion of the use of delay differential equations in modelling cocoa agriculture, with a particular focus on capturing the effects of varying rainfall within these models. Then, skipping over a large number of processing steps (which naturally I hope to one day return to mathematically), we will explore the fluid dynamics of molten chocolate during the enrobing process, which is used to apply a chocolate coating to e.g. biscuits. Expect extensional flows, non-Newtonian behaviours and some light asymptotics, (hopefully!) presented in an applied and accessible fashion.
-
In this PSS, I will introduce you to the card game of Poker, focusing on the no-limit Texas Hold'em variety. I will also provide a brief overview of basic poker strategy, supported by some mathematical insights. Finally, we’ll run a quick poker tournament to put this theory into practice. The goal is to conduct the tournament during the full 2-hour PSS session, but if you need to leave partway through, that’s completely fine.
N.B. The tournament is free to enter, and the prize is chocolates. So, feel free to come along and play. But most importantly, remember: you’ve got to know when to hold ’em, know when to fold ’em, know when to walk away, and know when to run.
-
Splitting methods are widely used for solving initial value problems (IVPs) due to their ability to simplify complicated evolutions into more manageable subproblems which can be solved efficiently and accurately. Traditionally, these methods are derived using analytic and algebraic techniques from numerical analysis, including truncated Taylor series and their Lie algebraic analogue, the Baker-Campbell-Hausdorff formula. These tools enable the development of high-order numerical methods that provide exceptional accuracy for small timesteps. Moreover, these methods often (nearly) conserve important physical invariants, such as mass, unitarity, and energy. However, in many practical applications the computational resources are limited. Thus, it is crucial to identify methods that achieve the best accuracy within a fixed computational budget, which might require taking relatively large timesteps. In this regime, high-order methods derived with traditional methods often exhibit large errors since they are only designed to be asymptotically optimal. Machine Learning techniques offer a potential solution since they can be trained to efficiently solve a given IVP with less computational resources. However, they are often purely data-driven, come with limited convergence guarantees in the small-timestep regime and do not necessarily conserve physical invariants. In this work, we propose a framework for finding machine learned splitting methods that are computationally efficient for large timesteps and have provable convergence and conservation guarantees in the small-timestep limit. We demonstrate numerically that the learned methods, which by construction converge quadratically in the timestep size, can be significantly more efficient than established methods for the Schrödinger equation if the computational budget is limited.
-
In this talk, we will discuss the regularity theory of the Navier-Stokes equations connected to the millennium prize problem: it is not known whether solutions to the 3D incompressible Navier-Stokes equations starting from "nice" initial data stay "nice" forever. Supposing that there is a solution that does not, certain norms/quantities of the solution approach infinity as the solution loses its smoothness. I will introduce you to a certain new, interesting, local quantity that does this, talk about some methods in the proof, and explain why we care about it!
-
The majority of research in changepoint detection has a focus on a change in the mean of data. However, for financial data for example, it is often necessary to detect changes in variance and in the tail distribution. This talk will cover some underlying extreme value theory and will apply these methods to changepoint detection algorithms with an application in finance.
-
In this talk, we'll embark on a journey through the world of rough path theory—a framework for understanding complex signals driven by irregular, noisy inputs. At the heart of our exploration is the inverse problem of noise recovery: can we reconstruct the hidden noise driving a system, using only observations of the resulting path?
There's something for everyone on this journey: numerical analysis for the approximation aficionados, probability for the lovers of randomness, and a dash of algebra for those curious about compact ways to represent these rough paths. Join us as we navigate the noise to uncover the hidden forces guiding rough dynamics!
-
In 1637, Fermat posed a seemingly simple problem in the margin of his copy of arithmetica, and added the following note: "I have discovered a truly marvellous proof of this, which this margin is too narrow to contain", a sentence that would haunt (Halloween pun intended) mathematicians for centuries.
For 358 years an army of number theorists worked together desperately trying to prove this infamous problem. A proof was first announced in 1994 (then corrected and thus completed in 1995). In this talk we give an overview of how this proof was constructed over the years for this elusive problem. In particular we study a very special class of curves named elliptic curves, which have uses in all sorts of areas of pure mathematics but are also known for their use in cryptography.
As nearly all the number theory department is absent, I have written this talk for non-number theorists and so I will not assume much number theory specific knowledge.
-
Underdamped Langevin dynamics (ULD) is a stochastic differential equation of great interest to those in the molecular dynamics and machine learning communities. Interest from the latter is due to ergodic properties of ULD, meaning that under assumptions on the potential function, the process can be used as proxy for generating samples from an unnormalized (log-concave) target density. Due to the presence of a non-linear term, one cannot simulate ULD exactly and must apply a numerical discretisation.
In this talk, we introduce the numerical method “QUICSORT”, which is a third order convergent method for simulating ULD. Our method is based off discretising the “Shifted ODE” method of Foster, Lyons and Oberhauser, where Brownian motion is replaced by a piecewise linear path matching higher order stochastic integrals. During the talk we will motivate the construction of QUICSORT, discuss the key ideas in the proof of third order convergence (in the 2-Wasserstein metric) and highlight the key role played by numerical contractivity. To the best of our knowledge, this is the first method to achieve third order convergence over the infinite-time horizon for strongly convex MCMC problems. We conclude with a numerical experiment (courtesy of Andraž Jelinčič) illustrating the performance of our method against the current state of the art sampler for high dimensional problems.
Joint work with James Foster.
-
In this talk, I go over my passion project on the surprisingly non-trivial strategy of the kid's card game Top Trumps, featuring very accessible levels of probability, analysis and programming. The main focus is on the Lord of the Rings collectors edition decks and I will have a thing or two to say about some of the questionable choices they have made!
-
With an ulterior motive of trying to get into the flow ahead of an upcoming confirmation viva, I will be talking about my work from the first year of my PhD. I am looking for steady solutions to the free boundary Navier—Stokes equations on an inclined plane, by use of bifurcation theory techniques... i.e. yes, I'm literally trying to make waves with my research. Come along to hear about what a Fredholm operator is, how bifurcation can (hopefully!) be used to find rigorous solutions to PDEs, and possibly another bad fluid mechanics pun.
-
Have you ever wanted to go on a treasure hunt? Then come to PSS this week as I go through a list of clues that led me to find the elusive treasure (vector) that I require. We will traverse through black holes (singularities, no physics involved!), steer clear of stacks, and blow up some obstacles in the way. I will explain quiver varieties, and what I've been thinking about for the last couple of years, no knowledge of algebraic geometry required!