Local linear viscoelasticity of confined fluids

In this paper the authors propose a novel method to study the local linear viscoelasticity of fluids confined between two walls. The method is based on the linear constitutive equation and provides details about the real and imaginary parts of the local complex viscosity. They apply the method to a simple atomic fluid undergoing zero mean oscillatory flow using nonequilibrium molecular dynamics simulations. The method shows that the viscoelastic properties of the fluid exhibit dramatic spatial changes near the wall-fluid boundary due to the high density in this region. It is also shown that the real part of the viscosity converges to the frequency dependent local shear viscosity sufficiently far away from the wall. This also provides valuable information about the transport properties in the fluid, in general. The viscosity is compared with predictions from the local average density model. The two methods disagree in that the local average density model predicts larger viscosity variations near the wall-fluid boundary than what is observed through the method presented here.     

The structure of fluids confined in crystalline slitlike nanoscopic pores

Grand canonical and canonical ensemble Monte Carlo simulation methods are used to study the structure and phase behavior of Lennard-Jones fluids confined between the parallel (100) planes of the face centered cubic crystal. Thin slit pores with a width allowing for the formation of only up to five atomic layers are considered. The phase diagrams of the systems characterized by different pore width as well as by different strength of the fluid-pore walls potential are determined. It is shown that an enormously large number of different phase diagram topologies can occur, depending on the parameters of the problem (pore width, strength of the fluid-pore walls potential, etc).     

Structure of Lennard-Jones fluids confined in square nanoscale channels from density functional theory

The density distribution of Lennard-Jones fluids confined in square nanoscale channels with Lennard-Jones walls has been studied using the nonlocal density functional theory (DFT) based on the Tarazona model. The effect of channel lengths on the density profiles with various chemical potentials was discussed. It was found that there is an apparent layering phenomenon for the confined fluids due to the combining influences of the enhancing solid-fluid interaction and the excluded volume effect. The pronounced density peaks were observed at the corners of square channels due to the strong fluid-solid interactions. The grand canonical ensemble Monte Carlo simulation (GCEMC) was applied to test the nonlocal DFT results. The DFT calculations are in relatively good agreement with the GCEMC simulations. The adsorption isotherms in a series of square channels were evaluated based on the obtained density distributions. The adsorption mechanism within the square pores was investigated. A comparison between the adsorptions of the square pores with those of the corresponding slit-size pores has been given.     

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.     

The structure of fluids confined in crystalline slitlike nanoscopic pores: bilayers

Grand canonical and canonical ensemble Monte Carlo simulation methods are used to study the structure and phase behavior of Lennard-Jones fluids confined between the parallel (100) planes of the face centered cubic crystal. Ultra thin slit pores of the width allowing for the formation of only two adsorbate layers are considered. It is demonstrated that the structure of adsorbed phases is very sensitive to the wall-wall separation and to the strength of the fluid-wall potential. It is also shown that the structure of low temperature (solid) phases strongly depends on the fluid density. In particular, when the surface field is sufficiently strong, then the high density phases may exhibit a domain wall structure, quite the same as found in monolayer films adsorbed at a single substrate wall. On the other hand, the weakening of the surface potential leads to the regime in which only the hexagonally ordered bilayer structure is stable. The phase diagrams for a series of systems are estimated. It is shown that, depending on the pore width and the temperature, the condensation leads to the formation of the commensurate or incommensurate phases. The incommensurate phases may have the domain-wall or the hexagonal structure depending on the pore width and the strength of the fluid-wall potential.     

Using gamma distributions to predict self-diffusivities and density of states of fluids confined in carbon nanotubes

The density of states of rare gas atoms confined in carbon nanotubes is analyzed using a recently proposed model based on gamma distributions [Krishnan and Ayappa, J. Chem. Phys., 124 144503 (2006)]. The inputs into the model are the 2nd and 4th frequency moments that are obtained from molecular dynamics simulations. The predicted density of states, velocity autocorrelation functions and self-diffusivities are compared with those obtained from molecular dynamics simulations, for different nanotube loadings and temperatures. All results are reported for argon confined in a (16,16) carbon nanotube. The model predictions are extremely accurate at intermediate reduced densities of rhosigma(3) = 0.3, 0.4, where the majority of the self-diffusivity predictions lie within 10% of the simulation results. Since the frequency moments can be also obtained from Monte Carlo simulations, the study suggests an alternate route to the system dynamics of strongly confined fluids.     

Molecular simulation of condensation process of Lennard-Jones fluids confined in nanospace with jungle-gym structure

We report grand canonical Monte Carlo simulations for a Lennard-Jones (LJ) fluid modeled on methane confined in nanospace with jungle-gym-like (JG) cubic structure, which is typically found in porous coordination polymers. Pillars composing the cubic structure were modeled as structureless smooth solid rods made of LJ carbon. We examined the effects of pore size, pore geometry, rod thickness, and rod potential onto the condensation phenomena in the JG pore structure. The simulations clarified that the condensation pressure and adsorption amount in the JG structure were influenced by pore size and rod potential, while the transition type was determined by rod thickness. The characteristics of the JG structure lie in the sensitivity to the slight changes in pore size, rod thickness, and rod potential owing to the combination of the packing effect of molecules and the superposition effect of rod potentials.     

Viscosity of confined inhomogeneous nonequilibrium fluids

We use the nonlocal linear hydrodynamic constitutive model, proposed by Evans and Morriss [Statistical Mechanics of Nonequilibrium Liquids (Academic, London, 1990)], for computing an effective spatially dependent shear viscosity of inhomogeneous nonequilibrium fluids. The model is applied to a simple atomic fluid undergoing planar Poiseuille flow in a confined channel of several atomic diameters width. We compare the spatially dependent viscosity with a local generalization of Newton's law of viscosity and the Navier-Stokes viscosity, both of which are known to suffer extreme inaccuracies for highly inhomogeneous systems. The nonlocal constitutive model calculates effective position dependent viscosities that are free from the notorious singularities experienced by applying the commonly used local constitutive model. It is simple, general, and has widespread applicability in nanofluidics where experimental measurement of position dependent transport coefficients is currently inaccessible. In principle the method can be used to predict approximate flow profiles of any arbitrary inhomogeneous system. We demonstrate this by predicting the flow profile for a simple fluid undergoing planar Couette flow in a confined channel of several atomic diameters width.     

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.     

Flow and anchoring effects on nematic fluctuations in confined geometry

Light scattering observation of nematic director fluctuations in confined geometries can be used to obtain interaction parameters of liquid crystals with surfaces. We present the basics of the method and some examples of the results in planar and cylindrical geometries. These results were obtained after neglecting the coupling of the director motion to flow. We give analytical and numerical results of flow effects on director fluctuations in a slab. The backflow contribution to the effective viscosity is strongly suppressed so that the results for the anchoring energy remain valid. Modal dispersion relations show an interesting behaviour of avoiding crossings.     

Thermodynamic and structural properties of confined discrete-potential fluids

The thermodynamic and structural behaviors of confined discrete-potential fluids are analyzed by computer simulations, studying in a systematic way the effects observed by varying the density, temperature, and parameters of the potentials that characterize the molecule-molecule interactions. The Gibbs ensemble simulation technique for confined fluids [A. Z. Panagiotopoulos, Mol. Phys. 62, 701 (1987)] is applied to a fluid confined between two parallel hard walls. Two different systems have been considered, both formed by spherical particles that differ by the interparticle pair potential: a square well plus square shoulder or a square shoulder plus square well interaction. These model interactions can describe in an effective way pair potentials of real molecular and colloidal systems. Results are compared with the simpler reference systems of square-shoulder and square-well fluids, both under confinement. From the adsorption characterization through the use of density profiles, it is possible to obtain specific values of the interparticle potential parameters that result in a positive to negative adsorption transition.     

Local fluctuations in solution mixtures

An extension of the traditional Kirkwood-Buff (KB) theory of solutions is outlined which provides additional fluctuating quantities that can be used to characterize and probe the behavior of solution mixtures. Particle-energy and energy-energy fluctuations for local regions of any multicomponent solution are expressed in terms of experimentally obtainable quantities, thereby supplementing the usual particle-particle fluctuations provided by the established KB inversion approach. The expressions are then used to analyze experimental data for pure water over a range of temperatures and pressures, a variety of pure liquids, and three binary solution mixtures - methanol and water, benzene and methanol, and aqueous sodium chloride. In addition to providing information on local properties of solutions it is argued that the particle-energy and energy-energy fluctuations can also be used to test and refine solute and solvent force fields for use in computer simulation studies.     

Freezing of fluids confined between mica surfaces

Using grand ensemble simulations, we show that octamethyl-cyclo-tetra-siloxane (OMCTS) confined between two mica surfaces can form a variety of frozen phases which undergo solid-solid transitions as a function of the separation between the surfaces. For atomically smooth mica surfaces, the following sequence of transitions 1[triangle up] --> 1[triangle up]b --> 2B --> 2 square --> 2[triangle up] are observed in the one- and two-layered regimes, where n[triangle up], n[square], and nB denote triangular, square, and buckled phases, respectively, with the prefix n denoting the number of confined layers. The presence of potassium on mica is seen to have a strong influence on the degree of order induced in the fluid. The sequence of solid-solid transitions that occurs with the smooth mica surface is no longer observed. When equilibrated with a state point near the liquid-solid transition, a counterintuitive freezing scenario is observed in the presence of potassium. Potassium disrupts in-plane ordering in the fluid in contact with the mica surface, and freezing is observed only in the inner confined layers. The largest mica separations at which frozen phases were observed ranged from separations that could accommodate six to seven fluid layers. The extent of freezing and the square-to-triangular lattice transition was found to be sensitive to the presence of potassium as well as the thermodynamic conditions of the bulk fluid. The implications of our results on interpretation of surface force experiments as well as the generic phase behavior of confined soft spheres is discussed.     

Dielectric permittivity profiles of confined polar fluids

The dielectric response of a simple model of a polar fluid near neutral interfaces is examined by a combination of linear response theory and extensive molecular dynamics simulations. Fluctuation expressions for a local permittivity tensor epsilon(r) are derived for planar and spherical geometries, based on the assumption of a purely local relationship between polarization and electric field. While the longitudinal component of epsilon exhibits strong oscillations on the molecular scale near interfaces, the transverse component becomes ill defined and unphysical, indicating nonlocality in the dielectric response. Both components go over to the correct bulk permittivity beyond a few molecular diameters. Upon approaching interfaces from the bulk, the permittivity tends to increase, rather than decrease as commonly assumed, and this behavior is confirmed for a simple model of water near a hydrophobic surface. An unexpected finding of the present analysis is the formation of "electrostatic double layers" signaled by a dramatic overscreening of an externally applied field inside the polar fluid close to an interface. The local electric field is of opposite sign to the external field and of significantly larger amplitude within the first layer of polar molecules.     

A diffusion model for the fluids confined in micropores

The self-diffusion coefficients were calculated by molecular dynamics simulations and the effects of pore width, temperature, and fluid density on diffusion behavior of simple fluid argon and polar fluid water confined in micropores were analyzed and studied. A mathematical model describing diffusion behavior of fluids confined in micropores was proposed from the theories of molecular dynamics and molecular kinematics, and validated on the basis of the simulation results at various conditions. The model indicates that the diffusion coefficient is proportional to the square root of the pore width and to the temperature divided by the density squared. It is applicable to either liquid or gas states and only two parameters are required.     

Probing memory effects in confined fluids via diffusion measurements

Confinement of fluids in porous materials is widely exploited in a variety of technologies, including chemical conversion by heterogeneous catalysis and adsorption separations. Important fundamental phenomena associated with many-molecule interactions occur in such systems, including a remarkably long "memory" of the past when the actual amount of molecules in the pores dramatically depends on the history of how the external conditions have been changed. We demonstrate that the intrinsic diffusivity as measured by NMR serves as an excellent probe of the history-dependent states of the confined fluid. A remarkable feature of our results are differences in diffusivity between out-of-equilibrium states with the same density within the hysteresis loop. This reflects different spatial distributions of the confined fluid that accompany the arrested equilibration of the system in this region.     

Distribution of transverse chain fluctuations in harmonically confined semiflexible polymers

Two different experimental studies of polymer dynamics based on single-molecule fluorescence imaging have recently found evidence of heterogeneities in the widths of the putative tubes that surround filaments of F-actin during their motion in concentrated solution. In one [J. Glaser, D. Chakraborty, K. Kroy, I. Lauter, M. Degawa, N. Kirchesner, B. Hoffmann, R. Merkel, and M. Giesen, Phys. Rev. Lett. 105, 037801 (2010)], the observations were explained in terms of the statistics of a worm-like chain confined to a potential determined self-consistently by a binary collision approximation, and in the other [B. Wang, J. Guan, S. M. Anthony, S. C. Bae, K. S. Schweizer, and S. Granick, Phys. Rev. Lett. 104, 118301 (2010)], they were explained in terms of the scaling properties of a random fluid of thin rods. In this paper, we show, using an exact path integral calculation, that the distribution of the length-averaged transverse fluctuations of a harmonically confined weakly bendable rod (one possible realization of a semiflexible chain in a tube), is in good qualitative agreement with the experimental data, although it is qualitatively different in analytic structure from the earlier theoretical predictions. We also show that similar path integral techniques can be used to obtain an exact expression for the time correlation function of fluctuations in the tube cross section.     

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.     

Concentration fluctuations in binary mixtures of model polar fluids

We present calculations of the mean square concentration fluctuations, S_α(0), for binary mixtures of model polar fluids. The assumed pair interactions are taken to be of the forms of a hard core plus either dipole-dipole, quadrupole-quadrupole or dipole- quadrupole terms. The calculations are carried out within a mean field approximation. We have considered in some detail the interplay between size differences and the difference in the strength and range of the potentials in deciding how S_α(0) deviates from the ideal behaviour.