Theory of activated glassy relaxation, mobility gradients, surface diffusion, and vitrification in free standing thin films

Importance of Experiments That Can Test Theories Critically

General dynamic and thermodynamic properties of complex materials, including amorphous polymers and molecular glass-formers, have been established from the wealth of experimental data accumulated over the years. Naturally, these general properties attract researchers to construct theories and models to address and explain them. Often more than one theory with contrasting or even conflicting theoretical bases can equally explain a general property rather well. The correct explanation becomes unclear, and progress is stopped. The resolution of the problem comes when an innovative experiment is performed with insightful results that can critically test the premise and assumptions of each theory. This important role played by experimentalists is exemplified by the contributions of Mark Ediger in several general properties considered in this paper: (1) dynamics of the components in binary polymer blends; (2) breakdown of the Stokes-Einstein and the Debye-Stokes-Einstein relations; (3) enhancement of surface mobility and in relation to formation of ultrastable glasses; and (4) the Johari-Goldstein β-relaxation in ultrastable glasses. Different theories proposed to explain these properties are discussed, including the Coupling Model of the author.

Effects of Polymers on Cocrystal Growth in a Drug-Drug Coamorphous System: Relations between Glass-to-Crystal Growth and Surface-Enhanced Crystal Growth

Coamorphous and cocrystal drug delivery systems provide attractive crystal engineering strategies for improving the solubilities, dissolution rates, and oral bioavailabilities of poorly water-soluble drugs. Polymeric additives have often been used to inhibit the unwanted crystallization of amorphous drugs. However, the transformation of a coamorphous phase to a cocrystal phase in the presence of polymers has not been fully elucidated. Herein, we investigated the effects of low concentrations of the polymeric excipients poly(ethylene oxide) (PEO) and poly(vinylpyrrolidone) (PVP) on the growth of carbamazepine-celecoxib (CBZ-CEL) cocrystals from the corresponding coamorphous phase. PEO accelerated the growth rate of the cocrystals by increasing the molecular mobility of the coamorphous system, while PVP had the opposite effect. The coamorphous CBZ-CEL system exhibited two anomalously fast crystal growth modes: glass-to-crystal (GC) growth in the bulk and accelerated crystal growth at the free surface. These two fast growth modes both disappeared after doping with PEO (1-3% w/w) but were retained in the presence of PVP, indicating a potential correlation between the two fast crystal growth modes. We propose that the different effects of PEO and PVP on the crystal growth modes arose from weaker effects of the polymers on cocrystallization at the surface than in the bulk. This work provides a deep understanding of the mechanisms by which polymers influence the cocrystallization kinetics of a multicomponent amorphous phase and highlights the importance of polymer selection in stabilizing coamorphous systems or preparing cocrystals via solid-based methods.

Surface mobility gradient and emergent facilitation in glassy films

Confining glassy polymers into films can substantially modify their local and film-averaged properties. We present a lattice model of film geometry with void-mediated facilitation behaviors but free from any elasticity effect. We analyze the spatially varying viscosity to delineate the transport properties of glassy films. The film mobility measurements reported by Yang et al., Science, 2010, 328, 1676 are successfully reproduced. The flow exhibits a crossover from a simple viscous flow to a surface-dominated regime as the temperature decreases. The propagation of a highly mobile front induced by the free surface is visualized in real space. Our approach provides a microscopic treatment of the observed glassy phenomena.

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.

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.

Understanding creep suppression mechanisms in polymer nanocomposites through machine learning

While recent efforts have shown how local structure plays an essential role in the dynamic heterogeneity of homogeneous glass-forming materials, systems containing interfaces such as thin films or composite materials remain poorly understood. It is known that interfaces perturb the molecular packing nearby, however, numerous studies show the dynamics are modified over a much larger range. Here, we examine the dynamics in polymer nanocomposites (PNCs) using a combination of simulations and experiments and quantitatively separate the role of polymer packing from other effects on the dynamics, as a function of distance from the nanoparticle surfaces. After showing good qualitative agreement between the simulations and experiments in glassy structure and creep compliance, we use a machine-learned structure indicator, softness, to decompose polymer dynamics in our simulated PNCs into structure-dependent and structure-independent processes. With this decomposition, the free energy barrier for polymer rearrangement can be described as a combination of packing-dependent and packing-independent barriers. We find both barriers are higher near nanoparticles and decrease with applied stress, quantitatively demonstrating that the slow interfacial dynamics is not solely due to polymer packing differences, but also the change of structure-dynamics relationships. Finally, we present how this decomposition can be used to accurately predict strain-time creep curves for PNCs from their static configuration, providing additional insights into the effects of polymer-nanoparticle interfaces on creep suppression in PNCs.

Crystallization of Amorphous Nimesulide: The Relationship between Crystal Growth Kinetics and Liquid Dynamics

Understanding crystallization and its correlations with liquid dynamics is relevant for developing robust amorphous pharmaceutical solids. Herein, nimesulide, a classical anti-inflammatory agent, was used as a model system for studying the correlations between crystallization kinetics and molecular dynamics. Kinetic parts of crystal growth (u_kin) of nimesulide exhibited a power law dependence upon the liquid viscosity (η) as u_kin~η^-0.61. Bulk molecular diffusivities (D_Bulk) of nimesulide were predicted by a force-level statistical-mechanical model from the α-relaxation times, which revealed the relationship as u_kin~D_bulk^0.65. Bulk crystal growth kinetics of nimesulide in deeply supercooled liquid exhibited a fragility-dependent decoupling from τ_α. The correlations between growth kinetics and α-relaxation times predicted by the Adam-Gibbs-Vogel equation in a glassy state were also explored, for both the freshly made and fully equilibrated glass. These findings are relevant for the in-depth understanding and prediction of the physical stability of amorphous pharmaceutical solids.

Physical stability and dissolution behaviors of amorphous pharmaceutical solids: Role of surface and interface effects

Amorphous pharmaceutical solids (APS) are single- or multi-component systems in which drugs exist in high-energy states with long-range disordered molecular packing. APSs have become one of the most effective and widely used pharmaceutical delivery approaches for poorly water-soluble drugs in the last several decades. Considerable efforts have been made to investigate the physical stability and dissolution behaviors of APSs, however, the underlying mechanisms remain imperfectly understood. Recent studies reveal that surface and interface properties of APSs could strongly affect the physical stability and dissolution behaviors. This paper provides a comprehensive overview of recent studies focusing on the physical stability and dissolution behaviors of APSs from both surface and interface perspectives. We highlight the role of surface or interface properties in nucleation, crystal growth, phase separation, dissolution, and supersaturation. Meanwhile, the challenges and scope of research on surface and interface properties in the future are also briefly discussed. This review contributes to a better understanding of the surface- and interface-facilitated processes, which will provide more efficient and rational guidance for the design of APSs.

Surface diffusion of a glassy discotic organic semiconductor and the surface mobility gradient of molecular glasses

Surface diffusion has been measured in the glass of an organic semiconductor, MTDATA, using the method of surface grating decay. The decay rate was measured as a function of temperature and grating wavelength, and the results indicate that the decay mechanism is viscous flow at high temperatures and surface diffusion at low temperatures. Surface diffusion in MTDATA is enhanced by 4 orders of magnitude relative to bulk diffusion when compared at the glass transition temperature T_g. The result on MTDATA has been analyzed along with the results on other molecular glasses without extensive hydrogen bonds. In total, these systems cover a wide range of molecular geometries from rod-like to quasi-spherical to discotic and their surface diffusion coefficients vary by 9 orders of magnitude. We find that the variation is well explained by the existence of a steep surface mobility gradient and the anchoring of surface molecules at different depths. Quantitative analysis of these results supports a recently proposed double-exponential form for the mobility gradient: log D(T, z) = log D_v(T) + [log D₀ - log D_v(T)]exp(-z/ξ), where D(T, z) is the depth-dependent diffusion coefficient, D_v(T) is the bulk diffusion coefficient, D₀ ≈ 10^-8 m²/s, and ξ ≈ 1.5 nm. Assuming representative bulk diffusion coefficients for these fragile glass formers, the model reproduces the presently known surface diffusion rates within 0.6 decade. Our result provides a general way to predict the surface diffusion rates in molecular glasses.

Surface Diffusion Is Controlled by Bulk Fragility across All Glass Types

Surface diffusion is vastly faster than bulk diffusion in some glasses, but only moderately enhanced in others. We show that this variation is closely linked to bulk fragility, a common measure of how quickly dynamics is excited when a glass is heated to become a liquid. In fragile molecular glasses, surface diffusion can be a factor of 10^{8} faster than bulk diffusion at the glass transition temperature, while in the strong system SiO_{2}, the enhancement is a factor of 10. Between these two extremes lie systems of intermediate fragility, including metallic glasses and amorphous selenium and silicon. This indicates that stronger liquids have greater resistance to dynamic excitation from bulk to surface and enables prediction of surface diffusion, surface crystallization, and formation of stable glasses by vapor deposition.

Nature of dynamic gradients, glass formation, and collective effects in ultrathin freestanding films

Molecular, polymeric, colloidal, and other classes of liquids can exhibit very large, spatially heterogeneous alterations of their dynamics and glass transition temperature when confined to nanoscale domains. Considerable progress has been made in understanding the related problem of near-interface relaxation and diffusion in thick films. However, the origin of "nanoconfinement effects" on the glassy dynamics of thin films, where gradients from different interfaces interact and genuine collective finite size effects may emerge, remains a longstanding open question. Here, we combine molecular dynamics simulations, probing 5 decades of relaxation, and the Elastically Cooperative Nonlinear Langevin Equation (ECNLE) theory, addressing 14 decades in timescale, to establish a microscopic and mechanistic understanding of the key features of altered dynamics in freestanding films spanning the full range from ultrathin to thick films. Simulations and theory are in qualitative and near-quantitative agreement without use of any adjustable parameters. For films of intermediate thickness, the dynamical behavior is well predicted to leading order using a simple linear superposition of thick-film exponential barrier gradients, including a remarkable suppression and flattening of various dynamical gradients in thin films. However, in sufficiently thin films the superposition approximation breaks down due to the emergence of genuine finite size confinement effects. ECNLE theory extended to treat thin films captures the phenomenology found in simulation, without invocation of any critical-like phenomena, on the basis of interface-nucleated gradients of local caging constraints, combined with interfacial and finite size-induced alterations of the collective elastic component of the structural relaxation process.

Using Deposition Rate and Substrate Temperature to Manipulate Liquid Crystal-Like Order in a Vapor-Deposited Hexagonal Columnar Glass

We investigate vapor-deposited glasses of a phenanthroperylene ester, known to form an equilibrium hexagonal columnar phase, and show that liquid crystal-like order can be manipulated by the choice of deposition rate and substrate temperature during deposition. We find that rate-temperature superposition (RTS)-the equivalence of lowering the deposition rate and increasing the substrate temperature-can be used to predict and control the molecular orientation in vapor-deposited glasses over a wide range of substrate temperatures (0.75-1.0 T_g). This work extends RTS to a new structural motif, hexagonal columnar liquid crystal order, which is being explored for organic electronic applications. By several metrics, including the apparent average face-to-face nearest-neighbor distance, physical vapor deposition (PVD) glasses of the phenanthroperylene ester are as ordered as the glass prepared by cooling the equilibrium liquid crystal. By other measures, the PVD glasses are less ordered than the cooled liquid crystal. We explain the difference in the maximum attainable order with the existence of a gradient in molecular mobility at the free surface of a liquid crystal and its impact upon different mechanisms of structural rearrangement. This free surface equilibration mechanism explains the success of the RTS principle and provides guidance regarding the types of order most readily enhanced by vapor deposition. This work extends the applicability of RTS to include molecular systems with a diverse range of higher-order liquid-crystalline morphologies that could be useful for new organic electronic applications.

Diffusion-controlled and `diffusionless' crystal growth: relation between liquid dynamics and growth kinetics of griseofulvin

Understanding how liquid dynamics govern crystallization is critical for maintaining the physical stability of amorphous pharmaceutical formulations. In the present study, griseofulvin (GSF), a classic antifungal drug, was used as the model system to investigate the correlations between crystal growth kinetics and liquid dynamics. The temperature dependence of the kinetic part of the bulk crystal growth in a supercooled liquid of GSF was weaker than that of the structural relaxation time τα and scaled as τα −0.69. In the glassy state, GSF exhibited the glass-to-crystal (GC) growth behavior, whose growth rate was too fast to be under the control of the α-relaxation process. Moreover, from the perspective of τα, the GC growth of GSF also satisfied the general condition for GC growth to exist: D/u < 7 pm, where D is the diffusion coefficient and u the speed of crystal growth. Also compared were the fast surface crystal growth rates u s and surface relaxation times τsurface predicted by the random first-order transition theory. Here, the surface crystal growth rate u s of GSF exhibited a power-law dependence upon the surface structural relaxation time: u s ∝ τsurface −0.71, which was similar to that of the bulk growth rate and τα. These findings are important for understanding and predicting the crystallization of amorphous pharmaceutical solids both in the bulk and at the surface.

Theory of microstructure-dependent glassy shear elasticity and dynamic localization in melt polymer nanocomposites

We present an integrated theoretical study of the structure, thermodynamic properties, dynamic localization, and glassy shear modulus of melt polymer nanocomposites (PNCs) that spans the three microstructural regimes of entropic depletion induced nanoparticle (NP) clustering, discrete adsorbed layer driven NP dispersion, and polymer-mediated bridging network. The evolution of equilibrium and dynamic properties with NP loading, total packing fraction, and strength of interfacial attraction is systematically studied based on a minimalist model. Structural predictions of polymer reference interaction site model integral equation theory are employed to establish the rich behavior of the interfacial cohesive force density, surface excess, and a measure of free volume as a function of PNC variables. The glassy dynamic shear modulus is predicted to be softened, reinforced, or hardly changed relative to the pure polymer melt depending on system parameters, as a result of the competing and qualitatively different influences of interfacial cohesion (physical bonding), free volume, and entropic depletion on dynamic localization and shear elasticity. The localization of polymer segments is the dominant factor in determining bulk PNC softening and reinforcement effects for moderate to strong interfacial attractions, respectively. While in the athermal entropy-dominated regime, the primary origin of mechanical reinforcement is the stress stored in the aggregated NP subsystem. The PNC shear modulus is often qualitatively correlated with the segment localization length but with notable exceptions. The present work provides the foundation for developing a theory of segmental relaxation, T_g changes, and collective NP dynamics in PNCs based on a self-consistent treatment of the cooperative activated motions of segments and NPs.

Surface diffusion in glasses of rod-like molecules posaconazole and itraconazole: effect of interfacial molecular alignment and bulk penetration

The method of surface grating decay has been used to measure surface diffusion in the glasses of two rod-like molecules posaconazole (POS) and itraconazole (ITZ). Although structurally similar antifungal medicines, ITZ forms liquid-crystalline phases while POS does not. Surface diffusion in these systems is significantly slower than in the glasses of quasi-spherical molecules of similar volume when compared at the glass transition temperature T_g. Between the two systems, ITZ has slower surface diffusion. These results are explained on the basis of the near-vertical orientation of the rod-like molecules at the surface and their deep penetration into the bulk where mobility is low. For molecular glasses without extensive hydrogen bonds, we find that the surface diffusion coefficient at T_g decreases smoothly with the penetration depth of surface molecules and the trend has the double-exponential form for the surface mobility gradient observed in simulations. This supports the view that these molecular glasses have a similar mobility vs. depth profile and their different surface diffusion rates arise simply from the different depths at which molecules are anchored. Our results also provide support for a previously observed correlation between the rate of surface diffusion and the fragility of the bulk liquid.

Temperature-Dependent Spectroscopic Ellipsometry of Thin Polymer Films

Thin polymer films have found many important applications in organic electronics, such as active layers, protective layers, or antistatic layers. Among the various experimental methods suitable for studying the thermo-optical properties of thin polymer films, temperature-dependent spectroscopic ellipsometry plays a special role as a nondestructive and very sensitive optical technique. In this Review Article, issues related to the physical origin of the dependence of ellipsometric angles on temperature are surveyed. In addition, the Review Article discusses the use of temperature-dependent spectroscopic ellipsometry for studying phase transitions in thin polymer films. The benefits of studying thermal transitions using different cooling/heating speeds are also discussed. Furthermore, it is shown how the analysis and modeling of raw ellipsometric data can be used to determine the thermal properties of thin polymer films.

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.

A simulation study on the glass transition behavior and relevant segmental dynamics in free-standing polymer nanocomposite films

In polymer/nanoparticle composite (PNC) thin films, polymer chains experience strong confinement effects not only at the free surface area but also from nanoparticles (NPs). In this work, the influence of NP-polymer interaction and NP distribution on the polymer segmental dynamics and the glass transition behavior of PNC free-standing films are investigated through molecular dynamics simulations. We demonstrate that NPs will migrate to the film surface area and form an NP-concentrated layer when NP-polymer interactions are weak, while NPs are well dispersed in the bulk region when NP-polymer interactions are strong. In both cases, we find increases in the glass transition temperature Tg compared with the pure film without NPs, although with a different degree. The weakly interacting system has the same Tg as the pure bulk system without NPs. The NP layer formed at the surface area reduces both the mobility of the surface polymer beads and the mobility gradient in the film normal direction (MGFND), therefore resulting in an increase in the Tg which highlights the vital role of the mobile surface layer. In contrast, the NPs in the bulk region enlarge the MGFND. NPs have opposite influences on the polymer bead dynamic anisotropy when they interact weakly or strongly with polymers, weakened for the former and enhanced for the latter. These findings offer a clear picture of the segmental dynamics and glass transition behavior in free-standing PNC films with different NP-polymer interaction strengths. We hope these results will be helpful for the property design of related materials.

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.

Temperature-Independent Rescaling of the Local Activation Barrier Drives Free Surface Nanoconfinement Effects on Segmental-Scale Translational Dynamics near Tg

Near-interface alterations in dynamics and glass formation behavior have been the subject of extensive study for the past two decades, both because of their practical importance and in the hope of revealing underlying correlation lengths underpinning glass transition more generally. Here we employ molecular dynamics simulations of thick films to demonstrate that these effects emerge, for segmental-scale translational dynamics at low temperature, from a temperature-independent rescaling of the local activation barrier. This rescaling manifests as a fractional power law decoupling relationship of local dynamics relative to the bulk, with a transition from a regime of weak decoupling at high temperatures to a regime of strong decoupling at low temperatures. The range of this effect saturates at low temperatures, with 90% of the surface perturbation in the barrier lost over a range of 12 segmental diameters. These findings reduce the phenomenology of T_g nanoconfinement effects to two properties-a position-dependent, temperature independent, barrier rescaling factor and an onset time scale-while substantially constraining the predictions required from any theoretical explanation of this phenomenon.

Structural Origin of Enhanced Dynamics at the Surface of a Glassy Alloy

The enhancement of mobility at the surface of an amorphous alloy is studied using a combination of molecular dynamic simulations and normal mode analysis of the nonuniform distribution of Debye-Waller factors. The increased mobility at the surface is found to be associated with the appearance of Arrhenius temperature dependence. We show that the transverse Debye-Waller factor exhibits a peak at the surface. Over the accessible temperature range, we find that the bulk and surface diffusion coefficients obey the same empirical relationship with the respective Debye-Waller factors. Extrapolating this relationship to lower T, we argue that the observed decrease in the constraint at the surface is sufficient to account for the experimentally observed surface enhancement of mobility.

Disconnecting structure and dynamics in glassy thin films

Nanometrically thin glassy films depart strikingly from the behavior of their bulk counterparts. We investigate whether the dynamical differences between a bulk and thin film polymeric glass former can be understood by differences in local microscopic structure. Machine learning methods have shown that local structure can serve as the foundation for successful, predictive models of particle rearrangement dynamics in bulk systems. By contrast, in thin glassy films, we find that particles at the center of the film and those near the surface are structurally indistinguishable despite exhibiting very different dynamics. Next, we show that structure-independent processes, already present in bulk systems and demonstrably different from simple facilitated dynamics, are crucial for understanding glassy dynamics in thin films. Our analysis suggests a picture of glassy dynamics in which two dynamical processes coexist, with relative strengths that depend on the distance from an interface. One of these processes depends on local structure and is unchanged throughout most of the film, while the other is purely Arrhenius, does not depend on local structure, and is strongly enhanced near the free surface of a film.

Influence of Hydrogen Bonding on the Surface Diffusion of Molecular Glasses: Comparison of Three Triazines

Surface grating decay measurements have been performed on three closely related molecular glasses to study the effect of intermolecular hydrogen bonds on surface diffusion. The three molecules are derivatives of bis(3,5-dimethyl-phenylamino)-1,3,5-triazine and differ only in the functional group R at the 2-position, with R being C₂H₅, OCH₃, and NHCH₃, and referred to as "Et", "OMe", and "NHMe", respectively. Of the three molecules, NHMe forms more extensive intermolecular hydrogen bonds than Et and OMe and was found to have slower surface diffusion. For Et and OMe, surface diffusion is so fast that it replaces viscous flow as the mechanism of surface grating decay as temperature is lowered. In contrast, no such transition was observed for NHMe under the same conditions, indicating significantly slower surface diffusion. This result is consistent with the previous finding that extensive intermolecular hydrogen bonds slow down surface diffusion in molecular glasses and is attributed to the persistence of hydrogen bonds even in the surface environment. This result is also consistent with the lower stability of the vapor-deposited glass of NHMe relative to those of Et and OMe and supports the view that surface mobility controls the stability of vapor-deposited glasses.

How Do Organic Vapors Swell Ultrathin Films of Polymer of Intrinsic Microporosity PIM-1?

Dynamic sorption of ethanol and toluene vapor into ultrathin supported films of polymer of intrinsic microporosity PIM-1 down to a thickness of 6 nm are studied with a combination of in situ spectroscopic ellipsometry and in situ X-ray reflectivity. Both ethanol and toluene significantly swell the PIM-1 matrix and, at the same time, induce persistent structural relaxations of the frozen-in glassy PIM-1 morphology. For ethanol below 20 nm, three effects were identified. First, the swelling magnitude at high vapor pressures is reduced by about 30% as compared to that of thicker films. Second, at low penetrant activities (below 0.3p/p₀), films below 20 nm are able to absorb slightly more penetrant as compared with thicker films despite a similar swelling magnitude. Third, for the ultrathin films, the onset of the dynamic penetrant-induced glass transition P_g has been found to shift to higher values, indicating higher resistance to plasticization. All of these effects are consistent with a view where immobilization of the superglassy PIM-1 at the substrate surface leads to an arrested, even more rigid, and plasticization-resistant, yet still very open, microporous structure. PIM-1 in contact with the larger and more condensable toluene shows very complex, heterogeneous swelling dynamics, and two distinct penetrant-induced relaxation phenomena, probably associated with the film outer surface and the bulk, are detected. Following the direction of the penetrant's diffusion, the surface seems to plasticize earlier than the bulk, and the two relaxations remain well separated down to 6 nm film thickness, where they remarkably merge to form just a single relaxation.

Dynamics of poly(vinyl methyl ketone) thin films studied by local dielectric spectroscopy

Local dielectric spectroscopy, which entails measuring the change in resonance frequency of the conducting tip of an atomic force microscope to determine the complex permittivity of a sample with high spatial (lateral) resolution, was employed to characterize the dynamics of thin films of poly(vinyl methyl ketone) (PVMK) having different substrate and top surface layers. A free surface yields the usual speeding up of the segmental dynamics, corresponding to a glass transition suppression of 6.5° for 18 nm film thickness. This result is unaffected by the presence of a glassy, compatible polymer, poly-4-vinyl phenol (PVPh), between the metal substrate and the PVMK. However, covering the top surface with a thin layer of the PVPh suppresses the dynamics. The speeding up of PVMK segmental motions observed for a free surface is absent due to interfacial interactions of the PVMK with the glass layer, an effect not seen when the top layer is an incompatible polymer.

Decoupling of surface diffusion and relaxation dynamics of molecular glasses

Tobacco mosaic virus is used as a probe to measure surface diffusion of ultrathin films of N,N'-Bis(3-methylphenyl)-N,N'-diphenylbenzidine (TPD) (12 nm [Formula: see text] 53 nm, where [Formula: see text] is the film thickness) at various temperatures below the glass transition temperature, [Formula: see text], of all films. As the film thickness is decreased, [Formula: see text] decreases rapidly and the average film dynamics are enhanced by 6-14 orders of magnitude. We show that the surface diffusion is invariant of the film thickness decrease and the resulting enhanced overall mobility. The values of the surface diffusion coefficient and its temperature dependence are invariant of film thickness and are the same as the corresponding bulk values ([Formula: see text]400 nm). For the thinnest films ([Formula: see text]20 nm), the effective activation energy for rearrangement (temperature dependence of relaxation times) becomes smaller than the activation energy for surface diffusion. These results suggest that the fast surface diffusion is decoupled from film relaxation dynamics and is a solely free surface property.

Invariant Fast Diffusion on the Surfaces of Ultrastable and Aged Molecular Glasses

Surface diffusion of molecular glasses is found to be orders of magnitude faster than bulk diffusion, with a stronger dependence on the molecular size and intermolecular interactions. In this study, we investigate the effect of variations in bulk dynamics on the surface diffusion of molecular glasses. Using the tobacco mosaic virus as a probe particle, we measure the surface diffusion on glasses of the same composition but with orders of magnitude of variations in bulk relaxation dynamics, produced by physical vapor deposition, physical aging, and liquid quenching. The bulk fictive temperatures of these glasses span over 35 K, indicating 13 to 20 orders of magnitude changes in bulk relaxation times. However, the surface diffusion coefficients on these glasses are measured to be identical at two temperatures below the bulk glass transition temperature T_{g}. These results suggest that surface diffusion has no dependence on the bulk relaxation dynamics when measured below T_{g}.

Why is surface diffusion the same in ultrastable, ordinary, aged, and ultrathin molecular glasses?

The primitive/JG relaxation explains the same surface diffusion coefficient in ordinary, ultrastable and thin film glasses of OTP and TPD.

Bound Layers "Cloak" Nanoparticles in Strongly Interacting Polymer Nanocomposites

Polymer-nanoparticle (NP) interfacial interactions are expected to strongly influence the properties of nanocomposites, but surprisingly, experiments often report small or no changes in the glass transition temperature, T_g. To understand this paradoxical situation, we simulate nanocomposites over a broad range of polymer-NP interaction strengths, ε. When ε is stronger than the polymer-polymer interaction, a distinct relaxation that is slower than the main α-relaxation emerges, arising from an adsorbed "bound" polymer layer near the NP surface. This bound layer "cloaks" the NPs, so that the dynamics of the matrix polymer are largely unaffected. Consequently, T_g defined from the temperature dependence of the routinely measured thermodynamics or the polymer matrix relaxation is nearly independent of ε, in accord with many experiments. Apparently, quasi-thermodynamic measurements do not reliably reflect dynamical changes in the bound layer, which alter the overall composite dynamics. These findings clarify the relation between quasi-thermodynamic T_g measurements and nanocomposite dynamics, and should also apply to thin polymer films.