Influence of Pressure on Glass Formation in a Simulated Polymer Melt

Influence of Density and Pressure on Glass Formation in the Kob-Andersen Model

We systematically study glass formation in the well-known Kob-Andersen model over a wide range of densities and pressures as a basis for judging the "universality" of glass formation through a comparison to a recent systematic study of model polymeric glass-forming liquids. Our purpose is to establish general characteristics of glass formation, to identify new relations, and to discern which properties of glass-forming liquids are material-specific and which are "universal." To this end, we analyze a number of characteristic properties of glass formation, such as the structural relaxation time, self-diffusion coefficient, viscosity, characteristic temperatures, and fragility. We also consider a suite of properties presumably related to dynamic heterogeneity in an attempt to better understand its relation to structural relaxation. We demonstrate that the glassy dynamics in the Kob-Andersen model exhibit many of the essential trends observed in polymeric glass-forming liquids, pointing to a remarkable "universality" of many aspects of glass formation.

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.

Mesoscopic two-point collective dynamics of glass-forming liquids

The collective density-density and hydrostatic pressure-pressure correlations of glass-forming liquids are spatiotemporally mapped out using molecular dynamics simulations. It is shown that the sharp rise of structural relaxation time below the Arrhenius temperature coincides with the emergence of slow, nonhydrodynamic collective dynamics on mesoscopic scales. The observed long-range, nonhydrodynamic mode is independent of wave numbers and closely coupled to the local structural dynamics. Below the Arrhenius temperature, it dominates the slow collective dynamics on length scales immediately beyond the first structural peak in contrast to the well-known behavior at high temperatures. These results highlight a key connection between the qualitative change in mesoscopic two-point collective dynamics and the dynamic crossover phenomenon.

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.

Understanding the role of cross-link density in the segmental dynamics and elastic properties of cross-linked thermosets

Cross-linking is known to play a pivotal role in the relaxation dynamics and mechanical properties of thermoset polymers, which are commonly used in structural applications because of their light weight and inherently strong nature. Here, we employ a coarse-grained (CG) polymer model to systematically explore the effect of cross-link density on basic thermodynamic properties as well as corresponding changes in the segmental dynamics and elastic properties of these network materials upon approaching their glass transition temperatures (T_g). Increasing the cross-link density unsurprisingly leads to a significant slowing down of the segmental dynamics, and the fragility K of glass formation shifts in lockstep with T_g, as often found in linear polymer melts when the polymer mass is varied. As a consequence, the segmental relaxation time τ_α becomes almost a universal function of reduced temperature, (T - T_g)/T_g, a phenomenon that underlies the applicability of the "universal" Williams-Landel-Ferry (WLF) relation to many polymer materials. We also test a mathematical model of the temperature dependence of the linear elastic moduli based on a simple rigidity percolation theory and quantify the fluctuations in the local stiffness of the network material. The moduli and distribution of the local stiffness likewise exhibit a universal scaling behavior for materials having different cross-link densities but fixed (T - T_g)/T_g. Evidently, T_g dominates both τ_α and the mechanical properties of our model cross-linked polymer materials. Our work provides physical insights into how the cross-link density affects glass formation, aiding in the design of cross-linked thermosets and other structurally complex glass-forming materials.

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.

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.

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.

Packing and dynamics of a protein solution approaching the jammed state

The packing of proteins and their collective behavior in crowded media is crucial for the understanding of biological processes. Here we study the structural and dynamical evolution of solutions of the globular protein bovine serum albumin with increasing concentration via drying using small angle X-ray scattering and dynamic light scattering. We probe an evolving correlation peak on the scattering profile, corresponding to the inter-protein distance, ξ, which decreases following a power law of the protein volume fraction, φ. The rate of decrease in ξ becomes faster above a protein concentration of ∼200 mg ml-1 (φ = 0.15). The power law exponent changes from 0.33, which is typical of colloidal or protein solutions, to 0.41. During the entire drying process, we observe the development and the growth of two-step relaxation dynamics with increasing φ as revealed by dynamic light scattering. We find three different regimes of the dependence of ξ as a function of φ. In the dilute regime (φ < 0.22), protein molecules are far apart from each other compared to their size. In this case, the dynamics mainly corresponds to Brownian motion. At an intermediate concentration (0.22 < φ < 0.47), inter-protein distances become comparable to the size of protein molecules, leading to a preferential orientation of the ellipsoidal protein molecules along with a possible deformation. In this regime, the dynamics shows two distinct relaxation times. At a very high concentration (φ > 0.47), the system reaches a jammed state. Subsequently, the secondary relaxation time in this state becomes extremely slow. In this state, the protein molecules have approximately one hydration layer. This study contributes to the understanding of protein molecular packing in crowded environments and the phenomenon of density-driven jamming for soft matter systems.

Universal nature of dynamic heterogeneity in glass-forming liquids: A comparative study of metallic and polymeric glass-forming liquids

Glass-formation is a ubiquitous phenomenon that is often observed in a broad class of materials ranging from biological matter to commonly encountered synthetic polymer, as well as metallic and inorganic glass-forming (GF) materials. Despite the many regularities in the dynamical properties of GF materials, the structural origin of the universal dynamical properties of these materials has not yet been identified. Recent simulations of coarse-grained polymeric GF liquids have indicated the coexistence of clusters of mobile and immobile particles that appear to be directly linked, respectively, to the rate of molecular diffusion and structural relaxation. The present work examines the extent to which these distinct types of "dynamic heterogeneity" (DH) arise in metallic GF liquids (Cu-Zr, Ni-Nb, and Pd-Si alloys) having a vastly different molecular structure and chemistry. We first identified mobile and immobile particles and their transient clusters and found the DH in the metallic alloys to be remarkably similar in form to polymeric GF liquids, confirming the "universality" of the DH phenomenon. Furthermore, the lifetime of the mobile particle clusters was found to be directly related to the rate of diffusion in these materials, while the lifetime of immobile particles was found to be proportional to the structural relaxation time, providing some insight into the origin of decoupling in GF liquids. An examination of particles having a locally preferred atomic packing, and clusters of such particles, suggests that there is no one-to-one relation between these populations of particles so that an understanding of the origin of DH in terms of static fluid structure remains elusive.

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.

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.