“…These include three cases consisting of N hex = n × n hexagons and N = 3 n 2 + 4 n − 1 sites, and two hexagonal-shaped samples containing N hex = 7 and 19 hexagons (N = 30 and 72 spins, respectively). As in previous work [41] for periodic boundary conditions, we started by enumerating all allowed n -loop configurations { i }. We then make use of the fact that, apart from a global 2In the case of periodic boundary conditions, there are actually linear relations between the hard-hexagon states on the torus, but from the point of view of counting, this deficit is compensated by the state with ω − ( ì 0) = ω 0 , see also figure 1 [36,41].…”