BZPEERPapers, people, evidence and research pathsDiscover
Research tools
Network
Opportunities
Account
Source author record
Martin Hairer appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.
Catalog footprint
Research graph
Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.
BZPEER is loading the nearby papers, people, topics and institutions for this page.
Published work
We consider the Allen-Cahn equation $\partial_t u- Δu=u-u^3$ with a rapidly mixing Gaussian field as initial condition. We show that provided that the amplitude of the initial condition is not too large, the equation generates fronts described by nodal sets of the Bargmann-Fock Gaussian field, which then evolve according to mean curvature flow.
The past decade has seen tremendous progress in our understanding of the behaviour of many probabilistic models at or near their "critical point". On the 5th of July 2022, Hugo Duminil-Copin was awarded the Fields medal for the crucial role he played in many of these developments. In this short review article, we will try to put his work into context and present a small selection of his results.
We introduce a $1+1$-dimensional temperature-dependent model such that the classical ballistic deposition model is recovered as its zero-temperature limit. Its $\infty$-temperature version, which we refer to as the $0$-Ballistic Deposition ($0$-BD) model, is a randomly evolving interface which, surprisingly enough, does {\it not} belong to either the Edwards--Wilkinson (EW) or the Kardar--Parisi--Zhang (KPZ) universality class. We show that $0$-BD has a scaling limit, a new stochastic process that we call {\it Brownian Castle} (BC) which, although it is "free", is distinct from EW and, like any other renormalisation fixed point, is scale-invariant, in this case under the $1:1:2$ scaling (as opposed to $1:2:3$ for KPZ and $1:2:4$ for EW). In the present article, we not only derive its finite-dimensional distributions, but also provide a "global" construction of the Brownian Castle which has the advantage of highlighting the fact that it admits backward characteristics given by the (backward) Brownian Web (see [Tóth B., Werner W., Probab. Theory Related Fields, '98] and [L. R. G. Fontes, M. Isopi, C. M. Newman, and K. Ravishankar, Ann. Probab., '04]). Among others, this characterisation enables us to establish fine pathwise properties of BC and to relate these to special points of the Web. We prove that the Brownian Castle is a (strong) Markov and Feller process on a suitable space of càdlàg functions and determine its long-time behaviour. At last, we give a glimpse to its universality by proving the convergence of $0$-BD to BC in a rather strong sense.
We consider the KPZ equation in one space dimension driven by a stationary centred space-time random field, which is sufficiently integrable and mixing, but not necessarily Gaussian. We show that, in the weakly asymmetric regime, the solution to this equation considered at a suitable large scale and in a suitable reference frame converges to the Hopf-Cole solution to the KPZ equation driven by space-time Gaussian white noise. While the limiting process depends only on the integrated variance of the driving field, the diverging constants appearing in the definition of the reference frame also depend on higher order moments.
This paper studies the intermediate time behaviour of a small random perturbation of a periodic cellular flow. Our main result shows that on time scales shorter than the diffusive time scale, the limiting behaviour of trajectories that start close enough to cell boundaries is a fractional kinetic process: A Brownian motion time changed by the local time of an independent Brownian motion. Our proof uses the Freidlin-Wentzell framework, and the key step is to establish an analogous averaging principle on shorter time scales. As a consequence of our main theorem, we obtain a homogenization result for the associated advection-diffusion equation. We show that on intermediate time scales the effective equation is a fractional time PDE that arises in modelling anomalous diffusion.
Recently, a solution theory for one-dimensional stochastic PDEs of Burgers type driven by space-time white noise was developed. In particular, it was shown that natural numerical approximations of these equations converge and that their convergence rate in the uniform topology is arbitrarily close to $\frac{1}{6}$. In the present article we improve this result in the case of additive noise by proving that the optimal rate of convergence is arbitrarily close to $\frac{1}{2}$.
We review a series of forthcoming results leading to the construction of a natural evolution on the space of loops with values in a Riemannian manifold. In particular, this clarifies the algebraic structure of the renormalisation procedures appearing in the context of the theory of regularity structures.
We consider stochastic semi-linear evolution equations which are driven by additive, spatially correlated, Wiener noise, and in particular consider problems of heat equation (analytic semigroup) and damped-driven wave equations (bounded semigroup) type. We discretize these equations by means of a spectral Galerkin projection, and we study the approximation of the probability distribution of the trajectories: test functions are regular, but depend on the values of the process on the interval $[0,T]$. We introduce a new approach in the context of quantative weak error analysis for discretization of SPDEs. The weak error is formulated using a deterministic function (Itô map) of the stochastic convolution found when the nonlinear term is dropped. The regularity properties of the Itô map are exploited, and in particular second-order Taylor expansions employed, to transfer the error from spectral approximation of the stochastic convolution into the weak error of interest. We prove that the weak rate of convergence is twice the strong rate of convergence in two situations. First, we assume that the covariance operator commutes with the generator of the semigroup: the first order term in the weak error expansion cancels out thanks to an independence property. Second, we remove the commuting assumption, and extend the previous result, thanks to the analysis of a new error term depending on a commutator.
We consider a large class of $1+1$-dimensional continuous interface growth models and we show that, in both the weakly asymmetric and the intermediate disorder regimes, these models converge to Hopf-Cole solutions to the KPZ equation.
We prove a version of the Wong-Zakai theorem for one-dimensional parabolic nonlinear stochastic PDEs driven by space-time white noise. As a corollary, we obtain a detailed local description of solutions. Dedicated to the memory of Kiyosi Itô on the occasion of the 100th anniversary of his birth.
We study a Markov process with two components: the first component evolves according to one of finitely many underlying Markovian dynamics, with a choice of dynamics that changes at the jump times of the second component. The second component is discrete and its jump rates may depend on the position of the whole process. Under regularity assumptions on the jump rates and Wasserstein contraction conditions for the underlying dynamics, we provide a concrete criterion for the convergence to equilibrium in terms of Wasserstein distance. The proof is based on a coupling argument and a weak form of the Harris theorem. In particular, we obtain exponential ergodicity in situations which do not verify any hypoellipticity assumption, but are not uniformly contracting either. We also obtain a bound in total variation distance under a suitable regularising assumption. Some examples are given to illustrate our result, including a class of piecewise deterministic Markov processes.
The celebrated Hörmander condition is a sufficient (and nearly necessary) condition for a second-order linear Kolmogorov partial differential equation (PDE) with smooth coefficients to be hypoelliptic. As a consequence, the solutions of Kolmogorov PDEs are smooth at all positive times if the coefficients of the PDE are smooth and satisfy Hörmander's condition even if the initial function is only continuous but not differentiable. First-order linear Kolmogorov PDEs with smooth coefficients do not have this smoothing effect but at least preserve regularity in the sense that solutions are smooth if their initial functions are smooth. In this article, we consider the intermediate regime of nonhypoelliptic second-order Kolmogorov PDEs with smooth coefficients. The main observation of this article is that there exist counterexamples to regularity preservation in that case. More precisely, we give an example of a second-order linear Kolmogorov PDE with globally bounded and smooth coefficients and a smooth initial function with compact support such that the unique globally bounded viscosity solution of the PDE is not even locally Hölder continuous. From the perspective of probability theory, the existence of this example PDE has the consequence that there exists a stochastic differential equation (SDE) with globally bounded and smooth coefficients and a smooth function with compact support which is mapped by the corresponding transition semigroup to a function which is not locally Hölder continuous. In other words, degenerate noise can have a roughening effect. A further implication of this loss of regularity phenomenon is that numerical approximations may converge without any arbitrarily small polynomial rate of convergence to the true solution of the SDE. More precisely, we prove for an example SDE with globally bounded and smooth coefficients that the standard Euler approximations converge to the exact solution of the SDE in the strong and numerically weak sense, but at a rate that is slower then any power law.
We give a concise overview of the theory of regularity structures as first exposed in [Hai14]. In order to allow to focus on the conceptual aspects of the theory, many proofs are omitted and statements are simplified. In order to provide both motivation and focus, we concentrate on the study of solutions to the stochastic quantisation equations for the Euclidean $Φ^4_3$ quantum field theory which can be obtained with the help of this theory. In particular, we sketch the proofs of how one can show that this model arises quite naturally as an idealised limiting object for several classes of smooth models.
We consider stochastic differential equations of the form $dY_t=V(Y_t)\,dX_t+V_0(Y_t)\,dt$ driven by a multi-dimensional Gaussian process. Under the assumption that the vector fields $V_0$ and $V=(V_1,\ldots,V_d)$ satisfy Hörmander's bracket condition, we demonstrate that $Y_t$ admits a smooth density for any $t\in(0,T]$, provided the driving noise satisfies certain nondegeneracy assumptions. Our analysis relies on relies on an interplay of rough path theory, Malliavin calculus and the theory of Gaussian processes. Our result applies to a broad range of examples including fractional Brownian motion with Hurst parameter $H>1/4$, the Ornstein-Uhlenbeck process and the Brownian bridge returning after time $T$.
We introduce a new notion of "regularity structure" that provides an algebraic framework allowing to describe functions and / or distributions via a kind of "jet" or local Taylor expansion around each point. The main novel idea is to replace the classical polynomial model which is suitable for describing smooth functions by arbitrary models that are purpose-built for the problem at hand. In particular, this allows to describe the local behaviour not only of functions but also of large classes of distributions. We then build a calculus allowing to perform the various operations (multiplication, composition with smooth functions, integration against singular kernels) necessary to formulate fixed point equations for a very large class of semilinear PDEs driven by some very singular (typically random) input. This allows, for the first time, to give a mathematically rigorous meaning to many interesting stochastic PDEs arising in physics. The theory comes with convergence results that allow to interpret the solutions obtained in this way as limits of classical solutions to regularised problems, possibly modified by the addition of diverging counterterms. These counterterms arise naturally through the action of a "renormalisation group" which is defined canonically in terms of the regularity structure associated to the given class of PDEs. As an example of a novel application, we solve the long-standing problem of building a natural Markov process that is symmetric with respect to the (finite volume) measure describing the $Φ^4_3$ Euclidean quantum field theory. It is natural to conjecture that the Markov process built in this way describes the Glauber dynamic of 3-dimensional ferromagnets near their critical temperature.
We consider a two-parameter averaging-homogenization type elliptic problem together with the stochastic representation of the solution. A limit theorem is derived for the corresponding diffusion process and a precise description of the two-parameter limit behavior for the solution of the PDE is obtained.
In this article we consider rough differential equations (RDEs) driven by non-geometric rough paths, using the concept of branched rough paths introduced in Gubinelli (2004). We first show that branched rough paths can equivalently be defined as $γ$-Hölder continuous paths in some Lie group, akin to geometric rough paths. We then show that every branched rough path can be encoded in a geometric rough path. More precisely, for every branched rough path $\mathbf{X}$ lying above a path $X$, there exists a geometric rough path $\bar{\mathbf{X}}$ lying above an extended path $\bar X$, such that $\bar{\mathbf{X}}$ contains all the information of $\mathbf{X}$. As a corollary of this result, we show that every RDE driven by a non-geometric rough path $\mathbf{X}$ can be rewritten as an extended RDE driven by a geometric rough path $\bar{\mathbf{X}}$. One could think of this as a generalisation of the Itô-Stratonovich correction formula.
We propose a modification, based on the RESTART (repetitive simulation trials after reaching thresholds) and DPR (dynamics probability redistribution) rare event simulation algorithms, of the standard diffusion Monte Carlo (DMC) algorithm. The new algorithm has a lower variance per workload, regardless of the regime considered. In particular, it makes it feasible to use DMC in situations where the "naïve" generalisation of the standard algorithm would be impractical, due to an exponential explosion of its variance. We numerically demonstrate the effectiveness of the new algorithm on a standard rare event simulation problem (probability of an unlikely transition in a Lennard-Jones cluster), as well as a high-frequency data assimilation problem.
These are short notes from a series of lectures given at the University of Rennes in June 2013, at the University of Bonn in July 2013, at the XVIIth Brazilian School of Probability in Mambucaba in August 2013, and at ETH Zurich in September 2013. We give a concise overview of the theory of regularity structures as exposed in Hairer (2014). In order to allow to focus on the conceptual aspects of the theory, many proofs are omitted and statements are simplified. We focus on applying the theory to the problem of giving a solution theory to the stochastic quantisation equations for the Euclidean $Φ^4_3$ quantum field theory.
We study the stochastic Allen-Cahn equation driven by a noise term with intensity $\sqrt{\varepsilon}$ and correlation length $δ$ in two and three spatial dimensions. We study diagonal limits $δ, \varepsilon \to 0$ and describe fully the large deviation behaviour depending on the relationship between $δ$ and $\varepsilon$. The recently developed theory of regularity structures allows to fully analyse the behaviour of solutions for vanishing correlation length $δ$ and fixed noise intensity $\varepsilon$. One key fact is that in order to get non-trivial limits as $δ\to 0$, it is necessary to introduce diverging counterterms. The theory of regularity structures allows to rigorously analyse this renormalisation procedure for a number of interesting equations. Our main result is a large deviation principle for these renormalised solutions. One interesting feature of this result is that the diverging renormalisation constants disappear at the level of the large deviations rate function. We apply this result to derive a sharp condition on $δ, \varepsilon$ that guarantees a large deviation principle for diagonal schemes $\varepsilon, δ\to 0$ for the equation without renormalisation.
We present a series of recent results on the well-posedness of very singular parabolic stochastic partial differential equations. These equations are such that the question of what it even means to be a solution is highly non-trivial. This problem can be addressed within the framework of the recently developed theory of "regularity structures", which allows to describe candidate solutions locally by a "jet", but where the usual Taylor polynomials are replaced by a sequence of custom-built objects. In order to illustrate the theory, we focus on the particular example of the Kardar-Parisi-Zhang equation, a popular model for interface propagation.
We study the problem of sampling high and infinite dimensional target measures arising in applications such as conditioned diffusions and inverse problems. We focus on those that arise from approximating measures on Hilbert spaces defined via a density with respect to a Gaussian reference measure. We consider the Metropolis-Hastings algorithm that adds an accept-reject mechanism to a Markov chain proposal in order to make the chain reversible with respect to the target measure. We focus on cases where the proposal is either a Gaussian random walk (RWM) with covariance equal to that of the reference measure or an Ornstein-Uhlenbeck proposal (pCN) for which the reference measure is invariant. Previous results in terms of scaling and diffusion limits suggested that the pCN has a convergence rate that is independent of the dimension while the RWM method has undesirable dimension-dependent behaviour. We confirm this claim by exhibiting a dimension-independent Wasserstein spectral gap for pCN algorithm for a large class of target measures. In our setting this Wasserstein spectral gap implies an $L^2$-spectral gap. We use both spectral gaps to show that the ergodic average satisfies a strong law of large numbers, the central limit theorem and nonasymptotic bounds on the mean square error, all dimension independent. In contrast we show that the spectral gap of the RWM algorithm applied to the reference measures degenerates as the dimension tends to infinity.
We provide a mathematical study of the modified Diffusion Monte Carlo (DMC) algorithm introduced in the companion article \cite{DMC}. DMC is a simulation technique that uses branching particle systems to represent expectations associated with Feynman-Kac formulae. We provide a detailed heuristic explanation of why, in cases in which a stochastic integral appears in the Feynman-Kac formula (e.g. in rare event simulation, continuous time filtering, and other settings), the new algorithm is expected to converge in a suitable sense to a limiting process as the time interval between branching steps goes to 0. The situation studied here stands in stark contrast to the "naïve" generalisation of the DMC algorithm which would lead to an exponential explosion of the number of particles, thus precluding the existence of any finite limiting object. Convergence is shown rigorously in the simplest possible situation of a random walk, biased by a linear potential. The resulting limiting object, which we call the "Brownian fan", is a very natural new mathematical object of independent interest.
We introduce the dynamical sine-Gordon equation in two space dimensions with parameter $β$, which is the natural dynamic associated to the usual quantum sine-Gordon model. It is shown that when $β^2 \in (0,\frac{16π}{3})$ the Wick renormalised equation is well-posed. In the regime $β^2 \in (0,4π)$, the Da Prato-Debussche method applies, while for $β^2 \in [4π,\frac{16π}{3})$, the solution theory is provided via the theory of regularity structures (Hairer 2013). We also show that this model arises naturally from a class of $2+1$-dimensional equilibrium interface fluctuation models with periodic nonlinearities. The main mathematical difficulty arises in the construction of the model for the associated regularity structure where the role of the noise is played by a non-Gaussian random distribution similar to the complex multiplicative Gaussian chaos recently analysed by Lacoin, Rhodes and Vargas (2013).
We study approximations to a class of vector-valued equations of Burgers type driven by a multiplicative space-time white noise. A solution theory for this class of equations has been developed recently in [Hairer, Weber, Probab. Theory Related Fields, to appear]. The key idea was to use the theory of controlled rough paths to give definitions of weak / mild solutions and to set up a Picard iteration argument. In this article the limiting behaviour of a rather large class of (spatial) approximations to these equations is studied. These approximations are shown to converge and convergence rates are given, but the limit may depend on the particular choice of approximation. This effect is a spatial analogue to the Itô-Stratonovich correction in the theory of stochastic ordinary differential equations, where it is well known that different approximation schemes may converge to different solutions.
In this article, we consider the problem of homogenising the linear heat equation perturbed by a rapidly oscillating random potential. We consider the situation where the space-time scaling of the potential's oscillations is \textit{not} given by the diffusion scaling that leaves the heat equation invariant. Instead, we treat the case where spatial oscillations are much faster than temporal oscillations. Under suitable scaling of the amplitude of the potential, we prove convergence to a deterministic heat equation with constant potential, thus completing the results previously obtained in \cite{MR2962093}.
We consider differential equations driven by rough paths and study the regularity of the laws and their long time behavior. In particular, we focus on the case when the driving noise is a rough path valued fractional Brownian motion with Hurst parameter $H\in(\frac{1}{3},\frac{1}{2}]$. Our contribution in this work is twofold. First, when the driving vector fields satisfy Hörmander's celebrated "Lie bracket condition," we derive explicit quantitative bounds on the inverse of the Malliavin matrix. En route to this, we provide a novel "deterministic" version of Norris's lemma for differential equations driven by rough paths. This result, with the added assumption that the linearized equation has moments, will then yield that the transition laws have a smooth density with respect to Lebesgue measure. Our second main result states that under Hörmander's condition, the solutions to rough differential equations driven by fractional Brownian motion with $H\in(\frac{1}{3},\frac{1}{2}]$ enjoy a suitable version of the strong Feller property. Under a standard controllability condition, this implies that they admit a unique stationary solution that is physical in the sense that it does not "look into the future."
We construct solutions to Burgers type equations perturbed by a multiplicative space-time white noise in one space dimension. Due to the roughness of the driving noise, solutions are not regular enough to be amenable to classical methods. We use the theory of controlled rough paths to give a meaning to the spatial integrals involved in the definition of a weak solution. Subject to the choice of the correct reference rough path, we prove unique solvability for the equation and we show that our solutions are stable under smooth approximations of the driving noise.
We consider a class of stochastic PDEs of Burgers type in spatial dimension 1, driven by space-time white noise. Even though it is well known that these equations are well posed, it turns out that if one performs a spatial discretization of the nonlinearity in the "wrong" way, then the sequence of approximate equations does converge to a limit, but this limit exhibits an additional correction term. This correction term is proportional to the local quadratic cross-variation (in space) of the gradient of the conserved quantity with the solution itself. This can be understood as a consequence of the fact that for any fixed time, the law of the solution is locally equivalent to Wiener measure, where space plays the role of time. In this sense, the correction term is similar to the usual Itô-Stratonovich correction term that arises when one considers different temporal discretizations of stochastic ODEs.
We introduce a new concept of solution to the KPZ equation which is shown to extend the classical Cole-Hopf solution. This notion provides a factorisation of the Cole-Hopf solution map into a "universal" measurable map from the probability space into an explicitly described auxiliary metric space, composed with a new solution map that has very good continuity properties. The advantage of such a formulation is that it essentially provides a pathwise notion of a solution, together with a very detailed approximation theory. In particular, our construction completely bypasses the Cole-Hopf transform, thus laying the groundwork for proving that the KPZ equation describes the fluctuations of systems in the KPZ universality class. As a corollary of our construction, we obtain very detailed new regularity results about the solution, as well as its derivative with respect to the initial condition. Other byproducts of the proof include an explicit approximation to the stationary solution of the KPZ equation, a well-posedness result for the Fokker-Planck equation associated to a particle diffusing in a rough space-time dependent potential, and a new periodic homogenisation result for the heat equation with a space-time periodic potential. One ingredient in our construction is an example of a non-Gaussian rough path such that the area process of its natural approximations needs to be renormalised by a diverging term for the approximations to converge.
We study the spatial homogenisation of parabolic linear stochastic PDEs exhibiting a two-scale structure both at the level of the linear operator and at the level of the Gaussian driving noise. We show that in some cases, in particular when the forcing is given by space-time white noise, it may happen that the homogenised SPDE is not what one would expect from existing results for PDEs with more regular forcing terms.
We consider the stochastic Allen-Cahn equation driven by mollified space-time white noise. We show that, as the mollifier is removed, the solutions converge weakly to 0, independently of the initial condition. If the intensity of the noise simultaneously converges to 0 at a sufficiently fast rate, then the solutions converge to those of the deterministic equation. At the critical rate, the limiting solution is still deterministic, but it exhibits an additional damping term.
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.
The aim of this note is to provide a short and self-contained proof of Hörmander's theorem about the smoothness of transition probabilities for a diffusion under Hörmander's "brackets condition". While both the result and the technique of proof are well-known, the exposition given here is novel in two aspects. First, we introduce Malliavin calculus in an "intuitive" way, without using Wiener's chaos decomposition. While this may make it difficult to prove some of the standard results in Malliavin calculus (boundedness of the derivative operator in $L^p$ spaces for example), we are able to bypass these and to replace them by weaker results that are still sufficient for our purpose. Second, we introduce a notion of "almost implication" and "almost truth" (somewhat similar to what is done in fuzzy logic) which allows, once the foundations of Malliavin calculus are laid out, to give a very short and streamlined proof of Hörmader's theorem that focuses on the main ideas without clouding it by technical details.
We consider a diffusion process with coefficients that are periodic outside of an "interface region" of finite thickness. The question investigated in this article is the limiting long time/large scale behavior of such a process under diffusive rescaling. It is clear that outside of the interface, the limiting process must behave like Brownian motion, with diffusion matrices given by the standard theory of homogenization. The interesting behavior therefore occurs on the interface. Our main result is that the limiting process is a semimartingale whose bounded variation part is proportional to the local time spent on the interface. The proportionality vector can have nonzero components parallel to the interface, so that the limiting diffusion is not necessarily reversible. We also exhibit an explicit way of identifying its parameters in terms of the coefficients of the original diffusion. Similarly to the one-dimensional case, our method of proof relies on the framework provided by Freidlin and Wentzell [Ann. Probab. 21 (1993) 2215--2245] for diffusion processes on a graph in order to identify the generator of the limiting process.
A series of recent articles introduced a method to construct stochastic partial differential equations (SPDEs) which are invariant with respect to the distribution of a given conditioned diffusion. These works are restricted to the case of elliptic diffusions where the drift has a gradient structure and the resulting SPDE is of second-order parabolic type. The present article extends this methodology to allow the construction of SPDEs which are invariant with respect to the distribution of a class of hypoelliptic diffusion processes, subject to a bridge conditioning, leading to SPDEs which are of fourth-order parabolic type. This allows the treatment of more realistic physical models, for example, one can use the resulting SPDE to study transitions between meta-stable states in mechanical systems with friction and noise. In this situation the restriction of the drift being a gradient can also be lifted.
We investigate the behaviour of a chain of interacting Brownian particles with one end fixed and the other moving away at slow speed, in the limit of small noise. The interaction between particles is through a pairwise potential with finite range. We consider both overdamped and underdamped dynamics.
This article is devoted to the numerical study of various finite difference approximations to the stochastic Burgers equation. Of particular interest in the one-dimensional case is the situation where the driving noise is white both in space and in time. We demonstrate that in this case, different finite difference schemes converge to different limiting processes as the mesh size tends to zero. A theoretical explanation of this phenomenon is given and we formulate a number of conjectures for more general classes of equations, supported by numerical evidence.
We demonstrate that stochastic differential equations (SDEs) driven by fractional Brownian motion with Hurst parameter H > 1/2 have similar ergodic properties as SDEs driven by standard Brownian motion. The focus in this article is on hypoelliptic systems satisfying Hörmander's condition. We show that such systems satisfy a suitable version of the strong Feller property and we conclude that they admit a unique stationary solution that is physical in the sense that it does not "look into the future". The main technical result required for the analysis is a bound on the moments of the inverse of the Malliavin covariance matrix, conditional on the past of the driving noise.
The Metropolis-Adjusted Langevin Algorithm (MALA), originally introduced to sample exactly the invariant measure of certain stochastic differential equations (SDE) on infinitely long time intervals, can also be used to approximate pathwise the solution of these SDEs on finite time intervals. However, when applied to an SDE with a nonglobally Lipschitz drift coefficient, the algorithm may not have a spectral gap even when the SDE does. This paper reconciles MALA's lack of a spectral gap with its ergodicity to the invariant measure of the SDE and finite time accuracy. In particular, the paper shows that its convergence to equilibrium happens at exponential rate up to terms exponentially small in time-stepsize. This quantification relies on MALA's ability to exactly preserve the SDE's invariant measure and accurately represent the SDE's transition probability on finite time intervals.
We consider a one-dimensional diffusion process with coefficients that are periodic outside of a finite 'interface region'. The question investigated in this article is the limiting long time / large scale behaviour of such a process under diffusive rescaling. Our main result is that it converges weakly to a rescaled version of skew Brownian motion, with parameters that can be given explicitly in terms of the coefficients of the original diffusion. Our method of proof relies on the framework provided by Freidlin and Wentzell for diffusion processes on a graph in order to identify the generator of the limiting process. The graph in question consists of one vertex representing the interface region and two infinite segments corresponding to the regions on either side.
In this article, we show how the theory of rough paths can be used to provide a notion of solution to a class of nonlinear stochastic PDEs of Burgers type that exhibit too high spatial roughness for classical analytical methods to apply. In fact, the class of SPDEs that we consider is genuinely ill-posed in the sense that different approximations to the nonlinearity may converge to different limits. Using rough paths theory, a pathwise notion of solution to these SPDEs is formulated, and we show that this yields a well-posed problem, which is stable under a large class of perturbations, including the approximation of the rough driving noise by a mollified version and the addition of hyperviscosity. We also show that under certain structural assumptions on the coefficients, the SPDEs under consideration generate a reversible Markov semigroup with respect to a diffusion measure that can be given explicitly.
We consider a class of singular perturbations to the stochastic heat equation or semilinear variations thereof. The interesting feature of these perturbations is that, as the small parameter epsilon tends to zero, their solutions converge to the 'wrong' limit, i.e. they do not converge to the solution obtained by simply setting epsilon = 0. A similar effect is also observed for some (formally) small stochastic perturbations of a deterministic semilinear parabolic PDE. Our proofs are based on a detailed analysis of the spatially rough component of the equations, combined with a judicious use of Gaussian concentration inequalities.
We study a model of two interacting Hamiltonian particles subject to a common potential in contact with two Langevin heat reservoirs: one at finite and one at infinite temperature. This is a toy model for 'extreme' non-equilibrium statistical mechanics. We provide a full picture of the long-time behaviour of such a system, including the existence / non-existence of a non-equilibrium steady state, the precise tail behaviour of the energy in such a state, as well as the speed of convergence toward the steady state. Despite its apparent simplicity, this model exhibits a surprisingly rich variety of long time behaviours, depending on the parameter regime: if the surrounding potential is 'too stiff', then no stationary state can exist. In the softer regimes, the tails of the energy in the stationary state can be either algebraic, fractional exponential, or exponential. Correspondingly, the speed of convergence to the stationary state can be either algebraic, stretched exponential, or exponential. Regarding both types of claims, we obtain matching upper and lower bounds.
It is well-known under the name of `periodic homogenization' that, under a centering condition of the drift, a periodic diffusion process on R^d converges, under diffusive rescaling, to a d-dimensional Brownian motion. Existing proofs of this result all rely on uniform ellipticity or hypoellipticity assumptions on the diffusion. In this paper, we considerably weaken these assumptions in order to allow for the diffusion coefficient to even vanish on an open set. As a consequence, it is no longer the case that the effective diffusivity matrix is necessarily non-degenerate. It turns out that, provided that some very weak regularity conditions are met, the range of the effective diffusivity matrix can be read off the shape of the support of the invariant measure for the periodic diffusion. In particular, this gives some easily verifiable conditions for the effective diffusivity matrix to be of full rank. We also discuss the application of our results to the homogenization of a class of elliptic and parabolic PDEs.