On the existence of a hysteresis loop in open and closed end pores

Abstract

We have studied the adsorption of argon at 87 K in slit pores of finite length with a smooth graphitic potential, open at both ends or closed at one end. Simulations were carried out using conventional GCMC (grand canonical Monte Carlo) or kMC (kinetic Monte Carlo) in the canonical ensemble with extremely long Markov chain, of at least 2 × 108 configurations; selected simulations with much longer Markov chains do not show any change in the results. When the pore width is in the micropore range (0.65 nm), type I isotherms are obtained for both pore models and for both simulation methods. However, wider pores (1, 2 and 3 nm in width) all exhibit hysteresis loops in the GCMC simulations, while in the canonical ensemble simulations, the isotherms pass through a sigmoid van der Waals type loop in the transition region. This loop locates the true equilibrium transition. For the pores with one closed end, this transition is close to, or coincides with, the adsorption branch of the GCMC hysteresis loop, but for the open-ended pores, it is more closely associated with the desorption branch. In a separate study of adsorption hysteresis in an infinitely long slit pore, using both simulation techniques, the van der Waals loop follows the adsorption branch of the GCMC isotherm to the transition, then reverts to a long vertical section that falls midway between the two hysteresis branches and finally moves to the desorption transition close to the evaporation pressure. An examination of molecular distributions inside the pores reveals two coexisting phases in the canonical simulations, whereas in the grand canonical simulations, the molecules are uniformly distributed along the length of the pores.