Freezing of Lennard-Jones fluid in cylindrical nanopores under tensile conditions

We have shown the Lennard-Jones (LJ) phase diagram for a slit-shaped nanopore by molecular simulations and thermodynamically predicted the results with no adjustable parameter. With this success, LJ phase diagrams are predictable. In this study, the freezing of an LJ CH₄ capillary condensate under a tensile condition in a nonstructural carbon nanopore with a cylindrical geometry was examined using molecular dynamics (MD) simulation. We employ a unit cell in contact with a bulk vapor phase, which is useful for the determination of the bulk vapor pressure in equilibrium with the molecules in a pore. The MD simulation results show liquid-solid (amorphous) phase transitions with a variation in the bulk vapor pressure. The frozen particles are arranged in concentric circular regions along the wall similar to those reported by Maddox and Gubbins. The freezing points are determined from the variations in density, enthalpy, arrangement, and structural functions. The obtained liquid-solid coexistence points are found to exhibit a significant dependence of the freezing point on the equilibrium bulk vapor pressure, forming an extraordinarily skewed curve on the p-T diagram, in contrast to the bulk phase coexistence that is represented by an almost vertical line. The origin of the significant dependence is considered to be the Laplace effect on the capillary condensate similar to the case with a slit-shaped pore. A simple model, which the authors earlier presented for slit-shaped nanopores, successfully predicted the p-T relation of the freezing point for cylindrical nanopores as well.