Influence of Backbone Rigidity on Nanoscale Confinement Effects in Model Glass-Forming Polymers

Mechanism of Density Evolution of Polystyrene Adsorbed Layers on the Substrate

The density evolution of polystyrene (PS) adsorbed layers on phenyl-modified SiO2-Si substrates was investigated. The thickness and density of flattened layer on substrates with above 75% phenyl content increased over annealing time and could approach 4.7 nm and 1.37 g/cm3 at equilibrium, respectively, which were much higher than those on SiO2-Si. The annealing time for flattened chains to reach equilibrium increased with an increasing phenyl content on the substrate. The interface sensitive sum frequency generation vibrational spectroscopy (SFG) technique revealed that both the amount and the strength of the interfacial π-π interaction between the phenyl groups of substrates and in PS chains increased with annealing time. This resulted in more stretched chains perpendicularly, leading to a denser and thicker adsorbed layer with a closest-packing structure, driven by favorable enthalpy processes. Our work provides important insight into the densification mechanism of adsorbed flattened layers.

Role of additive size in the segmental dynamics and mechanical properties of cross-linked polymers

Thermoset materials often involve the addition of molecular and nanoparticle additives to alter various chemo-physical properties of importance in their ultimate applications. The resulting compositional heterogeneities can lead to either enhancement or degradation of thermoset properties, depending on the additive chemical structure and concentration. We tentatively explore this complex physical phenomenon through the consideration of a model polymeric additive to our coarse-grained (CG) thermoset investigated in previous works by simply varying the size of additive segments compared to those of polymer melt. We find that the additive modified thermoset material becomes chemically heterogeneous from additive aggregation when the additive segments become much smaller than those of the thermoset molecules, and a clear evidence is observed in the spatial distribution of local molecular stiffness estimated from Debye-Waller factor 〈u2〉. Despite the non-monotonic variation trends observed in dynamical and mechanical properties with decreasing additive segmental size, both the structural relaxation time and moduli (i.e., shear modulus and bulk modulus) exhibit scaling laws with 〈u2〉. The present work highlights the complex role of additive size played in the dynamical and mechanical properties of thermoset polymers, which should provide a better understanding for the glass formation process of cross-linked polymer composites.

Unveiling the Nanoconfinement Effect on Crystallization of Semicrystalline Polymers Using Coarse-Grained Molecular Dynamics Simulations

Semicrystalline polymers under nanoconfinement show distinct structural and thermomechanical properties compared to their bulk counterparts. Despite extensive research on semicrystalline polymers under nanoconfinement, the nanoconfinement effect on the local crystallization process and the unique structural evolution of such polymers have not been fully understood. In this study, we unveil such effects by using coarse-grained molecular dynamics simulations to study the crystallization process of a model semicrystalline polymer-polyvinyl alcohol (PVA)-under different levels of nanoconfinement induced by nanoparticles that are represented implicitly. We quantify in detail the evolution of the degree of crystallinity (XC) of PVA and examine distinct crystalline regions from simulation results. The results show that nanoconfinement can promote the crystallization process, especially at the early stage, and the interfaces between nanoparticles and polymer can function as crystallite nucleation sites. In general, the final XC of PVA increases with the levels of nanoconfinement. Further, nanoconfined cases show region-dependent XC with higher and earlier increase of XC in regions closer to the interfaces. By tracking region-dependent XC evolution, our results indicate that nanoconfinement can lead to a heterogenous crystallization process with a second-stage crystallite nucleation in regions further away from the interfaces. In addition, our results show that even under very high cooling rates, the nanoconfinement still promotes the crystallization of PVA. This study provides important insights into the underlying mechanisms for the intricate interplay between nanoconfinement and the crystallization behaviors of semicrystalline polymer, with the potential to guide the design and characterization of semicrystalline polymer-based nanocomposites.

Effect of penetrant-polymer interactions and shape on the motion of molecular penetrants in dense polymer networks

The diffusion of dilute molecular penetrants within polymers plays a crucial role in the advancement of material engineering for applications such as coatings and membrane separations. The potential of highly cross-linked polymer networks in these applications stems from their capacity to adjust the size and shape selectivity through subtle changes in network structures. In this paper, we use molecular dynamics simulation to understand the role of penetrant shape (aspect ratios) and its interaction with polymer networks on its diffusivity. We characterize both local penetrant hopping and the long-time diffusive motion for penetrants and consider different aspect ratios and penetrant-network interaction strengths at a variety of cross-link densities and temperatures. The shape affects the coupling of penetrant motion to the cross-link density- and temperature-dependent structural relaxation of networks and also affects the way a penetrant experiences the confinement from the network meshes. The attractive interaction between the penetrant and network primarily affects the former since only the system of dilute limit is of present interest. These results offer fundamental insights into the intricate interplay between penetrant characteristics and polymer network properties and also suggest future directions for manipulating polymer design to enhance the separation efficiency.

Fundamentals of Crystalline Evolution and Properties of Carbon Nanotube-Reinforced Polyether Ether Ketone Nanocomposites in Fused Filament Fabrication

As fused filament fabrication (FFF) continues to gain popularity, many studies are turning to nanomaterials or optimization of printing parameters to improve the materials' properties; however, many overlook how materials formulation and additive manufacturing (AM) processes cooperatively engineer the evolution of properties across length scales. Evaluating the in-process evolution of the nanocomposite using AM will provide a fundamental understanding of the material's microstructure, which can be tailored to create unique characteristics in functionality and performance. In this study, the crystallinity behavior of polyetheretherketone (PEEK) was studied in the presence of carbon nanotubes (CNTs) as a nucleation aid for improved crystallization during FFF processing. Using various characterization techniques and molecular dynamics simulations, it was discovered that the crystallization behavior of extruded filaments is very different from that of 3D printed roads. Additionally, the printed material exhibited cold crystallization, and the CNT addition increased the crystallization of printed roads, which were amorphous without CNT addition. Tensile strength and modulus were increased by as much as 42 and 51%, respectively, due to higher crystallinity during printing. Detailed knowledge on the morphology of PEEK-CNT used in FFF allows gaining a fundamental understanding of the morphological evolution occurring during the AM process that in turn enables formulating materials for the AM process to achieve tailored mechanical and functional properties, such as crystallinity or conductivity.

Effect of the polar group content on the glass transition temperature of ROMP copolymers

Polar groups have long been recognized to greatly influence the glass transition temperature (T_g) of polymers, but understanding the underlying physical mechanism remains a challenge. Here, we study the glass formation of ring-opening metathesis polymerization (ROMP) copolymers containing polar groups by employing all-atom molecular dynamics simulations. We show that although the number of hydrogen bonds (N_HB) and the cohesive energy density increase linearly as the content of polar groups (f_pol) increases, the T_g of ROMP copolymers increases with the increase of f_pol in a nonlinear fashion, and tends to plateau for sufficiently high f_pol. Importantly, we find that the increase rate of Gibbs free energy for HB breaking gradually slows down with the increase of f_pol, indicating that the HB is gradually stabilized. Therefore, T_g is jointly determined by N_HB and the strength of HBs in the system, while the latter dominates. Although N_HB increases linearly with increasing f_pol, the HB strength increases slowly with increasing f_pol, which leads to a decreasing rate of increase in T_g.

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.

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.

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

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

Ensemble fluctuations matter for variances of macroscopic variables

Extending recent work on stress fluctuations in complex fluids and amorphous solids we describe in general terms the ensemble average [Formula: see text] and the standard deviation [Formula: see text] of the variance [Formula: see text] of time series [Formula: see text] of a stochastic process x(t) measured over a finite sampling time [Formula: see text]. Assuming a stationary, Gaussian and ergodic process, [Formula: see text] is given by a functional [Formula: see text] of the autocorrelation function h(t). [Formula: see text] is shown to become large and similar to [Formula: see text] if [Formula: see text] corresponds to a fast relaxation process. Albeit [Formula: see text] does not hold in general for non-ergodic systems, the deviations for common systems with many microstates are merely finite-size corrections. Various issues are illustrated for shear-stress fluctuations in simple coarse-grained model systems.

The influence of additives on polymer matrix mobility and the glass transition

In the region near an interface, the microscopic properties of a glass forming liquid may be perturbed from their equilibrium bulk values. In this work, we probe how the interfacial effects of additive particles dispersed in a matrix can influence the local mobility of the material and its glass transition temperature, Tg. Experimental measurements and simulation results indicate that additives, such as nanoparticles, gas molecules, and oligomers, can shift the mobility and Tg of a surrounding polymer matrix (even for relatively small concentrations of additive; e.g., 5-10% by volume) relative to the pure bulk matrix, thus leading to Tg enhancement or suppression. Additives thus provide a potential route for modifying the properties of a polymer material without significantly changing its chemical composition. Here we apply the Limited Mobility (LM) model to simulate a matrix containing additive species. We show that both additive concentration, as well as the strength of its very local influence on the surrounding matrix material, will determine whether the Tg of the system is raised or lowered, relative to the pure matrix. We demonstrate that incorporation of additives into the simple LM simulation method, which has successfully described the behavior of bulk and thin film glassy solids, leads to direct connections with available experimental and simulation results for a broad range of polymer/additive systems.

Confinement effects on the mechanical heterogeneity of polymer fiber glasses

Both polymer fiber glasses and bulk polymer glasses exhibit nonlinear mechanical responses under uniaxial deformation. In polymer fibers, however, polymer chains are confined strongly and the surface area is relatively large compared to their volume. The confinement and the surface may lead to the spatially heterogeneous relaxation of chains in polymer fibers. In this work we perform molecular dynamics simulations and investigate the relation between the heterogeneous dynamics and the nonlinear mechanical responses at a molecular level. Our molecular simulations capture successfully not only the nonlinear mechanical response but also the dependence of mechanical properties on the strain rate of typical polymer glasses as in experiments. We find that the local elastic modulus and the nonaffine displacement are spatially heterogeneous in the pre-yield regime, which results in a lower elastic modulus for polymer fibers than bulk polymer glasses. In the post-yield regime, those mechanical properties become relatively homogeneous. Monomers with large nonaffine displacement are localized mainly at the interfacial region in the pre-yield regime while highly nonaffine monomers are distributed throughout the fibers in the post-yield regime. We show that the nonaffine displacement during deformation relates closely to the mechanical response of the polymer fibers. We also find that in the strain-hardening regime there is a significant difference in the energetic contribution to the stress between polymer fibers and bulk polymers, for which the modulus of the strain-hardening regime of the polymer fibers is smaller than that of bulk polymers.

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

We employ a new force-level statistical mechanical theory to predict spatial gradients of the structural relaxation time and T_g of polymer liquids near microscopically rough and smooth hard surfaces and contrast the results with vapor interface systems. Repulsive rough (smooth) surfaces induce large slowing down (modest speeding up) of the relaxation time compared to the bulk. Nevertheless, a remarkable degree of universality of distinctive dynamical behaviors is predicted for different polymer chemistries and all interfaces, including a double exponential form of the alpha time gradient, power law decoupling of the relaxation time from its bulk value with exponential spatial variation of the exponent, exponential spatial gradient of T_g, weak dependence of normalized T_g gradients on vitrification criterion, and near linear growth with cooling of the slowed down layer thickness near a rough hard interface. The results appear consistent with simulations and experiments, and multiple testable predictions are made.

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.

Design rules for glass formation from model molecules designed by a neural-network-biased genetic algorithm

The glass transition - an apparent amorphous solidification process - is a central feature of the physical properties of soft materials such as polymers and colloids. A key element of this phenomenon is the observation of a broad spectrum of deviations from an Arrhenius temperature of dynamics in glass-forming liquids, with the extent of deviation quantified by the "fragility" of glass formation. The underlying origin of "fragile" glass formation and its dependence on molecular structure remain major open questions in condensed matter physics and soft materials science. Here we employ molecular dynamics simulations, together with a neural-network-biased genetic algorithm, to design and study model rigid molecules spanning a broad range of fragilities of glass formation. Results indicate that fragility of glass formation can be controlled by tuning molecular asphericity, with extended molecules tending to exhibit low fragilities and compact molecules tending toward higher fragilities. The glass transition temperature itself, on the other hand, correlates well with high-temperature activation behavior and with density. These results point the way towards rational design of glass-forming liquids spanning a range of dynamical behavior, both via these physical insights and via future extensions of this evolutionary design strategy to real chemistries. Finally, we show that results compare well with predictions of the nonlinear Langevin theory of liquid dynamics, which is a precursor of the more recently developed elastically collective nonlinear Langevin equation theory of Mirigian and Schweizer, identifying this framework as a promising basis for molecular design of the glass transition.

Distinguishing failure modes in oligomeric polymer nanopillars

Brittle failure is ubiquitous in amorphous materials that are sufficiently cooled below their glass transition temperature, T_g. This catastrophic failure mode is limiting for amorphous materials in many applications, and many fundamental questions surrounding it remain poorly understood. Two challenges that prevent a more fundamental understanding of the transition between a ductile response at temperatures near T_g to brittle failure at lower temperatures are (i) a lack of computationally inexpensive molecular models that capture the failure modes observed in experiments and (ii) the lack of quantitative metrics that can distinguish various failure mechanisms. In this work, we use molecular dynamics simulations to capture ductile-to-brittle transition in glass-forming oligomeric polymer systems where we systematically vary both the temperature relative to T_g and the form of the interaction potential to induce a variety of failure modes. We characterized the effects of this new potential on macroscopic mechanical properties as well as microscopic structural and dynamical response during deformation. Finally, we develop several quantitative metrics to distinguish between different failure modes, and we find that the transition between catastrophic brittle failure, necking, and homogeneous plastic flow is gradual as the temperature is increased.

Influence of knot complexity on glass-formation in low molecular mass ring polymer melts

We perform molecular dynamics simulations on a coarse-grained polymer melt to study the dynamics of glass-formation in ring polymer melts of variable knot complexity. After generating melts of non-concatenated polymeric rings having a range of minimum crossing number values, mc, we compute the coherent intermediate scattering function, the segmental α-relaxation time, fragility, and the glass transition temperature as a function of mc. Variation of knot complexity is found to have a pronounced effect on the dynamics of polymer melts since both molecular rigidity and packing are altered, primary physical factors governing glass-formation in polymeric materials.

Polymer-Based Synthetic Routes to Carbon-Based Metal-Free Catalysts

Carbons are increasingly important as possible alternatives to expensive metal catalysts owing to the wide range of chemical properties they can exhibit and the growing set of synthetic routes available to produce them. This progress report discusses the process of making catalytic carbons from polymeric precursors, focusing on mechanisms of carbonization and how the polymer structures and synthetic procedures affect the resulting carbons. In considering what is necessary to move laboratory catalytic carbons to industrial and commercial applications, the cost and complexity to produce them are a considerable challenge to overcome. Industrially produced carbons are typically made from biopolymers such as lignin while many of the catalytic carbons studied in literature are from synthetic polymers. Thus, studying polymer-derived carbons can provide insights into the carbonization process and the properties of catalytic carbons, which can subsequently be translated to improve biopolymer-derived carbons in an economical way. Aspects of polymer carbonization discussed include carbonization mechanisms, effects of crosslinkers, polymer microstructure, heteroatom control, and effects of nanostructuring.

Different metrics for connecting mobility and glassiness in thin films

Data continue to accrue indicating that experimental techniques may differ in their sensitivity to mobility and glassiness. In this work the Limited Mobility (LM) kinetic model is used to show that two metrics for tracking sample mobility yield quantitatively different results for the glass transition and mobile layer thickness in systems where free surfaces are present. Both LM metrics track the fraction of material that embodies mobile free volume; in one it is relative to that portion of the sample containing any kind (mobile and dormant) of free volume, and in the other it is relative to the overall sample. Without any kind of optimization, use of the latter metric leads to semi-quantitative agreement with experimental film results, both for the mobile layer thickness and the dependence of sample glass transition temperature on film thickness. Connecting the LM predictions with experiment also produces a semi-quantitative mapping between LM model length and temperature scales, and those of real systems.

Universal localization transition accompanying glass formation: insights from efficient molecular dynamics simulations of diverse supercooled liquids

The origin of the precipitous dynamic arrest known as the glass transition is a grand open question of soft condensed matter physics. It has long been suspected that this transition is driven by an onset of particle localization and associated emergence of a glassy modulus. However, progress towards an accepted understanding of glass formation has been impeded by an inability to obtain data sufficient in chemical diversity, relaxation timescales, and spatial and temporal resolution to validate or falsify proposed theories for its physics. Here we first describe a strategy enabling facile high-throughput simulation of glass-forming liquids to nearly unprecedented relaxation times. We then perform simulations of 51 glass-forming liquids, spanning polymers, small organic molecules, inorganics, and metallic glass-formers, with longest relaxation times exceeding one microsecond. Results identify a universal particle-localization transition accompanying glass formation across all classes of glass-forming liquid. The onset temperature of non-Arrhenius dynamics is found to serve as a normalizing condition leading to a master collapse of localization data. This transition exhibits a non-universal relationship with dynamic arrest, suggesting that the nonuniversality of supercooled liquid dynamics enters via the dependence of relaxation times on local cage scale. These results suggest that a universal particle-localization transition may underpin the glass transition, and they emphasize the potential for recent theoretical developments connecting relaxation to localization and emergent elasticity to finally explain the origin of this phenomenon. More broadly, the capacity for high-throughput prediction of glass formation behavior may open the door to computational inverse design of glass-forming 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.

Experimental Test of the Cooperative Free Volume Rate Model under 1D Confinement: The Interplay of Free Volume, Temperature, and Polymer Film Thickness in Driving Segmental Mobility

We show that the Cooperative Free Volume (CFV) rate model, successful at modeling pressure-dependent dynamics, can be employed to describe the temperature and thickness dependence of the segmental time of polymers confined in thin films (1D confinement). The CFV model is based on an activation free energy that increases with the number of cooperating segments, which is determined by the system's free volume. Here, we apply the CFV model to new experimental results on the segmental relaxation of 1D confined poly(4-chlorostyrene), P4ClS, and find remarkable agreement over the whole temperature and thickness ranges investigated. This work further validates the robustness of the CFV model, which relates the effects of confinement on dynamics to pressure changes in the bulk, and supports the idea that confinement effects originate from local perturbations in density.

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.

Effect of substrate interactions on the glass transition and length-scale of correlated dynamics in ultra-thin molecular glass films

Interfacial interactions can play an important role in the glass transition temperature (T _g ) and relaxation dynamics of ultra-thin glass polymer films. We have recently shown that similar to the polymeric systems in ultra-thin molecular glass films of N, N'-Bis(3-methylphenyl)-N, N'-diphenylbenzidine (TPD), the T _g is reduced and dynamics are enhanced. Furthermore, in molecular glass systems, as the film thickness is reduced below ∼30 nm, the dynamics at the two interfaces correlate such that the range of the gradients in the dynamics induced by the free surface narrows compared to thicker films. These observations indicate that the dynamics of the glassy thin films are strongly correlated and cannot be explained by a simple two-layer model consisting of a bulk and a thin interfacial layer with a constant thickness and constant range of dynamical gradients. Here, we investigate the effect of film/substrate interactions on the film dynamics by varying the TPD/substrate interfacial interactions. We show that thin TPD films with thicknesses below ∼60 nm show a smaller extent of T _g reduction and enhanced dynamics when supported on a near-neutral substrate (wetting) compared to a weakly interacting (dewetting) substrate. However, the ∼30 nm length scale, where the activation energy significantly reduces from its bulk value as measured by the onset of the glass transition remains unchanged. Coarse-grained molecular dynamics simulation also shows a narrowing in the range of relaxation times once the thickness is sufficiently reduced for the two interfaces to dynamically correlate, consistent with previous work. These results suggest that the length-scale for the correlated dynamics is independent of interfacial interactions and the polymeric nature of the film and may originate from the bulk glass properties.

Monte Carlo simulation on the dynamics of a semi-flexible polymer in the presence of nanoparticles

The dynamics of a semi-flexible polymer chain in the presence of periodically distributed nanoparticles is simulated by using off-lattice Monte Carlo simulations. For repulsive or weak attractive nanoparticles, the dynamics are slowed down monotonically by increasing the chain stiffness k_θ or decreasing the inter-particle distance d. For strong attractive nanoparticles, however, the dynamics show nonmonotonic behaviors with k_θ and d. An interesting result is that a stiff polymer may move faster than a flexible one. The underlying mechanism is that the nanoparticle's attraction is weakened by the chain stiffness. The nonmonotonic behavior of the polymer's dynamics with k_θ is explained by the competition between the weakening effect of the chain stiffness on the nanoparticle's attraction and the intrinsic effect of chain stiffness which reduces the dynamics of the polymer. In addition, the nonmonotonic behavior of the polymer's dynamics with d is explained by the competition between the nanoparticle-exchange motion of the polymer dominated at small d and the desorption-and-adsorption motion at large d. The excluded volume effect of the nanoparticles plays a more important role for stiffer polymers as the attraction of the nanoparticles is weakened by the chain stiffness.

Investigation of nematic to smectic phase transition and dynamical properties of strongly confined semiflexible polymers using Langevin dynamics

We investigated the nematic to smectic phase transition for strongly confined semiflexible polymer solutions in slit-like confinements using GPU-accelerated Langevin dynamics. We characterized the phase transitions from the nematic to smectic phases for semi-flexible polymer solutions as the polymer density increased. The dependence for the lyotropic nematic to smectic transition can be collapsed by scaling exponents between 0.2 and 0.3. The smectic C phase is found for all the cases with the polymer orientation director tilted with respect to smectic layer lateral alignment. As the chain rigidity increases, the transition density decreases for systems in which the polymer persistence length (P) to slit height (H) ratios are 1.25, 2.5, 3.75, 5 and 25. We also characterized the polymer dynamics for the isotropic-nematic-smectic transitions. The overall polymer diffusivity decreased steadily as the polymer density increased. We observed anomalous polymer diffusion along the nematic director near the isotropic-nematic transition, similar to previously reported behavior for nematic-forming ellipsoids. Polymer diffusivity decreased sharply by two orders of magnitude upon the nematic-smectic transition.

Relaxation processes and glass transition of confined polymer melts: A molecular dynamics simulation of 1,4-polybutadiene between graphite walls

Molecular dynamics simulations of a chemically realistic model for 1,4-polybutadiene in a thin film geometry confined by two graphite walls are presented. Previous work on melts in the bulk has shown that the model faithfully reproduces static and dynamic properties of the real material over a wide temperature range. The present work studies how these properties change due to nano-confinement. The focus is on orientational correlations observable in nuclear magnetic resonance experiments and on the local intermediate incoherent neutron scattering function, F_s(q_z, z, t), for distances z from the graphite walls in the range of a few nanometers. Temperatures from about 2T_g down to about 1.15T_g, where T_g is the glass transition temperature in the bulk, are studied. It is shown that weakly attractive forces between the wall atoms and the monomers suffice to effectively bind a polymer coil that is near the wall. For a wide regime of temperatures, the Arrhenius-like adsorption/desorption kinetics of the monomers is the slowest process, while very close to T_g the Vogel-Fulcher-Tammann-like α-relaxation takes over. The α-process is modified only for z≤1.2 nm due to the density changes near the walls, less than expected from studies of coarse-grained (bead-spring-type) models. The weakness of the surface effects on the glass transition in this case is attributed to the interplay of density changes near the wall with the torsional potential. A brief discussion of pertinent experiments is given.

Side-group size effects on interfaces and glass formation in supported polymer thin films

Recent studies on glass-forming polymers near interfaces have emphasized the importance of molecular features such as chain stiffness, side-groups, molecular packing, and associated changes in fragility as key factors that govern the magnitude of T_g changes with respect to the bulk in polymer thin films. However, how such molecular features are coupled with substrate and free surface effects on T_g in thin films remains to be fully understood. Here, we employ a chemically specific coarse-grained polymer model for methacrylates to investigate the role of side-group volume on glass formation in bulk polymers and supported thin films. Our results show that bulkier side-groups lead to higher bulk T_g and fragility and are associated with a pronounced free surface effect on overall T_g depression. By probing local T_g within the films, however, we find that the polymers with bulkier side-groups experience a reduced confinement-induced increase in local T_g near a strongly interacting substrate. Further analyses indicate that this is due to the packing frustration of chains near the substrate interface, which lowers the attractive interactions with the substrate and thus lessens the surface-induced reduction in segmental mobility. Our results reveal that the size of the polymer side-group may be a design element that controls the confinement effects induced by the free surface and substrates in supported polymer thin films. Our analyses provide new insights into the factors governing polymer dynamics in bulk and confined environments.

Exploring the broadening and the existence of two glass transitions due to competing interfacial effects in thin, supported polymer films

In this report, we use ellipsometry to characterize the glass transition in ultra-thin films of poly(2-vinyl pyridine) (P2VP) supported on a silicon substrate. P2VP is known to have attractive substrate interactions, which can increase the T_g of ultra-thin films compared to the bulk value. Here, we use an extended temperature range to show that the glass transition can be extremely broad, indicating that a large gradient of the dynamics exists through the film with slow dynamics near the substrate and enhanced dynamics at the free surface. To observe the effect of these two interfaces on the average thin film dynamics, cooling rate-dependent T_g (CR-T_g) measurements were used to indirectly probe the average relaxation times of the films. We demonstrate that ultra-thin films have lower fragility compared to bulk films, and, when cooled at slow cooling rates (<1 K/min), exhibit extreme broadening of the dynamics (<70 nm) and eventually complete decoupling between the free surface and substrate regions to produce films with two distinct T_g's (<16 nm). T_g,high increases with decreasing thickness in a similar manner to what has been observed in previous studies on P2VP, and T_g,low decreases with decreasing film thickness in a similar manner to what has been observed in polymer films with enhanced free surfaces and neutral substrate interactions. These observations indicate that the dynamics in thin films of P2VP can be strongly coupled over a length scale of ∼10-20 nm, resulting in two co-existing layers with two distinct glass transitions when the range of the dynamical gradients become too large to sustain (breadth of the transition > 50 K).

Structure-property relationships from universal signatures of plasticity in disordered solids

When deformed beyond their elastic limits, crystalline solids flow plastically via particle rearrangements localized around structural defects. Disordered solids also flow, but without obvious structural defects. We link structure to plasticity in disordered solids via a microscopic structural quantity, "softness," designed by machine learning to be maximally predictive of rearrangements. Experimental results and computations enabled us to measure the spatial correlations and strain response of softness, as well as two measures of plasticity: the size of rearrangements and the yield strain. All four quantities maintained remarkable commonality in their values for disordered packings of objects ranging from atoms to grains, spanning seven orders of magnitude in diameter and 13 orders of magnitude in elastic modulus. These commonalities link the spatial correlations and strain response of softness to rearrangement size and yield strain, respectively.