Hubbard-type models on the hexagonal lattice are of great interest, as they provide realistic descriptions of graphene and other related materials. Hybrid Monte Carlo simulations offer a first-principles approach to study their phase structure. Here, we review the present status of our work in this direction.