A thermodynamic study of the mid-density scheme to determine the equilibrium phase transition in cylindrical pores

Computational Simulations of Nanoconfined Argon Film through Adsorption–Desorption in a Uniform Slit Pore

We performed hybrid grand canonical Monte Carlo/molecular dynamics (GCMC/MD) simulations to investigate the adsorption-desorption isotherms of argon molecules confined between commensurate and incommensurate contacts in nanoscale thickness. The recently proposed mid-density scheme was applied to the obtained hysteresis loops to produce a realistic equilibrium phase of nanoconfined fluids. The appropriate chemical potentials can be determined if the equilibrium structures predicted by GCMC/MD simulations are consistent with those observed in previously developed liquid-vapor molecular dynamics (LVMD) simulations. With the chemical potential as input, the equilibrium structures obtained by GCMC/MD simulations can be used as reasonable initial configurations for future metadynamics free energy calculations.

Determination of phase equilibria in confined systems by open pore cell Monte Carlo method

We present a modification of the molecular dynamics simulation method with a unit pore cell with imaginary gas phase [M. Miyahara, T. Yoshioka, and M. Okazaki, J. Chem. Phys. 106, 8124 (1997)] designed for determination of phase equilibria in nanopores. This new method is based on a Monte Carlo technique and it combines the pore cell, opened to the imaginary gas phase (open pore cell), with a gas cell to measure the equilibrium chemical potential of the confined system. The most striking feature of our new method is that the confined system is steadily led to a thermodynamically stable state by forming concave menisci in the open pore cell. This feature of the open pore cell makes it possible to obtain the equilibrium chemical potential with only a single simulation run, unlike existing simulation methods, which need a number of additional runs. We apply the method to evaluate the equilibrium chemical potentials of confined nitrogen in carbon slit pores and silica cylindrical pores at 77 K, and show that the results are in good agreement with those obtained by two conventional thermodynamic integration methods. Moreover, we also show that the proposed method can be particularly useful for determining vapor-liquid and vapor-solid coexistence curves and the triple point of the confined system.