In light of the coronavirus outbreak, we are following university adviceLink opens in a new window and suspending all SPAAM talks for term 3. If and when we learn more, we will be in touch.
- Last updated: 17 March 20
Term 2: Tuesdays 2-3pm, B3.01
Anna Skorobogatova | Limiting Sobolev-type inequalities for vector fields
Abstract: The classical Gagliardo-Nirenberg-Sobolev inequality on $n$ dimensions guarantees that for any $1 \leq p < n$, we can estimate the $L^{\frac{np}{n-p}}$ norm of any smooth compactly supported map $u$ by the $L^p$ norm of $Du$. Now consider a more general first order differential operator $\mathbb{B}$ instead of the full gradient. Can we get an analogous inequality when we replace $D$ by $\mathbb{B}$? Or do we really need all of the information that the full gradient provides? Which are the ‘correct’ conditions that need be imposed on the operator $\mathbb{B}$? In particular, we will see that the limiting case $p=1$ is special and needs to be handled separately.
Peter Mühlbacher | Random walks on the symmetric group and induced integer partitions
Abstract: Random walks on the symmetric group and induced integer partitions Abstract: How often do you have to swap (iid) elements to observe large permutation cycles? Does this mean we're close to a uniform permutation? Can we directly describe how these permutation cycles evolve in time? Heuristics and partial results suggesting some kind of universality will be presented.
Jaro Sant | Ergodicity & Speed of Convergence for 1D Diffusions
Abstract: Ergodicity is a central theme in the study of diffusion processes, and ensures that the long time behaviour of a process can be studied meaningfully. I'll introduce the notions of ergodicity and Harris recurrence for general state space stochastic processes, show how the two notions are linked in a general setting, and illustrate how this allows one to relate time averages to state space averages via Birkhoff's ergodic theorem. I will then focus on obtaining rates of convergence in this theorem, specifically for the case of 1D diffusions where the scenario is much easier to deal with. Through the use of speed and scale one can control the rate of convergence in the ergodic theorem by means of moments of hitting times of singletons. I will briefly comment on the extra problems that need to be tackled in higher dimensions time permitting.
Cathie Wells, University of Reading Reducing Aviation Emissions and Fuel Burn by Re-routing Transatlantic Flights
Introduction: Commercial long haul flight is at the brink of major change. After decades of limited communications across the North Atlantic, full satellite coverage will soon be available. At the same time, airlines are facing increasing pressure to ensure their flights are more fuel efficient and thus less polluting. Here we show that changes to current fixed track trajectories could significantly reduce fuel use for aircraft making the journey between London and New York.
Methods: Each simulation is completed for a constant air speed and at FL 350, an altitude of 35 000 feet. Optimal Control Theory is used to minimise time of flight in routing aircraft through actual wind fields obtained from a global atmospheric re-analysis dataset. New formulae are then applied to calculate fuel burn rate at the given air speed and the fuel used for the trajectory is found. Applying this method to routes obtained at a variety of air speeds, the most fuel efficient simulation can be chosen.
Results:Total fuel burned in simulations was compared with actual flight data. Three different models of fixed wing aircraft were considered, the Boeing 747-436, 777-236 and 787-9. Each has a unique physical profile leading to different fuel burn rates. This is reflected in the varied levels of savings possible. The importance of the change to air speed, compared with the re-routing of the flight, was also considered.
Conclusions: By making use of the available wind fields, this low cost method not only allows punctual arrivals, but also reduces fuel use and emissions, ensuring a solution that is for once both good for business and for the environment.
Matthew King | Sound and critical layers in sheared flow over acoustic liners
Abstract: Acoustics within mean flows, governed by the linearised Euler equations, has applications within various real world settings such as aeroacoustics. One area of study is the modelling of acoustic liners within turbo-fan jet engines. For this, a uniform flow profile is often considered within a cylindrical duct. More recently however, the need to consider a sheared flow profile has become apparent. We make use of the Pridmore-Brown equation in order to study the acoustics of such a system. This talk will look at solving the Pridmore-Brown equation near a regular singularity known as the critical layer, the behaviour of these solutions, and the difficulties induced by the presence of the critical layer.
Marco Palma | Quantifying uncertainty in brain-predicted age using scalar-on-image quantile regression
Abstract: Prediction of subject age from brain anatomical MRI has the potential to provide a sensitive summary of brain changes, indicative of different neurodegenerative diseases. However, existing studies typically neglect the uncertainty of these predictions.
In this work we take into account this uncertainty by applying methods of functional data analysis. In particular, we propose a penalised functional quantile regression model of age on brain structure with cognitively normal (CN) subjects in the Alzheimer’s Disease Neuroimaging Initiative (ADNI), and use it to predict brain age in Mild Cognitive Impairment (MCI) and Alzheimer’s Disease (AD) subjects. Unlike the machine learning approaches available in the literature of brain age prediction, which provide only point predictions, the outcome of our model is a prediction interval for each subject.
In this talk an overview of all the theoretical aspects (functional data analysis, multidimensional smoothing, quantile regression) will be provided, followed by a discussion of the main results and of their meaning in the field of neuroscience of ageing.
Carla Groenland, University of Oxford | Cyclically Covering Subspaces
Abstract: How 'small' can you take a linear subspace of F_2^n such that it contains at least one cyclic shift of all vectors of F_2^n, where the cyclic shift rotates the coordinates? I will explain why Artin's conjecture about primes has implications for this question, and how this led us to conjecture a generalisation of the so-called handshake lemma to Cayley graphs. All the definitions in this abstract will be explained within the talk.
Suzie Brown | Resampling in Sequential Monte Carlo
Abstract: Sequential Monte Carlo refers to a wide class of stochastic algorithms that generate samples from a sequence of posterior distributions by simulating a population of particles evolving through time. A key step in the algorithm is resampling, where a new set of particles is generated from the weighted particles of the previous generation. There are many ways to do resampling, leading to different properties that can affect the behaviour of the resulting Monte Carlo estimates. I will introduce and discuss a range of commonly used resampling schemes, and compare their properties. Time permitting, I will explore the link between resampling and particle genealogies.
Ian Hamilton | Predicting the past: A retrodictive model for modern rugby union
Abstract: It is not uncommon for competitors in sports tournaments to have results related to schedules of varying strength, consisting of stronger (weaker) opposition, more (fewer) matches, and a higher (lower) proportion of matches at home. One popular means of producing a rating model in such circumstances is the Bradley-Terry model. This talk introduces that model, discusses the extensions in the literature, and proposes a further one to handle modern rugby union, where teams gain points for wins and draws, but also bonus points for losing within a certain score, and for the number of tries scored. Additionally it will touch on the differences between predictive and retrodictive models, present a principled entropy-based motivation for the model, and argue for a pragmatic understanding of the use of a prior in this setting. The proposed model is applied to a national schools rugby tournament in order to investigate their current ranking method and to suggest improvements.
Trishen Gunaratnam, University of Bath | Phase transition for $\phi^4$
Abstract: $\phi^4$ is a continuum statistical mechanics model with unbounded spin. It arises in the study of (Euclidean) quantum field theory, as the continuum limit of Ising-type models near criticality, as the invariant measure of singular stochastic PDEs, and as the invariant measure of nonlinear Schrodinger equations. At large scales this model has a remarkable resemblance to the famous Ising model: It undergoes phase transition in 2 and 3 dimensions as the temperature is varied.
In this talk we describe how contour arguments can be used to prove phase transition for $\phi^4$. Precisely, we look at the sophisticated modification of the classical Peierls' argument (used to prove phase transition for the Ising model) adapted to the case of $\phi^4$ in 2 dimensions by Glimm, Jaffe and Spencer (CMP 1975). Time permitting, we discuss the extension to 3 dimensions and phase separation results that have recently been obtained in joint work with Ajay Chandra and Hendrik Weber.
Eleanor Archer | Edge-reinforced random walk and the localisation phenomenon
Abstract: Edge-reinforced random walk (ERRW) is a classical example of a self-interacting process. Focusing mainly on the case of ERRW on trees, we will introduce the ERRW and discuss some fundamental connections with random urn models which enable us to describe it as a random "mixture" of non-interacting random walks. This representation allows us to construct the scaling limit of ERRW on trees. Time permitting we will also discuss some properties of the scaling limit, including a quenched localisation result showing that the walk moves logarithmically slowly.
Dr. Alice Corbella | Towards automatic Zig Zag sapling and its use for epidemic inference
Abstract: Zig-Zag sampling, introduced by Bierkens et al. 2019, is based on the simulation of a piecewise deterministic Markov process (PDMP) whose switching rate $\lambda(t)$ is governed by the derivative of the log-target density. To our knowledge, Zig-Zag sampling has been used mainly on simple targets for which derivatives can be computed manually in a reasonable time.
To expand the applicability of this method, we incorporate Automatic Differentiation (AD) tools in the Zig-Zag algorithm, computing $ \lambda(t) $ automatically from the functional form of the log-target density. Moreover, to allow the simulation of the PDMP via thinning, we use standard optimization routines to find a local upper bound for the rate.
We present several implementations of our automatic Zig-Zag sampling and we measure the potential loss in computational time caused by AD and optimization routines. Among the examples, we consider the case of data arising from an epidemic which can be approximated by a deterministic system of equations; here manual derivation of the posterior density is practically infeasible due to the recursive relationships contained the likelihood function. Automatic Zig-Zag sampling successfully explores the parameter space and samples efficiently from the posterior distribution.
(Joint work with Gareth O. Roberts and Simon E. F. Spencer)
Abstract: We will talk about percolation, loops and the interplay in-between the two. A number of open questions will be introduced as well as potential routes of attacking them.
Samuel Forbes | MathSys 3rd Year | Wealth Distribution Models
Abstract: I will introduce the wealth distribution, discuss empirical data for the UK and summarise some wealth distribution models.
Finally, the SPAAM poem puzzle, we state, Just a bit of fun - for you to procrastinate! They’ll sometimes be easy, sometimes fascinate, But we do hope you’ll try and investigate, For answering correctly will surely elate! Don’t worry if you fail, you can always wait, With a new week comes the next iterate. And come along for the answer – don’t be late!
Ten coins are in a two-by-five grid, Exactly four of them are to be slid, Behold! Five straight lines of four coins now, But the question is: can you tell me how?
Jaromir Sant | MASDOC 3rd Year Mathematical Genetics - A Song of Wright & Fisher
Abstract:
Have you ever wondered why your eyes are blue, And not some other beautiful hue? Or perhaps you’d like to know Why it is that your dad’s got such a prominent brow. Or why your brother has a chin That could easily fit an industrial bin. Then worry not for I shall tell a story About genetics in all it’s shining glory. From Darwin to Fisher, from Mendel to Wright, Their theories started to shed some light. Then came Kimura claiming selection was dead, And the issue has been all but put to bed. This and some more I shall present, To those who choose the SPAAM seminar to attend, On Tuesday afternoon when the clock strikes three, Just an hour before we serve coffee and tea.
Hanson Bharth | MASDOC 3rd Year | Rolls of Castamere
Abstract:
For this talk I shall Introduce a mathematical model to Simulate a cylinder rolling along an elastic Half-space. By taking a Fourier transform we Obtain a general solution but require a Wiener-Hopf technique to Learn more about the solution.
SPAAM poem puzzle #2:
“I’ve found a really nice puzzle”, said Flo,
“Oh that's great!” exclaimed Mo, “Come on, let’s have a go!”
“C is a circle with centre O”,
“This is exciting!” shouted Mo “I like geometry - didn’t you know?”
“Corners chosen uniformly from C, form a random triangle - let’s call it T”
“Probability!? This I did not foresee, And is almost definitely not for me!”
“The probability that O is in T, That’s exactly one half - don't you agree?”
“No!” replied Mo, “on the contrary, It’s closer to 1 over 3!”
“You are all wrong!” interrupts Lee, "But the actual answer alludes me, Can you solve the mystery?" "And how about a bonus round?" "Can the expected area of T be found?"
And as the question began to hound From the 3 came no more sound.
Mattia Sanna | 3rd Year PhD | Computation on Galois representations
Abstract: Galois representations have a special role in algebraic number theory. They arise from different theories for example arithmetic geometry and particluar holomorphic functions called modular forms. It is quite natural then to ask whether two such Galois representations are isomorphic, and therefore provide a connection among several areas of study. One of the most famous result in this direction is the proof of Fermat's last theorem due to Andrew Wiles in which he proves that all the Galois representations attached to elliptic curves defined over comes from the ones associated to particular modular forms. We would like to extend Wile's result in the most general setting but unfortunately it is a very hard task. However, there exist some explicit methods to check whether two given Galois representations are isomorphic. In this talk we will present a brief introduction on the subject and the ideas behind the methods.
Emma Southall | 3rd Year MathSys Early Warning Signals of Disease Elimination
Abstract: TBC.
SPAAM Poem Puzzle 3:
A conspiracy theorist named Lee, Frank, a risk mitigation analyst was he, Jane, an inventor with robots she had three,
The 3 friends and 3 robots on a riverbank, All stood still and looking blank, ”How would we cross the river?” thought Frank,
A small boat and an oar is all they had to hand, The opposite riverbank, they wanted to land, Luckily, the robots could operate the boat - this was well planned,
Unluckily, the boat was small and could only carry two people/robots at a time, This line is only here so that I can rhyme, Let’s go back to the story of maritime,
The conspiracy theorist didn’t trust the robots at all, He was convinced that they would cause a brawl,
The risk analyst more sceptical now, Thought "I’m sure we can sort something - but how?"
Lee insisted that no one should, Ever be outnumbered by robots - that would be good! The inventor wondered why her robots were so misunderstood, But carried on anyway and quickly solved it for good.
Can you find a solution? How many times did Jane have to cross the river?
Michael Luya | MathSys The Mathematics behind Rotation Invariant Basis Representations for Machine Learning Models within Computational Chemistry
Abstract: How much mathematics is present in chemistry? In this presentation, Michael hopes to answer this question from a background of mathematics from the perspective of his PhD project, which is based in the design of machine learning infrastructures in computational chemistry. We will be looking at a variety of mathematical models that provide us with symmetry invariant parameters, that specifically satisfy rotation invariance, deriving them from quantum mechanical principles. Quantum mechanics itself provides us with a high level of accuracy, which we are currently interfacing with machine learning methods, which will allow us to make a variety of calculations. This talk however, will focus on the mathematics that is integral in getting us to such a stage.
Connah Johnson | MathSys Spatial inhomogeneities in biochemical reaction networks
Abstract: Living systems may be thought of as perturbations in a chemical space. Of the set of possible chemical combinations a select few occur within biology, held within specific ranges of environments. These environments themselves being host to maelstrom of chemical interactions which may be modeled as a reaction network. Some ways in which these physical and chemical aspects combine and subsequent emergent phenomena will be the subject of this talk.
SPAAM Poem Puzzle 4:
I spent a £5 note on 2 ice creams, They were Chocolate and Caramel Supremes, At £1.20 each, this was certainly the dream, The number of coins in change I got was 19, At least one 2p, 5p, 20p, 50p and £1 coins could be seen, And exactly one 1p and one 10p I got, How many 2p coins could be in the lot?
Nicolò Paviato | MASDOC 2nd Year | Convergence to Brownian motion in a Deterministic System
Abstract: In this talk we will discuss statistical properties of dynamical systems, i.e. an application of probability theory to a deterministic case. If a system displays a "random behaviour", then both the Central Limit Theorem and Donsker's Invariance Principle are valid for regular observables; these follow from an intrinsic martingale property of the dynamic. Finally we will provide a gist on how to improve rates of convergence for a deterministic functional CLT and to find some multidimensional estimates.
SPAAM Poem Puzzle 5:
I used to play noughts and crosses, But all I got was a number of losses, So all the rules I decided to change, And you might find this strange, When players play uniformly at random, "No skill!" you cry and cannot fathom, But see, the skill is in choosing whether to go second or first, Can you figure out which one is best and which one is worst?
Diogo Caetano | MASDOC 2nd Year The Cahn-Hilliard equation on an evolving surface
Abstract: In this talk, we develop some tools to pose PDEs on evolving hypersurfaces and derive the Cahn-Hilliard equation with a logarithmic potential on a 2-dimensional moving surface. We then compare this problem with the corresponding case of a fixed domain, and conclude by discussing some results regarding existence and uniqueness of (weak) solutions.
Mohammad Noorbakhsh | MathSys Causal network discovery from climate time series
Abstract: TBC.
SPAAM Poem Puzzle #6:
There once was a frog called Etienne, Who used to sit on a riverbank in the Glen, Lilly pads on water there was one shy off N, Formed a straight line just in front of Etienne, He chooses a number uniformly from 1 to N, And jumps to this numbered lillypad then, But if he chooses N, he crosses the river in the Glen. And if he lands on a lillypad, the process repeats again, Until on the other side of the riverbank is our dear Etienne
Given he jumps forwards on every jump, what is the expected number of jumps for Etienne to cross the river?
Letizia Angeli | MASDOC | Interacting Particle Systems Approximations of Feynman-Kac Formulae in continuous time
Abstract: In this talk, I will present a class of numerical algorithms - based on the evolution of interacting particle systems - for estimating large deviation quantities of additive path functionals of stochastic processes. Adapting already established results from the literature of particle filters and sequential Monte Carlo methods, we can study the convergence of the algorithms and provide a rigorous framework to evaluate and improve their efficiency.
Osian Shelley | MASDOC | Transactions costs, boon or bust?
Abstract: The financial crisis of 2008 has helped further the debate on the desirability of transaction taxes in financial markets. Indeed, recurring questions and themes relating to such taxation often gain relevance after periods of economic instability. Unfortunately, arguments for and against these taxes are often based on ad hoc heuristics. This is due in part to a lack of general equilibrium models, from which one can gather quantitative and qualitative guidance.
In the talk I aim to show the effects of a quadratic tax rate levied against two agents with heterogeneous risk aversions in a continuous-time, risk-sharing equilibrium model. In particular, I’ll talk about the utility loss incurred by the agents upon the introduction of the tax and whether this tax can benefit the agents upon redistribution.
Yani Pehova | Mathematics PhD Quasirandomness in large permutations via the flag algebra method
Abstract: A sequence (Pi) of permutations is said to be quasirandom if the pattern density of every 4-point permutation in (Pi) converges to 1/24. We show that there is a set S of eight 4-point permutations such that the sum of the pattern densities of the permutations from S in (Pi) converges to |S|/24 if and only if the sequence is quasirandom. In this talk I will go over the proof methods we use to show that S has the above property. In particular, I will give an outline of the philosophy of the flag algebra method, and some details on how to apply it.
This is joint work with Timothy Chan, Jon Noel, Maryam Sharifzadeh, Daniel Kral and Jan Volec.
Phil Herbert | MASDOC | The membrane mediated force on point attachments with application to a near spherical biomembrane
Abstract: We consider a hybrid model of a biomembrane with attached particles. The membrane is represented by a near spherical continuous surface, attached proteins are described by rigid bodies which are free to move tangentially and rotate in the axis normal to a reference point. As the standard energy for a membrane is highly non-linear, we consider a quadratic energy which may be shown to be an approximation of the Canham-Helfrich energy with a volume constraint and the deformations due to the attached proteins are imposed by point Dirichlet constraints. We show differentiability of the membrane energy with respect to parameterisation of these embedded particles and provide illustrative numerical examples.
Dom Brockington | MASDOC | Sticky Flows and KPZ Universality
Abstract: We shall introduce stochastic flows, and stochastic flows of kernels, before moving to the special case of sticky flows and their fluctuations. On the large deviation scale these fluctuations turn out to be Tracey-Widom GUE distributed, thus the sticky flows lie in the KPZ universality class.
We're happy to announce and formally invite everyone to the SPAAM seminar series, a new forum where Maths & Stats PhD students as well as early researchers can present and discuss their work. The series is being organised by the SIAM Student Chapter here at Warwick.
Tea & Coffee? Yes – afterwards in the Maths common room
The seminar series aims to bring together CDT streams as well as straight PhD students across Mathematics and Statistics.
If you're interested in giving a talk, please let us know!