Constrained molecular dynamics simulation of the quark-gluon plasma

We calculate the equation of state of a quark system interacting through a phenomenological potential, the Richardson's potential, at finite baryon density and zero temperature. In particular we study three different cases with different quark masses ($u$ and $d$), and different assumptions for the potential at large distances. We solve molecular dynamics with a constraint due to Pauli blocking and find evidence of a phase transition from ``nuclear'' to ``quark matter.'' The phase transition is analyzed also through the behavior of the $J∕\ensuremath{\Psi}$ embedded in the quark system. We show that the $J∕\ensuremath{\Psi}$ particle behaves as an order parameter.