Simulating startup shear of entangled polymer melts

Start-up shear rheology is a standard experiment used for characterizing polymer flow, and to test various models of polymer dynamics. A rich phenomenology is developed for behavior of entangled monodisperse linear polymers in such tests, documenting shear stress overshoots as a function of shear rates and molecular weights. A tube theory does a reasonable qualitative job at describing these phenomena, although it involves several drastic approximations and the agreement can be fortuitous. Recently, Lu and coworkers published several papers [e.g. Lu {\it et al.} {\it ACS Macro Lett}. 2014, 3, 569-573] reporting results from molecular dynamics simulations of linear entangled polymers, which contradict both theory and experiment. Based on these observations, they made very serious conclusions about the tube theory, which seem to be premature. In this letter, we repeat simulations of Lu {\it et al.} and systematically show that neither their simulation results, nor their comparison with theory are confirmed.

Interactions of an anionic surfactant with poly(oxyalkylene) copolymers in aqueous solution

The interactions of sodium dodecyl sulfate (SDS) with poly(ethylene oxide)/poly(alkylene oxide) (E/A) block copolymers are explored in this study: With respect to the specific compositional characteristics of the copolymer, introduction of SDS can induce fundamentally different effects to the self-assembly behavior of E/A copolymer solutions. In the case of the E18B10-SDS system (E = poly(ethylene oxide) and B = poly(butylene oxide)) development of large surfactant-polymer aggregates was observed. In the case of B20E610-SDS, B12E227B12-SDS, E40B10E40-SDS, E19P43E19-SDS (P = poly(propylene oxide)), the formation of smaller particles compared to pure polymeric micelles points to micellar suppression induced by the ionic surfactant. This effect can be ascribed to a physical binding between the hydrophobic block of unassociated macromolecules and the non-polar tail of the surfactant. Analysis of critical micelle concentrations (cmc*) of polymer-surfactant aqueous solutions within the framework of regular solution theory for binary surfactants revealed negative deviations from ideal behavior for E40B10E40-SDS and E19P43E19-SDS, but positive deviations for E18B10-SDS. Ultrasonic studies performed for the E19P43E19-SDS system enabled the identification of three distinct regions, corresponding to three main steps of the complexation; SDS absorption to the hydrophobic backbone of polymer, development of polymer-surfactant complexes and gradual breakdown of the mixed aggregates. (C) 2008 Elsevier Inc. All rights reserved.

Chapter 16: Electrospinning of polymer solutions

Electrospinning is a technique that involves the production of nanoscale to microscale sized polymer fibres through the application of an electric field to a droplet of polymer solution passed through a spinneret tip. This chapter considers the optimisisation of the electrospinning process and in particular the variation with solution concentration. We show the strong connection between overlapping chains and the successful spinning of fibres. We use small-angle neutron scattering to evaluate the molecular conformations in the solutions and in the fibres.

Shear banding in molecular dynamics of polymer melts

In order to establish constitutive equations for a viscoelastic fluid uniform shear flow is usually required. However, in the last 10 years S. Q. Wang and co-workers have demonstrated that some entangled polymers do not flow with the uniform shear rate as usually assumed, but instead choose to separate into fast and slow flowing regions. This phenomenon, known as shear banding, causes flow instabilities and in principle invalidates all rheological measurements when it occurs. In this Letter we report the first observation of shear banding in molecular dynamics simulations of entangled polymer melts. We show that our observations are in a very good agreement with the phenomenology developed by Fielding and Olmsted. Our findings provide a simple way of validating the empirical macroscopic phenomenology of shear banding. © 2012 American Physical Society

Modelling small angle neutron scattering data from electrospun fibres

Electrospinning is a technique employed to produce nanoscale to microscale sized fibres by the application of a high voltage to a spinneret containing a polymer solution. Here we examine how small angle neutron scattering data can be modelled to analyse the polymer chain conformation. We prepared 1:1 blends of deuterated and hydrogenated atactic-polystyrene fibres from solutions in N, N-Dimethylformamide and Methyl Ethyl Ketone. The fibres themselves often contain pores or voiding within the internal structure on the length scales that can interfere with scattering experiments. A model to fit the scattering data in order to obtain values for the radius of gyration of the polymer molecules within the fibres has been developed, that includes in the scattering from the voids. Using this model we find that the radius of gyration is 20% larger than in the bulk state and the chains are slightly extended parallel to the fibre axis.

Effect of water-soluble polymers, polyethylene glycol and poly(vinylpyrrolidone),on the gelation of aqueous micellar solutions of Pluronic copolymer F127

The micellization of F127 (E98P67E98) in dilute aqueous solutions of polyethylene glycol (PEG6000 and PEG35000) and poly(vinylpyrrolidone) (PVP K30 and PVP K90) is studied. The average hydrodynamic radius (rh,app) obtained from the dynamic light scattering technique increased with increase in PEG concentration but decreased on addition of PVP, results which are consistent with interaction of the micelles with PEG and the formation of micelles clusters, but no such interaction occurs with PVP. Tube inversion was used to determine the onset of gelation. The critical concentration of F127 for gelation increased on addition of PEG and of PVP K30 but decreased on addition of PVP K90. Small-angle X-ray scattering (SAXS) was used to show that the 30 wt% F127 gel structure (fcc) was independent of polymer type and concentration, as was the d-spacing and so the micelle hard-sphere radius. The maximum elastic modulus (G0 max) of 30 wt% F127 decreased from its value for water alone as PEG was added, but was little changed by adding PVP. These results are consistent with the packed-micelles in the 30 wt% F127 gel being effectively isolated from the polymer solution on the microscale while, especially for the PEG, being mixed on the macroscale.

Segmental dynamics in entangled linear polymer melts

We present molecular dynamics (MD) and slip-springs model simulations of the chain segmental dynamics in entangled linear polymer melts. The time-dependent behavior of the segmental orientation autocorrelation functions and mean-square segmental displacements are analyzed for both flexible and semiflexible chains, with particular attention paid to the scaling relations among these dynamic quantities. Effective combination of the two simulation methods at different coarse-graining levels allows us to explore the chain dynamics for chain lengths ranging from Z ≈ 2 to 90 entanglements. For a given chain length of Z ≈ 15, the time scales accessed span for more than 10 decades, covering all of the interesting relaxation regimes. The obtained time dependence of the monomer mean square displacements, g1(t), is in good agreement with the tube theory predictions. Results on the first- and second-order segmental orientation autocorrelation functions, C1(t) and C2(t), demonstrate a clear power law relationship of C2(t) C1(t)m with m = 3, 2, and 1 in the initial, free Rouse, and entangled (constrained Rouse) regimes, respectively. The return-to-origin hypothesis, which leads to inverse proportionality between the segmental orientation autocorrelation functions and g1(t) in the entangled regime, is convincingly verified by the simulation result of C1(t) g1(t)−1 t–1/4 in the constrained Rouse regime, where for well-entangled chains both C1(t) and g1(t) are rather insensitive to the constraint release effects. However, the second-order correlation function, C2(t), shows much stronger sensitivity to the constraint release effects and experiences a protracted crossover from the free Rouse to entangled regime. This crossover region extends for at least one decade in time longer than that of C1(t). The predicted time scaling behavior of C2(t) t–1/4 is observed in slip-springs simulations only at chain length of 90 entanglements, whereas shorter chains show higher scaling exponents. The reported simulation work can be applied to understand the observations of the NMR experiments.

Finite-stretching corrections to the Milner-Witten-Cates theory for polymer brushes

This paper investigates finite-stretching corrections to the classical Milner-Witten-Cates theory for semi-dilute polymer brushes in a good solvent. The dominant correction to the free energy originates from an entropic repulsion caused by the impenetrability of the grafting surface, which produces a depletion of segments extending a distance $\mu \propto L^{-1}$ from the substrate, where $L$ is the classical brush height. The next most important correction is associated with the translational entropy of the chain ends, which creates the well-known tail where a small population of chains extend beyond the classical brush height by a distance $\xi \propto L^{-1/3}$. The validity of these corrections is confirmed by quantitative comparison with numerical self-consistent field theory.

Significance of cross correlations in the stress relaxation of polymer melts

According to linear response theory, all relaxation functions in the linear regime can be obtained using time correlation functions calculated under equilibrium. In this paper, we demonstrate that the cross correlations make a significant contribution to the partial stress relaxation functions in polymer melts. We present two illustrations in the context of polymer rheology using (1) Brownian dynamics simulations of a single chain model for entangled polymers, the slip-spring model, and (2) molecular dynamics simulations of a multichain model. Using the single chain model, we analyze the contribution of the confining potential to the stress relaxation and the plateau modulus. Although the idea is illustrated with a particular model, it applies to any single chain model that uses a potential to confine the motion of the chains. This leads us to question some of the assumptions behind the tube theory, especially the meaning of the entanglement molecular weight obtained from the plateau modulus. To shed some light on this issue, we study the contribution of the nonbonded excluded-volume interactions to the stress relaxation using the multichain model. The proportionality of the bonded/nonbonded contributions to the total stress relaxation (after a density dependent "colloidal" relaxation time) provides some insight into the success of the tube theory in spite of using questionable assumptions. The proportionality indicates that the shape of the relaxation spectrum can indeed be reproduced using the tube theory and the problem is reduced to that of finding the correct prefactor. (c) 2007 American Institute of Physics

Peptide mediated formation of hierarchically organized solution and solid state polymer nanostructures

Biologically-inspired peptide sequences have been explored as auxiliaries to mediate self-assembly of synthetic macromolecules into hierarchically organized solution and solid state nanostructures. Peptide sequences inspired by the coiled coil motif and "switch" peptides, which can adopt both amphiphilic alpha-helical and beta-strand conformations, were conjugated to poly(ethylene glycol) (PEG). The solution and solid state self-assembly of these materials was investigated using a variety of spectroscopic, scattering and microscopic techniques. These experiments revealed that the folding and organization properties of the peptide sequences are retained upon conjugation of PEG and that they provide the driving force for the formation of the different nanoscale structures which were observed. The possibility of using defined peptide sequences to direct structure formation of synthetic polymers together with the potential of peptide sequences to induce a specific biological response offers interesting prospects for the development of novel self-assembled and biologically active materials.

Compression of polymer brushes: quantitative comparison of self-consistent field theory with experiment

The self-consistent field theory (SCFT) prediction for the compression force between two semi-dilute polymer brushes is compared to the benchmark experiments of Taunton et al. [Nature, 1988, 332, 712]. The comparison is done with previously established parameters, and without any fitting parameters whatsoever. The SCFT provides a significant quantitative improvement over the classical strong-stretching theory (SST), yielding excellent quantitative agreement with the experiment. Contrary to earlier suggestions, chain fluctuations cannot be ignored for normal experimental conditions. Although the analytical expressions of SST provide invaluable aids to understanding the qualitative behavior of polymeric brushes, the numerical SCFT is necessary in order to provide quantitatively accurate predictions.

Comparing tube models for predicting the linear rheology of branched polymer melts

The hierarchical and "bob" (or branch-on-branch) models are tube-based computational models recently developed for predicting the linear rheology of general mixtures of polydisperse branched polymers. These two models are based on a similar tube-theory framework but differ in their numerical implementation and details of relaxation mechanisms. We present a detailed overview of the similarities and differences of these models and examine the effects of these differences on the predictions of the linear viscoelastic properties of a set of representative branched polymer samples in order to give a general picture of the performance of these models. Our analysis confirms that the hierarchical and bob models quantitatively predict the linear rheology of a wide range of branched polymer melts but also indicate that there is still no unique solution to cover all types of branched polymers without case-by-case adjustment of parameters such as the dilution exponent alpha and the factor p(2) which defines the hopping distance of a branch point relative to the tube diameter. An updated version of the hierarchical model, which shows improved computational efficiency and refined relaxation mechanisms, is introduced and used in these analyses.

Linear viscoelasticity from molecular dynamics simulation of entangled polymers

The linear viscoelastic (LVE) spectrum is one of the primary fingerprints of polymer solutions and melts, carrying information about most relaxation processes in the system. Many single chain theories and models start with predicting the LVE spectrum to validate their assumptions. However, until now, no reliable linear stress relaxation data were available from simulations of multichain systems. In this work, we propose a new efficient way to calculate a wide variety of correlation functions and mean-square displacements during simulations without significant additional CPU cost. Using this method, we calculate stress−stress autocorrelation functions for a simple bead−spring model of polymer melt for a wide range of chain lengths, densities, temperatures, and chain stiffnesses. The obtained stress−stress autocorrelation functions were compared with the single chain slip−spring model in order to obtain entanglement related parameters, such as the plateau modulus or the molecular weight between entanglements. Then, the dependence of the plateau modulus on the packing length is discussed. We have also identified three different contributions to the stress relaxation:  bond length relaxation, colloidal and polymeric. Their dependence on the density and the temperature is demonstrated for short unentangled systems without inertia.