A note on the Fredrickson-Andersen one spin facilitated model in stationarity
Abstract.
This note discusses three problems related to the Fredrickson-Andersen one spin facilitated model in stationarity. The first, considered in 2008 in a paper of Cancrini, Martinelli, Roberto and Toninelli, is the spectral gap of the model’s infinitesimal generator. They study the decay of this spectral gap when the density is large, but in dimensions and higher, they do not find the exact exponent. They also show that the persistence function of the model has exponential tail, but the typical decay time is not analyzed. We will see that the correct exponent for the decay of the spectral gap in dimension and higher is , and discover how the time over which the persistence function decays diverges in high densities. We also discuss the scaling of the spectral gap in finite graphs.
1. Introduction and results
The purpose of this note is to present three small results on the Fredrickson-Andersen one spin facilitated model (FA1f), following [2], addressing two problems that have not been solved there and one tightly related problem in a slightly different setting studied in [4, 5, 6]. Since it is, in a sense, an extension of [2], the reader is referred to [2] for the relevant background, references, and complete introduction of the model and notation.
We will only briefly remind here that sites in could be either occupied or empty, with equilibrium probabilities and respectively (where is thought of as small). When a site has at least one empty neighbor, it is being resampled from equilibrium with rate , and otherwise its occupation cannot change. The process is reversible with respect to the invariant measure , given by an independent product of Bernoulli random variables with parameter . Probabilities and expected values with respect to the stochastic process are denoted and , where the subscript indicates that the initial configuration is drawn from equilibrium.
The first result here completes Theorem 6.4 of [2], which bounds the spectral gap of the FA1f model. It is shown there that the gap decays polynomially as the parameter tends to , and for dimensions and the exact exponent is identified, up to a logarithmic correction in dimension . For dimension , however, the exponent is bounded between and , and its exact value is not determined. The following theorem shows that the correct scaling is –
Theorem 1.1.
Consider the setting of [2, Theorem 6.4], in dimension . Then there exists a positive constant (possibly depending on ) such that
The second result presented here concerns with the persistence function. Recall that
In general, when the spectral gap is positive, , where could be chosen to be equal . We will see that for the FAf model on this choice is not optimal, and that the typical time to empty the origin scales (perhaps with lower order corrections) like the inverse of the spectral gap, without the additional factor of .
Theorem 1.2.
Consider the FAf model. Then there exists a positive constant such that
Moreover, for a different positive constant ,
The last theorem that will be presented is in a slightly different setting – the Fredrickson-Andersen model on a finite graph . A particularly interesting case, studied in [5, 6] and more recently in [4], is when for some positive constant , where denotes the set of vertices of . A lower bound on the spectral gap is given in [5, 6] and later on in [4] by suggesting a relaxation mechanism in which vacancies travel as random walkers on . The next theorem will bound the spectral gap of this model from above, showing that this mechanism has a leading contribution. Consider two independent continuous time random walks on , namely, each of the two walkers moves to each neighboring site with rate . For two vertices , we define to be the expected time that it takes for two such random walkers starting at and to reach distance at most . Let be its expected value, starting at two random positions, i.e., .
Theorem 1.3.
Consider the FAf model on a finite graph , and assume that for some positive constant . Let denote its spectral gap with respect to the product measure conditioned on having at least one vacancy. Then
for some positive constant .
Remark 1.4.
[4] give a lower bound on the spectral gap, which in various graphs is of the same order as the upper bound given in Theorem 1.3. In particular, on the two dimensional torus , , both the upper and lower bounds scale like . In view of this scaling, and the relaxation mechanism reflected in its proof, it seems that the correct scaling of the spectral gap is also in , coinciding with the lower bound of [2]. Unfortunately, the ideas in the proof of Theorem 1.3 do not seem to be easily adapted for the model on , and the problem remains open.
2. Notation
We will now recall some of the notation in [2] that will be used in this note.
-
•
For , and .
-
•
For and we denote (where the -dependence of is implicit).
-
•
The configuration space is , and the measure on this space is a product measure of Bernoulli random variables with parameter .
-
•
For and , the configuration which equals outside and different from at is denoted .
-
•
The FAf critical length is denoted , where is given in [2, Theorem 4.1], and does not depend on .
-
•
The constraint of the FAf dynamics, for a configuration and , is
-
•
The Dirichlet form of FAf operating on a local function is given by
For the FAf model on a finite graph with vertex set denote ; and the product measure of Bernoulli random variables with parameter , conditioned on having at least one empty site. The constraint is defined in the same manner as in , except that should be understood as the graph distance between and , that we denote . The Dirichlet form operating on is given by
Throughout the proof will denote a generic positive constant, and is assumed to be small enough.
3. Proof of Theorem 1.1
In order to bound the spectral gap from above, we need to find an appropriate test function , such that
Consider the box for . For a configuration and a site , the connected cluster of , denoted , is defined as the set of sites that are connected to via a path of empty sites in . If , its connected cluster is the empty set. This way, the set of empty sites in is partitioned in connected clusters, and we define:
(3.1) |
Proposition 3.1.
For the test function defined in equation (3.1),
(3.2) |
Proof.
This result is shown in [3] for the case of Bernoulli bond percolation. We will repeat their argument applied to our case for completeness.
First, note that we may write
where, when (and therefore ), we define .
Let , and for , define to be the indicator of the event, that the set is equal . Note that . For such a set , let , and define
When is such that ,
In order to use this identity, we split the variance over the different choices of :
Consider one of the summands in the above expression –
The first term is positive, and we will simply bound it by . In order to find the second term, we note that the variables and are independent, and therefore
Finally, since under the event the variables and are independent, the third term vanishes. Therefore,
This establishes inequality (3.2). ∎
Proposition 3.2.
For the test function defined in equation (3.1),
(3.3) |
Proof.
Recall first that
Consider a single term in that sum. First, note that by flipping a single site could change by at most . If is outside , flipping it could not change the number of clusters in and its contribution would be . If is on the boundary of (i.e., it is in and has a neighbor outside ), then
Finally, if is in but has no neighbors outside , the number of open clusters could only change if it has at least two empty neighbors –
The proof is now concluded by summing these options –
4. Proof of Theorem 1.2
4.1. Upper bound
The basic tool for the proof of the upper bounds on is the following result of [1] (see also [8, Section 4]):
Lemma 4.1.
Assume that, for some and any local function which vanishes on the event ,
(4.1) |
Then .
We will use a path argument, similar to [2, proof of Proposition 6.6], proving inequality (4.1) with the appropriate .
We start by defining a canonical geometric path, which is a discrete approximation of a straight segment. More precisely, for any , we will construct a nearest neighbor path with and whose distance from the line segment is small. The exact definition is rather cumbersome, and a reader who accepts that such a path could be constructed satisfying Observation 4.3 and Claim 4.6 (see Definition 4.5) may skip the technicalities involved in their proofs.
Definition 4.2.
Fix with . The canonical geometric path connecting to the origin is the path constructed as follows – consider the set defined as
For each there are values of for which , hence . We will order according to the lexicographic order, , so that for all , and in case of equality . Then
Observation 4.3.
Fix and . Then , i.e., the sites of the path are indexed by their norm.
Claim 4.4.
Proof.
We show this by induction. Start with , and consider the vector . By the construction of ,
and indeed .
For , there are two options – either and , or . In the first case, by induction and letting ,
Since the coordinates between and are not integer, we can replace by , and since is integer we may replace by . That is, .
Let us now consider the second case, where . First, by induction, noting that the coordinates after of are not integer,
On the other hand, we know that
so and the proof is complete. ∎
Definition 4.5.
Fix and with . The the -cone of is the set
Claim 4.6.
Fix and such that . Then .
Proof.
First, since for the cone consists of the points in of norm smaller than , its size is smaller than , so in what follows we may assume ; and in this case we will show that the stronger inequality holds.
Let , i.e., , so by Claim 4.4 there exist and such that and
Assume first , so in particular by the construction of the geometric path. must be contained in ; so we fix and let , such that
for some . That is, for all ,
allowing at most integer choices of . Finally, since , necessarily , so, still for , the number of possibilities for is bounded by
Finally, summing over all possible values of ,
For any , let
Claim 4.7.
Fix and such that . Then there exists a path of configurations and a sequence of sites such that:
-
(1)
and .
-
(2)
For any , , and
-
(3)
The sites all belong to the geometric path . Moreover, each site of appears at most twice in the sequence , and in particular .
-
(4)
For all , the number of sites in which are empty for is at most the number of sites in which are empty for .
-
(5)
Fix , , . Then there exist at most one configuration and one index such that , and . We write .
Proof.
The path is constructed in the same manner as [2, proof of Proposition 6.6] – let be an empty site in with minimal -norm, and denote . Then set, for ,
This sequence defines a path that indeed satisfied the conditions of the claim. ∎
Proposition 4.8.
For any local function that vanishes on the event and for any ,
Proof.
First, since is local, we may restrict ourselves to proving the inequality for FAf on a large finite set , so the configuration space is . This allows us to write the Dirichlet form as
Consider, for any , the path constructed in Claim 4.7, and for set
Note that we can bound, uniformly in ,
and by Claim 4.6, for every ,
By the Cauchy-Schwarz inequality and the properties of the path,
(4.2) | ||||
Note that we are allowed to divide by since , and hence it is non-zero. We can estimate more precisely:
so
By property 4 of the path we obtain
We now conclude by continuing the estimate (4.2) –
Remark 4.9.
If, rather than in inequality (4.1), we would like to bound , we could use Proposition 4.8 directly. By the Cauchy-Schwarz inequality
Choosing with small enough, is bounded below , so
This inequality is not entirely worthless, and it does bound from above by (see [8, Section 4] and equation (4.3) in the following section). However, in order to obtain the exponential tail in Theorem 1.2 a more sophisticated approach is required.
From now on we set . The lower bound on the spectral gap of [2, Theorem 6.4] is proven by introducing an auxiliary dynamics with large spectral gap and then comparing it to the FAf dynamics. For that objective they define the constraints , stating that none of the boxes is entirely occupied for ; with the associated Dirichlet form
The following Lemma is given in [2, equation (5.1)]:
Lemma 4.10.
The spectral gap associated with is at least .
Then, [2, Theorem 6.4] is proved by showing:
Proposition 4.11.
We will use Lemma 4.10 in order to prove the following claim:
Claim 4.12.
Assume that is small enough, and let be a function vanishing on the event . Then
Proof.
We are now ready to prove the upper bound on the persistence function – consider which vanishes on . Then, by Claim 4.12 and the fact that does not depend on the occupation in for ,
4.2. Lower bound
In order to bound the expected value of from below, we will use the following variational principle (see [8, Proposition 4.7]):
Lemma 4.13.
Let be the space of local functions that vanish on the event . Then
Remark 4.14.
It would be more convenient to use a homogeneous version of that variational principle –
and since for fixed the expression is maximized for ,
(4.3) |
We will now treat separately three different cases: , , and .
4.2.1.
For high dimensions, we will use the test function Its expected value is , and its Dirichlet form is given by Equation (4.3) now concludes the proof of this case. ∎
4.2.2.
In the one dimensional case we will use a test function similar to [2, proof of Theorem 6.4]. Let ,
and
(4.4) |
Proposition 4.15.
Consider defined in equation (4.4). Then
Proof.
First, note that is a geometric random variable with parameter , so we can calculate explicitly
and since is positive . ∎
Proposition 4.16.
Consider defined in equation (4.4). Then
Proof.
In order to bound we make the following observations –
-
(1)
For fixed , if then either , or .
-
(2)
For fixed , if then .
With these observations in mind,
Using these two propositions and equation (4.3), the case is concluded.∎
4.2.3.
Let , recalling , and . The test function we will use is
(4.5) |
Note that it vanishes on the event , and that it depends only on the occupation in .
Remark 4.17.
The function is used in [7] in a different context, in order to bound the relaxation time of the simple random walk on a certain graph that consists of two copies of (for some ). Though presented differently, the proof there is based on the fact that this function serves as a test function for the hitting time at of the random walk on ; and that for the dynamics to relax a random walk in one of the two copies of must first hit . Indeed, the bound obtained scales as the expected hitting time at the origin for a random walk in with jump rate .
Proposition 4.18.
Consider defined in equation (4.5). Then
Proof.
The proof is based on the fact, that the probability that is entirely occupied, given by , is bounded away from uniformly in (thanks to the choice ). In this case, equals , which is greater than . ∎
Proposition 4.19.
Consider defined in equation (4.5). Then
Proof.
The proof is based on the following observation:
Observation 4.20.
Fix and such that , , and . Then and .
Proof.
Since , there can be no vertex with and , so in particular . Moreover, since , it must have an empty neighbor , and since this neighbor has norm greater than , necessarily . Since, in addition, no empty site for has norm strictly smaller than , we conclude that . ∎
This observation implies in particular that for all and such that
Using this estimate,
5. Proof of Theorem 1.3
As in the proof of Theorem 1.1, we look for such that
where the variance is understood with respect to the measure .
The test function that we use is
(5.1) |
Before analyzing the function , note that solves the following Poisson problem:
where is the infinitesimal generator of two independent random walks on . Multiplying both sides by and averaging over and we obtain
(5.2) |
where the Dirichlet form is given for every by
Proposition 5.1.
For defined in equation (5.1),
Proof.
Under the probability that exactly one site is empty is of order (i.e., bounded away from uniformly in ). When this happens, , so
Moreover, the probability that there are exactly two vacancies is also of order . Under this event,
Proposition 5.2.
For defined in equation (5.1),
Proof.
We start the proof with an observation:
Observation 5.3.
Let and such that , , and . Then there exist such that , , , and
Proof.
First, recalling equation (5.1), when filling an empty site could only decrease, and since necessarily . Moreover, could only change if the maximum is attained at the pair for some , i.e., . Note that is non-zero, hence . Finally, means that has an empty neighbor ; and since in the configuration both and are empty . ∎
As a consequence of this observation, for all and such that ,
We can now use this estimate and calculate the Dirichlet form:
and the proposition follows form equation (5.2). ∎
Theorem 1.3 is a consequence of Propositions 5.1 and 5.2, using the variational characterization of the spectral gap. ∎
Remark 5.4.
Acknowledgments
I wish to thank Ivailo Hartarsky and Fabio Martinelli for the discussions and comments.
References
- [1] Amine Asselah and Paolo Dai Pra, Quasi-stationary measures for conservative dynamics in the infinite lattice, Ann. Probab. 29 (2001), no. 4, 1733–1754. MR 1880240
- [2] Nicoletta Cancrini, Fabio Martinelli, Cyril Roberto, and Cristina Toninelli, Kinetically constrained spin models, Probab. Theory Related Fields 140 (2008), no. 3-4, 459–504. MR 2365481
- [3] J. Theodore Cox and Geoffrey Grimmett, Central limit theorems for associated random variables and the percolation model, Ann. Probab. 12 (1984), no. 2, 514–528. MR 735851
- [4] Ivailo Hartarsky, Fabio Martinelli, and Cristina Toninelli, Coalescing and branching simple symmetric exclusion process, arXiv preprint arXiv:2006.01426 (2020).
- [5] Natesh S. Pillai and Aaron Smith, Mixing times for a constrained Ising process on the torus at low density, Ann. Probab. 45 (2017), no. 2, 1003–1070.
- [6] by same author, Mixing times for a constrained Ising process on the two-dimensional torus at low density, Ann. Inst. H. Poincaré Probab. Statist. 55 (2019), no. 3, 1649–1678.
- [7] Laurent Saloff-Coste, Lectures on finite markov chains, Lecture Notes in Mathematics, vol. 1665, pp. 301–413, Springer, 1996.
- [8] Assaf Shapira, Kinetically constrained models with random constraints, Ann. Appl. Probab. 30 (2020), no. 2, 987–1006.