Monte Carlo simulations of the NJL model near the nonzero temperature   phase transition
  We present results from numerical simulations of the Nambu-Jona-Lasinio model with an SU(2)xSU(2) chiral symmetry and N_c=4,8, and 16 quark colors at nonzero temperature. We performed the simulations by utilizing the hybrid Monte Carlo and hybrid Molecular Dynamics algorithms. We show that the model undergoes a second order phase transition. The critical exponents measured are consistent with the classical 3d O(4) universality class and hence in accordance with the dimensional reduction scenario. We also show that the Ginzburg region is suppressed by a factor of 1/N_c in accordance with previous analytical predictions. 
Costas Strouthos and Stavros Christofi