Direct Molecular Dynamics Simulation of Branch Point Motion in Asymmetric Star Polymer Melts

Capillary filling of star polymer melts in nanopores

The topology of a polymer profoundly influences its behavior. However, its effect on imbibition dynamics remains poorly understood. In the present work, capillary filling (during imbibition and following full imbibition) of star polymer melts was investigated by molecular dynamics simulations with a coarse-grained model. The reversal of imbibition dynamics observed for linear-chain systems was also present for star polymers. Star polymers with short arms penetrate slower than the prediction of the Lucas-Washburn equation, while systems with long arms penetrate faster. The radius of gyration increases during confined flow, indicating the orientation and disentanglement of arms. In addition, the higher the functionality of the star polymer, the more entanglement points are retained. Besides, a stiff region near the core segments of the stars is observed, which increases in size with functionality. The proportion of different configurations of the arms (e.g., loops, trains, tails) changes dramatically with the arm length and degree of confinement but is only influenced by the functionality when the arms are short. Following full imbibition, the different decay rates of the self-correlation function of the core-to-end vector illustrate that arms take a longer time to reach the equilibrium state as the functionality, arm length, and degree of confinement increase, in agreement with recent experimental findings. Furthermore, the star topology induces a stronger effect of adsorption and friction, which becomes more pronounced with increasing functionality.

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.

Dynamics of asymmetric star polymers under coarse grain simulations

The asymmetric star polymers are studied by coarse grain simulations. Each polymer chain is represented by number of consecutive soft blobs and additional uncrossability constraints are added to prevent chain crossings. In this work two types of asymmetric star polymers with different backbone lengths are structured. Their dynamical properties are discussed by comparisons with corresponding linear chains, the one covers chain length along with the asymmetric arm through the branch point to one of the symmetric arm, or the backbone chain between two symmetric arm ends, or the largest linear possesses the same molecular weight of the entire star. To reveal the influence of the asymmetric arm length on their relaxation decay times, the autocorrelation function of the vectors from each branching point to corresponding asymmetric arm end are calculated, results are compared with the symmetric star having the same backbone chain.

Determining Tube Theory Parameters by Slip-Spring Model Simulations of Entangled Star Polymers in Fixed Networks

Dynamical properties of branched polymer melts are determined by the polymer molecular weights and architectures containing junction points. Relaxation of entangled symmetric star polymers proceeds via arm-retraction and constraint release (CR). In this work, we investigate arm-retraction dynamics in the framework of a single-chain slip-spring model without CR effect where entanglements are treated as binary contacts, conveniently modeled as virtual "slip-links", each involving two neighboring strands. The model systems are analogous to isolated star polymers confined in a permanent network or a melt of very long linear polymers. We find that the distributions of the effective primitive path lengths are Gaussian, from which the entanglement molecular weight N e , a key tube theory parameter, can be extracted. The procured N e value is in good agreement with that obtained from mapping the middle monomer mean-square displacements of entangled linear chains in slip-spring model to the tube model prediction. Furthermore, the mean first-passage (FP) times of destruction of original tube segments by the retracting arm end are collected in simulations and examined quantitatively using a theory recently developed in our group for describing FP problems of one-dimensional Rouse chains with improbable extensions. The asymptotic values of N e as obtained from the static (primitive path length) and dynamical (FP time) analysis are consistent with each other. Additionally, we manage to determine the tube survival function of star arms μ ( t ) , or equivalently arm end-to-end vector relaxation function ϕ ( t ) , through the mean FP time spectrum τ ( s ) of the tube segments after careful consideration of the inner-most entanglements, which shows reasonably good agreement with experimental data on dielectric relaxation.

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.

A comparative study of coarse-graining methods for polymeric fluids: Mori-Zwanzig vs. iterative Boltzmann inversion vs. stochastic parametric optimization

We construct effective coarse-grained (CG) models for polymeric fluids by employing two coarse-graining strategies. The first one is a forward-coarse-graining procedure by the Mori-Zwanzig (MZ) projection while the other one applies a reverse-coarse-graining procedure, such as the iterative Boltzmann inversion (IBI) and the stochastic parametric optimization (SPO). More specifically, we perform molecular dynamics (MD) simulations of star polymer melts to provide the atomistic fields to be coarse-grained. Each molecule of a star polymer with internal degrees of freedom is coarsened into a single CG particle and the effective interactions between CG particles can be either evaluated directly from microscopic dynamics based on the MZ formalism, or obtained by the reverse methods, i.e., IBI and SPO. The forward procedure has no free parameters to tune and recovers the MD system faithfully. For the reverse procedure, we find that the parameters in CG models cannot be selected arbitrarily. If the free parameters are properly defined, the reverse CG procedure also yields an accurate effective potential. Moreover, we explain how an aggressive coarse-graining procedure introduces the many-body effect, which makes the pairwise potential invalid for the same system at densities away from the training point. From this work, general guidelines for coarse-graining of polymeric fluids can be drawn.

DNA as a Model for Probing Polymer Entanglements: Circular Polymers and Non-Classical Dynamics

Double-stranded DNA offers a robust platform for investigating fundamental questions regarding the dynamics of entangled polymer solutions. The exceptional monodispersity and multiple naturally occurring topologies of DNA, as well as a wide range of tunable lengths and concentrations that encompass the entanglement regime, enable direct testing of molecular-level entanglement theories and corresponding scaling laws. DNA is also amenable to a wide range of techniques from passive to nonlinear measurements and from single-molecule to bulk macroscopic experiments. Over the past two decades, researchers have developed methods to directly visualize and manipulate single entangled DNA molecules in steady-state and stressed conditions using fluorescence microscopy, particle tracking and optical tweezers. Developments in microfluidics, microrheology and bulk rheology have also enabled characterization of the viscoelastic response of entangled DNA from molecular levels to macroscopic scales and over timescales that span from linear to nonlinear regimes. Experiments using DNA have uniquely elucidated the debated entanglement properties of circular polymers and blends of linear and circular polymers. Experiments have also revealed important lengthscale and timescale dependent entanglement dynamics not predicted by classical tube models, both validating and refuting new proposed extensions and alternatives to tube theory and motivating further theoretical work to describe the rich dynamics exhibited in entangled polymer systems.

Communication: Role of short chain branching in polymer structure and dynamics

A comprehensive understanding of chain-branching effects, essential for establishing general knowledge of the structure-property-phenomenon relationship in polymer science, has not yet been found, due to a critical lack of knowledge on the role of short-chain branches, the effects of which have mostly been neglected in favor of the standard entropic-based concepts of long polymers. Here, we show a significant effect of short-chain branching on the structural and dynamical properties of polymeric materials, and reveal the molecular origins behind the fundamental role of short branches, via atomistic nonequilibrium molecular dynamics and mesoscopic Brownian dynamics by systematically varying the strength of the mobility of short branches. We demonstrate that the fast random Brownian kinetics inherent to short branches plays a key role in governing the overall structure and dynamics of polymers, leading to a compact molecular structure and, under external fields, to a lesser degree of structural deformation of polymer, to a reduced shear-thinning behavior, and to a smaller elastic stress, compared with their linear analogues. Their fast dynamical nature being unaffected by practical flow fields owing to their very short characteristic time scale, short branches would substantially influence (i.e., facilitate) the overall relaxation behavior of polymeric materials under various flowing conditions.

Incorporation of memory effects in coarse-grained modeling via the Mori-Zwanzig formalism

The Mori-Zwanzig formalism for coarse-graining a complex dynamical system typically introduces memory effects. The Markovian assumption of delta-correlated fluctuating forces is often employed to simplify the formulation of coarse-grained (CG) models and numerical implementations. However, when the time scales of a system are not clearly separated, the memory effects become strong and the Markovian assumption becomes inaccurate. To this end, we incorporate memory effects into CG modeling by preserving non-Markovian interactions between CG variables, and the memory kernel is evaluated directly from microscopic dynamics. For a specific example, molecular dynamics (MD) simulations of star polymer melts are performed while the corresponding CG system is defined by grouping many bonded atoms into single clusters. Then, the effective interactions between CG clusters as well as the memory kernel are obtained from the MD simulations. The constructed CG force field with a memory kernel leads to a non-Markovian dissipative particle dynamics (NM-DPD). Quantitative comparisons between the CG models with Markovian and non-Markovian approximations indicate that including the memory effects using NM-DPD yields similar results as the Markovian-based DPD if the system has clear time scale separation. However, for systems with small separation of time scales, NM-DPD can reproduce correct short-time properties that are related to how the system responds to high-frequency disturbances, which cannot be captured by the Markovian-based DPD model.

Perspectives on the viscoelasticity and flow behavior of entangled linear and branched polymers

We briefly review the recent advances in the rheology of entangled polymers and identify emerging research trends and outstanding challenges, especially with respect to branched polymers. Emphasis is placed on the role of well-characterized model systems, as well as the synergy of synthesis-characterization, rheometry and modeling/simulations. The theoretical framework for understanding the observed linear and nonlinear rheological phenomena is the tube model, which is critically assessed in view of its successes and shortcomings, and alternative approaches are briefly discussed. Finally, intriguing experimental findings and controversial issues that merit consistent explanation, such as shear banding instabilities, multiple stress overshoots in transient simple shear and enhanced steady-state elongational viscosity in polymer solutions, are discussed, and future directions such as branch point dynamics and anisotropic monomeric friction are outlined.

Construction of dissipative particle dynamics models for complex fluids via the Mori-Zwanzig formulation

We present a bottom-up coarse-graining procedure to construct mesoscopic force fields directly from microscopic dynamics. By grouping many bonded atoms in the molecular dynamics (MD) system into a single cluster, we compute both the conservative and non-conservative interactions between neighboring clusters. In particular, we perform MD simulations of polymer melts to provide microscopic trajectories for evaluating coarse-grained (CG) interactions. Subsequently, dissipative particle dynamics (DPD) is considered as the effective dynamics resulting from the Mori-Zwanzig (MZ) projection of the underlying atomistic dynamics. The forces between finite-size clusters have, in general, both radial and transverse components and hence we employ four different DPD models to account differently for such interactions. Quantitative comparisons between these DPD models indicate that the DPD models with MZ-guided force fields yield much better static and dynamics properties, which are consistent with the underlying MD system, compared to standard DPD with empirical formulae. When the rotational motion of the particle is properly taken into account, the entire velocity autocorrelation function of the MD system as well as the pair correlation function can be accurately reproduced by the MD-informed DPD model. Since this coarse-graining procedure is performed on an unconstrained MD system, our framework is general and can be used in other soft matter systems in which the clusters can be faithfully defined as CG particles.

Analytical Rheology of Polymer Melts: State of the Art

The extreme sensitivity of rheology to the microstructure of polymer melts has prompted the development of “analytical rheology,” which seeks inferring the structure and composition of an unknown sample based on rheological measurements. Typically, this involves the inversion of a model, which may be mathematical, computational, or completely empirical. Despite the imperfect state of existing models, analytical rheology remains a practically useful enterprise. I review its successes and failures in inferring the molecular weight distribution of linear polymers and the branching content in branched polymers.

In Silico Molecular Design, Synthesis, Characterization, and Rheology of Dendritically Branched Polymers: Closing the Design Loop

It has been a long held ambition of both industry and academia to understand the relationship between the often complex molecular architecture of polymer chains and their melt flow properties, with the goal of building robust theoretical models to predict their rheology. The established key to this is the use of well-defined, model polymers, homogeneous in chain length and architecture. We describe here for the first time, the in silico design, synthesis, and characterization of an architecturally complex, branched polymer with the optimal rheological properties for such structure-property correlation studies. Moreover, we demonstrate unequivocally the need for accurate characterization using temperature gradient interaction chromatography (TGIC), which reveals the presence of heterogeneities in the molecular structure that are undetectable by size exclusion chromatography (SEC). Experimental rheology exposes the rich pattern of relaxation dynamics associated with branched polymers, but the ultimate test is, of course, did the theoretical (design) model accurately predict the rheological properties of the synthesized model branched polymer? Rarely, if ever before, has such a combination of theory, synthesis, characterization, and analysis resulted in a "yes", expressed without doubt or qualification.

Synthesis, characterizations and biocompatibility of novel biodegradable star block copolymers based on poly[(R)-3-hydroxybutyrate] and poly(epsilon-caprolactone)

Star block copolymers based on poly[(R)-3-hydroxybutyrate] (PHB) and poly(epsilon-caprolactone) (PCL), termed SPHBCL, were successfully synthesized with structural variation on arm numbers and lengths via coupling reactions and ring opening polymerizations. Arm numbers 3, 4 and 6 of SPHBCL were synthesized by using different multifunctional cores, such as trimethyol propane, pentaerythiritol and dipentaerthritol, respectively. Gel permeation chromatography (GPC) and (1)H and (13)C nuclear magnetic resonance were used to characterize the structure of SPHBCL. GPC failed to produce accurate molecular weights of the SPHBCL due to the discrepancy of star copolymer structures. The melting temperature of SPHBCL decreased with increasing degree of branching. Thermal decomposition temperature was revealed to be lower than that of linear block copolymer LPHBCL counterparts based on PHB and PCL. Films made from various SPHBCL copolymers had different porous or networking surface morphology, and all possessed improved biocompatibility in terms of less blood clotting and more osteoblast cell growth compared with their corresponding homopolymers PHB and PCL. Among them, it was found, however, that the 4-arm star block copolymer 4SPHBCL-25 showed unique surface properties, i.e. a regular nanoravine structure was observed by scanning electron microscopy and atomic force microscopy. This 4-arm star copolymer also showed the best biocompatibility.

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.