Theory of Spatial Gradients of Relaxation, Vitrification Temperature and Fragility of Glass-Forming Polymer Liquids Near Solid Substrates

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.

Modulus alteration of thin polystyrene films by their neighboring PDMS: Soft and hard confinement

It is highly demanded to understand the confinement effect on nanoconfined polymers. Recent studies reported a strong perturbation of local dynamics and substantial alteration of glass transition temperature Tg at nanoscale. However, how confinement affects the mechanical properties of polymers is not fully understood. Here, we show that the modulus of thin polymer films could be remarkedly altered through a polymer-polymer interface. The modulus of a thin polystyrene (PS) film next to a polydimethylsiloxane (PDMS) was determined from the PS-PDMS bilayer bulging test. A series of experiments show that the modulus of PS can be increased up to 37%, when the modulus of the neighboring PDMS varies from 1.04 to 4.88 MPa. The results demonstrate a strong sensitivity of mechanical properties of thin polymers to the hard/soft environment, which we attribute to the change of high-mobility layer by the polymer-polymer interface.

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.

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.

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.

Confinement Effects on the Spatially Inhomogeneous Dynamics in Metallic Glass Films

This work develops the elastically collective nonlinear Langevin equation theory to investigate, for the first time, the glassy dynamics in capped metallic glass thin films. Finite-size effects on the spatial gradient of structural relaxation time and glass transition temperature (T_g) are calculated at different temperatures and vitrification criteria. Molecular dynamics is significantly slowed down near rough solid surfaces, and the dynamics at location far from the interfaces is sped up. In thick films, the mobility gradient normalized by the bulk value obeys the double-exponential form since interference effects between two surfaces are weak. Reducing the film thickness induces a strong dynamic coupling between two surfaces and flattens the relaxation gradient. The normalized gradient of the glass transition temperature is independent of vitrification time scale criterion and can be fitted by a superposition function as the films are not ultrathin. The local fragility is found to remain unchanged with location. This finding suggests that one can use Angell plots of bulk relaxation time and the T_g spatial gradient to characterize glassy dynamics in metallic glass films. Our computational results agree well with experimental data and simulation.

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.

Interfacial Dynamics in Supported Ultrathin Polymer Films-From the Solid to the Free Interface

Molecular dynamics in ultrathin layers is investigated using nanostructured electrodes to perform broadband dielectric spectroscopy measurements, and by atomistic molecular dynamics simulations. Using poly(vinyl acetate) as the model system and taking advantage of access to the distribution of relaxation times in an extended temperature range above the glass transition temperature, T_g, we demonstrate that while the mean rates of the segmental relaxation remain bulklike down to 12 nm film thickness, modified molecular mobilities arise in the interfacial zones. Combining results from simulations and experiments, we show unambiguously that both the slow relaxations arising from adsorbed polymer segments and the faster modes attributed to segments in the vicinity of the free interface have non-Arrhenius temperature activation. These interfacial regions span thicknesses of ∼1.5 nm each just above the calorimetric T_g independent of molecular weight and film thickness. These deviations at interfaces are relevant for applications of polymers in adhesion, coatings, and polymer nanocomposites.

Direct observation of morphological transition for an adsorbed single polymer chain

A better understanding of the structure of polymers at solid interfaces is crucial for designing various polymer nano-composite materials from structural materials to nanomaterials for use in industry. To this end, the first step is to obtain information on how synthetic polymer chains adsorb onto a solid surface. We closely followed the trajectory of a single polymer chain on the surface as a function of temperature using atomic force microscopy. Combining the results with a full-atomistic molecular dynamics simulation revealed that the chain became more rigid on the way to reaching a pseudo-equilibrium state, accompanied by a change in its local conformation from mainly loops to trains. This information will be useful for regulating the physical properties of polymers at the interface.

Local Glass Transition Temperature Tg(z) Within Polystyrene Is Strongly Impacted by the Modulus of the Neighboring PDMS Domain

Profiles in the local glass transition temperature T_g(z) within polystyrene (PS) next to polydimethylsiloxane (PDMS) domains were determined using a localized fluorescence method. By changing the base to cross-linker ratio, we varied the cross-link density and, hence, the Young's modulus of PDMS (Sylgard 184). The local T_g(z) in PS at a distance of z = 50 nm away from the PS/PDMS interface was found to shift by 40 K as the PDMS modulus was varied from 0.9 to 2.6 MPa, demonstrating a strong sensitivity of this phenomenon to the rigidity of the neighboring domain. The extent the T_g(z) perturbation persists away from the PS/PDMS interface, z ≈ 65-90 nm before bulk T_g is recovered, is much shorter for this strongly immiscible system compared with the weakly immiscible systems studied previously, which we attribute to a smaller interfacial width, as the χ parameter for PS/PDMS is an order of magnitude larger.

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

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