I will discuss Agol's proof of the Virtually Fibred Conjecture of
Thurston, focusing on the role played by the `RFRS' property. I will
then show how one can modify parts of Agol's proof by replacing some
topological considerations with a group theoretic statement about
virtual fibring of RFRS groups.
 

Atom Probe Tomography is a powerful 3D mass spectrometry technique. By pulsing the sample apex with an electric field, surface atoms are ionised and collected by a detector. A 3D image of estimated initial ion positions is constructed via an image reconstruction protocol. Current protocols assume ion trajectories follow a stereographic projection. However, this method assumes a hemispherical sample apex that fails to account for varying material ionisation rates and introduces severe distortions into atomic distributions for complex material systems.

We aim to develop continuum models and use this to derive a time-dependent mapping describing how ion initial positions on the sample surface correspond to final impact positions on the detector. When correctly calibrated with experiment, such a mapping could be used for performing reconstruction.

Currently we track the sample surface using a level set method, while the electric field is solved via BEM or a FEM-BEM coupling. These field calculations must remain accurate close to the boundary. Calibrating unknown evaporation parameters with experiment requires an ensemble of models per experiment. Therefore, we are also looking to maximise model efficiency via BEM compression methods i.e. fast multipole BEM. Efficiently constructing and reliably interpolating the non-bijective trajectory mapping, while accounting for ion trajectory overlap and instabilities (at sample surface corners), also presents intriguing problems.

This project is in collaboration with Cameca, the leading manufacturer of commercial atom probe instruments. If successful in minimising distortions such a technique could become valuable within the semiconductor industry.

This work aims at developing new approach for parameterizing mesoscale eddy effects for use in non-eddy-resolving ocean circulation models. These effects are often modelled as some diffusion process or a stochastic forcing, and the proposed approach is implicitly related to the latter category. The idea is to approximate transient eddy flux divergence in a simple way, to find its actual dynamical footprints by solving a simplified but dynamically relevant problem, and to relate the ensemble of footprints to the large-scale flow properties.

The Kelvin wave is perhaps the most important of the equatorially trapped waves in the terrestrial atmosphere and ocean, and plays a role in various phenomena such as tropical convection and El Nino. Theoretically, it can be understood from the linear dynamics of a stratified fluid on an equatorial beta plane, which, with simple assumptions about the disturbance structure, leads to wavelike solutions propagating along the equator, with exponential decay in latitude. However, when the simplest possible background flow is added (with uniform latitudinal shear), the Kelvin wave (but not the other equatorial waves) becomes unstable. This happens in an extremely unusual way: there is instability for arbitrarily small nondimensional shear p, and the growth rate is proportional to exp(-1/p^2) as p->0. This in contrast to most hydrodynamic instabilities, in which the growth rate typically scales as a positive power of p-p_c as the control parameter p passes through a critical value p_c.

This Kelvin wave instability has been established numerically by Natarov and Boyd, who also speculated as to the underlying mathematical cause. Here we show how the growth rate and full spatial structure of the instability may be derived using matched asymptotic expansions applied to the (linear) equations of motion. This involves an adventure with Whittaker functions in the exponentially-decaying tails of the Kelvin waves, and a trick to reveal the exponentially small growth rate from a formulation that only uses regular perturbation expansions. Numerical verification of the analysis is also interesting and challenging, since special high-precision solutions of the governing ODE are required even when the nondimensional shear is not that small (circa 0.5).

Jacob Bernoulli is known for his studies of the curves, infinitesimal math- ematics and statistics. However, before being a professor in mathematics, he taught experimental physics at the University of Basel. This explains his high interest in solving physical problems with newly developed Leibnizian calculus. In his scientific notebook, Meditationes, there are more than thirty notes about various mechanical problems for solving of which Bernoulli has applied Leibnizian calculus and has advanced this method along the way. A discussion with a craftsman brought Bernoulli’s attention to the problem of the strength of a beam early in his career and occupied his mind until his death. The craftsman’s narration based on his experience highlighted the flaws in Galilean-Leibnizian theory of the strength of a beam. This was the starting point of Bernoulli’s quest to mathematically find the profile of a bent beam (the Elastica Problem) and the physical laws governing it. He started a challenge to encourage other mathematicians of the time to study the problem, providing a hint hidden in an anagram. Although he published his solution of the Elastica Problem in 1694, that was not the end of the quest for him. Studying his unpublished notes in Meditationes reveals that over the last decade of his life, Bernoulli has reconsidered the problem. In my project, I demonstrate that he has found remarkable concepts such as mean tensile stress, and the notion of local stress-strain relation, etc.