There is a simple connection between ground states of antiferromagnetic Potts models and colorings of the plane: if the unit distance graph of the plane ($G=(\mathbb R^2,\{\{x,y\},d_2(x,y)=1\})$) is admissibly colored with $q$ colors, then this is a ground state of the $q$-state Potts model, and conversely, a ground state without frustration (therefore with trivially computed energy and ground state entropy) yields a coloring of $G$.
But the recent advances in our understanding of the ferromagnetic Potts models let us hope that a related statistical physics connection could yield insight into the chromatic number of the plane. Duminil-Copin et al. have proved that the phase transition for the Potts model is continuous up to and including $q=4$ states. And it is discontinuous/1st order for $q\ge 26$ and conjectured to be for 5 and more states. These results are for the square lattice but they have obtained related results on isoradial subgraphs of the plane, many such graphs are subgraphs of $G$. It seems on page 5 they conjecture that for all planar graphs the nature of the transition is the same given a number of states $q$. It seems to me that this should be interpreted as a hint that the chromatic number of the plane is 5. I think that Erdos and other experts never really knew what to believe about that number, which is 4, 5, 6, or 7. The usual example for a non-3-colorable subgraph is Moser's spindle. While a hexagonal tiling with different colors on adjacent tiles give an admissible 7-coloring.
The rationale for phase transitions at $q>4$ yielding admissible colorings could be, I am guessing, that if $G$ is $q$-colorable it is possible with $q$ states to have an ordered phase, with an infinite cluster (thus not admissible), at $T=T_t$ (say adiabatically continued raising temperature from $T<T_t$) coexist with a lower energy disordered phase, i.e. with small finite clusters (i.e. avoiding distance 1, of diameter<1), at the same temperature. This is a characteristic of discontinuous transitions (like water/gas below critical point) that there are 2 phases with different energies at the same temperature at the transition $T=T_t$.
So what do probabilists and statistical physicists think of this possible link? When you prove discontinuous phase transition what is the measure of the different phases at transition temperature? Can you have systems with big measure for only one side (ordered or disordered) of the transition?
Thanks a lot for any comments.
EDIT: I fixed a couple of my statements and added references. A quick description of the Potts model is that it is a generalization of the Ising model where the "spins" can take more than 2 values on the unit circle and the repulsion energy is proportional to the scalar product of those times a constant coefficient $J$, positive=ferromagnetic (adjacent spins tend to oppose), negative=antiferromagnetic. The energy can be simplified to just a delta function between spins instead of scalar product. The partition function at temp $1/\beta$ is $Z(\beta)=\sum_\sigma e^{-\beta\sum H(\sigma_{ij})}$ with $H(\sigma_{ij})$ the energy contribution from the $ij$ sites interaction in state $\sigma$.
2ND EDIT: The rationale for considering Potts models is trying to find a "local" interaction which would create proper colorings of $G$ by favoring adjacent spins to be aligned but far enough spins to be unaligned with high enough probability. The Potts model at high temperature favors some random flipping at mid distance while keeping adjacent spins clustering on the same color. So the clusters should not be too large and yield proper colorings "usually". Another justification I would bring is that proper colorings should have a complicated structure which would look random and that the only way of yielding that seems to be the "probabilistic method" and apart from local laws for constraining the desired properties on coloring I do not see what would. And local interactions should have all been considered by physicists, and finally because of the universality of behaviors for many classes of interactions looking at Potts models with various parameters should be hopefully enough to find the right universality class to build at least on some "visible" part of phase space (i.e. with positive propability, or for typical states in simulations with Boltzmann/Gibbs measure) proper colorings, at least on some part of the plane and then we can tile the whole plane with that part. Well I just try to write what comes up to me. Sorry for likers of clear thoughts.
Question of calculating the chromatic number of the plane
https://arxiv.org/pdf/1505.04159
Duminil-Copin's review at CDM conference
Some of Sokal's papers for a combinatorialist view on antiferromagnetic Potts models: https://arxiv.org/abs/cond-mat/9910503 https://arxiv.org/abs/cond-mat/0004231 https://arxiv.org/abs/math/0503607
A physicist's survey of Potts models http://www.yaroslavvb.com/papers/wu-potts.pdf
Actually I found out the Duminil-Copin et al. have proved their conjecture that the phase transition is discontinuous (as long "physically" proved by physicists). http://www.ihes.fr/~duminil/publi/2016discontinuousPotts.pdf