Relationships between self-diffusivity, packing fraction, and excess entropy in simple bulk and confined fluids

Static measures such as density and entropy, which are intimately connected to structure, have featured prominently in modern thinking about the dynamics of the liquid state. Here, we explore the connections between self-diffusivity, density, and excess entropy for two of the most widely used model "simple" liquids, the equilibrium Lennard-Jones and square-well fluids, in both bulk and confined environments. We find that the self-diffusivity data of the Lennard-Jones fluid can be approximately collapsed onto a single curve (i) versus effective packing fraction and (ii) in appropriately reduced form versus excess entropy, as suggested by two well-known scaling laws. Similar data collapse does not occur for the square-well fluid, a fact that can be understood on the basis of the nontrivial effects that temperature has on its static structure. Nonetheless, we show that the implications of confinement for the self-diffusivity of both of these model fluids, over a broad range of equilibrium conditions, can be predicted on the basis of knowledge of the bulk fluid behavior and either the effective packing fraction or the excess entropy of the confined fluid. Excess entropy is perhaps the most preferable route due to its superior predictive ability and because it is a standard, unambiguous thermodynamic quantity that can be readily predicted via classical density functional theories of inhomogeneous fluids.     

Does confining the hard-sphere fluid between hard walls change its average properties?

We use grand canonical transition-matrix Monte Carlo and discontinuous molecular dynamics simulations to generate precise thermodynamic and kinetic data for the equilibrium hard-sphere fluid confined between smooth hard walls. These simulations show that the pronounced inhomogeneous structuring of the fluid normal to the confining walls, often the primary focus of density functional theory studies, has a negligible effect on many of its average properties over a surprisingly broad range of conditions. We present one consequence of this insensitivity to confinement: a simple analytical equation relating the average density of the confined fluid to that of the bulk fluid with equal activity. Nontrivial implications of confinement for average fluid properties do emerge in this system, but only when the fluid is both (i) dense and (ii) confined to a gap smaller than approximately three particle diameters. For this limited set of conditions, we find that "in-phase" oscillatory deviations in excess entropy and self-diffusivity (relative to the behavior of the bulk fluid at the same average density) occur as a function of gap size. These paired thermodynamic/kinetic deviations from bulk behavior appear to reflect the geometric packing frustration that arises when the confined space cannot naturally accommodate an integer number of particle layers.     

Determination of surface tension in binary mixtures using transition-matrix Monte Carlo

We present a methodology based on grand-canonical transition-matrix Monte Carlo and finite-size scaling analysis to calculate surface tensions in binary mixtures. In particular, mixture transition-matrix Monte Carlo is first used to calculate apparent, system-size-dependent free-energy barriers separating coexisting fluid phases. Finite-size scaling is then used to extrapolate these values to the infinitely large system limit to determine the true thermodynamic surface tension. A key distinction of the methodology is that it yields the entire isothermal surface-tension curve for a binary mixture in a relatively small number of simulations. We demonstrate the utility of the method by calculating surface-tension curves for three binary Lennard-Jones mixtures. While we have only examined the surface tension of simple fluids in this work, the method is general and can be extended to molecular fluids as well as to determine interfacial tensions of liquid-liquid interfaces.     

Transport in nanoporous zeolites: relationships between sorbate size, entropy, and diffusivity

Molecular dynamics simulations have been performed on monatomic sorbates confined within zeolite NaY to obtain the dependence of entropy and self-diffusivity on the sorbate diameter. Previously, molecular dynamics simulations by Santikary and Yashonath [J. Phys. Chem. 98, 6368 (1994)], theoretical analysis by Derouane et al. [J. Catal. 110, 58 (1988)] as well as experiments by Kemball [Adv. Catal. 2, 233 (1950)] found that certain sorbates in certain adsorbents exhibit unusually high self-diffusivity. Experiments showed that the loss of entropy for certain sorbates in specific adsorbents was minimum. Kemball suggested that such sorbates will have high self-diffusivity in these adsorbents. Entropy of the adsorbed phase has been evaluated from the trajectory information by two alternative methods: two-phase and multiparticle expansion. The results show that anomalous maximum in entropy is also seen as a function of the sorbate diameter. Further, the experimental observation of Kemball that minimum loss of entropy is associated with maximum in self-diffusivity is found to be true for the system studied here. A suitably scaled dimensionless self-diffusivity shows an exponential dependence on the excess entropy of the adsorbed phase, analogous to excess entropy scaling rules seen in many bulk and confined fluids. The two trajectory-based estimators for the entropy show good semiquantitative agreement and provide some interesting microscopic insights into entropy changes associated with confinement.     

Determination of fluid-phase behavior using transition-matrix Monte Carlo: binary Lennard-Jones mixtures

We present a novel computational methodology for determining fluid-phase equilibria in binary mixtures. The method is based on a combination of highly efficient transition-matrix Monte Carlo and histogram reweighting. In particular, a directed grand-canonical transition-matrix Monte Carlo scheme is used to calculate the particle-number probability distribution, after which histogram reweighting is used as a postprocessing procedure to determine the conditions of phase equilibria. To validate the methodology, we have applied it to a number of model binary Lennard-Jones systems known to exhibit nontrivial fluid-phase behavior. Although we have focused on monatomic fluids in this work, the method presented here is general and can be easily extended to more complex molecular fluids. Finally, an important feature of this method is the capability to predict the entire fluid-phase diagram of a binary mixture at fixed temperature in a single simulation.     

Surface tension and vapor-liquid phase coexistence of confined square-well fluid

Phase equilibria of a square-well fluid in planar slit pores with varying slit width are investigated by applying the grand-canonical transition-matrix Monte Carlo (GC-TMMC) with the histogram-reweighting method. The wall-fluid interaction strength was varied from repulsive to attractive such that it is greater than the fluid-fluid interaction strength. The nature of the phase coexistence envelope is in agreement with that given in literature. The surface tension of the vapor-liquid interface is calculated via molecular dynamics simulations. GC-TMMC with finite size scaling is also used to calculate the surface tension. The results from molecular dynamics and GC-TMMC methods are in very good mutual agreement. The vapor-liquid surface tension, under confinement, was found to be lower than the bulk surface tension. However, with the increase of the slit width the surface tension increases. For the case of a square-well fluid in an attractive planar slit pore, the vapor-liquid surface tension exhibits a maximum with respect to wall-fluid interaction energy. We also report estimates of critical properties of confined fluids via the rectilinear diameter approach.     

Hard sphere fluids at a soft repulsive wall: a comparative study using Monte Carlo and density functional methods

Hard-sphere fluids confined between parallel plates at a distance D apart are studied for a wide range of packing fractions including also the onset of crystallization, applying Monte Carlo simulation techniques and density functional theory. The walls repel the hard spheres (of diameter σ) with a Weeks-Chandler-Andersen (WCA) potential V(WCA)(z) = 4ε[(σ(w)/z)(12) - (σ(w)/z)(6) + 1/4], with range σ(w) = σ/2. We vary the strength ε over a wide range and the case of simple hard walls is also treated for comparison. By the variation of ε one can change both the surface excess packing fraction and the wall-fluid (γ(wf)) and wall-crystal (γ(wc)) surface free energies. Several different methods to extract γ(wf) and γ(wc) from Monte Carlo (MC) simulations are implemented, and their accuracy and efficiency is comparatively discussed. The density functional theory (DFT) using fundamental measure functionals is found to be quantitatively accurate over a wide range of packing fractions; small deviations between DFT and MC near the fluid to crystal transition need to be studied further. Our results on density profiles near soft walls could be useful to interpret corresponding experiments with suitable colloidal dispersions.     

Confinement Induced Ordering in Fluid of Hard Ellipsoids

The ordering configurations of a fluid of anisotropic ellipsoids under the confinement of two apposing impenetrable walls are studied by Monte Carlo simulations. The excess adsorption of the fluid on the walls with respect to the aspect ratio has a maximum at the critical aspect ratio of 2.9 in high-density ellipsoid fluids, indicating an orientational ordering in the adjacent region of the walls, which is confirmed by probing into the density configurations and the orientational order parameter in the adjacent region of the walls for varying aspect ratios. In addition, the orientational order parameter in the bulk fluid at the same density is calculated, and it indicates an isotropic state as the bulk density is still below the bulk isotropic-to-nematic transition. Therefore, it can be concluded that the anisotropic ordering near the walls in the ellipsoid fluid that exhibits isotropic in the bulk is induced by the confinement effect of the walls.     

Molecular simulation study of effect of molecular association on vapor-liquid interfacial properties

Vapor-liquid interfacial properties of square-well associating fluids are studied via transition-matrix Monte Carlo simulation. Results for one-site and two-site association models are presented. Coexistence properties, surface tension, cluster distribution, density profile, and orientation profile are presented. Molecular association affects the interfacial properties and cluster fractions more than it affects the bulk densities. We observe that the surface tension exhibits a maximum with respect to association strength. This behavior is in agreement with the recent study of Peery and Evans for one site system using a square-gradient approach.     

Relaxation and short time dynamics of bulk liquids and fluids confined in spherical cavities and slit pores

The density of states for bulk and confined fluids have been modeled using a recently proposed gamma distribution (Krishnan, S. H.; Ayappa, K. G. J. Chem. Phys. 2004, 121, 3197). The gamma distribution results in a closed form analytical expression for the velocity autocorrelation function and the relaxation time of the fluid. The two parameters of the gamma distribution are related analytically to the second and fourth frequency moments of the fluid using short time expansions. The predictions by the proposed gamma model are compared with the velocity autocorrelation functions obtained using the theory of instantaneous normal modes (INMs) and from molecular dynamics simulations. The model is applied to a bulk soft sphere liquid and fluids confined in a spherical cavity and slit-shaped pores. The gamma model is able to capture the resulting changes in relaxation time due to changes in density and temperature extremely well for both the bulk liquid and confined inhomogeneous fluid situations. In all cases, the predictions by the gamma model are superior to those obtained from the INM theory. In the case of the fluid confined in a slit pore, the loadings were obtained from a grand canonical Monte Carlo simulation where the pore is equilibrated with a bulk fluid. This is similar to a confinement situation in a surface force apparatus. The predicted relaxation times vs pore widths from the gamma model are seen to accurately capture the oscillations due to formation and disruption of layers within the slit pore.     

Phase separation of the Widom-Rowlinson mixture in confined geometry

The behavior of binary Widom-Rowlinson [J. Chem. Phys. 52, 1670 (1970)] mixture in confined geometry is investigated. The influence of confinement on phase separation is examined. Coexistence curves for the mixture in slitlike pores and pores of complex geometry are calculated in Monte Carlo simulations. Finite size scaling analysis is used to determine precisely the location of critical point for the bulk mixture.     

Phase equilibria and interfacial tension of fluids confined in narrow pores

Correlation between phase behaviors of a Lennard-Jones fluid in and outside a pore is examined over wide thermodynamic conditions by grand canonical Monte Carlo simulations. A pressure tensor component of the confined fluid, a variable controllable in simulation but usually uncontrollable in experiment, is related with the pressure of a bulk homogeneous system in equilibrium with the confined system. Effects of the pore dimensionality, size, and attractive potential on the correlations between thermodynamic properties of the confined and bulk systems are clarified. A fluid-wall interfacial tension defined as an excess grand potential is evaluated as a function of the pore size. It is found that the tension decreases linearly with the inverse of the pore diameter or width.     

Combining reactive and configurational-bias Monte Carlo: confinement influence on the propene metathesis reaction system in various zeolites

In order to efficiently calculate chemical equilibria of large molecules in a confined environment the reactive Monte Carlo technique is combined with the configurational-bias Monte Carlo approach. To prove that detailed balance is fulfilled the acceptance rule for this combination of particular Monte Carlo techniques is derived in detail. Notably, by using this derivation all other acceptance rules of any Monte Carlo trial moves usually carried out in combination with the configurational-bias Monte Carlo approach can be deduced from it. As an application of the combination of reactive and configurational-bias Monte Carlo the influence of different zeolitic confinements (MFI, TON, LTL, and FER) on the reaction equilibrium and the selectivity of the propene metathesis reaction system was investigated. Compared to the bulk phase the conversion is increased significantly. The authors study this reaction system in the temperature range between 300 and 600 K, and the pressure range from 1 to 7 bars. In contrast to the bulk phase, pressure and temperature have a strong influence on the composition of the reaction mixture in confinement. At low pressures and temperatures both conversion and selectivity are highest. Furthermore, the equilibrium composition is strongly dependent on the type of zeolite. This demonstrates the important role of the host structure in catalytic systems.     

Chemical reaction equilibrium in nanoporous materials: NO dimerization reaction in carbon slit nanopores

We present a molecular-level simulation study of the effects of confinement on chemical reaction equilibrium in nanoporous materials. We use the reaction ensemble Monte Carlo (RxMC) method to investigate the effects of temperature, nanopore size, bulk pressure, and capillary condensation on the nitric oxide dimerization reaction in a model carbon slit nanopore in equilibrium with a bulk reservoir. In addition to the RxMC simulations, we also utilize the molecular-dynamics method to determine self-diffusion coefficients for confined nonreactive mixtures of nitric oxide monomers and dimers at compositions obtained from the RxMC simulations. We analyze the effects of the temperature, nanopore width, bulk pressure, and capillary condensation on the reaction equilibrium with respect to the reaction conversion, fluid structure, and self-diffusion coefficients. We show that the influence of the temperature, nanopore size, and capillary condensation on the confined reaction equilibrium is quite dramatic while the effect of the bulk pressure on the reaction equilibrium in the carbon slit nanopore is only moderate. This work is an extension of previous work by Turner et al. [J. Chem. Phys. 114, 1851 (2001)] on the confined reactive nitric oxide system.     

Excess free energy,enthalpy and entropy of surfactant-surfactant mixed micelle formation

Enthalpies of dilution and osmotic coefficients of sodium decylsulfate (NaDeS)-dodecyldimethylamine oxide (DDAO) mixtures in water were determined at 298 and 310 K, respectively. From the enthalpies of dilution, the apparent and then the partial molar relative enthalpies of the surfactant mixtures were calculated. From the osmotic coefficients, calculated at 298 K, the non-ideal free energies were derived. The latter were combined with the partial molar relative enthalpies to obtain the non-ideal entropies. From the apparent molar properties, using a previously reported approach, the excess thermodynamic properties for the surfactant-surfactant mixed micelle formation in water were evaluated as functions of the mixture composition at some total micellized concentration. In the whole range of the mixture composition, the excess free energy is negative, indicating that the mixed micelle formation is favoured with respect to that of pure micelles. This process is governed by the enthalpy and/or the entropy, depending on the mixture composition. The effect of the alkyl chain length was also studied by comparing the present results to those of the sodium dodecylsulfate-DDAO mixture.     

Structure, phase behavior, and inhomogeneous fluid properties of binary dendrimer mixtures

The effective pair potentials between different kinds of dendrimers in solution can be well approximated by appropriate Gaussian functions. We find that in binary dendrimer mixtures the range and strength of the effective interactions depend strongly upon the specific dendrimer architecture. We consider two different types of dendrimer mixtures, employing the Gaussian effective pair potentials, to determine the bulk fluid structure and phase behavior. Using a simple mean field density functional theory (DFT) we find good agreement between theory and simulation results for the bulk fluid structure. Depending on the mixture, we find bulk fluid-fluid phase separation (macrophase separation) or microphase separation, i.e., a transition to a state characterized by undamped periodic concentration fluctuations. We also determine the inhomogeneous fluid structure for confinement in spherical cavities. Again, we find good agreement between the DFT and simulation results. For the dendrimer mixture exhibiting microphase separation, we observe a rather striking pattern formation under confinement.     

Volumes and heat capacities of mixtures of dimethylformamide with several dialkylacetamides

Densities and heat capacities per unit volume of binary mixtures of dimethylformamide and a series of di-n-alkylacetamides have been measured and converted into excess molar volumes and heat capacities of the mixtures. In addition, the apparent and partial molar volumes and heat capacities of the various components have been evaluated. They vary smoothly with the mole fraction. The apparent molar heat capacities in the mixtures depend linearly on volume fraction, so that the partial molar heat capacities can be described using only one parameter for each mixture.     

Dynamics of fluid mixtures in nanospaces

A multicomponent extension of our recent theory of simple fluids [U. M. B. Marconi and S. Melchionna, J. Chem. Phys. 131, 014105 (2009)] is proposed to describe miscible and immiscible liquid mixtures under inhomogeneous, nonsteady conditions typical of confined fluid flows. We first derive from a microscopic level the evolution equations of the phase space distribution function of each component in terms of a set of self-consistent fields, representing both body forces and viscous forces (forces dependent on the density distributions in the fluid and on the velocity distributions). Second, we numerically solve the resulting governing equations by means of the lattice Boltzmann method, whose implementation contains novel features with respect to existing approaches. Our model incorporates hydrodynamic flow, diffusion, surface tension, and the possibility for global and local viscosity variations. We validate our model by studying the bulk viscosity dependence of the mixture on concentration, packing fraction, and size ratio. Finally, we consider inhomogeneous systems and study the dynamics of mixtures in slits of molecular thickness and relate structural and flow properties.     

Theory of gelation, vitrification, and activated barrier hopping in mixtures of hard and sticky spheres

Naive mode coupling theory (NMCT) and the nonlinear stochastic Langevin equation theory of activated dynamics have been generalized to mixtures of spherical particles. Two types of ideal nonergodicity transitions are predicted corresponding to localization of both, or only one, species. The NMCT transition signals a dynamical crossover to activated barrier hopping dynamics. For binary mixtures of equal diameter hard and attractive spheres, a mixture composition sensitive "glass-melting" type of phenomenon is predicted at high total packing fractions and weak attractions. As the total packing fraction decreases, a transition to partial localization occurs corresponding to the coexistence of a tightly localized sticky species in a gel-like state with a fluid of hard spheres. Complex behavior of the localization lengths and shear moduli exist because of the competition between excluded volume caging forces and attraction-induced physical bond formation between sticky particles. Beyond the NMCT transition, a two-dimensional nonequilibrium free energy surface emerges, which quantifies cooperative activated motions. The barrier locations and heights are sensitive to the relative amplitude of the cooperative displacements of the different species.