Branched polymers stress relaxation is at the center to their function as viscosity modifiers, though the fundamentals that underlie the correlation between the polymer topology and their impact on viscosity remains an open question. Here, the stress relaxation of short, branched polyethylene comb polymer melts is studied by molecular dynamics simulations. A coarse-grained model where four methylene groups constitute one bead is used, and the results are transposed to the atomistic level. For arms of length comparable to entanglement length ne of the linear polymer, we show that while increasing the number of branches with the same arm length decreases the plateau modulus, the terminal diffusive time does not change significantly. Increasing the arm length decreases the plateau modulus and increases the terminal time. As arms shorter than ne relax by the entanglement time, both the chain mobility and stress relaxation can be described by reptation of the backbone with an increased tube diameter and an increased friction coefficient; or in other words, the branches act as a solvent.
more »
« less
Viscoelastic Response of Dispersed Entangled Polymer Melts
Polymer synthesis routes result in macromolecules with molecular weight dispersity ĐM that depends on the polymerization mechanism. The lowest dispersity polymers are those made by anionic and atom-transfer radical polymerization, which exhibit narrow distributions ĐM = Mw/Mn ∼ 1.02–1.04. Even for small dispersity, the chain length can vary by a factor of two from the average. The impact of chain length dispersity on the viscoelastic response remains an open question. Here, the effects of dispersity on stress relaxation and shear viscosity of entangled polyethylene melts are studied using molecular dynamics simulations. Melts with chain length dispersity, which follow a Schulz–Zimm (SZ) distribution with ĐM = 1.0–1.16, are studied for times up to 800 μs, longer than the terminal time. These systems are compared to those with binary and ternary distributions. The stress relaxation functions are extracted from the Green–Kubo relation and from stress relaxation following a uniaxial extension. At short and intermediate time scales, both the mean squared displacement and the stress relaxation function G(t) are independent of ĐM. At longer times, the terminal relaxation time decreases with increasing ĐM. In this time range, the faster motion of the shorter chains results in constraint release for the longer chains.
more »
« less
- Award ID(s):
- 1905407
- NSF-PAR ID:
- 10350227
- Date Published:
- Journal Name:
- Macromolecules
- ISSN:
- 1539-8870
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
more » « less
Dissipative particle dynamics (DPD) simulations are performed on coarse-grained replicas of linear, monodisperse entangled polyethylene melts [Formula: see text] and [Formula: see text] undergoing both steady-state and transient planar elongational flow (PEF). The fidelity of the DPD simulations is verified by direct comparison of flow topological and rheological properties of a 334-particle chain liquid against the united-atom [Formula: see text] liquid, simulated using nonequilibrium molecular dynamics (NEMD). These DPD simulations demonstrate that a flow-induced coil-stretch transition (CST) and its associated hysteresis caused by configurational microphase separation, as observed in previous NEMD simulations of PEF, can be replicated using a more computationally efficient coarse-grained system. Results indicate that the breadth of the CST hysteresis loop is enlarged for the longer molecule liquid relative to the shorter one. Furthermore, relaxation simulations reveal that reducing the applied flow Deborah number ([Formula: see text]) from a high value corresponding to a homogeneous phase of highly stretched molecules to a [Formula: see text] within the biphasic region results in a two-stage relaxation process. There is a fast initial stratification of the kinetically trapped highly stretched chains into regions of highly extended and less extended chains, which displays similar behavior to a system undergoing a spinodal decomposition caused by spatial configurational free energy fluctuations. After a short induction period of apparently random duration, the less extended chain regions experience a stochastic nucleation event that induces configurational relaxation to domains composed of coiled molecules over a much longer time scale, leaving the more highly extended chains in surrounding sheetlike domains. The time scales of these two relaxation processes are of the same order of magnitude as the Rouse and disengagement times of the equilibrium liquids.
-
A series of poly( N -acryloyl glycinamide) (pNAGA) polymers were synthesized and studied as capture agents for surface-enhanced Raman scattering (SERS) detection of aflatoxin B1 (AFB1), a highly carcinogenic food-borne toxin. Four molecular weights of pNAGA were synthesized by reversible addition-fragmentation chain-transfer (RAFT) polymerization to study the dependence of affinity agent efficacy on chain length for this AFB1 sensing platform. Isothermal titration calorimetry (ITC) was used to verify the sign and magnitude of the enthalpic effects involved in polymer–AFB1 interactions in solution and to understand the effects of pNAGA chain length on AFB1 noncovalent binding. pNAGA–AFB1 interactions were found to be exothermic, and longer pNAGA chains generally resulted in smaller enthalpy decreases per repeat unit. With pNAGA 22 being thermodynamically the strongest affinity agent, we hypothesize that AFB1 affinity is determined by a balance between the configurational restrictions in pNAGA chains and the enthalpic advantage of binding AFB1. SERS spectral changes observed following AFB1 exposure were used to evaluate the influence of polymer molecular weight (2.0–5.2 kDa), order of attachment (pre- vs. post- functionalization of the substrate) and attachment chemistry (thiol vs. trithiocarbonate) on the sensitivity of AFB1 detection. The method by which target, polymer affinity agent, and signal transduction mechanism are combined was found to have significant impacts on the achieved sensitivity. The most effective polymer chain length (pNAGA 22 ), anchoring chemistry (thiol), and polymer/toxin assembly scheme (in-solution) allowed detection of 10 ppb AFB1 in water (below the FDA regulatory limit of 20 ppb), a hundred-fold improvement over SERS sensing without the pNAGA affinity agent.more » « less
-
Polymer molecular weight, or chain length distributions, are a core characteristic of a polymer system, with the distribution being intimately tied to the properties and performance of the polymer material. A model is developed for the ideal distribution of polymers made using reversible activation/deactivation of chain ends, with monomer added to the active form of the chain end. The ideal distribution focuses on living chains, with the system having minimal impact from irreversible termination or transfer. This model was applied to ATRP, RAFT, and cationic polymerizations, and was also used to describe complex systems such as blended polymers and block copolymers. The model can easily and accurately be fitted to molecular weight distributions, giving information on the ratio of propagation to deactivation, as well as the mean number of times a chain is activated/deactivated under the polymerization conditions. The mean number of activation cycles per chain is otherwise difficult to assess from conversion data or molecular weight distributions. Since this model can be applied to wide range of polymerizations, giving useful information on the underlying polymerization process, it can be used to give fundamental insights into macromolecular synthesis and reaction outcomes.more » « less
-
more » « less
Abstract Dielectric spectroscopy is extremely powerful to study molecular dynamics, because of the very broad frequency range. Often multiple processes superimpose resulting in spectra that expand over several orders of magnitude, with some of the contributions partially hidden. For illustration, we selected two examples, (i) normal mode of high molar mass polymers partially hidden by conductivity and polarization and (ii) contour length fluctuations partially hidden by reptation using the well-studied polyisoprene melts as example. The intuitive approach to describe experimental spectra and to extract relaxation times is the addition of two or more model functions. Here, we use the empirical Havriliak-Negami function to illustrate the ambiguity of the extracted relaxation time, despite an excellent agreement of the fit with experimental data. We show that there are an infinite number of solutions for which a perfect description of experimental data can be achieved. However, a simple mathematical relationship indicates uniqueness of the pairs of the relaxation strength and relaxation time. Sacrificing the absolute value of the relaxation time enables to find the temperature dependence of the parameters with a high accuracy. For the specific cases studied here, the time temperature superposition (TTS) is very useful to confirm the principle. However, the derivation is not based on a specific temperature dependence, hence, independent from the TTS. We compare new and traditional approaches and find the same trend for the temperature dependence. The important advantage of the new technology is the knowledge of the accuracy of the relaxation times. Relaxation times determined from data for which the peak is clearly visible are the same within the experimental accuracy for traditional and new technology. However, for data where a dominant process hides the peak, substantial deviations can be observed. We conclude that the new approach is particularly helpful for cases in which relaxation times need to be determined without having access to the associated peak position.
- Free Publicly Accessible Full Text
-
Accepted Manuscript1.0
- Journal Article:
- The DOI is not currently available.
