BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Jonathan C. Mattingly

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

ResearcherUnclaimed source record
math.PRmath.DSmath-phmath.MPmath.NAmath.STStatistics TheoryComputationmath.APmath.MGcond-mat.mtrl-scimath.CANumerical Analysisphysics.flu-dynphysics.soc-ph

Catalog footprint

What is connected

24works
15topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

Research graph

See the researcher in context

Open full explorer

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.

Published work

24 published item(s)

preprint2023arXiv

Gibbsian dynamics and the generalized Langevin equation

We study the statistically invariant structures of the nonlinear generalized Langevin equation (GLE) with a power-law memory kernel. For a broad class of memory kernels, including those in the subdiffusive regime, we construct solutions of the GLE using a Gibbsian framework, which does not rely on existing Markovian approximations. Moreover, we provide conditions on the decay of the memory to ensure uniqueness of statistically steady states, generalizing previous known results for the GLE under particular kernels as a sum of exponentials.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Multi-Scale Merge-Split Markov Chain Monte Carlo for Redistricting

We develop a Multi-Scale Merge-Split Markov chain on redistricting plans. The chain is designed to be usable as the proposal in a Markov Chain Monte Carlo (MCMC) algorithm. Sampling the space of plans amounts to dividing a graph into a partition with a specified number of elements which each correspond to a different district. The districts satisfy a collection of hard constraints and the measure may be weighted with regard to a number of other criteria. The multi-scale algorithm is similar to our previously developed Merge-Split proposal, however, this algorithm provides improved scaling properties and may also be used to preserve nested communities of interest such as counties and precincts. Both works use a proposal which extends the ReCom algorithm which leveraged spanning trees merge and split districts. In this work we extend the state space so that each district is defined by a hierarchy of trees. In this sense, the proposal step in both algorithms can be seen as a "Forest ReCom." We also expand the state space to include edges that link specified districts, which further improves the computational efficiency of our algorithm. The collection of plans sampled by the MCMC algorithm can serve as a baseline against which a particular plan of interest is compared. If a given plan has different racial or partisan qualities than what is typical of the collection of plans, the given plan may have been gerrymandered and is labeled as an outlier.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Non-reversible Markov chain Monte Carlo for sampling of districting maps

Evaluating the degree of partisan districting (Gerrymandering) in a statistical framework typically requires an ensemble of districting plans which are drawn from a prescribed probability distribution that adheres to a realistic and non-partisan criteria. In this article we introduce novel non-reversible Markov chain Monte-Carlo (MCMC) methods for the sampling of such districting plans which have improved mixing properties in comparison to previously used (reversible) MCMC algorithms. In doing so we extend the current framework for construction of non-reversible Markov chains on discrete sampling spaces by considering a generalization of skew detailed balance. We provide a detailed description of the proposed algorithms and evaluate their performance in numerical experiments.

0 citations0 reviewsNo code linkedOpen access
preprint2020arXiv

Singularities of invariant densities for random switching between two linear ODEs in 2D

We consider a planar dynamical system generated by two stable linear vector fields with distinct fixed points and random switching between them. We characterize singularities of the invariant density in terms of the switching rates and contraction rates. We prove boundedness away from those singularities. We also discuss some motivating biological examples.

0 citations0 reviewsNo code linkedOpen access
preprint2019arXiv

Geometric ergodicity of Langevin dynamics with Coulomb interactions

This paper is concerned with the long time behavior of Langevin dynamics of {\em Coulomb gases} in $\mathbf{R}^d$ with $d\geq 2$, that is a second order system of Brownian particles driven by an external force and a pairwise repulsive Coulomb force. We prove that the system converges exponentially to the unique Boltzmann-Gibbs invariant measure under a weighted total variation distance. The proof relies on a novel construction of Lyapunov function for the Coulomb system.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Sensitivity to switching rates in stochastically switched ODEs

We consider a stochastic process driven by a linear ordinary differential equation whose right-hand side switches at exponential times between a collection of different matrices. We construct planar examples that switch between two matrices where the individual matrices and the average of the two matrices are all Hurwitz (all eigenvalues have strictly negative real part), but nonetheless the process goes to infinity at large time for certain values of the switching rate. We further construct examples in higher dimensions where again the two individual matrices and their averages are all Hurwitz, but the process has arbitrarily many transitions between going to zero and going to infinity at large time as the switching rate varies. In order to construct these examples, we first prove in general that if each of the individual matrices is Hurwitz, then the process goes to zero at large time for sufficiently slow switching rate and if the average matrix is Hurwitz, then the process goes to zero at large time for sufficiently fast switching rate. We also give simple conditions that ensure the process goes to zero at large time for all switching rates.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Stochastic switching in infinite dimensions with applications to random parabolic PDEs

We consider parabolic PDEs with randomly switching boundary conditions. In order to analyze these random PDEs, we consider more general stochastic hybrid systems and prove convergence to, and properties of, a stationary distribution. Applying these general results to the heat equation with randomly switching boundary conditions, we find explicit formulae for various statistics of the solution and obtain almost sure results about its regularity and structure. These results are of particular interest for biological applications as well as for their significant departure from behavior seen in PDEs forced by disparate Gaussian noise. Our general results also have applications to other types of stochastic hybrid systems, such as ODEs with randomly switching right-hand sides.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Noise-Induced Stabilization of Planar Flows I

We show that the complex-valued ODE \begin{equation*} \dot z_t = a_{n+1} z^{n+1} + a_n z^n+\cdots+a_0, \end{equation*} which necessarily has trajectories along which the dynamics blows up in finite time, can be stabilized by the addition of an arbitrarily small elliptic, additive Brownian stochastic term. We also show that the stochastic perturbation has a unique invariant measure which is heavy-tailed yet is uniformly, exponentially attracting. The methods turn on the construction of Lyapunov functions. The techniques used in the construction are general and can likely be used in other settings where a Lyapunov function is needed. This is a two-part paper. This paper, Part I, focuses on general Lyapunov methods as applied to a special, simplified version of the problem. Part II of this paper extends the main results to the general setting.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Noise-Induced Stabilization of Planar Flows II

We continue the work started in Part I of this article, showing how the addition of noise can stabilize an otherwise unstable system. The analysis makes use of nearly optimal Lyapunov functions. In this continuation, we remove the main limiting assumption of Part I by an inductive procedure as well as establish a lower bound which shows that our construction is radially sharp. We also prove a version of Peskir's \cite{Peskir_07} generalized Tanaka formula adapted to patching together Lyapunov functions. This greatly simplifies the analysis used in previous works.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

On Unique Ergodicity in Nonlinear Stochastic Partial Differential Equations

We illustrate how the notion of asymptotic coupling provides a flexible and intuitive framework for proving the uniqueness of invariant measures for a variety of stochastic partial differential equations whose deterministic counterpart possesses a finite number of determining modes. Examples exhibiting parabolic and hyperbolic structure are studied in detail. In the later situation we also present a simple framework for establishing the existence of invariant measures when the usual approach relying on the Krylov-Bogolyubov procedure and compactness fails.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Sticky central limit theorems at isolated hyperbolic planar singularities

We derive the limiting distribution of the barycenter $b_n$ of an i.i.d. sample of $n$ random points on a planar cone with angular spread larger than $2π$. There are three mutually exclusive possibilities: (i) (fully sticky case) after a finite random time the barycenter is almost surely at the origin; (ii) (partly sticky case) the limiting distribution of $\sqrt{n} b_n$ comprises a point mass at the origin, an open sector of a Gaussian, and the projection of a Gaussian to the sector's bounding rays; or (iii) (nonsticky case) the barycenter stays away from the origin and the renormalized fluctuations have a fully supported limit distribution---usually Gaussian but not always. We conclude with an alternative, topological definition of stickiness that generalizes readily to measures on general metric spaces.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Invariant measure selection by noise: An Example

We consider a deterministic system with two conserved quantities and infinity many invariant measures. However the systems possess a unique invariant measure when enough stochastic forcing and balancing dissipation are added. We then show that as the forcing and dissipation are removed a unique limit of the deterministic system is selected. The exact structure of the limiting measure depends on the specifics of the stochastic forcing.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Redistricting and the Will of the People

We introduce a non-partisan probability distribution on congressional redistricting of North Carolina which emphasizes the equal partition of the population and the compactness of districts. When random districts are drawn and the results of the 2012 election were re-tabulated under the drawn districtings, we find that an average of 7.6 democratic representatives are elected. 95% of the randomly sampled redistrictings produced between 6 and 9 Democrats. Both of these facts are in stark contrast with the 4 Democrats elected in the 2012 elections with the same vote counts. This brings into serious question the idea that such elections represent the "will of the people." It underlines the ability of redistricting to undermine the democratic process, while on the face allowing democracy to proceed.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Sticky central limit theorems on open books

Given a probability distribution on an open book (a metric space obtained by gluing a disjoint union of copies of a half-space along their boundary hyperplanes), we define a precise concept of when the Fréchet mean (barycenter) is sticky. This nonclassical phenomenon is quantified by a law of large numbers (LLN) stating that the empirical mean eventually almost surely lies on the (codimension $1$ and hence measure $0$) spine that is the glued hyperplane, and a central limit theorem (CLT) stating that the limiting distribution is Gaussian and supported on the spine. We also state versions of the LLN and CLT for the cases where the mean is nonsticky (i.e., not lying on the spine) and partly sticky (i.e., is, on the spine but not sticky).

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Diffusion limits of the random walk Metropolis algorithm in high dimensions

Diffusion limits of MCMC methods in high dimensions provide a useful theoretical tool for studying computational complexity. In particular, they lead directly to precise estimates of the number of steps required to explore the target measure, in stationarity, as a function of the dimension of the state space. However, to date such results have mainly been proved for target measures with a product structure, severely limiting their applicability. The purpose of this paper is to study diffusion limits for a class of naturally occurring high-dimensional measures found from the approximation of measures on a Hilbert space which are absolutely continuous with respect to a Gaussian reference measure. The diffusion limit of a random walk Metropolis algorithm to an infinite-dimensional Hilbert space valued SDE (or SPDE) is proved, facilitating understanding of the computational complexity of the algorithm.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Geometric ergodicity of a bead-spring pair with stochastic Stokes forcing

We consider a simple model for the fluctuating hydrodynamics of a flexible polymer in dilute solution, demonstrating geometric ergodicity for a pair of particles that interact with each other through a nonlinear spring potential while being advected by a stochastic Stokes fluid velocity field. This is a generalization of previous models which have used linear spring forces as well as white-in-time fluid velocity fields. We follow previous work combining control theoretic arguments, Lyapunov functions, and hypo-elliptic diffusion theory to prove exponential convergence via a Harris chain argument. In addition we allow the possibility of excluding certain "bad" sets in phase space in which the assumptions are violated but from which the system leaves with a controllable probability. This allows for the treatment of singular drifts, such as those derived from the Lennard-Jones potential, which is a novel feature of this work.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Propagating Lyapunov Functions to Prove Noise--induced Stabilization

We investigate an example of noise-induced stabilization in the plane that was also considered in (Gawedzki, Herzog, Wehr 2010) and (Birrell, Herzog, Wehr 2011). We show that despite the deterministic system not being globally stable, the addition of additive noise in the vertical direction leads to a unique invariant probability measure to which the system converges at a uniform, exponential rate. These facts are established primarily through the construction of a Lyapunov function which we generate as the solution to a sequence of Poisson equations. Unlike a number of other works, however, our Lyapunov function is constructed in a systematic way, and we present a meta-algorithm we hope will be applicable to other problems. We conclude by proving positivity properties of the transition density by using Malliavin calculus via some unusually explicit calculations.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

A Theory of Hypoellipticity and Unique Ergodicity for Semilinear Stochastic PDEs

We present a theory of hypoellipticity and unique ergodicity for semilinear parabolic stochastic PDEs with "polynomial" nonlinearities and additive noise, considered as abstract evolution equations in some Hilbert space. It is shown that if Hormander's bracket condition holds at every point of this Hilbert space, then a lower bound on the Malliavin covariance operator $M_t$ can be obtained. Informally, this bound can be read as "Fix any finite-dimensional projection $Π$ on a subspace of sufficiently regular functions. Then the eigenfunctions of $M_t$ with small eigenvalues have only a very small component in the image of $Π$". We also show how to use a priori bounds on the solutions to the equation to obtain good control on the dependency of the bounds on the Malliavin matrix on the initial condition. These bounds are sufficient in many cases to obtain the asymptotic strong Feller property introduced in [HairerMattingly06]. One of the main novel technical tools is an almost sure bound from below on the size of "Wiener polynomials", where the coefficients are possibly non-adapted stochastic processes satisfying a Lipschitz condition. By exploiting the polynomial structure of the equations, this result can be used to replace Norris' lemma, which is unavailable in the present context. We conclude by showing that the two-dimensional stochastic Navier-Stokes equations and a large class of reaction-diffusion equations fit the framework of our theory.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

A weak trapezoidal method for a class of stochastic differential equations

We present a numerical method for the approximation of solutions for the class of stochastic differential equations driven by Brownian motions which induce stochastic variation in fixed directions. This class of equations arises naturally in the study of population processes and chemical reaction kinetics. We show that the method constructs paths that are second order accurate in the weak sense. The method is simpler than many second order methods in that it neither requires the construction of iterated Ito integrals nor the evaluation of any derivatives. The method consists of two steps. In the first an explicit Euler step is used to take a fractional step. This fractional point is then combined with the initial point to obtain a higher order, trapezoidal like, approximation. The higher order of accuracy stems from the fact that both the drift and the quadratic variation of the underlying SDE are approximated to second order.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Convergence of Numerical Time-Averaging and Stationary Measures via Poisson Equations

Numerical approximation of the long time behavior of a stochastic differential equation (SDE) is considered. Error estimates for time-averaging estimators are obtained and then used to show that the stationary behavior of the numerical method converges to that of the SDE. The error analysis is based on using an associated Poisson equation for the underlying SDE. The main advantage of this approach is its simplicity and universality. It works equally well for a range of explicit and implicit schemes including those with simple simulation of random variables, and for hypoelliptic SDEs. To simplify the exposition, we consider only the case where the state space of the SDE is a torus and we study only smooth test functions. However we anticipate that the approach can be applied more widely. An analogy between our approach and Stein's method is indicated. Some practical implications of the results are discussed.

0 citations0 reviewsNo code linkedOpen access
preprint2010arXiv

Rare Transition Events in Nonequilibrium Systems with State-Dependent Noise: Application to Stochastic Current Switching in Semiconductor Superlattices

Using recent mathematical advances, a geometric approach to rare noise-driven transition events in nonequilibrium systems is given, and an algorithm for computing the maximum likelihood transition curve is generalized to the case of state-dependent noise. It is applied to a model of electronic transport in semiconductor superlattices to investigate transitions between metastable electric field distributions. When the applied voltage $V$ is varied near a saddle-node bifurcation at $V_th$, the mean life time $<T>$ of the initial metastable state is shown to scale like $log<T> \propto |V_th - V|^{3/2}$ as $V\to V_th$.

0 citations0 reviewsNo code linkedOpen access
preprint2005arXiv

Stationary Solutions of Stochastic Differential Equation with Memory and Stochastic Partial Differential Equations

We explore Ito stochastic differential equations where the drift term possibly depends on the infinite past. Assuming the existence of a Lyapunov function, we prove the existence of a stationary solution assuming only minimal continuity of the coefficients. Uniqueness of the stationary solution is proven if the dependence on the past decays sufficiently fast. The results of this paper are then applied to stochastically forced dissipative partial differential equations such as the stochastic Navier-Stokes equation and stochastic Ginsburg-Landau equation.

0 citations0 reviewsNo code linkedOpen access