Disorder-assisted melting and the glass transition in amorphous solids

Fragility and thermal expansion control crystal melting and the glass transition

Analytical relations for the glass transition temperature, Tg, and the crystal melting temperature, Tm, are developed on the basis of nonaffine lattice dynamics. The proposed relations explain the following: (i) the seemingly universal factor of ≈2/3 difference between the glass transition temperature and the melting temperature of the corresponding crystal, and (ii) the recent empirical discovery that both Tg and Tm are proportional to the liquid fragility m divided by the thermal expansion coefficient α of the solid.

Hyperuniform disordered solids with crystal-like stability

Hyperuniform disordered solids, characterised by unusually suppressed density fluctuations at low wavenumbers (q), are of great interest due to their potentially distinct properties as a unique glass state. From the jamming perspective, there is ongoing debate about the relationship between hyperuniformity and the jamming transition, as well as whether hyperuniformity persists above the jamming point. Here, we successfully generate over-jammed disordered solids exhibiting the strongest class of hyperuniformity, characterised by a power-law density spectrum (qα with α = 4). By decompressing both hyperuniform and conventional over-jammed packings to their respective marginally jammed states, we identify protocol-independent exponents: α ≈ 0.25 for density hyperuniformity and α ≈ 2 for contact-number hyperuniformity, both associated with the jamming transition. Although both marginally jammed and conventional over-jammed packings exhibit marginal stability, we demonstrate that hyperuniform over-jammed packings possess exceptional stability across vibrational, kinetic, thermodynamic, and mechanical properties-similar to crystals. These findings suggest that hyperuniform over-jammed packings offer crucial insights into the ideal disordered solid state and stand out as promising candidates for disordered metamaterials, uniquely combining hyperuniformity with ultrastability.

Learning glass transition temperatures via dimensionality reduction with data from computer simulations: Polymers as the pilot case

Machine learning methods provide an advanced means for understanding inherent patterns within large and complex datasets. Here, we employ the principal component analysis (PCA) and the diffusion map (DM) techniques to evaluate the glass transition temperature (Tg) from low-dimensional representations of all-atom molecular dynamic simulations of polylactide (PLA) and poly(3-hydroxybutyrate) (PHB). Four molecular descriptors were considered: radial distribution functions (RDFs), mean square displacements (MSDs), relative square displacements (RSDs), and dihedral angles (DAs). By applying Gaussian Mixture Models (GMMs) to analyze the PCA and DM projections and by quantifying their log-likelihoods as a density-based metric, a distinct separation into two populations corresponding to melt and glass states was revealed. This separation enabled the Tg evaluation from a cooling-induced sharp increase in the overlap between log-likelihood distributions at different temperatures. Tg values derived from the RDF and MSD descriptors using DM closely matched the standard computer simulation-based dilatometric and dynamic Tg values for both PLA and PHB models. This was not the case for PCA. The DM-transformed DA and RSD data resulted in Tg values in agreement with experimental ones. Overall, the fusion of atomistic simulations and DMs complemented with the GMMs presents a promising framework for computing Tg and studying the glass transition in a unified way across various molecular descriptors for glass-forming materials.

Structural evolution during inverse vulcanization

Inverse vulcanization exploits S8 to synthesize polysulfides. However, evolution of products and its mechanism during inverse vulcanization remains elusive. Herein, inverse vulcanization curves are obtained to describe the inverse vulcanization process in terms of three stages: induction, curing and over-cure. The typical curves exhibit a moduli increment before declining or plateauing, reflecting the process of polysulfide network formation and loosing depending on monomers. For aromatic alkenes, in the over-cure, the crosslinked polysulfide evolves significantly into a sparse network with accelerated relaxation, due to the degradation of alkenyl moieties into thiocarbonyls. The inverse vulcanization product of olefins degrades slowly with fluctuated relaxation time and modulus because of the generation of thiophene moieties, while the inverse vulcanization curve of dicyclopentadiene has a plateau following curing stage. Confirmed by calculations, the mechanisms reveal the alkenyl groups react spontaneously into thiocarbonyls or thiophenes via similar sulfur-substituted alkenyl intermediates but with different energy barriers.

Unraveling the rheology of inverse vulcanized polymers

Multiple relaxation times are used to capture the numerous stress relaxation modes found in bulk polymer melts. Herein, inverse vulcanization is used to synthesize high sulfur content (≥50 wt%) polymers that only need a single relaxation time to describe their stress relaxation. The S-S bonds in these organopolysulfides undergo dissociative bond exchange when exposed to elevated temperatures, making the bond exchange dominate the stress relaxation. Through the introduction of a dimeric norbornadiene crosslinker that improves thermomechanical properties, we show that it is possible for the Maxwell model of viscoelasticity to describe both dissociative covalent adaptable networks and living polymers, which is one of the few experimental realizations of a Maxwellian material. Rheological master curves utilizing time-temperature superposition were constructed using relaxation times as nonarbitrary horizontal shift factors. Despite advances in inverse vulcanization, this is the first complete characterization of the rheological properties of this class of unique polymeric material.

Molecular-Level Relation between the Intraparticle Glass Transition Temperature and the Stability of Colloidal Suspensions

In many colloidal suspensions, the dispersed colloidal particles are amorphous solids, resulting from vitrification. A crucial open problem is understanding how colloidal stability is affected by the intraparticle glass transition. By dealing with the latter process from a solid-state perspective, we estabilish a proportionality relation between the intraparticle glass transition temperature, Tg, and the Hamaker constant, AH, of a generic suspension of nanoparticles. It follows that Tg can be used as a convenient parameter (alternative to AH) for controlling the stability of colloidal systems. Within the Derjaguin-Landau-Verwey-Overbeek theory, we show that the novel relationship, connecting Tg to AH, implies the critical coagulation ionic strength to be a monotonically decreasing function of Tg. We connect our predictions to recent experimental findings.

Bond Orientation-Determined Enthalpic Stress in Polymer Glasses upon Deformation

The mechanical properties of polymer glass are determined by both intermolecular local packing structures and aligned intrachain configurations. These configurations involve multiple space scales, and the underlying mechanism is not well understood yet. By applying mechanical stimulation to cold-drawn polymer glasses, the present simulation work shows a one-to-one correspondence between arising retractive stress and the segment orientation parameter on the length scale of the intrachain connecting bond. Such retractive stress is a newly produced enthalpic stress when segment orientation on the length scale of bonds and particle mobility coexist. This reveals a potential mechanism of how the intrachain orientation on the length scale of bonds influences the mechanical behaviors of polymer glasses.

A double-well potential model for glass transition in a glassy hydrogel undergoing bi-stable interactions with water

The key mechanisms for achieving ultra-high mechanical properties of glassy hydrogels have not been fully understood, and it is commonly believed that their glass transitions are the crucial reasons due to the existence of significant bi-stable interactions between polymer macromolecules and water molecules. In this study, a double-well potential model is formulated to describe the mechanical properties of glassy hydrogels undergoing glass transition, by combining phase evolution theory and a rubber elasticity model. Bi-stable interactions between polymer macromolecules and water molecules (for both the trapped and free water) have been characterized using this double-well potential model, and various parameters are studied, including depth of well (for elasticity), distance between two wells (for yielding), and energy difference between two wells (for transition probability). Furthermore, constitutive stress-strain relationships are developed to explore the working principles for achieving ultra-high mechanical properties of these glassy hydrogels. Finally, the effectiveness of the proposed models is verified using finite element analysis (FEA) and also the experimental results reported in the literature, thus providing physical and mechanical insights into the ultra-high mechanical properties of glassy hydrogels.

Probing the Free Volume in Polymers by Means of Positron Annihilation Lifetime Spectroscopy

Positron annihilation lifetime spectroscopy (PALS) is a valuable technique to investigate defects in solids, such as vacancy clusters and grain boundaries in metals and alloys, as well as lattice imperfections in semiconductors. Positron spectroscopy is able to reveal the size, structure and concentration of vacancies with a sensitivity of 10^-7. In the field of porous and amorphous systems, PALS can probe cavities in the range from a few tenths up to several tens of nm. In the case of polymers, PALS is one of the few techniques able to give information on the holes forming the free volume. This quantity, which cannot be measured with macroscopic techniques, is correlated to important mechanical, thermal, and transport properties of polymers. It can be deduced theoretically by applying suitable equations of state derived by cell models, and PALS supplies a quantitative measure of the free volume by probing the corresponding sub-nanometric holes. The system used is positronium (Ps), an unstable atom formed by a positron and an electron, whose lifetime can be related to the typical size of the holes. When analyzed in terms of continuous lifetimes, the positron annihilation spectrum allows one to gain insight into the distribution of the free volume holes, an almost unique feature of this technique. The present paper is an overview of PALS, addressed in particular to readers not familiar with this technique, with emphasis on the experimental aspects. After a general introduction on free volume, positronium, and the experimental apparatus needed to acquire the corresponding lifetime, some of the recent results obtained by various groups will be shown, highlighting the connections between the free volume as probed by PALS and structural properties of the investigated materials.

Combined description of pressure-volume-temperature and dielectric relaxation of several polymeric and low-molecular-weight organic glass-formers using SL-TS2 approach

We apply our recently-developed mean-field "SL-TS2" (two-state Sanchez-Lacombe) model to simultaneously describe dielectric α-relaxation time, τ_α, and pressure-volume-temperature (PVT) data in four polymers (polystyrene, poly(methylmethacrylate), poly(vinyl acetate) and poly(cyclohexane methyl acrylate)) and four organic molecular glass formers (ortho-terphenyl, glycerol, PCB-62, and PDE). Previously, it has been shown that for all eight materials, the Casalini-Roland thermodynamical scaling, τ_α = f(Tvγsp) (where T is temperature and v_sp is specific volume) is satisfied (R. Casalini and C. M. Roland, Phys. Rev. E, 2004, 69(6), 62501). It has also been previously shown that the same scaling emerges naturally (for sufficiently low pressures) within the "SL-TS2" framework (V. V. Ginzburg, Soft Matter, 2021, 17, 9094-9106). Here, we fit the ambient pressure curves for the relaxation time and the specific volume as functions of temperature for the eight materials and observe a good agreement between theory and experiment. We then use the Casalini-Roland scaling to convert those results into "master curves", thus enabling predictions of relaxation times and specific volumes at elevated pressures. The proposed approach can be used to describe other glass-forming materials, both low-molecular-weight and polymeric.

Kinetics of inherent processes counteracting crystallization in supercooled monatomic liquid

Crystallization of supercooled liquids is mainly determined by two competing processes associated with the transition of particles (atoms) from liquid phase to crystalline one and, vice versa, with the return of particles from crystalline phase to liquid one. The quantitative characteristics of these processes are the so-called attachment rateg+and the detachment rateg-, which determine how particles change their belonging from one phase to another. In the present study, acorrespondence rulebetween the ratesg+andg-as functions of the sizeNof growing crystalline nuclei is defined for the first time. In contrast to the well-known detailed balance condition, which relatesg+(N)andg-(N)atN=nc(wheren_cis the critical nucleus size) and is satisfied only at the beginning of the nucleation regime, the foundcorrespondence ruleis fulfilled at all the main stages of crystallization kinetics (crystal nucleation, growth and coalescence). On the example of crystallizing supercooled Lennard-Jones liquid, the rateg-was calculated for the first time at different supercooling levels and for the wide range of nucleus sizesN∈[nc;40nc]. It was found that for the whole range of nucleus sizes, the detachment rateg-is only≈2% less than the attachment rateg+. This is direct evidence that the role of the processes that counteract crystallization remains significant at all the stages of crystallization. Based on the obtained results, a kinetic equation was formulated for the time-dependent distribution function of the nucleus sizes, that is an alternative to the well-known kinetic Becker-Döring-Zeldovich-Frenkel equation.

Decoding polymer self-dynamics using a two-step approach

The self-correlation function and corresponding self-intermediate scattering function in Fourier space are important quantities for describing the molecular motions of liquids. This work draws attention to a largely overlooked issue concerning the analysis of these space-time density-density correlation functions of polymers. We show that the interpretation of non-Gaussian behavior of polymers is generally complicated by intrachain averaging of distinct self-dynamics of different segments. By the very nature of the mathematics involved, the averaging process not only conceals critical dynamical information, but also contributes to the observed non-Gaussian dynamics. To fully expose this issue and provide a thorough benchmark of polymer self-dynamics, we perform analyses of coarse-grained molecular dynamics simulations of linear and ring polymer melts as well as several theoretical models using a "two-step" approach, where interchain and intrachain averagings of segmental self-dynamics are separated. While past investigations primarily focused on the average behavior, our results indicate that a more nuanced approach to polymer self-dynamics is clearly required.

Molecular signatures of the glass transition in polymers

The glass transition temperature (T_{g}) is one of the most fundamental properties of polymers. T_{g} is predicted by some theories as a sudden change in a "macroscopic" quantity (e.g., compressibility). However, for systems with "soft" glass transitions where the change is gradual it becomes hard to pinpoint precisely the transition temperature as well as the set of molecular changes occurring during this transition. Here, we introduce two new molecular signatures for the glass transition of polymers that exhibit clear changes as one approaches T_{g}: (i) differential change of the probability distribution of dihedral angles as a function of temperature and (ii) the distribution of fractional of the time spent in the different torsional states. These new signatures provide insights into the glass transition in polymers by directly exhibiting the concept of spatial heterogeneity and dynamical ergodicity breaking in such systems, as well as provide a key step to quantitatively obtain the transition temperature from molecular characteristics of the polymeric systems.

Structural modification enhances the optoelectronic properties of defect blue phosphorene thin films

Active enhancement of the optical absorption coefficient to improve the light converting efficiency of thin-film solar cell materials is crucial to develop the next-generation solar cell devices. Here we report first-principles calculations with generalized gradient approximation to study the optoelectronic properties of pristine and divacancy (DV) blue phosphorene (BlueP) thin films under structural deformation. We show that instead of formingsp-like covalent bonds as in the pristine BlueP layer, a DV introduces two particular dangling bonds between the voids. Using a microscopic (non-) affine deformation model, we reveal that the orbital hybridization of these dangling bonds is strongly modified in both the velocity and vorticity directions depending on the type of deformation, creating an effective light trap to enhance the material absorption efficiency. Furthermore, this successful light trap is complemented by a clear signature ofσ+πplasmon when a DV BlueP layer is slightly compressive. These results demonstrate a practical approach to tailor the optoelectronic properties of low-dimensional materials and to pave a novel strategy to design functionalized solar cell devices from the bottom-up with selective defects.

Molecular rotors to probe the local viscosity of a polymer glass

We investigate the local viscosity of a polymer glass around its glass transition temperature using environment-sensitive fluorescent molecular rotors embedded in the polymer matrix. The rotors' fluorescence depends on the local viscosity, and measuring the fluorescence intensity and lifetime of the probe therefore allows to measure the local free volume in the polymer glass when going through the glass transition. This also allows us to study the local viscosity and free volume when the polymer film is put under an external stress. We find that the film does not flow homogeneously, but undergoes shear banding that is visible as a spatially varying free volume and viscosity.

Combined description of polymer PVT and relaxation data using a dynamic "SL-TS2" mean-field lattice model

We develop a combined model to describe the pressure-volume-temperature (PVT) thermodynamics and the α- and β-relaxation time dynamics in glass-forming amorphous materials. The PVT results are described using a two-state modification of the Sanchez-Lacombe equation of state (SL-EoS). The minimization of the Sanchez-Lacombe free energy expression allows one to calculate the density (or specific volume) and the "solid fraction" as a function of temperature and pressure. The solid fraction, ψ(T,P), is then substituted into the TS2 mean-field model (V. V. Ginzburg, Soft Matter, 2020, 16, 810), and the equilibrium relaxation times and the glass transition temperature are calculated. Finally, we use a dynamic model of a Tool-Narayanaswami-Moynihan (TNM)-type to describe the density relaxation as a function of the cooling rate. We applied the new framework to describe experimental data for polystyrene and poly(methylmethacrylate) and found a good qualitative and quantitative agreement.

LCST polymers with UCST behavior

In this study, temperature dependent behavior of dense dispersions of core crosslinked flower-like micelles is investigated. Micelles were prepared by mixing aqueous solutions of two ABA block copolymers with PEG B-blocks and thermosensitive A-blocks containing PNIPAM and crosslinkable moieties. At a temperature above the lower critical solution temperature (LCST), self-assembly of the polymers resulted in the formation of flower-like micelles with a hydrophilic PEG shell and a hydrophobic core. The micellar core was stabilized by native chemical ligation (NCL). Above the LCST, micelles displayed a radius of ∼35 nm, while a radius of ∼48 nm was found below the LCST due to hydration of the PNIPAM core. Concentrated dispersions of these micelles (≥7.5 wt%) showed glassy state behavior below a critical temperature (Tc: 28 °C) which is close to the LCST of the polymers. Below this Tc, the increase in the micelle volume resulted in compression of micelles together above a certain concentration and formation of a glass. We quantified and compared micelle packing at different concentrations and temperatures. The storage moduli (G') of the dispersions showed a universal dependence on the effective volume fraction, which increased substantially above a certain effective volume fraction of φ = 1.2. Furthermore, a disordered lattice model describing this behavior fitted the experimental data and revealed a critical volume fraction of φc = 1.31 close to the experimental value of φ = 1.2. The findings reported provide insights for the molecular design of novel thermosensitive PNIPAM nanoparticles with tunable structural and mechanical properties.

Determination of Young's Modulus of Active Pharmaceutical Ingredients by Relaxation Dynamics at Elevated Pressures

A new approach is theoretically proposed to study the glass transition of active pharmaceutical ingredients and a glass-forming anisotropic molecular liquid at high pressures. We describe amorphous materials as a fluid of hard spheres. Effects of nearest neighbor interactions and cooperative motions of particles on glassy dynamics are quantified through a local and collective elastic barrier calculated using the elastically collective nonlinear Langevin equation theory. Inserting two barriers into Kramer's theory gives the structural relaxation time. Then, we formulate a new mapping based on the thermal expansion process under pressure to intercorrelate particle density, temperature, and pressure. This analysis allows us to determine the pressure and temperature dependence of α relaxation. From this, we estimate the effective elastic modulus of amorphous materials and capture the effects of conformation on the relaxation process. Remarkably, our theoretical results agree well with experiments.

Microscopic ergodicity breaking governs the emergence and evolution of elasticity in glass-forming nanoclay suspensions

We report a study combining x-ray photon correlation spectroscopy (XPCS) with in situ rheology to investigate the microscopic dynamics and mechanical properties of aqueous suspensions of the synthetic hectorite clay Laponite, which is composed of charged, nanometer-scale, disk-shaped particles. The suspensions, with particle concentrations ranging from 3.25 to 3.75 wt %, evolve over time from a fluid to a soft glass that displays aging behavior. The XPCS measurements characterize the localization of the particles during the formation and aging of the soft-glass state. The fraction of localized particles, f_{0}, increases rapidly during the early formation stage and grows more slowly during subsequent aging, while the characteristic localization length r_{loc} steadily decreases. Despite the strongly varying rates of aging at different concentrations, both f_{0} and r_{loc} scale with the elastic shear modulus G^{'} in a manner independent of concentration. During the later aging stage, the scaling between r_{loc} and G^{'} agrees quantitatively with a prediction of naive mode coupling theory. Breakdown of agreement with the theory during the early formation stage indicates the prevalence of dynamic heterogeneity, suggesting the soft solid forms through precursors of dynamically localized clusters.

Analytical prediction of logarithmic Rayleigh scattering in amorphous solids from tensorial heterogeneous elasticity with power-law disorder

The damping or attenuation coefficient of sound waves in solids due to impurities scales with the wavevector to the fourth power, also known as Rayleigh scattering. In amorphous solids, Rayleigh scattering may be enhanced by a logarithmic factor although computer simulations offer conflicting conclusions regarding this enhancement and its microscopic origin. We present a tensorial replica field-theoretic derivation based on heterogeneous or fluctuating elasticity (HE), which shows that long-range (power-law) spatial correlations of the elastic constants, is the origin of the logarithmic enhancement to Rayleigh scattering of phonons in amorphous solids. We also consider the case of zero spatial fluctuations in the elastic constants, and of power-law decaying fluctuations in the internal stresses. Also in this case the logarithmic enhancement to the Rayleigh scattering law can be derived from the proposed tensorial HE framework.

Growth and arrest of topological cycles in small physical networks

The chordless cycle sizes of spatially embedded networks are demonstrated to follow an exponential growth law similar to random graphs if the number of nodes [Formula: see text] is below a critical value [Formula: see text] For covalent polymer networks, increasing the network size, as measured by the number of cross-link nodes, beyond [Formula: see text] results in a crossover to a new regime in which the characteristic size of the chordless cycles [Formula: see text] no longer increases. From this result, the onset and intensity of finite-size effects can be predicted from measurement of [Formula: see text] in large networks. Although such information is largely inaccessible with experiments, the agreement of simulation results from molecular dynamics, Metropolis Monte Carlo, and kinetic Monte Carlo suggests the crossover is a fundamental physical feature which is insensitive to the details of the network generation. These results show random graphs as a promising model to capture structural differences in confined physical networks.

Rheology of hard glassy materials

Glassy solids may undergo a fluidization (yielding) transition upon deformation whereby the material starts to flow plastically. It has been a matter of debate whether this process is controlled by a specific time scale, from among different competing relaxation/kinetic processes. Here, two constitutive models of cage relaxation are examined within the microscopic model of nonaffine elasto-plasticity. One (widely used) constitutive model implies that the overall relaxation rate is dominated by the fastest between the structural (α) relaxation rate and the shear-induced relaxation rate. A different model is formulated here which, instead, assumes that the slowest (global) relaxation process controls the overall relaxation. We show that the first model is not compatible with the existence of finite elastic shear modulus for quasistatic (low-frequency) deformation, while the second model is able to describe all key features of deformation of 'hard' glassy solids, including the yielding transition, the nonaffine-to-affine plateau crossover, and the rate-stiffening of the modulus. The proposed framework provides an operational way to distinguish between 'soft' glasses and 'hard' glasses based on the shear-rate dependence of the structural relaxation time.

Electroplasticization of Liquid Crystal Polymer Networks

Shape-shifting liquid crystal networks (LCNs) can transform their morphology and properties in response to external stimuli. These active and adaptive polymer materials can have impact in a diversity of fields, including haptic displays, energy harvesting, biomedicine, and soft robotics. Electrically driven transformations in LCN coatings are particularly promising for application in electronic devices, in which electrodes are easily integrated and allow for patterning of the functional response. The morphing of these coatings, which are glassy in the absence of an electric field, relies on a complex interplay between polymer viscoelasticity, liquid crystal order, and electric field properties. Morphological transformations require the material to undergo a glass transition that plasticizes the polymer sufficiently to enable volumetric and shape changes. Understanding how an alternating current can plasticize very stiff, densely cross-linked networks remains an unresolved challenge. Here, we use a nanoscale strain detection method to elucidate this electric-field-induced devitrification of LCNs. We find how a high-frequency alternating field gives rise to pronounced nanomechanical changes at a critical frequency, which signals the electrical glass transition. Across this transition, collective motion of the liquid crystal molecules causes the network to yield from within, leading to network weakening and subsequent nonlinear expansion. These results unambiguously prove the existence of electroplasticization. Fine-tuning the induced emergence of plasticity will not only enhance the surface functionality but also enable more efficient conversion of electrical energy into mechanical work.

Boson peak, elasticity, and glass transition temperature in polymer glasses: Effects of the rigidity of chain bending

The excess low-frequency vibrational spectrum, called boson peak, and non-affine elastic response are the most important particularities of glasses. Herein, the vibrational and mechanical properties of polymeric glasses are examined by using coarse-grained molecular dynamics simulations, with particular attention to the effects of the bending rigidity of the polymer chains. As the rigidity increases, the system undergoes a glass transition at a higher temperature (under a constant pressure), which decreases the density of the glass phase. The elastic moduli, which are controlled by the decrease of the density and the increase of the rigidity, show a non-monotonic dependence on the rigidity of the polymer chain that arises from the non-affine component. Moreover, a clear boson peak is observed in the vibrational density of states, which depends on the macroscopic shear modulus G. In particular, the boson peak frequency ω_BP is proportional to \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{G}$$\end{document}G. These results provide a positive correlation between the boson peak, shear elasticity, and the glass transition temperature.

Theory of applying shear strains from boundary walls: Linear response in glasses

We construct a linear response theory of applying shear deformations from boundary walls in the film geometry in Kubo's theoretical scheme. Our method is applicable to any solids and fluids. For glasses, we assume quasi-equilibrium around a fixed inherent state. Then, we obtain linear-response expressions for any variables including the stress and the particle displacements, even though the glass interior is elastically inhomogeneous. In particular, the shear modulus can be expressed in terms of the correlations between the interior stress and the forces from the walls. It can also be expressed in terms of the inter-particle correlations, as has been shown in the previous literature. Our stress relaxation function includes the effect of the boundary walls and can be used for inhomogeneous flow response. We show the presence of long-ranged, long-lived correlations among the fluctuations of the forces from the walls and the displacements of all the particles in the cell. We confirm these theoretical results numerically in a two-dimensional model glass. As an application, we describe emission and propagation of transverse sounds after boundary wall motions using these time-correlation functions. We also find resonant sound amplification when the frequency of an oscillatory shear approaches that of the first transverse sound mode.

Parameter-free predictions of the viscoelastic response of glassy polymers from non-affine lattice dynamics

We study the viscoelastic response of amorphous polymers using theory and simulations. By accounting for internal stresses and considering instantaneous normal modes (INMs) within athermal non-affine theory, we make parameter-free predictions of the dynamic viscoelastic moduli obtained in coarse-grained simulations of polymer glasses at non-zero temperatures. The theoretical results show very good correspondence with rheology data collected from molecular dynamics simulations over five orders of magnitude in frequency, with some instabilities that accumulate in the low-frequency part on approach to the glass transition. These results provide evidence that the mechanical glass transition itself is continuous and thus represents a crossover rather than a true phase transition. The relatively sharp drop of the low-frequency storage modulus across the glass transition temperature can be explained mechanistically within the proposed theory: the proliferation of low-eigenfrequency vibrational excitations (boson peak and nearly-zero energy excitations) is directly responsible for the rapid growth of a negative non-affine contribution to the storage modulus.

Linking slow dynamics and microscopic connectivity in dense suspensions of charged colloids

The quest to unravel the nature of the glass transition, where the viscosity of a liquid increases by many orders of magnitude, while its static structure remains largely unaffected, remains unresolved. While various structural and dynamical precursors to vitrification have been identified, a predictive and quantitative description of how subtle changes at the microscopic scale give rise to the steep growth in macroscopic viscosity is missing. It was recently proposed that the presence of long-lived bonded structures within the liquid may provide the long-sought connection between local structure and global dynamics. Here we directly observe and quantify the connectivity dynamics in liquids of charged colloids en route to vitrification using three-dimensional confocal microscopy. We determine the dynamic structure from the real-space van Hove correlation function and from the particle trajectories, providing upper and lower bounds on connectivity dynamics. Based on these data, we extend Dyre's model for the glass transition to account for particle-level structural dynamics; this results in a microscopic expression for the slowing down of relaxations in the liquid that is in quantitative agreement with our experiments. These results indicate how vitrification may be understood as a dynamical connectivity transition with features that are strongly reminiscent of rigidity percolation scenarios.

Nonmonotonic dependence of polymer-glass mechanical response on chain bending stiffness

We investigate the mechanical properties of amorphous polymers by means of coarse-grained simulations and nonaffine lattice dynamics theory. A small increase of polymer chain bending stiffness leads first to softening of the material, while hardening happens only upon further strengthening of the backbones. This nonmonotonic variation of the storage modulus G^{'} with bending stiffness is caused by a competition between additional resistance to deformation offered by stiffer backbones and decreased density of the material due to a necessary decrease in monomer-monomer coordination. This counterintuitive finding suggests that the strength of polymer glasses may in some circumstances be enhanced by softening the bending of constituent chains.

Shear-stress fluctuations and relaxation in polymer glasses

We investigate by means of molecular dynamics simulation a coarse-grained polymer glass model focusing on (quasistatic and dynamical) shear-stress fluctuations as a function of temperature T and sampling time Δt. The linear response is characterized using (ensemble-averaged) expectation values of the contributions (time averaged for each shear plane) to the stress-fluctuation relation μ_{sf} for the shear modulus and the shear-stress relaxation modulus G(t). Using 100 independent configurations, we pay attention to the respective standard deviations. While the ensemble-averaged modulus μ_{sf}(T) decreases continuously with increasing T for all Δt sampled, its standard deviation δμ_{sf}(T) is nonmonotonic with a striking peak at the glass transition. The question of whether the shear modulus is continuous or has a jump singularity at the glass transition is thus ill posed. Confirming the effective time-translational invariance of our systems, the Δt dependence of μ_{sf} and related quantities can be understood using a weighted integral over G(t).

Isostaticity and the solidification of semiflexible polymer melts

Using molecular dynamics simulations of a tangent-soft-sphere bead-spring polymer model, we examine the degree to which semiflexible polymer melts solidify at isostaticity. Flexible and stiff chains crystallize when they are isostatic as defined by appropriate degree-of-freedom-counting arguments. Semiflexible chains also solidify when isostatic if a generalized isostaticity criterion that accounts for the slow freezing out of configurational freedom as chain stiffness increases is employed. The configurational freedom associated with bond angles (θ) can be associated with the characteristic ratio C_∞ = (1 + 〈cos(θ)〉)/(1 - 〈cos(θ)〉). We find that the dependence of the average coordination number at solidification [Z(T_s)] on chains' characteristic ratio C_∞ has the same functional form [Z ≃ a - b ln(C_∞)] as the dependence of the average coordination number at jamming [Z(ϕ_J)] on C_∞ in athermal systems, suggesting that jamming-related phenomena play a significant role in thermal polymer solidification.

Shear Modulus and Shear-Stress Fluctuations in Polymer Glasses

Using molecular dynamics simulation of a standard coarse-grained polymer glass model, we investigate by means of the stress-fluctuation formalism the shear modulus μ as a function of temperature T and sampling time Δt. While the ensemble-averaged modulus μ(T) is found to decrease continuously for all Δt sampled, its standard deviation δμ(T) is nonmonotonic, with a striking peak at the glass transition. Confirming the effective time-translational invariance of our systems, μ(Δt) can be understood using a weighted integral over the shear-stress relaxation modulus G(t). While the crossover of μ(T) gets sharper with an increasing Δt, the peak of δμ(T) becomes more singular. It is thus elusive to predict the modulus of a single configuration at the glass transition.

Thermodynamic scaling of relaxation: insights from anharmonic elasticity

Using molecular dynamics simulations of a molecular liquid, we investigate the thermodynamic scaling (TS) of the structural relaxation time [Formula: see text] in terms of the quantity [Formula: see text], where T and ρ are the temperature and density, respectively. The liquid does not exhibit strong virial-energy correlations. We propose a method for evaluating both the characteristic exponent [Formula: see text] and the TS master curve that uses experimentally accessible quantities that characterise the anharmonic elasticity and does not use details about the microscopic interactions. In particular, we express the TS characteristic exponent [Formula: see text] in terms of the lattice Grüneisen parameter [Formula: see text] and the isochoric anharmonicity [Formula: see text]. An analytic expression of the TS master curve of [Formula: see text] with [Formula: see text] as the key adjustable parameter is found. The comparison with the experimental TS master curves and the isochoric fragilities of 34 glassformers is satisfying. In a few cases, where thermodynamic data are available, we test (i) the predicted characteristic exponent [Formula: see text] and (ii) the isochoric anharmonicity [Formula: see text], as drawn by the best fit of the TS of the structural relaxation, against the available thermodynamic data. A linear relation between the isochoric fragility and the isochoric anharmonicity [Formula: see text] is found and compared favourably with the results of experiments with no adjustable parameters. A relation between the increase of the isochoric vibrational heat capacity due to anharmonicity and the isochoric fragility is derived.

Direct Observation of Entropic Stabilization of bcc Crystals Near Melting

Crystals with low latent heat are predicted to melt from an entropically stabilized body-centered cubic symmetry. At this weakly first-order transition, strongly correlated fluctuations are expected to emerge, which could change the nature of the transition. Here we show how large fluctuations stabilize bcc crystals formed from charged colloids, giving rise to strongly power-law correlated heterogeneous dynamics. Moreover, we find that significant nonaffine particle displacements lead to a vanishing of the nonaffine shear modulus at the transition. We interpret these observations by reformulating the Born-Huang theory to account for nonaffinity, illustrating a scenario of ordered solids reaching a state where classical lattice dynamics fail.

Jamming of Semiflexible Polymers

We study jamming in model freely rotating polymers as a function of chain length N and bond angle θ_{0}. The volume fraction at jamming ϕ_{J}(θ_{0}) is minimal for rigid-rodlike chains (θ_{0}=0), and increases monotonically with increasing θ_{0}≤π/2. In contrast to flexible polymers, marginally jammed states of freely rotating polymers are highly hypostatic, even when bond and angle constraints are accounted for. Large-aspect-ratio (small θ_{0}) chains behave comparably to stiff fibers: resistance to large-scale bending plays a major role in their jamming phenomenology. Low-aspect-ratio (large θ_{0}) chains behave more like flexible polymers, but still jam at much lower densities due to the presence of frozen-in three-body correlations corresponding to the fixed bond angles. Long-chain systems jam at lower ϕ and are more hypostatic at jamming than short-chain systems. Implications of these findings for polymer solidification are discussed.

Direct link between boson-peak modes and dielectric α-relaxation in glasses

We compute the dielectric response of glasses starting from a microscopic system-bath Hamiltonian of the Zwanzig-Caldeira-Leggett type and using an ansatz from kinetic theory for the memory function in the resulting generalized Langevin equation. The resulting framework requires the knowledge of the vibrational density of states (DOS) as input, which we take from numerical evaluation of a marginally stable harmonic disordered lattice, featuring a strong boson peak (excess of soft modes over Debye ∼ω_{p}^{2} law). The dielectric function calculated based on this ansatz is compared with experimental data for the paradigmatic case of glycerol at T≲T_{g}. Good agreement is found for both the reactive (real) part of the response and for the α-relaxation peak in the imaginary part, with a significant improvement over earlier theoretical approaches. On the low-frequency side of the α peak, the fitting supports the presence of ∼ω_{p}^{4} modes at vanishing eigenfrequency as recently shown [E. Lerner, G. During, and E. Bouchbinder, Phys. Rev. Lett. 117, 035501 (2016)PRLTAO0031-900710.1103/PhysRevLett.117.035501]. α-wing asymmetry and stretched-exponential behavior are recovered by our framework, which shows that these features are, to a large extent, caused by the soft boson-peak modes in the DOS.

Polymer glass transition occurs at the marginal rigidity point with connectivity z* = 4

We re-examine the physical origin of the polymer glass transition from the point of view of marginal rigidity, which is achieved at a certain average number of mechanically active intermolecular contacts per monomer. In the case of polymer chains in a melt/poor solvent, each monomer has two neighbors bound by covalent bonds and also a number of central-force contacts modelled by the Lennard-Jones (LJ) potential. We find that when the average number of contacts per monomer (covalent and non-covalent) exceeds the critical value z* ≈ 4, the system becomes solid and the dynamics arrested - a state that we declare the glass. Coarse-grained Brownian dynamics simulations show that at sufficient strength of LJ attraction (which effectively represents the depth of quenching, or the quality of solvent) the polymer globule indeed crosses the threshold of z*, and becomes a glass with a finite zero-frequency shear modulus, G∝ (z-z*). We verify this by showing the distinction between the 'liquid' polymer droplet at z < z*, which changes shape and adopts the spherical conformation in equilibrium, and the glassy 'solid' droplet at z > z*, which retains its shape frozen at the moment of z* crossover. These results provide a robust microscopic criterion to tell the liquid apart from the glass for the linear polymers.

Elastic moduli and vibrational modes in jammed particulate packings

When we elastically impose a homogeneous, affine deformation on amorphous solids, they also undergo an inhomogeneous, nonaffine deformation, which can have a crucial impact on the overall elastic response. To correctly understand the elastic modulus M, it is therefore necessary to take into account not only the affine modulus M_{A}, but also the nonaffine modulus M_{N} that arises from the nonaffine deformation. In the present work, we study the bulk (M=K) and shear (M=G) moduli in static jammed particulate packings over a range of packing fractions φ. The affine M_{A} is determined essentially by the static structural arrangement of particles, whereas the nonaffine M_{N} is related to the vibrational eigenmodes. We elucidate the contribution of each vibrational mode to the nonaffine M_{N} through a modal decomposition of the displacement and force fields. In the vicinity of the (un)jamming transition φ_{c}, the vibrational density of states g(ω) shows a plateau in the intermediate-frequency regime above a characteristic frequency ω^{*}. We illustrate that this unusual feature apparent in g(ω) is reflected in the behavior of M_{N}: As φ→φ_{c}, where ω^{*}→0, those modes for ω<ω^{*} contribute less and less, while contributions from those for ω>ω^{*} approach a constant value which results in M_{N} to approach a critical value M_{Nc}, as M_{N}-M_{Nc}∼ω^{*}. At φ_{c} itself, the bulk modulus attains a finite value K_{c}=K_{Ac}-K_{Nc}>0, such that K_{Nc} has a value that remains below K_{Ac}. In contrast, for the critical shear modulus G_{c}, G_{Nc} and G_{Ac} approach the same value so that the total value becomes exactly zero, G_{c}=G_{Ac}-G_{Nc}=0. We explore what features of the configurational and vibrational properties cause such a distinction between K and G, allowing us to validate analytical expressions for their critical values.

Importance of non-affine viscoelastic response in disordered fibre networks

Disordered fibre networks are ubiquitous in nature and have a wide range of industrial applications as novel biomaterials. Predicting their viscoelastic response is straightforward for affine deformations that are uniform over all length scales, but when affinity fails, as has been observed experimentally, modelling becomes challenging. Here we present a numerical methodology, related to an existing framework for amorphous packings, to predict the steady-state viscoelastic spectra and degree of affinity for disordered fibre networks driven at arbitrary frequencies. Applying this method to a peptide gel model reveals a monotonic increase of the shear modulus as the soft, non-affine normal modes are successively suppressed as the driving frequency increases. In addition to being dominated by fibril bending, these low frequency network modes are also shown to be delocalised. The presented methodology provides insights into the importance of non-affinity in the viscoelastic response of peptide gels, and is easily extendible to all types of fibre networks.

Reversibility and hysteresis of the sharp yielding transition of a colloidal glass under oscillatory shear

The mechanical response of glasses remains challenging to understand. Recent results indicate that the oscillatory rheology of soft glasses is accompanied by a sharp non-equilibrium transition in the microscopic dynamics. Here, we use simultaneous x-ray scattering and rheology to investigate the reversibility and hysteresis of the sharp symmetry change from anisotropic solid to isotropic liquid dynamics observed in the oscillatory shear of colloidal glasses (D. Denisov, M.T. Dang, B. Struth, A. Zaccone, P. Schall, Sci. Rep. 5 14359 (2015)). We use strain sweeps with increasing and decreasing strain amplitude to show that, in analogy with equilibrium transitions, this sharp symmetry change is reversible and exhibits systematic frequency-dependent hysteresis. Using the non-affine response formalism of amorphous solids, we show that these hysteresis effects arise from frequency-dependent non-affine structural cage rearrangements at large strain. These results consolidate the first-order-like nature of the oscillatory shear transition and quantify related hysteresis effects both via measurements and theoretical modelling.

Rigidity in Condensed Matter and Its Origin in Configurational Constraint

Motivated by the formal argument that a nonzero shear modulus is the result of averaging over a constrained configuration space, we demonstrate that the shear modulus calculated over a range of temperatures and averaging times can be expressed (relative to its infinite frequency value) as a single function of the mean squared displacement. This result is shown to hold for both a glass-liquid and a crystal-liquid system.

Anomalous energy cascades in dense granular materials yielding under simple shear deformations

By using molecular dynamics (MD) simulations of dense granular particles in two dimensions, we study turbulent-like structures of their non-affine velocities under simple shear deformations. We find that the spectrum of non-affine velocities, introduced as an analog of the energy spectrum for turbulent flows, exhibits the power-law decay if the system is yielding in a quasi-static regime, where large-scale collective motions and inelastic interactions of granular particles are crucial for the anomalous cascade of kinetic energy. Based on hydrodynamic equations of dense granular materials, which include both kinetic and contact contributions in constitutive relations, we derive a theoretical expression for the spectrum, where a good agreement between the result of MD simulations and theoretical prediction is established over a wide range of length scales.

Glass transition of two-dimensional 80-20 Kob-Andersen model at constant pressure

We reconsider numerically the two-dimensional version of the Kob-Andersen model (KA2d) with a fraction of 80% of large spheres. A constant moderate pressure is imposed while the temperature T is systematically quenched from the liquid limit through the glass transition at [Formula: see text] down to very low temperatures. Monodisperse Lennard-Jones (mdLJ) bead systems, forming a crystal phase at low temperatures, are used to highlight several features of the KA2d model. As can be seen, e.g. from the elastic shear modulus G(T), determined using the stress-fluctuation formalism, our KA2d model is a good glass-former. A continuous cusp-singularity, [Formula: see text] with [Formula: see text], is observed in qualitative agreement with other recent numerical and theoretical work, however in striking conflict with the additive jump discontinuity predicted by mode-coupling theory.

Influence of molecular-weight polydispersity on the glass transition of polymers

It is well known that the polymer glass transition temperature T_{g} is dependent on molecular weight, but the role of molecular-weight polydispersity on T_{g} is unclear. Using molecular-dynamics simulations, we clarify that for polymers with the same number-average molecular weight, the molecular-weight distribution profile (either in Schulz-Zimm form or in bimodal form) has very little influence on the glass transition temperature T_{g}, the average segment dynamics (monomer motion, bond orientation relaxation, and torsion transition), and the relaxation-time spectrum, which are related to the local nature of the glass transition. By analyzing monomer motions in different chains, we find that the motion distribution of monomers is altered by molecular-weight polydispersity. Molecular-weight polydispersity dramatically enhances the dynamic heterogeneity of monomer diffusive motions after breaking out of the "cage," but it has a weak influence on the dynamic heterogeneity of the short time scales and the transient spatial correlation between temporarily localized monomers. The stringlike cooperative motion is also not influenced by molecular-weight polydispersity, supporting the idea that stringlike collective motion is not strongly correlated with chain connectivity.

Interatomic repulsion softness directly controls the fragility of supercooled metallic melts

We present an analytic scheme to connect the fragility and viscoelasticity of metallic glasses to the effective ion-ion interaction in the metal. This is achieved by an approximation of the short-range repulsive part of the interaction, combined with nonaffine lattice dynamics to obtain analytical expressions for the shear modulus, viscosity, and fragility in terms of the ion-ion interaction. By fitting the theoretical model to experimental data, we are able to link the steepness of the interionic repulsion to the Thomas-Fermi screened Coulomb repulsion and to the Born-Mayer valence electron overlap repulsion for various alloys. The result is a simple closed-form expression for the fragility of the supercooled liquid metal in terms of few crucial atomic-scale interaction and anharmonicity parameters. In particular, a linear relationship is found between the fragility and the energy scales of both the screened Coulomb and the electron overlap repulsions. This relationship opens up opportunities to fabricate alloys with tailored thermoelasticity and fragility by rationally tuning the chemical composition of the alloy according to general principles. The analysis presented here brings a new way of looking at the link between the outer shell electronic structure of metals and metalloids and the viscoelasticity and fragility thereof.