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 evidencies of a phase transition from "nuclear" to "quark matter", which is analyzed also through the behaviour of the $J/\Psi$ embedded in the quark system. We show that the $J/\Psi$ particle behaves as an order parameter.