String-like collective motion in the α- and β-relaxation of a coarse-grained polymer melt

Influence of swelling on the elasticity of polymer networks cross-linked in the melt state: Test of the localization model of rubber elasticity

The elasticity of polymer networks, formed by cross-linking high molecular mass polymers in the melt state and then swollen by a solvent, involves contributions from both the presence of cross-link network junctions and the interchain interactions associated with the combined effect of excluded volume interactions and topological constraints that become modified when the network is swollen. We test the capacity of the previously developed localization model of rubber elasticity, a mean field "tube model," to describe changes in elasticity observed in classical experimental studies of the mechanical properties of this type of network. In order to obtain a satisfactory comparison to the experiments, it was found to be necessary to account for the independently observed tendency of the network junctions to become localized in the network with appreciable swelling, as well as the interchain interactions emphasized in previous discussions of the localization model.

Stringlet excitation model of the boson peak

The boson peak (BP), a low-energy excess in the vibrational density of states over the Debye contribution, is often identified as a characteristic of amorphous solid materials. Despite decades of efforts, its microscopic origin still remains a mystery. Recently, it has been proposed, and corroborated with simulations, that the BP might stem from intrinsic localized modes involving one-dimensional (1D) string-like excitations ("stringlets"). We build on a theory originally proposed by Lund that describes the localized modes as 1D vibrating strings, but we specify the stringlet size distribution to be exponential, as observed in simulations. We provide an analytical prediction for the BP frequency ωBP in the temperature regime well below the observed glass transition temperature Tg. The prediction involves no free parameters and accords quantitatively with prior simulation observations in 2D and 3D model glasses based on inverse power law potentials. The comparison of the string model to observations is more uncertain when compared to simulations of an Al-Sm metallic glass material at temperatures well above Tg. Nonetheless, our stringlet model of the BP naturally reproduces the softening of the BP frequency upon heating and offers an analytical explanation for the experimentally observed scaling with the shear modulus in the glass state and changes in this scaling in simulations of glass-forming liquids. Finally, the theoretical analysis highlights the existence of a strong damping for the stringlet modes above Tg, which leads to a large low-frequency contribution to the 3D vibrational density of states, observed in both experiments and simulations.

A closer examination of the nature of atomic motion in the interfacial region of crystals upon approaching melting

Although crystalline materials are often conceptualized as involving a static lattice configuration of particles, it has recently become appreciated that string-like collective particle exchange motion is a ubiquitous and physically important phenomenon in both the melting and interfacial dynamics of crystals. This type of collective motion has been evidenced in melting since early simulations of hard disc melting by Alder et al. [Phys. Rev. Lett. 11(6), 241-243 (1963)], but a general understanding of its origin, along with its impact on melting and the dynamics of crystalline materials, has been rather slow to develop. We explore this phenomenon further by focusing on the interfacial dynamics of a model crystalline Cu material using molecular dynamics simulations where we emphasize the geometrical nature and spatial extent of the atomic trajectories over the timescale that they are caged, and we also quantify string-like collective motion on the timescale of the fast β-relaxation time, τf, i.e., "stringlets." Direct visualization of the atomic trajectories in their cages over the timescale over which the cage persists indicates that they become progressively more anisotropic upon approaching the melting temperature Tm. The stringlets, dominating the large amplitude atomic motion in the fast dynamics regime, are largely localized to the crystal interfacial region and correspond to "excess" modes in the density of states that give rise to a "boson peak." Moreover, interstitial point defects occur in direct association with the stringlets, demonstrating a link between classical defect models of melting and more recent studies of melting emphasizing the role of this kind of collective motion.

Solid-that-Flows Picture of Glass-Forming Liquids

This perspective article reviews arguments that glass-forming liquids are different from those of standard liquid-state theory, which typically have a viscosity in the mPa·s range and relaxation times on the order of picoseconds. These numbers grow dramatically and become 1012 - 1015 times larger for liquids cooled toward the glass transition. This translates into a qualitative difference, and below the "solidity length" which is roughly one micron at the glass transition, a glass-forming liquid behaves much like a solid. Recent numerical evidence for the solidity of ultraviscous liquids is reviewed, and experimental consequences are discussed in relation to dynamic heterogeneity, frequency-dependent linear-response functions, and the temperature dependence of the average relaxation time.

Confinement effect of inter-arm interactions on glass formation in star polymer melts

We utilized molecular dynamic simulation to investigate the glass formation of star polymer melts in which the topological complexity is varied by altering the number of star arms (f). Emphasis was placed on how the "confinement effect" of repulsive inter-arm interactions within star polymers influences the thermodynamics and dynamics of star polymer melts. All the characteristic temperatures of glass formation were found to progressively increase with increasing f, but unexpectedly the fragility parameter KVFT was found to decrease with increasing f. As previously observed, stars having more than 5 or 6 arms adopt an average particle-like structure that is more contracted relative to the linear polymer size having the same mass and exhibit a strong tendency for intermolecular and intramolecular segregation. We systematically analyzed how varying f alters collective particle motion, dynamic heterogeneity, the decoupling exponent ζ phenomenologically linking the slow β- and α-relaxation times, and the thermodynamic scaling index γt. Consistent with our hypothesis that the segmental dynamics of many-arm star melts and thin supported polymer films should exhibit similar trends arising from the common feature of high local segmental confinement, we found that ζ increases considerably with increasing f, as found in supported polymer films with decreasing thickness. Furthermore, increasing f led to greatly enhanced elastic heterogeneity, and this phenomenon correlates strongly with changes in ζ and γt. Our observations should be helpful in building a more rational theoretical framework for understanding how molecular topology and geometrical confinement influence the dynamics of glass-forming materials more broadly.

Competing Effects of Molecular Additives and Cross-Link Density on the Segmental Dynamics and Mechanical Properties of Cross-Linked Polymers

The introduction of molecular additives into thermosets often results in changes in their dynamics and mechanical properties that can have significant ramifications for diverse applications of this broad class of materials such as coatings, high-performance composites, etc. Currently, there is limited fundamental understanding of how such additives influence glass formation in these materials, a problem of broader significance in glass-forming materials. To address this fundamental problem, here, we employ a simplified coarse-grained (CG) model of a polymer network as a model of thermoset materials and then introduce a polymer additive having the same inherent rigidity and polymer-polymer interaction strength as the cross-linked polymer matrix. This energetically "neutral" or "self-plasticizing" additive model gives rise to non-trivial changes in the dynamics of glass formation and provides an important theoretical reference point for the technologically more important case of interacting additives. Based on this rather idealized model, we systematically explore the combined effect of varying the additive mass percentage (m) and cross-link density (c) on the segmental relaxation dynamics and mechanical properties of a model thermoset material with additives. We find that increasing the additive mass percentage m progressively decreases both the glass-transition temperature Tg and the fragility of glass formation, a trend opposite to increasing c so that these thermoset variables clearly have a competing effect on glass formation in these model materials. Moreover, basic mechanical properties (i.e., bulk, shear, and tensile moduli) likewise exhibit a competitive variation with the increase of m and c, which are strongly correlated with the Debye-Waller parameter ⟨u2⟩, a measure of material stiffness at a molecular scale. Our findings prove beneficial in the development of structure-property relationships for the cross-linked polymers, which could help guide the design of such network materials with tailored physical properties.

Approach to hyperuniformity in a metallic glass-forming material exhibiting a fragile to strong glass transition

We investigate a metallic glass-forming (GF) material (Al₉₀Sm₁₀) exhibiting a fragile-strong (FS) glass-formation by molecular dynamics simulation to better understand this highly distinctive pattern of glass-formation in which many of the usual phenomenological relations describing relaxation times and diffusion of ordinary GF liquids no longer apply, and where instead genuine thermodynamic features are observed in response functions and little thermodynamic signature is exhibited at the glass transition temperature, T_g. Given the many unexpected similarities between the thermodynamics and dynamics of this metallic GF material with water, we first focus on the anomalous static scattering in this liquid, following recent studies on water, silicon and other FS GF liquids. We quantify the "hyperuniformity index" H of our liquid, which provides a quantitative measure of molecular "jamming". To gain insight into the T-dependence and magnitude of H, we also estimate another more familiar measure of particle localization, the Debye-Waller parameter 〈u²〉 describing the mean-square particle displacement on a timescale on the order of the fast relaxation time, and we also calculate H and 〈u²〉 for heated crystalline Cu. This comparative analysis between H and 〈u²〉 for crystalline and metallic glass materials allows us to understand the critical value of H on the order of 10^-3 as being analogous to the Lindemann criterion for both the melting of crystals and the "softening" of glasses. We further interpret the emergence of FS GF and liquid-liquid phase separation in this class of liquids to arise from a cooperative self-assembly process in the GF liquid.

The Phonon Theory of Liquids and Biological Fluids: Developments and Applications

Among the three basic states of matter (solid, liquid, and gas), the liquid state has always eluded general theoretical approaches for describing liquid energy and heat capacity. In this Viewpoint, we derive the phonon theory of liquids and biological fluids stemming from Frenkel's microscopic picture of the liquid state. Specifically, the theory predicts the existence of phonon gaps in vibrational spectra of liquids and a thermodynamic boundary in the supercritical state. Direct experimental evidence reaffirming these theoretical predictions was achieved through a combination of techniques using static compression X-ray diffraction and inelastic X-ray scattering on deeply supercritical argon in a diamond anvil cell. Furthermore, these findings inspired and then led to the discovery of phonon gaps in liquid crystals (mesogens), block copolymers, and biological membranes. Importantly, phonon gaps define viscoelastic crossovers in cellular membranes responsible for lipid self-diffusion, lateral molecular-level stress propagation, and passive transmembrane transport of small molecules and solutes. Finally, molecular interactions mediated by external stimuli result in synaptic activity controlling biological membranes' plasticity resulting in learning and memory. Therefore, we also discuss learning and memory effects─equally important for neuroscience as well as for the development of neuromorphic devices─facilitated in biological membranes by external stimuli.

Construction of a quantitative relation between structural relaxation and dynamic heterogeneity by vibrational dynamics in glass-forming liquids and polymers

The structural relaxation slows down drastically upon approaching the glass transition, accompanied by the significant growth of dynamic heterogeneity. The fundamental question of elusiveness and interest is whether there exists an underlying quantitative relationship between structural relaxation and dynamic heterogeneity. Here, we reveal that b̃ which is related to the reduced mean square displacements to overcome the energy barriers of activated jumps, instead of the kinetic fragility m, is the genuine key parameter connecting dynamic heterogeneity with structural relaxation for varying types of glass formers. Furthermore, based on the dependence of dynamic heterogeneity on the Debye-Waller factor we obtained a direct quantitative relation between dynamic heterogeneity and structural relaxation is built for different glass-forming liquids. More importantly, a scaling collapse of structural relaxation and dynamic heterogeneity is achieved by the important parameter b̃. These results are of fundamental and critical importance for developing a unified theory of glassy dynamics.

The initiation of shear band formation in deformed metallic glasses from soft localized domains

It has long been thought that shear band (SB) formation in amorphous solids initiates from relatively "soft" regions in the material in which large-scale non-affine deformations become localized. The test of this hypothesis requires an effective means of identifying "soft" regions and their evolution as the material is deformed to varying degrees, where the metric of "softness" must also account for the effect of temperature on local material stiffness. We show that the mean square atomic displacement on a caging timescale ⟨u²⟩, the "Debye-Waller factor," provides a useful method for estimating the shear modulus of the entire material and, by extension, the material stiffness at an atomic scale. Based on this "softness" metrology, we observe that SB formation indeed occurs through the strain-induced formation of localized soft regions in our deformed metallic glass free-standing films. Unexpectedly, the critical strain condition for SB formation occurs when the softness (⟨u²⟩) distribution within the emerging soft regions approaches that of the interfacial region in its undeformed state, initiating an instability with similarities to the transition to turbulence. Correspondingly, no SBs arise when the material is so thin that the entire material can be approximately described as being "interfacial" in nature. We also quantify relaxation in the glass and the nature and origin of highly non-Gaussian particle displacements in the dynamically heterogeneous SB regions at times longer than the caging time.

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.

A Dynamically Correlated Network Model for the Collective Dynamics in Glass-Forming Molecular Liquids and Polymers

The non-Arrhenius behavior of segmental dynamics in glass-forming liquids is one of the most profound mysteries in soft matter physics. In this article, we propose a dynamically correlated network (DCN) model to understand the growing behavior of dynamically correlated regions during cooling, which leads to the viscous slowdown of supercooled liquids. The fundamental concept of the model is that the cooperative region of collective motions has a network structure that consists of string-like parts, and networks of various sizes interpenetrate each other. Each segment undergoes dynamical coupling with its neighboring segments via a finite binding energy. Monte Carlo simulations showed that the fractal dimension of the DCNs generated at different temperatures increased and their size distribution became broader with decreasing temperature. The segmental relaxation time was evaluated based on a power law with four different exponents for the activation energy of rearrangement with respect to the DCN size. The results of the present DCN model are consistent with the experimental results for various materials of molecular and polymeric liquids.

The Relationship between Free Volume and Cooperative Rearrangement: From the Temperature-Dependent Neutron Total Scattering Experiment of Polystyrene

Although many theories have been proposed to describe the nature of glass formation, its microscopic picture is still missing. Here, by a combination of neutron scattering and molecular dynamics simulation, we present the temperature-dependent atomic structure variation of polystyrene at the glass formation, free volume and cooperative rearrangement. When it is close to glass formation, the polymer is confined in tubes, whose diameter is the main chain–main chain distance, in a “static cage” from its neighbors. This definition can not only account for the kinetic pathway dependence of Williams-Landel-Ferry (WLF) free volume, but also be testified in a set of six polymers. However, the free volume which allows a monomer to move cannot be found in any frame of its real-space image. Monomers, thus, have to move cooperatively to be out of the cage. During glass formation, dynamic heterogeneity develops, and string-like cooperative rearrangement region (CRR) grows over a long range of time and length scales. All of these CRRs tend to walk through loose “static cages”. Our observation unifies the concepts of free volume and cooperative rearrangement. The former is a statistical average leading to a polydisperse “static cage” formation; while a loose “static cage” provides the way that CRRs move.

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

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

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

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

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

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

The polarizability response of a glass-forming liquid reveals intrabasin motion and interbasin transitions on a potential energy landscape

Potential energy landscape (PEL) concepts have been useful in conceptualizing the effects of intermolecular interactions on dynamic and thermodynamic properties of liquids and glasses. "Basins", or regions of reduced potential energy associated with locally preferred molecular packing are important PEL features. The molecular configurations at the bottom of these basins are referred to as inherent structures (ISs). Experimental methods for directly characterizing PEL features such as these are rare, largely relegating PEL concepts to theory and simulation studies, and impeding their exploration in real systems. Recently, we showed that quasielastic neutron scattering (QENS) data from propylene carbonate (PC) exhibit signatures of picosecond timescale motion that are consistent with intrabasin motion and interbasin transitions [Cicerone et al., J. Chem. Phys., 2017, 146, 054502]. Here we present optically-heterodyne-detected optical Kerr effect (OHD-OKE) spectroscopy studies on PC. The data exhibit signatures of motion within and transitions between basins that agree quantitatively with and extend the QENS results. We show that the librational component of the OKE response corresponds to intrabasin dynamics, and the enigmatic intermediate OKE response corresponds to interbasin transition events. The OKE data extend the measurement range of these parameters and reveal their utility in characterizing PEL features of real systems.

Steady state dynamic dependence between local mobility and non-affine fluctuations in two-dimensional aggregates

Motivated by qualitative experimental observations in collective behavior of self-propelled camphor particles at air-water interfaces, we study a generic aggregate forming system in two dimensions using canonical ensemble constant temperature molecular dynamics simulation. The aggregates form due to the competition between short-range attraction and long-range repulsion of pair-wise interactions as a generic proxy for the specific case of short-range capillary attraction competing with long-range Marangoni-assisted repulsion in camphor boat systems. Choosing the appropriate set of interaction parameters, we focus on characterising the local dynamics in two specific limiting morphologies, viz. compact and string-like aggregates. We focus on the temporal evolution of the mobility of an individual particle and the dynamic change in its nearest neighbourhood, measured in terms of the Debye-Waller factor ([Formula: see text]) and the non-affine parameter ([Formula: see text]), respectively (both defined in the text), and their interrelation over several lengths of observation time [Formula: see text]. The distribution for both measures are found to follow the relation: [Formula: see text] for the measured quantity x. The exponent [Formula: see text] is equal to two and one respectively, for the compact and string-like morphologies following the respective ideal fractal dimension of these aggregates. A functional dependence between these two observables is determined from a detailed statistical analysis of their joint and conditional distributions. The results obtained can readily be used and verified by experiments on aggregate forming systems more generic than the specific camphor boat system that motivated us, such as globular proteins, nanoparticle self-assembly etc. Further, the insights gained from this study might be useful to understand the evolution of collective dynamics in diverse glass-forming systems.

Predictive relation for the α-relaxation time of a coarse-grained polymer melt under steady shear

We examine the influence of steady shear on structural relaxation in a simulated coarse-grained unentangled polymer melt over a wide range of temperature and shear rates. Shear is found to progressively suppress the α-relaxation process observed in the intermediate scattering function, leading ultimately to a purely inertially dominated β-relaxation at high shear rates, a trend similar to increasing temperature. On the basis of a scaling argument emphasizing dynamic heterogeneity in cooled liquids and its alteration under material deformation, we deduce and validate a parameter-free scaling relation for both the structural relaxation time τα from the intermediate scattering function and the "stretching exponent" β quantifying the extent of dynamic heterogeneity over the entire range of temperatures and shear rates that we can simulate.

Do String-like Cooperative Motions Predict Relaxation Times in Glass-Forming Liquids?

The Adam-Gibbs theory of glass formation posits that the growth in the activation barrier of fragile liquids on cooling emerges from a loss of configurational entropy and concomitant growth in "cooperatively rearranging regions" (CRRs). A body of literature over 2 decades has suggested that "string-like" cooperatively rearranging clusters observed in molecular simulations may be these CRRs-a scenario that would have profound implications for the understanding of the glass transition. The central element of this postulate is the report of an apparent zero-parameter relationship between the mass of string-like CRRs and the relaxation time. Here, we show, based on molecular dynamics simulations of multiple glass-forming liquids, that this finding is the result of an implicit adjustable parameter-a "replacement distance". This parameter is equivalent to an adjustable exponent within a generalized Adam-Gibbs relation, such that it tunes the entire functional form of the relation. Moreover, we are unable to find any objective criterion, based on the radial distribution function or the cluster fractal dimension, for selecting this replacement distance across multiple systems. We conclude that the present data do not establish that string-like cooperative rearrangements, as presently defined, are predictive of segmental relaxation via an Adam-Gibbs-like physical model.

Collective Motion in the Interfacial and Interior Regions of Supported Polymer Films and Its Relation to Relaxation

To understand the role of collective motion in the often large changes in interfacial molecular mobility observed in polymer films, we investigate the extent of collective motion in the interfacial regions of a thin supported polymer film and within the film interior by molecular dynamics simulation. Contrary to commonly stated expectations, we find that the extent of collective motion, as quantified by string-like molecular exchange motion, is similar in magnitude in the polymer-air interfacial layer as the film interior and distinct from the bulk material. This finding is consistent with Adam-Gibbs description of the segmental dynamics within mesoscopic film regions, where the extent of collective motion is related to the configurational entropy of the film as a whole rather than a locally defined extent of collective motion or configurational entropy.

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

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

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.