Departmental Seminars & Events

We will decide on speakers for Trinity term 2024.

We develop a comprehensive geometrically-exact theory for an end-loaded elastic rod constrained to deform on a cylindrical surface. By viewing the rod-cylinder system as a special case of an elastic braid, we are able to obtain all forces and moments imparted by the deforming rod to the cylinder as well as all contact reactions. This framework allows us to give a complete treatment of static friction consistent with force and moment balance. In addition to the commonly considered model of hard frictionless contact, we analyse two friction models in which the rod, possibly with intrinsic curvature, experiences either lateral or tangential friction. As applications of the theory we study buckling of the constrained rod under compressive and torsional loads, finding critical loads to depend on Coulomb-like friction parameters, as well as the tendency of the rod to lift off the cylinder under further loading. The cylinder can also have arbitrary orientation relative to the direction of gravity. The cases of a horizontal and vertical cylinder, with gravity having only a lateral or axial component, are amenable to exact analysis, while numerical results map out the transition in buckling mechanism between the two extremes. Weight has a stabilising effect for near-horizontal cylinders, while for near-vertical cylinders it introduces the possibility of buckling purely due to self-weight. Our results are relevant for many engineering and medical applications in which a slender structure winds inside or outside a cylindrical boundary.

We introduce a new method - ESPIRA (Estimation of Signal Parameters via Iterative Rational Approximation) \cite{DP22,  DPP21} - for the recovery of complex exponential  sums
$$
f(t)=\sum_{j=1}^{M} \gamma_j \mathrm{e}^{\lambda_j t},
$$
that are determined by a finite number of parameters: the order $M$, weights $\gamma_j \in \mathbb{C} \setminus \{0\}$ and nodes  $\mathrm{e}^{\lambda_j} \in \mathbb{C}$ for $j=1,...,M$.  Our new recovery procedure is based on the observation that Fourier coefficients (or DFT coefficients) of exponential sums have a special rational structure.  To  reconstruct this structure in a stable way we use the AAA algorithm  proposed by Nakatsukasa et al.   We show that ESPIRA can be interpreted as a matrix pencil method applied to Loewner matrices. 

During the talk we will demonstrate that ESPIRA outperforms Prony-like methods such as ESPRIT and MPM for noisy data and for signal approximation by short exponential sums.  

Bibliography
N. Derevianko,  G.  Plonka, 
Exact reconstruction of extended exponential sums using rational approximation of their Fourier coefficients, Anal.  Appl.,  20(3),  2022,  543-577.

N. Derevianko,  G. Plonka,  M. Petz, 
From ESPRIT to ESPIRA: Estimation of signal parameters by iterative rational approximation,   IMA J. Numer. Anal.,  43(2),  2023, 789--827.  

Y. Nakatsukasa, O. Sète,   L.N. Trefethen,  The AAA algorithm for rational approximation.
SIAM J. Sci. Comput., 40(3),   2018,  A1494–A1522.  

A conjecture of Malle predicts an asymptotic formula for the number of number fields with given Galois group and bounded discriminant. Malle conjectured the shape of the formula but not the leading constant. We present a new conjecture on the leading constant motivated by a version for algebraic stacks of Peyre's constant from Manin's conjecture. This is joint work with Tim Santens.

We consider the reinforcement learning problem of controlling an unknown dynamical system to maximise the long-term average reward along a single trajectory. Most of the literature considers system interactions that occur in discrete time and discrete state-action spaces. Although this standpoint is suitable for games, it is often inadequate for systems in which interactions occur at a high frequency, if not in continuous time, or those whose state spaces are large if not inherently continuous. Perhaps the only exception is the linear quadratic framework for which results exist both in discrete and continuous time. However, its ability to handle continuous states comes with the drawback of a rigid dynamic and reward structure.

        This work aims to overcome these shortcomings by modelling interaction times with a Poisson clock of frequency $\varepsilon^{-1}$ which captures arbitrary time scales from discrete ($\varepsilon=1$) to continuous time ($\varepsilon\downarrow0$). In addition, we consider a generic reward function and model the state dynamics according to a jump process with an arbitrary transition kernel on $\mathbb{R}^d$. We show that the celebrated optimism protocol applies when the sub-tasks (learning and planning) can be performed effectively. We tackle learning by extending the eluder dimension framework and propose an approximate planning method based on a diffusive limit ($\varepsilon\downarrow0$) approximation of the jump process.

        Overall, our algorithm enjoys a regret of order $\tilde{\mathcal{O}}(\sqrt{T})$ or $\tilde{\mathcal{O}}(\varepsilon^{1/2} T+\sqrt{T})$ with the approximate planning. As the frequency of interactions blows up, the approximation error $\varepsilon^{1/2} T$ vanishes, showing that $\tilde{\mathcal{O}}(\sqrt{T})$ is attainable in near-continuous time.

We prove that any adjoint Chevalley group over an arbitrary commutative ring is regularly bi-interpretable with this ring. The same results hold for central quotients of arbitrary Chevalley groups and for Chevalley groups with bounded generation.
Also, we show that the corresponding classes of Chevalley groups (or their central quotients) are elementarily definable and even finitely axiomatizable.

We study supersymmetric Berry connections of 2d N = (2,2) gauged linear sigma models (GLSMs) quantized on a circle, which are periodic monopoles, with the aim to provide a fruitful physical arena for recent mathematical constructions related to the latter. These are difference modules encoding monopole solutions via a Hitchin-Kobayashi correspondence established by Mochizuki. We demonstrate how the difference modules arises naturally by studying the ground states as the cohomology of a one-parameter family of supercharges. In particular, we show how they are related to one kind of monopole spectral data, a deformation of the Cherkis–Kapustin spectral curve, and relate them to the physics of the GLSM. By considering states generated by D-branes and leveraging the difference modules, we derive novel difference equations for brane amplitudes. We then show that in the conformal limit, these degenerate into novel difference equations for hemisphere partition functions, which are exactly calculable. When the GLSM flows to a nonlinear sigma model with Kähler target X, we show that the difference modules are related to deformations of the equivariant quantum cohomology of X.

The Junior Algebra and Representation Theory Seminar will kick-off the start of Hilary term with a social event in the common room. Come to catch up with your fellow students and maybe play a board game or two. Afterwards we'll have lunch together.

From genetic studies and astrophysics simulations to statistical modelling and AI, scientific computing has enabled amazing discoveries and there is no doubt it will continue to do so. However, the corresponding environmental impact is a growing concern in light of the urgency of the climate crisis, so what can we all do about it? Tackling this issue and making it easier for scientists to engage with sustainable computing is what motivated the Green Algorithms project. Through the prism of the GREENER principles for environmentally sustainable science, we will discuss what we learned along the way, how to estimate the impact of our work and what levers scientists and institutions have to make their research more sustainable. We will also debate what hurdles exist and what is still needed moving forward.

Deep Gaussian processes have proved remarkably successful as a tool for various statistical inference tasks. This success relates in part to the flexibility of these processes and their ability to capture complex, non-stationary behaviours. 

In this talk, we will introduce the general framework of deep Gaussian processes, in which many examples can be constructed, and demonstrate their superiority in inverse problems including computational imaging and regression.

 We will discuss recent algorithmic developments for efficient sampling, as well as recent theoretical results which give crucial insight into the behaviour of the methodology.

 

We consider the sharp interface limit problem for 1D stochastic Allen-Cahn equation, and extend a classic result by Funaki to the full small noise regime. One interesting point is that the notion of "small noise" turns out to depend on the topology one uses. The main new idea in the proof is the construction of a series of functional correctors, which are designed to recursively cancel out potential divergences. At a technical level, in order to show these correctors are well behaved, we also develop a systematic decomposition of functional derivatives of the deterministic Allen-Cahn flow of all orders, which might have its own interest.
Based on a joint work with Wenhao Zhao (EPFL) and Shuhan Zhou (PKU).

A cap set is a subset of $\mathbb{F}_3^n$ with no solutions to $x + y + z = 0$ other than when $x = y = z$, or equivalently no non-trivial $3$-term arithmetic progressions. The cap set problem asks how large a cap set can be, and is an important problem in additive combinatorics and combinatorial number theory. In this talk, I will introduce the problem, give some background and motivation, and describe how I was able to provide the first progress in 20 years on the lower bound for the size of a maximal cap set. Building on a construction of Edel, we use improved computational methods and new theoretical ideas to show that, for large enough $n$, there is always a cap set in $\mathbb{F}_3^n$ of size at least $2.218^n$. I will then also discuss recent developments, including an extension of this result by Google DeepMind.

We study the stochastic quantisation for the fractional $\varphi^4$ theory. The model has been studied by Brydges, Mitter and Scopola in 2003 as a natural extension of $\phi^4$ theories to fractional sub-critical dimensions. The stochastic quantisation equation is given by the (formal) SPDE 

\[

(\partial_t + (-\Delta)^{s}) \varphi = - \lambda \varphi^3 + \xi

\]

where $\xi$ is a space-time white noise over the three dimensional torus. The equation is sub-critical for $s > \frac{3}{4}$.

We derive a priori estimates in the full sub-critical regime $s>\frac{3}{4}$. These estimates rule out explosion in finite time and they imply the existence of an invariant measure with a standard Krylov-Bogoliubov argument. 

Our proof is based on the strategy developed for the parabolic case $s=1$ in [Chandra, Moinat, Weber, ARMA 2023]. In order to implement this strategy here, a new Schauder estimate for the fractional heat operator is developed. Additionally, several algebraic arguments from [Chandra, Moinat, Weber, ARMA 2023] are streamlined significantly. 

This is joint work with Hendrik Weber (Münster). 

One of the oldest open problems in representation theory is to classify the irreducible unitary representations of a semisimple Lie group G_R. Such representations play a fundamental role in harmonic analysis and the Langlands program and arise in physics as the state space of quantum mechanical systems in the presence of G_R-symmetry. Most unitary representations of G_R are realized, via some kind of induction, from unitary representations of proper Levi subgroups. Thus, the major obstacle to understanding the unitary dual of G_R is identifying the "non-induced" unitary representations of G_R. In previous joint work with Losev and Matvieievskyi, we have proposed a general construction of these non-induced representations, which we call "unipotent" representations of G_R. Unfortunately, the methods we employ do not provide a proof that these representations are unitary. In this talk, I will explain how one can apply Saito's theory of mixed Hodge modules to overcome this difficulty, giving a uniform proof of the unitarity of all unipotent representations. This is joint work in progress with Dougal Davis

The saturation number of a graph is a famous and well-studied counterpoint to the Turán number, and the rainbow saturation number is a generalisation of the saturation number to the setting of coloured graphs. Specifically, for a given graph $F$, an edge-coloured graph is $F$-rainbow saturated if it does not contain a rainbow copy of $F$, but the addition of any non-edge in any colour creates a rainbow copy of $F$. The rainbow saturation number of $F$ is the minimum number of edges in an $F$-rainbow saturated graph on $n$ vertices. Girão, Lewis, and Popielarz conjectured that, like the saturation number, for all $F$ the rainbow saturation number is linear in $n$. I will present our attractive and elementary proof of this conjecture, and finish with a discussion of related results and open questions.

Measure equivalence was introduced by Gromov as a measure-theoretic analogue to quasi-isometry between finitely generated groups. In this talk I will present measure equivalence classification results for right-angled Artin groups, and more generally graph products. This is based on joint works with Jingyin Huang and with Amandine Escalier. 

Please note unusual day and room. 

Motivated by the SYZ conjecture, it is expected that $G_2$ and Spin(7)-manifolds also admit calibrated fibrations. One potential way to construct examples is via gluing of complex fibrations, as in the program of Kovalev. For this to succeed we need that the fibration property is stable under deformation of the ambient Spin(7)-structure. Here the main difficulty lies in the analysis of the singular fibres. In this talk I will present a stability result for fibrations with conically singular Cayleys modeled on the complex cone $\{x^2 + y^2 + z^2 = 0\}$ in ${\mathbb C}^3$.

Let X and H be large, and consider n ranging from 1 to X. For an arithmetic function f(n), what is the best mean square approximation of f(n) by a restricted divisor sum (a function of the sort sum_{d|n, d < H} a_d)? I hope to explain how for a wide variety of arithmetic functions, when X grows and H grows like a power of X, a solution of this problem is connected to the evaluation of random matrix integrals. The problem is connected to some combinatorial formula for computing high moments of traces of random unitary matrices and I hope to discuss this also.

One possible version of the Kirchberg—Phillips theorem states that simple, separable, nuclear, purely infinite C*-algebras are classified by KK-theory. In order to generalize this result to non-simple C*-algebras, Kirchberg first restricted his attention to those that absorb the Cuntz algebra O2 tensorially. C*-algebras in this class carry no KK-theoretical information in a strong sense, and they are classified by their ideal structure alone. It should be mentioned that, although this result is in Kirchberg’s work, its full proof was first published by Gabe. In joint work with Gábor Szabó, we showed a generalization of Kirchberg's O2-stable theorem that classifies G-C*-algebras up to cocycle conjugacy, where G is any second-countable, locally compact group. In our main result, we assume that actions are amenable, sufficiently outer, and absorb the trivial action on O2 up to cocycle conjugacy. In very recent work, I moreover show that the range of the classification invariant, consisting of a topological dynamical system over primitive ideals, is exhausted for any second-countable, locally compact group.

In this talk, I will recall the classification of O2-stable C*-algebras, and describe their classification invariant. Subsequently, I will give a short introduction to the C*-dynamical working framework and present the classification result for equivariant O2-stable actions. Time permitting, I will give an idea of how one can build a C*-dynamical system in the scope of our classification with a prescribed invariant. 

TBA

Traditionally, stochastic processes are modelled one of two ways: a continuum Fokker-Planck approach, where a PDE is solved to determine the time evolution of the probability density, or a Langevin approach, where the SDE describing the system is sampled, and multiple simulations are used to collect statistics. There is also a third way: the functional or path integral. Originally developed by Wiener in the 1920s to model Brownian motion, path integrals were famously applied to quantum mechanics by Feynman in the 1950s. However, they also have much to offer classical stochastic processes (and statistical physics).  

In this talk I will introduce the formalism at a physicist’s level of rigour, and focus on determining the dominant contribution to the path integral when the noise is weak. There exists a remarkable correspondence between the most-probable stochastic paths and Hamiltonian dynamics in an effective potential [1,2,3]. I will then discuss some applications, including reaction pathways conditioned on finite time [2]. We demonstrate that the most probable pathway at a finite time may be very different from the usual minimum energy path used to calculate the average reaction rate. If time permits, I will also discuss the extremely nonlinear crystal dislocation response to applied stress [4].  

[1] Ge, Hao, and Hong Qian. Int. J. Mod. Phys. B 26.24 1230012 (2012)     

[2] Fitzgerald, Steve, et al. J. Chem. Phys. 158.12 (2023).

[3] Honour, Tom and Fitzgerald, Steve. in press J. Phys. A (2024)

[4] Fitzgerald, Steve. Sci. Rep. 6 (1) 39708 (2016)

Artificial Intelligence (AI) will strongly determine our future prosperity and well-being. Due to its generic nature, AI will have an impact on all sciences and business sectors, our private lives and society as a whole. AI is pre-eminently a multidisciplinary technology that connects scientists from a wide variety of research areas, from behavioural science and ethics to mathematics and computer science.

Without downplaying the importance of that variety, it is apparent that mathematics can and should play an active role. All the more so as, alongside the successes of AI, also critical voices are increasingly heard. As Robert Dijkgraaf (former director of the Princeton Institute of Advanced Studies) observed in May 2019: ”Artificial intelligence is in its adolescent phase, characterised by trial and error, self-aggrandisement, credulity and lack of systematic understanding.” Mathematics can contribute to the much-needed systematic understanding of AI, for example, greatly improving reliability and robustness of AI algorithms, understanding the operation and sensitivity of networks, reducing the need for abundant data sets, or incorporating physical properties into neural networks needed for super-fast and accurate simulations in the context of digital twinning.

Mathematicians absolutely recognize the potential of artificial intelligence, machine learning and (deep) neural networks for future developments in science, technology and industry. At the same time, a sound mathematical treatment is essential for all aspects of artificial intelligence, including imaging, speech recognition, analysis of texts or autonomous driving, implying it is essential to involve mathematicians in all these areas. In this presentation, we highlight the role of mathematics as a key enabling technology within the emerging field of scientific machine learning. Or, as I always say: ”Real intelligence is needed to make artificial intelligence work.”

We show the super-hedging duality for multi-action options which generalise American options to a larger space of actions (possibly uncountable) than {stop, continue}. We put ourselves in the framework of Bouchard & Nutz model relying on analytic measurable selection theorem. Finally we consider information delay on the action component of the product space. Information delay is expressed as a possibility to look into the future in the dual formulation. This is a joint work with Ivan Guo, Shidan Liu and Zhou Zhou.

The talk is based on joint work with Lasse Grimmelt. We prove a theorem that allows one to count solutions to determinant equations twisted by a periodic weight with high uniformity in the modulus. It is obtained by using spectral methods of $\operatorname{SL}_2(\mathbb{R})$ automorphic forms to study Poincaré series over congruence subgroups while keeping track of interactions between multiple orbits. This approach offers increased flexibility over the widely used sums of Kloosterman sums techniques. We give applications to correlations of the divisor function twisted by periodic functions and the fourth moment of Dirichlet $L$-functions on the critical line.

Periodic point clouds naturally arise when modelling large homogenous structures like crystals. They are naturally attributed with a map to a d-dimensional torus given by the quotient of translational symmetries, however there are many surprisingly subtle problems one encounters when studying their (persistent) homology. It turns out that bisheaves are a useful tool to study periodic data sets, as they unify several different approaches to study such spaces. The theory of bisheaves and persistent local systems was recently introduced by MacPherson and Patel as a method to study data with an attributed map to a manifold through the fibres of this map. The theory allows one to study the data locally, while also naturally being able to appeal to local systems of (co)sheaves to study the global behaviour of this data. It is particularly useful, as it permits a persistence theory which generalises the notion of persistent homology. In this talk I will present recent work on the theory and implementation of bisheaves and local systems to study 1-periodic simplicial complexes. Finally, I will outline current work on generalising this theory to study more general periodic systems for d-periodic simplicial complexes for d>1.