The phase diagram of the Lennard-Jones fluid using temperature dependent interaction parameters

The forces of interaction between argon atoms can be described by the Lennard-Jones potential model. It is hypothesised that the use of temperature dependent interaction parameters, instead of using temperature independent interaction parameters, may lead to improvement in the prediction of the vapour–liquid coexistence curve. Published second virial coefficient data were used to fit a simple two-parameter temperature dependent model for the collision diameter and well depth. Vapour–liquid coexistence curve for argon was simulated in the NVT Gibbs ensemble Monte Carlo technique. The simulations were carried out using each of the temperature independent and temperature dependent parameters in the temperature range: 110–148 K. The critical temperature and density were determined using the Ising-scaling model. The results using temperature dependent parameters produce, overall, a more accurate phase diagram compared to the diagram generated using temperature independent interaction parameters. The root mean square deviation is reduced by 42.1% using temperature dependent interaction parameters. Also, there was no significant difference between the results obtained using temperature dependent interaction parameters and the highly accurate and computationally demanding phase diagrams based on three body contributions.