The retention-probability conjecture for Gibbs point processes

Let S{\mathcal S} be the state space, let λ(u;)\lambda(u;\emptyset) denote the conditional intensity at uu given the empty configuration, and let α\alpha be a finite upper bound for the conditional intensity. Let pλ,α(u)p_{\lambda,\alpha}(u) be the retention probability appearing in the thinning representation, and let GG and V{\mathcal V} denote the generating functional and void probability, respectively.

Retention-probability conjecture. The function pλ,α()p_{\lambda,\alpha}(\cdot) in the thinning representation is

pλ,α(u)=λ(u;)α.p_{\lambda,\alpha}(u)=\frac{\lambda(u;\emptyset)}{\alpha}.

Consequently,

G(f)=G(1g)=exp{Sg(u)λ(u;)du}G(f)=G(1-g)=\exp\left\{-\int_{\mathcal S}g(u)\lambda(u;\emptyset)\,\mathrm{d}u\right\}

and

V(B)=exp{Bλ(u;)du}.{\mathcal V}(B)=\exp\left\{-\int_B\lambda(u;\emptyset)\,\mathrm{d}u\right\}.

The conjecture would give an explicit retention probability and therefore explicit forms for the generating functional and void probability. It is motivated by the expected dependence of the retention probability on the conditional intensity and by the invariance of the representation under the choice of any admissible upper bound α\alpha; the supplied text does not establish the formula.

Progress summary

Solved

An unverified posted counterexample claims the conjecture is false even for a simple hard-core process, while no independent mathematical confirmation has appeared.

Cronie formulated the conjecture in 2024: the thinning retention probability should equal the empty-configuration conditional intensity divided by the stability bound, yielding explicit generating and void probabilities. The paper labels this as unproved and reports only that ongoing work suggests it may hold.

Known results

  • Last and Otto (2023): a thinning representation with an unspecified retention probability.
  • Stucki and Schuhmacher (2014): upper and lower bounds for the generating functional of locally stable Gibbs point processes.
  • Cronie (2024): formulas involving the unknown function pλ,αp_{\lambda,\alpha}, but no proof of the conjectured identity.

Posted attempt

A reader-written argument claims a counterexample on S=[0,1]S=[0,1] with R>1R>1: the normalized hard-core process has void probability 1/21/2, whereas the conjecture predicts e1e^{-1}. It also identifies an alleged invalid factorization of dependent retention marks. This claimed counterexample has not been independently verified.

Current status (as of August 2026): The conjecture is unproved in the published source, while a posted counterexample claims it is false; neither the refutation nor the original formula has independent verification.

Sources
Sources & referencesView supporting material

Primary source

Ottmar Cronie, “New density/likelihood representations for Gibbs models based on generating functionals of point processes”, arXiv:2406.07075 (2024).

Solutions 1

Counterexample

The conjecture is false even for the source's own hard-core Gibbs process. Its asserted generating-functional theorem is also false; in fact, the proposed identity holds exactly for Poisson processes.

Take

S=[0,1],β=α=1,R>1.S=[0,1],\qquad \beta=\alpha=1,\qquad R>1.

The hard-core Papangelou conditional intensity is

λ(u;ξ)=1{ξ=}.\lambda(u;\xi) = \mathbf1_{\{\xi=\varnothing\}}.

This satisfies the stated local-stability hypothesis

λ(u;ξ)α=1,λ(u;)=1.\lambda(u;\xi)\le\alpha=1, \qquad \lambda(u;\varnothing)=1.

Since RR exceeds the diameter of SS, admissible configurations contain at most one point. Their correctly normalized Janossy densities are

j0=12,j1(u)=121[0,1](u),jr=0(r2).j_0=\frac12,\qquad j_1(u)=\frac12\mathbf1_{[0,1]}(u),\qquad j_r=0\quad(r\ge2).

Therefore the actual void probability and probability generating functional are

V(S)=12,G(f)=12(1+01f(u)du).\mathcal V(S)=\frac12,\qquad G(f)=\frac12\left(1+\int_0^1f(u)\,du\right).

The conjecture instead predicts

V(S)=exp(01λ(u;)du)=e1,\mathcal V(S) = \exp\left(-\int_0^1\lambda(u;\varnothing)\,du\right) =e^{-1},

and

G(f)=exp(01(1f(u))du).G(f)= \exp\left(-\int_0^1(1-f(u))\,du\right).

Hence already f0f\equiv0 gives

12e1.\boxed{\frac12\ne e^{-1}.}

More generally, every activity β>0\beta>0 gives true void probability 1/(1+β)1/(1+\beta), while the conjecture predicts eβe^{-\beta}.

The source's own displayed hard-core Janossy formula also fails normalization: for these parameters it assigns density e1e^{-1} to both the empty configuration and the singleton stratum, giving total mass

e1+01e1du=2e1.e^{-1}+\int_0^1e^{-1}\,du=\frac2e\ne1.

The obstruction persists even with strictly positive densities and finite interaction energies. For 0<γ<10<\gamma<1, define

Zγ=r0γ(r2)r!,jr(u1,,ur)=γ(r2)Zγ.Z_\gamma=\sum_{r\ge0}\frac{\gamma^{\binom r2}}{r!}, \qquad j_r(u_1,\ldots,u_r) = \frac{\gamma^{\binom r2}}{Z_\gamma}.

Then

λ(u;ξ)=γξ1,λ(u;)=1,\lambda(u;\xi)=\gamma^{|\xi|}\le1, \qquad \lambda(u;\varnothing)=1,

but

Zγ<e,Vγ(S)=Zγ1>e1.Z_\gamma<e,\qquad \mathcal V_\gamma(S)=Z_\gamma^{-1}>e^{-1}.

Thus a continuum of full-support locally stable Gibbs processes also contradicts the conjecture.

The source's preceding claimed theorem fails because it replaces

E[xY(1M(x)g(x))|Y]\mathbb E\left[ \prod_{x\in Y}(1-M(x)g(x)) \,\middle|\,Y \right]

by

xY(1E[M(x)Y]g(x))\prod_{x\in Y} \left(1-\mathbb E[M(x)\mid Y]g(x)\right)

for dependent retention marks. Condition on Y={x,y}Y=\{x,y\}, retain exactly one point with equal probabilities, and take g(x)=g(y)=1g(x)=g(y)=1. The actual expression is 00, whereas the incorrectly factorized expression is 1/41/4.

Finally, the sharp corrected characterization is

G(f)=exp[S(1f(u))λ(u;)du] for every fX is Poisson with intensity λ(u;).G(f)= \exp\left[-\int_S(1-f(u))\lambda(u;\varnothing)\,du\right] \text{ for every }f \quad\Longleftrightarrow\quad X\text{ is Poisson with intensity }\lambda(u;\varnothing).

Indeed, the right side is exactly the Poisson generating functional and uniquely determines the point-process law.

Source: arXiv:2406.07075, Conjecture 1, Theorem 1, and the displayed hard-core example.

0 endorsements
Shivam Patel ·