Fast generation of ultrastable computer glasses by minimization of an augmented potential energy

Irreversible Monte Carlo Algorithms for Hard Disk Glasses: From Event-Chain to Collective Swaps

Equilibrium sampling of the configuration space in disordered systems requires algorithms that bypass the glassy slowing down of the physical dynamics. Irreversible Monte Carlo algorithms breaking detailed balance successfully accelerate sampling in some systems. We first implement an irreversible event-chain Monte Carlo algorithm in a model of continuously polydisperse hard disks. The effect of collective translational moves marginally affects the dynamics and results in a modest speedup that decreases with density. We then propose an irreversible algorithm performing collective particle swaps which outperforms all known Monte Carlo algorithms. We show that these collective swaps can also be used to prepare very dense jammed packings of disks.

Elasticity of self-organized frustrated disordered spring networks

There have been some interesting recent advances in understanding the notion of mechanical disorder in structural glasses and the statistical mechanics of these systems' low-energy excitations. Here we contribute to these advances by studying a minimal model for structural glasses' elasticity in which the degree of mechanical disorder-as characterized by recently introduced dimensionless quantifiers-is readily tunable over a very large range. We comprehensively investigate a number of scaling laws observed for various macro, meso and microscopic elastic properties, and rationalize them using scaling arguments. Interestingly, we demonstrate that the model features the universal quartic glassy vibrational density of states as seen in many atomistic and molecular models of structural glasses formed by cooling a melt. The emergence of this universal glassy spectrum highlights the role of self-organization (toward mechanical equilibrium) in its formation, and elucidates why models featuring structural frustration alone do not feature the same universal glassy spectrum. Finally, we discuss relations to existing work in the context of strain stiffening of elastic networks and of low-energy excitations in structural glasses, in addition to future research directions.

Effects of polydispersity on the plastic behaviors of dense two-dimensional granular systems under shear

We study particle-scale motion in sheared highly polydisperse amorphous materials, in which the largest particles are as much as ten times the size of the smallest. We find strikingly different behavior from the more commonly studied amorphous systems with low polydispersity. In particular, an analysis of the nonaffine motion of particles reveals qualitative differences between large and small particles: The smaller particles have dramatically more nonaffine motion, which is induced by the presence of the large particles. We characterize how the nonaffine motion changes from the low- to high-polydispersity regimes. We further demonstrate a quantitative way to distinguish between "large" and "small" particles in systems with broad distributions of particle sizes. A macroscopic consequence of the nonaffine motion is a decrease in the energy dissipation rate for highly polydisperse samples, which is due both to a geometric consequence of the changing jamming conditions for higher polydispersity and to the changing character of nonaffine motion.

Creating bulk ultrastable glasses by random particle bonding

A recent breakthrough in glass science has been the synthesis of ultrastable glasses via physical vapor deposition techniques. These samples display enhanced thermodynamic, kinetic and mechanical stability, with important implications for fundamental science and technological applications. However, the vapor deposition technique is limited to atomic, polymer and organic glass-formers and is only able to produce thin film samples. Here, we propose a novel approach to generate ultrastable glassy configurations in the bulk, via random particle bonding, and using computer simulations we show that this method does indeed allow for the production of ultrastable glasses. Our technique is in principle applicable to any molecular or soft matter system, such as colloidal particles with tunable bonding interactions, thus opening the way to the design of a large class of ultrastable glasses.

Active responsive colloids driven by intrinsic dichotomous noise

We study the influence of intrinsic noise on the structure and dynamics of responsive colloids (RCs), which actively change their size and mutual interactions. The colloidal size is explicitly resolved in our RC model as an internal degree of freedom (DOF) in addition to the particle translation. A Hertzian pair potential between the RCs leads to repulsion and shrinking of the particles, resulting in an explicit responsiveness of the system to self-crowding. To render the colloids active, their size is internally driven by a dichotomous noise, randomly switching ("breathing") between growing and shrinking states with a predefined rate, as motivated by recent experiments on synthetic active colloids. The polydispersity of this dichotomous active responsive colloid (D-ARC) model can be tuned by the parameters of the noise. Utilizing stochastic computer simulations, we study crowding effects on the spatial distributions, relaxation times, and self-diffusion of dense suspensions of the D-ARCs. We find a substantial influence of the "built-in" intrinsic noise on the system's behavior, in particular, transitions from unimodal to bimodal size distributions for an increasing colloid density as well as intrinsic noise-modified diffusive translational dynamics. We conclude that controlling the noise of internal DOFs of a macromolecule or cell is a powerful tool for active colloidal materials to enable autonomous changes in the system's collective structure and dynamics towards the adaptation of macroscopic properties to external perturbations.

Predicting Defects in Soft Sphere Packings near Jamming Using the Force Network Ensemble

Amorphous systems of soft particles above jamming have more contacts than are needed to achieve mechanical equilibrium. The force network of a granular system with a fixed contact network is thus underdetermined and can be characterized as a random instantiation within the space of the force network ensemble. In this Letter, we show that defect contacts that are not necessary for stability of the system can be uniquely identified by examining the boundaries of this space of allowed force networks. We further show that, for simulations in the near jamming limit, this identification is nearly always correct and that defect contacts are broken under decompression of the system.

Structural Measures as Guides to Ultrastable States in Overjammed Packings

Jammed, disordered packings of given sets of particles possess a multitude of equilibrium states with different mechanical properties. Identifying and constructing desired states, e.g., of superior stability, is a complex task. Here, we show that in two-dimensional particle packings the energy of all metastable states (inherent structures) is reliably classified by simple scalar measures of local steric packing. These structural measures are insensitive to the particle interaction potential and so robust that they can be used to guide a modified swap algorithm that anneals polydisperse packings toward low-energy metastable states exceptionally fast. The low-energy states are extraordinarily stable against applied shear, so that the approach also efficiently identifies ultrastable packings.

Active responsive colloids coupled to different thermostats

We introduce a model of active responsive colloids (ARCs) in which an internal degree of freedom (DoF) of a single colloidal particle is "activated" by coupling it to a different thermostat than for the translational DoFs. As for the responsive internal DoF, we consider specifically the size (diameter) of the spherical particles, which is confined by a harmonic parent potential being either entropic or energetic in nature. The ARCs interact via a repulsive Hertzian pair potential, appropriate to model hydrogels or elastic colloids, and are studied for various densities using Brownian dynamics simulations. We tune the internal activity in the nonequilibrium steady state by scanning through a wide range of internal temperatures, both smaller ("colder") and larger ("hotter") than the translational temperature. The results show a rich and intriguing behavior for the emergent property distributions, colloidal pair structure, and the diffusive translational dynamics controlled by the internal activity, substantially depending on whether the internal DoF moves in an entropic or energetic potential. We discuss theoretical thermal limits and phenomenological models which can explain some of the nonequilibrium trends qualitatively. Our study indicates that particle dynamical polydispersity as well as the structure and dynamics of dense macromolecular suspensions can be vastly tuned by internal activity in terms of internal "hot" or "cold" fluctuating states.

Transient learning degrees of freedom for introducing function in materials

SignificanceMany protocols used in material design and training have a common theme: they introduce new degrees of freedom, often by relaxing away existing constraints, and then evolve these degrees of freedom based on a rule that leads the material to a desired state at which point these new degrees of freedom are frozen out. By creating a unifying framework for these protocols, we can now understand that some protocols work better than others because the choice of new degrees of freedom matters. For instance, introducing particle sizes as degrees of freedom to the minimization of a jammed particle packing can lead to a highly stable state, whereas particle stiffnesses do not have nearly the same impact.

Mean-field description for the architecture of low-energy excitations in glasses

In amorphous materials, groups of particles can rearrange locally into a new stable configuration. Such elementary excitations are key as they determine the response to external stresses, as well as to thermal and quantum fluctuations. Yet, understanding what controls their geometry remains a challenge. Here we build a scaling description of the geometry and energy of low-energy excitations in terms of the distance to an instability, as predicted, for instance, at the dynamical transition in mean-field approaches of supercooled liquids. We successfully test our predictions in ultrastable computer glasses, with a gapped spectrum and an ungapped (regular) spectrum. Overall, our approach explains why excitations become less extended, with a higher energy and displacement scale upon cooling.

Modulating internal transition kinetics of responsive macromolecules by collective crowding

Packing and crowding are used in biology as mechanisms to (self-)regulate internal molecular or cellular processes based on collective signaling. Here, we study how the transition kinetics of an internal "switch" of responsive macromolecules is modified collectively by their spatial packing. We employ Brownian dynamics simulations of a model of Responsive Colloids, in which an explicit internal degree of freedom-here, the particle size-moving in a bimodal energy landscape self-consistently responds to the density fluctuations of the crowded environment. We demonstrate that populations and transition times for the two-state switching kinetics can be tuned over one order of magnitude by "self-crowding." An exponential scaling law derived from a combination of Kramers' and liquid state perturbation theory is in very good agreement with the simulations.

Marginal stability in memory training of jammed solids

Memory encoding by cyclic shear is a reliable process to store information in jammed solids, yet its underlying mechanism and its connection to the amorphous structure are not fully understood. When a jammed sphere packing is repeatedly sheared with cycles of the same strain amplitude, it optimizes its mechanical response to the cyclic driving and stores a memory of it. We study memory by cyclic shear training as a function of the underlying stability of the amorphous structure in marginally stable and highly stable packings, the latter produced by minimizing the potential energy using both positional and radial degrees of freedom. We find that jammed solids need to be marginally stable in order to store a memory by cyclic shear. In particular, highly stable packings store memories only after overcoming brittle yielding and the cyclic shear training takes place in the shear band, a region which we show to be marginally stable.

Disorder perturbation expansion for athermal crystals

We introduce a perturbation expansion for athermal systems that allows an exact determination of displacement fields away from the crystalline state as a response to disorder. We show that the displacement fields in energy-minimized configurations of particles interacting through central potentials with microscopic disorder can be obtained as a series expansion in the strength of the disorder. We introduce a hierarchy of force-balance equations that allows an order-by-order determination of the displacement fields, with the solutions at lower orders providing sources for the higher-order solutions. This allows the simultaneous force-balance equations to be solved, within a hierarchical perturbation expansion to arbitrary accuracy. We present exact results for an isotropic defect introduced into the crystalline ground state at linear order and second order in our expansion. We show that the displacement fields produced by the defect display interesting self-similar properties at every order. We derive a |δr|∼1/r and |δf|∼1/r^{2} decay for the displacement fields and excess interparticle forces at large distances r away from the defect. Finally, we derive nonlinear corrections introduced by the interactions between defects at second order in our expansion. We verify our exact results with displacement fields obtained from energy-minimized configurations of soft disks.

Soft Particles at Liquid Interfaces: From Molecular Particle Architecture to Collective Phase Behavior

Soft particles such as microgels can undergo significant and anisotropic deformations when adsorbed to a liquid interface. This, in turn, leads to a complex phase behavior upon compression. To date, experimental efforts have predominantly provided phenomenological links between microgel structure and resulting interfacial behavior, while simulations have not been entirely successful in reproducing experiments or predicting the minimal requirements for the desired phase behavior. Here, we develop a multiscale framework to link the molecular particle architecture to the resulting interfacial morphology and, ultimately, to the collective interfacial phase behavior. To this end, we investigate interfacial morphologies of different poly(N-isopropylacrylamide) particle systems using phase-contrast atomic force microscopy and correlate the distinct interfacial morphology with their bulk molecular architecture. We subsequently introduce a new coarse-grained simulation method that uses augmented potentials to translate this interfacial morphology into the resulting phase behavior upon compression. The main novelty of this method is the possibility to efficiently encode multibody interactions, the effects of which are key to distinguishing between heterostructural (anisotropic collapse) and isostructural (isotropic collapse) phase transitions. Our approach allows us to qualitatively resolve existing discrepancies between experiments and simulations. Notably, we demonstrate the first in silico account of the two-dimensional isostructural transition, which is frequently found in experiments but elusive in simulations. In addition, we provide the first experimental demonstration of a heterostructural transition to a chain phase in a single-component system, which has been theoretically predicted decades ago. Overall, our multiscale framework provides a phenomenological bridge between physicochemical soft-particle characteristics at the molecular scale and nanoscale and the collective self-assembly phenomenology at the macroscale, serving as a stepping stone toward an ultimately more quantitative and predictive design approach.

Internally Stressed and Positionally Disordered Minimal Complexes Yield Glasslike Nonphononic Excitations

Glasses, unlike their crystalline counterparts, exhibit low-frequency nonphononic excitations whose frequencies ω follow a universal D(ω)∼ω^{4} density of states. The process of glass formation generates positional disorder intertwined with mechanical frustration, posing fundamental challenges in understanding the origins of glassy nonphononic excitations. Here we suggest that minimal complexes-mechanically frustrated and positionally disordered local structures-embody the minimal physical ingredients needed to generate glasslike excitations. We investigate the individual effects of mechanical frustration and positional disorder on the vibrational spectrum of isolated minimal complexes, and demonstrate that ensembles of marginally stable minimal complexes yield D(ω)∼ω^{4}. Furthermore, glasslike excitations emerge by embedding a single minimal complex within a perfect lattice. Consequently, minimal complexes offer a conceptual framework to understand glasslike excitations from first principles, as well as a practical computational method for introducing them into solids.

Mechanical disorder of sticky-sphere glasses. II. Thermomechanical inannealability

Many structural glasses feature static and dynamic mechanical properties that can depend strongly on glass formation history. The degree of universality of this history dependence and what it is possibly affected by are largely unexplored. Here we show that the variability of elastic properties of simple computer glasses under thermal annealing depends strongly on the strength of attractive interactions between the glasses' constituent particles-referred to here as glass "stickiness." We find that in stickier glasses the stiffening of the shear modulus with thermal annealing is strongly suppressed, while the thermal-annealing-induced softening of the bulk modulus is enhanced. Our key finding is that the characteristic frequency and density per frequency of soft quasilocalized modes becomes effectively invariant to annealing in very sticky glasses; the latter are therefore deemed "thermomechanically inannealable." The implications of our findings and future research directions are discussed.

Mechanical disorder of sticky-sphere glasses. I. Effect of attractive interactions

Recent literature indicates that attractive interactions between particles of a dense liquid play a secondary role in determining its bulk mechanical properties. Here we show that, in contrast with their apparent unimportance to the bulk mechanics of dense liquids, attractive interactions can have a major effect on macro- and microscopic elastic properties of glassy solids. We study several broadly applicable dimensionless measures of stability and mechanical disorder in simple computer glasses, in which the relative strength of attractive interactions-referred to as "glass stickiness"-can be readily tuned. We show that increasing glass stickiness can result in the decrease of various quantifiers of mechanical disorder, on both macro- and microscopic scales, with a pair of intriguing exceptions to this rule. Interestingly, in some cases strong attractions can lead to a reduction of the number density of soft, quasilocalized modes, by up to an order of magnitude, and to a substantial decrease in their core size, similar to the effects of thermal annealing on elasticity observed in recent works. Contrary to the behavior of canonical glass models, we provide compelling evidence indicating that the stabilization mechanism in our sticky-sphere glasses stems predominantly from the self-organized depletion of interactions featuring large, negative stiffnesses. Finally, we establish a fundamental link between macroscopic and microscopic quantifiers of mechanical disorder, which we motivate via scaling arguments. Future research directions are discussed.

Thermal origin of quasilocalized excitations in glasses

Key aspects of glasses are controlled by the presence of excitations in which a group of particles can rearrange. Surprisingly, recent observations indicate that their density is dramatically reduced and their size decreases as the temperature of the supercooled liquid is lowered. Some theories predict these excitations to cause a gap in the spectrum of quasilocalized modes of the Hessian that grows upon cooling, while others predict a pseudogap D_{L}(ω)∼ω^{α}. To unify these views and observations, we generate glassy configurations of controlled gap magnitude ω_{c} at temperature T=0, using so-called breathing particles, and study how such gapped states respond to thermal fluctuations. We find that (i) the gap always fills up at finite T with D_{L}(ω)≈A_{4}(T)ω^{4} and A_{4}∼exp(-E_{a}/T) at low T, (ii) E_{a} rapidly grows with ω_{c}, in reasonable agreement with a simple scaling prediction E_{a}∼ω_{c}^{4} and (iii) at larger ω_{c} excitations involve fewer particles, as we rationalize, and eventually become stringlike. We propose an interpretation of mean-field theories of the glass transition, in which the modes beyond the gap act as an excitation reservoir, from which a pseudogap distribution is populated with its magnitude rapidly decreasing at lower T. We discuss how this picture unifies the rarefaction as well as the decreasing size of excitations upon cooling, together with a stringlike relaxation occurring near the glass transition.

Simple and Broadly Applicable Definition of Shear Transformation Zones

Plastic deformation in amorphous solids is known to be carried by stress-induced localized rearrangements of a few tens of particles, accompanied by the conversion of elastic energy to heat. Despite their central role in determining how glasses yield and break, the search for a simple and generally applicable definition of the precursors of those plastic rearrangements-the so-called shear transformation zones (STZs)-is still ongoing. Here we present a simple definition of STZs-based solely on the harmonic approximation of a glass's energy. We explain why and demonstrate directly that our proposed definition of plasticity carriers in amorphous solids is more broadly applicable compared to anharmonic definitions put forward previously. Finally, we offer an open-source library that analyzes low-lying STZs in computer glasses and in laboratory materials such as dense colloidal suspensions for which the harmonic approximation is accessible. Our results constitute a physically motivated methodological advancement towards characterizing mechanical disorder in glasses, and understanding how they yield.

Stable glassy configurations of the Kob-Andersen model using swap Monte Carlo

The swap Monte Carlo algorithm allows the preparation of highly stable glassy configurations for a number of glass-formers but is inefficient for some models, such as the much studied binary Kob-Andersen (KA) mixture. We have recently developed generalizations to the KA model where swap can be very effective. Here, we show that these models can, in turn, be used to considerably enhance the stability of glassy configurations in the original KA model at no computational cost. We successfully develop several numerical strategies both in and out of equilibrium to achieve this goal and show how to optimize them. We provide several physical measurements indicating that the proposed algorithms considerably enhance mechanical and thermodynamic stability in the KA model, including a transition toward brittle yielding behavior. Our results thus pave the way for future studies of stable glasses using the KA model.

Flow heterogeneities in supercooled liquids and glasses under shear

Using extensive nonequilibrium molecular dynamics simulations, we investigate a glass-forming binary Lennard-Jones mixture under shear. Both supercooled liquids and glasses are considered. Our focus is on the characterization of inhomogeneous flow patterns such as shear bands that appear as a transient response to the external shear. For the supercooled liquids, we analyze the crossover from Newtonian to non-Newtonian behavior with increasing shear rate γ[over ̇]. Above a critical shear rate γ[over ̇]_{c} where a non-Newtonian response sets in, the transient dynamics are associated with the occurrence of short-lived vertical shear bands, i.e., bands of high mobility that form perpendicular to the flow direction. In the glass states, long-lived horizontal shear bands, i.e., bands of high mobility parallel to the flow direction, are observed in addition to vertical ones. The systems with shear bands are characterized in terms of mobility maps, stress-strain relations, mean-squared displacements, and (local) potential energies. The initial formation of a horizontal shear band provides an efficient stress release, corresponds to a local minimum of the potential energy, and is followed by a slow broadening of the band towards the homogeneously flowing fluid in the steady state. Whether a horizontal or a vertical shear band forms cannot be predicted from the initial undeformed sample. Furthermore, we show that with increasing system size, the probability for the occurrence of horizontal shear bands increases.

Non-Markovian Methods in Glass Transition

A model for the heterogeneity of local dynamics in polymer and other glass-forming materials is provided here. The fundamental characteristics of the glass transition phenomenology emerge when simulating a condensed matter open cluster that has a strong interaction with its heterogeneous environment. General glass transition features, such as non-exponential structural relaxations, the slowing down of relaxation times with temperature and specific off-equilibrium glassy dynamics can be reproduced by non-Markovian dynamics simulations with the minimum computer resources. Non-Markovian models are shown to be useful tools for obtaining insights into the complex dynamics involved in the glass transition phenomenon, including whether or not there is a need for a growing correlation length or the relationship between the non-exponentiality of structural relaxations and dynamic heterogeneity.

Ultrastable Metallic Glasses In Silico

We develop a generic strategy and simple numerical models for multicomponent metallic glasses for which the swap Monte Carlo algorithm can produce highly stable equilibrium configurations equivalent to experimental systems cooled more than 10^{7} times slower than in conventional simulations. This paves the way for a deeper understanding of the thermodynamic, dynamic, and mechanical properties of metallic glasses. As first applications, we considerably extend configurational entropy measurements down to the experimental glass temperature, and demonstrate a qualitative change of the mechanical response of metallic glasses of increasing stability toward brittleness.

Universality of the Nonphononic Vibrational Spectrum across Different Classes of Computer Glasses

It has been recently established that the low-frequency spectrum of simple computer glass models is populated by soft, quasilocalized nonphononic vibrational modes whose frequencies ω follow a gapless, universal distribution D(ω)∼ω^{4}. While this universal nonphononic spectrum has been shown to be robust to varying the glass history and spatial dimension, it has so far only been observed in simple computer glasses featuring radially symmetric, pairwise interaction potentials. Consequently, the relevance of the universality of nonphononic spectra seen in simple computer glasses to realistic laboratory glasses remains unclear. Here, we demonstrate the emergence of the universal ω^{4} nonphononic spectrum in a broad variety of realistic computer glass models, ranging from tetrahedral network glasses with three-body interactions, through molecular glasses and glassy polymers, to bulk metallic glasses. Taken together with previous observations, our results indicate that the low-frequency nonphononic vibrational spectrum of any glassy solid quenched from a melt features the universal ω^{4} law, independently of the nature of its microscopic interactions.

Thermodynamic equivalence of cyclic shear and deep cooling in glass formers

The extreme slowing down associated with glass formation in experiments and in simulations results in serious difficulties to prepare deeply quenched, well annealed, glassy material. Recently, methods to achieve such deep quenching were proposed, including vapor deposition on the experimental side and "swap Monte Carlo" and oscillatory shearing on the simulation side. The relation between the resulting glasses under different protocols remains unclear. Here we show that oscillatory shear and swap Monte Carlo result in thermodynamically equivalent glasses sharing the same statistical mechanics and similar mechanical responses under external strain.

Glass Stability Changes the Nature of Yielding under Oscillatory Shear

We perform molecular dynamics simulations to investigate the effect of a glass preparation on its yielding transition under oscillatory shear. We use swap Monte Carlo to investigate a broad range of glass stabilities from poorly annealed to highly stable systems. We observe a qualitative change in the nature of yielding, which evolves from ductile to brittle as glass stability increases. Our results disentangle the relative role of mechanical and thermal annealing on the mechanical properties of amorphous solids, which is relevant for various experimental situations from the rheology of soft materials to fatigue failure in metallic glasses.

Athermal Fluctuations in Disordered Crystals

We analyze the fluctuations in particle positions and interparticle forces in disordered crystals composed of jammed soft particles in the limit of weak disorder. We demonstrate that such athermal systems are fundamentally different from their thermal counterparts, characterized by constrained fluctuations of forces perpendicular to the lattice directions. We develop a disorder perturbation expansion in polydispersity about the crystalline state, which we use to derive exact results to linear order. We show that constrained fluctuations result as a consequence of local force balance conditions, and are characterized by non-Gaussian distributions, which we derive exactly. We analytically predict several properties of such systems, including the scaling of the average coordination with polydispersity and packing fraction, which we verify with numerical simulations using soft disks with one-sided harmonic interactions.

Pinching a glass reveals key properties of its soft spots

It is now well established that glasses feature quasilocalized nonphononic excitations-coined "soft spots"-, which follow a universal [Formula: see text] density of states in the limit of low frequencies ω. All glass-specific properties, such as the dependence on the preparation protocol or composition, are encapsulated in the nonuniversal prefactor of the universal [Formula: see text] law. The prefactor, however, is a composite quantity that incorporates information both about the number of quasilocalized nonphononic excitations and their characteristic stiffness, in an apparently inseparable manner. We show that by pinching a glass-i.e., by probing its response to force dipoles-one can disentangle and independently extract these two fundamental pieces of physical information. This analysis reveals that the number of quasilocalized nonphononic excitations follows a Boltzmann-like law in terms of the parent temperature from which the glass is quenched. The latter, sometimes termed the fictive (or effective) temperature, plays important roles in nonequilibrium thermodynamic approaches to the relaxation, flow, and deformation of glasses. The analysis also shows that the characteristic stiffness of quasilocalized nonphononic excitations can be related to their characteristic size, a long sought-for length scale. These results show that important physical information, which is relevant for various key questions in glass physics, can be obtained through pinching a glass.

Nonlinear quasilocalized excitations in glasses: True representatives of soft spots

Structural glasses formed by quenching a melt possess a population of soft quasilocalized excitations-often called "soft spots"-that are believed to play a key role in various thermodynamic, transport, and mechanical phenomena. Under a narrow set of circumstances, quasilocalized excitations assume the form of vibrational (normal) modes, that are readily obtained by a harmonic analysis of the multidimensional potential energy. In general, however, direct access to the population of quasilocalized modes via harmonic analysis is hindered by hybridizations with other low-energy excitations, e.g., phonons. In this series of papers we reintroduce and investigate the statistical-mechanical properties of a class of low-energy quasilocalized modes-coined here nonlinear quasilocalized excitations (NQEs)-that are defined via an anharmonic (nonlinear) analysis of the potential-energy landscape of a glass, and do not hybridize with other low-energy excitations. In this paper, we review the theoretical framework that embeds a micromechanical definition of NQEs. We demonstrate how harmonic quasilocalized modes hybridize with other soft excitations, whereas NQEs properly represent soft spots without hybridization. We show that NQEs' energies converge to the energies of the softest, nonhybridized harmonic quasilocalized modes, cementing their status as true representatives of soft spots in structural glasses. Finally, we perform a statistical analysis of the mechanical properties of NQEs, which results in a prediction for the distribution of potential-energy barriers that surround typical inherent states of structural glasses, as well as a prediction for the distribution of local strain thresholds to plastic instability.

Universal non-mean-field scaling in the density of states of amorphous solids

Amorphous solids have excess soft modes in addition to the phonon modes described by the Debye theory. Recent numerical results show that if the phonon modes are carefully removed, the density of state of the excess soft modes exhibit universal quartic scaling, independent of the interaction potential, preparation protocol, and spatial dimensions. We hereby provide a theoretical framework to describe this universal scaling behavior. For this purpose, we extend the mean-field theory to include the effects of finite-dimensional fluctuation. Based on a semiphenomenological argument, we show that mean-field quadratic scaling is replaced by the quartic scaling in finite dimensions. Furthermore, we apply our formalism to explain the pressure and protocol dependence of the excess soft modes.

Theory for the density of interacting quasilocalized modes in amorphous solids

Quasilocalized modes appear in the vibrational spectrum of amorphous solids at low frequency. Though never formalized, these modes are believed to have a close relationship with other important local excitations, including shear transformations and two-level systems. We provide a theory for their frequency density, D_{L}(ω)∼ω^{α}, that establishes this link for systems at zero temperature under quasistatic loading. It predicts two regimes depending on the density of shear transformations P(x)∼x^{θ} (with x the additional stress needed to trigger a shear transformation). If θ>1/4, then α=4 and a finite fraction of quasilocalized modes form shear transformations, whose amplitudes vanish at low frequencies. If θ<1/4, then α=3+4θ and all quasilocalized modes form shear transformations with a finite amplitude at vanishing frequencies. We confirm our predictions numerically.