Designing Phononic Band Gaps with Sticky Potentials

Spectral gaps in the vibrational modes of disordered solids are key design elements in the synthesis and control of phononic metamaterials that exhibit a plethora of novel elastic and mechanical properties. However, reliably producing these gaps often require a high degree of network specificity through complex control optimization procedures. In this work, we present as an additional tool to the existing repertoire, a numerical scheme that rapidly generates sizeable spectral gaps in absence of any fine tuning of the network structure or elastic parameters. These gaps occur even in disordered polydisperse systems consisting of relatively few particles ($N \sim 10^2-10^3$). Our proposed procedure exploits sticky potentials that have recently been shown to suppress the formation of soft modes, thus effectively recovering the linear elastic regime where band structures appear, at much shorter length scales than in conventional models of disordered solids. Our approach is relevant to design and realization of gapped spectra in a variety of physical setups ranging from colloidal suspensions to 3D-printed elastic networks.

Spectral gaps in the vibrational modes of disordered solids are key design elements in the synthesis and control of phononic metamaterials that exhibit a plethora of novel elastic and mechanical properties. However, reliably producing these gaps often require a high degree of network specificity through complex control optimization procedures. In this work, we present as an additional tool to the existing repertoire, a numerical scheme that rapidly generates sizeable spectral gaps in absence of any fine tuning of the network structure or elastic parameters. These gaps occur even in disordered polydisperse systems consisting of relatively few particles (N ∼ 10 2 − 10 3 ). Our proposed procedure exploits sticky potentials that have recently been shown to suppress the formation of soft modes, thus effectively recovering the linear elastic regime where band structures appear, at much shorter length scales than in conventional models of disordered solids. Our approach is relevant to design and realization of gapped spectra in a variety of physical setups ranging from colloidal suspensions to 3D-printed elastic networks.

I. INTRODUCTION
Low energy excitations of vibrational modes in disordered solids have been a subject of recent interest, for it presents an exciting approach to the construction of tunable metamaterials with a wide range of applications. From its collective contributions to bulk acoustic or sound attenuation properties [1], to the remarkable precision in conformational changes made available with the advent of directed mechanical response [2][3][4], these soft modes have offered an unprecedented level of control over the design of modern metamaterials due to its remarkable localization and topological properties [5][6][7][8]. In particular, recent phononic applications originating from these concepts include the design of auxetic materials with negative possion ratio [9][10][11], allosteric inspired nano-or macro-scale levers capable of eliciting a targeted mechanical response at a distance [12,13], and the experimental realization of materials with negative effective stiffness [14][15][16].
To robustly implement these novel features in a sustained and controllable manner requires the existence of spectral band gaps that isolate the specific desired functionality of low energy excitations from the remainder of its spectrum that constitute noise to the intended response. As such, understanding and controlling the properties of spectral gaps forms an important aspect to metamaterial design that goes beyond simply optimizing for the displacement field of the low frequency eigenmodes. This concern for the spectral architecture has motivated studies on the formation and tunability of band gaps in disordered photonic systems, that have recently also been incorporating hyperuniform point distributions that display peculiar behavior in Fourier space and their resulting structure factor [17][18][19]. Previous studies addressing phononic or acoustic band gaps in discrete elastic systems have largely focused on ordered or quasi-periodic systems [20][21][22][23], with far less attention paid to spectral gaps of disordered solids [24].
In this work, we present a numerical scheme capable of rapidly generating phononic spectral band gaps that dictate the robustness of programmable mechanical behavior in disordered elastic networks. Unlike existing strategies that depend on optimizing linear response [24] or some secondary cost functions [12,13], this scheme robustly generates spectral gaps with no fine tuning of the elastic and topological properties. As such, our approach serves not only as a standalone algorithm, but can also be used to generate precursors used in combination with the various existing optimization techniques to accentuate particular functionalities, while having the assurance of sizeable spectral gaps at the vicinity of a pre-determined frequency.
Phononic band gaps are naturally present in all systems as long as their size is larger than their characteristic disorder length scale. In this limit, the vibrational spectrum exhibits well defined phonon gaps that characterizes linear elasticity. However, for conventional systems of particles interacting via purely repulsive [25] or Lennard-Jones(LJ)-like [42] potentials, this disorder length scale is large. For instance, in a study of two-dimensional LJ systems, convergence to the elastic limit was only observed for systems comprising millions of particles [26]. From a design perspective, devising disordered and gapped mass-spring networks of a small manageable size that one might conceivably construct for example via 3d printing techniques, thus becomes the key issue of concern. Here, we introduce an approach to do so.
This approach in essence involves first obtaining energy minimal configurations of pairwise interacting particles mediated through a family of sticky potentials that are short range repulsive, but attractive at intermediate distances of separation [27,28]. Recent work suggest that such interactions result in the drastic reduction in the availability of soft vibrational modes in its jammed state [29,30], and thus these networks exhibits gapped spectra reminiscent to that of linear elasticity at much shorter lengths scale compared to their WCA or LJ-like counterparts, despite being inherently disordered from a structural perspective [29]. The arrived energy minimal configuration can then be mapped onto a harmonically interacting linear elastic network that serves as the blue print for construction of the target phononic metamaterial.
In the following, we first introduce in Sec. (II) the general class of sticky potentials and energy minimization protocols specific to this work, before briefly discussing elastic properties of the networks derived from energy minimal configurations. In Sec. (III), we quantify the occurrence statistics of the spectral gap and their corresponding frequencies, showing that they remain robust and highly predictable over a range of (III A) densities and (III B) parameters of the potential. Moreover, we show that these spectral gaps remain pronounced even in small systems (III C), which are often of particular interest to practical realizations. In Sec. (IV), we demonstrate how removal of the pre-stress contribution can further enhance performance of the spectral gaps. Lastly we conclude our work by summarizing in a schematic (V), the step-wise procedures involved in obtaining robust gapped spectra in linear elastic networks.

II. STICKY POTENTIALS AND ELASTIC NETWORKS
The key towards reliable generation of phononic band gaps in small disordered systems lies in the emergent collective properties of particles interacting via a class of "sticky" potentials [27]. In this section, we review the functional form of this potential that has been adapted specifically for this work, and briefly discuss its various properties.
Traditional atomistic models of disordered solids have largely focused on the purely repulsive Weeks-Chandler-Andersen (WCA) [25] potential, or on the Lennard-Jones (LJ) one [42] which possesses a slowly decaying attractive tail ∝ 1/r 6 , that captures qualitative behavior of physical Van-der-Waals interactions. Sticky potentials seek to vary this effective range of the attractive tail at intermediate length scales. This seemingly modest modification to the local potential has far reaching consequences to the dynamical behavior of its relaxation dynamics [27] and the elastic properties [31] of its jammed configuration at sufficiently high densities [28] and is also known to result in systems with negative thermal expansion [32].
More importantly, the energy minimal state in solids or otherwise inherent state [33], where pressure (p) and shear modulus (µ) remains positive definite, also exhibits phononic spectra {ω n } and density of states D(ω) that differ significantly from the typical behavior of jammed or LJ-like disordered solids [29,30]. As observed in all disordered solids [34][35][36][37][38], the vibrational density of states of these systems is the superposition of a quartic scaling behavior regime D(ω) ∼ A 4 ω 4 comprising quasilocalized modes, and of a Debye-like contribution, D(ω) ∼ A 4 ω d−1 in d spatial dimensions. However, quasi-localized soft modes in sticky potentials are severely suppressed [29,30], e.g. the parameter A 4 is unusually small. This work thus focuses on how one might exploit the natural suppression of these quasi-localized low frequency modes to produce robust gapped spectra that adheres closely to linear elasticity even in small disordered systems.
We note that several variants of these sticky potentials exists in the current literature [27][28][29][30][31], with specifics of its parameterization playing a non-trivial role in the behavior of physical observables [29,30]. These variabilities seem preliminarily related to the degree of "stickiness", an emergent measure of the relative attractive to repulsive forces present in the energy minimal state [29]. Here, instead of discussing the various trade-offs involved, we focus on a particular version [27] of sticky potentials, and suggest corresponding parameter values that reliably produce gapped spectra with high fidelity. In this potential (1), particles at pair separations r i j ≤ r min i j ≡ 2 1/6 σ i j interact repulsively via Kob-Andersen / Lennard-Joneslike interactions (LJ) [42] similar to purely repulsive WCA potentials [25] where σ i j ≡ (σ i + σ j )/2 is the effective length scale that characterizes polydispersity σ i of interacting particles. In addition to this contact repulsion at extremely short length scales, particles also attract each other at intermediate distances r min i j ≤ r i j ≤ r c , in which r c is a cut-off parameter where the potential vanishes. The remaining coefficients a 0 , a 1 and c 2l are then uniquely obtained by solving equations that ensure continuity in U(r) and its first ∂ r U(r) and second derivatives ∂ 2 r U(r) at distances r min i j and r c . Such a parameterization allows for the repulsive ( f (r) > 0) part of the interaction to be unchanged for different values of r c while varying the domain in r i j where interactions are attractive ( f (r) < 0) (see Fig. 1a). In the limit of r c → r min , U(r) becomes discontinuous at r min such that non-zero region of U(r) is purely repulsive and resembles qualitatively the WCA potential, while r c → ∞ recovers the long range interacting LJ potential. Previous numerical studies indicate that r c ≥ 2.4 essentially recovers the LJ limit in which the phononic spectra is known to be continuous [27]. As such, interesting parameter regimes of the sticky potentials are restricted to intermediate values of r c ∼ 1.2.
Unlike purely repulsive systems, the resulting bond tensions f ≡ −∂ r U and stiffnesses k ≡ ∂ 2 r U between particles derived from sticky potentials can assume negative values as seen in Fig. 1(b) and (c) respectively. We note that while r c → ∞ reduces the magnitude of extreme values in the local bond stiffnesses which may in principle ease its construction, the broadening of regime in r i j where k is negative as r c increases, renders such systems to be impractical especially for the purpose of constructing elastic networks, and thus may instead be detrimental to its design. Moreover, a relatively large r c as compared to the average inter-particle distance ρ −1/d results in extended interactions with other particles beyond its immediate neighbors. This further  For a given r c and number density ρ , the sticky potential can be used to deterministically generate particle configurations in arbitrary dimensions using a variety of energy minimization algorithms that ranges from gradient descent to monte carlo based methods [43,44]. In this work, we implement a conjugate gradient algorithm [45] to obtain particle configurations starting from initial states generated from a homogeneous Poisson point process, with all data presented unless otherwise stated, derived from systems of N = 2000 particles in two dimensions. σ i selected from a uniform distribution within [0.8, 1.2] such that the length scale of the system is given by σ i = 1 as conventionally adopted in previous studies of sticky potentials [27]. The system is highly polydisperse, the standard deviation of the σ distribution being 12%σ, and hence lack both translational and positional order. ε ≡ 1 sets the energy scale of the system, and all particles have mass m = 1.
Termination of the minimization is imposed upon mechanical equilibrium based on a threshold defined by magnitude of the largest net force acting on each particle -max i || j i f i j || ≤ 10 −12 .
The linear response regime of this energy minimal configuration is investigated by Taylor expanding to second order the energy of the system around the equilibrium configuration [? ]. In doing this operation, one maps the positions of particles in the energy minimal configurations { R i } as shown in Fig. 2(a) on to an equivalent elastic network Fig.(2(b) where { R i } now represent the position of nodes, and pairwise interactions (r i j ≤ r c ) are replaced by effective harmonic springs of stiffness k i j = ∂ r U(r i j ) and rest lengths l 0,i j = r i j + ( f i j /k i j ). The response of this effective elastic network informs on vibrational properties of the energy minimal configuration in the linear response regime, and thus in turn fully determine the phononic spectrum of the material. In Fig. 3, we show the statistical features of elastic properties of derived networks for an ensemble of N samples = 2500 configurations for various ρ evaluated at r c = 1.2. While potentials generically permit the existence of negative stiffness k < 0, a good choice of r c that restricts the domain at which this happens can significantly reduce the probability of encountering unconducive environments for a wide range of densities ( Fig.3(a)). Indeed, k < 0 is undesired for material design, as it is difficult to fabricate unit springs or linear elastic segments with negative stiffness. Moreover, we note that the distribution of k can exhibit bi-modalities at sufficiently high ρ where the average inter-particle length scale enters deep into the repulsive regime and is thus related to the existence of contacts with particles in the first and second coordination shell as mentioned above. On the other hand, the distribution of bond tensions (Fig.3b) remains uni-modal at all densities considered but attractive interactions ( f < 0) diminish with increasing ρ, indicating the reduction of stickiness driven by a change in density that was also observed in [29]. These properties collectively may provide additional avenues of tunability for multi-functionalities that depend on local stiffnesses of the material [13] that are no doubt interesting possibilities that lies beyond the scope of this work .

III. SPECTRAL PROPERTIES
Having obtained the elastic networks representing energy minimal configurations of sticky potentials in (II), we evaluate and discuss the novel spectral properties of sticky solids in this section. In general, the vibrational modes represent the elementary excitations of a (meta)stable configuration and can be obtained from the dynamical matrix or otherwise hessian of its potential: where α, β labels the indices of orthogonal axes and i, j indexes the particles. Upon diagonalization, the hessian yields normal modes |n with eigenfrequencies ω n = √ −λ n where λ n represents eigenvalues of the n-th mode.
For a typical disordered system interacting via repulsive WCA potential, the length scale in two dimensions and corresponding system sizes at which continuum elasticity sets in is tremendously large (≥ 4 × 10 4 ) [26] such that phononic excitations are very low in frequency and thus do not appear in significantly smaller systems (N ∼ 10 2 − 10 3 ). As such, the low frequency branch is dominated by quasi-localized soft modes which results in a continuous spectra in which no spectral gaps exists as we show in Fig.4, by plotting the average frequency ω n of a given mode number n calculated for systems of size N = 2000. Now, unlike their WCA or LJ counterparts, solids of sufficient stickiness suppresses the density of quasi-localized soft modes [29] such that the spectrum is not continuous and gapped as seen in Fig. 4. More importantly, the typical occurrence of these gaps and frequencies of the quasi-degenerate branches in {ω n }, all closely resemble predictions made based on the fundamental frequency of the transverse mode ω e suggested by linear elasticity, which we evaluate and plot as blue horizontal lines in Fig. 4 for non-trivial modes in the first two branches of the spectrum. From a design perspective, the question of interest is how to robustly generate and maximize the observed gaps.
To investigate this issue and formally characterize properties and assess robustness of these spectral gaps, we define a dimensionless characteristic gap frequency (ω i ) for the i-th linear elastic branch: that marks the onset (lower) frequency of the gap and a corresponding dimensionless gap parameter (γ i ) between the i and i + 1 linear elastic branch: where ∆ω n = ω n+1 − ω n and . . . Ω i denotes the average over modes in the i-th branch of the vibrational spectrum. Visual representation of these quantities are shown schematically in Fig.4. In the following subsections, we examine the properties and consistency of these spectral parameters derived for a "clean" ensemble consisting of post-selected samples that strongly adheres to linear elasticity under various conditions. The clean ensemble consists of samples that do not contain rattlers or bonds with negative stiffness, and also do not include samples where either of the gaps occurs at a mode number not in accordance to linear elasticity. This is done for the ease of consistency in physical realizations that do not require fine tuning. Specifically we investigate the ρ and r c dependence of this ensemble in III A and III B respectively, and show that these gaps can persists even in systems as small as N = 500 in section III C.
A. ρ dependence In Fig. 5(a,c) we show the distribution in γ i for various ρ with r c = 1.2 for the i = 1, 2 linear elastic branches respectively. We see that with increasing density, the profiles of gap distributions comprising of clean samples shift towards higher values, indicating that the non-dimensional features of these gaps are accentuated at higher frequencies. This observation is also supported by the largely increasing trend in the average of this distribution as shown in Fig. 6(b).
Furthermore, we note that the low frequency tail of P(γ i ) decreases with increasing ρ suggesting that the proportion of systems with γ i ∼ O(0) is vanishing at sufficiently high density. At the same time, the distributions of characteristic frequencies at which the gap occurs Fig. 5(b,d), and their respective mean values Fig. 6(c), also suggest that ω i typically increases with ρ, while the fluctuations about its mean values decreases, resulting in a strongly peaking distribution in ω i for high density configurations. For the intended purpose of robustly producing physically implementable spectral gaps in materials with high fidelity, a consistently larger γ i with smaller fluctuations in ω i is intuitively desired. This might suggest that higher parameter values of ρ may be conducive to the production of gapped systems. However we show in Fig.6(a) that the fraction of samples considered clean and thus constituting the ensemble corresponding to the observed distributions, is a non-monotonic function of the density. In fact, it dramatically decreases at higher densities ( f clean ∼ 0.4 for ρ = 1.3), which implies that the intuitive notion of preferring higher ρ may not be entirely productive in recovering linear elasticity in small systems. Indeed, with increasing densities, the typical inter-particle length scale ∼ ρ −1/d shifts towards the purely repulsive regime in U(r), which thus suggest behavior that rapidly approaches ungapped WCA driven systems with increasing density.
This perspective on the role of ρ in gap generation is further confirmed by studying the average fraction of sticky f s and super-sticky f ss particles across all samples generated, (i.e inclusive of "unclean" configurations) which we show in the inset of Fig.6(a). Sticky particles refer to nodes with at least one attractive interaction, while super sticky particles are those that only interact attractively (i.e no repulsive interaction). Note that no super sticky particles are present in our data for the two highest densities considered, indicating that f ss < 10 −6 for ρ ≥ 1.2. Now, we see that although the behavior of γ i and ω i is desired for higher ρ within the ensemble of clean samples, the emergent stickiness monotonically decreases, resulting in the presence of fewer clean samples. This nuance in the ρ dependence may be of concern to material design, and its relevance is dependent on whether post selection is part of the protocol in question. Densities within the range -0.9 ρ 1.2 are generally recommended.

B. r c dependence
Analogous to the analysis performed in the previous subsection that had focused on the role of density, we investigate here the dependence of gap properties on the parameter r c that controls the width of attractive interactions in the potential (1). In particular, we consider configurations at fixed density ρ = 1, and evaluate the probability distributions of γ i (Fig.7(a,b)) and ω i (Fig.7(c,d)) for r c = {1.15, 1.20, 1.25}. Corresponding ensemble averaged quantities are shown in Fig.8.
As with the scenario encountered with increasing ρ in sec. III A, the distribution profiles (Fig.7) and mean values of γ i (Fig.8) shifts in the positive direction with increasing r c , and this again is accompanied by a decrease in fluctuations of the characteristic frequencies at which the respective gaps occur 7(c,d). However, the limit r c → ∞ recovers qualitatively, the long range characteristics of the LJ potential that results in configurations that does not exhibit gapped spectra. To resolve this corundum , we point towards the dramatic reduction in the fraction of clean samples ( f clean ∼ 0.1) present in r c = 1.25 that we highlight is significantly lower than even its value in the highest density scenario considered in III A. This behavior stems from the broadening regime of U(r) that yields negative bond stiffness with increase in r c Fig. 1(c), that are discarded when constructing the clean ensemble. These results suggest that perhaps only a very narrow region of parameter values in the vicinity of r c ∼ 1.2 are ideally suitable for the robust generation phononic band gaps.
Furthermore, we note that unlike the scenario encountered when varying ρ the fraction of sticky and super sticky particles, do not seem to change significantly with r c even when f clean decreases dramatically (see inset of Fig.8(a)). Hence, direct association between the sample stickiness of individual configurations and their corresponding gap properties seem to be tenuous despite the critical nature of the attractive component in U(r) in recovering linear elasticity. In fact, we do not find any meaningful correlations between sample stickiness s as defined in [29], with the gap properties for all combinations of r c and ρ considered. With that said, the numerical values of the Pearson coefficients (C s,γ i ) indicate a weak positive correlation between sample stickiness and the gap size -0.047 < C s,γ 1 < 0.058 and 0.054 < C s,γ 2 < 0.110, that may be of some limited use in identifying configurations with a given gap size based on its stickiness.

C. System size dependence
We investigate the effects of system size on the behavior and statistical properties of spectral gaps. As mentioned earlier, the predictions of linear elasticity is made increasingly precise by an increase in system size of disordered materials [26]. The use of sticky potentials in this work act to suppress the density of soft modes, such that linear elastic behavior and the occurrence of gapped spectra sets in at significantly smaller length scales. Here we show in Fig.9, that these gaps persist and their statistical properties remain robust even in small systems for ρ = 1.2 and r c = 1.2, (i.e parameter values that have been previously determined to be in the ideal range for gapped spectra creation). While we see that the fluctuations in ω 1 decrease and typical values of γ 1 increase with larger system sizes as expected also for purely repulsive potentials [26], the spectral gap in systems as small as N = 500 remains well-defined and pronounced. This N is significantly smaller than the N ∼ 4 × 10 5 found by [26] needed to recover continuum elasticity in two dimensions. Lastly, we note that while the gap is unambiguously present even at N = 500, the degeneracy of frequencies within each linear elastic branch are increasingly lifted with smaller system sizes (see Fig. 9(a)). This may aid in selecting or exciting particular eigenmodes for practical applications so long as the broadening of bands do not destroy the primary features of the spectral gap.

IV. ELIMINATING PRE-STRESS CONTRIBUTIONS
In this section, we discuss and show how eliminating pre-stress contributions to the dynamical matrix (2) can be used to obtain stress-free normal modes that we also show serve to enhance spectral gap features of the elastic network. To begin with, we write down explicitly the off-diagonal terms i j of the hessian where k i j and l 0,i j are respectively the stiffness and rest lengths of the harmonic interaction between particles i and j that constitute the effective elastic network. Now, in the absence of pre-stress, the terms involving f i j , l 0,i j vanishes, and we are thus able to write down a stress-free hessianM that is completely characterized by the set of stiffnesses k i j : Upon diagonalizingM, we then obtain the vibrational spectrum of the stress-free network with corresponding eigenfrequencies ω n that now depends only on a set of stiffness parameters {k i j }. We note that existing approaches in designing allosteric metamaterial that tune for specific functionalities, also eliminate pre-stress contributions by construction, and in this regard is unlike natural occurring proteins it is trying to simulate [12]. However, this reduction in the number of local elastic parameters of the network have practical advantages when it comes to the actual synthesis of the material notwithstanding its positive effects on the spectrum that we discuss subsequently.
In Fig.10(a) we plot the ordered set of eigenfrequencies for elastic networks with (orange) and without pre-stress(green) contributions derived respectively from M andM originating from the same energy minimal configuration of a sticky solid at r c = 1.2 and ρ = 1. The frequencies are scaled by their respective gap frequency of the first linear elastic branch ω 1 and ω 1 , such that the enhancement to the first gap parameter is made visually clear. We note that this accentuation of linear elastic behavior upon elimination of pre-stress is not unique to sticky potentials, but a generic feature also present in conventional models [47].  Furthermore, in Fig.10(b) we see that the probability distribution of γ i is shifted in the positive direction for stress-free systems indicating that the enlargement in gap size is statistically significant. Fluctuations in the corresponding gap frequencies ω i is also observed to decrease in Fig.10(c). These results indicate that the elimination of pre-stress contributions in the energy minimal configuration, typically enhance the gap features in the derived elastic network. Here we reiterate for clarity that unlike WCA driven systems [12], these gaps are reliably present in sticky networks even without this explicit removal of pre-stress, and hence do not require any fine tuning to generate. The procedure detailed in this section serves as an optional amplification of the gaps.
We remark that "unclean" samples due to the appearance of bonds with negative stiffness previously discussed in section III A and III B can be made experimentally relevant by simply discarding bonds with negative stiffness (i.e replacing all k < 0 by k = 0) when considering the stress-free hessianM [7]. For all densities considered at r c = 1.2, we find that the typical deviation of γ i upon removing these bonds are less than 1% from their corresponding values in the stress-free system.
Lastly, we note that imprecision in printing or random variability involved in experimental realization of our approach may also result in deviations in the bond stiffness k used in practise. However, we find that random perturbations in the values of k do not significantly affect the properties of the spectral gaps. In particular, by imposing random perturbations selected from a uniform distribution δk ∼ [−kδ max , +kδ max ], we find that the typical deviation of γ i upon imposing these perturbations are less than 3% for δ max as large as 0.05.

V. CONCLUSION
In this work, we enumerated the spectral properties of disordered solids consisting of particles mediated by a class of sticky potentials and detailed its various parameter and size dependencies. For a wide range of parameters, we show that such systems adhere to continuum elasticity even in small systems and that sizeable spectral gaps are consistently and predictably present in absence of any fine tuning. In addition, we discuss appropriate parameter selection and various strategies including the elimination of pre-stress that collectively enhance gap features for practical applications. To succinctly summarize these findings and formalize our approach, we present our results as a step-wise algorithm in generating elastic networks with robust spectral gaps that is depicted by the schematic flow chart in Fig.11.
Lastly, we emphasize that while this approach in itself, reliably generates spectral gaps without fine tuning, the adaptation of these sticky potentials in constructing elastic networks is meant to be used in conjunction with the various exisiting protocols that optimize for specific functionalities in mechanical response [12,13,24] or novel bulk [9,10] and topological properties [5,6]. The potential enhancement of various functionalities beyond the accentuation of spectral gaps brought about by the use of sticky potentials is thus a very interesting avenue for investigation.

III. SPECTRAL PROPERTIES
Vibrational modes with frequencies {! n } are obtained by diagonalizing the dynamical matrix given by: Characteristic gap frequency (! ? i ) of the i-th linear elastic branch: Dimensionless gap parameter ( i ) between i and i + 1 linear elastic branch : where ! n = ! n+1 ! n and h. . . i ⌦i denotes the average over modes in the i-th branch of the vibrational spectrum. Linear elastic frequency: where µ is the shear modulus of the system.

IV. ELIMINATING PRE-STRESS CONTRIBUTIONS
The o↵-diagonal terms of the dynamical matrix is explictly given by: hi ji 2 6 6 6 6 6 6 6 4 X ↵= 8 > > > < > > > : 3  7 7 7 7 7 7 7 5 where k i j and l 0,i j are respectively the sti↵ness and rest lengths of the harmonic interaction between particles i and j that constit the e↵ective elastic network. By neglecting the pre-stress contributions (i.e setting f i j = l 0,i j = 0), we are able to write down a stress-free hessian tha completely characterized by the set of sti↵nesses n k i j o : Upon diagonalizingM, we obtain the vibrational spectrum of the stress-free network with frequencies ! 0 where k i j and l 0,i j are respectively the sti↵ness and rest lengths of the harmonic interaction between particles i and j that constitute the e↵ective elastic network. By neglecting the pre-stress contributions (i.e setting f i j = l 0,i j = 0), we are able to write down a stress-free hessian that is completely characterized by the set of sti↵nesses n k i j o : Upon diagonalizingM, we obtain the vibrational spectrum of the stress-free network with frequencies ! 0