The interactions between particles are described by a virial expansion up to the fourth order in our model. Note, the expansion at least up to the third order is required to guarantee a physically relevant solution for the negative second virial coefficient. The second viral coefficient B ij 2 between the particles i and j can be easily evaluated from the Lennard-Jones pair potential