Stochastic entangled chain dynamics of dense polymer solutions

We propose an adjustable-parameter-free, entangled chain dynamics model of dense polymer solutions. The model includes the self-consistent dynamics of molecular chains and solvent by describing the former via coarse-grained polymer dynamics that incorporate hydrodynamic interaction effects, and the latter via the forced Stokes equation. Real chain elasticity is modeled via the inclusion of a Pincus regime in the polymer's force-extension curve. Excluded volume effects are taken into account via the combined action of coarse-grained intermolecular potentials and explicit geometric tracking of chain entanglements. We demonstrate that entanglements are responsible for a new (compared to phantom chain dynamics), slow relaxation mode whose characteristic time scale agrees very well with experiment. Similarly good agreement between theory and experiment is also obtained for the equilibrium chain size. We develop methods for the solution of the model in periodic flow domains and apply them to the computation of entangled polymer solutions in equilibrium. We show that the number of entanglements Π agrees well with the number of entanglements expected on the basis of tube theory, satisfactorily reproducing the latter's scaling of Π with the polymer volume fraction φ. Our model predicts diminishing chain size with concentration, thus vindicating Flory's suggestion of excluded volume effects screening in dense solutions. The predicted scaling of chain size with φ is consistent with the heuristic, Flory theory based value.     

Dynamic light scattering from non-entangled wormlike micellar solutions

The fuzzy cylinder theory, originally proposed for conventional polymer solutions, was applied to wormlike micellar solutions to take into account effects of the intermicellar collision and hydrodynamic interaction on the self-diffusion of wormlike micelles in solution at finite concentrations. Previously reported apparent hydrodynamic radius data obtained by dynamic light scattering for non-entangled wormlike micelles formed in aqueous solution by non-ionic surfactants, polyoxyethylene monoalkyl ethers C(i)E(j), were analyzed by this theory to estimate the persistence length q of the wormlike micelles. The results of q estimated were consistent with those obtained from radius of gyration data obtained by static light scattering.     

高分子动力学的单链模型

高分子单链模型是高分子稀溶液理论研究的基本模型.对其进行深入地分析,不仅有助于解决高分子稀溶液体系中溶液黏度和分子链扩散等基本问题,而且能够增进人们对高分子链结构与溶液性质间关联性的理解.虽然基于经典连续性介质力学的流体动力学理论可以定性,甚至半定量地获得稀溶液的一些重要性质,但是,随着科学技术的发展,人们从分子水平上建立了许多描述高分子稀溶液性质的模型和理论,期望能够定量地描述高分子稀溶液的性质.本文以高分子稀溶液中3个典型的单链模型为例(包括:不含流体力学相互作用的Rouse模型、含二体流体力学相互作用的Zimm模型和含多体流体力学相互作用的部分穿透球模型),综述高分子稀溶液的重要性质,并详细地给出其动力学方程的推导过程及其重要的研究进展.特别是,对于Rouse模型,本文还将其预言结果拓展到了短链高分子流体体系;此外,还介绍了这一领域的关键科学问题、发展前景和研究方向.     

COMPUTER SIMULATION OF POLYMER SOLUTION THERMODYNAMICS

The statistical counting method for the computer simulation of the ther-modynamic quantities of polymer solution has been reviewed. The calculating results fora single athermal chain confirm the theory of the renormalization group. The results forthe athermal solution are consistent with the scaling law of the osmotic pressure with theexponent 2.25. The results for a single chain with the segmental interaction are in a goodagreement with the exact results obtained by the direct counting method. The results forthe polymer solution show us that the Flory-Huggins parameter is strongly dependent onboth the polymer concentration and the interaction energy between segments. Monte carlo simulation; Polymer solution; Thermodynamic quantities;Translational entropy; Conformational entropy; Scaling law     

Conformation and diffusion behavior of ring polymers in solution: a comparison between molecular dynamics, multiparticle collision dynamics, and lattice Boltzmann simulations

We have studied the effect of chain topology on the structural properties and diffusion of polymers in a dilute solution in a good solvent. Specifically, we have used three different simulation techniques to compare the chain size and diffusion coefficient of linear and ring polymers in solution. The polymer chain is modeled using a bead-spring representation. The solvent is modeled using three different techniques: molecular dynamics (MD) simulations with a particulate solvent in which hydrodynamic interactions are accounted through the intermolecular interactions, multiparticle collision dynamics (MPCD) with a point particle solvent which has stochastic interactions with the polymer, and the lattice Boltzmann method in which the polymer chains are coupled to the lattice fluid through friction. Our results show that the three methods give quantitatively similar results for the effect of chain topology on the conformation and diffusion behavior of the polymer chain in a good solvent. The ratio of diffusivities of ring and linear polymers is observed to be close to that predicted by perturbation calculations based on the Kirkwood hydrodynamic theory.     

Chain overlap and entanglements in dilute polymer solutions: Brownian dynamics simulation

We have analyzed chain conformations and the existence — or otherwise — of chain overlaps and entanglements in dilute polymer solutions (at concentrations c < C^*, c^* = critical concentration). The fundamental problem of existence of chain overlaps in dilute solutions is also related to the drag reduction phenomenon (DR). Some experimental results pertinent to DR are explained in terms of entanglements even for solutions at concentrations defined in ppm. We report results of Brownian dynamics simulations of polymer solutions in which the equations of motion of the chains are solved by using the Langevin equation. Chains move according to actions of a systematic frictional force and a randomly fluctuating force w(t), where t is time. In addition, a shear flow field can be introduced into the model. To evaluate the structure of polymer chains in solution we have devised a measure of interchain contacts and two different measures of entanglements. The results for c = 0.3 c^* demonstrate that both chain entanglements and overlaps take place even in dilute solution. They also confirm predictions from an earlier combinatorial model.     

Electrostatic relaxation and hydrodynamic interactions for self-diffusion of ions in electrolyte solutions

The concentration dependence of self-diffusion of ions in solutions at large concentrations has remained an interesting yet unsolved problem. Here we develop a self-consistent microscopic approach based on the ideas of mode-coupling theory. It allows us to calculate both contributions which influence the friction of a moving ion: the ion atmosphere relaxation and hydrodynamic interactions. The resulting theory provides an excellent agreement with known experimental results over a wide concentration range. Interestingly, the mode-coupling self-consistent calculation of friction reveal a nonlinear coupling between the hydrodynamic interactions and the ion atmosphere relaxation which enhances ion diffusion by reducing friction, particularly at intermediate ion concentrations. This rather striking result has its origin in the similar time scales of the relaxation of the ion atmosphere relaxation and the hydrodynamic term, which are essentially given by the Debye relaxation time. The results are also in agreement with computer simulations, with and without hydrodynamic interactions.     

Hydrodynamic interaction in polymer solutions simulated with dissipative particle dynamics

The authors analyzed extensively the dynamics of polymer chains in solutions simulated with dissipative particle dynamics (DPD), with a special focus on the potential influence of a low Schmidt number of a typical DPD fluid on the simulated polymer dynamics. It has been argued that a low Schmidt number in a DPD fluid can lead to underdevelopment of the hydrodynamic interaction in polymer solutions. The authors' analyses reveal that equilibrium polymer dynamics in dilute solution, under typical DPD simulation conditions, obey the Zimm [J. Chem. Phys. 24, 269 (1956)] model very well. With a further reduction in the Schmidt number, a deviation from the Zimm model to the Rouse model is observed. This implies that the hydrodynamic interaction between monomers is reasonably developed under typical conditions of a DPD simulation. Only when the Schmidt number is further reduced, the hydrodynamic interaction within the chains becomes underdeveloped. The screening of the hydrodynamic interaction and the excluded volume interaction as the polymer volume fraction is increased are well reproduced by the DPD simulations. The use of soft interaction between polymer beads and a low Schmidt number do not produce noticeable problems for the simulated dynamics at high concentrations, except for the entanglement effect which is not captured in the simulations.     

Effects of ion atmosphere relaxation on dipole isomerization reactions in electrolyte solutions

A theoretical study of the effects of ion atmosphere relaxation on the rate of a model dipole isomerization reaction in electrolyte solutions is presented. The time-dependent ion atmosphere friction is calculated by using a molecular hydrodynamic theory which properly includes the static and dynamic interionic correlations through ionic structure factors and van Hove functions. The rate constant is determined by employing the well-known Grote-Hynes theory. Numerical results are obtained for the time-dependent ion atmosphere friction and for the rate of isomerization reaction in electrolyte solutions of varying ion concentration and dielectric constant. It is found that the ion atmosphere friction can have significant effects in reducing the rate of isomerization below the prediction of equilibrium solvation transition state theory.     

Dynamical problems in the theory of polymer solutions

Qualitative discrepancies are found between what is predicted by available theory and what is actually observed, for several concentration regimes of the dynamical properties of polymer solutions. The difficulties are most severe, from the standpoint of experiment or simulation as well as theory, for the entanglement concentration regime. However, the classical problems of chain polymers in dilute solution are not fully understood. For example, the constants of proportionality that relate hydrodynamic radii to the radius of gyration, in the nondraining limit and in theta solvents, may not be universal constants. That is, the proportionality constants may vary with polymer and solvent species. Discrepancies between theory and experiment are discussed for the two different systems, dilute chains and semidilute rods. Speculation is offered on the resolution of these difficulties.     

Fourth‐order hydrodynamic contribution to the polymer self‐diffusion coefficient

A hydrodynamic scattering treatment of interacting polymer chains is extended to obtain the five‐point chain–chain–chain–chain–chain hydrodynamic interaction tensor. The tensor is used to calculate the second‐order concentration correction to the self‐diffusion coefficient of a polymer in solution. The self‐similarity assumption of the hydrodynamic scaling model of polymer dynamics is tested against these calculations. © 2004 Wiley Periodicals, Inc. J Polym Sci Part B: Polym Phys 42: 1663–1670, 2004     

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.     

改进的聚合物溶液理论

在高分子溶液理论中引入Gibbs分布 ,用统计物理学方法重新推导出了聚合物溶液的热力学公式 .将高分子溶液的自由能和熵分三部分进行了计算 ,无热平动部分 ,无热构象部分和构象有热部分 .无热平动自由能和无热构象自由能分别等于Flory Huggins混合自由能公式的前两项 ,构象有热部分引入了Gibbs分布 ,考虑了链段 溶剂分子相互作用对高分子构象的影响 .在分子间的相互作用足够小时 ,又回到了FH公式     

Time-dependent Diffusion Coefficient and Conventional Diffusion Constant of Nanoparticles in Polymer Melts by Mode-coupling Theory (cited: 2)

Time-dependent diffusion coefficient and conventional diffusion constant are calculated and analyzed to study diffusion of nanoparticles in polymer melts. A generalized Langevin equa-tion is adopted to describe the diffusion dynamics. Mode-coupling theory is employed to calculate the memory kernel of friction. For simplicity, only microscopic terms arising from binary collision and coupling to the solvent density fluctuation are included in the formalism. The equilibrium structural information functions of the polymer nanocomposites required by mode-coupling theory are calculated on the basis of polymer reference interaction site modelwith Percus-Yevick closure. The effect of nanoparticle size and that of the polymer size are clarified explicitly. The structural functions, the friction kernel, as well as the diffusion coefficient show a rich variety with varying nanoparticle radius and polymer chain length. We find that for small nanoparticles or short chain polymers, the characteristic short time non-Markov diffusion dynamics becomes more prominent, and the diffusion coefficient takes longer time to approach asymptotically the conventional diffusion constant. This constant due to the microscopic contributions will decrease with the increase of nanoparticle size, while increase with polymer size. Furthermore, our result of diffusion constant from mode-coupling theory is compared with the value predicted from the Stokes-Einstein relation. It shows that the microscopic contributions to the diffusion constant are dominant for small nanoparticles or long chain polymers. Inversely, when nanonparticle is big, or polymer chain is short, the hydrodynamic contribution might play a significant role.     

Length-scale dependent relaxation shear modulus and viscoelastic hydrodynamic interactions in polymer liquids

A quantitative theory of hydrodynamic interactions in unentangled polymer melts and concentrated solutions is presented. The study is focussed on the pre-Rouse transient time regimes (t < τ(R), the Rouse relaxation time) where the hydrodynamic response is governed mainly by the viscoelastic effects. It is shown that transient viscoelastic hydrodynamic interactions are not suppressed (screened) at large distances and are virtually independent of polymer molecular mass. A number of transient regimes of unusual and qualitatively different behavior of isotropic and anisotropic hydrodynamic response functions are elucidated. The regimes are characterized in terms of two main length-scale dependent characteristic times: momentum spreading time τ(i) ∝ r(4∕3) and viscoelastic time τ(?) ∝ r(4). It is shown that for t > τ(i) the viscoelastic hydrodynamic interactions can be described in terms of the time or length scale dependent effective viscosity which, for t < τ(R) and/or for r < R(coil), turns out to be much lower than the macroscopic "polymer" viscosity η(m). The theory also involves a quantitative analysis of the length-scale dependent stress relaxation in polymer melts. The general predictions for hydrodynamic interactions in thermostated systems with Langevin friction are obtained as well.     

Coagulation of tobacco mosaic virus in alcohol-water-LiCl solutions

The coagulation and colloidal stability of tobacco mosaic virus (TMV) in alcohol-water-LiCl solutions were studied. Without the addition of LiCl salt, the coagulation was promoted by the increase of hydrophobicity of the alcohols that is proportional to their alkyl chain length and concentration. Addition of the LiCl salt reduced the electrostatic repulsion between TMV particles resulting in coagulation in methanol-water and ethanol-water solutions. In water-alcohol-LiCl mixture, the coagulation of TMV was driven by both the hydrophobic interaction of the solution and the screening effect of the salt simultaneously. To understand the particle-particle interaction during the coagulation, the interaction energy was calculated using DLVO theory. Considering the electrostatic repulsive energy, van der Waals attractive energy, and hydrophobic interaction energy, the total energy profiles were obtained. The experiment and model calculation results indicated that the increase of alcohol concentration would increase hydrophobic attraction energy so that the coagulation is promoted. These results provide the fundamental understanding on the coagulation of biomolecular macromolecules.     

On the theory of dilute polyelectrolyte solutions: Extensions,refinements, and experimental tests

We synthesize a quantitative theory for the radius of gyration, second virial coefficient, intrinsic viscosity, and friction coefficient for polyelectroytes in dilute solution from existing treatments of electrostatic and hydrodynamic interactions within and among wormlike chains. Comparison with data for K-PSS demonstrates the importance of accounting for nonlinearities in the electrostatics and the finite diameter of the polymer backbone.     

Viscoelastic properties of an isolated polymer chain with arbitrary flexibility,chain length and hydrodynamic interactions from 4‐dimensional Dirac propagator

Kholodenko's theory of semiflexible polymer chains, the conformation and properties of which are obtained from the Dirac propagator, shows applicability to dilute solutions of semiflexible polymers of arbitrary persistence and contour lengths by calculating the static scattering function and the squared end‐to‐end distance of the polymer chain. In the present work, the theory is extended and applied to obtain the intrinsic viscosity with consideration of hydrodynamic interactions. The intrinsic viscosity formula is derived as function of chain length and persistence length. The hydrodynamic interactions are also taken into account following the Kirkwood and Riseman scheme. From this calculation, we obtain the general expression for the intrinsic viscosity and diffusion coefficients covering the whole range of chain flexibilities without confusion with the excluded volume effects. Calculated limiting values of hydrodynamical observables are in complete agreement with those known for random coils and rigid rods.     

Transport coefficients of aqueous dodecyltrimethylammonium bromide solutions: comparison between experiments, analytical calculations and numerical simulations

We study dynamical properties of ionic species in aqueous solutions of dodecyltrimethylammonium bromide, for several concentrations below and above the critical micellar concentration (cmc). New experimental determinations of the electrical conductivity are given which are compared to results obtained from an analytical transport theory; transport coefficients of ions in these solutions above the cmc are also computed from Brownian dynamics simulations. Analytical calculations as well as the simulation treat the solution within the framework of the continuous solvent model. Above the cmc, three ionic species are considered: the monomer surfactant, the micelle and the counterion. The analytical transport theory describes the structural properties of the electrolyte solution within the mean spherical approximation and assumes that the dominant forces which determine the deviations of transport processes from the ideal behavior (i.e., without any interactions between ions) are hydrodynamic interactions and electrostatic relaxation forces. In the simulations, both direct interactions and hydrodynamic interactions between solutes are taken into account. The interaction potential is modeled by pairwise repulsive 1/r(12) interactions and Coulomb interactions. The input parameters of the simulation (radii and self-diffusion coefficients of ions at infinite dilution) are partially obtained from the analytical transport theory which fits the experimental determinations of the electrical conductivity. Both the electrical conductivity of the solution and the self-diffusion coefficients of each species computed from Brownian dynamics are compared to available experimental data. In every case, the influence of hydrodynamic interactions (HIs) on the transport coefficients is investigated. It is shown that HIs are crucial to obtain agreement with experiments. In particular, the self-diffusion coefficient of the micelle, which is the largest and most charged species in the present system, is enhanced when HIs are included whereas the diffusion coefficients of the monomer and the counterion are roughly not influenced by HIs.     

A NOVEL INTERPRETATION OF CONCENTRATION DEPENDENCE OF VISCOSITY OF DILUTE POLYMER SOLUTION*

The concentration dependence of the reduced viscosity of dilute polymer solution is interpreted in the light of anew concept of the self-association of polymer chains in dilute solution. The apparent self-association constant is defined asthe molar association constant divided by the molar mass of individual polymer chain and is numerically interconvertiblewith the Huggins coefficient. The molar association constant is directly proportional to the effective hydrodynamic volume ofthe polymer chain in solution and is irrespective of the chain architecture. The effective hydrodynamic volume accounts forthe non-spherical conformation of a short polymer chain in solution and is a product of a shape factor and hydrodynamicvolume. The observed enhancement of Huggins coefficient for short chain and branched polymer is satisfactorily interpretedby the concept of self-association. The concept of self-association allows us to predict the existence of a boundaryconcentration C_s (dynamic contact concentration) which divides the dilute polymer solution into two regions.