{'fraction_non_alphanumeric': 0.043071397063457056, 'fraction_numerical': 0.01370596760980528, 'mean_word_length': 4.609579356663132, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 6, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
Quasiperiodic granular chains and Hofstadter butterflies

We study quasiperiodicity-induced localization of waves in strongly precompressed granular chains. We propose three different setups, inspired by the Aubry--Andr\'e (AA) model, of quasiperiodic chains; and we use these models to compare the effects of on-site and off-site quasiperiodicity in nonlinear lattices. When there is purely on-site quasiperiodicity, which we implement in two different ways, we show for a chain of spherical particles that there is a localization transition (as in the original AA model). However, we observe no localization transition in a chain of cylindrical particles in which we incorporate quasiperiodicity in the distribution of contact angles between adjacent cylinders by making the angle periodicity incommensurate with that of the chain. For each of our three models, we compute the Hofstadter spectrum and the associated Minkowski--Bouligand fractal dimension, and we demonstrate that the fractal dimension decreases as one approaches the localization transition (when it exists). Finally, in a suite of numerical computations, we demonstrate both localization and that there exist regimes of ballistic, superdiffusive, diffusive, and subdiffusive transport. Our models provide a flexible set of systems to study quasiperiodicity-induced analogs of Anderson phenomena in granular chains that one can tune controllably from weakly to strongly nonlinear regimes.


Introduction
Quasicrystals are solids whose structure is ordered but not periodic [1]. For many years, it was thought that quasicrystals could only be produced artificially. However, almost a decade ago, the first natural quasicrystal was found in Russia [2]. A common type of quasicrystal arises when atoms are arranged so that they possess symmetries, such as 5-fold symmetry, that are forbidden to periodic crystals 1 . A famous two-dimensional (2D) example is a Penrose tiling [4]. More generally, quasiperiodic structures can exist in any number of dimensions as structures with a broken translational symmetry. In 1D, the most common models used in the study of quasiperiodic systems are a Fibonacci quasicrystal [5] and the Aubry-André (AA) model [6]. These models are topologically equivalent to each other [7], in the sense that it is possible to continuously deform one into another without closing any bulk gaps.
A key feature of quasiperiodic potentials, arising prominently in the study of Schrödinger equations [6], is the transition from a "metallic" phase, in which all eigenstates are delocalized, to an "insulating" phase, in which they are localized. See, e.g., the analysis in [8,9] and the experiments in [11]. It is of considerable interest to extend these studies in various ways, including to nonlinear systems (see, e.g., the work of [10] and references therein), many-body systems, discrete systems, and settings with controllable interactions. For example, there have been relevant investigations in both bosonic and fermionic settings [12,13].
In the present paper, we study the effects of quasiperiodicity in strongly precompressed granular chains [14,15], in which neighboring particles interact with each other according to a Hertzian potential. Our aim is to illustrate that localization of eigenmodes can occur in quasiperiodic granular chains and to explore the conditions -with respect to both models and experimental setups -in which it occurs. We illuminate these conditions by comparing a variety of different models with one or both of on-site and off-site quasiperiodic structures. As demonstrated in a wealth of research over more than three decades [16], granular chains are extremely versatile, as one can adjust interaction potentials; readily tune them between almost linear, weakly nonlinear, and strongly nonlinear regimes by applying different precompression strengths; construct them using particles of different sizes, shapes, and material properties; and so on. This has yielded a wealth of insights about a diverse set of physical phenomena, including acoustic realizations of many concepts from condensed-matter physics [17]. Most centrally for our discussion, this includes the dynamics of wave transport and localization in disordered nonlinear systems in both theoretical [18,19] and experimental [20] studies. Other phenomena (both theoretical and applied) from condensed-matter and quantum physics that have been realized in granular crystals include an analog of a Ramsauer-Townsend resonance in the form of a square well [21], switching and rectification [22], and others. More broadly, granular chains provide a wonderful playground that enables systematic exploration of the role of lattice structure (e.g., material heterogeneity [23,24]) and fundamental dynamic (e.g., rogue waves [25]) and thermodynamic (e.g., equipartition [26]) phenomena.
The remainder of our paper is organized as follows. In Section 2, we briefly review the AA model. In Section 3, we present three models of 1D quasiperiodic lattices: two with on-site quasiperiodicity and one with off-site quasiperiodicity. In Section 4, we linearize the governing equations of our three models. In Section 5, we demonstrate that such models may possess a Hofstadter butterfly structure. In Section 6, we examine energy transport and localization in our models. We conclude and suggest some interesting directions for future research in Section 7.

A brief review of the Aubry-André model
The prototypical form of the Aubry-André (AA) model at the level of a tight-binding model is where Ψ is a wavefunction, n indexes the lattice site, E is energy, and V is a potential. We suppose that the on-site energy is modulated by a lattice distortion with wave vector q ∈ R, which is incommensurate with 2π and has phase r. We also suppose that V (x) = V (x + 2πn), where x ∈ R and n ∈ Z.
In [6], Aubry and André proved several fundamental properties of the eigenmodes of (2.1). They showed that a ground state exists and that it undergoes a transition from analyticity for λ < λc to nonanalyticity for λ > λc for certain λc(q) and when q is not a Liouville number. That is, q ∈ R \ Q and there exist γ and r such that is satisfied for any rational number p 1 /p 2 . Aubry and André also showed, using a perturbative approach, that the analyticity-breaking causes the eigenmodes of (2.1) to have very rich spatial properties. Specifically, one can write the eigenmodes of the modulated system as where vm are the coefficients of the Fourier expansion of V . The eigenmodes and eigenvalues for λ = 0 are given by e ikn and 2 cos(k), respectively. For the series in (2.3) to be convergent, one needs to satisfy a Diophantine condition. That is, there exist two positive constants K and β such that for any integers m and n.
For the sake of simplicity, consider the special case in which the phase r = 0 and V (ξn) = cos(2πξx), where ξ equals the golden ratio (1 + √ 5)/2. So, Eq. (2.1) takes the form: EΨn = Ψ n+1 + Ψ n−1 + λ cos(2πξn)Ψn. In this case, Aubry and André showed that for λ > λc = 2, all of the eigenmodes of (2.1) are exponentially localized, as in the Anderson model (in which the potential V is disordered rather than quasiperiodic) [27], with the same characteristic localization length (2.5) That is, Ψn decays asymptotically as e −n/ζ as n → ∞. However, when λ < 2, most eigenmodes are given by extended, modulated plane waves. Interestingly, this implies that the loss of analyticity is also associated with a transition at λc to spatial localization of the eigenmodes, where λc = 2 in this specific example. This transition is called a localization transition or Aubry-André transition. This result is a generic phenomenon in Schrödinger lattices, and it is thus relevant for a diverse variety of physical systems [11,[28][29][30], including photonic lattices, Bose-Einstein condensates, and many others. Additionally, the spectra of the corresponding Schrödinger operators (with two or even more frequencies) is a topic of intense mathematical interest; see, e.g., [31] and references therein.

Implementing the AA model in granular chains
Several recent studies have generalized conventional granular chains in various ways. They have yielded several interesting insights, and they promise to result in a host of others in the coming years [16,17]. One type of a generalized granular chain is a cradle system [32], in which particles are attached to linear oscillators, enabling the use of on-site potentials in a way that is independent of particle-particle interactions. Several potential realizations of such a setting were given in [33] (although they have yet to be implemented experimentally, to the best of our knowledge). Another fascinating variant arises from examining a chain of particles with internal resonators, such as by coupling a secondary particle inside a principal one. This leads to a locally-resonant granular chain, which is sometimes called a mass-in-mass (or mass-withmass, if the secondary particle is external) chain [34][35][36]. Additionally, the use of particles with non-spherical geometries can drastically modify particle-particle interactions. For instance, with cylindrical particles, although one has the same functional relation between the force and the displacements as with spherical particles, one can tune the magnitude of such interactions by changing the contact angle between adjacent cylinders [37,38]. We will consider granular chains with all of the above types of variations. The equations of motion in our general setting arë where un is the displacement of the nth particle (where n ∈ {1, 2, . . ., N }) measured from its equilibrium position in the initially compressed chain, vn is the displacement of the nth interior mass (when one particle is located inside another), and is a static displacement for each particle that arises from the static load F 0 = const. There is a Hertzian interaction between a pair of particles only when they are in contact, so each particle is affected directly only by its nearest neighbors and experiences a force from a neighbor only when it overlaps with it. This yields In our subsequent discussions (see Section 3), we will consider various special cases of (3.1) and (3.2), depending on the type of particle that we use to construct chains. Specifically, we work with three different models: two of them have spherical particles, and one has cylindrical particles. For our analysis and computations, we assume that αn, βn, and γn vary sinusoidally in space according to the following formulas: where ξ is the golden mean √ 5+1 2 (unless we explicitly state otherwise),ᾱ 1 >ᾱ 2 ≥ 0, andβ,γ ≥ 0. For simplicity, we separately examine the effects of (3.5)-(3.7). This leads to three different models in which it may be possible to observe an AA transition. In Eqs. (3.5)-(3.7),ᾱ 2 ,β, andγ are the quasiperiodic parameters that determine the strengths of the modulations for the different terms in Eqs.   so the coupling parameter αn is given by αn = An/mn, where , (3.8) and the nth particle has elastic modulus En, Poisson ratio νn, radius Rn, and mass mn. We assume that the particles are identical, so En = E, νn = ν, Rn = R, and mn = m. This, in turn, implies that αn =ᾱ 1 , and we letᾱ 1 = 1 without loss of generality.
(i) Model Ia:β = 0 andᾱ 2 =γ = 0 Suppose that the particles in the chain are attached to a mechanical restitution mechanism, such that there is a linear force in the equations of motion (3.1,3.2). In the limit of small angles, this can describe the well-known Newton's cradle, a system in which particles are aligned in one dimension and are suspended from a ceiling by strings so that the particles collide with each other along one dimension and oscillate. In the top panel of Fig. 1, we show a schematic of this system. Studies of this system have focused primarily on waves that arise by releasing one of the particles at one end with some velocity [33]. This produces a transfer of energy across the chain in the form of traveling waves [32]. Mulansky and Pikovsky studied disorder in closely related (nonlinearly coupled, and locally linear or nonlinear) oscillator systems [39], showing numerically and using a fractional-nonlinear-diffusion approach that energy transport is subdiffusive. This helps further motivate the study of modulated systems, such as the AA model, in granular chains. In our case, the equation of motions arë whereβ > 0. In other words, the Hookean spring constants are positive for all n. (ii) Model Ib:γ = 0 andᾱ 2 =β = 0 We now consider chains composed of particles that include an internal degree of freedom (i.e., mass-in-mass particles) [34]. Previous studies have focused on the generation [35,36], of such lattices, their traveling-wave solutions [40,41], and their (bright and dark) breather-like excitations [42,43]. These works illustrate that coherent structures and their dynamics are enriched significantly by the presence of the internal degree of freedom (DOF). To give just one example, incorporating an internal DOF in the particles can lead to nonlocal solitary waves with nonvanishing tails (so-called "nanoptera"), which have been observed experimentally in woodpile granular chains [44]. In our setting, we envision embedding a particle in the interior of each host particle, such that a particle and its interior mass are coupled via a linear restitution mechanism (such as a Hookean spring). The equations of motion, upon quasiperiodic modulation of the mass-in-mass (MiM) resonator, arë whereγ > 0.
(b) Off-site quasiperiodicity: The AA model with cylindrical particles Another approach for incorporating quasiperiodicity in a granular chain is by tuning particleparticle interactions. In existing experimental setups, to implement an AA model, one can use chains of cylindrical particles (rather than spherical ones). We are motivated by recent experiments [37], where it was demonstrated that cylindrical particles offer more flexibility than spherical particles for tuning particle-particle interactions, as one can change the contact angle between contiguous cylinders. Moreover, spatial and even temporal (periodic) variation of contacts between cylindrical particles has been proposed as an efficient way for implementing various functionalities, including that of an acoustic transistor in [38].
(i) Model II:ᾱ 2 = 0 andβ =γ = 0 To give equations of motion for a chain of cylindrical particles, we first need to know the form of the Hertzian coefficient in this case. For identical cylinders, the interaction coefficients are [46] αn(φn) = Y m g 1 (φn) g 2 (φn) g 3 (φn) where Y depends on the physical parameters of the particles, m is the mass of a particle, and φn is the contact angle between cylinders n − 1 and n and it is defined mod π/2. Explicit forms of Y and g i are where E is the elastic modulus, ν is the Poisson ratio, and R is the radius of the circular cross section of the cylinders. The functions K and E, respectively, are the complete elliptical integrals In panels (b) and (c), we show the contact-angle distributions for two cases: (i)ᾱ 1 = 3 andᾱ 2 = 1 and (ii)ᾱ 1 = 3 and α 2 = 3. We use arrows to represent the mapping process that we describe in the text.
of the first and second kinds [45]. They are where e = 1 − (b/a) 2 is the eccentricity of the contact area, a is the semi-major axis length of the ellipse, and b is its semi-minor axis length. One can approximate the quotient b/a by [(1 − cos φn)/(1 + cos φn)] 2/3 . This yields the following equations of motion: One can then control the interactions between particles by changing φn [38]. This raises the following question of what the distribution of contact angles {φn} has to be to obtain αn(φn) = α 1 +ᾱ 2 cos(2πξn). We address this issue by numerically inverting Eq. (3.12), so that quasiperiodic variation of αn yields a quasiperiodic variation of angles in the interval (φ min , φmax). In Fig. 2, we show the contact-angle distributions generated in two cases: (i) whenᾱ 1 = 3 andᾱ 2 = 1 and (ii) α 1 = 3 andᾱ 2 = 3. We also note that α(φn) → ∞ as φn → 0 and that αn(φn) has a lower bound at φn = π/2.

Linear approximation
Depending on the strength of the precompression that we apply to a granular chain and the magnitude of the strains that arise in (or are exerted on) the chain, one can expand the Hertzian force in a power series about the equilibrium state. This process reduces the equations of motion for the granular chain to ones that resemble those for a Fermi-Pasta-Ulam-Tsingou (FPUT) chain [17]. In particular, if the precompression is strong enough -specifically, if δn |u n−1 − un| for By considering (without loss of generality at the level of this linear approximation) a complex representation of the wavefunctions, un = φne iωt and vn = ψne iωt , and defining E = −ω 2 , we obtain Eφn = Bnφ n−1 + B n+1 φ n+1 − (Bn + B n+1 + βn + γn)φn + γnψn , (4.3) an eigenvalue problem that we can solve numerically by diagonalization. Using this linear description, we can now address the issue of a localization ("metal-insulator") transition for suitable incommensurate periodic coefficient variations of different types.

(a) Linear spectrum and localization transition
In addition to the linear spectrum, which we obtain by solving (4.3), we also compute the inverse participation ratio (IPR)  qualitative understanding of the nature of the linear modes, and a transition in the IPR also gives a way to quantitatively describe the AA transition.
In Fig. 3, we show the linear spectrum and IPR as a function of the quasiperiodicity parameter (which isβ in model Ia andγ in model Ib) for a chain of N = 200 particles for models Ia and Ib. We observe that these two models have a complex structure of bands and gaps, with some frequencies that appear isolated in the gaps and others that form bands that appear to cluster. Isolated frequencies are associated with modes that are similar to impurity modes. Similar structures of bands and gaps have been observed in other physical systems, such as in optics (see, e.g., [11,28]). Interestingly, we observe from the IPR that AA transitions occur in granular chains, most prominently in model Ia, where the transition is effectively identical to that in the original AA model. This is a consequence of modulating only the on-site potential with an external mechanism, so linearizations of the two systems yield the same equations. If Bn = 1 and γn = 0 for all n, the transition occurs atβ = 2. This differs starkly from the localization properties of the linear modes in the Anderson model, where low-frequency linear modes remain extended independently of the amount of disorder [18]. In Fig. 4(a), we show the transition to localization in the fundamental mode of model Ia. For model Ib, we double the number of modes in the system, because we double the number of DOF by incorporating the internal particles. This system has a very rich spectrum, where the upper part (half of the modes) has the same structure as in model Ia, but there is also a bottom part (the other half of the modes) associated with modes that do not undergo the localization transition and consequently are extended independently of the modulation. This is straightforward to explain by writing the system (4.1) in terms of in-phase (xn = un + vn) and out-of-phase (yn = un − vn) variables. This yieldsẋn = f (xn, x n±1 , yn, y n±1 ) andẏn = g(xn, x n±1 , yn, y n±1 ) − 2γnyn. Only the equations forẏn are affected explicitly by the quasiperiodicity. In this case, 2γn enters as a prefactor of yn, instead of γn as in model Ia. This explains why the localization transition occurs atγ = 1 instead of atγ = 2. Note that modes in the upper part of the spectrum also correspond to out-of-phase modes (between un and vn), whereas the bottom part of the spectrum is associated with in-phase modes. The latter do not see the quasiperiodicity in practice (because they effectively satisfy the original Hertzian dynamics without the MiM contribution), so they are generically extended.
In contrast to models Ia and Ib, model II does not have a localization transition; instead, we observe that all modes are extended, except for the ones that are associated with isolated frequencies in the gaps. In Fig. 4(b), we show the IPR as a function ofᾱ 2 for model II withᾱ 1 = 3. This suggests that, without an on-site potential, one cannot observe this sort of transition in granular chains of cylindrical particles. In the future, it will be particularly worthwhile to explore the generality of this conclusion. Specifically, a relevant question is whether it is generically the case that it is impossible for inter-site interactions, modulated by one or more frequencies, to induce a localization transition in a granular chain.

Hofstadter butterfly
Another property of the AA model's spectrum is its fractal nature. To explore this, we compute the spectrum as a function of ξ (see Fig. 5), and we observe a structure that is known as a Hofstadter butterfly. The butterfly is a footprint of the spectrum's fractality, and one can see its statistical self-similarity in the figure.
The Hofstadter butterfly was first predicted in 1976 [47] in a completely different system: Bloch electrons on 2D lattices and in the presence of an orthogonal magnetic fields. In typical crystals, one needs to use magnetic fields that are at least of the order of several thousands of tesla to observe a Hofstadter butterfly. As a result, it took until 1997 -in a microwave system [48] -that a Hofstadter butterfly was observed experimentally. A Hofstadter butterfly was then observed in graphene in 2013 [49] and using interacting photons in superconducting qubits in 2017 [50]. The possibility to also observe Hofstadter butterflies in granular chains is very exciting, given the simplicity and controllability of the latter system.   To characterize the self-similarity of the spectrum in the different cases, we compute the Minkowski-Bouligand fractal dimension (Dm), which (by Moran's theorem) is the same as the Hausdorff dimension (D h ) for strictly self-similar fractals [51,52]. The procedure to numerically compute D M is known as box-counting, and we proceed as follows. First, we map the Hofstadter spectrum into a square of 480 × 480 pixels, we then partition the square into boxes of characteristic size (side length) l B , and finally we count the number N B of boxes that include at least one point of the spectrum.  (b) Minkowski-Bouligand fractal dimension D M as a function of the quasiperiodic parameterβ,γ, andᾱ 2 for models Ia, Ib, and II, respectively. We use the parameter valueᾱ 1 = 3 for model II.
of the logarithm of the data using gradient descent. In Fig. 6(a), we show an example of the box counting for the Hofstadter spectrum of model Ia at the localization transition (i.e., whenβ = 2).
To show how the band gaps are filled with boxes, we have superimposed the N B boxes over the spectrum for different values of l B . We now compute the fractal dimension D M as a function of the quasiperiodic parametersᾱ 2 , β, andγ for our three models. We expect D M to be between 1 and 2, because D M = 1 for a line and D M = 2 for a plane. In Fig. 6(b), we show the results of our computations. We observe that the minimum fractal dimension occurs at the same point as the localization transition for models Ia (atβ = 2) and Ib (atγ = 1). We calculate that D M 1.69 for model Ia and D M 1.82 for model Ib. It is interesting to note the similar non-monotonic dependence of the fractal dimension on the model parameters in models Ia and Ib. Presumably, this arises from the aforementioned similarity of the former model and the out-of-phase excitations of the latter model. In contrast, in model II, forᾱ 1 = 3 andᾱ 2 ∈ (0, 3), we observe that the fractal dimension decreases monotonically asᾱ 2 increases.

Energy transport and localization
Another important issue, which one can examine in several ways, is how energy transport is affected by the quasiperiodicity [53,54]. For instance, one can compute a second moment m 2 of the energy distribution as a function of time to quantitatively characterize the temporal evolution of the energy distribution's width [18][19][20][55][56][57][58][59][60]. Following recent studies in disordered granular chains [20], we study the evolution of the energy distribution {Hn} N n=1 immediately after the impact of a striker against the first particle (n = 1). We write the particles' energy as The total energy is H = n Hn, which is a conserved quantity. Note that, in the context of experiments, one should also consider dissipation which is neglected here; see a relevant summary of models thereof in [17]. Importantly, we expect that our results will be robust enough to be observable experimentally even in the presence of a small amount of dissipation. This claim is supported by recent experimental results on the observation of other linear phenomena (e.g., Anderson-like localization [20] and an analog of a Ramsauer-Townsend resonance [21]) in granular chains. We then compute the second moment m 2 (t) = n (n − 1) 2 Hn n Hn , (6.2) and we estimate a scaling relationship between m 2 and t. When the scaling is approximated reasonably as a power law -for which, in exact form, m 2 (t) ∼ t η as t → ∞ for some exponent η -one can categorize as ballistic the case when η = 2, superdiffusive the one when η ∈ (1, 2), diffusive when η = 1, and subdiffusive when η ∈ (0, 1). If η = 0, we say that there is no diffusion; in other words, all energy remains localized. In a perfectly homogeneous granular chain, it is known that energy transport is ballistic. However, when disorder is added to the chain, the dynamics can change drastically, and one can observe different energy transport regimes [18][19][20].
These previous studies have focused on the interplay between disorder and nonlinearity. Here, in contrast, we show that even in a strongly precompressed (i.e., almost linear) chain, one can obtain any desired exponent diffusion η ∈ [0, 2] for the energy transport. However, we find that on-site potentials are essential to have localization. One advantage of working with quasiperiodic chains instead of disordered ones is that we do not need to compute averages over a large number of realizations to obtain robust insights. Our quasiperiodic chains are produced in a deterministic way, so given a parameter ξ, one gets one specific chain. This enables us to cover the whole parameter space with considerably fewer computations than when studying disordered chains. We are also interested in characterizing energy transport in realistic frameworks, so we set N = 21, which gives a long enough chain to qualitatively capture the nature of transport, at least in several recent experimental and theoretical explorations [18,20].
To integrate (3.1), we use a fifth-order explicit Runge-Kutta (RK5) method with a step size of dt = 0.01. We set δn = 1 for all n. We also set u 0 (t) = u N +1 (t) = 0 for all t, so we have fixed boundary conditions on both sides. We use the stopping criterion for our simulations that either T = 20 or that energy reaches the boundary opposite to the one impacted by the striker. We use the former condition to stop the code in cases in which all of the energy is trapped in the form of localized states. In other words, there is no diffusion. This is expected, for instance, in model Ia for β > 2 (i.e., after the localization transition occurs). However, it is an uncommon scenario in model II, which does not have a localization transition. In Fig. 7, we show our results for energy transport in our three models. In models Ia and Ib, we explore the parameters rangesβ ∈ [0, 4] andγ ∈ [0, 3], respectively, so we can compare energy transport on both sides of the localization transition. In model II, we considerᾱ 1 = 3 andᾱ 2 ∈ [0, 3]. For all three models, and for ξ ∈ [0, 1], we can tune the energy transport over from subdiffusive to ballistic behaviors. This allows a great deal of control of the energy-transport properties, and it is remarkable that we are able to do this using a deterministic model. In model Ia, we also observe localization (for which η = 0) in contrast with observations in disordered granular crystals [18][19][20]. This suggests that the inclusion of on-site potentials is crucial for this localization phenomenon.

Conclusions
In this article, we introduced different types of quasiperiodic granular chains that were inspired by the work by Aubry and André in condensed-matter physics and by recent developments (cradle and mass-in-mass systems) in the context of granular lattice systems. We studied the localization and spectral properties of such chains. To achieve each type of quasiperiodic chain, we incorporated spatial modulation (which is incommensurate with the chain's period) into one of the physical parameters. We proposed three models: models Ia and Ib use spherical particles, in which quasiperiodicity enters via an on-site potential -either in the form of a local oscillator  (as in a cradle) or in the form of a local resonator (for a mass in mass system) -and model II uses cylindrical particles. In models Ia and Ib, we demonstrated the existence of an analog of the wellknown AA transition. However, in model II, we showed that, without an on-site potential and with quasiperiodicity affecting only inter-site interactions, a localization transition cannot occur. We also computed the Hofstadter spectrum for each of the models and studied their fractal properties by computing the Minkowski-Bouligand fractal dimension. In models Ia and Ib, we showed that the minimum fractal dimension D M of the spectrum coincides with the point at which the localization transition occurs. For model II, we observed that the spectrum's fractal dimension decreases monotonically as a function of the quasiperiodic parameter.
Finally, we numerically studied energy transport by exciting the granular chains with a striker. We demonstrated the existence of different regimes -ranging from ballistic to subdiffusive transport -as well as localization. In contrast to prior work, which achieved such control using a combination of disorder and adjusting the nonlinearity strength [18][19][20], we were able to control the energy transport using a deterministic model in a strongly precompressed (and almost linear) granular chain.
Naturally, it will be particularly valuable to implement some of these ideas in laboratory experiments, as one can then further explore the role that granular systems can play in the study of quasiperiodic operators [48][49][50]. Achieving an experimental realization of a Hofstadter butterfly in a granular chain would also be very exciting in its own right. Among the settings that we have proposed in this paper, model Ib (i.e., the chain of mass-in-mass particles) is the clearest candidate for observing a localization transition (in the out-of-phase variables), given that the cradle system has yet to be experimentally realized. Arguably, the woodpile setup of [44] may also constitute an excellent playground for such studies. However, a key consideration, given experimental limitations, is that one seeks to build a chain with as few particles as possible such that one can (still) observe the relevant phenomenology.
There are also numerous open issues to explore computationally and theoretically. Examples include the effect of nonlinearity (e.g., through larger excitation amplitudes) on these modes and their localization, how these phenomena differ for granular crystals in different numbers of dimensions, and others. Relevant extensions will be considered in future studies.