Role of Hydrodynamic Interactions in the Deformation of Star Polymers in Poiseuille Flow

Stretching polymer in fluid flow is a vital process for studying and utilizing the physical properties of these molecules,such as DNA linearization in nanofluidic channels.We studied the role of hydrodynamic interactions(His)in stretching a free star polymer in Poiseuille flow through a tube using mesoscale hydrodynamic simulations.As increasing the flow strength,star polymers migrate toward the centerline of tube due to His,whereas toward the tube wall in the absence of His.By analyzing the end monomer distribution and the perturbed flow around the star polymer,we found that the polymer acts like a shield against the flow,leading to additional hydrodynamic drag forces that compress the arm chains in the front of the star center toward the tube axis and lift the arm chai ns at the back toward the tube wall.The balanced hydrodynamic forces freeze the polymer into a trumpet structure,where the arm chains maintain a steady strongly stretched state at high flow strength.In contrast,the polymer displays remarkably large conformational change when switching off His.Our simulation results explained the coupling between His and the structure of star polymers in Poiseuille flow.     

Cross-streamline migration of a semiflexible polymer in a pressure driven flow

Experiments and simulations on single α-actin filaments in the Poiseuille flow through a microchannel show that the center-of-mass probability density across the channel assumes a bimodal shape as a result of pronounced cross-streamline migration. We reexamine the problem and perform Brownian dynamics simulations for a bead-spring chain with bending elasticity. Hydrodynamic interactions between the pointlike beads are taken into account by the two-wall Green tensor of the Stokes equations. Our simulations reproduce the bimodal distribution only when hydrodynamic interactions are taken into account. Numerical results on the orientational order of the end-to-end vector of the model polymer are also presented together with analytical hard-needle expressions at zero flow velocity. We derive a Smoluchowski equation for the center-of-mass distribution and carefully analyze the different contributions to the probability current that causes the bimodal distribution. As for flexible polymers, hydrodynamic repulsion explains the depletion at the wall. However, in contrast to flexible polymers, the deterministic drift current mainly determines migration away from the centerline and thereby depletion at the center. Diffusional currents due to a position-dependent diffusivity become less important with increasing polymer stiffness.     

Flow-induced translocation of polymers through a fluidic channel: a dissipative particle dynamics simulation study

The dynamics of flow-induced translocation of polymers through a fluidic channel has been studied by dissipative particle dynamics (DPD) approach. Unlike implicit solvent models, the many-body energetic and hydrodynamic interactions are preserved naturally by incorporating explicit solvent particles in this approach. The no-slip wall boundary and the adaptive boundary conditions have been implemented in the modified DPD approach to model the hydrodynamic flow within a specific wall structure of fluidic channel and control the particles' density fluctuations. The results show that the average translocation time versus polymer chain length satisfies a power-law scaling of τ ～N(1.152). The conformational changes and translocation dynamics of polymers through the fluidic channel have also been investigated in our simulations, and two different translocation processes, i.e., the single-file and double-folded translocation events, have been observed in detail. These findings may be helpful in understanding the conformational and dynamic behaviors of such polymer and/or DNA molecules during the translocation processes.     

良溶剂中环形高分子单链的结构与动力学性质的模拟研究

本文采用多粒子碰撞动力学与分子动力学耦合的模拟方法研究了环形高分子单链在良溶剂中的静态与动态性质,并与线形分子进行了对比.研究发现,环形高分子链内粒子之间的平均距离小于线形链,即粒子排列得更加紧密;相应的均方回转半径也小于线形链,线形链与环形链的均方回转半径的比值为1.77;同时,环形链扩散的速度也比线形链快,两者比值为1.10.模拟结果揭示了扩散行为是排斥体积作用和流体力学相互作用耦合的结果,在扩散过程中,流体力学相互作用消减了排斥体积作用对扩散行为的贡献.此外,通过对有和没有流体力学相互作用的多粒子碰撞动力学得到的结果作对比,研究了流体力学相互作用对高分子静态和动态行为的影响,结果表明,流体力学相互作用使高分子链在极稀溶液中的扩散速度变快.     

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.     

Screening of hydrodynamic interactions in Brownian rod suspensions

We present the details and results of a simulation study addressing the dynamics and rheology of rod suspensions over a wide regime of concentrations ranging from dilute to concentrated systems. Our study compares the results of two complementary simulation methods. The first method adapts a recently proposed explicit solvent simulation strategy and incorporates both hydrodynamical effects and steric interactions between the rod units. We compare the results of such a method with those obtained from a Brownian dynamics simulation approach which retains the steric interactions but neglects the effects of hydrodynamic interactions. Overall, our results in the context of the translational and rotational diffusivities are in agreement with the hydrodynamical predictions in the dilute regime and the corresponding results of the tube model and its extensions thereof in the semidilute regimes. The latter results suggest that effects of hydrodynamic interactions on the translational and rotational diffusivities are secondary relative to the steric interactions and at best lead only to a small correction to the results of the classical tube model. Our results in the context of linear viscoelasticity also broadly confirms the predictions of the tube model for the storage and loss moduli and allows us to extract for the first time the independent hydrodynamic and Brownian contributions to the zero shear viscosity. While the relative magnitudes of these contributions are consistent with the theoretical predictions, the quantitative magnitudes are quite different from the theoretical predictions. Overall, these results confirm the validity of the hydrodynamic "screening" hypothesis and ratify the neglect of hydrodynamical stresses in quantifying the linear rheology of Brownian rod suspensions.     

Lattice-Boltzmann simulations of the dynamics of polymer solutions in periodic and confined geometries

A numerical method to simulate the dynamics of polymer solutions in confined geometries has been implemented and tested. The method combines a fluctuating lattice-Boltzmann model of the solvent [Ladd, Phys. Rev. Lett. 70, 1339 (1993)] with a point-particle model of the polymer chains. A friction term couples the monomers to the fluid [Ahlrichs and Dunweg, J. Chem. Phys. 111, 8225 (1999)], providing both the hydrodynamic interactions between the monomers and the correlated random forces. The coupled equations for particles and fluid are solved on an inertial time scale, which proves to be surprisingly simple and efficient, avoiding the costly linear algebra associated with Brownian dynamics. Complex confined geometries can be represented by a straightforward mapping of the boundary surfaces onto a regular three-dimensional grid. The hydrodynamic interactions between monomers are shown to compare well with solutions of the Stokes equations down to distances of the order of the grid spacing. Numerical results are presented for the radius of gyration, end-to-end distance, and diffusion coefficient of an isolated polymer chain, ranging from 16 to 1024 monomers in length. The simulations are in excellent agreement with renormalization group calculations for an excluded volume chain. We show that hydrodynamic interactions in large polymers can be systematically coarse-grained to substantially reduce the computational cost of the simulation. Finally, we examine the effects of confinement and flow on the polymer distribution and diffusion constant in a narrow channel. Our results support the qualitative conclusions of recent Brownian dynamics simulations of confined polymers [Jendrejack et al., J. Chem. Phys. 119, 1165 (2003) and Jendrejack et al., J. Chem. Phys. 120, 2513 (2004)].     

Alternative theory of hydrodynamic interactions in polymer solutions

The model is presented for coarse grained dynamics of macromolecules in dilute solutions. The coarse graining is achieved by dividing the polymer chain into subchains, consisting of many monomers, and spatial averaging over lengths that are large compared to the mean-square end-to-end distance of subchains and small compared to macromolecule size. Kinetic equations of the model are derived from first principles of statistical mechanics under the assumption that subchain center of mass positions and solvent flow velocity field are the only slow variables of the system. In this approach hydrodynamic interactions result from the intercomponent friction forces between polymer and solvent instead of boundary conditions on the bead surfaces as in traditional theories. The integrodifferential diffusion equation is obtained for steady flows with the kernel involving the Oseen tensor multiplied by equilibrium distribution in the space of the subchain center of mass positions.     

Dynamic light scattering and optical absorption study of poly(monosubstituted)acetylene polymers and copolymers

Polymers and copolymers obtained from the polymerization of substituted phenylacetylenes, namely PMOPPA and copolymers PPA/PMOPPA, PMOPPA/PNPPA, PPA/PNPPA, have been characterized by means of dynamic light scattering and UV–Vis absorption measurements. The analysis of the light intensity correlation function shows a similar behaviour for all the systems investigated, with three decay modes correlated to the molecular structures. The steric hyndrance and the electron polarization modify the hydrodynamic radius of the polymers and cause polymer/solvent interactions, whose strength depends on the solvent. The polymers structure does not change appreciably during the time lasting for the experiments, as assessed by UV–Vis optical tests.     

Effect of polymer adsorption on PEO-coated latex particles during salt-induced aggregation

The salt-induced aggregation of polystyrene particles in dilute aqueous solutions has been studied by means of dynamic light scattering measurements and the hydrodynamic radius of the resulting aggregates has been evaluated during the time evolution of the whole process. Poly(ethylene oxide) (PEO) polymer adsorbed on the particle surface at different amounts has been used to modify the inter-particle interactions resulting in the formation of clusters of increasing size or in the stabilization of the suspension, depending on the polymer molecular weight. The aggregation regime, i.e. a diffusion limited cluster aggregation (DLCA) occurring in the polymer-free latex suspension, is partially modified according to the polymer percentage adsorbed on the particle surface. At high polymer content, the polystyrene latex undergoes a complete steric stabilization. The deviation from a DLCA regime has been observed for different polymer contents and for polymers of different molecular weights, from 1.5 to 2000 kD. The alterations of the aggregation rates, induced by the polymer interactions, are presented and briefly discussed.     

Polymer dynamics in repton model at large fields

Polymer dynamics at large fields in Rubinstein-Duke repton model is investigated theoretically. Simple diagrammatic approach and analogy with asymmetric simple exclusion models are used to analyze the reptation dynamics of polymers. It is found that for polyelectrolytes the drift velocity decreases exponentially as a function of the external field with an exponent depending on polymer size and parity, while for polyampholytes the drift velocity is independent of polymer chain size. However, for polymers, consisting of charged and neutral blocks, the drift velocity approaches the constant limit which is determined by the size of the neutral block. The theoretical arguments are supported by extensive numerical calculations by means of density-matrix renormalization group techniques.     

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.     

高分子动力学的单链模型

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

Influence of ionic strength, sample size, and flow conditions on the retention behavior of pullulan in flow field-flow fractionation

Polymer molecular parameters such as hydrodynamic size are expected to be invariant regardless of the technique used to measure them, and to vary only, to some extent, with the solvent power and the polymer structure and properties as predicted from polymer chemistry. The hydrodynamic size of five pullulan standards derived from FlFFF in solutions of different ionic strength appears to correlate well to molecular mass as expected for neutral polymers for all fractions except that of lower mass. The correlation also holds for large amounts of injected sample even though with a slope which increases with rising polymer load. The evidence that the same result is obtained also for low sample amounts but with a higher cross-flow rate is interpreted as the manifestation of the presence of hydrodynamic interactions in concentrated polymer systems.     

Mesoscopic description of solvent effects on polymer dynamics

Solvent effects on polymer dynamics and structure are investigated using a mesoscopic solvent model that accounts for hydrodynamic interactions among the polymer beads. The simulation method combines molecular dynamics of the polymer chain, interacting with the solvent molecules through intermolecular forces, with mesoscopic multiparticle collision dynamics for the solvent molecules. Changes in the intermolecular forces between the polymer beads and mesoscopic solvent molecules are used to vary the solvent conditions from those for good to poor solvents. Polymer collapse and expansion dynamics following changes in solvent conditions are studied for homopolymer and block copolymer solutions. The frictional properties of polymers are also investigated.     

Conformational Properties of Comb-shaped Polyelectrolytes with Negatively Charged Backbone and Neutral Side Chains Studied by a Generic Coarse-grained Bead-and-Spring Model

A generic coarse-grained bead-and-spring model,mapped onto comb-shaped polycarboxylate-based(PCE)superplasticizers,is developed and studied by Langevin molecular dynamics simulations with implicit solvent and explicit counterions.The agreement on the radius of gyration of the PCEs with experiments shows that our model can be useful in studying the equilibrium sizes of PCEs in solution.The effects of ionic strength,side-chain number,and side-chain length on the conformational behavior of PCEs in solution are explored.Single-chain equilibrium properties,including the radius of gyration,end-to-end distance and persistenee length of the polymer backbone,shape-asphericity parameter,and the mean span dimension,are determined.It is found that with the increase of ionic strength,the equilibrium sizes of the polymers decrease only slightly,and a linear dependenew of the persistence length of backbone on the Debye screening length is found,in good agreement with the theory developed by Dobrynin.Increasing side-chain numbers and/or side-chain lengths increases not only the equilibrium sizes(radius of gyration and mean span)of the polymer as a whole,but also the persistence length of the backbone due to excluded volume interactions.     

聚羧酸系高性能减水剂分子在稀溶液中尺寸及其在体积排除色谱表征中表观分子量的模型研究

针对被称为"第一代聚羧酸高性能减水剂"(以下简称为MPEG-type PCE)的甲基丙烯酸(MAA)/烯酸甲酯(MAA-MPEG)梳状共聚物分子,从高分子物理基础理论出发,构建等效自由连接链模型,结合前人的理论结果和实验数据,得到了MPEG-type PCE分子的回转半径、流体力学半径及其相应的支化参数的数学表达式.在此基础上,报道了以下三方面的工作:首先,将计算结果与文献中的实验结果进行比较,检验模型的合理性;其次,利用所建立的数学模型考察主链分子量、侧链分子量和侧链接枝密度对PCE分子的回转半径和流体力学半径的影响;最后,结合近年来发展的体积排除色谱分离理论,对PCE分子的真实分子量与其常规体积排除色谱"表观分子量"(又被称为GPC分子量)两者之间的差异进行了分析.本文所提出的计算模型和数学表达式没有不确定的指前因子,可用来估算MPEG-type PCE分子在稀水溶液中的尺寸以及根据其GPC分子量估算真实分子量.     

Mean square radius of gyration and hydrodynamic radius of jointed star (dumbbell) and H-comb polymers

Equations for the mean square radius of gyration and the hydrodynamic radius for jointed stars (dumbbells) and H-combs are derived, based on random flight statistics for each subchain. Comparision with literature data on computer simulations and experimental data for H-combs show good agreement for the g-value of the mean square radius of gyration even in good solvents. This suggests that for the mean square radius of gyration the relative dimension of a H-comb relative to the linear molecule of the same degree of polymerization is not altered significantly by long range interactions, as in the case of star polymers. For the hydrodynamic radius the situation is different. Fair agreement is found when comparing our results with viscosity measurements under θ conditions, while in good solvents the viscosity of the branched molecule is more reduced than predicted by our results.     

Solution properties of ion-containing polymers in polar solvents

The placement of ionic groups within the molecular structure of a polymer produces marked modification in physical properties. A large number of studies have been performed on these ion-containing polymers, but few have focused on the effects of anion–cation interactions (i.e., counterion binding or ionization) on hydrodynamic volume, especially as the molecular structure of the solvent and nature of counterion are varied. In this study changes in hydrodynamic volume are followed through reduced viscosity measurements as a function of the abovementioned molecular parameters. The dilute solution properties of various polyelectrolytes that contain sulfonate and carboxylate groups were investigated as a function of the counterion structure, charge density, molecular weight, and solvent structure. The polymeric materials were selected because of their specific chemical structure and physical properties. In the first instance a (2-acrylamide-2 methylpropanesulfonic acid)-acrylamide-sodium vinyl sulfonate terpolymer was synthesized and subsequently neutralized with a series of bases. Viscometric measurements on these materials indicate that the nature of the cation affects the ability of the polyelectrolyte to expand its hydrodynamic volume at low polymer levels. The magnitude of the molecular expansion is shown to be due in part to the ability of the counterion to dissociate from the backbone chain, which, in turn, is directly related to the solvent structure. The changes in solution behaviour of these inomers lend support for the existence of ion pairs (i.e., site binding) and ionized moieties on the polymer chains. Measurements performed in a variety of solvent systems further confirm this interpretation. In addition, and acrylamide-sodium vinyl sulfonate copolymer was partially hydrolyzed with sodium hydroxide to study the effect of varying the charge density at a constant degree of polymerization and counterion structure. The results show that the charge density has a significant effect on the magnitude of the reduced viscosity and dilute solution behaviour. These observations, made in aqueous and nonaqueous solvents, are related to the interrelation of hydrodynamic volume, counterion concentration, and site binding. Again the controlling factor is the degree of site binding of the counterion onto the polymer backbone. Finally, we observe that the increased hydrodynamic volume affects viscosity behavior beyond the polyelectrolyte effect regime. If the average charge density on the macromolecule is relative high and/or the molecular weight is large (≥ 10⁶) sufficient intermolecular interactions will occur to produce rapid changes in reduced viscosity.     

Simulation of GPC‐Distribution Coefficients of Linear and Star‐Shaped Molecules in Spherical Pores

Simulations of the distribution coefficients of linear and star‐shaped polymers in spherical pores were performed in order to predict the GPC‐elution behavior of star‐shaped polymers relative to that of linear polymers. Self avoiding walks were generated on a tetrahedral lattice to simulate good solvent conditions. It was found that neither the molecular weight nor the mean squared radius of gyration of the polymer serves as a universal factor to determine the distribution coefficient. However, the calculated distribution coefficients correlate well with the calculated hydrodynamic radii even for different topologies. For molecules at same elution volume the ratios of molecular weights of star and linear polymer agree well with exact calculations for Gaussian chains. These ratios are nearly independent of pore geometry (spherical or cylindrical).