Simulation studies on architecture dependence of unentangled polymer melts

Non-Rouse behavior of short ring polymers in melts by molecular dynamics simulations

Short ring polymers are expected to behave nearly Rouse-like due to the little effect of topological constraints of non-knot and non-concatenation. However, this notion is questioned because of several simulation and experiment findings in recent times, which requires a further more quantitative study. Therefore, we perform a deep investigation of statics and dynamics of flexible short ring polymers (N < 2Ne) in melts via molecular dynamics simulations by further taking linear analogues as well as all-crossing ring and linear polymers with switched off topological constraints for comparisons and demonstrate the noticeable deviations from the Rouse model in terms of local and global scales. Although the overall size is compact, the subchains are swollen, which is traced back to the deeper "segmental correlation hole" effect. The same scaling relationship of the non-Gaussian deviation of the static structure factor holds, but the deviation magnitude of rings is larger than that of linear analogues. By checking the non-Gaussian parameter and autocorrelation function of center-of-mass velocity, the physical origin of anomalous sub-diffusions of short rings is identified as unscreened viscoelastic hydrodynamic interactions and not correlation hole effects, like linear analogues.

Distinguishing between Linear and Star Polystyrenes with Unentangled Arms by Dynamic Oscillatory Shear Tests

Linear and nonlinear viscoelastic properties of star polystyrene (PS) melts with unentangled arms were measured by using small-amplitude and medium-amplitude oscillatory shear (SAOS and MAOS) tests. For comparison purposes, such tests were also conducted on entangled linear and star PS melts. Interestingly, the linear viscoelastic properties of unentangled star PS were quantitatively described using the Lihktman-McLeish model for entangled linear chains, indicating that unentangled stars were indistinguishable from linear chains by using relaxation spectra. By contrast, the relative intrinsic nonlinearity (Q0), one of the MAOS material functions, exhibited a difference between unentangled star and linear PS. When the maximum Q0 value (Q0,max) was plotted against the entanglement number of span molecules (Zs), unentangled star PS exhibited larger Q0,max values than linear PS, which was quantitatively predicted via the multimode K-BKZ model. Therefore, in the unentangled regime, star PS was concluded to be characterized by intrinsically higher relative nonlinearity than linear PS.

Slip-Spring Hybrid Particle-Field Molecular Dynamics for Coarse-Graining Branched Polymer Melts: Polystyrene Melts as an Example

The topology of chains significantly modifies the dynamical properties of polymer melts. Here, we extend a recently developed efficient simulation method, namely the slip-spring hybrid particle-field (SS-hPF) model, to study the structural and dynamical properties of branched polymer melts over large spatial-temporal scales. In the coarse-grained SS-hPF simulation of polymers, the bonded potentials are derived by iterative Boltzmann inversion from the underlying fine-grained model. The nonbonded potentials are computed from a density functional field instead of pairwise interactions used in standard molecular dynamics simulations, which increases the computational efficiency by a factor of 10-20. The entangled dynamics is lost due to the soft-core nature of density functional field interactions. It is recovered by a multichain slip-spring model that is rigorously parametrized from existing experimental or simulation data. To quantitatively predict the relaxation and diffusion of branched polymers, which are dominated by arm retraction rather than chain reptation, the slip-spring algorithm is augmented to improve the polymer dynamics near the branch point. Multiple dynamical observables, e.g., diffusion coefficients, arm relaxations, and tube survival probabilities, are characterized in an example coarse-grained model of symmetric and asymmetric star-shaped polystyrene melts. Consistent dynamical behaviors are identified and compared with theoretical predictions. With a single rescaling factor, the prediction of diffusion coefficients agrees well with the available experimental measurements. In this work, an efficient approach is provided to build chemistry-specific coarse-grained models for predicting the dynamics of branched polymers.

Effects of Slip-Spring Parameters and Rouse Bead Density on Polymer Dynamics in Multichain Slip-Spring Simulations

The multichain slip-spring (MCSS) model is one of the coarse-grained models of polymers developed in the niche between bead-spring models and tube type descriptions. In this model, polymers are represented by Rouse chains connected by virtual springs that temporally connect the chains, hop along the chain, and are constructed and annihilated at the chain ends. Earlier studies have shown that MCSS simulations can nicely reproduce entangled and unentangled polymer dynamics. However, the model parameters have been chosen arbitrarily, and their effects have not been reported. In this study, for the first time, we systematically investigated the effects of model parameters: fugacity of virtual springs, its intensity, and the Rouse bead density. We validated the employed simulation code by confirming that the statistics of the system follow the theoretical setup. Namely, the virtual spring density is correctly controlled, and polymer chains exhibit ideal chain statistics irrespective of the chosen parameter values. For diffusion and linear viscoelasticity, simulation results obtained for different parameters can be superposed with each other by conversion factors for the bead number per chain and units of length, time, and modulus. These conversion factors follow scaling laws concerning the number of Rouse segments between two consecutive anchoring points of virtual springs along the polymer chain. Besides, diffusion and viscoelasticity excellently agree with literature data for the standard bead-spring simulation. These results imply that the coarse-graining level for the MCSS model can be arbitrarily chosen and controlled by model parameters.

Prediction of crossover in the molecular weight dependence of polyethylene viscosity using a polymer free volume theory

Based upon the Doolittle concept that viscosity and free volume are inversely related, we used the Boltzmann equation along with a polymer free volume theory to calculate the viscosity (η) of polyethylene with three different molecular structures - linear, ring and four-arm symmetrical star - over a molecular weight (M) range of 420-14 000 g mol-1. Free volume parameters were estimated using the Polymer Reference Interaction Site Model (PRISM) and generic van der Waals (GvdW) equation. The polymer free volume theory was able to describe the crossovers in the molecular weight dependence of the viscosity of the aforementioned molecular structures. In particular, the crossover for the linear structure was predicted at about 3000 g mole-1 with η ∼ M1.5 in the unentangled regime and η ∼ M3.3 in the entangled regime that agree reasonably well with experiment. The predictions for the other two structures also agree with the available simulation data. We also demonstrated that the accuracy of the viscosity prediction was sensitive to the difference between two free volume parameters (i.e., (φ+ - F)). Here, F signifies the probability of a bead finding free volume greater than the critical free volume while the fraction of such beads (φ+) can be used to calculate the activation energy. The theory also reproduced the temperature dependence of η for the linear structure at different M, giving apparent activation energy (Eappa) values in the range of 5.30-7.70 kcal mole-1 that are in good agreement with experimental values of 5.50-6.75 kcal mole-1. This work demonstrates for the first time that viscosity of polymer melts can be determined from the polymer free volume theory.

Influence of Branching on the Configurational and Dynamical Properties of Entangled Polymer Melts

We probe the influence of branching on the configurational, packing, and density correlation function properties of polymer melts of linear and star polymers, with emphasis on molecular masses larger than the entanglement molecular mass of linear chains. In particular, we calculate the conformational properties of these polymers, such as the hydrodynamic radius R h , packing length p, pair correlation function g ( r ) , and polymer center of mass self-diffusion coefficient, D, with the use of coarse-grained molecular dynamics simulations. Our simulation results reproduce the phenomenology of simulated linear and branched polymers, and we attempt to understand our observations based on a combination of hydrodynamic and thermodynamic modeling. We introduce a model of "entanglement" phenomenon in high molecular mass polymers that assumes polymers can viewed in a coarse-grained sense as "soft" particles and, correspondingly, we model the emergence of heterogeneous dynamics in polymeric glass-forming liquids to occur in a fashion similar to glass-forming liquids in which the molecules have soft repulsive interactions. Based on this novel perspective of polymer melt dynamics, we propose a functional form for D that can describe our simulation results for both star and linear polymers, covering both the unentangled to entangled polymer melt regimes.

Modeling of Entangled Polymer Diffusion in Melts and Nanocomposites: A Review

This review concerns modeling studies of the fundamental problem of entangled (reptational) homopolymer diffusion in melts and nanocomposite materials in comparison to experiments. In polymer melts, the developed united atom and multibead spring models predict an exponent of the molecular weight dependence to the polymer diffusion very similar to experiments and the tube reptation model. There are rather unexplored parameters that can influence polymer diffusion such as polymer semiflexibility or polydispersity, leading to a different exponent. Models with soft potentials or slip-springs can estimate accurately the tube model predictions in polymer melts enabling us to reach larger length scales and simulate well entangled polymers. However, in polymer nanocomposites, reptational polymer diffusion is more complicated due to nanoparticle fillers size, loading, geometry and polymer-nanoparticle interactions.

Individual circular polyelectrolytes under shear flow

Individual circular polyelectrolytes in simple shear flow are studied by means of mesoscale hydrodynamic simulations, revealing the complex coupling effects of shear rate, electrostatic interaction, and circular architecture on their conformational and dynamical properties. Shear flow deforms the polyelectrolyte and strips condensed counterions from its backbone. A decrease in condensed counterions alters electrostatic interactions among charged particles, affecting shear-induced polymer deformation and orientation. Circular architecture determines the features of deformation and orientation. At weak electrostatic interaction strengths, the polyelectrolyte changes its shape from an oblate ring at small shear rates to a prolate ring at large shear rates, whereas strong electrostatic interaction strengths are associated with a transition from a prolate coil to a prolate ring. Circular polyelectrolytes exhibit tumbling and tank-treading motions in the range of large shear rates. Further study reveals a similarity between the roles of intramolecular electrostatic repulsion and chain rigidity in shear-induced dynamics.

Comparison among multi-chain models for entangled polymer dynamics

Although lots of coarse-grained models have been proposed to trace the long-term behaviors of entangled polymers, compatibility among the different models has not been frequently discussed. In this study, some dynamical and static quantities, such as diffusion, relaxation modulus, chain dimension, and entanglement density, were examined for the multi-chain slip-link model (primitive chain network model) and the multi-chain slip-spring model, and the results were compared with those reported for the standard bead-spring model. For the diffusion, three models are compatible with scale-conversion parameters for units of length, time and bead (segment) number (or the molecular weight). The relaxation modulus is also compatible given that the model dependence can be accommodated by the entanglement density and the additional scale-conversion for the unit of modulus. The chain dimension is reasonably coincident with small deviations due to the weak non-Gaussianity of the models. Apart from these plausible compatibilities, significant discrepancies have been found for the inter-chain cross-correlations in the relaxation modulus.

Structural Mechanism for Viscosity of Semiflexible Polymer Melts in Shear Flow

The viscosities of semiflexible polymers with different chain stiffnesses in shear flow are studied via nonequilibrium molecular dynamics techniques. The simulation reproduces the experimentally observed results, giving a complete picture of viscosity as chain stiffness increases. Analysis of flow-induced changes in chain conformation and local structure indicates two distinct mechanisms behind a variety of viscosity curves. For polymers of small stiffnesses, it is related to flow-induced changes in chain conformation and, for those of large stiffnesses, to flow-induced instabilities of nematic structures. The four-region flow curve is confirmed for polymers of contour length close to persistence length and understood by combining the two structural mechanisms. Thus, these findings clarify the microscopic structures indicated by the macroscopic viscosity.