Imperial College

We introduce constrained variational principles for fluid dynamics on rough paths. The advection of the fluid is constrained to be the sum of a vector field which represents coarse-scale motion and a rough (in time) vector field which parametrizes fine-scale motion. The rough vector field is regarded as fixed and the rough partial differential equation for the coarse-scale velocity is derived as a consequence of being a critical point of the action functional.

The action functional is perturbative in the sense that if the rough vector f ield is set to zero, then the corresponding variational principle agrees with the reduced (to the vector fields) Euler-Poincare variational principle introduced in Holm, Marsden and Ratiu (1998). More precisely, the Lagrangian in the action functional encodes the physics of the fluid and is a function of only the coarse-scale velocity. 

By parametrizing the fine-scales of fluid motion with a rough vector field, we preserve the pathwise nature of deterministic fluid dynamics and establish a flexible framework for stochastic parametrization schemes. The main benefit afforded by our approach is that the system of rough partial differential equations we derive satisfy essential conservation laws, including Kelvin’s circulation theorem. This talk is based on recent joint work with Dan Crisan, Darryl Holm, and Torstein Nilssen.

The notion of emergence is at the core of many of the most challenging open scientific questions, being so much a cause of wonder as a perennial source of philosophical headaches. Two classes of emergent phenomena are usually distinguished: strong emergence, which corresponds to supervenient properties with irreducible causal power; and weak emergence, which are properties generated by the lower levels in such "complicated" ways that they can only be derived by exhaustive simulation. While weak emergence is generally accepted, a large portion of the scientific community considers causal emergence to be either impossible, logically inconsistent, or scientifically irrelevant.

In this talk we present a novel, quantitative framework that assesses emergence by studying the high-order interactions of the system's dynamics. By leveraging the Integrated Information Decomposition (ΦID) framework [1], our approach distinguishes two types of emergent phenomena: downward causation, where macroscopic variables determine the future of microscopic degrees of freedom; and causal decoupling, where macroscopic variables influence other macroscopic variables without affecting their corresponding microscopic constituents. Our framework also provides practical tools that are applicable on a range of scenarios of practical interest, enabling to test -- and possibly reject -- hypotheses about emergence in a data-driven fashion. We illustrate our findings by discussing minimal examples of emergent behaviour, and present a few case studies of systems with emergent dynamics, including Conway’s Game of Life, neural population coding, and flocking models.
[1] Mediano, Pedro AM, Fernando Rosas, Robin L. Carhart-Harris, Anil K. Seth, and Adam B. Barrett. "Beyond integrated information: A taxonomy of information dynamics phenomena." arXiv preprint arXiv:1909.02297 (2019).
 

Tautological bundles on Hilbert schemes of points often enter into enumerative and physical computations. I will explain how to use the Donaldson-Thomas theory of toric threefolds to produce combinatorial identities that are expressed geometrically using tautological bundles on the Hilbert scheme of points on a surface. I'll also explain how these identities can be used to study Euler characteristics of tautological bundles over Hilbert schemes of points on general surfaces.

The Asai (or twisted tensor) L-function attached to a Bianchi modular form is the 'restriction to the rationals' of the standard L-function. Introduced by Asai in 1977, subsequent study has linked its special values to the arithmetic of the corresponding form. In this talk, I will discuss joint work with David Loeffler in which we construct a p-adic Asai L-function -- that is, a measure on Z_p* that interpolates the critical values L^As(f,chi,1) -- for ordinary weight 2 Bianchi modular forms. We use a new method for constructing p-adic L-functions, using Kato's system of Siegel units to build a 'Betti analogue' of an Euler system, building on algebraicity results of Ghate. I will start by giving a brief introduction to p-adic L-functions and Bianchi modular forms, and if time permits, I will briefly mention another case where the method should apply, that of non-self-dual automorphic representations for GL(3).

By viewing a stochastic process as a random variable taking values in a path space, the support of its law describes the set of all attainable paths. In this talk, we show that the support of the law of a solution to a path-dependent stochastic differential equation is given by the image of the Cameron-Martin space under the flow of mild solutions to path-dependent ordinary differential equations, constructed by means of the vertical derivative of the diffusion coefficient. This result is based on joint work with Rama Cont and extends the Stroock-Varadhan support theorem for diffusion processes to the path-dependent case.

I will give an overview of some recent progress on potential automorphy results over CM fields, that is joint work with Allen, Calegari, Gee, Helm, Le Hung, Newton, Scholze, Taylor, and Thorne. I will focus on explaining an application to the generalized Ramanujan-Petersson conjecture. 

Convex and in particular semidefinite relaxations (SDP) for non-convex continuous quadratic optimisation can provide tighter bounds than traditional linear relaxations. However, using SDP relaxations directly in Branch&Cut is impeded by lack of warm starting and inefficiency when combined with other cut classes, i.e. the reformulation-linearization technique. We present a general framework based on machine learning for a strong linear outer-approximation that can retain most tightness of such SDP relaxations, in the form of few strong low dimensional linear cuts selected offline. The cut selection complexity is taken offline by using a neural network estimator (trained before installing solver software) as a selection device for the strongest cuts. Lastly, we present results of our method on QP/QCQP problem instances.

We analyse under which conditions equilibration between two competing effects, repulsion modelled by nonlinear diffusion and attraction modelled by nonlocal interaction, occurs. I will discuss several regimes that appear in aggregation diffusion problems with homogeneous kernels. I will first concentrate in the fair competition case distinguishing among porous medium like cases and fast diffusion like ones. I will discuss the main qualitative properties in terms of stationary states and minimizers of the free energies. In particular, all the porous medium cases are critical while the fast diffusion are not. In the second part, I will discuss the diffusion dominated case in which this balance leads to continuous compactly supported radially decreasing equilibrium configurations for all masses. All stationary states with suitable regularity are shown to be radially symmetric by means of continuous Steiner symmetrisation techniques. Calculus of variations tools allow us to show the existence of global minimizers among these equilibria. Finally, in the particular case of Newtonian interaction in two dimensions they lead to uniqueness of equilibria for any given mass up to translation and to the convergence of solutions of the associated nonlinear aggregation-diffusion equations towards this unique equilibrium profile up to translations as time tends to infinity. This talk is based on works in collaboration with S. Hittmeir, B. Volzone and Y. Yao and with V. Calvez and F. Hoffmann.

In this talk we address the numerical approximation of Mean Field Games with local couplings. For finite difference discretizations of the Mean Field Game system, we follow a variational approach, proving that the schemes can be obtained as the optimality system of suitably defined optimization problems. In order to prove the existence of solutions of the scheme with a variational argument, the monotonicity of the coupling term is not used, which allow us to recover general existence results. Next, assuming next that the coupling term is monotone, the variational problem is cast as a convex optimization problem for which we study and compare several proximal type methods. These algorithms have several interesting features, such as global convergence and stability with respect to the viscosity parameter. We conclude by presenting numerical experiments assessing the performance of the proposed methods. In collaboration with L. Briceno-Arias (Valparaiso, CL) and F. J. Silva (Limoges, FR).