Past Partial Differential Equations Seminar

Systems of high-contrast resonators can be used to control and manipulate wave-matter interactions at scales that are much smaller than the operating wavelengths. The aim of this talk is to review recent studies of ordered and disordered systems of subwavelength resonators and to explain some of their topologically protected localization properties. Both reciprocal and non-reciprocal systems will be considered.
 

In the first half of the talk I will review the theory of nuclear magnetic resonance (NMR), leading to the Bloch-Torrey PDE. I will then describe the pulsed-gradient spin-echo method for measuring the Fourier transform of the voxel-averaged propagator of the Bloch-Torrey equation.  This technique permits one to compute the diffusion coefficient in a voxel. For complex biological tissue, as in the brain, the standard model represents spin-echo as a multiexponential signal, whose exponents and coefficients describe the diffusion coefficients and volume fractions of isolated tissue compartments, respectively. The question of identifying these parameters from experimental measurements leads us to investigate the degree of well-posedness of this problem that I will discuss in the second half of the talk. We show that the parameter reconstruction problem exhibits power law transition to ill-posedness, and derive the explicit formula for the exponent by reformulating the problem in terms of the integral equation that can be solved explicitly. This is a joint work with my Ph.D. student Henry J. Brown.

In this talk we consider the four-waves spatially homogeneous kinetic equation arising in weak wave turbulence theory from the microscopic Fermi-Pasta-Ulam-Tsingou (FPUT) oscillator chains.  This equation is sometimes referred to as the Phonon Boltzmann Equation. I will discuss the global existence and stability of solutions of the kinetic equation near the Rayleigh-Jeans (RJ) thermodynamic equilibrium solutions. This is a joint work with Pierre Germain (Imperial College London) and Joonhyun La (KIAS).

The Stein Variational Gradient Descent method is a variational inference method in statistics that has recently received a lot of attention. The method provides a deterministic approximation of the target distribution, by introducing a nonlocal interaction with a kernel. Despite the significant interest, the exponential rate of convergence for the continuous method has remained an open problem, due to the difficulty of establishing the related so-called Stein-log-Sobolev inequality. Here, we prove that the inequality is satisfied for each space dimension and every kernel whose Fourier transform has a quadratic decay at infinity and is locally bounded away from zero and infinity. Moreover, we construct weak solutions to the related PDE satisfying exponential rate of decay towards the equilibrium. The main novelty in our approach is to interpret the Stein-Fisher information as a duality pairing between $H^{-1}$ and $H^{1}$, which allows us to employ the Fourier transform. We also provide several examples of kernels for which the Stein-log-Sobolev inequality fails, partially showing the necessity of our assumptions. This is a joint work with J. A. Carrillo and J. Warnett. 

When it comes to the nonlinear heat equation u_t - \Delta u = u^p, a sharp condition for the stability of positive radial steady states was derived in the classical paper by Gui, Ni and Wang.  In this talk, I will present some recent joint work with Daniel Devine that focuses on a more general system of reaction-diffusion equations (which is also also known as the parabolic Henon-Lane-Emden system).  We obtain a sharp condition that determines the stability of positive radial steady states, and we also study the separation property of these solutions along with their asymptotic behaviour at infinity.

The Camassa–Holm equation, which is nonlinear one-dimensional nonlinear PDE which is completely integrable and has  applications in several areas, has received considerable attention. We will discuss recent work regarding the Camassa—Holm equation with transport noise, more precisely, the equation $u_t+uu_x+P_x+\sigma u_x \circ dW=0$ and $P-P_{xx}=u^2+u_x^2/2$. În particular, we will show existence of a weak, global, dissipative solution of the Cauchy initial-value problem on the torus.  This is joint work with L. Galimberti (King’s College), K.H. Karlsen (Oslo), and P.H.C. Pang (NTNU/Oslo).

The Riemann problem is an IVP having simple piecewise constant initial data that is invariant under scaling. In 1D, the problem was originally considered by Riemann during the 19th century in the context of gas dynamics, and the general theory was more or less completed by Lax and Glimm in the mid-20th century. In 2D and MD, the situation is much more complicated, and very few analytic results are available. We discuss a shock reflection problem for the Euler equations for potential flow, with initial data that generates four interacting shockwaves. After reformulating the problem as a free boundary problem for a nonlinear PDE of mixed hyperbolic-elliptic type, the problem is solved via a sophisticated iteration procedure. The talk is based on joint work with G-Q Chen (Oxford) et. al. arXiv:2305.15224, to appear in JEMS (2025).

Minimal discrete energy problems arise in a variety of scientific contexts – such as crystallography, nanotechnology, information theory, and viral morphology, to name but a few.     Our goal is to analyze the structure of configurations generated by optimal (and near optimal)-point configurations that minimize the Riesz s-energy over a sphere in Euclidean space R^d and, more generally, over a bounded manifold. The Riesz s-energy potential, which is a generalization of the Coulomb potential, is simply given by 1/r^s, where r denotes the distance between pairs of points. We show how such potentials for s>d and their minimizing point configurations are ideal for use in sampling surfaces.

Connections to the results by Field's medalist M. Viazovska and her collaborators on best-packing and universal optimality in 8 and 24 dimensions will be discussed. Finally we analyze the minimization of a "k-nearest neighbor" truncated version of Riesz energy that reduces the order N^2 computation for energy minimization to order N log N , while preserving global and local properties.

In this talk I will give a short introduction to moderately interacting particle systems and the general notion of fluctuations around the mean-field limit. We will see how a central limit theorem can be shown for moderately interacting particles on the whole space for certain types of interaction potentials. The interaction potential approximates singular attractive potentials of sub-Coulomb type and we can show that the fluctuations become asymptotically Gaussians. The methodology is inspired by the classical work of Oelschläger in the 1980s on fluctuations for the porous-medium equation. To allow for attractive potentials we use a new approach of quantitative mean-field convergence in probability in order to include aggregation effects. 

I will introduce the class of causally-null-compactifiable spacetimes that can be canonically converted into compact timed-metric spaces using the cosmological time function of Andersson-Galloway-Howard and the null distance of Sormani-Vega. This class of space-times includes future developments of compact initial data sets and regions exhausting asymptotically flat space-times. I will discuss various intrinsic notions of distance between such space-times and show that some of them are definite in the sense that they are equal to zero if and only if there is a time-oriented Lorentzian isometry between the space-times. These definite distances allow us to define notions of convergence of space-times to limit space-times that are not necessarily smooth. This is joint work with Christina Sormani.

The Dyson Brownian motion (DMB) is a system of interacting Brownian motions with logarithmic interaction potential, which was introduced by Freeman Dyson '62 in relation to the random matrix theory. In this talk, we discuss the case where the number of particles is infinite and show that the DBM induces a diffusion structure on the configuration space having the Bakry-Émery lower Ricci curvature bound. As an application, we show that the DBM can be realised as the unique Benamou-Brenier-type gradient flow of the Boltzmann-Shannon entropy associated with the sine_beta point process. 

Understanding the mechanisms behind the spatial distribution, self-organisation and aggregation of organisms is a central issue in both ecology and cell biology. Since self-organisation at the population level is the cumulative effect of behaviours at the individual level, it requires a mathematical approach to be elucidated.
In nature, every individual, be it a cell or an animal, inspects its territory before moving. The process of acquiring information from the environment is typically non-local, i.e. individuals have the ability to inspect a portion of their territory. In recent years, a growing body of empirical research has shown that non-locality is a key aspect of movement processes, while mathematical models incorporating non-local interactions have received increasing attention for their ability to accurately describe how interactions between individuals and their environment can affect their movement, reproduction rate and well-being. In this talk, I will present a study of a class of advection-diffusion equations that model population movements generated by non-local species interactions. Using a combination of analytical and numerical tools, I will show that these models support a wide variety of spatio-temporal patterns that are able to reproduce segregation, aggregation and time-periodic behaviours commonly observed in real systems. I will also show the existence of parameter regions where multiple stable solutions coexist and hysteresis phenomena.
Overall, I will describe various methods for analysing bifurcations and pattern formation properties of these models, which represent an essential mathematical tool for addressing fundamental questions about the many aggregation phenomena observed in nature.

The "index of hypocoercivity" is defined via a coercivity-type estimate for the self-adjoint/skew-adjoint parts of the generator, and it quantifies `how degenerate' a hypocoercive evolution equation is, both for ODEs and for evolutions equations in a Hilbert space. We show that this index characterizes the polynomial decay of the propagator norm for short time and illustrate these concepts for the Lorentz kinetic equation on a torus. Discrete time analogues of the above systems (obtained via the mid-point rule) are contractive, but typically not strictly contractive. For this setting we introduce "hypocontractivity" and an "index of hypocontractivity" and discuss their close connection to the continuous time evolution equations.

This talk is based on joint work with F. Achleitner, E. Carlen, E. Nigsch, and V. Mehrmann.

References:
1) F. Achleitner, A. Arnold, E. Carlen, The Hypocoercivity Index for the short time behavior of linear time-invariant ODE systems, J. of Differential Equations (2023).
2) A. Arnold, B. Signorello, Optimal non-symmetric Fokker-Planck equation for the convergence to a given equilibrium, Kinetic and Related Models (2022).
3) F. Achleitner, A. Arnold, V. Mehrmann, E. Nigsch, Hypocoercivity in Hilbert spaces, J. of Functional Analysis (2025).
 

On arbitrary Carnot groups, the only hypoelliptic Hodge-Laplacians on forms that have been introduced are 0-order pseudodifferential operators constructed using the Rumin complex.  However, to address questions where one needs sharp estimates, this 0-order operator is not suitable. Indeed, this is a rather difficult problem to tackle in full generality, the main issue being that the Rumin exterior differential is not homogeneous on arbitrary Carnot groups. In this talk, I will focus on the specific example of the free Carnot group of step 3 with 2 generators, where it is possible to introduce different hypoelliptic Hodge-Laplacians on forms. Such Laplacians can be used to obtain sharp div-curl type inequalities akin to those considered by Bourgain & Brezis and Lanzani & Stein for the de Rham complex, or their subelliptic counterparts obtained by Baldi, Franchi & Pansu for the Rumin complex on Heisenberg groups

We prove the local Lipschitz continuity of energy minimizing harmonic maps between singular spaces, more specifically from the n-dimensional Heisenberg group into CAT(0) spaces. The present result paves the way for a general regularity theory of sub-elliptic harmonic maps, providing a versatile approach applicable beyond the Heisenberg group.  Joint work with Yaoting Gui and Jürgen Jost.

We propose a density functional theory of Thomas-Fermi-(von Weizsacker) type to describe the response of a single layer of graphene to a charge some distance away from the layer. We formulate a variational setting in which the proposed energy functional admits minimizers. We further provide conditions under which those minimizers are unique. The associated Euler-Lagrange equation for the charge density is also obtained, and uniqueness, regularity and decay of the minimizers are proved under general conditions. For a class of special potentials, we also establish a precise universal asymptotic decay rate, as well as an exact charge cancellation by the graphene sheet. In addition, we discuss the existence of nodal minimizers which leads to multiple local minimizers in the TFW model. This is a joint work with Cyrill Muratov (University of Pisa).

When studying interacting particle systems, two distinct categories emerge: indistinguishable systems, where particle identity does not influence system dynamics, and non-exchangeable systems, where particle identity plays a significant role. One way to conceptualize these second systems is to see them as particle systems on weighted graphs. In this talk, we focus on the latter category. Recent developments in graph theory have raised renewed interest in understanding largepopulation limits in these systems. Two main approaches have emerged: graph limits and mean-field limits. While mean-field limits were traditionally introduced for indistinguishable particles, they have been extended to the case of non-exchangeable particles recently. In this presentation, we introduce several models, mainly from the field of opinion dynamics, for which rigorous convergence results as N tends to infinity have been obtained. We also clarify the connection between the graph limit approach and the mean-field limit one. The works discussed draw from several papers, some co-authored with Nastassia Pouradier Duteil and David Poyato.

We study the mathematical properties of time-dependent flows of incompressible fluids that respond as an Euler fluid until the modulus of the symmetric part of the velocity gradient exceeds a certain, a-priori given but arbitrarily large, critical value. Once the velocity gradient exceeds this threshold, a dissipation mechanism is activated. Assuming that the fluid, after such an activation, dissipates the energy in a specific manner, we prove that the corresponding initial-boundary-value problem is well-posed in the sense of Hadamard. In particular, we show that for an arbitrary, sufficiently regular, initial velocity there is a global-in-time unique weak solution to the spatially-periodic problem. This is a joint result with Miroslav Bulíček. 

We present a general concept that is suitable for studying the stability of equilibria for open systems in continuum thermodynamics. We apply such concept to a generalized Newtonian incompressible heat conducting fluid with prescribed nonuniform temperature on the boundary and with the no-slip boundary conditions for the velocity in three dimensional domain. For large class of constitutive relation for the Cauchy stress, we identify a class of proper solutions converging to the equilibria exponentially in a suitable metric and independently of the distance to equilibria at the initial time. Consequently, the equilibrium is nonlinearly stable and attracts all weak solutions from that class. The proper solutions exist and satisfy entropy (in)equality.

We study fully discrete approximation of the 2D Euler equations for ideal homogeneous fluids. We focus on spectral methods and  discuss rates of convergence of velocity and vorticity under different assumptions on the smoothness of the data.

The curvature-dimension condition, CD(K,N) for short, and the (weaker) measure contraction property, or MCP(K,N), are two synthetic notions for a metric measure space to have Ricci curvature bounded from below by K and dimension bounded from above by N. In this talk, we investigate the validity of these conditions in sub-Finsler geometry, which is a wide generalization of Finsler and sub-Riemannian geometry. Firstly, we show that sub-Finsler manifolds equipped with a smooth strongly convex norm and with a positive smooth measure can not satisfy the CD(K,N) condition for any K and N. Secondly, we focus on the sub-Finsler Heisenberg group, where we show that, on the one hand, the CD(K,N) condition can not hold for any reference norm and, on the other hand, the MCP(K,N) may hold or fail depending on the regularity of the reference norm. 

First we will show that the continuum counterpart of the discrete individual-based mechanical model that describes the dynamics of two contiguous cell populations is given by a free-boundary problem for the cell densities.  Then, in addition to interactions, we will consider the microscopic movement of cells and derive a fractional cross-diffusion system as the many-particle limit of a multi-species system of moderately interacting particles.