No Arabic abstract
We analyze the fluctuations in particle positions and inter-particle forces in disordered jammed crystals 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.
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 $|delta r| sim 1/r$ and $|delta f| sim 1/r^2$ decay for the displacement fields and excess forces at large distances $r$ away from the defect. Finally we derive non-linear 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.
The properties of crystals consisting of several components can be widely tuned. Often solid solutions are produced, where substitutional or interstitional disorder determines the crystal thermodynamic and mechanical properties. The chemical and structural disorder impedes the study of the elasticity of such solid solutions, since standard procedures like potential expansions cannot be applied. We present a generalization of a density-functional based approach recently developed for one-component crystals to multi-component crystals. It yields expressions for the elastic constants valid in solid solutions with arbitrary amounts of point defects and up to the melting temperature. Further, both acoustic and optical phonon eigenfrequencies can be computed in linear response from the equilibrium particle densities and established classical density functionals. As a proof of principle, dispersion relations are computed for two different binary crystals: A random fcc crystal as an example for a substitutional, and a disordered sodium chloride structure as an example of an interstitial solid solution. In cases where one of the components couples only weakly to the others, the dispersion relations develop characteristic signatures. The acoustic branches become flat in much of the first Brillouin zone, and a crossover between acoustic and optic branches takes place at a wavelength which can far exceed the lattice spacing.*
We derive exact results for displacement fields that develop as a response to external pinning forces in two dimensional athermal networks. For a triangular lattice arrangement of particles interacting through soft potentials, we develop a Greens function formalism which we use to derive exact results for displacement fields produced by localized external forces. We show that in the continuum limit the displacement fields decay as $1/r$ at large distances $r$ away from a force dipole. Finally, we extend our formulation to study correlations in the displacement fields produced by the external pinning forces. We show that uncorrelated pinned forces at each vertex give rise to long-range correlations in displacements in athermal systems, with a non-trivial system size dependence. We verify our predictions with numerical simulations of athermal networks in two dimensions.
Athermal systems across a large range of length scales, ranging from foams and granular bead packings to crumpled metallic sheets, exhibit slow stress relaxation when compressed. Experimentally they show a non-monotonic stress response when decompressed somewhat after an initial compression, i.e. under a two-step, Kovacs-like protocol. It turns out that from this response one can tell the age of the system, suggesting an interpretation as a memory effect. In this work we use a model of an athermal jammed solid, specifically a binary mixture of soft harmonic spheres, to explore this phenomenon through in-silico experiments. Using extensive simulations under conditions analogous to those in experiment, we observe identical phenomenology in the stress response under a two--step protocol. Our model system also recovers the behaviour under a more recently studied three-step protocol, which consists of a compression followed by a decompression and then a final compression. We show that the observed response in both two-step and three-step protocols can be understood using Linear Response Theory. In particular, a linear scaling with age for the two-step protocol arises generically for slow linear responses with power law or logarithmic decay and does not in itself point to any underlying aging dynamics.
We develop an elasto-plastic description for the transient dynamics prior to steady flow of athermally yielding materials. Our mean-field model not only reproduces the experimentally observed non-linear time dependence of the shear-rate response to an external shear-stress, but also allows for the determination of the different physical processes involved in the onset of the re-acceleration phase after the initial critical slowing down and a distinct well defined fluidization phase. The evidenced power-law dependence of the fluidization time on the distance of the applied to an age dependent static yield stress is not universal but strongly dependent on initial conditions.