Progress towards a phenomenological picture and theoretical understanding of glassy dynamics and vitrification near interfaces and under nanoconfinement

Thickness-Dependent Segmental Dynamics in Supported Thin Films: Insights from a Dynamically Correlated Network Model

A large body of experimental studies shows that the local dynamics in supercooled liquids are significantly altered by spatial nanoconfinement. In a previous study, we proposed a concept of a dynamically correlated network (DCN) model, which assumes that segments in a supercooled liquid undergo cooperative rearrangements within a network-like cluster. We further demonstrated that a model modified for freestanding thin films can predict for the glass transition dynamics in atactic polystyrene (PS) films consistent with experimental results. In this study, we adapted the model to apply it to supported thin films by introducing a layer of virtual vacant segments at the free surface and virtual anchoring segments at the liquid/substrate interface. The latter segments, carrying a finite number of virtual segments, reduce mobility at the interface. We evaluated the cooperative cluster size and distribution with respect to temperature and film thickness, along with the average relaxation time and glass transition temperature Tg for supported thin films of PS. The model predicted that the thickness dependence of Tg for PS becomes stronger with increasing time scale, and this result agreed well with experimental data across different timescales from pseudothermodynamic and dynamic measurements. The results provide insights into the origin of the dynamical decoupling between pseudothermodynamic and dynamic glass transition behaviors.

Quantitative Understanding of Liquid Dynamics at Interfaces from a Free-Energy Landscape Perspective

On the basis of molecular dynamics simulations of water and ethanol in nanopores, we devise a methodology to determine the free-energy landscape (FEL) imposed by an interface on an adjoining liquid directly from the particle trajectories. The methodology merely uses the statistical mechanical relation between occupancy and energy and, hence, is particularly suitable in complex situations, e.g., for disordered or rough atomistic interfaces and molecular liquids, as encountered in many biological, geological, and technological situations. Moreover, we show that the thus-obtained FEL enables a quantitative understanding of interface effects on liquid dynamics. Specifically, by determining the local minima and barriers of the FEL and using an Arrhenius-like relation, we reproduce the very strong spatial variation of the structural relaxation time of water and ethanol across nanopores over a broad temperature range. We anticipate that the proposed FEL approach is transferable to various other liquids and interfaces.

Visualization of Nanocrystallites in Organic Semiconducting Blends Using Cryo-Electron Microscopy

The efficiency of an organic solar cell is highly dependent on the complex, interpenetrating morphology, and molecular order within the composite phases of the bulk heterojunction (BHJ) blend. Both these microstructural aspects are strongly influenced by the processing conditions and chemical design of donor/acceptor materials. To establish improved structure-function relationships, it is vital to visualize the local microstructural order to provide specific local information about donor/acceptor interfaces and crystalline texture in BHJ blend films. The visualization of nanocrystallites, however, is difficult due to the complex semi-crystalline structure with few characterization techniques capable of visualizing the molecular ordering of soft materials at the nanoscale. Here, it is demonstrated how cryo-electron microscopy can be utilized to visualize local nanoscale order. This method is used to understand the distribution/orientation of crystallites in a BHJ blend. Long-range (>300 nm) texturing of IEICO-4F crystallites oriented in an edge-on fashion is observed, which has not previously been observed for spin-coated materials. This approach provides a wealth of quantitative information about the texture and size of nanocrystallites, which can be utilized to understand charge generation and transport in organic film. This study guides tailoring the material design and processing conditions for high-performance organic optoelectronic devices.

Intramolecular dynamic coupling slows surface relaxation of polymer glasses

Over the past three decades, studies have indicated a mobile surface layer with steep gradients on glass surfaces. Among various glasses, polymers are unique because intramolecular interactions - combined with chain connectivity - can alter surface dynamics, but their fundamental role has remained elusive. By devising polymer surfaces occupied by chain loops of various penetration depths, combined with surface dissipation experiments and Monte Carlo simulations, we demonstrate that the intramolecular dynamic coupling along surface chains causes the sluggish bulk polymers to suppress the fast surface dynamics. Such effect leads to that accelerated segmental relaxation on polymer glass surfaces markedly slows when the surface polymers extend chain loops deeper into the film interior. The surface mobility suppression due to the intramolecular coupling reduces the magnitude of the reduction in glass transition temperature commonly observed in thin films, enabling new opportunities for tailoring polymer properties at interfaces and under confinement and producing glasses with enhanced thermal stability.

Predicting compatibilized polymer blend toughness

Polymer blends can yield superior materials by merging the unique properties of their components. However, these mixtures often phase separate, leading to brittleness. While compatibilizers can toughen these blends, their vast design space makes optimization difficult. Here, we develop a model to predict the toughness of compatibilized glassy polymer mixtures. This theory reveals that compatibilizers increase blend toughness by creating molecular bridges that stitch the interface together. We validate this theory by directly comparing its predictions to extensive molecular dynamics simulations in which we vary polymer incompatibility, chain stiffness, compatibilizer areal density, and blockiness of copolymer compatibilizers. We then parameterize the model using self-consistent field theory and confirm its ability to make predictions for practical applications through comparison with simulations and experiments. These results suggest that the theory can optimize compatibilizer design for industrial glassy polymer blends in silico while providing microscopic insight, allowing for the development of next-generation mixtures.

Local dynamics and failure of inhomogeneous polymer networks

Inhomogeneous crosslinked polymers are powerful platforms for materials design, because they can be synthesized from materials that provide complimentary properties to the resulting gel. For example, a membrane with both glassy and rubbery domains will be mechanically robust while enabling transport. The dynamics, and mechanical and failure properties of rubbery/glassy conetworks are only beginning to be studied, and there is likely to be strong heterogeneities in the dynamics and mechanical response. In this study, we use coarse-grained molecular dynamics simulations to generate microphase separated rubbery/glassy polymer networks with a bicontinuous morphology via in silico crosslinking. We study the effect of phase boundary on the local mobility gradient, and our simulation results reveal an asymmetric shift in the local mobility gradient across the interface that extends deeper into the phase with a lower Tg when the system temperature is between the glass transition temperatures of the two phases. Moreover, by employing a model that allows bond breaking, we examine the microscopic mechanism for failure in these networks as a function of the molecular weight of polymer strands between crosslinks and the number fraction of the glassy domain. Under uniaxial extension, we find that the stress is initially larger in the glassy domain. As the deformation proceeds, the segmental dynamics of the two phases homogenize, and subsequently bond breaking begins.

A Molecular Dynamics Study of Mechanical and Conformational Properties of Conjugated Polymer Thin Films

Understanding and predicting the mechanical and conformational properties of conjugated polymer (CP) thin films are a central focus in flexible electronic device research. Employing molecular dynamics simulations with an architecture-transferable chemistry-specific coarse-grained (CG) model of poly(3-alkylthiophene)s (P3ATs), developed by using an energy renormalization approach, we investigate the mechanical and conformational behavior of P3AT thin films during deformation. The density profiles and measures of local mobility identify a softer interfacial layer for all films, the thickness of which does not depend on M w or side-chain length. Remarkably, Young's modulus measured via nanoindentation is more sensitive to M w than for tensile tests, which we attribute to distinct deformation mechanisms. High-M w thin films show increased toughness, whereas longer side-chain lengths of P3AT resulted in lower Young's modulus. Fractures in low-M w thin films occur through chain pullout due to insufficient chain entanglement and crazing in the plastic region. Importantly, stretching promoted both chain alignment and longer conjugation lengths of P3AT, potentially enhancing its electronic properties. For instance, at room temperature, stretching P3HT thin films to 150% increases the conjugated length of P3HT thin films from 2.7 nm to 4.7 nm, aligning with previous experimental findings and all-atom simulation results. Furthermore, high-M w thin films display elevated friction forces due to the chain accumulation on the indenter, with negligible variations in the friction coefficient across all thin film systems. These findings offer valuable insights that enhance our understanding and guide the rational design of CP thin films in flexible electronics.

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.

Dynamics of polylactic acid under ultrafine nanoconfinement: The collective interface effect and the spatial gradient

Polymers under nanoconfinement can exhibit large alterations in dynamics from their bulk values due to an interface effect. However, understanding the interface effect remains a challenge, especially in the ultrafine nanoconfinement region. In this work, we prepare new geometries with ultrafine nanoconfinement ∼10nm through controlled distributions of the crystalline phases and the amorphous phases of a model semi-crystalline polymer, i.e., the polylactic acid. The broadband dielectric spectroscopy measurements show that ultrafine nanoconfinement leads to a large elevation in the glass transition temperature and a strong increment in the polymer fragility index. Moreover, new relaxation time profile analyses demonstrate a spatial gradient that can be well described by either a single-exponential decay or a double-exponential decay functional form near the middle of the film with a collective interface effect. However, the dynamics at the 1-2 nm vicinity of the interface exhibit a power-law decay that is different from the single-exponential decay or double-exponential decay functional forms as predicted by theories. Thus, these results call for further investigations of the interface effect on polymer dynamics, especially for interfaces with perturbed chain packing.

Surface Glass Transition Temperature Region of Diarylethene Films Determined by Nano-Marangoni Effect

For the last two decades, research has addressed whether the glass transition temperature and the molecular motions on the surface of organic films are significantly different from those inside the bulk glasses. It is reported that the surface of the photochromic diarylethene film prepared by vacuum deposition has fluidity and the vacuum deposition of small amount of rubrene molecules induces surface tension fluctuations, generating dents due to the Marangoni flow in nanoscale. The depth of the dents increases in proportion to these radii for the colorless diarylethene film with a bulk glass transition temperature (Tg) close to room temperature. On the other hand, in the colored diarylethene obtained by UV irradiation to the colorless film, the depth becomes constant at a certain level. The Tg distribution in the depth direction is clarified based on an analysis of the dent depth. By approximating the obtained Tg depth distribution with an exponential function, the outermost surface Tg is about 100 K lower than the bulk Tg in the case of photoisomerized diarylethene.

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.

Eliminating the Tg-confinement and fragility-confinement effects in poly(4-methylstyrene) films by incorporation of 3 mol % 2-ethylheyxl acrylate comonomer

Nanoconfined poly(4-methylstyrene) [P(4-MS)] films exhibit reductions in glass transition temperature (Tg) relative to bulk Tg (Tg,bulk). Ellipsometry reveals that 15-nm-thick P(4-MS) films supported on silicon exhibit Tg - Tg,bulk = - 15 °C. P(4-MS) films also exhibit fragility-confinement effects; fragility decreases ∼60% in going from bulk to a 20-nm-thick film. Previous research found that incorporating 2-6 mol  % 2-ethylhexyl acrylate (EHA) comonomer in styrene-based random copolymers eliminates Tg- and fragility-confinement effects in polystyrene. Here, we demonstrate that incorporating 3 mol  % EHA in a 4-MS-based random copolymer, 97/3 P(4-MS/EHA), eliminates the Tg- and fragility-confinement effects. The invariance of fragility with nanoconfinement of 97/3 P(4-MS/EHA) films, hypothesized to originate from the interdigitation of ethylhexyl groups, indicates that the presence of EHA prevents the free surface from perturbing chain packing and the cooperative mobility associated with Tg. This method of eliminating confinement effects is advantageous as it relies on the simplest of polymerization methods and neat copolymer only slightly altered in composition from homopolymer. We also investigated whether we could eliminate the Tg-confinement effect with low levels of 2-ethylhexyl methacrylate (EHMA) in 4-MS-based or styrene-based copolymers. Although EHMA is structurally nearly identical to EHA, 4-MS-based and styrene-based copolymers incorporating 4 mol  % EHMA exhibit Tg-confinement effects similar to P(4-MS) and polystyrene. These results support the special character of EHA in eliminating confinement effects originating at free surfaces.

Development of broad modulus profile upon polymer-polymer interface formation between immiscible glassy-rubbery domains

Interfaces of glassy materials such as thin films, blends, and composites create strong unidirectional gradients to the local heterogeneous dynamics that can be used to elucidate the length scales and mechanisms associated with the dynamic heterogeneity of glasses. We focus on bilayer films of two different polymers with very different glass transition temperatures ([Formula: see text]) where previous work has demonstrated a long-range (∼200 nm) profile in local [Formula: see text] is established between immiscible glassy and rubbery polymer domains when the polymer-polymer interface is formed to equilibrium. Here, we demonstrate that an equally long-ranged gradient in local modulus [Formula: see text] is established when the polymer-polymer interface ([Formula: see text]5 nm) is formed between domains of glassy polystyrene (PS) and rubbery poly(butadiene) (PB), consistent with previous reports of a broad [Formula: see text] profile in this system. A continuum physics model for the shear wave propagation caused by a quartz crystal microbalance across a PB/PS bilayer film is used to measure the viscoelastic properties of the bilayer during the evolution of the PB/PS interface showing the development of a broad gradient in local modulus [Formula: see text] spanning [Formula: see text]180 nm between the glassy and rubbery domains of PS and PB. We suggest these broad profiles in [Formula: see text] and [Formula: see text] arise from a coupling of the spectrum of vibrational modes across the polymer-polymer interface as a result of acoustic impedance matching of sound waves with [Formula: see text] nm during interface broadening that can then trigger density fluctuations in the neighboring domain.

Absolute local conformation of poly(methyl methacrylate) chains adsorbed on a quartz surface

Polymer chains at a buried interface with an inorganic solid play a critical role in the performance of polymer nanocomposites and adhesives. Sum frequency generation (SFG) vibrational spectroscopy with a sub-nanometer depth resolution provides valuable information regarding the orientation angle of functional groups at interfaces. However, in the case of conventional SFG, since the signal intensity is proportional to the square of the second-order nonlinear optical susceptibility and thereby loses phase information, it cannot be unambiguously determined whether the functional groups face upward or downward. This problem can be solved by phase-sensitive SFG (ps-SFG). We here applied ps-SFG to poly(methyl methacrylate) (PMMA) chains in direct contact with a quartz surface, shedding light on the local conformation of chains adsorbed onto the solid surface. The measurements made it possible to determine the absolute orientation of the ester methyl groups of PMMA, which were oriented toward the quartz interface. Combining ps-SFG with all-atomistic molecular dynamics simulation, the distribution of the local conformation and the driving force are also discussed.

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.

Glass formation and dynamics of model polymer films with one versus two active interfaces

Polymers and other glass-forming liquids can exhibit profound alterations in dynamics in the nanoscale vicinity of interfaces, over a range appreciably exceeding that of typical interfacial thermodynamic gradients. The understanding of these dynamical gradients is particularly complicated in systems with internal or external nanoscale dimensions, where a gradient nucleated at one interface can impinge on a second, potentially distinct, interface. To better understand the interactions that govern system dynamics and glass formation in these cases, here we simulate the baseline case of a glass-forming polymer film, over a wide range of thickness, supported on a dynamically neutral substrate that has little effect on nearby dynamics. We compare these results to our prior simulations of freestanding films. Results indicate that dynamical gradients in our simulated systems, as measured based upon translational relaxation, are simply truncated when they impinge on a secondary surface that is locally dynamically neutral. Altered film behavior can be described almost entirely by gradient effects down to the thinnest films probed, with no evidence for finite-size effects sometimes posited to play a role in these systems. Finally, our simulations predict that linear gradient overlap effects in the presence of symmetric dynamically active interfaces yield a non-monotonic variation of the whole free standing film stretching exponent (relaxation time distribution breadth). The maximum relaxation time distribution breadth in simulation is found at a film thickness of 4-5 times the interfacial gradient range. Observation of this maximum in experiment would provide an important validation that the gradient behavior observed in simulation persists to experimental timescales. If validated, observation of this maximum would potentially also enable determination of the dynamic gradient range from experimental mean-film measurements of film dynamics.

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.

Soft pinning: Experimental validation of static correlations in supercooled molecular glass-forming liquids

Enormous enhancement in the viscosity of a liquid near its glass transition is a hallmark of glass transition. Within a class of theoretical frameworks, it is connected to growing many-body static correlations near the transition, often called "amorphous ordering." At the same time, some theories do not invoke the existence of such a static length scale in the problem. Thus, proving the existence and possible estimation of the static length scales of amorphous order in different glass-forming liquids is very important to validate or falsify the predictions of these theories and unravel the true physics of glass formation. Experiments on molecular glass-forming liquids become pivotal in this scenario as the viscosity grows several folds (∼ 10 14 ), and simulations or colloidal glass experiments fail to access these required long-time scales. Here we design an experiment to extract the static length scales in molecular liquids using dilute amounts of another large molecule as a pinning site. Results from dielectric relaxation experiments on supercooled Glycerol with different pinning concentrations of Sorbitol and Glucose, as well as the simulations on a few model glass-forming liquids with pinning sites, indicate the versatility of the proposed method, opening possible new avenues to study the physics of glass transition in other molecular liquids.

Surface premelting and melting of colloidal glasses

The nature of liquid-to-glass transition is a major puzzle in science. A similar challenge exists in glass-to-liquid transition, i.e., glass melting, especially for the poorly investigated surface effects. Here, we assemble colloidal glasses by vapor deposition and melt them by tuning particle attractions. The structural and dynamic parameters saturate at different depths, which define a surface liquid layer and an intermediate glassy layer. The power-law growth of both layers and melting front behaviors at different heating rates are similar to crystal premelting and melting, suggesting that premelting and melting can be generalized to amorphous solids. The measured single-particle kinetics reveal various features and confirm theoretical predictions for glass surface layer.

Influence of the Graft Length on Nanocomposite Structure and Interfacial Dynamics

Both the dispersion state of nanoparticles (NPs) within polymer nanocomposites (PNCs) and the dynamical state of the polymer altered by the presence of the NP/polymer interfaces have a strong impact on the macroscopic properties of PNCs. In particular, mechanical properties are strongly affected by percolation of hard phases, which may be NP networks, dynamically modified polymer regions, or combinations of both. In this article, the impact on dispersion and dynamics of surface modification of the NPs by short monomethoxysilanes with eight carbons in the alkyl part (C₈) is studied. As a function of grafting density and particle content, polymer dynamics is followed by broadband dielectric spectroscopy and analyzed by an interfacial layer model, whereas the particle dispersion is investigated by small-angle X-ray scattering and analyzed by reverse Monte Carlo simulations. NP dispersions are found to be destabilized only at the highest grafting. The interfacial layer formalism allows the clear identification of the volume fraction of interfacial polymer, with its characteristic time. The strongest dynamical slow-down in the polymer is found for unmodified NPs, while grafting weakens this effect progressively. The combination of all three techniques enables a unique measurement of the true thickness of the interfacial layer, which is ca. 5 nm. Finally, the comparison between longer (C₁₈) and shorter (C₈) grafts provides unprecedented insight into the efficacy and tunability of surface modification. It is shown that C₈-grafting allows for a more progressive tuning, which goes beyond a pure mass effect.

How Tuning Interfaces Impacts the Dynamics and Structure of Polymer Nanocomposites Simultaneously

Fundamental understanding of the macroscopic properties of polymer nanocomposites (PNCs) remains difficult due to the complex interplay of microscopic dynamics and structure, namely interfacial layer relaxations and three-dimensional nanoparticle (NP) arrangements. The effect of surface modification by alkyl methoxysilanes at different grafting densities has been studied in PNCs made of poly(2-vinylpyridine) and spherical 20 nm silica NPs. The segmental dynamics has been probed by broadband dielectric spectroscopy and the filler structure by small-angle X-ray scattering and reverse Monte Carlo simulations. By combining the particle configurations with the interfacial layer properties, it is shown how surface modification tunes the attractive polymer-particle interactions: bare NPs slow down the polymer interfacial layer dynamics over a thickness of ca. 5 nm, while grafting screens these interactions. Our analysis of interparticle spacings and segmental dynamics provides unprecedented insights into the effect of surface modification on the main characteristics of PNCs: particle interactions and polymer interfacial layers.

End-Tethered Chains Increase the Local Glass Transition Temperature of Matrix Chains by 45 K Next to Solid Substrates Independent of Chain Length

The local glass transition temperature T_g of pyrene-labeled polystyrene (PS) chains intermixed with end-tethered PS chains grafted to a neutral silica substrate was measured by fluorescence spectroscopy. To isolate the impact of the grafted chains, the films were capped with bulk neat PS layers eliminating competing effects of the free surface. Results demonstrate that end-grafted chains strongly increase the local T_g of matrix chains by ≈45 K relative to bulk T_g, independent of grafted chain molecular weight from M_n = 8.6 to 212 kg/mol and chemical end-group, over a wide range of grafting densities σ = 0.003 to 0.33 chains/nm² spanning the mushroom-to-brush transition regime. The tens-of-degree increase in local T_g resulting from immobilization of the chain ends by covalent bonding in this athermal system suggests a mechanism that substantially increases the local activation energy required for cooperative rearrangements.

Relaxation dynamics of deformed polymer nanocomposites as revealed by small-angle scattering and rheology

The relaxation dynamics of polystyrene (PS)/silica nanocomposites after a large step deformation are studied by a combination of small-angle scattering techniques and rheology. Small-angle X-ray scattering measurements and rheology show clear signatures of nanoparticle aggregation that enhances the mechanical properties of the polymer nanocomposites (PNCs) in the linear viscoelastic regime and during the initial phase of stress relaxation along with accelerated relaxation dynamics. Small-angle neutron scattering experiments under the zero-average-contrast condition reveal, however, smaller structural anisotropy in the PNCs than that in the neat polymer matrix, as well as accelerated anisotropy relaxation. In addition, the degrees of anisotropy reduction and relaxation dynamics acceleration increase with increasing nanoparticle loading. These results are in sharp contrast to the prevailing viewpoint of enhanced molecular deformation as the main mechanism for the mechanical enhancement in PNCs. Furthermore, the observed acceleration of stress relaxation and reduction in structural anisotropy point to two types of nonlinear effects in the relaxation dynamics of PNCs at large deformation.

Nonmonotonic Dynamical Correlations beneath the Surface of Glass-Forming Liquids

Collective motion over increasing length scales is a signature of the vitrification process of liquids. We demonstrate how distinct static and dynamic length scales govern the dynamics of vitrifying films. In contrast to a monotonically growing static correlation length, the dynamical correlation length that measures the extent of surface-dynamics acceleration into the bulk displays a striking nonmonotonic temperature evolution that is robust also against changes in detailed interatomic interaction. This nonmonotonic change defines a crossover temperature T_{*} that is distinct from the critical temperature T_{c} of mode-coupling theory. We connect this nonmonotonic change to a morphological change of cooperative rearrangement regions of fast particles, and to the point where the decoupling of fast-particle motion from the bulk relaxation is most sensitive to fluctuations. We propose a rigorous definition of this new crossover temperature T_{*} within a recent extension of mode-coupling theory, the stochastic β-relaxation theory.

Microscale mobile surface double layer in a glassy polymer

This study examines the origin of the widely different length scales, h_t-nanometers to micrometers-that have been observed for the propagation of the near-surface enhanced mobility in glassy polymers. Mechanical relaxations of polystyrene films with thicknesses, h, from 5 nm to 186 μm have been studied. For h < ~1 μm, the films relaxed faster than the bulk and the relaxation time decreased with decreasing h below ~100 nm, consistent with the enhanced dynamics originating from a near-surface nanolayer. For h > ~1 μm, a bulk-like relaxation mode emerged, while the fast mode changed to one that extended over ~1 μm from the free surface. These findings evidence that the mobile surface region is inhomogeneous, comprising a nanoscale outer layer and a slower microscale sublayer that relax by different mechanisms. Consequently, measurements probing the enhanced mobility of different mechanisms may find vastly different h_t's as shown by the literature.

Decoupling of Glassy Dynamics from Viscosity in Thin Supported Poly(n-butyl methacrylate) Films

We utilized fast scanning calorimetry to characterize the glass transition temperature (T g) and intrinsic molecular mobility of low-molecular-weight poly(n-butyl methacrylate) thin films of varying thicknesses. We found that the T g and intrinsic molecular mobility were coupled, showing no film thickness-dependent variation. We further employed a unique noncontact capillary nanoshearing technique to directly probe layer-resolved gradients in the rheological response of these films. We found that layer-resolved shear mobility was enhanced with a reduction in film thickness, whereas the effective viscosity decreased. Our results highlight the importance of polymer-substrate attractive interactions and free surface-promoted enhanced mobility, establishing a competitive nanoconfinement effect in poly(n-butyl methacrylate) thin films. Moreover, the findings indicate a decoupling in the thickness-dependent variation of T g and intrinsic molecular mobility with the mechanical responses (shear mobility and effective viscosity).

Screening and collective effects in randomly pinned fluids: a new theoretical framework

We propose a theoretical framework for the dynamics of bulk isotropic hard-sphere systems in the presence of randomly pinned particles and apply this theory to supercooled water to validate it. Structural relaxation is mainly governed by local and non-local activated process. As the pinned fraction grows, a local caging constraint becomes stronger and the long range collective aspect of relaxation is screened by immobile obstacles. Different responses of the local and cooperative motions results in subtle predictions for how the alpha relaxation time varies with pinning and density. Our theoretical analysis for the relaxation time of water with pinned molecules quantitatively well describe previous simulations. In addition, the thermal dependence of relaxation for unpinned bulk water is also consistent with prior computational and experimental data.

Orientational wetting and dynamical correlations toward glass transition on the surface of imidazolium-based ionic liquids

Surface induces many fascinating physical phenomena, such as dynamic acceleration, surface anchoring, and orientational wetting, and, thus, is of great interest to study. Here, we report classic molecular dynamics simulations on the free-standing surface of imidazolium-based ionic liquids (ILs) [C4mim][PF6] and [C10mim][PF6]. On [C10mim][PF6] surface, a significant orientational wetting is observed, with the wetting strength showing a diverging tendency. Depth of the wetting was captured from the density and orientational order profile by a static length, which remarkably increases below the temperature Tstat upon cooling down. The dynamical correlation length that measures the distance of surface-dynamics acceleration into the bulk was characterized via the spatial-dependent mobility. The translational correlation exhibits a similar drastic increment at Tstat, while the rotational correlation drastically increases at a lower temperature Trot. We connect these results to the dynamics in bulk liquids, by finding Tstat and Trot that correspond to the onset temperatures where the liquids become cooperative for translational and rotational relaxation, respectively. This signifies the importance of collective dynamics in the bulk on the orientational wetting and surface dynamics in the ILs.

Surface Mobility of Amorphous Indomethacin Containing Moisture and a Surfactant: A Concentration-Temperature Superposition Principle

An amorphous material can have vastly higher mobility on the surface than in the bulk and, as a result, shows fast surface crystallization. Most amorphous materials contain multiple components, but the effect of composition on surface dynamics remains poorly understood. In this study, the surface mobility of amorphous indomethacin was measured using the method of surface-grating decay in the presence of moisture and the surfactant Tween 20. It is found that both components significantly enhance the surface mobility, and their effects are well described by the principle of concentration-temperature superposition (CTS); that is, the same surface dynamics is observed at the same T_g-normalized temperature T/T_g, where T_g is the composition-dependent glass transition temperature. For doped indomethacin showing CTS, the mechanism of surface evolution for a 1000 nm wavelength surface grating transitions from viscous flow at high temperatures to surface diffusion at low temperatures at 1.04 T_g. For the surfactant-doped system, the T_g used is the value for the surface layer that reflects the surface enrichment of the surfactant (measured by X-ray photoelectron spectroscopy). At a high surfactant concentration (>10% by weight), the surface-grating decay rate in the surface-diffusion regime is limited by the large, slow-diffusing surfactant molecules; in this case, CTS holds only for the viscous-flow regime. The CTS principle allows the prediction of the surface dynamics of multicomponent amorphous materials.

Role of torsional potential in chain conformation, thermodynamics, and glass formation of simulated polybutadiene melts

For polymer chains, the torsional potential is an important intramolecular energy influencing chain flexibility and segmental dynamics. Through molecular dynamics simulations of an atomistic model for melts of cis-trans-1,4-polybutadiene (PBD), we explore the effect of the torsions on conformational properties (bond vector correlations and mean-square internal distances), fundamental thermodynamic quantities (density, compressibility, internal energy, and specific heat), and glass transition temperature T_g. This is achieved by systematically reducing the strength of the torsional potential, starting from the chemically realistic chain (CRC) model with the full potential toward the freely rotating chain (FRC) model without the torsional potential. For the equilibrium liquid, we find that the effect of the torsions on polymer conformations is very weak. Still weaker is the influence on the monomer density ρ and isothermal compressibility κ_T of the polymer liquid, both of which can be considered as independent of the torsional potential. We show that a van der Waals-like model proposed by Long and Lequeux [Eur. Phys. J. E 4, 371 (2001)] allows us to describe very well the temperature (T) dependence of ρ and κ_T. We also find that our data obey the linear relation between 1/k_BTρκ_T and 1/T (with the Boltzmann constant k_B) that has recently been predicted and verified on the experiment by Mirigian and Schweizer [J. Chem. Phys. 140, 194507 (2014)]. For the equilibrium liquid, simulations result in a specific heat, at constant pressure and at constant volume, which increases on cooling. This T dependence is opposite to the one found experimentally for many polymer liquids, including PBD. We suggest that this difference between simulation and experiment may be attributed to quantum effects due to hydrogen atoms and backbone vibrations, which, by construction, are not included in the classical united-atom model employed here. Finally, we also determine T_g from the density-temperature curve monitored in a finite-rate cooling process. While the influence of the torsional potential on ρ(T) is vanishingly small in the equilibrium liquid, the effect of the torsions on T_g is large. We find that T_g decreases by about 150 K when going from the CRC to the FRC model.

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.

A β-NMR study of the depth, temperature, and molecular-weight dependence of secondary dynamics in polystyrene: Entropy–enthalpy compensation and dynamic gradients near the free surface

We investigated the depth, temperature, and molecular-weight (MW) dependence of the γ-relaxation in polystyrene glasses using implanted 8Li+ and β-detected nuclear magnetic resonance. Measurements were performed on thin films with MW ranging from 1.1 to 641 kg/mol. The temperature dependence of the average 8Li spin–lattice relaxation time [Formula: see text] was measured near the free surface and in the bulk. Spin–lattice relaxation is caused by phenyl ring flips, which involve transitions between local minima over free-energy barriers with enthalpic and entropic contributions. We used transition state theory to model the temperature dependence of the γ-relaxation, and hence [Formula: see text]. There is no clear correlation of the average entropy of activation [Formula: see text] and enthalpy of activation [Formula: see text] with MW, but there is a clear correlation between [Formula: see text] and [Formula: see text], i.e., entropy–enthalpy compensation. This results in the average Gibbs energy of activation, [Formula: see text], being approximately independent of MW. Measurements of the temperature dependence of [Formula: see text] as a function of depth below the free surface indicate the inherent entropic barrier, i.e., the entropy of activation corresponding to [Formula: see text] = 0, has an exponential dependence on the distance from the free surface before reaching the bulk value. This results in [Formula: see text] near the free surface being lower than the bulk. Combining these observations results in a model where the average fluctuation rate of the γ-relaxation has a “double-exponential” depth dependence. This model can explain the depth dependence of [Formula: see text] in polystyrene films. The characteristic length of enhanced dynamics is ∼6 nm and approximately independent of MW near room temperature.

Theory for the Elementary Time Scale of Stress Relaxation in Polymer Nanocomposites

We construct a microscopic theory for the elementary time scale of stress relaxation in dense polymer nanocomposites. The key dynamical event is proposed to involve the rearrangement of cohesive segment-nanoparticle (NP) tight bridging complexes via an activated small NP dilational motion, which allows the confined segments to relax. The corresponding activation energy is determined by the NP bridge coordination number and potential of mean force barrier. The activation energy varies nonlinearly with interfacial cohesion strength and NP concentration, and a universal master curve is predicted. The theory is in very good agreement with experiments. The underlying ideas are relevant to a variety of other hybrid macromolecular materials involving hard particles and soft macromolecules.

Characterization of the Interfacial Orientation and Molecular Conformation in a Glass-Forming Organic Semiconductor

The ability to control structure in molecular glasses has enabled them to play a key role in modern technology; in particular, they are ubiquitous in organic light-emitting diodes. While the interplay between bulk structure and optoelectronic properties has been extensively investigated, few studies have examined molecular orientation near buried interfaces despite its critical role in emergent functionality. Direct, quantitative measurements of buried molecular orientation are inherently challenging, and many methods are insensitive to orientation in amorphous soft matter or lack the necessary spatial resolution. To overcome these challenges, we use polarized resonant soft X-ray reflectivity (p-RSoXR) to measure nanometer-resolved, molecular orientation depth profiles of vapor-deposited thin films of an organic semiconductor Tris(4-carbazoyl-9-ylphenyl)amine (TCTA). Our depth profiling approach characterizes the vertical distribution of molecular orientation and reveals that molecules near the inorganic substrate and free surface have a different, nearly isotropic orientation compared to those of the anisotropic bulk. Comparison of p-RSoXR results with near-edge X-ray absorption fine structure spectroscopy and optical spectroscopies reveals that TCTA molecules away from the interfaces are predominantly planar, which may contribute to their attractive charge transport qualities. Buried interfaces are further investigated in a TCTA bilayer (each layer deposited under separate conditions resulting in different orientations) in which we find a narrow interface between orientationally distinct layers extending across ≈1 nm. Coupling this result with molecular dynamics simulations provides additional insight into the formation of interfacial structure. This study characterizes the local molecular orientation at various types of buried interfaces in vapor-deposited glasses and provides a foundation for future studies to develop critical structure-function relationships.

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.

The dynamics of freestanding films: predictions for poly(2-chlorostyrene) based on bulk pressure dependence and thoughtful sample averaging

In this paper we model the segmental relaxation in poly(2-chlorostyrene) 18 nm freestanding films, using only data on bulk samples to characterize the system, and predict film relaxation times (τ) as a function of temperature that are in semi-quantitative agreement with film data. The ability to translate bulk characterization into film predictions is a direct result of our previous work connecting the effects of free surfaces in films with those of changing pressure in the bulk. Our approach combines the Locally Correlated Lattice (LCL) equation of state for prediction of free volume values (V_free) at any given density (ρ), which are then used in the Cooperative Free Volume (CFV) rate model to predict τ(T, V_free). A key feature of this work is that we calculate the locally averaged density profile as a function of distance from the surface, ρ_av(z), using the CFV-predicted lengthscale, L_coop(z), over which rearranging molecular segments cooperate. As we have shown in the past, ρ_av(z) is significantly broader than the localized profile, ρ(z), which translates into a relaxation profile, τ(z), exhibiting a breadth that mirrors experimental and simulated results. In addition, we discuss the importance of averaging the log of position dependent relaxation times across a film sample (〈log τ(z)〉), as opposed to averaging the relaxation times, themselves, in order to best approximate a whole sample-averaged value that can be directly compared to experiment.

Theory of Transient Localization, Activated Dynamics, and a Macromolecular Glass Transition in Ring Polymer Liquids

We construct a segmental scale force level theory for the center-of-mass diffusion constant and corresponding relaxation time for globally compact unconcatenated ring polymer solutions and melts (degree of polymerization N). The approach is based on slowly decaying macromolecular scale intermolecular force dynamic correlations as the origin of their unusual dynamics. Unentangled Rouse, weakly caged, and activated regimes are predicted. The barrier of the activated regime scales linearly with N and as a power law of concentration, which drives a kinetic glass transition on the radius-of-gyration scale. The values of N at the two dynamic crossovers (Rouse to weakly caged, weakly caged to activated) are proportional, with nonuniversality entering mainly via macromolecular volume fraction and dimensionless compressibility. Quantitative comparisons with simulation data reveal good agreement. Aspects of intermediate time dynamics are analyzed, and predictions are made for the conditions required to observe a macromolecular glass transition in the laboratory and on the computer.

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.

Glass transition of ion-containing polymer melts in bulk and thin films

Ion-containing polymers often are good glass formers, and the glass transition temperature is an important parameter to consider for practical applications, which prescribes the working temperature range for different mechanical and dynamic properties. In this work, we present a systematic molecular dynamics simulation study on the coupling of ionic correlations with the glass transition, based on a generic coarse-grained model of ionic polymers. The variation of the glass transition temperature is examined concerning the influence of the electrostatic interaction strength, charge fraction, and charge sequence. The interplay with the film thickness effect is also discussed. Our results reveal a few typical features about the glass transition process that are in qualitative agreement with previous studies, further highlighting the effects of counterion entropy at weak ionic correlations and physical crosslinking of ionic aggregates at strong ionic correlations. Detailed parametric dependencies are displayed, which demonstrate that introducing strong ionic correlations promotes vitrification while adopting a precise charge sequence and applying strong confinement with weak surface affinity reduce the glass transition temperature. Overall, our investigation provides an improved picture towards a comprehensive understanding of the glass transition in ion-containing polymeric systems from a molecular simulation perspective.

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.

Mobility gradients yield rubbery surfaces on top of polymer glasses

Many emerging materials, such as ultrastable glasses^1,2 of interest for phone displays and OLED television screens, owe their properties to a gradient of enhanced mobility at the surface of glass-forming liquids. The discovery of this surface mobility enhancement^3-5 has reshaped our understanding of the behaviour of glass formers and of how to fashion them into improved materials. In polymeric glasses, these interfacial modifications are complicated by the existence of a second length scale-the size of the polymer chain-as well as the length scale of the interfacial mobility gradient^6-9. Here we present simulations, theory and time-resolved surface nano-creep experiments to reveal that this two-scale nature of glassy polymer surfaces drives the emergence of a transient rubbery, entangled-like surface behaviour even in polymers comprised of short, subentangled chains. We find that this effect emerges from superposed gradients in segmental dynamics and chain conformational statistics. The lifetime of this rubbery behaviour, which will have broad implications in constraining surface relaxations central to applications including tribology, adhesion, and surface healing of polymeric glasses, extends as the material is cooled. The surface layers suffer a general breakdown in time-temperature superposition (TTS), a fundamental tenet of polymer physics and rheology. This finding may require a reevaluation of strategies for the prediction of long-time properties in polymeric glasses with high interfacial areas. We expect that this interfacial transient elastomer effect and TTS breakdown should normally occur in macromolecular systems ranging from nanocomposites to thin films, where interfaces dominate material properties^5,10.