By adopting the isothermal–isochoric ensemble instead of the isothermal–isobaric one, the density of the modified Lennard-Jones system could be specified unambiguously even under extremely low-pressure conditions. The Gibbs free energy of the gaseous state could thereby be calculated with sufficient accuracy. For practical applications of the free energy, the lowest correction in the virial expansion was found to be sufficient under pressure conditions below 0.001e/σ3, where e and σ are the energy and length parameters, respectively, of the modified Lennard-Jones potential. Comparison of this free energy with those for the condensed phases obtained hitherto gave well-defined sublimation points. These points then defined the equilibrium solid–gas phase boundary. The present study thus confirmed the location of not only a sublimation point found previously under marginally stable conditions of an isobaric ensemble, but also the triple point; its temperature and pressure were speculated to be 0.61e/kB and 0....