Defect Modes in One-Diensional Granular Crystals N. Boechler, G. Theocharis, Y. Man,2, P. G. Kevrekidis 3, and C. Daraio Graduate Aerospace Laboratories (GALCIT), California Institute of Technology, Pasadena, CA 925, USA 2 Keble College, University of Oxford, Oxford OX 3PG, UK 3 Departent of Matheatics and Statistics, University of Massachusetts, Aherst, MA 3-455, USA arxiv:3.2483v [cond-at.other] 2 Mar 2 We study the vibrational spectra of one-diensional statically copressed granular crystals (arrays of elastic particles in contact) containing defects. We focus on the prototypical settings of one or two spherical defects (particles of saller radii) interspersed in a chain of larger unifor spherical particles. We easure the near-linear frequency spectru within the spatial vicinity of the defects, and identify the frequencies of the localized defect odes. We copare the experientally deterined frequencies with those obtained by nuerical eigen-analysis and by analytical expressions based on few-site considerations. We also present a brief nuerical and experiental exaple of the nonlinear generalization of a single-defect localized ode. PACS nubers: 63.2.Ry, 63.2.Pw, 45.7.-n, 46.4.-f INTRODUCTION Defect odes in crystals have long been studied in the real of solid state physics [, 2]. The presence of defects or disorder is known to enable localized lattice vibrations, whose associated frequencies have been easured in the spectra of real crystals (see [, 3] and references therein). More recently, this study has been extended to include other exaples, including superconductors [4] and electron-phonon interactions [5]. Siilar phenoena have also been observed in nonlinear systes, including photonic crystals [6], optical waveguide arrays [7 9], dielectric superlattices (with ebedded defect layers) [], icroechanical cantilever arrays [], and Bose- Einstein condensates of atoic vapors [2, 3]. Granular crystals are nonlinear systes coposed of densely-packed particles interacting through Hertzian contacts [4 7]. These systes present a rearkable ability to tune their dynaic response fro linear to strongly nonlinear regies [5]. This has allowed the exploration of fundaental nonlinear wavefors such as traveling waves[5 8] and discrete breathers[9]. Granular crystals have also been proposed for several engineering applications, such as energy absorbing layers [2 23], actuating devices [24], and sound scrablers [25, 26]. The presence of defects in statically uncopressed (or weakly copressed, as copared to the relative dynaic displaceents) granular chains excited by ipulsive loading has been studied in a nuber of previous works that have reported the existence of interesting dynaic responses such as the fragentation of waves, anoalous reflections, and energy trapping [2 23, 27 32]. In this paper, we study the response of strongly copressed granular crystals, with one or two defects (extending our earlier theoretical work [33]), excited by continuous signals. We easure the frequency response of the syste and reveal localized odes due to the presence of defects. We report that the nuber of localized odes irrors that of the defects, and note that the frequencies of such odes depend on: (i) the ratio of the defect ass to the ass of the particles in the unifor chain, (ii) the relative proxiity of ultiple defects, (iii) the geoetric and aterial properties of the particles coposing the crystal, and (iv) the static load. We copare our experiental findings with nuerical coputations and with theoretical analysis approxiating the behavior of a few sites in the vicinity of the defect(s). Finally, we deonstrate that as we go fro the linear to the nonlinear regie, nonlinear deforations of the linear defect odes (with appropriately downshifted frequencies) are sustained by the syste. EXPERIMENTAL SETUP We asseble D granular crystals, siilar to those described in [9, 34], coposed of N = 2 statically copressed stainless steel spherical particles (36 type, with elastic odulus E = 93 GPa and Poisson ratio ν b =.3 [35]), as shown in Fig. (a). The chain is coposed of unifor particles of (easured) radius R = 9.53 and ass M = 28.84 g, except for one (or two) light-ass stainless steel defect particles. The spheres are held in a D configuration using four polycarbonate bars(2.7 diaeter) that are aligned by polycarbonate guide plates spaced at approxiately 2 c intervals along the axis of the crystal. The defect particles, which are of saller radii than the rest of the particles of the chain, are aligned with the axis of the crystal using polycarbonate support rings. Dynaic perturbations are applied to the chain by a piezoelectric actuator ounted on a steel cube (which acts as a rigid wall). The particles are statically copressed by a load of F = 2 N. The static load is applied using a soft spring (of stiffness.24 kn/), which is copressed between the last particle in the chain and a second steel cube bolted to the table. The applied static load is easured by a calibrated load cell placed between the spring and the steel cube.
2 We easure the dynaic force signals of the propagating waves with custo-ade force sensors consisting of a piezo-electric disk ebedded inside two halves of a stainless steel particle with radius R = 9.53. The sensorparticles are carefully constructed to reseble the ass, shape, and contact properties of the other spherical particlescoposingthe restofthe crystal[2, 25,26, 37,38]. displaceent obtained fro a static load F is n = (F /A n ) 2/3 [4, 5], where R n is the radius of the nth particle. In order to study the linear spectru of the inhoogeneous granular crystal, we linearize Eq. () about the equilibriu state under the presence of the static load. This yields the following linear syste [9, 33, 36]: a) Piezoelectric actuator Large steel spheres Polycarbonate guide rails Dynaic force sensor 2 Sall steel defect sphere and polycarbonate alignent ring Soft spring Wall Static load cell n ü n = K n ( ) K n+ ( + ), (2) where K n = 3 2 A2/3 n F /3. Following [33], we siplify Eq.(2) to the eigensyste: ω 2 Mu = Λu, (3) b) Transfer Function 5 4 6 8 2 4 FIG. : [Color online] a) Scheatic diagra of the experiental setup for the hoogeneous chain with a single defect configuration. b) Experiental transfer functions (as defined in the single-defect: near linear regie section) for a granular crystal with a static load of F = 2 N and a defect-bead of ass = 5.73 g located at site n def = 2. Blue (dark-grey) [red (light-grey)] curves corresponds to transfer function obtainedfro theforcesignal ofasensor particleplacedatn = 4 [n = 2]. The diaond arker is the defect ode. The triangle arker is the upper acoustic cutoff ode. The vertical black dashed line is the theoretically predicted defect ode frequency, and the vertical solid black line is the theoretically predicted upper acoustic cutoff frequency. THEORETICAL MODEL We consider the D inhoogeneous crystal of N beads as a chain of nonlinear oscillators [5]: n ü n = A n [ n + ] p + A n+ [ n+ + + ] p +, () where[y] + denotesthe positivepartofy (which signifies that adjacent particles interact only when they are in contact), is the displaceent of the nth sphere (where n {,,N}) around the static equilibriu, n is the ass of the nth particle, and the coefficients A n depend on the exponent p and the geoetry/aterial properties of adjacent beads. The exponent p = 3/2 represents the Hertz law potential between adjacent spheres [4]. In this case, A n = 2E 3( ν 2 ) ( Rn R n R n +R n ) /2, and the static where M is a N N diagonal atrix with eleents M nn = n, and u is the displaceent vector. Λ is a N N triagonal atrix with eleents Λ n = [K n + ( δ nn )K n+ ]δ n +K n+ (δ n +δ n+ ),whereδ isthe Kronecker delta and we consider left-fixed and right-free boundary conditions. The right-free boundary assuption derives fro the low stiffness of the static copression spring (Fig. (a)) as copared to the stiffness of the particles in contact. SINGLE-DEFECT: NEAR-LINEAR REGIME In this section, we study D granular crystals that are hoogeneous except for one light-ass defect bead at site n def, as shown in Fig. (a). Solving the eigenvalue proble of Eq. (3), for such a granular crystal, we obtain the eigenfrequencies and the corresponding spatial profile of the odes of the syste. The presence of the single light-ass defect generates a localized ode (see also [32, 33]), centered at the defect site, which we will refer to as the defect ode. The defect ode aplitude decays exponentially away fro the defect site and its 4K RR M frequency f d is such that f d > f c, where f c = is the upper cutoff frequency of the acoustic band of the hoogeneous host crystal (where K RR = 3 2 A2/3 RR F/3 is the linear stiffness of the contact between two beads with radius R). The spatial profile of this ode consists of adjacent particles oscillating out of phase (see inset in Fig. (2)). As the radius of the defect bead becoes saller, the difference between f d and f c becoes larger, while the defect ode becoes ore spatially localized. We observe that for the granular crystals studied here, with r radii ratios of R <.7, the defect ode involves the otion of up to approxiately three beads, i.e., the displaceents of the beads at n n def +2 and n n def 2 arenegligible. Becausein thisrangeofradiiratiothe otion of the particles can be accurately approxiated by three beads, we consider the particles at n = n def ± 2 as fixed walls, in order to find an analytical approxiation for the frequency of the defect ode. Solving for
3 the eigenfrequencies of this reduced three-bead syste, we find that the ode corresponding to the out of phase otion can be analytically approxiated by Eq. (4) f 3bead = 2K Rr M +K RR +K Rr + 8K Rr K RR M +[2K Rr M +(K RR +K Rr )] 2 2M (4) where K Rr = 3 2 A2/3 Rr F/3 is the linear stiffness of the contact between a defect-bead and a bead of radius R. We conduct experients to identify the frequency of the defect ode in granular crystals with a single lightass defect as shown in Fig. (a). We place the defect particle at site n def = 2 (close to the actuator) so that the energy applied by the actuator, at the defect ode frequency, will not be copletely attenuated by the unifor crystal, which acts as a echanical frequency filter, before it arrives at the defect site. Because of the localized nature of the defect ode, placing a defect particle (of radius r 7.4 ) at site n def = 2 or further into the chain akes nearly no difference on the frequency of the defect ode. For instance, for a defect particle of radius r = 7.4, we nuerically calculate (using Eq. 3) the difference in the defect ode frequency for the cases where a defect particle is placed at site n def = 2 or n def =, to be 3 Hz. Conversely, because of the presence of the fixed boundary and the larger localization length of the defect ode, for a defect particle of r = 8.73, we calculate the difference in defect ode frequency, between sites n def = 2 and n def =, to be 68 Hz. The defect particles are stainless steel spheres of saller radii, r = [3.97, 4.76, 5.56, 6.35, 7.4, 7.94], and easured asses of = [2.8, 3.6, 5.73, 8.54, 2.9, 6.65] g, respectively. We experientally characterize the linear spectru of this syste by applying low aplitude (approxiately 2 N) bandwidth liited noise(3 25 khz for the two sallest defect particles, and 3 5 khz otherwise) via the piezoelectric actuator. We calculate the transfer functions, specific to the sensor location, by averaging the Power Spectral Densities(PSD[39]) of 6 forcetie histories, easured with the ebedded sensors, and dividing by the average PSD level in the 3 8 khz range (corresponding to the acoustic band). We ebed sensors in particles at sites n = 4 and n = 2. In Fig. (b) we show the transfer functions for the granular crystal with defect radius r = 5.56. The red (light-grey) and blue (dark-grey) curves are the transfer functions for the sensors at sites n = 4 and n = 2, respectively. We denote the experiental cutoff frequency by the triangular arker (found by identifying the last peak in the acoustic band) and defect frequency as the diaond arker on the n = 4 transfer function. The vertical lines denote the theoretically deterined upper cutoff frequency of the acoustic band and the defect frequency (Eq. 4). The presence of the defect ode can be clearly identified in the vicinity of the defect (at n = 4), but is not visible far fro the defect (at n = 2). We repeat the process of easuring the transfer function and identifying the defect ode frequency 6 ties, re-assebling the crystal after each repetition. In Fig. 2, we plot the average frequency of the 6 experientally identified defect odes as a function of the ass ratio M (blue [dark-grey] solid line connecting the closed diaonds). We also plot, for coparison, the defect frequency predicted by the analytical expression of Eq. (4) (green [light-grey] dashed line connecting the crosses), and the nuerical eigenanalysis of Eq. (3) corresponding to the experiental setup (black solid line connecting the open diaonds). The error bars on the experiental data are ±2σ where σ is the standard deviation of the identified defect frequencies over the 6 repetitions. Coparing the analytical three-bead approxiation with the nuerical eigenfrequencies, we find an excellent agree- ent for ass ratios of M <.6. Coparing the experiental data with the nuerics, we find an upshift of 5 %, siilar to the upshift observed in [9, 34]. For the r = 5.56 defect, the average experiental defect frequency is f exp d = 3.59 khz and the average experiental cutoff frequency is fc exp = 8.36 khz. In coparison, the theoretical 3-bead approxiation gives a defect frequency of fd 3bead = 2.84 khz and the eigenproble of Eq. (3) gives a defect frequency of fd nu = 2.85 khz, while the analytically calculated cutoff frequency was f c = 8.2 khz. Possible reasons for these upshifts have been identified in [9, 34] and the references therein, such as error in the aterial paraeters, nonlinear elasticity, surface roughness, dissipative echaniss and isalignent of the particles. We note that a systeatic error in the easureent of the static load could also cause such an upshift. Nevertheless, it is clear fro Fig. 2 that the functional dependence of the relevant frequencies on the ass ratio (of defect to regular beads) is accurately captured by our analytical and nuerical results.
4 22 2 8 6 4 2 8 (noralized).5 /M=.2 5 5 2 Bead Nuber.2.4.6.8 /M FIG. 2: [Color online] Frequency of the defect ode, with defect-bead placed at n def = 2, as a function of ass ratio /M. Solid blue line (dark-grey, closed diaonds) corresponds to experients, solid black line (open diaonds) to nuerically obtained eigenfrequencies (see Eq. (3)), and green dashed line (light-grey, x arkers) to the analytical prediction of the three-beads approxiation (see Eq. (4)). The error bars account for statistical errors on the easured frequencies and are ±2σ. Inset: The noralized defect ode for M =.2. TWO-DEFECTS: NEAR-LINEAR REGIME We study granular crystal configurations with two identical light-ass defects to better understand the effects of increasing heterogeneity on the spectral response of the syste. The localized ode due to the presence of a light-ass defect (for ass ratios M <.6) has a spatial localization length of about three particles (larger particles have a greater localization length, and saller particles have a shorter localization length), as described in the one-defect case and shown in Fig. 2. We can thus expect that two light-ass defects placed far fro each other in a granular crystal (sufficiently outside this localization length) would have siilar frequencies and ode shapes independent of the presence of the other. However, as the two defect particles are brought closer together (within the localization length), each ode influences the other. For a sufficiently sall ass ratio, this results in the creation of two defect odes at different frequencies; one with the defect particles oving out of phase, and the other with the defect particles oving in phase. For the case of nearestneighbor identical defects, our theoretical analysis can be extended by using a four-particle analogy. In this case, using the notation s = K Rr (M + ) + K RR, s 2 = 4K Rr K RR M+(K RR +K Rr (M +)) 2, s 3 = s + 2K rr M and s 4 = 4(2K rr K RR + K Rr (2K rr + K RR ))M + (2K rr M + K RR + K Rr (M + )) 2, we obtain the following frequencies f () 4bead = f (2) 4bead = 2M (s ± s 2 ), (5) 2M (s 3 ± s 4 ). (6) The two highest frequencies correspond to the linear defect ode frequencies. Naturally, this analytical approach can be extended to ore distant defects, although we do not present such algebraically intensive cases here. In Fig. 3, we show the behavior of two r = 5.56 defects in a N = 2 particle granular crystal under F = 2Nstaticload(siiliartotheconfigurationshown infig.(a)), wherethefirstdefectisatsitek = n def = 2 and the second defect is at a variable position between site l = n def2 = 3 and l = n def2 = 6. We use the sae experiental ethod as in the single defect case except now we use a noise range between 3 2 khz and we place the first sensor at n = l +. We show the experientally deterined PSD transfer function for the case of l k = in Fig. 3(a), with sensors at site n = 4 (blue [dark-grey]) and n = 2 (red [light-grey]). As described in [33], the existence of two separate defect odes for the case where the defect particles are adjacent to each other (l k = ), depends on the ass ratio of the defect particles to those of the rest of the crystal. Here the ass ratio is such that two odes are present, as can be seen in the blue (dark-grey) curve in Fig. 3(a). The two distinct odes, which we denote by the open square and closed circular arkers, have frequencies above the acoustic band. The square arkers denote the ode with defect beads oving out of phase, and the closed circular arker corresponds to the ode with defect particles oving in phase, as shown by the nuerically calculated eigenodes in Fig. 3(c) and(d) respectively[33]. In Fig. 3(b) we plot the experientally deterined frequencies of both odes as a function of the inter-defect particle distance (l k). The solid blue (dark-grey) lines are the experiental data, and the dashed black lines are the frequencies obtained fro solving the eigenvalue proble of Eq. (3). The green (light-grey) x-arkers denote the frequencies calculated with Eqs. (5)-(6), for the l k = case. It is evident that the analytical results agree closely with the nuerically calculated eigen-frequencies. The error bars on the experiental data correspond to the ±2σ standard deviation as calculated in the single-defect case. We see close qualitative agreeent between the experiental data and the nuerical predictions, but also the sae systeatic upshift as observed in the single defect case and [9, 34]. Fro Fig. 3(b) we can see that as the defects are placed three or ore particles apart, the frequencies of the defect odes converged to approxiately the sae value, suggesting the defects respond independently of each other.
5 Transfer Function (nor.) 5 5 5 l k= a) c) 5 5 2 Bead Nuber (nor.) 6 4 2 b) 2 3 4 Defect Distance (beads) d) l k= 5 5 2 Bead Nuber FIG. 3: [Color online] (a) Experiental transfer functions for a granular crystal with two defect-beads of ass ratio M =.2 at n def = 2 and n def = 3 (in contact). Blue (darkgrey) [red (light-grey)] curve corresponds to transfer function obtained fro the force signal of a custo sensor placed at n = 4 [n = 2]. (b) Frequencies of the defect odes as a function of the distance between the. The solid line denotes experiental data, the dashed line the nuerically obtained eigenfrequncies, and the x arkers the frequencies fro the analytical expresssions of Eqs. (5)- (6). (c),(d) The noralized defect ode shapes corresponding to the defect odes identified in (a) with frequency of the sae arker-type. SINGLE DEFECT: NONLINEAR LOCALIZED MODES As shown in [33], the interplay of the inherent nonlinearity of the granular crystal with the linear localization due to the defect results in the presence of robust nonlinear localized odes (NLMs). The frequency of these odes depends not only on the static load and the aterial values of the beads but also on the aplitude of the oscillations. In order to find this dependence, we apply Netwon s ethod (see [33] and references therein) for the experiental, single-defect, configuration of Fig. (a). For the nuerical calculations in this section, we calculate experiental contact coefficients following a procedure siilar to the one described in [9]. The experiental contact coefficients obtained are A exp RR =.79 N/µ3/2 for the contact between two R = 9.53 beads and A exp rr = 9.95 N/ µ3/2 for the contact between the R = 9.53 and the r = 5.56 beads. In coparison, the values of the coefficient A, as calculated by the aterial values and used for the previous sections of the paper, are A RR = 9.76 N/µ 3/2 and A Rr = 8.38 N/ µ 3/2. In Fig. 4(a), we show the frequency of the nuerically deterined NLM as a function of the averaged dynaic forcefortheparticleatsiten = 3. Thelattercorresponds to the average of the two dynaic contact forces adjacent to the particle, which is analogous to what is easured experientally by the dynaic force sensor [25]. In Fig. 4(b), we plot the nuerically deterined noralized NLM shape at f b = 3.28 khz. Coparing this NLM shape to the linear ode shape of the sae frequency (see inset of Fig. 2), we can see that the NLM has a slightly odified (ore asyetric) spatial profile. TheexperientalsetupusedforthestudyoftheNLMs is the sae asin the case ofthe linear singledefect experients (as shown in Fig. (a)) except we place sensors in particles at sites n = 3, n = 5, and n = 2. Additionally, we replace the n = particle with an ebedded actuator particle, so as to apply high aplitude (approxiately N), short tie pulse (approxiately µs) perturbations directly to the defect particle. Exciting such a pulse creates an initial condition in the crystal that resebles the predicted defect NLM shape. The ebedded actuator particle is siilar in construction to the sensors but with a piezoceraic construction/geoetry ore appropriate for high force aplitude actuation (Piezoechanik PCh 5/5x5/2 Piezo-chip). The force-tie history of the dynaic force easured by the sensor at site n = 3 is shown in Fig. 4(c). A sharp excitation is evident at tie t =, followed by periodic oscillations with a decaying envelope, due to the inherent dissipation in the syste. As shown by the paraetric continuation in Fig. 4(a), NLMs corresponding to the defect ode at higher aplitudes have a frequency deeper into the gap than its linear counterpart. However, for the aplitudes observed here this is only a slight shift (up to 2 Hz over 7 N). We study, in ore detail, the response of two selected tie regions of the force-tie history shown in Fig. 4(c), to experientally deonstrate the frequency shift characteristic of higher aplitude NLMs. The two nonoverlapping tie regions are of length T = 5. s. The red (light-grey) tie region begins iediately following the arrival of the initial actuated pulse, and presents a axiu aplitude of 7 N. The blue (dark-grey) tie region starts T = 6 s after the begining of the previous tie region, and presents a axiu aplitude of.3 N. We calculate PSDs for both tie regions (frequency resolution δf = 95 Hz) as shown in Fig. 4(d). The PSDs shown in Fig. 4(d) correspond to the tie regions of the sae color shown in Fig. 4(c). Here, the PSDs are noralized by dividing the PSD by the peak PSD aplitude of the identified defect ode. It is evident that the peak in the PSD spectru corresponding to the tie region with larger force aplitude presents a lower characteristic frequency (i.e., it is further into the gap) with respect to the peak representing the tie region with lower force aplitudes. This is in agreeent with the shift predicted by the paraetric continuation analysis shown in Fig. 4(a). The peak frequency of the PSD of the high force aplitude tie region is
6 f def = 3.28 khz, and the peak frequency of the PSD of the low aplitude tie region is f def = 3.48 khz, where f def = 3.48 khz is closer to the ean experientally deterined linear defect ode frequency (shown by the dashed line in Fig. 4(d)). Force (N) 3.6 3.5 3.4 5 Force (N) a) c) 2 Tie (s) (noralized) Nor. PSD.5 d) 3 5 5 2 Bead Nuber 2 4 FIG. 4: [Color online] (a) Nuerical frequency continuation of the nonlinear defect odes corresponding to the experiental setup in Fig. (a). (b) Nuerically calculated spatial profile of the nonlinear localized ode with frequency f def = 3.28 khz. (c) Measured force-tie history of a sensor at site n = 3, where a high aplitude, short width, force pulse is applied to the granular crystal. (d) Noralized PSD for the easured tie regions of the sae color in (c); closed and open diaonds correspond to the high and low aplitude tie regions respectively. The vertical dashed line is the ean experientally deterined linear defect ode frequency. CONCLUSIONS We studied the response of statically copressed granular crystals containing light-ass defects, and characterized their near-linear spectra by applying continuous excitation. We deonstrated that such chains support localized odes with frequencies above that of their acoustic band cutoff, using approxiate few-bead analytical calculations, nuerics, and experients. The nuber of supported localized odes depends on the nuber of defects, while their frequencies depend on the inter-defect distance, on the ratio M of defect to regular asses (and the geoetric/elastic properties of the beads), and on the static load. We also briefly described the nonlinear generalizations of such odes, departing fro the nearlinear regie, and showed a downshift of the corresponding defect ode frequencies with increasing aplitude. This study is iportant for understanding the interplay of disorder and nonlinearity in discrete systes, and the results reported ay be relevant in the design of applications involving vibrational energy trapping. b) We thank Stéphane Job for help with the experiental setup. GT and PGK acknowledge support fro the A.S. Onassis Public Benefit Foundation through RZG 3/2-2 and PGK also through NSF-CMMI- 337. CD acknowledges support fro NSF-CMMI- 96954 and NSF-CMMI-84454 (CAREER). [] A. A. Maradudin, E. W. Montroll, and G. H. Weiss, Theory of Lattice Dynaics in the Haronic Approxiation, Acadeic, New York (963). [2] I. M. Lifschitz, Nuovo Ciento, Suppl. 3, 76 (956); I. M. Lifschitz and A. M. Kosevich, Rep. Progr. Phys. 29, 27 (966). [3] G. Lucovsky, M.H. Brodsky and E. Burstein, Phys. Rev. B 2, 3295 (97). [4] A. F. Andreev, JETP Lett. 46, 584 (987); A. V. Balatsky, Nature (London) 43 77 (2). [5] M. I. Molina and G. P. Tsironis, Phys. Rev. B 47, 533 (993); G. P. Tsironis, M. I. Molina, and D. Hennig, Phys. Rev. E 5, 2365 (994). [6] S. Y. Jin et al., Science 282, 274 (998); M. G. Khazhinsky and A. R. McGurn, Phys. Lett. A 237, 75 (998). [7] U. Peschel et al., Appl. Phys. Lett., 75, 348 (999). [8] R. Morandotti et al., Opt. Lett. 28, 834 (23). [9] E. Sirnov et al., Opt. Express 4, 248 (26). [] E. Lidorikis et al., Phys. Rev. B 56, 59 (997). [] M. Sato et al., Europhys. Lett. 66, 38 (24). [2] P. Engels and C. Atherton, Phys. Rev. Lett. 99, 645 (27). [3] T.W. Neely et al., Phys. Rev. Lett. 4, 64 (2). [4] K. L. Johnson, Contact Mechanics, (Cabridge University Press, 985). [5] V. F. Nesterenko, Dynaics of Heterogeneous Materials (Springer-Verlag, New York, NY, 2). [6] S. Sen et al., Phys. Rep. 462, 2 (28). [7] C. Coste, E. Falcon, and S. Fauve, Phys. Rev. E. 56, 64 (997). [8] M. A. Porter et al., Phys. Rev. E 77, 56(R) (28). [9] N. Boechler et al., Phys. Rev. Lett. 4, 24432 (2); see also Phys. Rev. E 82, 5664 (2). [2] C. Daraio et al., Phys. Rev. Lett. 96, 582 (26). [2] J. Hong, Phys. Rev. Lett. 94, 8 (25). [22] F. Fraternali, M. A. Porter, and C. Daraio, Mech. Adv. Mat. Struct. 7, (2). [23] R. Doney and S. Sen, Phys. Rev. Lett. 97, 5552 (26). [24] D. Khatri, C. Daraio, and P. Rizzo, SPIE 6934, 6934U (28). [25] C. Daraio et al., Phys. Rev. E 72, 663 (25). [26] V. F. Nesterenko et al., Phys. Rev. Lett. 95, 5872 (25). [27] E. Hascoet and H. J. Herann, Eur. Phys. J. B 4, 83 (2). [28] E. J. Hinch and S. Saint-Jean, Proc. R. Soc. Lond. A 455, 32 (999) [29] J. B. Hong and A. G. Xu, App. Phys. Lett. 8, 4868 (22). [3] S. Sen, M. Manciu, and J. D. Wright, Phys. Rev. E 57, 2386 (998). [3] M. Manciu, S. Sen, and A. J. Hurd, Physica A 274, 588 (999); M. Manciu, S. Sen, and A. J. Hurd, Physica A
7 274, 67 (999). [32] S. Job et al., Phys. Rev. E 8, 2562(R) (29). [33] G. Theocharis et al., Phys. Rev. E 8, 666 (29). [34] N. Boechler et al., arxiv:.586, J. Appl. Phys. in press (2). [35] Metals Handbook, th Ed., ASM International, Materials Park, OH (99). [36] E. B. Herbold et al., Acta Mechanica 25, 85 (29). [37] S. Job et al., Phys. Rev. Lett. 94, 782 (25). [38] V. F. Nesterenko, A. N. Lazaridi, and E. B. Sibiryakov, Prikl. Mekh. Tekh. Fiz. 2, 9 (995). [39] A.V. Oppenhei and R.W. Schafer, Discrete-Tie Signal Processing (Prentice-Hall, 989).