Mass dependence of the activation enthalpy and entropy of unentangled linear alkane chains

Understanding Relaxation in the Kob-Andersen Liquid Based on Entropy, String, Shoving, Localization, and Parabolic Models

We assess the validity of a range of models of glass formation based on molecular dynamics simulation results of the Kob-Andersen (KA) model system under a wide range of constant volume and constant pressure conditions. These models include the Adam-Gibbs model emphasizing configurational entropy, the string model emphasizing collective particle exchange motion, the shoving model emphasizing material elasticity, the localization model emphasizing dynamical free volume, and parabolic models based on the ideas of dynamic facilitation and, alternatively, the hypothesis that glass formation involves an avoided critical point. We demonstrate that these seemingly disparate models all provide a reasonable description of structural relaxation and diffusion data for the KA model system under all simulation conditions considered. Hence, the present study points to some unity in our understanding of the relationship between leading models of glass formation, supporting inferences drawn from previous studies of polymeric glass-forming liquids.

Diffusion Coefficients of Variable-Size Amphiphilic Additives in a Glass-Forming Polyethylene Matrix

Diffusion of additives in polymers is an important issue in the plastics industry since migratory-type molecules are widely used to tune the properties of polymeric composites. Predicting the diffusional behavior of new additives can minimize the need for repetitive experiments. This work presents molecular dynamics simulations at the microsecond time scale and uses the MARTINI force field to estimate self-diffusion coefficients, D, of six monounsaturated amides and their analogs carboxylic acids in polyethylene matrices (PE, MW = 5600 Da). The results are strongly influenced by the glass-forming properties of the PE matrix, which we characterize by three distinct temperatures. The metastability region (T < 325 K), the glass transition temperature (Tg = 256-260 K), and the end of the transition (T ≅ 200 K). Self-diffusion mechanisms are inferred from the results of the dependence of D on the molecular mass of the additive, observing a Rouse-like behavior at high temperatures and deviations from it within the metastability region of the matrix. Interestingly, D values are nonsensitive to the nature of the considered polar head for additives of similar size. The temperature-dependent behavior of D follows, at fixed additive size, a linear Arrhenius pattern at high temperatures and a super Arrhenius trend at lower temperatures, which is well represented with a power law equation as predicted by the Mode Coupling Theory (MCT). We offer a conceptual explanation for the observed super-Arrhenius behavior. This explanation draws on Truhlar and Kohen's interpretation of the available energies at both the initial and the transition states along the diffusion pathway. The matrix's mobility significantly affects solute self-diffusion, yielding equal activation enthalpies for the Arrhenius region or the same power law parameters for the super-Arrhenius regime. Finally, we establish a one-to-one time-equivalence of the self-diffusion processes between CG and all-atom systems for the largest additives and the PE matrix in the high-temperature regime.

Coexistence of Crumpling and Flat Sheet Conformations in Two-Dimensional Polymer Networks: An Understanding of Aggrecan Self-Assembly

We investigate the conformational properties of self-avoiding two-dimensional (2D) ideal polymer networks with tunable mesh sizes as a model of self-assembled structures formed by aggrecan. Polymer networks having few branching points and large enough mesh tend to crumple, resulting in a fractal dimension of d_{f}≈2.7. The flat sheet behavior (d_{f}=2) emerges in 2D polymer networks having more branching points at large length scales; however, it coexists with crumpling conformations at intermediate length scales, a feature found in scattering profiles of aggrecan solutions. Our findings bridge the long-standing gap between theories and simulations of polymer sheets.

A β-NMR study of the depth, temperature, and molecular-weight dependence of secondary dynamics in polystyrene: Entropy–enthalpy compensation and dynamic gradients near the free surface

We investigated the depth, temperature, and molecular-weight (MW) dependence of the γ-relaxation in polystyrene glasses using implanted 8Li+ and β-detected nuclear magnetic resonance. Measurements were performed on thin films with MW ranging from 1.1 to 641 kg/mol. The temperature dependence of the average 8Li spin–lattice relaxation time [Formula: see text] was measured near the free surface and in the bulk. Spin–lattice relaxation is caused by phenyl ring flips, which involve transitions between local minima over free-energy barriers with enthalpic and entropic contributions. We used transition state theory to model the temperature dependence of the γ-relaxation, and hence [Formula: see text]. There is no clear correlation of the average entropy of activation [Formula: see text] and enthalpy of activation [Formula: see text] with MW, but there is a clear correlation between [Formula: see text] and [Formula: see text], i.e., entropy–enthalpy compensation. This results in the average Gibbs energy of activation, [Formula: see text], being approximately independent of MW. Measurements of the temperature dependence of [Formula: see text] as a function of depth below the free surface indicate the inherent entropic barrier, i.e., the entropy of activation corresponding to [Formula: see text] = 0, has an exponential dependence on the distance from the free surface before reaching the bulk value. This results in [Formula: see text] near the free surface being lower than the bulk. Combining these observations results in a model where the average fluctuation rate of the γ-relaxation has a “double-exponential” depth dependence. This model can explain the depth dependence of [Formula: see text] in polystyrene films. The characteristic length of enhanced dynamics is ∼6 nm and approximately independent of MW near room temperature.

Simulation of the Coronal Dynamics of Polymer-Grafted Nanoparticles

Polymer-grafted nanoparticles (PGNPs) are an important component of many advanced materials. The interplay between the nanoparticle surface curvature and spatial confinement by neighboring chains produces a complex set of structural and dynamical behaviors in the polymer corona surrounding the nanoparticle. For example, experiments have shown that the inner portion of the corona is more stretched and relaxes more slowly than the outer region. Here, we perform systematic core-modified dissipative particle dynamics (CM-DPD) simulations and analyze the relaxation dynamics using proper orthogonal decomposition (POD) of the monomer coordinates. We find that grafted chains relax more slowly than free chains and that the relaxation time of the grafted chains scales inversely with the confinement strength. For PGNPs in a polymer melt, the relaxation processes are always Rouse-like. However, we observe either Zimm-like or Rouse-like dynamics for PGNPs in solution depending on the confinement strength.

Dynamically consistent coarse-grain simulation model of chemically specific polymer melts via friction parameterization

Coarse-grained (CG) models of polymers involve grouping many atoms in an all-atom (AA) representation into single sites to reduce computational effort yet retain the hierarchy of length and time scales inherent to macromolecules. Parameterization of such models is often via "bottom-up" methods, which preserve chemical specificity but suffer from artificially accelerated dynamics with respect to the AA model from which they were derived. Here, we study the combination of a bottom-up CG model with a dissipative potential as a means to obtain a chemically specific and dynamically correct model. We generate the conservative part of the force-field using the iterative Boltzmann inversion (IBI) method, which seeks to recover the AA structure. This is augmented with the dissipative Langevin thermostat, which introduces a single parameterizable friction factor to correct the unphysically fast dynamics of the IBI-generated force-field. We study this approach for linear polystyrene oligomer melts for three separate systems with 11, 21, and 41 monomers per chain and a mapping of one monomer per CG site. To parameterize the friction factor, target values are extracted from the AA dynamics using translational monomer diffusion, translational chain diffusion, and rotational chain motion to test the consistency of the parameterization across different modes of motion. We find that the value of the friction parameter needed to bring the CG dynamics in line with AA target values varies based on the mode of parameterization with short-time monomer translational dynamics requiring the highest values, long-time chain translational dynamics requiring the lowest values, and rotational dynamics falling in between. The friction ranges most widely for the shortest chains, and the span narrows with increasing chain length. For longer chains, a practical working value of the friction parameter may be derived from the rotational dynamics, owing to the contribution of multiple relaxation modes to chain rotation and a lack of sensitivity of the translational dynamics at these intermediate levels of friction. A study of equilibrium chain structure reveals that all chains studied are non-Gaussian. However, longer chains better approximate ideal chain dimensions than more rod-like shorter chains and thus are most closely described by a single friction parameter. We also find that the separability of the conservative and dissipative potentials is preserved.

Observation of General Entropy-Enthalpy Compensation Effect in the Relaxation of Wrinkled Polymer Nanocomposite Films

Surface-textured polymer nanocomposite (PNC) films are utilized in many device applications, and therefore understanding the relaxation behavior of such films is important. By extending an in situ wrinkle relaxation method, we observed that the thermal stability of wrinkled PNC films, both above and below the glass transition temperature (T_g), is proportional to a film's nanoparticle (polymer grafted and bare) concentration, with a slope that changes sign at a compensation temperature (T_comp) that is determined to be in the vicinity of the film's T_g. This provides unambiguous confirmation of entropy-enthalpy compensation (EEC) as a general feature of PNC films, implying that the stability of PNC films changes from being enhanced to becoming diminished by simply passing through this characteristic temperature, a phenomenon having evident practical ramifications. We suggest EEC will also arise in films where residual stresses are associated with the film fabrication process, which is relevant to nanotech device applications.

Conformation and dynamics of ring polymers under symmetric thin film confinement

Understanding the structure and dynamics of polymers under confinement has been of widespread interest, and one class of polymers that have received comparatively little attention under confinement is that of ring polymers. The properties of non-concatenated ring polymers can also be important in biological fields because ring polymers have been proven to be a good model to study DNA organization in the cell nucleus. From our previous study, linear polymers in a cylindrically confined polymer melt were found to segregate from each other as a result of the strong correlation hole effect that is enhanced by the confining surfaces. By comparison, our subsequent study of linear polymers in confined thin films at similar levels of confinements found only the onset of segregation. In this study, we use molecular dynamics simulation to investigate the chain conformations and dynamics of ring polymers under planar (1D) confinement as a function of film thickness. Our results show that conformations of ring polymers are similar to the linear polymers under planar confinement, except that ring polymers are less compressed in the direction normal to the walls. While we find that the correlation hole effect is enhanced under confinement, it is not as pronounced as the linear polymers under 2D confinement. Finally, we show that chain dynamics far above Tg are primarily affected by the friction from walls based on the monomeric friction coefficient we get from the Rouse mode analysis.

Predictive relation for the α-relaxation time of a coarse-grained polymer melt under steady shear

We examine the influence of steady shear on structural relaxation in a simulated coarse-grained unentangled polymer melt over a wide range of temperature and shear rates. Shear is found to progressively suppress the α-relaxation process observed in the intermediate scattering function, leading ultimately to a purely inertially dominated β-relaxation at high shear rates, a trend similar to increasing temperature. On the basis of a scaling argument emphasizing dynamic heterogeneity in cooled liquids and its alteration under material deformation, we deduce and validate a parameter-free scaling relation for both the structural relaxation time τα from the intermediate scattering function and the "stretching exponent" β quantifying the extent of dynamic heterogeneity over the entire range of temperatures and shear rates that we can simulate.

A blob model to parameterize polymer hole free volumes and solute diffusion

Solute diffusion in solid polymers has tremendous applications in packaging, reservoir, and biomedical technologies but remains poorly understood. Diffusion of non-entangled linear solutes with chemically identical patterns (blobs) deviates dramatically in polymers in the solid-state (α_lin > 1, Macromolecules 2013, 46, 874) from their behaviors in the molten state (α_lin = 1, Macromolecules, 2007, 40, 3970). This work uses the scale invariance of the diffusivities, D, of linear probes D(N·M_blob + M_anchor,T,T_g) = N^{-α_lin(T,T_g)}D(M_blob + M_anchor,T,T_g) comprising N identical blobs of mass M_blob and possibly one different terminal pattern (anchor of mass M_anchor) to evaluate the amounts of hole-free volume in seven polymers (aliphatic, semi-aromatic and aromatic) over a broad range of temperatures (-70 K ≤T-T_g≤ 160 K). The new parameterization of the concept of hole-free volumes opens the application of the free-volume theory (FVT) developed by Vrentas and Duda to practically any polymer, regardless of the availability of free-volume parameters. The quality of the estimations was tested with various probes including n-alkanes, 1-alcohols, n-alkyl acetates, and n-alkylbenzene. The effects of enthalpic and entropic effects of the blobs and the anchor were analyzed and quantified. Blind validation of the reformulated FVT was tested successfully by predicting from first principles the diffusivities of water and toluene in amorphous polyethylene terephthalate from 4 °C to 180 °C and in various other polymers. The new blob model would open the rational design of additives with controlled diffusivities in thermoplastics.

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.

Influence of knot complexity on glass-formation in low molecular mass ring polymer melts

We perform molecular dynamics simulations on a coarse-grained polymer melt to study the dynamics of glass-formation in ring polymer melts of variable knot complexity. After generating melts of non-concatenated polymeric rings having a range of minimum crossing number values, mc, we compute the coherent intermediate scattering function, the segmental α-relaxation time, fragility, and the glass transition temperature as a function of mc. Variation of knot complexity is found to have a pronounced effect on the dynamics of polymer melts since both molecular rigidity and packing are altered, primary physical factors governing glass-formation in polymeric materials.

Weak and Strong Gels and the Emergence of the Amorphous Solid State

Gels are amorphous solids whose macroscopic viscoelastic response derives from constraints in the material that serve to localize the constituent molecules or particles about their average positions in space. These constraints may either be local in nature, as in chemical cross-linking and direct physical associations, or non-local, as in case of topological "entanglement" interactions between highly extended fiber or sheet structures in the fluid. Either of these interactions, or both combined, can lead to "gelation" or "amorphous solidification". While gels are often considered to be inherently non-equilibrium materials, and correspondingly termed "soft glassy matter", this is not generally the case. For example, the formation of vulcanized rubbers by cross-linking macromolecules can be exactly described as a second order phase transition from an equilibrium fluid to an equilibrium solid state, and amorphous solidification also arises in diverse physical gels in which molecular and particle localization occurs predominantly through transient molecuar associations, or even topological interactions. As equilibrium, or near equilibrium systems, such gels can be expected to exhibit universal linear and non-linear viscoelastic properties, especially near the "critical" conditions at which the gel state first emerges. In particular, a power-law viscoelastic response is frequently observed in gel materials near their "gelation" or "amorphous solidification" transition. Another basic property of physical gels of both theoretical and practical interest is their response to large stresses at constant shear rate or under a fixed macrocopic strain. In particular, these materials are often quite sensitive to applied stresses that can cause the self-assembled structure to progressively break down under flow or deformation. This disintegration of gel structure can lead to "yield" of the gel material, i.e., a fluidization transition, followed by shear thinning of the resulting heterogeneous "jelly-like" fluid. When the stress is removed, however, the material can relax back to its former equilibrium gel state, i.e., gel rejuvenation. In constrast, a non-equilibrium material will simply change its form and properties in a way that depends on processing history. Physical gels are thus unique self-healing materials in which the existence of equilibrium ensures their eventual recovery. The existence of equilibrium also has implications for the nature of both the linear and non-linear rheological response of gel materials, and the present paper explores this phenomenon based on simple scaling arguments of the kind frequently used in describing phase transitions and the properties of polymer solutions.

Energy Renormalization Method for the Coarse-Graining of Polymer Viscoelasticity

Developing temperature transferable coarse-grained (CG) models is essential for the computational prediction of polymeric glass-forming (GF) material behavior, but their dynamics are often greatly altered from those of all-atom (AA) models mainly because of the reduced fluid configurational entropy under coarse-graining. To address this issue, we have recently introduced an energy renormalization (ER) strategy that corrects the activation free energy of the CG polymer model by renormalizing the cohesive interaction strength ε as a function of temperature T, i.e., ε(T), thus semiempirically preserving the T-dependent dynamics of the AA model. Here we apply our ER method to consider-in addition to T-dependency-the frequency f-dependent polymer viscoelasticity. Through smallamplitude oscillatory shear molecular dynamics simulations, we show that changing the imposed oscillation f on the CG systems requires changes in ε values (i.e., ε(T, f)) to reproduce the AA viscoelasticity. By accounting for the dynamic fragility of polymers as a material parameter, we are able to predict ε(T, f) under coarse-graining in order to capture the AA viscoelasticity, and consequently the activation energy, across a wide range of T and f in the GF regime. Specifically, we showcase our achievements on two representative polymers of distinct fragilities, polybutadiene (PB) and polystyrene (PS), and show that our CG models are able to sample viscoelasticity up to the megahertz regime, which approaches state-of-the-art experimental resolutions, and capture results sampled via AA simulations and prior experiments.

Influence of polymer architectures on diffusion in unentangled polymer melts

Recent simulations have indicated that the thermodynamic properties and the glassy dynamics of polymer melts are strongly influenced by the average molecular shape, as quantified by the radius of gyration tensor of the polymer molecules, and that the average molecular shape can be tuned by varying the molecular topology (e.g., ring, star, linear chain, etc.). In the present work, we investigate if the molecular shape is similarly a predominant factor in understanding the polymer center of mass diffusion D in the melt, as already established for polymer solutions. We find that all our D data for unentangled polymer melts having a range of topologies can be reasonably described as a power law of the polymer hydrodynamic radius, Rh. In particular, this scaling is similar to the scaling of D for a tracer sphere having a radius on the order of the chain radius of gyration, Rg. We conclude that the chain topology influences the molecular dynamics in as much as the polymer topology influences the average molecular shape. Experimental evidence seems to suggest that this situation is also true for entangled polymer melts.

Buckling a Semiflexible Polymer Chain under Compression

Instability and structural transitions arise in many important problems involving dynamics at molecular length scales. Buckling of an elastic rod under a compressive load offers a useful general picture of such a transition. However, the existing theoretical description of buckling is applicable in the load response of macroscopic structures, only when fluctuations can be neglected, whereas membranes, polymer brushes, filaments, and macromolecular chains undergo considerable Brownian fluctuations. We analyze here the buckling of a fluctuating semiflexible polymer experiencing a compressive load. Previous works rely on approximations to the polymer statistics, resulting in a range of predictions for the buckling transition that disagree on whether fluctuations elevate or depress the critical buckling force. In contrast, our theory exploits exact results for the statistical behavior of the worm-like chain model yielding unambiguous predictions about the buckling conditions and nature of the buckling transition. We find that a fluctuating polymer under compressive load requires a larger force to buckle than an elastic rod in the absence of fluctuations. The nature of the buckling transition exhibits a marked change from being distinctly second order in the absence of fluctuations to being a more gradual, compliant transition in the presence of fluctuations. We analyze the thermodynamic contributions throughout the buckling transition to demonstrate that the chain entropy favors the extended state over the buckled state, providing a thermodynamic justification of the elevated buckling force.

Energy-Renormalization for Achieving Temperature Transferable Coarse-Graining of Polymer Dynamics

The bottom-up prediction of the properties of polymeric materials based on molecular dynamics simulation is a major challenge in soft matter physics. Coarse-grained (CG) models are often employed to access greater spatiotemporal scales required for many applications, but these models normally experience significantly altered thermodynamics and highly accelerated dynamics due to the reduced number of degrees of freedom upon coarse-graining. While CG models can be calibrated to meet certain properties at particular state points, there is unfortunately no temperature transferable and chemically specific coarse-graining method that allows for modeling of polymer dynamics over a wide temperature range. Here, we pragmatically address this problem by "correcting" for deviations in activation free energies that occur upon coarse-graining the dynamics of a model polymeric material (polystyrene). In particular, we propose a new strategy based on concepts drawn from the Adam-Gibbs (AG) theory of glass formation. Namely we renormalize the cohesive interaction strength and effective interaction length-scale parameters to modify the activation free energy. We show that this energy-renormalization method for CG modeling allows accurate prediction of atomistic dynamics over the Arrhenius regime, the non-Arrhenius regime of glass formation, and even the non-equilibrium glassy regime, thus allowing for the predictive modeling of dynamic properties of polymer over the entire range of glass formation. Our work provides a practical scheme for establishing temperature transferable coarse-grained models for predicting and designing the properties of polymeric materials.

Impact of Monovalent Counter-ions on the Conformation of Flexible Polyelectrolytes Having Different Molecular Architectures

We explore the impact of monovalent counter-ions on the molecular conformation of highly charged flexible polyelectrolytes for a range of molecular topologies (linear chains, stars, and unknotted and trefoil rings) by molecular dynamics simulations that include an explicit solvent having short range interaction with the polyelectrolyte. In particular, we investigate how the counter-ions near the polyelectrolytes with variable mass influence the average molecular shape. We also characterize the interfatially "bound" counter-ions by calculating the time-averaged number of interfacial counter-ions, as well as the degree to which the polyelectrolytes wrap around the counter-ions by calculating the number of contacts between the counter-ions and the polyelectrolyte.