In physics, with the advent of topological materials, and in chemistry, within the scope of chemical graph theory, topological invariants and/or indices have been considered to successfully characterize innumerous systems. In particular, strong links have been identified (both numerically and analytically) between properties of the Ising model on a lattice L and features of the so-called spanning trees (STs) of L. Nontheless, studies exploring this connection tend to address only a handful of cases given the demands of the necessary calculations. But examining only a few instances prevents one from looking for general trends across numerous L's, which could eventually reveal universal traits. In this contribution, we present the most comprehensive investigation to date, analyzing the Ising-ST relation for all the L's belonging to the families
of
–uniform periodic tiling of the plane, in a total of 1248 lattices. With this goal, we develop optimized protocols (taking advantage of a recently proposed
representation for
) to compute for each L its ST constant λ and the Kac–Ward matrix. The determinant of the latter yields the Ising model free energy and consequently the critical temperature
. Then, considering the relatively large sample generated, we use machine learning techniques, which disclose a general correlation between the Ising critical temperature and the ST constant, described by a simple quadratic polynomial function P. As a benchmark, we test P for some arbitrary lattices (outside
), finding rather satisfactory fittings. These results point to a useful classification scheme for Ising
in 2D, demonstrating that λ can be a relevant topological concept to investigate lattice models. Finally, as a positive 'side-effect' of our computations, for these
's we confirm (and even improve) a conjectured inequality associating λ and the effective coordinator number κ of a lattice.