Dynamical Theory of Segmental Relaxation and Emergent Elasticity in Supercooled Polymer Melts

Polymer-Polymer Interfacial Perturbation on the Glass Transition of Supported Low Molecular Weight Polystyrene Thin Films

Clarifying interfacial perturbation on polymer relaxation is important for polymer material development. Herein we investigated polymer-polymer interfacial perturbation on low molecular weight (MW) polystyrene (PS) thin film (15-180 nm) glass transition by depositing various polymers atop PS films. Overall, rubbery topcoats induced T_g depression of PS thin film (below 60 nm), while glassy topcoats induced T_g elevation of PS thin film (below 30 nm). Importantly, for the rubbery topcoat, T_g perturbation strength is largely dependent on the T_g difference between interfacial polymers and a larger T_g difference would induce stronger perturbation, while for the glassy topcoat this dependence is inconspicuous. Meanwhile, the interfacial perturbation length during PS glass transition by rubbery topcoats is estimated to be around 8 nm, while it is considered to be about 3.5 nm for glassy topcoats. The different interfacial perturbation length induced by disparate topcoats was accounted for by their different perturbation strength on adjacent PS molecules and disparate interfacial roughness. The results can promote the understanding of polymer interfacial perturbation and benefit the design and development of polymer-based materials.

Linear and nonlinear rheology and structural relaxation in dense glassy and jammed soft repulsive pNIPAM microgel suspensions

We present an integrated experimental and quantitative theoretical study of the mechanics of self-crosslinked, slightly charged, repulsive pNIPAM microgel suspensions over a very wide range of concentrations (c) that span the fluid, glassy and putative "soft jammed" regimes. In the glassy regime we measure a linear elastic dynamic shear modulus over 3 decades which follows an apparent power law concentration dependence G' ∼ c5.64, a variation that appears distinct from prior studies of crosslinked ionic microgel suspensions. At very high concentrations there is a sharp crossover to a nearly linear growth of the modulus. To theoretically understand these observations, we formulate an approach to address all three regimes within a single conceptual Brownian dynamics framework. A minimalist single particle description is constructed that allows microgel size to vary with concentration due to steric de-swelling effects. Using a Hertzian repulsion interparticle potential and a suite of statistical mechanical theories, quantitative predictions under quiescent conditions of microgel collective structure, dynamic localization length, elastic modulus, and the structural relaxation time are made. Based on a constant inter-particle repulsion strength parameter which is determined by requiring the theory to reproduce the linear elastic shear modulus over the entire concentration regime, we demonstrate good agreement between theory and experiment. Testable predictions are then made. We also measured nonlinear rheological properties with a focus on the yield stress and strain. A theoretical analysis with no adjustable parameters predicts how the quiescent structural relaxation time changes under deformation, and how the yield stress and strain change as a function of concentration. Reasonable agreement with our observations is obtained. To the best of our knowledge, this is the first attempt to quantitatively understand structure, quiescent relaxation and shear elasticity, and nonlinear yielding of dense microgel suspensions using microscopic force based theoretical methods that include activated hopping processes. We expect our approach will be useful for other soft polymeric particle suspensions in the core-shell family.

Theory of the spatial transfer of interface-nucleated changes of dynamical constraints and its consequences in glass-forming films

We formulate a new theory for how caging constraints in glass-forming liquids at a surface or interface are modified and then spatially transferred, in a layer-by-layer bootstrapped manner, into the film interior in the context of the dynamic free energy concept of the Nonlinear Langevin Equation (NLE) theory approach. The dynamic free energy at any mean location (cage center) involves contributions from two adjacent layers where confining forces are not the same. At the most fundamental level of the theory, the caging component of the dynamic free energy varies essentially exponentially with distance from the interface, saturating deep enough into the film with a correlation length of modest size and weak sensitivity to the thermodynamic state. This imparts a roughly exponential spatial variation of all the key features of the dynamic free energy required to compute gradients of dynamical quantities including the localization length, jump distance, cage barrier, collective elastic barrier, and alpha relaxation time. The spatial gradients are entirely of dynamical, not structural or thermodynamic, origin. The theory is implemented for the hard sphere fluid and diverse interfaces which can be a vapor, a rough pinned particle solid, a vibrating (softened) pinned particle solid, or a smooth hard wall. Their basic description at the level of the spatially heterogeneous dynamic free energy is identical, with the crucial difference arising from the first layer where dynamical constraints can be weakened, softened, or hardly changed depending on the specific interface. Numerical calculations establish the spatial dependence and fluid volume fraction sensitivity of the key dynamical property gradients for five different model interfaces. A comparison of the theoretical predictions for the dynamic localization length and glassy modulus with simulations and experiments for systems with a vapor interface reveals good agreement. The present advance sets the stage for using the Elastically Collective NLE theory to make quantitative predictions for the alpha relaxation time gradient, decoupling phenomena, T_g gradient, and many film-averaged properties of both model and experimental (colloids, molecules, and polymers) systems with diverse interfaces and chemical makeup.

Optimizing the Grafting Density of Tethered Chains to Alter the Local Glass Transition Temperature of Polystyrene near Silica Substrates: The Advantage of Mushrooms over Brushes

We measured the local glass transition temperature T_g(z) of polystyrene (PS) as a function of distance z from a silica substrate with end-grafted chains using fluorescence, where competing effects from the free surface have been avoided to focus only on the influence of the tethered interface. The local T_g(z) increase next to the chain-grafted substrate is found to exhibit a maximum increase of 49 ± 2 K relative to bulk at an optimum grafting density that corresponds to the mushroom-to-brush transition regime. This perturbation to the local T_g(z) dynamics of the matrix is observed to persist out to a distance z ≈ 100-125 nm for this optimum grafting density before bulk T_g is recovered, a distance comparable to that previously observed by Baglay and Roth [J. Chem. Phys. 2017, 146, 203307] for PS next to the higher-T_g polymer polysulfone.

Correlated matrix-fluctuation-mediated activated transport of dilute penetrants in glass-forming liquids and suspensions

We formulate a microscopic, force-level statistical mechanical theory for the activated diffusion of dilute penetrants in dense liquids, colloidal suspensions, and glasses. The approach explicitly and self-consistently accounts for coupling between penetrant hopping and matrix dynamic displacements that actively facilitate the hopping event. The key new ideas involve two mechanistically (at a stochastic trajectory level) coupled dynamic free energy functions for the matrix and spherical penetrant particles. A single dynamic coupling parameter quantifies how much the matrix displaces relative to the penetrant when the latter reaches its transition state which is determined via the enforcement of a temporal causality or coincidence condition. The theory is implemented for dilute penetrants smaller than the matrix particles, with or without penetrant-matrix attractive forces. Model calculations reveal a rich dependence of the penetrant diffusion constant and degree of dynamic coupling on size ratio, volume fraction, and attraction strength. In the absence of attractions, a near exponential decrease of penetrant diffusivity with size ratio over an intermediate range is predicted, in contrast to the much steeper, non-exponential variation if one assumes local matrix dynamical fluctuations are not correlated with penetrant motion. For sticky penetrants, the relative and absolute influence of caging versus physical bond formation is studied. The conditions for a dynamic crossover from the case where a time scale separation between penetrant and matrix activated hopping exists to a "slaved" or "constraint release" fully coupled regime are determined. The particle mixture model is mapped to treat experimental thermal systems and applied to make predictions for the diffusivity of water, toluene, methanol, and oxygen in polyvinylacetate liquids and glasses. The theory agrees well with experiment with values of the penetrant-matrix size ratio close to their chemically intuitive values.

Segment-scale, force-level theory of mesoscopic dynamic localization and entropic elasticity in entangled chain polymer liquids

We develop a segment-scale, force-based theory for the breakdown of the unentangled Rouse model and subsequent emergence of isotropic mesoscopic localization and entropic elasticity in chain polymer liquids in the absence of ergodicity-restoring anisotropic reptation or activated hopping motion. The theory is formulated in terms of a conformational N-dynamic-order-parameter generalized Langevin equation approach. It is implemented using a universal field-theoretic Gaussian thread model of polymer structure and closed at the level of the chain dynamic second moment matrix. The physical idea is that the isotropic Rouse model fails due to the dynamical emergence, with increasing chain length, of time-persistent intermolecular contacts determined by the combined influence of local uncrossability, long range polymer connectivity, and a self-consistent treatment of chain motion and the dynamic forces that hinder it. For long chain melts, the mesoscopic localization length (identified as the tube diameter) and emergent entropic elasticity predictions are in near quantitative agreement with experiment. Moreover, the onset chain length scales with the semi-dilute crossover concentration with a realistic numerical prefactor. Distinctive novel predictions are made for various off-diagonal correlation functions that quantify the full spatial structure of the dynamically localized polymer conformation. As the local excluded volume constraint and/or intrachain bonding spring are softened to allow chain crossability, the tube diameter is predicted to swell until it reaches the radius-of-gyration at which point mesoscopic localization vanishes in a discontinuous manner. A dynamic phase diagram for such a delocalization transition is constructed, which is qualitatively consistent with simulations and the classical concept of a critical entanglement degree of polymerization.

Glass transition of polymers in bulk, confined geometries, and near interfaces

When cooled or pressurized, polymer melts exhibit a tremendous reduction in molecular mobility. If the process is performed at a constant rate, the structural relaxation time of the liquid eventually exceeds the time allowed for equilibration. This brings the system out of equilibrium, and the liquid is operationally defined as a glass-a solid lacking long-range order. Despite almost 100 years of research on the (liquid/)glass transition, it is not yet clear which molecular mechanisms are responsible for the unique slow-down in molecular dynamics. In this review, we first introduce the reader to experimental methodologies, theories, and simulations of glassy polymer dynamics and vitrification. We then analyse the impact of connectivity, structure, and chain environment on molecular motion at the length scale of a few monomers, as well as how macromolecular architecture affects the glass transition of non-linear polymers. We then discuss a revised picture of nanoconfinement, going beyond a simple picture based on interfacial interactions and surface/volume ratio. Analysis of a large body of experimental evidence, results from molecular simulations, and predictions from theory supports, instead, a more complex framework where other parameters are relevant. We focus discussion specifically on local order, free volume, irreversible chain adsorption, the Debye-Waller factor of confined and confining media, chain rigidity, and the absolute value of the vitrification temperature. We end by highlighting the molecular origin of distributions in relaxation times and glass transition temperatures which exceed, by far, the size of a chain. Fast relaxation modes, almost universally present at the free surface between polymer and air, are also remarked upon. These modes relax at rates far larger than those characteristic of glassy dynamics in bulk. We speculate on how these may be a signature of unique relaxation processes occurring in confined or heterogeneous polymeric systems.

Big Effect of Small Nanoparticles: A Shift in Paradigm for Polymer Nanocomposites

Polymer nanocomposites (PNCs) are important materials that are widely used in many current technologies and potentially have broader applications in the future due to their excellent property tunability, light weight, and low cost. However, expanding the limits in property enhancement remains a fundamental scientific challenge. Here, we demonstrate that well-dispersed, small (diameter ∼1.8 nm) nanoparticles with attractive interactions lead to unexpectedly large and qualitatively different changes in PNC structural dynamics in comparison to conventional nanocomposites based on particles of diameters ∼10-50 nm. At the same time, the zero-shear viscosity at high temperatures remains comparable to that of the neat polymer, thereby retaining good processability and resolving a major challenge in PNC applications. Our results suggest that the nanoparticle mobility and relatively short lifetimes of nanoparticle-polymer associations open qualitatively different horizons in the tunability of macroscopic properties in nanocomposites with a high potential for the development of advanced functional materials.

Microscopic Theory for the Role of Attractive Forces in the Dynamics of Supercooled Liquids

We formulate a microscopic, no adjustable parameter, theory of activated relaxation in supercooled liquids directly in terms of the repulsive and attractive forces within the framework of pair correlations. Under isochoric conditions, attractive forces can nonperturbatively modify slow dynamics, but at high enough density their influence vanishes. Under isobaric conditions, attractive forces play a minor role. High temperature apparent Arrhenius behavior and density-temperature scaling are predicted. Our results are consistent with recent isochoric simulations and isobaric experiments on a deeply supercooled molecular liquid. The approach can be generalized to treat colloidal gelation and glass melting, and other soft matter slow dynamics problems.

Theory of activated penetrant diffusion in viscous fluids and colloidal suspensions

We heuristically formulate a microscopic, force level, self-consistent nonlinear Langevin equation theory for activated barrier hopping and non-hydrodynamic diffusion of a hard sphere penetrant in very dense hard sphere fluid matrices. Penetrant dynamics is controlled by a rich competition between force relaxation due to penetrant self-motion and collective matrix structural (alpha) relaxation. In the absence of penetrant-matrix attraction, three activated dynamical regimes are predicted as a function of penetrant-matrix size ratio which are physically distinguished by penetrant jump distance and the nature of matrix motion required to facilitate its hopping. The penetrant diffusion constant decreases the fastest with size ratio for relatively small penetrants where the matrix effectively acts as a vibrating amorphous solid. Increasing penetrant-matrix attraction strength reduces penetrant diffusivity due to physical bonding. For size ratios approaching unity, a distinct dynamical regime emerges associated with strong slaving of penetrant hopping to matrix structural relaxation. A crossover regime at intermediate penetrant-matrix size ratio connects the two limiting behaviors for hard penetrants, but essentially disappears if there are strong attractions with the matrix. Activated penetrant diffusivity decreases strongly with matrix volume fraction in a manner that intensifies as the size ratio increases. We propose and implement a quasi-universal approach for activated diffusion of a rigid atomic/molecular penetrant in a supercooled liquid based on a mapping between the hard sphere system and thermal liquids. Calculations for specific systems agree reasonably well with experiments over a wide range of temperature, covering more than 10 orders of magnitude of variation of the penetrant diffusion constant.