Dynamical heterogeneity in a vapor-deposited polymer glass

Correlating Ultrastability with Fragility and Surface Mobility in Vapor Deposited Tetrahedral Glasses

Several experiments on molecular and metallic glasses have shown that the ability of vapor deposition to produce ultrastable glasses is correlated to their structural and thermodynamic properties. Here we investigate the vapor deposition of a class of tetrahedral materials (including silicon and water) via molecular dynamics simulations of the generalized Stillinger-Weber potential. By changing a single parameter that controls the local tetrahedrality, we show that the emergence of ultrastable behavior is correlated with an increase in the fragility of the model. At the same time, while the mobility of the surface compared to the bulk shows only slight changes at low temperature, with increasing the tetrahedrality, it displays a significant enhancement toward the glass transition temperature. Our results point toward a strong connection between bulk dynamics, surface dynamics, and glass-ultrastability ability in this class of materials.

Entropic timescales of dynamic heterogeneity in supercooled liquid

Non-Gaussian displacement distributions are universal predictors of dynamic heterogeneity in slowly varying environments. Here, we explore heterogeneous dynamics in supercooled liquid using molecular dynamics simulations and show the efficiency of the relative-entropy based measure, negentropy, in quantifying dynamic heterogeneity over the widely used non-Gaussian parameter. Our analysis shows that the heterogeneity quantified by the negentropy is significantly different from the one obtained using the conventional moment-based definition that considers deviation from Gaussianity up to lower-order moments. We extract the timescales of dynamic heterogeneity using the two methods and show that the differential changes diverge as the system experiences strong intermittency near the glass transition. Further, we interpret the entropic timescales and discuss the general implications of our work.

A simple mean-field model of glassy dynamics and glass transition

We propose a phenomenological model to describe the equilibrium dynamic behavior of amorphous glassy materials. It is assumed that a material can be represented by a lattice of cooperatively re-arranging regions (CRRs), with each CRR having two states, the low-temperature "solid" and the high-temperature "liquid". At low temperatures, the material exhibits two characteristic relaxation times, corresponding to the slow large-scale motion between the "solid" CRRs (α-relaxation) and the faster local motion within individual CRRs (β-relaxation). At high temperatures, the α- and β-relaxation times merge, as observed experimentally and suggested by the "Coupling Model" framework. Our new approach is labeled "Two-state, two (time)scale model" or TS2. It is shown that the TS2 treatment can successfully describe the "two-Arrhenius" relaxation time behavior described in several recent experiments. We also apply TS2 to describe the pressure- and molecular-weight dependence of the glass transition temperature in bulk polymers, as well as its dependence on film thickness in thin films.

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.

Spatial Dependence of Non-Gaussian Diffusion of Nanoparticles in Free-Standing Thin Polymer Films

The addition of nanoparticles (NPs) to a free-standing polymer film affects the properties of the film such as viscosity and glass transition temperature. Recent experiments, for example, showed that the glass transition temperature of thin polymer films was dependent on how NPs were distributed within the polymer films. However, the spatial arrangement of NPs in free-standing polymer films and its effect on the diffusion of NPs and polymers remain elusive at a molecular level. In this study, we employ generic coarse-grained models for polymers and NPs and perform extensive molecular dynamics simulations to investigate the diffusion of polymers and NPs in free-standing thin polymer films. We find that small NPs are likely to stay at the interfacial region of the polymer film, while large NPs tend to stay at the center of the film. On the other hand, as the interaction between a NP and a monomer becomes more attractive, the NP is more likely to be placed at the film center. The diffusion of monomers slows down slightly as more NPs are added to the film. Interestingly, the NP diffusion is dependent strongly on the spatial arrangement of the NPs: NPs at the interfacial region diffuse faster and undergo more non-Gaussian diffusion than NPs at the film center, which implies that the interfacial region would be more mobile and dynamically heterogeneous than the film center. We also find that the mechanism for non-Gaussian diffusion of NPs at the film center differs from that at the interfacial region and that the NP diffusion would reflect the local viscosity of the polymer films.

A comparative study on the dynamic heterogeneity of supercooled polymers under nanoconfinement

Dynamic heterogeneity (DH) is a universal property of glass transition phenomena. In this work, we perform a comparative analysis of DH for pure polymer and polymer/nanoparticle composite systems in both film and bulk states via molecular dynamics simulations. We find that the dynamic gradient and the faster average dynamics due to the presence of a free surface are two leading factors, resulting from a nanoconfinement effect, which influence different parts of DH in a film system. The dynamic gradient results from differences in dynamics at different distances from the mobile surface, which induces a large deviation from the Gaussian distribution for the displacement distribution in the film. At the same time, the maximum string size which describes the region size for cooperative motion (dynamic correlation) can also be influenced by the dynamic gradient, although this influence is much weaker than that on the displacement distribution. On the other hand, reflecting temporal fluctuations of dynamics or temporal parts of DH, characteristic peak times of the non-Gaussian parameter and string size, and the ratio between persistent times and exchange times which describe the dynamic exchange properties, are mainly influenced by the faster dynamics on average. Our results demonstrate that measuring different properties (dynamic distribution, dynamic correlation or dynamic exchange) place an emphasis on distinct temporal and spatial parts of DH. It is necessary to use combinational measurements of these properties to give a complete picture of DH in nanoconfinement environments.

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

To understand the role of collective motion in the often large changes in interfacial molecular mobility observed in polymer films, we investigate the extent of collective motion in the interfacial regions of a thin supported polymer film and within the film interior by molecular dynamics simulation. Contrary to commonly stated expectations, we find that the extent of collective motion, as quantified by string-like molecular exchange motion, is similar in magnitude in the polymer-air interfacial layer as the film interior and distinct from the bulk material. This finding is consistent with Adam-Gibbs description of the segmental dynamics within mesoscopic film regions, where the extent of collective motion is related to the configurational entropy of the film as a whole rather than a locally defined extent of collective motion or configurational entropy.

Heterogeneous dynamics of unentangled chains in polymer nanocomposites

We present a systematic investigation on the effect of adding nanoparticles on the dynamics of polymer chains by using coarse-grained molecular dynamics simulation. The dynamics is characterized by three aspects: molecular motion, relaxation at different length scales, and dynamical heterogeneity. It is found that the motion of polymer chains slows down and the deviation from Gaussian distribution becomes more pronounced with increasing nanoparticle volume fractions. For polymer nanocomposites with R ≤ R_g, the relaxation at the wave vector q = 7.0 displays multistep decay, consistent with the previous reports in strongly interacting polymer nanocomposites. Moreover, a qualitatively universal law is established that dynamic heterogeneity at whole chain's scale follows a nonmonotonic increase with increasing nanoparticle loadings, where the volume fraction of the maximum dynamic heterogeneity corresponds to the particle loading when the average distance between nanoparticles is equal to the Kuhn length of polymer chains. We show that the decoupling between whole chain's dynamics and segment dynamics is responsible for the nonmonotonic behavior of dynamic heterogeneity of whole chains.

Comparison of single particle dynamics at the center and on the surface of equilibrium glassy films

Glasses prepared by vapor depositing molecules onto a properly prepared substrate can have enhanced kinetic stability when compared with glasses prepared by cooling from the liquid state. The enhanced stability is due to the high mobility of particles at the surface, which allows them to find lower energy configurations than for liquid cooled glasses. Here we use molecular dynamics simulations to examine the temperature dependence of the single particle dynamics in the bulk of the film and at the surface of the film. First, we examine the temperature dependence of the self-intermediate scattering functions for particles in the bulk and at the surface. We then examine the temperature dependence of the probability of the logarithm of single particle displacements for bulk and surface particles. Both bulk and surface particle displacements indicate populations of slow and fast particles, i.e., heterogeneous dynamics. We find that the temperature dependence of the surface dynamics mirrors the bulk despite being several orders of magnitude faster.

Why we need to look beyond the glass transition temperature to characterize the dynamics of thin supported polymer films

There is significant variation in the reported magnitude and even the sign of [Formula: see text] shifts in thin polymer films with nominally the same chemistry, film thickness, and supporting substrate. The implicit assumption is that methods used to estimate [Formula: see text] in bulk materials are relevant for inferring dynamic changes in thin films. To test the validity of this assumption, we perform molecular simulations of a coarse-grained polymer melt supported on an attractive substrate. As observed in many experiments, we find that [Formula: see text] based on thermodynamic criteria (temperature dependence of film height or enthalpy) decreases with decreasing film thickness, regardless of the polymer-substrate interaction strength ε. In contrast, we find that [Formula: see text] based on a dynamic criterion (relaxation of the dynamic structure factor) also decreases with decreasing thickness when ε is relatively weak, but [Formula: see text] increases when ε exceeds the polymer-polymer interaction strength. We show that these qualitatively different trends in [Formula: see text] reflect differing sensitivities to the mobility gradient across the film. Apparently, the slowly relaxing polymer segments in the substrate region make the largest contribution to the shift of [Formula: see text] in the dynamic measurement, but this part of the film contributes less to the thermodynamic estimate of [Formula: see text] Our results emphasize the limitations of using [Formula: see text] to infer changes in the dynamics of polymer thin films. However, we show that the thermodynamic and dynamic estimates of [Formula: see text] can be combined to predict local changes in [Formula: see text] near the substrate, providing a simple method to infer information about the mobility gradient.

Origin of Ultrastability in Vapor-Deposited Glasses

Glass films created by vapor-depositing molecules onto a substrate can exhibit properties similar to those of ordinary glasses aged for thousands of years. It is believed that enhanced surface mobility is the mechanism that allows vapor deposition to create such exceptional glasses, but it is unclear how this effect is related to the final state of the film. Here we use molecular dynamics simulations to model vapor deposition and an efficient Monte Carlo algorithm to determine the deposition rate needed to create ultrastable glassy films. We obtain a scaling relation that quantitatively captures the efficiency gain of vapor deposition over bulk annealing, and demonstrates that surface relaxation plays the same role in the formation of vapor-deposited glasses as bulk relaxation does in ordinary glass formation.

Relating Ultrastable Glass Formation to Enhanced Surface Diffusion via the Johari-Goldstein β-Relaxation in Molecular Glasses

Glasses are materials essential for modern technology; they are usually prepared by cooling liquids. Recently, novel ultrastable glasses (SGs) with extraordinary thermodynamic and kinetic stability have been created by vapor deposition at appropriate substrate temperatures. However, the underlying mechanism for the formation of SGs is still not established. For most of the molecular SGs created so far, we demonstrate that the formation of SGs is closely related to the Johari-Goldstein β-relaxation from the fact that the lowest substrate temperatures possible for the formation of SGs match the secondary glass-transition temperatures, where the β-relaxation time reaches 10³ s. Theoretically the β-relaxation time via the primitive relaxation time of the coupling model has proven capable of accounting for the enhancement of molecular mobility at the surface. Thus our findings provide evidence to support that the immense enhancement of molecular diffusion at the surface is critical for the formation of SGs. The result has implications in the design and fabrication of SGs.