This paper deals with stationary Gibbsian point processes on the plane with an interaction that depends on the tiles of the Delaunay triangulation of points via a bounded triangle potential. It is shown that the class of these Gibbs processes includes all minimisers of the associated free energy density and is therefore nonempty. Conversely, each such Gibbs process minimises the free energy density, provided the potential satisfies a weak long-range assumption.
"Variational Characterisation of Gibbs Measures with Delaunay Triangle Interaction." Electron. J. Probab. 14 2438 - 2462, 2009. https://doi.org/10.1214/EJP.v14-713