Comparison of several classical density functional theories for the adsorption of flexible chain molecules into cylindrical nanopores

Heat of adsorption, adsorption stress, and optimal storage of methane in slit and cylindrical carbon pores predicted by classical density functional theory

Temperature, pressure and pore-size dependences of the heat of adsorption, adsorption stress, and adsorption capacity of methane in simple models of slit and cylindrical carbon pores are studied using classical density functional theory (CDFT) and grand-canonical Monte-Carlo (MC) simulation. Studied properties depend nontrivially on the bulk pressure and the size of the pores. Heat of adsorption increases with loading, but only for sufficiently narrow pores. While the increase is advantageous for gas storage applications, it is less significant for cylindrical pores than for slits. Adsorption stress and the average adsorbed fluid density show oscillatory dependence on the pore size and increase with bulk pressure. Slit pores exhibit larger amplitude of oscillations of the normal adsorption stress with pore size increase than cylindrical pores. However, the increase of the magnitude of the adsorption stress with bulk pressure increase is more significant for cylindrical than for slit pores. Adsorption stress appears to be negative for a wide range of pore sizes and external conditions. The pore size dependence of the average delivered density of the gas is analyzed and the optimal pore sizes for storage applications are estimated. The optimal width of slit pore appears to be almost independent of storage pressure at room temperature and pressures above 10 bar. Similarly to the case of slit pores, the optimal radius of cylindrical pores does not exhibit much dependence on the storage pressure above 15 bar. Both optimal width and optimal radii of slit and cylindrical pores increase as the temperature decreases. A comparison of the results of CDFT theory and MC simulations reveals subtle but important differences in the underlying fluid models employed by the approaches. The differences in the high-pressure behaviour between the hard-sphere 2-Yukawa and Lennard-Jones models of methane, employed by the CDFT and MC approaches, respectively, result in an overestimation of the heat of adsorption by the CDFT theory at higher loadings. However, both adsorption stress and adsorption capacity appear to be much less sensitive to the differences between the models and demonstrate excellent agreement between the theory and the computer experiment.

Adsorption and Phase Behavior of Pure/Mixed Alkanes in Nanoslit Graphite Pores: An iSAFT Application

The prediction of fluid phase behavior in nanoscale pores is critical for shale gas/oil development. In this work, we use a molecular density functional theory (DFT) to study the effect of molecular size and shape on partitioning to graphite nanopores as a model of shale. Here, interfacial statistical associating fluid theory (iSAFT) is applied to model alkane (C₁ - C₈) adsorption/desorption/phase behavior in graphite slit pores for both pure fluids and mixtures. The pure component parameters were fit to the bulk saturated liquid density and vapor pressure data in selected temperature ranges. The potential of interaction between the fluid and graphite is modeled with a Steele 10-4-3 potential that is fit to the potential of mean force from single-molecule simulations. Good agreement is found between theory and molecular simulation for the density distributions of pure components in slit pores. The critical properties of methane, ethane, and their mixtures as well as the shift in bubble point and dew point densities were studied, showing good agreement with simulation. The competitive adsorption of mixtures of normal and branched alkanes in graphite pores was also studied. Heavier components more strongly adsorb up to the point that the entropic penalty due to confinement reduces adsorption.

Exponential approximation for one-component Yukawa plasma

A theory based on the exponential approximation of the liquid-state theory is applied to study properties of several models of one-component Yukawa plasma characterized by different values of the screening parameter z. The results of the new theory are compared to the results of a conventional theory, which is based on the first-order mean spherical approximation, and to the results of a Monte Carlo simulation. The new theory shows improvements in the predictions for the thermodynamic and structural properties of Yukawa plasmas with high and intermediate values of the screening parameter, z, and coupling parameter, Γ. For low values of z and Γ, the new theory is comparable in accuracy to the conventional theory, which in turn agrees well with the results of the Monte Carlo simulation.