Source author record

Jonathan Mattingly

Jonathan Mattingly appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

math.PR Computational Geometry math-ph math.AP math.AT math.MP

Catalog footprint

What is connected

4works

6topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2022arXiv

Convergence of Stratified MCMC Sampling of Non-Reversible Dynamics

We present a form of stratified MCMC algorithm built with non-reversible stochastic dynamics in mind. It can also be viewed as a generalization of the exact milestoning method, or form of NEUS. We prove convergence of the method under certain assumptions, with expressions for the convergence rate in terms of the process's behavior within each stratum and large scale behavior between strata. We show that the algorithm has a unique fixed point which corresponds to the invariant measure of the process without stratification. We will show how the speeds of two versions of the new algorithm, one with an extra eigenvalue problem step and one without, relate to the mixing rate of a discrete process on the strata, and the mixing probability of the process being sampled within each stratum. The eigenvalue problem version also relates to local and global perturbation results of discrete Markov chains, such as those given by Van Koten, Weare et. al.

preprint2020arXiv

Limiting Behaviors of High Dimensional Stochastic Spin Ensembles

Lattice spin models in statistical physics are used to understand magnetism. Their Hamiltonians are a discrete form of a version of a Dirichlet energy, signifying a relationship to the Harmonic map heat flow equation. The Gibbs distribution, defined with this Hamiltonian, is used in the Metropolis-Hastings (M-H) algorithm to generate dynamics tending towards an equilibrium state. In the limiting situation when the inverse temperature is large, we establish the relationship between the discrete M-H dynamics and the continuous Harmonic map heat flow associated with the Hamiltonian. We show the convergence of the M-H dynamics to the Harmonic map heat flow equation in two steps: First, with fixed lattice size and proper choice of proposal size in one M-H step, the M-H dynamics acts as gradient descent and will be shown to converge to a system of Langevin stochastic differential equations (SDE). Second, with proper scaling of the inverse temperature in the Gibbs distribution and taking the lattice size to infinity, it will be shown that this SDE system converges to the deterministic Harmonic map heat flow equation. Our results are not unexpected, but show remarkable connections between the M-H steps and the SDE Stratonovich formulation, as well as reveal trajectory-wise out of equilibrium dynamics to be related to a canonical PDE system with geometric constraints.

preprint2014arXiv

Probabilistic Fréchet Means for Time Varying Persistence Diagrams

In order to use persistence diagrams as a true statistical tool, it would be very useful to have a good notion of mean and variance for a set of diagrams. In 2011, Mileyko and his collaborators made the first study of the properties of the Fréchet mean in $(\mathcal{D}_p,W_p)$, the space of persistence diagrams equipped with the p-th Wasserstein metric. In particular, they showed that the Fréchet mean of a finite set of diagrams always exists, but is not necessarily unique. The means of a continuously-varying set of diagrams do not themselves (necessarily) vary continuously, which presents obvious problems when trying to extend the Fréchet mean definition to the realm of vineyards. We fix this problem by altering the original definition of Fréchet mean so that it now becomes a probability measure on the set of persistence diagrams; in a nutshell, the mean of a set of diagrams will be a weighted sum of atomic measures, where each atom is itself a persistence diagram determined using a perturbation of the input diagrams. This definition gives for each $N$ a map $(\mathcal{D}_p)^N \to \mathbb{P}(\mathcal{D}_p)$. We show that this map is Hölder continuous on finite diagrams and thus can be used to build a useful statistic on time-varying persistence diagrams, better known as vineyards.

preprint2008arXiv

A stochastic-Lagrangian particle system for the Navier-Stokes equations

This paper is based on a formulation of the Navier-Stokes equations developed by P. Constantin and the first author (\texttt{arxiv:math.PR/0511067}, to appear), where the velocity field of a viscous incompressible fluid is written as the expected value of a stochastic process. In this paper, we take $N$ copies of the above process (each based on independent Wiener processes), and replace the expected value with $\frac{1}{N}$ times the sum over these $N$ copies. (We remark that our formulation requires one to keep track of $N$ stochastic flows of diffeomorphisms, and not just the motion of $N$ particles.) We prove that in two dimensions, this system of interacting diffeomorphisms has (time) global solutions with initial data in the space $\holderspace{1}α$ which consists of differentiable functions whose first derivative is $α$ Hölder continuous (see Section \ref{sGexist} for the precise definition). Further, we show that as $N \to \infty$ the system converges to the solution of Navier-Stokes equations on any finite interval $[0,T]$. However for fixed $N$, we prove that this system retains roughly $O(\frac{1}{N})$ times its original energy as $t \to \infty$. Hence the limit $N \to \infty$ and $T\to \infty$ do not commute. For general flows, we only provide a lower bound to this effect. In the special case of shear flows, we compute the behaviour as $t \to \infty$ explicitly.