Asymptotics of pressure tensor in wedge-shaped film between solids characterized by dispersion forces

Asymptotic equations are derived for all four (three diagonal and one off-diagonal) components of the Irving-Kirkwood pressure tensor in a wedge-shaped film of a fluid located between solid phases characterized by dispersion forces. The equations describe the dependence of the aforementioned components on the position inside the film at large distances from the solid surfaces. It is shown that, when the walls of a slit are of different natures, the pressure tensor depends not only on the opening angle of the wedge-shaped slit, bus also on the position of the interface between the solids in contact outside of the slit. As a consequence, the symmetry of the spatial dependence of the tensor components with respect to the median plane is violated.     

Asymptotics of distribution functions in wedge-shaped film

Asymptotic formulas are derived for one- and two-particle distribution functions in a liquid wedge-shaped film. These formulas are valid at fairly large distances from the solid walls that confine the film. Along with the dependence on distances to solid walls, these formulas also determine the dependence on the opening angle of the wedge-shaped pore. In the limit of a straight angle, formulas are transformed into asymptotic formulas for the plane surface layer; in the limit of small angles, they are transformed into corresponding formulas for the plane film.     

Fluid in a closed narrow slit

The behavior of a fluid inside a closed narrow slit between solid walls is examined on the basis of the density functional theory. It is shown that the constraint of constant number of molecules leads to interesting effects which are absent when the slit is open and in contact with a reservoir. If the slit walls are identical, the density profiles at low temperatures or at high average densities rhoav of the fluid molecules in the slit have a sharp maximum in the middle of the slit, the value of the density at maximum being comparable to that of a liquid. The density of fluid at the walls is in this case comparable to the density of a vapor phase. At high temperatures or at low rhoav the fluid density in the middle of the slit is of the same order of magnitude as at the walls. For nonidentical walls the density maximum is shifted towards the wall with a stronger wall-fluid interaction. The transition between the two types (with and without the sharp maximum) of density profiles with the change of temperature in the slit occurs in a narrow range of temperatures, this range being larger for narrower slits. The pressures which the fluid exerts on the walls as well as the forces per unit area arising due to stresses in the sidewalls of the system can decrease with increasing rhoav. Such a behavior is not possible for homogeneous systems and can be explained by analyzing the fluid density at the walls when rhoav increases. The normal and transversal components of the pressure tensor were calculated as functions of the distance from the wall using the equation of hydrostatic equilibrium and direct calculation of the forces between molecules, respectively. The normal component decreases with increasing distance near the wall in contrast to the normal component near the liquid-vapor interface reported previously in the literature. The behavior of the transverse component does not depend on the fluid-solid interaction and is comparable to that for a liquid-vapor interface.     

Adsorption of a binary mixture of adhesive fluids in planar pores: a Monte Carlo study

Grand canonical Monte Carlo simulation is used to study the adsorption of a binary sticky hard-sphere fluid mixture in planar pores. The wall-component 1 and wall-component 2 contact densities are determined to calculate the pressure as a function of the composition of the mixture and the separation between the walls. From these data dependence of the solvation force between the plates on pore width is estimated. The simulation results are compared with the predictions of the Percus-Yevick approximation for planar pores. The density profiles of particular components show interesting shapes stemming from the interplay between the steric effects and the competitive adhesion among all possible species pairs. It is shown that narrowing of the pore causes selective partitioning of individual components of the mixture between the bulk phase and the interior of the pore. The agreement between the two methods is better at wider pores and for the component comprised of weakly adhesive particles.     

Computer simulation of Lennard-Jones fluid adsorption in wedge-shaped pores with smooth walls

The adsorption isotherms and local density distributions in wedge-shaped pores with smooth solid walls are obtained for a molecular fluid by the grand canonical Monte Carlo and multiple histogram reweighting methods. The effect of wedge opening angle on the phase behavior of the fluid is considered, and the peculiarities of the mechanism of wedge-shaped pore filling are established.     

An asymmetric slit: Pressure on larger wall

Local normal components of pressure tensor have been compared for the walls of a plane-parallel slit between a cylinder face and a planar surface. It has been shown that the normal pressure on the infinite surface is nonzero even outside of the slit. The correction for the finite size of the slit to the local normal pressure is more noticeable on this surface than on the surface of the cylinder.     

Grand canonical-like molecular dynamics simulations: application to anisotropic mass diffusion in a nanoporous medium

In this work, we describe two grand canonical-like molecular dynamics approaches to investigate mass diffusion phenomenon of a simple Lennard-Jones fluid confined between solid surfaces and in direct contact with reservoirs. In the first method, the density is used as the control variable in the reservoir whereas it is the pressure in the second method. Both methods provide consistent results, however, the constant density approach is the most efficient with respect to the computational time and implementation. Then, employing the constant density approach, we have studied the transient behavior of the diffusion process associated with the migration of one fluid into another one confined between parallel solid walls. Results have shown that the evolution of molar fraction of the invading fluid follows roughly a 1D diffusion model when the solid phase is weakly or moderately adsorbent with a characteristic time increasing when the pore width decreases. However, when the adsorption is high and the pore width small (i.e., below ten molecular sizes), the apparent mass diffusion in the adsorbed layer is reduced compared to that in the center of the slit pore. Hence, this mass diffusion process becomes a two-dimension phenomenon that must take into account an effective mass diffusion coefficient varying locally.     

Complex phase behavior of a fluid in slits with semipermeable walls modified with tethered chains

We study the phase behavior of a two-component fluid in a pore with the walls modified by tethered chains. The walls are completely permeable for one component of the fluid and completely impenetrable for the second component. The fluid is perfectly mixed in a bulk phase. We have found that depending on the details of the model the fluid undergoes capillary condensation inside the pore and wetting and layering transitions at the outer walls. Moreover, we have found transitions connected with the change of symmetry of the distribution of chains and fluid inside the pore.     

Pressure tensor in wedge-shaped cavity of solid with dispersion forces

The pressure tensor in a wedge-shaped cavity of a solid with dispersion forces is calculated. All four (three diagonal and one nondiagonal) components of the pressure tensor are calculated as functions of their positions inside the cavity. It is shown that these components satisfy the mechanical equilibrium condition. The dependence of the pressure tensor on the opening angle of the wedge-shaped cavity is determined. Disjoining pressure of vacuum film for its two alternative definitions is calculated.     

Approximation of contacts for calculating the velocity of the thermal motion of rodlike molecules in narrow slitlike pores

Expressions for calculating the thermal velocities of the motion of rodlike molecules in slitlike pores were obtained within the framework of a lattice gas model valid over a wide range of fluid densities (from rarefied gases to liquids) and temperatures, including the critical region. The translational and rotational motion of molecules was described using the transition-state theory for nonideal reaction systems, which takes into account the effect of the neighboring molecules on the activation barrier height. The local distributions of the components of the mixture under equilibrium conditions were calculated by describing lateral interactions in the approximation of isolated contacts. The equations of the model reflect the fact that the distributions of the components in the direction perpendicular to the pore walls (due to the effect of adsorption forces) and along the pore axis (if capillary condensation occurs) exhibit a strong anisotropy.     

On the calculation of pressure tensor in wedge-shaped cavity

A new simplified expression for the components of a pressure tensor in a wedge-shaped cavity of a solid with dispersion forces is derived. All four (three diagonal and one off-diagonal) components of the pressure tensor are calculated as functions of their positions inside the cavity. It is verified that these functions satisfy the conditions for mechanical equilibrium.     

Self-Diffusion Coefficients for Pure and Mixed Adsorbate Fluids in Narrow Pores

The equilibrium distribution and the concentration dependence of the local and average self-diffusion coefficients for pure fluid and binary mixture components in narrow slitlike pores were analyzed. The coefficients were calculated using the lattice gas model in the quasi-chemical approximation on the assumption of a spherical shape and approximately equal sizes of the components. For the pure adsorbate, these calculations were compared with molecular dynamics simulations. Both methods gave similar concentration profile changes and dynamic characteristics of interlayer particle redistributions in strong nonuniform adsorption fields for dense fluids. A satisfactory agreement was obtained for the temperature dependences of the self-diffusion coefficients along the pore axis. The influence of the molecule–wall potential and of intermolecular interaction were considered. The self-diffusion coefficients of the adsorbate were shown to strongly depend on the density of the mixture and the distance from pore walls.     

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.     

Molecular simulation of pressure-driven fluid flow in nanoporous membranes

An extended nonequilibrium molecular dynamics technique has been developed to investigate the transport properties of pressure-driven fluid flow in thin nanoporous membranes. Our simulation technique allows the simulation of the pressure-driven permeation of liquids through membranes while keeping a constant driving pressure using fluctuating walls. The flow of argon in the liquid state was simulated on applying an external pressure difference of 2.4x10(6) Pa through the slitlike and cylindrical pores. The volume flux and velocity distribution in the membrane pores were examined as a function of pore size, along with the interaction with the pore walls, and these were compared with values estimated using the Hagen-Poiseuille flow. The calculated velocity strongly depends on the strength of the interaction between the fluid and the atoms in the wall when the pore size is approximately<20sigma. The calculated volume flux also shows a dependence on the interaction between the fluid and the atoms in the wall. The Hagen-Poiseuille law overestimates or underestimates the flux depending on the interaction. From the analysis of calculated results, a good linear correlation between the density of the fluid in the membrane pores and the deviation of the flux estimated from the Hagen-Poiseuille flow was found. This suggests that the flux deviation in nanopore from the Hagen-Poiseuille flow can be predicted based on the fluid density in the pores.     

Phase behavior of a fluid confined in slitlike pores with walls modified by preadsorbed chain molecules

We have studied the microscopic structure, thermodynamics of adsorption, and phase behavior of Lennard-Jones fluid in slitlike pores with walls modified due to preadsorption of chain molecules. The chain species are grafted at the walls by terminating segments. Our theoretical considerations are based on a density functional approach in the semigrand canonical ensemble. The applied constraint refers to the constant number of grafted chain molecules in the pore without restriction of the number of chains at each of the walls. We have observed capillary condensation of Lennard-Jones fluid combined with the change of the distribution of chains from nonsymmetric to symmetric with respect to the pore walls. The phase diagrams of the model are analyzed in detail, dependent on the pore width, length of chains, and grafted density.     

Asymptotics of pressure tensor in the wedge-shaped film of simple fluid: 1. Calculation of Irving-Kirkwood tensor

The spatial behavior of the components of a pressure tensor inside the filled wedge-shaped cavity of a solid, which was calculated in a previous communication, is analyzed. The approximation for introducing and studying the disjoining pressure in a wedge-shaped film is developed.     

Asymptotics of pressure tensor in wedge-shaped film of simple fluid: II. Spatial dependence of pressure tensor and disjoining pressure

The spatial behavior of the components of a pressure tensor inside the filled wedge-shaped cavity of a solid, which was calculated in a previous communication, is analyzed. The approximation for introducing and studying the disjoining pressure in a wedge-shaped film is developed.     

Pressure driven flow of polymer solutions in nanoscale slit pores

Polymer solutions subject to pressure driven flow and in nanoscale slit pores are systematically investigated using the dissipative particle dynamics approach. The authors investigated the effect of molecular weight, polymer concentration, and flow rate on the profiles across the channel of the fluid and polymer velocities, polymer density, and the three components of the polymers radius of gyration. They found that the mean streaming fluid velocity decreases as the polymer molecular weight and/or polymer concentration is increased, and that the deviation of the velocity profile from the parabolic profile is accentuated with increase in polymer molecular weight or concentration. They also found that the distribution of polymers conformation is highly anisotropic and nonuniform across the channel. The polymer density profile is also found to be nonuniform, exhibiting a local minimum in the center plane followed by two symmetric peaks. They found a migration of the polymer chains either from or toward the walls. For relatively long chains, as compared to the thickness of the slit, a migration toward the walls is observed. However, for relatively short chains, a migration away from the walls is observed.     

Diffusion and velocity relaxation of a Brownian particle immersed in a viscous compressible fluid confined between two parallel plane walls

The diffusion tensor and velocity correlation function of a Brownian particle immersed in a viscous compressible fluid confined between two parallel plane walls are calculated in point approximation. The fluid is assumed to satisfy stick boundary conditions at the walls. It is found that the velocity correlation function decays asymptotically with a negative t(-2) long-time tail due to coupling to overdamped sound waves. The coefficient of the long-time tail is calculated and shown to be independent of fluid viscosity.