Particle rearrangement and softening contributions to the nonlinear mechanical response of glasses

Shear softening and hardening of a two-dimensional Yukawa solid

Langevin dynamical simulations are performed to study the elastic behaviors of two-dimensional (2D) solid dusty plasmas under the periodic shear deformation. The frequency- and strain-dependent shear moduli G(ω,γ) of our simulated 2D Yukawa solid are calculated from the ratio of the shear stress to strain in different orientations. The shear-softening and -hardening properties in different lattice orientations are discovered from the obtained G(ω,γ). The component of the elastic constant tensor corresponding to the shear deformation is also calculated, whose variation trend exactly agrees with the discovered shear-softening and -hardening features in different shear directions. It is also found that the shear modulus of the 2D Yukawa solid always increases monotonically with the frequency.

Mean-Field Predictions of Scaling Prefactors Match Low-Dimensional Jammed Packings

No known analytic framework precisely explains all the phenomena observed in jamming. The replica theory for glasses and jamming is a mean-field theory which attempts to do so by working in the limit of infinite dimensions, such that correlations between neighbors are negligible. As such, results from this mean-field theory are not guaranteed to be observed in finite dimensions. However, many results in mean field for jamming have been shown to be exact or nearly exact in low dimensions. This suggests that the infinite dimensional limit is not necessary to obtain these results. In this Letter, we perform precision measurements of jamming scaling relationships between pressure, excess packing fraction, and number of excess contacts from dimensions 2-10 in order to extract the prefactors to these scalings. While these prefactors should be highly sensitive to finite dimensional corrections, we find the mean-field predictions for these prefactors to be exact in low dimensions. Thus the mean-field approximation is not necessary for deriving these prefactors. We present an exact, first-principles derivation for one, leaving the other as an open question. Our results suggest that mean-field theories of critical phenomena may compute more for d≥d_{u} than has been previously appreciated.

Contact network changes in ordered and disordered disk packings

We investigate the mechanical response of packings of purely repulsive, frictionless disks to quasistatic deformations. The deformations include simple shear strain at constant packing fraction and at constant pressure, "polydispersity" strain (in which we change the particle size distribution) at constant packing fraction and at constant pressure, and isotropic compression. For each deformation, we show that there are two classes of changes in the interparticle contact networks: jump changes and point changes. Jump changes occur when a contact network becomes mechanically unstable, particles "rearrange", and the potential energy (when the strain is applied at constant packing fraction) or enthalpy (when the strain is applied at constant pressure) and all derivatives are discontinuous. During point changes, a single contact is either added to or removed from the contact network. For repulsive linear spring interactions, second- and higher-order derivatives of the potential energy/enthalpy are discontinuous at a point change, while for Hertzian interactions, third- and higher-order derivatives of the potential energy/enthalpy are discontinuous. We illustrate the importance of point changes by studying the transition from a hexagonal crystal to a disordered crystal induced by applying polydispersity strain. During this transition, the system only undergoes point changes, with no jump changes. We emphasize that one must understand point changes, as well as jump changes, to predict the mechanical properties of jammed packings.

Shear jamming, discontinuous shear thickening, and fragile states in dry granular materials under oscillatory shear

We numerically study the linear response of two-dimensional frictional granular materials under oscillatory shear. The storage modulus G^{'} and the loss modulus G^{''} in the zero strain rate limit depend on the initial strain amplitude of the oscillatory shear before measurement. The shear jammed state (satisfying G^{'}>0) can be observed at an amplitude greater than a critical initial strain amplitude. The fragile state is defined by the emergence of liquid-like and solid-like states depending on the form of the initial shear. In this state, the observed G^{'} after the reduction of the strain amplitude depends on the phase of the external shear strain. The loss modulus G^{''} exhibits a discontinuous jump corresponding to discontinuous shear thickening in the fragile state.

Intrinsic dissipation mechanisms in metallic glass resonators

Micro- and nanoresonators have important applications including sensing, navigation, and biochemical detection. Their performance is quantified using the quality factor Q, which gives the ratio of the energy stored to the energy dissipated per cycle. Metallic glasses are a promising material class for micro- and nanoscale resonators since they are amorphous and can be fabricated precisely into complex shapes on these length scales. To understand the intrinsic dissipation mechanisms that ultimately limit large Q-values in metallic glasses, we perform molecular dynamics simulations to model metallic glass resonators subjected to bending vibrations at low temperatures. We calculate the power spectrum of the kinetic energy, redistribution of energy from the fundamental mode of vibration, and Q vs the kinetic energy per atom K of the excitation. In the harmonic and anharmonic response regimes where there are no atomic rearrangements, we find that Q → ∞ over the time periods we consider (since we do not consider coupling to the environment). We identify a characteristic K_r above which atomic rearrangements occur, and there is significant energy leakage from the fundamental mode to higher frequencies, causing finite Q. Thus, K_r is a critical parameter determining resonator performance. We show that K_r decreases as a power-law, K_r ∼ N^-k, with increasing system size N, where k ≈ 1.3. We estimate the critical strain ⟨γ_r⟩∼ 10^-8 for micrometer-sized resonators below which atomic rearrangements do not occur in the millikelvin temperature range, and thus, large Q-values can be obtained when they are operated below γ_r. We also find that K_r for amorphous resonators is comparable to that for resonators with crystalline order.