  We study an associating lattice gas (ALG) using Monte Carlo simulation and
solutions on Husimi lattices. In this model, the molecules have an
orientational degree of freedom and the interactions depend on the relative
orientations of nearest-neighbor molecules. We focus on the transition between
the high-density liquid (HDL) phase and the isotropic gas phase in the limit of
full occupancy ($\mu \to \infty$). Simulation results show a continuous phase
transition at $\tau_c=k_BT_c/\gamma=0.4763(1)$ (where $-\gamma$ is the bond
energy) between the low-temperature HDL phase, with a non-vanishing mean
orientation of the molecules, and the high-temperature isotropic phase. Results
for critical exponents and the Binder cumulant indicate that the transition
belongs to the three-state Potts model universality class, even though the ALG
Hamiltonian does not have the full permutation symmetry of the Potts model. In
contrast with simulation, the Husimi lattice results furnish a discontinuous
phase transition, characterized by a discontinuity of the nematic order
parameter. The transition temperatures ($\tau_c=0.51403$ and $0.51207$ for
trees built with triangles and hexagons, respectively) are slightly higher than
the one found via simulation. Since the Husimi lattice studies show that the
ALG phase diagram features a discontinuous gas-HDL line for finite $\mu$, three
possible scenarios arise. The first is that in the limit $\mu \to \infty$ the
first-order line ends in a critical point; the second is a change in the nature
of the transition at some finite chemical potential; the third is that the
entire line is one of continuous phase transitions. Results from other ALG
models and the fact that mean-field approximations show a discontinuous phase
transition for the three-state Potts model (known to possess a continuous
transition) lends some weight to the third alternative.
