Search

Disordered solids respond to quasistatic shear with intermittent avalanches of plastic activity, an example of the crackling noise observed in many nonequilibrium critical systems. The temporal power spectrum of activity within disordered solids consists of three distinct domains: a novel power-law rise with frequency at low frequencies indicating anticorrelation, white-noise at intermediate frequencies, and a power-law decay at high frequencies. As the strain rate increases, the white-noise regime shrinks and ultimately disappears as the finite strain rate restricts the maximum size of an avalanche. A new strain-rate- and system-size-dependent scaling theory is derived for power spectra in both the quasistatic and finite-strain-rate regimes. This theory is validated using data from overdamped two- and three-dimensional molecular dynamics simulations. We identify important exponents in the yielding transition including the dynamic exponent $z$ which relates the size of an avalanche to its duration, the fractal dimension of avalanches, and the exponent characterizing the divergence in correlations with strain rate. Results are related to temporal correlations within a single avalanche and between multiple avalanches.

Rate-effects in sheared disordered solids are studied using molecular dynamics simulations of binary Lennard-Jones glasses in two and three dimensions. In the quasistatic (QS) regime, systems exhibit critical behavior: the magnitudes of avalanches are power-law distributed with a maximum cutoff that diverges with increasing system size $L$. With increasing rate, systems move away from the critical yielding point and the average flow stress rises as a power of the strain rate with exponent $1/\beta$, the Herschel-Bulkley exponent. Finite-size scaling collapses of the stress are used to measure $\beta$ as well as the exponent $\nu$ which characterizes the divergence of the correlation length. The stress and kinetic energy per particle experience fluctuations with strain that scale as $L^{-d/2}$. As the largest avalanche in a system scales as $L^\alpha$, this implies $\alpha < d/2$. The diffusion rate of particles diverges as a power of decreasing rate before saturating in the QS regime. A scaling theory for the diffusion is derived using the QS avalanche rate distribution and generalized to the finite strain rate regime. This theory is used to collapse curves for different system sizes and confirm $\beta/\nu$.

The packing and flow of aspherical frictional particles are studied using discrete element simulations. Particles are superballs with shape $|x|^{s}+|y|^{s}+|z|^{s} = 1$ that varies from sphere ($s=2$) to cube ($s=\infty$), constructed with an overlapping-sphere model. Both packing fraction, $\phi$, and coordination number, $z$, decrease monotonically with microscopic friction $\mu$, for all shapes. However, this decrease is more dramatic for larger $s$ due to a reduction in the fraction of face-face contacts with increasing friction. For flowing grains, the dynamic friction $\tilde{\mu}$ - the ratio of shear to normal stresses - depends on shape, microscopic friction and inertial number $I.$ For all shapes, $\tilde{\mu}$ grows from its quasi-static value $\tilde{\mu}_0$ as $(\tilde{\mu}-\tilde{\mu}_0) = dI^\alpha,$ with different universal behavior for frictional and frictionless shapes. For frictionless shapes the exponent $\alpha \approx 0.5$ and prefactor $d \approx 5\tilde{\mu}_0$ while for frictional shapes $\alpha \approx 1$ and $d$ varies only slightly. The results highlight that the flow exponents are universal and are consistent for all the shapes simulated here., Comment: 6 pages, 5 figures

While nearly all theoretical and computational studies of entangled polymer melts have focused on uniform samples, polymer synthesis routes always result in some dispersity, albeit narrow, of distribution of molecular weights (D_M=M_w/M_n ~ 1.02-1.04). Here the effects of dispersity on chain mobility are studied for entangled, disperse melts using a coarse-grained model for polyethylene. Polymer melts with chain lengths set to follow a Schulz-Zimm distribution for the same average M_w = 36 kg/mol with D_M = 1.0 to 1.16, were studied for times of $600-800$ $\mu$s using molecular dynamics simulations. This time frame is longer than the time required to reach the diffusive regime. We find that dispersity in this range does not affect the entanglement time or tube diameter. However, while there is negligible difference in the average mobility of chains for the uniform distribution D_M=1.0 and D_M = 1.02, the shortest chains move significantly faster than the longest ones offering a constraint release pathway for the melts for larger D_M.

The nested sampling algorithm has been shown to be a general method for calculating the pressure-temperature-composition phase diagrams of materials. While the previous implementation used single-particle Monte Carlo moves, these are inefficient for condensed systems with general interactions where single-particle moves cannot be evaluated faster than the energy of the whole system. Here we enhance the method by using all-particle moves: either Galilean Monte Carlo or a total enthalpy Hamiltonian Monte Carlo algorithm, introduced in this paper. We show that these algorithms enable the determination of phase transition temperatures with equivalent accuracy to the previous method at $1/N$ of the cost for an $N$-particle system with general interactions, or at equal cost when single particle moves can be done in $1/N$ of the cost of a full $N$-particle energy evaluation., Comment: 15 pages, 11 figures, 10 supplementary material files

Coupled length and time scales determine the dynamic behavior of polymers and underlie their unique viscoelastic properties. To resolve the long-time dynamics it is imperative to determine which time and length scales must be correctly modeled. Here we probe the degree of coarse graining required to simultaneously retain significant atomistic details and access large length and time scales. The degree of coarse graining in turn sets the minimum length scale instrumental in defining polymer properties and dynamics. Using linear polyethylene as a model system, we probe how coarse graining scale affects the measured dynamics. Iterative Boltzmann inversion is used to derive coarse-grained potentials with 2-6 methylene groups per coarse-grained bead from a fully atomistic melt simulation. We show that atomistic detail is critical to capturing large scale dynamics. Using these models we simulate polyethylene melts for times over 500 {\mu}s to study the viscoelastic properties of well-entangled polymer chains., Comment: 6 pages, 6 figures

The structure and interactions of coated silica nanoparticles have been studied in water using molecular dynamics simulations. For 5 nm diameter amorphous silica nanoparticles we studied the effects of varying the chain length and grafting density of polyethylene oxide (PEO) on the nanoparticle coating's shape and on nanoparticle-nanoparticle effective forces. For short ligands of length $n=6$ and $n=20$ repeat units, the coatings are radially symmetric while for longer chains ($n=100$) the coatings are highly anisotropic. This anisotropy appears to be governed primarily by chain length, with coverage playing a secondary role. For the largest chain lengths considered, the strongly anisotropic shape makes fitting to a simple radial force model impossible. For shorter ligands, where the coatings are isotropic, we found that the force between pairs of nanoparticles is purely repulsive and can be fit to the form $(R/2r_\text{core}-1)^{-b}$ where $R$ is the separation between the center of the nanoparticles, $r_\text{core}$ is the radius of the silica core, and $b$ is measured to be between 2.3 and 4.1., Comment: 20 pages, 6 figures

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 thousands to millions of 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 damage 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., Comment: 17 pages, 13 figures

Simulations are used to determine the effect of inertia on athermal shear of a two-dimensional binary Lennard-Jones glass. In the quasistatic limit, shear occurs through a series of rapid avalanches. The distribution of avalanches is analyzed using finite-size scaling with thousands to millions of particles. Inertia takes the system to a new underdamped universality class rather than driving the system away from criticality. Scaling exponents are determined for the underdamped and overdamped limits and a critical damping that separates the two regimes. Systems are in the overdamped universality class even when most vibrational modes are underdamped., Comment: 5 pages, 3 figures