Heterogeneous Rouse Model Predicts Polymer Chain Translational Normal Mode Decoupling

Generalized entropy theory investigation of the relatively high segmental fragility of many glass-forming polymers

We utilize the generalized entropy theory (GET) of glass formation to address one of the most singular and least understood properties of polymer glass-forming liquids in comparison to atomic and small molecule liquids-the often relatively high fragility of the polymer dynamics on a segmental scale, ms. Based on this highly predictive framework of both the thermodynamics and segmental dynamics in terms of molecular structure, polymer backbone and side-group rigidities, and intermolecular interaction strength, we first analyze the relation between ms and the ratio, , where Sc is the configurational entropy density of the polymer fluid, equals Sc at the onset temperature TA for non-Arrhenius relaxation, and Tg is the glass transition temperature at which the structural relaxation time τα equals 100 s. While the reduced activation energy estimated from an Arrhenius plot (i.e., differential activation energy) normalized by kBTg is determined to be not equal to the actual activation energy, we do find that an apparently general nonlinear relation between ms and holds to a good approximation for a large class of polymer models, . The predicted ranges of ms and are consistent with experimental estimates for high molecular-mass polymer, oligomeric, small molecule, and atomic glass-forming liquids. In particular, relatively high values of ms are found for polymers having complex monomer structures and significant chain stiffness. The variation of ms with molecular mass, chain stiffness, and intermolecular interaction strength can be traced to the variation of , which is shown to provide a measure of packing frustration defined in terms of the dimensionless thermal expansion coefficient and isothermal compressibility. The often relatively high fragility and large extent of cooperative motion are found in the GET to derive from the often relatively large packing frustration in this class of polymer glass-forming liquids. Finally, we also develop a tentative model of the "dynamical segmental relaxation time" based on the GET, in which the polymers on a coarse-grained scale are modeled as strings of structureless "beads", as assumed in the Rouse and reptation models of polymer dynamics.

Rouse mode analysis of chain relaxation in reversibly crosslinked polymer melts

Polymer melts with chains undergoing reversible crosslinking have distinctively favorable dynamic properties, e.g., self-healing and reprocessability. In these situations, there are two relevant elementary time scales: the segmental and the sticker association times. A convenient framework to model these situations is the sticky Rouse model and here we perform hybrid molecular dynamics (MD)-Monte Carlo (MC) simulations to examine its relevance. In agreement with the underpinning idea discussed above, we find that reversibly crosslinked chains show two distinct modes of relaxation behavior depending on the magnitude of bond lifetimes. For bond lifetimes shorter than the chain end-to-end relaxation time, the polymers exhibit essentially Rouse-like dynamics, but with an apparently increased local friction relative to the non-sticky analog. For longer bond lifetimes, the chains exhibit two modes of relaxation: the faster mode is independent of the bond lifetime, but the slower mode is controlled by it. However, these slower mode results are not consistent with the predictions of the sticky Rouse model. Our Rouse mode analysis as a function of chain length, N, implies that this is likely a result of the relatively short N's employed, but it nevertheless suggests that theories need to include these short chain effects if they are to be relevant to experimental systems with short chains following Rouse dynamics.

Does the Naı̈ve Mode-Coupling Power Law Divergence Provide an Objective Determination of the Crossover Temperature in Glass Formation Behavior?

The glass formation temperature range is commonly divided into a weakly supercooled regime at higher temperatures and a deeply supercooled regime at lower temperatures, with a change in the physical mechanisms that govern dynamics often postulated to occur at the crossover between these regimes. This crossover temperature Tc is widely determined based on a fit of relaxation time vs temperature data to a power law divergence form predicted by the naı̈ve mode coupling theory (MCT). Here, we show, based on simulation data spanning polymeric, small molecule organic, metallic, and inorganic glass formers, that this approach does not yield an objective measure of a crossover temperature. Instead, the value of Tc is determined by the lowest temperature Tmin employed in the fit, and no regime of stationary or convergent Tc value is generally observed as Tmin is varied. Nor does the coefficient of determination R2 provide any robust means of selecting a fit range and thus a value of Tc. These results may require a re-evaluation of published results that have employed the fit MCT Tc value as a metric of temperature-dependent dynamics or a benchmark for depth of supercooling, and they highlight a need for the field to converge on a more objective determination of any posited crossover temperature between high and low temperature regimes of glass formation.

Interplay between dynamic heterogeneity and interfacial gradients in a model polymer film

Glass-forming liquids exhibit long-lived, spatially correlated dynamical heterogeneity, in which some nm-scale regions in the fluid relax more slowly than others. In the nanoscale vicinity of an interface, glass-formers also exhibit the emergence of massive interfacial gradients in glass transition temperature Tg and relaxation time τ. Both of these forms of heterogeneity have a major impact on material properties. Nevertheless, their interplay has remained poorly understood. Here, we employ molecular dynamics simulations of polymer thin films in the isoconfigurational ensemble in order to probe how bulk dynamic heterogeneity alters and is altered by the large gradient in dynamics at the surface of a glass-forming liquid. Results indicate that the τ spectrum at the surface is broader than in the bulk despite being shifted to shorter times, and yet it is less spatially correlated. This is distinct from the bulk, where the τ distribution becomes broader and more spatially organized as the mean τ increases. We also find that surface gradients in slow dynamics extend further into the film than those in fast dynamics-a result with implications for how distinct properties are perturbed near an interface. None of these features track locally with changes in the heterogeneity of caging scale, emphasizing the local disconnect between these quantities near interfaces. These results are at odds with conceptions of the surface as reflecting simply a higher "rheological temperature" than the bulk, instead pointing to a complex interplay between bulk dynamic heterogeneity and spatially organized dynamical gradients at interfaces in glass-forming liquids.

Glass formation and dynamics of model polymer films with one versus two active interfaces

Polymers and other glass-forming liquids can exhibit profound alterations in dynamics in the nanoscale vicinity of interfaces, over a range appreciably exceeding that of typical interfacial thermodynamic gradients. The understanding of these dynamical gradients is particularly complicated in systems with internal or external nanoscale dimensions, where a gradient nucleated at one interface can impinge on a second, potentially distinct, interface. To better understand the interactions that govern system dynamics and glass formation in these cases, here we simulate the baseline case of a glass-forming polymer film, over a wide range of thickness, supported on a dynamically neutral substrate that has little effect on nearby dynamics. We compare these results to our prior simulations of freestanding films. Results indicate that dynamical gradients in our simulated systems, as measured based upon translational relaxation, are simply truncated when they impinge on a secondary surface that is locally dynamically neutral. Altered film behavior can be described almost entirely by gradient effects down to the thinnest films probed, with no evidence for finite-size effects sometimes posited to play a role in these systems. Finally, our simulations predict that linear gradient overlap effects in the presence of symmetric dynamically active interfaces yield a non-monotonic variation of the whole free standing film stretching exponent (relaxation time distribution breadth). The maximum relaxation time distribution breadth in simulation is found at a film thickness of 4-5 times the interfacial gradient range. Observation of this maximum in experiment would provide an important validation that the gradient behavior observed in simulation persists to experimental timescales. If validated, observation of this maximum would potentially also enable determination of the dynamic gradient range from experimental mean-film measurements of film dynamics.

Effect of the polar group content on the glass transition temperature of ROMP copolymers

Polar groups have long been recognized to greatly influence the glass transition temperature (T_g) of polymers, but understanding the underlying physical mechanism remains a challenge. Here, we study the glass formation of ring-opening metathesis polymerization (ROMP) copolymers containing polar groups by employing all-atom molecular dynamics simulations. We show that although the number of hydrogen bonds (N_HB) and the cohesive energy density increase linearly as the content of polar groups (f_pol) increases, the T_g of ROMP copolymers increases with the increase of f_pol in a nonlinear fashion, and tends to plateau for sufficiently high f_pol. Importantly, we find that the increase rate of Gibbs free energy for HB breaking gradually slows down with the increase of f_pol, indicating that the HB is gradually stabilized. Therefore, T_g is jointly determined by N_HB and the strength of HBs in the system, while the latter dominates. Although N_HB increases linearly with increasing f_pol, the HB strength increases slowly with increasing f_pol, which leads to a decreasing rate of increase in T_g.

Mobility gradients yield rubbery surfaces on top of polymer glasses

Many emerging materials, such as ultrastable glasses^1,2 of interest for phone displays and OLED television screens, owe their properties to a gradient of enhanced mobility at the surface of glass-forming liquids. The discovery of this surface mobility enhancement^3-5 has reshaped our understanding of the behaviour of glass formers and of how to fashion them into improved materials. In polymeric glasses, these interfacial modifications are complicated by the existence of a second length scale-the size of the polymer chain-as well as the length scale of the interfacial mobility gradient^6-9. Here we present simulations, theory and time-resolved surface nano-creep experiments to reveal that this two-scale nature of glassy polymer surfaces drives the emergence of a transient rubbery, entangled-like surface behaviour even in polymers comprised of short, subentangled chains. We find that this effect emerges from superposed gradients in segmental dynamics and chain conformational statistics. The lifetime of this rubbery behaviour, which will have broad implications in constraining surface relaxations central to applications including tribology, adhesion, and surface healing of polymeric glasses, extends as the material is cooled. The surface layers suffer a general breakdown in time-temperature superposition (TTS), a fundamental tenet of polymer physics and rheology. This finding may require a reevaluation of strategies for the prediction of long-time properties in polymeric glasses with high interfacial areas. We expect that this interfacial transient elastomer effect and TTS breakdown should normally occur in macromolecular systems ranging from nanocomposites to thin films, where interfaces dominate material properties^5,10.

The microscopic origins of stretched exponential relaxation in two model glass-forming liquids as probed by simulations in the isoconfigurational ensemble

The origin of stretched exponential relaxation in supercooled glass-forming liquids is one of the central questions regarding the anomalous dynamics of these fluids. The dominant explanation for this phenomenon has long been the proposition that spatial averaging over a heterogeneous distribution of locally exponential relaxation processes leads to stretching. Here, we perform simulations of model polymeric and small-molecule glass-formers in the isoconfigurational ensemble to show that stretching instead emerges from a combination of spatial averaging and locally nonexponential relaxation. The results indicate that localities in the fluid exhibiting faster-than-average relaxation tend to exhibit locally stretched relaxation, whereas slower-than-average relaxing domains exhibit more compressed relaxation. We show that local stretching is predicted by loose local caging, as measured by the Debye-Waller factor, and vice versa. This phenomenology in the local relaxation of in-equilibrium glasses parallels the dynamics of out of equilibrium under-dense and over-dense glasses, which likewise exhibit an asymmetry in their degree of stretching vs compression. On the basis of these results, we hypothesize that local stretching and compression in equilibrium glass-forming liquids results from evolution of particle mobilities over a single local relaxation time, with slower particles tending toward acceleration and vice versa. In addition to providing new insight into the origins of stretched relaxation, these results have implications for the interpretation of stretching exponents as measured via metrologies such as dielectric spectroscopy: measured stretching exponents cannot universally be interpreted as a direct measure of the breadth of an underlying distribution of relaxation times.

Breakdown of the Stokes-Einstein Equation for Solutions of Water in Oil Reverse Micelles

An experimental study is presented for the reverse micellar system of 15% by mass polydisperse hexaethylene glycol monodecylether (C₁₀E₆) in cyclohexane with varying amounts of added water up to 4% by mass. Measurements of viscosity and self-diffusion coefficients were taken as a function of temperature between 10 and 45 °C at varying sample water loads but fixed C₁₀E₆/cyclohexane composition. The results were used to inspect the validity of the Stokes-Einstein equation for this system. Unreasonably small reverse average micelle radii and aggregation numbers were obtained with the Stokes-Einstein equation, but reasonable values for these quantities were obtained using the ratio of surfactant-to-cyclohexane self-diffusion coefficients. While bulk viscosity increased with increasing water load, a concurrent expected decrease of self-diffusion coefficient was only observed for the surfactant and water but not for cyclohexane, which showed independence of water load. Moreover, a spread of self-diffusion coefficients was observed for the protons associated with the ethylene oxide repeat unit in samples with polydisperse C₁₀E₆ but not in a sample with monodisperse C₁₀E₆. These findings were interpreted by the presence of reverse micelle to reverse micelle hopping motions that with higher water load become increasingly selective toward C₁₀E₆ molecules with short ethylene oxide repeat units, while those with long ethylene oxide repeat units remain trapped within the reverse micelle because of the increased hydrogen bonding interactions with the water inside the growing core of the reverse micelle. Despite the observed breakdown of the Stokes-Einstein equation, the temperature dependence of the viscosities and self-diffusion coefficients was found to follow Arrhenius behavior over the investigated range of temperatures.

Dynamic Properties of Supercooled Chlorinated Biphenyls

A study of the dynamics of a series of biphenyl compounds having varying chlorine levels was carried out. Increasing the chlorine content increases the glass transition temperature and makes the dynamics substantially more sensitive to density changes. Nonetheless, in the vicinity of their respective glass transitions, the different liquids display very similar extents of dynamic correlation and dynamic heterogeneity. The slight narrowing of the relaxation peak with increasing chlorine follows the general trend of the effect of increasing molecular polarity. This relationship between the peak breadth and dipole moment was reproduced in molecular dynamics simulations of a simplified model of the Aroclor molecule.

Progress towards a phenomenological picture and theoretical understanding of glassy dynamics and vitrification near interfaces and under nanoconfinement

The nature of alterations to dynamics and vitrification in the nanoscale vicinity of interfaces-commonly referred to as "nanoconfinement" effects on the glass transition-has been an open question for a quarter century. We first analyze experimental and simulation results over the last decade to construct an overall phenomenological picture. Key features include the following: after a metrology- and chemistry-dependent onset, near-interface relaxation times obey a fractional power law decoupling relation with bulk relaxation; relaxation times vary in a double-exponential manner with distance from the interface, with an intrinsic dynamical length scale appearing to saturate at low temperatures; the activation barrier and vitrification temperature T_g approach bulk behavior in a spatially exponential manner; and all these behaviors depend quantitatively on the nature of the interface. We demonstrate that the thickness dependence of film-averaged T_g for individual systems provides a poor basis for discrimination between different theories, and thus we assess their merits based on the above dynamical gradient properties. Entropy-based theories appear to exhibit significant inconsistencies with the phenomenology. Diverse free-volume-motivated theories vary in their agreement with observations, with approaches invoking cooperative motion exhibiting the most promise. The elastically cooperative nonlinear Langevin equation theory appears to capture the largest portion of the phenomenology, although important aspects remain to be addressed. A full theoretical understanding requires improved confrontation with simulations and experiments that probe spatially heterogeneous dynamics within the accessible 1-ps to 1-year time window, minimal use of adjustable parameters, and recognition of the rich quantitative dependence on chemistry and interface.

Universal localization transition accompanying glass formation: insights from efficient molecular dynamics simulations of diverse supercooled liquids

The origin of the precipitous dynamic arrest known as the glass transition is a grand open question of soft condensed matter physics. It has long been suspected that this transition is driven by an onset of particle localization and associated emergence of a glassy modulus. However, progress towards an accepted understanding of glass formation has been impeded by an inability to obtain data sufficient in chemical diversity, relaxation timescales, and spatial and temporal resolution to validate or falsify proposed theories for its physics. Here we first describe a strategy enabling facile high-throughput simulation of glass-forming liquids to nearly unprecedented relaxation times. We then perform simulations of 51 glass-forming liquids, spanning polymers, small organic molecules, inorganics, and metallic glass-formers, with longest relaxation times exceeding one microsecond. Results identify a universal particle-localization transition accompanying glass formation across all classes of glass-forming liquid. The onset temperature of non-Arrhenius dynamics is found to serve as a normalizing condition leading to a master collapse of localization data. This transition exhibits a non-universal relationship with dynamic arrest, suggesting that the nonuniversality of supercooled liquid dynamics enters via the dependence of relaxation times on local cage scale. These results suggest that a universal particle-localization transition may underpin the glass transition, and they emphasize the potential for recent theoretical developments connecting relaxation to localization and emergent elasticity to finally explain the origin of this phenomenon. More broadly, the capacity for high-throughput prediction of glass formation behavior may open the door to computational inverse design of glass-forming materials.

Influence of wall heterogeneity on nanoscopically confined polymers

We investigate via molecular dynamics simulations the behaviour of a polymer melt confined between surfaces with increasing spatial correlation (patchiness) of weakly and strongly interacting sites. Beyond a critical patchiness, we find a dramatic dynamic decoupling, characterized by a steep growth of the longest relaxation time and a constant diffusion coefficient. This arises from dynamic heterogeneities induced by the walls in the adjacent polymer layers, leading to the coexistence of fast and slow chain populations. Structural variations are also present, but they are not easy to detect. Our work opens the way to a better understanding of adhesion, friction, rubber reinforcement by fillers, and many other open issues involving the dynamics of polymeric materials on rough, chemically heterogeneous and possibly "dirty" surfaces.