Geometric Formula for 2d Ising Zeros: Examples & Numerics

A geometric formula for the zeros of the partition function of the inhomogeneous 2d Ising model was recently proposed in terms of the angles of 2d triangulations embedded in the flat 3d space. Here we proceed to an analytical check of this formula on the cubic graph, dual to a double pyramid, and provide a thorough numerical check by generating random 2d planar triangulations. Our method is to generate Delaunay triangulations of the 2-sphere then performing random local rescalings. For every 2d triangulations, we compute the corresponding Ising couplings from the triangle angles and the dihedral angles, and check directly that the Ising partition function vanishes for these couplings (and grows in modulus in their neighborhood). In particular, we lift an ambiguity of the original formula on the sign of the dihedral angles and establish a convention in terms of convexity/concavity. Finally, we extend our numerical analysis to 2d toroidal triangulations and show that the geometric formula does not work and will need to be generalized, as originally expected, in order to accommodate for non-trivial topologies.