1. Introduction
Point processes form the main building block of stochastic geometry. They serve to construct a broad spectrum of geometric models used for analysing spatial data, for instance in material science, particle physics, telecommunications, and biology. The theory of point processes is now available for very abstract topological settings (see e.g. [Reference Daley and Vere-Jones5, Reference Schneider and Weil26]), but in the current paper, we restrict ourselves to point configurations randomly sprinkled in $\mathbb{R}^d$ . This restriction is natural and still allows us to build many random geometric structures by using certain connections among points (e.g. random graphs or tessellations), by adding certain marks to each point (a number, set, etc.), or by a combination of these methods.
A stationary (translation-invariant) point process $\Gamma$ in $\mathbb{R}^d$ is called hyperuniform if it exhibits small density fluctuations. Formally, $\Gamma$ is hyperuniform if the number of points of the process in a bounded domain $N_\Lambda\;:\!=\; N_\Lambda(\Gamma)$ fluctuates at a lower order then the volume of the set, i.e.,
where by $\Lambda \nearrow \mathbb{R}^d$ we understand any sequence of increasing sets $\Lambda_1 \subseteq \Lambda_2 \subseteq \ldots$ tending to the whole space $\mathbb{R}^d$ . Historically, hyperuniformity is related to the study of the compressibility of matter. In statistical mechanics, compressibility is a measure of the relative volume change of a fluid or solid in response to a pressure, which can be understood as the relative size of fluctuations in particle density. Nonetheless, knowing the order of the variance is of interdisciplinary interest. In material science, the concept of hyperuniformity enables the characterization of naturally organized structures such as crystals and quasicrystals (see [Reference Björklund and Hartnick2, 20]). The quantitative characterization of fluctuations in the number of particles has a long history in statistical physics as well. Under certain constraints on the two-point correlation function (in the physics literature referred to as sum rules), the authors of [Reference Martin and Yalcin18] (with an extension in [Reference Lebowitz17]) showed for infinite classical systems of particles with long-range interactions that the variance of $N_\Lambda$ should increase as the surface area. A rigorous result for one-dimensional Coulomb systems is shown in [Reference Kunz15]. In the theory of random matrices, one-dimensional point patterns associated with the eigenvalues have been characterized by their density fluctuations (see e.g. [Reference Metha19]). The measurement of galaxy density fluctuations is a standard approach to studying the structure of the universe (see [Reference Peebles21]). The concepts of hyperuniformity and density fluctuations have been identified across many other areas of fundamental science, such as computer science, number theory, and the biological sciences. Yet the theoretical understanding of such systems is still limited. The first attempt to rigorously handle the concept of hyperuniformity in the physics of matter was established in the seminal paper [Reference Torquato and Stillinger29] and the subsequent papers of the authors. The current state of the art is summarized in the survey [Reference Torquato28].
It seems to be a fundamental question of great practical interest to determine whether a given point process is hyperuniform or not. Unsurprisingly, this question has become a popular topic of study for researchers in stochastic geometry and related fields. In particular, the study of Coulomb and Riesz gases in the context of the theory of point processes is of considerable interest in both mathematics and physics. Some answers have been given for $d=2,3$ in [Reference Chatterjee4, Reference Leblé16]. By estimating the structure factor, the authors of [Reference Hawat, Gautier, Bardenet and Lachièze-Rey12, Reference Klatt, Last and Henze14] provide tests of hyperuniformity based on point samples. It is also of great value in stochastic geometry and spatial statistics to be able not only to prove hyperuniformity, but also to verify when it does not hold. Nondegeneracy of the asymptotic variance (1.1) is one of the key assumptions for geometric central limit theorems (as in [Reference Błaszczyszyn, Yogeshwaran and Yukich3, Reference Penrose and Yukich22, Reference Schreiber and Yukich27]).
The question of hyperuniformity is immediate for some special point processes, especially if the exact distribution of $N_{\Lambda}$ is tractable. This is the case for a stationary Poisson point process with intensity $\lambda>0$ , which obviously is not hyperuniform. Here, the number of points in any set is independent of the outside configurations. In general, it is unclear what happens if interactions between the points inside the set and outside are introduced into the model. A standard approach to generating point patterns with interactions among points is to consider a Gibbs modification of some underlying measure, typically a Poisson point process through energy. The simplest model takes into account only interactions between pairs of points; this is known as a Gibbs point process with pairwise interactions. Interactions among k-tuples of points and other more complicated types of interactions are also widely studied in stochastic geometry and spatial statistics.
In this paper, we show that short-range Gibbs point processes are not hyperuniform. By a short-range process, we simply mean a point process such that interactions among points weaken with the distance. This phenomenon can be interpreted through the Papangelou intensity. A point process $\Gamma$ has a Papangelou intensity $\lambda^*$ if for any non-negative function f we can write
Intuitively, we interpret $\lambda^*(x,\gamma)\mathrm{d}\, x$ as the conditional probability of observing a point in the infinitesimally small neighbourhood of x given that $\Gamma$ agrees with a configuration $\gamma$ outside this neighbourhood. For example, if $\varphi\;:\; \mathbb{R}^d\times \mathbb{R}^d\to \mathbb{R} \cup \{+\infty\}$ measures the interaction between two points, then a Gibbs point process with pairwise interactions has a Papangelou intensity of the form
where $z>0, \beta \geq 0$ are usually called activity and inverse temperature. Now, a point process is short-range if
fast enough as $\|x-y\| \to \infty$ for almost all configurations $\gamma$ . The interpretation is straightforward. Any point y in the configuration $\gamma$ plays only a negligible role when a new point x is introduced far away from y. For a Gibbs point process with pair potential $\varphi$ , (1.2) simply translates to the condition that $\varphi(x,y)\to 0$ fast enough.
In this paper, we prove a kind of compressibility of any scale in the bulk of the interacting particle systems. Our main theorem is formulated for infinite-volume Gibbs point processes with a general Papangelou intensity $\lambda^*$ . We claim that if the Papangelou intensity satisfies some moment conditions, it is enough to verify that
for some $\alpha\geq 1$ to ensure that $\Gamma$ is not hyperuniform. Specially for pairwise interactions with potential $\Phi$ , (1.3) translates into the integrability condition
which, combined with superstability and lower-regularity as understood in [Reference Ruelle25], generates a non-hyperuniform point process. In the context of pair potentials, the authors of [Reference Ginibre11, Reference Ruelle25] also provide some variance estimates, but only for the finite-volume Gibbs measures. A more general question than the initial one is addressed in [Reference Xia and Yukich30], where the authors study volume-order fluctuations of score functions (including the number of points). The setting is for short-range Gibbs processes that are dominated by a Poisson point process. It corresponds to purely repulsive interactions, which can be translated to the Papangelou intensity by assuming $\lambda^*\geq C>0$ everywhere. Such a setting is also discussed here in Section 3.2. For our main theorem, we consider more general interactions that can be arbitrary within a short range and vanishing at a long range. The theorem also covers processes that cannot be coupled with a Poisson point process. The techniques in the two papers are completely different. To show our results, we use an approach based on the Georgii–Nguyen–Zessin (GNZ) formalism, rather than the algorithmic construction of the Gibbs point processes used in [Reference Xia and Yukich30].
The organization of the paper is as follows. In Section 2, we introduce all the necessary notation covering the theory of point processes with marks, and we present two versions of our main result. First, we state a non-hyperuniformity result for Gibbs point processes with pairwise interactions or interactions that become deterministic for distant points. Next, we state a theorem for Gibbs point processes such that the interactions are random everywhere, but also weak at long distances, which is suitable for a large variety of geometrically based interactions. Section 3 provides some further investigation and a couple of hints for verifying Assumption 1.3. It is followed by Section 4, where we provide a comprehensive application of our main results, discussing first pairwise interactions and then more general geometrical interactions. Finally, in Section 5, we give rigorous proofs of our main statements.
2. Main results
2.1. Notation
To describe marked Gibbs point processes in the space $\mathbb{R}^d$ , we proceed as follows. Let $\mathcal{B}$ be the Borel $\sigma$ -algebra on $\mathbb{R}^d$ and $\lambda$ the Lebesgue measure on $(\mathbb{R}^d, \mathcal{B})$ . Furthermore, we introduce a complete, separable mark space $\mathbb{M}$ , equipped with the associated Borel $\sigma$ -algebra $\mathcal{B}_{\mathbb{M}}$ and a finite measure $\lambda_{\mathbb{M}}$ .
We consider point configurations on the space $E\;:\!=\; \mathbb{R}^d \times \mathbb{M}$ with $\sigma$ -field $\mathcal{E}\;:\!=\; \mathcal{B}\times \mathcal{B}_{\mathbb{M}}$ , where each point in $\mathbb{R}^d$ is associated with a mark belonging to $\mathbb{M}$ . The reference measure on E is the product measure $\nu\;:\!=\; \lambda\otimes \lambda_{\mathbb{M}}$ . The marks are mutually independent random variables whose distribution $\mathbb{Q}_{\mathbb{M}}$ does not depend on the location.
Let $\textbf{N}$ denote the set of all locally finite configurations in E, i.e.,
Moreover, we denote by $\textbf{N}_f$ the subset of $\textbf{N}$ consisting of finite configurations. We endow $\textbf{N}$ with $\mathcal{N}$ , which is the smallest $\sigma$ -field such that all projections $\gamma\mapsto \gamma\cap B$ are measurable for all $B \in \mathcal{E}$ . For a point configuration $\gamma \in \textbf{N}$ and a fixed set $\Lambda \subset \mathbb{R}^d$ , we denote by $\gamma_\Lambda$ the restriction of $\gamma$ to the set $\Lambda\times \mathbb{M}$ , i.e. $\gamma_\Lambda \;:\!=\; \gamma \cap (\Lambda \times \mathbb{M})$ .
By a point process we mean a probability measure $\Gamma$ on $(\textbf{N},\mathcal{N})$ . For $u \in \mathbb{R}^d$ , we let $\tau_u\;:\;(x,m)\mapsto (x+u,m)$ be the shift in the position coordinate. If $\gamma=\{\textbf{x}_1,\textbf{x}_2,\ldots\}\in \textbf{N}$ , we write $\tau_u \gamma = \{\tau_u \textbf{x}_1,\tau_u\textbf{x}_2,\ldots\}$ . Moreover, we call a point process stationary if its distribution is invariant with respect to $\tau_u$ , i.e., $\Gamma \stackrel{D}{=}\tau_u\Gamma$ for all $u \in \mathbb{R}^d$ .
For a bounded set $\Lambda\subset \mathbb{R}^d$ and $\gamma \in \textbf{N}$ , we denote the number of points of $\gamma$ occurring in $\Lambda$ by $N_\Lambda\;:\!=\;N_\Lambda (\gamma)=\sum_{\textbf{x}\in \gamma}\textbf{1}_{\Lambda\times\mathbb{M}}(\textbf{x})$ . If $\Gamma$ is a marked point process, then $N_\Lambda(\Gamma)$ is a random variable with values in $\mathbb{N}\cup\{0\}$ .
2.2. Gibbs point process
Definition 2.1. Let $\lambda^*\;:\;E\times\textbf{N} \to \mathbb{R}$ be some measurable function. We call $\Gamma$ a Gibbs point process associated with the Papangelou intensity $\lambda^*$ if for all positive measurable $f\;:\;E\times\textbf{N}\to \mathbb{R}$ it solves the Georgii–Nguyen–Zessin equations
Remark 2.1. (The form of the Papangelou intensity.) Usually, the Papangelou intensity is given in the form
where $h\;:\;E\times\textbf{N} \to \mathbb{R}$ is called local energy, $z>0$ is the activity parameter, and $\beta\geq 0$ is the inverse temperature.
Here we consider only stationary Gibbs point processes. In this case, the related Papangelou intensity $\lambda^*$ is necessarily translation-invariant simultaneously with respect to both coordinates, meaning that
Remark 2.2. The existence of a point process satisfying Definition 2.1 is not discussed in the present paper; it is always implicitly assumed that we are given a well-defined Gibbs point process. However, note that (1.2) is related to the condition for existence as in [Reference Dereudre and Vasseur9]. We mention some relevant existence results in the course of the text. In addition, the point process $\Gamma$ may not be uniquely determined by (2.1). If we make a claim about a Gibbs process with Papangelou intensity $\lambda^*$ , we mean that it holds for all processes corresponding to this Papangelou intensity.
2.3. Main results
Our first result concerns point processes whose interactions are short-range and for which the ratio (1.2) has an upper bound which is uniform in $\textbf{N}$ . An important class of short-range Gibbs point processes is given by those determined by an integrable pair potential (see Section 4.1 for details). For multibody potentials, we replace the integrability condition by a suitable deterministic bound, which is integrable in the sense of [Reference Ginibre11, Reference Ruelle25], except that we do not force the potentials to be integrable around the origin.
Theorem 2.1. Let $\Gamma$ be a stationary Gibbs point process on $\mathbb{R}^d$ with Papangelou intensity $\lambda^*$ , and let the following assumptions be satisfied:
-
(a1) For all $m \in \mathbb{M}$ , $\mathbb{E}\lambda^*((0,m),\Gamma)^2<\infty$ .
-
(a2) There exist a function $\phi\;:\;\mathbb{R}^d\to \mathbb{R}$ and some $\delta\geq 0$ such that $\int_{\mathbb{R}^d\setminus B(0,\delta)}\phi(x)\mathrm{d}\, x<\infty$ and, for any $\|x-y\|>\delta$ ,
\begin{align*} \left|1-\frac{\lambda^*(\textbf{x},\gamma\cup \{\textbf{y}\})}{\lambda^*(\textbf{x},\gamma)}\right|\leq\phi(x-y). \end{align*}
Then there is a constant $\mathcal{C}_{nhyp}>0$ , not depending on $\Lambda$ , such that
In particular, $\Gamma$ is not hyperuniform.
Theorem 2.1 applies to point processes where the interaction among distant points is deterministically bounded. This, however, is often too restrictive for processes with a random range of interaction (e.g. processes with energies based on Voronoi tessellations). For such processes, we provide another result.
Theorem 2.2. Let $\Gamma$ be a stationary Gibbs point process on $\mathbb{R}^d$ with Papangelou intensity $\lambda^*$ . Assume that there are some $\alpha_1, \alpha_2 > 1$ , $\frac{1}{\alpha_1}+\frac{1}{\alpha_2} = 1$ , and $\delta\geq 0$ such that the following hold:
-
(A1) $ \int_{\mathbb{M}} \mathbb{E} |\lambda^*((0,m),\Gamma)|^{2\alpha_1}\lambda_{\mathbb{M}}(\mathrm{d}\, m) < +\infty$ , and
-
(A2) $\int_{\mathbb{R}^d\setminus B(0,\delta)}\left(\int_{\mathbb{M}^2}\mathbb{E} \left|1-\frac{\lambda^*((0,m_1),\Gamma \cup \{(y,m_2)\})}{\lambda^*((0,m_1),\Gamma)}\right|^{\alpha_2}\mathrm{d}\,\lambda_{\mathbb{M}}(m_1,m_2)\right)^{1/\alpha_2} \mathrm{d}\, y <\infty.$
Then there exists $\mathcal{C}_{nhyp}>0$ , not depending on $\Lambda$ , such that
In particular, $\Gamma$ is not hyperuniform.
The proofs of Theorems 2.2 and 2.2 are postponed to Section 5. The constant $\mathcal{C}_{nhyp}$ can be given in a closed form directly from the proof. We provide some lower bounds for these constants in special cases in Section 4.
Remark 2.3. (Assumptions (A1)–(A2).) Note that $\mathbb{E}\lambda^*(\textbf{0},\Gamma) \;:\!=\; \lambda$ is the intensity of $\Gamma$ . Loosely speaking, Assumption (A1) of Theorem 2.2 prevents the point process $\Gamma$ from having too many points in a unit window. In other words, it forces some stability assumptions (e.g. superstability for pair potentials). On the other hand, Assumption (A2) states that the interactions among points of $\Gamma$ become negligible with the distance. Again, we do not force the interactions among points close to each other to be bounded. A similar interpretation applies for Assumptions (a1) and (a2) of Theorem 2.1.
Remark 2.4. (Shape of $\Lambda$ .) Note that the proof and the lower bounds (2.2)–(2.7) do not depend on the shape of the window $\Lambda$ .
Remark 2.5. (Unmarked case.) In the unmarked case, one can choose the mark space $\mathbb{M}$ with just one atom m and set $\lambda_{\mathbb{M}}(m)=1$ . The Papangelou intensity obviously does not depend on the choice of m, so we may set $\Gamma'= \{x_i;\; \,(x_i,m_i)\in \Gamma\}$ to be the unmarked point process and write $\lambda^*((x,m),\Gamma)=\lambda^*(x,\Gamma')$ almost surely (a.s.). Then the assumptions of Theorem 2.2 become
-
(A1) $ \mathbb{E} |\lambda^* (0, \Gamma')|^{2\alpha_1}< +\infty$ ,
-
(A2) $\int_{\mathbb{R}^d\setminus B(0,\delta)}\left(\mathbb{E} \left|1-\frac{\lambda^*(0,\Gamma' \cup \{y\})}{\lambda^*(0,\Gamma')}\right|^{\alpha_2}\right)^{1/\alpha_2} \mathrm{d}\, {y} <\infty.$
3. Conditions for Assumptions (a1)–(a2) (resp. (A1)–(A2))
In this section, we provide a short cookbook of conditions that can be imposed on the Papangelou intensity $\lambda^*$ to guarantee the validity of Assumptions (a1)–(a2) of Theorem 2.1 (resp. Assumptions (A1)–(A2) of Theorem 2.2).
3.1. Stability and range of interaction
Definition 3.1. (Local stability.) The Papangelou intensity $\lambda^*$ is called
-
locally stable from above if there is a constant $C_1 <\infty$ such that $\lambda^*(\textbf{x},\gamma) \leq C_1$ for all $\textbf{x}\in E$ and $\gamma \in \textbf{N}$ ,
-
locally stable from below if there is a constant $C_2> 0$ such that $\lambda^*(\textbf{x},\gamma) \geq C_2$ for all $\textbf{x}\in E$ and $\gamma \in \textbf{N}$ , and
-
double locally stable if it is simultaneously locally stable from above and from below.
Note that a point process whose Papangelou intensity is locally stable from above trivially satisfies Assumption (a1) of Theorem 2.1 (resp. (A1) of Theorem 2.2). Yet this is a very frequent assumption in the literature.
Definition 3.2. (Range of interaction.) Let $\Gamma$ be a Gibbs point process with Papangelou intensity $\lambda^*$ . Then $\Gamma$ has
-
finite range of interaction if there exists $R>0$ such that for all $\gamma \in \mathcal{N}$
\begin{align*}\lambda^*(\textbf{x},\gamma)= \lambda^*(\textbf{x}, \gamma \cap B(x,R)\times \mathbb{M}), \quad \forall \textbf{x}\;:\!=\;(x,m)\in E;\end{align*} -
random finite range of interaction if for all $\textbf{x}\;:\!=\;(x,m) \in E$ there is an a.s. finite random variable $R_{\textbf{x}}\;:\!=\; R_{\textbf{x}}(\gamma)$ such that
\begin{align*}\lambda^*(\textbf{x},\gamma) = \lambda^*(\textbf{x}, \gamma \cap B(x,R_{\textbf{x}})\times\mathbb{M}), \quad \text{ for $\Gamma$-almost all }\Gamma;\;\end{align*} -
decreasing range of interaction if
(3.1) \begin{align}R_\textbf{x}(\gamma) \geq R_\textbf{x}(\gamma \cup \{\textbf{y}\})\quad \forall \textbf{x}, \textbf{y} \in E, \gamma \in \textbf{N}.\end{align}
Remark 3.1. Note that (3.1) is in fact a natural condition arising from many geometric models, including Voronoi tessellations, Delaunay triangulations, and k-nearest-neighbours graphs.
Proposition 3.1. Let $\Gamma$ be a Gibbs point process associated with Papangelou intensity $\lambda^*$ . Then the following hold:
-
(i) A Gibbs point process whose Papangelou intensity is locally stable from above satisfies Assumption (a1) of Theorem 2.1 (resp. (A1) of Theorem 2.2 with any $\alpha_1>0$ ).
-
(ii) Suppose $\Gamma$ has a random finite range of interaction and there exists a function $f\;:\;\mathbb{R}^+\to \mathbb{R}$ such that $\lambda^*(\textbf{x},\gamma)\leq f(R_{\textbf{x}}(\gamma))$ for all $\textbf{x} \in E$ and $\Gamma$ -almost all $\gamma$ . If $\mathbb{E} f(R_\textbf{0})^{2\alpha_1}<\infty$ then $\Gamma$ satisfies Assumption (A1) of Theorem 2.2.
-
(iii) If $\Gamma$ has a finite range of interaction, then it satisfies Assumption (a2) of Theorem 2.1 and Assumption (A2) of Theorem 2.2 for any $\alpha_2>1$ .
-
(iv) If $\lambda^*$ is double locally stable and $\Gamma$ has a decreasing random finite range of interaction such that $\mathbb{E} R^\alpha_{\textbf{0}} <\infty$ for some $\alpha> d$ , then Assumption (A2) of Theorem 2.2 is satisfied for any $\alpha_2 \leq \alpha/d$ .
Proof.
-
(i) Since $\lambda^*$ is uniformly bounded, any moment is finite.
-
(ii) For any $m \in \mathbb{M}$ , we have
\begin{align*}\mathbb{E}\lambda^*((0,m),\Gamma)^{2\alpha_1}\leq\mathbb{E} f(R_{(0,m)})^{2\alpha_1}<\infty. \end{align*}Since $\lambda_{\mathbb{M}}$ was assumed to be finite, (A1) is satisfied for any mark distribution. -
(iii) Consider $\delta = R$ and $\phi \equiv 0$ in Theorem 2.1, resp. $\delta = R$ in Theorem 2.2.
-
(iv) Take $\alpha_2 \in (1,\alpha/d)$ . By using the Markov inequality, we arrive at
\begin{align*} \int_{\mathbb{R}^d}&\left(\int_{\mathbb{M}^2}\mathbb{E} \left|1-\frac{\lambda^*((0,m_1),\Gamma \cup \{(y,m_2)\})}{\lambda^*((0,m_1),\Gamma)}\right|^{\alpha_2}\mathrm{d}\, \lambda^2_{\mathbb{M}}(m_1,m_2) \right)^{1/\alpha_2} \mathrm{d}\, y\\ &\leq \left(1+\frac{C_2}{C_1}\right)\int_{\mathbb{R}^d} \left[ \int_{\mathbb{M}^2}\mathbb{P} \left(|y|\leq R_{(0,m_1)}(\Gamma\cup\{(y,m_2)\})\right)\mathrm{d}\, \lambda^2_{\mathbb{M}}(m_1,m_2)\right]^{1/\alpha_2}\mathrm{d}\, y\\ &\leq \left(1+\frac{C_2}{C_1}\right)\lambda_{\mathbb{M}}(\mathbb{M})^{1/\alpha_2}\int_{\mathbb{R}^d} \left[ \int_{\mathbb{M}}\mathbb{P} \left(|y|\leq R_{(0,m_1)}(\Gamma)\right)\mathrm{d}\, \lambda_{\mathbb{M}}(m_1)\right]^{1/\alpha_2}\mathrm{d}\, y\\ &\leq \left(1+\frac{C_2}{C_1}\right)\lambda_{\mathbb{M}}(\mathbb{M})^{1/\alpha_2}\int_{\mathbb{R}^d} \frac{1}{|y|^{\alpha/\alpha_2}} \mathrm{d}\, y \left( \int_{\mathbb{M}}\mathbb{E} R_{(0,m_1)}^{\alpha}(\Gamma)\mathrm{d}\, \lambda_{\mathbb{M}}(m_1)\right)^{1/\alpha_2}<\infty.\end{align*}
In conclusion, the task of verifying the assumptions of Theorems 2.1 and 2.2 often translates to estimating the moments, resp. the tail probabilities of the radius of interactions. In more than a few situations, one can benefit from the literature on stabilization and stochastic comparison to the Poisson point process or other random structures. The latter is described in the following section.
3.2. Stochastic comparison tools
In this section, we explore tools to estimate the moments $\mathbb{E} R_{\textbf{0}}(\Gamma)^\alpha$ , $\alpha>0$ , of the range of interaction needed in the previous section. Usually, for an infinite-volume Gibbs point process, this is not a straightforward task, since we do not possess the local distribution of the number of points. We use stochastic comparison with some random object which is easier to handle.
For this purpose, we consider the usual order on $\textbf{N}$ . For $\gamma_1, \gamma_2 \in \textbf{N}$ , we write $\gamma_1\leq \gamma_2$ if $\gamma_1(B)\leq \gamma_2(B)$ for all $B \in \mathcal{E}$ . In the language of point sets, this means that $\gamma_1$ has fewer points than $\gamma_2$ . A function $f\;:\;\textbf{N}\to \mathbb{R}$ is called increasing if $f(\gamma_1)\leq f(\gamma_2)$ whenever $\gamma_1\leq \gamma_2$ . It is decreasing if $-f$ is increasing.
Definition 3.3. We say that a point process $\Gamma_1$ is stochastically dominated by a point process $\Gamma_2$ (and we write $\Gamma_1\ll\Gamma_2)$ if $\int h \mathrm{d}\, \Gamma_1 \leq \int h \mathrm{d}\, \Gamma_2$ for all increasing functions h. Vice versa, we say that $\Gamma_1$ stochastically minorates $\Gamma_2$ .
By the famous Strassen theorem, $\Gamma_1\ll\Gamma_2$ if and only if there is a coupling of $\Gamma_1$ and $\Gamma_2$ that is supported in the set $\{(\gamma_1,\gamma_2), \gamma_1\leq \gamma_2\}$ .
First, we recall the ‘Poisson sandwich inequality’ from [Reference Georgii and Küneth10] which stochastically connects a Gibbs point process with a stationary Poisson point process.
Proposition 3.2. Let $\Gamma$ be a Gibbs point process associated with Papangelou intensity $\lambda^*$ , and let $\Pi_\rho$ denote a stationary Poisson point process with intensity $\rho$ . Then the following hold:
-
(i) If $\lambda^*$ is locally stable from below, then $\Pi_{{C_1}}\ll\Gamma$ .
-
(ii) If $\lambda^*$ is locally stable from above, then $\Gamma\ll\Pi_{{C_2}}$ ,
Here, $C_1$ and $C_2$ are the constants from Definition 3.1.
Corollary 3.1. Let $\Gamma$ be associated with the Papangelou intensity $\lambda^*$ and assume that the corresponding range of interaction $R_{\textbf{x}}$ is decreasing. If $\lambda^*$ is locally stable from below, then
where $\Pi_{C_1}$ is as in Proposition 3.2.
Remark 3.2. (Increasing range of interaction.) If the range of interaction is increasing with respect to the order on $\textbf{N}$ , then Corollary 3.1 provides the inequality in the opposite direction. In this situation, the local stability of $\lambda^*$ from above produces similar upper bounds for the tail probabilities of $R_\textbf{x}(\Gamma)$ .
Consequently, for $\Gamma$ as in Corollary 3.1 and $f\;:\;\mathbb{R} \to \mathbb{R}$ an increasing function, we may verify the assumptions of Proposition 3.1(ii) and Proposition 3.1(iv) by computing
if $f\;:\;\mathbb{R} \to \mathbb{R}$ is some increasing function.
Theorems 2.1 and 2.2 provide non-hyperuniformity results for Gibbs processes with interactions that are weak at long distances. As mentioned before, we assume nothing about the interactions among close points. Such points can generate an unboundedly large amount of energy (imagine Coulomb interaction), and therefore the value $\lambda^*(\textbf{x},\gamma)$ can be arbitrarily close to zero when $\textbf{x}$ is at a small distance from $\gamma$ , i.e. when $d(\textbf{x},\gamma)\;:\!=\; \inf\{d(\textbf{x},\textbf{y});\; \textbf{y} \in \gamma\}$ is small, where $d(\textbf{x},\textbf{y})$ is some distance between two marked points. This situation does not allow one to use minoration by a Poisson point process. However, in this situation, we are able to construct a coupling with a Bernoulli field.
For this purpose, we introduce the mapping $I^{s}\;:\; \textbf{N} \to \{0,1\}^{\mathbb{Z}^d}$ as follows. We split the space $\mathbb{R}^d$ into a collection of disjoint cubes $\mathcal{D}\;:\!=\;\{\mathcal{D}_i;\; i\in \mathbb{Z}^d\}$ of a common side length $s>0$ such that there exists $k \in \mathbb{Z}^d$ with the origin being one of the vertices of $\mathcal{D}_k$ . Then, we define
In fact, in many situations, we do not need to know the exact positions of the points of the process to estimate the range of interaction. The only necessary information is often that there is at least one point in a given region. That is exactly the meaning of $I^{\delta}$ .
In order to formulate our next result, we need some partial order on $\{0,1\}^{Z^d}$ . We write that $l_1 \leq l_2$ if $l_1,l_2 \in \{0,1\}^{\mathbb{Z}^d}$ are such that $l_{1,i}=1$ implies $l_{2,i}=1$ for all $i \in \mathbb{Z}^d$ . As usual, we say that $R\;:\;\{0,1\}^{Z^d}\to \mathbb{R}$ is increasing if $R(l_1)\leq R(l_2)$ whenever $l_1\leq l_2$ and decreasing if $-R$ is increasing. We write $X\ll Y$ for two random variables X, Y with values in $\{0,1\}^{\otimes \mathbb{Z}^d}$ if $\mathbb{P}(X_{i}=1;\; i \in J)\leq \mathbb{P}(Y_{i}=1, i \in J)$ for any $J \subset \mathbb{Z}^d$ .
Proposition 3.3. Assume there are constants $C_1, C_2, \delta>0$ such that $\lambda^*(\textbf{x},\gamma)\leq C_1$ everywhere and $\lambda^*(\textbf{x},\gamma)\geq C_2$ whenever $d(\textbf{x},\gamma)\geq \delta$ . If $\Gamma$ is the corresponding Gibbs point process, then for any $\varepsilon>0$ there exists $p>0$ such that
where B is a random variable with values in $\{0,1\}^{\otimes \mathbb{Z}^d}$ such that $B_i$ , $i \in \mathbb{Z}^d$ , are independent with Bernoulli distribution B(p).
See Section 5.2 for the proof.
Remark 3.3. (On the constant $\varepsilon$ .) The constant $\varepsilon$ in Proposition 3.3 is merely an auxiliary for the proof. We need the side length of $\mathcal{D}_k$ to be slightly bigger than $2\delta$ so that we can fit another box $\mathcal{C}_k$ of side length $\varepsilon$ inside $\mathcal{D}_k$ in such a way that it has distance exactly $\delta$ from the complement of $\mathcal{D}_k$ . Depending on $\varepsilon$ , though, one can optimize the value of p for further applications.
Corollary 3.2. Let the range of interaction $R_{\textbf{x}}$ be decreasing with respect to the order in $\textbf{N}$ , and suppose there exists a decreasing function $R'\;:\;\{0,1\}^{Z^d}\to \mathbb{R}$ such that $R_{\textbf{x}}(\gamma)\leq R'(I^{\delta}(\gamma))$ for all $\gamma \in \textbf{N}$ and some $\delta>0$ . Then, under the assumptions and notation of Proposition 3.3, we have
Note that $R_{\textbf{x}}$ and R’ are defined on different spaces, and so the assumption of being decreasing has a slightly different meaning.
We know exactly the distribution of B. Therefore, it is usually a simple task to compute $\mathbb{P}(R'(B)> r)$ . Then, using Corollary 3.2, we can estimate the moments of $R_{\textbf{0}}(\Gamma)$ similarly as in (3.2) for any increasing function $f\;:\;\mathbb{R}\to \mathbb{R}$ by
Some applications of both Corollary 3.1 and Corollary 3.2 will be presented in Sections 4.3 and 4.4.
4. Examples
4.1. Pair potentials
A classical example of a Gibbs point process in $\mathbb{R}^d$ is the model with pairwise interactions. In this section, we omit the marks.
Definition 4.1. (Pair potential.) A Gibbs point process has a pair potential if there is a measurable, symmetric function $\Phi\;:\; \mathbb{R}^d \to \mathbb{R}\cup\{+ \infty\}$ such that the Papangelou intensity has the form
In the following definition, we recall a classical stability assumption from the pairwise model originating in [Reference Ruelle25]. Note that these conditions jointly guarantee the existence of a Gibbs point process with pair potential $\Phi$ .
Definition 4.2. (Superstable pairwise interactions.) Let $\Gamma$ be a Gibbs point process with pair potential $\Phi$ . For any finite configuration $\gamma \in \textbf{N}_f$ , we define the energy of this configuration by
We say that $\Phi$ is
-
superstable if for any bounded $\Lambda \subset \mathbb{R}^d$ there exist constants $A>0, B\geq 0$ such that
\begin{align*}H(\gamma_\Lambda)\geq A N_{\Lambda}(\gamma)^2 - B N_\Lambda(\gamma),\quad \forall \gamma \in \textbf{N},\end{align*} -
lower-regular if there is a positive, decreasing function $\varphi\;:\; [0,\infty)\to \mathbb{R}$ such that
\begin{align*}\int_0^{\infty} x^{d-1}\varphi(x) \mathrm{d}\, x <\infty\end{align*}and $\Phi(x)\geq - \varphi(|x|)$ for all $x \in \mathbb{R}^d$ , and -
integrable if
(4.1) \begin{align} \int_{\mathbb{R}^d}\left|1-e^{-\beta\Phi(y)}\right|\mathrm{d}\, y <\infty.\end{align}
Remark 4.1. (Assumptions (a1)–(a2) for superstable interactions.) A Gibbs point process with an integrable pair potential $\Phi$ automatically satisfies Assumption (a2) of Theorem 2.1 with $\phi=\Phi$ , since
Therefore, Assumption (a2) of Theorem 2.1 reduces to the standard integrability condition for the pair potential $\Phi$ . On the other hand, the existence of the second moment stems from the assumption of superstability and lower-regularity (see [Reference Ruelle25, Corollary 5.3]).
Corollary 4.1. A Gibbs point process with superstable, lower-regular, and integrable pair potential is not hyperuniform.
Remark 4.2 (Bound on the asymptotic variance.) Following the proof of Theorem 2.1, one can derive the following lower bound for the asymptotic variance of Theorem 2.1 for the Gibbs point process with integrable pair potential $\Phi$ :
Example 4.1. The class of Gibbs point processes with superstable, integrable pairwise interactions is large and covers many standard examples, including the following:
-
1. The Strauss process, i.e. a Gibbs point process $\Gamma$ with pair potential
\begin{align*}\Phi(x-y)=\begin{cases}1 \quad \mathrm{ if } |x-y|\leq R,\\0 \quad \mathrm{ otherwise,}\end{cases}\end{align*}for some $R\in [0,\infty)$ . Let $\lambda\;:\!=\;\mathbb{E}\lambda^*(0,\Gamma)$ be the intensity of the process; then, directly from (4.2), we get the lower bound\begin{align*}\mathcal{C}_{nhyp}\geq \frac{\lambda^2}{z+z^2|B(0,R)|(1-e^{-\beta})}.\end{align*} -
2. Riesz gases with $s>d$ , i.e. processes with a pair potential of the form $\Phi(x)=\|x\|^{-s}$ , $x \in \mathbb{R}^d$ , $s \in \mathbb{R}$ . The case $s\leq d$ (including that of a Coulomb gas, $s=d-2$ ) determines a non-integrable pair potential.
-
3. The Lennard–Jones pair potential, given by
\begin{align*}\Phi(x)= A\|x\|^{-\alpha_1}- B\|x\|^{-\alpha_2}\end{align*}for some $A, B >0$ and $\alpha_1> \alpha_2> d$ .
4.2. Widom–Rowlinson models
We start with a simple model in which the size of the balls is the same deterministic constant for all the points. Thus, for $R>0$ , we define
The Widom–Rowlinson point process is a Gibbs point process $\Gamma$ with an energy function
The Papangelou intensity has the form
for some $z>0$ , $\beta\geq 0$ . Clearly, $z e^{-\beta |B(0,R)|}\leq \lambda^*(x,\gamma)\leq z$ ; hence all moments of $\lambda^*$ are finite. Moreover, $\Gamma$ has a finite range of interaction that produces non-hyperuniformity by (i) and (iii) of Proposition 3.1. Alternatively, we can directly apply Theorem 2.1 with $\delta=0$ and $\phi(x)=I_{B(0,2R)}(x) e^{\beta |B(0,R)|}$ . Similarly, a Gibbs point process with quermass interaction among balls with a deterministic size is non-hyperuniform, i.e. a Gibbs point process with energy
where Per is the perimeter and $\chi$ the Euler–Poincaré characteristic (the number of connected components minus the number of holes).
Moreover, using the above-mentioned estimates for the Papangelou intensity and following the steps in the proof of Theorem 2.1, we obtain the following lower bound for the asymptotic variance:
Next, assume each point $x\in \Gamma$ is equipped with a non-negative random variable $R_x$ independently and with the same distribution $\mathbb{Q}$ . We are now in the setting of a marked Gibbs point process with $\mathbb{M}= \mathbb{R}_+$ and $\lambda_{\mathbb{M}}=\mathbb{Q}$ . As in the previous example, we define
and
Note that a Gibbs point process associated with Papangelou intensity (4.3) is well defined. The existence of the process with interactions also involving other Minkowski functionals is proved in [Reference Dereudre6].
Corollary 4.2. Assume that $\mathbb{E}_{\mathbb{Q}} e^{\alpha \beta |B(0,1)| R^d}<\infty$ for some $\alpha>1$ ; then a Gibbs point process defined by the Papangelou intensity (4.3) is non-hyperuniform.
Proof. As in Example 4.2, $\lambda^*\leq z$ . Hence, Assumption (A1) of Theorem 2.2 is trivially satisfied by Proposition 3.1(ii). It remains to show that Assumption (A2) is true for some $\alpha_2>1$ . Furthermore, for any $\textbf{x}\;:\!=\;(x,R_x), \textbf{y}\;:\!=\;(y,R_y) \in E\;:\!=\;\mathbb{R}^d\times\mathbb{R}_+$ , we have
which does not depend on the configuration $\gamma$ . Therefore, for any $\alpha_2>1$ ,
Now, take $\zeta_1, \zeta_2>1$ such that $\zeta_1<\alpha$ and $1/\zeta_1+1/\zeta_2=1$ . Furthermore, set $\alpha_2 = \alpha/\zeta_1$ and take $\zeta>\alpha_2 \zeta_2 d$ . We successively use (4.4), Hölder’s inequality with $\zeta_1,\zeta_2$ , and Markov’s inequality with $\zeta$ to get
Remark 4.3. (Gibbs particle process.) Generally, we may replace balls of random radius by random compact sets to define the Gibbs particle process. Assuming that the diameters of these sets have suitable exponential moments, the resulting process is again non-hyperuniform. Consequently, the crucial assumption of positive asymptotic variance in [Reference Beneš, Hofer-Temmel, Last and Večeřa1] is satisfied, enabling us to study the limit behaviour of U-statistics of Gibbs particle processes.
4.3. Voronoi interactions
In the previous example, we had direct control over the sizes of the discs through marks and hence over the range of interaction. Here, we may encounter cells that are very large, that interfere with distant neighbourhoods, and whose distribution is generally unknown. To define a Gibbs point process with interactions among Voronoi cells, for non-empty $\gamma\in \textbf{N}$ and $x \in \gamma$ , we denote the cell around x by
The set $C(x, \gamma)$ represents those points in $\mathbb{R}^d$ such that x is the nearest point to them out of all points in $\gamma$ . If $C(x,\gamma \cup\{x,y\})\cap C(y,\gamma\cup\{x,y\})\neq \emptyset$ , we write $x\stackrel{\gamma}{\sim} y$ . Finally, if $\gamma\in \textbf{N}$ and $x\notin \gamma$ , we write $C(x,\gamma)$ for $C(x,\gamma\cup \{x\}).$ The set $C(x, \gamma)$ is in fact a closed convex set, since it can be written as an intersection of closed half-spaces. The set of all closed convex sets is denoted by $\mathcal{C}^d$ .
We consider a function $\Phi\;:\;\mathcal{C}^d\to \mathbb{R}\cup \{\infty\}$ and write $\Phi(x,\gamma)$ for $\Phi(C(x,\gamma))$ . For a finite configuration $\gamma \in \textbf{N}$ , we consider energy
The Gibbs point process associated to H is a point process satisfying (GNZ) with Papangelou intensity of the form
In order to analyse Assumptions (A1) and (A2) of Theorem 2.2, we provide a list of assumptions imposed on the function $\Phi$ .
Definition 4.3. We say that the function $\Phi$ is
-
sub-additive if $\Phi(C)\leq \Phi(C_1)+\Phi(C_2)$ whenever $C =C_1\cup C_2$ , where $C,C_1,C_2 \in \mathcal{C}^d$ ,
-
increasing if $\Phi(C)\leq \Phi(C')$ whenever $C,C' \in \mathcal{C}^d$ and $C \subseteq C'$ , and
-
controlled by the volume if there exists a constant $K>0$ such that $|\Phi(C)|\leq \min\{|C|,K\}$ for all $C \in \mathcal{C}^d$ .
Note that, by [Reference Dereudre, Drouilhet and Georgii8], a Gibbs point process associated with $\lambda^*$ in (4.6) with $\Phi$ as above exists. The conditions on $\Phi$ stated above allow us to use stochastic minoration by the Poisson point process in order to gain control over the size of the typical cell.
Proposition 4.1. If $\Phi$ is sub-additive, increasing, and controlled by the volume, then
Proof. For any $y \in \gamma$ , we have $C(y,\gamma \cup \{x\})\subseteq C(y,\gamma)$ , while both sides are elements of $\mathcal{C}^d$ and hence $\Phi(y,\gamma \cup \{x\}) \leq \Phi(y,\gamma)$ since $\Phi$ is increasing. Consequently,
Moreover, $\Phi$ is uniformly bounded from above by K, yielding $-\beta \Phi(x,\gamma)\geq -\beta K$ . This proves the lower bound.
To obtain the upper bound, for $\gamma \in \textbf{N}$ and $y\stackrel{\gamma}{\sim} x$ we define the set $K_y\;:\!=\; C(y,\gamma)\setminus C(y,\gamma\cup\{x\})\in \mathcal{C}^d.$ Then
but also
since $|K_y\cap K_{y'}|=0$ for any $y\neq y'$ .
Note that it is enough to consider only the neighbouring points of x in (4.6), that is, to sum over $\{y \in \gamma;\; x\stackrel{\gamma}{\sim} y\}$ . Otherwise, adding a point x to the configuration $\gamma$ (or removing it) does not affect the shape of $C(y,\gamma)$ (or, consequently, the value of $\Phi(y,\gamma)$ ). Using the sub-additivity of $\Phi$ , control by the volume, and (4.8), we arrive at
The latter combined with the fact that $-\Phi\leq K$ produces the desired upper bound.
Corollary 4.3. Let $\Phi$ be some function on $\mathcal{C}^d$ that is sub-additive, increasing, and controlled by the volume. Then a Gibbs point process $\Gamma$ associated with the Papangelou intensity (4.6) is non-hyperuniform for any $\beta < \beta_c$ , where $\beta_c$ is the unique solution of the equation
with
The proof is postponed to Section 5.
Remark 4.4. (Values of $C_d$ and $\beta_c$ .) The constant $C_d$ can easily be evaluated. It can be seen that it decreases with the dimension d; for instance $C_1 = \frac{1}{6}$ , $C_2= \frac{1}{36}$ , $C_3\approx 0.006$ , etc. As a consequence, the exact value of $\beta_c$ can be also given in terms of z and K. If $d=2$ and $z=K=1$ , then $\beta_c \approx 0.03$ .
Remark 4.5. (General tessellations.) We can also formulate Corollary 4.3 for interactions based on more general tessellations, such as the Laguerre and Johnson–Mehl types. In this context, the shape of the cells is determined additionally by marks, which justifies our consideration of marks in Theorems 2.1 and 2.2.
4.4. Interactions based on k-nearest-neighbours graph
For $k\in \mathbb{N}$ and $\gamma \in \textbf{N}$ , denote by $v^i(x,\gamma)$ the ith nearest neighbour of x in $\gamma$ , $i=1,\ldots,k$ , and let $V^k(x,\gamma)\;:\!=\;\{v^i(x,\gamma);\; i=1,\ldots,\min\{k,N(\gamma)-1\}\}$ be the set of the first k neighbours of x in $\gamma$ . Here, $N(\gamma)$ is the cardinality of $\gamma$ . If there are two or more points at the same distance from a given point, we use the lexicographic ordering as a tie-breaker to determine the k-nearest-neighbours structure. However, such ties have zero probability for the random point sets considered here.
For a finite configuration $\gamma$ , we consider an energy of the form
where $\Phi\;:\;\mathbb{R}^d\to \mathbb{R}\cup \{\infty\}$ is some measurable function. Adding a new point x to the configuration will change the nearest-neighbour structure, but only locally. If $x \in V^k(y,\gamma\cup\{x\})$ for some $y \in \gamma$ , then $V^k(y,\gamma\cup\{x\})= V^{k-1}(y,\gamma)\cup\{x\}.$ Otherwise, it is obvious that $V^k(y,\gamma\cup\{x\})=V^k(y,\gamma)$ . Because of this fact, the corresponding Papangelou intensity takes the form (2.1) with
Proposition 4.2. There exists a constant $N_d$ , depending only on the dimension d, such that
Proof. In (4.9), the indicator $\textbf{1}_{V^k(y,\gamma\cup\{x\})}(x)$ takes the value 1 only for finitely many points $y \in \gamma$ . The number of such points is random, yet bounded by $k N_d$ , where $N_d$ depends only on the dimension (see [Reference Penrose23, Lemma 4.3] for the proof).
Corollary 4.4. Let $\Phi\;:\;\mathbb{R}^d\to \mathbb{R}$ be such that $\|\Phi\|_{\infty}<\infty$ . Then the Gibbs point process defined by the local energy (4.9) is non-hyperuniform.
The Papangelou intensity in Corollary 4.4 is double locally stable. Therefore, it is enough to verify the moment condition of Proposition 3.1(iv). Alternatively, if we assume that $\Phi$ is decreasing and non-negative, we are allowed to have an explosion around the origin.
Corollary 4.5. If $\Phi(x)\;:\!=\; \Phi(\|x\|)$ is a decreasing and non-negative function on $\mathbb{R}_+$ , then the Gibbs point process defined by (4.9) is non-hyperuniform.
Example 4.2. (Coulomb interaction.) For $d\geq 3$ , let $\Phi(x)= \frac{1}{\|x\|^{d-2}}$ . Then by Corollary 4.5, the Gibbs point process defined by (4.9) is non-hyperuniform.
For the proofs, see Section 5.
5. Proofs of the main results
This section aims to present the proof of the main result of this paper, together with the proofs of our examples and auxiliary propositions that involve more technical details.
5.1. Proofs of Theorem 2.1 and Theorem 2.2
We start by proving Theorem 2.2. A slight modification of the proof then yields the statement of Theorem 2.1.
Proof. We aim to show that there is a constant $F>-1$ not depending on $\Lambda$ such that
By the stationarity of $\Gamma$ , this already implies (2.3).
In the spirit of the proof of [Reference Ruelle25, Proposition 4.1], we define $R_\Lambda\;:\!=\; R_\Lambda (\gamma)\;:\!=\; \int_{\Lambda\times\mathbb{M}} \lambda^*(\textbf{x},\gamma)\nu(\mathrm{d}\, \textbf{x})$ for $\gamma \in \mathcal{N}$ . One has directly by (GNZ) that $\mathbb{E} N_\Lambda = \mathbb{E} R_\Lambda$ . Moreover, note that
where for the latter equality we apply (GNZ) twice, first with $f(\textbf{x},\gamma)\;:\!=\; \textbf{1}_{\Lambda\times\mathbb{M}} (\textbf{x}) \sum_{\textbf{y}\in \gamma} \textbf{1}_{\Lambda\times\mathbb{M}}(\textbf{y})$ and then with $g(\textbf{y},\gamma)\;:\!=\; \textbf{1}_{\Lambda\times\mathbb{M}} (\textbf{y}) \int_{\Lambda\times\mathbb{M}} \lambda^*(\textbf{x},\gamma\cup \{\textbf{y}\})\nu(\mathrm{d}\, \textbf{x}).$ Similarly, by the definition of $R_\Lambda$ and (GNZ), we have
Suppose for a moment that
for some $F\in (-1,\infty)$ . We then have that
Here, we have used consecutively (GNZ), the Cauchy–Schwartz inequality, (5.4), and the assumption (5.5). From this, the assertion (2.3) immediately follows.
It remains to check the validity of (5.5). First, we write $\mathcal{M}=\mathcal{M}_1 + \mathcal{M}_2+ \mathcal{M}_3$ , where
Note that
To proceed further, we define the quantities
which are both assumed to be finite. Using Fubini’s theorem, Hölder’s inequality with respect to the product measure $P_{\Gamma}\otimes \lambda_{\mathbb{M}}^2$ ( $P_{\Gamma}$ is the distribution of $\Gamma$ ), and finally the stationarity of $\Gamma$ , we arrive at
By Assumptions (A1) and (A2), $F_1 <\infty$ . By the same arguments,
is finite by (A1). Finally, $\mathcal{M}_3\leq 0$ and hence (5.5) holds true with $F=F_1 + F_2$ .
We continue with a proof of Theorem 2.1.
Proof. Clearly, for the constants $\mathcal{M}_2, \mathcal{M}_3$ from the proof of Theorem 2.2, one can use the same arguments in order to find finite upper bounds. For $\mathcal{M}_1$ , we use the Cauchy–Schwartz inequality to see that
The rest of the proof follows the same path as the proof of Theorem 2.2.
5.2. Proof of Proposition 3.3
Proof. For all $k\in \mathbb{Z}^d$ we denote by $\mathcal{C}_k$ a cube of a side length $\varepsilon$ lying in the centre of $\mathcal{D}_k$ , with the two cubes having parallel edges. Now, let $k \in \mathbb{Z}^d$ and $\gamma \in \textbf{N}$ be chosen to be arbitrary, but fixed. Let $\Pi^z_{\mathcal{D}_k}$ be the marked Poisson point process in $\mathcal{D}_k\times \mathbb{M}$ with intensity z and $\textbf{N}_{\mathcal{D}_k}\subset \textbf{N}$ the set of configurations restricted to $\mathcal{D}_k\times\mathbb{M}$ . It can be shown that $\Gamma$ also satisfies the Dobrushin–Lanford–Ruelle (DLR) equations (see [Reference Dereudre7, Theorem 1]). Then, the distribution of points in $\mathcal{D}_k$ given the configuration in $\mathcal{D}_k^C$ is precisely given by
where
Since $\lambda^*$ is locally stable from above,
The constant $Z^{\beta,z}$ depends on $\delta$ , $\varepsilon$ , $C_1$ , and d, but not on k or on the configuration $\gamma_{\mathcal{D}_k^C}$ . Using the fact that $d(\textbf{x}, \gamma_{\mathcal{D}_k^C})\geq \delta$ for any $\textbf{x}\in \mathcal{C}_k\times\mathbb{M}$ , we have that
where p depends neither on k nor on the boundary condition $\gamma_{\mathcal{D}_k^C}$ . Note that the quantity $\lambda^*(\gamma_{\mathcal{C}_k},\gamma_{\mathcal{D}_k^C})$ is well defined, since $\gamma_{\mathcal{C}_k}$ is assumed to be a.s. a one-point set.
To prove the statement, we construct a disagreement coupling of $I^{2\delta+\varepsilon}(\Gamma)$ and $B \sim B(p)^{\otimes \mathbb{Z}^d}$ . For this purpose, let I be a finite index set and $X\;:\!=\;(X_i)_{i \in I}$ random variables, not necessarily independent or identically distributed, with values in $\{0,1\}$ . Assume that
Then also
and similarly,
for any $J\subseteq I\setminus\{i\}$ .
Let $U\;:\!=\;(U_i)_{i \in I}$ be a vector of independent and identically distributed (i.i.d.) uniform random variables on [0,1] and define $B_{i}\;:\!=\; \textbf{1}[U_{i}<p], i \in I$ . Clearly, $B_I\;:\!=\;(B_i)_{i \in I}$ is a vector of i.i.d. Bernoulli variables with parameter p.
For the coupling, we define $Z_{i_1}= \textbf{1}[U_{i_1} < \mathbb{P}(X_{i_1}=1)]$ , $i_1 \in I$ . It is easy to check that $\mathbb{P}(Z_{i_1}= z)=\mathbb{P}(X_{i_1}= z)$ for $z \in \{0,1\}$ and $\mathbb{P}(B_{i_1}\leq Z_{i_1})=1$ . Inductively, for $k\geq 1$ , let $((Z_{i_1},\ldots Z_{i_k}), (B_{i_1},\ldots B_{i_k}))$ be a coupling of $(X_{i_1},\ldots X_{i_k})$ and $(B_{i_1},\ldots B_{i_k})$ , such that $\mathbb{P}(Z_{i_j}\geq B_{i_j})=1$ for any $j=1,\ldots,k$ . Then we define
Again, $\mathbb{P}(Z_{i_{k+1}}\geq B_{i_{k+1}})=1$ , and for $z_1,\ldots z_k \in \{0,1\}$ we compute
Let $Z\;:\!=\; (Z_i)_{i \in I}$ . We have shown that $(Z,B_I)$ is a coupling of X and $B_I$ yielding
This remains true if we turn to the limit and take $I_n \nearrow \mathbb{Z}^d$ . The choice $X = (I^{2\delta+\varepsilon}(\Gamma))_{i \in \mathbb{Z}^d}$ satisfies the assumption (5.7), and hence, by the arguments above,
5.3. Proof of Corollary 4.3
Proof. As before, one needs to verify Assumptions (A1)–(A2) of Theorem 2.2. Here, it is enough to do so in the form without marks.
By Proposition 4.1, the Papangelou intensity of $\Gamma$ is locally stable from below, and hence by Proposition 3.2 there exists a Poisson point process $\Pi$ such that $\Pi \ll \Gamma$ and $\Pi$ has the intensity $z e^{-K \beta}$ , where K is as in (4.7) and $z, \beta$ are the parameters defining the Papangelou intensity (2.1) of the Gibbs point process $\Gamma$ . This stochastic minoration can be interpreted via coupling of $\Pi$ and $\Gamma$ such that $\gamma'\subseteq \gamma$ whenever $\gamma'\sim\Pi$ and $\gamma \sim \Gamma$ . For Voronoi tessellations, it translates to $C(x,\gamma)\subseteq C(x,\gamma')$ for any $x \in \gamma\cap\gamma'$ . Combining this with the upper bound from Proposition 4.1, we may write
Now we aim to construct a ball B(0,R) such that $R=R(\Pi)$ is a random variable depending on $\Pi$ and $C(0,\Pi)\subseteq B(0,R)$ a.s. We use exactly the approach in the proof of [Reference Penrose24, Lemma 5.1], the only difference being that we are interested in more precise estimates of the tail probabilities $\mathbb{P}(R>t)$ , $t>0$ .
Let $K_1,\ldots, K_J$ be set of circular cones with apices in the origin with angular radii $\pi/6$ . We do not expect the cones to have zero-volume intersections, yet we choose J to be the minimum value such that $\cup_{j=1}^J K_j =\mathbb{R}^d$ . Note that J depends on the dimension and is finite. For each $j=1,\ldots,J$ choose $x_j \in \gamma' \cup K_j$ to be the closest point to the origin.
Denote by $H_x(y)\;:\!=\; \{z \in \mathbb{R}^d;\; \|z-y\|\leq \|z-x\|\}$ the closed half-space induced by points that are closer to y than to x. By the definition of a Voronoi cell,
Set $R\;:\!=\;\max\{\|x_j\|; j=1,\ldots,J\}$ . We need to verify that $C(0,\Pi) \subseteq B(0,R)$ a.s. To do so, we choose $y \in H_{x_j}(0) \cap K_j$ and show that $\|y\|\leq \|x_j\|$ for any $j\in \{1,\ldots,J\}$ . It follows from a simple computation that
where the first inequality holds since $y \in H_{x_j}(0)$ and the last one from the fact that $y \in K_j$ (note that $z_1,z_2 \in K_j$ implies $\langle z_1, z_2\rangle \geq 1/2\|z_1\|\|z_2\|$ ). Finally, we have that $C(0,\Gamma)\subseteq B(0,R(\Pi))$ a.s.
Using void probabilities of the Poisson point process $\Pi$ , we arrive at the estimate
where $c_d \;:\!=\; |K_1\cap B(0,1)|_d$ . We leave it to the reader to check that
Ultimately, for (A1), we have that
The latter integral converges as long as
Similarly, for (A2), we have from the assumptions on $\Phi$ that
Therefore, by an additional application of Fubini’s theorem and the fact that $1/\alpha_2\leq 1$ ,
A computation in the same spirit as in (5.8) shows that the expectation above is finite as long as
The two inequalities (5.9) and (5.10) are optimal for $\alpha_1=\frac{3}{2}$ and $\alpha_2=3$ .
Remark 5.1. We highlight two facts regarding the construction of the estimates of the tail probabilities of R. First, the construction of the estimates is not hindered by the fact that the cones $K_1,\ldots, K_J$ may intersect and we may possibly choose some point $x \in \gamma'$ multiple times. Second, it can be seen from the exponential form of the tail probabilities that the circular cones are the optimal choice for this construction. They give us more precise estimates than any other solids would do (e.g. non-circular and non-intersecting).
5.4. Proofs of Corollary 4.4 and Corollary 4.5
We begin with a proof of Corollary 4.4.
Proof. The Papangelou intensity is double locally stable by Proposition 4.2; ergo Assumption (A1) of Theorem 2.2 is justified by Proposition 3.1(i).
It remains to validate (A2). For that, we construct $R\;:\!=\; R(\gamma)$ such that R is the range of interaction, i.e.
The construction is almost the same as in the proof of [Reference Penrose and Yukich22, Lemma 6.1], where the authors prove stabilization for a planar undirected k-nearest-neighbours graph (there is an edge between x, y whenever $x \in V^k(y,\gamma)$ or $y \in V^k(x,\gamma)$ ). Here, we consider a directed graph, i.e. such that there is an edge pointing from x to y whenever $y\in V^k(x,\gamma)$ . The idea, however, remains the same.
Let J be the smallest integer such that $K_1,\ldots,K_J$ are cones with apex at the origin and angular radius at most $\pi/6$ , such that $\cup_{j=1}^J K_j= \mathbb{R}^d$ . Note that, unlike in the proof of Corollary 4.3, we do not need to optimize the shape of the cones. In fact, by Proposition 3.1(iv), all we need to show is that the radius $R$ has finite $\alpha$ th moment for some $\alpha>d$ and set $\alpha_2=\alpha/d$ . Let $D\;:\!=\;D(\gamma)$ be the smallest $t>0$ such that there are at least $k+1$ points in each $K_j\cap B(0,t)$ , $j=1,\ldots,J$ , and set $R=2D$ .
First, for any $y\in \gamma$ such that $0\in V^k(y,\gamma)$ there exists $j\in\{1,\ldots,J\}$ such that $y \in K_j\cap B(0,D)$ . Otherwise, y would have at least $k+1$ points that are closer to y than the origin, by the construction of D, and that contradicts $0\in V^k(y,\gamma)$ . In addition, $V^k(y,\gamma)\in B(0,R)$ . No point from $B(0,R)^C$ can be among the k nearest neighbours of y, because $y\in K_j\cap B(0,D)$ implies that there are at least k points closer to y than a potential neighbour outside B(0,R). We conclude that R is a decreasing range of interaction for the Papangelou intensity $\lambda^*$ as in Definition 3.2.
By Proposition 4.2 and Proposition 3.2, there is a Poisson point process $\Pi$ with intensity $\lambda\;:\!=\;z e^{-\beta (1+2 N_d) k N_d\|\Phi\|_{\infty}}$ such that $\Pi \ll \Gamma$ . By Corollary 3.1, for any $t>0$ ,
Let $\alpha>d$ . Then
The latter term is finite, as it is a finite sum of converging integrals. Thus, Assumption (A2) of Theorem 2.2 is satisfied for $\alpha_2=\alpha/d$ by Proposition 3.1(iv).
Next we prove Corollary 4.5.
Proof. Since $\Phi\geq 0$ everywhere, the Papangelou intensity is locally stable from above; hence Assumption (A1) of Theorem 2.2 is trivially satisfied by Proposition 3.1(i). For (A2), an easy computation leads to
where we used the fact that $\Phi(x)$ is non-negative and decreasing. Let $\delta>0$ and $\zeta <\infty$ be such that $\Phi(x)<\zeta$ whenever $\|x\|\geq \delta/2$ . Then, for $\alpha_2>1$ ,
where $A\;:\!=\; \exp\{\beta N_d \zeta \}<\infty$ . In the third line of the latter expression, we used that $\|x\|\geq \|y\|/2 \geq \delta/2$ . Otherwise, $\|x-0\|\leq \|x-y\|$ implies that if $0 = v^{k+1}(x,\Gamma \cup\{0,y\})$ , then y cannot be among the k nearest neighbours of x.
Again, we let J be the smallest integer such that $K_1,\ldots,K_J$ are cones with apex at the origin and angular radius at most $\pi/6$ such that $\cup_{j=1}^J K_j= \mathbb{R}^d$ . For $\gamma \in \textbf{N}$ , we define $D\;:\!=\;D(\gamma)$ as the smallest $t>0$ such that there are at least $k+1$ points of $\gamma$ in each $K_j\cap B(0,t)$ , $j=1,\ldots,J$ , and set $R=2D$ . Then
If $\|y\|> R$ , then any point $x \in \gamma$ with $y \in V^k(x,\gamma\cup \{0,y\})$ and $0=v^{k+1}(x, \gamma\cup \{0,y\})$ satisfies $x \in B(0, D)^C$ . But then there is $j \in \{1,\ldots, J\}$ with $x \in K_j$ such that at least $k+1$ points in the cone $K_j$ are closer to x than the origin, contradicting the fact that $0 = v^{k+1}(x,\gamma\cup\{0,y\})$ .
We are now in a position to apply Corollary 3.2. It can be seen from (4.9) that $\lambda^*(x,\gamma)\leq z$ . This is due to the fact that $\Phi$ is non-negative and decreasing. Furthermore, if $d(x,\gamma)\geq \delta/2$ , then $\lambda^*(x,\gamma)\geq z \exp\{- (1+N_d)k \zeta\}$ . Here, we also used the fact that the number of summands is bounded by a deterministic constant (see Proposition 4.2). By Proposition 3.3, for every $\varepsilon>0$ there is some $p>0$ such that $I^{2\delta+\varepsilon}(\Gamma)$ is minorated by a Bernoulli field B with parameter p.
For $l \in \mathbb{L}^d\;:\!=\; \{0,1\}^{\mathbb{Z}^d}$ , we define R’(l) by taking the smallest $t>0$ such that, for all $j=1,\ldots,J$ , the cone $K_j$ fully contains at least $k+1$ cubes $\mathcal{D}_{i_1},\ldots, \mathcal{D}_{i_{k+1}}\in \mathcal{D}$ such that $l_{i_1},\ldots, l_{i_{k+1}}=1$ . Then R’ is decreasing. By the construction of $I^{2\delta+\varepsilon}(\Gamma)$ and Proposition 3.3 we have that, a.s.,
and hence, by Corollary 3.2,
For $j \in \{1,\ldots,J\}$ , take the cube $\mathcal{D}_{i_1}\subset K_j$ that is closest to the origin, and denote by $c_j\;:\!=\; \inf_{x \in \mathcal{D}_{i_1}} d(x,0)$ its distance to the origin. Then there exists another cube $\mathcal{D}_{i_2}\subset K_j$ sharing exactly one vertex with $\mathcal{D}_{i_1}$ at a distance $c_1+c_2$ from the origin, where $c_2$ is the body diagonal length of $\mathcal{D}_{i_1}$ . Note that $c^j_1$ and $c_2$ depend on d, $\delta$ , and $\varepsilon$ , and $c^j_1$ moreover depends on j. Inductively, we construct a chain of cubes $\mathcal{D}_{i_j}$ , $j \in \mathbb{N}$ , all fully included in the cone $K_j$ . Let $B=(B_i)_{i \in \mathbb{Z}^d}$ be distributed according to $B(p)^{\otimes \mathbb{Z}^d}$ . Define $R^j_{diag}(B)$ as the smallest $t= c^j_1+c_2 q$ such that $\sum_{j=1}^q B_{i_j}=k+1$ . To simplify the notation, for the rest of the proof let C be a universal finite constant depending on d, k, $\delta$ , $\varepsilon$ , $\alpha_2$ , $c_2$ , and $c^j_1$ , $j=1,\ldots,J$ . It can be seen that
Finally, since $p \in (0,1)$ ,
Acknowledgement
The authors would like to thank the anonymous referees and the editor handling our paper for their careful reading and comments, which improved the paper.
Funding information
This work was supported in part by the Labex CEMPI (ANR-11-LABX-0007-01), the ANR project RANDOM (ANR-19-CE24-0014), the CNRS GdR 3477 GeoSto, and the Czech Science Foundation, project no. 22-15763S.
Competing interests
There were no competing interests to declare which arose during the preparation or publication process of this article.