Molecular dynamics simulation study of friction and diffusion of a tracer in a Lennard–Jones solvent

The friction and diffusion coefficients of a tracer in a Lennard–Jones (LJ) solvent are evaluated by equilibrium molecular dynamics simulations in a microcanonical ensemble. The solvent molecules interact through a repulsive LJ force each other and the tracer of diameter σ₂ interacts with the solvent molecules through the same repulsive LJ force with a different LJ parameter σ. Positive deviation of the diffusion coefficient D of the tracer from a Stokes–Einstein behavior is observed and the plot of 1/D versus σ₂ shows a linear behavior. It is also observed that the friction coefficient ζ of the tracer varies linearly with σ₂ in accord with the prediction of the Stokes formula but shows a smaller slope than the Stokes prediction. When the values of ratios of sizes between the tracer and solvent molecules are higher than 5 approximately, the behavior of the friction and diffusion coefficients is well described by the Einstein relation D = k _B T/ζ, from which the tracer is considered as a Brownian particle.     

Diffusion coefficient of Brownian particle in rough micro-channel

This study examines the feasibility of using of the lattice Boltzmann method to determine how the surface roughness of a quadrate channel affects the diffusion coefficient of Brownian particle(s). The surface was represented by a regular array of spheres. Surface roughness reduced the diffusion coefficient of the Brownian particle(s) because of a change in the velocity autocorrelation function decay and in pressure. Additionally, the neighboring particles increased the diffusion coefficient of Brownian particle.     

Calculation of the diffusion coefficient in a liquid based on the force covariance function

Methods of optimal use of statistical information on the Brownian motion of molecules in a liquid for calculating the diffusion coefficient are considered. Relationships between covariance functions for the displacement of a molecule, its velocity, and the force acting on the molecule are established. It is proposed the method of calculation of the diffusion coefficient estimates should be made using the force covariance function, which enables to reduce the time of simulation of the molecular motion. The proposed method is used for callculating diffusion coefficients for some systems.     

Diffusion of nanoparticles in dense fluids

The diffusion process of a single spherical nanoparticle immersed in a fluid solvent is studied by molecular dynamics simulations. When the nanoparticle mass stays constant, it is shown that, at short times, the decay of the nanoparticle velocity autocorrelation function is strongly modified when the particle diameter increases. It is also shown that, at large times, the characteristic algebraic decay induced by the hydrodynamic correlations between the solvated particle and the solvent presents a scaling behavior depending on the particle diameter.     

Dynamics of swelling/contracting hard spheres surmised by an irreversible Langevin equation

The diffusion of molecules through uniform homogeneous materials can readily be described by Brownian motion or generalizations thereof. The further generalization of these models to describe molecular diffusion through heterogeneous and nonstationary solvents is much less understood. Phenomenological nonstationary generalizations of the generalized Langevin equation (GLE) have earlier been developed satisfying the fluctuation-dissipation relationship in quasi-equilibrium limits while exhibiting somewhat complex behavior away from equilibrium. This reduced-dimensional representation should be capable of describing the diffusion of a particle through a colloidal suspension whose average particle size is tuned by an external driving force such as pH. A simple particle model of such a process involves the motion of a hard-sphere particle in an explicit environment of swelling hard spheres. The velocity autocorrelation functions observed in a large number of simulations of the particle model under various swelling rates agree precisely with those of a single form of the nonstationary phenomenological model. Though this procedure is not an explicit projection of the mechanical system onto the nonstationary GLE, it does show that the latter correctly describes the dynamics of the projected coordinate--namely, diffusion of the solute--under nonequilibrium conditions. Both nonequilibrium solvent models lead to behavior reminiscent of beta-relaxation processes at packing fractions substantially below that of the glass transition.     

Scaling of the memory function and Brownian motion

It has been recently shown that the velocity autocorrelation function of a tracer particle immersed in a simple liquid scales approximately with the inverse of its mass. With increasing mass the amplitude is systematically reduced and the velocity autocorrelation function tends to a slowly decaying exponential, which is characteristic for Brownian motion. We give here an analytical proof for this behavior and comment on the usual explanation for Brownian dynamics which is based on the assumption that the memory function is proportional to a Dirac distribution. We also derive conditions for Brownian dynamics of a tracer particle which are entirely based on properties of its memory function.     

Communication: Turnover behavior of effective mobility in a tube with periodic entropy potential

Using Brownian dynamics simulations, we study the effective mobility and diffusion coefficient of a point particle in a tube formed from identical compartments of varying diameter, as functions of the driving force applied along the tube axis. Our primary focus is on how the driving force dependences of these transport coefficients are modified by the changes in the compartment shape. In addition to monotonically increasing or decreasing behavior of the effective mobility in periodic entropy potentials reported earlier, we now show that the effective mobility can even be nonmonotonic in the driving force.     

Estimating the viscoelastic moduli of complex fluids from observation of Brownian motion of a particle confined to a harmonic trap

A procedure is proposed to estimate the viscoelastic properties of a complex fluid from the behavior of the velocity autocorrelation function of a suspended Brownian particle, trapped in a harmonic potential. The procedure is tested for a model complex fluid with a given frequency-dependent shear viscosity. The analysis shows that the procedure can provide a rather accurate prediction of the viscoelastic properties of the fluid on the basis of experimental data on the velocity autocorrelation function of the trapped Brownian particle in a limited range of time.     

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.     

Spectrum of position fluctuations of a Brownian particle bound in a harmonic trap near a plane wall

The spectrum of position fluctuations of a Brownian particle bound in a harmonic trap near a plane wall is calculated from an approximate result for the Fourier transform of the velocity autocorrelation function. Both a no-slip and a perfect slip boundary condition at the wall are considered. In both cases at low frequency the calculated spectrum differs markedly from recent experimental data. It is suggested that a partial slip boundary condition with a frequency-dependent slip coefficient may explain the experimental results.     

Theory of nanoparticle diffusion in unentangled and entangled polymer melts

We propose a statistical dynamical theory for the violation of the hydrodynamic Stokes-Einstein (SE) diffusion law for a spherical nanoparticle in entangled and unentangled polymer melts based on a combination of mode coupling, Brownian motion, and polymer physics ideas. The non-hydrodynamic friction coefficient is related to microscopic equilibrium structure and the length-scale-dependent polymer melt collective density fluctuation relaxation time. When local packing correlations are neglected, analytic scaling laws (with numerical prefactors) in various regimes are derived for the non-hydrodynamic diffusivity as a function of particle size, polymer radius-of-gyration, tube diameter, degree of entanglement, melt density, and temperature. Entanglement effects are the origin of large SE violations (orders of magnitude mobility enhancement) which smoothly increase as the ratio of particle radius to tube diameter decreases. Various crossover conditions for the recovery of the SE law are derived, which are qualitatively distinct for unentangled and entangled melts. The dynamical influence of packing correlations due to both repulsive and interfacial attractive forces is investigated. A central finding is that melt packing fraction, temperature, and interfacial attraction strength all influence the SE violation in qualitatively different directions depending on whether the polymers are entangled or not. Entangled systems exhibit seemingly anomalous trends as a function of these variables as a consequence of the non-diffusive nature of collective density fluctuation relaxation and the different response of polymer-particle structural correlations to adsorption on the mesoscopic entanglement length scale. The theory is in surprisingly good agreement with recent melt experiments, and new parametric studies are suggested.     

Statistical mechanical theory for steady state systems. VIII. General theory for a Brownian particle driven by a time- and space-varying force

A Brownian particle subject to a time- and space-varying force is studied with the second entropy theory for nonequilibrium statistical mechanics. A fluctuation expression is obtained for the second entropy of the path, and this is maximized to obtain the most likely path of the particle. Two approaches are used, one based on the velocity correlation function and one based on the position correlation function. The approaches are a perturbation about the free particle result and are exact for weak external forces. They provide a particularly simple way of including memory effects in time-varying driven diffusion. The theories are tested against computer simulation data for a Brownian particle trapped in an oscillating parabolic well. They accurately predict the phase lag and amplitude as a function of drive frequency, and they account quantitatively for the memory effects that are important at high frequencies and that are missing in the simplest Langevin equation.     

Combining molecular dynamics with Lattice Boltzmann: a hybrid method for the simulation of (charged) colloidal systems

We present a hybrid method for the simulation of colloidal systems that combines molecular dynamics (MD) with the Lattice Boltzmann (LB) scheme. The LB method is used as a model for the solvent in order to take into account the hydrodynamic mass and momentum transport through the solvent. The colloidal particles are propagated via MD and they are coupled to the LB fluid by viscous forces. With respect to the LB fluid, the colloids are represented by uniformly distributed points on a sphere. Each such point [with a velocity V(r) at any off-lattice position r] is interacting with the neighboring eight LB nodes by a frictional force F = xi0(V(r)-u(r)), with xi0 being a friction coefficient and u(r) being the velocity of the fluid at the position r. Thermal fluctuations are introduced in the framework of fluctuating hydrodynamics. This coupling scheme has been proposed recently for polymer systems by Ahlrichs and Dunweg [J. Chem. Phys. 111, 8225 (1999)]. We investigate several properties of a single colloidal particle in a LB fluid, namely, the effective Stokes friction and long-time tails in the autocorrelation functions for the translational and rotational velocity. Moreover, a charged colloidal system is considered consisting of a macroion, counterions, and coions that are coupled to a LB fluid. We study the behavior of the ions in a constant electric field. In particular, an estimate of the effective charge of the macroion is yielded from the number of counterions that move with the macroion in the direction of the electric field.     

Transport of a heated granular gas in a washboard potential

We study numerically the motion of a one dimensional array of Brownian particles in a washboard potential, driven by an external stochastic force and interacting via short range repulsive forces. In particular, we investigate the role of instantaneous elastic and inelastic collisions on the system dynamics and transport. The system displays a locked regime, where particles may move only via activated processes and a running regime where particles drift along the direction of the applied field. By tuning the value of the friction parameter controlling the Brownian motion we explore both the overdamped dynamics and the underdamped dynamics. In the two regimes we considered the mobility and the diffusivity of the system as functions of the tilt and other relevant control parameters such as coefficient of restitution, particle size, and total number of particles. We find that while in the overdamped regime the results for the interacting systems present similarities with the known noninteracting case, in the underdamped regime the inelastic collisions determine a rich variety of behaviors among which is an unexpected enhancement of the inelastic diffusion.     

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.     

Thermodiffusion in model nanofluids by molecular dynamics simulations

In this work, a new algorithm is proposed to compute single particle (infinite dilution) thermodiffusion using nonequilibrium molecular dynamics simulations through the estimation of the thermophoretic force that applies on a solute particle. This scheme is shown to provide consistent results for model nanofluids in the liquid state (spherical nonmetallic nanoparticles+Lennard-Jones fluid) where it appears that thermodiffusion amplitude, as well as thermal conductivity, decreases with nanoparticle concentration. Then, by changing the nature of the nanoparticle (size, mass, and internal stiffness) and that of the solvent (quality and viscosity), various trends are exhibited. In all cases, the single particle thermodiffusion is positive, i.e., the nanoparticle tends to migrate toward the cold area. The single particle thermal diffusion coefficient is shown to be independent of the size of the nanoparticle (diameter of 0.8-4 nm), whereas it increases with the quality of the solvent and is inversely proportional to the viscosity of the fluid. In addition, this coefficient is shown to be independent of the mass of the nanoparticle and to increase with the stiffness of the nanoparticle internal bonds. Besides, for these configurations, the mass diffusion coefficient behavior appears to be consistent with a Stokes-Einstein-like law.     

Design of a microfabricated magnetic cell separator.

A new magnetic separation idea utilizing several ideas from microfabrication and nanomagnetics is presented. The basic idea comes from our earlier work using asymmetry in obstacles and Brownian motion to effect separation of objetcs [10] by moving them in streams whose angle to the hydrodynamic average velocity is a function of the diffusion coefficient of the object. The device we propose here is not technically a Brownian ratchet device but uses the idea of force which acts at angle to the hydrodynamic flow. In our case, the force is generated by a magnetic field gradient which comes from an array of magnetized wires which lie at an angle 0 to a hydrodynamic field flow. The sum of the hydrodynamic force and the magnetic force create a new vector which as in the case of the Brownian ratchet moves the cell out of the main stream direction.     

Communication: drift velocity of Brownian particle in a periodically tapered tube induced by a time-periodic force with zero mean: dependence on the force period

We study the drift of a Brownian particle in a periodically tapered tube, induced by a longitudinal time-periodic force of amplitude ∣F∣ that alternates in sign every half-period. The focus is on the velocity dependence on the force period, which is usually considered not tractable analytically. For large ∣F∣ we derive an analytical solution that gives the velocity as a function of the amplitude and the period of the force as well as the geometric parameters of the tube. The solution shows how the velocity decreases from its maximum value to zero as the force period decreases from infinity (adiabatic regime) to zero. Our analytical results are in excellent agreement with those obtained from 3D Brownian dynamics simulations.     

Dynamics of liquid–liquid demixing: mesoscopic simulations of polymer solutions

In this article, the mesoscopic simulation method dissipative particle dynamics (DPD) is applied to study the dynamics of polymer–solvent liquid–liquid phase separation. It will be shown that the degree of branching has a pronounced effect on the radius of gyration and the centre of mass diffusion of the polymer. Based on the simulation results it can be concluded that the difference in chemical potential between the mixed and the demixed state is the main driving force behind the centre of mass diffusion (and thus phase separation), rather than the reduced radius of gyration due to to polymer chain collapse.     

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.