How to set system in simulation so get the accurate density

GROMACS version:
GROMACS modification: Yes/No
I’m working on md simulation of ionic liquid.the density got from equlibrition is not equal to reported one.how to set the system so get the correct density.i crossed checked all bonded and non bonded parameters.all are correct.