Local elasticity map and plasticity in a model Lennard-Jones glass

Acoustic excitations and elastic heterogeneities in disordered solids

In the recent years, much attention has been devoted to the inhomogeneous nature of the mechanical response at the nanoscale in disordered solids. Clearly, the elastic heterogeneities that have been characterized in this context are expected to strongly affect the nature of the sound waves which, in contrast to the case of perfect crystals, cannot be completely rationalized in terms of phonons. Building on previous work on a toy model showing an amorphization transition, we investigate the relationship between sound waves and elastic heterogeneities in a unified framework by continuously interpolating from the perfect crystal, through increasingly defective phases, to fully developed glasses. We provide strong evidence of a direct correlation between sound wave features and the extent of the heterogeneous mechanical response at the nanoscale.

Surface shear-transformation zones in amorphous solids

We perform a systematic study of the characteristics of shear transformation zones (STZs) that nucleate at free surfaces of two-dimensional amorphous solids subject to tensile loading using two different atomistic simulation methods, the standard athermal, quasistatic (AQ) approach and our recently developed self-learning metabasin escape (SLME) method, to account for the finite temperature and strain-rate effects. In the AQ, or strain-driven limit, the nonaffine displacement fields of surface STZs decay exponentially away from their centers at similar decay rates as their bulk counterparts, though the direction of maximum nonaffine displacement is tilted away from the tensile axis due to surface effects. Using the SLME method at room temperature and at the high strain rates that are seen in classical molecular dynamics simulations, the characteristics for both bulk and surface STZs are found to be identical to those seen in the AQ simulations. However, using the SLME method at room temperature and experimentally relevant strain rates, we find a transition in the surface STZ characteristics where a loss in the characteristic angular tensile-compression symmetry is observed. Finally, the thermally activated surface STZs exhibit a slower decay rate in the nonaffine displacement field than do strain-driven surface STZs, which is characterized by a larger drop in potential energy resulting from STZ nucleation that is enabled by the relative compliance of the surface as compared to the bulk.

Order parameter for structural heterogeneity in disordered solids

We construct a structural order parameter from the energy equipartition of normal modes of vibration to quantify the structural heterogeneity in disordered solids. The order parameter exhibits strong spatial correlations with low-temperature dynamics and local structural entropy. To characterize the role of particles with the most defective local structures identified by the order parameter, we pin them and measure the system response. It turns out that particles with the largest value of the order parameter are responsible for the quasilocalized low-frequency vibration, instability, softening, and nonaffinity of disordered solids. The order parameter thus crucially links the heterogeneous structure to low-temperature dynamics and mechanical properties of disordered solids.

Microrheology of colloidal systems

Microrheology was proposed almost twenty years ago as a technique to obtain rheological properties in soft matter from the microscopic motion of colloidal tracers used as probes, either freely diffusing in the host medium, or subjected to external forces. The former case is known as passive microrheology, and is based on generalizations of the Stokes-Einstein relation between the friction experienced by the probe and the host-fluid viscosity. The latter is termed active microrheology, and extends the measurement of the friction coefficient to the nonlinear-response regime of strongly driven probes. In this review article, we discuss theoretical models available in the literature for both passive and active microrheology, focusing on the case of single-probe motion in model colloidal host media. A brief overview of the theory of passive microrheology is given, starting from the work of Mason and Weitz. Further developments include refined models of the host suspension beyond that of a Newtonian-fluid continuum, and the investigation of probe-size effects. Active microrheology is described starting from microscopic equations of motion for the whole system including both the host-fluid particles and the tracer; the many-body Smoluchowski equation for the case of colloidal suspensions. At low fluid densities, this can be simplified to a two-particle equation that allows the calculation of the friction coefficient with the input of the density distribution around the tracer, as shown by Brady and coworkers. The results need to be upscaled to agree with simulations at moderate density, in both the case of pulling the tracer with a constant force or dragging it at a constant velocity. The full many-particle equation has been tackled by Fuchs and coworkers, using a mode-coupling approximation and the scheme of integration through transients, valid at high densities. A localization transition is predicted for a probe embedded in a glass-forming host suspension. The nonlinear probe-friction coefficient is calculated from the tracer's position correlation function. Computer simulations show qualitative agreement with the theory, but also some unexpected features, such as superdiffusive motion of the probe related to the breaking of nearest-neighbor cages. We conclude with some perspectives and future directions of theoretical models of microrheology.

Experimental signatures of a nonequilibrium phase transition governing the yielding of a soft glass

We present direct experimental signatures of a nonequilibrium phase transition associated with the yield point of a prototypical soft solid-a binary colloidal glass. By simultaneously quantifying single-particle dynamics and bulk mechanical response, we identified the threshold for the onset of irreversibility with the yield strain. We extracted the relaxation time from the transient behavior of the loss modulus and found that it diverges in the vicinity of the yield strain. This critical slowing down is accompanied by a growing correlation length associated with the size of regions of high Debye-Waller factor, which are precursors to yield events in glasses. Our results affirm that the paradigm of nonequilibrium critical phenomena is instrumental in achieving a holistic understanding of yielding in soft solids.

Evolution of elastic heterogeneity during aging in metallic glasses

The properties of glasses vary widely depending on the way they are prepared, even though their structures appear similar. We show that the local potential energy landscape (PEL) sensitively reflects the stability differences through simulation of local structural excitation in a model metallic glass. It is observed that the spectrum of local structural excitation develops a pseudogap at low energies as the glass becomes more stable. We also demonstrate that the spatial variation of the atomic level shear modulus, rather than the distribution of the magnitude of the single atom shear modulus, is more closely related to the nature of the PEL and the stabilities of glasses. In particular, local aggregation of atoms with low shear modulus greatly contributes to instability of the system.

Spatiotemporal correlations between plastic events in the shear flow of athermal amorphous solids

The slow flow of amorphous solids exhibits striking heterogeneities: swift localised particle rearrangements take place in the midst of a more or less homogeneously deforming medium. Recently, experimental as well as numerical work has revealed spatial correlations between these flow heterogeneities. Here, we use molecular dynamics (MD) simulations to characterise the rearrangements and systematically probe their correlations both in time and in space. In particular, these correlations display a four-fold azimuthal symmetry characteristic of shear stress redistribution in an elastic medium and we unambiguously detect their increase in range with time. With increasing shear rate, correlations become shorter-ranged. In addition, we study a coarse-grained model motivated by the observed flow characteristics and challenge its predictions directly with the MD simulations. While the model captures both macroscopic and local properties rather satisfactorily, the agreement with respect to the spatiotemporal correlations is at most qualitative. The discrepancies provide important insight into relevant physics that is missing in all related coarse-grained models that have been developed for the flow of amorphous materials so far, namely the finite shear wave velocity and the impact of elastic heterogeneities on stress redistribution.

Crossover from random three-dimensional avalanches to correlated nano shear bands in metallic glasses

When applying mechanical stress to a bulk metallic glass it responds with elastic and/or plastic deformations. A comprehensive microscopic theory for the plasticity of amorphous solids remains an open task. Shear transformation zones consisting of dozens of atoms have been identified as smallest units of deformation. The connexion between local formation of shear transformations zones and the creation of macroscopic shear bands can be made using statistical analysis of stress/energy drops or strain slips during mechanical loading. Numerical work has proposed a power law dependence of those energy drops. Here we present an approach to circumvent the experimental resolution problem using a waiting time analysis. We report on the power law-distributed deformation behaviour and the observation of a crossover in the waiting times statistics. This crossover indicates a transition in the plastic deformation behaviour from three-dimensional random activity to a two-dimensional nano shear band sliding.

Time-dependent elastic response to a local shear transformation in amorphous solids

The elastic response of a two-dimensional amorphous solid to induced local shear transformations, which mimic the elementary plastic events occurring in deformed glasses, is investigated via molecular-dynamics simulations. We show that for different spatial realizations of the transformation, despite relative fluctuations of order one, the long-time equilibrium response averages out to the prediction of the Eshelby inclusion problem for a continuum elastic medium. We characterize the effects of the underlying dynamics on the propagation of the elastic signal. A crossover from a propagative transmission in the case of weakly damped dynamics to a diffusive transmission for strong damping is evidenced. In the latter case, the full time-dependent elastic response is in agreement with the theoretical prediction, obtained by solving the diffusion equation for the displacement field in an elastic medium.

Predicting plasticity with soft vibrational modes: from dislocations to glasses

We show that quasilocalized low-frequency modes in the vibrational spectrum can be used to construct soft spots, or regions vulnerable to rearrangement, which serve as a universal tool for the identification of flow defects in solids. We show that soft spots not only encode spatial information, via their location, but also directional information, via directors for particles within each soft spot. Single crystals with isolated dislocations exhibit low-frequency phonon modes that localize at the core, and their polarization pattern predicts the motion of atoms during elementary dislocation glide in two and three dimensions in exquisite detail. Even in polycrystals and disordered solids, we find that the directors associated with particles in soft spots are highly correlated with the direction of particle displacements in rearrangements.

Soft modes and nonaffine rearrangements in the inherent structures of supercooled liquids

We find that the hierarchical organization of the potential energy landscape in a model supercooled liquid can be related to a change in the spatial distribution of soft normal modes. For groups of nearby minima, between which fast relaxation processes typically occur, the localization of the soft modes is very similar. The spatial distribution of soft regions changes, instead, for minima between which transitions relevant to structural relaxation occur. This may be the reason why the soft modes are able to predict spatial heterogeneities in the dynamics. Nevertheless, the very softest modes are only weakly correlated with dynamical heterogeneities and instead show higher statistical overlap with regions in the local minima that would undergo nonaffine rearrangements if subjected to a shear deformation. This feature of the supercooled liquid is reminiscent of the behavior of nonaffine deformations in amorphous solids, where the very softest modes identify the loci of plastic instabilities.

Visualizing the strain evolution during the indentation of colloidal glasses

We use an analog of nanoindentation on a colloidal glass to elucidate the incipient plastic deformation of glasses. By tracking the motion of the individual particles in three dimensions, we visualize the strain field and glass structure during the emerging deformation. At the onset of flow, we observe a power-law distribution of strain indicating strongly correlated deformation, and reflecting a critical state of the glass. At later stages, the strain acquires a Gaussian distribution, indicating that plastic events become uncorrelated. Investigation of the glass structure using both static and dynamic measures shows a weak correlation between the structure and the emerging strain distribution. These results indicate that the onset of plasticity is governed by strong power-law correlations of strain, weakly biased by the heterogeneous glass structure.

Effect of inertia on sheared disordered solids: critical scaling of avalanches in two and three dimensions

Molecular dynamics simulations with varying damping are used to examine the effects of inertia and spatial dimension on sheared disordered solids in the athermal quasistatic limit. In all cases the distribution of avalanche sizes follows a power law over at least three orders of magnitude in dissipated energy or stress drop. Scaling exponents are determined using finite-size scaling for systems with 10(3)-10(6) particles. Three distinct universality classes are identified corresponding to overdamped and underdamped limits, as well as a crossover damping that separates the two regimes. For each universality class, the exponent describing the avalanche distributions is the same in two and three dimensions. The spatial extent of plastic deformation is proportional to the energy dissipated in an avalanche. Both rise much more rapidly with system size in the underdamped limit where inertia is important. Inertia also lowers the mean energy of configurations sampled by the system and leads to an excess of large events like that seen in earthquake distributions for individual faults. The distribution of stress values during shear narrows to zero with increasing system size and may provide useful information about the size of elemental events in experimental systems. For overdamped and crossover systems the stress variation scales inversely with the square root of the system size. For underdamped systems the variation is determined by the size of the largest events.

Strain-rate and temperature-driven transition in the shear transformation zone for two-dimensional amorphous solids

We couple the recently developed self-learning metabasin escape algorithm, which enables efficient exploration of the potential energy surface (PES), with shear deformation to elucidate strain-rate and temperature effects on the shear transformation zone (STZ) characteristics in two-dimensional amorphous solids. In doing so, we report a transition in the STZ characteristics that can be obtained through either increasing the temperature or decreasing the strain rate. The transition separates regions having two distinct STZ characteristics. Specifically, at high temperatures and high strain rates, we show that the STZs have characteristics identical to those that emerge from purely strain-driven, athermal quasistatic atomistic calculations. At lower temperatures and experimentally relevant strain rates, we use the newly coupled PES + shear deformation method to show that the STZs have characteristics identical to those that emerge from a purely thermally activated state. The specific changes in STZ characteristics that occur in moving from the strain-driven to thermally activated STZ regime include a 33% increase in STZ size, faster spatial decay of the displacement field, a change in the deformation mechanism inside the STZ from shear to tension, a reduction in the stress needed to nucleate the first STZ, and finally a notable loss in characteristic quadrupolar symmetry of the surrounding elastic matrix that has previously been seen in athermal, quasistatic shear studies of STZs.

Heterogeneous shear elasticity of glasses: the origin of the boson peak

The local elasticity of glasses is known to be inhomogeneous on a microscopic scale compared to that of crystalline materials. Their vibrational spectrum strongly deviates from that expected from Debye's elasticity theory: The density of states deviates from Debye's law, the sound velocity shows a negative dispersion in the boson-peak frequency regime and there is a strong increase of the sound attenuation near the boson-peak frequency. By comparing a mean-field theory of shear-elastic heterogeneity with a large-scale simulation of a soft-sphere glass we demonstrate that the observed anomalies in glasses are caused by elastic heterogeneity. By observing that the macroscopic bulk modulus is frequency independent we show that the boson-peak-related vibrational anomalies are predominantly due to the spatially fluctuating microscopic shear stresses. It is demonstrated that the boson-peak arises from the steep increase of the sound attenuation at a frequency which marks the transition from wave-like excitations to disorder-dominated ones.

Measuring spatial distribution of the local elastic modulus in glasses

Glasses exhibit spatially inhomogeneous elastic properties, which can be investigated by measuring their elastic moduli at a local scale. Various methods to evaluate the local elastic modulus have been proposed in the literature. A first possibility is to measure the local stress-local strain curve and to obtain the local elastic modulus from the slope of the curve or, equivalently, to use a local fluctuation formula. Another possible route is to assume an affine strain and to use the applied global strain instead of the local strain for the calculation of the local modulus. Most recently, a third technique has been introduced, which is easy to be implemented and has the advantage of low computational cost. In this contribution, we compare these three approaches by using the same model glass and reveal the differences among them caused by the nonaffine deformations.

Analysis of local properties during a scratch test on a polymeric surface using molecular dynamics simulations

This work demonstrates a possible route to connect a particle (chain) based understanding with continuum mechanical questions about contact mechanics. The bond orientation, chain conformation and stress field of a polymer film were analyzed during scratch tests (tangential contact) using a molecular dynamics (MD) simulation approach. Scratch tests with a conical tip at constant scratching velocity were simulated on linear amorphous polymer surfaces at various temperatures and roughnesses of the tip and for various interactions between the tip and the particles of the polymer chains. The second Legendre polynomial (computed for small domains around the tip) gave the bond orientation inside the polymer film during sliding of the tip. The gyration tensor (layer-resolved in the direction of the polymer film thickness) provided information about the conformation of the polymer chains. These results allowed us to argue in favor of Briscoe's hypothesis (thin film sheared vs. "bulk" compressive behavior) concerning the friction properties of the polymer surfaces. Finally, the first stress measurements of the virial stress tensor (in sub-boxes placed in the MD cell) revealed a complex combination between compressive hydrostatic pressure and shear stress, which may be interpreted as a complex sheared domain at the interface.

Relationship between bond-breakage correlations and four-point correlations in heterogeneous glassy dynamics: configuration changes and vibration modes

We investigate the dynamic heterogeneities of glassy particle systems in the theoretical schemes of bond breakage and four-point correlation functions. In the bond-breakage scheme, we introduce the structure factor S(b)(q,t) and the susceptibility χ(b)(t) to detect the spatial correlations of configuration changes. Here χ(b)(t) attains a maximum at t=t(b)(max) as a function of time t, where the fraction of the particles with broken bonds φ(b)(t) is about 1/2. In the four-point scheme, treating the structure factor S(4)(q,t) and the susceptibility χ(4)(t), we detect superpositions of the heterogeneity of bond breakage and that of thermal low-frequency vibration modes. While the former grows slowly, the latter emerges quickly to exhibit complex space-time behavior. In two dimensions, the vibration modes extending over the system yield significant contributions to the four-point correlations, which depend on the system size logarithmically. A maximum of χ(4)(t) is attained at t=t(4)(max), where these two contributions become of the same order. As a result, t(4)(max) is considerably shorter than t(b)(max).

Impulsive correction to the elastic moduli obtained using the stress-fluctuation formalism in systems with truncated pair potential

The truncation of a pair potential at a distance rc is well known to imply, in general, an impulsive correction to the pressure and other moments of the first derivatives of the potential. That, depending on rc, the truncation may also be of relevance to higher derivatives is shown theoretically for the Born contributions to the elastic moduli obtained using the stress-fluctuation formalism in d dimensions. Focusing on isotropic liquids for which the shear modulus G must vanish by construction, the predicted corrections are tested numerically for binary mixtures and polydisperse Lennard-Jones beads in, respectively, d=3 and 2 dimensions. Both models being glass formers, we comment briefly on the temperature (T) dependence of the (corrected) shear modulus G(T) around the glass transition temperature Tg.

Mechanical properties of thin confined polymer films close to the glass transition in the linear regime of deformation: theory and simulations

Over the past twenty years experiments performed on thin polymer films deposited on substrates have shown that the glass transition temperature T(g) can either decrease or increase depending on the strength of the interactions. Over the same period, experiments have also demonstrated that the dynamics in liquids close to the glass transition temperature is strongly heterogeneous, on the scale of a few nanometers. A model for the dynamics of non-polar polymers, based on percolation of slow subunits, has been proposed and developed over the past ten years. It proposes a unified mechanism regarding these two features. By extending this model, we have developed a 3D model, solved by numerical simulations, in order to describe and calculate the mechanical properties of polymers close to the glass transition in the linear regime of deformation, with a spatial resolution corresponding to the subunit size. We focus on the case of polymers confined between two substrates with non-negligible interactions between the polymer and the substrates, a situation which may be compared to filled elastomers. We calculate the evolution of the elastic modulus as a function of temperature, for different film thicknesses and polymer-substrate interactions. In particular, this allows to calculate the corresponding increase of glass transition temperature, up to 20 K in the considered situations. Moreover, between the bulk T(g) and T(g) + 50 K the modulus of the confined layers is found to decrease very slowly in some cases, with moduli more than ten times larger than that of the pure matrix at temperatures up to T(g) + 50 K. This is consistent with what is observed in reinforced elastomers. This slow decrease of the modulus is accompanied by huge fluctuations of the stress at the scale of a few tens of nanometers that may even be negative as compared to the solicitation, in a way that may be analogous to mechanical heterogeneities observed recently in molecular dynamics simulations. As a consequence, confinement may result not only in an increase of the glass transition temperature, but in a huge broadening of the glass transition.

Mechanical behavior of linear amorphous polymers: comparison between molecular dynamics and finite-element simulations

This paper studies the rheology of weakly entangled polymer melts and films in the glassy domain and near the rubbery domain using two different methods: molecular dynamics (MD) and finite element (FE) simulations. In a first step, the uniaxial mechanical behavior of a bulk polymer sample is studied by means of particle-based MD simulations. The results are in good agreement with experimental data, and mechanical properties may be computed from the simulations. This uniaxial mechanical behavior is then implemented in FE simulations using an elasto-viscoelasto-viscoplastic constitutive law in a continuum mechanics (CM) approach. In a second step, the mechanical response of a polymer film during an indentation test is modeled with the MD method and with the FE simulations using the same constitutive law. Good agreement is found between the MD and CM results. This work provides evidence in favor of using MD simulations to investigate the local physics of contact mechanics, since the volume elements studied are representative and thus contain enough information about the microstructure of the polymer model, while surface phenomena (adhesion and surface tension) are naturally included in the MD approach.

Local anisotropy in globally isotropic granular packings

We report on two-dimensional computer simulations of frictionless granular packings at various area fractions φ above the jamming point φ(c). We measure the anisotropy in coarse-grained stress ε(s) and shear modulus ε(m) as functions of coarse-graining scale, R. ε(s) can be collapsed onto a master curve after rescaling R by a characteristic length scale ξ and ε(s) by an anisotropy magnitude A. Both A and ξ accelerate as φ→φ(c) from above, consistent with a divergence at φ(c). ε(m) shows no characteristic length scale and has a nontrivial power-law form, ε(m)~R(-0.62), over almost the entire range of R at all φ. These results suggest that the force chains present in the spatial structure of the quenched stress may be governed by different physics than the anomalous elastic response near jamming.

Towards ultrastrong glasses

The development of new glassy materials is key for addressing major global challenges in energy, medicine, and advanced communications systems. For example, thin, flexible, and large-area glass substrates will play an enabling role in the development of flexible displays, roll-to-roll processing of solar cells, next-generation touch-screen devices, and encapsulation of organic semiconductors. The main drawback of glass and its limitation for these applications is its brittle fracture behavior, especially in the presence of surface flaws, which can significantly reduce the practical strength of a glass product. Hence, the design of new ultrastrong glassy materials and strengthening techniques is of crucial importance. The main issues regarding glass strength are discussed, with an emphasis on the underlying microscopic mechanisms that are responsible for mechanical properties. The relationship among elastic properties and fracture behavior is also addressed, focusing on both oxide and metallic glasses. From a theoretical perspective, atomistic modeling of mechanical properties of glassy materials is considered. The topological origin of these properties is also discussed, including its relation to structural and chemical heterogeneities. Finally, comments are given on several toughening strategies for increasing the damage resistance of glass products.

Measurement of correlations between low-frequency vibrational modes and particle rearrangements in quasi-two-dimensional colloidal glasses

We investigate correlations between low-frequency vibrational modes and rearrangements in two-dimensional colloidal glasses composed of thermosensitive microgel particles, which readily permit variation of the sample packing fraction. At each packing fraction, the particle displacement covariance matrix is measured and used to extract the vibrational spectrum of the "shadow" colloidal glass (i.e., the particle network with the same geometry and interactions as the sample colloid but absent damping). Rearrangements are induced by successive, small reductions in the packing fraction. The experimental results suggest that low-frequency quasilocalized phonon modes in colloidal glasses, i.e., modes that present low energy barriers for system rearrangements, are spatially correlated with rearrangements in this thermal system.

Vibrational modes identify soft spots in a sheared disordered packing

We analyze low-frequency vibrational modes in a two-dimensional, zero-temperature, quasistatically sheared model glass to identify a population of structural "soft spots" where particle rearrangements are initiated. The population of spots evolves slowly compared to the interval between particle rearrangements, and the soft spots are structurally different from the rest of the system. Our results suggest that disordered solids flow via localized rearrangements that tend to occur at soft spots, which are analogous to dislocations in crystalline solids.

Simulated glass-forming polymer melts: glass transition temperature and elastic constants of the glassy state

By means of molecular-dynamics simulation we study a flexible and a semiflexible bead-spring model for a polymer melt on cooling through the glass transition. Results for the glass transition temperature T(g) and for the elastic properties of the glassy state are presented. We find that T(g) increases with chain length N and is for all N larger for the semiflexible model. The N dependence of T(g) is compared to experimental results from the literature. Furthermore, we characterize the polymer glass below T(g) via its elastic properties, i.e., via the Lamé coefficients λ and μ. The Lamé coefficients are determined from the fluctuation formalism which allows to split λ and μ into affine (Born term) and nonaffine (fluctuation term) contributions. We find that the fluctuation term represents a substantial correction to the Born term. Since the Born terms for λ and μ are identical, the fluctuation terms are responsible for the different temperature dependence of the Lamé coefficients. While λ decreases linearly on approaching T(g) from below, the shear modulus μ displays a much stronger decrease near T(g). From the present simulation data it is not possible to decide whether μ takes a finite value at T(g), as would be expected from mode-coupling theory, or vanishes continuously, as suggested by recent work from replica theory.

Local elastic properties of a metallic glass

The nature of non-crystalline materials causes the local potential energy of a cluster of atoms or molecules to vary significantly in space. Different configurations of an ensemble of atoms in a metallic glass lead therefore to a distribution of elastic constants which also changes in space. This is totally different to their crystalline counterparts, where a long-range order exists in space and therefore a much more unified elastic modulus is expected. Using atomic force acoustic microscopy, we present data which show that the local so-called indentation modulus M indeed exhibits a wide distribution on a scale below 10 nm in amorphous PdCuSi, with ΔM/M≈30%. About 10(4) atoms are probed in an individual measurement. Crystallized PdCuSi shows a variation that is 10-30 times smaller and which is determined by the resolution of the microscope and by the polycrystalline structure of the material.

Role of local order in the small-scale plasticity of model amorphous materials

The mechanism of plastic flow in amorphous solids involves nucleation-controlled shear transformations, triggered under stress from fertile sites. However, the origin of these sites is still a matter of debate. In this paper, we show that the connection between local plastic activity and coordination defects in amorphous systems depends on the nature of the interatomic interactions. In particular, the directionality of the bonds, as quantified by the three-body term in Stillinger-Weber-like interactions, affects not only the role of local defects, but also the size of the plastic rearrangements, and the global stress-strain behavior. We study the effect of structure changes due to different quenching rates as well. We conclude the paper by a comparison between amorphous plasticity and the Peierls-Nabarro theory of plasticity in crystals.

Continuum modeling of bulk metallic glasses and composites

At low temperatures, monolithic bulk metallic glasses (BMGs) exhibit high strength and large elasticity limits. On the other hand, BMGs lack overall ductility due to highly localized deformation mechanisms. Recent experimental findings suggest that the problem of catastrophic failure by shear band propagation in BMGs can be mitigated by tailoring microstructural features at different length scales to promote more homogeneous plastic deformation. Herein, based on a continuum approach, we present a quantitative analysis of the effects of microstructure on the deformation behavior of monolithic BMGs and BMG composites. In particular, simulations highlight the importance of short-ranged structural correlations on ductility in monolithic BMGs and demonstrate that particle size controls the ductility of BMG composites. In broader terms, our results provide new avenues for further improvements to the mechanical properties of BMGs.

Plasticity and dynamical heterogeneity in driven glassy materials

Many amorphous glassy materials exhibit complex spatio-temporal mechanical response and rheology, characterized by an intermittent stress strain response and a fluctuating velocity profile. Under quasistatic and athermal deformation protocols this heterogeneous plastic flow was shown to be composed of plastic events of various sizes, ranging from local quadrupolar plastic rearrangements to system spanning shear bands. In this paper, through numerical study of a 2D Lennard-Jones amorphous solid, we generalize the study of the heterogeneous dynamics of glassy materials to the finite shear rate (gamma not equal to 0) and temperature case (T not equal to 0). In practice, we choose an effectively athermal limit (T approximately 0) and focus on the influence of shear rate on the rheology of the glass. In line with previous works we find that the model Lennard-Jones glass follows the rheological behavior of a yield stress fluid with a Herschel-Bulkley response of the form, sigma = sigmaY + c1gamma(beta). The global mechanical response obtained through the use of Molecular Dynamics is shown to converge in the limit gamma --> 0 to the quasistatic limit obtained with an energy minimization protocol. The detailed analysis of the plastic deformation at different shear rates shows that the glass follows different flow regimes. At sufficiently low shear rates the mechanical response reaches a shear-rate-independent regime that exhibits all the characteristics of the quasistatic response (finite-size effects, cascades of plastic rearrangements, yield stress, ...). At intermediate shear rates the rheological properties are determined by the externally applied shear rate and the response deviates from the quasistatic limit. Finally at higher shear the system reaches a shear-rate-independent homogeneous regime. The existence of these three regimes is also confirmed by the detailed analysis of the atomic motion. The computation of the four-point correlation function shows that the transition from the shear-rate-dominated to the quasistatic regime is accompanied by the growth of a dynamical cooperativity length scale xi that is shown to diverge with shear rate as xi is proportional to gamma(-nu), with nu approximately 0.2 -0.3. This scaling is compared with the prediction of a simple model that assumes the diffusive propagation of plastic events.

Predicting plastic flow events in athermal shear-strained amorphous solids

We propose a method to predict the value of the external strain where a generic amorphous solid will fail by a plastic response (i.e., an irreversible deformation), solely on the basis of measurements of the nonlinear elastic moduli. While usually considered fundamentally different, with the elastic properties describing reversible phenomena and plastic failure epitomizing irreversible behavior, we show that the knowledge of some nonlinear elastic moduli is enough to predict where plasticity sets in.

Plastic deformations in crystal, polycrystal, and glass in binary mixtures under shear: collective yielding

Using molecular dynamics simulation, we examine the dynamics of crystal, polycrystal, and glass in a Lennard-Jones binary mixture composed of small and large particles in two dimensions. The crossovers occur among these states as the composition c is varied at fixed size ratio. Shear is applied to a system of 9000 particles in contact with moving boundary layers composed of 1800 particles. The particle configurations are visualized with a sixfold orientation angle αj(t) and a disorder variable Dj(t) defined for particle j, where the latter represents the deviation from hexagonal order. Fundamental plastic elements are classified into dislocation gliding and grain boundary sliding. At any c, large-scale yielding events occur on the acoustic time scale. Moreover, they multiply occur in narrow fragile areas, forming shear bands. The dynamics of plastic flow is highly hierarchical with a wide range of time scales for slow shearing. We also clarify the relationship between the shear stress averaged in the bulk region and the wall stress applied at the boundaries.

From macroscopic yield criteria to atomic stresses in polymer glasses

The relationship between macroscopic shear yield criteria and local stress distributions in deformed polymer glasses is investigated via molecular dynamics simulations on different scales of coarse-graining. Macroscopic shear stresses at the yield point obey a pressure-modified von Mises (pmvM) criterion for many different loading conditions and strain rates. Average local stresses in small volume elements obey the same yield criterion for volumes containing approx. 100 atoms or more. Qualitatively different behavior is observed on smaller scales: the average octahedral atomic shear stress has a simple linear relationship to hydrostatic pressure regardless of macroscopic stress state and failure mode. Local plastic events are identified through a threshold in the mean-squared nonaffine displacement and compared to the local stress state. We find that the pmvM criterion only predicts local yield events when stress and displacements are averaged over at least 100 atoms. By contrast, macroscopic shear yield criteria appear to lose their ability to predict plastic activity on the atomic scale.

Central role of thermal collective strain in the relaxation of structure in a supercooled liquid

The spatial distribution of structural relaxation in a supercooled liquid is studied using molecular dynamics simulations of a two-dimensional binary mixture. It is shown that the spatial heterogeneity of the relaxation along with the time scale of the relaxation is determined, not by the frequency with which particles move a distance pi/2kBragg, but by the frequency with which particles can achieve persistent displacements. We show that these persistent displacements are achieved through the coupled action of local reorganizations and unrecoverable thermal strains.