I am trying to generate a phase diagram of pressure versus temperature in a NPT ensemble using molecular dynamics method. I am using the Lennard Jones potential with 108 particles and Berendsen thermostat and barostat to fix the parameter.
In order to observe phase transition I plot the radial distribution function (RDF) and also density versus temperature graph for various fixed temperature. However it is hard to differentiate liquid and gaseous phase from RDF.
Another peculiar thing is that there should be a discontinuity in density as the transition temperature is reached, which I observe only for solid liquid transition but not for the liquid vapor transition. The change in density should be huge for liquid vapor transition. Any idea on what possibly goes wrong?
cut off radius = 2.5
Truncated and shifted Lennard Jones potential
number of particle = 108
number of steps = 2000
timestep = 0.004
Berendsen barostat of scale coeefficient of 0.1 is used to scale the size of simulation box
The density versus temperature graph for pressure of 1.0MPa:
The solid liquid transition can be seen but not the liquid vapor transition.