Researcher profile

Jonathan Mattingly

Jonathan Mattingly contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 15 - UnverifiedVerification L1Unclaimed author
3works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

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

Published work

3 published item(s)

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.

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.