Search

We report a multi-chain approach for dissipative particle dynamics where the uncrossability constraints of polymer chains are mimicked by temporary cross-links, so-called slip-springs. The conformational statistics of the chains are not affected by the introduction of slip-springs. Dynamical properties such as mean square displacements, diffusion coefficient, and longest relaxation time are in good agreement with the results of reptation theory. According to our analysis, the present formalism is 500 times faster and requires 7 times fewer beads than conventional generic polymer models employing Newtonian dynamics and excluded-volume potentials. [ABSTRACT FROM AUTHOR]

For branched polymers, the curvilinear motion of the branch point along the backbone is a significant relaxation source but details of this motion have not been well understood. This study conducts multi-chain sliplink simulations to examine effects of the spatial fluctuation and curvilinear hopping of the branch point on the viscoelastic relaxation. The simulation is based on the primitive chain network model that allows the spatial fluctuations of sliplink and branch point and the chain sliding along the backbone according to the subchain tension, chemical potential gradients, drag force against medium, and random force. The sliplinks are created and/or disrupted through the motion of chain ends. The curvilinear hopping of the branch point along the backbone is allowed to occur when all sliplinks on a branched arm are lost. The simulations considering the fluctuation and the hopping of the branch point described well the viscoelastic data for symmetric and asymmetric star polymers with a parameter set common to the linear polymer. On the other hand, the simulations without the branch point motion predicted unreasonably slow relaxation for asymmetric star polymers. For asymmetric star polymers, further tests with and without the branch point hopping revealed that the hopping is much less important compared to the branch point fluctuation when the lengths of the short and long backbone arms are not very different and the waiting time for the branch point hopping (time for removal of all sliplinks on the short arm) is larger than the backbone relaxation time. Although this waiting time changes with the hopping condition, the above results suggest a significance of the branch point fluctuation in the actual relaxation of branch polymers. [ABSTRACT FROM AUTHOR]

Doi and Edwards (DE) proposed that the relaxation of entangled linear polymers under large deformation occurs in two steps: the fast chain contraction (via the longitudinal Rouse mode of the chain backbone) and the slow orientational relaxation (due to reptation). The DE model assumes these relaxation processes to be independent and decoupled. However, this decoupling is invalid for a generalized convective constraint release (CCR) mechanism that releases the entanglement on every occasion of the contraction of surrounding chains. Indeed, the decoupling does not occur in the sliplink models where the entanglement is represented by the binary interaction (hooking) of chains. Thus, we conducted primitive chain network simulations based on a multichain sliplink model to investigate the chain contraction under step shear. The simulation quantitatively reproduced experimental features of the nonlinear relaxation modulus G(t,γ). Namely, G(t,γ) was cast in the time-strain separable form, G(t,γ)=h(γ)G(t) with h(γ)=damping function and G(t)=linear modulus, but this rigorous separability was valid only at times t comparable to the terminal relaxation time, although a deviation from this form was rather small (within ±10%) at t>τR (longest Rouse relaxation time). A molecular origin of this delicate failure of time-strain separability at t∼τR was examined for the chain contour length, subchain length, and subchain stretch. These quantities were found to relax in three steps, the fast, intermediate, and terminal steps, governed by the local force balance between the subchains, the longitudinal Rouse relaxation, and the reptation, respectively. The contributions of the terminal reptative mode to the chain length relaxation as well as the subchain length/stretch relaxation, not considered in the original DE model, emerged because the sliplinks (entanglement) were removed via the generalized CCR mechanism explained above and the reformation of the sliplinks was slow at around the chain center compared to the more rapidly fluctuating chain end. The number of monomers in the subchain were kept larger at the chain center than at the chain end because of the slow entanglement reformation at the center, thereby reducing the tension of the stretched subchain at the chain center compared to the DE prediction. This reduction of the tension at the chain center prevented completion of the length equilibration of subchains at t∼τR (which contradicts to the DE prediction), and it forces the equilibration to complete through the reptative mode at t>τR. The delicate failure of time-strain separability seen for G(t,γ) at t∼τR reflects this retarded length equilibration. [ABSTRACT FROM AUTHOR]

The primitive chain network (PCN) model successfully employed to simulate the rheology of entangled polymers is here tested versus less coarse-grained (lattice or atomistic) models for what concerns the structure of the network at equilibrium (i.e., in the absence of flow). By network structure, we mean the distributions of some relevant quantities such as subchain length in space or in monomer number. Indeed, lattice and atomistic simulations are obviously more accurate, but are also more difficult to use in nonequilibrium flow situations, especially for long entangled polymers. Conversely, the coarse-grained PCN model that deals more easily with rheology lacks, strictly speaking, a rigorous foundation. It is therefore important to verify whether or not the equilibrium structure of the network predicted by the PCN model is consistent with the results recently obtained by using lattice and atomistic simulations. In this work, we focus on single chain properties of the entangled network. Considering the significant differences in modeling the polymer molecules, the results here obtained appear encouraging, thus providing a more solid foundation to Brownian simulations based on the PCN model. Comparison with the existing theories also proves favorable. [ABSTRACT FROM AUTHOR]

Molecular theories for polymer rheology are based on conformational dynamics of the polymeric chain. Hence, measurements directly related to molecular conformations appear more appealing than indirect ones obtained from rheology. In this study, primitive chain network simulations are compared to experimental data of entangled DNA solutions [Teixeira et al., Macromolecules 40, 2461 (2007)]. In addition to rheological comparisons of both linear and nonlinear viscoelasticities, a molecular extension measure obtained by Teixeira et al. through fluorescent microscopy is compared to simulations, in terms of both averages and distributions. The influence of flow on conformational distributions has never been simulated for the case of entangled polymers, and how DNA molecular individualism extends to the entangled regime is not known. The linear viscoelastic response and the viscosity growth curve in the nonlinear regime are found in good agreement with data for various DNA concentrations. Conversely, the molecular extension measure shows significant departures, even under equilibrium conditions. The reason for such discrepancies remains unknown. [ABSTRACT FROM AUTHOR]

In condensed polymeric liquids confined in slit channels, the movement of chains is constrained by two factors: entanglement among the chains and the excluded volume between the chains and the wall. In this study, we propose a wall boundary (WB) model for the primitive chain network (PCN) model, which describes the dynamics of polymer chains in bulk based on coarse graining upon the characteristic molecular weight of the entanglement. The proposed WB model is based on the assumptions that (i) polymers are not stuck but simply reflected randomly by the wall, and (ii) subchains below the entanglement length scale behave like those in bulk even near the wall. Using the WB model, we simulate the dynamics of entangled polymer chains confined in slit channels. The results show that as the slit narrows, the chains are compressed in the direction normal to the wall, while they are expanded in the parallel direction. In addition, the relaxation time of the end-to-end vector increases, and the diffusivity of the center of mass decreases. The compression in the normal direction is a natural effect of confinement, while the expansion is introduced by a hooking process near the wall. The trends revealed that the relaxation time and diffusivity depend on the increase in friction due to an increased number of entanglements near the wall, which is also associated with the hooking process in the PCN model. These results are expected within the assumptions of the PCN model. Thus, the proposed WB model can successfully reproduce the effects of wall confinement on chains. [ABSTRACT FROM AUTHOR]

Simulation results of the primitive chain network (PCN) model for entangled polymers are compared here to existing data of diffusion coefficient, linear and nonlinear shear and elongational rheology of monodisperse polystyrene melts. Since the plateau modulus of polystyrene is well known from the literature, the quantitative comparison between the whole set of data and simulations only requires a single adjustable parameter, namely, a basic time. The latter, however, must be consistent with the known rheology of unentangled polystyrene melts, i.e., with Rouse behavior, and is therefore not really an adjustable parameter. The PCN model adopted here is a refined version of the original model, incorporating among other things a more accurate description of chain end dynamics as well as finite extensibility effects. In the new version, we find good agreement with linear rheology, virtually without adjustable parameters. It is also shown that, at equilibrium, Gaussian statistics are well obeyed in the simulated network. In the nonlinear range, excellent agreement with data is found in shear, whereas discrepancies and possible inadequacies of the model emerge in fast uniaxial elongational flows, even when accounting for finite extensibility of the network strands. [ABSTRACT FROM AUTHOR]

The differences in formation and structural properties of polymer networks consisting of end-linked flexible or rigid chains were studied by molecular dynamics simulation. Networks were formed from monodisperse, linear, short, flexible or rigid chains with functional end groups and a stoichiometric ratio of trifunctional cross-linkers. The rigid chains had a rodlike shape defined by an angle potential, while the flexible chains had no angle potential. In order to understand the influence of chain rigidity, all parameters of precursor chains (length, reactivity, bond potential, nonbonding potential) were the same, with the exception of the angle potential. The system density ρ, corresponding to the concentration of monomer in solvent, was varied from 0.01 to 0.11. Different network structures resulting from the different processes of network formation were observed. Simulations showed that the flexible chains created an inhomogeneous network on a large scale via microgel cluster formation, in agreement with experimental observations, whereas the rigid chains rapidly created a homogeneous network in the entire system volume without first generating microgel clusters, with the additional difference that they gave rise to mutually interpenetrating networks at the local scale. [ABSTRACT FROM AUTHOR]

More than 100 μs dissipative particle dynamics simulations were carried out to investigate the spontaneous formation process of threadlike micelles from the random configuration for surfactant molecules. Stable spherical micelles were formed during the earlier stage. These spherical micelles fused to each other and grew into rodlike and threadlike micelles during the later stage. The length and radius of a micelle were estimated by tracing the backbone positions and the distance between the head group particles and the backbone of the micelles, respectively. The ratio of the largest to the smallest principal moments of inertia for each micelle was calculated as the micelle shape. [ABSTRACT FROM AUTHOR]

The concept of dynamic tube dilation (DTD) is here used to formulate a new simulation scheme to obtain the linear viscoelastic response of long chains with a large number of entanglements. The new scheme is based on the primitive chain network model previously proposed by some of the authors, and successfully employed to simulate linear and nonlinear behavior of moderately entangled polymers. Scaling laws are generated by the DTD concept, and allow for prediction of the linear response of very long chains on the basis of suitable simulations performed on shorter ones, without introducing adjustable parameters. Tests of the method against existing data for linear monodisperse polyisoprene and polystyrene show good quantitative agreement.© 2004 American Institute of Physics. [ABSTRACT FROM AUTHOR]

A new model for Brownian dynamics simulations of entangled polymeric liquids is proposed here. Chains are coarse grained at the level of segments between consecutive entanglements; hence, the system is in fact a network of primitive chains. The model incorporates not only the “individual” mechanisms of reptation and tube length fluctuation, but also collective contributions arising from the 3D network structure of the entangled system, such as constraint release. Chain coupling is achieved by fulfilling force balance on the entanglement nodes. The Langevin equation for the nodes contains both the tension in the chain segments emanating from the node and an osmotic force arising from density fluctuations. Entanglements are modeled as slip links, each connecting two chain strands. The motion of monomers through slip links, which ultimately generates reptation as well as tube length fluctuations, is also described by a suitable Langevin equation. Creation and release of entanglements is controlled by the number of monomers at the chain ends. In a creation event, the partner chain segment is chosen randomly among those spatially close to the advancing chain end. To validate the model, equilibrium dynamics simulations were run for monodisperse linear chains containing up to Z=40 entanglements. The results show, in agreement with experiments, (i) a Z[sup 3.5±0.1] dependence of the longest relaxation time, (ii) a Z[sup -2.4±0.2] dependence of the self-diffusion coefficient, and (iii) a relaxation modulus proportional to the square of the end-to-end vector correlation function, consistently with the dynamic tube dilation concept. © 2001 American Institute of Physics. [ABSTRACT FROM AUTHOR]