My system comprises of oppositely charged polymers, counter-ions and monovalent salts ions.All charges are monovalent(+1/-1). I am done with the basic simulations for the dilute phase calculating thE RDFs For the process. In order to calculate phase change, what should be implemented in the system. I have implemented overlapping, bondconstraint criteria, long range electrostatic energy, short range and bending potentials of polymer chain. My system is off lattice and simulation box is 15*!5*15. I am working in C++