arXiv:1709.10424v2 [math.PR] 25 Mar 2018

Tulasi Ram Reddy · Sreekar Vadlamani · D. Yogeshwaran

Central limit theorem for exponentially quasi-local statistics of
spin models on Cayley graphs
Abstract Central limit theorems for linear statistics of lattice random fields (including spin models) are usually
proven under suitable mixing conditions or quasi-associativity. Many interesting examples of spin models do not satisfy mixing conditions, and on the other hand, it does not seem easy to show central limit theorem for local statistics
via quasi-associativity. In this work, we prove general central limit theorems for local statistics and exponentially
quasi-local statistics of spin models on discrete Cayley graphs with polynomial growth. Further, we supplement
these results by proving similar central limit theorems for random fields on discrete Cayley graphs taking values in
a countable space, but under the stronger assumptions of α -mixing (for local statistics) and exponential α -mixing
(for exponentially quasi-local statistics). All our central limit theorems assume a suitable variance lower bound like
many others in the literature. We illustrate our general central limit theorem with specific examples of lattice spin
models and statistics arising in computational topology, statistical physics and random networks. Examples of clustering spin models include quasi-associated spin models with fast decaying covariances like the off-critical Ising
model, level sets of Gaussian random fields with fast decaying covariances like the massive Gaussian free field and
determinantal point processes with fast decaying kernels. Examples of local statistics include intrinsic volumes, face
counts, component counts of random cubical complexes while exponentially quasi-local statistics include nearest
neighbour distances in spin models and Betti numbers of sub-critical random cubical complexes.
Keywords Clustering spin models · mixing random fields · central limit theorem · Cayley graphs · fast decaying
covariance · Gaussian random field · determinantal point process · exponentially quasi-local statistics · cubical
complexes · nearest-neighbour graphs
Mathematics Subject Classification (2000) 82B20 Lattice systems (Ising, dimer, Potts, etc.) and systems on
graphs · 60G60 Random fields · 60F05 Central limit theorems and other weak theorems · 60D05 Geometric
probability and stochastic geometry.
DY’s research was supported by DST-INSPIRE Faculty fellowship and CPDA from the Indian Statistical Institute.
Tulasi Ram Reddy
Division of Sciences, New York University Abu Dhabi
E-mail: tulasi@nyu.edu
Sreekar Vadlamani
TIFR Center for Applicable Mathematics, Bangalore, and
Department of Statistics, Lund University
E-mail: sreekar@tifrbng.res.in
D. Yogeshwaran
Theoretical Statistics and Mathematics Unit, Indian Statistical Institute, Bangalore.
E-mail: d.yogesh@isibang.ac.in

2

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

1 Introduction:
Our results in full generality are on finitely generated countably infinite Cayley graphs but we shall restrict ourselves
to the simple case of integer lattices i.e., Zd , to motivate and illustrate our results in the introduction. A simple model
of random lattice networks is to consider a random subset of Zd (the d-dimensional integer lattice) as nodes and
define edges or other features depending on the geometry of the nodes. Under such a set-up, various performance
measures/functionals/statistics of the network reduce to geometric functionals/statistics of the underlying set of rand
dom nodes. Mathematically, the set of random nodes is a random element in {0, 1}Z with 1 denoting the presence
of a node at the corresponding location in Zd and 0 denoting absence of a node. Denoting the random node set
as P, various geometric and topological features of the random network are encapsulated in the random subset
(called as cubical complex) C(P) = ∪x∈P (Q + x) with Q = [−1/2, 1/2]d being the centred unit cube and + denoting the Minkowski sum. A common way to understand C(P) is to investigate the asymptotics of C(Pn ) where
Pn = P ∩ [−n, n]d . Though the mathematical model described above is possibly one of the simpler ones, it appears
in various avatars in diverse areas ranging from statistical physics, digital geometry, cosmology, stereology etc.
One of the oft-used model for lattice random networks (cf. [43,32]) is the percolation model arising from statistical
physics ([39]). We note here that percolation models encode only pairwise interaction between nodes and indeed,
it is very common for various models of complex networks to be modelled as graphs using pairwise interactions
between the nodes (see [13, Section 6]). However, in many of these models the interactions are not just pairwise, but
happen between subsets of nodes, referred to as higher-order interactions. Hypergraphs represent a natural choice
for modelling such higher-order interactions, and so do cubical complexes. The additional topological structure of
cubical complexes also makes them suitable for modelling of surfaces or other topological structures as well. A hypergraph can be constructed easily from a cubical complex by making {z1 , . . . , zk } a k-hyperedge if ∩ki=1 (Q+zi ) 6= 0.
/
Due to the choice of Q, we shall have at most 2d-hyperedges on the hypergraph constructed from C(P). However,
it is possible to consider any other compact subset of Rd instead of Q to build such hypergraphs, and obtain more
general hypergraphs. But for ease of illustrating our main theorems, we shall restrict ourselves to this simple model
of cubical complex or the corresponding hypergraph. Both cubical complexes and hypergraphs are very useful models of complex networks. For example, [2,3] proposed simplicial complexes to model social networks, and we refer
the reader to [59] for a survey on this line of research using simplicial complexes for network analysis. Though
simplicial complexes are different from cubical complexes, our methods are applicable to simplicial complexes
built on spin models as well. In [54], hypergraphs have also been proposed as models for cellular networks such as
protein interaction, reaction and metabolic networks. Apart from these examples, further examples are discussed in
[29,69] including food webs and collaboration networks as well. Analysis of such networks involve understanding
of the corresponding cubical complex or its generalisation (see Section 2.2). Though we shall not explicitly state
applications of our results to such performance metrics, we hope that it will be conspicuous to the reader that our
applications in Section 2.2 can be extended naturally to handle other hypergraph statistics as well.
In the community of statistical physicists, where P is termed as a spin model, there is considerable interest in
understanding the connectivity properties of spin models ([27,39,33]). This is all the more important when one
considers random cubical complexes as models of discrete random surfaces as in [34,35]. With this perspective, it
is natural to study connectivity properties of surfaces. In Sections 2.2.2 and 2.2.3, we shall mention asymptotics of
nearest neighbour distances and Betti numbers of C(Pn ). The latter are a measure of high-dimensional connectivity
of surfaces.
Besides being considered as models of complex networks or discrete surfaces, cubical or cell complexes have often
been used in digital image analysis. A simple digital image is nothing but black and white values (i.e., 0, 1-values)

CLT for statistics of spin models

3

assigned to lattice points. Geometry and topology of digital images are of interest in image processing ([38,57,56,
81]), morphology ([83]), cosmology and stereology ([46,88]). We shall explicitly mention some of the statistics (for
example, subgraph counts, component counts) in the Section 2.2.1 and these are motivated by those that arise in the
above referenced literature. For example, recently Minkowski tensors of random cubical complexes have received
attention in astronomy (see [37,55]). Minkowski tensors are a generalisation of intrinsic volumes and asymptotics
of the latter shall concern us in Section 2.2.1.
It is pertinent to mention here that in [21, Section 6], it was said about random fields on Rd that “It would be desirable to prove limit theorems for joint distributions of various surface characteristics of different classes of random
fields . . . it might be of interest to prove limit theorems involving other Minkowski functionals for level sets such as
the boundary length or the Euler characteristics.” Such a question remains unproven even for random fields parameterised by Zd . As level-sets of random fields are an example of spin models (see below), we can safely remark
that our general central limit theorems have reduced the task of proving a central limit theorem for many statistics
(including those mentioned above) of various spin models on many Cayley graphs to that of proving variance lower
bounds. Indeed, this represents one of our major contributions.
We remark here that other statistics of interest arising naturally in signal-to-interference-plus-noise networks (see
[43,32]), k-nearest neighbour edges, local topological numbers as in [81, Section VI], discrete Morse critical points
as in [31] and local porosity as in [46], also fall within our framework.
Even though the afore-cited applications are important contributions of our article as well as crucial motivations for
our work, our main contribution can be said to lie within the realm of limit theorems for asymptotically independent
random fields as in [87,25,18,45,20,8]. We shall elaborately survey this literature to contextualise our results. We
shall make more specific remarks on the relation between this literature (see Remark 1.14) and our contribution after
stating our results.
A lattice random field X = {Xz : z ∈ Zd } is a collection of S -valued random variables (S is countable) indexed by
the lattice points. Various statistics of the random field can be expressed as sums of score functions ξ (x, X ) which
encode the interaction of Xz with X . Setting Xn = {Xz : z ∈ [−n, n]d }, one is interested in the asymptotics for
Hnξ :=

∑

ξ (z, Xn ).

(1.1)

z∈[−n,n]d

Not surprisingly, a large section of literature on random fields is devoted to the simplest of score functions i.e,
ξ
ξ (z, Xn ) = Xz . In what follows, we shall call the corresponding Hn as linear statistic. In contrast to linear statistic, two other general statistics - local statistics and exponentially quasi-local statistic - have received very little
attention. The two are defined rigorously in Section 1.2. Briefly, ξ is a local statistic if there exists a (deterministic)
r ∈ N such that ξ (x, {Xz }z∈x+[−n,n]d ) = ξ (x, {Xz }z∈x+[−r,r]d ) for all n > r. Such statistics are also U-statistics. In
vague terms, ξ is exponentially quasi-local if r is allowed to be random with suitably decaying tail probability. It
is to be understood that we refer to the central limit theorem (often abbreviated as CLT) whenever we use the word
asymptotics or limit theorems below.
A spin model P can be naturally viewed as a random field by setting Xz = 1[z ∈ P]. Under such a consideration,
all the aforementioned statistics are either local or exponentially quasi-local statistics. Spin models also arise as
level-sets of random fields i.e., set P := {z ∈ Zd : Xz ≥ u} for some u ∈ R. Such level-sets are of interest as well in
the applications mentioned above and represent yet another use for our results.

4

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

While asymptotics for linear statistics of X follow naturally when X is an i.i.d. random field, asymptotics of
local or exponentially quasi-local statistics can be deduced from the powerful martingale-difference based central
limit theorem of [76]. For example, we refer the reader to recent application of this general central limit theorem
to random topology in [47]. Asymptotics for local statistics can also be deduced from the results for U-statistics of
i.i.d. random fields. Once we drop the assumption of independence, the classical methods or even those in [76] fail.
A natural heuristic is that stationary ‘asymptotically’ independent random fields shall display the same asymptotic
behaviour as i.i.d. random fields. However it is a challenge to theorise a notion of ‘asymptotic independence’ that is
powerful enough to prove various limit theorems and yet accommodative enough to include many interesting examples. To the best of our knowledge, two successful approaches to capture the notion of ‘asymptotic independence’
in random fields are mixing and (BL, Θ )-dependence.
In short, mixing conditions require that the ‘distance’ between sigma-algebras of {Xz : z ∈ A} and {Xz : z ∈ B} for
A, B ⊂ Zd decay suitably fast as a function of the distance between A and B. We shall mostly focus on α -mixing in
this article and refer to Section 1.3 for more details.
(BL, Θ )-dependence, on the other hand, requires ‘suitable decay’ of covariance of f (Xz : z ∈ A) and g(Xz : z ∈ B) for
A, B ⊂ Zd far apart, and bounded Lipschitz functions f , g. The ‘suitable decay’ allows for dependence on cardinalities of A, B, but in a specific manner, and this indeed is the obstacle in relating mixing to (BL, Θ )-dependence. We
refer to Section 1.4 for details.
However, we also shall use the weaker notion of clustering arising in statistical physics (cf. [67,68]) which roughly
states that, for disjoint sets A, B ⊂ Zd
|P(A ∪ B ⊂ P) − P(A ⊂ P)P(B ⊂ P)|
decays exponentially fast as a function of the distance between A and B with the rate of decay allowed to depend on
cardinalities of A, B with a lot more flexibility (see Definition 1.3). These different notions of ‘asymptotic independence’ and their relations are summarised in the following figure.
ρ -mixing

Quasi-association
with summable covariances

α -mixing

(BL, Θ ) dependence
under + condition

under + condition

clustering of spin models
Fig. 1 Relationships between mixing, quasi-associativity, (BL, Θ ) and clustering. Here + condition stands for the assumption that the
random field is a spin model and coefficient is fast decaying as in the Definition 1.3. See Section 2 for precise statements and proofs of
these relations.

While asymptotics for linear statistics are proven under all three conditions - α -mixing ([16]), (BL, Θ )-dependence
([20]) and clustering ([67,68]) - the non-degeneracy of the limit requires additional assumptions. However, the limit
theorems as stated in any of these papers do not apply to local or exponentially quasi-local statistics.
ξ

Thus, in order to use α -mixing or (BL, Θ )-dependence or clustering to prove a CLT for Hn , one needs to show that
α -mixing or (BL, Θ )-dependence or clustering holds for the random field {ξ (z, Xn )}z∈Wn under appropriate assump-

CLT for statistics of spin models

5

tions on the random field X and the score functions ξ . If we consider ξ generating local statistics, such an approach
can be carried out successfully for α -mixing random field X (see Theorem 1.12). The (BL, Θ )-dependence argument does not apply to local statistics i.e., the random field {ξ (x, Xn )} need not be (BL, Θ )-dependent even if X is
so. If we consider ξ generating exponentially quasi-local statistics, assuming exponential α -mixing of X and using
ideas involving clustering, we prove a central limit theorem (see Theorem 1.13).
Further, even when these two notions have been used to prove central limit theorems, the attention has been restricted
to either random fields on Zd or Rd but rarely on more general spaces. The notions of asymptotic independence and
exponentially quasi-local can be naturally defined on any metric space and so it is obvious to wonder whether such
extensions to more general spaces and their consequences have been considered. While the ergodic theorem (or
strong law in our language) has been studied extensively in the ergodic theory or dynamical systems literature, the
extension of central limit theorem to more general spaces (mainly to ‘nice’ groups or their Cayley graphs) has received some attention but still many questions persist. For example, we refer the reader to the survey of [24], and to
the recent articles [23, Theorems 3.4 and 3.5], [8, Theorems 1.1 and 1.5]. Though the language is vastly different,
ξ
the crux of these results is again that the CLT holds for Hn under suitable weak mixing conditions (which are more
closer to our clustering condition) on the random fields {ξ (x, Pn )}z∈Wn where Wn is now a ball of radius n centred
at the identity element in a group or its Cayley graph. We note here that CLTs have been established on general
metric measure spaces for statistics of Poisson point processes (see [63]), however to the best of our knowledge, we
are not aware of any work under α -mixing or (BL, Θ ) assumptions.
While all the aforementioned results are important precursors to our article, to the best of our knowledge, we are not
aware of a CLT with easy-to-use geometric conditions on ξ , and simple mixing conditions on X . Such conditions
are more common in the literature on point processes in continuous settings such as Euclidean spaces or some nice
compact manifolds (see [91,77,11,63]). In this article, we prove similar generic central limit theorems for ‘nice’
geometric statistics of the form (1.1) for suitably mixing or clustering random fields on Cayley graphs of finitely
generated infinite groups. We make extensive remarks following our main theorems (see Remark 1.14) comparing
our results with those in literature, and also indicating that the “problem of establishing sufficient conditions on a
function defined on an ergodic dynamical system in order that the central limit theorem holds” (see Derriennic [24]) is far from being answered satisfactorily.
Given the above context, we remark that while many of the interesting models (for example, massive Gaussian
free field, off-critical Ising model, determinantal random fields, et al.) satisfy the (BL, Θ )-dependence condition, not
many satisfy α -mixing. In other words, there is a trade-off between (BL, Θ )-dependence and α -mixing in that the
former admits more examples but the latter is powerful enough in proving asymptotics. To our advantage, clustering
condition manages to retain the benefits of both α -mixing and (BL, Θ )-dependence. We summarise our central limit
theorem in relation to those in the literature in the following table.
Linear statistics

Local statistics

Quasi-local statistics

Clustering
[67] *+
Theorem 1.10 *+
Theorem 1.10 *+
(BL, Θ )
[20] *
Theorem 1.10 *+
Theorem 1.10 *+
α -mixing
[19] *
Theorem 1.12 *
Theorem 1.13*
Table 1 Central Limit Theorems : The above table summarises the central limit theorems proved in this paper, and its precursors in the
literature. Here, the + refers to the + condition as in Figure 1 and ∗ refers to assumption on suitable variance lower bounds. ∗+ means
that both ∗ and + conditions hold. We are not explicitly referring to other moment bounds as these can be shown to hold in many of our
examples.

6

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

The direct precursor to this article is [11], where asymptotics for geometric statistics of clustering point processes in
Rd were proven. While the multivariate central limit theorem (Theorem 1.11) and α -mixing central limit theorems
(Theorems 1.12 and 1.13) are new, analogues of our Theorems 1.7 - 1.10 in Rd are known (see [11, Theorems 1.11
- 1.14]). While one expects CLT results for point processes in Rd to hold for Zd , for more general spaces it is apriori
not obvious the inter-relation between growth and amenability of the spaces on one hand and correlation decay
rates of the spin models and quasi-locality of the statistics on the other hand. This paper generalises [11] in this
different direction by giving sufficient conditions for the aforementioned inter-relations at least on discrete spaces.
This is a step in the direction of CLTs of geometric statistics of mixing random fields on more general spaces.
As mentioned in Remark 1.14, our work indicates a way to unify the two frameworks as well as extend similar
generic limit theorems to more general spaces. As with many other papers using clustering condition (see, for
instance, [67,68,5,71,8]), we also derive suitable bounds on mixed moments of the random field {ξ (z, Pn )}z∈Wn
and then use the cumulant method to prove the central limit theorem. The idea of bounding mixed moments via
factorial moment expansion of [9,10] is borrowed from [11]. However, the statements and proofs of theorems
here are simpler exploiting the structure of discrete Cayley graphs, and also lead to applications with minimal
assumptions (see Section 2.2). This also increases the scope of our applications. For example, the application to
k-nearest neighbour graphs in the continuum set-up (see [11, Section 2.3.4]) needs the assumption of determinantal
point processes while our application here (Theorem 2.9) works for any clustering spin model due to Lemma 4.1.
Also, our Lemma 4.1 will be very useful in verifying exponentially quasi-locality for many statistics of clustering
spin models. One distinct advantage of discrete point processes over their continuum versions is that usage of Palm
theory can be avoided completely. This combined with discreteness of the space yields very easy-to-verify moment
conditions compared to the Rd case. For example, local statistics on discrete spaces will be U-statistics and hence
trivially satisfy the moment conditions unlike their continuum counterparts. Comparing our applications in Section
2.2 with those of [11, Section 2.3] will emphasize this point. Further, we are able to compare clustering with mixing
conditions more directly on discrete Cayley graphs as well as furnish central limit theorems for exponentially quasilocal statistics of exponentially α -mixing random fields.
Organisation of the paper : We shall first introduce all the preliminaries – Cayley graphs, random fields, exponentially quasi-local statistics, clustering spin models, mixing coefficients, quasi-association and (BL, Θ )-dependence
– in Sections 1.1 - 1.4. After introducing the necessary notions, we shall state all our results in Section 1.5 with a
detailed discussion of our results in Remark 1.14. In Section 2, we shall primarily discuss examples and applications
of our results on clustering spin models. In addition, we shall also elaborate on the connections between different
notions of ‘asymptotic independence’. Examples of spin models satisfying our assumptions are mentioned in Section 2.1, and applications of our general results to random cubical complexes are described in detail in Section 2.2.
All our proofs are in Section 3 including the crucial factorial moment expansion in Section 3.1. Finally, we conclude
with a void probability bound needed in our applications in the Section 4.
1.1 Preliminary notions and notation
Cayley graphs : We shall briefly define the necessary notions related to Cayley graphs here and we point the reader
to the rich literature [62,65,79] for more details. Let (V, +) be a countably infinite group with a finite symmetric set
of generators S i.e., h ∈ S iff −h ∈ S, where −h denotes the inverse of h in V . By calling S as generators, we assume
that V is the smallest subgroup containing S. Further, denoting the identity by O, we shall assume that O ∈
/ S. The
Cayley graph on (V, +) is the graph G = (V, E) whose vertex set is V and (g1 , g2 ) is an edge if −g1 + g2 ∈ S. In
other words, for all g ∈ G and h ∈ S, (g, g + h) ∈ E. By symmetry of S, it is easy to see that (g, h) ∈ E iff (h, g) ∈ E
and so G is an undirected graph. Since O ∈
/ S and S is the generator of V , G is a simple, connected, regular graph
with vertex degree |S|. For any h ∈ V , the group operation g 7→ g + h is called as translation of G. We emphasise

CLT for statistics of spin models

7

that we have not assumed commutativity of the group, thus g + h need not be same as h + g for all g, h. We shall
now make the blanket assumption that all our Cayley graphs are countably infinite but finitely generated. We shall
call such Cayley graphs as discrete Cayley graphs in the rest of the article.
We shall use d := dG to denote the usual graph distance. Observe that d(x, y) = d(O, −x + y) which we shall also
denote as | − x + y|. For a set A ⊂ V , let |A| denote its cardinality, Wr (y) := {x : d(y, x) ≤ r} denote the ball of radius r
centred at x and Wr := Wr (O). Since the underlying graph is a Cayley graph, Wr (y) ∼
= Wr (x) for all x, y ∈ V where ∼
=
denotes the graph isomorphism. Setting wr := |Wr |, there are two trivial bounds for growth of wn for discrete Cayley
graphs. In particular, there exists d1 ≥ 1 and a constant D > 0 such that for all n ≥ 1,
Dnd1 ≤ wn ≤ |S|(|S| − 1)n.

(1.2)

The lower bound is a result of wn being a strictly increasing N-valued function (See [62, Section VI.A.3]), and the
upper bound follows by comparing with a regular tree of vertex degree |S|. For a subset of vertices, A ⊂ V , we define
/ A such that (x, y) ∈ E}.
its inner vertex boundary as ∂ A := {x ∈ A : there exists y ∈
Definition 1.1 (b-Amenable Cayley graphs) Let G be a discrete Cayley graph as defined above. We say that a
Cayley graph is b-amenable (amenable with respect to the sequence of balls) if
lim

n→∞

|∂ Wn|
→ 0.
|Wn |

We say that a Cayley graph has polynomial growth if for some d2 > 0, we have that |Wn | = O(nd2 )1 .
Polynomial growth is stronger than b-amenability, which in turn is stronger than amenability (see [62, VII.34]).
b-amenability is equivalent to the condition that |Wn |/|Wn−1| → 1. On a discrete Cayley graph, it is easy to note
that |Wn+m | ≤ |Wn ||Wm | and hence by Fekete’s sub-additive lemma limn→∞ |Wn |1/n = wG ∈ [1, ∞) exists. Further, if
wG > 1 (in particular, includes discrete Cayley graphs with exponential growth) then G is not b-amenable (see [62,
VII.34]). We discuss the use of b-amenability assumption in Remark 1.14[(8)].
A large class of examples of discrete Cayley graphs with polynomial growth can be constructed via Gromov’s
famous characterisation of groups with polynomial growth ([40] and see also [79, Theorem 10.1]). We shall quickly
mention this result. For two subgroups A, B, define [A, B] to be the commutator subgroup generated by all a−1 b−1 ab
for a ∈ A, b ∈ B. A group V is said to be nilpotent if the descending series of subgroups (called as lower central series)
defined by Vn+1 = [Vn ,V ] terminates in the trivial subgroup {O}. With this quick definition, Gromov’s theorem can
be stated as follows : A finitely generated group has polynomial growth if and only if it is virtually nilpotent (i.e.,
contains a nilpotent subgroup of finite index). From the definitions, it is easy to see that abelian groups are nilpotent
(since V1 = {O}) and hence any group containing an abelian subgroup of finite index is virtually nilpotent. The class
of abelian groups includes not only our motivating example of integer lattices Zd , d ≥ 1 but other lattices on the
Euclidean plane as well. Further, for a virtually nilpotent group, it was shown in [72, Theorem 51] that there exists
a d ∈ N (called the degree of polynomial growth) such that
lim n−d wn ∈ (0, ∞).

n→∞

(1.3)

Without defining the group, we mention that discrete Heisenberg groups of all dimensions (i.e., Hn (Z) over the ring
Z, see [79, Section 4.1]) is also an example of a nilpotent group and its degree of polynomial growth is 4. For more
examples refer to the previous references, especially [80, Chapter 3], [62, Chapters VI and VII].
1 Here,

we have used the Big O notation of Bachmann–Landau.

8

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

1.2 Random fields and score functions
Let X := {Xx }x∈V , be an S -valued random field, where S is a countable set. We say that the random field X is
d

stationary if {Xx }x∈V = {Xx+y }x∈V for all y ∈ V . The configuration space of such a random field is Ω := S V , which
is equipped with its canonical σ -algebra.

We shall now introduce score functions which are defined on the configuration space. Specifically, a score function
ξ is a function ξ : V × Ω → R, which is measurable with respect to product σ -algebra on V × Ω . We shall assume
all our score functions ξ to be translation-invariant i.e., for all y ∈ V , ξ (· + y, ϑy (·)) = ξ (·, ·), where ϑy : Ω → Ω is
translation by y in the configuration space defined by (ϑy ω )x = ωx+y for all x ∈ V . Next, for ω ∈ Ω , set ω ∩ S :=
{ωx : x ∈ S} for S ⊂ V which is another random field but indexed by S. Also, for any r ∈ N, we set ωr := ω ∩ Wr .
We say that a score function is local if there exists r > 0 such that for any ω ∈ Ω it holds that

ξ (O, ω ) = ξ0 (O, ωr ),

(1.4)

for a (measurable) score function ξ0 : S Wr → R. In statistical physics, the aforementioned local statistics are referred
to as local observables (see [85, Section 2]) or local function (see [33, Definition 3.11]). Further, due to discreteness
of the underlying space, the total mass ((1.1)) induced by local score functions are nothing but local U-statistics and
this is important in simplifying some of our proofs.
More generally, for a fixed score function ξ , not necessary local, we define its radius of stabilisation Rξ (x, ω ) as
follows: for any fixed ω ∈ Ω and r ∈ N, define Ωr (ω ) = {ω ′ ∈ Ω : ωr′ = ωr }, then
R(O, ω ) := Rξ (O, ω ) = inf{r > 0 : ξ (O, ω ) = ξ (O, ω ′ ), ∀ ω ′ ∈ Ωr (ω )},

(1.5)

and R(x, ω ) = R(O, ϑ−x (w)). Measurability of R(O, .) follows because for r ∈ N, we have that
{w : R(O, w) ≤ r} = ∪w∈Ω Ωrξ (w),
ξ

ξ

ξ

/ Now, since Ωr (w) is measurable
where Ωr (w) = Ωr (w) if ξ (O, .) : Ωr (w) → R is a constant and else Ωr (w) = 0.
and there are at most countably many distinct Ωr (w)’s, measurability of R(O, .) follows.
We shall adopt the convention that W∞ = V and Xn = X ∩Wn . Trivially, X∞ = X .
Definition 1.2 (Stabilising score function) We say that ξ is a stabilising score function on X if there is a constant
A > 0 and a function ϕ (t) with ϕ (t) ↓ 0 as t → ∞, such that
sup sup P(R(x, Xn ) ≥ t) ≤ Aϕ (t).

1≤n≤∞ x∈Wn

We say that ξ is a exponentially quasi-local score function on X if ξ is a stabilising score function as in Definition
1.2 and
log ϕ (t)
lim sup
< 0,
(1.6)
tc
t→∞
for some c ∈ (0, ∞). If ξ is a local score function, then R(x, X ) ≤ r a.s., i.e., ϕ (t) = 0 for all t > r and hence it is
trivially exponentially quasi-local. Our definition 1.2 is a quantified version of of [33, Lemma 6.21], where there is
no assumption on the rate of decay. Since such functions are termed as quasi-local in [33, Lemma 6.21], we feel it
is apt to call our functions as exponentially quasi-local. However, we wish to point out that in stochastic geometry
literature, such score functions are also termed as exponentially stabilizing.

CLT for statistics of spin models

9

As is to be expected, we shall need a suitable moment condition as well on the pair (ξ , X ). We say that (ξ , X )
satisfies the p-moment condition if
sup sup E(max{|ξ (x, Xn )|, 1} p ) ≤ M p < ∞,

(1.7)

1≤n≤∞ x∈Wn

where we assume without loss of generality that M p is non-decreasing in p. Though the finiteness above does not
depend upon taking maximum, the definition of M p as above will be notationally convenient in our proofs.
Spin models: A specific class of random fields when S = {0, 1} are called spin models. Let us denote N = {0, 1}V
as the space of spin configurations2 on V . Alternatively, we can think of µ = {µx }x∈V ∈ N as the space of simple
point measures in V by setting µ (·) = ∑x∈A µx . We can also identify µ with its support {x ∈ V : µx = 1}. We shall
use either the measure-theoretic notation µ (·) or the set-theoretic notation µ ⊂ V .
Definition 1.3 (Clustering spin models) By a spin model or point process P on V , we refer to a N -valued
d
random variable. Further, let P be stationary (i.e., P + x = P for all x ∈ V ) and that P is non-degenerate (i.e.,
P(O ∈ P) ∈ (0, 1)). We say that a non-degenerate and stationary P is a clustering spin model, if for all p, q ≥ 1
there exists constants Cp+q , c p+q and a fast decreasing function φ (i.e., φ ≤ 1, decreasing and lims→∞ sm φ (s) = 0 for
all m ≥ 1) such that for all distinct x1 , . . . , x p+q ∈ V , we have that
|P({x1 , . . . , x p+q } ⊂ P) − P({x1 , . . . , x p } ⊂ P)P({x p+1 , . . . , x p+q } ⊂ P)| ≤ Cp+q φ (c p+q s).

(1.8)

where s = d({x1 , . . . , x p }, {x p+1 , . . . , x p+q }) := min1≤i≤p,1≤ j≤q | − xi + x p+ j |. Without loss of generality, we assume
that Ck is non-decreasing in k and ck is non-increasing in k.
We shall say that φ is summable in G if for all c > 0, p ≥ 1,

∑ φ 1/p (c|x|) = ∑ |∂ Wn|φ 1/p (cn) < ∞.

x∈V

(1.9)

n≥0

Trivially, compactly supported φ are summable on any discrete Cayley graph G and fast decreasing φ are summable
on any discrete Cayley graph with polynomial growth. If lim infn→∞ n−1−ε log φ (n) < 0 for some ε > 0, then φ is
clearly summable on any discrete Cayley graph because of the upper bound on wn in (1.2).
Remark 1.4 Though seemingly simple, clustering has various implications for spin models (see Lemma 4.1, for
instance). Another important consequence is that if P clusters, so does P c := V \ P, where we have identified P
with supp(P) (see (4.2) for a proof). If Ck , ck are the clustering constants and φ the clustering function for P, then
2kCk , ck are the clustering constants and φ the clustering function for P c .
Next, setting the joint density or k-correlation functions of P as

ρ (k) (x1 , . . . , xk ) = P({x1 , . . . , xk } ⊂ P),

where x1 , . . . , xk are assumed to be distinct and else ρ (k) (x1 , . . . , xk ) = 0.

A spin model P is said to be exponentially clustering if P is a clustering spin model as in Definition 1.3 with
ck ≥ c > 0, and Ck = O(kak ) for some a ∈ [0, 1), and the clustering function satisfies the growth condition3
lim sup
t→∞

log φ (t)
< 0,
tb

(1.10)

2 Usually, lattice spin configurations are defined as elements of {−1, +1}V but this is trivially equivalent to our definition. Further,
these spin-configurations are also referred more specifically as two-state spin configurations to emphasise that spins here can take only
two values instead of multiple values.
3 We have referred stretched exponential or super-exponential also as exponential for convenience.

10

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

for some b ∈ (0, ∞).
Since spin models are also random fields, we can define score functions as we did for random fields and all the related
definitions carry forward for spin models as well with a few changes. For a spin model, we assume ξ (x, µ ) = 0
if x ∈
/ µ . We say that ξ satisfies the exponential growth condition if for some κ ∈ [0, ∞) and for all t > 0 and
µ ∈ N (with µt = µ ∩Wt ), we have that
κ
|ξ (O, µt )| ≤ ec∗t .
(1.11)
Remark 1.5 Local score functions (say with R(·, P) ≤ r a.s.) satisfy the growth condition (1.11) and p-moment
condition (1.7). To see this, define
kξt k∞ := sup |ξ (O, µ ∩Wt )| =
µ ∈N

sup
µ ∈N (Wt )

|ξ (O, µ )|,

(1.12)

where again µ ∩ Wt denotes the restriction of the spin model to the ball Wt and N (Wt ) denotes the space of spin
configurations on Wt . Since N (Wt ) is finite for all t > 0, we trivially have that kξt k∞ < ∞ for all t > 0. Thus by
‘locality’ of ξ , we have that |ξ (O, µ )| ≤ kξr k∞ for any µ ∈ N and implying the p-moment condition and the
exponential growth condition trivially.
Remark 1.6 Let ξ be a exponentially quasi-local statistic as in (1.6) and satisfy polynomial growth condition i.e.,
there exists κ ∈ [0, ∞) such that for all µ ∈ N and t > 0, we have that
|ξ (O, µ ∩Wt )| ≤ c∗t κ .

(1.13)

Then, using the exponential tail decay of the radius of stabilisation along with the polynomial growth condition,
we have that for all p > 1,
E(max{|ξ (x, Pn )| p , 1}) ≤ E(max{c∗p Rξ (x, Pn ) pκ , 1}) < ∞.
Thus, we have shown that ξ satisfies the p-moment condition (1.7) for all p > 1 if it is exponentially quasi-local and
satisfies the polynomial growth condition. As we shall see, our examples of exponentially quasi-local statistics shall
satisfy this weaker condition of polynomial growth and thereby removing the need to prove the p-moment condition
as well.

1.3 Mixing random fields :
While the notion of clustering suffices for capturing ‘asymptotic independence’ of spin models, we shall need
stronger notions to do so for general random fields. Let X = {Xx : x ∈ V } be a stationary random field as in Section
1.2, and let σS := σS (X ) = σ (Xs , s ∈ S) for a subset S ⊂ V . We shall often omit the reference to X , when the
underlying random field is clear. We define the following two mixing coefficients for subsets S, T ⊂ V

α (S, T ) :=

sup
A∈σS ,B∈σT

|P(A ∩ B) − P(A)P(B)|,

ρ (S, T ) :=

sup
X∈L2 (S),Y ∈L2 (T )

|Corr (X,Y ) |,

where L2 (S) denotes all square-integrable, σS measurable functions. Further we define

α (s) := sup{α (S, T ) : d(S, T ) ≥ s} ; ρ (s) := sup{ρ (S, T ) : d(S, T ) ≥ s}.

CLT for statistics of spin models

11

A random field X is said to be α -mixing (resp. ρ -mixing) if α (s) → 0 (resp. ρ (s) → 0) as s → ∞. Using the above
notation and by choosing appropriate indicator functions, one can trivially obtain the following relationship for any
fixed s ∈ (0, ∞),
1
α (s) ≤ ρ (s).
4

(1.14)

In light of (1.14), it is clear that ρ -mixing implies α -mixing. However, it is well known that the two mixing properties
are equivalent for any dimension. In fact for d ≥ 2, one can show that

ρ (s) ≤ 2π α (s).

(1.15)

We refer the reader to [16, Theorem 1] for the complete statement and proof when the underlying parameter space
for the random field is assumed to be either Zd , or Rd . However, it is not difficult to observe that the arguments set
forth in proving (1.15) can be used verbatim to conclude the same result even when the underlying parameter space
is assumed to be a discrete Cayley graph as in this paper, with the understanding that generators play the same role
as eigen bases of Zd or Rd .
Standard examples of mixing random fields include stationary Gaussian random fields on Zd with spectral density
bounded away from zero (see Section 2.1 [25, Theorem 2]). In fact, one can obtain precise decay rates of the
mixing coefficients of stationary Gaussian random field on Zd if in addition to the condition on the spectral density
stated above, the covariance of the random field decays appropriately as a function of the distance between random
variables (see Section 2.1 [25, Corollary 2]). We refer the reader to [18] for related discussion and many more
examples of processes satisfying various mixing conditions.

1.4 Quasi-associative and (BL, Θ )
As stated earlier, we shall now introduce (BL, Θ )-dependence and quasi-association to provide a complete account
of our results for various dependence structures. A square integrable random field {Xx : x ∈ V } is quasi associated
if for any disjoint sets U,W ⊂ V , and any bounded Lipschitz functions f ∈ R|U| and g ∈ R|W | , we have
|Cov ( f (XU ), g(XW )) | ≤ L f Lg

∑

x∈U;y∈W

|Cov (Xx , Xy ) |,

where L f and Lg are Lipschitz constants of f and g, respectively. This definition is further refined in [20] by restricting the class of test functions to bounded Lipschitz functions on appropriate domain, and employing the Bernstein
blocks technique to introduce a sequence Θ = {θr }r∈Z+ of nonnegative numbers satisfying θr → 0 as r → ∞. For a
given Θ , a random field {Xx , x ∈ V } is said to be (BL, Θ ) dependent if for any disjoint sets I, J ⊂ V , and f ∈ BL(R|I| ),
g ∈ BL(R|J| ), with r = d(I, J), we have
Cov ( f (XI ), g(XJ )) ≤ L f Lg min (|I|, |J|) θr .

(1.16)

A stationary quasi-associated random field {Xx , x ∈ V } can be tested to be (BL, Θ ) by checking the validity of (1.16)
for θr = ∑|x|≥r |Cov (X0 , Xx )|. Please refer to [20] for more details.

12

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

1.5 Our results
As mentioned in the introduction, our main statistic of interest is the sum of score functions of a spin model, i.e.,
our interest lies in the asymptotics of
Hnξ = H ξ (Pn ) :=

∑

ξ (x, Pn ),

x∈Pn

ξ

as n → ∞. We shall now state our abstract theorems for Hn and mention examples of scores and spin models in the
ξ
next section. For distinct x1 , . . . , x p and k1 , . . . , k p ≥ 1, define the mixed moments of Hn (·) for 1 ≤ n ≤ ∞ as
k ,...,k p

mn1

p

(x1 , . . . , x p ) := E(∏ ξ (xi , Pn )ki ).

(1.17)

i=1

The above quantity is well-defined if ξ ’s are non-negative or if we assume max{k1 , . . . , k p }th moment of ξ (x, Pn )
exists for all x ∈ Wn . In our results the later shall hold and hence we shall not have to bother ourselves about existence
k ,...,k
of mn1 p (x1 , . . . , x p ).
These play a crucial role in the proof analogous to that of moments for a usual random variable. The following
theorem is the key step in our proofs of the weak law of large numbers and the central limit theorem.
Theorem 1.7 (Clustering of mixed moments) Let G be a discrete Cayley graph which together with (ξ , P) satisfy
one of the following two conditions :
1. G is b-amenable; P is a clustering spin model and ξ is a local score function as in (1.4).
2. G has polynomial growth; P is an exponentially clustering spin model as in (1.10), ξ is a exponentially quasilocal score function as in (1.6) satisfying the exponential growth condition (1.11) and the p-moment condition
(1.7) for all p ≥ 1.

Then, the mixed moments of {ξ (x, Pn )}x∈Pn for 1 ≤ n ≤ ∞ satisfies clustering, i.e., there exists constants C̃K , c̃K
p+q
such that for all z1 , . . . , z p+q with s = d({z1 , . . . , z p }, {z p+1 , . . . , z p+q }) and K = ∑i=1
ki , ki ≥ 1 for i = 1, . . . , p + q,
we have that
k ,...,k p+q

|mn1

k ,...,k p

(z1 , . . . , z p+q ) − mn1

k

(z1 , . . . , z p )mnp+1

,...,k p+q

(z p+1 , . . . , z p+q )| ≤ C̃K φ̃ (c̃K s),

(1.18)

for a fast decreasing function φ̃ . Further, under Assumption 1 if φ is summable as in (1.9), so is φ̃ .
We shall first state the weak law of large numbers and then the central limit theorem, for which we recall that
wn = |Wn |, where Wn is a ball of radius n in G, defined earlier in Section 1.1.
Theorem 1.8 (Weak law of large numbers for exponentially quasi-local statistics of clustering spin models)
Let G be a discrete Cayley graph which together with (ξ , P) satisfy one of the following two conditions :
1. G is b-amenable; P is a clustering spin model with summable φ as in (1.9) and ξ is a local score function as in
(1.4).
2. G has polynomial growth; P is an exponentially clustering spin model as in (1.10), ξ is a exponentially quasilocal score function as in (1.6) satisfying the exponential growth condition (1.11) and the p-moment condition
(1.7) for some p > 2.
Then as n → ∞,

ξ
2
w−1
n Var(Hn ) → σ (ξ , P) :=

∑ Cov (ξ (O, P), ξ (z, P)) ∈ [0, ∞),

z∈V

and further
P

ξ
w−1
n Hn → E(ξ (O, P)) ∈ R.

CLT for statistics of spin models

13

The following abstract central limit theorem is an extension of [67, Theorem 1] to general random fields (See also
[67, Remark 1]). Since we are unable to find a general statement of the form below, we shall sketch the proof of this
later.
Theorem 1.9 (CLT for clustering random fields.) Let Xn := {Xn,x }x∈Wn , n ≥ 1 be a sequence of random fields such
that supn≥1 supx∈Wn E(|Xn,x | p ) < ∞ for all p ≥ 1. Further, we assume that Xn satisfy
clustering of mixed moments i.e., there exists constants CKX , cXK such that for all z1 , . . . , z p+q with
p+q
s = d({z1 , . . . , z p }, {z p+1 , . . . , z p+q }) and K = ∑i=1
ki , ki ≥ 1 for i = 1, . . . , p + q, we have that
p+q

p

q

ki
ki
ki
|E( ∏ Xn,z
) − E(∏ Xn,z
)E(∏ Xn,z
)| ≤ CKX φX (cXK s),
i
i
i
i=1

i=1

(1.19)

i=1

where φX is a fast-decreasing function and satisfies the summability condition as in (1.9). Set Hn := ∑x∈Wn Xn,x .
Further, if for some ν > 0, it holds that
Var(Hn ) = Ω (wn ν ),
then we have that as n → ∞,

Hn − E(Hn ) d
p
⇒ N(0, 1).
Var(Hn )

Combining Theorems 1.7 and 1.9, we obtain easily the following result that is very convenient for use in applications
as we shall see in Section 2.2.
Theorem 1.10 (CLT for exponentially quasi-local statistics of clustering spin models) Let G be a discrete Cayley graph and together with (ξ , P) satisfy one of the following two conditions :
1. G is b-amenable; P is a clustering spin model with summable φ as in (1.9) and ξ is a local score function as in
(1.4).
2. G has polynomial growth; P is an exponentially clustering spin model as in (1.10), ξ is a exponentially quasilocal score function as in (1.6) satisfying the exponential growth condition (1.11) and the p-moment condition
(1.7) for all p ≥ 1.
Further if for some ν > 0 it holds that
then we have that as n → ∞,

Var(Hnξ ) = Ω (wn ν ),
 
ξ
ξ
Hn − E Hn d
q
⇒ N(0, 1).
ξ
Var(Hn )

ξ

ξ

ξ

In applications, one is also interested in a joint distributional limit of the vector (Hn 1 , . . . , Hn k ) where Hn i is the total
mass corresponding to the score function ξi . We shall combine Theorems 1.8 and 1.10 along with the Cramér-Wold
theorem to derive the following multivariate central limit theorem.
Theorem 1.11 (Multivariate CLT for exponentially quasi-local statistics of clustering spin models) Let G be a
discrete Cayley graph and together with (ξi , P), i = 1, . . . , k satisfy one of the following two conditions :
1. G is b-amenable; P is a clustering spin model with summable φ as in (1.9) and ξi , i = 1, . . . , k are all local
score functions as in (1.4).
2. G has polynomial growth; P is an exponentially clustering spin model as in (1.10), ξi , i = 1, . . . , k are all
exponentially quasi-local score functions as in (1.6) satisfying the exponential growth condition (1.11) and the
p-moment condition (1.7) for all p ≥ 1.

14

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran
ξ

ξ

Set H̄n := (Hn 1 , . . . , Hn k ). We have that as n → ∞,
H̄n − E(H̄n ) d
⇒ N(0, Σ ),
√
wn
where Σ := (Σ (i, j))1≤i, j≤k is the covariance matrix defined by

Σ (i, j) :=

∑ Cov (ξi(O, P), ξ j (z, P)) ∈ R.

z∈V

Though there is no variance lower bound assumption in our multivariate CLT, it is implicit in the result. Note that
the limiting Gaussian vector is non-degenerate iff σ 2 (ξi , P) > 0 for some i ∈ {1, . . . , k}.
We shall now move to next class of results which are based on different set of assumptions involving mixing coefficients. Broadly, the results are the same as stated above in the case of clustering random fields, but with little leeway
allowing more general random fields X . The interest again lies in the asymptotic distribution of appropriately
normalised and scaled sums of
Hnξ = H ξ (Xn ) := ∑ ξ (x, Xn ).
x∈Xn

We note here that when the context is clear, we often omit reference to the underlying random field.
Theorem 1.12 (CLT for local statistics of α -mixing random fields) Let G be a discrete Cayley graph with finite
symmetric generator, and X be stationary, α -mixing random field defined on G. Then, writing ξ for a local statistic
as defined in (1.4), we have as n → ∞,
 
ξ
ξ
Hn − E Hn d
⇒ N(0, σ 2 ),
√
wn
where σ 2 = ∑z∈G Cov (ξ (O, X ), ξ (z, X )) .
Next, we shall state the analogous limit theorem stated for the sums of exponentially quasi-local score functions
evaluated on random field satisfying certain assumptions on the rate of decay of the mixing coefficients.
Theorem 1.13 (CLT for exponentially quasi-local statistics of exponential α -mixing random fields) Let G be
a discrete Cayley graph and together with (ξ , X ) such that G has polynomial growth as in Definition 1.1, X is an
exponential α -mixing spin model (i.e., lim sups→∞ s− log α (s) < 0), ξ is a exponentially quasi-local score function
as in (1.6) and the p-moment condition (1.7) for all p ≥ 1. Additionally, let us assume that
Var(Hnξ ) = Ω (wn ν ),
for some ν > 0. Then as n → ∞,

 
ξ
ξ
Hn − E Hn d
q
⇒ N(0, 1).
ξ
Var(Hn )

A multivariate version of the above theorem can also be concluded using precisely the same set of arguments as put
forth in the proof of Theorem 1.11.
Remark 1.14 (Remarks on our results and future directions)

CLT for statistics of spin models

15

1. Before comparing our results specifically with CLTs in the literature on mixing, (BL, Θ )-dependence or ergodic theory, we wish to comment on the general points of likeness and unlikeness between our results and
those. Firstly, our CLT does not require volume-order variance growth unlike the CLTs available in the aforementioned literature. Secondly, as mentioned in the introduction, we provide what we believe as easy-to-use
geometric conditions on ξ and mixing/clustering conditions on P for CLTs to hold. All the CLTs, including
ours, require a non-trivial variance lower bound assumption for the limit to be non-degenerate.
2. The assumption of S being countable is mainly required to ensure measurability of the radius of stabilization
in Theorem 1.13 but the proof of theorem 1.12 will work even with S being a Polish space.
3. Comparison with CLTs under mixing conditions : CLTs under α -mixing condition are known for linear statistics
([74]) and here with suitable additional assumptions, we have extended it to local and exponentially quasi-local
statistics in Theorems 1.12 and 1.13. But for spin models clustering is a simpler condition to check than α mixing as attested by the many spin models (see Section 2) that satisfy clustering condition.
4. Comparison with CLTs under (BL, Θ )-dependence : Again, CLTs under (BL, Θ )-dependence are proven for linear statistics ([20]) and many spin models do satisfy this condition. But as mentioned before, it is far from clear
whether local or exponentially quasi-local statistics of (BL, Θ )-dependent random fields are (BL, Θ )-dependent.
This problem arises mainly due to the specific structure of the covariance decay required in the (BL, Θ ) dependence condition.
5. Comparison with CLTs on general spaces : Though the underlying spaces considered in [24,8] are far more
general than ours, the exponential mixing conditions assumed on the random field {ξ (z, P)}z∈V (see [8, Definition 2.1]) is similar to our clustering of mixed moments as in (3.32) and the condition [8, (1.7)] plays the
role of our summability condition (1.9). Further, for (strictly) exponential clustering (i.e., b ≥ 1 in (1.10)) and
Cayley graphs with sub-exponential growth, our summability condition (1.9) holds and so does [8, (1.7)] (see
[8, Section 3]). However, we also allow for sub-exponential clustering (see (1.10)) as well provided it is suitably
fast-decreasing depending on the score function ξ and the growth of the Cayley graph.
6. Normal approximation : While our focus has only been on central limit theorems, it is not uncommon to ask
for rates of convergence in central limit theorems. The well-known Stein’s method has often been used to derive such rates. For example, rates of normal convergence for linear statistics (and also some local and global
statistics) for Ising model and some other specific particle systems have been derived recently in [36]. For some
models, the clustering property of spin models play a crucial role. [36] exploits the positive association property
of spin models and hence applies for increasing statistics of spin models but it is not clear if it applies to clustering spin models or non-linear statistics like in our CLTs. Another way to obtain rates of normal convergence is
by obtaining suitable bounds on the growth of cumulants (see [42, Lemma 4.2] and [82]). This method of normal
approximation necessitates a more precise quantification of our cumulant bounds. Further, [42, Lemma 4.2] also
gives cumulant bounds needed for moderate deviations. This would be a worthwhile direction to pursue in the
future.
7. Cumulant Bounds : Another use of cumulant bounds to derive CLT is in [26,30] where such bounds are crucially used in the weighted dependency graph method to prove CLTs. In [26], CLTs for local and some global
statistics of the Ising model (see Section 2.1.2) are proved using bounds on cumulants and a generalisation of
the dependency graph method. Again, it requires restrictions on regimes for the Ising model and variance lower

16

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

bounds but it is not obvious if the methods can be adapted to other similar spin models.
8. Scaling limits : Suppose that V = Zd , the integer lattice with the generators of the group being {z = (z1 , . . . , zd ) ∈
Zd : kzk1 = 1} with k · k1 denoting the ℓ1 distance. The corresponding Cayley graph on Zd is nothing but
V = Zd ; E = {(z1 , z2 ) : kz1 − z2 k1 = 1}.

(1.20)

Then it is possible to consider a suitably scaled version of the random field {ξ (x, Pn )}x∈Wn and study its scaling
limit. Two possible choices for scaling are either to consider the random field
ξ
Xn1 := {(Var(Hn ))−1/2 (ξ (n−1 x, Pn ) − E(ξ (n−1 x, Pn )))}x∈Wn or the random field
Xn2 := {(Var(∑nx≤y≤(n+1)x ξ (x, P)))−1/2 ∑nx≤y≤(n+1)x (ξ (x, P)−E(ξ (x, P)))}x∈Zd where ≤ is the co-ordinate
wise ordering of points in Zd . The former scaling is more in the spirit of the scaling considered in [11, (1.3)]
and the latter is in spirit of the scaling considered in [67, Theorem 2]. In both the cases, the limit is expected to
be a Gaussian random field with a white noise like structure i.e., the covariance matrix of all finite dimensional
marginals converging to a diagonal matrix (see [11, (1.21)] and [67, Theorem 2]).
9. We are not aware of any examples of b-amenable Cayley graphs exhibiting non polynomial growth. However,
our proof methods for b-amenable graphs should allow our results to be proven under the assumption that wn
grows sub-exponentially i.e., lim supn→∞ n−1 log wn = 0. For such groups, there exists a subsequence nk → ∞ as
k → ∞ such that w−1
nk |∂ Wnk | → 0 as k → ∞. Thus, if we take sums over Wnk for the subsequence nk chosen as
above instead of taking sums over Wn in (1.1), one would expect our results to hold under such asymptotics as
well.
10. We have restricted ourselves to a class of amenable Cayley graphs but it is natural to ask whether one can consider a more general class of graphs. Deriving the motivation from various probabilistic studies (see [1,7,65,
79]), two possible classes of graphs that are suitable to such a study are unimodular random graphs and vertex
transitive graphs. To further emphasise the need for such a study, even graphs on stationary point processes as
studied in [11] can be considered as unimodular random graphs (see [4, Section 5]). Thus a study on unimodular
random graphs can unify the framework in this article and that of [11] apart from considerably extending the
scope of applications.
11. Another persistent but unavoidable assumption not only in our CLTs but in various such generic CLTs in the
literature (including those cited here) is the variance lower bound condition. Such lower bounds are usually
shown by ad-hoc methods. Primarily, for sums of stationary sequences, variance lower bounds can be obtained
under conditions involving the spectral density of the random variables (see Theorem 2 in Chapter 1.5 of [25]).
Alternatively, under the stationary and strong mixing condition, together with appropriate summability condition
of the covariance, the necessary and sufficient condition for meaningful variance lower bounds of partial sums,
is that the variance must grow to infinity (see [17, Lemma 1], or [74, Theorem 2.1]). However, for sums of nonstationary sequences of random variables, the condition ρ ′ (X , 1) < 1 provides suitable variance lower bounds
(see [19, Theorem 2.2]). In this context, a very simple and natural question follows: for an α -mixing random
field X satisfying ρ ′ (X , 1) < 1, writing ρ ′ (X ξ , k) as the ρ -mixing coefficient of the field X ξ , is it possible
to conclude ρ ′ (X ξ , 1) < 1, even for a local statistic ξ ?
2 Examples and applications
In this section, we illustrate our main theorem (Theorem 1.10) by providing examples of exponentially quasi-local
statistics and clustering spin models. Though we shall mainly focus on a variety of applications to random cubical

CLT for statistics of spin models

17

complexes, we shall also hint at others. Also, we shall not mention applications of our mixing CLTs (Theorems
1.12 and 1.13) but we hope the discussion in the introduction and our applications for clustering spin models will
convince the reader that such results are feasible as well.
Prior to discussing the examples, we shall detail relationship between clustering and other measures of of association, specifically, mixing and (BL, Θ ), which are going to be used in this section. Observe that the various mixing
coefficients merely extract the dependence structure without delving into the decay rates of the inherent dependence
structure whereas other dependence structures like the clustering and (BL, Θ ) do, depend and, provide certain information regarding the decay rates of the inherent dependence structure. Heuristically, appropriate decay assumptions
on mixing coefficients may establish a link from mixing to clustering, and possibly to (BL, Θ ). In this direction, we
state the following two propositions establishing connection between mixing and clustering conditions and between
(BL, Θ ) dependence structure and the clustering.
Proposition 2.1 Let P be a stationary, (strong) α -mixing spin random field indexed by G, a discrete Cayley graph.
Assume that the mixing coefficient α (s) is a fast decreasing function. Then, such X also satisfies the clustering
condition with Ck = ck = 1 and φ (s) = α (s).
The proof is rather straightforward, thus we leave it to the reader.
Proposition 2.2 Let P be a stationary, spin random field defined on a discrete Cayley graph G, and for any vertex
u ∈ V denote Zu = 1{u ∈ P} Let the covariance function of Z, given by r(u) = Cov (Z0 , Zu ), satisfy the (BL, Θ )
condition stated in (1.16), with θ (k) = ∑kuk≥k |r(u)|, such that θ (k) is a fast decreasing function. Then, the random
field Z is also a clustering random field with Ck = k, ck = 1 and φ (s) = θ (s).
Proof For any regions P, Q ⊂ G such that |P| = p, |Q| = q, and dist(P, Q) = s, define f : {0, 1} p → R and g : {0, 1}q →
p
R as f (t1 , . . . ,t p ) = ∏i=1
1{ti =xi } and g(s1 , . . . , sq ) = ∏qj=1 1{si =x p+i } for some fixed (x1 , . . . , x p+q ). Next, setting P =
(u1 , . . . , u p )
and
Q = (v1 , . . . , vq ),
clearly,
E( f (ZP )) = P(Zu1 = x1 , . . . , Zu p = x p )
and E(g(ZQ )) = P(Zv1 = x p+1 , . . . , Zvq = x p+q ). Therefore,


|P Zu1 , . . . , Zu p , Zv1 , . . . , Zvq = (x1 , . . . , x p+q ) − P(Zu1 = x1 , . . . , Zu p = x p ) P(Zv1 = x p+1 , . . . , Zvq = x p+q )
= |Cov ( f (ZP ), g(ZQ ))| ≤ (p ∧ q) L f Lg

Now the result follows easily as L f = Lg = 1.

∑

kuk≥s

|Cov (Z0 , Zu ) |.

✷

2.1 Examples of clustering spin models
The simplest example of a clustering spin model is one with i.i.d. spins. The clustering property, which captures
asymptotic independence in a strong way, is a natural condition that is expected to hold in statistical physical models
which have weak dependences. Without delving into the details, we shall mention a few illustrative examples in
this part of the section, and specifically restrict our attention to spin models on the lattice Zd , unless mentioned
otherwise. Notice that by the commutative property of Zd , the distance d(x, y) = | − y + x| = |x − y|, which matches
the ℓ1 distance on Zd .
2.1.1 Level sets of Gaussian fields
Let X = {X(x)}x∈V be a stationary Gaussian random field whose covariance kernel is exponentially decaying i.e.,
ω (x, y) = Cov (X(x), X(y)) is such that lim inf|y|→∞ log ω|y|(O,y) < 0. Further for simplicity, assume that ω (x, y) is a

18

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

function of |x − y| alone. The super-level sets at level u of this Gaussian field, defined as Pu := {x ∈ V : X(x) ≥ u},
is a spin model. To show clustering of Pu , we shall use the following total-variation distance bound between Gaussian random vectors from [6].
We recall the definition of total variation distance dTV between two probability measures µ and ν on a sigma-algebra
F is given by
dTV (µ , ν ) = sup | µ (A) − ν (A)|.
A∈F

Theorem 2.3 (Theorem 4.3 in [6])
 Let X =(X1 , X2 ) andY = (Y1 ,Y2 ) be random Gaussian vectors (not necessarily
Σ11 Σ12
Σ11 0
and ΣY =
. Assume that the size of vectors X1 and Y1 be
centred) with covariances ΣX =
T Σ
0 Σ22
Σ12
22
m, while that of X2 and Y2 be n. Let µX and µY be the laws of the vectors whose entries are 1 or 0, whether the
corresponding entries in the vectors X and Y are positive or not, moreover Σ11 and Σ22 has 1 on diagonal. Then
14

8

1

dTV (µX , µY ) ≤ 2 5 (m + n) 5 (max |Σ12 (i, j)|) 5 .
i, j

Remark 2.4 Theorem 2.3 of Beffara and Gayet was originally stated for centred Gaussian vectors in [6]. By using
the trivial bound P(|X| < ε ) ≤ ε for any Gaussian random variable X with arbitrary mean, where ever necessary,
in the proof of Theorem 2.3 it trivially extends to any non-centred Gaussian vector. Hence, Theorem 2.3 is true for
arbitrary level sets of a Gaussian vector.
Thus trivially from Theorem 2.3, and the corresponding definitions, we have the following corollary.
Corollary 2.5 Let X be a stationary Gaussian field, on a discrete Cayley graph G, with X(O) having unit variance.
14 8
The level sets Pu satisfy exponential clustering with clustering constants Ck = 2 5 k 5 , ck = 1 and clustering function
φ (t) = ω (O, x), where |x| = t.
Massive Gaussian free field : Massive Gaussian free field, on the lattice Zd , is defined as the Gaussian field Gm (·)
′
whose covariance kernel is gm (x, y) = am e−am |x−y| for all x, y ∈ Zd where am and a′m are some constants depending on
m (see [33, Proposition 8.30]). The function gm is nothing but the discrete Green’s function for the operator m2 + ∆ ,
where m > 0 is considered as the mass and ∆ is the discrete Laplacian operator on Zd . The super-level set of the
massive Gaussian free field at level u, Pu := {z : Gm (z) ≥ u} is a spin model. Invoking Theorem 2.3 along with
14 8
Corollary 2.5, Pu is an exponentially clustering spin model as in (1.10) with clustering constants Ck = 2 5 k 5 , ck = 1
′
and clustering function φ (t) = am e−amt .
2.1.2 Ising model
In this example, we consider the simplest form of the Ising model with interactions from the nearest neighbours
alone. For a finite set Λ ⊂ Zd , the Ising model is defined as the probability distribution on ΩΛ := {±1}Λ such that
for σ = {σi }i∈Λ ∈ ΩΛ , the probability distribution is given by








1
P(σ ) =
exp −β ∑ σi σ j − h ∑ σi ,


ZΛ

i, j∈Λ
i∈Zd 


|i− j|=1

CLT for statistics of spin models

19

where β > 0 is the inverse temperature, h ∈ R is the external field, ZΛ is a suitable normalising constant, and the
term in the exponential is referred to as interaction potential. By making the transformation σi → σi2+1 , the Ising
model is a spin model in the framework of this article. For any d ≥ 1 and h = 0, there exists βc (d), such that when
0 < β < βc (d) a unique measure on Zd is obtained by suitably taking the weak limit as Λ approaches Zd . For
any d ≥ 1 and h 6= 0, the weak limit exists whenever β > 0. For d = 1, the weak limit exists for any h ∈ R, for
any β ∈ (0, ∞). We refer the reader to [33] (Chapters 3 and 6, Theorem 3.25 and Exercise 6.17) for more on the
characterisation of these limits. In [78] it was established that the Ursell functions decay exponentially, and that the
infinite volume limit of these does not depend on the sequence in which the limit is taken. Therefore, combining
results from [33] and [78] for the Ising model on any dimension d ≥ 1 with 0 < β < βc (d), in the absence of
external magnetic field (i.e., h = 0), and for any β > 0 in the presence of external magnetic field (i.e., h 6= 0), the
Ursell functions decay exponentially. Ursell functions (mixed cumulants) of a random field X = {X(x)}x∈V are
defined as follows : For n ≥ 1, x1 , . . . , xn ∈ V,
un (x1 , . . . , xn ) =

∂
∂
...
log Eez1 X(x1 )+···+zn X(xn ) |z1 =···=zn =0 .
∂ z1
∂ zn

Joint moments can be expressed using the relation,
E(Xx1 Xx2 . . . Xxn ) = ∑ ∏ u|B| (Xxi : i ∈ B),
π B∈π

where the sum is over all partitions π of {1, 2, . . . , n} and the product is over all blocks B in the partition π .
P(Xx1 = Xx2 = · · · = Xx p+q = 1)−P(Xx1 = · · · = Xx p = 1)P(Xx p+1 = · · · = Xx p+q = 1) = ∑ ∏ u|B′ | (Xxi : i ∈ B′ ), (2.1)
π ′ B′

where the sum is over the partitions π ′ of {1, 2, . . . , p + q}, whose blocks B′ has at-least one term from each of the
sets {1, 2, . . . , p} and {p + 1, . . . , p + q}. From the exponential decay of Ursell functions, each term in the right hand
side of 2.1 is decaying exponentially in the distance between the sets {x1 , . . . , x p } and {x p+1 , . . . , x p+q }.
Therefore in the specified regimes ({β > 0}\{h = 0; β ≥ βc (d)}), the Ising model is an exponentially clustering
spin model. Referring to Definition 1.3, we have Ck = O(k2 ), ck = 1 and the clustering function is the two point
function φk (·) = u2 (·), where u2 (t) = |P(Xx1 = Xx2 = 1) − P(Xx1 = 1)P(Xx2 = 1)| for some |x1 − x2 | = t and is
exponentially decreasing in t.
2.1.3 Determinantal point processes
Determinantal point processes have been of considerable interest in probability and statistical physics literature. We
refer the reader to [49, Chapter 4] for an introduction, [64] for these processes on discrete structures and [66] for
stationary determinantal point processes on Zd .
A stationary determinantal point process P on any discrete Cayley graph G is defined to be a spin model whose
probabilities are determined by the relations P({x1 , . . . , xk } ⊂ P) = det(K(xi , x j ))i, j , where K is a suitable nonnegative definite real-valued kernel. From these relations it follows that the kernel K has to be a contraction and
invariant under group action. Consider a stationary determinantal point process on the graph G whose kernel K
decays exponentially with graph distance. It is known that (see [15, Theorem 3.4]) determinantal point processes
are negatively associated, which implies that they are quasi-associated (for related discussion see [22, Definitions
1.1, 1.2 and 1.3]). Hence, from the definition of quasi-association (see Section 1.4) and by using arguments as in
the proof of Proposition 2.2, we have that the stationary determinantal point process is exponentially clustering with
Ck = O(k2 ), ck = 1 and the clustering function φk (t) = |K(O, z)| for some |z| = t.

20

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

2.1.4 Other possible clustering spin models:
The phenomenon of exponential decay of covariances (two point functions) of the fields in distance is sometimes
termed as ‘massive’, and it is a common feature in many statistical physics models. To mention a few, δ -pinned
models for Gaussian fields studied in [50,89,14], stochastic interface models considered in [34,89], Gross-Neveu
model considered in [58] exhibit this behaviour in suitable regimes. These are only indicative references for the
literature and by no means exhaustive. One can define spin models using these fields, for e.g. choosing a level set or
determining the spins by value of the field in a finite neighbourhood. It is natural to expect that the corresponding
spin models arising from these fields satisfy exponential clustering conditions.
We mentioned in Section 2.1.2 that the Ising model satisfies clustering condition for high temperatures i.e, β small
enough. One generalisation of Ising model are Gibbs spin models by taking more general interaction potentials.
Various Gibbs spin models with finite-range (or suitably local) potentials are also expected to satisfy clustering
condition for high temperatures (see [48,41,61]).

2.2 Applications to random cubical complexes.
We shall illustrate our central limit theorem for local and exponentially quasi-local statistics of clustering spin
models using random cubical complexes. We shall now re-introduce cubical complexes a little more formally before detailing our results about the same. For more details on cubical complexes, see [52,90,47] or [56, Section 6.4].
Set Q = W1/2 = [− 21 , 12 ]d , the unit cube centred at origin and Qx = Q+x, the shifted cube for any x ∈ Rd . Given a spin
model µ , define C(µ ) := ∪x∈µ Qx . C(µ ) can be viewed as a random subset of Rd as is done in various applications
in image analysis, stereology and mathematical morphology. Alternatively, define Fk (µ ) = {(x0 , . . . , xk ) ∈ µ (k) :
∩ki=0 Qxi 6= 0}
/ as the k-faces of the cubical complex. By default, set F0 (µ ) = µ ⊂ Zd . Note that Fk (µ ) = 0/ for
2d−1
k ≥ 2d. We shall denote elements of Fk as [x0 , . . . , xk ]. The cubical complex is defined as K (µ ) := ∪k=0
Fk (µ ).
The k-skeleton of K is defined as K k (µ ) := ∪kj=0 Fk (µ ). Trivially, note that the 1-skeleton is the graph with vertex
set µ and edge-set {(x, y) ∈ µ 2 : Qx ∩ Qy 6= 0}.
/ This is same as the Cayley graph on Zd defined in (1.20). We
shall not give more details here on cubical complexes apart and will rather refer the reader to [52]. Given two
cubical complexes K , L , a cubical homomorphism is a map f : K 0 → L 0 such that whenever [x0 , . . . , xk ] ∈ K ,
then [ f (x0 ), . . . , f (xk )] ∈ L . We say that f is a cubical isomorphism if it is a bijection and f −1 is also a cubical
homomorphism.

2.2.1 Local counts and intrinsic volumes:
We first define an abstract class of local statistics that shall include various statistics of interest about cubical complexes and then state asymptotics for this local statistic using our main theorem (Theorem 1.10). Let h : {0, 1}Wk → R
be a function of spin models on Wk . The local score function ξ defined for µ ∈ N , x ∈ µ ,

ξ (x, µ ) := h(µ ∩Wk (x))

(2.2)

Further, set the total mass H h (µ ) = H ξ (µ ) := ∑x∈µ ξ (x, µ ). We shall now give some examples of h and hence of ξ
as well. Let z1 , . . . , zk ∈ Zd be such that Γk := K ({z1 , . . . , zk }) is connected as a set (or equivalently the 1-skeleton

CLT for statistics of spin models

21

is connected). Define the following two h functions on µ ∈ {0, 1}Wk :
hΓk (µ ) =
ĥΓk (µ ) =

1
k!
1
k!

∑

1[K (O, x2 . . . , xk ) ∼
= Γk ]

(2.3)

∑

1[K (O, x2 , . . . , xk ) ∼
= Γk ]1[∪ki=1 (µ ∩W1 (xi )) = {O, x2 , . . . , xk }]

(2.4)

(x2 ,...,xk )∈ µ (k−1)

(x2 ,...,xk

)∈ µ (k−1)

Observe that the the total mass H hΓk (µ ) counts the number of isomorphic copies of Γk in K (µ ) and H ĥΓk (µ ) counts
the number of components in K (µ ) isomorphic to Γk . These are two basic statistics of interest in combinatorial
topology.
Let ∆k := {∪z∈µ Qz : µ ∈ {0, 1}Wk } be the space of all cubical complexes in Wk , and g : ∆k → R be a real valued
functional defined on cubical complexes. Then, setting h(µ ) := g(∪z∈µ Qz ), and by appropriately choosing g, we
can get ξ j (for 0 ≤ j ≤ d) to be a score function such that
H j (µ ) := H ξ j (µ ) = V j (C(µ )),

(2.5)

where V j is the jth intrinsic volume. See [11, Section 2.3.3] or [90, Section 4] for the details of precise definition of
g and ξ j ’s. A precise definition of ξ j can be provided by using the finite additivity of intrinsic volumes. In fact, by
the famed Hadwiger’s characterisation theorem ([84, Theorem 14.4.6]) the additive property on convex ring, motion
invariance and appropriate continuity along with the assumption that V j (rA) = r j V j (A) uniquely characterises the
jth intrinsic volume. V0 is the famed Euler-Poincaré characteristic, and for a convex set A, Vd (A) & Vd−1 (A) are
valuations of A which equal the volume and surface measure of A, respectively. For more details on intrinsic volumes,
we refer the reader to [84].
While first-order asymptotics of the above three statistics can be derived easily using ergodic theory arguments,
results akin to CLTs are proven only in special cases. For example, if the spins are i.i.d., a central limit theorem for
H j ’s (i.e., V j ’s) for j = 0, . . . , d is proven in [90, Theorem 9]. A similar approach can be used to prove central limit
theorems for H Γk and H ξ̂ in the case of i.i.d. spins. As an application of our main theorem, we now reduce the proof
of such theorems for general spin models to that of deriving suitable variance lower bounds.
Theorem 2.6 Let P be a clustering spin model as in Definition 1.3, h a general local function as defined in (2.2)
and ξ the score function induced by h. Set Hn := H ξ (Pn ). Then if for some ν > 0, it holds that
Var(Hn ) = Ω (wνn ),
d

then we have that Var(Hn )−1/2 (Hn − E(Hn )) ⇒ N(0, 1).
We had earlier alluded to the use of integral geometric statistics in morpho-metric analysis of digital images and
the simplest of such statistics being the intrinsic volumes. We shall now show an interesting application of Theorem
1.11 that shall give the joint distribution of the intrinsic volumes of random cubical complexes.
Corollary 2.7 Let P be a clustering spin model as in Definition 1.3 and ξi , i = 0, . . . , d be score functions defined
such that H ξi (µ ) = H j (µ ) (see (2.5)) for any spin model µ . Set H̄n := (H0 (Pn ), . . . , Hd (Pn )). Then we have that
H̄n − E(H̄n ) d
⇒ N(0, Σ ),
√
wn
where Σ := (Σ (i, j))0≤i, j≤d is as defined in Theorem 1.11.

22

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran
h

In general, our multivariate central limit theorem can be applied to (Hnh1 , . . . , Hn k ) for local functions h1 , . . . , hk :
{0, 1}Wr → R. For example, hi = hΓi or ĥΓi where Γi , 1 = 1, . . . , k are k distinct connected complexes. Similar multivariate central limit theorems for subgraph counts have been proven for (Euclidean) Poisson point processes in [75,
Chapter 3]. In the case of i.i.d. spins, one can derive multivariate CLTs for local statistics with rates can be derived
using results in [60]. For example, such a multivariate CLT for intrinsic volumes can be found in [60, Section 4.2].
2.2.2 Nearest neighbour graphs:
We shall now illustrate CLT for exponentially quasi-local score functions via the following graph, which often arises
in computational geometry. It is called the nearest-neighbour graph, and is defined as follows: the vertices are points
of µ and (x, y) is a directed edge (called the nearest neighbour edge (NNE)) if µ ∩ Wr (x) = {x} for all r < |x − y|.
Then the nearest neighbour distance score function is given by

ξNN (x, µ ) =

∑
y∈Zd

|x − y|1[(x, y) an NNE],

(2.6)

where |.| denotes the L1 Euclidean distance. The distribution of {ξNN (x, µ )}x∈µ is important in understanding the
structure of spacing of the points in µ . In this case HNN (µ ) is the total edge length of the ‘weighted’ nearest
neighbour graph on µ i.e., we count edges that occur as a NNE for both the points twice. It is possible to redefine the score function such that HNN (µ ) yields the total edge-length of the nearest neighbour graph but this shall
complicate our analysis a little more. Hence, we avoid it here.
Theorem 2.8 Let P be an exponential clustering spin model as in (1.10), ξ := ξNN (·) be score function of the
nearest neighbour distance as defined in (2.6). Set Hn := H ξ (Pn ). Then if for some ν > 0, it holds that
Var(Hn ) = Ω (wνn ),
d

then we have that Var(Hn )−1/2 (Hn − E(Hn )) ⇒ N(0, 1).
Proof By definition of ξNN , we have that RNN (x, µ ) := Rξ (x, µ ) = inf{r : µ ∩Wr (x) 6= 0}.
/ Thus, by Lemma 4.1, we
obtain that for some A, a′ , ν > 0,
′ ν

/ ≤ Ae−a r .
P(RNN (x, µ ) > r) = P(Pn ∩Wr (x) = 0)
This proves exponentially quasi-locality of the score function as in (1.6) with c = ν . As for the exponential growth
condition (1.11), note that
|ξNN (x, µ ∩Wt (x))| ≤ |∂ Wt (x)|t ≤ Ct d ,
i.e., we have the polynomial growth condition and thus by Remark 1.6, ξNN satisfies the p-moment condition (1.7)
for all p > 1 as well. Thus we have verified all the conditions in (2) of Theorem 1.10 and hence the central limit
theorem follows.
✷
As mentioned in the introduction, one can also study k-nearest neighbour edges i.e., (x, y) form a k-NNE if y is
among the first k nearest neighbours of x i.e., µ ∩ Wr (x) ≤ k for all r < |x − y|. We remark here that in the case of
nearest neighbour functional defined on Zd with i.i.d. spins, one may use [76, Theorem 2.1] to estimate asymptotic
variance. Further, the void probability bound in Lemma 4.1 can be used to show exponentially quasi-locality of
various score functions such as statistics (intrinsic volumes of faces, in-radii of cells) of Voronoi tessellations ([84,
Section 10.2]) and other proximity graphs.

CLT for statistics of spin models

23

2.2.3 Topology of random cubical complexes:
We shall now introduce a slightly more complicated score function from algebraic topology. With increased interest in applied topology, there have been many studies on topological invariants, like Betti numbers, of random
complexes (see [12,53]). The non-trivial Betti numbers of a “nice” set A ⊂ Rd (or a cubical complex K ) denoted
by β0 (A), . . . , βd−1 (A) are a measure of connectivity of the set A. For example, β0 counts the number of connected
components and βd−1 (A) counts the number of connected components of Ac minus one. The other Betti numbers are
harder to describe and in simple terms, βk counts the number of ‘distinct’ (k + 1)-dimensional holes enclosed in the
set A. A little more formally, one constructs a group called the ‘kth homology group’ with non-trivial k-dimensional
cycles or equivalently (k + 1)-dimensional holes and we define βk as the rank of the kth homology group. We assume
that all the homology groups here are vector spaces (over an arbitrary field) and hence the ranks are well-defined. For
details, we refer the reader to [47, Section 2.1], [56, Section 6.4], or [28,52]. Topological properties of random cubical complexes and slightly more general models were investigated recently in [47]. While first-order asymptotics
([47, Section 2.4]) were proven for general ergodic spin models, a central limit theorem ([47, Theorem 2.11 and
2.12]) again required the assumption of i.i.d. spins. We shall now state a central limit theorem for Betti numbers of
random cubical complexes on subcritical spin models. Let C(x, µ ) denote the component (i.e., maximal connected
set or sub-complex) containing x in C(µ ) or K (µ ). We shall use C(x, µ ) to denote a subset of Zd as well as the
corresponding union of unit cubes. We shall say that P is a subcritical spin model (in terms of percolation) if for
all x, y ∈ Zd ,
∗
P(y ∈ C(x, P)) ≤ C∗ e−c |x−y| ,
for some positive constants C∗ , c∗ .

Theorem 2.9 Let P be an exponential clustering spin model as in (1.10) and further, let P or P c be a subcritical
spin model. Then if, Var(βk (C(Pn ))) = Ω (wνn ), for some ν > 0, then
d

Var(βk (C(Pn )))−1/2 (βk (C(Pn )) − E(βk (C(Pn )))) ⇒ N(0, 1).
Further, we have that
d

(wn )−1/2 (β0 (C(Pn )) − E(β0 (C(Pn ))) , . . . , βd−1 (C(Pn )) − E(βd−1 (C(Pn )))) ⇒ N(0, Σ ),
where Σ is as defined in Theorem 1.11 with ξk as defined below in (2.7).
Proof We shall give the proof for the case when P is subcritical and then by duality arguments in algebraic topology
extend it to the case when P c is subcritical. We shall now define an appropriate score function and then prove
exponentially quasi-locality of the same. For a fixed k ∈ {0, . . . , d − 1}, define the score function ξk as

ξk (x, µ ) :=

βk (C(x, µ ))
,
|C(x, µ )|

(2.7)

where |C(x, µ )| stands for the (vertex) cardinality of the connected component of x. Since all the cubes are of unit
volume and disjoint, |C(x, µ )| = Vd (C(x, µ )). By the property of Betti numbers of cubical complexes βk (C(x, µ )) ≤
|C(x, µ )| (see [47, Lemma 3.1]), we have that |ξk (x, µ )| ≤ 1 and hence ξk satisfies the power-growth condition (1.11)
and p-moment condition (1.7) for all p ≥ 1. Further, as βk (A ∪ B) = βk (A) + βk (B) for disjoint A and B, we have

βk (C(µ )) =

∑ ξk (x, µ ).

x∈ µ

Let d(x, µ ) := inf{s : C(x, µ ) ⊂ Ws−1 (x)} be the L1 distance from x to the furthest vertex on the boundary of its
connected component. Suppose r > d(x, µ ) and µ ′ is a spin model such that µ ′ ∩ Wr (x) = µ ∩ Wr (x). Trivially, we

24

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

have that C(x, µ ) = C(x, µ ′ ) and hence ξk (x, µ ) = ξk (x, µ ′ ). In other words, configurations which differ from µ
outside Wr (x) keep C(x, µ ) unchanged. Thus, we derive that R(x, µ ) ≤ d(x, µ ). Now, observe that
P(d(x, Pn ) > t) ≤

∑
y:|y−x|=t

∗

P(y ∈ C(x, P)) ≤ C′t d−1 e−c t ,

and so ξk (x, P) is exponentially quasi-local as in (1.6), implying (ξk , P) satisfies the assumptions in (2) of Theorem 1.10 and hence the CLT follows.
Now suppose that P c is subcritical. By Remark 1.4, we have that P c is also exponentially clustering as in (1.6). Set
µn∗ := Wn+1 \ Wn . By the universal coefficient theorem for simplicial homology ([70, Theorem 45.8]) and Alexander duality ([86, Theorem 16])4 , we derive that β̃k (C(Pn )) = β̃d−k−1 (C((Pn )c )). Further from the homotopy
equivalence of C((µ ∩ Wn ) ∪ Wnc ) and C((µ ∩ Wn ) ∪ µn∗ ) for any spin model µ , we have that β̃d−k−1 (C((Pn )c )) =
β̃d−k−1 (C((P c )n ∪ µn∗ )), where β̃k (·) := βk (·) − 1[k = 0], and we point out that (Pn )c 6= (P c )n .
Now, we can again define appropriate score functions ξk′ as above for β̃d−k−1 (C((P c )n ∪ µn )) and by subcriticality of P c show that (ξk′ , P c ) satisfy the assumption (2) of Theorem 1.10. Thus, the central limit theorem follows
even in the case when P c is subcritical. The multivariate central limit theorem shall follow as we have shown that
✷
(ξk , P), (ξk′ , P c ), k = 0, . . . , d − 1 satisfy the assumptions in (2) of Theorem 1.11.
As we noted earlier, the random cubical complexes considered in [47] are more general and in particular, allow for
lower-dimensional cubical subsets of Q = [−1/2, 1/2]d to be present in the cubical complex without the entire cube
Q being present. To incorporate such a model into our set-up, one can consider a more general spin-model wherein
the underlying lattice L is generated by integer-valued translations of {(a1 , . . . , ad ) : ai ∈ {−1/2, 0, 1/2}} which
are nothing but vertices of the cube Q or mid-points of lower-dimensional cubes in Q. A spin being positive at a
site in L is equivalent to the particular lower-dimensional cube being present where the vertices are nothing but
0-dimensional cubes.

3 Proofs
3.1 Factorial Moment Expansion and clustering of random measures
ξ

The key tool in our proofs is the clustering of Hn derived using the factorial moment expansion introduced in [9,
10]. We shall now first recall the factorial moment expansion and then prove the clustering result. Recall that V is
the vertex set of the Cayley graph G, and N is the space of all spin configurations, or point measures, on V . Let ≺
be a total order on V such that if |u| < |v| then u ≺ v. One can construct such orders by starting with an arbitrary
order on the generators S and then constructing the corresponding lexicographic or co-lexicographic order on V .
For µ ∈ N on V and x ∈ V , set µ|x (·) := µ (· ∩ {y ∈ V : y ≺ x}). By local-finiteness of µ and the property of the
ordering, µ|x (V ) ≤ µ (W|x| ) and so µ|x is a finite measure for all x ∈ V . We denote the null measure by o i.e., o(B) = 0
for all subsets B ⊂ V . For a measurable function ψ : N → R, non-negative integer l and x1 , . . . , xl ∈ V , we define
the factorial moment expansion (FME) kernels as follows : for l = 0, D0 ψ (µ ) := ψ (o). For l ≥ 1,
Dlx1 ,...,xl ψ (µ ) :=

∑ (−1)l−|J|ψ (µ|x∗ + ∑ δx j ),

J⊂[l]

(3.1)

j∈J

4
Alexander duality shows isomorphism of kth homology group with (d − k − 1)th cohomology group and the universal coefficient
theorem gives equality of the rank of (d − k − 1)th cohomology group with the corresponding homology group.

CLT for statistics of spin models

25

where [l] = {1, . . . , l} and x∗ := min{x1 , . . . , xl }. Note that Dlx1 ,...,xl ψ (µ ) is a symmetric function. We say that ψ is
≺-continuous at ∞ if for all µ ∈ N , we have that
lim ψ (µ|x ) = ψ (µ ).

(3.2)

|x|↑∞

Theorem 3.1 ([10, Theorem 3.1]) Let P be a simple point process on V and ψ : N → R be ≺-continuous at ∞
and assume that for all l ≥ 1
6=

∑
(y1 ,...,yl

and

1
l! (y

)∈V l

6=

∑
l
1 ,...,yl )∈V

E[|Dly1 ,...,yl ψ (P \ {y1 , . . . , yl })|1[{y1 , . . . , yl } ⊂ P]] < ∞

(3.3)

E[Dly1 ,...,yl ψ (P \ {y1 , . . . , yl })1[{y1 , . . . , yl } ⊂ P]] → 0 as l → ∞,

(3.4)

6=

where ∑ denotes that the summation is over only distinct y1 , . . . , yl . Then E[ψ (P)] has the following factorial
moment expansion
∞

1
l=1 l! (y

E[ψ (P)] = ψ (o) + ∑

6=

Dly1 ,...,yl ψ (o)ρ (l) (y1 , . . . , yl ).

∑
1 ,...,yl

(3.5)

)∈V l

Though the actual statement of [10, Theorem 3.1] involves Palm expectations and probabilities, all of them vanish in
this discrete case due to the simple form of conditional probabilities. The proof of Theorem 1.7 shall use the above
FME for E (ψ (Pn )), where ψ (µ ) is the following product of the score functions
p

ψ (µ ) = ψk1 ,...,k p (x1 , . . . , x p ; µ ) := ∏ ξ (xi , µ )ki ,

(3.6)

i=1

with k1 , . . . , k p ≥ 1. For ease of stating the FME, we shall consider the following modified functional ψ !
p
p
ki

ψ ! (µ ) = ψk!1 ,...,k p (x1 , . . . , x p ; µ ) := ψ (µ + ∑ δx j ) = ∏ ξ xi , µ + ∑ pj=1 δx j .
j=1

(3.7)

i=1

Proposition 3.2 Let P be a clustering spin model on V as in Definition 1.3 and ξ be a local statistic (i.e., R(z, P) ≤
r for all z a.s.) as in (1.4). Then for distinct x1 , . . . , x p ∈ V , non-negative integers k1 , . . . , k p and 1 ≤ n ≤ ∞, the
functional ψ ! at (3.7) admits the FME
k ,...,k p

mn1

L

(x1 , . . . , x p ) = E(ψk1 ,...,k p (x1 , . . . , x p ; Pn ))

p 1
= ψk!1 ,...,k p (x1 , . . . , x p ; o) + ∑l=1
l!

6=

∑ Dly1 ,...,yl ψk!1 ,...,k p (x1 , . . . , x p ; o)ρ (p+l) (x1 , . . . , x p , y1 , . . . , yl ),

(y1 ,...,yl )∈
p
(∪i=1Wr,n (xi ))(l)

p
where L p = | ∪i=1
Wr (xi )| and Wr,n (x) := Wr (x) ∩Wn .

(3.8)

26

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

p
/ ∪i=1
Wr (xi ) for some k ≤ l,
Proof Firstly, by the local property of ξ , we have that if y1 , . . . , yl ∈ V such that yk ∈
then for any spin model µ on V ,
Dly1 ,...,yl ψk1 ,...,k p (x1 , . . . , x p ; µ ) = 0.
(3.9)

The proof follows easily by noting that in (3.1), the terms of the form J ∪ {k} and J with J ⊂ [l] \ {k} cancel out
each other due to the locality property of ξ (see [11, (3.11)] for a detailed proof).
Thus, from (3.9) we have that the FME kernels for ψk!1 ,...,k p (x1 , . . . , x p ; µ ) (i.e., Dy1 ,...,yl ψk!1 ,...,k p (x1 , . . . , x p ; µ ) ) are
p
non-trivial only if y1 , . . . , yl ∈ ∪i=1
Wr (xi ). Since the correlation functions are zero if the coordinates repeat, we have
p
(l)
that ρ (y1 , . . . , yl ) is non-trivial for y1 , . . . , yl ∈ ∪i=1
Wr (xi ) only if l ≤ L p where L p is as defined in the proposition.
Thus, for l > L p
6=

∑
(y1 ,...,yl

)∈V l

p
Wr (xi ))]] = 0.
E[Dly1 ,...,yl ψk1 ,...,k p (x1 , . . . , x p ; Pn )1[{y1 , . . . , yl } ⊂ Pn ∩ (∪i=1

Since local property of ξ guarantees ≺-continuity of ψ , we only need to prove (3.3) to verify the conditions of
p
Theorem 3.1 and then we shall show that it reduces to (3.8). Setting Kp = ∑i=1
ki , note that by definition of the FME
kernels in (3.1) and from the definition of kξ k∞ at the beginning of the proof we trivially have that
K

|Dly1 ,...,yl ψk1 ,...,k p (x1 , . . . , x p ; Pn )| ≤ 2l kξ k∞p ,
where kξ k∞ = kξr k∞ is as defined in Remark 1.5. Now using triviality of FME kernels as in (3.9), we get that for
any l ≥ 1,
6=

K

∑
(y1 ,...,yl )∈V l

E[|Dly1 ,...,yl ψk1 ,...,k p (x1 , . . . , x p ; Pn )|] ≤ 2l kξ k∞p Llp < ∞,

and thus (3.3) is valid for ψk1 ,...,k p (x1 , . . . , x p ) i.e.,
k ,...,k p

mn1

Lp

1
l!
l=1

6=

Dly1 ,...,yl ψk1 ,...,k p (x1 , . . . , x p ; o)ρ (l) (y1 , . . . , yl ).

∑

(x1 , . . . , x p ) = ψk1 ,...,k p (x1 , . . . , x p ; o) + ∑

(y1 ,...,yl )
p
∈(∪i=1Wr,n (xi ))l

But by definition of ξ , we have that if xi ∈
/ µ for some 1 ≤ i ≤ p, then ψk1 ,...,k p (x1 , . . . , x p ; µ ) = 0. Thus, Dly1 ,...,yl ψk1 ,...,k p (x1 , . . . , x p ; o) =
0 unless {x1 , . . . , x p } ⊂ {y1 , . . . , yl }. Rearranging and relabelling the coordinates, we shall assume that the summand
in the above FME expansion is non-trivial only for {x1 , . . . , x p , y1 , . . . , yl } ⊂ V for some l ≥ 0. Thus, we get that
k ,...,k p

mn1

Lp

1
(p
+
l)!
l=0

(x1 , . . . , x p ) = ∑

6=

∑

(z1 ,...,z p+l )
p
∈(∪i=1Wr,n (xi )) p+l

ψ (o)ρ (p+l) (z1 , . . . , z p+l )1[{x1 , . . . , x p } ⊂ {z1 , . . . , z p+l }].
Dzp+l
1 ,...,z p+l

Now fixing an l ≥ 0. Again combining the definition of ξ with the definition of FME kernels as in (3.1), we have
that
p

ψ (o) =
Dxp+l
1 ,...,x p ,y1 ,...,yl

∑ (−1)l−|J| ψ ( ∑ δxi + ∑ δy j ) =
J⊂[l]

i=1

j∈J

∑ (−1)l−|J| ψ ! ( ∑ δy j ) = Dly1,...,yl ψ ! (o).
J⊂[l]

j∈J

By the invariance of the RHS under permutations of {x1 , . . . , x p , y1 , . . . , yl }, we have that

ψ (o)ρ (p+l) (z1 , . . . , z p+l )1[{x1 , . . . , x p } ⊂ {z1 , . . . , z p+l }] = Dly1 ,...,yl ψ ! (o),
Dzp+l
1 ,...,z p+l

CLT for statistics of spin models

27

where we have assumed that {z1 , . . . , z p+l } = {x1 , . . . , x p , y1 , . . . , yl }. Now the proof of (3.8) follows by correctly
counting {z1 , . . . , z p+l } = {x1 , . . . , x p , y1 , . . . , yl }.
✷
Proof (Theorem 1.7) Case (1) - ‘Local’ score functions : Firstly, we shall prove the bounded stabilisation radius case
i.e., assume a.s. for all x ∈ V , R(x, P) ≤ r. Choose s > 2r. Fix x1 , . . . , x p+q and assume that s = d({x1 , . . . , x p }, {x p+1 , . . . , x p+q }).
Setting
p
p+q
L p = | ∪i=1
Wr (xi )|, Lq = | ∪qi=1 Wr (x p+i )|, L p+q = | ∪i=1
Wr (xi )|,
(3.10)
p
we see that trivially L p+q = L p + Lq as the sets ∪i=1
Wr (xi ), ∪qi=1Wr (x p+i ) are disjoint. Note that L p ≤ pwr and
p
p+q
similarly for Lq , L p+q . Further set Kp = ∑i=1 ki , Kq = ∑qi=1 k p+i , K = ∑i=1
ki and let ψ ! be as defined in (3.7). Thus
applying Proposition 3.2 we obtain that

mk1 ,...,k p+q (x1 , . . . , x p+q )
L p+q

=

∑
l=0
L p+q

=

∑
l=0

1
l!
1
l!

(y1 ,...yl )∈(∪i=1 Wr,n (xi ))(l)
l

∑

j=0

L p+q l

=

Dly1 ,...,yl ψk!1 ,...,k p+q (x1 , . . . , x p+q ; o)ρ (l+p+q) (x1 , . . . , x p+q , y1 , . . . , yl ),

∑
p+q

∑∑

l=0 j=0

∗
l!
Dly1 ,...,yl ψk!1 ,...,k p+q (x1 , . . . , x p+q ; o) × ρ (l+p+q) (x1 , . . . , x p+q , y1 , . . . , yl ),
j!(l − j)! ∑
y

∗
1
∑ (−1)l−|J| ψk!1,...,k p+q (x1 , . . . , x p+q; ∑ δy j ) ρ (l+p+q) (x1 , . . . , x p+q, y1 , . . . , yl ), (3.11)
j!(l − j)! ∑
y J⊂[l]
j∈J

∗

p
where, ∑ denotes summation over the set {(y1 , . . . yl ) ∈ (∪i=1
Wr,n (xi ))( j) × (∪qi=1Wr,n (x p+i ))(l− j) }. Proceeding again
y

∗

as above, and writing ∑ as summation over the set {(y1 , . . . , yl1 , z1 , . . . , zl2 ) ∈
y,z

p
(∪i=1
Wr,n (xi ))(l1 ) × (∪qi=1Wr,n (x p+i ))(l2 ) }, we derive that

mk1 ,...,k p (x1 , . . . , x p )mk p+1 ,...,k p+q (x p+1 , . . . , x p+q )
Lp

=

Lq

1 ∗
∑
=0 l1 ! l2 ! y,z J

l1 =0 l2

(−1)l1 +l2 −|J1 |−|J2 | ψk!1 ,...,k p (x1 , . . . , x p ;

∑

∑ ∑

1 ⊂[l1 ],J2 ⊂[l2 ]

×ψk! p+1 ,...,k p (x p+1 , . . . , x p+q ;
L p +Lq l

=

∑ ∑

l=0 j=0

∗

1
j!(l − j)! ∑
y J

∑ δz j ) ρ (l1 +p) (x1 , . . . , x p+q, y1 , . . . , yl1 )ρ (l2 +q) (x p+1 , . . . , x p+q, y1 , . . . , yl2 ),

j∈J2

∑

1 ⊂[ j],J2 ⊂[l]\[ j]

×ψk! p+1 ,...,k p+q (x p+1 , . . . , x p+q ;
L p+q l

=

∑∑

l=0 j=0

∑ δy j )

j∈J1

(−1)l−|J1 |−|J2 | ψk!1 ,...,k p (x1 , . . . , x p ; ∑ δyi )

∑ δyi ) ρ

i∈J1

( j+p)

(x1 , . . . , x p+q , y1 , . . . , y j )ρ (l− j+q) (x p+1 , . . . , x p+q , y j+1 , . . . , yl ),

i∈J2

∗

1
∑ (−1)l−|J|ψk!1 ,...,k p+q (x1 , . . . , x p+q ; ∑ δyi )
j!(l − j)! ∑
y J⊂[l]
i∈J

×ρ ( j+p) (x1 , . . . , x p+q , y1 , . . . , y j )ρ (l− j+q) (x p+1 , . . . , x p+q , y j+1 , . . . , yl ),

(3.12)

where in the last equality we have used the fact that for all J ⊂ [l] with J1 = J ∩ [ j], J2 = J \ J1 , we have that

ψk!1 ,...,k p+q (x1 , . . . , x p+q ; ∑ δyi ) = ψk!1 ,...,k p (x1 , . . . , x p ; ∑ δyi )ψk! p+1 ,...,k p (x p+1 , . . . , x p+q ; ∑ δyi ).
i∈J

i∈J1

i∈J2

28

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

This follows from the definitions of ψ ! (see (3.7)), R(·, P) and that R(xi , P) ≤ r for all i ∈ [p + q]. We now set
L̃ p+q = L p+q + p + q and note that L̃ p+q + p + q ≤ (p + q)(1 + wr ) ≤ K(1 + wr ). Comparing (3.11) with (3.12), and
using Definition 1.3 and Remark 1.5, we have that
|mk1 ,...,k p+q (x1 , . . . , x p+q ) − mk1 ,...,k p (x1 , . . . , x p )mk p+1 ,...,k p+q (x p+1 , . . . , x p+q )|
L p+q l

≤

Cl+p+q φ (cl+p+q s)
j!(l − j)!
j=0

∑∑
l=0

∑ |ψk!1,...,k p+q (x1 , . . . , x p+q ; ∑ δyi )|,
i∈J

(y1 ,...,yl )∈
J⊂[l]
p
q
(∪i=1Wr,n (xi ))( j) ×(∪i=1Wr,n (x p+i ))(l− j)

L p+q l

≤ φ (cL̃ p+q s)

∑

∑∑

l=0 j=0

Kwr
Cl+p+q l
Cl+K (4Kwr )l
2 kξ kK∞ Llp+q ≤ φ (cK(1+wr ) s)kξ kK∞ ∑
.
j!(l − j)!
l!
l=0

(3.13)

Thus, we have proven clustering of mixed moment functions for local score functions with φ̃ (·) = φ (·), C̃K =
r
kξ kK∞ ∑Kw
l=0
as well.

Cl+K (4Kwr )l
, c̃K
l!

= cK(1+wr ) . Also trivially note that summability of φ as in (1.9) implies summability of φ̃

p
Case (2) - ‘Quasi-local’ score functions : Let us fix x1 , . . . , x p+q , and assume that s = d({xi }i=1
, {x p+ j }qj=1 ). Further
γ
without loss of generality, let s > 4 and n ≫ s. Set t = t(s) := (s/4) for a γ ∈ (0, 1) to be chosen later. Define for
all xi ,
ξ̃ (xi , Pn ) = ξ (x, Pn )1[R(xi , Pn ) ≤ t],

By definition of radius of stabilisation, we have that ξ̃ (z, P) = ξ (z, P ∩ Wt (z)) for all z ∈ V and so ξ̃ is a local
p
ξ̃ (xi , Pn )ki ]. Now, by
statistic with radius of stabilisation R̃(x, Pn ) ≤ t. Further set m̃k1 ,...,k p (x1 , . . . , x p ; n) := E[∏i=1
Hölder’s inequality, the moment condition (1.7) and exponentially quasi-locality (1.6), we have that
p

p

|mk1 ,...,k p (x1 , . . . , x p ; n) − m̃k1 ,...,k p (x1 , . . . , x p ; n)| ≤ |E[∏ ξ (xi , Pn )ki ] − E[∏ ξ̃ (xi , Pn )ki ]|
≤

K /(K +1)
M p p p (pAϕ (t))1/(Kp +1)

≤

i=1
K/(K+1)
1/(K+1)
MK
(AK ϕ (t))
.

i=1

Further, for any reals A′ , B′ , A′′ , B′′ with |B′′ | ≤ |B′ | we have that |A′′ B′′ − A′ B′ | ≤ (|A′ | + |B′ |)(|A′′ − A′ | + |B′′ − B′ |).
Using this bound and Hölder’s inequality as above, we have that
|mk1 ,...,k p+q (x1 , . . . , x p+q ; n) − mk1 ,...,k p (x1 , . . . , x p ; n)mk p+1 ,...,k p+q (x p+1 , . . . , x p+q ; n)|
≤ |mk1 ,...,k p+q (x1 , . . . , x p+q ; n) − m̃k1 ,...,k p+q (x1 , . . . , x p+q ; n)|

+|mk1 ,...,k p (x1 , . . . , x p ; n)mk p+1 ,...,k p+q (x p+1 , . . . , x p+q ; n) − m̃k1 ,...,k p (x1 , . . . , x p )m̃k p+1 ,...,k p+q (x p+1 , . . . , x p+q )|

+|m̃k1 ,...,k p+q (x1 , . . . , x p+q ; n) − m̃k1 ,...,k p (x1 , . . . , x p ; n)m̃k p+1 ,...,k p+q (x p+1 , . . . , x p+q ; n)|,
2K/(K+1)

≤ 5MK

(AK ϕ (t))1/(K+1) + |m̃k1 ,...,k p+q (x1 , . . . , x p+q ) − m̃k1 ,...,k p (x1 , . . . , x p ; n)m̃k p+1 ,...,k p+q (x p+1 , . . . , x p+q ; n)|,

(3.14)

p
Set L p (t) = | ∪i=1
Wt (xi )| and similarly Lq (t), L p+q (t). Now using (3.13) for ξ̃ , we have that
Kwt

Cl+p+q (4Kwt )l
.
l!
l=0

|m̃k1 ,...,k p+q (x1 , . . . , x p+q ) − m̃k1 ,...,k p (x1 , . . . , x p )m̃k p+1 ,...,k p+q (x p+1 , . . . , x p+q )| ≤ φ (cL̃ p+q s)kξt kK∞ ∑

CLT for statistics of spin models

29

Now using the polynomial growth of wt (see Definition (1.1)) exponential growth condition (1.11) and exponential
clustering condition (1.10), we derive that
|m̃k1 ,...,k p+q (x1 , . . . , x p+q ) − m̃k1 ,...,k p (x1 , . . . , x p )m̃k p+1 ,...,k p+q (x p+1 , . . . , x p+q )|

≤ c′1 exp(−csb ) × exp(CKt κ ) × c′3 exp(c′4t d/a ) ≤ c′5 exp(−csb ) × exp(c′6 sγ (κ +d/a) ),
where all the c′i are constants depending on a, d and K. Thus, if we choose γ such that γ (κ + d/a) < b/2, we have
that
b/2
|m̃k1 ,...,k p+q (x1 , . . . , x p+q ) − m̃k1 ,...,k p (x1 , . . . , x p )m̃k p+1 ,...,k p+q (x p+1 , . . . , x p+q )| ≤ C̃K e−c̃K s ,
for two constants Ck′ , c′k ∈ (0, ∞) as required to complete the proof of the theorem.

✷

3.2 Proof of the general central limit theorem - Theorem 1.9
As stated earlier, we shall only be giving a sketch of the proof here as a similar theorem or variant has appeared in
all the papers using the cumulant method such as[67,68,51,5,71,11,26,8]. For more details, see the proof in [11,
Section 4.4]. We shall use the arguments therein adapted to the discrete Cayley graph case by essentially setting
f ≡ 1 and defining the moment measures appropriately.
We shall first define cumulants and state their relations to moments. For a random variable Y with all moments
being finite, the cumulants Sk (Y ), k ≥ 1 are formally defined as the coefficients in the power series expansion of the
log Laplace transform of Y in the negative domain i.e.,
log E(etY ) = log(1 + ∑ Mk (Y )t k ) =
k≥1

∑ Sk (Y )t k ,
k≥1

E(Y k )

where Mk (Y ) :=
is the kth moment of Y . Further, by formally manipulating the above series expansion, we
can derive the following useful relation between moments and cumulants (see [73, Proposition 3.2.1]) :
l

Sk =

∑

γ ={γ (1),...,γ (l)}∈Π [k]

(−1)l−1 (l − 1)! ∏ M|γ (l)| ,

(3.15)

i=1

where Π [k] denotes the set of all unordered partitions of the set [k] = {1, . . . , k}, a partition γ ∈ Π [k] is represented
as γ = {γ (1), . . . , γ (l)} with l representing the number of classes and γ (i), the elements in the ith class. |γ (i)| denotes
the cardinality of γ (i) and we shall also denote l by |γ |. Equipped with this relation, cumulants can be defined as
long as all moments exist without being concerned about the existence of the log Laplace transform. We refer the
reader to [73, Chapter 3] for more details on cumulants.
Proof (Proof of Theorem 1.9) For convenience, we shall drop the superscript X in our notations referring to the
underlying random field Xn . Let Mk,n , Sk,n denote the moments and cumulants of Hn respectively and the moments
and cumulants of H̄n = H√n −E(Hn ) are denoted by M̄k,n , S̄k,n respectively. The existence of the above moments and
Var(Hn )

cumulants follow from the assumption of moment condition in the theorem and Hölder’s inequality. By (3.15),
we have that S̄1,n = 0, S̄2,n = 1 for all n ≥ 1. Further, for a random variable Y and a constant c ∈ R, we have that
Sk (cY ) = ck Sk (Y ) for all k ≥ 1 and Sk (Y + c) = Sk (Y ) for all k ≥ 2. Thus we obtain that for all k ≥ 2,
S̄k,n = Var(Hn )−k/2 Sk,n .

30

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

Now by the above relation, the variance lower bound assumption and an extension of the classical cumulant method
using Marcinkiewicz’s theorem (see [51, Theorem 1]), we have the required normal convergence if we show that for
all k ≥ 3
Sk,n = O(wn ).
(3.16)
The rest of the proof will consist of showing the above bound. First define mixed moment functions of the random
ki
field Xn as m(k1 ,...,k p ) (x1 , . . . , x p ; n) = E(∏ni=1 Xn,x
i ) for 1 ≤ n < ∞, x1 , . . . , x p ∈ Wn and k1 , . . . , k p ≥ 1. Now define
(k ,...,k p )

the truncated mixed moment functions or Ursell functions m⊤ 1
for 1 ≤ i ≤ p and inductively,
(k ,...,k p )

m⊤ 1

(x1 , . . . , x p ; n) := m(k1 ,...,k p ) (x1 , . . . , x p ; n) −

(1)

(x1 , . . . , x p ; n) as follows : Set m⊤ (·) = m(1) (·)
l

∑

γ ={γ (1),...,γ (l)}∈Π [p]
l>1

(k : j∈γ (i))

∏ m⊤ j
i=1

(x j ; j ∈ γ (i); n).

(3.17)

The mixed moment functions and Ursell functions exist because of the moment condition assumed in the theorem
and Hölder’s inequality. The Ursell functions are crucial to our analysis because by using the above definition and
that of cumulant, one can derive that (see [44, Section 2] or [11, Section 4.4])
(k ,...,k p )

∑

∑

Sk,n =

k1 +···+k p =k x1 ,...,x p ∈Wn
k1 ,··· ,k p ≥1

m⊤ 1

(x1 , . . . , x p ; n).

Observe that for all k ≥ 1, Π [k] is finite and so the first summand is over finitely many terms for all n ≥ 1. Hence
the proof of (3.16) follows if for any k ≥ 1 and all k1 , . . . , k p ≥ 1 such that k1 + · · · + k p = k, we show that

∑

x1 ,...,x p ∈Wn

(k ,...,k p )

|m⊤ 1

(x1 , . . . , x p ; n)| = O(wn ).

(3.18)

Now onwards, we fix a 1 ≤ n < ∞ and drop the reference to n in the notation for Ursell functions and mixed
moment functions. A partition γ = {γ (1), . . . , γ (l)} is said to refine a partition σ = {σ (1), . . . , σ (l1 )} if for all i ∈ [l],
γ (i) ⊂ σ ( j) for some j ∈ [l1 ]. Else, we say that γ mixes σ . Using this definition and the definition of Ursell functions,
we derive that for any I ( [p]
l

c

m(k j : j∈I) (x j ; j ∈ I)m(k j : j∈I ) (x j ; j ∈ I c ) =

∑

γ ={γ (1),...,γ (l)}∈Π [p]
γ refines {I,I c }

(k : j∈γ (i))

∏ m⊤ j
i=1

(x j ; j ∈ γ (i)),

and consequently we have that
(k ,...,k p )

m⊤ 1

c

(x1 , . . . , x p ; n) = m(k1 ,...,k p ) (x1 , . . . , x p ; n) − m(k j : j∈I) (x j ; j ∈ I)m(k j : j∈I ) (x j ; j ∈ I c )
l

∑

+

γ ={γ (1),...,γ (l)}∈Π [p]
l>1,γ mixes {I,I c }

(k : j∈γ (i))

∏ m⊤ j
i=1

(x j ; j ∈ γ (i)).

(3.19)

Now, by induction one can show that for all p ≥ 1 and for any configuration x1 , . . . , x p ∈ V , there exists a partition
I, I c of [p] with d({x j ; j ∈ I}, {x j ; j ∈ I c }) ≥ diam(x1 , . . . , x p )/(p − 1) where diam is the graph diameter of the
set and defined as diam(x1 , . . . , x p ) := supi6= j k − xi + x j k. Thus using (3.19), the above partition and the clustering
condition (1.19), we can inductively show that
sup

sup

1≤n≤∞ x1 ,...,x p ∈Wn

(k ,...,k p )

|m⊤ 1

⊤
(x1 , . . . , x p ; n)| ≤ C⊤
p φ (c p diam(x1 , . . . , x p )),

(3.20)

CLT for statistics of spin models

31

⊤
where C⊤
p , c are finite constants. Now using (3.20), we derive that

sup

∑

(k ,...,k p )

x1 ∈Wn x2 ,...,x p ∈Wn

|m⊤ 1

(x1 , . . . , x p ; n)| ≤ sup C⊤
p
x1 ∈Wn

= sup
x1 ∈Wn

p

∑

∏ φ (c⊤p d(x1 , xi))1/(p−1) ,

x2 ,...,x p ∈Wn i=2

C⊤
p(

∑

1/(p−1) p−1
φ (c⊤
)
< ∞,
p d(x1 , x))

(3.21)

x∈Wn

where the finiteness is due to the summability of φ as in (1.9). The above bound trivially implies (3.18) and thus we
have completed the proof.

3.3 Proof of the weak law of large numbers and the multivariate central limit theorem -Theorems 1.8 and
1.11
ξ

Armed with the powerful clustering result for Hn , we can now give the proofs of the weak law of large numbers and
the central limit theorems.
Proof (Proof of Theorem 1.8) We will prove the results under Assumption 2. The proof can be trivially adapted
under Assumption 1 by setting ϕ (t) = 0 for large enough t and using summability of φ̃ as in (1.9) due to Theorem
1.7. We shall first show that
ξ
1
w−1
(3.22)
n E(Hn ) → m (O),
and this along with the variance asymptotics to be proved and Chebyshev’s inequality suffices to prove the weak law
of large numbers. By stationarity of P, we have that E(ξ (x, P)) = E(ξ (O, P)) and hence
−1
ξ
|w−1
n E(Hn ) − E(ξ (O, P))| ≤ wn

∑

x∈Wn

E(|ξ (x, Pn ) − ξ (x, P)|1[R(x, Pn) ≥ d(x, ∂ Wn )]),

1/p
≤ 2w−1
n (M p )

∑

x∈Wn

P(R(x, Pn ) ≥ d(x, ∂ Wn ))1/q ,

(3.23)

where in the last inequality we have used Hölder’s inequality with q ≥ 1 such that 1p + 1q = 1 (p > 1 is chosen as in
the theorem assumption) and the moment condition (1.7). By the property of radius of stabilisation (1.6), we have
that for any r > 0,
w−1
n

∑

x∈Wn

P(R(x, Pn ) ≥ d(x, ∂ Wn ))1/q ≤ w−1
n A

∑

x∈Wn

ϕ (n − |x|)1/q + wn−1 |∂ Wn |
r

1/q
−1
1/q
+ Aw−1
+ A ∑ w−1
≤ w−1
n |Wn \Wn−r | ≤ wn |Wn−r |Aϕ (r)
n |Wn−r |Aϕ (r)
n− j |∂ Wn− j |.
j=0

Thus, by the b-amenability of G we conclude that the second term above converges to 0, and we have
lim sup w−1
n
n→∞

∑

x∈Wn

P(R(x, Pn ) ≥ d(x, ∂ Wn ))1/q ≤ Aϕ (r)1/q ,

for any r > 0. Now letting r → ∞, the proof of expectation asymptotics (3.22) is complete.
Now moving onto variance asymptotics, we have that
Var(Hnξ ) =

∑

x,y∈Wn


E(ξ (x, Pn )ξ (y, Pn )) − E(ξ (x, Pn ))E(ξ (y, Pn )) .

(3.24)

32

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

For x ∈ Wn , we set Pnx = P ∩ (Wn − x) and Wn′ (z) = Wn ∩ (−z +Wn )c . Now, by change of variables, we have

E(ξ (x, Pn )ξ (z + x, Pn )) − E(ξ (x, Pn))E(ξ (z + x, Pn )) ,
Var(Hnξ ) =
∑
=

x∈Wn ,z∈Wn −x

1,1

∑

z∈Zd ,x∈Wn

=

∑
z∈Zd ,x∈Wn

[m (0, z; Pnx ) − m1 (O; Pnx )m1 (z; Pnx )]1[x ∈ (−z +Wn)],

[m1,1 (0, z; Pnx ) − m1 (O; Pnx )m1 (z; Pnx )] −

∑
z,x∈Zd

[m1,1 (0, z; Pnx ) − m1 (O; Pnx )m1 (z; Pnx )]1[x ∈ Wn′ (z)].
(3.25)

By Theorem 1.7, we have that
w−1
n

∑
z,x∈Zd

|m1,1 (0, z; Pnx ) − m1 (O; Pnx )m1 (z; Pnx )|1[x ∈ Wn′ (z)] ≤ w−1
n C̃2

∑ φ̃ (c̃2 |z|)|Wn′ (z)|.

z∈Zd

′
Since ∑z∈Zd φ̃ (c̃2 |z|) is summable as φ̃ is fast-decreasing, w−1
n |Wn (z)| ≤ 1 and further, for all z ∈ V , by b-amenability
of the discrete Cayley graph, we have that
′
−1
w−1
n |Wn (z)| ≤ wn

n

∑
j=n−|z|

|∂ W j | → 0.

So, we can use dominated convergence theorem to conclude that

∑
z,x∈Zd

[m1,1 (0, z; Pnx ) − m1 (O; Pnx )m1 (z; Pnx )]1[x ∈ Wn′ (z)] → 0,

as n → ∞ as required. Now, the proof is complete if we show that
w−1
n

∑
z∈Zd ,x∈Wn

[m1,1 (0, z; Pnx ) − m1 (O; Pnx )m1 (z; Pnx )] →

∑ (m1,1 (O, z) − m1 (O)2 ).

(3.26)

z∈Zd

The above convergence result follows if we show that

∑ w−1
n ∑

x∈Wn

z∈Zd

|m1,1 (0, z; Pnx ) − m1 (O; Pnx )m1 (z; Pnx ) − m1,1 (O, z) + m1 (O)2 | → 0.

From Theorem 1.7, we have that the modulus is bounded above by 2C̃2 φ̃ (c̃2 |z|) and hence we have that
w−1
n

∑

x∈Wn

|m1,1 (0, z; Pnx ) − m1 (O; Pnx )m1 (z; Pnx ) − m1,1 (O, z) + m1 (O)2 | ≤ 2C̃2 φ̃ (c̃2 |z|),

and since the RHS is summable over z ∈ V , we can apply dominated convergence to conclude (3.26) as required
provided we prove that for all z ∈ V as n → ∞,
w−1
n

∑

x∈Wn

|m1,1 (0, z; Pnx ) − m1 (O; Pnx )m1 (z; Pnx ) − m1,1 (O, z) + m1 (O)2 | → 0.

(3.27)

This can be shown term-wise. By using the arguments in the proof of expectation asymptotics (3.22), we derive that

1
x
1
1
x
1
1
x
1
x
1
2
−1
w−1
n ∑ |m (O; Pn )m (z; Pn ) − m (O) | ≤ 2wn M1 ∑ |m (O; Pn ) − m (O)| + |m (z; Pn ) − m (O)| .
x∈Wn

x∈Wn

Now note that by Cauchy-Schwarz inequality and definition of radius of stabilisation,
1/2

|m1 (z; Pnx ) − m1 (z)| ≤ 2M2 P(R(z, Pn ) ≥ d(x, ∂ Wn ))1/2 ,

CLT for statistics of spin models

33

and similarly for the other term |m1 (O; Pnx ) − m1 (O)|. Now we can argue exactly like in the proof of (3.22) to obtain
that
1
x
1
x
1
2
w−1
n ∑ |m (O; Pn )m (z; Pn ) − m (O) | → 0.
x∈Wn

Now, again using triangle inequality, we have that
|E (ξ (O, Pnx )ξ (z, Pnx)) − E(ξ (O, P)ξ (z, P))|


≤ (E(ξ (O, Pnx )2 )1/2 + E(ξ (z, P)2)1/2 ) × E(|ξ (O, Pnx ) − ξ (O, P)|2)1/2 + E(|ξ (z, Pnx) − ξ (z, P)|2)1/2 ,

1/2
≤ 2M2 E(|ξ (O, Pnx ) − ξ (O, P)|2)1/2 + E(|ξ (z, Pnx) − ξ (z, P)|2)1/2 .

Now, choose q such that 2/p + 1/q = 1 for p > 2 as assumed in the Theorem 1.8. Then we can derive that


E |ξ (O, Pnx ) − ξ (O, P)|2 ≤ E |ξ (O, Pnx) − ξ (O, P)|21[R(O, Pn ) ≥ d(x, ∂ Wn )]


≤ E ξ (O, Pnx )2 1[R(O, Pn ) ≥ d(x, ∂ Wn )] + E ξ (O, P)2 1[R(O, Pn ) ≥ d(x, ∂ Wn )]

+2E |ξ (O, Pnx )ξ (O, P)|1[R(O, Pn) ≥ d(x, ∂ Wn )] ,

≤ E(|ξ (O, Pnx | p )2/p P(R(O, Pn ) ≥ d(x, ∂ Wn ))1/q + E(|ξ (O, P)| p)2/p P(R(O, Pn ) ≥ d(x, ∂ Wn ))1/q
+2E(|ξ (O, Pnx |)| p )1/p E(|ξ (O, P)| p)1/p P(R(O, Pn ) ≥ d(x, ∂ Wn ))1/q ,
2/p

≤ 4M p P(R(O, Pn ) ≥ d(x, ∂ Wn ))1/q ,
and similarly,

2/P

E(|ξ (z, Pnx ) − ξ (z, P)|2) ≤ 4M p P(R(z, Pn ) ≥ d(x, ∂ Wn ))1/q .
Thus, we have that
1,1
x
1,1
w−1
n ∑x∈Wn |m (0, z; Pn ) − m (O, z)|

≤

(2+p)/2p −1
wn ×

4M p

∑ (P(R(O, Pn) ≥ d(x, ∂ Wn ))1/2q + P(R(z, Pn) ≥ d(x, ∂ Wn))1/2q )

x∈Wn

→ 0,

and this proves (3.27) and consequently (3.26) as well. This completes the proof of variance asymptotics.

✷

Proof (Proof of Theorem 1.11) The proof of this theorem is via Cramér-Wold theorem. Let N = (N1 , . . . , Nk ) be the
Gaussian vector with mean zero and covariance matrix Σ := (Σi, j ) as defined in the Theorem. By the Cramér–Wold
theorem, it suffices to show that for all t ∈ Rk , we have
 

H̄n − E(H̄n )
d
t,
→ ht, Ni,
(3.28)
√
wn
where ht, si := ∑ki=1 ti si stands for the usual dot product between two k-dimensional real vectors. Fix t ∈ Rk . For
ξ
x ∈ V and spin model µ , set ξ (x, µ ) := ∑i ti ξi (x, µ ). Note that ht, H̄n i = Hn . Further, we observe that (ξ , P) satisfies
Assumption 1 (resp. Assumption 2) of Theorem 1.10 if all the (ξi , P), i = 1, . . . , k satisfies Assumption 1 (resp.
Assumption 2) of Theorem 1.10. Further, as with the variance asymptotics in Theorem 1.8, we can show that for all
1 ≤ i, j ≤ k,


ξ

j
ξi
w−1
→ Σi, j .
n Cov Hn , Hn

34

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

One can also deduce the above asymptotic directly by noting that

 1
ξ
ξ +ξ
ξ 
Cov Hnξi , Hn j = Var(Hn i j ) − Var(Hnξi ) − Var(Hn j ) ,
2

and the asymptotics for each of the three terms on the RHS follow from Theorem 1.8. After appropriate cancellations, we obtain the desired covariance asymptotics. Thus, we derive that
ξ
w−1
n Var(Hn ) →

∑

tit j Σi, j = Var(ht, Ni).

1≤i, j≤k

ξ

If Var(ht, Ni) = 0, then ht, Ni = 0 a.s. i.e., a degenerate random variable. Then, since w−1
n Var(Hn ) → 0, we derive
(3.28) via Chebyshev’s inequality trivially in this case. Alternatively, if Var(ht, Ni) > 0, then since (ξ , P) satisfies
ξ
either Assumptions 1 or 2 in Theorem 1.10 and we also have that Var(Hn ) = Ω (wn ), we can conclude (3.28) from
Theorem 1.10 and thereby completing the proof of our multivariate central limit theorem.

3.4 Proof of Theorems 1.12 and 1.13
Before we start documenting the proofs of related to central limit theorems for sums of score functions of random
fields satisfying certain mixing conditions, we remind the reader that we have so far stayed clear of an important
aspect, that of non-degeneracy of the limiting distribution. We refer the reader to [25], and references therein, for a
detailed discussion on various conditions / assumptions which ensure positivity of the limiting variance. We shall
not delve any deeper into this very delicate issue of the non-degeneracy of the limiting distribution.
Proof (Proof of Theorem 1.12) We begin the proof by first noting that we shall follow the proof of similar limit
theorem mentioned in [19, Theorem 1.1]. Observe that the strong α -mixing condition as stated in [19] is weaker
than the strong α -mixing we defined in Section 1.3. Therefore, if X satisfies the strong α -mixing condition stated
in our Section 1.3, then X also satisfies the strong α -mixing condition stated in [19].
Consider the total mass,
(3.29)
Hnξ = ∑ ξ (x, Xn ).
x∈Wn

Note that since ξ is a local functional of the point process with with finite radius, say r, the random variables
ξ (x, Xn ) are not identically distributed for x ∈ Xn . However, as observed earlier, the random field {ξ (x, X )}x∈G
enξ = ∑x∈W ξ (x, X ), and show that the
is indeed stationary. We, therefore shall prove the invariance principle for H
n
ξ
ξ
en ) is negligible under the volume scaling. We shall break the proof into two steps: first we shall
difference (Hhn − H
i
ξ
enξ )2 is o(wn ), thereafter we shall prove the invariance principle for H
enξ .
prove that E (Hn − H
S TEP 1: Let r (fixed) be the range of the local functional ξ , we have


enξ =
Hnξ − H
∑ [ξ (x, Xn) − ξ (x, X )] .

(3.30)

x∈Wn \Wn−r

Writing ζx,n = (ξ (x, Xn ) − ξ (x, X )), consider the triangular array {ζx,n : x ∈ Wn \ Wn−r , n ≥ 1}. Clearly, ζx,n is
another local functional with range r. Therefore, Cov (ζx,n , ζy,n ) ≤ ρ (|x − y| + 2r). Now let j be the index for which
ρ ( j) < 1, then following the same arguments as set forth in [19], we conclude that

2
enξ ≤ C ∑ E (ξ (x, Xn ) − ξ (x, X ))2 ,
E Hnξ − H
x∈Wn \Wn−r

CLT for statistics of spin models

35

where the constant C is given by 2 ( j − 2r)d / (1 − ρ ( j)). Using now the regularity assumption concerning the moments of the local score function ξ , we conclude that,

2
enξ ≤ C1 |Wn \Wn−r | .
E Hnξ − H
Next, using the estimate obtained earlier in the proof of Theorem 1.8, we obtain the desired conclusion
2

enξ = o(wn ).
E Hnξ − H
S TEP 2: We shall now move to the second step of our proof, which involves the invariance principle
 
enξ − E H
enξ
H
d
⇒ N(0, σ 2 ).
√
wn

(3.31)

First we note that the Lindeberg condition for X is automatically satisfied due to the stationarity of the field X ,
and by the same argument, the random field ξ also does satisfy the Lindeberg condition. Then, we observe that if
X is an α -mixing random field then so is ξ (·, X ), which in turn implies that ξ satisfies the mixing condition stated
in [19, Theorem 1.1]. More precisely, using the same notation as set forth earlier, and writing α and αξ for the
mixing coefficients corresponding to the field X and ξ , respectively, we have αξ (s) ≤ α (s + 2r), implying that the
random field ξ is α -mixing if the original field X is α -mixing. Now we invoke Theorem 1.1 of [19] to conclude
✷
the statement (3.31), which in turn, concludes the proof.
Next, in order to prove Theorem 1.13, we shall use the following technical observation.
Proposition 3.3 Let G be a discrete Cayley graph and together with (ξ , X ) such that G has polynomial growth
as in Definition 1.1, X is an exponentially mixing spin model, ξ is a exponentially quasi-local score function as
in (1.6) satisfying the p-moment condition (1.7) for all p ≥ 1. Then, the random field {ξ (x, Xn )}x∈Wn , 1 ≤ n ≤ ∞
satisfies clustering in terms of mixed moments, i.e., there exists constants C̃K , c̃K such that for all x1 , . . . , x p+q with
p+q
s = d({x1 , . . . , x p }, {x p+1 , . . . , x p+q }) and K = ∑i=1
ki , ki ≥ 1 for i = 1, . . . , p + q, such that writing
"
#
p

mk1 ,...,k p (x1 , . . . , x p ) := E

∏ ξ (xi, Xn )ki

,

i=1

we have that
k ,...,k p+q

|mn1

k ,...,k p

(z1 , . . . , z p+q ) − mn1

k

(z1 , . . . , z p )mnp+1

,...,k p+q

(z p+1 , . . . , z p+q )| ≤ C̃K φ̃ (c̃K s),

(3.32)

for a fast decreasing function φ̃ . Further, under Assumption 1 if φ is summable as in (1.9), so is φ̃ .
Proof We shall now try to mimic the arguments used in the proof of Theorem 1.7. Let us fix x1 , . . . , x p+q and
assume that s = d({x1 , . . . , x p }, {x p+1 , . . . , x p+q }). Further, without loss of generality, let s > 4 and n ≥ 10s. Set
t = t(s) := (s/4)γ for a γ ∈ (0, 1) to be chosen later. Define for all xi

ξ̃ (xi , Xn ) = ξ (x, Xn )1[R(xi , Xn ) ≤ t].

p
Clearly, ξ̃ is a local statistic with radius of stabilisation R̃(x, Xn ) ≤ t. Further set m̃k1 ,...,k p (x1 , . . . , x p ; n) := E[∏i=1
ξ̃ (xi , Xn )ki ].
Now, proceeding as in the proof of (3.14) and again using Hölder’s inequality, the moment condition (1.7) and exponentially quasi-locality (1.6), we have that

|mk1 ,...,k p+q (x1 , . . . , x p+q ; n) − mk1 ,...,k p (x1 , . . . , x p ; n)mk p+1 ,...,k p+q (x p+1 , . . . , x p+q ; n)|
2K/(K+1)

≤ 5MK

(AK ϕ (t))1/(K+1)

k1 ,...,k p+q

+|m̃

(x1 , . . . , x p+q ) − m̃k1 ,...,k p (x1 , . . . , x p ; n)m̃k p+1 ,...,k p+q (x p+1 , . . . , x p+q ; n)|.

(3.33)

36

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

Note, however, by the choice of t and ϕ , we can replace (AK ϕ (t))1/(K+1) by an appropriate fast decaying functional
φ̃1 (cK s). Then, using Theorem 3 in Section 1.2.2 of [25] for the local statistics ξ̃ (xi , Xn ), we obtain
k ,...,k p+q

|m̃n1

k ,...,k p

k

≤ 2 [α (s − 2t)]1/3 E

p

∏
i=1

j=1

Applying (1.7) and using Hölder’s inequality, we can bound

! 1/3 
3

p

E

∏

ξ ki (xi , Xn )

i=1

q



Finally, setting φ̃2 (s) = [α (s − 2t)]1/3 , we have
k ,...,k p+q

|m̃n1

,...,k p+q

(x1 , . . . , x p )m̃np+1
(x p+1 , . . . , x p+q )|
!3 1/3 
!3 1/3
q
 E ∏ ξ k p+ j (x p+ j , Xn )
 .
ξ ki (xi , Xn )

(x1 , . . . , x p+q ) − m̃n1


k ,...,k p

(z1 , . . . , z p+q ) − m̃n1

E

∏

j=1

!3 1/3
 ≤ C̃1,K .
ξ k p+ j (x p+ j , Xn )
k

(z1 , . . . , z p )m̃np+1

,...,k p+q

(z p+1 , . . . , z p+q )| ≤ C̃1,K φ̃2 (s).
2K/(K+1)

Collating (3.33) and (3.34), now we set φ̃ (cK s) = max(φ̃1 (cK s), φ̃2 (s)), and C̃K = 5MK
the exponential decay of the α -mixing coefficient, we conclude this theorem.

(3.34)

+ C̃1,K . Finally, using
✷

Proof (Proof of Theorem 1.13) The proof follows from Proposition 3.3 and Theorem 1.9.

✷

4 Appendix: Void probability estimate
Lemma 4.1 Let P be an exponential clustering spin model as in (1.10). Then, there exists constants A, a′ , ν ∈ (0, ∞)
such that for any t > 0,
′ ν
sup sup P(Pn (Wt (z)) = 0) ≤ Ae−a t .
(4.1)
n≥1 z∈Wn

Proof Without loss of generality, we shall assume that c = 1 in (1.10). Define the random field Zx = 1[x ∈ P] for
x ∈ V . Thus, P(Pn (S) = 0) = E(∏x∈S∩Wn (1 − Zx )). Firstly, let P, Q ⊂ V such that s := d(P, Q) and |P| = p, |Q| = q.
Then we have by the clustering property of the spin model that
E(

∏

(1 − Zx )) − E( ∏ (1 − Zx ))E( ∏ (1 − Zx )) =

x∈P∪Q

x∈P

x∈Q

≤

∑

(−1)|A|+|B| E(

A⊂P,B⊂Q,
A∪B6=0/

∑

A⊂P,B⊂Q,
A∪B6=0/

∏

x∈A∪B


Zx ) − E( ∏ Zx )E( ∏ Zx ) ,
x∈A

C|A|+|B| φ (c|A|+|B| s) ≤ 2 p+qCp+q φ (s).

x∈B

(4.2)

Further, for any S ⊂ V such that s := s(S) = minx6=y∈S | − y + x|, |S| = l. Using the triangle inequality recursively, we
can conclude by induction that
|E(∏(1 − Zx )) − ∏ E(1 − Zx )| ≤ (l − 1)2l Cl φ (s).
x∈S

x∈S

Now let a∗ ∈ (0, ∞] be such that E(Zx ) = ρ (1) (O) = 1 − e−a∗ . Then combining with the above inequalities and that
Cl ≤ C∗ l al by exponential clustering assumption, we have that for S as above,
P(P(S) = 0) ≤ C∗ (l − 1)(2l a )l φ (s) + e−la∗ .

(4.3)

CLT for statistics of spin models

37
′

′

For any t large and any z ∈ Wn for n > 10t, we can find a subset S′ ⊂ Wt (z) ∩Wn such that |S′ | ∈ [A′t α , A′t α + 2] and
s(S′ ) = w∗ t β for 0 < α ′ , β < ∞. Such a subset exists because for all z ∈ Wn , there is at least a (non-intersecting) path
of length t in Wt (z) ∩Wn . Hence such a choice of S′ can be made with α , β = 1/2. Then, we further choose a subset
S ⊂ S′ such that |S| ∈ [A∗t α , A∗t α + 2] where 2α < bβ with b as in (1.10). As s(S) ≥ s(S′ ), we can derive from (4.3)
that
α
α
b bβ
sup P(Pn (Wt (z)) = 0) ≤ C∗∗ t α (2t aα )A∗t e−(w∗ ) t + e−a∗t ,
z∈Wn

where C∗∗ is a product of the various constants involved. Thus Pn satisfies(4.1) for all large enough n.
Acknowledgements
DY is thankful for the discussions with Matthew Wright which led to his interest in this question and especially
the applications to random cubical complexes. The authors are also thankful to numerous comments by anonymous
referees that has lead to an improved presentation.
References
1.
2.
3.
4.
5.
6.
7.
8.
9.
10.
11.
12.
13.
14.
15.
16.
17.
18.
19.
20.
21.
22.
23.
24.

Aldous, D., Lyons, R.: Processes on unimodular random networks. Electron. J. Probab. 12, 1454–1508 (2007)
Atkin, R.: An algebra for patterns on a complex, I. International Journal of Man-Machine Studies 6(3), 285 – 307 (1974)
Atkin, R.: An algebra for patterns on a complex, II. International Journal of Man-Machine Studies 8(5), 483 – 498 (1976)
Baccelli, F., Haji-Mirsadeghi, M.O., Khezeli, A.: Dynamics on unimodular random graphs. arXiv:1608.05940 (2016)
Baryshnikov, Y., Yukich, J., et al.: Gaussian limits for random measures in geometric probability. Ann. Appl. Prob. 15(1A), 213–253
(2005)
Beffara, V., Gayet, D.: Percolation of random nodal lines. arXiv:1605.08605 (2016)
Benjamini, I.: Coarse geometry and randomness, École d’Été de Probabilités de Saint-Flour, vol. 2100. Springer (2013)
Björklund, M., Gorodnik, A.: Central limit theorems for group actions which are exponentially mixing of all orders.
arXiv:1706.09167 (2017)
Błaszczyszyn, B.: Factorial moment expansion for stochastic systems. Stoch. Proc. Appl 56(2), 321–335 (1995)
Błaszczyszyn, B., Merzbach, E., Schmidt, V.: A note on expansion for functionals of spatial marked point processes. Stat. Probab.
Letters 36(3), 299–306 (1997)
Blaszczyszyn, B., Yogeshwaran, D., Yukich, J.E.: Limit theory for geometric statistics of point processes having fast decay of
correlations. arXiv:1606.03988 (2018)
Bobrowski, O., Kahle, M.: Topology of random geometric complexes: a survey. arXiv:1409.4734 (2017)
Boccaletti, S., Latora, V., Moreno, Y., Chavez, M., Hwang, D.U.: Complex networks: structure and dynamics. Physics Reports
424(4), 175–308 (2006)
Bolthausen, E., Cipriani, A., Kurt, N.: Exponential decay of covariances for the supercritical membrane model. Comm. Math. Phys.
353(3), 1217–1240 (2017)
Borcea, J., Brändén, P., Liggett, T.M.: Negative dependence and the geometry of polynomials. J. Amer. Math. Soc. 22(2), 521–567
(2009)
Bradley, R.: Equivalent mixing conditions for random fields. Ann. Probab 21(4), 1921–1926 (1993)
Bradley, R.: On quantiles and the central limit question for strongly mixing sequences. J. Theor. Probab. 10(2), 1921–1926 (1997)
Bradley, R.: Basic properties of strong mixing conditions : a survey and some open questions. Probab. Surveys 2, 107–144 (2005)
Bradley, R., Tone, C.: A central limit theorem for non-stationary strongly mixing random fields. J. Theor. Probab. 2, 107–144
(2015)
Bulinski, A., Spodarev, E.: Central limit theorems for weakly dependent random fields. In: Stochastic geometry, spatial statistics
and random fields, Lecture Notes in Math., vol. 2068, pp. 337–383. Springer, Heidelberg (2013)
Bulinski, A., Spodarev, E., Timmermann, F.: Central limit theorems for the excursion set volumes of weakly dependent random
fields. Bernoulli 18(1), 100–118 (2012)
Bulinski, A., Suquet, C.: Normal approximation for quasi-associated random fields. Statist. Probab. Lett. 54(2), 215–226 (2001)
d
Cohen, G., Conze, J.P.: Almost mixing of all orders and clt for some Zd actions on subgroups of F̥
p . arXiv:1609.06484 (2016)
Derriennic, Y.: Some aspects of recent works on limit theorems in ergodic theory with special emphasis on the central limit theorem.
Discrete and continuous dynamical systems 15(1) (2006)

38

Tulasi Ram Reddy, Sreekar Vadlamani, D. Yogeshwaran

25. Doukhan, P.: Properties and examples. In: Mixing, Lecture Notes in Statistics, vol. 85. Springer-Verlag, New York (1994)
26. Dousse, J., Féray, V.: Weighted dependency graphs and the Ising model. arXiv:1610.05082 (2016)
27. Duminil-Copin, H.: Graphical representations of lattice spin models. Lecture notes of Cours Peccot du Collège de France. Spartacus
(2015). URL http://www.ihes.fr/~ duminil/publi/2016Peccot.pdf
28. Edelsbrunner, H., Harer, J.: Computational Topology, An Introduction. American Mathematical Society, Providence, RI (2010)
29. Estrada, E., Rodriguez-Velazquez, J.A.: Complex networks as hypergraphs. arXiv:physics/0505137 (2005)
30. Féray, V.: Weighted dependency graphs. arXiv:1605.03836 (2016)
31. Forman, R.: A user’s guide to discrete Morse theory. Sém. Lothar. Combin. 48, Art. B48c, 35 (2002)
32. Franceschetti, M., Meester, R.: Random networks for communication: from statistical physics to information systems, vol. 24.
Cambridge University Press (2008)
33. Friedli, S., Velenik, Y.: Statistical mechanics of lattice systems: a concrete mathematical introduction. to be published by Cambridge
University Press (2017). URL http://www.unige.ch/math/folks/velenik/smbook/
34. Funaki, T.: Stochastic interface models. In: J. Picard (ed.) Lectures on Probability Theory and Statistics., pp. 103–274. Springer
(2005). Lecture Notes in Mathematics, vol. 1869
35. Giacomin, G.: Aspects of statistical mechanics of random surfaces.
IHP Lecture notes (2001).
URL
https://www.lpma-paris.fr/modsto/_media/users/giacomin/ihp.pdf
36. Goldstein, L., Wiroonsri, N.: Stein’s method for positively associated random variables with applications to the Ising and voter
models, bond percolation, and contact process. arXiv:1603.05322 (2016)
37. Göring, D., Klatt, M., Stegmann, C., Mecke, K.: Morphometric analysis in gamma-ray astronomy using Minkowski functionalssource detection via structure quantification. Astronomy & Astrophysics 555, A38 (2013)
38. Gray, S.B.: Local properties of binary images in two dimensions. IEEE Transactions on Computers 20(5), 551–561 (1971)
39. Grimmett, G.: Probability on graphs: random processes on graphs and lattices, vol. 1. Cambridge University Press (2010)
40. Gromov, M.: Groups of polynomial growth and expanding maps. Appendix by J. Tits. Publ. Math. de l’I.H.E.S. 53, 53–78 (1981)
41. Gross, L.: Decay of correlations in classical lattice models at high temperature. Comm. Math. Phys. 68(1), 9–27 (1979)
42. Grote, J., Thäle, C.: Gaussian polytopes: a cumulant-based approach. arXiv:1602.06148 (2016)
43. Haenggi, M.: Interference in lattice networks. arXiv:1004.0027 (2010)
44. Hegerfeldt, G.C.: Noncommutative analogs of probabilistic notions and results. J. Funct. Anal. 64(3), 436–456 (1985)
45. Heinrich, L.: Asymptotic methods in statistics of random point processes. In: Stochastic geometry, spatial statistics and random
fields, Lecture Notes in Math., vol. 2068, pp. 115–150. Springer, Heidelberg (2013)
46. Hilfer, R.: Local porosity theory and stochastic reconstruction for porous media. In: Statistical Physics and Spatial Statistics, pp.
203–241. Springer (2000)
47. Hiraoka, Y., Tsunoda, K.: Limit theorems on random cubical homology. arXiv:1612.08485 (2016)
48. Holley, R.A., Stroock, D.W.: Applications of the stochastic Ising model to the Gibbs states. Comm. Math. Phys. 48(3), 249–265
(1976)
49. Hough, J.B., Krishnapur, M., Peres, Y., Virág, B.: Zeros of Gaussian analytic functions and determinantal point processes, University Lecture Series, vol. 51. American Mathematical Society, Providence, RI (2009)
50. Ioffe, D., Velenik, Y.: A note on the decay of correlations under δ -pinning. Probab. Theory Related Fields 116(3), 379–389 (2000)
51. Janson, S.: Normal convergence by higher semi-invariants with applications to sums of dependent random variables and random
graphs. Ann. Prob. 16(1), 305–312 (1988)
52. Kaczynski, T., Mischaikow, K., Mrozek, M.: Computational homology. New York (2004)
53. Kahle, M.: Topology of random simplicial complexes: a survey. AMS Contemp. Math 620, 201–222 (2014)
54. Klamt, S., Haus, U.U., Theis, F.: Hypergraphs and cellular networks. PLoS computational biology 5(5), e1000,385 (2009)
55. Klatt, M.A.: Morphometry of random spatial structures in physics.
Ph.D. thesis (2016).
URL
https://opus4.kobv.de/opus4-fau/frontdoor/index/index/docId/7654 .
Friedrich-Alexander-Universität
Erlangen-Nürnberg
56. Klette, R., Rosenfeld, A.: Digital geometry: geometric methods for digital picture analysis. Elsevier (2004)
57. Kong, T.Y., Rosenfeld, A.: Digital topology: introduction and survey. Computer Vision, Graphics, and Image Processing 48(3),
357–393 (1989)
58. Kopper, C., Magnen, J., Rivasseau, V.: Mass generation in the large N Gross-Neveu-model. Communications in Mathematical
Physics 169(1), 121–180 (1995)
59. Kraetzl, M., Laubenbacher, R., Gaston, M.E.: Combinatorial and algebraic approaches to network analysis (2001). DSTO Internal
Report
60. Krokowski, K., Thäle, C., et al.: Multivariate central limit theorems for rademacher functionals with applications. Elec. J. Prob. 22
(2017)
61. Künsch, H.: Decay of correlations under Dobrushin’s uniqueness condition and its applications. Comm. Math. Phy. 84(2), 207–222
62. de La Harpe, P.: Topics in geometric group theory. University of Chicago Press (2000)
63. Lachieze-Rey, R., Schulte, M., Yukich, J.E.: Normal approximation for stabilizing functionals. arXiv:1702.00726 (2017)

CLT for statistics of spin models

39

64. Lyons, R.: Determinantal probability measures. Publ. Math. Inst. Hautes Études Sci. (98), 167–212 (2003)
65. Lyons, R., Peres, Y.: Probability on trees and networks, Cambridge Series in Statistical and Probabilistic Mathematics, vol. 42.
Cambridge University Press, New York (2016)
66. Lyons, R., Steif, J.E.: Stationary determinantal processes: phase multiplicity, Bernoullicity, entropy, and domination. Duke Math.
J. 120(3), 515–575 (2003)
67. Malyshev, V.A.: The central limit theorem for Gibbsian random fields. Soviet Mathematics Doklady 16, 1141–1145 (1975)
68. Martin, P.A., Yalcin, T.: The charge fluctuations in classical Coulomb systems. J. Stat. Phys 22(4), 435–463 (1980)
69. Michoel, T., Nachtergaele, B.: Alignment and integration of complex networks by hypergraph-based spectral clustering. Physical
Review E 86(5), 056,111 (2012)
70. Munkres, J.: Elements of algebraic topology. Addison-Wesley (1984)
71. Nazarov, F., Sodin, M.: Correlation functions for random complex zeroes: strong clustering and local universality. Comm. Math.
Phys. 310(1), 75–98 (2012)
72. Pansu, P.: Croissance des boules et des géodésiques fermées dans les nilvariétés. Ergodic Theory and Dynamical Systems 3(3),
415–445 (1983)
73. Peccati, G., Taqqu, M.S.: Wiener chaos: moments, cumulants and diagrams, Bocconi & Springer Series, vol. 1. Springer, Milan;
Bocconi University Press, Milan (2011)
74. Peligrad, M.: Maximum of partial sums and in invariance principle for a class of weak dependent random variables. Proc. AMS
126(4), 1181–1189 (1998)
75. Penrose, M.: Random geometric graphs, Oxford Studies in Probability, vol. 5. Oxford University Press, Oxford (2003)
76. Penrose, M.D.: A central limit theorem with applications to percolation, epidemics and Boolean models. Ann. Probab. 29(4),
1515–1546 (2001)
77. Penrose, M.D., Yukich, J.E.: Limit theory for point processes in manifolds. Ann. Appl. Prob. 23(6), 2161–2211 (2013)
78. Penrose, O., Lebowitz, J.L.: On the exponential decay of correlation functions. Comm. Math. Phys. 39(3), 165–184 (1974)
79. Pete, G.: Probability and geometry on groups.
Lecture notes for a graduate course (2017).
URL
http://math.bme.hu/~ gabor/PGG.pdf
80. Roe, J.: Lectures on coarse geometry, vol. 31. American Mathematical Soc. (2003)
81. Saha, P.K., Strand, R., Borgefors, G.: Digital topology and geometry in medical imaging: a survey. IEEE transactions on medical
imaging 34(9), 1940–1964 (2015)
82. Saulis, L., Statulevicius, V.: Limit theorems for large deviations. Kluwer Academic, Dordrecht (1991)
83. Schladitz, K., Ohser, J., Nagel, W.: Measurement of intrinsic volumes of sets observed on lattices. Discrete Geometry for Computer
Imagery pp. 247–258 (2006)
84. Schneider, R., Weil, W.: Stochastic and integral geometry. Probability and its Applications. Springer-Verlag, Berlin (2008)
85. Schonmann, R.H.: Theorems and conjectures on the droplet-driven relaxation of stochastic Ising models. In: Probability and Phase
Transition, pp. 265–301. Springer (1994)
86. Spanier, E.H.: Algebraic topology. McGaw-Hill Book Co., New York (1966)
87. Sunklodas, J.: Approximation of distributions of sums of weakly dependent random variables by the normal distribution, pp. 113–
165. Springer-Verlag Berlin Heidelberg (1991)
88. Svane, A.M.: Valuations in image analysis, pp. 435–454. Springer International Publishing, Cham (2017)
89. Velenik, Y.: Localization and delocalization of random interfaces. Probab. Surv 3, 112–169 (2006)
90. Werman, M., Wright, M.: Intrinsic volumes of random cubical complexes. Discrete Comput. Geom. 56, 93–113 (2016)
91. Yukich, J.: Limit theorems in discrete stochastic geometry. In: Stochastic geometry, spatial statistics and random fields, Lecture
Notes in Math., vol. 2068, pp. 239–275. Springer, Heidelberg (2013)

