Past events

Causal optimal transport and the related adapted Wasserstein distance have recently been popularized as a more appropriate alternative to the classical Wasserstein distance in the context of stochastic analysis and mathematical finance. In this talk, we establish some interesting consequences of causality for transports on the space of continuous functions between the laws of stochastic differential equations.
 

We first characterize bicausal transport plans and maps between the laws of stochastic differential equations. As an application, we are able to provide necessary and sufficient conditions for bicausal transport plans to be induced by bi-causal maps. Analogous to the classical case, we show that bicausal Monge transports are dense in the set of bicausal couplings between laws of SDEs with unique strong solutions and regular coefficients.

 This is a joint work with Rama Cont.

This talk explores recent advancements in stress and flux-based finite element methods. It focuses on addressing the limitations of traditional finite elements, in order to describe complex material behavior and engineer new metamaterials.

Stress and flux-based finite element methods are particularly useful in error estimation, laying the groundwork for adaptive refinement strategies. This concept builds upon the hypercircle theorem [1], which states that in a specific energy space, both the exact solution and any admissible stress field lie on a hypercircle. However, the construction of finite element spaces that satisfy admissible states for complex material behavior is not straightforward. It often requires a relaxation of specific properties, especially when dealing with non-symmetric stress tensors [2] or hyperelastic materials.

Alternatively, methods that directly approximate stresses can be employed, offering high accuracy of the stress fields and adherence to physical conservation laws. However, when approximating eigenvalues, this significant benefit for the solution's accuracy implies that the solution operator cannot be compact. To address this, the solution operator must be confined to a subset of the solution that excludes the stresses. Yet, due to compatibility conditions, the trial space for the other solution components typically does not yield the desired accuracy. The second part of this talk will therefore explore the Least-Squares method as a remedy to these challenges [3].

To conclude this talk, we will emphasize the integration of those methods within global solution strategies, with a particular focus on the challenges regarding model order reduction methods [4].

[1] W. Prager, J. Synge. Approximations in elasticity based on the concept of function space.

Quarterly of Applied Mathematics 5(3), 1947.

[2] FB, K. Bernhard, M. Moldenhauer, G. Starke. Weakly symmetric stress equilibration and a posteriori error estimation for linear elasticity, Numerical Methods for Partial Differential Equations 37(4), 2021.

[3] FB, D. Boffi. First order least-squares formulations for eigenvalue problems, IMA Journal of Numerical Analysis 42(2), 2023.

[4] FB, D. Boffi, A. Halim. A reduced order model for the finite element approximation of eigenvalue problems,Computer Methods in Applied Mechanics and Engineering 404, 2023.

We study “bottom-up” models for the collective motion of large groups of animals. Similar models can be encoded into (micro)robotic matter, capable of sensing light and processing information. Agents are endowed only with visual sensing and information processing. We study a model in which moving agents reorientate to maximise the path-entropy of their visual environment over paths into the future. There are general arguments that principles like this that are based on retaining freedom in the future may confer fitness in an uncertain world. Alternative “top-down” models are more common in the literature. These typically encode coalignment and/or cohesion directly and are often motivated by models drawn from physics, e.g. describing spin systems. However, such models can usually give little insight into how co-alignment and cohesion emerge because these properties are encoded in the model at the outset, in a top-down manner. We discuss how our model leads to dynamics with striking similarities with animal systems, including the emergence of coalignment, cohesion, a characteristic density scaling anddifferent behavioural phenotypes. The dynamics also supports a very unusual order-disorder transition in which the order (coalignment) initially increases upon the addition of sensory or behavioural noise, before decreasing as the noise becomes larger.

This is a pre-seminar meeting for Margaret Bilu's talk "A motivic circle method", which takes place later in the day at 5PM in L3.

I will present novel freeness results in ultraproducts of tracial von Neumann algebras. As a particular case, I will show that if a and b are the generators of the free group F_2, then the relative commutants of a and b in the ultraproduct of the free group factor are free with respect to the ultraproduct trace. The proof is based on a surprising application of Lp-boundedness results of Fourier multipliers in free group factors for p > 2. I will describe applications of these results to absorption and model theory of II_1 factors. This is joint work with Adrian Ioana.

We will look at a number of interesting examples — some proven, others merely conjectured — of Hamburger moment sequences in combinatorics. We will consider ways in which this positivity may be expected: for instance, in different types of combinatorial statistics on perfect matchings that encode moments of noncommutative analogues of the classical Central Limit Theorem. We will also consider situations in which this positivity may be surprising, and where proving it would open up new approaches to a class of very hard open problems in combinatorics.

Random temporal graphs are a version of the classical Erdős-Rényi random graph $G(n,p)$ where additionally, each edge has a distinct random time stamp, and connectivity is constrained to sequences of edges with increasing time stamps. We are interested in the asymptotics for the distances in such graphs, mostly in the regime of interest where the average degree $np$ is of order $\log n$ ('near' the phase transition).

More specifically, we will discuss the asymptotic lengths of increasing paths: the lengths of the shortest and longest paths between typical vertices, as well as the maxima between any two vertices; this also covers the (temporal) diameter. In the regime $np \gg \log n$, longest increasing paths were studied by Angel, Ferber, Sudakov and Tassion.

The talk contains joint work with Nicolas Broutin and Gábor Lugosi.

The classification of measure preserving actions up to orbit equivalence has attracted a lot of interest in the last 25 years. The goal of this talk is to survey the major discoveries in the field, including Popa's cocycle and orbit equivalence superrigidity theorem and discuss some recent superrigidity results for dense subgroups of Lie groups acting by translation.

Optimization problems constrained to a smooth manifold can be solved via the framework of Riemannian optimization. To that end, a geometrical description of the constraining manifold, e.g., tangent spaces, retractions, and cost function gradients, is required. In this talk, we present a novel approach that allows performing approximate Riemannian optimization based on a manifold learning technique, in cases where only a noiseless sample set of the cost function and the manifold’s intrinsic dimension are available.

In this talk we will review the problem of sharpness in percolation, tracing its origins back to the seminal works of Menshikov, Grimmett-Marstrand and Aizenman-Barsky, which successfully settled the question in the context of Bernoulli independent percolation. Then we will present some recent advancements on the field, which have opened up the possibility of investigating dependent percolation models. Special emphasis will be given to the Interpolation technique, which has proved itself very effective. In particular, it has been used to establish the sharpness for Interlacements Percolation, a model introduced by Sznitman with notoriously intricate dependencies.

This talk is based on a joint work with Duminil-Copin, Goswami, Rodriguez and Severo

Let G be a reductive group defined over a local field of characteristic 0 (real or p-adic). By Harish-Chandra’s regularity theorem, the character Θ_π of an irreducible representation π of G is given by a locally integrable function f_π on G. It turns out that f_π has even better integrability properties, namely, it is locally L^{1+r}-integrable for some r>0. This gives rise to a new singularity invariant of representations \e_π by considering the largest such r.

We explore \e_π, show it is bounded below only in terms of the group G, and calculate it in the case of a p-adic GL(n). To do so, we relate \e_π to the integrability of Fourier transforms of nilpotent orbital integrals appearing in the local character expansion of Θ_π. As a main technical tool, we use explicit resolutions of singularities of certain hyperplane arrangements. We obtain bounds on the multiplicities of K-types in irreducible representations of G for a p-adic G and a compact open subgroup K.

Based on a joint work with Itay Glazer and Julia Gordon.

Gaussian processes provide a powerful probabilistic kernel learning framework, which allows high-quality nonparametric learning via methods such as Gaussian process regression. Nevertheless, its learning phase requires unrealistic massive computations for large datasets. In this talk, we present a quadrature-based approach for scaling up Gaussian process regression via a low-rank approximation of the kernel matrix. The low-rank structure is utilized to achieve effective hyperparameter learning, training, and prediction. Our Gauss-Legendre features method is inspired by the well-known random Fourier features approach, which also builds low-rank approximations via numerical integration. However, our method is capable of generating high-quality kernel approximation using a number of features that is poly-logarithmic in the number of training points, while similar guarantees will require an amount that is at the very least linear in the number of training points when using random Fourier features. The utility of our method for learning with low-dimensional datasets is demonstrated using numerical experiments.

The Bethe-Gauge Correspondence (BGC) of Nekrasov and Shatashvili, linking quantum integrable spin chains to two-dimensional supersymmetric gauge theories with N=2 supersymmetry, stands out as a significant instance of the deep connection between supersymmetric gauge theories and integrable models. In this talk, I will delve into this correspondence and its origins for superspin chains. To achieve this, I will first elucidate the Bethe Side and its corresponding Gauge Side of the BGC. Subsequently, it becomes evident that the BGC can be naturally realized within String Theory. I will initially outline the brane configuration for the realization of the Gauge Side. Through the use of string dualities, this brane configuration will be mapped to another, embodying the Bethe Side of the correspondence. The 4D Chern-Simons theory plays a crucial role in this latter duality frame, elucidating the integrability of the Bethe Side. Lastly, I will elaborate on computing the main object of interest for integrable superspin chains—the R-matrix—from the BGC. While this provides a rather comprehensive picture of the correspondence, some important questions remain for further clarification. I will summarize some of the most interesting ones at the end of the talk.

In a recent preprint, together with Bailleul and Chevyrev we introduced a class of random fields which try to model the basic properties of quantum fields. I will try to explain the basic ideas and some of the many open problems.

To read the preprint, please click here.

Busenberg and Travis suggested in 1983 a population system that exhibits complete segregation of the species. This system can be rigorously derived from interacting particle systems in a mean-field-type limit. It consists of parabolic cross-diffusion equations with an indefinite diffusion matrix. It is known that this system can be formulated in terms of so-called entropy variables such that the transformed equations possess a positive semidefinite diffusion matrix. We consider in this talk the case of incomplete diffusion, which means that the diffusion matrix has zero eigenvalues, and the problem is not parabolic in the sense of Petrovskii. 

We show that the cross-diffusion equations can be written as a normal form of symmetric hyperbolic-parabolic type beyond the Kawashima-Shizuta theory. Using results for symmetric hyperbolic systems, we prove the existence of a unique local classical solution. As solutions may become discontinuous in finite time, only global solutions with very low regularity can be expected. We prove the existence of global dissipative measure-valued solutions satisfying a weak-strong uniqueness property. The proof is based on entropy methods and a finite-volume approximation with a mesh-dependent artificial diffusion. 

Many congruences between modular forms (or at least their q-expansions) can be explained by the theory of $p$-adic families of modular forms. In this talk, I will discuss properties of eigenvarieties, a geometric interpretation of the idea of $p$-adic families. In particular, focusing initially on the well-understood case of (elliptic) modular forms, before delving into the considerably murkier world of Bianchi modular forms. In this second case, this work gives numerical verification of a couple of conjectures, including BSD by work of Loeffler and Zerbes.

We propose a novel generative model for time series based on Schrödinger bridge (SB) approach. This consists in the entropic interpolation via optimal transport between a reference probability measure on path space and a target measure consistent with the joint data distribution of the time series. The solution is characterized by a stochastic differential equation on finite horizon with a path-dependent drift function, hence respecting  the temporal dynamics of the time series distribution. We  estimate the drift function from data samples by nonparametric, e.g. kernel regression methods,  and the simulation of the SB diffusion  yields new synthetic data samples of the time series. The performance of our generative model is evaluated through a series of numerical experiments.  First, we test with autoregressive models, a GARCH Model, and the example of fractional Brownian motion,  and measure the accuracy of our algorithm with marginal, temporal dependencies metrics, and predictive scores. Next, we use our SB generated synthetic samples for the application to deep hedging on real-data sets. 

Let $(M, g)$ be a Riemannian manifold equipped with a calibration $k$-form $\alpha$. In earlier work with Cheng and Madnick (AJM 2021), we studied the analytic properties of a special class of $k$-harmonic maps into $M$ satisfying a first order nonlinear PDE, whose images (away from a critical set) are $\alpha$-calibrated submanifolds of $M$. We call these maps Smith immersions, as they were originally introduced in an unpublished preprint of Aaron Smith. They have nice properties related to conformal geometry, and are higher-dimensional analogues of the $J$-holomorphic map equation. In new joint work (arXiv:2311.14074) with my PhD student Anton Iliashenko, we have obtained analogous results for maps out of $M$. Slightly more precisely, we define a special class of $k$-harmonic maps out of $M$, satisfying a first order nonlinear PDE, whose fibres (away from a critical set) are $\alpha$-calibrated submanifolds of $M$. We call these maps Smith submersions. I will give an introduction to both of these sets of equations, and discuss many future questions.

Mathematical methods are developed to characterize the asymptotics of recurrent neural networks (RNN) as the number of hidden units, data samples in the sequence, hidden state updates, and training steps simultaneously grow to infinity. In the case of an RNN with a simplified weight matrix, we prove the convergence of the RNN to the solution of an infinite-dimensional ODE coupled with the fixed point of a random algebraic equation. 
The analysis requires addressing several challenges which are unique to RNNs. In typical mean-field applications (e.g., feedforward neural networks), discrete updates are of magnitude O(1/N ) and the number of updates is O(N). Therefore, the system can be represented as an Euler approximation of an appropriate ODE/PDE, which it will converge to as N → ∞. However, the RNN hidden layer updates are O(1). Therefore, RNNs cannot be represented as a discretization of an ODE/PDE and standard mean-field techniques cannot be applied. Instead, we develop a fixed point analysis for the evolution of the RNN memory state, with convergence estimates in terms of the number of update steps and the number of hidden units. The RNN hidden layer is studied as a function in a Sobolev space, whose evolution is governed by the data sequence (a Markov chain), the parameter updates, and its dependence on the RNN hidden layer at the previous time step. Due to the strong correlation between updates, a Poisson equation must be used to bound the fluctuations of the RNN around its limit equation. These mathematical methods allow us to prove a neural tangent kernel (NTK) limit for RNNs trained on data sequences as the number of data samples and size of the neural network grow to infinity.

Say we want to view the function-Rips multifiltration as an estimator. Then, what is the target? And what kind of consistency, bias, or convergence rate, should we expect? In this talk I will present on-going joint work with Ethan André (Ecole Normale Supérieure) that aims at laying the algebro-topological ground to start answering these questions.

Motile cells inside living tissues often encounter junctions, where their path branches into several alternative directions of migration. We present a theoretical model of cellular polarization for cells migrating along one-dimensional lines, exhibiting diverse migration modes. When arriving at a symmetric Y-junction and extending protrusions along the different paths that emanate from the junction. The model predicts the spontaneous emergence of deterministic oscillations between competing protrusions, whereby the cellular polarization and growth alternates between the competing protrusions. These predicted oscillations are found experimentally for two different cell types, noncancerous endothelial and cancerous glioma cells, migrating on patterned network of thin adhesive lanes with junctions. Finally we present an analysis of the migration modes of multicellular "trains" along one-dimensional tracks.

The Junior Algebra and Representation Theory Seminar will kick-off the start of Hilary term with a social event in the common room. Come to catch up with your fellow students and maybe play a board game or two. Afterwards we'll have lunch together.