Hydrodynamic interaction solutions

In dilute solutions, tire dependence of tire diffusion coefficient on tire molecular weight is different from tliat found in melts, eitlier entangled or not. This difference is due to tire presence of hydrodynamic interactions among tire solvent molecules. Such interactions arise from tire necessity to transfer solvent molecules from tire front to tire back of a moving particle. The motion of tire solvent gives rise to a flow field which couples all molecules over a... 

In dilute polymer solutions, hydrodynamic interactions lead to a concerted motion of tire whole polymer chain and tire surrounding solvent. The folded chains can essentially be considered as impenneable objects whose hydrodynamic radius is / / is tire gyration radius defined as... 

Equation (23) predicts a dependence of xR on M2. Experimentally, it was found that the relaxation time for flexible polymer chains in dilute solutions obeys a different scaling law, i.e. t M3/2. The Rouse model does not consider excluded volume effects or polymer-solvent interactions, it assumes a Gaussian behavior for the chain conformation even when distorted by the flow. Its domain of validity is therefore limited to modest deformations under 0-conditions. The weakest point, however, was neglecting hydrodynamic interaction which will now be discussed. 

The Zimm model predicts correctly the experimental scaling exponent xx ss M3/2 determined in dilute solutions under 0-conditions. In concentrated solution and melts, the hydrodynamic interaction between the polymer segments of the same chain is screened by the host molecules (Eq. 28) and a flexible polymer coil behaves much like a free-draining chain with a Rouse spectrum in the relaxation times. 

First approaches to approximating the relaxation time on the basis of molecular parameters can be traced back to Rouse [33]. The model is based on a number of boundary assumptions (1) the solution is ideally dilute, i.e. intermolecular interactions are negligible (2) hydrodynamic interactions due to disturbance of the medium velocity by segments of the same chain are negligible and (3) the connector tension F(r) obeys an ideal Hookean force law. 

There are several attractive features of such a mesoscopic description. Because the dynamics is simple, it is both easy and efficient to simulate. The equations of motion are easily written and the techniques of nonequilibriun statistical mechanics can be used to derive macroscopic laws and correlation function expressions for the transport properties. Accurate analytical expressions for the transport coefficient can be derived. The mesoscopic description can be combined with full molecular dynamics in order to describe the properties of solute species, such as polymers or colloids, in solution. Because all of the conservation laws are satisfied, hydrodynamic interactions, which play an important role in the dynamical properties of such systems, are automatically taken into account. 

Multiparticle collision dynamics can be combined with full molecular dynamics in order to describe the behavior of solute molecules in solution. Such hybrid MPC-MD schemes are especially useful for treating polymer and colloid dynamics since they incorporate hydrodynamic interactions. They are also useful for describing reactive systems where diffusive coupling among solute species is important. 

Most descriptions of the dynamics of molecular or particle motion in solution require a knowledge of the frictional properties of the system. This is especially true for polymer solutions, colloidal suspensions, molecular transport processes, and biomolecular conformational changes. Particle friction also plays an important role in the calculation of diffusion-influenced reaction rates, which will be discussed later. Solvent multiparticle collision dynamics, in conjunction with molecular dynamics of solute particles, provides a means to study such systems. In this section we show how the frictional properties and hydrodynamic interactions among solute or colloidal particles can be studied using hybrid MPC-MD schemes. 

Brownian motion theory may be generalized to treat systems with many interacting B particles. Such many-particle Langevin equations have been investigated at a molecular level by Deutch and Oppenheim [58], A simple system in which to study hydrodynamic interactions is two particles fixed in solution at a distance Rn- The Langevin equations for the momenta P, (i = 1,2)... 

It is known that polymer dynamics is strongly influenced by hydrodynamic interactions. When viewed on a microscopic level, a polymer is made from molecular groups with dimensions in the angstrom range. Many of these monomer units are in close proximity both because of the connectivity of the chain and the fact that the polymer may adopt complicated conformations in solution. Polymers are solvated by a large number of solvent molecules whose molecular dimensions are comparable to those of the monomer units. These features make the full treatment of hydrodynamic interactions for polymer solutions very difficult. 

The dynamical properties of polymer molecules in solution have been investigated using MPC dynamics [75-77]. Polymer transport properties are strongly influenced by hydrodynamic interactions. These effects manifest themselves in both the center-of-mass diffusion coefficients and the dynamic structure factors of polymer molecules in solution. For example, if hydrodynamic interactions are neglected, the diffusion coefficient scales with the number of monomers as D Dq /Nb, where Do is the diffusion coefficient of a polymer bead and N), is the number of beads in the polymer. If hydrodynamic interactions are included, the diffusion coefficient adopts a Stokes-Einstein formD kltT/cnr NlJ2, where c is a factor that depends on the polymer chain model. This scaling has been confirmed in MPC simulations of the polymer dynamics [75]. 

The nature of the solvent influences both the structure of the polymer in solution and its dynamics. In good solvents the polymer adopts an expanded configuration and in poor solvents it takes on a compact form. If the polymer solution is suddenly changed from good to poor solvent conditions, polymer collapse from the expanded to compact forms will occur [78], A number of models have been suggested for the mechanism of the collapse [79-82], Hydrodynamic interactions are expected to play an important part in the dynamics of the collapse and we show how MPC simulations have been used to investigate this problem. Hybrid MD-MPC simulations of the collapse dynamics have been carried out for systems where bead-solvent interactions are either explicitly included [83] or accounted for implicitly in the multiparticle collision events [84, 85]. 

Multiparticle collision dynamics describes the interactions in a many-body system in terms of effective collisions that occur at discrete time intervals. Although the dynamics is a simplified representation of real dynamics, it conserves mass, momentum, and energy and preserves phase space volumes. Consequently, it retains many of the basic characteristics of classical Newtonian dynamics. The statistical mechanical basis of multiparticle collision dynamics is well established. Starting with the specification of the dynamics and the collision model, one may verify its dynamical properties, derive macroscopic laws, and, perhaps most importantly, obtain expressions for the transport coefficients. These features distinguish MPC dynamics from a number of other mesoscopic schemes. In order to describe solute motion in solution, MPC dynamics may be combined with molecular dynamics to construct hybrid schemes that can be used to explore a variety of phenomena. The fact that hydrodynamic interactions are properly accounted for in hybrid MPC-MD dynamics makes it a useful tool for the investigation of polymer and colloid dynamics. Since it is a particle-based scheme it incorporates fluctuations so that the reactive and nonreactive dynamics in small systems where such effects are important can be studied. 

The earliest and simplest approach in this direction starts from Langevin equations with solutions comprising a spectrum of relaxation modes [1-4], Special features are the incorporation of entropic forces (Rouse model, [6]) which relax fluctuations of reduced entropy, and of hydrodynamic interactions (Zimm model, [7]) which couple segmental motions via long-range backflow fields in polymer solutions, and the inclusion of topological constraints or entanglements (reptation or tube model, [8-10]) which are mutually imposed within a dense ensemble of chains. 

The dynamics of polymer solution is governed by the hydrodynamic interaction between moving segments mediated by the solvent [90,91]. This interaction is long range and couples the motions of the different segments strongly. In this... 

In fact, the diffusion constant in solutions has the form of an Einstein diffusion of hard spheres with radius Re. For a diffusing chain the solvent within the coil is apparently also set in motion and does not contribute to the friction. Thus, the long-range hydrodynamic interactions lead, in comparison to the Rouse model, to qualitatively different results for both the center-of-mass diffusion—which is not proportional to the number of monomers exerting friction - as well as for the segment diffusion - which is considerably accelerated and follows a modified time law t2/3 instead of t1/2. 

A comparison with Burchard s first cumulant calculations shows qualitative agreement, in particular with respect to the position of the minimum. Quantitatively, however, important differences are obvious. Both the sharpness as well as the amplitude of the phenomenon are underestimated. These deviations may originate from an overestimation of the hydrodynamic interaction between segments. Since a star of high f internally compromises a semi-dilute solution, the back-flow field of solvent molecules will be partly screened [40,117]. Thus, the effects of hydrodynamic interaction, which in general eases the renormalization effects owing to S(Q) [152], are expected to be weaker than assumed in the cumulant calculations and thus the minimum should be more pronounced than calculated. Furthermore, since for Gaussian chains the relaxation rate decreases... 

It is generally accepted that in semi-dilute solutions under good solvent conditions both the excluded volume interactions and the hydrodynamic interactions are screened owing to the presence of other chains [4,5,103], With respect to the correlation lengths (c) and H(c) there is no consensus as to whether these quantities have to be equal [11] or in general would be different [160],... 

The hydrodynamic interaction is introduced in the Zimm model as a pure intrachain effect. The molecular treatment of its screening owing to presence of other chains requires the solution of a complicated many-body problem [11, 160-164], In some cases, this problem can be overcome by a phenomenological approach [40,117], based on the Zimm model and on the additional assumption that the average hydrodynamic interaction in semi-dilute solutions is still of the same form as in the dilute case. 

Fig. 59. Incomplete screening of hydrodynamic interactions in semi-dilute polymer solutions. Presentation of different regimes which are passed with increasing concentration. A,C Unscreened and screened Zimm relaxation, respectively, B enhanced Rouse relaxation. (Reprinted with permission from [12]. Copyright 1987 Vieweg and Sohn Verlagsgemeinschaft, Wiesbaden)...

Fig. 64. Single-chain behavior in semi-dilute PDMS/d-chlorbenzene solutions. Line-shape parameter (3 as a function of Q at the concentration c = 0.18 and c = 0.45, indicating the occurance of two crossover effects, as predicted by the concept of incompletely screened hydrodynamic interactions. (----), (---) asymptotic Zimm and Rouse behavior, respectively. (Reprinted with per-...

In the general case, when arbitrary interaction profiles prevail, the particle deposition rate must be obtained by solving the complete transport equations. The first numerical solution of the complete convective diffusional transport equations, including London-van der Waals attraction, gravity, Brownian diffusion and the complete hydrodynamical interactions, was obtained for a spherical collector [89]. Soon after, numerical solutions were obtained for a panoplea of other collector geometries... 

All measurements, of course, have to be made at a finite concentration. This implies that interparticle interactions cannot be fully neglected. However, in very dilute solutions we can safely assume that more than two particles have only an extremely small chance to meet [72]. Thus only the interaction between two particles has to be considered. There are two types of interaction between particles in solution. One results from thermodynamic interactions (repulsion or attraction), and the other is caused by the distortion of the laminar fiow due to the presence of the macromolecules. If the particles are isolated only the laminar flow field is perturbed, and this determines the intrinsic viscosity but when the particles come closer together the distorted flow fields start to overlap and cause a further increase of the viscosity. The latter is called the hydrodynamic interaction and was calculated by Oseen to various approximations [3,73]. Figure 7 elucidates the effect. 

In order to resolve these challenges, it is essential to account for chain connectivity, hydrodynamic interactions, electrostatic interactions, and distribution of counterions and their dynamics. It is possible to identify three distinct scenarios (a) polyelectrolyte solutions with high concentrations of added salt, (b) dilute polyelectrolyte solutions without added salt, and (c) polyelectrolyte solutions above overlap concentration and without added salt. If the salt concentration is high and if there is no macrophase separation, the polyelectrolyte solution behaves as a solution of neutral polymers in a good solvent, due to the screening of electrostatic interaction. Therefore for scenario... 

Now, the decay rate of the incoherent dynamic structure factor is proportional to fc(i+2v)/v xherefore, the -dependence of the decay rate for salt-free solutions is independent of whether the hydrodynamic interaction is present or not. 

This reduces to the result of Eqs. (119) and (175) for dilute solutions ( Rg). In the Rouse regime where hydrodynamic interaction is screened, D becomes... 

Therefore in this Rouse regime of unentangled semidilute solutions where hydrodynamic interaction is screened, both the reduced viscosity and reduced modulus decrease with increase in polymer concentration in salt free solutions... 

In inhnitely dilute solutions where full hydrodynamic interaction between segments is present as in the Kirkwood-Riseman-Zimm model, the electrophoretic mobility is independent of N in both low (kR 1) and high (kR 1) salt concentrations ... 

But p decreases with salt concentration with an apparent exponent of k which changes from 0 at low salt concentration to — at high salt concentrations. The N-independence of p arises from a cancellation between hydrodynamic interaction and electrostatic coupling between the polyelectrolyte and other ions in the solution. It is to be noted that the self-translational diffusion coefficient D is proportional to as in the Zimm model with full... 

---