Hybrid method coupling fluctuating hydrodynamics and molecular dynamics for the simulation of macromolecules

We present a hybrid computational method for simulating the dynamics of macromolecules in solution which couples a mesoscale solver for the fluctuating hydrodynamics (FH) equations with molecular dynamics to describe the macromolecule. The two models interact through a dissipative Stokesian term first introduced by Ahlrichs and Dunweg [J. Chem. Phys. 111, 8225 (1999)]. We show that our method correctly captures the static and dynamical properties of polymer chains as predicted by the Zimm model. In particular, we show that the static conformations are best described when the ratio sigma/b=0.6, where sigma is the Lennard-Jones length parameter and b is the monomer bond length. We also find that the decay of the Rouse modes' autocorrelation function is better described with an analytical correction suggested by Ahlrichs and Dunweg. Our FH solver permits us to treat the fluid equation of state and transport parameters as direct simulation parameters. The expected independence of the chain dynamics on various choices of fluid equation of state and bulk viscosity is recovered, while excellent agreement is found for the temperature and shear viscosity dependence of center of mass diffusion between simulation results and predictions of the Zimm model. We find that Zimm model approximations start to fail when the Schmidt number Sc < or approximately 30. Finally, we investigate the importance of fluid fluctuations and show that using the preaveraged approximation for the hydrodynamic tensor leads to around 3% error in the diffusion coefficient for a polymer chain when the fluid discretization size is greater than 50 A.     

高分子动力学的单链模型

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

Reduction of the effective shear viscosity in polymer solutions due to crossflow migration in microchannels: Effective viscosity models based on DPD simulations

Molecular dynamics simulations (dissipative particle dynamics–DPD) were developed and used to quantify wall-normal migration of polymer chains in microchannel Poseuille flow. Crossflow migration due to viscous interaction with the walls results in lowered polymer concentration near the channel walls. A larger fraction of the total flow volume becomes depleted of polymer when the channel width h decreases into the submicron range, significantly reducing the effective viscosity. The effective viscosity was quantified in terms of channel width and Weissenberg number Wi, for 5% polymer volume fraction in water. Algebraic models for the depletion width δ(Wi, h) and effective viscosity μ_e(δ/h, Wi) were developed, based on the hydrodynamic theory of Ma and Graham and our simulation results. The depletion width model can be applied to longer polymer chains after a retuning of the polymer persistence length and the corresponding potential/thermal energy ratio.     

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.     

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.     

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)].     

End-monomer Dynamics in Semiflexible Polymers

Spurred by an experimental controversy in the literature, we investigate the end-monomer dynamics of semiflexible polymers through Brownian hydrodynamic simulations and dynamic mean-field theory. Precise experimental observations over the last few years of end-monomer dynamics in the diffusion of double-stranded DNA have given conflicting results: one study indicated an unexpected Rouse-like scaling of the mean squared displacement (MSD) ?r(2)(t)? ~ t(1/2) at intermediate times, corresponding to fluctuations at length scales larger than the persistence length but smaller than the coil size; another study claimed the more conventional Zimm scaling ?r(2)(t)? ~ t(2/3) in the same time range. Using hydrodynamic simulations, analytical and scaling theories, we find a novel intermediate dynamical regime where the effective local exponent of the end-monomer MSD, α(t) = d log?r(2)(t)?/d log t, drops below the Zimm value of 2/3 for sufficiently long chains. The deviation from the Zimm prediction increases with chain length, though it does not reach the Rouse limit of 1/2. The qualitative features of this intermediate regime, found in simulations and in an improved mean-field theory for semiflexible polymers, in particular the variation of α(t) with chain and persistence lengths, can be reproduced through a heuristic scaling argument. Anomalously low values of the effective exponent α are explained by hydrodynamic effects related to the slow crossover from dynamics on length scales smaller than the persistence length to dynamics on larger length scales.     

Long-time dynamics of Rouse-Zimm polymers in dilute solutions with hydrodynamic memory

The dynamics of flexible polymers in dilute solutions is studied taking into account the hydrodynamic memory, as a consequence of fluid inertia. As distinct from the Rouse-Zimm (RZ) theory, the Boussinesq friction force acts on the monomers (beads) instead of the Stokes force, and the motion of the solvent is governed by the nonstationary Navier-Stokes equations. The obtained generalized RZ equation is solved approximately using the preaveraging of the Oseen tensor. It is shown that the time correlation functions describing the polymer motion essentially differ from those in the RZ model. The mean-square displacement (MSD) of the polymer coil is at short times approximately t(2) (instead of approximately t). At long times the MSD contains additional (to the Einstein term) contributions, the leading of which is approximately t. The relaxation of the internal normal modes of the polymer differs from the traditional exponential decay. It is displayed in the long-time tails of their correlation functions, the longest lived being approximately t(-3/2) in the Rouse limit and t(-5/2) in the Zimm case, when the hydrodynamic interaction is strong. It is discussed that the found peculiarities, in particular, an effectively slower diffusion of the polymer coil, should be observable in dynamic scattering experiments.     

A practical method to avoid bond crossing in two-dimensional dissipative particle dynamics simulations

Dissipative particle dynamics (DPD) simulation technique is an effective method targeted on mesoscopic simulations in which the interactions between particles are soft. As a result, it inevitably causes bond crossing and interpenetration between particles. Here we develop a practical method based on the two-dimensional DPD model which can extremely reduce the possibility of bond crossing. A rigid core is added to each particle by modifying the form of the conservative force in DPD so that the particles cannot penetrate each other. Then by adjusting the spring constant of the bond, we can impose a simple geometry constraint so that the bond crossing can hardly take place. Furthermore, we take into account an analytic geometry constraint in the polymerization model of DPD by which we can successfully avoid the severe bond crossing problem during bond generation in two dimensions. A parabola fitting between the pressure and the particle number density shows that our modified DPD model with small rigid cores can still be mapped onto the Flory-Huggins model, and the mesoscopic length scale of our simulations does not change. By analyzing the mean-square displacement of the innermost monomer and the center of mass of the chains, we find a t(8/15) power law of the polymer dynamics in our model instead of the Rouse prediction supporting the recent results in literature.     

From molecular dynamics to hydrodynamics: a novel Galilean invariant thermostat

This paper proposes a novel thermostat applicable to any particle-based dynamic simulation. Each pair of particles is thermostated either (with probability P) with a pairwise Lowe-Andersen thermostat [C. P. Lowe, Europhys. Lett. 47, 145 (1999)] or (with probability 1-P) with a thermostat that is introduced here, which is based on a pairwise interaction similar to the Nosé-Hoover thermostat. When the pairwise Nosé-Hoover thermostat dominates (low P), the liquid has a high diffusion coefficient and low viscosity, but when the Lowe-Andersen thermostat dominates, the diffusion coefficient is low and viscosity is high. This novel Nosé-Hoover-Lowe-Andersen thermostat is Galilean invariant and preserves both total linear and angular momentum of the system, due to the fact that the thermostatic forces between each pair of the particles are pairwise additive and central. We show by simulation that this thermostat also preserves hydrodynamics. For the (noninteracting) ideal gas at P = 0, the diffusion coefficient diverges and viscosity is zero, while for P > 0 it has a finite value. By adjusting probability P, the Schmidt number can be varied by orders of magnitude. The temperature deviation from the required value is at least an order of magnitude smaller than in dissipative particle dynamics (DPD), while the equilibrium properties of the system are very well reproduced. The thermostat is easy to implement and offers a computational efficiency better than (DPD), with better temperature control and greater flexibility in terms of adjusting the diffusion coefficient and viscosity of the simulated system. Applications of this thermostat include all standard molecular dynamic simulations of dense liquids and solids with any type of force field, as well as hydrodynamic simulation of multiphase systems with largely different bulk viscosities, including surface viscosity, and of dilute gases and plasmas.     

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.     

Polymer chain dynamics at interfaces: role of boundary conditions at solid interface

Using classical molecular dynamics simulations, we study the dynamical properties of a single polymer chain dissolved in an explicit solvent and strongly adsorbed at solid-liquid interface. To circumvent a serious challenge posed by finite size effects due to long-range hydrodynamic effects, we developed a correction procedure that substantially limits the finite size effects. Concurrently, we provide an analysis of distinctly different size effects in the directions transverse and normal to the interface. We find that on analytically smooth interfaces, corresponding to the slip boundary condition, the motions of the polymer chain and the surrounding solvent are hydrodynamically coupled. This leads to the chain diffusion coefficient D scaling with the chain degree of polymerization N as D approximately N(-3/4), consistent with the Zimm dynamics for strongly adsorbed chains. Introduction of transverse forces at the interface results in loss of correlation between the motion of the polymer chain and the solvent. Consequently, D approximately N(-1), which is a characteristic of the Rouse dynamics.     

Computer simulation of linkage of two ring chains

We performed off-lattice Monte Carlo simulations of links of two model ring chains with chain length N up to 32,768 in the theta solution or amorphous bulk state by using a random walk model (Model I), and molecular dynamics simulations of two model ring chains in solution with excluded volume interaction (Model II) to investigate topological effects on the geometry of link and ring conformation. In the case of Model I, the mean squared linking number, its distribution, and the size of two chains with fixed linking number are investigated. Our simulation results confirm the previous theoretical prediction that the mean squared linking number decays as pe(-qs(2)) with the distance of centers of chain mass s, where p and q are found to be chain length dependent and q asymptotically approaches to 0.75 as chain length increases. The linking number distribution of two chains has a universal form for long chains, but our simulation results clearly show that the distribution function deviates from the Gaussian distribution, a fact not predicted by any previous theoretical work. A scaling prediction is proposed to predict the link size, and is checked for our simulations for the Model II. The simulation results confirmed the scaling prediction of the blob picture that the link with linking number m occupies a compact volume of m blobs, and the size of the link is asymptotic to R(L) ≈ bN(ν)m(1/3-ν), where N is the chain length, and v is the Flory exponent of polymer in solutions.     

Nonequilibrium simulations of filled polymer networks: searching for the origins of reinforcement and nonlinearity

We present nonequilibrium dissipative particle dynamics (DPD) simulations of cross-linked elastomers containing solid filler particles at 30% volume fraction. We study systematically the effect of the morphology (dispersed or aggregated particles) and of the effective particle-particle interactions. In addition, we have experimented by replacing the standard harmonic DPD bonds with other potential functions, conceived to deal with the finite extensibility of the polymer chains and the possibility of a slow equilibrium between strongly and weakly adsorbed chains at the rubber-filler interface. The simulation results shed some light on the basic mechanisms of rubber reinforcement, including the nonlinearity and history dependence commonly known as "Payne effect" and "Mullins effect."     

Comparison of Brownian dynamics algorithms with hydrodynamic interaction

The hydrodynamic interaction is an essential effect to consider in Brownian dynamics simulations of polymer and nanoparticle dilute solutions. Several mathematical approaches can be used to build Brownian dynamics algorithms with hydrodynamic interaction, the most common of them being the exact but time demanding Cholesky decomposition and the Chebyshev polynomial expansion. Recently, Geyer and Winter [J. Chem. Phys. 130, 1149051 (2009)] have proposed a new approximation to treat the hydrodynamic interaction that seems quite efficient and is increasingly used. So far, a systematic comparison among those approaches has not been clearly made. In this paper, several features and the efficiency of typical implementations of those approaches are evaluated by using bead-and-spring chain models. The different sensitivity to the bead overlap detected for the different implementations may be of interest to select the suitable algorithm for a given simulation.     

Modeling of polyethylene and poly (L-lactide) polymer blends and diblock copolymer: chain length and volume fraction effects on structural arrangement

Dissipative particle dynamics (DPD), a mesoscopic simulation approach, has been used to investigate the chain length effect on the structural property of the immiscible polyethylene (PE)/poly(L-lactide) (PLLA) polymer in a polymer blend and in a system with their diblock copolymer. In this work, the interaction parameter in DPD simulation, related to the Flory-Huggins interaction parameter chi, is estimated by the calculation of mixing energy for each pair of components in molecular dynamics simulation. The immiscibility property of PE and PLLA polymers induces the phase separation and exhibits different architectures at different volume fractions. In order to observe the structural property, the radius of gyration is used to observe the detailed arrangement of the polymer chains. It shows that the structure arrangement of a polymer chain is dependent on the phase structure and has a significantly different structural arrangement character for the very short chains in the homopolymer and copolymers. The chain length effect on the degree of stretching or extension of polymers has also been observed. As the chain length increases, the chain exhibits more stretching behavior at lamellae, perforated lamellae, and cylindrical configurations, whereas the chain exhibits a similar degree of stretching or extension at the cluster configuration.     

流场驱动高分子链迁移穿过微通道的耗散粒子动力学模拟

利用耗散粒子动力学模拟方法研究了高分子链在流场驱动作用下迁移穿过微通道过程中的链构象变化和动力学行为.在足够大的流场力驱动作用下,高分子链在沿着流场方向逐渐被拉伸,从而能够穿过管径小于其自身尺寸的微通道.耗散粒子动力学模拟结果表明高分子链的迁移过程主要分为3个步骤:(1)在流场驱动作用下,高分子链漂移并逐渐靠近微通道入口;(2)高分子链逐渐调整自身构象,并使其部分进入微通道;(3)高分子链成功穿过微通道.同时,模拟还发现当高分子链尺寸大于微通道细管道管径时,高分子链穿过微通道所需的平均迁移时间随着流量的增加而逐渐减小.此外,为了研究高分子链刚性对高分子链穿过微通道的影响,模型中还引入了蠕虫状高分子链模型.模拟结果发现,高分子链的链刚性越强,其迁移穿过微通道的时间越长.     

Monte carlo simulation of polymers in steady potential flows

This work presents a study of the configurational properties of bead-spring chains in steady-state potential flows. These properties are obtained from off-lattice Monte Carlo simulations. The results of our simulations compare favourably with theoretical results for phantom chains. We investigate the effects of excluded volume and hydrodynamic interaction on the configurations of the chains, and compare our results to the predictions of scaling arguments. The scaling law observed here for the latter case deviates from theoretical predictions but is in agreement with findings reported from molecular dynamics simulations.     

Dissipative Particle Dynamics Simulations of Polymer Brushes: Comparison with Molecular Dynamics Simulations

Summary: The structure of polymer brushes is investigated by dissipative particle dynamics (DPD) simulations that include explicit solvent particles. With an appropriate choice of the DPD interaction parameters , we obtain good agreement with previous molecular dynamics (MD) results where the good solvent behavior has been modeled by an effective Lennard–Jones potential. The present results confirm that DPD simulation techniques can be applied for large length scale simulations of polymer brushes. A relation between the different length scales and is established.

Polymer brush at a solid–liquid interface.     

Molecular dynamics simulations of aqueous solutions of ethanolamines

We report on molecular dynamics simulations performed at constant temperature and pressure to study ethanolamines as pure components and in aqueous solutions. A new geometric integration algorithm that preserves the correct phase space volume is employed to study molecules having up to three ethanol chains. The most stable geometry, rotational barriers, and atomic charges were obtained by ab initio calculations in the gas phase. The calculated dipole moments agree well with available experimental data. The most stable conformation, due to intramolecular hydrogen bonding interactions, has a ringlike structure in one of the ethanol chains, leading to high molecular stability. All molecular dynamics simulations were performed in the liquid phase. The interaction parameters are the same for the atoms in the ethanol chains, reducing the number of variables in the potential model. Intermolecular hydrogen bonding is also analyzed, and it is shown that water associates at low water mole fractions. The force field reproduced (within 1%) the experimental liquid densities at different temperatures of pure components and aqueous solutions at 313 K. The excess and partial molar volumes are analyzed as a function of ethanolamine concentration.