Dereudre David
email: david.dereudre@univ-lille1.fr

Introduction to the theory of Gibbs point processes

The Gibbs point processes (GPP) constitute a large class of point processes with interaction between the points. The interaction can be attractive, repulsive, depending on geometrical features whereas the null interaction is associated with the so-called Poisson point process. In a first part of this mini-course, we present several aspects of finite volume GPP defined on a bounded window in R d . In a second part, we introduce the more complicated formalism of infinite volume GPP defined on the full space R d . Existence, uniqueness and non-uniqueness of GPP are non-trivial questions which we treat here with completely self-contained proofs. The DLR equations, the GNZ equations and the variational principle are presented as well. Finally we investigate the estimation of parameters. The main standard estimators (MLE, MPLE, Takacs-Fiksel and variational estimators) are presented and we prove their consistency. For sake of simplicity, during all the mini-course, we consider only the case of finite range interaction and the setting of marked points is not presented.

Introduction

The spatial point processes are well studied objects in probability theory and statistics for modelling and analysing spatial data which appear in several disciplines as statistical mechanics, material science, astronomy, epidemiology, plant ecology, seismology, telecommunication, and others [START_REF] Baddeley | Case studies in spatial point process models[END_REF][START_REF] Daley | An introduction to the theory of point processes[END_REF]. There exist many models of such random points configurations in space and the most popular one is surely the Poisson point process. It corresponds to the natural way of producing independent locations of points in space without interaction. For dependent random structures, we can mention for instance the Cox processes, determinantal point processes, Gibbs point processes, etc. None of them is established as the most relevant model for ap-plications. In fact the choice of the model depends on the nature of the dataset, the knowledge of (physical or biological) mechanisms producing the pattern, the aim of the study (theoretical, applied or numerical).

In this mini-course, we focus on Gibbs point processes (GPP) which constitute a large class of points processes, able to fit several kinds of patterns and which provide a clear interpretation of the interaction between the points, such as attraction or repulsion depending on their relative position. Note that this class is particularly large since several point processes can be represented as GPP (see [START_REF] Georgii | Conditional intensity and Gibbsianness of determinantal point processes[END_REF][START_REF] Kozlov | Description of a point random field by means of the Gibbs potential[END_REF] for instance). The main disadvantage of GPP is the complexity of the model due to an intractable normalizing constant which appears in the local conditional densities. Therefore their analytical studies are in general based on implicit equilibrium equations which lead to complicated and delicate analysis. Moreover, the theoretical results which are needed to investigate the Gibbs point process theory are scattered across several publications or books. The aim of this mini-course is to provide a solid and selfcontained theoretical basis for understanding deeply the Gibbs point process theory. The results are in general not exhaustive but the main ideas and tools are presented in accordance with modern and recent developments. The main strong restriction here involves the range of the interaction, which is assumed to be finite. The infinite range interaction requires the introduction of tempered configuration spaces and for sake of simplicity we decided to avoid this level of complexity. The mini-course is addressed for Master and Phd students and also for researchers who want to discover or investigate the domain. The manuscript is based on a mini-course given during the conference of GDR 3477 géométrie stochastique, at university of Nantes in April 2016.

In a first section, we introduce the finite volume GPP on a bounded window Λ ⊂ R d . They are simply defined as point processes in Λ whose the distributions are absolutely continuous with respect to the Poisson point process distribution. The unnormalized densities are of form z N e -β H , where z and β are positive parameters (called respectively activity and inverse temperature), N is the number of points and H an energy function. Clearly, these distributions favour (or penalize) configurations with low (or high) energy E. This distortion strengthens as β is large. The parameter z allows to tune the mean number of points. This setting is relatively simple since all the objects are defined explicitly. However, the intractable normalization constant is ever a problem and most of quantities are not computable. Several standard notions (DLR and GNZ equations, Ruelle's estimates, etc.) are treated in this first section as a preparation for the more complicated setting of infinite volume GPP developed in the second section. Note that we do not present the setting of marked Gibbs point processes in order to keep the notations as simple as possible. However, all the results can be easily extended in this case.

In a second section, we present the theory of infinite volume GPP in R d . There are several motivations for studying such infinite volume regime. Firstly, the GPP are the standard models in statistical physics for modelling systems with a large number of interacting particles (around 10 23 according to the Avogadro's number). Therefore, the case where the number of particles is infinite is an idealization of this setting and furnishes microscopic descriptions of gas, liquid or solid. Macroscopic quantities like the density of particles, the pressure and the mean energy are consequently easily defined by mean values or laws of large numbers. Secondly, in the spatial statistic context, the asymptotic properties of estimators or tests are obtained when the observation window tends to the full space R d . This strategy requires the existence of infinite volume models. Finally, since the infinite volume GPP are stationary (shift invariant) in R d , several powerful tools, as the ergodic theorem or the central limit Theorem for mixing field, are available in this infinite volume regime.

The infinite volume Gibbs measures are defined by a collection of implicit DLR equations (Dobrushin, Lanford and Ruelle). The existence, uniqueness and nonuniqueness are non trivial questions which we treat in depth with self-contained proofs in this second section. The phase transition between uniqueness and non uniqueness is one of the most difficult conjectures in statistical physics. This phenomenon is expected to occur for all standard interactions although it is proved rigorously only for few models. The area interaction is one of such models and the complete proof of its phase transition is given here. The GNZ equations, the variational principle are discussed as well.

In the last section, we investigate the estimation of parameters which appear in the distribution of GPP. For sake of simplicity we deal only with the activity parameter z and the inverse temperature β . We present several standard procedures (MLE, MPLE, Takacs-Fiksel procedure) and a new variational procedure. We show the consistency of estimators, which highlights that many theoretical results are possible in spite of lack of explicit computations. We will see that the GNZ equations play a crucial role in this task. For sake of simplicity the asymptotic normality is not presented but some references are given.

Let us finish this introduction by giving standard references. Historically, the GPP have been introduced for statistical mechanics considerations and an unovoidable reference is the book by Ruelle [START_REF] Ruelle | Statistical mechanics: Rigorous results[END_REF]. Important theoretical contributions are also developed in two Lecture Notes [START_REF] Georgii | Canonical Gibbs measures[END_REF][START_REF] Preston | Random fields[END_REF] by Georgii and Preston. For the relations between GPP and stochastic geometry, we can mention the book [START_REF] Chiu | Stochastic geometry and its applications[END_REF] by Chiu et al. and for spatial statistic and numerical considerations, the book by Møller and Waagepetersen [START_REF] Møller | Statistical inference and simulation for spatial point processes[END_REF] is the standard reference. Let us mention also the book [START_REF] Van Lieshout | Markov point processes and their applications[END_REF] by van Lieshout on the applications of GPP.

Finite volume Gibbs point processes

In this first section we present the theory of Gibbs point process on a bounded set Λ ⊂ R d . A Gibbs point process (GPP) is a point process with interactions between the points defined via an energy functional on the space of configurations. Roughly speaking, the GPP produces random configurations for which the configurations with low energy have more chance to appear than the configurations with high energy (see Definition 2). In Section 1.1 we recall succinctly some definitions of point process theory and we introduce the reference Poisson point process. The energy functions are discussed in Section 1.2 and the definiton of finite volume GPP is given in Section 1.3. Some first properties are presented as well. The central DLR equations and GNZ equations are treated in Sections 1.4 and 1.5. Finally we finish the first section by giving Ruelle estimates in the setting of superstable and lower regular energy functions.

Poisson point process

In this first section, we describe briefly the setting of point process theory and we introduce the reference Poisson point process. We only give the main definitions and concepts and we suggest [START_REF] Daley | An introduction to the theory of point processes[END_REF][START_REF] Matthes | Infinitely divisible point processes[END_REF] for a general presentation.

The space of configurations C is defined as the set of locally finite subsets in R d :

C = {γ ⊂ R d , γ Λ := γ ∩ Λ is finite for any bounded set Λ ⊂ R d }.
Note that we consider only the simple point configurations, which means that the points do not overlap. We denote by C f the space of finite configurations in C and by C Λ the space of finite configurations inside Λ ⊂ R d .

The space C is equipped with the sigma-field F C generated by the counting functions N Λ for all bounded measurable Λ ⊂ R d , where N Λ : γ → #γ Λ . A point process Γ is then simply a measurable function from any probability space (Ω , F , P) to (C , F C ). As usual, the distribution (or the law) of a point process Γ is defined by the image of P to (C , F C ) by the application Γ . We say that Γ has finite intensity if, for any bounded set Λ , the expectation µ(Λ ) := E(N Λ (Γ )) is finite. In this case, µ is a sigma-finite measure called intensity measure of Γ . When µ = ζ λ d , where λ d is the Lebesgue measure on R d and ζ ≥ 0 a positive real, we simply say that Γ has finite intensity ζ .

The main class of point processes is the family of Poisson point processes, which furnish the natural way of producing independent points in space. Let µ be a sigmafinite measure in R d . A Poisson point process with intensity µ is a point process Γ such that, for any bounded Λ in R d , these properties both occur

• The random variable N Λ (Γ ) is distributed following a Poisson distribution with parameter µ(Λ ).

• Given the event {N Λ (Γ ) = n}, the n points in Γ Λ are independent and distributed following the distribution µ Λ /µ(Λ ).

The distribution of such a Poisson point process is denoted by π µ . When the intensity is µ = ζ λ d , we say that the Poisson point process is stationary (or homogeneous) with intensity ζ > 0, and denote its distribution π ζ . For any measurable set Λ ⊂ R d , we denote by π

Energy functions

In this section, we present the energy functions with the standard assumptions which we assume in this mini-course. The choices of energy functions come from two main motivations. First, the GPP are natural models in statistical physics for modelling continuum interacting particles systems. In general, in this setting the energy function is a sum of the energy contribution of all pairs of points (see expression [START_REF] Baddeley | Variational estimators for the parameters of Gibbs point process models[END_REF]). The GPP are also used in spatial statistics to fit as best as possible the real datasets. So, in a first step, the energy function is chosen by the user with respect to the characteristics of the dataset. Then the parameters are estimated in a second step.

Definition 1. An energy function is a measurable function

H : C f → R ∪ {+∞}
such that the following assumptions hold

• H is non-degenerate: H( / 0) < +∞.
• H is hereditary: for any γ ∈ C f and x ∈ γ then

H(γ) < +∞ ⇒ H(γ\{x}) < +∞.
• H is stable: there exists a constant A such that for any

γ ∈ C f H(γ) ≥ AN R d (γ).
The stability implies that the energy is superlinear. If the energy function H is positive then the choice A = 0 works but in the interesting cases, the constant A is negative. The hereditary means that the set of allowed configurations (configurations with finite energy) is stable when points are removed. The non-degeneracy is very natural. Without this assumption, the energy would be equal to infinity everywhere (by hereditary).

1) Pairwise interaction. Let us start with the most popular energy function which is based on a function (called pair potential)

ϕ : R + → R ∪ {+∞}.
The pairwise energy function is defined for any γ ∈ C f by

H(γ) = ∑ {x,y}⊂γ ϕ(|x -y|). (1)
Note that such an energy function is trivially hereditary and non-degenerate. The stability is more delicate and we refer to general results in [START_REF] Ruelle | Statistical mechanics: Rigorous results[END_REF]. However if ϕ is positive the result is obvious.

A standard example coming from statistical physics is the so-called Lennard-Jones pair potential where ϕ(r) = ar -12 + br -6 with a > 0 and b ∈ R. In the interesting case b < 0, the pair potential ϕ(r) is positive (repulsive) for small r and negative (attractive) for large r. The stability is not obvious and is proved in Proposition 3.2.8 in [START_REF] Ruelle | Statistical mechanics: Rigorous results[END_REF].

The Strauss interaction corresponds to the pair potential ϕ(r) = 1 [0,R] (r) where R > 0 is a support parameter. This interaction exhibits a constant repulsion between the particles at distance smaller than R. This simple model is very popular in spatial statistics.

The multi-Strauss interaction corresponds to the pair potential

ϕ(r) = k ∑ i=1 a i 1 ]R i-1 ,R i ] ,
where (a i ) 1≤i≤k is a sequence of real numbers and 0 = R 0 < R 1 < . . . < R k a sequence of increasing real numbers. Clearly, the pair potential exhibits a constant attraction or repulsion at different scales. The stability occurs provided that the parameter a 1 is large enough (see Section 3.2 in [START_REF] Ruelle | Statistical mechanics: Rigorous results[END_REF]).

2) Energy functions coming from geometrical objects. Several energy functions are based on local geometrical characteristics. The main motivation is to provide random configurations such that special geometrical features appear with higher probability under the Gibbs processes than the original Poisson point process. In this paragraph we give examples related to the Delaunay-Voronoi diagram. Obviously other geometrical graph structures could be considered.

Let us recall that for any

x ∈ γ ∈ C f the Voronoi cell C(x, γ) is defined by C(x, γ) = w ∈ R d , such that ∀y ∈ γ |x -w| ≤ |x -y| .
The Delaunay graph with vertices γ is defined by considering the edges

D(γ) = {x, y} ⊂ γ such that C(x, γ) ∩C(y, γ) = / 0 .
See [START_REF] Møller | Lectures on random Voronoȋ tessellations[END_REF] for a general presentation on the Delauany-Voronoi tessellations.

A first geometric energy function can be defined by

H(γ) = ∑ x∈γ 1 C(x,γ) is bounded ϕ(C(x, γ)), (2)
where ϕ is any function from the space of polytopes in R d to R. Examples of such functions ϕ are the Area, the (d -1)-Hausdorff measure of the boundary, the number of faces, etc... Clearly these energy functions are non-degenerate and hereditary.

The stability holds as soon as the function ϕ is bounded from below. Another kind of geometric energy function can be constructed via a pairwise interaction along the edges of the Delaunay graph. Let us consider a finite pair potential ϕ : R + → R. Then the energy function is defined by

H(γ) = ∑ {x,y}⊂D(γ) ϕ(|x -y|) (3)
which is again clearly non-degenerate and hereditary. The stability occurs in dimension d = 2 thanks to Euler's formula. Indeed the number of edges in the Delaunay graph is linear with respect to the number of vertices. Therefore the energy function is stable as soon as the pair potential ϕ is bounded from below. In higher dimension d > 2, the stability is more complicated and not really understood. Obviously, if ϕ is positive, the stability occurs.

Let us give a last example of geometric energy function which is not based on the Delaunay-Voronoi diagram but on a germ-grain structure. For any radius R > 0 we define the germ-grain structure of γ ∈ C by

L R (γ) = x∈γ B(x, R),
where B(x, R) is the closed ball centred at x with radius R. Several interesting energy functions are built from this germ-grain structure. First the Widom-Rowlinson interaction is simply defined by

H(γ) = Area(L R (γ)), (4)
where the "Area" is simply the Lebesgue measure λ d . This model is very popular since it is one of a few models for which the phase transition result is proved (see Section 2.8). This energy function is sometimes called Area-interaction [START_REF] Baddeley | Area-interaction point processes[END_REF][START_REF] Widom | New model for the study of liquid-vapor phase transitions[END_REF]. If the Area functional is replaced by any linear combination of the Minkowski functionals we obtain the Quermass interaction [START_REF] Dereudre | The existence of quermass-interaction processes for nonlocally stable interaction and nonbounded convex grains[END_REF]. Another example is the random cluster interaction defined by

H(γ) = Ncc(L R (γ)), (5)
where Ncc denotes the functional which counts the number of connected components. This energy function is introduced first in [START_REF] Chayes | The analysis of the Widom-Rowlinson model by stochastic geometric methods[END_REF] for its relations with the Widom-Rowlinson model. See also [START_REF] Dereudre | Infinite volume continuum random cluster model[END_REF] for a general study in the infinite volume regime.

Contents

Finite Volume GPP

Let Λ ⊂ R d such that 0 < λ d (Λ ) < +∞. In this section we define the finite volume GPP on Λ and we give its first properties.

Definition 2. The finite volume Gibbs measure on Λ with activity z > 0, inverse temperature β ≥ 0 and energy function H is the distribution

P z,β Λ = 1 Z z,β Λ z N Λ e -β H π Λ , (6) where Z
Λ ≥ π Λ ( / 0)e -β H({ / 0}) = e -λ d (Λ ) e -β H({ / 0}) > 0
and thanks to the stability of

H Z z,β Λ ≤ e -λ d (Λ ) +∞ ∑ n=0 (ze -β A λ d (Λ )) n n! = e λ d (Λ )(ze -β A -1) < +∞.
In the case β = 0, we recover that P z,β Λ is the Poisson point process π z Λ . So the activity parameter z is the mean number of points per unit volume when the interaction is null. When the interaction is active (β > 0), P z,β Λ favours the configurations with low energy and penalizes the configurations with high energy. This distortion strengthens as β is large.

There are many motivations for the exponential form of the density in (6). Historically, it is due to the fact that the finite volume GPP solves the variational principle of statistical physics. Indeed, P z,β Λ is the unique probability measure which realizes the minimum of the free excess energy, equal to the mean energy plus the entropy. It expresses the common idea that the equilibrium states in statistical physics minimize the energy and maximize the "disorder". This result is presented in the following proposition. Recall first that the relative entropy of a probability measure P on C Λ with respect to the Poisson point process π ζ Λ is defined by

I(P|π ζ Λ ) = log( f )dP if P π z Λ with f = dP dπ ζ Λ +∞ otherwise. ( 7
)
Proposition 1 (Variational Principle). Let H be an energy function, z > 0, β ≥ 0.

Then {P z,β Λ } = argmin P∈P Λ β E P (H) -log(z)E P (N Λ ) + I(P|π Λ )
, where P Λ is the space of probability measures on C Λ with finite intensity and E P (H) is the expectation of H under P, which is always defined (maybe equal to infinity) since H is stable.

Proof. First we note that

β E P z,β Λ (H) -log(z)E P z,β Λ (N Λ ) + I(P z,β Λ |π Λ ) = β HdP z,β Λ -log(z)E P z,β Λ (N Λ ) + log z N Λ e -β H Z z,β Λ dP z,β Λ = -log(Z z,β Λ ). (8)
Λ ) ≥ log { f >0} z N Λ e -β H dπ Λ = log z N Λ e -β H f -1 dP ≥ log z N Λ e -β H f -1 dP = -β E P (H) -log(z)E P (N Λ ) -log( f )dP.
The second inequality, due to the Jensen's inequality, is an equality if and only if z N Λ e -β H f -1 is P a.s. constant which is equivalent to

P = P z,β
Λ . The proposition is proved.

The parameters z and β allow to fit the mean number of points and the mean value of the energy under the GPP. Indeed when z increases, the mean number of points increases as well and similarly when β increases, the mean energy decreases. This phenomenon is expressed in the following proposition. The proof is a simple computation of derivatives.

Let us note that it is not easy to tune both parameters simultaneously since the mean number of points changes when β is modified (and vice versa). The estimation of the parameters z and β is discussed in the last Section 3. Let us finish this section by explaining succinctly how to simulate such finite volume GPP. There are essentially two algorithms. The first one is based on a MCMC procedure where GPP are viewed as equilibrium states of Markov chains. The simulation is obtained by letting run for a long enough time the Markov chain. The simulation is not exact and the error is essentially controlled via a monitoring approach (see [START_REF] Møller | Statistical inference and simulation for spatial point processes[END_REF]). The second one is a coupling from the past algorithm which provided exact simulations. However, the computation time is often very long and these algorithms are not really that used in practice (see [START_REF] Kendall | Perfect simulation using dominating processes on ordered spaces, with application to locally stable point processes[END_REF]).

DLR equations

The DLR equations are due to Dobrushin, Lanford and Ruelle and give the local conditional distributions of GPP in any bounded window ∆ given the configuration outside ∆ . We need to define a family of local energy functions (H ∆ ) ∆ ⊂R d . Definition 3. For any bounded set ∆ and any finite configuration γ ∈ C f we define

H ∆ (γ) := H(γ) -H(γ ∆ c ), with the the convention ∞ -∞ = 0.
The quantity H ∆ (γ) gives the energetic contribution of points in γ ∆ towards the computation of the energy of γ. As an example, let us compute these quantities in the setting of pairwise interaction introduced in (1);

H ∆ (γ) = ∑ {x,y}⊂γ ϕ(|x -y|) -∑ {x,y}⊂γ ∆ c ϕ(|x -y|) = ∑ {x, y} ⊂ γ {x, y} ∩ ∆ = / 0 ϕ(|x -y|).
Note that H ∆ (γ) does not depend only on points in ∆ . However, trivially we have H(γ) = H ∆ (γ) + H(γ ∆ c ), which shows that the energy of γ is the sum of the energy H ∆ (γ) plus something which does not depends on γ ∆ . Proposition 3 (DLR equations for finite volume GPP). Let ∆ ⊂ Λ be two bounded sets in R d with λ d (∆ ) > 0. Then for P z,β

Λ -a.s. all γ ∆ c P z,β Λ (dγ ∆ |γ ∆ c ) = 1 Z z,β ∆ (γ ∆ c ) z N ∆ (γ) e -β H ∆ (γ) π ∆ (dγ ∆ ), (9)
where Z z,β

∆ (γ ∆ c ) is the normalizing constant z N ∆ (γ) e -β H ∆ (γ) π ∆ (dγ ∆ ).
In particular the right term in (9) does not depend on Λ .

Proof. From the definition of H ∆ and the stochastic properties of the Poisson point process we have

P z,β Λ (dγ) = 1 Z z,β Λ z N Λ (γ) e -β H(γ) π Λ (dγ) = 1 Z z,β Λ z N ∆ (γ) e -β H ∆ (γ) z N Λ \∆ (γ) e -β H(γ Λ \∆ ) π ∆ (dγ ∆ )π Λ \∆ (dγ Λ \∆ ).
This expression ensures that the unnormalized conditional density of P z,β

Λ (dγ ∆ |γ ∆ c ) with respect to π ∆ (dγ ∆ ) is γ ∆ → z N ∆ (γ) e -β H ∆ (γ) . The normalization is necessary Z z,β ∆ (γ ∆ c
) and the proposition is proved.

The DLR equations give the local conditional marginal distributions of GPP. They are the main tool to understand the local description of P z,β Λ , in particular when Λ is large. Note that the local marginal distributions (not conditional) are in general not accessible. It is a difficult point of the theory of GPP. This fact will be reinforced in the infinite volume regime, where the local distributions can be non-unique.

The DLR equations have a major issue due the the intractable normalization constant Z z,β ∆ (γ ∆ c ). In the next section the problem is partially solved via the GNZ equations.

GNZ equations

The GNZ equations are due to Georgii, Nguyen and Zessin and have been introduced first in [START_REF] Nguyen | Integral and differential characterizations Gibbs processes[END_REF]. They generalize the Slivnyak-Mecke formulas for Poisson point processes. In this section we present and prove these equations. We need first to define the energy of a point inside a configuration. Definition 4. Let γ ∈ C f be a finite configuration and x ∈ R d . Then the local energy of x in γ is defined by

h(x, γ) = H({x} ∪ γ) -H(γ), with the convention +∞ -(+∞) = 0. Note that if x ∈ γ then h(x, γ) = 0.
Proposition 4 (GNZ equations). For any positive measurable function f from

R d × C f to R, ∑ x∈γ f (x, γ\{x})P z,β Λ (dγ) = z Λ f (x, γ)e -β h(x,γ) dxP z,β Λ (dγ). (10)
Proof. Let us decompose the left term in [START_REF] Daley | An introduction to the theory of point processes[END_REF].

∑ x∈γ f (x, γ\{x})P z,β Λ (dγ) = 1 Z z,β Λ ∑ x∈γ f (x, γ\{x})z N Λ (γ) e -β H(γ) π Λ (dγ) = e -λ d (Λ ) Z z,β Λ +∞ ∑ n=1 z n n! n ∑ k=1 Λ k f (x k , {x 1 , . . . , x n }\{x k })e -β H({x 1 ,...,x n }) dx 1 . . . dx n = e -λ d (Λ ) Z z,β Λ +∞ ∑ n=1 z n (n -1)! Λ k f (x, {x 1 , . . . , x n-1 })e -β H({x 1 ,...,x n-1 }) e -β h(x,{x 1 ,...,x n-1 }) dx 1 . . . dx n-1 dx = z Z z,β Λ Λ f (x, γ)z N Λ (γ) e -β H(γ) e -β h(x,γ) π Λ (dγ)dx = z Λ f (x, γ)e -β h(x,γ) dxP z,β Λ (dγ).
As usual the function f in ( 10) can be chosen without a constant sign. We just need to check that both terms in [START_REF] Daley | An introduction to the theory of point processes[END_REF] are integrable.

In the following proposition we show that the equations GNZ (10) characterize the probability measure

P z,β Λ . Proposition 5. Let Λ ⊂ R d bounded such that λ d (Λ ) > 0. Let P be a probability measure on C Λ such that for any positive measurable function f from R d × C f to R ∑ x∈γ f (x, γ\{x})P(dγ) = z Λ f (x, γ)e -β h(x,γ) dxP(dγ).
Then it holds that P = P z,β Λ .

Proof. Let us consider the measure

Q = 1 {H<+∞} z -N Λ e β H P. Then ∑ x∈γ f (x, γ\{x})Q(dγ) = ∑ x∈γ f (x, γ\{x})1 {H(γ)<+∞} z -N Λ (γ) e β H(γ) P(dγ) = z -1 ∑ x∈γ f (x, γ\{x})1 {H(γ\{x})<+∞} 1 {h(x,γ\{x})<+∞} z -N Λ (γ\{x}) e β H(γ\{x}) e β h(x,γ\{x}) P(dγ) = Λ f (x, γ)1 {H(γ)<+∞} 1 {h(x,γ)<+∞} e -β h(x,γ) z -N Λ (γ) e β H(γ) e β h(x,γ) dxP(dγ) = Λ f (x, γ)1 {h(x,γ)<+∞} dxQ(dγ).
We deduce that Q satisfies the Slivnyak-Mecke formula on {γ ∈ C Λ , H(γ) < +∞}. It is well-known (see [START_REF] Matthes | Infinitely divisible point processes[END_REF] for instance) that it implies that the measure Q (after normalization) is the Poisson point process π Λ restricted to {γ ∈ C Λ , H(γ) < +∞}. The proposition is proved.

These last two propositions show that the GNZ equations contain completely the informations on P z,β Λ . Note again that the normalization constant Z z,β Λ is not present in the equations.

Ruelle estimates

In this section we present Ruelle estimates in the context of superstable and lower regular energy functions. These estimates are technical and we refer to the original paper [START_REF] Ruelle | Superstable interactions in classical statistical mechanics[END_REF] for the proofs.

Definition 5. An energy function H is said superstable if H = H 1 + H 2 where H 1 is an energy function (see Definition (1)
) and H 2 is a pairwise energy function defined in (1) with a non-negative continuous pair potential ϕ such that ϕ(0) > 0. The energy function H is said lower regular if there exists a summable decreasing sequence of positive reals (ψ k ) k≥0 (i.e. ∑ +∞ k=0 ψ k < +∞) such that for any finite configurations

γ 1 and γ 2 H(γ 1 ∪ γ 2 ) -H(γ 1 ) -H(γ 2 ) ≥ -∑ k,k ∈Z d ψ k-k N 2 [k+[0,1] d ] (γ 1 ) + N 2 [k +[0,1] d ] (γ 2 ) . (11)
Let us give the main example of superstable and lower regular energy function.

Proposition 6 (Proposition 1.3 [START_REF] Ruelle | Superstable interactions in classical statistical mechanics[END_REF]). Let H be a pairwise energy function with a pair potential ϕ = ϕ 1 +ϕ 2 where ϕ 1 is stable and ϕ 2 is non-negative continuous with ϕ 2 (0) > 0. Moreover, we assume that there exists a positive decreasing function

ψ from R + to R such that +∞ 0 r d-1 ψ(r)dr < +∞
and such that for any x ∈ R, ϕ(x) ≥ -ψ( x ). Then the energy function H is superstable and lower regular.

In particular, the Lennard-Jones pair potential or the Strauss pair potential defined in Section 1.2 are superstable and lower regular. Note also that all geometric energy functions presented in Section 1.2 are not superstable.

Proposition 7 (corollary 2.9 [START_REF] Ruelle | Superstable interactions in classical statistical mechanics[END_REF]). Let H be a superstable and lower regular energy function. Let z > 0 and β > 0 be fixed. Then for any bounded subset ∆ ⊂ R d with λ d (∆ ) > 0 there exist two positive constants c 1 , c 2 such that for any bounded set Λ and k ≥ 0

P z,β Λ (N ∆ ≥ k) ≤ c 1 e -c 2 k 2 . ( 12
)
In particular, Ruelle estimates [START_REF] Dereudre | The existence of quermass-interaction processes for nonlocally stable interaction and nonbounded convex grains[END_REF] ensure that the random variable N ∆ admits exponential moments for all orders under P z,β Λ . Surprisingly, the variate N 2 ∆ admits exponential moments for small orders. This last fact is not true under the Poisson point process π z Λ = P z,0 Λ . The interaction between the points improves the integrability properties of the GPP with respect to the Poisson point process.

Infinite volume Gibbs point processes

In this section we present the theory of infinite volume GPP corresponding to the case "Λ = R d " of the previous section. Obviously, a definition inspired by (6) does not work since the energy of an infinite configuration γ is meaningless. A natural construction would be to consider a sequence of finite volume GPP (P

z,β Λ n ) n≥1 on bounded windows Λ n = [-n
, n] d and let n tend to infinity. It is more or less what we do in the following Sections 2.1 and 2.2, except that the convergence occurs only for a subsequence and that the field is stationarized (see equation ( 14)). As far as we know, there does not exist a general proof of the convergence of the sequence (P z,β Λ n ) n≥1 without extracted a subsequence. The stationarization is a convenient setting here in order to use the tightness entropy tools. In Sections 2.3 and 2.4 we prove that the accumulation points P z,β satisfy the DLR equations which is the standard definition of infinite volume GPP (see Definition 8). We make precise that the main new assumption in this section is the finite range property (see Definition 7). It means that the points interact with each other only if their distance is smaller than a fixed constant R > 0. The GNZ equations in the infinite volume regime are discussed in Section 2.5. The varitional characterisation of GPP, in the spirit of Proposition 1, is presented in Section 2.6. Uniqueness and non-uniqueness results of infinite volume GPP are treated in Sections 2.7 and 2.8. These results, whose proofs are completely self contained here, ensure the existence of a phase transition for the Area energy function presented in [START_REF] Baddeley | Area-interaction point processes[END_REF]. It means that the associated infinite volume Gibbs measures are unique for some parameters (z, β ) and non-unique for other parameters.

The local convergence setting

In this section we define the topology of local convergence which is the setting we use to prove the existence of an accumulation point for the sequence of finite volume Gibbs measures.

First, we say that a function from C to R is local if there exists a bounded set

∆ ⊂ R d such that for all γ ∈ C , f (γ) = f (γ ∆ ).
Definition 6. The local convergence topology on the space of probability measures on C is the smallest topology such that for any local bounded function f from C to R the function P → f dP is continuous. We denote by τ L this topology.

Let us note that the continuity of functions f in the previous definition is not required. For instance the function γ → f (γ) = 1 N ∆ (γ)≥k , where ∆ is a bounded set in R d and k any integer, is a bounded local function. For any vector u ∈ R d we denote by τ u the translation by the vector u acting on R d or C . A probability P on C is said stationary (or shift invariant) if for any vector u ∈ R d P = P • τ -1 u . Our tightness tool is based on the specific entropy which is defined for any stationary probability P on C by

I ζ (P) = lim n→+∞ 1 λ d (Λ n ) I(P Λ n |π ζ Λ n ), (13)
where

I(P Λ n |π ζ Λ n )
is the relative entropy of P Λ n , the projection of P on Λ n , with respect to π ζ Λ n (see Definition [START_REF] Chayes | The analysis of the Widom-Rowlinson model by stochastic geometric methods[END_REF]). Note that the specific entropy I ζ (P) always exists (i.e. the limit in (13) exists); see chapter 15 in [START_REF] Georgii | Gibbs measures and phase transitions[END_REF]. The tightness tool presented in Lemma 1 below is a consequence of the following proposition.

Proposition 8 (Proposition 15.14 [START_REF] Georgii | Gibbs measures and phase transitions[END_REF]). For any ζ > 0 and any value K ≥ 0, the set {P ∈ P such that I ζ (P) ≤ K} is sequentially compact for the topology τ L , where P is the space of stationary probability measures on C with finite intensity.

An accumulation point P z,β

In this section we prove the existence of an accumulation point for a sequence of stationarized finite volume GPP. To the end we consider the Gibbs measures (P

z,β Λ n ) n≥1 on Λ n := [-n, n] d , where (P z,β
Λ ) is defined in (6) for any z > 0, β ≥ 0 and energy function H. We assume that H is stationary, which means that for any vector u ∈ R d and any finite configuration

γ ∈ C f H(τ u (γ)) = H(γ).
For any n ≥ 1, the empirical field Pz,β Λ n is defined by the probability measure on C such that for any test function

f Contents f (γ) Pz,β Λ n (dγ) = 1 λ d (Λ n ) Λ n f (τ u (γ))P z,β Λ n (dγ)du. (14)
The probability measure Pz,β Λ n can be interpreted as the Gibbs measure P z,β

Λ n where the origin of the space R d (i.e. the point {0}) is replaced by a random point chosen uniformly inside Λ n . It is a kind of stationarization of P z,β Λ n and any accumulation point of the sequence ( Pz,β Λ n ) n≥1 is necessary stationary.

Proposition 9. The sequence ( Pz,β Λ n ) n≥1 is tight for the τ L topology. We denote by P z,β any of its accumulation points.

Proof. Our tightness tool is the following lemma whose the proof is a consequence of Proposition 8 (See also Proposition 15.52 in [START_REF] Georgii | Gibbs measures and phase transitions[END_REF]).

Lemma 1. The sequence ( Pz,β Λ n ) n≥1 is tight for the τ L topology if there exits ζ > 0 such that sup n≥1 1 λ d (Λ n ) I(P z,β Λ n |π ζ Λ n ) < +∞. (15)
So, let us compute

I(P z,β Λ n |π ζ Λ n )
and check that we can find ζ > 0 such that (15) holds.

I(P z,β Λ n |π ζ Λ n ) = log dP z,β Λ n dπ ζ Λ n dP z,β Λ n = log dP z,β Λ n dπ Λ n + log dπ Λ n dπ ζ Λ n dP z,β Λ n = log z N Λn e -β H Z z,β Λ n + log e (ζ -1)λ d (Λ n ) 1 ζ N Λn dP z,β Λ n = -β H + log z ζ N Λ n dP z,β Λ n + (ζ -1)λ d (Λ n ) -log(Z z,β Λ n ).
Thanks to the non degeneracy and the stability of H we find that

I(P z,β Λ n |π ζ Λ n ) ≤ -Aβ + log z ζ N Λ n dP z,β Λ n + λ d (Λ n ) (ζ -1) + 1 + β H({ / 0}) .
Choosing ζ > 0 such that -Aβ + log(z/ζ ) ≤ 0 we obtain

I(P z,β Λ n |π ζ Λ n ) ≤ λ d (Λ n )(ζ + β H({ / 0})
and ( 15) holds. Proposition 9 is proved.

In the following, for sake of simplicity, we say that Pz,β Λ n converges to P z,β although it occurs only for a subsequence.

Note that the existence of an accumulation points holds under very weak assumptions on the energy function H. Indeed the two major assumptions are the stability and the stationarity. The superstability or the lower regularity presented in Definition 5 are not required here. However, if the energy function H is superstable and lower regular, then the accumularion points P z,β inherits Ruelle estimates [START_REF] Dereudre | The existence of quermass-interaction processes for nonlocally stable interaction and nonbounded convex grains[END_REF]. This fact is obvious since the function γ → 1 {N ∆ (γ)≥k} is locally bounded.

Corollary 1. Let H be a superstable and lower regular energy function (see Definition 5). Let z > 0 and β > 0 be fixed. Then for any bounded subset ∆ ⊂ R d with λ d (∆ ) > 0, there exists c 1 and c 2 two positive constants such that for any k ≥ 0

P z,β (N ∆ ≥ k) ≤ c 1 e -c 2 k 2 . ( 16
)
The important point now is to prove that P z,β satisfies good stochastic properties as for instance the DLR or GNZ equations. At this stage, without extra assumptions, these equations are not necessarily satisfied. Indeed it is possible to build energy functions H such that the accumulation point P z,β is degenerated and charges only the empty configuration. In this mini-course our extra assumption is the finite range property presented in the following section. More general settings have been investigated for instance in [START_REF] Dereudre | Existence of Gibbsian point processes with geometry-dependent interactions[END_REF] or [START_REF] Ruelle | Statistical mechanics: Rigorous results[END_REF].

The finite range property

The finite range property expresses that further a certain distance distance R > 0 the points do not interact each other. Let us recall the Minkoswki ⊕ operator acting on sets in R d . For any two sets A, B ⊂ R d , the set A ⊕ B is defined by {x + y, x ∈ A and y ∈ B}.

∈ C f H ∆ (γ) := H(γ) -H(γ ∆ c ) = H(γ ∆ ⊕B(0,R) ) -H(γ ∆ ⊕B(0,R)\∆ c ).
Let us illustrate the finite range property in the setting of pairwise interaction defined in [START_REF] Baddeley | Variational estimators for the parameters of Gibbs point process models[END_REF]. Assume that the interaction potential ϕ : R + → R ∪ {+∞} has a support included in [0, R]. Then the associated energy function has a finite R;

Contents H ∆ (γ) = ∑ {x, y} ⊂ γ {x, y} ∩ ∆ = / 0 |x -y| ≤ R ϕ(|x -y|) = ∑ {x, y} ⊂ γ ∆ ⊕B(0,R) {x, y} ∩ ∆ = / 0 ϕ(|x -y|).
Also the area energy function (4) inherits the finite range property. A simple computation gives

H ∆ (γ) = Area   x∈γ ∆ B(x, R)\ x∈γ ∆ ⊕B(0,2R)\∆ B(x, R)   (17)
which provides a range of interaction equals to 2R. Let us note that the energy functions defined in ( 2),( 3) and ( 5) do not have the finite range property. Similarly the pairwise energy function [START_REF] Baddeley | Variational estimators for the parameters of Gibbs point process models[END_REF] with the Lennard-Jones potential is not finite range since the support of the pair potential is not bounded. A truncated version of such potential is sometimes considered.

Let us finish this section by noting that the finite range property allows to extend the domain of definition of H ∆ from the space C f to the set C . Indeed, since

H ∆ (γ) = H ∆ (γ ∆ ⊕B(0,R) )
, this equality provides a definition of H ∆ (γ) when γ is in C . This point is crucial in order to correctly define the DLR equations in the infinite volume regime.

DLR equations

In section 1 on the finite volume GPP, the DLR equations are presented as properties for P z,β Λ (see Section 1.4). In the setting of infinite volume GPP, the DLR equations are the main points of the definition of GPP. Definition 8 (infinite volume GPP). Let H be a stationary and finite range energy function. A stationary probability P on C is an infinite volume Gibbs measure with activity z > 0, inverse temperature β ≥ 0 and energy function H if for any bounded

∆ ⊂ R d such that λ d (∆ ) > 0 then for P-a.s. all γ ∆ c P(dγ ∆ |γ ∆ c ) = 1 Z z,β ∆ (γ ∆ c ) z N ∆ (γ) e -β H ∆ (γ) π ∆ (dγ ∆ ), (18)
where

Z z,β ∆ (γ ∆ c ) is the normalizing constant z N ∆ (γ) e -β H ∆ (γ) π ∆ (dγ ∆ ).
As usual, an infinite volume GPP is a point process whose distribution is an infinite volume Gibbs measure.

Note that the DLR equations ( 18) make sense since H ∆ (γ) is well defined for any configuration γ ∈ C (see the end of Section 2.3). Note also that the DLR equations ( 18) can be reformulated in an integral form. Indeed P satisfies [START_REF] Dobrushin | A criterion of the uniqueness of Gibbsian fields in the noncompact case[END_REF] if and only if for any local bounded function

f from C to R f dP = f (γ ∆ ∪ γ ∆ c ) 1 Z z,β ∆ (γ ∆ c ) z N ∆ (γ ∆ ) e -β H ∆ (γ ∆ ∪γ ∆ c ) π ∆ (dγ ∆ )P(dγ). (19)
The term "equation" is now highlighted by the formulation [START_REF] Fiksel | Estimation of parametrized pair potentials of marked and nonmarked Gibbsian point processes[END_REF] since the unknown variate P appears in both left and right sides. The existence, uniqueness and non-uniqueness of solutions of such DLR equations are non trivial questions. In the next theorem, we show that the accumulation point P z,β obtained in Section 2.2 is such a solution. Infinite volume Gibbs measure exist and the question of existence is solved. The uniqueness and non-uniqueness are discussed in Sections 2.7 and 2.8.

Theorem 1. Let H be a stationary and finite range energy function. Then for any z > 0 and β ≥ 0 the probability measure P z,β defined in Proposition 9 is an infinite volume Gibbs measure.

Proof. We have just to check that P z,β satisfies, for any bounded ∆ and any positive local bounded function f , the equation [START_REF] Fiksel | Estimation of parametrized pair potentials of marked and nonmarked Gibbsian point processes[END_REF]. Let us define the function f ∆ by

f ∆ : γ → f (γ ∆ ∪ γ ∆ c ) 1 Z z,β ∆ (γ ∆ c ) z N ∆ (γ ∆ ) e -β H ∆ (γ ∆ ∪γ ∆ c ) π ∆ (dγ ∆ ).
Since f is local and bounded and since H is finite range, the function f ∆ is bounded and local as well. From the convergence of the sequence ( Pz,β Λ n ) n≥1 to P z,β with respect to the τ L topology, we have

f ∆ dP z,β = lim n→∞ f ∆ d Pz,β Λ n = lim n→∞ 1 λ d (Λ n ) Λ n f ∆ (τ u (γ))P z,β Λ n (dγ)du. = lim n→∞ 1 λ d (Λ n ) Λ n f (γ ∆ ∪ τ u (γ) ∆ c ) z N ∆ (γ ∆ ) Z z,β ∆ (τ u (γ) ∆ c ) e -β H ∆ (γ ∆ ∪τ u (γ) ∆ c ) π ∆ (dγ ∆ )P z,β Λ n (dγ)du = lim n→∞ 1 λ d (Λ n ) Λ n f τ u γ τ -u (∆ ) ∪ γ τ -u (∆ ) c z N τ -u (∆ ) (γ τ -u (∆ ) ) Z z,β τ -u (∆ ) (γ τ -u (∆ ) c ) e -β H τ -u (∆ ) γ τ -u (∆ ) ∪γ τ -u (∆ ) c π τ -u (∆ ) (dγ τ -u (∆ ) )P z,β Λ n (dγ)du. ( 20
)
Denoting by Λ * n the set of u ∈ Λ n such that τ -u (∆ ) ⊂ Λ n , by Proposition 3, P z,β

Λ n satisfies the DLR equation on τ -u (∆ ) as soon as τ -u (∆ ) ⊂ Λ n (i.e. u ∈ Λ * n ). It follows that for any u ∈ Λ * n f (τ u γ)P z,β Λ n (dγ) = f τ u γ τ -u (∆ ) ∪ γ τ -u (∆ ) c z N τ -u (∆ ) (γ τ -u (∆ ) ) Z z,β τ -u (∆ ) (γ τ -u (∆ ) c ) e -β H τ -u (∆ ) γ τ -u (∆ ) ∪γ τ -u (∆ ) c π τ -u (∆ ) (dγ τ -u (∆ ) )P z,β Λ n (dγ). (21)
By noting that

λ d (Λ * n ) is equivalent to λ d (Λ n )
when n goes to infinity, we obtain in compiling ( 20) and ( 21)

f ∆ dP z,β = lim n→∞ 1 λ d (Λ n ) Λ * n f (τ u γ)P z,β Λ n (dγ)du = lim n→∞ f (γ) Pz,β Λ n (dγ) = f dP z,β
which gives the expected integral DLR equation on ∆ with test function f .

GNZ equations

In this section we deal with the GNZ equations in the infinite volume regime. As in the finite volume case, the main advantage of such equations is that the intractable normalization factor Z z,β Λ is not present. Note first that, in the setting of finite range interaction R > 0, the local energy h(x, γ) defined in Definition 4 is well-defined for any configuration γ ∈ C even if γ is infinite. Indeed, we clearly have h(x, γ) = h(x, γ B(x,R) ).

Theorem 2. Let P be a probability measure on C . Let H be a finite range energy function and z > 0, β ≥ 0 be two parameters. Then P is an infinite volume Gibbs measure with energy function H, activity z > 0 and inverse temperature β if and only if for any positive measurable function f from

R d × C to R ∑ x∈γ f (x, γ\{x})P(dγ) = z R d f (x, γ)e -β h(x,γ) dxP(dγ). (22)
Proof. Let us start with the proof of the "only if" part. Let P be an infinite volume Gibbs measure. By standard monotonicity arguments it is sufficient to prove [START_REF] Georgii | Gibbs measures and phase transitions[END_REF] for any local positive measurable function f . So let ∆ ⊂ R d be a bounded set such that f (x, γ) = 1 ∆ (x) f (x, γ ∆ ). Applying now the DLR equation ( 19) on the set ∆ we find

∑ x∈γ f (x, γ\{x})P(dγ) = ∑ x∈γ ∆ f (x, γ ∆ \{x}) 1 Z z,β ∆ (γ ∆ c ) z N ∆ (γ ∆ ) e -β H ∆ (γ ∆ ∪γ ∆ c ) π ∆ (dγ ∆ )P(dγ).
By computations similar to those developed in the proof of Proposition 4, we obtain

∑ x∈γ f (x, γ\{x})P(dγ) = z ∆ f (x, γ ∆ ) 1 Z z,β ∆ (γ ∆ c ) e -β h(x,γ ∆ ∪γ ∆ c ) z N ∆ (γ ∆ ) e -β H ∆ (γ ∆ ∪γ ∆ c ) π ∆ (dγ ∆ )dxP(dγ) = z R d f (x, γ)e -β h(x,γ) dxP(dγ).
Let us now turn to the "if part". Applying equation ( 22) to the function f (x, γ) = ψ(γ ∆ c ) f (x, γ) where f is a local positive function with support ∆ and ψ a positive test function we find

ψ(γ ∆ c ) ∑ x∈γ ∆ f (x, γ\{x})P(dγ) = z ψ(γ ∆ c ) R d f (x, γ)e -β h(x,γ) dxP(dγ).
This implies that for P almost all γ ∆ c the conditional probability measure P(dγ ∆ |γ ∆ c ) solves the GNZ equations on ∆ with local energy function γ ∆ → h(x, γ ∆ ∪ γ ∆ c ). Following an adaptation of the proof of Proposition 5, we get that

P(dγ ∆ |γ ∆ c ) = 1 Z z,β ∆ (γ ∆ c ) z N ∆ (γ) e -β H ∆ (γ) π ∆ (dγ ∆ ),
which is exactly the DLR equation ( 18) on ∆ . The theorem is proved.

Let us finish this section with an application of the GNZ equations which highlights that some properties of infinite volume GPP can be extracted from the implicit GNZ equations.

Proposition 10. Let Γ be a infinite volume GPP for the hardcore pairwise interaction ϕ(r) = +∞1 [0,R] (r) (see Definition (1)) and the activity z > 0. Then

z 1 + zv d R d ≤ E N [0,1] d (Γ ) ≤ z, ( 23
)
where v d is the volume of the unit ball in R d .

Note that the inverse temperature β does not play any role here and that E P N [0,1] d (Γ ) is simply the intensity of Γ .

Contents

Proof. The local energy of such harcore pairwise interaction is given by

h(x, γ) = ∑ y∈γ B(x,R) ϕ(|x -y|) = +∞1 γ B(x,R) = / 0 .
So the GNZ equation ( 22) with the function

f (x, γ) = 1 [0,1] d (x) gives E N [0,1] d (Γ ) = z [0,1] d P(Γ B(x,R) = / 0)dx = z P(Γ B(0,R) = / 0),
which provides a relation between the intensity and the spherical contact distribution of Γ . The upper bound in [START_REF] Georgii | Stochastic comparison of point random fields[END_REF] follows. For the lower bound we have

E P N [0,1] d (Γ ) = z P(Γ B(0,R) = / 0) ≥ z 1 -E P N B(0,R) (Γ ) = z 1 -v d R r E P N [0,1] d (Γ ) .
Note also that a natural upper bound for E P N [0,1] d is obtained via the closed packing configuration. For instance, in dimension d = 2, it gives the upper bound π/(2 √ 3R 2 ).

Variational principle

In this section, we extend the variational principle for finite volume GPP presented in Proposition 1 to the setting of infinite volume GPP. For brevity we present only the result without the proof which can be found in [START_REF] Dereudre | Variational principle for Gibbs point processes with finite range interaction[END_REF].

The variational principle claims that the Gibbs measures are the minimizers of the free excess energy defined by the sum of the the mean energy and the specific entropy. Moreover, the minimum is equal to minus the pressure. Let us first define all these macroscopic quantities.

Let us start by introducing the pressure with free boundary condition. It is defined as the following limit

p z,β := lim n→+∞ 1 |Λ n | ln(Z z,β Λ n ), (24)
The existence of such limit is proved for instance in Lemma 1 in [START_REF] Dereudre | Variational principle for Gibbs point processes with finite range interaction[END_REF].

The second macroscopic quantity involves the mean energy of a stationary probability measure P. It is also defined by a limit but, in opposition to the pressure, we have to assume that it exists. The proof of such existence is generally based on stationary arguments and nice representations of the energy contribution per unit volume. It depends strongly on the expression of the energy function H. Examples are given below. So for any stationary probability measure P on C we assume that the following limit exists in R ∪ {+∞},

H(P) := lim n→∞ 1 |Λ n | H(γ Λ n )dP(γ), (25)
and we call the limit mean energy of P.

We need to introduce a technical assumption on the boundary effects of H. We assume that for any infinite volume Gibbs measure P

lim n→∞ 1 |Λ n | ∂ H Λ n (γ)dP(γ) = 0, (26)
where

∂ H Λ n (γ) = H Λ n (γ) -H(γ Λ n ).
Theorem 3 (Variational Principle, Theorem 1, [START_REF] Dereudre | Variational principle for Gibbs point processes with finite range interaction[END_REF]).

We assume that H is stationary and finite range. Moreover, we assume that the mean energy exists for any stationary probability measure P (i.e. the limit (25) exists) and that the boundary effects assumption (26) holds. Let z > 0 and β ≥ 0 two parameters. Then for any stationary probability measure P on C with finite intensity

I 1 (P) + β H(P) -log(z)E P (N [0,1] d ) ≥ -p z,β , (27)
with equality if and only if P is a Gibbs measure with activity z > 0, inverse temperature β and energy function H.

Let us finish this section by presenting the two fundamental examples of energy functions satisfying the assumptions of Theorem 3.

Proposition 11. Let H be the Area energy function defined in (4). Then both limits (25) and (26) exist. In particular, the assumptions of Theorem 3 are satisfied and the variational principle holds.

Proof. Let us prove only that the limit (25) exists. The existence of limit ( 26) can be shown in the same way. By definition of H and the stationarity of P,

H(γ Λ n )P(dγ) = Area(L R (γ Λ n ))P(dγ) = λ d (Λ n ) Area(L R (γ) ∩ [0, 1] d )P(dγ) + Area(L R (γ Λ n )) -Area(L R (γ) ∩ [-n, n] d ) P(dγ). ( 28
)
By geometric arguments, we get that

Area(L R (γ Λ n )) -Area(L R (γ) ∩ [-n, n] d ) ≤ Cn d-1 ,
for some constant C > 0. We deduce that the limit (25) exists with

H(P) = Area(L R (γ) ∩ [0, 1] d )P(dγ).
Proposition 12. Let H be the pairwise energy function defined in (1) with a superstable, lower regular pair potential with compact support. Then the both limits (25) and (26) exist. In particular the assumptions of Theorem 3 are satisfied and the variational principle holds.

Proof. Since the potential ϕ is stable with compact support, we deduce that ϕ ≥ 2A and H is finite range and lower regular. In this setting, the existence of the limit ( 25) is proved in [START_REF] Georgii | Large deviations and the equivalence of ensembles for Gibbsian particle systems with superstable interaction[END_REF], Theorem 1 with

H(P) = 1 2 ∑ 0 =x∈γ ϕ(x)P 0 (dγ) if E P (N 2 [0,1] d ) < ∞ +∞ otherwise ( 29
)
where P 0 is the Palm measure of P. Recall that P 0 can be viewed as the natural version of the conditional probability P(.|0 ∈ γ) (see [START_REF] Matthes | Infinitely divisible point processes[END_REF] for more details). It remains to prove the existence of the limit [START_REF] Geyer | Simulation procedures and likelihood inference for spatial point processes[END_REF] for any Gibbs measure P on C . A simple computation gives that, for any γ ∈ C ,

∂ H Λ n (γ) = ∑ x∈γ Λ ⊕ n \Λn ∑ y∈γ Λn\Λ n ϕ(x -y),
where Λ ⊕ n = Λ n+R 0 and Λ n = Λ n-R 0 with R 0 an integer larger than the range of the interaction R.

Therefore thanks to the stationarity of P and the GNZ equations ( 22), we obtain

∂ H Λ n (γ)dP(γ) ≤ ∑ x∈γ Λ ⊕ n \Λn ∑ y∈γ\{x} |ϕ(x -y)|dP(γ) = z Λ ⊕ n \Λ n e -β ∑ y∈γ ϕ(x-y) ∑ y∈γ |φ (x -y)|dxdP(γ) = z|Λ ⊕ n \Λ n | e -β ∑ y∈γ B(0,R 0 ) ϕ(y) ∑ y∈γ B(0,R 0 ) |ϕ(y)|dP(γ). Since ϕ ≥ 2A, denoting by C := sup c∈[2A;+∞) |c|e -β c < ∞ we find that ∂ H Λ n (γ)dP(γ) ≤ zC|Λ ⊕ n \Λ n | N B(0,R 0 ) (γ)e -2β AN B(0,R 0 ) (γ) dP(dγ). ( 30
)
Using Ruelle estimates ( 16), the integral in the right term of ( 30) is finite. The boundary assumption (26) follows.

A uniqueness result

In this section we investigate the uniqueness of infinite volume Gibbs measures. The common belief claims that the Gibbs measures are unique when the activity z or (and) the inverse temperature β are small enough (low activity, high temperature regime). The non-uniqueness phenomenon (discussed in the next section) are in general related to some issues with the energy part in the variational principle (see Theorem 3). Indeed, either the mean energy has several minimizers or there is a conflict between the energy and the entropy. Therefore it is natural to expect the the Gibbs measures are unique when β is small enough. When z is small, the mean number of points per unit volume is low and so the energy is in general low as well.

As far as we know, there do not exist general results which prove the uniqueness for small β or small z. In the case of pairwise energy functions (1), the uniqueness for any β > 0 and z > 0 small enough is proved via the Kirkwood-Salsburg equations (see Theorem 5.7 [START_REF] Ruelle | Superstable interactions in classical statistical mechanics[END_REF]). An extension of the Dobrushin uniqueness criterium in the continuum is developed as well [START_REF] Dobrushin | A criterion of the uniqueness of Gibbsian fields in the noncompact case[END_REF]. The uniqueness of GPP can also be obtained via the cluster expansion machinery which provides a power series expansion of the partition function when z and β are small enough. This approach has been introduced first by Mayer and Montroll [START_REF] Mayer | Molecular distributions[END_REF] and we refer to [START_REF] Poghosyan | Abstract cluster expansion with applications to statistical mechanical systems[END_REF] for a general presentation.

In this section we give a simple and self-contained proof of the uniqueness of GPP for all β ≥ 0 and any z > 0 small enough. We just assume that the energy function H has a local energy h uniformly bounded from below. This setting covers for instance the case of pairwise energy function [START_REF] Baddeley | Variational estimators for the parameters of Gibbs point process models[END_REF] with non-negative pair potential or the Area energy function [START_REF] Baddeley | Area-interaction point processes[END_REF].

Let us start by recalling the existence of a percolation threshold for the Poisson Boolean model. For any configuration γ ∈ C the percolation of L R (γ) = ∪ x∈γ B(x, R) means the existence of an unbounded connected component in L R (γ) .

Proposition 13 (Theorem 1 [START_REF] Hall | On continuum percolation[END_REF]). For any d ≥ 2, there exists 0 < z d < +∞ such that for z < z d , π z (L 1/2 percolates ) = 0 and for z > z d , π z (L 1/2 percolates ) = 1.

The value z d is called the percolation threshold of the Poisson Boolean model with radius 1/2. By scale invariance, the percolation threshold for any other radius R is simply z d /(2R) d . The exact value of z d is unknown but numerical studies provide for instance the approximation z 2 1.4 in dimension d = 2. Theorem 4. Let H be an energy function with finite range R > 0 such that the local energy h is uniformly bounded from below by a constant C. Then for any β ≥ 0 and z < z d e Cβ /R d , there exists an unique Gibbs measure with energy function H, activity z > 0 and inverse temperature β .

Proof. The proof is based on two main ingredients. The first one is the stochastic domination of Gibbs measures, with uniformly bounded from below local energy function h, by Poisson processes. This result is given in the following lemma, whose proof can be found in [START_REF] Georgii | Stochastic comparison of point random fields[END_REF]. The second ingredient is a disagreement percolation result presented in Lemma 3 below. Lemma 2. Let H be an energy function such that the local energy h is uniformly bounded from below by a constant C. Then for a any bounded set ∆ and any outside configuration γ ∆ c the Gibbs distribution inside ∆ given by

P z,β (dγ ∆ |γ ∆ c ) = 1 Z z,β ∆ (γ ∆ c ) z N ∆ (γ) e -β H ∆ (γ) π ∆ (dγ ∆ )
is stochastically dominated by the Poisson point distribution π ze -Cβ ∆ (dγ ∆ ).

Thanks to Strassen's Theorem, this stochastic domination can be interpreted via the following coupling (which could be the definition of the stochastic domination): There exist two point processes Γ and Γ on ∆ such that Γ ⊂ Γ , Γ ∼ P z,β (dγ ∆ |γ ∆ c ) and Γ ∼ π ze -Cβ ∆ (dγ ∆ ). Now the rest of the proof of Theorem 4 consists in showing that the Gibbs measure is unique as soon as π ze -Cβ (L R/2 percolates) = 0. Roughly speaking, if the dominating process does not percolate, the information coming from the boundary condition does not propagate in the heart of the model and the Gibbs measure is unique. To prove rigorously this phenomenon, we need a disagreement percolation argument introduced first in [START_REF] Van Den Berg | Disagreement percolation in the study of Markov fields[END_REF]. For any sets A, B ∈ R d , we denote by A B the set

(A c ⊕ B) c . Lemma 3. Let γ 1
∆ c and γ 2 ∆ c be two configurations on ∆ c . For any R > R, there exist three point processes

Γ 1 , Γ 2 and Γ on ∆ such that Γ 1 ⊂ Γ , Γ 2 ⊂ Γ , Γ 1 ∼ P z,β (dγ ∆ |γ 1 ∆ c ), Γ 2 ∼ P z,β (dγ ∆ |γ 2 ∆ c ) and Γ ∼ π ze -Cβ ∆ (dγ ∆ ). Moreover, denoting by L ∆ R /2 (Γ ) the connected components of L R /2 (Γ ) which are inside ∆ B(0, R /2), then Γ 1 = Γ 2 on the set L ∆ R /2 (Γ ).
Proof. Let us note first that, by Lemma 2, there exist three point processes Γ 1 , Γ 2 and Γ on ∆ such that

Γ 1 ⊂ Γ , Γ 2 ⊂ Γ , Γ 1 ∼ P z,β (dγ ∆ |γ 1 ∆ c ), Γ 2 ∼ P z,β (dγ ∆ |γ 2 ∆ c ) and Γ ∼ π ze -Cβ ∆ (dγ ∆ ).
The main difficulty is now to show that we can build Γ 1 and Γ 2 such that Γ 1 = Γ 2 on the set L ∆ R /2 (Γ ). Let us decompose ∆ via a grid of small cubes where each cube has a diameter smaller than ε = (R -R)/2. We define an arbitrary numeration of these cubes (C i ) 1≤i≤m and we construct progressively the processes Γ 1 , Γ 2 and Γ on each cube C i . Assume that they are already constructed on C I := ∪ i∈I C i with all the expected properties:Γ 1

C I ⊂ Γ C I , Γ 2 C I ⊂ Γ C I , Γ 1 C I ∼ P z,β (dγ C I |γ 1 ∆ c ), Γ 2 C I ∼ P z,β (dγ C I |γ 2 ∆ c ), Γ C I ∼ π ze -Cβ C I (dγ C I ) and Γ 1 C I = Γ 2 C I on the set L ∆ R /2 (Γ C I ). Let us consider the smaller index j ∈ {1, . . . m}\I such that either the distances d(C j , γ 1 ∆ c ) or d(C j , γ 2 ∆ c ) or d(C j ,Γ C I ) is smaller than R -ε.
• If such an index j does not exist, by the finite range property the following Gibbs distributions coincide on ∆ I = ∆ \C I ;

P z,β (dγ ∆ I |γ 1 ∆ c ∪ Γ 1 C I ) = P z,β (dγ ∆ I |γ 2 ∆ c ∪ Γ 2 C I ).
Therefore we define Γ 1 , Γ 2 and Γ on ∆ I by considering Γ 1 ∆ I and Γ ∆ I as in Lemma 2 and by putting Γ 2 ∆ I = Γ 1 ∆ I . We can easily check that all expected properties hold and the full construction of Γ 1 , Γ 2 and Γ is over.

• If such an index j does exist, we consider the double coupling construction of Γ 1 , Γ 2 and Γ on ∆ I . It means that

Γ 1 ∆ I ⊂ Γ ∆ I , Γ 2 ∆ I ⊂ Γ ∆ I , Γ 1 ∆ I ∼ P z,β (dγ ∆ I |γ 1 ∆ c ∪ Γ 1 C I ), Γ 2 ∆ I ∼ P z,β (dγ ∆ I |γ 2 ∆ c ∪ Γ 2 C I ) and Γ ∆ I ∼ π ze Cβ ∆ I (dγ ∆
). Now we keep these processes Γ 1 ∆ I , Γ 2 ∆ I and Γ ∆ I only on the window C j . The construction of the processes Γ 1 , Γ 2 and Γ is now over C I ∪ C j and we can check again that all expected properties hold. We go on to the construction of the processes on a new cube in (C i ) i∈{1,...n}\{I, j} and so on.

Let us now finish the proof of Theorem 4 by considering two infinite volume GPP Γ 1 and Γ 2 with distribution P 1 and P 2 . We have to show that for any local event A P 1 (A) = P 2 (A). We denote by ∆ 0 the support of such an event A. Let us consider a bounded subset ∆ ⊃ ∆ 0 and three new processes Γ 1 ∆ , Γ 2 ∆ and Γ ∆ on ∆ constructed as in Lemma 3. Precisely, for any i = 1, 2

Γ i ∆ ⊂ Γ ∆ , Γ ∼ π ze -Cβ ∆ , the conditional distribution of Γ i ∆ given Γ i ∆ c is P z,β (| Γ i ∆ c ) and Γ 1 ∆ = Γ 2 ∆ on the set L ∆ R /2 (Γ ∆ ). The parameter R > R is chosen such that ze -Cβ R d < z d (31)
which is possible by assumption on z.

Thanks to the DLR equations ( 18), for any i = 1, 2 the processes Γ i ∆ and Γ i ∆ have the same distributions and therefore P i (A) = P(Γ i ∆ ∈ A). Denoting by {∆ ↔ ∆ 0 } the event that there exists a connected component in L R /2 (Γ ) which intersects (∆ B(0, R /2)) c and ∆ 0 , we obtain that

|P 1 (A) -P 2 (A)| = |P(Γ 1 ∆ ∈ A) -P(Γ 2 ∆ ∈ A)| ≤ E 1 {∆ ↔∆ 0 } 1 Γ 1 ∆ ∈A -1 Γ 2 ∆ ∈A + E 1 {∆ ↔∆ 0 } c 1 Γ 1 ∆ ∈A -1 Γ 2 ∆ ∈A ≤ P({∆ ↔ ∆ 0 }) + E 1 {∆ ↔∆ 0 } c 1 Γ 1 ∆ ∈A -1 Γ 1 ∆ ∈A = P({∆ ↔ ∆ 0 }). ( 32
)
By the choice of R in inequality [START_REF] Jensen | Pseudolikelihood for exponential family models of spatial point processes[END_REF] and Proposition 13, it follows that

π ze -Cβ L R /2 percolates = 0
and we deduce, by a monotonicity argument, the probability P({∆ ↔ ∆ 0 }) tends to 0 when ∆ tends to R d (see [START_REF] Meester | Continuum percolation, volume 119 of Cambridge Tracts in Mathematics[END_REF] for details on equivalent characterizations of continuum percolation). The left term in [START_REF] Kendall | Perfect simulation using dominating processes on ordered spaces, with application to locally stable point processes[END_REF] does not depend on ∆ and therefore it is null. Theorem 4 is proved.

Contents

A non-uniqueness result

In this section we discuss the non-uniqueness phenomenon of infinite volume Gibbs measures. It is believed to occur for almost all models provided that the activity z or the inverse temperature β is large enough. However, in the present continuous setting without spin, it is only proved for few models and several old conjectures are still valid. For instance, for the pairwise Lennard-Jones interaction defined in (1), it is conjectured that for β large (but not too large) there exists an unique z such that the Gibbs measures are not unique. It would correspond to a liquid-vapour phase transition. Similarly for β very large, it is conjectured that the non-uniqueness occurs as soon as z is larger than a threshold z β . It would correspond to a crystallization phenomenon for which a symmetry breaking may occur. Indeed, it is expected, but not proved at all, that some continuum Gibbs measures would be not invariant under symmetries like translations, rotations, etc. This conjecture is probably one of the most important and difficult challenges in statistical physics. In all cases, the nonuniqueness appear when the local distribution of infinite volume Gibbs measures depend on the boundary conditions "at infinity".

In this section we give a complete proof of such non-uniqueness result for the Area energy interaction presented in (4). This result has been first proved in [START_REF] Widom | New model for the study of liquid-vapor phase transitions[END_REF] but our proof is inspired by the one given in [START_REF] Chayes | The analysis of the Widom-Rowlinson model by stochastic geometric methods[END_REF]. Roughly speaking, we build two different Gibbs measures which depend, via a percolation phenomenon, on the boundary conditions "at infinity". In one case, the boundary condition "at infinity" is empty and in the other case the boundary condition is full of particles. We show that the intensity of both infinite volume Gibbs measures are different.

Let us cite another famous non-uniqueness result for attractive pair and repulsive four-body potentials [START_REF] Lebowitz | Liquid-vapor phase transitions for systems with finite-range interactions[END_REF]. As far as we know, this result and the one presented below on the Area interaction, are the only rigorous proofs of non-uniqueness results for continuum particles systems without spin.

Theorem 5. For z = β large enough, the infinite volume Gibbs measures for the Area energy function H presented in (4), the activity z and the inverse temperature β are not unique.

Proof. In all the proof we fix z = β . Let us consider following finite volume Gibbs measures on Λ n = [-n, n] d with different boundary conditions:

dP Λ n (γ) = 1 Z Λ n 1 γ Λn\Λ n = / 0 z N Λn (γ) e -zArea Λ n ∩L R (γ) dπ Λ n (γ),
and

dQ Λ n (γ) = 1 Z Λ n z N Λn (γ) e -zArea Λ n ∩L R (γ) dπ Λ n (γ),
where

Λ n = Λ n B(0, R/2). Recall that R is the radius of balls in L R (γ) = ∪ x∈γ B(x, R
) and that the range of the interaction is 2R. As in Section 2.2 we consider the associated empirical fields PΛ n and QΛ n defined by

f (γ)d PΛ n (γ) = 1 λ d (Λ n ) Λ n f (τ u (γ))dP Λ n (γ)du and f (γ)d QΛ n (γ) = 1 λ d (Λ n ) Λ n f (τ u (γ))dQ Λ n (γ)du,
where f is any measurable bounded test function. Following the proof of Proposition 9 we get the existence of an accumulation point P (respectively Q) for ( PΛ n ) (respectively ( QΛ n )). As in Theorem 1, we show that P and Q satisfy the DLR equations and therefore they are both infinite volume Gibbs measures for the Area energy function, the activity z and the inverse temperature β = z. Now it remains to prove that P and Q are different when z is large enough. Note that the difference between P and Q comes only from their boundary conditions "at infinity" (i.e. the boundary conditions of P Λ n and Q Λ n when n goes to infinity).

Let us start with a representation of P Λ n and Q Λ n via the two type Widom-Rowlinson model on Λ n . Consider the following event of allowed configurations on C 2

Λ n A = (γ 1 , γ 2 ) ∈ C 2 Λ n , s.t. a) L R/2 (γ 1 ) ∩ L R/2 (γ 2 ) = / 0 b) L R/2 (γ 1 ) ∩ Λ c n = / 0 ( 33
)
which assumes first that the balls with radii R/2 centred at γ 1 and γ 2 do not overlap and secondly that the balls centred at γ 1 are completely inside Λ n . The two type Widom-Rowlinson model on Λ n with boundary condition b) is the probability measure PΛ n on C 2 Λ n which is absolutely continuous with respect to the product

(π z Λ n ) ⊗2 with density 1 Zn 1 A (γ 1 , γ 2 )z N Λn (γ 1 ) z N Λn (γ 2 ) dπ Λ n (γ 1 )dπ Λ n (γ 2 ),
where ZΛ n is a normalization factor.

Lemma 4. The first marginal (respectively the second marginal) distribution of

PΛ n is P Λ n (respectively Q Λ n ).
Proof. By definition of PΛ n , its first marginal admits the following unnormalized density with respect to π Λ n (dγ 1 )

f (γ 1 ) = 1 A (γ 1 , γ 2 )z N Λn (γ 1 ) z N Λn (γ 2 ) dπ Λ n (γ 2 ) = e (z-1)λ d (Λ n ) z N Λn (γ 1 ) 1 A (γ 1 , γ 2 )dπ z Λ n (γ 2 ) = e (z-1)λ d (Λ n ) z N Λn (γ 1 ) 1 γ 1 Λn\Λ n = / 0 e -zArea Λ n ∩L R (γ 1 )
which is proportional to the density of P Λ n . A similar computation gives the same result for Q Λ n .

Contents

Now let us give a representation of the two type Widom-Rowlinson model via the random cluster model. The random cluster process R Λ n is a point process on Λ n distributed by

1 Ẑn z N Λn (γ) 2 N Λn cc (γ) dπ Λ n (γ),
where N Λ n cc (γ) is the number of connected components of L R/2 (γ) which are completely included in Λ n . Then we build two new point processes Γ 1

Λ n and Γ 2 Λ n by splitting randomly and uniformly the connected component of

R Λ n . Each connected component inside Λ n is given to Γ 1 Λ n or Γ 2
Λ n with probability an half each. The connected components hitting Λ c n are given to Γ 2 Λ n . Rigorously this construction is done by the following way. Let us consider (C i (γ)) 1≤i≤N Λn cc (γ) the collection of connected components of L R/2 (γ) inside Λ n . Let (ε i ) i≥1 be a sequence of independent Bernoulli random variables with parameter 1/2. The processes Γ 1

Λ n and Γ 2 Λ n are defined by

Γ 1 Λ n = 1≤i≤N Λn cc (R Λn ), ε i =1 R Λ n ∩C i (R Λ n ) and Γ 2 Λ n = R Λ n \ Γ 1 Λ n . Lemma 5. The distribution of ( Γ 1 Λ n , Γ 2 Λ n ) is the two-type Widom-Rowlinson model with boundary condition b). In particular, Γ 1 Λ n ∼ P Λ n and Γ 2 Λ n ∼ Q Λ n . Proof.
For any bounded measurable test function f we have

E( f ( Γ 1 Λ n , Γ 2 Λ n )) = E    f    1≤i≤N Λn cc (R Λn ), ε i =1 R Λ n ∩C i (R Λ n ), R Λ n ∩ 1≤i≤N Λn cc (R Λn ), ε i =1 C i (R Λ n ) c       = 1 Ẑn ∑ (ε i )∈{0,1} N Λn cc (γ) 1 2 N Λn cc (γ) f   1≤i≤N Λn cc (γ), ε i =1 γ ∩C i (γ), γ ∩ 1≤i≤N Λn cc (γ), ε i =1 C i (γ) c   z N Λn (γ) 2 N Λn cc (γ) dπ Λ n (γ) = 1 Ẑn ∑ (ε x )∈{0,1} γ (1 A f ) x∈γ, ε x =1 {x}, γ\ x∈γ, ε x =1 {x} z N Λn (γ) dπ Λ n (γ) = 1 Ẑn (1 A f )   (x,ε x )∈ γ, ε x =1 {x}, (x,ε x )∈ γ, ε x =0 {x}   (2z) N Λn (γ) d πΛ n ( γ)
where πΛ n is a marked Poisson point process on Λ n × {0, 1}. It means that the points are distributed by π Λ n and that each point x is marked independently by a Bernoulli variable ε x with parameter 1/2. We obtain

E( f ( Γ 1 Λ n , Γ 2 Λ n )) = e |Λ n | Ẑn (1 A f )   (x,ε x )∈ γ, ε x =1 {x}, (x,ε x )∈ γ, ε x =0 {x}   z N Λn (γ) d π2 Λ n ( γ) = e |Λ n | Ẑn (1 A f ) γ 1 , γ 2 z N Λn (γ 1 ) z N Λn (γ 2 ) dπ Λ n (γ 1 )dπ Λ n (γ 2 ),
which proves the Lemma.

Note that the random cluster process R Λ n is a finite volume GPP with energy function Ĥ = -N Λ n cc , activity z and inverse temperature log [START_REF] Baddeley | Case studies in spatial point process models[END_REF]. Its local energy ĥ is defined by ĥ

(x, γ) = N Λ n cc (γ) -N Λ n cc (γ ∪ {x}).
Thanks to a geometrical argument, it is not difficult to note that ĥ is uniformly bounded from above by a constant c d (depending only on the dimension d). For instance, in the case d = 2, a ball with radius R/2 can overlap at most 5 disjoints balls with radius R/2 and therefore c 2 = 5 -1 = 4 is suitable.

By Lemma 2, we deduce that the distribution of

R Λ n dominates the Poisson point distribution π 2ze -c d Λ n . So we choose z > z d e c d 2R d
which implies that the Boolean model with intensity 2ze -c d and radii R/2 percolates with probability one (see Proposition 13). For any γ ∈ C , we denote by C ∞ (γ) the unbounded connected components in L R/2 (γ) (if it exists) and we define by α the intensity of points in C ∞ (γ) under the distribution π 2ze -c d ;

α := N [0,1] d γ ∩C ∞ (γ) dπ 2ze -c d (γ) > 0. (34)
We are now in position to finish the proof of Theorem 5 by proving that the difference in intensities between Q and P is larger than α.

The local convergence topology τ L ensures that, for any local bounded function f , the evaluation P → f dP is continuous. Actually, the continuity of such evaluation holds for the larger class of functions f satisfying: i) f is local on some bounded set ∆ ii) there exists

A > 0 such that | f (γ)| ≤ A(1 + #(γ)).
In particular, the application P → i(P) := N [0,1] d (γ)P(dγ) is continuous (see [START_REF] Georgii | Large deviations and the maximum entropy principle for marked point random fields[END_REF] for details). We deduce that

i( Q) -i( P) = N [0,1] d (γ)d Q(γ) -N [0,1] d (γ)d P(γ) = lim n→∞ N [0,1] d (γ)d QΛ n (γ) -N [0,1] d (γ)d PΛ n (γ) = lim n→∞ 1 λ d (Λ n ) Λ n N [0,1] d (τ u γ)dQ Λ n (γ) -N [0,1] d (τ u γ)dP Λ n (γ) du.
By the representation of P Λ n and Q λ n given in Lemma 5, we find

i( Q) -i( P) = lim n→∞ 1 λ d (Λ n ) Λ n E N [0,1] d (τ u Γ 2 Λ n ) -N [0,1] d (τ u Γ 1 Λ n ) du = lim n→∞ 1 λ d (Λ n ) Λ n E N τ u [0,1] d (R Λ n ∩C b (R Λ n )) du,
where

C b (γ) are the connected components of L R/2 (γ) hitting Λ c n . Since the distri- bution of R Λ n dominates π 2ze -c d Λ n , i( Q) -i( P) ≥ lim n→∞ 1 λ d (Λ n ) [-n,n-1] d N τ u [0,1] d γ ∩C(γ) ∞ dπ 2ze -c d Λ n (γ)du, ≥ lim n→∞ 1 λ d (Λ n ) [-n,n-1] d αdu = α > 0.
The theorem is proved.

3 Estimation of parameters.

In this section we investigate the parametric estimation of the activity z * and the inverse temperature β * of an infinite volume Gibbs point process Γ . As usual the star specifies that the parameters z * , β * are unknown whereas the variable z and β are used for the optimization procedures. Here the dataset is the observation of Γ trough the bounded window Λ n = [-n, n] d (i.e. the process Γ Λ n ). The asymptotic means that the window Λ n increases to the whole space R d (i.e. n goes to infinity) without changing the realization of Γ . For sake of simplicity, we decide to treat only the case of two parameters (z, β ) but it would be possible to consider energy functions depending on an extra parameter θ ∈ R p . The case where H depends linearly on θ can be treated exactly as z and β . For the non linear case the setting is much more complicated and each procedure has to be adapted. References are given in each section.

In all the section, we assume that the energy function H is stationary and has a finite range R > 0. The existence of Γ is therefore guaranteed by Theorem 1. The procedures presented below are not affected by the uniqueness or non-uniqueness of the distribution of such GPP.

In Section 3.1, we start by presenting the natural maximum likelihood estimator. Afterwards, in Section 3.2, we introduce the general Takacs-Fiksel estimator which is a mean-square procedure based on the GNZ equations. The standard maximum pseudo-likelihood estimator is a particular case of such estimator and is presented in Section 3.3. An application to an unobservable issue is treated in Section 3.4. The last Section 3.5 is devoted to a new estimator based on a variational GNZ equation.

Maximum likelihood estimator

The natural method to estimate the parameters is the likelihood inference. However a practical issue is that the likelihood depends on the intractable partition function. In the case of sparse data, approximations were first proposed in [START_REF] Ogata | Likelihood analysis of spatial point patterns[END_REF], before simulation-based methods have been developed [START_REF] Geyer | Simulation procedures and likelihood inference for spatial point processes[END_REF]. Here, we treat only the theoretical aspects of the MLE and these practical issues are not investigated. Definition 9. The maximum likelihood estimator of (z * , β * ) is given for any n ≥ 1 by

(ẑ n , βn ) = argmax z>0,β ≥0 1 Z z,β Λ n z N Λn (Γ ) e -β H(Γ Λn ) . (35)
Note that the argmax is not necessarily unique and that the boundary effects are not considered in this version of MLE. Other choices could be considered.

In this section we show the consistency of such estimators. The next natural question concerns the asymptotic distribution of the MLE but this problem is more arduous and is still partially unsolved today. Indeed, Mase [START_REF] Mase | Uniform LAN condition of planar Gibbsian point processes and optimality of maximum likelihood estimators of soft-core potential functions[END_REF] and Jensen [START_REF] Jensen | Asymptotic normality of estimates in spatial point processes[END_REF] proved that the MLE is asymptotically normal when the parameters z and β are small enough. Without these conditions, phase transition may occur and some longrange dependence phenomenon can appear. The MLE might then exhibit a non standard asymptotic behavior, in the sense that the rate of convergence might differ from the standard square root of the size of the window and the limiting law might be non-gaussian.

The next theorem is based on a preprint by Mase [START_REF] Mase | Asymptotic properties of MLEs of Gibbs models on rd[END_REF]. See alse [START_REF] Dereudre | Consistency of likelihood estimation for Gibbs point processes[END_REF] for general results on consistency. Theorem 6. We assume that the energy function H is stationary, finite range and not almost surely constant (i.e. there exists a subset Λ ⊂ R d such that H(γ Λ ) is not π Λ (dγ Λ ) almost surely constant). We assume also that the mean energy exists for any stationary probability measure P (i.e. the limit (25) exists) and that the boundary effects assumption [START_REF] Geyer | Simulation procedures and likelihood inference for spatial point processes[END_REF] holds. Moreover we assume that for any ergodic Gibbs measure P, the following limit holds for P-almost every γ

lim n →∞ 1 λ d (Λ n ) H(γ Λ n ) = H(P). (36)
Then, almost surely the parameters (ẑ n , βn ) converge to (z * , β * ) when n goes to infinity.

Proof. Let us assume that the Gibbs distribution P of Γ is ergodic. Ortherwise P can be represented as a mixture of ergodic stationary Gibbs measures (see [START_REF] Preston | Random fields[END_REF], Theorem 2.2 and 4.1). Therefore the proof of the consistency of the MLE reduces to the case when P is ergodic, which is assumed henceforth.

Let us consider the log-likelihood contrast function

K n (θ , β ) = -log(Z e -θ ,β Λ n ) -θ N Λ n (Γ ) -β H(Γ Λ n )
related to the parametrization θ =log(z). It is clear that (ẑ n , βn ) = (e -θn , βn ) where ( θn , βn ) is the argmax of (θ , β ) → K n (θ , β ). So it is sufficient to show that ( θn , βn ) converges almost surely to (log(z * ), β * ). The limit [START_REF] Georgii | Conditional intensity and Gibbsianness of determinantal point processes[END_REF], the ergodic Theorem and the assumption [START_REF] Mase | Asymptotic properties of MLEs of Gibbs models on rd[END_REF] imply the existence of the following limit contrast function

K(θ , β ) := -p e -θ ,β -θ E P (N [0,1] d (Γ )) -β H(P) = lim n→∞ K n (θ , β ) λ d (Λ n ) .
The variational principle (Theorem 3) ensures that (θ , β ) → K(θ , β ) is lower than I 1 (P) with equality if and only if P is a Gibbs measure with energy function H, activity z and inverse temperature β . Since H is not almost surely constant, it is easy to see that two Gibbs measures with different parameters z, β are different (this fact can be viewed used the DLR equations in a very large box Λ ). Therefore K(θ , β ) is maximal, equal to I 1 (P), if and only if (θ , β ) = (θ * , β * ).

Therefore it remains to prove that the maximizers of (θ , β ) → K n (θ , β ) converge to the unique maximizer of (θ , β ) → K(θ , β ). First note that the functions K n are concave. Indeed, the Hessian of K n is negative since

∂ 2 K n (θ , β ) ∂ 2 θ = -Var P e -θ ,β Λn (N Λ n ), ∂ 2 K n (θ , β ) ∂ 2 β = -Var P e -θ ,β Λn (H) and ∂ 2 K n (θ , β ) ∂ θ ∂ β = -Cov P e -θ ,β Λn (N Λ n , H).
The convergence result for the argmax follows since the function

(θ , β ) → K(θ , β ) is necessarily strictly concave at (θ * , β * ) because K(θ , β ) is maximal uniquely at (θ * , β * ).
Let us finish this section with a discussion on the extra assumption [START_REF] Mase | Asymptotic properties of MLEs of Gibbs models on rd[END_REF] which claims that the empirical mean energy converges to the expected value energy. This assumption is in general proved via the ergodic theorem or a law of large numbers. In the case of the Area energy function H defined in (4), it is a direct consequence of a decomposition as in [START_REF] Hall | On continuum percolation[END_REF] and the ergodic Theorem. In the case of pairwise interaction, the verification follows essentially the proof of Proposition 12.

Takacs-Fiksel estimator

In this section we present an estimator introduced in the eighties by Takacs and Fiksel [START_REF] Fiksel | Estimation of parametrized pair potentials of marked and nonmarked Gibbsian point processes[END_REF][START_REF] Takacs | Estimator for the pair-potential of a Gibbsian point process[END_REF]. It is based on the GNZ equations presented in Section 2.5. Let us start by explaining briefly the procedure. Let f be a test function from R d × C to R. We define the following quantity for any z > 0, β > 0 and γ ∈ C

C z,β Λ n ( f , γ) = ∑ x∈γ Λn f (x, γ\{x}) -z Λ n e -β h(x,γ) f (x, γ)dx. ( 37
)
By the GNZ equation ( 22) we obtain

E C z * ,β *
Λ n ( f ,Γ ) = 0 where Γ is a GPP with parameter z * and β * . Thanks to the ergodic Theorem it follows that for n large enough

C z * ,β * Λ n ( f ,Γ ) λ d (Λ n ) ≈ 0.
Then the Takacs-Fiksel estimator is defined as a mean-square method based on functions C z * ,β * Λ n ( f k ,Γ ) for a collection of test functions ( f k ) 1≤k≤K . Definition 10. Let K ≥ 2 be an integer and ( f k ) 1≤k≤K a family of K functions from R d × C to R. The Takacs-Fiksel estimator (ẑ n , βn ) of (z * , β * ) is defined by

(ẑ n , βn ) = argmin (z,β )∈D K ∑ k=1 C z,β Λ n ( f k ,Γ ) 2 ,
where D ⊂ (0, +∞) × [0, +∞) is a bounded domain containing (z * , β * ).

In opposition to the MLE procedure, the contrast function does not depend on the partition function. This estimator is explicit except for the computation of integrals and the optimization procedure. In [START_REF] Coeurjolly | Takacs-Fiksel method for stationary marked Gibbs point processes[END_REF] the Takacs-Fiksel procedure is presented in a more general setting including the case where the functions f k depend on parameters z and β . This generalization may lead to a simpler procedure in choosing f k such that the integral term in (37) is explicitly computable.

In the rest of the section, we prove the consistency of the estimator. General results on consistency and asymptotic normality are developed in [START_REF] Coeurjolly | Takacs-Fiksel method for stationary marked Gibbs point processes[END_REF]. Theorem 7 (Consistency). We make the following integrability assumption: for any

1 ≤ k ≤ K E | f k (0,Γ )|(1 + |h(0,Γ )|) sup (z,β )∈D e -β h(0,Γ ) < +∞. (38)
We assume also the following identifiability condition: the equality

K ∑ k=1 E f k (0,Γ ) ze -β h(0,Γ ) -z * e -β * h(0,Γ ) 2 = 0 ( 39
)
holds if and only (z, β ) = (z * , β * ). Then the Takacs-Fiksel estimator (ẑ n , βn ) presented in Definition 10 converges almost surely to (z * , β * ).

Proof. As in the proof of Theorem 6, without loss of generality, we assume that the Gibbs distribution of Γ is ergodic. Therefore, thanks to the ergodic Theorem, almost surely for any 1

≤ k ≤ K lim n →∞ C z,β Λ n ( f k ,Γ ) λ d (Λ n ) = E   ∑ x∈Γ [0,1] d f k (x,Γ \x)   -zE [0,1] d e -β h(x,Γ ) f k (x,Γ )dx . ( 40
)
By the GNZ equation ( 22)

E   ∑ x∈Γ [0,1] d f k (x,Γ \x)   = z * E [0,1] d e -β * h(x,Γ ) f k (x,Γ )dx . (41)
Using the stationarity and compiling ( 40) and ( 41), we obtain that the contrast function

K n (z, β ) = K ∑ k=1 C z,β Λ n ( f k ,Γ ) λ d (Λ n ) 2 admits almost surely the limit lim n →∞ K n (z, β ) = K(z, β ) := K ∑ k=1 E f k (0,Γ ) ze -β h(0,Γ ) -z * e -β * h(0,Γ ) 2 ,
which is null if and only if (z, β ) = (z * , β * ). Therefore it remains to prove that the minimizers of the contrast function converge to the minimizer of the limit contrast function. In the previous section we solved a similar issue for the MLE procedure using the convexity of contrast functions. This argument does not work here and we need more sophisticated tools. We define by W n (.) the modulus of continuity of the contrast function K n ; let η be a positive real

W n (η) = sup |K n (z, β )-K n (z , β )|, with (z, β ), (z , β ) ∈ D, (z-z , β -β ) ≤ η .
Lemma 6 (Theorem 3.4.3 [START_REF] Guyon | Random fields on a network[END_REF]). Assuming that there exists a sequence (ε l ) l≥1 , which goes to zero when l goes to infinity, such that for any l ≥ 1

P lim sup n →+∞ W n 1 l ≥ ε l = 0 ( 42
)
then almost surely the minimizers of (z, β ) → K n (z, β ) converges to the minimizer of (z, β ) → K(z, β ).

Let us show that the assertion (42) holds. Thanks to equalities [START_REF] Molchanov | Consistent estimation of the parameters of Boolean models of random closed sets[END_REF], [START_REF] Møller | Lectures on random Voronoȋ tessellations[END_REF] and assumption [START_REF] Mayer | Molecular distributions[END_REF], there exists a constant C 1 such that for n large enough, any 1 ≤ k ≤ K and any (z,

β ) ∈ D |C z,β Λ n ( f k ,Γ )| λ d (Λ n ) ≤ C 1 . (43)
We deduce that for n large enough

|K n (z, β ) -K n (z , β )| ≤ C 1 λ d (Λ n ) K ∑ k=1 Λ n | f k (x,Γ )| ze -β h(x,Γ ) -z e -β h(x,Γ ) dx ≤ C 1 |β -β | λ d (Λ n ) max 1≤k≤K Λ n | f k (x,Γ )h(x,Γ )| sup (z,β )∈D ze -β h(x,Γ ) dx + C 1 |z -z | λ d (Λ n ) max 1≤k≤K Λ n | f k (x,Γ )| sup (z,β )∈D e -β h(x,Γ ) dx.
By the ergodic Theorem, the following convergences hold almost surely

lim n →+∞ 1 λ d (Λ n ) Λ n | f k (x,Γ )h(x,Γ )| sup (z,β )∈D ze -β h(x,Γ ) dx = E | f k (0,Γ )h(0,Γ )| sup (z,β )∈D ze -β h(0,Γ ) < +∞, and lim n →+∞ 1 λ d (Λ n ) Λ n | f k (x,Γ )| sup (z,β )∈D e -β h(x,Γ ) dx = E | f k (0,Γ )| sup (z,β )∈D e -β h(0,Γ ) < +∞.
This implies the existence of a constant C 2 > 0 such that for n large enough, any 1 ≤ k ≤ K and any (z,

β ) ∈ D |K n (z, β ) -K n (z , β )| < C 2 (z -z , β -β ) .
The assumption [START_REF] Møller | Likelihood inference for unions of interacting discs[END_REF] occurs with the sequence ε l = C 2 /l and Theorem 7 is proved.

Remark 1 (On the integrability assumption). The integrability assumption [START_REF] Mayer | Molecular distributions[END_REF] is sometimes difficult to check, especially when the local energy h(0, γ) is not bounded from below. For instance in the setting of pairwise energy function H defined in (1) with a pair potential ϕ having negative values, Ruelle estimates ( 16) are very useful. Indeed, by stability of the energy function, the potential ϕ is necessary bounded from below by 2A and therefore

E e -β h(0,Γ ) < E e -2Aβ N B(0,R) (Γ )) < +∞,
where R is the range of the interaction.

Remark 2 (On the identifiability assumption). In the identifiability assumption (39), the sum is null if and only if each term is null. Assuming that the functions are regular enough, each term is null as soon as (z, β ) belongs to a 1-dimensional manifold embedded in R 2 containing (z * , β * ). Therefore, assumption [START_REF] Meester | Continuum percolation, volume 119 of Cambridge Tracts in Mathematics[END_REF] claims that (z * , β * ) is the unique element of these K manifolds. If K ≤ 2, there is no special geometric argument to ensure that K 1-dimensional manifolds in R 2 have an unique intersection point. For this reason, it is recommended to choose K ≥ 3. See Section 5 in [START_REF] Coeurjolly | Takacs-Fiksel method for stationary marked Gibbs point processes[END_REF] for more details and complements on this identifiability assumption.

Maximum pseudo-likelihood estimator

In this section we present the maximum pseudo-likelihood estimator, which is a particular case of the Takacs-Fiksel estimator. This procedure has been first introduced by Besag in [START_REF] Besag | Spatial interaction and the statistical analysis of lattice systems[END_REF] and popularized by Jensen and Moller in [START_REF] Jensen | Pseudolikelihood for exponential family models of spatial point processes[END_REF] and Baddeley and Turner in [START_REF] Baddeley | Practical maximum pseudolikelihood for spatial point patterns (with discussion)[END_REF]. Definition 11. The maximum pseudo-likelihood estimator (ẑ n , βn ) is defined as a Takacs-Fiksel estimator (see Definition 10)

with K = 2, f 1 (x, γ) = 1 and f 2 (x, γ) = h(x, γ).
This particular choice of functions f 1 , f 2 simplifies the identifiability assumption [START_REF] Meester | Continuum percolation, volume 119 of Cambridge Tracts in Mathematics[END_REF]. The following theorem is an adaptation of Theorem 7 in the present setting of MPLE. The asymptotic normality is investigated first in [START_REF] Jensen | On asymptotic normality of pseudo likelihood estimates for pairwise interaction processes[END_REF] (see also [6] for more general results).

Theorem 8 (Consistency). Assuming

E (1 + h(0,Γ ) 2 ) sup (z,β )∈D e -β h(0,Γ ) < +∞ (44)
and P h(0,Γ ) = h(0, / 0) < 1,

then the maximum pseudo-likelihood estimator (ẑ n , βn ) converges almost surely to (z * , β * ).

Proof. Let us check the assumptions of Theorem 7. Clearly, the integrability assumption (44) ensures the integrability assumptions [START_REF] Mayer | Molecular distributions[END_REF] with f 1 = 1 and f 2 = h. So it remains to show that assumption [START_REF] Ogata | Likelihood analysis of spatial point patterns[END_REF] implies the identifiability assumption [START_REF] Meester | Continuum percolation, volume 119 of Cambridge Tracts in Mathematics[END_REF]. Consider the parametrization z = e -θ and ψ the function

ψ(θ , β ) = E e -θ * -β * h(0,Γ ) (e U -U -1) , with U = β * h(0,Γ ) + θ * -β h(0,Γ ) -θ .
The function ψ is convex, non negative and equal to zero if and only if U is almost surely equal to zero. By assumption [START_REF] Ogata | Likelihood analysis of spatial point patterns[END_REF] this fact occurs when (z, β ) = (z * , β * ).

Therefore the gradient ∇ψ = 0 if and only (z, β ) = (z * , β * ). Noting that

∂ ψ(θ , β ) ∂ θ = E z * e -β * h(0,Γ ) -ze -β h(0,Γ ) and ∂ ψ(θ , β ) ∂ β = E h(0,Γ ) z * e -β * h(0,Γ ) -ze -β h(0,Γ ) ,
the identification assumption (39) holds. The theorem is proved.

Solving an unobservable issue

In this section we give an application of the Takacs-Fiksel procedure in a setting of partially observable dataset. Let us consider a Gibbs point process Γ for which we observe only L R (Γ ) in place of Γ . This setting appears when Gibbs point processes are used for producing random surfaces via germ-grain structures (see [START_REF] Møller | Likelihood inference for unions of interacting discs[END_REF] for instance). Applications for modelling micro-structure in materials or micro-emulsion in statistical physics are developed in [START_REF] Chiu | Stochastic geometry and its applications[END_REF].

The goal is to furnish an estimator of z * and β * in spite of this unobservable issue. Note that the number of points (or balls) is not observable from L R (Γ ) and therefore the MLE procedure is not achievable, since the likelihood is not computable. When β is known and fixed to zero, it corresponds to the estimation of the intensity of the Boolean model from its germ-grain structure (see [START_REF] Molchanov | Consistent estimation of the parameters of Boolean models of random closed sets[END_REF] for instance).

In the following we assume that Γ a Gibbs point process for the Area energy function defined in (4), the activity z * and the inverse temperature β * . This choice is natural since the energy function depends on the observations L R (Γ ). The more gen-eral setting of Quermass interaction is presented in [START_REF] Dereudre | Estimation of the intensity parameter of the germ-grain quermass-interaction model when the number of germs is not observed[END_REF] but for sake of simplicity, we treat only here the simpler case of Area interaction.

We opt for a Takacs-Fiksel estimator but the main problem is that the function

C z,β Λ n ( f , γ) = ∑ x∈γ Λn f (x, γ\{x}) -z
Λ n e -β h(x,γ) f (x, γ)dx, which appears in the procedure, is not computable since the positions of points are not observable. The main idea is to choose the function f properly such that the sum is observable although each term of the sum is not. To this end, we define

f 1 (x, γ) = Surface ∂ B(x, R) ∩ L c R (γ)
and f 2 (x, γ) = 1 {B(x,R)∩L R (γ)= / 0} , where ∂ B(x, R) is the boundary of the ball B(x, R) (i.e. the sphere S(x, R)) and the "Surface" means the (d -1)-dimensional Hausdorff measure in R d . Clearly the function f 1 gives the surface of the portion of the sphere S(x, R) outside the germgrain structure L R (γ). The function f 2 indicates if the ball B(x, R) hits the germ-grain structure L R (γ). Therefore we obtain that

∑ x∈γ Λn f 1 (x, γ\{x}) = Surface ∂ L R (γ Λ n ) and ∑ x∈γ Λn f 2 (x, γ\{x}) = N iso L R (γ Λ n ) ,
where N iso (L R (γ Λ n )) is the number of isolated balls in the germ-grain structure L R (γ Λ n ). Let us note that these quantities are not exactly observable since, in practice, we observe L R (γ) ∩ Λ n rather than L R (γ Λ n ). However, if we omit this boundary effect, the values C z,β Λ n ( f 1 ,Γ ) and C z,β Λ n ( f 2 ,Γ ) are observable and the Takacs-Fiksel procedure is achievable. The consistency of the estimator is guaranteed by Theorem 7. The integrability assumption [START_REF] Mayer | Molecular distributions[END_REF] is trivially satisfied since the functions f 1 , f 2 and h are uniformly bounded. The verification of the identifiability assumption [START_REF] Meester | Continuum percolation, volume 119 of Cambridge Tracts in Mathematics[END_REF] is more delicate and we refer to [START_REF] Coeurjolly | Takacs-Fiksel method for stationary marked Gibbs point processes[END_REF], example 2 for a proof. Numerical estimations on simulated and real datasets can be found in [START_REF] Dereudre | Estimation of the intensity parameter of the germ-grain quermass-interaction model when the number of germs is not observed[END_REF].

A variational estimator

In this last section, we present a new estimator based on a variational GNZ equation which is a mix between the standard GNZ equation and an integration by parts formula. This equation has been first introduced in [START_REF] Dereudre | Diffusion infini-dimensionnelles et champs de Gibbs sur l'espace des trajectoires continues[END_REF] for statistical mechanics issues and used recently in [START_REF] Baddeley | Variational estimators for the parameters of Gibbs point process models[END_REF] for spatial statistic considerations. In the following, we present first this variational equation and afterwards we introduce its associated estimator of β * . The estimation of z * is not considered here. Theorem 9. Let Γ be a GPP for the energy function H, the activity z and the inverse temperature β . We assume that, for any γ ∈ C , the function x → h(x, γ) is differentiable on R d \γ. Let f be a function from R d × C to R which is differentiable and with compact support with respect to the first variable. Moreover we assume the integrability of both terms below. Then

E ∑ x∈Γ ∇ x f (x,Γ \{x}) = β E ∑ x∈Γ f (x,Γ \{x})∇ x h(x,Γ \{x}) . ( 46
)
Proof. By the standard GNZ equation ( 22) applied to the function ∇ x f , we obtain

E ∑ x∈Γ ∇ x f (x,Γ \{x}) = zE R d
e -β h(x,Γ ) ∇ x f (x,Γ )dx .

By a standard integration by part formula with respect to the first variable x, we find that

E ∑ x∈Γ ∇ x f (x,Γ \{x}) = zβ E R d ∇ x h(x,Γ )e -β h(x,Γ ) f (x,Γ )dx .
Using again the GNZ equation we finally obtain [START_REF] Poghosyan | Abstract cluster expansion with applications to statistical mechanical systems[END_REF]. Note that equation ( 46) is a vectorial equation. For convenience it is possible to obtain a real equation by summing each coordinate of the vectorial equation. The gradient operator is simply replaced by the divergence operator.

Remark 3 (on the activity parameter z). The parameter z does not appear in the variational GNZ equation [START_REF] Poghosyan | Abstract cluster expansion with applications to statistical mechanical systems[END_REF]. Therefore these equations do not characterize the Gibbs measures as in Proposition 5. Actually these variational GNZ equations characterize the mixing of Gibbs measures with random activity (See [START_REF] Dereudre | Diffusion infini-dimensionnelles et champs de Gibbs sur l'espace des trajectoires continues[END_REF] for details).

Let us now explain how to estimate β * from these variational equations. When the observation window Λ n is large enough we identify the expectations of sums in [START_REF] Poghosyan | Abstract cluster expansion with applications to statistical mechanical systems[END_REF] by the sums. Then the estimator of β * is simply defined by βn = ∑ x∈Γ Λn div x f (x,Γ \{x})

∑ x∈Γ Λn f (x,Γ \{x})div x h(x,Γ \{x}) . ( 47
)
Note that this estimator is very simple and quick to compute in comparison to the MLE, MPLE or the general Takacs-Fiksel estimators. Indeed, in [START_REF] Preston | Random fields[END_REF], there are only elementary operations (no optimization procedure, no integral to compute).

Let us now finish this section with a consistency result. More general results for consistency, asymptotic normality and practical estimations are available in [START_REF] Baddeley | Variational estimators for the parameters of Gibbs point process models[END_REF]. Contents Theorem 10. Let Γ be a GPP for a stationary and finite range energy function H, activity z * and inverse temperature β * . We assume that, for any γ ∈ C , the function x → h(x, γ) is differentiable on R d \γ. Let f be a stationary function from R d × C to R, differentiable with respect to the first variable and such that E (| f (0,Γ | + |∇ x f (0,Γ )| + | f (0,Γ )∇ x h(0,Γ )|)e -β * h(0,Γ ) < +∞ [START_REF] Ruelle | Statistical mechanics: Rigorous results[END_REF] and E f (0,Γ )div x h(0,Γ )e -β * h(0,Γ ) = 0. (

) 49
Then the estimator βn converges almost surely to β * .

Proof. As usual, without loss of generality, we assume that the Gibbs distribution of Γ is ergodic. Then by the ergodic theorem the following limits both hold almost surely , where the denominator is not null thanks to assumption [START_REF] Ruelle | Superstable interactions in classical statistical mechanics[END_REF]. Therefore it remains to prove the following variational GNZ equation

E   ∑ x∈Γ [0,1] d ∇ x f (x,Γ \{x})   = β * E   ∑ x∈Γ [0,1] d f (x,Γ \{x})∇ x h(x,Γ \{x})   . (52)
Note that this equation is not a direct consequence of the variational GNZ equation [START_REF] Poghosyan | Abstract cluster expansion with applications to statistical mechanical systems[END_REF] since the function x → f (x, γ) does not have a compact support. We need the following cut-off approximation. Let us consider (ψ n ) n≥1 any sequence of functions from R d to R such that ψ n is differentiable, equal to 1 on Λ n , 0 on Λ c n+1 and such that |∇ψ n | and |ψ n | are uniformly bounded by a constant C (which does not depend on n). It is not difficult to build such a sequence of functions. Let us now apply the variational GNZ equation [START_REF] Poghosyan | Abstract cluster expansion with applications to statistical mechanical systems[END_REF] to the function (x, γ) → ψ n (x) f (x, γ), we obtain

E ∑ x∈Γ ψ n (x)∇ x f (x,Γ \{x}) + E ∑ x∈Γ ∇ x ψ n (x) f (x,Γ \{x}) = β * E ∑ x∈Γ ψ n (x) f (x,Γ \{x})∇ x h(x,Γ \{x}) . (53)

ζΛ

  the distribution of a Poisson point process with intensity ζ λ d Λ which is also the distribution of a stationary Poisson point process with intensity ζ restricted to Λ . For sake of brevity, π and π Λ denote the distribution of Poisson point processes with intensity ζ = 1.

Proposition 2 .

 2 The function z → E P z,β Λ (N Λ ) is continuous and differentiable, with derivative z → Var P z,β Λ (N Λ )/z on (0, +∞). Similarly the function β → E P z,β Λ (H) is continuous and differentiable with derivative β → -Var P z,β Λ (H) on R + .

Definition 7 .

 7 The energy function H has a finite range R > 0 if for every bounded ∆ , the local energy H ∆ (see Definition 3) is a local function on ∆ ⊕ B(0, R). It means that for any finite configuration γ

lim n →+∞ 1 λ 1 λ|

 11 d (Λ n ) ∑ x∈Γ Λn div x f (x,Γ \{x}) = E   ∑ x∈Γ [0,1] d div x f (x,Γ \{x}) d (Λ n ) ∑ x∈Γ Λn f (x,Γ \{x})div x h(x,Γ \{x}) ,Γ \{x})div x h(x,Γ \{x}) expectations in[START_REF] Takacs | Estimator for the pair-potential of a Gibbsian point process[END_REF] and (51) are finite since by the GNZ equations, the stationarity and assumption (48)E   ∑ x∈Γ [0,1] d |div f (x,Γ \{x})|   = E |div f (0,Γ )|e -β * h(0,Γ ) < +∞ and f (x,Γ \{x})divh(x,Γ \{x})|   = E | f (0,Γ )divh(0,Γ )|e -β * h(0,Γ ) < +∞.We deduce that almost surelylim n →+∞ βn = E ∑ x∈Γ [0,1] d div f (x,Γ \{x}) E ∑ x∈Γ [0,1] d f (x,Γ \{x})divh(x,Γ \{x})

≤≤

  Thanks to the GNZ equation and the stationarity we getE ∑ x∈Γ ψ n (x)∇ x f (x,Γ \{x})λ d (Λ n )E   ∑ x∈Γ [0,1] d ∇ x f (x,Γ \{x})   Cz * λ d (Λ n+1 \Λ n )E |∇ x f (0,Γ )|e -β * h(0,Γ ) , and E ∑ x∈Γ ψ n (x) f (x,Γ \{x})∇ x h(x,Γ \{x}) -λ d (Λ n )E   ∑ x∈Γ [0,1] d f (x,Γ \{x})∇ x h(x,Γ \{x})  Cz * λ d (Λ n+1 \Λ n )E | f (0,Γ )∇ x h(0,Γ )|e -β * h(0,Γ ) ,

  So for any P ∈ P Λ such that E P (H) < +∞ and I(P|π Λ ) < +∞ let us show that β E P (H)log(z)E P (N Λ ) + I(P|π Λ ) ≥log(Z

	This equality implies that the minimum of β E P (H) -log(z)E P (N Λ ) + I(P|π Λ )
	should be equal to -log(Z Λ ). z,β z,β Λ ) with equality if and only if P = P z,β Λ . Let f be the density of P with respect to π Λ .
	log(Z	z,β

Acknowledgement: The author thanks P. Houdebert, A. Zass and the anonymous referees for the careful reading and the interesting comments. This work was supported in part by the Labex CEMPI (ANR-11-LABX-0007-01), the CNRS GdR 3477 GeoSto and the ANR project PPP (ANR-16-CE40-0016).

Contents

and finally

Therefore, dividing equation ( 53) by λ d (Λ n ), using the previous approximations and letting n go to infinity, we find exactly the variational equation ( 52). The theorem is proved.