Equilibrium Properties of a Model Electrolyte Adsorbed in Quenched Disordered Charged Media: the ROZ Theory and GCMC Simulations

Fluids in porous media. IV. Quench effect on chemical potential

It appears to be a common sense to measure the crowdedness of a fluid system by the densities of the species constituting it. In the present work, we show that this ceases to be valid for confined fluids under some conditions. A quite thorough investigation is made for a hard sphere (HS) fluid adsorbed in a hard sphere matrix (a quench-annealed system) and its corresponding equilibrium binary mixture. When fluid particles are larger than matrix particles, the quench-annealed system can appear much more crowded than its corresponding equilibrium binary mixture, i.e., having a much higher fluid chemical potential, even when the density of each species is strictly the same in both systems, respectively. We believe that the insight gained from this study should be useful for the design of functionalized porous materials.

A three dimensional integral equation approach for fluids under confinement: Argon in zeolites

In this work, we explore the ability of an inhomogeneous integral equation approach to provide a full three dimensional description of simple fluids under conditions of confinement in porous media. Explicitly, we will consider the case of argon adsorbed into silicalite-1, silicalite-2, and an all-silica analogue of faujasite, with a porous structure composed of linear (and zig-zag in the case of silicalite-1) channels of 5-8 Å diameter. The equation is based on the three dimensional Ornstein-Zernike approximation proposed by Beglov and Roux [J. Chem. Phys. 103, 360 (1995)] in combination with the use of an approximate fluid-fluid direct correlation function furnished by the replica Ornstein-Zernike equation with a hypernetted chain closure. Comparison with the results of grand canonical Monte Carlo/molecular dynamics simulations evidences that the theory provides an accurate description for the three dimensional density distribution of the adsorbed fluid, both at the level of density profiles and bidimensional density maps across representative sections of the porous material. In the case of very tight confinement (silicalite-1 and silicalite-2), solutions at low temperatures could not be found due to convergence difficulties, but for faujasite, which presents substantially larger channels, temperatures as low as 77 K are accessible to the integral equation. The overall results indicate that the theoretical approximation can be an excellent tool to characterize the microscopic adsorption behavior of porous materials.

A computational study of electrolyte adsorption in a simple model for intercalated clays

A pillared interlayered clay is represented by a two-dimensional quenched charged disordered medium, in which the pillar configuration is produced by the quench of a two-dimensional electrolyte and the subsequent removal of the anions (that act as a template). The cation charge is counterbalanced by a neutralizing background that is an ideal representation of the layer's negative charge in the experimental system. In this paper we investigate the adsorption of electrolyte particles in this charged disordered medium resorting both to the use of the replica Ornstein-Zernike equation in the hypernetted chain approximation and grand canonical Monte Carlo simulations. The theoretical approach qualitatively reproduces the simulated behavior of the adsorbed fluids. Theoretical estimates of the material porosities obtained for various types of pillar distributions are in good agreement with the simulation. We investigate the influence of the matrix on correlation functions and adsorption isotherms.

Electrolyte exclusion from charged adsorbent: replica Ornstein-Zernike theory and simulations

Structural and thermodynamic properties of the restrictive primitive model +1:-1 electrolyte solution adsorbed in a disordered charged media were studied by means of the Grand Canonical Monte Carlo simulation and the replica Ornstein-Zernike theory. Disordered media (adsorbent, matrix) was represented by a distribution of negatively charged hard spheres frozen in a particular equilibrium distribution. The annealed counterions and co-ions were assumed to be distributed within the nanoporous adsorbent in thermodynamic equilibrium with an external reservoir of the same electrolyte. In accordance with the primitive model of electrolyte solutions, the solvent was treated as a dielectric continuum. The simulations were performed for a set of model parameters, varying the net charge of the matrix (i.e., concentrations of matrix ions) and of annealed electrolyte, in addition to the dielectric constant of the invading solution. The concentration of adsorbed electrolyte was found to be lower than the corresponding concentration of the equilibrium bulk solution. This electrolyte "exclusion" depends strongly on the dielectric constant of the invading solution, as also on concentrations of all components. The most important parameter is the net charge of the matrix. Interestingly, the electrolyte rejection decreases with increasing Bjerrum length for the range of parameters studied here. The latter finding can be ascribed to strong inter-ionic correlation in cases where the Bjerumm length is high enough. To a minor extent, the adsorption also depends on the spacial distribution of fixed charges in adsorbent material. The replica Ornstein-Zernike theory was modified to cater for this model and tested against the computer simulations. For the range of parameters explored in this work, the agreement between the two methods is very good. These calculations were also compared with the results of the classical Donnan theory for electrolyte exclusion.

A stepwise approximation for modeling of the wall–fluid potential of a mesoscopic pore

A common computational method for the characterization of porous materials is to calculate the adsorption isotherm of fluids in the materials from the pre-assumed wall-fluid potential. If the wall-fluid potential is unknown, the common computational method becomes invalid. In a realistic experiment, however, it is common to know the experimental adsorption isotherm of nitrogen and not to know the wall-fluid potential. Here we propose a stepwise approximation for modeling wall-fluid potential under conditions where only the adsorption isotherm of nitrogen is measured experimentally. Based on the modeled wall-fluid potential, we can characterize the porous materials and predict the adsorption of other adsorbates on the materials. It is expected that the approach would provide a powerful means for the characterization of novel materials under conditions where only the experimental adsorption isotherm is available.

Fluids in porous media. I. A hard sponge model

The morphology of many porous materials is spongelike. Despite the abundance of such materials, simple models which allow for a theoretical description of these materials are still lacking. Here, we propose a hard sponge model which is made by digging spherical cavities in a solid continuum. We found an analytical expression for describing the interaction potential between fluid particles and the spongelike porous matrix. The diagrammatic expansions of different correlation functions are derived as well as that of grand potential. We derived also the Ornstein-Zernike (OZ) equations for this model. In contrast to Madden-Glandt model of random porous media [W. G. Madden and E. D. Glandt, J. Stat. Phys. 51, 537 (1988)], the OZ equations for a fluid confined in our hard sponge model have some similarity to the OZ equations of a three-component fluid mixture. We show also how the replica method can be extended to study our sponge model and that the same OZ equations can be derived also from the extended replica method.

Integral equation study of a Stockmayer fluid adsorbed in polar disordered matrices

Based on replica integral equations in the (reference) hypernetted chain approximation we investigate the structural features and phase properties of a dipolar Stockmayer fluid confined to a disordered dipolar matrix. The integral equations are applied to the homogeneous high-temperature phase where the system is globally isotropic. At low densities we find the influence of dipolar interactions between fluid (f) and matrix (m) particles to be surprisingly similar to the previously investigated effect of attractive isotropic (fm) interactions: the critical temperature of the vapor-liquid transition decreases with increasing (fm) coupling, while the critical density increases. The anisotropic nature of the dipolar (fm) interactions turns out to play a more dominant role at high fluid densities where we observe a pronounced sensitivity in the dielectric constant and a strong degree of local orientational ordering of the fluid particles along the local fields generated by the matrix. Moreover, an instability of the dielectric constant, which is a precursor of ferroelectric ordering occurring both in bulk Stockmayer fluids and in fluids in nonpolar matrices, is observed only for very small dipolar (fm) couplings.