Collective diffusion in charge-stabilized suspensions: concentration and salt effects

The authors present a joint experimental-theoretical study of collective diffusion properties in aqueous suspensions of charge-stabilized fluorinated latex spheres. Small-angle x-ray scattering and x-ray photon correlation spectroscopy have been used to explore the concentration and ionic-strength dependence of the static and short-time dynamic properties including the hydrodynamic function H(q), the wave-number-dependent collective diffusion coefficient D(q), and the intermediate scattering function over the entire accessible range. They show that all experimental data can be quantitatively described and explained by means of a recently developed accelerated Stokesian dynamics simulation method, in combination with a modified hydrodynamic many-body theory. In particular, the behavior of H(q) for de-ionized and dense suspensions can be attributed to the influence of many-body hydrodynamics, without any need for postulating hydrodynamic screening to be present, as it was done in earlier work. Upper and lower boundaries are provided for the peak height of the hydrodynamic function and for the short-time self-diffusion coefficient over the entire range of added salt concentrations.     

Sedimentation of concentrated monodisperse colloidal suspensions: role of collective particle interaction forces

The sedimentation velocities and concentration profiles of low-charge, monodisperse hydroxylate latex particle suspensions were investigated experimentally as a function of the particle concentration to study the effects of the collective particle interactions on suspension stability. We used the Kossel diffraction technique to measure the particle concentration profile and sedimentation rate. We conducted the sedimentation experiments using three different particle sizes. Collective hydrodynamic interactions dominate the particle-particle interactions at particle concentrations up to 6.5 vol%. However, at higher particle concentrations, additional collective particle-particle interactions resulting from the self-depletion attraction cause particle aggregation inside the suspension. The collective particle-particle interaction forces play a much more important role when relatively small particles (500 nm in diameter or less) are used. We developed a theoretical model based on the statistical particle dynamics simulation method to examine the role of the collective particle interactions in concentrated suspensions in the colloidal microstructure formation and sedimentation rates. The theoretical results agree with the experimentally-measured values of the settling velocities and concentration profiles.     

Fluctuating lattice-Boltzmann model for complex fluids

We develop and test numerically a lattice-Boltzmann (LB) model for nonideal fluids that incorporates thermal fluctuations. The fluid model is a momentum-conserving thermostat, for which we demonstrate how the temperature can be made equal at all length scales present in the system by having noise both locally in the stress tensor and by shaking the whole system in accord with the local temperature. The validity of the model is extended to a broad range of sound velocities. Our model features a consistent coupling scheme between the fluid and solid molecular dynamics objects, allowing us to use the LB fluid as a heat bath for solutes evolving in time without external Langevin noise added to the solute. This property expands the applicability of LB models to dense, strongly correlated systems with thermal fluctuations and potentially nonideal equations of state. Tests on the fluid itself and on static and dynamic properties of a coarse-grained polymer chain under strong hydrodynamic interactions are used to benchmark the model. The model produces results for single-chain diffusion that are in quantitative agreement with theory.     

Short-time transport properties in dense suspensions: from neutral to charge-stabilized colloidal spheres

We present a detailed study of short-time dynamic properties in concentrated suspensions of charge-stabilized and of neutral colloidal spheres. The particles in many of these systems are subject to significant many-body hydrodynamic interactions. A recently developed accelerated Stokesian dynamics (ASD) simulation method is used to calculate hydrodynamic functions, wave-number-dependent collective diffusion coefficients, self-diffusion and sedimentation coefficients, and high-frequency limiting viscosities. The dynamic properties are discussed in dependence on the particle concentration and salt content. Our ASD simulation results are compared with existing theoretical predictions, notably those of the renormalized density fluctuation expansion method of Beenakker and Mazur [Physica A 126, 349 (1984)], and earlier simulation data on hard spheres. The range of applicability and the accuracy of various theoretical expressions for short-time properties are explored through comparison with the simulation data. We analyze, in particular, the validity of generalized Stokes-Einstein relations relating short-time diffusion properties to the high-frequency limiting viscosity, and we point to the distinctly different behavior of de-ionized charge-stabilized systems in comparison to hard spheres.     

How protein surfaces induce anomalous dynamics of hydration water

Water around biomolecules slows down with respect to pure water, and both rotation and translation exhibit anomalous time dependence in the hydration shell. The origin of such behavior remains elusive. We use molecular dynamics simulations of water dynamics around several designed protein models to establish the connection between the appearance of the anomalous dynamics and water-protein interactions. For the first time we quantify the separate effect of protein topological and energetic disorder on the hydration water dynamics. When a static protein structure is simulated, we show that both types of disorder contribute to slow down water diffusion, and that allowing for protein motion, increasing the spatial dimensionality of the interface, reduces the anomalous character of hydration water. The rotation of water is, instead, altered by the energetic disorder only; indeed, when electrostatic interactions between the protein and water are switched off, water reorients even faster than in the bulk. The dynamics of water is also related to the collective structure--à voir the hydrogen bond (H-bond) network--formed by the solvent enclosing the protein surface. We show that, as expected for a full hydrated protein, when the protein surface offers pinning sites (charged or polar sites), the superficial water-water H-bond network percolates throughout the whole surface, hindering the water diffusion, whereas it does not when the protein surface lacks electrostatic interactions with water and the water diffusion is enhanced.     

Slow dynamics of a colloidal lamellar phase

We used x-ray photon correlation spectroscopy to study the dynamics in the lamellar phase of a platelet suspension as a function of the particle concentration. We measured the collective diffusion coefficient along the director of the phase, over length scales down to the interparticle distance, and quantified the hydrodynamic interaction between the particles. This interaction sets in with increasing concentration and can be described qualitatively by a simplified model. No change in the microscopic structure or dynamics is observed at the transition between the fluid and the gel-like lamellar phases.     

Prediction of collective diffusion coefficient of bovine serum albumin in aqueous electrolyte solution with hard-core two-Yukawa potential

A new method to predict concentration dependence of collective diffusion coefficient of bovine serum albumin (BSA) in aqueous electrolyte solution is developed based on the generalized Stokes-Einstein equation which relates the diffusion coefficient to the osmotic pressure. The concentration dependence of osmotic pressure is evaluated using the solution of the mean spherical approximation for the two-Yukawa model fluid. The two empirical correlations of sedimentation coefficient are tested in this work. One is for a disordered suspension of hard spheres, and another is for an ordered suspension of hard spheres. The concentration dependence of the collective diffusion coefficient of BSA under different solution conditions, such as pH and ionic strength is predicted. From the comparison between the predicted and experimental values we found that the sedimentation coefficient for the disordered suspension of hard spheres is more suitable for the prediction of the collective diffusion coefficients of charged BSA in aqueous electrolyte solution. The theoretical predictions from the hard-core two-Yukawa model coupled with the sedimentation coefficient for a suspension of hard spheres are in good agreement with available experimental data, while the hard sphere model is unable to describe the behavior of diffusion due to its neglect of the double-layer repulsive charge-charge interaction between BSA molecules.     

Langevin dynamics for rigid bodies of arbitrary shape

We present an algorithm for carrying out Langevin dynamics simulations on complex rigid bodies by incorporating the hydrodynamic resistance tensors for arbitrary shapes into an advanced rotational integration scheme. The integrator gives quantitative agreement with both analytic and approximate hydrodynamic theories for a number of model rigid bodies and works well at reproducing the solute dynamical properties (diffusion constants and orientational relaxation times) obtained from explicitly solvated simulations.     

Self-diffusion and collective diffusion of charged colloids studied by dynamic light scattering

A microemulsion of decane droplets stabilized by a nonionic surfactant film is progressively charged by substitution of a nonionic surfactant molecule by a cationic surfactant. We check that the microemulsion droplets remain identical within the explored range of volume fraction (0.02-0.18) and of the number of charges per droplet (0-40). We probe the dynamics of these microemulsions by dynamic light scattering. Despite the similar structures of the uncharged and charged microemulsions, the dynamics are very different. In the neutral microemulsion, the fluctuations of polarization relax, as is well-known, via the collective diffusion of the droplets. In the charged microemulsions, two modes of relaxation are observed. The fast one is ascribed classically to the collective diffusion of the charged droplets coupled to the diffusion of the counterions. The slow one has, to our knowledge, not been observed previously neither in similar microemulsions nor in charged spherical colloids. We show that the slow mode is also diffusive and suggest that its possible origin is the relaxation of local charge fluctuations via the local exchange of droplets bearing different numbers of charges. The diffusion coefficient associated with this mode is then the self-diffusion coefficient of the droplets.     

Structures and dynamics of thermosensitive microgel suspensions studied with three-dimensional cross-correlated light scattering

The structure factors, short- and long-time diffusion coefficients, and hydrodynamic interactions of concentrated poly(N-isopropylacryamide) microgel suspensions were measured with simultaneous static and dynamic three-dimensional cross-correlated light scattering. The data are interpreted through comparison to hard sphere theory. The structure factors are known to be described well by the hard sphere approximation. When the structure factor is fit to an effective hard sphere volume fraction and radius, the diffusion and hydrodynamic interactions are also well described by the hard sphere model. We demonstrate that one single hard sphere volume fraction is sufficient to describe the microgel structures, hydrodynamic interactions, and long- and short-time collective diffusion coefficients. This result is surprising because the particle form of the microgels at these temperatures is not rigid, but rather "fuzzy" spheres with dangling polymer chains.     

Pair diffusion, hydrodynamic interactions, and available volume in dense fluids

We calculate the pair diffusion coefficient D(r) as a function of the distance r between two hard sphere particles in a dense monodisperse fluid. The distance-dependent pair diffusion coefficient describes the hydrodynamic interactions between particles in a fluid that are central to theories of polymer and colloid dynamics. We determine D(r) from the propagators (Green's functions) of particle pairs obtained from molecular dynamics simulations. At distances exceeding ~3 molecular diameters, the calculated pair diffusion coefficients are in excellent agreement with predictions from exact macroscopic hydrodynamic theory for large Brownian particles suspended in a solvent bath, as well as the Oseen approximation. However, the asymptotic 1/r distance dependence of D(r) associated with hydrodynamic effects emerges only after the pair distance dynamics has been followed for relatively long times, indicating non-negligible memory effects in the pair diffusion at short times. Deviations of the calculated D(r) from the hydrodynamic models at short distances r reflect the underlying many-body fluid structure, and are found to be correlated to differences in the local available volume. The procedure used here to determine the pair diffusion coefficients can also be used for single-particle diffusion in confinement with spherical symmetry.     

Dynamics of proteins: light scattering study of dilute and dense colloidal suspensions of eye lens homogenates

We report a dynamic light scattering study on protein suspensions of bovine lens homogenates at conditions (pH and ionic strength) similar to the physiological ones. Light scattering data were collected at two temperatures, 20 and 37 degrees C, over a wide range of concentrations from the very dilute limit up to the dense regime approaching the physiological lens concentration. A comparison with experimental data from intact bovine lenses was advanced, revealing differences between dispersions and lenses at similar concentrations. In the dilute regime, two scattering entities were detected and identified with the long-time self-diffusion modes of alpha-crystallins and their aggregates, which naturally exist in lens nucleus. Upon increasing protein concentration, significant changes in time correlation function were observed starting at approximately 75 mg ml(-1), where a new mode originating from collective diffusive motions becomes visible. Self-diffusion coefficients are temperature insensitive, whereas the collective diffusion coefficient depends strongly on temperature revealing a reduction of the net repulsive interparticle forces with decreasing temperature. While there are no rigorous theoretical approaches on particle diffusion properties for multicomponent, nonideal hard sphere polydispersed systems, as the suspensions studied here, a discussion of the volume fraction dependence of the long-time self-diffusion coefficient in the context of existing theoretical approaches was undertaken. This study is purported to provide some insight into the complex light scattering pattern of intact lenses and the interactions between the constituent proteins that are responsible for lens transparency. This would lead to understand basic mechanisms of specific protein interactions that lead to lens opacification (cataract) under pathological conditions.     

Diffusion of spheres in isotropic and nematic networks of rods: electrostatic interactions and hydrodynamic screening

Translational diffusion of a small charged tracer sphere in isotropic and nematic suspensions of long and thin charged rods is investigated as a function of ionic strength and rod concentration. A theory for the diffusive properties of a small sphere is developed, where both (screened) hydrodynamic interactions and charge interactions between the tracer sphere and the rod network are analyzed. Hydrodynamic interactions are formulated in terms of the hydrodynamic screening length. As yet, there are no independent theoretical predictions for the hydrodynamic screening length for rod networks. Experimental tracer-diffusion data are presented for various ionic strengths as a function of the rod concentration, both in the isotropic and nematic states. Orientational order parameters are measured for the same ionic strengths as a function of the rod concentration. The hydrodynamic screening length is determined from these experimental data and scaling relations obtained from the above mentioned theory. For the isotropic networks, a master curve is found for the hydrodynamic screening length as a function of the rod concentration. For the nematic networks the screening length turns out to be a very sensitive function of the orientational order parameter.     

Anomalous nanoparticle diffusion in polymer solutions and melts: a mode-coupling theory study

Mode-coupling theory is employed to study diffusion of nanoparticles in polymer melts and solutions. Theoretical results are directly compared with molecular dynamics simulation data for a similar model. The theory correctly reproduces the effects of the nanoparticle size, mass, particle-polymer interaction strength, and polymer chain length on the nanoparticle diffusion coefficient. In accord with earlier experimental, simulation, and theoretical work, it is found that when the polymer radius of gyration exceeds the nanoparticle radius, the Stokes-Einstein relation underestimates the particle diffusion coefficient by as much as an order of magnitude. Within the mode-coupling theory framework, a microscopic interpretation of this phenomenon is given, whereby the total diffusion coefficient is decomposed into microscopic and hydrodynamic contributions, with the former dominant in the small particle limit, and the latter dominant in the large particle limit. This interpretation is in agreement with previous mode-coupling theory studies of anomalous diffusion of solutes in simple dense fluids.     

Molecular dynamics simulation studies of the transport and adsorption of a charged macromolecule onto a charged adsorbent solid surface immersed in an electrolytic solution

Molecular dynamics simulations were performed in order to study the transport and adsorption of a charged macromolecule (desmopressin) onto a charged solid surface in an electrolytic solution. The strong Coulombic interaction from the charged solid surface represents the major force for accelerating, orienting, entrapping in the electrical double layer, and adsorbing the macromolecule onto the charged solid surface. The macromolecule is flattened as it approaches the charged surface, giving rise to a stronger surface exclusion effect that shields surface sites. When adsorbed, the macromolecule is restrained by a surface interaction more than one hundred times stronger than the thermal energy, of which 99.8% results from the strong dominant Coulombic interaction, and trapped by a hydration layer adjacent to the surface. This leads to zero lateral displacement of the adsorbed macromolecule and indicates that surface diffusion is a physically implausible mechanism in similar systems. Explicit solvent is required for realistic representation of the macromolecular structure and the surface interaction energy. The adsorbed macromolecule also decreased the electrostatic potential gradient perpendicular to the charged solid surface and introduced additional electrostatic potential gradients laterally. The results obtained from the molecular dynamics simulations confirm the importance of electrophoretic migration and support the physical mechanisms used in a macroscopic continuum model that predicts an overshoot in the concentration of a charged macromolecule in the adsorbed phase under certain conditions of pH and ionic strength.     

Simulations of polyelectrolyte dynamics in an externally applied electric field in confined geometry

We consider the dynamics of charged polymers in free solution in a slit geometry under the influence of an electrical field, applied at an angle to the plane parallel walls of the confinement. The simulations are carried out using the Brownian dynamics method with explicit counterions and implicit hydrodynamics. The hydrodynamic interactions between all the particles and the plane parallel walls are taken into account using a diffusion matrix which depends on slit geometry and the actual polyelectrolyte-solute conformations. We observe a selective transport of the charged polymers, as a function of the degree of polymerization and slit height.     

Onset of slow dynamics in difluorotetrachloroethane glassy crystal

Complementary neutron spin-echo and x-ray experiments and molecular-dynamics simulations have been performed on difluorotetrachloroethane (CFCl2-CFCl2) glassy crystal. Static, single-molecule reorientational dynamics and collective dynamics properties are investigated. Our results confirm the strong analogy between molecular liquids and plastic crystals. The orientational disorder is characterized at different temperatures and a change in the nature of rotational dynamics is observed. A careful check of the rotational diffusion model is performed using self-angular correlation functions Cl with high l values and compared to results obtained on molecular liquids composed of A-B dumbbells. Below the crossover temperature at which slow dynamics emerge, we show that some scaling predictions of the mode coupling theory hold and that alpha-relaxation times and nonergodicity parameters are controlled by the nontrivial static correlations.     

Thermodynamics and diffusion in size-symmetric and asymmetric dense electrolytes

MD simulation results for model size-symmetric and asymmetric electrolytes at high densities and temperatures (well outside the liquid-gas coexistence region) are generated and analyzed focusing on thermodynamic and diffusion properties. An extension of the mean spherical approximation for electrolytes originally derived for charged hard sphere fluids is adapted to these systems by exploiting the separation of short range and Coulomb interaction contributions intrinsic to these theoretical models and is found to perform well for predicting equation of state quantities. The diffusion coefficients of these electrolytes can also be reasonably well predicted using entropy scaling ideas suitably adapted to charged systems and mixtures. Thus, this approach may provide an avenue for studying dense electrolytes or complex molecular systems containing charged species at high pressures and temperatures.     

Mesoscopic model for diffusion-influenced reaction dynamics

A hybrid mesoscopic multiparticle collision model is used to study diffusion-influenced reaction kinetics. The mesoscopic particle dynamics conserves mass, momentum, and energy so that hydrodynamic effects are fully taken into account. Reactive and nonreactive interactions with catalytic solute particles are described by full molecular dynamics. Results are presented for large-scale, three-dimensional simulations to study the influence of diffusion on the rate constants of the A + C <==> B + C reaction. In the limit of a dilute solution of catalytic C particles, the simulation results are compared with diffusion equation approaches for both the irreversible and reversible reaction cases. Simulation results for systems where the volume fraction phi of catalytic spheres is high are also presented, and collective interactions among reactions on catalytic spheres that introduce volume fraction dependence in the rate constants are studied.     

Toward realistic modeling of dynamic processes in cell signaling: quantification of macromolecular crowding effects

One of the major factors distinguishing molecular processes in vivo from biochemical experiments in vitro is the effect of the environment produced by macromolecular crowding in the cell. To achieve a realistic modeling of processes in the living cell based on biochemical data, it becomes necessary, therefore, to consider such effects. We describe a protocol based on Brownian dynamics simulation to characterize and quantify the effect of various forms of crowding on diffusion and bimolecular association in a simple model of interacting hard spheres. We show that by combining the elastic collision method for hard spheres and the mean field approach for hydrodynamic interaction (HI), our simulations capture the correct dynamics of a monodisperse system. The contributions from excluded volume effect and HI to the crowding effect are thus quantified. The dependence of the results on size distribution of each component in the system is illustrated, and the approach is applied as well to the crowding effect on electrostatic-driven association in both neutral and charged environments; values for effective diffusion constants and association rates are obtained for the specific conditions. The results from our simulation approach can be used to improve the modeling of cell signaling processes without additional computational burdens.