Notes on Hydrodynamic Limits and Related Topics
Sunder Sethuraman
University of Arizona
In these lecture notes, we discuss various ‘hydrodynamic LLN’ and ‘CLT’ scaling limits, among others, in types of stochastic interacting particle systems, connecting ‘microscopic’ behaviors to continuum laws. Via ‘short stories’, the aim is to present some of the ‘basics’ for students and those entering the field, as a complement to books such as [130], [133], [147], [148]. To be concrete, attention is restricted to a few ‘mass conservative’ systems on discrete spaces, namely exclusion and zero-range processes, that have proved robust in the study of different phenomena.
After preliminaries, we discuss the ‘entropy’ and ‘relative entropy’ methods to prove hydrodynamic limits of the bulk mass in finite volume, as well as other items such as construction of systems in infinite volume and the structure of their invariant measures, and scaling limits of local functionals, such as occupation times of sites and the motion of a tagged particle. In the last part, we also discuss equilibrium fluctuations of the bulk mass when the process starts from an invariant measure.
While such a list forms a brief introduction, there are of course several results not covered. For instance, those with respect to large deviations, mixing times, metastability and condensation phenomena, non-equilibrium fluctuations, the ‘KPZ’ equation and ‘KPZ’ fixed point, and integrable probability, though referenced, are not much discussed.
Acknowledgements. These notes were begun for a topics course given in Spring 2012 at the University of Arizona, with input from research papers and texts, and revised during a 2026 sabbatical, supported in part by the Simons Foundation. They have also served as background for recent minicourses in Crete 2023 (Stochastic Methods in Finance and Physics IV), Rio 2024 (IMPA), and Hermosillo 2025 (XV Symposium on Probability and Stochastic Processes). Many thanks to the participants in these venues for their comments and interest.
Contents
- 1 Section
- 2 Section
- 3 Section
- 4 Section
- 5 Section
- 6 Section
- 7 Section
- 8 Section
- 9 Section
- 10 Section
- 11 Section
- 12 Section
- Bibliography
Section 1 Preliminaries and hydrodynamics of independent random walks
Before discussing a basic example, illustrating possibilities in the study of ‘hydrodynamics of stochastic particle systems’, we recall some basic notions in Markov chains to set the stage.
1.1. Markov chains
1.1.1. Construction
A family of random variables taking values on a countable state space is called a ‘discrete time Markov chain’ if the ‘stationary Markov property’ is satisfied:
for all and . When and , the last quantity on the right-side represents a ‘transition probability’ of the Markov chain, denoted . The -step probability satisfies a recurrence, .
We now construct a ‘continuous time Markov chain’ on with ‘skeleton’ and a transition probability vanishing on the diagonal, that is for all . Let be a collection of positive numbers, and let be a collection of independent identically distributed exponential random variables with rate , independent of the skeleton discrete time chain.
Define the process as follows: Initially, . After time , the process jumps to value , and after a subsequent time , the process jumps to value , and so on. Let for . Then,
where . Sufficient conditions for include the cases that the state space is finite, or that . We will assume from now on that the process is ‘regular’, that is , so that the process is defined for all time .
One can show that the ‘continuous time’ chain satisfies the stationary Markov property, which in this context is equivalent to
| (1.1.1) | |||||
for all , and .
Conversely, given a process satisfying (1.1.1) and also the ‘jump property’ that there exists a sequence of strictly increasing stopping times such that and is constant on intervals and for , one can determine a unique skeleton discrete time Markov chain and positive jump parameters such that , , and are independent exponentials with rates for . Chains with the same skeleton and jump parameters have the same joint distributions. The condition ensures vanishes on the diagonal.
1.1.2. Generators and Kolmogorov equations
In the discrete time setting, we can define the transition operator , which is a matrix when is finite. Then, the th powers give the th step probabilities, .
Computing the -time probabilities for the continuous time Markov chain is more complicated. Define the transition probability
Then, by the Markov property we have
or in terms of operators , the semigroup property.
Given regularity of the process, the transition functions are differentiable in time, and satisfy the ‘first jump’ relation
from which the backward equation is seen:
Similarly, from a decomposition of the jump time just before time , one has the forward equation
Exercise 1.1.1.
Review and perform these derivations.
Define the operator
Then, neatly expressed, the backward and forward equations become
Also, and as .
When the space is finite, is a ‘generator’ matrix, that is for and , and by solving the ODE’s, one obtains which can be computed in some cases. More generally, the semigroup can be understood in terms of the Hille-Yosida formulation.
Let be a bounded function on the state space. In the discrete time case, define , which is the conditional expectation of given . Then, is the conditional expectation of given , where is the -fold convolution, or th step transition probability.
In the continuous case, define which is the conditional expectation of given . In this framework, the generator is often expressed in terms of its action on compactly supported :
1.1.3. Invariant measures
Let be a probability measure on . In the discrete time situation, define . Hence, we see that is the distribution at time when the initial state is distributed according to .
In the continuous-time model, define
We say that is an ‘invariant measure’ for discrete time chains if , and for continuous time chains if for all . We also say that is a ‘reversible’ invariant measure in discrete time chains if for all . In continuous time chains, is ‘reversible’ when is self-adjoint, that is , or in terms of the inner product, for all compactly supported .
One can verify that in the continuous time setting that being invariant is equivalent to or in terms of expectations for all compactly supported . Also, being reversible is equivalent to for all , or in terms of inner products for all compactly supported .
In the finite state space case, invariant measures always exist, and if the skeleton chain can reach every state from any state in finite time, that is is irreducible, the invariant measure is unique. However, in the countable state case there may be no invariant measures.
Reversibility has the following interesting implication. Fix a time , and consider for . Suppose that initially is distributed according to an invariant measure . Then, it can be seen that is a continuous time Markov chain with transition probability . In particular, when is reversible, and in this case the ‘forward in time’ and ‘backward in time’ chains have the same distribution. We remark that it may be easier to find directly a reversible measure and therefore an invariant measure by checking the reversibility conditions.
1.1.4. Examples
We will content ourselves for the moment with two basic continuous time examples, the two-state Markov chain, and random walk.
Example 1.1.2.
Let correspond to states ‘on’ and ‘off’, or sometimes ‘empty’ and ‘occupied’. Here, state can transition to state and vice versa. Let and be the corresponding jump rates. The skeleton chain probabilities are . The generator matrix is
and correspondingly, it is an exercise in diagonalization to compute and to find the unique invariant measure.
Exercise 1.1.3.
Compute the invariant measure and . Is it reversible?
Example 1.1.4.
We define the Poisson process with parameter . Let , the whole numbers. Let . The skeleton chain is deterministic where transitions are to the nearest right site: for . Then, with , the continuous time chain counts the number steps made up to time . The trajectories are step functions, which are right-continuous, with left limits.
Exercise 1.1.5.
Compute the backward/forward equation, and show that the transition probability is in ‘Poisson’ form .
Example 1.1.6.
Let the -dimensional torus where , or integers modulo . Let also , and be a finite-range, translation-invariant transition probability: For all , if some , and . We will also assume that is irreducible. For instance, the nearest-neighbor, symmetric case is one possibility.
Then, for and . In particular, the uniform distribution is the unique invariant measure. Moreover, is reversible exactly when for all .
Now let be the number of jumps before time . Since the jump rates are all , we see that is a Poisson process with rate . In particular, the transition probability satisfies
We now consider a sequence of chains on a sequence of torii as increases. Define to be the mean displacement of the position. Then, it is not difficult to extablish the following law of large numbers,
Here, where is the unit circle.
Exercise 1.1.7.
Show this LLN (law of large numbes) by say variance computations. Noting is Poisson with variance is useful. One can do it also by regeneration, and other methods.
Also, if , let be the matrix of covariances, for . We have the central limit theorem,
Exercise 1.1.8.
There are a few ways to show this CLT (central limit theorem). One way is to compute the moment generating function or characteristic function, noting independence of and the displacements.
Finally, when , we say the walk is asymmetric, and when and is not symmetric, we say the walk is mean-zero asymmetric, and when is symmetric, we say the walk is symmetric.
1.2. Evolution of the mass of independent random walks
We would like to understand the space-time evolution of the mass in a system of particles with a conservation law. Perhaps the simplest model is that of unlabeled non-interacting random walks on a -dimensional torus with locations. When is large, and one looks at the system from afar, after long times, one can discern better the motion of the bulk of the mass rather than individual components. The goal is to make precise the evolution of the mass in this scale in terms of a continuum equation.
We will be working with a Markov chain on , where we recall , which governs the motion of independent random walks on specified in Example 1.1.6. Since we are interested in the ‘mass’ of particles, we will consider the occupation numbers at each location on the lattice . That is, let be the position of the th particle at time . Define
We now observe that the process is a Markov chain. Indeed, given independence and the Markov property of the individual particle movement, by splitting over all possibilities, the Markov property of can be deduced.
Exercise 1.2.1.
Show this Markov property.
1.2.1. Invariant measures and distribution at time
What are the invariant measures for the process? Since the process , corresponding to particles, is irreducible, there is a unique invariant measure. It is not so easy to characterize it immediately.
Indeed, let us relax the assumption that there are particles in the system. If we do not specify the initial number of particles, then is no longer irreducible, since there is no birth or death possible: For instance, a system with particles cannot evolve to one with random walks. Nevertheless, we may specify in good form several invariant measures for this ‘relaxed’ reducible system.
Recall the Poisson distribution with parameter : for . Its moment generating function is given by
For a nonnegative function , define the product measure on so that the means
When is constant, we denote .
The process belongs to the space of right-continuous paths with left limits in . We will denote by and the probability measure and expectation with respect to the evolution of the process when initially is distributed according to . On the other hand, refers to the expectation with respect to on .
Now, starting from , we observe that we may calculate the distribution at later times .
Proposition 1.2.2.
Under , the distribution of is the inhomogeneous product of Poisson measures where is the expectation with respect to the position of a random walk with rates starting from the origin at time .
That the later time -distribution would still be a product measure is a feature of the independence of the particles–not a generic feature in more general interacting particle systems!
When is constant, we immediately arrive at the following characterization.
Corollary 1.2.3.
The measures are invariant for the Markov chain .
Proof of Proposition 1.2.2. We need only compute the moment generating function of . Write
and, for ,
where denotes the position at time of the th particle initially at location .
Since particles move independently, and initially there are a Poisson number of particles on each site of the lattice,
where is the position of a random walk on starting at the origin.
Now,
and
Hence, after a calculation,
Since
the proof concludes. ∎
To come back to our initial question, we remark that the invariant measure can be decomposed in terms of its restrictions to the sets for . These are closed and irreducible for the motion. Then, each is the unique invariant measure on the restriction, and does not depend on . In physics terminology, is the ‘grand canonical’ measure and is the ‘canonical’ one.
Since the mean of under equals , it makes sense to call the ‘mass density’ of the process. In this way, is a family of invariant measures indexed by density .
When is not constant, the measures are examples of ‘local equilibrium’ measures, as near a continuity of point of , they distribute particles near near that of the invariant measure . One may consider other non-product ‘local equilibrium’ and ‘nonequilibrium’ measures, but, as we have seen, these local equilibrium measures allow for some computations.
1.2.2. Hydrodynamics in mean-value
If we start with a local equilibrium measure , then initially, the means are a discretization of the density ‘profile’ defined on the continuous space . At a later time , the means have evolved to . But, how to understand a macroscopic picture?
The answers depend on the particular time and space scales chosen in the problem. We will think of as embedded in where grid points on are separated by distance . In this way, a ‘macroscopic’ point on corresponds to the ‘microscopic’ point . As we will see, time should now be appropriately speeded up to see movement of the system. How fast this speed up should be will depend on the structure of the underlying jump probability .
When is fixed, one can see that converges weakly to . Hence, at a continuity point for , we have . So, if we do not speed up time at all, the system does not move.
Asymmetric motions
In the asymmetric setting, let . Since in probability, we have, for ,
In this case, when the initial profile is continuous,
Therefore, if we speed up time by a factor of , we see that the density profile has translated by . In this sense is referred to as the ‘microscopic’ time, and as the ‘macroscopic’ time. Moreover, the density satisfies
| (1.2.1) |
This makes sense as individual particles displace an order microscopic locations at microscopic time .
Mean-zero motions
However, when , particles do not displace as much, but follow the ‘square root’ law, when say the underlying jump rates are finite-range. In other words, displacements are of order at time , or alternatively of order at times . The latter version fits in nicely with our space scaling of .
By the central limit theorem for random walks in this case, and, when again is continuous, we have
Here, is the periodic extension of with period , and is the Gaussian density with covariance . It follows that , being a convolution with the , satisfies the heat equation on :
| (1.2.2) |
1.2.3. Conclusion
What we have done so far is to derive a macroscopic ‘hydrodynamic limit’ for the mean values over . In the next Section 2, we will view the hydrodynamic limit as a full fledged law of large numbers of the empirical measure of particles. We also give some physical motivation for the name ‘hydrodynamics’.
1.3. Notes
The material on Markov chains is a standard treatment based on the development in [127], [178] and [130, Appendix 1]. The discussion on hydrodynamics of independent walks follows that in [130, Chapter 1] and [135]. See also [157] for a more general treatment. An early rigorous work on hydrodynamics of independent particle systems is [62].
Section 2 Empirical measures, and a first look at hydrodynamics of exclusion processes
We introduce a notion of ‘hydrodynamics’ via empirical measures, which will be used in the analysis of more general interacting particle systems. Then, after a preliminary discussion of certain martingales in the Markov chain context, we discuss the hydrodynamics of exclusion processes. At the end, we comment on the physical motivation behind the name ‘hydrodynamics’.
2.1. Empirical measures and view of ‘hydrodynamics’ as a LLN
We may recast the derivation of ‘hydrodynamics’ in systems of independent random walks in terms of the random space-time scaled empirical distribution
Here, is a member of , the nonnegative measures on . Recall our definition of the independent particle process in Section 1: We derived the distribution of the particle numbers at later times, starting from a local equilibrium measure . However, it will be easier in what follows to consider the weaker notion of the asymptotic behavior of the empirical distribution. In this sense, what is meant by ‘hydrodynamics’ is a law of large numbers, characterizing ‘first-order’ behavior.
Let be a smooth, bounded function on , and write
| (2.1.1) |
Here, the pairing is another way to write the expectation . Recall, from Section 1, the distribution of starting from is a product of Poisson measures with intensity , where the expectation is with respect to the random walk . Then, in mean-value, starting in , we have
Depending on whether or , in which case or , the last quantity as before tends respectively to
| (2.1.2) |
However, the variance of (2.1.1), when starting in , since the occupation numbers are independent Poisson variables with intensities , vanishes:
as and are bounded. Hence, converges variously to the expressions in (2.1.2) in probability, depending on the drift .
In particular, we have shown, with respect to the initial distribution for the process , the empirical measure , a random element of , converges in probability to the deterministic measure corresponding to the macroscopic space-time mass evolution.
Here, the topology on used is as follows: Consider , the space of real continuous functions on endowed with the sup-metric. Let be a dense, countable family of continuous functions in . Then, define the distance on by
Therefore, capturing the limit behavior of for each continuous function is enough to compute the limits of .
2.2. Exclusion processes
The exclusion process is one of the ‘canonical’ interacting particle systems, introduced in [205]. Consider particles on with the minimal interaction that no particle can jump onto another. Accordingly, a configuration of occupation numbers belongs to the finite state space where or depending on whether is empty or occupied at time . Informally, updates in that each particle is a continuous time random walk carrying an exponential clock. When a clock rings, the particle tries to move with skeleton jump probability . However, if the site chosen is already occupied, the jump is suppressed, and all clocks reset. We will restrict attention to jump probabilities which are tranlsation-invariant: .
More formally, infinitesimally can change to when a particle at displaces by increment where
with rate . The ‘exclusion’ factor ‘’ is exactly when is occupied and the destination is unoccupied; otherwise, it is . The generator of the process is given by
| (2.2.1) |
for functions .
In the following, we will also assume that is finite-range, that is for for some . Of course, a case is when is nearest-neighbor, that is when the range .
When is symmetric, the process is called the ‘symmetric simple exclusion process’. When is asymmetric, is termed the ‘asymmetric exclusion process’. When the range , the label ‘simple’ is sometimes added to the names.
There is a simplification of the form of when is symmetric. Namely, since and , we have
| (2.2.2) | |||||
The last line follows as the term in curly braces equals exactly when the difference in square brackets vanishes.
Let be the product measure on with Bernoulli marginals with success probability . Recall from Subsection 1.1.3 that stands for the expectation under , and and for the process measure and expectation when starting in .
Proposition 2.2.1.
The measures are invariant measures for , for both symmetric and asymmetric . In the symmetric case, is also reversible.
Proof.
First note that vanishes if . Hence, we may write
where we dropped the factor ‘’.
Note also, under the change of measure , which exchanges values and , the measure remains the same. Hence, for functions , the term
Then, by collecting terms, with simple manipulation,
where the -adjoint is seen as the exclusion generator with single particle jump probability .
Now, since , by inspection, we have that for all . This shows is invariant.
Moreover, as when is symmetric, is reversible. ∎
From this result, we can deduce when there are exactly particles in the system that is the unique ‘canonical’ invariant measure for the motion on .
Exercise 2.2.2 (Duality).
Suppose is symmetric. Show that the space of linear combinations of occupation variables for remains invariant under action by generator . In fact, the space of linear combinations of for is closed. We comment that this property is sometimes referred to as ‘duality’. Hint: A case of the property is straightforwardly seen by computing the action on the variable .
2.3. Martingales and Markov chains
Recall that a martingale corresponding to sigma-fields is a random process, adapted to the filtration , which satisfies
for all .
Exercise 2.3.1.
Let be a Poisson process with rate . Then, is a martingale with respect to the ‘natural’ sigma-fields . Also, is also a martingale with respect to .
Let be a Markov process on a countable state space . Let be a twice continuously differentiable function whose first and second partial time derivatives are uniformly bounded. Define
These two processes, which we show below are martingales, will be very useful in our stochastic analysis of Markov systems. The term,
is the ‘(predictable) quadratic variation’ of the martingale .
Proposition 2.3.2.
With respect to natural sigma-fields , both and are martingales.
Proof.
We will show that is a martingale when does not depend on time. Generalizations and verification of as a martingale are left to the reader. We need only show that
Now, and from the Markov property. From the forward equation, the derivative of the left-side of the above display in equals
At time , the left-side also vanishes. This concludes the proof. ∎
Exercise 2.3.3.
Complete the proof of Proposition 2.3.2. Hint: With respect to , we need to show
When , the relation holds. Hence, if one shows the derivatives with respect to match, that is
2.4. Sketch of hydrodynamics for exclusion processes
Let be a continuous function. We will denote by the product measure with Bernoulli marginal at site with success probability .
Our goal is to analyze the asymptotic behavior of the empirical measure with respect simple exclusion process ,
in a time scale to be chosen later. We will start the process from local equilibrium measures .
In general, one doesn’t expect to consist of independent random variables, even if initially at they are independent. Instead of computing the mean and variance, as with independent particles, we will use the martingale formulation.
To understand main ideas, let be a smooth function, not depending on time. Treating
as a function of the Markov process on the state space , we obtain a microscopic evolution equation:
where is a martingale (cf. Proposition 2.3.2) with quadratic variation
Now, the martingale is negligible in the limit: We compute that
| (2.4.1) | ||||
In the last line, we have used that the occupation variables are bounded by and that the jump probabilities are finite-range.
A calculation shows that
| (2.4.2) | ||||
where .
Exercise 2.4.1.
Verify the derivation of the quadratic variation given in (2.4.1).
2.4.1. Symmetric case.
When is symmetric, recalling the simpler form of the generator (2.2.2), a further summation by parts is possible and we obtain
where .
Now, where we recall
Hence, if , we have that
Putting these estimates together, along with (2.4.1), for symmetric , we have ‘closed’ the equation:
This suggests in the limit that the empirical measures converge in the weak sense to a solution of the heat equation
The goal of Chapter 3 is to make precise this statement for symmetric simple exclusion.
The reader will have noticed that the above display is virtually the same equation derived for independent particles. This is due to the definition of the empirical measure , which uses the occupation variables as the masses at locations . The ‘duality’ property of such functions with respect to symmetric exclusion mentioned earlier (cf. Exercise 2.2.2) allows to recover the heat equation. There are differences however in the associated fluctuations, for instance as seen in Sections 11, and 12 (specialized to independent particles).
2.4.2. Drift case.
When is asymmetric, say , as in the independent particle model, we should choose . But, one cannot ‘close’ the equation. One has to deal with the term
composed of ‘two-point’ functions . In the limit, such a term due to ‘local averaging’ should be replaced by a quadratic function of the empirical density.
Formally, one would obtain a form of Burger’s equation
There is much work on such hyperbolic mass conservation laws. In particular, there may be several solutions even starting from smooth initial data. It will turn out that the solution found by hydrodynamics is the unique ‘entropy’ or ‘vanishing viscosity’ solution. We will discuss the case in the context of ‘TASEP’, that is in totally asymmetric simple exclusion setting in when and f or in Section 6.
2.4.3. Asymmetric, mean-zero case.
In the final case when is asymmetric, but mean-zero, that is for some , but , there are other complications. The system is referred to as a ‘non-gradient’ model. Although, we should speed up time by , a second summation-by-parts as in the symmetric situation cannot be done. Indeed, multiplying (2.4.2) by , we have
Unless is symmetric, one cannot rewrite the expression in curly braces as the exact difference of a function and its translate.
Let be the standard basis in . It can be shown that the sum over dotted with can be approximated by a difference , a ‘gradient’. The function depends on . With respect to a certain ’homogenization’ of , the hydrodynamic equation can then be written down as a nonlinear Heat equation,
We refer to [15], [130], [213] for more discussion and details.
2.5. Why is it called ‘hydrodynamics’?
Finally, we comment on why the scaling limit is called a ‘hydrodynamic limit’. The origins go back to the study of identical mass particles with certain positions and momenta moving on a torus , with width , according to Newtonian dynamics:
where the energy is given by
and in terms of an even, nonnegative, smooth, compactly supported function .
In the system, the total mass, the three components of momentum, and energy are conserved. The question is how to understand in a scaling limit, as time is speeded up by , and space is rescaled by , the evolution of the local density , momenta , and energy , that is the ‘hydrodynamic flow’ of the system.
One can form the scaled empirical measures corresponding to these quantities: Writing for the point mass ,
where energy of the th particle is
On the torus , one can define a parameter family of canonical point measures , corresponding to density , velocity , and inverse temperature , which are invariant for the dynamics: Here, points are distributed in according to joint density, with respect to Lebesgue measure,
where is a normalization. The momenta are i.i.d. vectors with independent Gaussian components , independent of the positions . An infinite volume limit of these measures, as , might be taken under some conditions.
The goal is to show in some sense the limits for where satisfies an ‘Euler’ equation
Here is a matrix. The rigorous passage to such a limit in general is open! In computing , one has to ‘close’ the expressions in terms of functions of the empirical measures. Formally, one can do this and derive the form of , subject to a certain ansatz.
For instance, let us derive the equation for the formal limiting ‘local average’ density in terms of the formal limiting ‘local average’ momentum . With respect to a test function ,
Then,
To derive an equation for , however, we will need to understand how to average time-dependent quantities involving nonlinear terms such as
One would expect that these quantities could be replaced by space averages with respect to an infinite volume local equilibrium measure.
Such an ergodic theorem for the purely deterministic coupled ODE dynamics to ‘close’ the equation is however difficult to obtain. One of the best results is in [164] where a small amount of noise is added to the dynamics to make it into a reasonable Markov process so that local averaging can be done, and a rigorous limit can be proved. See also the recent developments in the solution of Hilbert’s sixth problem with respect to Boltzmann’s equation [59].
One motivation, among others, for our study of stochastic interacting particle systems is that the randomness of the microscopic motions may allow for suitable ergodic theorems. In a sense, to quote S.R.S. Varadhan, instead of ‘approximating an exact problem’, we will try to solve ‘exactly an approximate problem’. Of course, in the modeling of various phenomena, adding noise to the dynamics is often natural.
Exercise 2.5.1.
Consider the ansatz, with respect to the canonical measures , that
and the ‘pressure’
Note that
as is anti-symmetric. Derive formally from the ansatz that the equation for is given by
We note that the last equation for the energy should be
2.6. Notes
The material on martingales can be found in [66] for instance, and other places. Similar treatments of the hydrodynamics of exclusion can be found in [130] and [213].
The exclusion process–see [100] for a retrospective–has many properties which make it amenable to calculation. It has proved to be a versatile model, which can be defined on very general graphs. See [58], [147], [148], [149], [150], [190] for detailed studies using different properties.
We comment one may formulate other empirical measures of interest, those of ‘higher order’ [46], or with respect to ‘local functions’ [130][Chapter III, Lemma V.5.5].
As one can see in the sketch of the hydrodynamics for symmetric exclusion processes, requirements may be weakened on the initial measure. In fact, what is needed is a guarantee of a law of large numbers at time . The concept of ‘very weak local equilibrium’ given below (cf. Chapter III in [130]) is sufficient and somewhat general.
Let be a function. We say that a sequence of probability measures on is a ‘very weak local equilibrium’ according to profile if
for all bounded, continuous .
We remark that may be degenerate, that is supported on a single configuration, and that the sequence may consist of deterministic configurations which satisfy the law of large numbers in the definition above.
Section 3 Proof of hydrodynamics for symmetric exclusion processes
After stating the main theorem, we provide an outline of the proof and provide details on the associated steps. Such an outline can be used to prove ‘hydrodynamic limits’ in other models, in particular in zero-range processes in Chapter 4.
3.1. Statement of hydrodynamics
Recall the notation from the last subsection, with respect to -dimensional symmetric exclusion processes with finite-range symmetric jump rate , in particular the definition of the operator with , and the measure . Recall also that and stand for the process measure and expectation starting from .
Theorem 3.1.1.
Consider the symmetric exclusion process with finite-range translation-invariant jump probabilities. Then, starting from the local equilibrium measure , associated with continuous profile , the empirical measure converges in probability to the measure where satisfies the hydrodynamic equation
| (3.1.1) |
The proof is given through the following rough steps, which are then explained in more detail in later subsections.
Step 1. Consider the trajectory of empirical measures indexed in an interval of time, . Here, is time length, fixed throughout. Let be the law of the trajectory. The first step is to show that the laws are tight in the space of measure-valued right-continuous trajectories with left limits, the Skorohod space . We will show this tightness in the stronger uniform topology.
Step 2. Given tightness, we will show that any subsequential limit of must be supported on trajectories satisfying
| (3.1.2) |
This equation is similar to what was derived in the last subsection when did not depend on time. However, to input into PDE uniqueness results, we will need to derive the equation with respect to this wider class of functions . Moreover, after also showing in Step 3 below that trajectories under a limit point are absolutely continuous, we will be able to conclude the associated densities satisfy a weak formulation of (3.1.1).
Step 3. We will show that supports trajectories such that for each , is absolutely continuous with respect to Lebesgue measure on , and hence can be written as , which is a priori random. Therefore, from Step 2, is a weak solution of the hydrodynamic equation (3.1.1). By uniqueness of weak solutions to the heat equation in the class of bounded solutions, we see that is actually deterministic. In particular, all subsequential limits of converge to the point mass supported on the trajectory . Since tightness was proved in the uniform topology, this trajectory is continuous in time. [This can also be inferred from regularity results in PDE.] Hence, it can be concluded, at a fixed time , that converges weakly to a constant, the measure , and hence in fact converges in probability.
3.1.1. Proof of Step 2: Identification
We recall that we have almost shown Step 2 in the last subsection. To this end, let be a smooth function. By the development in Subsection 2.4.1, we obtain
Therefore, if is a limit point of , we obtain (3.1.2).
We now discuss more carefully some topological considerations and argue Step 1, the most difficult part. Afterwards, we concentrate on Step 3 in a subsequent subsection.
3.2. Topology and compactness
Before making calculations with respect to Step 1, we first recall some definitions and results on weak convergence, and the spaces and with a view toward characterizing when a sequence of probability measures on these spaces is tight. More details can be found in [27], [130][Section 4.1].
First, we recall that a family of probability measures on a metric space is relatively compact if any subsequence of the family has a weakly convergent subsequence. We say that the family is tight if for each there is a compact set such that all measures give at least weight to . Recall that is a complete metric space if all Cauchy sequences converge in , and is a separable space if it contains a countable dense set of points.
Proposition 3.2.1 (Prokhorov’s theorem).
Let be a complete, separable metric space. Then, a family of probability measures on is relatively compact exactly when the family is tight.
Moreover, if is a complete metric space, then a family of probability measures on is relatively compact when the family is tight.
We will consider in the following a generic complete, separable space with metric . Often, will be the space of finite measures on equipped with the metric
where is a countable dense set of functions in , the space of continuous functions on the torus endowed with the ‘sup’ metric between and . We may take . Here, is a complete, separable metric space. Moreover, a set is relatively compact exactly when .
At other times, may be with usual Euclidean metric.
Now, we consider the space with the ‘uniform’ distance,
It is known that this space is a complete separable space, and therefore amenable to Prokhorov’s theorem.
Since our basic building blocks in our study of hydrodynamics involve jump processes, the space of continuous trajectories is not sufficient for our purposes. We will often focus on the space of right-continuous trajectories with left limits. Unfortunately, the uniform distance will not make this space a complete, separable metric space. Define to be the set of strictly increasing continuous functions of into itself, and
Then, define the Skorohod distance between elements in as
In some sense, the Skorohod distance compares two trajectories allowing small variation both in space and in time. To contrast, the uniform distance only compares variation in space. With the Skorohod distance, is a complete, separable metric space.
How to characterize compact sets in these spaces? Consider the following moduli of continuity:
where refers to a partition such that for .
Then, belongs to exactly when , and belongs to exactly when .
Exercise 3.2.2.
Relate the two moduli by showing that
Exercise 3.2.3.
Show, when and , that but .
Compact sets in these spaces, since they are complete, are characterized as follows.
Proposition 3.2.4.
A set belonging to or is compact exactly when
- (a)
is relatively compact in for each .
- (b)
where on and on .
We remark, when , the above characterization reduces to the Ascoli-Arzela condition with respect to equicontinuous families of trajectories on the finite interval , and (a) can be replaced by ‘ is relatively compact in ’.
Recall now Exercise 3.2.2. We now have the following claim.
Proposition 3.2.5.
A family of probability measures on is tight exactly when
- (1)
For each , the distributions of under are tight.
- (2)
For every , .
Moreover, a sufficient condition for is
- (2’)
In condition , replace with .
Exercise 3.2.6.
Observe that any limit point of satisfying is supported on continuous paths. This is a well-known property; see [27].
It is not so easy to work with directly. However, when , one can understand a family of probability measures on by their actions on smooth functions on .
Proposition 3.2.7.
A family of probability measures on is tight if the distributions of under for are tight in , that is satisfying (1) and (2’) in Proposition 3.2.5, for each .
3.3. Proof of Step 1: Tightness
We are now back to considering the tightness of for which are elements of . From Proposition 3.2.7, we need only show for a smooth function , tightness of the distributions of for which are elements of .
From Proposition 3.2.5, we need to show conditions and . Condition is the simplest, and follows straightforwardly as for each ,
Then, is a tight sequence in as the sequence is uniformly bounded in with full probability.
Verifying Condition , and therefore tightness in the uniform topology, is a little more involved. Since
we need only show condition (2’) for each term separately. The initial term does not contribute in this respect.
We now bound the second term. Condition (2’) follows as
where is the range of the probability .
For the third term, to treat condition (2’), we first recall Doob’s inequality (cf. [63]): For a martingale , and , we have
3.4. Proof of Step 3: Absolute continuity
To show is absolutely continuous, observe for any continuous function that
since at most one particle is allowed per site. Hence, since the function
is continuous with respect to the Skorohod topology, any limit point satisfies, with full probability (see Portmanteau theorem [27][Chapter 1]),
Hence, any limit point of is supported on trajectories with the property that is absolutely continuous for all .
Exercise 3.4.1.
Detail the use of the Portmanteau theorem in the previous paragraph.
At this point, as noted in the initial strategy, for each , can be written as where may be random. However, by Step 2, is a weak solution to the hydrodynamic equation, which has a unique bounded solution. Therefore, is deterministic, and where is the hydrodynamic density. We remark when then is actually a ‘classical’ solution given in terms of a Gaussian kernel convolution with .
What we have shown is that the law of , where initial configurations are distributed according to , converges to the point mass on trajectory . To conclude convergence at a fixed time , we note that in Step 1, tightness of was obtained in the uniform topology. Hence, the trajectory in the support of the limit is continuous for all . For , let be the projection function,
Now, is not continuous on . However, since the limit is supported on a continuous trajectory, is continuous on the support of . Now it is known that if and is a continuous function almost surely on the support of , then (see [27][Chapter 1]). In other words, the projection converges in law to the point mass at , and therefore also in probability.
3.5. Notes
There are other proofs of hydrodynamics for symmetric exclusion processes, for instance using correlation functions and ‘duality’ as in [58] or superexponential estimates as in [213]. We have followed mostly the treatment with some caveats in [130], following the strategy, in the symmetric exclusion context, expounded in [106], a classic in the field. It is worth noting by this method we have shown existence of weak solutions to the PDE (3.1.1).
Section 4 Entropy, Dirichlet form, and averaging principle in zero-range processes
We now consider another well-studied ‘mass conservative’ interacting particle system, the so called ‘zero-range’ model. Unlike for symmetric simple exclusion, the evolution equation for the mass empirical measure does not close, and one must invoke an ‘averaging principle’ to approximate a nonlinear term in terms of a function of the empirical measure. Here, we develop the ‘entropy method’ for this rigorous replacement in the associated hydrodynamics. To this end, notions of entropy and Dirichlet form, as well as the ‘basic coupling’, are introduced before defining the zero-range process. At the end, we give an outline to derive the associated hydrodynamical equation.
The following Chapter 5 is devoted to proving steps in the outline, including the famous so called ‘1-block’ and ‘2-block’ lemmas introduced in [106].
4.1. Relative entropy
Let be a countable space. For probability measures and on , define the relative entropy of with respect to as
| (4.1.1) |
where the supremum is over bounded functions . Recall, as before, stands for the expectation under .
Since the expression
is invariant for all constants , the above supremum can be taken over bounded nonnegative functions .
Exercise 4.1.1.
Show that is nonnegative, convex and lower semi-continuous in the argument (where if converges weakly to ). Hint: For nonnegativity, substitute equal to a constant.
Lemma 4.1.2 (Entropy inequality).
For all bounded functions and , we have
Also, when , for , one has
Exercise 4.1.3.
Prove this lemma. Hint: Multiply and divide by .
Lemma 4.1.4.
We have implies that . Also, if , then
Moreover, when is finite, we have .
Proof.
If , there is an such that but . In the variational definition of , we can insert the function to see . Hence, .
We now suppose and prove the equality in the display. By simple truncations, approximate in (4.1.1) by the supremum over functions which vanish except for a finite number of points in . We now find the maximum of
over functions supported on where and . Indeed, the functional is concave and takes its maximum where its gradient vanishes. Computing the gradient, we obtain the maximum occurs when
for each . We can add a constant to , noting and for , and not change this equation–the mapping is extended to functions which are constant off . Hence, we may choose so that
With this choice, we obtain for , with convention . As and , since vanishes, we have the upper bound
We obtain the lower bound by inserting into the variational formula for , and then taking . Hence, implies that .
The second statement follows finiteness of . ∎
4.2. Entropy with respect to Markov chains
Let be a continuous time Markov chain on a finite space . Let be an invariant measure. Recall and stand for the semigroup and generator of the process. We write in the following .
Lemma 4.2.1.
For probability measures such that , we have , and for .
Proof.
We show first that is absolutely continuous with respect to , and therefore by Lemma 4.1.4. As , and , we have . Then, by the invariance , we conclude
showing the absolute continuity.
Now, write in form , where is convex. Since
we have to finish,
Exercise 4.2.2.
One can update Lemma 4.2.1 to countable state space by showing that if . Perform this deduction when is countable.
We now find the rate of decrease in the entropy.
Lemma 4.2.3.
Under the assumptions of Lemma 4.2.1, we have
Proof.
Recall from the forward and backward equations that . Then, as there is no problem to interchange limits and sums in finite state space,
The first term on the right-hand side equals
after interchanging the sum on and .
On the other hand, the second term vanishes:
Here, is a function on for each , and so given is invariant. ∎
4.2.1. Dirichlet forms
For a continuous time Markov process on a finite state space with invariant measure , define the Dirichlet form on functions by
We now define the adjoint by the relation for . Then, a simple computation shows
The operators and may also be defined. If is reversible, . In particular, is a reversible operator, , and is anti-symmetric in that for all . Recall, also as we have seen before,
Then, by straightforward calculation, as also ,
Since is reversible, by interchanging and in the above sum and averaging the two expressions, we obtain the first of the following formulas. The second follows by anti-symmetry of .
Lemma 4.2.4.
We have
| (4.2.1) | |||||
We have the following properties of the Dirichlet form. For a nonnegative function, let .
Lemma 4.2.5.
We have is nonnegative, convex, and lower semi-continuous. When the chain is irreducible and , then is a constant function.
Proof.
Nonnegativity follows from the expression (4.2.1). After squaring out terms, convexity follows by Jensen inequality applied to the crossterm. Clearly, the Dirichlet form is continuous, and therefore lower semi-continuous, in its argument as the space is finite.
When , we have for all where . If the chain is irreducible, all values of must be the same. ∎
Remark 4.2.6.
Naturally, the above notions of Dirichlet form extend to countable state Markov chains, with respect to compactly supported functions , and to those in in the domain of the generator .
4.2.2. Connection between entropy and Dirichlet form
Let be a continuous time finite space Markov chain with invariant measure . Recall in Lemma 4.2.3, that the derivative of equals
We now bound this expression in terms of the form .
Lemma 4.2.7.
We have
Proof.
The argument follows from an application of the inequality for , which we leave to the reader. ∎
4.3. Zero-range models
We now discuss the ‘zero-range’ system of interacting random walks on the -dimensional torus with locations. Informally, particles interact infinitesimally only with those on their sites through their jump times, hence the name ‘zero-range’.
More specifically, at a vertex with particles, a particle displaces by with rate where again we assume is a finite-range translation-invariant jump probability on and is a function such that and for . An alternate description is that each vertex has its own exponential clock which rings at rate when there are particles at the vertex. When rung, one of the particle is selected at random and then it displaces according to .
Formally, consider where denotes the number of particles at at time . The process on the countable configuration space is a continuous time Markov chain with semigroup and generator
| (4.3.1) |
where is the configuration obtained from by moving a particle from to :
4.3.1. Invariant measures
As with the exclusion process, there is an associated family of invariant measures depending on particle density. Let be the product measure with common marginal on given by
Here, is the normalization which converges for .
An interesting point is that when , this is the model of ‘independent’ random walks considered in Section 1 when there is no interaction. In this case, of course are Poisson measures.
We now index the family in terms of ‘density’. Let be the density, . One can check that the function is a strictly increasing function which may or may not diverge when , depending on the structure of . In any case, for we can invert and define
the product measure which places a mean number of particles at each location in . The term is sometimes called a ‘fugacity’ or ‘chemical potential’. One may calculate that is also the mean of the rate function: . In the following, we drop the superscript ‘’ and write to simplify notation.
Exercise 4.3.1.
For bounded, local , that is supported only on a finite number of variables , show
where and is the configuration with exactly one particle at .
Lemma 4.3.2.
Let be such that . Then, is an invariant measure for process on .
Moreover, does not depend on and is the unique invariant measure on when is irreducible.
Proof.
Recall the expression for in (4.3.1). We need only show that for all bounded local functions . Since is translation-invariant, it is doubly stochastic, .
The second statement for , the restriction of to , as for independent particles, follows by noting is a closed, irreducible set when is irreducible. ∎
Local equilibrium. In the following, with respect to a continuous nonnegative function , define
to be the (inhomogeneous) product measure with marginals over sites . Sometimes is referred to as a ‘local equilibrium’ measure. Again, we will drop the superscript ‘’ and write to simplify notation.
Remark 4.3.3.
We comment that divergence of as allows to define and for any and , useful in the formulation of the hydrodynamic limit.
A sufficient condition is that diverges as . This is the case when is an increasing function. See Lemma II.3.3 and Section II.3 in [130] for an argument and more discussion.
Exercise 4.3.4.
Suppose . Then, . By explicit calculation, show that diverges as .
4.3.2. Connection with entropy
The next results will be used to bound a certain ‘entropy production’, discussed in Section 5, when starting from . Since here is countable but not finite, we now update the results in Subsection 4.2 to this context.
Lemma 4.3.5.
Consider and such that . Then, , and for .
Proof.
Exercise 4.3.6.
Extend the definition of with respect to (4.2.1) in the current setting:
| (4.3.2) |
Lemma 4.3.7.
Under the assumptions of Lemma 4.3.5, we have, for ,
Proof.
By Lemma 4.3.5, and therefore by Lemma 4.1.4, we may write
where for each . Decomposing on the number of particles in , the entropy equals
Recall, here is a finite set, and . Let also and .
On , note by mass conservation that . Then, on , . Therefore, inputting these relations,
Since is itself a relative entropy of the measures and on , it is nonnegative. Also, the relative entropy for each . Note also that , as is a finite space, is differentiable.
Hence, given finiteness of , we may write it as a sum of two finite, nonnegative terms, one of which, , does not depend on . Therefore, we are able to write the difference
The right-hand side, by the arguments of Lemmas 4.2.3 and 4.2.7, since is an invariant measure on , is evaluated as
Since the summands, as , are all negative, we may interchange the sum on and the integral. The right-hand side above equals . We obtain the desired statement as . ∎
Remark 4.3.8.
Since and so is continuous in , by Fatou’s lemma, one may bound the upper derivative:
See also discussion in Appendix I.1.9 in [130].
4.4. Basic coupling
We discuss a useful coupling in the context of zero-range processes when is an increasing function. First, we define the notion of ‘stochastic domination’ on the partially ordered set . We say two probability measures , on a countable set are ordered, , if for all coordinatewise increasing functions .
There is an interesting characterization of ordered measures for whose proof we refer to [147][Theorem II.2.4].
Proposition 4.4.1.
We have exactly when there is a bivariate distribution on such that marginally variables and are distributed according to and and .
In this next result, does not have to be increasing.
Lemma 4.4.2.
Let . The marginals exactly when .
Proof.
If , since is increasing, we have from which one deduces .
For the converse, since linear combinations of are dense in the set of increasing functions, it is enough to show that
for all . This is equivalent, after cancellation, to showing
which is equivalent to
The last inequality will be shown if term by term the inequality is true. But, since .∎
We now show the existence of the so-called ‘basic coupling’ for the zero-range process when is increasing. Another name for such an existence is that the zero-range process is ‘attractive’.
Proposition 4.4.3.
When is increasing, there is a joint process on , starting from where , such that at all later times , and marginally are the zero-range processes starting from and respectively and a.s.
Proof.
Let be the joint probability on such that a.s. (Proposition 4.4.1). Consider the Markov process on with initial distribution and generator
One can see, by inserting a function of coordinate only or of coordinate only, that the marginal processes are as desired.
To check that a.s., recall the joint process is a continuous time Markov chain on a countable state space. At time , we have arranged that . At the next jump time , by the specification of the rates, we see that still . Hence, the joint process remains ordered for all time . ∎
Exercise 4.4.4.
There is another way to show that for all a.s. (cf. Theorem II.5.2 in [130]). Compute that . Then, with denoting the process expectation starting from , is increasing in : . Hence,
The next exercise will be useful in the later proof of hydrodynamics. Recall that and are well-defined when is increasing by Remark 4.3.3.
Exercise 4.4.5.
Let . Show that , and consequently , when is increasing.
We will also need later an estimate on if is a Lipschitz function, for all .
Lemma 4.4.6.
If is Lipschitz, then is also Lipschitz.
Proof.
Let . Then, by the basic coupling,
As we have seen, the zero-range model allows several ‘closed form’ calculations. The following is another one of these. We will not use this exercise in the sequel.
Exercise 4.4.7.
Calculate that where is the variance of under .
4.5. What is the hydrodynamical equation?
We will start the process from initial configurations distributed according to . Denote, as before, for fixed , that and are the distribution of and its expectation, when is governed by .
We now informally compute the generator action to guess the hydrodynamic behavior. Recall stands for the empirical measure
For a smooth , consider the equation
where is a martingale. One may compute the quadratic variation,
which shows, easily if is bounded, that the martingale is negligible in the limit.
We now compute
Here, again when is symmetric (or mean-zero), and otherwise.
When is symmetric, similar to symmetric exclusion, one can change to and average. The right-side equals
where we recall .
The trouble now is that this weighted average of functions does not close in terms of the empirical measure. The main point of ‘hydrodynamics’ is to approximate this average by a function of the empirical measure.
What should it be? In the microscopic scale, near point there is lots of local particle movement. One expects the distribution of particles in an neighborhood of at time to be in ‘local’ equilibrium given by the nearby density, such as where
| (4.5.1) |
and is a small parameter.
Then, one might expect that
where we recall .
Now, this ‘local’ density, for , can be re-expressed in terms of the empirical measure:
| (4.5.2) |
which if is further approximated by
Putting this together, one ‘closes’ the equation and obtains, after first and then , that
which is the weak formulation of
where is as defined before, and .
4.6. Assumptions and statement of hydrodynamics
We will make the following assumptions to help simplify the proof and introduce main ideas. We remark however that hydrodynamics has been shown for much more general zero-range processes; see the later Notes for details on extensions.
We will assume that is a Lipschitz function, bounded and increasing, with invariant measures at all densities :
- •
for
- •
,
- •
.
Note by Remark 4.3.3 that as , and so is defined for all densities .
A standard example is , under which is Geometric, that is for and .
Also, to reduce notation, we assume is nearest-neighbor and symmetric:
- •
for coordinate unit vectors , so that .
Recall that stands for the process measure when starting in .
Theorem 4.6.1.
We have for all and that
where is the unique weak solution of the hydrodynamic equation
4.6.1. Strategy
We now separate the proof into two rough steps.
Step 1. We will again consider the measures which govern the trajectories of empirical distributions on . The first task is to show that is tight, and all limit points are concentrated on trajectories of measures with densities .
Step 2. We show that all limit points are supported on weak solutions to the hydrodynamical equation in . It is a result in PDE that such weak solutions are unique. Hence, the measures converge to a single limit point. One concludes convergence at a fixed time, as for simple exclusion processes. The development on entropy, Dirichlet forms, and coupling will be useful in proving Step 2 in Section 5.
4.7. Notes
The zero-range process was introduced in [205]; see also [5], [68]. It has been a good vehicle to study a variety of phenomena, beyond hydrodynamics, such as condensation and metastability (cf. [8], [17], [136], [137], [154]).
Much of the development of entropy and its use in various applications, including hydrodynamics, stems from the work of Guo, Papanicolaou, and Varadhan (GPV) [106]. On the other hand the basic coupling, and its use in particle systems analysis was advanced by Liggett [147].
In terms of hydrodynamics, one may weaken the assumptions on . In particular, the increasing or boundedness assumptions are not needed when is mean-zero; see [130] and the Notes of Section 5. See the Notes of Section 6 for remarks when is asymmetric and not mean-zero.
There are at least five different proofs of the hydrodynamic behavior when say is symmetric. One way, the ‘entropy’ method [106], to be outlined in Section 5, shows that the limiting hydrodynamic density satisfies a weak form of hydrodynamical PDE. Another way, the ‘relative entropy’ method [221], assumes that a classical, smooth solution exists, and then shows that the relative entropy of with respect to vanishes; see Section 6. One may use a Hopf-Lax formula, and the basic coupling, to prove in the hydrodynamical limit for systems with drift [189], [6], [10]. A method using a logarithmic Sobolev inequality and compensated compactness ideas can be employed in [80], [81]. Also, the hydrodynamic limit may be seen as a ‘gradient flow’, and found via Gamma convergence methods [61], [71], [101]. The first two methods are also discussed in [130]. See also [58] for a treatment of the GPV ‘entropy’ method, and [83], [213] for treatments of the ‘relative entropy’ method.
Section 5 Entropy method for hydrodynamics of zero-range processes
A rigorous proof of hydrodynamics (Theorem 4.6.1) is given for a class of zero-range models through the ‘entropy’ method of [106]. The main idea is that the drift computed from the generator action, an average of a nonlinear function of the occupation variables, may be understood in the scaling limit as a function of the limiting empirical density.
This sort of ergodic theorem, which proceeds in two replacement steps, the ‘1-block’ and ‘2-block’ lemmas which introduce additional scales, is facilitated by estimates on the Dirichlet form of a nonstationary Radon-Nikodym density function, and local central limit theorem asymptotics.
5.1. Proof of Step 1: Tightness and absolute continuity
Recall the notation from the Section 4, in particular the development and assumptions from Subsections 4.5 and 4.6. As before stands for the expectation under , and and denote the process measure and expectation when starting in .
Analogously as for simple exclusion, the first step can be divided into the tasks:
- •
Show that the measures on which govern starting from local equilibrium are tight.
- •
Show that all limit points of are supported on trajectories with densities .
Lemma 5.1.1.
is tight.
Proof.
Considering Proposition 3.2.7, we need only show items (1) for each that is tight on , and
Now, by the basic coupling and Exercise 4.4.5, we have that where . Hence, (1) follows as the first moments, , are uniformly bounded:
To establish (2’), we consider, as we did for the exclusion process, the drift and martingale terms separately. For the drift term, since is bounded, we have the uniform bound on ,
Noting the form of the quadratic variation with in Subsection 4.5, the proof is similar to that as for exclusion processes. ∎
We now establish the following characterization of limit points.
Lemma 5.1.2.
All limit points of are supported on trajectories with densities such that, for a constant ,
Before beginning the proof, we remark that absolute continuity of the limit trajectories means that particles cannot pile up to form point masses. Given the process is ‘attractive’, that is the basic coupling holds, one can bound the probability of large piles in terms of estimates where is an absolute bound on the hydrodynamic density. Without attractiveness, the proof is harder and we refer to [130] for a general argument.
Proof.
We show first that the trajectories are absolutely continuous under a limit point . Let be a smooth function. Note, for and for all large , by attractiveness and truncation bounds, that
The last term is an increasing function of , which we will show is small. The role of is to introduce a truncation. Let . Therefore, by the basic coupling, Chebychev’s inequality and invariance of , for , we have
Hence, by simple estimates,
and by weak convergence, since the function is continuous in the Skorohod topology and for closed sets,
Since is arbitrary, we have a.s. under that all trajectories satisfy
Now, note by the tightness estimate (2’) already proven in the previous lemma that all trajectories under are continuous in time in the underlying topology (e.g. is continuous for ). Then, by choosing to approximate for , we obtain , that is is absolutely continuous for each .
To show the display in the lemma, recall (4.5.1) and consider the bound
via the basic coupling and Schwarz inequality
Hence, as is a limit point, by Fatou’s lemma again, we have
where is the ball of radius around . Taking limit on , and another application of Fatou’s lemma along with Lebesgue’s differentiation theorem (a.e. is a Lebesgue point), we finish the proof. ∎
5.2. Proof of Step 2: Identification given replacement homogenization
We now supply, modulo a replacement estimate, the proof of Theorem 4.6.1. Note that is a bounded function as is bounded. Recall also the notation in (4.5.1). Let be the shift operator, and .
Theorem 5.2.1 (Replacement).
For , we have that
where
Notice that is a bounded function as is bounded.
Let us now see how this replacement allows to finish the proof of hydrodynamics.
5.2.1. Proof of Theorem 4.6.1
First, we establish an equation that the densities under a limit point must satisfy. Let . Then, following the derivation when did not depend on time in Subsection 4.5, we obtain
where by Doob’s inequality,
since is bounded.
On the other hand, by Theorem 5.2.1, since is continuous, we have
Now, again the quantity in absolute value in (5.2.1) is a continuous function of in the Skorohod topology. Hence, any limit point , by Fatou’s lemma, satisfies
where .
But, since trajectories have densities under , and is bounded and continuous, by Lebesgue’s differentiation theorem and dominated convergence,
Finally, we obtain therefore is supported on trajectories satisfying
which are weak solutions to
5.3. Proof of the replacement homogenization
After preliminaries and some reductions, we prove Theorem 5.2.1 at the end of this subsection. Since we start out of a stationary distribution, we have to understand the nonstationary contribution. Controlling the entropy and Dirichlet form of the associated Radon-Nikodym probability density is the first step.
For , let be a reference measure throughout the rest of Section 5. Let also be the Radon-Nikodym density at macroscopic time . Let . By Lemma 4.3.5, we know for some constant . Since the time scaling is , we have by Lemma 4.3.7 for each that
Let . Then, by convexity of the entropy and Dirichlet form, we have
| (5.3.1) |
For later reference, the following estimate will allow to introduce a truncation corresponding to the number of particles in a region.
Lemma 5.3.1.
For , we have
Proof.
By Markov’s inequality, the basic coupling as is an increasing function, and Schwarz inequality, we have the left-hand side is bounded by
5.3.1. Reductions
Consider the following bound of the integrand in Theorem 5.2.1, by adding and subtracting terms. Write
| (5.3.2) |
The first term on the right-hand side introduces the scale and more averaging: As is bounded and is uniformly continuous, it vanishes when and . If were , it would be of order .
The second term, bringing an absolute value inside the sum, is bounded by where we observe
Moreover, we may introduce a truncation by Lemma 5.3.1, so that what remains to bound is for the expectation of
In the third term, as is Lipschitz, we bound it by
We may further write the -window term in terms of an average of disjoint -window terms for as
where are centers of the cubes in the decomposition and the last sum is over at most variables in possibly uncompleted -windows. One may also add into this sum the contributions from the -windows neighboring , say those with centers within of . Then, noting the -overcount of these ‘edge’ items, the third term is further bounded by
The -process expectation of the last term is of order since by the basic coupling . We may also introduce truncations, for by Lemma 5.3.1, so that we need only bound the expectation of
| (5.3.3) |
5.3.2. Statements of and -block lemmas
Recall the definition of the density near (5.3.1). By the development in the previous subsection, to bound the second term in (5.3.1), it is enough to estimate
Since we know from (5.3.1), choosing , it will be sufficient to show the following.
Proposition 5.3.2 (1-block lemma).
We have for all that
where the supremum is over nonnegative densities with .
The third term in (5.3.1) is handled similarly. Noting the formulation in (5.3.3), it is enough to show the following limit.
Proposition 5.3.3 (2-block lemma).
We have for that
5.4. Proof of 1-block lemma
To prove Proposition 5.3.2, the main idea is that the Dirichlet form of the density vanishes as . In some sense, the optimal density is almost constant. Things then reduce to a standard ergodic theorem or law of large numbers associated with i.i.d. random variables distributed according to .
To make this strategy precise, define
Then, by translation-invariance of ,
| (5.4.1) |
Observe that is a translation-invariant density.
Let now
Denote . Since depends only on variables where , we have the right-side of (5.4.1) equals
where is the product measure over sites in with marginal .
Now, recall that the Dirichlet form has explicit form (4.3.2). For such that , let
Then, ; note there are edges in .
Define also the Dirichlet form corresponding to the dynamics restricted to with invariant measure .
note there are edges in .
Lemma 5.4.1.
We have that
Proof.
Note and are averages: For instance, is a conditional expectation. Then, the first and second inequalities follow as the Dirichlet form is convex. The equality follows as and the measure are translation-invariant: Namely, for and . The last line follows from the bound . ∎
Now, taking into account (5.4.1) and the display after, and Lemma 5.4.1, we need to show
| (5.4.2) |
where the supremum is over densities with respect to .
In fact, since is uniformly bounded below for all large because , in (5.4.2), we may replace the integrating measure by which supports only a finite number of configurations. Also, as
we may view the supremum as over the collection of nonnegative with respect to which are uniformly bounded . The corresponding collection of subprobablity measures , being on a compact space, is tight and has a converging subsequence by a form of Prokhorov’s theorem.
Hence, to evaluate (5.4.2), for fixed , we consider a sequence in which approaches the limit supremum. Densities can be found on which the supremum value is well approximated. By tightness, we can find a subsequence where converges to for which necessarily . Hence, it is enough to show
Now, given and , we know that is constant on the configurations such that for . Therefore, as , decomposing along such ‘hyperplanes’, it is enough to show that
| (5.4.3) |
The left-hand side expectation is the same as
Here, we removed the restriction of the measure to since it does not matter. Notice also that the integrating measure is the canonical measure on with particles. Hence, by the form , the parameter of the underlying grand-canonical measure does not matter, and can be chosen as we like, say .
We may rewrite the above display as
| (5.4.4) | ||||
It is an exercise to see that the denominator is bounded away from .
Exercise 5.4.2.
Observe by a local central limit theorem, namely an expansion in the characteristic function, that
Now, the variables in (5.4.4) are mean-zero and in . Hence, by an application of Schwarz inequality and Exercise 5.4.2 again, we have that (5.4.4) is at most of order which vanishes.
This concludes the proof of Proposition 5.3.2. ∎
5.5. Proof of 2-block lemma
The argument for the proof of Proposition 5.3.3 is similar to the proof of the ‘1-block Lemma’ Proposition 5.3.2. Now, we have to control the differences of averages in two separated blocks of width . To compare these averages, we need again to show that the localized Dirichlet form of the density function is small. Since the jump probabilities are nearest-neighbor, to localize on the two blocks of width , we will need to extend the Dirichlet form and dynamics to include a long jump from one block to the other. In this way, the localized system can mix, and the difference in the averages will wash out. There is a cost for this extension which however can be overcome.
As before, with the ‘-block Lemma’, the first step is to estimate the display in Proposition 5.3.3 in terms of the averaged density over shifts . We need to show that
Here the supremum is over translation invariant densities .
Next, we localize to the union of the two blocks and . Let be the product measure over with marginal . Let also where . Then, as before, the above display reduces to
for say .
The question now is how to treat the Dirichlet form . If we use the argument for the ‘-block Lemma’, then a density in the limit would be constant on each hyperplane for where . But, the constants could be different on each block!
The trick is to add a (non-local) Dirichlet form bond corresponding to a jump say from to . Any such jump from a site in to would work. Then, if the Dirichlet form localized to , including this extra bond term, vanishes, the density would take the same constant value on both blocks. In terms of dynamics, the zero-range process associated to this modified Dirichlet form is irreducible on configurations in .
For , define
noting Exercise 4.3.1. By itself, is the Dirichlet form for the dynamics which moves a particle from to and back according to zero-range dynamics, and as such shares all the convexity and other Dirichlet form properties that we used in the proof of the ‘1-block Lemma’.
By the estimate and adding and subtracting several terms,
where corresponds to a nearest-neighbor path from to in steps.
We will apply the above estimate to the translation invariant density . Since , we have that . Now, if and , since the separation between and is of order , we conclude .
Let now where . By convexity, we have
| (5.5.1) |
for fixed and all large .
Hence, it is enough to show for each constant that
Here, the supremum is over densities with respect to .
But, at this point, the proof can follow the method as for the ‘-block Lemma’. We just point out that in the above display that, as and share no terms, with respect to , they are independent. This was the reason to avoid the nearest-neighbor cubes in the original decomposition of in Subsection 5.3.1. This finishes the proof of Proposition 5.3.3. ∎
5.6. Notes
We have followed the scheme of Chapter V in [130], although there are differences. The reader will have noted that translation invariance is much used in the above derivations. We remark that there are ways to handle non-translation invariant dynamics if the inhomogeneity is ‘slowly varying’; see [52], and [95], [69], [76], [138] when the particles move in a random environment.
Also, we comment that diffusive scaling was crucial to prove the ‘-block’ estimate; see (5.5.1). In asymmetric models, with hyperbolic scaling , although the ‘-block’ result still holds, the ‘-block’ estimate is not generally available. Such a ‘-block Lemma’ is useful for the ‘relative entropy’ method discussed in Section 6 to deduce hydrodynamics, at least for short times, in asymmetric processes.
Finally, we comment, beyond Zero-range models, the ‘entropy’ method works well in symmetric rate processes () with conserved quantities that is of ‘gradient’ type, e.g. allowing a twice sum-by-parts evaluation
for some local function (depending only on a finite number of variables ) to be homogenized, and invariant measures indexed to the conserved quantities with sufficient decorrelation properties (cf. for instance [85]). As alluded in Subsection 2.4.3, diffusively scaled ‘nongradient’ models can also be analyzed by versions of the ‘entropy’ method (cf. [15], [130], [213]).
Section 6 Relative entropy method and hydrodynamics of TASEP
We discuss the ‘relative entropy’ method in the context of hydrodynamics of the totally asymmetric simple exclusion process (TASEP) in . The underlying notion is to measure how far the distribution of the process at time is away from a specified product measure on . If distance, measured in terms of relative entropy, is not so far away, then hydrodynamics will follow from calculations with respect to independent variables governed by . Such methods and ideas have proved useful for other processes and in other problems.
6.1. Statement of hydrodynamics and sketch of the argument
Recall the notation in Chaper 2, with respect to asymmetric exclusion processes on . As before, denotes the expectation under , and and the process measure and expectation when starting in .
When and allows movement only to the nearest right site, that is and for otherwise, the process is known as the totally asymmetric simple exclusion process or TASEP for short. To reduce notation, we will concentrate in the following on TASEP, although calculations may be extended to finite-range asymmetric processes in with non-zero drift.
As mentioned in Section 2, the nontrivial time scaling is when , the ‘hyperbolic’ or ‘Euler’ scale, and the hydrodynamic density satisfies . In the context of TASEP, as , the equation simplifies to
| (6.1.1) |
Solutions to this one-dimensional Burgers-type hyperbolic conservation law on are not necessarily unique. However, when starting from smooth initial profile , bounded above and below on ,
it is known that a unique classical smooth solution, also bounded above and below by and , with continuous bounded derivatives, exists up to a short time , which we now fix. We mention one may specify uniquely a ‘physical’ (or by other names, ‘entropy’ or ‘viscosity’) solution for all times by imposing that the solution satisfies extra conditions; see [67] as a general reference and for a comprehensive discussion.
Recall the empirical measure on :
Define, for , with respect to the hydrodynamic density, the product measures on :
Suppose that initially the process is distributed according to . In terms of the process semigroup , let be the distribution of starting from .
We will show the following form of the hydrodynamic limit.
Theorem 6.1.1.
Let be a function. Then, for , we have the convergence in probability,
where satisfies (6.1.1) with .
Proof.
The sketch of the argument is in a few steps.
Step 1. Since for all , we have . Noting Lemma 4.1.4, we will show that the relative entropy between the distribution of and is . This step is the main part of the proof. It takes the -block lemma as input, and is discussed in the next subsection.
Let now and
An exponential bound may be found as the variables are independent under with means . Therefore, given the relative entropy estimate in Step 1, the theorem follows. ∎
Exercise 6.1.2.
Show the ‘Chernoff’ exponential bound in Step 2, by computing the moment-generating functions. The inequality will be useful.
6.2. Proof of Step 1: entropy estimate
Let denote the product invariant measure of TASEP with density . Recall the form of in (2.2.1) (in and , otherwise). One may compute that the -adjoint operator of is the generator of the asymmetric exclusion process, speeded up by , where particles jumps left, that is when and otherwise:
The probabilities satisfy the forward equation (cf. Exercise 6.2.1),
| (6.2.1) |
However, is no longer a product measure, the process having mixed things up from the initial condition . One feels though from the intuition given in Section 2 that the process may not be far from . The relative entropy is a convenient, and analytically tractable measure of how far apart they are.
Exercise 6.2.1.
Perform the computation of . Show also satisfies (6.2.1) by differentiating where the density and is a bounded function. The relation for all will be useful.
Step A. We first update the conclusion of Lemmas 4.2.3 and 4.2.7, one of the differences being that is not the invariant measure.
Lemma 6.2.2.
We have
Proof.
Write
Note that as . Also,
Then,
| (6.2.2) |
Now, consider the relation
for . Then, as , with and , we observe
which further equals
We may insert this expression into (6.2.2). Dividing/multiplying by again and moving to on the other side, we obtain the expression desired. ∎
We remark that for the ‘relative entropy’ method in the sequel, we will not use the first term, a ‘Dirichlet’ form expression, bounding it above by . However, the the full force of Lemma 6.2.2 has been useful in other problems; see the Notes.
Hence,
Observe
Hence, under the condition and , corresponding to , we have by Taylor approximation that
Therefore, cancelling the ‘’, we obtain
On the other hand,
Recall . Then, by Lemma 6.2.2, dropping the negative ‘Dirichlet’ term, and dividing by , we conclude
Step C. Now, we invoke a form of the -block lemma in Section 5. We state it in in the context of TASEP, although it can be generalized to and other processes. Recall that .
Proposition 6.2.3.
Let be continuous. Let be a local function and . Then, for ,
Sketch of proof.
Consider the -block shown in Section 5, namely Theorem 5.2.1 with time-dependent replacing , replacing , and limits replacing . The same scheme can be followed here: With respect to the first two terms on the right-hand side of (5.3.1), continuity of in the space variable handles the first term, and as it is Proposition 5.3.2 takes care of the second term.
A difference is that the time scaling previously was . Although this was important for the later -block estimate Proposition 5.3.3, for the -block estimate, the Euler scale suffices. Indeed, we will now obtain the Dirichlet form of the density is of order , instead of when . Such an estimate is enough to obtain the -block limit. These computations are left to the reader. ∎
Observe that
and that it vanishes when . Hence,
with respect to a constant .
Subtracting and adding terms in the summation (6.2) is negligible: Indeed,
Therefore, we have that
Step D. We now estimate the expectation above, via use of the entropy inequality (Lemma 4.1.2): Let . Then, with ,
The variables are -dependent. To gain some independence, divide the interval into blocks of width , where the possible last block may be of length less than . We may separate the sum over as a sum over such that and where refers to the set , and also a sum over the remaining indices in the possible overflow block .
Let . The sum corresponding to the overflow block can be bounded . We may also center where . Note for .
Then, for each , since are independent, by Hölder’s inequality,
Step E. We now apply a concentration inequality for subgaussian random variables. We say is -subgaussian if
for all .
Lemma 6.2.4.
Let be -subgaussian. Then, if ,
In our context, the variable , with respect to , is subgaussian with parameter
Exercise 6.2.5.
Since on , we have . Then,
when .
Hence, for such a ,
which vanishes as and then .
Step F. Putting things together, we observe that the entropy is bounded as
We have taken the initial distribution to be , in which case (an estimate with respect to the initial condition would also suffice). By Gronwall’s bound, we conclude as desired that
as and . ∎
6.3. Notes
We have followed the scheme in [213], although there are differences; see also Chapter VI in [130]. The hydrodynamic limit presented in Theorem 6.1.1 may be extended to all times , by different methods, where is the ‘physical’ solution mentioned earlier. See for instance [182], when starting from a ‘step’ initial profile , and [10], [179] for general initial data.
However, in general, hydrodynamics for all times has not been shown for processes with drift in Euler scale that do not satisfy the ‘basic coupling’ or are not ‘attractive’ in the sense given in Subsection 4.4. It is an open problem to show such hydrodynamics. Here, TASEP can be verified to be ‘attractive’, used in [10], [179]. See also [199] for hydrodynamics of long-range asymmetric exclusion and other processes.
The ‘relative entropy’ method is due to H.T. Yau [221]. Other treatments of the method in different contexts may be found in [130], [83]. Instead of using the concentration inequality for the last argument, one may invoke large deviations bounds as in [130] for instance.
While the ‘entropy’ method of Guo-Papanicolaou-Varadhan [106] presented in Section 5 yields existence of a weak solution to the hydrodynamic partial differential equation, the ‘relative entropy’ method shows uniqueness of classical solutions of the hydrodynamic equation: Although a requirement for the method is that a priori some smoothness of the hydrodynamic PDE solution is known.
Section 7 Construction of particle systems in infinite volume
We construct, using the method of Liggett and Spitzer [151] and Andjel [5], zero-range particle systems on . The approach is to take a limit of the processes restricted to large but finite sets. Such a construction would also hold for exclusion and related models. Other construction techniques are mentioned in the Notes subsection.
7.1. What does it mean to construct a process?
In Section 1, on finite or countable state spaces , ‘construction of the process’ meant determining the transition probabilities of a Markov chain, in terms of infinitesimal rates, from which probabilities could be found for . One could then write down then the semigroup operator , acting on bounded functions, and understand its connection to the generator in terms of backward and forward equations.
On more exotic spaces , to construct a Markov process usually means building a semigroup , that is an operator, with the ‘Chapman-Kolmogorov’ property for , acting on a collection of functions on , and relating to a generator via backward and forward evolution equations. In this context, one can usually associate, by Kolmogorov’s extension theorem, a probability , with initial condition , on the Borel sets in . Sometimes, we will want the semigroup to have certain properties, such as ‘strong continuity’ or the ‘Feller property’, but this is not guaranteed (cf. [145][p. 247]).
For the particle systems we have studied, exclusion and zero-range models, we would like to extend the background space to . Then, the configuration space would be for exclusion systems, and for zero-range models. With respect to exclusion systems, is compact, which is helpful and allows different ways to construct the process. However, for zero-range models, since is not compact, some care must be taken and assumptions on the parameters, namely the rate function and transition probability , should be made.
Interestingly, if is in general unbounded, but Lipschitz, not all configurations may be ‘allowed‘. In other words, we will construct the zero-range process on a strict subset , so that once begun in , the process will stay in . Part of the reason for this restriction is that if the rate function is large, particles jump faster, and the process may be influenced from particles at ‘infinity’; see the end of Section 2 in [5] for an example.
On the other hand, if is bounded, then the process can be constructed on the full by a different method [111].
7.2. Zero-range model and statement of results
We will assume in the following that satisfies and for , and is Lipschitz: There is a constant such that
- (LIP)
In particular, we do not assume in this Section 7 that is increasing.
Also, we will take the jump probability on ( and ) to be irreducible and such that for all . Note that we do not assume that is translation-invariant, e.g. , or that is finite-range, that is when for for some .
We now specify the allowed configuration space . For , let
where is the probability that a single particle reaches in steps from . Observe that
| (7.2.1) |
Exercise 7.2.1.
Show (7.2.1).
Define for the norm
The collection of allowed configurations will be
| (7.2.2) |
Examples of configurations in include those with only a finite number of particles, for instance: We say a configuration is finite if the number of particles .
We now define a class of ‘Lipschitz’ functions on which we construct the process. We say that is Lipschitz if there is a constant such that
for all . Let be the smallest such constant .
Denote by the collection of all Lipschitz functions on . Note that includes ‘simple’ functions, those say depending on a finite number of variables and taking on a finite number of values. Moreover, trivially, belongs to with .
Define the operator , which we will identify later as our generator, on functions in and by
Lemma 7.2.2.
The operator is well defined for and , and
Proof.
We now come to the main theorems.
Theorem 7.2.3.
There exists a semigroup on , which has specification for and finite configurations in terms of the countable-state process.
Moreover, for and , the semigroup satisfies
and . Also,
We observe, when there are only a finite number of particles in the system, and are the semigroup and generator of a countable state Markov chain.
More properties are given in the following result.
Theorem 7.2.4.
We have for and that
- (i)
and
- (ii)
- (iii)
.
7.2.1. The meaning of Theorem 7.2.3
Recall that ‘simple’ or cylinder functions, for and , , are Lipschitz functions. Since, for a given , one may approximate by finite configurations such that as , by Theorem 7.2.3, the semigroup may be computed from that of countable state Markov chains: . Therefore, the finite dimensional distributions of for a given initial configuration may be identified.
Then, by Kolmogorov’s extension theorem, there exists a probability measure on Borel sets in (not necessarily for the moment on !) for each and such that
| (7.2.3) |
for cylinder functions .
Lemma 7.2.5.
For , we have , and therefore the measure concentrates on .
Proof.
Lemma 7.2.6.
For , we may identify for or that
Proof.
When , one can approximate by nonnegative simple functions (belonging to ). Then, by monotone convergence, one may take a limit in (7.2.3), and define in this way.
7.3. Construction estimates on a finite cube
The strategy to construct an infinite volume semigroup is first to obtain estimates on the semigroup and generator which are well defined when the space is finite, and then use these estimates to define as the volume grows. Throughout this subsection, the underlying space will be a cube of width :
Let and be the (countable-state) zero-range process corresponding to a transition probability on , with no transitions from to .
We now bound the means of such a process at times .
Lemma 7.3.1.
For , and , we have
Here, is the -step transition probability from to .
Proof.
The argument is by coupling the zero-range process to a continuous-time multitype branching process on with generator
Here, is the configuration which adds a particle at site to . Note also the sum is over where . An inspection of the formula reveals that is such that each particle at gives birth to a new particle at rate . This new particle is then displaced by with probability . Alternatively, each particle gives birth to two new particles before it dies; one is kept at the birth location , and the other is displaced by .
The coupling of and as follows: Since by (LIP), , whenever a zero-range particle displaces from to , one may arrange also that the branching process gives birth at and creates a new particle at . In this way, if the two processes are started from the same initial configuration, then for all and .
In particular, when coordinatewise, we have for all .
To finish, we need only calculate . Observe that
Then, satisfies where is the transition matrix. Hence, . In particular,
Lemma 7.3.2.
For and , we have and .
Proof.
Consider the basic coupling in Subsection 4.4, the joint Markov process on generated by
Let be the coupled process semigroup. Both marginals are zero-range process, and if initially , then the coordinatewise ordering is preserved, .
Now write
We give now an estimate of the derivative of the right-hand side. Note
and . Then, adding together the positive and negative parts of , with an application of (7.2.1), we obtain
Hence, since , by comparison, we obtain
which finishes the proof. ∎
Let now and be zero-range processes on according to jump probabilities and that do not allow transitions from to .
Lemma 7.3.3.
For , we have
Lemma 7.3.4.
We have
Exercise 7.3.5.
Lemma 7.3.6.
For , we have
| (7.3.1) | ||||
7.4. Extension to and proofs of the construction theorems
Let increase to . For a transition probability on , define the modified transition probabilities:
Note that does not allow transitions from to .
It will be useful to observe that . Hence, satisfies (7.2.1) with constant instead of .
Let and be the zero-range semigroup and generator corresponding to on .
Proposition 7.4.1.
For and , we have converges uniformly on sets of bounded , sets of with bounded , and sets of functions with bounded .
Proof.
We show that forms a Cauchy sequence with the uniformity properties. Consider and , where . Both are transition probabilities on . We have that the integrand in (7.3.1) is bounded by
Here, we used (7.2.1) several times.
Hence, the integrand is dominated. Given vanishes, the integrand vanishes pointwise for each , as .
Plugging into (7.3.1), we see that the convergence is uniform as desired. ∎
7.4.1. Proof of Theorem 7.2.3.
By Proposition 7.4.1, we may now define, for and that
Moreover, by Lemma 7.3.2, since satisfies (7.2.1) with constant , we have
Consider now the two properties stated in Theorem 7.2.3. We first establish the semigroup property for : Namely, for . Already this property holds for . It is sufficient to show for and that
The first limit follows from the uniform convergence in Proposition 7.4.1, given uniformly in .
For the second limit, for fixed and , by construction, we may find by Kolmogorov’s extension theorem (as discussed in Subsection 7.2.1) a probability on such that and, for ,
Now observe
Hence, by dominated convergence,
To establish the integral property, note that it holds for the finite-volume semigroup and generator :
We can already pass to the limit on the left-hand side. To pass on the right-hand side we need to show that
| (7.4.1) |
and that the the integrand is dominated. Domination holds as
| (7.4.2) |
We leave to the reader to show (7.4.1). ∎
7.4.2. Proof of Theorem 7.2.4
We will use the previous estimates to show the items in Theorem 7.2.4.
Proof of (i). From Lemma 7.2.2 and Theorem 7.2.3 (see also (7.4.1) and (7.4.2)), we have that
Now, plugging into the integral expression in Theorem 7.2.3, and integrating, we complete (i).
Proof of (ii). If we can show that is continuous at , then (ii) follows from the integral formula in Theorem 7.2.3. Now, from part (i), is continuous at . Write
To show continuity of , we dominate
The right-hand side bound is summable:
Hence, the desired continuity follows from dominated convergence.
Proof of (iii). Write
The first equality is from part (ii). The second equation is from the semigroup property proved in Theorem 7.2.3. The third equality is from dominated convergence and Lemma 7.2.6, writing for . To justify these steps, we note the pointwise convergence is established already in part (ii). The domination and belonging to follows from part (i): For , using ,
This completes the argument. ∎
7.5. Notes
The material follows [5] which is almost always cited when working with zero-range processes. On the other hand, other construction methods exist. For instance, if is bounded, one can construct the semigroup by the Hille-Yosida theorem as in [111]. [145] gives another way, along the lines developed above, with different estimates, to construct the semigroup. See also [14] for yet another construction.
The exclusion process can also be constructed on in the way above, but can be done also using the Hille-Yosida theorem; see [147]. In this setting, local functions, those which depend on a finite number of variables , provide a natural core for the generator.
Section 8 Invariant measures: ergodicity and extremality
We consider the invariant measures in zero-range processes on . Given the mass-conservative dynamics, there are several invariant measures corresponding to different particle densities. To determine them, since the set of invariant measures is a convex set, it is useful to identify extreme points of this set. In particular, if an invariant measure is an extreme point, the process run under it has interesting ergodic properties; see Subsection 8.3.
Our focus here will be to show that a class of product invariant measures are extreme points. With respect to zero-range models, many things are known, but there are still several open questions.
8.1. Invariant measures for zero-range processes
To simplify the discussion, we will assume in the remainder of Section 8 that is translation-invariant, but not necessarily finite-range. However, a good theory exists when is not translation-invariant, remarked upon in the Notes subsection. We will also keep the assumptions on made in Section 7: Namely, , for , and the (LIP) condition, for .
Recall the discussion of invariant measures for the zero-range process on a torus in Subsection 4.3. Our initial goal will be to extend these notions to . Recall the marginal where
for , and form the measure
As before, is a strictly increasing function of , and hence the inverse exists. One then defines
for .
Recall , and from Section 7, the norm where , and subset .
Lemma 8.1.1.
We have and hence the measure fully charges , that is
| (8.1.1) |
Moreover, we have .
Proof.
Write
where , since by translation-invariance of .
Moreover, by Schwarz inequality,
We now update the finite-volume definition of an invariant measure to . Recall the Lipschitz functions and the semigroup of the process on from Section 7. We will say that a probability measure on is an invariant measure if
for all bounded Lipschitz functions . Let denote the convex set of invariant measures for the process on .
Our first main theorem is the following infinite volume invariance.
Theorem 8.1.2 (Invariance).
For , we have .
8.2. Proof of the invariance of in infinite volume
It will be convenient to give a generator characterization of invariance of a measure.
Proposition 8.2.1.
Let be a probability measure on such that . Then, the following are equivalent:
- (1)
for all bounded
- (2)
.
In what follows, we will only use that ‘(1) implies (2)’, and so we will prove this part of the result.
Proof of Proposition 8.2.1: (1) (2).
For bounded , we have by Theorem 7.2.3 that
Since , we have is supported on . For and , by Lemma 7.2.2 and Theorem 7.2.3 we have and
Moreover, as is bounded, is a bounded, Lipschitz function.
By the domination, we may interchange the integrals by Fubini’s theorem. Then, given (1), we have
and so conclude (2). ∎
Remark 8.2.2.
To prove the converse, ‘(2) implies (1)’, we refer the reader to [5][Lemma 2.9].
We now assert that
| (8.2.1) |
is an invariant measure of the process with translation-invariant transition probability restricted to the finite set . Indeed, this follows immediately from (the proof of) Lemma 4.3.2 in Section 4.
Proof of Theorem 8.1.2.
We observe the following properties of :
(a). First, is a transition probability: for all . If , it is trivial. If , noting that , we have . Moreover, has no transitions from to its complement.
(b). In addition, for all : If , the claim is trivial. If , write
The last quantity equals if is doubly-stochastic, the case if is translation-invariant.
(c). Also, for all . Eventually, . The claim follows as .
Now, by the comments near (8.2.1), as nothing moves away from when starting in , we have restricted to is invariant. Then, we have for bounded where is the generator for the zero-range process on according to transition probability . Therefore, to show , it is enough to show
One may write, noting (LIP), , when and , and for , that
For fixed ,
Hence, since and is doubly stochastic, the penultimate display vanishes as by dominated convergence. ∎
8.3. Extremality, harmonicity and ergodicity
We digress for the moment to an abstract setting. Let be a space with Borel sets . Let be a Markov process on with process semigroup , acting on bounded functions (and well-defined extensions). We say here that is an invariant measure if for all bounded functions on (and therefore functions) and , we have . Let be the process measure on the path space with initial distribution .
Definition 8.3.1.
We say is an extremal invariant measure if the the following property holds. When for and invariant probability measures and we have , then .
We note if the process has only one invariant measure , then of course it is extremal. The import of the definition comes when the process is reducible in some way. For instance, in finite state Markov chains, with exactly two irreducible components and , the extreme invariant measures are exactly the unique invariant measures supported on and respectively.
One might ask why is it useful to know when an invariant measure is extremal. It turns out there is an interesting connection with harmonic functions and shift-ergodicity.
Exercise 8.3.2.
With respect to invariant measure , noting the representation , show that is an contraction. Let be the adjoint of . Show also that is an contraction for each .
We will say that is harmonic, if for all .
Lemma 8.3.3.
The space is the direct sum of and the closure of .
Proof.
Let be perpendicular to all functions in , the closure of . Then, for all and . This means for all . Therefore, as is a contraction, . Hence, and or .
On the other hand, let . Then, . Since , by a similar argument as above, we conclude and so and . Hence, . ∎
We have the following ergodic theorem due to Von Neumann.
Proposition 8.3.4.
Let and define to be the projection onto the subspace . In other words, is the conditional expectation with respect to the time-shift invariant sets .
Then, we have the convergence as ,
Proof.
Decompose where . Clearly, . Since can be approximated by in , we need to understand the contribution of the term , in form for some and . We see
in from Hölder’s inequality as is an contraction. ∎
We come now to the main result of the subsection. Recall stands for the process measure when starting in .
Proposition 8.3.5.
Let be an invariant measure. All are equivalent:
- (a)
For sets , a.s. .
- (b)
is ergodic: For each , , a.s.
- (c)
is extremal.
Note that part (b) is the same as ‘time-shift ergodicity’ or in other words that the shift invariant -field is trivial: Sets satisfy or .
Proof.
‘bc’ Let be an invariant measure whose path measure is ergodic. Write for and invariant measures and . Let now be a bounded function. Then, by (b), as , converges to in and therefore also in . Moreover, by Proposition 8.3.4, converges to in . Hence, a.s. and taking expectation, . This gives for and therefore .
‘ab’ Let be an invariant measure and suppose that is not ergodic. Then there exists an such that is not constant -a.s. Let be such that , where .
Now, as -a.s. and is a positive contraction taking to , we have that -a.s.: Indeed, first, as is a positive operator, , so that, as is an contraction, we have . Note also . Then, we have , and so -a.s. . Therefore, is harmonic. Further, if are harmonic, then is harmonic. Correspondingly, is harmonic. Of course, is harmonic. All of this gives that is a sequence of uniformly bounded harmonic functions. The limit, as , is which is therefore harmonic by dominated convergence. Hence, by (a), is constant, a contradiction.
‘ca’ Let be such that -a.s. and for . By this relation, the process begun on stays in and, if begun in stays in , with -probability . Then, we have that and are distinct invariant probability measures such that . Therefore is not extremal. ∎
8.4. Extremality of
We now address the ‘extremality’ of in the convex set of invariant measures . To do this rigorously, we will need to extend the process so that the extended semigroup and generator act on functions, not just Lipschitz functions . In this way, we can fit into Subsection 8.3, and in particular will be an invariant measure in this setting.
However, to present the main ideas we assume this extension has been done, and that adjoints can be taken. In the next subsection, we detail the extension to a process. Recall is the symmetrized transition probability.
Theorem 8.4.1.
When is irreducible on , the invariant measure is extremal.
Proof.
The idea is simple. We will need to understand the Dirichlet form of the process. In the next subsection, we show for , the domain of the extended generator , that the associated Dirichlet form satisfies
Now, let be a bounded harmonic function, that is when , a.s. Hence, as the limit,
(trivially) exists a.s., we conclude and , a.s.
Therefore, , and so all summands vanish when . In particular, as exactly when ,
Consequently, is invariant to the motion of particles! So, by irreducibility of and countability of , one has
where is the configuration which exchanges values and .
Therefore, is finite permutation invariant. Since is a product measure with i.i.d. marginals, by Hewitt-Savage law, is constant a.s.-. Inputting into Proposition 8.3.5 shows extremality. ∎
8.5. Extension to an process
We first extend the semigroup defined on to . We define the concept of a Markov semigroup on .
Definition 8.5.1.
A family of linear operators on is a Markov semigroup if (a) ; (b) for , is right-continuous in on ; (c) satisfies the semigroup property for ; (d) for ; (e) for and nonnegative .
Lemma 8.5.2.
on extends by continuity to a Markov semigroup on .
Proof.
For , we have , and hence by Lemma 8.1.1. Also, for , we may write . Therefore, by Jensen inequality, . Then,
Given that simple functions, those supported on finitely many variables and taking finitely many values, are contained in , we observe is dense in . Therefore, for , we may define as the Cauchy limit of for a sequence of simple functions which converge to in . Such an extension, given Theorems 7.2.3 and 7.2.4, fulfills the definition of a Markov semigroup on , details left to the reader. ∎
Exercise 8.5.3.
Verify that satisfies the definition of an invariant measure given at the beginning of Subsection 8.3 with respect to .
There is a one-to-one correspondence between Markov semigroups and associated generators by the Hille-Yosida theorem; see [66].
Proposition 8.5.4 (Hille-Yosida Theorem).
For a Markov semigroup on , define
Then, if , we have and .
Let be the generator associated with semigroup with domain .
Lemma 8.5.5.
We have , on extends on , and is the closure of on , that is the graph of is the closure in of the graph of .
Proof.
Since by Lemma 8.1.1, and on extends on , we have by dominated convergence for that
in as . Hence, , , and is an extension of on .
We now give the formula for the Dirichlet form. Recall the symmetrized transition probability .
Lemma 8.5.6.
For , we have that
| (8.5.1) |
Proof.
We note, by the last line in the proof of Lemma 8.5.5, that bounded functions in form a ‘core’ for . First, the formula (8.5.1) may be obtained for bounded , given that is explicitly defined and , and is left as an exercise.
We now extend the representation to . Let be the right-hand side of the display in the lemma. For , take bounded so that and in . Then
by Fatou’s lemma. Therefore, and in particular, for . However also,
which vanishes as . Hence, . With the inequality , one may bound
Hence, writing , we may conclude to finish the proof. ∎
8.6. Notes
In these notes, we have followed [5] for the invariance of and [194] for the extension to (see also Section IV.4 in [147]) and extremality of ; see these papers for extensions to non translation-invariant transition probabilities. See also [181][Section 2] for more on the connections between extremality, harmonicity, and ergodicity.
A natural question is what are all the extremals of the process. In [5], it is shown in when is an increasing function that are all the extremals! When is not translation-invariant, also in when is increasing, all the extremals are found–in this case, the invariant measures are not necessarily translation-invariant [5]. When is positive-recurrent, the extremals concentrate on configurations with a finite number of particles [215], [5].
For simple exclusion on , since local functions, those depending on a finite number of variables , form a core, one may verify that the Bernoulli product measures for are invariant for the infinite volume process (as compared to the finite-volume one in Proposition 2.2.1); see Chapter VIII in [147]. One may also use the arguments given here for the zero-range process. Extremality of can also be seen from the arguments for Theorem 8.4.1.
However, for zero-range and exclusion models, there remain many open problems! In particular, determining all the extremal invariant measures in all model settings is still not fully resolved. See [147], [194], as well as [3], [16], [48], [31], [32], [91], [126], [152] for more references, and discussion.
Section 9 Additive functionals and central limit theorems
The occupation time of a set, for instance, is an additive functional of a process. One would like to know the first order and second order behaviors, namely the LLN and CLT for these objects. This is an old subject, and in stochastic particle systems where the configuration space is difficult, with interesting scalings and limits.
We focus here mostly on symmetric exclusion processes on , when starting from an invariant measure. After a general introduction, we discuss ergodic and fluctuation behaviors, when starting under an invariant measure. The Kipnis-Varadhan CLT, which has wide application, not only to particle systems, is stated and proved. Comments on the behaviors in asymmetric simple exclusion are also made.
9.1. Basic problem and ergodic theory
Consider a Markov process on a state space . Let be an invariant measure for the process. Let be an function. Recall that stands for the expectation under , and and for the process measure and expectation when starting in .
The basic problem is to determine the behavior as of
In the following, we call an ‘additive functional’ since under has the same distribution as under .
If is extremal, then, as we have seen in Proposition 8.3.5, the process started from initial distribution is ergodic with respect to time shifts, and hence (by Birkhoff’s theorem)
This covers the case when is unique, for instance, in positive-recurrent Markov chains. When is null-recurrent, different behaviors may emerge (see [168]).
The next order question is to ask about the fluctuations: Find the limits of
One expects a Gaussian limit, with sufficient mixing. In the following, we will assume is extremal and that is already centered, that is to simplify expressions.
Consider the example of finite-state irreducible Markov chains. We may view and as elements of . Note that for all where is the generator of the process. Since is the unique invariant measure, . The orthogonal complement of this null space is : For , we have . Then, as is invariant.
As a consequence, given , we have and hence belongs to the range , and for some function . That is, solves a ‘Poisson’ equation with respect to the generator.
We may now write
| (9.1.1) |
where
is a martingale with respect to natural sigma-fields. The quadratic variation, as we have observed in Subsection 2.3, equals
whose mean is
where is the Dirichlet form.
The second term in (9.1.1) is of order since is a function on finite state space, and hence bounded. The first term, however, can be treated with a martingale central limit theorem.
Proposition 9.1.1.
Let be a mean-zero martingale with stationary, ergodic increments such that . Then,
Exercise 9.1.2.
To complete the argument, note the martingale has stationary and ergodic increments when the process is started from . Hence, by Proposition 9.1.1, we conclude
9.2. Occupation functionals in exclusion systems
Consider the -dimensional exclusion process with semigroup and generator ,
lifting the finite-volume definition in Section 2. Such an infinite-volume system on may be well constructed; see the Notes in Section 7. In the following, we will assume the jump probability is finite-range and translation-invariant.
Let be the Bernoulli product measure with marginal success probability on . These measures can be shown invariant and extremal for the process, and when is symmetric they are also reversible; see the Notes in Section 8. Compare also with the finite-volume discussion in Section 2. In the following, we will work with the -exclusion processes.
Consider now the ‘additive functional’ question for occupation functions . Given is extremal and therefore ergodic with respect to time shifts (cf. Proposition 8.3.5), we see that the LLN behavior starting under is clear: For ,
both in and -a.s. by Proposition 8.3.4 and Birkhoff’s ergodic theorem.
However, the fluctuation behavior of , with respect to is different depending on the dimension , the symmetry/asymmetry of the jump probability , and the density . A more or less complete theory is known, except for some interesting asymmetric cases in which connect with ‘KPZ’ class phenomenon, and ‘mean-zero asymmetric’ cases.
We first consider the types of variances that might be obtained. What is the variance of , what are its limits, and how does the limit depend on ? For a general local function , that is one that depends only on a finite number of variables , let us try to get a formula for . Write, using stationarity, that
By reversibility of , the semigroup is self-adjoint, and we may express
Recall, .
Then, when the variance is scaled by , its limit exists, possibly infinite, and is in form
The last equality follows from monotone convergence. This formula is the familiar ‘sum of correlations’ expression with respect to sums of correlated random variables.
9.2.1. Symmetric and duality relations
In symmetric exclusion, any mean-zero, local function can be written as
Here, the decomposition is over ‘degrees’, that is, functions which are the product of centered, normalized occupation variables.
This basis is quite central to symmetric exclusion because of the ‘duality’ relation. Namely, for a function
where , we have
where is the transition probability of -particle simple exclusion from a configuration to configuration in time .
One way to prove this relation is to observe that
where is the transition rate of -particle symmetric exclusion; see Exercise 2.2.2.
At this point, we can use the duality relation, and independence of coordinates under , to evaluate the term
Now, when is the centered occupation function, for all , and . Then, the variance is calculated as
where is the return probability of a random walk according to jump probabilities !
From local limit theorems, when the transition probability is simple, for the standard basis , we have the following asymptotics.
Proposition 9.2.1.
For , under symmetric simple exclusion processes,
9.2.2. Asymmetric
By asymmetric, we mean that . The case that is not symmetric but , the ‘mean-zero’ asymmetric situation is not considered here; see the Notes for citations.
As in symmetric exclusion, one can understand the general asymptotics by calculating the term . However, the ‘duality relation’ is no longer valid in the form given, since the generator does not preserve the ‘degree’ of a function, that is applied to a function of coordinates is no longer a linear combination of coordinate functions.
However, for , write
The basic coupling, with respect to exclusion, couples two copies of the exclusion process starting from configurations where for and , . Then, has generator
The process always majorizes , with exactly one discrepancy, whose position we label .
The dynamics of is as follows: Infinitesimally, it displaces by with rate . The term corresponds to the discrepancy, or ‘second-class’ particle as it is sometimes known, moving by its own intention, and the term refers to when a particle at would move to location in which case by the basic coupling, accedes and takes the place .
Then, we have that
Hence, by monotone convergence the limit exists,
In mean value, substituting for , we see that a ‘mean’ infinitesimal drift is
This leads to the conjecture that
which has been proved; see the Notes.
Proposition 9.2.2.
For , with respect to asymmetric exclusion with drift , we have
when or when .
However, when , we have . In , we have .
Remark 9.2.3.
We remark the orders expected when in are and respectively. These orders connect with certain KPZ class phenomena, and are open to verify. See the Notes for more discussion.
9.3. Central limit theorems and norms
Our goal now is to prove asymptotic normality of when . Except for the two open cases in Proposition 9.2.2, a central limit theorem also holds when is normalized by the square root of its variance, although we do not discuss these results here; see the following Notes.
We focus on the symmetric case where the Kipnis-Varadhan CLT applies. The Kipnis-Varadhan theorem is a CLT for additive functionals of a Markov process, reversible and ergodic, on general state space .
Theorem 9.3.1.
Consider an Markov process begun with ergodic, reversible invariant measure . Let be an function such that . Then,
A functional CLT/invariance principle can also be proved under the assumption of the theorem, but we do not discuss this extension here.
The main tool to prove Theorem 9.3.1 is to approximate by a martingale, and then to use Proposition 9.1.1. The idea is to try to write . However, this is not possible in general. If it were possible, one could use the results in [26]. A resolvent type equation always holds however, that can be worked with. The following, which can be proved by this approach, is sufficient for the martingale approximation.
Proposition 9.3.2.
Under the assumptions of Theorem 9.3.1, there is a martingale with stationary and ergodic increments such that
where
9.3.1. and norms
Before proving Proposition 9.3.2, some definitions will be useful. The generator of exclusion is well defined on the core of local functions. Define, for local functions and , the semi-norm by
Note that, since is a nonnegative self-adjoint operator,
After modding out by functions with , define the space as the completion with respect to .
Now, we define for a local function that
Again, after modding out by functions , define as the completion with respect to . When , is a bounded operator (bounded by from resolvent formulas–see the next subsection), and
Both spaces , are Hilbert spaces where innerproducts are given by polarization:
These two norms are dual to each other: For local functions,
When , and are called the and norms respectively. Note that .
We note there is another formula for of note:
| (9.3.1) |
We remark, when is not self-adjoint, there are useful notions of and spaces; see [18][Lemma 2.1], [133] for instance.
Sometimes the norm is referred to as the ‘variance’ norm. Since is increasing as , and , we have
The last evaluation of the limit follows as is nonnegative since is self-adjoint.
In particular, we observe the following.
Lemma 9.3.4.
We have and so
9.3.2. Step 1: Resolvent calculations
Exercise 9.3.5.
Show that . This bound is part of the standard theory of ‘resolvents’.
Let us multiply the resolvent equation by and take expectation:
Then,
which gives immediately that
uniformly in .
Also, since , we observe is a bounded operator, with bound . Hence, for all ,
Now, by the uniform boundedness principle, one can find a subsequence which converges weakly to an element :
For local, and as local functions are dense in , we have
and since , we conclude weakly in . Therefore,
| (9.3.3) |
9.3.3. Step 2: Strong approximations
At this point, we would like to claim that
so that and satisfies a ‘weak’ Poisson equation. However, the weak convergence shown, as there are limits in in both components of , is not strong enough to conclude this relation immediately.
Nevertheless, we will show that the convergences can be strengthened as follows.
Proposition 9.3.6.
We have
- (i)
.
- (ii)
There is a such that strongly in .
Also .
Proof of Proposition 9.3.2.
Given (i) and (ii) above, to be shown later, we now finish the proof of the Proposition 9.3.2. Write
where
is a martingale and is the remainder.
We claim now . Indeed, we need only evaluate
| (9.3.4) |
First, as and , we have and so
Second, to finish,
As a consequence, we observe
as from (ii). Hence, being a Cauchy sequence, in . Therefore, (with respect to the natural filtration of ) is a martingale with ergodic and stationary increments such that
and
Now, the error is handled as follows: By the equation , we see that a limit holds as in . Hence, and
By Schwarz inequality, using that is invariant, we have by (i), choosing , that
as .
On the other hand,
as .
Finally, we need to identify the variance: By (ii) and Lemma 9.3.4,
9.3.4. Step 3: Proof of Theorem 9.3.1
We now have a representation
The error vanishes in and .
The martingale central limit theorem, Proposition 9.1.1, finishes the proof. ∎
9.3.5. Step 4: Proof of Proposition 9.3.6
We now recall Mazur’s theorem (see [148][Lemma 4.38]). A proof is given at the end.
Lemma 9.3.7.
In a Hilbert space, if weakly, there is a convex combination of the which converges strongly to , that is .
Now, let be a a convex combination of such that strongly in . Then, as weakly in and is linear, weakly in . In fact, since , as is self-adjoint, we see that is a Cauchy in , and hence converges strongly to in .
Then,
Now, returning to the subsequence , by lower semicontinuity and weak convergence of in , and (9.3.2), we have
The same calculation can be repeated with ‘’ replaced by ‘’.
Therefore, we conclude , which means strongly in , and also that .
By the same arguments, on any subsequence of , a further subsequence can be found so that . Hence, , showing part (i).
To show part (ii), we need to show that the limit is obtained on any subsequential limit of strongly in . To this end, as before, suppose is a subsequence converging strongly to in . We conclude as before that strongly in . We now show that which will finish the claim.
Write,
since noting (9.3.2). Now, weakly in (cf. before (9.3.3)). Hence, via the penultimate equation, and the relation , we may approximate by and by . Then, by the last equation, we see that the right-hand side vanishes. Therefore, in . This finishes the proof of Proposition 9.3.6. ∎
Proof of Mazur’s Theorem. Let us assume the limit is without loss of generality. From the weak convergence, one can find iteratively a subsequence with such that
Then,
which vanishes as . Hence, we can take for . ∎
9.4. Notes
Aside from the original paper [131], treatments of the Kipnis-Varadhan theorem, and applications to additive functionals and tagged particles, can be found in [133], [148]. There have been generalizations of the Kipnis-Varadhan theorem to nonreversible situations [202], [212], [159], some of which has been used here, albeit in the reversible setting. In the nonreversible case, as we have seen, it may be that but where the norm is with respect to the symmetrized operator , e.g. for in . However, a non-asymptotic bound holds; see Lemma 12.2.1, Remark 12.2.2. An open problem of interest is to show in the general nonreversible case, with a condition such as , that a CLT holds.
The asymptotic variances in Proposition 9.2.1 and associated CLT’s were first proved in [129]. For more general additive functionals, when is not necessarily the occupation function of a site, CLT’s and diffusive variance criteria are proved in [203] for reversible mass conservative systems.
For the asymmetric process with a non-zero drift, Proposition 9.2.2 is proved in a combination of papers [18], [191], [195], [196], and associated CLT’s are shown in [193]. In particular, the paper [195] shows that norms are comparable across finite-range exclusion processes with the same drift. For instance such comparisons would allow to lift variance calculations from asymmetric simple exclusion processes, e.g. ‘ASEP’ or ‘TASEP’, as in [191], to more general finite-range processes.
The ‘second-class’ particle , introduced in Subsection 9.2.2 in the context of ASEP is of interest by itself. When the process starts under , in , scales as [13], [173], and in scales as [222], both ‘super-diffusive’. In , the variance is diffusive of order ; see [143]. However, when the initial condition is a ‘step’ condition (to the left and right with different densities) in , the variance may have different orders [72], [75]; see also [1] and references therein. Moreover, among other properties, the second-class particle in tracks the location of microscopic ‘shocks’ in the system [72], [73]. See also [2] and references therein for a ‘speed’ process coupling between second-class particles with different priorities.
The superdiffusive conjectures in Remark 9.2.3 regarding , for and , rely on a ‘Gaussian ansatz’, that the local limit terms in and in scale like the power . Verification is an open problem; see however [12] for a weak form of the ansatz in .
In the mean-zero asymmetric setting, variance orders are the same as in the symmetrized process, and CLTs hold when the norm is finite. This would include occupation functionals in . See [212] for a Kipnis-Varadhan-type CLT using a ‘sector inequality’ for the generators of mean-zero processes. In mean-zero asymmetric processes, fractional Brownian motion limits for properly scaled, centered occupation variables have been shown [96]. For mean-zero asymmetric systems, it is open to show a fluctuation limit for the properly scaled, centered occupation functionals. See [21], [22] for more discussion of these problems in simple exclusion and zero-range processes, when starting from an invariant measure.
Section 10 A tagged particle in symmetric exclusion on
Understanding the motion of a tracer particle, as it interacts with others, is a basic applied concern. We consider here the problem in the finite range symmetric exclusion processes when started under an invariant measure . Except for the case in dimension when the jump law is nearest-neighbor, the tagged particle motion is diffusive and converges to a Brownian motion. In the exceptional case, it can be shown that the motion is subdiffusive and converges to a fractional Brownian motion of Hurst parameter .
First, we discuss LLNs. Then, the diffusive behavior is considered using the Kipnis-Varadhan theorem. Finally, we discuss an exceptional subdiffusive case.
10.1. Tagged problem setting
Consider the -dimensional exclusion process with finite-range translation-invariant jump probability , as discussed in Section 9. We will assume that the symmetrized jump probability is irreducible. The system consists typically of an infinite number particles. Let us identify one of them initially and let be its position at time . To fix things, suppose it starts at the origin initially, .
How to capture the evolution of , say LLN and CLT statements? With respect to its own history, it is not Markovian because of influence of the other particles. However, if we consider and together, then the joint process is Markovian with generator
It will be convenient to consider ‘Lagrangian’ coordinates, or those in the reference frame of the tagged motion. Define , that is for all where is shift by . Then, the joint process is also Markovian and has generator
where is the configuration which exchanges the values and , displacing the particle at the origin, namely the tagged particle by , and then shifts the reference frame to this new origin.
In particular, an explicit description of is given by
Both joint processes can be constructed with a core of local functions. Interestingly, the process by itself is a Markov process. This can be seen as the generator acting on functions of alone,
does not depend on , and hence the associated semigroup also does not depend on . The process ‘drives’ the process in that can be recovered in terms of reference frame shifts.
Define as the count of shifts of of displacement up to time . Then,
| (10.1.1) |
where the sum may be restricted to in the support of . For instance, in one dimension, when is nearest-neighbor, is the number of right shifts minus the number of left shifts.
Moreover, it is not difficult to find invariant measures for generated by .
Lemma 10.1.1.
The Bernoulli product measure , conditioned to have a particle at the origin, is invariant for . When is symmetric, is reversible.
The Dirichlet form can be computed also on local functions:
| (10.1.2) | |||||
where is the symmetrized probability, , assumed to be irreducible on .
Lemma 10.1.2.
is extremal for the process .
Exercise 10.1.3.
Exercise 10.1.4.
Now, as before, with zero-range processes, we can extend the process to an process. The problem now is to understand the LLN and fluctuations for in terms of the formula (10.1.1).
10.1.1. Martingales for
Each count , as with a Poisson process, can be compensated by its intensity, , to form a martingle
with quadratic variation
Then,
is a martingale.
How to verify the above statements? Technically, we should have formed the generator for the process . Then, with , we may compute and to see that
and
Although is not a bounded function, one may apply truncations to bring it in to the domain of the generator , while noting , a Poisson distributed r.v., to help remove the truncations.
Now, we may write
| (10.1.3) |
The first term, is another (vector valued) martingale with quadratic variation, for , given by
| (10.1.4) |
Remark 10.1.5.
We also note, since is extremal, by Proposition 8.3.5, has stationary and ergodic increments and, by (10.1.4),
| (10.1.5) |
The fourth moment of may also be bounded via a generalization of the Burkholder-Gundy-Davis inequalities:
where is the compensator of the jumps process and measures the size of a jump; see [116][Theorem I.3.17] and [110][Corollary 2.4]. It is also known that ; see [110][page 4] and [144][Lemma 4.1].
As the range of jumps , we have . Also, the integrand of is bounded. Then, we may bound
10.2. LLN for
The law of large numbers for is already interesting. In some sense should be a random walk, but it is impeded by the other particles. How much is the question.
Theorem 10.2.1.
Starting under , we have a.s. and in that
Proof.
Write
Since the fourth moment of is bounded by and is extremal and therefore ergodic to time shifts (Proposition 8.3.5), we have by the ergodic theorem (Birkhoff) that the right side converges to its mean, equal to , both a.s. and in . ∎
We comment the more simple estimate (cf. (10.1.5)) does not lead immediately to a.s. convergence, although convergence would hold.
10.3. General CLT for
We now discuss a central limit theorem for when the jump probability is symmetric, first proved in [131]. It can be extended to a functional CLT with respect to a Brownian motion limit, but that is not done here. We mostly follow the scheme in [131].
Theorem 10.3.1.
Let be symmetric. Starting under , we have
where is a covariance matrix.
The limiting covariance is not explicit, although there are physics conjectures, and some rigorous results for its behavior as a function of [139]. However, it is nondegenerate, except for a particular case.
Proposition 10.3.2.
Except in the case, and is nearest-neighbor (), the covariance is nondegenerate.
We first prove Theorem 10.3.1 and then later Proposition 10.3.2. Let be the local, bounded drift function. For , it will be useful to bound the norm, with respect to , of .
Lemma 10.3.3.
We have belongs to .
Proof.
Since is symmetric, is a linear combination of for in the support of . Consider, as is invariant with respect to the operator , for local functions that
Now, the last quantity is less than where . Hence, belongs to . ∎
Proof of Theorem 10.3.1.
From the martingale decomposition (10.1.3), we have
Since is symmetric, the process is reversible with respect to . We want to apply the Kipnis-Varadhan apparatus to the second term.
By Lemma 10.3.3, we have . Then, by Proposition 9.3.2, there is a square integrable martingale and error such that vanishes in ,
and
| (10.3.1) |
Hence,
At this point, one completes the argument by the martingale central limit Theorem 9.1.1, since has stationary and ergodic increments and, noting (10.3.1), (10.1.5),
Finally, we comment that the covariance satisfies
10.3.1. Proof of Proposition 10.3.2
In the exceptional case, the two martingales and cancel each other, leading to subdiffusive fluctuations for discussed in the next subsection. However, in all other situations, full cancellation does not occur.
The argument relies on the following bound. Recall the Dirichlet form from (10.1.2).
Lemma 10.3.4.
We have for and that
Proof.
From the form of , we need only show the result for where is in the support of . The key point is that, aside from the exceptional case, one can either go around, in terms of exchanges of coordinates, the origin in or hop over it in , without disturbing the tagged particle sitting at the origin.
One can rewrite, by irreducibility of , that
in terms of a sequence from to so that for . Then,
Proof of Proposition 10.3.2.
From the formula for the limiting variance, we need to show there is an such that
In the argument below, in the support of will suffice.
Step 1. Note that is bounded, and consider the resolvent equation
Multiply by and take expectation to obtain . Since by Lemma 10.3.3, we obtain that . Hence, we obtain the inequality
| (10.3.2) |
by Lemma 10.3.4.
Let . Then,
By a similar argument as for (9.3.4) (in the proof of the Kipnis-Varadhan Proposition 9.3.2), we have
An explicit computation gives that the right hand side equals
The first term is . Then, by dropping the second term, we have the lower bound
Step 2. Now, if vanishes as , we have
and so the limiting variance would be which is positive for in the support of .
On the other hand, if , then also does not vanish by (10.3.2). Hence, still the limiting variance is positive. ∎
10.4. Subdiffusive CLT in the exceptional case
Since particles are ordered in when the transition probability is nearest-neighbor, one may suspect that the expected displacement of a particle n time may be less than diffusive when is additionally symmetric. This is indeed the case, when starting under .
Theorem 10.4.1.
Starting under , we have
where .
A functional CLT is also known, but here instead of Brownian motion the limit will be a fractional BM with Hurst parameter [165].
There are a couple of ways to prove this subdiffusive CLT, first proved in [9]. Another proof was given in [183]. We present a method which makes a connection with the current through the bond based on [57].
Define
where is the number of particles which start on the left of and move across the bond up to time , and is the number of particles to the right of crossing to the left.
Then, we have the following relations between and :
| (10.4.1) |
We will show the following result which will imply Theorem 10.4.1. The proof of Theorem 10.4.2 is deferred to Subsection 10.4.3.
Theorem 10.4.2.
Starting under , we have
We now give the proof of Theorem 10.4.1.
Proof of Theorem 10.4.1.
Write, for , that
Since, under , are i.i.d. Bernoulli random variables,
Then, from Theorem 10.4.2, we have
The right-hand side is the same as .
A similar argument works for . ∎
10.4.1. Stirring representation
One way to prove Theorem 10.4.2 is to use a ‘graphical’ representation of symmetric exclusion process (cf. [147][Chapter 8). We present the representation in the nearest-neighbor setting in .
Recall that the generator for symmetric simple exclusion has form (cf. (2.2.2))
In other words, the process evolves by ‘exchanging’ values with neighboring coordinate variables.
Now, consider a collection of Poisson processes with intensity indexed by the bonds . Place a particle at each vertex in . At the event times of the Poisson process with index , the particles at and exchange positions. Let be the location of the particle intially at at time . Marginally, each has the statistics of a nearest-neighbor symmetric random walk, although jointly they are dependent.
We claim that
that is exactly when there is an such that and . Infinitesimally, the possible transitions are exchanges of nearest-neighbor values which is the same as given in the generator . Then,
Define now
which represent the number of stirring particles starting from the left of the origin which end up on the right at time , and vice versa.
Since, in the stirring process, a crossing of the bond in one direction corresponds to a crossing in the other direction, as all sites are occupied by stirring particles, we have is constant in time. But, , and hence for all . When , let be the random locations where , and be the random locations where . Then, we can write the current as
| (10.4.2) |
where ; Note the time dependence on comes from the random locations which depend on . Conditionally on , and the random locations, the variables are independent with respect to . Moreover, conditionally, only the variable , if , is different from which are i.i.d., taking values with mean and variance .
10.4.2. Estimates on
In the following, and refer to the stirring process probability and expectation when is distributed according to . We have the following properties of .
Lemma 10.4.3.
We have
Proof.
Write, by symmetry,
where is a simple random walk starting at the origin and is the positive part.
Now, is a uniformly integrable sequence of random variables which converges in distribution to . Hence,
We also have that the variables and are negatively correlated, which makes intuitive sense since the exchange mechanism is repulsive to an extent.
Lemma 10.4.4.
For , we have that
We refer the reader to Lemma 4.12 [147] or page 365 [9] for proofs of the lemma (in a more general context). See also the discussion in [150], which shows a very general form of negative association, a ‘strong Rayleigh’ condition for symmetric exclusion processes started from product measures.
A consequence of the Lemma 10.4.4 is that .
10.4.3. Proof of Theorem 10.4.2
Consider the following steps.
Step 1. To prove the theorem, noting (10.4.2), we show that the following difference converges in probability, as ,
We make explicit the contribution of the term which might reflect the tagged particle initially at the origin: Write, noting empty sums vanish,
with respect to variables which are i.i.d. mean and variance , conditionally on .
Then, the variance, noting the conditional mean vanishes,
Since is the sum of negatively correlated random variables, we have and also by Lemma 10.4.3. Hence, for , we can apply Chebychev’s inequality to get
Step 2. Now, since is bounded, a.s. Hence, to finish the argument, we will take the limit of
Note, by Lemma 10.4.3, that with high probability for each . Then, there are roughly summands in the sum above, with high probability.
Hence, the last display may be approximated by
where is the characteristic function of conditioned on and the random locations, which has an explicit distribution on values with mean and variance . One now follows the strategy for the usual CLT to obtain the last term converges to the characteristic function of a Normal r.v. with variance as desired. ∎
10.5. Notes
The tagged particle problem, starting from , is somewhat complete at the level of LLN and CLT for symmetric exclusion processes with finite range on . See [147], [148] and references therein for treatments. See also [117] for longer range symmetric exclusion processes. Large and moderate deviations have also been studied [155], [200], [201], [218]. See also [209] and references therein for related work on random walk in random environment.
For asymmetric simple exclusion, a CLT has been shown in diffusive scale when [202], or when and the jumps are nearest-neighbor [128]. Also, when the jumps are mean-zero, finite-range, but asymmetric in , the CLT also holds [212].
When the jumps are asymmetric with a nonzero drift, the variance of the tagged particle at time is shown to be diffusive, in the sense of Laplace transforms for all [197]. However, the associated CLT, when and when and the jumps are not nearest-neighbor, is open.
When starting out of the invariant measure, the LLN behavior has been studied in [180], but less is known for the fluctuations or other limits; see however [118], [47], [49], [50], [65], [211].
See also [114], [120], [121], [198] for discussions and some results for a tagged particle in zero-range models. On other spaces, such as trees, see [45], [93] for recent developments.
We comment that the tagged particle problem is an old one going back at least to Einstein’s 1905 opus. Recently, derivation of Brownian motion rigorously from deterministic hard sphere dynamics with elastic collisions, where initially positions of particles are in thermal equilibrium, has been addressed in [29]. See [90], [7], [206], [210] for reviews and discussions.
Section 11 Equilibrium fluctuations of symmetric exclusion
We motivate the study of ‘equilibrium’ fluctuations of -dimensional symmetric exclusion via another proof of the -scaling limit of the current, Theorem 10.4.2, in the one dimensional nearest-neighbor process. ‘Equilibrium’ fluctuations capture the CLT behavior of the empirical mass density around the constant solution of the hydrodynamic equation, when started from an invariant measure . The limiting equation is a type of linear stochastic heat equation, an SPDE whose solution is an infinite dimensional ‘Ornstein-Uhlenbeck’ process.
We will focus our discussion in the symmetric simple exclusion context when with respect to the standard basis .
11.1. Another derivation of the current fluctuations limit
Consider the nearest-neighbor, symmetric exclusion setting discussed in Subsection 10.4. Recall the current through the bond in is the number of particles crossing the bond from left to right minus the number crossing from right to left up to time .
A moment’s thought gives, reflecting the boundary values at , that
the difference being equal to , or . Then, formally,
However, the display does not make sense as typically there are an infinite number of particles in the system.
We may truncate however in the following way: Let
Write, in terms of a scaling parameter , that
At the same time, we have
Following our custom, since space is scaled by , we will speed up time by , and define the mass ‘fluctuation field’ with respect to by its action on :
| (11.1.1) |
here written for .
In , scaling the current by the fourth root of the time scaling, we have for all that
| (11.1.2) |
The idea now is that the first two terms on the right should be given in terms of a limit fluctuation field , which we must define, while the last term on the right hand side should vanish as and .
One can adjust the construction of the process (cf. Remark 10.1.5), to include counts and , keeping track of the numbers of particles crossing from left to right and vice versa. Then, and
are martingales. Since jumps are not simultaneous in the process, are orthogonal martingales.
Exercise 11.1.1.
Show by decomposing on the possible crossing times of bonds and , which are stopping times occurring a.s. at distinct times, and the martingale property. That is, write and a similar formula for where are the jumps on these bonds up to time . Then, we have
equals zero if since jumps are not simultaneous. But, if , then since is a martingale, the display also vanishes.
In the following, recall as before that stands for the probability and expectation under , and and for the process measure and expectation when starting in .
Lemma 11.1.2.
Starting in an invariant measure, , we have
11.2. Infinite dimensional Ornstein-Uhlenbeck process limit
We consider now the general setting. Consider the space , consisting of compactly supported, real functions on . Its dual is the space of distributions, that is the continuous linear functionals acting on .
Recall the definition of in (11.1.1). To capture the evolution of for a fixed , write, with respect to the finite-range symmetric exclusion generator (cf. (2.2.2)),
where, after the usual twice summation-by-parts under symmetry,
and is a martingale such that
is a martingale. The last integral in the display, the quadratic variation of , can be evaluated as
Here, similar to the hydrodynamics calculations (cf. Subsection 2.4), and are the scaled discrete gradient and Laplacian.
Since , we may subtract and write
| (11.2.1) |
Now, as before with respect to the hydrodynamics, we have two steps:
Step 1: Show tightness of in an appropriate space, and continuity of limit trajectories under limit points.
Step 2: Identify the limit points in terms of a unique ‘infinite dimensional Ornstein-Uhlenbeck’ process.
Putting it together, we will arrive at the following informal description: converges to which solves, in the nearest-neighbor symmetric setting,
| (11.2.2) |
more explained in the following subsections.
Spaces
A natural space will be , the distribution valued right-continuous with left limits trajectories on , with the strong dual topology.
We remark there are other natural spaces one could consider. For instance, could act on the Schwartz space of rapidly decreasing functions , in which case would be a member of , the space of ‘tempered distributions.’ One could work also with the Hermite basis based spaces containing . Then, would be a member of . In another direction, we could have also specified the lattice as , instead of . See Chapter 11 [130] for instance where spaces with respect to are used; on , the Hermite spaces would use definitions in [177][p. 142].
11.2.1. A precise statement of (11.2.2)
Let be the probability measure on governing , when the underlying nearest-neighbor, symmetric exclusion process starts from invariant measure .
Theorem 11.2.1.
We have converges to , concentrated on , governing a Gaussian Markov random field with mean and covariance,
for all , .
The process governed by is sometimes called a ‘generalized Ornstein-Uhlenbeck process’, as discussed in Holley-Stroock [112]. The scalings and the limit in Theorem 11.2.1, where space is scaled by , time by and the empirical measure by , are sometimes referred to as a ‘Edwards-Wilkinson’ fluctuation limit.
Holley-Stroock martingale problem
The existence/uniqueness of in Theorem 11.2.1 follows from a ‘martingale problem’ characterization.
Let be the nonnegative self-adjoint operator defined on domain and let be the associated heat semigroup. Define to be the linear gradient operator. Let also be the sigma-field in generated by a process for and .
Theorem 11.2.2.
Suppose is a probability measure governing , concentrating on , and for each ,
and
are , -martingales. Then, for all , and subsets , a.s.,
| (11.2.3) | |||
Therefore, the finite-dimensional distributions of , and hence the measure on , are determined by its restriction to .
Sketch of proof of Theorem 11.2.1
Assuming Steps 1,2, we outline the proof of Theorem 11.2.1. In our context, the restriction to is already known: Under the invariant measure , we have converges to a Gaussian field with mean zero and covariance
| (11.2.4) |
Exercise 11.2.3.
Show that the joint distribution of is Gaussian with covariance given via (11.2.4). One can use the Cramér-Wold device.
Also, from the martingale property in Theorem 11.2.2, the formula (11.2.3), and that , we can argue
| (11.2.5) | |||||
This covariance is exactly what is given in Theorem 11.2.1.
Now, since is a continuous martingale with quadratic variation , by Levy’s characterization, we have that is distributed as a Brownian motion. Hence, we have
| (11.2.6) |
where is the infinite dimensional Brownian motion (Gaussian Markov random field) with covariance
In this way, (11.2.6) gives a meaning to the integral form of (11.2.2).
Exercise 11.2.4.
Use polarization with the martingales and to show the formula in the last display.
11.2.2. Application: Proof of Theorem 10.4.2
We return to our motivating example with respect to current fluctuations. By (11.1.2) and Lemma 11.1.2, since does not depend on , we have uniformly in that
is a Cauchy sequence in . Note, by stationarity, that
Then, by Theorem 11.2.1, approximating by smooth compactly supported functions, we have for fixed as that
Since is Cauchy in , we denote its limit by , whose distribution is a mean-zero Gaussian.
11.3. Proof of Steps 1,2: Tightness and identification
We first restate Steps 1, 2 given informally before (11.2.2) in the setting of nearest-neighbor symmetric exclusion.
Step 1: Show tightness of governing , members in , in the uniform topology of . Hence, trajectories under limit points will be continuous.
Step 2: Identify the limit points in terms of the unique ‘infinite dimensional Ornstein-Uhlenbeck’ process given in Theorem 11.2.1.
In the following, recall that is nearest-neighbor, symmetric and translation-invariant.
11.3.1. Moments of
Recall the martingales and specified in (11.2.1). Before going to the proofs of Steps 1,2, we record some moment estimates of for .
By taking expectation of the quadratic variation, we obtain
| (11.3.1) |
To bound the fourth moment, we will take the following approach. Other arguments, via Burkholder-Gundy-Davis inequalities may also be used.
Lemma 11.3.1.
For all local functions , , and , we have
is a martingale.
Proof.
The lemma is a type of ‘Girsanov’ formula. See [66][around p. 175] which shows is a local martingale. Since we are dealing with the exclusion process, where occupation numbers are bounded, it is integrable and so a martingale. ∎
Lemma 11.3.2.
For , and , we have
Proof.
Let be the martingale with
By explicit calculation,
Morevoer,
| (11.3.2) |
Exercise 11.3.3.
Make the computation in the proof of the above lemma to match fourth powers of to obtain the lemma statement.
11.3.2. Proof of Step 2, given Step 1
We will apply results in the limit theory of martingales in the context of symmetric simple exclusion. References in this vein include [66], [116], [216], among others.
Suppose is a limit point found with respect to the uniform topology of , necessarily supported on continuous -valued trajectories. We will identify the subsequence by itself to streamline notation.
We now show that and are martingales, and thereby identify via Theorem 11.2.2.
Consider the formula for in (11.2.1). Given Step 1 and that the function is continuous in the uniform topology, the term would converge in distribution to a limit
Observe, by (11.3.1), that . Therefore, we conclude is uniformly integrable for each . Since the martingales , we conclude the limit is an martingale by Theorem IX.1.12 in [116].
Also, as
| (11.3.3) | ||||
we have converges to in probability.
By the continuous mapping theorem, . Moreover, by the convergence in probability implied by (11.3.3), we observe . To see that the limit is an martingale, we observe is uniformly integrable by the fourth moment estimate Lemma 11.3.2. Hence, as desired is a martingale, again by Theorem IX.1.12 in [116].
11.3.3. Proof of Step 1
To show tightness of with respect to the uniform topology, as is the ‘inductive limit’ of nuclear Frechét spaces, it is sufficient to verify tightness of for each ; see [78].
We remark, as in the study of hydrodynamics, tightness with respect to uniform topology implies that in the Skorohod topology. We have therefore the following tightness criterion. Denote the uniform modulus of continuity by
Lemma 11.3.4.
A family of probability measures on is tight if for each ,
- (a)
and,
- (b)
for all , .
Proof.
Prelimit, we have
| (11.3.4) | ||||
We need to show each of these terms satisfies the criteria.
The term . The mean square, starting from the invariant measure , is bounded in the limit by , and hence item (a) holds. Part (b) holds trivially.
The integral term. Part (a) holds trivially. For item (b), by Chebychev and Schwarz inequalities, we may bound
Starting from , the last quantity is bounded by , which vanishes as .
The term . Part (a) holds trivially. For item (b), We employ the standard ‘three argument, dividing the interval into subintervals (a similar scheme was used in Subsection 3.3). Noting when belong to a single subinterval or are in adjacent ones, we may bound
where the max is over the subintervals and is the left endpoint of the th one. Then, by stationarity, Chebychev and Doob’s inequality and Lemma 11.3.2,
The right-hand side of the display vanishes as and . ∎
11.4. Notes
The proof given here for the current fluctuations Theorem 10.4.2 follows [118], and the idea of truncation goes back to [183].
We comment, with respect to equilibrium fluctuations, in the context of symmetric exclusion models, as for the hydrodynamic limit, there is no ‘replacement’ estimate needed; see also the treatment in [130]. The theory of generalized Ornstein Uhlenbeck processes originates in [112].
Recently, a proof of fluctuations starting from ‘flat’ states via ‘discrete regularity structures’ has been given in symmetric exclusion [113].
‘Nonequilibrium’ fluctuations have also been shown in symmetric exclusion [175]. In asymmetric models, fluctuations starting from an invariant measure are also known in [42].
11.4.1. KPZ etc.
In , for the nearest-neighbor asymmetric process, including TASEP, there are many works describing the behaviors of currents, ‘height’ functions, and fluctuation fields starting from an invariant measure and other initial conditions. Instead of a generalized OU process, under time scaling , space scaling , and say scaling of the ‘height’ functions, the limits are different, those with respect to a ‘KPZ’ (Kardar-Parisi-Zhang) fixed point (cf. [158]).
On the other hand, if the asymmetry is weak in the sense the drift of the transition probability is , under diffusive time/space scaling, the limit of certain translated fluctuation fields solves a ‘stochastic Burgers’ or sometimes called ‘KPZ-Burgers’ equation.
These limits connect with integrable probability, last passage percolation and the ‘Directed Landscape’, polymers, singular KPZ SPDE, and random matrices. See [4], [11], [25], [51], [55], [56], [97], [98], [103], [99], [104], [105], [107], [115], [125], [153], [171], [207], [219], [220] for reviews and discussion, among other references.
In , the work [39] on scaled continuum stochastic Burgers equation shows Gaussian fluctuations. In , Gaussian behaviors have been shown with respect to scaled ‘subcritical’ continuum KPZ equations [44], [40], [102]. In , Gaussian behaviors have been seen starting from asymmetric exclusion [141]. As mentioned in [39], in , it is open to see such behavior starting from asymmetric exclusion, or other particle systems.
Section 12 Boltzmann-Gibbs principle for symmetric zero-range processes
We discuss the equilibrium fluctuations of -dimensional symmetric finite-range zero-range processes, that is the scaling limit of the empirical mass fluctuation field, starting from an invariant measure . Unlike for symmetric exclusion models, a ‘replacement’ at the fluctuation level must be made in order to identify the limit field as a generalized OU process. This replacement, sometimes called the ‘Boltzmann-Gibbs’ principle, is of its own interest. The proof we give makes use of a ‘spectral gap’ estimate.
12.1. Statement of equilibrium fluctuations
Recall, from Sections 4 and 7, that the rate function and the jump probability define the zero-range process with generator
To reduce notation, we will assume that is symmetric, translation-invariant, and nearest-neighbor: for the standard basis . Also, we remind that satisfies , for . We will also impose the following.
- (Lip)
for all .
- (M)
There exists and such that for all .
The first assumption is something we have already seen in the construction of the process on in Section 7, which includes the independent particle process when . The second assures a uniform ‘spectral gap’ that will be discussed later in Subsection 12.4.
We will also fix an invariant measure , a product of ‘Poisson’ like marginals, and the process will be assumed to begin under this extremal invariant measure. We will assume the process is an process on . Recall Section 8 for specifications and details. There is no difficulty in assuming if preferred that the process is on the torus where there are no construction issues. As before, denotes the probability and expectation under , and and the process measure and expectation when starting in .
Recall, as in Section 11, that for fixed smooth with compact support stands for the fluctuation field,
To derive the limit field, write as before
where, after the usual summation-by-parts,
Here, , and is a martingale such that
is a martingale. The last integral in the display, equal to , can be evaluated as
These calculations are analogous to those for hydrodynamics in Sections 4, 5.
Following the method described for exclusion processes in Section 11, we have two steps:
Step 1: Show tightness of in an appropriate space, and continuity of limit trajectories under limit points.
Step 2: Identify the limit points in terms of a unique ‘infinite dimensional Ornstein-Uhlenbeck’ process
As before, the space of trajectories is , where is the dual space of distributions with respect to . Step 1, tightness, is accomplished as before, and it is left to the reader to verify the proof in the zero-range setting.
The more interesting part, for zero-range processes, is Step 2 where instead of the occupation variable we have a function of it, namely in both martingales above. If the normalization were instead of , replacing the nonlinear function with a function of the mass empirical density is the ‘standard’ hydrodynamic replacement. Here, we have to work a little harder. However, starting in equilibrium helps.
The following ‘Boltzmann-Gibbs’ estimate allows the replacement in .
Theorem 12.1.1.
For smooth, compactly supported , we have
The replacement, however, in the square martingale will be a consequence of the following limit.
| (12.1.1) |
Hence, we arrive at the following result, specializing to the nearest-neighbor setting, when with respect to the standard basis .
Theorem 12.1.2.
We have that converges to where
| (12.1.2) |
The characterization of in terms of the generalized OU martingale problem of Holley and Stroock is as before with symmetric exclusion in Section 11. Though, the operators and differ in prefactor constants.
Analogous to before with symmetric exclusion,
where . Also, the covariance of and satisfies the formula (11.2.5) where is the semigroup associated to .
One way to look at (12.1.2) is to relate it to the hydrodynamic equation:
with initial condition . The quantity can be seen as an ‘error’ via a ‘linearization’ of the hydrodynamic equation about the equilibrium density (the constant solution when starting from ), where . See [206] for more physical intuition behind this intepretation.
We remark in passing, in the case , the setting of independent particles, the replacement Theorem 12.1.1 is not needed as, analogous to symmetric exclusion, the martingale is already ‘closed’ with respect to , that is fully expressed in terms of and .
12.2. Kipnis-Varadhan estimate
We begin with the following ‘Kipnis-Varadhan’ [131] non-asymptotic bound of independent interest, helpful to establish the Boltzmann-Gibbs principle. Recall the notions of and norms from Section 9. These extend to the zero-range context, where the local functions used to define spaces and are local functions.
Lemma 12.2.1.
For all local functions, we have for the symmetric zero-range process
Note that it may be that diverges for a given .
Proof.
Write the resolvent equation, for ,
Multiplying by and integrating, we have
where . Note that , and so . Now, consider the martingale
with quadratic variation .
Write
Since , from squaring the left hand side of the above display, using , by stationarity, we have
Choosing , we see that the left hand side is bounded by
12.3. Derivation of Boltzmann-Gibbs estimate
We need to show the variance of
vanishes in the limit. It will be enough to bound the variance of
| (12.3.1) |
for , from which Theorem 12.1.1 can be deduced.
Write (12.3.1) as
| (12.3.2) | ||||
Here, is a block of width centered at , and is another scaling parameter. Related to the study of the hydrodynamic limit, the idea is that considers the approximation of by its conditional expectation average with respect to a local density. The term gives the error with respect to the leading order terms in this conditional expectation.
The strategy will be to bound the norm of , and to use Schwarz inequality and Taylor expansions with , a sort of ‘equivalence of ensembles’ estimate. These are done in the next two subsections. In Subsection 12.3.3, we assemble these bounds to prove Theorem 12.1.1.
12.3.1. Bound on
To bound the variance of , we will bound the norm of its integrand, denoted . Since time has been sped up by , the generator of is . For local functions , we would like to demonstrate
where the Dirichlet form and is a bound of the maximum ratio . Then, the norm of is bounded by , which we show vanishes as .
To this end, write
given depends only on variables indexed by , where is the shift operator, is the restriction, and . We comment, as is a product measure, that .
We may further evaluate the right-hand side as
where is the ‘canonical’ measure, that is conditioned on there being particles in , and .
Recall (cf. (4.5.1) from Section 4). Underlying the above localization, that is for adding and subtracting the conditional expectation , is that we can now solve a certain ‘Poisson’ equation, no matter the particle number , since for all .
Indeed, note that is the unique invariant measure for the irreducible zero-range dynamics localized to with generator
and is mean-zero with respect to , which assigns probabilities to configurations in . Moreover, one may verify that is reversible with respect to . Also, the Dirichlet form with respect to may be computed,
So, in particular, is a nonnegative, reversible operator.
Exercise 12.3.1.
Verify the reversibility of and the form of the Dirichlet form .
Now, thinking of as a vector in , since it is in null space of the transpose as it is a stationary measure, and so that is orthogonal to , we conclude that belongs to the range of . See Subsection 9.1 for related decompositions.
Hence, we may solve for some function . In particular, as is nonnegative, and reversible with respect to ,
At this point, we state a spectral gap inequality to be discussed later. Namely, for mean-zero functions, we have that
where a bound of the inverse of the ‘spectral gap’ is independent of when (Lip) and (M) hold for the rate [142].
Hence, by the relation for , we have
Note, by unraveling the measures, that
Also, by estimating the overcount of terms over bonds , we have
where the full Dirichlet form is evaluated
In addition, .
Then, coming back to (12.3), summing over , we have
after optimizing on . Here, we used translation-invariance to yield and the constant may have changed line to line.
Hence, the norm of is bounded
and therefore the variance of by Lemma 12.2.1 is bounded
| (12.3.3) |
as for each fixed.
12.3.2. Bound on
Since is smooth, we may replace with in . By invariance of , adding and subtracting , and Schwarz inequality,
Then, by Schwarz inequality, using invariance and translation-invariance of , we have
| (12.3.4) | ||||
The factor arises since the summands over are independent only when they are separated by distance .
We now estimate the mean-square of
Outside truncation
We first truncate the number of particles in . Let . Recall . Then, by Schwarz inequality,
| (12.3.5) |
Here, we used Chebychev to bound .
Inserting into (12.3.4), we see that the multiplying factor is compensated, and the cost of truncation is small, uniformly over , for large .
Inside truncation
We now use a Taylor expansion and a local central limit theorem. For , let . Given the truncation, we have that is bounded away from and uniformly in .
Since is a product measure with marginals in the form given in Subsection 4.3, the conditional expectation does not depend on the density and can be changed to another, say ; see near (5.4.4) for a similar discussion with respect to the -block hydrodynamics lemma. Recall the calculation in Exercise 4.3.1,
where is the configuration with a single particle at . Then,
| (12.3.6) | ||||
Recall . By the local central limit Theorem VII.13 in [167] (applied with there), when and is bounded uniformly away from and , we have for and large ,
For , we see , whereas for , we have , in which case for large .
12.3.3. Proof of Theorem 12.1.1
12.4. Spectral gap
One can define the ‘spectral gap’ for a reversible process generator in terms of a Poincare inequality:
Then, the gap would be the reciprocal of the smallest for which the inequality is true for all .
Equivalently, the gap is the difference between the two largest eigenvalues of . Since is the largest eigenvalue, the gap would be the negative of the second largest eigenvalue, . Since is a matrix, one can in principle compute it exactly, but this is usually hard to do when the state space gets large.
The spectral gap has connections with the mixing time of the process. Recall that
The larger , the faster the approach of to its mean .
In zero-range processes, since the jump times are controlled by the rate , one expects the orders of the mixing time to depend on . This is in fact the case. When , the independent case, the mixing is rapid, and corresponds to the mixing time of a single random walk on which is . For instance, consider , then a nearest-neighbor walk, in moving from one end of the interval to the other, has to take steps which usually takes units of time.
If is not too different from the independent case, one expects similar behavior, not dependent on the number of particles, and this is the result quoted and used above, valid under our assumptions (Lip) and (M) [142].
If is sublinear, that is for , then it has been shown the spectral gap depends on the number of particles , namely the gap is where [163]. For instance, if , movement becomes difficult, and the gap vanishes!
If , then it can be seen that the gap is [162]. When , this is the case where there is a connection between the zero-range process and simple exclusion: The number of spaces between particles in simple exclusion correspond to the number of particles in the zero-range process.
It seems to be an open problem to characterize the gap when is bounded in general. One presumably expects the behavior as in the last case.
12.5. Notes
The Boltzmann-Gibbs principle, under a stationary measure, was named and proved first by Brox and Rost [33] which did not use a spectral gap assumption. See also [130] for a presentation along this line. The notion of using the spectral gap estimate to prove the Boltzmann-Gibbs principle seems to be newer; see also [97], [98].
While equilibrium fluctuations of particle systems with finite-range symmetric interactions is relatively well understood (cf. [130], [133], [206]), a main open problem is to understand ‘nonequilibrium’ fluctuations of the density field in a general class of particle systems. A main difficulty is to perform the needed ‘replacement’ estimates in the fluctuation scales. Results in dimension are known for finite-range symmetric models; see [43], [74]. Recently, these results have been generalized to , for a class of particle systems [123]; see also [53], [86], [119], and discussions therein.
Bibliography
- [1] Adams, D., Jay, J., Balázs (2026) Second class particle behaviour in blocking ASEP. ALEA Lat. Am. J. Probab. Stat. 23 367–416.
- [2] Aggrawal, A., Corwin, I., Ghosal, P. (2023) The ASEP speed process. Adv. in Math. 442 109004
- [3] Amir, G., Bahadoran, C., Busani, O., Saada, E. (2025) Invariant measures for multilane exclusion process. Ann. Inst. H., Poincaré (B) Probab. Stat. 61 2184–2234.
- [4] Amir, G., Corwin, I., Quastel, J. (2011) Probability distribution of the free energy of the continuum directed random polymer in 1+1 dimensions. Comm. Pure Appl. Math. 64 466–537.
- [5] Andjel, E.D. (1982) Invariant measures for the zero range process. Ann. Probab. 10 525-547.
- [6] Andjel, E. and Vares, M.E. (1987) Hydrodynamic equations for attractive particle systems on . J. Stat. Phys. 47 265-288.
- [7] Ardourel, V. (2022) Brownian motion from a deterministic system of particles Synthese 200 29
- [8] Armendáriz, I., Grosskinsky, S. and Loulakis, M. (2017) Metastability in a condensing zero-range process in the thermodynamic limit. Probab. Theory and Rel. Fields 169 105–175.
- [9] Arratia, R. (1983) The motion of a tagged particle in the simple symmetric exclusion system on . Ann. Probab. 11 362-373.
- [10] Bahadoran, C. Guiol, H., Ravishankar, K., and Saada, E. (2006) Euler hydrodynamics of one-dimensional attractive particle systems. Ann. Probab. 34 1339-1369.
- [11] Baik, J. (2023) KPZ limit theorems. In International Congress of Mathematicians 4190–4211. European Mathematical Society-EMS-Publishing House GmbH.
- [12] Baik J., Ferrari P.L., Péché S. (2014) Convergence of the two-point function of the stationary TASEP. Singular Phenomena and Scaling in Mathematical Models Springer, 91–110.
- [13] Balász, M., Seppääinen, T. (2010) Order of current variance and diffusivity in the asymmetric simple exclusion process. Ann. of Math. 171 1237–1265.
- [14] Balázs, M., Rassoul-Agha, F., Seppäläinen, T., Sethuraman, S. (2007) Existence of the zero range process and a deposition model with superlinear growth rates. Ann. Probab. 35 1201–1249T
- [15] Bannai, K., Sasada, M. (2026) Varadhan?s decomposition of shift-invariant closed L 2-forms for large scale interacting systems on Euclidean lattices. Elec. J. Probab. 31 1–49.
- [16] Basu, R., Sarkar, S., Sly, A. (2017) Invariant measures for TASEP with a slow bond. arXiv:1704.07799
- [17] Beltrán, J., Landim, C. (2015) A martingale approach to metastability. Probab. Theory Rel. Fields 161 267–307.
- [18] Bernardin, C. (2004) Fluctuations in the occupation time of a site in the asymmetric simple exclusion process. Ann. Probab. 32 855–879.
- [19] Bernardin, C., Funaki, T., Sethuraman, S. (2021) Derivation of coupled KPZ-Burgers equation from multi-species zero-range processes. Ann. Appl. Probab. 31 1966–2017.
- [20] Bernardin, C., Gonçalves, P. and Jiménez-Oviedo, B. (2021) A microscopic model for a one parameter class of fractional Laplacians with Dirichlet boundary conditions. Archive for Rational Mechanics and Analysis 239 1–48.
- [21] Bernardin, C., Gonçalves, P., Sethuraman, S. (2014) Equilibrium fluctuations of additive functionals of zero-range models. (2014) In Conference Proceedings, ‘Particle systems and PDE’ Dec. 5-7, 2012, U. Minho, Springer proceedings in Mathematics and Statistics, 75, 143–169.
- [22] Bernardin, C., Gonçalves, P., Sethuraman, S. (2016) Occupation times of long-range exclusion and connections KPZ class exponents. Prob. Theory Rel. Fields 166, 365–428.
- [23] Bertini, L., De Sole, A., Gabrielli, D., Jona-Lasinio, G., Landim, C. (2015) Macroscopic fluctuation theory. Reviews of Modern Physics 87 593–636.
- [24] Bertini, L., Gabrielli, D., Landim, C. (2023) Concurrent Donsker?Varadhan and hydrodynamical large deviations. Ann. Probab. 51, 1298–1341.
- [25] Bertini, L.; Giacomin, G. (1997) Stochastic Burgers and KPZ equations from particle systems. Comm. Math. Phys. 183 571–607.
- [26] Bhattacharya, R.N. (1982) On the functional central limit theorem and the law of the iterated logarithm for Markov processes. Z. Wahrsch. Verw. Gebiete 60 185–201.
- [27] Billingsley, P. (1968) Convergence of Probability Measures. John Wiley & Sons, New York.
- [28] Bodineau, T., Bauerschmidt,, R., Dagallier. B. (2025) A criterion on the free energy for log-Sobolev inequalities in mean-field particle systems. HAL-open science, hal-05014812
- [29] Bodineau,T., Gallagher, I., Saint-Raymond, L. (2016)The Brownian motion as the limit of a deterministic system of hard-spheres. Inventiones Mathematicae 203 493–553.
- [30] Bodineau, T. and Derrida, B., (2006) Current Large Deviations for Asymmetric Exclusion Processes with Open Boundaries. J. Stat. Phys. 123 277–300.
- [31] Bramson, M., Liggett, T.M. (2005) Exclusion processes in higher dimensions: Stationary measures and convergence. Ann. Probab. 33 2255–2313.
- [32] Bramson, M., Liggett, T.M. and Mountford, T. (2002) Characterization of stationary measures for one-dimensional exclusion processes. Ann. Probab. 30 1539–1575.
- [33] Brox, T., Rost, H. (1984) Equilibrium fluctuations of stochastic particle systems: The role of conserved quantities. Ann. Probab. 12 742-759.
- [34] Caputo, P. (2004) Spectral gap inequalities in product spaces with conservation laws. in Advanced Studies in Pure Mathematics, H. Osada and Funaki T eds., Japan.
- [35] Caputo, P., Liggett, T. M. and Richthammer, T. (2010). Proof of Aldous’ spectral gap conjecture. J. Amer. Math. Soc. 23 831–851.
- [36] Caputo, P., Quattropani, M. and Sau, F. (2025) A universal cutoff phenomenon for mean-field exchange models. arXiv:2506.12816.
- [37] Cancrini, N., Martinelli, F., Roberto, C. (2002) The logarithmic Sobolev constant of Kawasaki dynamics under a mixing condition revisited. Ann. IHP Prob. Stat. 38 385-436.
- [38] Cannizzaro, G., Gonçalves, P., Misturini, R., Occelli, A. (2025) From abc to kpz. Probab. Theory Rel. Fields 191 361–420.
- [39] Cannizzaro, G., Gubinelli, M., Toninelli, F. (2024) Gaussian Fluctuations for the Stochastic Burgers Equation in Dimension. Commun. Math. Phys. 405 89.
- [40] Caravenna, F., Sun, R., Zygouras, N. (2020) The two-dimensional KPZ equation in the entire subcritical regime. Ann. Probab. 48 1086–1127.
- [41] Carinci, G., Franceschini, C., Groenevelt, W. (2021) -orthogonal dualities for asymmetric particle systems. Elec. J. Probab. 26 1–38.
- [42] Chang, C.C., Landim, C., Olla, S. (2001) Equilibrium fluctuations of asymmetric simple exclusion processes in dimension . Probab. Theory Related Fields 119 381–409.
- [43] Chang, C.C.; Yau, H.T. (1992) Fluctuations of the one dimensional Ginzburg-Landau models in nonequilibrium. Commun. Math. Phys. 145 209-234.
- [44] Chatterjee, S., Dunlap, A. (2020) Constructing a solution of the (2 + 1)-dimensional KPZ equation. Ann. Probab. 48 1014–1055.
- [45] Chen, D., Chen, P., Gantert, N. and Schmid, D. (2019) Limit theorems for the tagged particle in exclusion processes on regular trees. Elec. Commun. Probab. 24 1–10.
- [46] Chen, J.P., Sau, F. (2021) Higher-order Hydrodynamics and Equilibrium Fluctuations of Interacting Particle Systems. Markov Proc. Rel. Fields. 27 339–380.
- [47] Conroy, M., Gonzalez Casanova, A., Sethuraman, S. (2025) Point process convergence of extremes in K-symmetric exclusion. arXiv:2506.12632
- [48] Conroy, M., Rastegar, R., Sethuraman, S. (2026) Remarks on exchangeability of invariant measures for symmetric exclusion. preprint. Available at https://sundersethuraman.github.io/papers/remarks_exclusion.pdf
- [49] Conroy, M., Sethuraman, S. (2022) Gumbel laws in the symmetric exclusion process. Commun. Math. Phys. 402 723–764.
- [50] Conroy, M., Sethuraman, S. (2025) Poisson statistics, vanishing correlations, and extremal particle limits for symmetric exclusion in . Elec. J. Probab. 30 1–38.
- [51] Corwin, I. (2016) Kardar-parisi-zhang universality. Notices of the AMS 63 230–239.
- [52] Covert, P., Rezakhanlou, F. (1997) Hydrodynamic limit for particle systems with nonconstant speed parameter. J. Stat. Phys. 88, 383–426.
- [53] Dagalier, B., Landim, C. (2025) Critical dynamical fluctuations in reaction-diffusion processes. arXiv:2505.20008
- [54] Dai Pra, P., Posta, G. (2005) Logarithmic Sobolev inequality for zero-range dynamics. Ann. Probab. 33 2355–2401.
- [55] Das, S. and Zhu, W. (2022) Upper-tail large deviation principle for the ASEP. Elec. J. Probab. 27 1–34.
- [56] Dauvergne, D., Ortmann, J. and Virág, B. (2022) The directed landscape. Acta Mathematica 229 201–285.
- [57] De Masi, A. and Ferrari, P.A. (2002) Flux fluctuations in the one dimensional nearest neighbors symmetric simple exclusion process. J. Stat. Phys. 107 677-683
- [58] De Masi, A.; Presutti, E. (1991) Mathematical Methods for Hydrodynamic Limits. Lecture Notes in Mathematics, 1501, Springer-Verlag, Berlin.
- [59] Deng, Y., Zaher, M., Ma, X. (2025) Hilbert’s sixth problem: derivation of fluid equations via Boltzmann’s kinetic theory. arXiv:2503.01800
- [60] Diaconis, P., Saloff-Coste, L. (1996) Logarithmic Sobolev inequalities for finite Markov chains. Ann. Appl. Probab. 6 695–750.
- [61] Dizdar, D., Menz, G., Otto, F., Wu, T. (2024) A quantitative hydrodynamic limit of the Kawasaki dynamics. Elec. J. Probab. 29, #192, 1–57.
- [62] Dobrushin, R., and Siegmund-Schultze, R. (1982) The hydrodynamic limit for systems of particles with independent evolution. Math. Nachr. 105 199-224.
- [63] Durrett, R. (2005) Probability: Theory and Examples Third Edition. Duxbury Advanced Series, Brooks/Cole, Belmont, CA.
- [64] El Kettani, P, Funaki, T., Hilhorst, D., Park, H.J., Sethuraman, S. (2022) Mean curvature interface limit from Glauber + Zero-range interacting particles. Commun. Math. Phys. 394, 1173-1223.
- [65] Erhad, D., Franco, T., Xu, T. (2024) Nonequilibrium joint fluctuations for current and occupation time in the symmetric exclusion process. Electron. J. Probab. 29 1–53.
- [66] Ethier, S. N., Kurtz, T. G. (1986). Markov processes: Characterization and convergence. New York: Wiley.
- [67] Evans, L. C. (2010). Partial differential equations (2nd ed.). American Mathematical Society.
- [68] Evans, M.R., Hanney, T. (2005) Nonequilibrium statistical mechanics of the zero-range process and related models. J. Phys. A: Math. Gen. 38 195-240.
- [69] Faggionato, A. (2025) An ergodic theorem with weights and applications to random measures, homogenization and hydrodynamics. Stoch. Proc. Appl. 180 104522.
- [70] Faggionato, A., Jara, M., Landim, C. (2009) Hydrodynamic limit of one dimensional subdiffusive exclusion processes with random conductances. Probab. Theory Related Fields 144 633-667.
- [71] Fathi, M., Simon, M. (2016). The gradient flow approach to hydrodynamic limits for the simple exclusion process. In: Gonçalves, P., Soares, A. (eds) From Particle Systems to Partial Differential Equations III. Springer Proceedings in Mathematics & Statistics, 162, Springer, Cham.
- [72] Ferrari, P.A. and Fontes, L.R.G. (1994) Shock fluctuations in the asymmetric simple exclusion process. Probab. Theory Rel. Fields 99 305–319.
- [73] Ferrari, P.A., Kipnis, C., Saada, E. (1991). Microscopic structure of travelling waves in the asymmetric simple exclusion process. Ann. Probab. 19 226–244.
- [74] Ferrari, P.A., Presutti, E., Vares, M.E. (1988) Nonequilibrium fluctuations for a zero-range process. Stoch. Proc. Appl. 26 31-45.
- [75] Ferrari, P.L., Nejjar, P. (2024) The second class particle process at shocks. Stoch. Proc. Appl. 170 104298.
- [76] Floreani, S., Redig, F., Sau, F. (2021) Hydrodynamics for the partial exclusion process in random environment. Stoch. Proc. Appl. 142, 124–158.
- [77] Fonseca-Mora, C.A. (2019) Tightness and weak convergence of probabilities on the Skorokhod space on the dual of a nuclear space and applications. arXiv:1806:10231
- [78] Fouque, J.P. (1984) La convergence en loi pour les processus á valeurs dans un espace nuclèaire. Ann. IHP Prob. Stat. 20, 225–245
- [79] Franco, T., Gonçalves, P., Neumann, A. (2015) Phase transition of a heat equation with Robin’s boundary conditions and exclusion process. Trans. Amer. Math. Soc. 367 6131–6158.
- [80] Fritz, J. (2001) An Introduction to the Theory of Hydrodynamic Limits Lectures in Mathematical Sciences 18 The University of Tokyo, Tokyo.
- [81] Fritz, J. and Nagy, K. (2006) On uniqueness of the Euler limit of one-component lattice gas models. ALEA Lat. Am. J. Probab. Stat. 1 367-392.
- [82] Fu, B., Funaki, T., Sethuraman, S., Venkataramani, S.C. (2025) Coupled KPZ equations and their decoupleability. arXiv:2508.04637
- [83] Funaki, T. (2018) Hydrodynamic Limit for Exclusion Processes. Commun. Math. Stat. 6, 417–480.
- [84] Funaki, T. (2026) Interface motion from Glauber-Kawasaki dynamics of non-gradient type. Probab. Theory Rel. Fields https://doi.org/10.1007/s00440-026-01473-1
- [85] Funaki, T., Handa, K, Uchiyama, K. (1991) Hydrodynamic limit of one-dimensional exclusion processes with speed change Ann. Probab. 19 245–265.
- [86] Funaki, T., Landim, C., Sethuraman, S. (2026) Linear fluctuation of interfaces in Glauber-Kawasaki dynamics. To appear Commun. Math. Phys.
- [87] Funaki, T., van Meurs, P., Sethuraman, S., Tsunoda, K. (2023) Motion by mean curvature from Glauber-Kawasaki dynamics with speed change. J. Stat. Phys. 190 45.
- [88] Funaki, T., van Meurs, P., Sethuraman, S., Tsunoda, K. (2023) Constant-speed interface flow from unbalanced Glauber-Kawasaki dynamics. Ensaios Math. 38 223–248.
- [89] Funaki, T., Tsunoda, K. (2019) Motion by mean curvature from Glauber-Kawasaki dynamics. J. Stat. Phys. 177 183–208.
- [90] Gallagher, I. (2019) From Newton to Navier-Stokes, or how to connect fluid mechanics equations from microscopic to macroscopic scales. Bulletin of the American Mathematical Society 56 65–85.
- [91] Gantert, N., Georgiou, N. and Schmid, D. (2021) The TASEP on Galton-Watson trees. Elec. J. Probab. 26 1–38. .
- [92] Gantert, N., Nestoridi, E., Schmid, D. (2023) Mixing times for the simple exclusion process with open boundaries. Ann. Appl. Probab. 33 1172–1212.
- [93] Gantert, N., Schmid, D. (2020) The speed of the tagged particle in the exclusion process on Galton-Watson trees. Elect. J. Probab. 25 1–27.
- [94] Giardiná, C., Redig, F. (2025) Duality for Markov Processes. Grundlehren der mathematischen Wissenschaften 365 Springer, Cham.
- [95] Gonçalves, P., Jara, M. (2008) Scaling limits for gradient systems in random environment. J. Stat. Phys. 131, 691–716.
- [96] Gonçalves, P., Jara, M. (2013) Scaling limits of additive functionals of interacting particle systems. Commun. Pure Appl. Math. 66 649–677.
- [97] Gonçalves, P. and Jara, M. (2014) Nonlinear fluctuations of weakly asymmetric interacting particle systems. Archive for Rational Mechanics and Analysis 212 597–644.
- [98] Gonçalves, P., Jara, M., Sethuraman, S. (2015) A stochastic Burgers equation from a class of microscopic interactions. Ann. Probab. 43 286–338.
- [99] Gonçalves, P., Perkowski, N., Simon, M. (2020) Derivation of the stochastic Burgers equation with Dirichlet boundary conditions from the WASEP. Annales Henri Lebesgue 3 87–167.
- [100] Griffeath, D. (1993) Frank Spitzer’s pioneering work on interacting particle system. Ann. Probab. 21 608-621.
- [101] Grunewald, N., Otto, F., Villani, C., and Westdickenberg, M. (2009) A two-scale approach to logarithmic Sobolev inequalities and the hydrodynamic limit. Ann. IHP Prob. Stat. 45, no. 2, 302–351.
- [102] Gu, Y. (2020) Gaussian fluctuations from the 2D KPZ equation. Stoch. Partial Differ. Equ. Anal. Comput. 8 150–185.
- [103] Gu, Y., Komorowski, T. (2023) Another look at the Balázs-Quastel-Sepp aläinen theorem. Trans. Amer. Math. Soc. 376 2947–2962.
- [104] Gubinelli, M., Perkowski, N. (2018) Energy solutions of KPZ are unique. Journal of the American Mathematical Society 31 427–471.
- [105] Gubinelli, M., Perkowski, N. (2017) KPZ reloaded. Commun. Math. Phys. 349 165–269.
- [106] Guo, M.Z., Papanicoloau, G., Varadhan, S.R.S. (1988) Nonlinear diffusion limit for a system with nearest-neighbor interactions. Commun. Pure Appl. Math. 118 31-59.
- [107] Hairer, M. (2014) A theory of regularity structures. Inventiones Mathematicae 198 269–504.
- [108] Helland, I.S. (1982) Central limit theorems for martingales with discrete or continuous time. Scand. J. Statist. 9 79–94
- [109] Hermon, J., Salez, J. (2019) A version of Aldous? spectral-gap conjecture for the zero range process. Ann. Appl. Probab. 29 2217–2229.
- [110] Hernández-Hernández, E. M., Jacka, S.D. (2022) A generalisation of the Burkholder-Davis-Gundy inequalities. Elec. Commun. Probab. 27 #50, 1–8.
- [111] Holley, R. (1970) A class of interaction in an infinite particle system. Adv. in Math. 5 291-309.
- [112] Holley, R.; Stroock, D. (1978) Generalized Ornstein-Uhlenbeck processes and infinite particle branching Brownian motions. RISMA Kyoto Publ. 14 741-788.
- [113] Huang, R., Matetski, K., Weber, H. (2025) Scaling limit of a weakly asymmetric simple exclusion process in the framework of regularity structures. arXiv:2505.00621.
- [114] Hudiani, M., Landim, C., Sethuraman, S. (2025) Singular diffusion limit of a tagged particle in zero range processes with Sinai-type random environment. arXiv:2502.17365
- [115] Gubinelli, M., Imkeller, P., Perkowski, N. (2015) Paracontrolled distributions and singular PDEs. Forum of Mathematics, Pi 3 e6.
- [116] Jacod, J., Shiryaev, A. (2013) Limit Theorems for Stochastic Processes 288 Springer Science and Business Media.
- [117] Jara, M.D. (2009) Nonequilibrium scaling limit for a tagged particle in the simple exclusion process with long jumps. Commun. Pure Appl. Math. 62 198–214.
- [118] Jara, M.; Landim, C. (2006) Nonequilibrium central limit theorem for a tagged particle in symmetric simple exclusion. Ann. Inst. Henri Poincaré Prob. et Statistiques 42 567-577.
- [119] Jara, M., Landim, C. (2023) The stochastic heat equation as the limit of a stirring dynamics perturbed by a voter model. Ann. Appl. Probab. 33 4163–4209.
- [120] Jara, M., Landim, C., Sethuraman, S. (2009) Nonequilibrium fluctuations for a tagged particle in mean-zero one dimensional zero-range processes. Probab. Theory Rel. Fields 145, 565–590.
- [121] Jara, M., Landim, C., Sethuraman, S. (2013) Nonequilibrium fluctuations for a tagged particle in one dimensional sublinear rate zero-range processes. Ann. IHP Prob. Stat. 49, 611–637.
- [122] Jara, M., Landim, C., and Teixeira, A. (2011) Quenched scaling limits of trap models. Ann. Probab. 39 176–223.
- [123] Jara, M., Menezes, O. (2018) Non-equilibrium fluctuations of interacting particle systems. arXiv:1810.09526
- [124] Jara, M., Peterson, J. (2017) Hydrodynamic limit for a system of independent, sub-ballistic random walks in a common random environment. Ann. IHP Prob. Stat. 53 1747–1792.
- [125] Johansson, K. (2000) Shape fluctuations and random matrices. Commun. Math. Phys. 209 437–476.
- [126] Jung, P. (2004) Perturbations of the symmetric exclusion process Markov Proc. Rel. Fields 10 565–584.
- [127] Karlin, S., and Taylor, H.M. (1975) A first course in stochastic processes. Second edition. Academic Press, New York-London.
- [128] Kipnis, C. (1986) Central limit theorems for infinite series of queues and applications to simple exclusion. Ann. Probab. 14 397-408.
- [129] Kipnis, C. (1987) Fluctuations des temps d’occupation d’un site dans l’exclusion simple symmetrique. Ann. I. H. Poincare Sect. B (N.S.) 23 21-35.
- [130] Kipnis, C.; Landim, C. (1999) Scaling limits of interacting particle systems. Grundlehren der Mathematischen Wissenschaften 320 Springer-Verlag, Berlin.
- [131] Kipnis, C., Varadhan, S. R. S. (1986) Central limit theorem for additive functionals of reversible markov processes. Commun. Math. Phys. 104 1-19.
- [132] Kipnis, C., Olla, S., Varadhan, S.R.S. (1989) Hydrodynamics and large deviations for simple exclusion processes. Comm. Pure Appl. Math. 42, 115–137.
- [133] Komorowski, T; Landim, C; Olla, S. (2012) Fluctuations in Markov Processes: Time Symmetry and Martingale Approximation. 345 Grundlehren der Mathematischen Wissenschaften Springer-Verlag, Berlin.
- [134] Labbé, C., Lacoin, H. (2019) Cutoff phenomenon for the asymmetric simple exclusion process and the biased card shuffling. Ann. Probab. 47 1541–1586.
- [135] Landim, C. (2004) Hydrodynamic limit of interacting particle systems. in ICTP Lecture notes. 57-100. School and Conference on Probability Theory, May 13-17, 2002. Ed. G. Lawler. Abdus Salam ICTP Publications, Trieste, Italy.
- [136] Landim, C. (2019) Metastable Markov chains. Probab. Surveys 16 143–227.
- [137] Landim, C., Lee, J., Seo, I. (2024) Metastability and time scales for parabolic equations with drift 1: the first time scale. Archive for Rational Mechanics and Analysis 248 78.
- [138] Landim, C., Pacheco, C.G., Sethuraman, S., Xue, J. (2023) On a nonlinear spde derived from a hydrodynamic limit in a Sinai-type random environment. Ann. Appl. Probab. 33, 200–237.
- [139] Landim, C; Olla, S; Varadhan, SRS (2001) Symmetric simple exclusion process: Regularity of the self-diffusion coefficient. Commun. Math. Phys. 224 307–321.
- [140] Landim, C., Olla, S., Yau, H.T. (1997) First order correction for the hydrodynamic limit of asymmetric simple exclusion processes in dimensions . Commun. Pure Appl. Math 50 149-203.
- [141] Landim, C., Olla, S., Yau, H. (2004) Diffusive behaviour of the equilibrium fluctuations in the asymmetric exclusion processes. Adv. Stud. Pure Math. 39 307–324.
- [142] Landim, C.; Sethuraman, S.; Varadhan, S.R.S. (1996) Spectral gap for zero range dynamics. Ann. Probab. 24 1871-1902.
- [143] Landim, C., Yau, H. T. (1997) Fluctuation-dissipation equation of asymmetric simple exclusion processes. Probab. Theory Rel. Fields 108 321–356.
- [144] Lenglart, E., Lépingle, D. and Pratelli, M. (1980) Présentation unifiée de certaines inégalités de la théorie des martingales. Séminaire de Probabilités 14 26–48 .
- [145] Liggett, T.M. (1973) An infinite particle system with zero-range interactions. Ann. Probab. 1 240-253.
- [146] Liggett, T.M. (1978) Random invariant measures for Markov chains and independent particle systems. Z. Wahr. Gebiete. 45, 297-313.
- [147] Liggett, T. M. (1985) Interacting Particle Systems Grundlehren der Mathematischen Wissenschaften 276 Springer-Verlag, New York.
- [148] Liggett, T.M. (1999) Stochastic Interacting Systems: Contact, Voter and Exclusion Processes. Grundlehren der Mathematischen Wissenschaften 324 Springer-Verlag, New York.
- [149] Liggett, T.M. (2004) Interacting particle systems–an introduction. in ICTP Lecture notes. 1-57. School and Conference on Probability Theory, May 13-17, 2002. Ed. G. Lawler. Abdus Salam ICTP Publications, Trieste, Italy.
- [150] Liggett, T.M. (2009) Distributional limits for the symmetric exclusion process. Stoch. Proc. Appl. 119 1–15.
- [151] Liggett, T.M. and Spitzer, F. (1981) Ergodic theorems for coupled random walks and other systems with locally interacting components. Prob. Theory Rel. Fields 56 443-468.
- [152] Lin, Y. (2023) Classification of stationary distributions for the stochastic vertex models. Elec. J. Probab. 28 1–40.
- [153] Liu, Z. (2022) Multipoint distribution of TASEP. Ann. Probab. 50 1255–1321.
- [154] Loulakis, M., Stamatakis, M.G. (2019) Generalized Young measures and the hydrodynamic limit of condensing zero range processes arXiv:1910.00493
- [155] Mallick, K., Moriya, H. and Sasamoto, T. (2022) Exact solution of the macroscopic fluctuation theory for the symmetric exclusion process. Physical Review Letters 129 040601.
- [156] Mariani, M. (2010) Large deviations principles for stochastic scalar conservation laws. Probab. Theory Rel. Fields 147 607–648.
- [157] Martin-Löf, A. (1976) Limit theorems for the motion of a Poisson system of independent Markovian particles with high density. Z. Wahr. Gebiete 34 205–223.
- [158] Matetski, K., Quastel, J., Remenik, D., (2021) The KPZ fixed point. Acta Math. 227 no. 1, 115–203.
- [159] Maxwell, M., Woodroofe, M. (2000) Central limit theorems for additive functionals of Markov chains. Ann. Probab. 28 713-724.
- [160] Mitoma, I. (1983) Tightness of probabilities on and . Ann. Probab. 11 989–999.
- [161] Morrey, C.B. (1955) On the derivation of the equation of hydrodynamics from statistical mechanics. Commun. Math. Pure. Appl. VIII 279–326.
- [162] Morris, B. (2006) Spectral gap for the zero range process with constant rate. Ann. Probab. 34 1645-1664.
- [163] Nagahata, Y. (2010) Spectral gap for zero-range processes with jump rate . Stoch. Proc. Appl. 120 949-958.
- [164] Olla, S., Varadhan, S.R.S., Yau, H.T. (1993) Hydrodynamical limit for a Hamiltonian system with weak noise. Commun. Math. Phys. 155 523-560.
- [165] Peligrad, M. and Sethuraman, S. (2008) On fractional Brownian motion limits in one dimensional nearest-neighbor symmetric simple exclusion. ALEA Lat. Am. J. Prob. Stat. 4 245-255.
- [166] Peterson, J. (2010) Systems of one-dimensional random walks in a common random environment. Electron. J. Probab. 15, 102–1040.
- [167] Petrov, V.V. (1975) Sums of Independent Random Variables, Springer Verlag, New York.
- [168] Papanicolaou, G. C.; Stroock, D.; Varadhan, S. R. S. (1977) Martingale approach to some limit theorems. Papers from the Duke Turbulence Conference (Duke Univ., Durham, N.C., 1976), Paper No. 6, ii+120 pp. Duke Univ. Math. Ser., Vol. III, Duke Univ., Durham, N.C.
- [169] Parthasarathy, K.R. (2005) Probability Measures on Metric Spaces 352 American Mathematical Soc., AMS Chelsea Publishing, Providence, R.I.
- [170] Popkov, V., Schadschneider, A., Schmidt, J., Schütz, G. (2015) Fibonacci family of dynamical universality classes. Proceedings of the National Academy of Sciences 112 12645–12650.
- [171] Quastel J. (2011) Introduction to KPZ. Current Developments in Mathematics 2011, 125–194.
- [172] Quastel J, Tsai LC. (2025) Hydrodynamic large deviations of TASEP. Commun. Pure Appl. Math. 78 913–994.
- [173] Quastel, J., Valkó, B. (2008) A note on the diffusivity of finite-range asymmetric exclusion processes on . Progress in Probability 60 54–549, Birkhauser Verlag Basel/Switzerland.
- [174] Rameh, O., (2024) Cutoff phenomenon for asymmetric zero range process with monotone rates. arXiv:2410.05220.
- [175] Ravishankar, K. (1992) Fluctuations from the hydrodynamical limit for the symmetric simple exclusion in . Stoc. Proc. Appl. 42 31-37.
- [176] Redig, F., Sau, F. (2018) Factorized duality, stationary product measures and generating functions. J. Stat. Phys. 172. 980–1008.
- [177] Reed, M. and Simon, B. (1980) Methods of Modern Mathematical Physics. I: Functional Analysis Revised and Enlarged Edition. Academic Press, New York.
- [178] Resnick, S. (1992) Adventures in stochastic processes. Birkhäuser Boston, Inc., Boston, MA.
- [179] Rezakhanlou, F. (1991) Hydrodynamic limit for attractive particle systems on . Commun.Math. Phys. 140, 417–448.
- [180] Rezakhanlou, F., (1994) Evolution of tagged particles in non-reversible particle systems. Commun. Math. Phys. 165 1–32.
- [181] Rosenblatt M. (1967) Transition probability operators, in: Proc. Fifth Berkeley Symp. Math. Statist. Prob., 2, 473–483.
- [182] Rost, H. (1981) Nonequilibrium behavior of many particle systems: density profile and local equilibria. Z. Wahrs. Verw. Gebiete, 58, 41–53.
- [183] Rost, H.; Vares, M.E. (1985) Hydrodynamics of a one dimensional nearest neighbor model. Contemp. Math. 41 329-342.
- [184] Roy, D., Dhar, A., Khanin, K., Kulkarni, M., Spohn, H. (2024) Universality in coupled stochastic burgers systems with degenerate flux jacobian. J. Stat. Mech.: Theory and Experiment 2024 033209.
- [185] Roy, D., Dhar, A., Kulkarni, M., Spohn, H. (2025) Fixed points and universality classes in coupled Kardar-Parisi-Zhang equations. J. Stat. Mech.: Theory and Experiment 2025 073209.
- [186] Saint-Raymond, L. (2009) Hydrodynamic Limits of the Boltzmann Equation Lecture Notes in Mathematics 1971 Springer, Berlin.
- [187] Salez, J. (2023) Universality of cutoff for exclusion with reservoirs. Ann. Probab. 51 478–494.
- [188] Schmidt, J., Krajnik, Z., Popkov, V. (2025) Universality in many-body driven systems with an umbilic point. J. Stat. Mech.: Theory and Experiment 2025 083202.
- [189] Seppäläinen, T. (1999) Existence of hydrodynamics for the totally asymmetric simple -exclusion process. Ann. Probab. 27 361–415.
- [190] Seppäläinen, T. (2008) Translation Invariant Simple Exclusion Processes. Available at http://www.math.wisc.edu/~seppalai/excl-book/etusivu.html
- [191] Seppäläinen, T. and Sethuraman, S. (2003) Transience of second-class particles and diffusive bounds for additive functionals in one-dimensional asymmetric exclusion processes. Ann. Probab. 31 148-169.
- [192] Sethuraman, S. (2000) A martingale central limit theorem. Available at https://sundersethuraman.github.io/notes/wi_mart1.pdf
- [193] Sethuraman, S. (2000) Central limit theorems for additive functionals of the simple exclusion process. Ann. Probab. 28 277-302; Correction (2006).
- [194] Sethuraman, S. (2001) On extremal measures for conservative particle systems. Ann. IHP Prob. Stat. 37 139-154.
- [195] Sethuraman, S., (2003) An equivalence of H(-1) norms for the simple exclusion process. Ann. Probab. 31 35–62.
- [196] Sethuraman, S. (2006) Superdiffusivity of occupation-time variance in 2-dimensional asymmetric exclusion processes with density 1/2. J. Stat. Phys. 123 787–802.
- [197] Sethuraman, S. (2006) Diffusive variance for a tagged particle in asymmetric simple exclusion. ALEA 1 305-332.
- [198] Sethuraman, S. (2007) On diffusivity of a tagged particle in asymmetric zero-range dynamics. Ann. IHP Prob. Stat. 43 215–232.
- [199] Sethuraman, S., Shahar, D. (2018) Hydrodynamic limits for long-range asymmetric interacting particle systems. Electron. J. Probab. 23 1–54.
- [200] Sethuraman, S., Varadhan, S.R.S. (2013) Large deviations for the current and tagged particle in 1D nearest-neighbor symmetric simple exclusion. Ann. Probab. 41 1461–1512.
- [201] Sethuraman, S., Varadhan, S.R.S. (2025) Atypical behaviors of a tagged particle in asymmetric simple exclusion. Ann. Probab. 53 1857–1896.
- [202] Sethuraman, S., Varadhan, S.R.S., Yau, H.T. (2000) Diffusive limit of a tagged particle in asymmetric simple exclusion processes. Commun. Pure and Appl. Math. 53 972-1006.
- [203] Sethuraman, S., Xu, L. (1996) A central limit theorem for reversible exclusion and zero-range particle systems. Ann. Probab. 24 1842-1870.
- [204] Sethuraman, S. and Xue, J. (2024) Condensation, boundary conditions, and effects of slow sites in zero-range systems. Ann. Probab. 52 1048–1092.
- [205] Spitzer, F. (1970) Interaction of Markov processes. Adv. in Math. 5 246-290.
- [206] Spohn, H. (1991) Large Scale Dynamics of Interacting Particles. Springer-Verlag, Berlin.
- [207] Spohn, H. (2020) The 1+ 1 dimensional Kardar-Parisi-Zhang equation: more surprises. J. Stat. Mech.: Theory and Experiment 4 04400
- [208] Spohn, H., Stoltz, G. (2015) Nonlinear fluctuating hydrodynamics in one dimension: The case of two conserved fields. J. Stat. Phys. 160 861–884.
- [209] Tóth, B. (2024) Central limit theorem for random walks in divergence free random drift field – revisited. Revue Roumaine de Mathématiques Pures et Appliquées 69 585–601.
- [210] Tóth, B. (2025) Semi-quenched invariance principle for the random Lorentz gas: Beyond the Boltzmann?Grad limit. Entropy 27 397.
- [211] Tracy, C.A., Widom, H. (2009) On ASEP with step Bernoulli initial condition. J. Stat. Phys. 137 825–838.
- [212] Varadhan, S. R. S. (1995) Self-diffusion of a tagged particle in equilibrium for asymmetric mean zero random walk with simple exclusion. Ann. I.H.P. Prob. et Stat. 31 273-285.
- [213] Varadhan, S.R.S. (2000) Lectures on hydrodynamic scaling. in Hydrodynamic Limits and Related Topics. Fields Institute Communications 27 3–42.
- [214] Varadhan, S.R.S. (2004) Large deviations for the asymmetric simple exclusion process. In Stochastic analysis on large scale interacting systems. 39 Mathematical Society of Japan 1–28.
- [215] Waymire, E. (1976) Zero-range interaction at Bose-Einstein speeds under a positive recurrent single particle law. Ann. Probab. 8 441–450.
- [216] Whitt, W. (2007) Proofs of the martingale FCLT. Probability Surveys 4 (2007) 268–302.
- [217] Xu, T., Zhao, L. (2025) Nonequilibrium fluctuations for the occupation time of the SSEP in . arXiv:2512.09424
- [218] Xue, X., Zhao, L. (2024) Sample path MDP for the current and the tagged particle in the SSEP. Elec.. J. Probab. 29 1–39.
- [219] Yang, K. (2023) Fluctuations for some nonstationary interacting particle systems via Boltzmann?Gibbs principle. Forum of Mathematics, Sigma 11, e32, 1–86.
- [220] Yang, K. (2025) Time-inhomogeneous KPZ equation from non-equilibrium Ginzburg-Landau SDEs. Elec. J. Probab. 30 1–155.
- [221] Yau, H.T. (1991) Relative entropy and hydrodynamics of Ginzburg-Landau models. Lett. Math. Phys. 22 63–80.
- [222] Yau, H.T. (2004) The law of the two dimensional asymmetric simple exclusion process. Annals of Mathematics 159 377–405.
- [223] Zhou, Z. (2020) Hydrodynamic limit for a d-dimensional open symmetric exclusion process Elec. Commun. Probab. 25 (2020) #76 1--8.