Hybrid Monte Carlo Simulation of Graphene on the Hexagonal Lattice
  We present a method for direct hybrid Monte Carlo simulation of graphene on the hexagonal lattice. We compare the results of the simulation with exact results for a unit hexagonal cell system, where the Hamiltonian can be solved analytically. 
R. C. Brower, C. Rebbi, D. Schaich
hep-lat cond-mat.str-el