HTML conversions sometimes display errors due to content that did not convert correctly from the source. This paper uses the following packages that are not yet supported by the HTML conversion tool. Feedback on these issues are not necessary; they are known and are being worked on.
failed: xpatch
failed: xpatch
Authors: achieve the best HTML results from your LaTeX submissions by following these best practices.
Dynamic Screening by Plasticity in Amorphous Solids
H. George E. Hentschel
Dept. of Physics, Emory University, Atlanta Ga. 30322
Anna Pomyalov
Dept. of Chemical Physics, The Weizmann Institute of Science, Rehovot 76100, Israel
Itamar Procaccia
Dept. of Chemical Physics, The Weizmann Institute of Science, Rehovot 76100, Israel
Sino-Europe Complex Science Center, School of Mathematics, North University of China, Shanxi, Taiyuan 030051, China
Oran Szachter
Racah Institute of Physics, The Hebrew University of Jerusalem, Jerusalem, Israel 9190
(December 20, 2023)
Abstract
In recent work it was shown that elasticity theory can break down in amorphous solids subjected to nonuniform static loads. The elastic fields are screened by geometric dipoles; these stem from gradients of the quadrupole field associated with plastic responses. Here we study the dynamical responses induced by oscillatory loads. The required modification to classical elasticity is described. Exact solutions for the displacement field in circular geometry are presented, demonstrating that dipole screening results in essential departures from the expected predictions of classical elasticity theory. Numerical simulations are conducted to validate the theoretical predictions and to delineate their range of validity.
I Introduction
The response of amorphous solids to oscillatory strain is a well studied subject, mainly in protocols employing simple shear, cf. for example [1, 2, 3]. While important and indicating a lot of interesting physics, systems under oscillatory shear do not succumb easily to analytic scrutiny, making the comparison of analytic predictions to experiments or simulations quite difficult to accomplish. In this work we turn to oscillatory forcing by another protocol, of an oscillatory inflation of an inner circular boundary, creating non-uniform oscillatory strain on an amorphous solids that is contained between this and a much larger outer circular boundary that is maintained stationary.
The reason for this choice is two-fold. First, in a series of recent works it was shown that classical elasticity theory needs to be reconsidered and amended to describe correctly the mechanical response of amorphous solids subjected to time-independent non-uniform strains [4, 5, 6, 7, 8, 9, 10]. A new phase is formed in which the elastic fields are being screened by emergent geometric dipoles in the resulting displacement field. The transition from a quadrupolar screening in the solid phase to a dipolar screening is reminiscent of the structural transition in 2D crystals from hexagonal solids to hexatics, except that here the transition is in the emergent displacement field and the structure of the solid is always amorphous [9].
Our theory provides a classical field theory for describing the mechanical state of a deformed amorphous solids, and it predicts anomalous behavior that is observed in both numerical simulations and experiments.
In this paper we extend the theory to a new direction, to describe mechanical responses of amorphous solids to dynamical loading.
The second reason for choosing the present geometry is that it allows considerable analytic progress. We consider the effects of inertia and plastic responses in amorphous matter. For the sake of obtaining analytic solutions we will focus on an amorphous solid contained in an annulus of outer radius and inner radius , such that the inner radius oscillates with a fixed frequency . Initially, before oscillations begin, the material will be brought to mechanical equilibrium in which the resultant force on each particle vanishes. Once oscillations start, we will be interested in the time-dependent displacement field, evaluated with respect to the initial equilibrated positions. We will show that classical elasticity fails to predict correctly this displacement field, and that an appropriate theory requires taking into account the screening introduced by quadrupolar and dipolar charges that form due to plastic responses.
The structure of the paper is as follows: in Sect. II we review the dynamics as expected for our configuration from the solution of the equations of motion dictated by classical elasticity. We then solve analytically for the displacement field of a purely elastic medium which is subjected to oscillatory inflation of an inner boundary. We find that the dynamics is already non-trivial, exhibiting interesting features. Section III introduces quadruple and then dipole screening, leading to predictions of a rich array of expected solutions in which new length scales emerge spontaneously, breaking down elasticity theory. In Sect. IV we describe simulation results to test the prediction of the theory. The comparison of theory to simulations calls for some careful considerations. First, it is important that the dynamics will describe oscillations around a well defined mechanical equilibrium state. Second, we need to deal with the issue of dissipation. In the numerics we introduce dissipation by damping terms in the particle collisions. In the theory there are dissipative terms proportional the rate of change of displacement fields. This requires careful discussion. Finally, nonlinear effects should be kept at bay. When all these are considered, we find indeed dynamical responses that are in accord with the novel theory.
Section V offers conclusion and a discussion of the road ahead.
II Dynamics of amorphous solids in two dimensions
II.1 The displacement field in purely elastic solids
We start by reviewing the classical approach to dynamics and dissipation á la Landau & Lifshitz [11], and then we develop the basic ideas that are called for to accommodate the physics of amorphous solids.
In a purely elastic medium of mass density , denote the stress field by , and the displacement field by . The equation of equilibrium is given by
(1)
i.e. the force per unit volume is zero in equilibrium . Once we are not in mechanical equilibrium the force does not vanish but is equal to the acceleration times the mass per unit volume
(2)
For an isotropic body, one can write in terms of the displacement to get [11]
(3)
where and are Lamé coefficients.
II.1.1 Dynamical response to oscillations of an inner circle
Elasticity theory holds equally well in two and in three dimensions. The theory described below can be easily extended to three dimension, and see for example [8]. In this paper the simulations are performed in two dimensions for computational efficiency, and therefore
we specialize the equations to the case of an oscillating inner circle of initial radius and time dependent radius , in a system bounded by a fixed, rigid outer circular boundary of radius . The simulations were done according to the following protocol. At the beginning of each periodic driving cycle, the inner circle was inflated to , and the system was equilibrated. Then we periodically inflate and deflate the central disk according to
(4)
where are parameters. After a given number of oscillations, the system is again allowed to equilibrate. Then the procedure is repeated for ten times. The reason for this initial set up is to guarantee that the subsequent oscillations are around a well defined equilibrated state and the resulting dynamics is reproducible. It is important that the oscillating displacement field is measured with respect to a well-defined equilibrated state.
In order to compare the numerical responses to the offered theory, we calculate and analyze the angle-averaged radial displacement.
Clearly, in a realistic amorphous solid such an oscillatory inflation exerts work, and without dissipation the system will not reach a stationary state. Thus in the numerical simulations presented below, cf. Sect. IV, we will add a small dissipative term to the dynamics of the disks that comprise the amorphous solids. This results in the system reaching an oscillatory steady state.
At this point we continue to focus on an ideal elastic medium to solve Eq. (3) as it stands, without adding dissipation that in an ideal elastic solid does not exist.
In the presence of circular symmetry we can focus on the radial component of the displacement field, i.e
(5)
The equation for becomes
(6)
This equation of motion is solved with the boundary conditions
(7)
The displacement field is measured from the equilibrated configuration obtained after the expansion of the inner boundary to . Accordingly we seek a solution of the form
(8)
where
(9)
The dynamic response solves the Bessel equation (9). Interestingly, this equation is of the same form as that found in anomalous static anomalous elasticity cf. Ref. [4]. Here the frequency related expression acts as dynamical screening. The solution of Eq. 9 together with boundary conditions (II.1.1) in terms of the first order
Bessel and Neumann functions is
(10)
where
An important feature of this solution is that is fully determined by the geometry of the system, the prescribed frequency, i.e. by and the mechanical properties of the media and . In particular, the number of nodes of the oscillating Bessel functions is determined by the prescribed frequency . In order to achieve presentation of the results which does not depend on the system size and material properties, we rewrite the equation and its solutions in dimensionless variables and . This is done using a characteristic length the outer radius and characteristic time where is the dilatational speed
Note that in this result , , , and the frequency is dimensionless.
To get familiarity with the type of functions involved we show in Fig. 1 a few examples of as given by Eq. (13). Later we will learn how dissipation and dipole screening change this function in realistic situations.
II.1.2 Adding dissipation
In reality, in any experimental or model system, dissipation will be
crucial to balance the work done by an oscillating pulsar at the center of the circular system. In the simulations, we employ disks that interact via Hertzian repulsive forces (see Appendix A), but we introduce dissipation by adding to the particle dynamics a dumping force proportional to their velocity , , which is applied to all the disks, to remove any excess energy. In the macroscopic modeling we need therefore to take into account an effective viscous term, that is added to Eq.(3):
(14)
At this point the function is kept unspecified, later it will be chosen to accommodate the results of the numerical simulations. For any choice of the appearance of will now mix sine and cosine functions. In dimensionless units we seek a solution for the radial component in the form
(15)
Next we solve for the functions and by substitution the ansatz Eq. (15) into the boundary value equation Eq. (14), matching together terms in and . We find
(16)
These two coupled real equations can be solved by introducing one complex variable
The solutions of Eq. (19) are a combination of a first order
Bessel and Neumann functions of the complex variable .
Thus we can write
(20)
To determine the coefficients and we need to fit the boundary conditions. They are
at the outer boundary. At the oscillating inner boundary , . Accordingly,
, and .
Using these boundary conditions, we define
These exact solutions are expected to
be relevant as long as dipole screening is absent. To assess the effects of the latter we now add quadrupole and dipole
contributions to the theoretical discussion.
III Anomalous Dynamics due to plastic Screening
In this section we derive the equations of motion in the presence of quadrupolar and dipolar screening. This derivation follows the ideas presented in Ref. [4] for the static problem, supplemented with inertial and dissipative contributions as required. In each case we need to write down the appropriate
Lagrangian which reflects the interactions that are taken into account. We denote the quadrupolar tensor field as . Physically, this is the field associated with the eigenstrains of the Eshelby quadrupoles that are formed by plastic responses. The dipolar field is simply the gradient of this density. The displacement field is the vector field and the strain tensor
. The Euclidean metric is denoted .
In all cases we start with the
Euler-Lagrange equations in the standard form [12], with playing
the role of a fundamental field, which in our case can be the displacement or the quadrupole density.
(23)
III.1 Quadrupole Screening
To derive the equations of motion under quadrupole screening we assume that the density
of quadrupoles is sufficiently low, or that the gradients of their density are negligible.
Thus to quadratic order we can write the Lagrangian for an elastic medium with quadrupoles as
(24)
Computing the derivative with respect to the quadrupole field we find
(25)
At this point we do not have dipoles in the Lagrangian, or
Here we denoted the inverse of as and defined a new (presently unknown) tensor of coefficients .
Eq. (27) is important, showing that the quadrupolar tensor is not some mysterious entity, but that it is induced by the strain field in the system, with the quartic tensor providing the link.
At this point we can define a normalized stress field and write
(30)
where the last equality follows from Eq. (27) and homogeneity and isotropy, with being a scalar number. The upshot of this calculation is that the dynamics of the displacement field is unchanged
in form compared to Eq. (2), except for a renormalization of the elastic moduli. Thus in the present case the introduction of the dissipative term and the analytic solutions for the displacement field follow verbatim the theory presented in the previous section.
III.2 Dipole Screening
The situation changes qualitatively with dipole screening, the equations of motion change their form. To derive these equations we will assert that the renormalization due to quadrupole screening is already included in our Lagrangian, in the form of a renormalized stress tensor (that will be again denoted as ). To quadratic order in the fields we therefore write
(31)
Computing the derivative with respect to the dipole field we find
(32)
The expression in the square brackets is a constant, that can be taken as zero using the translational invariance of the displacement field. From this point onward we lose the gauge freedom of the displacement field since we chose a gauge. As before, we use this equation to express the dipole field in terms of the fundamental displacement field:
(33)
Computing the derivatives with respect to the displacement field and its derivative we find
(34)
(35)
Combining these two equations together we write
(36)
Inverting for we find the modified equation of motion
(37)
In a homogeneous and isotropic medium we find
(38)
We note that Eq. (37) breaks translational symmetry together with the introduction of
a length scale, since has the dimension of an inverse scale. This is further discussed at
great length (and excuse the pun) in Sect. V below.
III.3 Solutions of Anomalous Dynamic dipole screening
As done before, we write in terms of the displacement, to get the explicit equation in the case of dipole screening. As before, we add the dissipative term to the resulting equation and end up with
(39)
We note that the addition of the screening term does not
change the linearity of the equation, and therefore we can use the solution presented in Subsect. II.1.2 with very little modification. Referring again to a solution in the form of Eq. (15), to solve for and we only need to change the variable to a new variable where
(40)
where .
The Eqs. (21) and (22) are unchanged except that the everywhere is replaced by .
To provide a feeling to the nature of the solutions of Eq. (39) we show in Figs. 2 and 3 the functions and for
various frequencies and different values of the screening parameter .
In all cases the dissipation function with and . The choice of this value of has been made purely by comparison to simulations. But one can rationalize it by noting that the geometric spreading of 2-dimensional waves originating from a point source is expected to decay like [13]. On the other hand, since the number of collisions between disks is increasing like , in balance one can expect that a dissipation that increases like should suffice to keep the energy budget finite.
IV Numerical Simulations
IV.1 Setup
The numerical setup consists of a 2D concentric circular enclosure
with an inner circle of radius , filled with a bi-disperse mixture of frictionless disks. The disks interact through normal Herzian forces only [see Appendix A]. The interaction with the outer wall is also Hertzian, with the interaction coefficient equal to the interaction coefficient between disks. Additional viscous dumping force is applied to all discs to remove excess of energy from the system. Here is the microscopic damping coefficient (to be distinguished from the macroscopic of Eq. (14)), and is the velocity of the ’th disc.
The simulations are carried out using LAMMPS [14] for the dynamics of the granular system. The driving is realized by periodic inflation of the radius of the inner enclosure wall which is placed at the center of the system.
The granular system was prepared at a desired packing fraction and pressure, starting with the
creation of a randomly distributed set of particles, half of which has a radius
and and the other half with , with . Next the system is equilibrated by energy minimization. Then the radius of the enclosing circle is decreased by small decrements from the initial radius to a value corresponding to the largest desired packing fraction according to . At each step, the system is again equilibrated. A set of configurations for various pressures , or equivalently packing fractions , was obtained.
Having a system at mechanical equilibrium we can start oscillating the radius of the inner circle. The procedure followed was explained in Subsect. II.1.1. As stated there, to see reproducible results it
is important to guarantee that the oscillations take place around a true
stationary state, in the sense that one is guaranteed that stopping the
oscillations at any point and performing energy minimization would result in the same mechanical equilibrium independent of when and where the oscillatory driving is stopped.
Another lesson from the numerical experiments is that for the present system of Hertzian disks one is limited in the frequencies of driving that can be simulated. Increasing the frequency too much results in the creation of a hole around the oscillating boundary - the system does not have enough time to return to contact with the boundary when this boundary recedes back too rapidly. We therefore limited the simulations, and the comparison with the theory, to relatively low frequencies as we will see next. It is likely that in systems in which there are attractive forces between the constituents the creation of the hole can be eliminated, and higher frequencies could be studied. On the other hand, it is possible that for any system higher frequencies would trigger nonlinear interactions, requiring additional scrutiny of the equations of motion. At this point we leave these interesting questions to future study.
IV.2 Results for the dominant radial mode
To compare the theory with the numerical experiments, we note that the measured functions are calculated by angle-averaging over a band of finite width and therefore we amend the boundary conditions at the inner boundary as
In Fig. 4 we present comparisons between the measured functions and and the prediction of the theory with dissipation and screening for a number of driving frequencies. The parameters and that were used to fit the analytical solution defined by equations (15),(20),(40), and (IV.2), are typical for our system. The quality of the fits is typical, as long as the frequency of oscillations is not too high, we find very consistent agreement between the analytical theory and the simulations. As can be seen in Appendix B, the transverse components of the displacement field (as well as the nonlinear contributions) are much smaller than the radial ones,
and thus we do not refer to these components in this paper.
The dependence of the inverse scale and the dissipative parameter on the frequency of oscillations appears linear. In Fig. 5 we plot an ensemble-averaged values (over ten independent configuration ) of the and . The observed linearly is a reflection of the linear regime of the system’s dynamics. It is likely that for higher frequencies, with nonlinear effects becoming relevant, also this observed linearity will be lost.
Our results indicate that as the frequency decreases, goes to zero. This means that the low frequency limit does not coincide with a single inflation of the inner circle [4]. The training of our system by repeated oscillation until a stable steady state is obtained removes the anomalous responses to a single inflation step. The system becomes quasi-elastic and not static-anomalous.
V Conclusions and the road ahead
It is important to stress that consistent and reproducible results for the measured displacement field under oscillatory forcing depend crucially on insisting upon the existence of true stationary state around which the oscillations take place. As explained, the requirement is that stopping the dynamics at any point in time, and applying energy minimization algorithms result in the same equilibrated configuration up to some severe tolerance requirements. This requires repeated training of the system until such a stable steady state configuration is attained. In terms of the energy picture such a stable steady state corresponds to a sufficiently deep local minimum around which the oscillations take place.
A second requirement for the successful correspondence between the analytic theory and the simulations is a synchronization between the oscillatory driving and the system’s dynamics. In the present case, when the driving frequency exceeds in dimensionless units, the system reaction cannot follow the driving, entering a dynamical regime that cannot be treated with the present linear theory. Needless to say, this regime is of interest, but is beyond the scope of the present paper.
Finally, the results of this study indicate very strongly that linear classical elasticity fails to provide a valid description of the response of amorphous solids to oscillatory driving, in much the same way as it fails to describe the response to static driving as has been shown in recent work. It is therefore worthwhile to study in the near future also the effect of screening on oscillatory driving by simple or pure shear, since these modes of forcing have attracted a lot of recent interest due to the presence of shear banding and fracture. Such a future direction will automatically necessitate a nonlinear extension of the screening theory, an effort that is under present active research.
Appendix A Hertzian interaction
The disks interact via normal Hertzian forces
(43)
Here is the overlap distance of two discs, is the elastic constant of the normal contact, is the viscoelactic constant for normal contact, and is the effective mass.
In our simulations all , .
Additional viscous dumping force is applied to all discs to remove excess of energy from the system. Here is the damping coefficient and is the velocity of the disc. The simulations
are carried out using LAMMPS [14] for the dynamics of the granular system using LJ units.
Appendix B The dominance of the driving mode over nonlinear and azimuthal contributions
Here we present the evidence that the driving mode is dominant in the considered range of parameters. Since the original displacement field has both radial and azimuthal components defined as . In principle we can solve our equations for both components, but it turns out that the azimuthal component is small compared to the radial one and we did not consider it for the present conditions. In Fig. 6 we compare the radial contributions at the largest driving frequency in our range with the azimuthal contribution at the same frequency, and . This justifies the present neglect of this component. In addition, we confirm that the radial components of the lowest order nonlinear contribution with double frequency, and , are also negligible at our range of frequencies, as is shown in Fig. 7.
Acknowledgements.
We thank Michael Moshe for useful discussions at the early stages of this project. This work has been supported by the joint grant between the Israel Science Foundation and the National Science Foundation of China, and by the Minerva Foundation, Munich, Germany.
References
[1]
Davide Fiocco, Giuseppe Foffi, and Srikanth Sastry.
Oscillatory athermal quasistatic deformation of a model glass.
Phys. Rev. E, 88:020301, Aug 2013.
[2]
Nikolai V. Priezjev.
Shear band formation in amorphous materials under oscillatory shear
deformation.
Metals, 10(3), 2020.
[3]
Ido Regev, Ido Attia, Karin Dahmen, Srikanth Sastry, and Muhittin Mungan.
Topology of the energy landscape of sheared amorphous solids and the
irreversibility transition.
Phys. Rev. E, 103:062614, Jun 2021.
[4]
Anaël Lemaître, Chandana Mondal, Michael Moshe, Itamar Procaccia, Saikat
Roy, and Keren Screiber-Re’em.
Anomalous elasticity and plastic screening in amorphous solids.
Phys. Rev. E, 104:024904, Aug 2021.
[5]
Bhanu Prasad Bhowmik, Michael Moshe, and Itamar Procaccia.
Direct measurement of dipoles in anomalous elasticity of amorphous
solids.
Phys. Rev. E, 105:L043001, Apr 2022.
[6]
Avanish Kumar, Michael Moshe, Itamar Procaccia, and Murari Singh.
Anomalous elasticity in classical glass formers.
Phys. Rev. E, 106:015001, Jul 2022.
[7]
Chandana Mondal, Michael Moshe, Itamar Procaccia, Saikat Roy, Jin Shang, and
Jie Zhang.
Experimental and numerical verification of anomalous screening theory
in granular matter.
Chaos, Solitons and Fractals, 164:112609, 2022.
[8]
Harish Charan, Michael Moshe, and Itamar Procaccia.
Anomalous elasticity and emergent dipole screening in
three-dimensional amorphous solids.
Phys. Rev. E, 107:055005, May 2023.
[9]
Yuliang Jin, Itamar Procaccia, and Tuhin Samanta.
An intermediate phase between jammed and un-jammed amorphous solids,
2023.
[10]
Chandana Mondal, Michael Moshe, Itamar Procaccia, and Saikat Roy.
Dipole screening in pure shear strain protocols of amorphous solids,
2023.
[11]
Lev Davidovich Landau and Eugin M Lifshitz.
Course of Theoretical Physics Vol 7: Theory of Elasticity.
Pergamon press, 1959.
[12]
Sean M Carroll.
Spacetime and geometry.
Cambridge University Press, 2019.
[13]
Sebastiano Foti, Carlo G Lai, Glenn J Rix, and Claudio Strobbia.
Surface wave methods for near-surface site characterization.
CRC press, 2014.
[14]
A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown,
P. S. Crozier, P. J. in ’t Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen,
R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton.
LAMMPS - a flexible simulation tool for particle-based materials
modeling at the atomic, meso, and continuum scales.
Comp. Phys. Comm., 271:108171, 2022.