Related papers: Sampling Biased Monotonic Surfaces using Exponenti…
Sampling from binary quadratic distributions (BQDs) is a fundamental but challenging problem in discrete optimization and probabilistic inference. Previous work established theoretical guarantees for stochastic localization (SL) in…
In this paper, we consider inverse time-harmonic acoustic and electromagnetic scattering from locally perturbed rough surfaces in three dimensions. The scattering interface is supposed to be the graph of a Lipschitz continuous function with…
For Markov chain Monte Carlo methods, one of the greatest discrepancies between theory and system is the scan order - while most theoretical development on the mixing time analysis deals with random updates, real-world systems are…
Mixing in viscous fluids is challenging, but chaotic advection in principle allows efficient mixing. In the best possible scenario,the decay rate of the concentration profile of a passive scalar should be exponential in time. In practice,…
Mixing of materials is fundamental to many natural phenomena and engineering applications. The presence of discontinuous deformations - such as shear banding or wall slip - creates new mechanisms for mixing and transport beyond those…
Markov chains are one of the well-known tools for modeling and analyzing stochastic systems. At the same time, they are used for constructing random walks that can achieve a given stationary distribution. This paper is concerned with…
Consider the barycentric subdivision which cuts a given triangle along its medians to produce six new triangles. Uniformly choosing one of them and iterating this procedure gives rise to a Markov chain. We show that almost surely, the…
A mixed basis approach based on density functional theory is employed for low dimensional systems. The basis functions are taken to be plane waves for the periodic direction multiplied by B-spline polynomials in the non-periodic direction.…
Uniform sampling of binary matrix with fixed margins is an important and difficult problem in statistics, computer science, ecology and so on. The well-known swap algorithm would be inefficient when the size of the matrix becomes large or…
We call "flippable tilings" of a constant curvature surface a tiling by "black" and "white" faces, so that each edge is adjacent to two black and two white faces (one of each on each side), the black face is forward on the right side and…
We consider the problem of generating uniformly random partitions of the vertex set of a graph such that every piece induces a connected subgraph. For the case where we want to have partitions with linearly many pieces of bounded size, we…
High-dimensional distributions, especially those with heavy tails, are notoriously difficult for off-the-shelf MCMC samplers: the combination of unbounded state spaces, diminishing gradient information, and local moves results in…
We investigate and prove the mathematical properties of a general class of one-dimensional unimodal smooth maps perturbed with a heteroscedastic noise. Specifically, we investigate the stability of the associated Markov chain, show the weak…
We present a novel approach to quantizing Markov chains. The approach is based on the Markov chain coupling method, which is frequently used to prove fast mixing. Given a particular coupling, e.g., a grand coupling, we construct a…
We use molecular dynamics (MD) to simulate an unstable homogeneous mixture of binary fluids (AB), confined in a slit pore of width $D$. The pore walls are assumed to be flat and structureless, and attract one component of the mixture (A)…
It is common practice in Markov chain Monte Carlo to update the simulation one variable (or sub-block of variables) at a time, rather than conduct a single full-dimensional update. When it is possible to draw from each full-conditional…
Sampling permutations from S_n is a fundamental problem from probability theory. The nearest neighbor transposition chain \cal{M}}_{nn} is known to converge in time \Theta(n^3 \log n) in the uniform case and time \Theta(n^2) in the constant…
We prove that square-tiled surfaces having fixed combinatorics of horizontal cylinder decomposition and tiled with smaller and smaller squares become asymptotically equidistributed in any ambient linear $GL(\mathbb R)$-invariant suborbifold…
We consider the advection equation on $\mathbb{T}^2$ with a real analytic and time-periodic velocity field that alternates between two Hamiltonian shears. Randomness is injected by alternating the vector field randomly in time between just…
A key task in Bayesian statistics is sampling from distributions that are only specified up to a partition function (i.e., constant of proportionality). However, without any assumptions, sampling (even approximately) can be #P-hard, and few…