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

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.     

Combining molecular dynamics with Lattice Boltzmann: a hybrid method for the simulation of (charged) colloidal systems

We present a hybrid method for the simulation of colloidal systems that combines molecular dynamics (MD) with the Lattice Boltzmann (LB) scheme. The LB method is used as a model for the solvent in order to take into account the hydrodynamic mass and momentum transport through the solvent. The colloidal particles are propagated via MD and they are coupled to the LB fluid by viscous forces. With respect to the LB fluid, the colloids are represented by uniformly distributed points on a sphere. Each such point [with a velocity V(r) at any off-lattice position r] is interacting with the neighboring eight LB nodes by a frictional force F = xi0(V(r)-u(r)), with xi0 being a friction coefficient and u(r) being the velocity of the fluid at the position r. Thermal fluctuations are introduced in the framework of fluctuating hydrodynamics. This coupling scheme has been proposed recently for polymer systems by Ahlrichs and Dunweg [J. Chem. Phys. 111, 8225 (1999)]. We investigate several properties of a single colloidal particle in a LB fluid, namely, the effective Stokes friction and long-time tails in the autocorrelation functions for the translational and rotational velocity. Moreover, a charged colloidal system is considered consisting of a macroion, counterions, and coions that are coupled to a LB fluid. We study the behavior of the ions in a constant electric field. In particular, an estimate of the effective charge of the macroion is yielded from the number of counterions that move with the macroion in the direction of the electric field.     

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.     

Dynamics and scaling of polymers in a dilute solution: analytical treatment in two and higher dimensions

We consider the dynamical scaling of a single polymer chain in good solvent. In the case of two-dimensional systems, Shannon and Choy [Phys. Rev. Lett. 79, 1455 (1997)] have suggested that the dynamical scaling for a dilute polymer solution breaks down. Using scaling arguments and analytical calculations based on the Zimm model, we show that the dynamical scaling of a dilute two-dimensional polymer system holds when the relevant dynamical quantities are properly extracted from finite systems. Most important, the polymer diffusion coefficient in two dimensions scales logarithmically with system size, in excellent agreement with our extensive computer simulations. This scaling is the reason for the failure of the previous attempts to resolve the dynamical scaling of dilute two-dimensional polymer systems. In three and higher dimensions our analytic calculations are in agreement with previous results in the literature.     

Rheology of gelling polymers in the Zimm model

In order to study rheological properties of gelling systems in dilute solution, we investigate the viscosity and the normal stresses in the Zimm model [B. H. Zimm, J. Chem. Phys. 24, 269 (1956)]. for randomly cross-linked monomers. The distribution of cluster topologies and sizes is assumed to be given either by Erdos-Renyi random graphs or three-dimensional bond percolation. Within this model the critical behavior of the viscosity and of the first normal stress coefficient is determined by the power-law scaling of their averages over clusters of a given size n with n. We investigate these scaling relations numerically and conclude that the scaling exponents are independent of the hydrodynamic interaction strength. The numerically determined exponents agree well with experimental data for branched polymers. However, we show that this traditional model of polymer physics is not able to yield a critical divergence at the gel point of the viscosity for a polydisperse dilute solution of gelation clusters. A generally accepted scaling relation for the Zimm exponent of the viscosity is thereby disproved.     

One-dimensional models of polymer dynamics. I. Dashpot-spring models applied to a single rotor

The proposion of Clark and Zimm that a dashpot and a spring can be used in place of a set of rotational barriers in polymer dynamics is studied. The simplest possible case is examined here, that of a single rotor. Reasons for altering the Clark–Zimm diffusion equation are presented and an alternative diffusion equation is proposed. The results of both diffusion equations for the correlation function 〈exp[?iθ(0)] exp[iθ(t)]〉 (θ is the angular position of the rotor) are compared with the correlation function for a rotor in an n-fold cosine potential. Although the two diffusion equations differ, both agree well with the n-fold cosine model for barriers above a few k_BT. This agreement is obtained in the absence of adjustable parameters, and motivates the application of these two diffusion equations to polymeric systems.     

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.     

Molecular dynamics simulations of polymer transport in nanocomposites

Molecular dynamics simulations on the Kremer-Grest bead-spring model of polymer melts are used to study the effect of spherical nanoparticles on chain diffusion. We find that chain diffusivity is enhanced relative to its bulk value when polymer-particle interactions are repulsive and is reduced when polymer-particle interactions are strongly attractive. In both cases chain diffusivity assumes its bulk value when the chain center of mass is about one radius of gyration R(g) away from the particle surface. This behavior echoes the behavior of polymer melts confined between two flat surfaces, except in the limit of severe confinement where the surface influence on polymer mobility is more pronounced for flat surfaces. A particularly interesting fact is that, even though chain motion is strongly speeded up in the presence of repulsive boundaries, this effect can be reversed by pinning one isolated monomer onto the surface. This result strongly stresses the importance of properly specifying boundary conditions when the near surface dynamics of chains are studied.     

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.     

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.     

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.     

The single chain limit of structural relaxation in a polyolefin blend

The influence of composition on component dynamics and relevant static properties in a miscible polymer blend is investigated using molecular dynamics simulation. Emphasis is placed on dynamics in the single chain dilution limit, as this limit isolates the role of inherent component mobility in the polymer's dynamic behavior when placed in a blend. For our systems, a biased local concentration affecting dynamics must arise primarily from chain connectivity, which is quantified by the self-concentration, because concentration fluctuations are minimized due to restraints on chain lengths arising from simulation considerations. The polyolefins simulated [poly(ethylene-propylene) (PEP) and poly(ethylene-butene) (PEB)] have similar structures and glass transition temperatures, and all interactions are dispersive in nature. We find that the dependence of dynamics upon composition differs between the two materials. Specifically, PEB (slower component) is more influenced by the environment than PEP. This is linked to a smaller self-concentration for PEB than PEP. We examine the accuracy of the Lodge-McLeish model (which is based on chain connectivity acting over the Kuhn segment length) in predicting simulation results for effective concentration. The model predicts the simulation results with high accuracy when the model's single parameter, the self-concentration, is calculated from simulation data. However, when utilizing the theoretical prediction of the self-concentration the model is not quantitatively accurate. The ability of the model to link the simulated self-concentration with biased local compositions at the Kuhn segment length provides strong support for the claim that chain connectivity is the leading cause of distinct mobility in polymer blends. Additionally, the direct link between the willingness of a polymer to be influenced by the environment and the value of the self-concentration emphasizes the importance of the chain connectivity. Furthermore, these findings are evidence that the Kuhn segment length is the relevant length scale controlling segmental dynamics.     

Influence of molecular topology on the static and dynamic properties of single polymer chain in solution

The influence of molecular topology on the structural and dynamic properties of polymer chain in solution with ring structure, three-arm branched structure, and linear structure are studied by molecular dynamics simulation. At the same degree of polymerization (N), the ring-shaped chain possesses the smallest size and largest diffusion coefficient. With increasing N, the difference of the radii of gyration between the three types of polymer chains increases, whereas the difference of the diffusion coefficients among them decreases. However, the influence of the molecular topology on the static and the dynamic scaling exponents is small. The static scaling exponents decrease slightly, and the dynamic scaling exponents increase slightly, when the topology of the polymer chain is changed from linear to ring-shaped or three-arm branched architecture. The dynamics of these three types of polymer chain in solution is Zimm-like according to the dynamic scaling exponents and the dynamic structure factors.     

高分子动力学的单链模型

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

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.     

Structure and dynamics of short chain molecules in disordered porous materials: a molecular dynamics simulation study

The static and dynamic properties of short polymer chains in disordered materials are studied using discontinuous molecular dynamics simulations. The polymers are modeled as chains of hard spheres and the matrix is a collection of fixed hard spheres. The simulations show that the chain size is a nonmonotonic function of the matrix concentration for all polymer concentrations. The dependence of polymer diffusion D on the degree of polymerization N becomes stronger as the matrix concentration is increased. At high matrix concentrations we observe a decoupling between translational and rotational diffusion, i.e., the rotational relaxation time becomes very large but the translational diffusion is not affected significantly. We attribute this to the trapping of a small number of polymers. Under these conditions the polymer chains diffuse via a hopping mechanism.     

Universal properties of a single polymer chain in slit: Scaling versus molecular dynamics simulations

We revisit the classical problem of a polymer confined in a slit in both of its static and dynamic aspects. We confirm a number of well known scaling predictions and analyze their range of validity by means of comprehensive molecular dynamics simulations using a coarse-grained bead-spring model of a flexible polymer chain. The normal and parallel components of the average end-to-end distance, mean radius of gyration and their distributions, the density profile, the force exerted on the slit walls, and the local bond orientation characteristics are obtained in slits of width D=4/10 (in units of the bead diameter) and for chain lengths N=50/300. We demonstrate that a wide range of static chain properties in normal direction can be described quantitatively by analytic model-independent expressions in perfect agreement with computer experiment. In particular, the observed profile of confinement-induced bond orientation is shown to closely match theory predictions. The anisotropy of confinement is found to be manifested most dramatically in the dynamic behavior of the polymer chain. We examine the relation between characteristic times for translational diffusion and lateral relaxation. It is demonstrated that the scaling predictions for lateral and normal relaxation times are in good agreement with our observations. A novel feature is the observed coupling of normal and lateral modes with two vastly different relaxation times. We show that the impact of grafting on lateral relaxation is equivalent to doubling the chain length.     

Interfacial friction between semiflexible polymers and crystalline surfaces

The results obtained from molecular dynamics simulations of the friction at an interface between polymer melts and weakly attractive crystalline surfaces are reported. We consider a coarse-grained bead-spring model of linear chains with adjustable intrinsic stiffness. The structure and relaxation dynamics of polymer chains near interfaces are quantified by the radius of gyration and decay of the time autocorrelation function of the first normal mode. We found that the friction coefficient at small slip velocities exhibits a distinct maximum which appears due to shear-induced alignment of semiflexible chain segments in contact with solid walls. At large slip velocities, the friction coefficient is independent of the chain stiffness. The data for the friction coefficient and shear viscosity are used to elucidate main trends in the nonlinear shear rate dependence of the slip length. The influence of chain stiffness on the relationship between the friction coefficient and the structure factor in the first fluid layer is discussed.