Effect of the wall on the velocity autocorrelation function and long-time tail of Brownian motion in a viscous compressible fluid

Brownian motion of a particle situated near a wall bounding the fluid in which it is immersed is affected by the wall. Specifically, it is assumed that a viscous compressible fluid fills a half space bounded by a plane wall, and that the fluid flow satisfies stick boundary conditions at the wall. The fluctuation-dissipation theorem shows that the velocity autocorrelation function of the Brownian particle can be calculated from the frequency-dependent admittance valid locally. The admittance can be found from the linearized Navier-Stokes equations. The t(-3/2) long-time tail of the velocity relaxation function, valid in bulk fluid, is obliterated by the wall and replaced by a t(-5/2) long-time tail of positive amplitude for motions parallel to the wall and by a t(-5/2) long-time tail of negative amplitude for motions perpendicular to the wall. In both cases the amplitude of the t(-5/2) long-time tail turns out to be independent of fluid compressibility and bulk viscosity.     

Effect of the wall on the velocity autocorrelation function and long-time tail of Brownian motion

Brownian motion of a particle situated near a wall bounding the fluid in which it is immersed is affected by the wall. Specifically, it is assumed that an incompressible viscous fluid fills a half-space bounded by a plane wall and that the fluid flow satisfies stick boundary conditions at the wall. The fluctuation-dissipation theorem shows that the velocity autocorrelation function of the Brownian particle can be calculated from the frequency-dependent admittance valid locally. It is shown that the t(-3/2) long-time tail of the velocity relaxation function, valid in bulk fluid, is obliterated and replaced by a t(-5/2) long-time tail of positive amplitude for motions parallel to the wall and by a t(-5/2) long-time tail of negative amplitude for motions perpendicular to the wall. The latter finding is at variance with an earlier calculation by Gotoh and Kaneda.     

Momentum and stress relaxation in fluids illustrating caging

The self diffusion coefficient, shear viscosity, and velocity time correlation function are calculated for a hard sphere fluid under a severe assumption, namely, the friction arises from uncorrelated binary collisions and from correlated backscattering (caging) collisions as represented in the memory function. Relaxation of the memory function from its zerotime caging value is described as a diffusion process. Derived diffusion coefficients and the shear viscosities, relative to their Enskog values decrease and increase with density, respectively, in a monotonic and gradual fashion in contrast with simulation values that show a precipitous change near the fluid-solid transition. In the present pair diffusion model, the velocity time correlation function vanishes at the proper time but its tail is overly damped relative to the simulation data. A weak breakdown of the Stokes-Einstein relation is also predicted.     

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.     

Polar solvation dynamics in supercritical fluids: a mode-coupling treatment

A mode-coupling treatment of polar solvation dynamics in supercritical fluids is presented. The equilibrium solvation time correlation function for the solute fluctuating transition frequency is obtained from the mode-coupling theory method and from molecular-dynamics simulations. The theory is shown to be in good agreement with the simulation. The solvation time correlation function exhibits three distinct time scales, with rapid initial decay, followed by a recurrence at intermediate times, and a slowly decaying long-time tail. Our theoretical analysis shows that the short-time decay arises from the coupling of the solute energy gap to the solvent polarization modes, the recurrence at intermediate times is due to the energy modes, while the slow long-time decay reflects the coupling to the number density modes.     

Dynamics of polymers in a particle-based mesoscopic solvent

We study the dynamics of flexible polymer chains in solution by combining multiparticle-collision dynamics (MPCD), a mesoscale simulation method, and molecular-dynamics simulations. Polymers with and without excluded-volume interactions are considered. With an appropriate choice of the collision time step for the MPCD solvent, hydrodynamic interactions build up properly. For the center-of-mass diffusion coefficient, scaling with respect to polymer length is found to hold already for rather short chains. The center-of-mass velocity autocorrelation function displays a long-time tail which decays algebraically as (Dt)(-3/2) as a function of time t, where D is the diffusion coefficient. The analysis of the intramolecular dynamics in terms of Rouse modes yields excellent agreement between simulation data and results of the Zimm model for the mode-number dependence of the mode-amplitude correlation functions.     

Time correlation functions and chemical reaction rates

The relation between time-dependent density correlation functions and chemical reaction rates is investigated. The case of a simple two-atom reaction in a fluid system is studied in detail. It is found that the reaction is obtainable from a three-density correlation function, taken in the long-time, long-wavelength limit. This result is analogous to the well known relation between the diffusion rate and the two-density correlation in a fluid. The present calculation is valid in an arbitrary fluid system, and is applicable in particular to reactions in a catalytic medium.     

Communication: exploring the reorientation of benzene in an ionic liquid via molecular dynamics: effect of temperature and solvent effective charge on the slow dynamics

The rotational time correlation function (RTCF) of solute benzene molecules in the ionic liquid (1-butyl-3-methylimidazolium chloride) has been studied using classical molecular dynamics simulation. The effect of solvent charge on the functional form of RTCF was investigated by comparing four force fields for the solvent where the total charge on the anion and the cation was set to ±1e, ±0.7e, ±0.5e, and 0, respectively. For all three charged solvent models, the RTCF exhibits a long-time tail where the relaxation rate exhibits a significant slowdown. This feature is strengthened by higher solvent charges as well as lower temperatures, indicating the influence of the strong Coulombic fields arising from the solvent charges. The long-time tail is caused by the extraordinarily slow solvent structural relaxation of ionic liquids compared to the time scale of their local vibrational and librational dynamics.     

Transport coefficients of square-well fluids

We analyze the analytical form of the velocity time correlation function of a hard sphere system obtained by employing generalized Langevin equation for a square-well fluid. The self-diffusion coefficient and shear viscosity have been calculated using this analytical form of velocity tcf for a square-well fluid. The addition of an attractive square-well potential in place of hard sphere leads to a substantial influence on transport coefficients. Unlike harmonic model diffusion coefficient no longer vanishes. A breakdown of the Stokes–Einstein relation is observed at low densities for a square-well fluid.     

Low-frequency velocity correlation spectrum of fluid in a rectangular microcapillary

In addition to the fast correlation for local stochastic motion, the velocity correlation function in a fluid enclosed within the pore boundaries features a slow long time-tail decay. At late times, the flow approaches that of an incompressible fluid. Here, we consider the motion of a viscous fluid, at constant temperature, in a rectangular semipermeable channel. The fluid is driven through the rectangular capillary by a uniform main pressure gradient. Tiny pressure gradients are allowed perpendicular to the main flux. We solve numerically the three-dimensional Navier-Stokes equations for the velocity field to obtain the steady solution. We then set and solve the Langevin equation for the fluid velocity. We report hydrodynamic fluctuations for the center-line velocity together with the corresponding relaxation times as a function of the size of the observing region and the Reynolds number. The effective diffusion coefficient for the fluid in the microchannel is also estimated (Deff = 1.43 x 10(-10) m2.s-1 for Re = 2), which is in accordance with measurements reported for a similar system (Stepisnik, J.; Callaghan, P. T. Physica B 2000, 292, 296-301; Stepisnik, J.; Callaghan, P. T. Magn. Reson. Imaging 2001, 19, 469-472).     

Generalized Langevin dynamics of a nanoparticle using a finite element approach: thermostating with correlated noise

A direct numerical simulation (DNS) procedure is employed to study the thermal motion of a nanoparticle in an incompressible Newtonian stationary fluid medium with the generalized Langevin approach. We consider both the Markovian (white noise) and non-Markovian (Ornstein-Uhlenbeck noise and Mittag-Leffler noise) processes. Initial locations of the particle are at various distances from the bounding wall to delineate wall effects. At thermal equilibrium, the numerical results are validated by comparing the calculated translational and rotational temperatures of the particle with those obtained from the equipartition theorem. The nature of the hydrodynamic interactions is verified by comparing the velocity autocorrelation functions and mean square displacements with analytical results. Numerical predictions of wall interactions with the particle in terms of mean square displacements are compared with analytical results. In the non-Markovian Langevin approach, an appropriate choice of colored noise is required to satisfy the power-law decay in the velocity autocorrelation function at long times. The results obtained by using non-Markovian Mittag-Leffler noise simultaneously satisfy the equipartition theorem and the long-time behavior of the hydrodynamic correlations for a range of memory correlation times. The Ornstein-Uhlenbeck process does not provide the appropriate hydrodynamic correlations. Comparing our DNS results to the solution of an one-dimensional generalized Langevin equation, it is observed that where the thermostat adheres to the equipartition theorem, the characteristic memory time in the noise is consistent with the inherent time scale of the memory kernel. The performance of the thermostat with respect to equilibrium and dynamic properties for various noise schemes is discussed.     

Molecular dynamics of some crotonates using infrared band shapes

The dipole correlation function of three esters of crotonic acid have been computed from the infrared absorption bands. Both short-time and long-time behaviour of the correlation functions are computed and discussed. The correlation times are also calculated. The short-time behaviour was compared with the correlation functions generated as an even-power series with second and fourth moments. They were also compared with the correlation functions of the freely rotating molecule.     

Steady-state hydrodynamics of a viscous incompressible fluid with spinning particles

The steady-state hydrodynamics of a viscous incompressible fluid with spinning particles is studied on the basis of extended Stokes equations. The profiles of flow velocity and spin velocity in simple flow situations may be used to determine the vortex viscosity and spin viscosity of the molecular liquid or fluid suspension. As an example, one situation studied is the flow generated by a uniform torque density in a planar layer of infinite fluid. The spinning particles drive a nearly uniform flow on either side of the layer, in opposite directions on the two sides. The Green function of the extended Stokes equations is derived. The translational and rotational friction coefficients of a sphere with no-slip boundary conditions, and the corresponding flow profiles, are calculated.     

Molecular dynamics of immiscible fluids in chemically patterned nanochannels

Molecular dynamics simulations of chain molecules are used to elucidate physical phenomena involved in flows of dense immiscible fluids in nanochannels. We first consider a force driven flow in which the channel walls are homogeneous and wetting to one fluid and nonwetting to the other fluid. The coating of the walls by the wetting fluid provides a fluctuating surface that confines the flow of the nonwetting fluid. The resulting dissipation yields stationary Poiseuille-like flows in contrast to the accelerating nature of flow in the absence of the coating. We then consider walls consisting of patches whose wetting preferences to a fluid alternate along the walls. In the resulting flow, the immiscible components exhibit periodic structures in their velocity fields such that the crests are located at the wettability steps in contrast to the behavior of a single fluid for which the crest occurs in the wetting region. We demonstrate that for a single fluid, the modulated velocity field scales with the size of the chain molecules.     

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.     

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.     

Structure of inhomogeneous polymer solutions: a density functional approach

The structure of polymer solutions confined between surfaces is studied using a density functional theory where the polymer molecules have been modeled as a pearl necklace of freely jointed hard spheres and the solvent as hard spheres. The present theory uses the concept of universality of the free energy density functional to obtain the first-order direct correlation function of the nonuniform system from that of the corresponding uniform system, calculated through the Verlet-modified type bridge function. The uniform bulk fluid direct correlation function required as input has been calculated from the reference interaction site model integral equation theory using the Percus-Yevick closure relation. The calculated results on the density profiles of the polymer as well as the solvent are shown to compare well with computer simulation results.     

The residence probability: single molecule fluorescence correlation spectroscopy and reversible geminate recombination

The residence probability of a freely diffusing particle within an open d-dimensional ball is calculated as a function of time, for an initial distribution that is either a spherical delta function or uniform within the sphere. The latter is equivalent to the autocorrelation function (ACF) of fluorescence correlation spectroscopy (FCS) when utilizing near-field scanning optical microscopy (NSOM) probes. Starting from the general equation for the Laplace transform of the residence probability, we solve it in Laplace space for any dimensionality, inverting it into the time domain in one- and three-dimensions. The short- and long-time asymptotic behaviors of the residence probability are derived and compared with the exact results. Approximations for the two-dimensional ACF are discussed, and a new approximation is derived for the NSOM-FCS ACF. Also of interest is an analytic expression for a three-dimensional ACF, which could be useful for two-photon FCS. Analogy with the binding probability for reversible geminate recombination suggests that more information could be extracted from the long-time tails in FCS experiments.     

Note: Aris-Taylor dispersion from single-particle point of view

When a point Brownian particle diffuses in a straight circular tube in the presence of a laminar stationary flow of the liquid, its effective diffusion coefficient along the tube axis increases compared to its value in the absence of flow. The effective diffusion coefficient as a function of the average fluid velocity and the tube radius is given by the Aris-Taylor formula. We give a new derivation of this formula, which is based on consideration of the axial displacement of the particle that moves in the plane normal to the tube axis along a given trajectory. The result is obtained by averaging the displacement and its square over different realizations of the particle trajectory and analyzing the long-time asymptotic behavior of the two moments.