University of Sussex

The Wasserstein metric originates in the theory of optimal transport, and among many other applications, it provides a natural way to measure how evenly distributed a finite point set is. We give a survey of classical and more recent results that describe the behaviour of some random point processes in Wasserstein metric, including the eigenvalues of some random matrix models, and explain the connection to the logarithm of the characteristic polynomial of a random unitary matrix. We also discuss a simple random walk model on the unit circle defined in terms of a quadratic irrational number, which turns out to be related to surprisingly deep arithmetic properties of real quadratic fields.

In this talk, I review some recent results obtained for black holes using
effective field theory methods applied to quantum gravity, in particular the
unique effective action. Black holes are complex thermodynamical objects
that not only have a temperature but also have a pressure. Furthermore, they
have quantum hair which provides a solution to the black hole information
paradox.

We study expectations of powers and correlations for characteristic polynomials of N x N non-Hermitian random matrices. This problem is related to the analysis of planar models (log-gases) where a Gaussian (or other) background measure is perturbed by a finite number of point charges in the plane. I will discuss the critical asymptotics, for example when a point charge collides with the boundary of the support, or when two point charges collide with each other (coalesce) in the bulk. In many of these situations, we are able to express the results in terms of Painlevé transcendents. The application to certain d-fold rotationally invariant models will be discussed. This is joint work with Alfredo Deaño (University of Kent).

We analyze a fully discrete numerical scheme for solving a parabolic PDE on a moving surface. The method is based on a diffuse interface approach that involves a level set description of the moving surface. Under suitable conditions on the spatial grid size, the time step and the interface width we obtain stability and error bounds with respect to natural norms. Test calculations are presented that confirm our analysis.

In this talk, I will present a novel solution strategy to efficiently and accurately compute approximate solutions to semilinear optimal control problems, focusing on the optimal control of phase field formulations of geometric evolution laws.
The optimal control of geometric evolution laws arises in a number of applications in fields including material science, image processing, tumour growth and cell motility.
Despite this, many open problems remain in the analysis and approximation of such problems.
In the current work we focus on a phase field formulation of the optimal control problem, hence exploiting the well developed mathematical theory for the optimal control of semilinear parabolic partial differential equations.
Approximation of the resulting optimal control problem is computationally challenging, requiring massive amounts of computational time and memory storage.
The main focus of this work is to propose, derive, implement and test an efficient solution method for such problems. The solver for the discretised partial differential equations is based upon a geometric multigrid method incorporating advanced techniques to deal with the nonlinearities in the problem and utilising adaptive mesh refinement.
An in-house two-grid solution strategy for the forward and adjoint problems, that significantly reduces memory requirements and CPU time, is proposed and investigated computationally.
Furthermore, parallelisation as well as an adaptive-step gradient update for the control are employed to further improve efficiency.
Along with a detailed description of our proposed solution method together with its implementation we present a number of computational results that demonstrate and evaluate our algorithms with respect to accuracy and efficiency.
A highlight of the present work is simulation results on the optimal control of phase field formulations of geometric evolution laws in 3-D which would be computationally infeasible without the solution strategies proposed in the present work.

Linear elasticity can be rigorously derived from finite elasticity in the case of small loadings in terms of \Gamma-convergence. This was first done by Dal Maso-Negri-Percivale in the case of one-well energies with super-quadratic growth. This has been later generalised to different settings, in particular to the case of multi-well energies where the distance between the wells is very small (comparable to the size of the load). I will discuss recent developments in the case when the distance between the wells is arbitrary. In this context linear elasticity can be derived by adding to the multi-well energy a singular higher order term which penalises jumps from one well to another. The size of the singular term has to satisfy certain scaling assumptions which turn out to be optimal. (This is joint work with Alicandro, Dal Maso and Lazzaroni.) 

In this seminar I will present a semi-langrangian discretisation of the Monge-Ampère operator, which is of interest in optimal transport 
and differential geometry as well as in related fields of application.

I will discuss the proof of convergence to viscosity solutions. To address the challenge of uniqueness and convexity we draw upon the classical relationship with Hamilton-Jacobi-Bellman equations, which we extend to the viscosity setting. I will explain that the monotonicity of semi-langrangian schemes implies that they possess large stencils, which in turn requires careful treatment of the boundary conditions.

The contents of the seminar is based on current work with X Feng from the University of Tennessee.