Past

The classical Turán problem asks: given a graph $H$, how many edges can an $4n$-vertex graph have while containing no isomorphic copy of $H$? By viewing $(k+1)$-uniform hypergraphs as $k$-dimensional simplicial complexes, we can ask a topological version (first posed by Nati Linial): given a $k$-dimensional simplicial complex $S$, how many facets can an $n$-vertex $k$-dimensional simplicial complex have while containing no homeomorphic copy of $S$? Until recently, little was known for $k > 2$. In this talk, we give an answer for general $k$, by way of dependent random choice and the combinatorial notion of a trace-bounded hypergraph. Joint work with Jason Long and Bhargav Narayanan.

Non-Hermitian random matrices with complex eigenvalues represent a truly two-dimensional (2D) Coulomb gas at inverse temperature beta=2. Compared to their Hermitian counter-parts they enjoy an enlarged bulk and edge universality. As an application to ecology we model large scale data of the approximately 2D distribution of buzzard nests in the Teutoburger forest observed over a period of 20 y. These birds of prey show a highly territorial behaviour. Their occupied nests are monitored annually and we compare these data with a one-component 2D Coulomb gas of repelling charges as a function of beta. The nearest neighbour spacing distribution of the nests is well described by fitting to beta as an effective repulsion parameter, that lies between the universal predictions of Poisson (beta=0) and random matrix statistics (beta=2). Using a time moving average and comparing with next-to-nearest neighbours we examine the effect of a population increase on beta and correlation length.
 

Often in scattering applications it is advantageous to reformulate the problem as an integral equation, discretize, and then solve the resulting linear system using a fast direct solver. The computational cost of this approach is typically dominated by the work needed to compress the coefficient matrix into a rank-structured format. In this talk, we present a novel technique which exploits the bandlimited-nature of solutions to the Helmholtz equation in order to accelerate this procedure in environments where multiple frequencies are of interest.

This talk is based on joint work with Gunnar Martinsson (UT Austin).

-

A link for this talk will be sent to our mailing list a day or two in advance.  If you are not on the list and wish to be sent a link, please contact @email.

Coadmissible modules over Frechet-Stein algebras arise naturally in p-adic representation theory, e.g. in the study of locally analytic representations of p-adic Lie groups or the function spaces of rigid analytic Stein spaces. We show that in many cases, the category of coadmissible modules admits an exact and fully faithful embedding into the category of complete bornological modules, also preserving tensor products. This allows us to introduce derived methods to the study of coadmissible modules without forsaking the analytic flavour of the theory. As an application, we introduce six functors for Ardakov-Wadsley's D-cap-modules and discuss some instances where coadmissibility (in a derived sense) is preserved.

Self-similar fragmentation processes are random models for particles that are subject to successive fragmentations. When the index of self-similarity is negative the fragmentations intensify as the masses of particles decrease. This leads to a shattering phenomenon, where the initial particle is entirely reduced to dust - a set of zero-mass particles - in finite time which is what we call the extinction time. Equivalently, these extinction times may be seen as heights of continuous compact rooted trees or scaling limits of heights of sequences of discrete trees. Our objective is to set up precise bounds for the large time asymptotics of the tail distributions of these extinction times.

We consider a nonlinear flow on simplicial complexes related to the simplicial Laplacian, and show that it is a generalization of various consensus and synchronization models commonly studied on networks. In particular, our model allows us to formulate flows on simplices of any dimension, so that it includes edge flows, triangle flows, etc. We show that the system can be represented as the gradient flow of an energy functional, and use this to deduce the stability of various steady states of the model. Finally, we demonstrate that our model contains higher-dimensional analogues of structures seen in related network models.

arXiv link: https://arxiv.org/abs/2010.07421

Motivated by the work of K.R. Rajagopal, the objective of the talk is to discuss the construction and analysis of numerical approximations to a class of models that fall outside the realm of classical Cauchy elasticity. The models under consideration are implicit and nonlinear, and are referred to as strain-limiting, because the linearised strain remains bounded even when the stress is very large, a property that cannot be guaranteed within the framework of classical elastic or nonlinear elastic models. Strain-limiting models can be used to describe, for example, the behavior of brittle materials in the vicinity of fracture tips, or elastic materials in the neighborhood of concentrated loads where there is concentration of stress even though the magnitude of the strain tensor is limited.

We construct a finite element approximation of a strain-limiting elastic model and discuss the theoretical difficulties that arise in proving the convergence of the numerical method. The analytical results are illustrated by numerical experiments.

The talk is based on joint work with Andrea Bonito (Texas A&M University) and Vivette Girault (Sorbonne Université, Paris).

--

A link for this talk will be sent to our mailing list a day or two in advance.  If you are not on the list and wish to be sent a link, please contact @email.

Coronary heart disease is characterised by the formation of plaque on artery walls, restricting blood flow. If a plaque deposit ruptures, blood clot formation (thrombosis) rapidly occurs with the potential to fatally occlude the vessel within minutes. Von Willebrand Factor (vWF) is a shear-sensitive protein which has a critical role in blood clot formation in arteries. At the high shear rates typical in arterial constrictions (stenoses), vWF undergoes a conformation change, unfolding and exposing binding sites and facilitating rapid platelet deposition. 

To understand the effect of  stenosis geometry and blood flow conditions on the unfolding of vWF and subsequent platelet binding, we developed a continuum model for the initiation of thrombus formation by vWF in an idealised arterial stenosis. In this talk I will discuss modelling proteins in flow using viscoelastic fluid models, the insight asymptotic reductions can offer into this complex system and some of the challenges of studying fast arterial blood flows. 

The propagation of high-frequency electromagnetic waves can be analyzed using the geometrical optics approximation. In the case of large but finite frequencies, the geometrical optics approximation is no longer accurate, and polarization-dependent corrections at first order in wavelength modify the propagation of light in an inhomogenous medium via a spin-orbit coupling mechanism. This effect, known as the spin Hall effect of light, has been experimentally observed. In this talk I will discuss recent work which generalizes the spin Hall effect to the propagation of light and gravitational waves in inhomogenous spacetimes. This is based on joint work with Marius Oancea and Jeremie Joudioux.

In 1966 Chen Jingrun showed that every large even integer can be written as the sum of two primes or the sum of a prime and a semiprime. To date, this weakened version of Goldbach's conjecture is one of the most remarkable results of sieve theory. I will talk about the big ideas which paved the way to this proof and the ingenious trick which led to Chen's success. No prior knowledge of sieve theory required – all necessary techniques will be introduced in the talk.

The binormal flow is a model for the dynamics of a vortex filament in a 3-D inviscid incompressible fluid. This flow is also related to the classical continuous Heisenberg model in ferromagnetism and to the 1-D cubic Schrödinger equation. In this lecture I will first talk about the state of the art of the binormal flow conjecture, as well as about mathematical methods and results for the binormal flow. Then I will introduce a class of solutions at the critical level of regularity that generate singularities in finite time and describe some of their properties. These results are joint work with Luis Vega.

We discuss a novel backward Ito-Ventzell formula and an extension of the Aleeksev-Gröbner interpolating formula to stochastic flows. We also present some natural spectral conditions that yield direct and simple proofs of time uniform estimates of the difference between the two stochastic flows when their drift and diffusion functions are not the same, yielding what seems to be the first results of this type for this class of  anticipative models.

We illustrate the impact of these results in the context of diffusion perturbation theory, interacting diffusions and discrete time approximations.

In rational conformal field theory, the sewing and factorization properties are probably the most important properties that conformal blocks satisfy. For special examples such as Weiss-Zumino-Witten models and minimal models, these two properties were proved decades ago (assuming the genus is ≤1 for the sewing theorem). But for general (strongly) rational vertex operator algebras (VOAs), their proofs were finished only very recently. In this talk, I will first motivate the definition of conformal blocks and VOAs using Segal's picture of CFT. I will then explain the importance of Sewing and Factorization Theorem in the construction of full rational conformal field theory.

The torus T of projective space also acts on the Hilbert
scheme of subschemes of projective space, and the T-graph of the
Hilbert scheme has vertices the fixed points of this action, and edges
the closures of one-dimensional orbits. In general this graph depends
on the underlying field. I will discuss joint work with Rob
Silversmith, in which we construct of a subgraph, which we call the
spine, of the T-graph of Hilb^N(A^2) that is independent of the choice
of field. The key technique is an understanding of the tropical ideal,
in the sense of tropical scheme theory, of the ideal of the universal
family of an edge in the spine.

Recent developments on the black hole information paradox have shown that Euclidean wormholes — so called “replica wormholes’’  — can dominate the von Neumann entropy as computed by a gravitational path integral, and that inclusion of these wormholes results in a unitary Page curve. This development raises some puzzles from the perspective of factorization, and has raised questions regarding what the gravitational path integral is computing. In this talk, I will focus on understanding the relationship between the gravitational path integral and the partition function via the gravitational free energy (more generally the generating functional). A proper computation of the free energy requires a replica trick distinct from the usual one used to compute the entropy. I will show that in JT gravity there is a regime where the free energy computed without replica wormholes is pathological. Interestingly, the inclusion of replica wormholes is not quite sufficient to resolve the pathology: an alternative analytic continuation is required. I will discuss the implications of this for various interpretations of the gravitational path integral (e.g. as computing an ensemble average) and also mention some parallels with spin glasses. 

 I will be reviewing our recent article arXiv:2101.02218 where we propose a simple method for detecting global (a.k.a. non-perturbative) anomalies for generic quantum field theories. The basic idea is to study how the symmetries are realized on the Hilbert space of the theory. I will present several elementary examples where everything can be solved explicitly. After that, we will use these results to make non-trivial predictions about strongly interacting theories.

In this session, William Durham from the Bank of England will give a presentation about working as a mathematician in the BoE, and will give advice on interviewing for non-academic jobs. He has previously provided mock interviews in our department for jobs aimed at mathematicians with PhDs, and is happy to conduct some mock interviews (remotely, of course) for individuals as well.

Please email Helen McGregor (@email) by Monday 22 February if you might be interested in having a mock interview with William Durham on 5 March.
 

For $G$ a finite group, $p$ a prime and $(K, \mathcal{O}_K, k)$ a $p$-modular system the group ring $\mathcal{O}_K G$ is an $\mathcal{O}_k$-order in the $K$-algebra $KG.$ Graduated $\mathcal{O}_K$-orders are a particularly nice class of $\mathcal{O}_K$-orders first introduced by Zassenhaus. In this talk will see that an $\mathcal{O}_K$-order $\Lambda$ in a split $K$-algebra $A$ is graduated if the decomposition numbers for the regular $A$-module are no greater than $1$. Furthermore will see that graduated orders can be described (not uniquely) by a tuple $n$ and a matrix $M$ called the exponant matrix. Finding a suitable $n$ and $M$ for a graduated order $\Lambda$ in the $K$-algebra $A$ provides a parameterisation of the $\Lambda$-lattices inside the regular $A$-module. Understanding the $\mathcal{O}_K G$-lattices inside representations of certain groups $G$ is of interest to those involved in the Langlands programme as well as of independent interest to algebraists.

We propose a mathematical model that unifies the psychiatric concepts of drug-induced incentive salience (IST), reward prediction error

(RPE) and opponent process theory (OPT) to describe the emergence of addiction within substance abuse. The biphasic reward response (initially

positive, then negative) of the OPT is activated by a drug-induced dopamine release, and evolves according to neuro-adaptative brain

processes.  Successive drug intakes enhance the negative component of the reward response, which the user compensates for by increasing the

drug dose.  Further neuroadaptive processes ensue, creating a positive feedback between physiological changes and user-controlled drug

intake. Our drug response model can give rise to qualitatively different pathways for an initially naive user to become fully addicted.  The

path to addiction is represented by trajectories in parameter space that depend on the RPE, drug intake, and neuroadaptive changes.

We will discuss how our model can be used to guide detoxification protocols using auxiliary substances such as methadone, to mitigate withdrawal symptoms.

If this is useful here are my co-authors:
Davide Maestrini, Tom Chou, Maria R. D'Orsogna

Low rank orthogonal tensor approximation (LROTA) is an important problem in tensor computations and their applications. A classical and widely used algorithm is the alternating polar decomposition method (APD). In this talk, I will first give very a brief introduction to tensors and their decompositions. After that, an improved version named iAPD of the classical APD will be proposed and all the following four fundamental properties of iAPD will be discussed : (i) the algorithm converges globally and the whole sequence converges to a KKT point without any assumption; (ii) it exhibits an overall sublinear convergence with an explicit rate which is sharper than the usual O(1/k) for first order methods in optimization; (iii) more importantly, it converges R-linearly for a generic tensor without any assumption; (iv) for almost all LROTA problems, iAPD reduces to APD after finitely many iterations if it converges to a local minimizer. If time permits, I will also present some numerical experiments.

Our understanding and ability to compute the solutions to nonlinear partial differential equations has been strongly curtailed by our inability to effectively parameterize the inertial manifold of their solutions.  I will discuss our ongoing efforts for using machine learning to advance the state of the art, both for developing a qualitative understanding of "turbulent" solutions and for efficient computational approaches.  We aim to learn parameterizations of the solutions that give more insight into the dynamics and/or increase computational efficiency. I will discuss our recent work using machine learning to develop models of the small scale behavior of spatio-temporal complex solutions, with the goal of maintaining accuracy albeit at a highly reduced computational cost relative to a full simulation.  References: https://www.pnas.org/content/116/31/15344 and https://arxiv.org/pdf/2102.01010.pdf 

We explore reinforcement learning methods for finding the optimal policy in the linear quadratic regulator (LQR) problem. In particular, we consider the convergence of policy gradient methods in the setting of known and unknown parameters. We are able to produce a global linear convergence guarantee for this approach in the setting of finite time horizon and stochastic state dynamics under weak assumptions. The convergence of a projected policy gradient method is also established in order to handle problems with constraints. We illustrate the performance of the algorithm with two examples. The first example is the optimal liquidation of a holding in an asset. We show results for the case where we assume a model for the underlying dynamics and where we apply the method to the data directly. The empirical evidence suggests that the policy gradient method can learn the global optimal solution for a larger class of stochastic systems containing the LQR framework and that it is more robust with respect to model mis-specification when compared to a model-based approach. The second example is an LQR system in a higher-dimensional setting with synthetic data.

This talk concerns applications of differential geometry in numerical optimization. They arise when the optimization problem can be formulated as finding an optimum of a real-valued cost function defined on a smooth nonlinear search space. Oftentimes, the search space is a "matrix manifold", in the sense that its points admit natural representations in the form of matrices. In most cases, the matrix manifold structure is due either to the presence of certain nonlinear constraints (such as orthogonality or rank constraints), or to invariance properties in the cost function that need to be factored out in order to obtain a nondegenerate optimization problem. Manifolds that come up in practical applications include the rotation group SO(3) (generation of rigid body motions from sample points), the set of fixed-rank matrices (low-rank models, e.g., in collaborative filtering), the set of 3x3 symmetric positive-definite matrices (interpolation of diffusion tensors), and the shape manifold (morphing).

In the recent years, the practical importance of optimization problems on manifolds has stimulated the development of geometric optimization algorithms that exploit the differential structure of the manifold search space. In this talk, we give an overview of geometric optimization algorithms and their applications, with an emphasis on the underlying geometric concepts and on the numerical efficiency of the algorithm implementations.

A link for this talk will be sent to our mailing list a day or two in advance.  If you are not on the list and wish to be sent a link, please contact @email.

The pandemic has had a deleterious influence on the Hollywood film
industry.  Fortunately,  however, the thin film industry continues to
flourish.  A host of effects are responsible for confined liquids
exhibiting properties that differ from their bulk counterparts. For
example, the dominant polarization and surface forces across a layered
system can control the material behavior on length scales vastly larger
than the film thickness.  This basic class of phenomena, wherein
volume-volume interactions create large pressures, are at play in,
amongst many other settings, wetting, biomaterials, ceramics, colloids,
and tribology.  When the films so created involve phase change and are
present in disequilibrium, the forces can be so large that they destroy
the setting that allowed them to form in the first place. I will
describe the connection between such films in a semi-traditional wetting
dynamics geometry and active brownian dynamics.  I then explore their
power to explain a wide range of processes from materials- to astro- to
geo-science.