Nonadditive Binary Hard Sphere Mixture in Disordered Hard Sphere Matrices: Integral Equations and Computer Simulation

Investigation of Fluid-Fluid and Solid-Solid Phase Separation of Symmetric Nonadditive Hard Spheres at High Density

We have calculated the values of the critical packing fractions for the mixtures of symmetric nonadditive hard spheres at high densities for small values of the nonadditivity parameter. Calculations have been performed for solid-solid and fluid-fluid demixing transitions. A cluster algorithm for Monte Carlo simulations in a semigrand ensemble was used, and the waste recycling method was applied to improve the accuracy of the calculations. The finite size scaling analysis was employed to compute the critical packing fractions for infinite systems with high accuracy.

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.

Phase separation of binary nonadditive hard sphere fluid mixture confined in random porous media

I analyze the fluid-fluid phase separation of nonadditive hard sphere fluid mixture absorbed in random porous media. An equation of state is derived by using the perturbation theory to this complex system with quenched disorders. The results of this theory are in good agreement with those obtained from semi-grand canonical ensemble Monte Carlo simulations. The contact value of the fluid-fluid radial distribution functions of the reference which is the key point of the perturbation process is derived as well, the comparison against Monte Carlo simulations shows that it has an excellent accuracy.

Fluids in porous media: the case of neutral walls

The bulk phase behavior of a fluid is typically altered when the fluid is brought into confinement by the walls of a random porous medium. Inside the porous medium, phase-transition points are shifted, or may disappear altogether. A crucial determinant is how the walls interact with the fluid particles. In this work, we consider the situation whereby the walls are neutral with respect to the liquid and vapor phases. In order to realize the condition of strict neutrality, we use a symmetric binary mixture inside a porous medium that interacts identically with mixture species. Monte Carlo simulations are then used to obtain the phase behavior. Our main finding is that, in the presence of the porous medium, a liquid-vapor critical point still exists. At the critical point, the distribution of the order parameter remains scale invariant, but self-averaging is violated. These findings provide further evidence that random confinement by neutral walls induces critical behavior of the random Ising model (i.e., Ising models with dilution type disorder, where the disorder couples to the energy).

Colloid-polymer mixtures in the presence of quenched disorder: a theoretical and computer simulation study

We use theory and computer simulation to study the structure and phase behavior of colloid-polymer mixtures in the presence of quenched disorder. The Asakura-Oosawa model (AO) (Asakura and Oosawa 1954 J. Chem. Phys. 22 1255) is used to describe the colloid-colloid, colloid-polymer, and polymer-polymer pair interactions. We then investigate the behavior of this model in the presence of frozen-in (quenched) obstacles. The obstacles will be placed according to two different scenarios, both of which are experimentally feasible. In the first scenario, polymers are distributed at positions drawn from an ideal gas configuration. In the second scenario, colloidal particles are distributed at positions drawn from an equilibrium hard sphere configuration. We investigate how the unmixing transition of the AO model is affected by the type of quenched disorder. The theoretical formalism is based on the replica method of Given and Stell (1994 Physica A 209 495). Our foremost aim is to test the accuracy of three common closures to the replica Ornstein-Zernike equations, namely the hypernetted chain, the Percus-Yevick, and the Martinov-Sarkisov equations. The accuracy is determined by comparison with grand canonical Monte Carlo simulations. We find that, for quenched polymer disorder, all three closures perform remarkably well. However, when quenched colloid disorder is considered, i.e. the second mentioned scenario, the predictions of all three closures worsen dramatically.

Hybrid Density Functional Theory for Homopolymer Mixtures Confined in a Selective Nanoslit

By integrating polymer density function theory (DFT) and single-chain molecular simulation, a hybrid DFT is developed for homopolymer mixtures confined in a selective nanoslit. Two weighting functions are adopted separately in the polymer DFT for repulsive and attractive contributions to the excess free energy functional. The theoretical results agree well with simulation data for the density profiles, configurations (tail, loop and train), adsorption amounts, layer thicknesses, and partition coefficients. The polymer-slit interaction is found to have a large effect on the density profiles and partition coefficients but is found to have a small effect on the average sizes and percentages of the configurations. Nearly half of the polymer segments form tails, and the other half form trains. In addition, bridges are observed to form for sufficiently long polymer chains. As the length difference between two polymers increases, the effect of chain connectivity becomes increasingly important.

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.