Modeling molecular transport in slit pores

We examine the transport of methane in microporous carbon by performing equilibrium and nonequilibrium molecular dynamics simulations over a range of pore sizes, densities, and temperatures. We interpret these simulation results using two models of the transport process. At low densities, we consider a molecular flow model, in which intermolecular interactions are neglected, and find excellent agreement between transport diffusion coefficients determined from simulation, and those predicted by the model. Simulation results indicate that the model can be applied up to fluid densities of the order to 0.1-1 nm(-3). Above these densities, we consider a slip flow model, combining hydrodynamic theory with a slip condition at the solid-fluid interface. As the diffusion coefficient at low densities can be accurately determined by the molecular flow model, we also consider a model where the slip condition is supplied by the molecular flow model. We find that both density-dependent models provide a useful means of estimating the transport coefficient that compares well with simulation.     

Adsorption of water molecules in slit pores

We report molecular dynamics (MD) simulations on the adsorption of water in attractive and repulsive slit pores, where the slit and a bulk region are in contact with each other. Water structure, surface force and adsorption behavior are investigated as a function of the overall density in the bulk region. The gas–liquid transition in both types of pores occurs at similar densities of the bulk region.     

Anisotropic dynamics of dipolar liquids in narrow slit pores

We report molecular dynamics simulation results for Stockmayer fluids confined to narrow slitlike pores with structureless, nonconducting walls. The translational and rotational dynamics of the dipolar particles have been investigated by calculating autocorrelation functions, diffusion coefficients, and relaxation times for various pore widths (five or less particle diameters) and directions parallel and perpendicular to the walls. The dynamic properties of the confined systems are compared to bulk properties, where corresponding bulk and pore states at the same temperature and chemical potential are determined in parallel grand canonical Monte Carlo simulations. We find that the dynamic behavior inside the pore depends on the distance from the walls and can be strongly anisotropic even in globally isotropic systems. This concerns especially the particles in the surface layers close to the walls, where the single particle and collective dipolar relaxation resemble that of true two-dimensional dipolar fluids with different in-plane and out-of-plane relaxations. On the other hand, bulklike relaxation is observed in the pore center of sufficiently wide pores.     

Dielectric response of polar liquids in narrow slit pores

Based on molecular dynamics (MD) simulations and a simple (Stockmayer) model we investigate the static and dynamic dielectric response of polar liquids confined to narrow slit pores. The MD simulations are used to calculate the time-dependent polarization fluctuations along directions parallel and perpendicular to the walls, from which the components of the frequency-dependent dielectric tensor can be derived via linear response theory. Our numerical results reveal that the system's response is strongly anisotropic. The parallel dielectric function, epsilonparallel(omega), has Debye-like character very similar to the corresponding isotropic bulk function, epsilonbulk(omega), at the same chemical potential. Indeed, the main confinement effect on epsilonparallel(omega) consists in a shift toward smaller values relative to the bulk function. On the other hand, in the perpendicular direction we observe a characteristic peak in the absorption part of the dielectric function, epsilonperpendicular(omega). This peak is absent in the bulk system and reflects strongly pronounced, damped oscillations in the polarization fluctuations normal to the walls. We discuss two possible origins of the oscillations (and the resulting absorption peak), that is collective oscillations of dipoles in clusters formed parallel to the walls, and the existence of a "dipolaron mode" previously observed in MD simulations of bulk polar fluids.     

Branched versus linear alkane adsorption in carbonaceous slit pores

The presence of carbonaceous deposits on the internal surfaces of a spark ignition engine has been linked in the literature to impaired vehicle performance, as manifested by increased knocking, higher fuel consumption, higher emissions and other adverse effects. One of the proposed mechanisms, in which the deposits affect the processes in the engine, is the adsorption and desorption of fuel components in the pores within the deposit. In this article we investigate this mechanism in more detail by considering single component adsorption of normal and branched alkanes in selected model slit pores representing the structure of the deposits. We further extend these studies to the binary mixture of normal heptane and isooctane, corresponding to a primary reference fuel blend. In particular, we show that in larger pores adsorption selectivity towards isooctane is about 1.2 on average throughout the pressure range. However, in the smaller 10 Å pore selectivity towards isooctane can be in excess of three as a result of packing effects. These results are then placed in the context of engine performance issues.     

Polymolecular adsorption and capillary condensation in narrow slit pores

Capillary condensation and polymolecular adsorption in narrow slits has been calculated, where the fields of surface forces overlap one another. The calculations were carried out on the basis of macroscopic theory of dispersion forces and the isotherms of lone adsorption layers at the free surface. It has been shown that under the effect of mutual attraction through a gap, polymolecular adsorption films lose their stability long before their thickness has approached the half-width of a flat slit. This results in hysteresis of the capillary condensation in an ensemble of plane-parallel slits.

In the case of systems having strong adsorbate-adsorbate interaction, there has been detected the existence of the lower limit of sizes of slit pores, wherein the capillary meniscus can coexist with adsorption films. With a slit width smaller than the critical one, the meniscus is likely to form a finite contact angle with “dry” surfaces of a slit. Thus an explanation has been given of the lower limit of the capillary condensation in an ensemble of flat-surface, slit pores. In the case of strong adsorbate-adsorbent interaction, the coexistence of meniscus with adsorption films within the scope of the approach used is possible in slits of any width.

The value of corrections for the surface forces effect to be entered in the calculations of slit pores dimensions has been analyzed on the basis of the capillary condensation data obtained.

In wedge-shaped slits there also exists, besides lower limit the upper limit of capillary hysteresis.     

Orientational orders of small anisotropic molecules confined in slit pores

Based on a constant-pressure Monte Carlo molecular simulation, we have studied orientationally ordered transitions of small anisotropic molecules confined in two parallel hard walls. These molecules are modeled by the hard Gaussian overlap model. The molecular elongations of the chosen molecules are so small that the molecules cannot form stable liquid-crystal (LC) phases in the bulk. But in the slit pores, we found, while the distance between two walls of the pores decreases to the molecular scale, an orientationally ordered phase can form. It shows that even hard confining surfaces favor the alignment of the small anisotropic molecules. Thus we conclude that the required molecular elongation for forming LC phases will decrease in confinement. Our results indicate that some non-LC small molecules may form stable LC phases due to the inducement of confining surfaces.     

Negative flow activation energy of polymer solutions

Using equations which describe the concentration dependence of the specific viscosity of polymer solutions and the temperature dependence of intrinsic viscosity and of the Huggins viscosity parameter, conditions were sought for which solution viscosity increases with temperature (flow activation energy is negative). It was found that such an effect might appear near to the θ-temperature (regardless of the θ-temperature being LCST or UCST) particularly for polymers with high molecular weight. The concentration range depends on the system being endo- or exothermal. The conclusions are in agreement with experimental results.     

Supermolecular structures and flow birefringence in polymer solutions

Experimental studies of supermolecular structures and localized flow birefringence in solutions of high-molecular weight polymer are described. Advantage is taken of poly(ethylene oxide) and polyisobutylene. Supermolecular structures are examined with the aid of optical microscopy using freeze-dried samples of the polymer solutions. Birefringence is investigated that arises in planar elongational flow in a cross-slot cell. Flow velocities at which the onset of the localized birefringence occurs are determined. Then these velocities are correlated with viscoelastic characteristics of the solutions. The presence of a liquid-crystalline fibrillar network in the polymer solutions exhibiting flow birefringence is ascertained. The fibrils are birefringent objects. The fibrils are birefringent objects. The localized birefringence phenomenon is explained in term of the orientation of the fibrils in elongational flow. It has been shown that the onset of localized birefringence occurs at a critical Weissenberg number, the value of which is close to unity.     

Cross-streamline migration of a semiflexible polymer in a pressure driven flow

Experiments and simulations on single α-actin filaments in the Poiseuille flow through a microchannel show that the center-of-mass probability density across the channel assumes a bimodal shape as a result of pronounced cross-streamline migration. We reexamine the problem and perform Brownian dynamics simulations for a bead-spring chain with bending elasticity. Hydrodynamic interactions between the pointlike beads are taken into account by the two-wall Green tensor of the Stokes equations. Our simulations reproduce the bimodal distribution only when hydrodynamic interactions are taken into account. Numerical results on the orientational order of the end-to-end vector of the model polymer are also presented together with analytical hard-needle expressions at zero flow velocity. We derive a Smoluchowski equation for the center-of-mass distribution and carefully analyze the different contributions to the probability current that causes the bimodal distribution. As for flexible polymers, hydrodynamic repulsion explains the depletion at the wall. However, in contrast to flexible polymers, the deterministic drift current mainly determines migration away from the centerline and thereby depletion at the center. Diffusional currents due to a position-dependent diffusivity become less important with increasing polymer stiffness.     

Monte Carlo simulation of n-alkane adsorption isotherms in carbon slit pores

The single component adsorption of alkanes in carbon slit pores was studied using configurational-biased grand canonical Monte Carlo simulations. Wide ranges of temperature, pressure, alkane chain length, and slit height were studied to evaluate their effects on adsorption. Adsorption isotherms and density and orientation profiles were calculated. The behavior of long alkanes at high temperatures was found to be similar to short alkanes at lower temperatures. This suggests that the isotherms may be related through the Polanyi potential theory.     

Surface energy gradient driven convection for generating nanoscale and microscale patterned polymer films using photosensitizers

The Marangoni effect describes how fluid flows in response to gradients in surface energy. This phenomenon could be broadly harnessed to pattern the surface topography of polymer films if generalizable techniques for programming surface energy gradients existed. Here, a near UV–visible light (NUV–vis) photosensitizer, 9,10‐dibromo‐anthracene (DBA), was doped into thin films of a model polymer, poly(isobutyl methacrylate). After exposure to light through a photomask and heating above the glass transition, thermolysis of photo‐oxidized DBA and grafting to the polymer promoted flow of the film material into the exposed regions. This mechanism did not significantly alter the molecular weight of PiBMA or the film's glass transition temperature, but resulted in an increase in film surface energy as indicated by a decrease in water contact angle. Film height variations of 580 nm were produced using a mask with 12.5 μm features; a mask with 800 nm features was also employed to generate topographic features of corresponding width without expensive contacting equipment. Due to the broad absorbance spectra of DBA, highly accessible and/or unconventional light sources may be employed in this process; this advantage was demonstrated by patterning with sunlight. The nonspecific radical‐mediated nature of the DBA grafting reaction makes this a promising approach for many classes of polymers. © 2014 Wiley Periodicals, Inc. J. Polym. Sci., Part B: Polym. Phys. 2014 , 52, 1195–1202     

Molecular dynamics simulation of benzene in graphite and amorphous carbon slit pores

It is well known that confining a liquid into a pore strongly alters the liquid behavior. Investigations of the effect of confinement are of great importance for many scientific and technological applications. Here, we present a study of the behavior of benzene confined in carbon slit pores. Two types of pores are considered–graphite and amorphous carbon ones. We show that the effect of different pore structure is of crucial importance for the benzene behavior. © 2013 Wiley Periodicals, Inc.     

Simplified crossover droplet model for adsorption of pure fluids in slit pores

We present a generalized crossover (GC) model for the excess adsorption of pure fluids at a flat solid-liquid interface, which reproduces scaling behavior of the excess adsorption in the critical region and is reduced to the classical, van der Waals-type analytical model far away from the bulk critical point. In developing this model, we used the density-functional theory (DFT) approach for the order parameter profile calculations with a generalized corresponding states model for the local free-energy density. The GC DFT model well represents the available experimental adsorption data for Kr/graphite, C2H4/graphite, C3H8/graphite, CO2/silica, and SF6/graphite systems in the entire density range 0 < rho < or = 3rhoc and temperatures up to 1.7Tc. In the critical region 0.5 rhoc < r < or = 1.5rhoc and T < or = 1.15Tc, the GC DFT model is consistent with the predictions of the asymptotic renormalization-group crossover model for the critical adsorption in a semi-infinite system developed earlier. For the excess adsorption on the critical isochore, both theories predict a scaling-law behavior Gamma proportional tau(-nu+beta), but fail to reproduce a "critical depletion" of the excess adsorption along the critical isochore of the SF6/graphite system near Tc. We show that an anomalous decrease of adsorption observed in this system at tau = T/Tc - 1 < 10(-2) can be explained by finite-size effect and develop a simplified crossover droplet (SCD) model for the excess adsorption in a slit pore. With the effective size of the pore of L = 50 nm, the SCD model reproduces all available experimental data for SF6/graphite, including the critical isochore data where tau-->0, within experimental accuracy. At L > xib (where xib is a bulk correlation length) the SCD model is transformed into the GC DFT model for semi-infinite systems. Application of the SCD model to the excess adsorption of carbon dioxide on the silica gel is also discussed.     

Determination of molecular parameters from polymer solutions in shear flow

A light scattering photometer with a double cylinder shear cell has been developed which allows detection not only in one plane as usual, but in a half-spherical range around the scattering volume. The anisotropic scattering behavior of oriented and deformed polymer molecules in dilute solution was investigated. All measurements were performed on polystyrene in several viscous solvents. From measurements in the plane of flow the average orientation angle was determined. By variation of detector position and wavelength of the primary beam, the determination of all three axes of the molecular gyration space of polymer coils was possible. Compared to predictions of the well-known dynamic theories by Kuhn, Rouse and Zimm, corresponding orientation data were found while the molecular deformation ratio shows much lower values than expected.     

Layering, condensation, and evaporation of short chains in narrow slit pores

The phase behavior of short-chain fluids in slit pores is investigated by using a nonlocal-density-functional theory that takes into account the effects of segment size, chain connectivity, and van der Waals attractions explicitly. The layering and capillary condensation/evaporation transitions are examined at different chain length, temperature, pore width, and surface energy. It is found that longer chains are more likely to show hysteresis loops and multilayer adsorptions along with the capillary condensation and evaporation. Decreasing temperature favors the inclusion of layering transitions into the condensation/evaporation hysteresis loops. For large pores, the surface energy has relatively small effect on the pressures of the capillary condensation and evaporation but affects significantly on the layering pressures. It is also observed that all phase transitions within the pore take place at pressures lower than the corresponding bulk saturation pressure. The critical temperature of condensation/evaporation is always smaller than that of the bulk fluid. All coexistence curves for confined phase transitions are contained within the corresponding bulk vapor-liquid coexistence curve. As in the bulk phase, the longer the chain length, the higher are the critical temperatures of phase transitions in the pore.     

Pressure dependence of the demixing of polymer solutions determined by viscometry

Summary From the break-down in the viscosity of a polymer solution, associated with the entrance into the two phase region, the pressure dependence of the demixing temperatures of solutions of polystyrene (M = 600.000) in cyclohexane, cyclopentane, diethylmalonate and 1-phenyldecane was measured up to 1000 bar. The application of pressure increases the solubility of polystyrene in cyclopentane and diethylmalonate, but decreases that in 1-phenyldecane; in the case of cyclohexane, a pressure of optimum miscibility is observed at ca. 120 bar.The discussion of these findings, together with further information from the literature, demonstrates that current theories cannot even predict the pressure influences qualitatively. For an easy forecast of the effects, the distance of the (upper) critical solution temperature,T _c from the melting point of the solvent,T _MP, can, however, serve as a guideline: From the present material it can be concluded thatT _c is increased by pressure if (T —T _MP)/T _MP (K/K) is less than 0,20 but decreased it if it exceeds 0,25.Dedicated to Prof. Dr. K. Ueberreiter in honor of his 70th birthday.     

受限狭缝中CO_2-CH_4体系的汽液相平衡性质

受限条件下CO2-CH4体系的相平衡性质对化工工艺过程的设计具有非常重要的意义.采用Gibbs系综Monte Carlo模拟,对220K下CO2-CH4体系在主体相和受限狭缝中的相平衡性质进行了系统地研究.通过主体相模拟与实验结果比较,验证了流体分子势能参数的合理性;通过改变狭缝壁面原子的能量参数,研究了受限环境对CO2-CH4体系汽液相平衡性质的影响.与主体相相比,在硬壁狭缝中,CO2-CH4体系的露点压力增加,泡点压力降低,压力-组成相图变窄,且体系更容易达到超临界状态;在吸引狭缝中,随壁面原子能量参数的增大,CO2-CH4体系的压力-组成相图上移,临界点处CH4的摩尔分数减小,相图变窄.在体系汽液相总组成相同情况下,硬壁狭缝内体系的汽液相密度均比主体相中小;随壁面原子能量参数增大,气相密度变大、液相密度在CH4的摩尔分数较小时变大而当CH4的摩尔分数达到一定值后反而减小.在体系汽液相总组成相同时,受限环境下的汽化热比主体相的汽化热小且随壁面吸引势的增强越来越小;在主体相和硬壁狭缝中体系的汽化热随CH4含量的增加单调减小,而当壁面势能参数较大时汽化热随CH4含量增加先增大后减小.     

Diffusioosmosis of electrolyte solutions in a fine capillary slit

The steady diffusioosmotic flows of an electrolyte solution along a charged plane wall and in a capillary channel between two identical parallel charged plates generated by an imposed tangential concentration gradient are theoretically investigated. The plane walls may have either a constant surface potential or a constant surface charge density. The electrical double layers adjacent to the charged walls may have an arbitrary thickness and their electrostatic potential distributions are determined by the Poisson-Boltzmann equation. Solving a modified Navier-Stokes equation with the constraint of no net electric current arising from the cocurrent diffusion, electric migration, and diffusioosmotic convection of the electrolyte ions, the macroscopic electric field and the fluid velocity along the tangential direction induced by the imposed electrolyte concentration gradient are obtained semianalytically as a function of the lateral position in a self-consistent way. The direction of the diffusioosmotic flow relative to the concentration gradient is determined by the combination of the zeta potential (or surface charge density) of the wall, the properties of the electrolyte solution, and other relevant factors. For a given concentration gradient of an electrolyte along a plane wall, the magnitude of fluid velocity at a position in general increases with an increase in its electrokinetic distance from the wall, but there are exceptions. The effect of the lateral distribution of the induced tangential electric field and the relaxation effect in the double layer on the diffusioosmotic flow are found to be very significant.