New method for atomistic modeling of the microstructure of activated carbons using hybrid reverse Monte Carlo simulation

We propose a new hybrid reverse Monte Carlo (HRMC) procedure for atomistic modeling of the microstructure of activated carbons whereby the guessed configuration for the HRMC construction simulation is generated using the characterization results (pore size and pore wall thickness distributions) obtained by the interpretation of argon adsorption at 87 K using our improved version of the slit-pore model, termed the finite wall thickness (FWT) model (Nguyen, T. X.; Bhatia, S. K. Langmuir 2004, 20, 3532). This procedure overcomes limitations arising from the use of short-range potentials in the conventional HRMC method, which make the latter unsuitable for carbons such as activated carbon fibers that are anisotropic with medium-range ordering induced by their complex pore structure. The newly proposed approach is applied specifically for the atomistic construction of an activated carbon fiber ACF15, provided by Kynol Corporation (Nguyen, T. X.; Bhatia, S. K. Carbon 2005, 43, 775). It is found that the PSD of the ACF15's constructed microstructure is in good agreement with that determined using argon adsorption at 87 K. Furthermore, we have also found that the use of the Lennard-Jones (LJ) carbon-fluid interaction well depth obtained from scaling the flat graphite surface-fluid interaction well depth taken from Steele (Steele, W. A. Surf. Sci. 1973, 36, 317) provides an excellent prediction of experimental adsorption data including the differential heat of adsorption of simple gases (Ar, N2, CH4, CO2) over a wide range of temperatures and pressures. This finding is in agreement with the enhancement of the LJ carbon-fluid well depth due to the curvature of the carbon surface, found by the use of ab initio calculations (Klauda, J. B.; Jiang, J.; Sandler, S. I. J. Phys. Chem. B 2004, 108, 9842)