Past

As parallel computers approach Exascale (10^18 floating point operations per second), processor failure and data corruption are of increasing concern. Numerical linear algebra solvers are at the heart of many scientific and engineering applications, and with the increasing failure rates, they may fail to compute a solution or produce an incorrect solution. It is therefore crucial to develop novel parallel linear algebra solvers capable of providing correct solutions on unreliable computing systems. The common way to mitigate failures in high performance computing systems consists of periodically saving data onto a reliable storage device such as a remote disk. But considering the increasing failure rate and the ever-growing volume of data involved in numerical simulations, the state-of-the-art fault-tolerant strategies are becoming time consuming, therefore unsuitable for large-scale simulations. In this talk, we will present a  novel class of fault-tolerant algorithms that do not require any additional resources. The key idea is to leverage the knowledge of numerical properties of solvers involved in a simulation to regenerate lost data due to system failures. We will also share the lessons learned and report on the numerical properties and the performance of the new resilience algorithms.

Function solutions to linear PDEs often carry rigidity properties directly associated to the equation they satsify. However, the realm of solutions covers a much larger sets of solutions. For instance, we can speak of measure solutions, as opposed to classical $C^\infty$ functions or even $L^p$ functions. It is only logical to expect that the “better” space the solution lives in, the more rigid its properties will be.

Measure solutions lie just at a comfortable half of this threshold: it is a sufficently large space which allows for a rich range of new structures; but is sufficiently rigid to preserve a meaningful geometrical pattern. For example, have you ever wondered how gradients look like in the space of measures? What about other PDE structures? In this talk I will discuss these general questions, a few examples of them, and a new theoretical approach to its understanding via PDE theory, harmonic analysis, and geometric measure theory methods.

In 1996 using techniques from model theory and intersection theory, Hrushovski obtained a generalisation of the Lang-Weil estimates. Subsequently the estimates have found applications in group theory, algebraic dynamics and algebraic geometry. We shall discuss a geometric proof of the non-uniform version of these estimates.

Property (FA) is one of the `rigidity properties’ defined for groups, concerning the way a group can act on trees. We’ll take a look at why you might be interested in an action on a tree, what the property is, and then investigate which automorphism groups of free products have it.

We study the problem of super-resolution using TV norm minimisation, where we recover the locations and weights of non-negative point sources from a few samples of their convolution with a Gaussian kernel. A practical approach is to solve the dual problem. In this talk, we study the stability of solutions with respect to the solutions to the dual problem. In particular, we establish a relationship between perturbations in the dual variable and the primal variables around the optimiser. This is achieved by applying a quantitative version of the implicit function theorem in a non-trivial way.

Deformation quantization modules or DQ-modules where introduced by M. Kontsevich to study the deformation quantization of complex Poisson varieties. It has been advocated that categories of DQ-modules should provide invariants of complex symplectic varieties and in particular a sort of complex analog of the Fukaya category. Hence, it is natural to aim at describing the functors between such categories and relate them with categories appearing naturally in algebraic geometry. Relying, on methods of homotopical algebra, we obtain an analog of Orlov representation theorem for functors between categories of DQ-modules and relate these categories to deformations of the category of quasi-coherent sheaves.
 

We consider a sparse multi-task regression framework for fitting a collection of related sparse models. Representing models as nodes in a graph with edges between related models, a framework that fuses lasso regressions with the total variation penalty is investigated. Under a form of generalised restricted eigenvalue assumption, bounds on prediction and squared error are given that depend upon the sparsity of each model and the differences between related models. This assumption relates to the smallest eigenvalue restricted to the intersection of two cone sets of the covariance matrix constructed from each of the agents' covariances. In the case of a grid topology high-probability bounds are given that match, up to log factors, the no-communication setting of fitting a lasso on each model, divided by the number of agents.  A decentralised dual method that exploits a convex-concave formulation of the penalised problem is proposed to fit the models and its effectiveness demonstrated on simulations. (Joint work with Sahand Negahban and Patrick Rebeschini)

Multiple scales analysis is a powerful asymptotic technique for problems where the solution depends on two scales of widely different sizes. Standard multiple scales involves the introduction of a macroscale and microscale which are assumed to be independent. A common (and usually acceptable) assumption is that when considering behaviour on the microscale, the macroscale variable can be taken as constant, however there are instances where this assumption is not valid. In this talk, I will explain one such situation, that is, when considering conductive-radiative thermal transfer within a solid matrix with spherical perforations and discuss the appropriate measures when converting the radiative boundary condition into multiple-scales form.
 

Quantum fields can sometimes have negative energy density.  In gravitational contexts, this threatens to permit both causality violations (such as traversable wormholes, warp drives, and time machines) and violations of the Second Law for black holes.  I will discuss the thermodynamic principles that rule out such pathological situations.  These principles have led us to an interesting lower bound on the energy flux, even for field theories in flat spacetime! This Quantum Null Energy Condition has now been proven for all relativistic field theories.  I will give an intuitive argument explaining why such ``quantum energy conditions'' ought to hold. 
 

Complexity Science aims to understand what is that makes some systems to be "more than the sum of their parts". A natural first step to address this issue is to study networks of pairwise interactions, which have been done with great success in many disciplines -- to the extend that many people today identify Complexity Science with network analysis. In contrast, multivariate complexity provides a vast and mostly unexplored territory. As a matter of fact, the "modes of interdependency" that can exist between three or more variables are often nontrivial, poorly understood and, yet, are paramount for our understanding of complex systems in general, and emergence in particular. 
In this talk we present an information-theoretic framework to analyse high-order correlations, i.e. statistical dependencies that exist between groups of variables that cannot be reduced to pairwise interactions. Following the spirit of information theory, our approach is data-driven and model-agnostic, being applicable to discrete, continuous, and categorical data. We review the evolution of related ideas in the context of theoretical neuroscience, and discuss the most prominent extensions of information-theoretic metrics to multivariate settings. Then, we introduce the O-information, a novel metric that quantify various structural (i.e. synchronous) high-order effects. Finally, we provide a critical discussion on the framework of Integrated Information Theory (IIT), which suggests an approach to extend the analysis to dynamical settings. To illustrate the presented methods, we show how the analysis of high-order correlations can reveal critical structures in various scenarios, including cellular automata, Baroque music scores, and various EEG datasets.

References:
[1] F. Rosas, P.A. Mediano, M. Gastpar and H.J. Jensen, ``Quantifying High-order Interdependencies via Multivariate Extensions of the Mutual Information'', submitted to PRE, under review.
https://arxiv.org/abs/1902.11239
[2] F. Rosas, P.A. Mediano, M. Ugarte and H.J. Jensen, ``An information-theoretic approach to self-organisation: Emergence of complex interdependencies in coupled dynamical systems'', in Entropy, vol. 20 no. 10: 793, pp.1-25, Sept. 2018.
https://www.mdpi.com/1099-4300/20/10/793

 

The principal ideal theorem (1930) ascertained that any number field K embeds into a finite extension, called the Hilbert class field of K, in which every ideal of the original field became principal -- however the Hilbert class field itself will not necessarily have class number 1. The class field tower problem asked whether iteratively taking Hilbert class fields must stabilize after finitely many steps. In 1964, it was finally answered in the negative by Golod and Shafarevich who produced infinitely many examples and pioneered the framework that is still the most common setting for deciding when a number field will have an infinite class field tower.

In this talk, I will sketch the proof of their cohomological result and explain how it settled the class field tower problem.

In this talk, I will present our recent progress collaborated with Prof. Gui-Qiang G. Chen and Prof. Paolo Secchi on two kinds of characteristic discontinuities: relativistic vortex sheets in three-dimensional Minkowski spacetime and multi-dimensional thermoelastic contact discontinuities.
 

A group is said to satisfy the Tits Alternative if its finitely generated subgroups exhibit a striking dichotomy: they are either "big" (they contain a non-abelian free subgroup) or "small" (they are virtually soluble). Many groups of geometric interest have been shown to satisfy the Tits Alternative: linear groups, mapping class groups of hyperbolic surfaces, etc. In this talk, I will explain how one can use ideas from group actions in negative curvature to prove such a dichotomy. In particular, I will show how one can prove a strengthening of the Tits Alternative for a large class of Artin groups. This is joint work with Piotr Przytycki.

Neural networks are undoubtedly successful in practical applications. However complete mathematical theory of why and when machine learning algorithms based on neural networks work has been elusive. Although various representation theorems ensures the existence of the ``perfect’’ parameters of the network, it has not been proved that these perfect parameters can be (efficiently) approximated by conventional algorithms, such as the stochastic gradient descent. This problem is well known, since the arising optimisation problem is non-convex. In this talk we show how the optimization problem becomes convex in the mean field limit for one-hidden layer networks and certain deep neural networks. Moreover we present optimality criteria for the distribution of the network parameters and show that the nonlinear Langevin dynamics converges to this optimal distribution. This is joint work with Kaitong Hu, Zhenjie Ren and Lukasz Szpruch. 

Partition Lie algebras are generalisations of rational differential graded Lie algebras which, by a recent result of Mathew and myself, govern the formal deformation theory of algebro-geometric objects in finite and mixed characteristic. In this talk, we will take a closer look at these new gadgets and discuss some of their applications in algebra and topology

Supergravity backgrounds are an essential ingredient in string theory or field theories via AdS/CFT. The simplest example of a 4d, N=1 background is the product of four-dimensional Minkowski space with a seven-dimensional manifold with G_2 holonomy in M-theory. For more complicated backgrounds where we allow non-zero fluxes, the supersymmetry conditions can be rephrased in terms of G-structure data. The geometry of these backgrounds is often complicated and their general features are not well understood.

In this talk, I will define the analogue of G_2 geometry for generic 4d, N=1 backgrounds with flux in both type II and eleven-dimensional supergravity. The geometry is characterised by a G-structure in 'exceptional generalised geometry' that includes G_2 structures and Hitchin's generalised geometry as subcases. Supersymmetry is then equivalent to integrability of the structures, which appears as an involutivity condition and a moment map for diffeomorphisms and gauge transformations. I will show how this works in a few simple examples and discuss how this can be used to understand general properties of supersymmetric backgrounds.

 

A panel discussion on non-academic careers for mathematicians with PhDs, featuring Katia Babbar (AI Wealth Technologies & QuantBright), Jara Imbers (Risk Management Solutions) and Tom Hawes (Smith Institute).
 

Cardiac fibrosis plays a significant role in the disruption of healthy electrical signalling in the heart, creating structural heterogeneities that induce and stabilise arrhythmia.  However, a proper understanding of the consequences of cardiac fibrosis must take into account the complex and highly variable patterns of its spatial localisation in the heart, which significantly affects the extent and manner of its impacts on cardiac wave propagation. In this work we present a methodology for the algorithmic generation of fibrotic patterns via Perlin noise, a technique for computationally efficient generation of textures in computer graphics.

Our approach works directly from image data to create populations of pattern realisations that all resemble the target image under a set of metrics. Our technique thus serves as a type of data enrichment, enabling analysis of how variability in the precise placement of fibrotic structures modulates their electrophysiological impact. We demonstrate our method, and the types of analysis it can enable, using a widely referenced histological image of four different types of microfibrotic structure. Our generator and Bayesian tuning method prove flexible enough to successfully capture each of these very distinct patterns.

We demonstrate the importance of this tool, by presenting 2D simulations overlayed on the generated images that highlight the effects of microscopic variability on the electrophysiological impact of fibrosis. Finally, we discuss the application of our methodology to the increasingly available imaging data of fibrotic patterning on a more macroscopic scale, and indeed to other areas of science underpinned by image based modelling and simulation.    

Dimension reduction is an overarching theme in data science: we enjoy finding informative patterns, features or substructures in large, complex data sets. Within the field of network science, an important problem of this nature is to identify core-periphery structure. Given a network, our task is to assign each node to either the core or periphery. Core nodes should be strongly connected across the whole network whereas peripheral nodes should be strongly connected only to core nodes. More generally, we may wish to assign a non-negative value to each node, with a larger value indicating greater "coreness." This type of problem is related to, but distinct from, commumnity detection (finding clusters) and centrality assignment (finding key players), and it arises naturally in the study of networks in social science and finance. We derive and analyse a new iterative algorithm for detecting network core-periphery structure.

Using techniques in nonlinear Perron-Frobenius theory we prove global convergence to the unique solution of a relaxed version of a natural discrete optimization problem. On sparse networks, the cost of each iteration scales linearly with the number of nodes, making the algorithm feasible for large-scale problems. We give an alternative interpretation of the algorithm from the perspective of maximum likelihood reordering of a new logistic core--periphery random graph model. This viewpoint also gives a new basis for quantitatively judging a core--periphery detection algorithm. We illustrate the algorithm on a range of synthetic and real networks, and show that it offers advantages over the current state-of-the-art.

This is joint work with Francesco Tudisco (Strathclyde)

In Formula 1 engineers strive to produce the fastest car possible for their drivers. A lap simulation provides an objective evaluation of the performance of the car and the subsequent lap time achieved. Using this information, engineers aim to test new car concepts, determine performance limitations or compromises, and identify the sensitivity of performance to car setup parameters.

The latest state of the art lap simulation techniques use optimal control approaches. Optimisation methods are employed to derive the optimal control inputs of the car that achieve the fastest lap time within the constraints of the system. The resulting state trajectories define the complete behaviour of the car. Such approaches aim to create more robust, realistic and powerful simulation output compared to traditional methods.

In this talk we discuss our latest work in this area. A dynamic vehicle model is used within a free-trajectory solver based on direct optimal control methods. We discuss the reasons behind our design choices, our progress to date, and the issues we have faced during development. Further, we look at the short and long term aims of our project and how we wish to develop our mathematical methods in the future.

I will discuss the arithmetic significance of Fourier expansions of modular forms at cusps. I will talk about joint work with F. Brunault, where we determine the number field generated by Fourier coefficients of newforms at a cusp. Then I will discuss joint work with A. Saha and K. Česnavičius where we find denominator bounds for Fourier expansions at cusps and apply these bounds to a conjecture on the Manin constants of elliptic curves.