Improved modeling of two-dimensional transitions in dense phases on crystalline surfaces. Krypton-graphite system

This paper presents a refined technique to describe two-dimensional phase transitions in dense fluids adsorbed on a crystalline surface. Prediction of parameters of 2D liquid-solid equilibrium is known to be an extremely challenging problem, which is mainly due to a small difference in thermodynamic functions of coexisting phases and lack of accuracy of numerical experiments in case of their high density. This is a serious limitation of various attempts to circumvent this problem. To improve this situation, a new methodology based on the kinetic Monte Carlo method was applied. The methodology involves analysis of equilibrium gas-liquid and gas-solid systems undergoing an external potential, which allows gradual shifting parameters of the phase coexistence. The interrelation of the chemical potential and tangential pressure for each system is then treated with the Gibbs-Duhem equation to obtain the point of intersection corresponding to the liquid/solid-solid equilibrium coexistence. The methodology is demonstrated on the krypton-graphite system below and above the 2D critical temperature. Using experimental data on the liquid-solid and the commensurate-incommensurate transitions in the krypton monolayer derived from adsorption isotherms, the Kr-graphite Lennard-Jones parameters have been corrected resulting in a higher periodic potential modulation.