Reconciling computational and experimental trends in the temperature dependence of the interfacial mobility of polymer films

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.

Gradient in refractive index reveals denser near free surface region in thin polymer films

A gradient in refractive index that is linear in magnitude with depth into the film is used to fit ellipsometric data for thin polymer films of poly(methyl methacrylate) (PMMA), polystyrene (PS), and poly(2-vinyl pyridine) (P2VP). We find that the linear gradient model fits provide more physically realistic refractive index values for thin films compared with the commonly used homogeneous Cauchy layer model, addressing recent reports of physically unrealistic density increases. Counter to common expectations of a simple free volume correlation between density and dynamics, we find that the direction of refractive index (density) gradient indicates a higher density near the free surface, which we rationalize based on the observed faster free surface dynamics needed to create vapor deposited stable glasses with optimized denser molecular packings. The magnitude of refractive index gradient is observed to be three times larger for PMMA than for PS films, while P2VP films exhibit a more muted response possibly reflective of a decoupling in free surface and substrate dynamics in systems with strong interfacial interactions.

Polymers under nanoconfinement: where are we now in understanding local property changes?

Polymers are increasingly being used in applications with nanostructured morphologies where almost all polymer molecules are within a few tens to hundreds of nanometers from some interface. From nearly three decades of study on polymers in simplified nanoconfined systems such as thin films, we have come to understand property changes in these systems as arising from interfacial effects where local dynamical perturbations are propagated deeper into the material. This review provides a summary of local glass transition temperature T_g changes near interfaces, comparing across different types of interfaces: free surface, substrate, liquid, and polymer-polymer. Local versus film-average properties in thin films are discussed, making connections to other related property changes, while highlighting several historically important studies. By experimental necessity, most studies are on high enough molecule weight chains to be well entangled, although aspects that connect to lower molecule weight materials are described. Emphasis is made to identify observations and open questions that have yet to be fully understood such as the evidence of long-ranged interfacial effects, finite domain size, interfacial breadth, and chain connectivity.

Dynamic heterogeneity, cooperative motion, and Johari-Goldstein [Formula: see text]-relaxation in a metallic glass-forming material exhibiting a fragile-to-strong transition

We investigate the Johari-Goldstein (JG) [Formula: see text]-relaxation process in a model metallic glass-forming (GF) material ([Formula: see text]), previously studied extensively by both frequency-dependent mechanical measurements and simulation studies devoted to equilibrium properties, by molecular dynamics simulations based on validated and optimized interatomic potentials with the primary aim of better understanding the nature of this universal relaxation process from a dynamic heterogeneity (DH) perspective. The present relatively low temperature and long-time simulations reveal a direct correspondence between the JG [Formula: see text]-relaxation time [Formula: see text] and the lifetime of the mobile particle clusters [Formula: see text], defined as in previous DH studies, a relationship dual to the corresponding previously observed relationship between the [Formula: see text]-relaxation time [Formula: see text] and the lifetime of immobile particle clusters [Formula: see text]. Moreover, we find that the average diffusion coefficient D nearly coincides with [Formula: see text] of the smaller atomic species (Al) and that the 'hopping time' associated with D coincides with [Formula: see text] to within numerical uncertainty, both trends being in accord with experimental studies. This indicates that the JG [Formula: see text]-relaxation is dominated by the smaller atomic species and the observation of a direct relation between this relaxation process and rate of molecular diffusion in GF materials at low temperatures where the JG [Formula: see text]-relaxation becomes the prevalent mode of structural relaxation. As an unanticipated aspect of our study, we find that [Formula: see text] exhibits fragile-to-strong (FS) glass formation, as found in many other metallic GF liquids, but this fact does not greatly alter the geometrical nature of DH in this material and the relation of DH to dynamical properties. On the other hand, the temperature dependence of the DH and dynamical properties, such as the structural relaxation time, can be significantly altered from 'ordinary' GF liquids.

Localization model description of the interfacial dynamics of crystalline Cu and [Formula: see text] metallic glass nanoparticles

Many of the special properties of nanoparticles (NPs) and nanomaterials broadly derive from the significant fraction of particles (atoms, molecules or segments of polymeric molecules) in the NP interfacial region in which the interparticle interactions are characteristically highly anharmonic in comparison to the bulk material. This leads to relatively large mean square particle displacements relative to the material interior, often resulting in a strong increase interfacial mobility and reactivity in both crystalline and glass NPs. The 'Debye-Waller factor', or the mean square particle displacement [Formula: see text] on a ps 'caging' timescale relative to the square of the average interparticle distance [Formula: see text], provides an often experimentally accessible measure of the strength of this anharmonic interaction. The Localization Model (LM) of the dynamics of condensed materials relates this thermodynamic property to the structural relaxation time [Formula: see text], determined from the intermediate scattering function, without any free parameters. Moreover, the LM allows for the prediction of the diffusion coefficient D when combined with the 'decoupling' or Fractional Stokes-Einstein relation linking [Formula: see text] to D. In the current study, we employed classical molecular dynamics simulation to investigate the structural relaxation and diffusion of model [Formula: see text] metallic glass and Cu crystalline NPs with different sizes. As with previous studies validating the LM on model bulk and crystalline materials, and for the interfacial dynamics of thin crystalline and metallic glass films, we find the LM model also describes the interfacial dynamics of model crystalline metal (Cu) and metallic glass ([Formula: see text] NPs to a good approximation, further confirming the generality of the model.

Factors correlating to enhanced surface diffusion in metallic glasses

The enhancement of surface diffusion (D_S) over the bulk (D_V) in metallic glasses (MGs) is well documented and likely to strongly influence the properties of glasses grown by vapor deposition. Here, we use classical molecular dynamics (MD) simulations to identify different factors influencing the enhancement of surface diffusion in MGs. MGs have a simple atomic structure and belong to the category of moderately fragile glasses that undergo pronounced slowdown of bulk dynamics with cooling close to the glass transition temperature (T_g). We observe that D_S exhibits a much more moderate slowdown compared to D_V when approaching T_g, and D_S/D_V at T_g varies by two orders of magnitude among the MGs investigated. We demonstrate that both the surface energy and the fraction of missing bonds for surface atoms show good correlation to D_S/D_V, implying that the loss of nearest neighbors at the surface directly translates into higher mobility, unlike the behavior of network-bonded and hydrogen-bonded organic glasses. Fragility, a measure of the slowdown of bulk dynamics close to T_g, also correlates to D_S/D_V, with more fragile systems having larger surface enhancement of mobility. The deviations observed in the fragility-D_S/D_V relationship are shown to be correlated to the extent of segregation or depletion of the mobile element at the surface. Finally, we explore the relationship between the diffusion pre-exponential factor (D₀) and the activation energy (Q) and compare it to a ln(D₀)-Q correlation previously established for bulk glasses, demonstrating similar correlations from MD as in the experiments and that the surface and bulk have very similar ln(D₀)-Q correlations.

Localization model description of the interfacial dynamics of crystalline Cu and Cu64Zr36 metallic glass films

Recent studies of structural relaxation in Cu-Zr metallic glass materials having a range of compositions and over a wide range of temperatures and in crystalline UO₂ under superionic conditions have indicated that the localization model (LM) can predict the structural relaxation time τ_α of these materials from the intermediate scattering function without any free parameters from the particle mean square displacement ⟨r²⟩ at a caging time on the order of ps, i.e., the "Debye-Waller factor" (DWF). In the present work, we test whether this remarkable relation between the "fast" picosecond dynamics and the rate of structural relaxation τ_α in these model amorphous and crystalline materials can be extended to the prediction of the local interfacial dynamics of model amorphous and crystalline films. Specifically, we simulate the free-standing amorphous Cu₆₄Zr₃₆ and crystalline Cu films and find that the LM provides an excellent parameter-free prediction for τ_α of the interfacial region. We also show that the Tammann temperature, defining the initial formation of a mobile interfacial layer, can be estimated precisely for both crystalline and glass-forming solid materials from the condition that the DWFs of the interfacial region and the material interior coincide.