Collective Motion in the Interfacial and Interior Regions of Supported Polymer Films and Its Relation to Relaxation

Effect of the nature of the solid substrate on spatially heterogeneous activated dynamics in glass forming supported films

We extend the force-level elastically collective nonlinear Langevin equation theory to treat the spatial gradients of the alpha relaxation time and glass transition temperature, and the corresponding film-averaged quantities, to the geometrically asymmetric case of finite thickness supported films with variable fluid-substrate coupling. The latter typically nonuniversally slows down motion near the solid-liquid interface as modeled via modification of the surface dynamic free energy caging constraints that are spatially transferred into the film and which compete with the accelerated relaxation gradient induced by the vapor interface. Quantitative applications to the foundational hard sphere fluid and a polymer melt are presented. The strength of the effective fluid-substrate coupling has very large consequences for the dynamical gradients and film-averaged quantities in a film thickness and thermodynamic state dependent manner. The interference of the dynamical gradients of opposite nature emanating from the vapor and solid interfaces is determined, including the conditions for the disappearance of a bulk-like region in the film center. The relative importance of surface-induced modification of local caging vs the generic truncation of the long range collective elastic component of the activation barrier is studied. The conditions for the accuracy and failure of a simple superposition approximation for dynamical gradients in thin films are also determined. The emergence of near substrate dead layers, large gradient effects on film-averaged response functions, and a weak non-monotonic evolution of dynamic gradients in thick and cold films are briefly discussed. The connection of our theoretical results to simulations and experiments is briefly discussed, as is the extension to treat more complex glass-forming systems under nanoconfinement.

Confinement effect of inter-arm interactions on glass formation in star polymer melts

We utilized molecular dynamic simulation to investigate the glass formation of star polymer melts in which the topological complexity is varied by altering the number of star arms (f). Emphasis was placed on how the "confinement effect" of repulsive inter-arm interactions within star polymers influences the thermodynamics and dynamics of star polymer melts. All the characteristic temperatures of glass formation were found to progressively increase with increasing f, but unexpectedly the fragility parameter KVFT was found to decrease with increasing f. As previously observed, stars having more than 5 or 6 arms adopt an average particle-like structure that is more contracted relative to the linear polymer size having the same mass and exhibit a strong tendency for intermolecular and intramolecular segregation. We systematically analyzed how varying f alters collective particle motion, dynamic heterogeneity, the decoupling exponent ζ phenomenologically linking the slow β- and α-relaxation times, and the thermodynamic scaling index γt. Consistent with our hypothesis that the segmental dynamics of many-arm star melts and thin supported polymer films should exhibit similar trends arising from the common feature of high local segmental confinement, we found that ζ increases considerably with increasing f, as found in supported polymer films with decreasing thickness. Furthermore, increasing f led to greatly enhanced elastic heterogeneity, and this phenomenon correlates strongly with changes in ζ and γt. Our observations should be helpful in building a more rational theoretical framework for understanding how molecular topology and geometrical confinement influence the dynamics of glass-forming materials more broadly.

Ejected microcrystals probe jammed states of droplets in cyclodextrin-based emulsions

The cyclodextrin (CD)-based emulsions exhibit complex instability behaviors such as rapid flocculation and creaming, and how to capture droplet dispersion states of the emulsions remains a great challenge. Here we prepare the CD-based emulsions with different oil-water volume ratios and CD concentrations by using high-pressure homogenization, and characterize the emulsion droplets by using optical microscopy and confocal laser scanning microscopy. We evaluate the effects of homogenization pressure on the stability of the emulsions, identify armored droplets with different surface features, measure interfacial concentrations of adsorbed ICs microcrystals, and observe ejection of the oil/CD inclusion complexes (ICs) microcrystals from the droplet surface. The droplet dispersion states are sensitive to the dynamic buildup and evolving morphologies of the interfacial microcrystals, and there are clear correlations between the properties of the ejected microcrystals and the characteristics of the emulsions. We ascribe the subsequent ejection of ICs microcrystals from the droplet surface to consolidation and deformation of the films formed between neighboring droplets. The ejection of the ICs microcrystals affords a simple method to detect the droplet-droplet interactions and phase transitions in the CD-based emulsions, which might be a generic feature in the broader context of the creaming processes of emulsions.

Explaining the Sensitivity of Polymer Segmental Relaxation to Additive Size Based on the Localization Model

We use molecular simulations to examine how the dynamics of a coarse-grained polymer melt are altered by additives of variable size and interaction strength with the polymer matrix. The effect of diluent size σ on polymer dynamics changes significantly when its size is comparable to the polymer segment size. For each σ, we show that the localization model (LM) quantitatively describes the dependence of the segmental relaxation time τ on temperature T in terms of dynamic free volume, quantified by the Debye-Waller factor ⟨u^{2}⟩. Within this model, we show that the additive size alone controls the functional form of the T dependence. The LM parameters reach asymptotic values when the diluent size exceeds the monomer size, converging to a limit applicable to macroscopic interfaces.

Activation free energy gradient controls interfacial mobility gradient in thin polymer films

We examine the mobility gradient in the interfacial region of substrate-supported polymer films using molecular dynamics simulations and interpret these gradients within the string model of glass-formation. No large gradients in the extent of collective motion exist in these simulated films, and an analysis of the mobility gradient on a layer-by-layer basis indicates that the string model provides a quantitative description of the relaxation time gradient. Consequently, the string model indicates that the interfacial mobility gradient derives mainly from a gradient in the high-temperature activation enthalpy ΔH₀ and entropy ΔS₀ as a function of depth z, an effect that exists even in the high-temperature Arrhenius relaxation regime far above the glass transition temperature. To gain insight into the interfacial mobility gradient, we examined various material properties suggested previously to influence ΔH₀ in condensed materials, including density, potential and cohesive energy density, and a local measure of stiffness or u²(z)^-3/2, where u²(z) is the average mean squared particle displacement at a caging time (on the order of a ps). We find that changes in local stiffness best correlate with changes in ΔH₀(z) and that ΔS₀(z) also contributes significantly to the interfacial mobility gradient, so it must not be neglected.

Energy barriers and cooperative motion at the surface of freestanding glassy polystyrene films

We investigate the near-surface relaxation of freestanding atactic polystyrene films with molecular dynamics simulations. As in previous coarse-grained simulations, relaxation times for backbone segments and phenyl rings are linked to their bulk relaxation times via a power-law coupling relation. Variation of the coupling exponent with distance from the surface is consistent with depth-dependent activation barriers. We also quantify a reduction in dynamical heterogeneity at the interface, which can be interpreted in the framework of cooperative models for glassy dynamics.

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.

What does the instantaneous normal mode spectrum tell us about dynamical heterogeneity in glass-forming fluids?

We examine the instantaneous normal mode spectrum of model metallic and polymeric glass-forming liquids. We focus on the localized modes in the unstable part of the spectrum [unstable localized (UL) modes] and find that the particles making the dominant contribution to the participation ratio form clusters that grow upon cooling in a fashion similar to the dynamical heterogeneity in glass-forming fluids, i.e., highly mobile (or immobile) particles form clusters that grow upon cooling; however, a comparison of the UL mode clusters to the mobile and immobile particle clusters indicates that they are distinct entities. We also show that the cluster size provides an alternate method to distinguish localized and delocalized modes, offering a significant practical advantage over the finite-size scaling approach. We examine the trajectories of particles contributing most to the UL modes and find that they have a slightly enhanced mobility compared to the average, and we determine a characteristic time quantifying the persistence time of this excess mobility. This time scale is proportional to the structural relaxation time τ_α of the fluid, consistent with a prediction by Zwanzig [Phys. Rev. 156, 190 (1967)] for the lifetime of collective excitations in cooled liquids. Evidently, these collective excitations serve to facilitate relaxation but do not actually participate in the motion associated with barrier crossing events governing activated transport. They also serve as a possible concrete realization of the "facilitation" clusters postulated in previous modeling of glass-forming liquids.

Universal nature of dynamic heterogeneity in glass-forming liquids: A comparative study of metallic and polymeric glass-forming liquids

Glass-formation is a ubiquitous phenomenon that is often observed in a broad class of materials ranging from biological matter to commonly encountered synthetic polymer, as well as metallic and inorganic glass-forming (GF) materials. Despite the many regularities in the dynamical properties of GF materials, the structural origin of the universal dynamical properties of these materials has not yet been identified. Recent simulations of coarse-grained polymeric GF liquids have indicated the coexistence of clusters of mobile and immobile particles that appear to be directly linked, respectively, to the rate of molecular diffusion and structural relaxation. The present work examines the extent to which these distinct types of "dynamic heterogeneity" (DH) arise in metallic GF liquids (Cu-Zr, Ni-Nb, and Pd-Si alloys) having a vastly different molecular structure and chemistry. We first identified mobile and immobile particles and their transient clusters and found the DH in the metallic alloys to be remarkably similar in form to polymeric GF liquids, confirming the "universality" of the DH phenomenon. Furthermore, the lifetime of the mobile particle clusters was found to be directly related to the rate of diffusion in these materials, while the lifetime of immobile particles was found to be proportional to the structural relaxation time, providing some insight into the origin of decoupling in GF liquids. An examination of particles having a locally preferred atomic packing, and clusters of such particles, suggests that there is no one-to-one relation between these populations of particles so that an understanding of the origin of DH in terms of static fluid structure remains elusive.