Binary fluids in planar nanopores: Adsorptive selectivity, heat capacity and self-diffusivity

Monte Carlo and molecular dynamics simulations are performed to study fluid adsorption of a two component fluid in slit pores of nanoscopic dimensions. The slit pores are immersed in a binary fluid bath, which is comprised of spherical molecules having a size ratio of 1.43, at constant temperature and composition. Pore width is varied to determine how the heat capacity and self-diffusion coefficient are linked to the composition and structure of the adsorbed fluid. In pores where the fluid structure is most pronounced, we observe: perfect (or near perfect) exclusion of one component by the other component, a heat capacity that rapidly oscillates and is of greater magnitude than in the fluid bath, and self-diffusion coefficients on the order of 10^–8 cm²/s. The behavior of the heat capacity and diffusion coefficients appears to arise from a near solid-like layering of OMCTS that occurs at certain favorable pore widths.     

Surface excess free energy of simple fluids confined in cylindrical pores by isothermal-isobaric Monte Carlo: influence of pore size

Confined fluid properties are mainly determined by interfacial phenomena characterized by surface quantities. Based on a simple model of Lennard-Jones particles confined in a cylindrical pore, this study introduces a grand potential surface quantity to quantify the difference in the thermodynamic pressure between the bulk and the confined fluids. The usual surface tension gamma defined as this grand potential difference for the same chemical potential in both confined and bulk states is generally strongly dependent on both the chemical potential and temperature. It is proposed here to introduce another surface quantity zeta which measures the thermodynamic pressure difference between confined and bulk states for identical densities. It is shown that this quantity is much less dependent on confined fluid density or chemical potential. It is actually constant along the gas-like and liquid-like adsorption/desorption branches for an irreversible isotherm (hysteresis), with a different value for each branch. For reversible supercritical isotherms, zeta is shown to remain constant in the low and high density parts of the isotherm. This independence on chemical potential (or equivalently fluid density) is believed to be of great interest for practical applications when one desires to calculate thermodynamic quantities such as the usual surface tension gamma or the thermodynamic pressure of a confined fluid for any given chemical potential and temperature. Such calculations are required to determine fundamental properties such as metastability or coexistence. The effects of temperature, fluid/substrate interaction strength, and pore size are studied.     

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.     

Torsion-induced phase transitions in fluids confined between chemically decorated substrates

In this paper we investigate the phase behavior of a "simple" fluid confined to a chemically heterogeneous slit pore of nanoscopic width s(z) by means of Monte Carlo simulations in the grand canonical ensemble. The fluid-substrate interaction is purely repulsive except for elliptic regions of semiaxes A and B attracting fluid molecules. On account of the interplay between confinement (i.e., s(z)) and chemical decoration, three fluid phases are thermodynamically permissible, namely, gaslike and liquidlike phases and a "bridge phase" where the molecules are preferentially adsorbed by the attractive elliptic patterns and span the gap between the opposite substrate surfaces. Because of their lack of cylindrical symmetry, bridge phases can be exposed to a torsional strain 0相似文献   

Relationships between three-body and two-body interactions in fluids and solids

Molecular dynamics data are reported for two-body and three-body interactions in noble gases at densities covering the gas, liquid, and solid phases. The data indicate that simple relationships exist between three- and two-body interactions in both fluid and solid phases. The relationship for liquids has a simple density dependence with only one external parameter. In contrast, the solid phase relationship depends both on density and on the square of density and requires the evaluation of two parameters. The relationships are tested for both system-size and temperature dependences. The values of the relationship parameters are only sensitive to system size when a small number of atoms are involved. For 500 or more atoms, they remain nearly constant. The relationships are valid for both subcritical and slightly supercritical temperatures. A practical benefit of the relationships is that they enable the use of two-body intermolecular potentials for the prediction of the properties of real systems without the computational expense of three-body calculations.     

Pseudocritical or hysteresis temperature versus pore size for simple fluids confined in cylindrical nanopores

The adsorption/desorption isotherms measured in nanoporous materials generally present a hysteresis. The hysteresis shrinks upon increasing the temperature (for a given pore size) or decreasing the pore size (for a given temperature), until it finally disappears at the so-called hysteresis (or pseudocritical) temperature T(h) or hysteresis (or pseudocritical) pore size R(h), not to be confused with a true critical point. In this paper, a Monte Carlo approach allowed calculating the surface free energy of confined fluid along the adsorption/desorption isotherms for various cylindrical pore sizes and temperatures. A simple phenomenological model then allowed exploiting these results to determine the relation between T(h) and R(h). The prediction is compared to various literature models and experimental data, showing agreement within uncertainties. On the other hand, the simulations cannot be used directly to predict T(h) and R(h) since they significantly overestimate the hysteresis width. The model predicts a nonlinear relation between the reduced hysteresis temperature and the inverse pore radius.     

Comparing landscape calculations with calorimetric data on ortho-terphenyl, and the question of the configurational fraction of the excess entropy

Mossa et al. [Phys. Rev. E 65, 041205 (2002)] have calculated the total and configurational entropies of supercooled ortho-terphenyl liquid using the potential-energy landscape formalism and a simplified model of the intermolecular potential. I show here that the agreement of their calculated configurational entropy with the experimental data depends on what is assumed about the configurational fraction of the excess entropy and its temperature dependence. In particular, if the configurational fraction is taken as 0.70 and independent of temperature the agreement is excellent; if a marked temperature dependence of that fraction inferred from calorimetric data is assumed the agreement is only fair at best. This marked temperature dependence of the configurational fraction also implies some implausible behavior of contributions to the excess entropy at the Kauzmann temperature, but no obvious reason for disregarding it presents itself.     

Fluidic assembly and packing of microspheres in confined channels

We study fluidic assembly and packing of spherical particles in rectilinear microchannels that are terminated by a flow constriction. First, we introduce a method for active assembly of particles in the confined microchannels by triggering a local constriction in the fluid channel using a partially closed membrane valve. This microfluidic valve allows active, on-demand particle assembly as opposed to previous passive assembly methods based on terminal channels and weirs. Second, we study the three-dimensional assembly and packing of particles against a weir in confined rectilinear microchannels. The packings result in achiral particle chains with alternating (zigzag) structure. This structure is characterized by a single, repeated bond angle whose components projected into the frame of the channel are quantified by confocal microscopy and image processing. Brownian dynamics simulation of the packing comprehensively delineates the range of bond angles possible in narrow, rectilinear microchannels as well as the complex dependence of these angles on the relative dimensions of the channel and particles. The simulations of the three-dimensional packings are accurately modeled by a compact theory based on trigonometric relationships. The experimentally measured bond angles show excellent agreement with the simulations, thereby validating the functional dependence of the achiral packing bond angles on channel dimensions. This functional relationship is immediately useful for the design of anisotropic particles by microfluidic synthesis.     

Cooperative effects, transport and entropy in simple liquids

We systematically investigate the cooperative effects in shear stress relaxation using equilibrium molecular-dynamics simulations in periodic boundary conditions containing a variable degree of strain. We show that, even in simple liquids, shear stress relaxation is a cooperative effect associated with a correlation length that increases with isobaric decrease in temperature. If the system size is less than the correlation length, shear stress in the system is determined by the boundary strain. Transport, however, does not depend on the boundary conditions. We relate these two effects to the number and properties of the configurations accessible to the system.     

Local structure and density fluctuations in confined fluids

Spatial confinement modifies the microscopic structure of dense fluids, thereby inducing for example structural forces between the confining walls. However, confinement also modifies the fluids' density fluctuations, resulting in more elusive but equally important effects. In this brief review it is shown that both of these phenomena are naturally analyzed using the confined fluid's pair densities, which have recently become also experimentally accessible. Two particular topics are discussed, namely, the mechanisms of oscillatory density profiles and ensuing solvation forces in dense confined fluids as well as the behavior of liquids in solvophobic confinement.     

Diffusivity, excess entropy, and the potential-energy landscape of monatomic liquids

The connection between thermodynamic, transport, and potential-energy landscape features is studied for liquids with Lennard-Jones-type pair interactions using both microcanonical molecular-dynamics and isothermal-isobaric ensemble Monte Carlo simulations. Instantaneous normal-mode and saddle-point analyses of two variants of the monatomic Lennard-Jones liquid have been performed. The diffusivity is shown to depend linearly on several key properties of instantaneous and saddle configurations-the energy, the fraction of negative curvature directions, and the mean, maximum, and minimum eigenvalues of the Hessian. Since the Dzugutov scaling relationship also holds for such systems [Nature (London) 381, 137 (1996)], the exponential of the excess entropy, within the two-particle approximation, displays the same linear dependence on energy landscape properties as the diffusivity.     

Thermodynamic, dynamic, structural, and excess entropy anomalies for core-softened potentials

Using molecular dynamic simulations, we study three families of continuous core-softened potentials consisting of two length scales: a shoulder scale and an attractive scale. All the families have the same slope between the two length scales but exhibit different potential energy gap between them. For each family three shoulder depths are analyzed. We show that all these systems exhibit a liquid-liquid phase transition between a high density liquid phase and a low density liquid phase ending at a critical point. The critical temperature is the same for all cases suggesting that the critical temperature is only dependent on the slope between the two scales. The critical pressure decreases with the decrease of the potential energy gap between the two scales suggesting that the pressure is responsible for forming the high density liquid. We also show, using the radial distribution function and the excess entropy analysis, that the density, the diffusion, and the structural anomalies are present if particles move from the attractive scale to the shoulder scale with the increase of the temperature indicating that the anomalous behavior depends only in what happens up to the second coordination shell.     

Longitudinal and volume viscosities of fluids confined in nanochannel

Longitudinal and volume viscosities of Lennard-Jones fluid, argon–krypton binary mixture and isotopic fluid mixture confined to nanochannels of different widths are calculated by employing theoretical technique based on Green–Kubo formula. A significant enhancement is observed in longitudinal and volume viscosities when width of the nanochannel is less than 10 nm. Effect of mass ratio of two species on longitudinal and volume viscosities is also studied for equimolar isotopic fluid mixture. It is found that enhancement in viscosity is more for larger mass ratios. It is also noted that enhancement in longitudinal and shear viscosities is more than volume viscosity.     

Solvation force, structure and thermodynamics of fluids confined in geometrically rough pores

The effect of periodic surface roughness on the behavior of confined soft sphere fluids is investigated using grand canonical Monte Carlo simulations. Rough pores are constructed by taking the prototypical slit-shaped pore and introducing unidirectional sinusoidal undulations on one wall. For the above geometry our study reveals that the solvation force response can be phase shifted in a controlled manner by varying the amplitude of roughness. At a fixed amplitude of roughness, a, the solvation force for pores with structured walls was relatively insensitive to the wavelength of the undulation, lambda for 2.3/=0.5. The predictions of the superposition approximation, where the solvation force response for the rough pores is deduced from the solvation force response of the slit-shaped pores, was in excellent agreement with simulation results for the structured pores and for lambda/sigma(ff)>/=7 in the case of smooth walled pores. Grand potential computations illustrate that interactions between the walls of the pore can alter the pore width corresponding to the thermodynamically stable state, with wall-wall interactions playing an important role at smaller pore widths and higher amplitudes of roughness.     

Molecular interaction of organic dyes in bulk and confined media

Molecular interactions of five thiazine dyes with increasing alkyl substitution have been studied in aqueous and microemulsion media at 303 K within a concentration range of (1.35–7.00) × 10^?4 M. The dimerization constant (K_d) values for the five dyes are ranged between 1.761 and 6.258 × 10³ l mol^?1 in bulk water media, where as in microemulsion media, K_d's are ranged between 1.760 and 4.110 × 10³ l mol^?1. Thionine (with no methyl substitution) and azure A (with two methyl substitution) displayed slightly larger K_d values in microemulsion water pools compared to bulk water while other dyes recorded significant drop in K_d values. The influence of microemulsion media on the molecular interaction of dyes has been explained in terms of electrostatic and hydrophobic factors. The monomer and the dimer spectra are explained in terms of molecular exciton model and the optical absorption parameters of both the species are reported in bulk and confined media.     

Waterlike structural and excess entropy anomalies in liquid beryllium fluoride

The relationship between structural order metrics and the excess entropy is studied using the transferable rigid ion model (TRIM) of beryllium fluoride melt, which is known to display waterlike thermodynamic anomalies. The order map for liquid BeF2, plotted between translational and tetrahedral order metrics, shows a structurally anomalous regime, similar to that seen in water and silica melt, corresponding to a band of state points for which average tetrahedral (q(tet)) and translational (tau) order are strongly correlated. The tetrahedral order parameter distributions further substantiate the analogous structural properties of BeF2, SiO2, and H2O. A region of excess entropy anomaly can be defined within which the pair correlation contribution to the excess entropy (S2) shows an anomalous rise with isothermal compression. Within this region of anomalous entropy behavior, q(tet) and S2 display a strong negative correlation, indicating the connection between the thermodynamic and the structural anomalies. The existence of this region of excess entropy anomaly must play an important role in determining the existence of diffusional and mobility anomalies, given the excess entropy scaling of transport properties observed in many liquids.     

An interpretative model for the anomalous behavior of some excess properties in mixed liquid systems: a relationship between excess molar volumes and excess compressibilities in strongly self-aggregated fluids

A simple theoretical model based on the hole theory of the liquid state has been developed in order to find a relationship between the excess molar volume and the excess compressibility of a two-component fluid. The model has been extended to investigate strongly associating fluids dissolved in an apolar solvent. Such a model nicely explains our recent results obtained from a Brillouin scattering study of the nonideal mixture between tert-butyl alcohol and 2-2(') dimethyl butane. Its validity, however, seems to be rather general and it could be usefully applied to rationalize the excess properties of several nonideal binary fluid mixtures.     

Phase behavior of magnetic nanoparticles dispersions in bulk and confined geometries

Dispersions of magnetic nanoparticles (ferrofluids) have been for a long time taken as examples of dipolar fluids. Theoretical papers conclude now that an important parameter is the ratio of the anisotropic attractions (dipolar ones) to the isotropic attractions (Van der Waals ones). It is confirmed by the recent experimental results concerning the behavior of small magnetic particles in bulk or confined 2D geometries.