1. First study of the infinite temperature dynamical phase diagram in an interacting chain with an aperiodic modulation of the potential
2. Very thorough numerical study of spectral and dynamical response

1. Hard to extract fine structure of the ETH and MBL phases in the model because of larger finite-size effects in the Fibonacci case as opposed to the random or Aubry-Andre case
2. No analytic estimate of exponent b (defined below Eq. 4) or critical value h_c (see report)

The authors present a very through numerical study of the infinite temperature dynamical phase diagram of the interacting Fibonacci chain. Even though the non-interacting model exhibits critically delocalized single-particle wave functions at all strengths of the potential, the interacting model localizes at larger strength of the potential. The authors present a qualitative argument for this localization and discuss various special features of the MBL phase that are a consequence of the correlated binary disorder.

The study is sound, interesting in the context of localization, and will likely motivate future work. For these reasons, the referee recommends publication. The referee has a few suggestions/questions though; addressing these would significantly strengthen the manuscript.

Big picture questions/suggestions:

1) The authors qualitatively argue that adding interactions smoothens the Fibonacci potential in a mean-field picture. Could this picture be made quantitative? For example, at large h and small Delta, one could start with a typical Slate determinant eigenstate at \Delta=0, compute the average occupation on each site and thus the mean-field, solve for the single-particle wavefunctions with this mean field, feed the new wave functions into the Slater determinant and iterate. The referee expects that this procedure should converge deep in the MBL phase because the eigenstates are approximately product states. If it works, this procedure would allow the authors to extract a modified b (that might be h-dependent) and an estimate for the critical h at small \Delta. It would also be interesting to see how the multi-fractal structure of the starting wave functions is forgotten and how a localization length emerges in this iteration process.

2) The claims about a “multi fractal ETH” phase seem dubious and inconsistent with the qualitative argument for MBL. If the effective exponent b is greater than 1, then the single-particle wave functions are either extended or localized (except at a critical value). So there is really no reason to expect the multi-fractal structure of the b=1 case to show up in the single-particle density matrix of the interacting ETH phase at large L. The structure that the authors are seeing must be a small-size effect.

3) The discussion about the t^z growth of entanglement in the ETH phase is very confusing. As best the referee understands, this type of growth is observed in the non-interacting Fibonacci chain, in random models in the Griffiths region, in the interacting Aubry-Andre model at finite size, as well as in the interacting Fibonacci chain. Do the authors understand the origin of the z>1 growth in the interacting Fibonacci chain? If it is related to the short-distance multi-fractal structure that they observe in the single-particle density matrix, can the exponent in Fig. 9 be related to z in Fig. 12?

Minor questions/suggestions:

1) As another referee suggested, it would be helpful to see the behavior of the saturation value of the entanglement entropy vs L in the ETH phase.

2) As the authors note in the conclusions, other aperiodic sequences have larger geometrical fluctuations (or positive wandering exponents). These sequences might be better for finite-size study of the ETH-MBL transition as the number of independent samples at size L would be larger than in the Fibonacci case.

3) In the top paragraph of Page 6, the authors say that the conductivity at the Fermi level decays as L^{-\alpha}. The authors are presumably talking about T=0, but the relevant properties for this study are at T=infinity. How does the conductivity of the free model behave as a function of L in the infinite T ensemble?

See above

1- MBL investigated in a new form quasi-periodic potential where the single particle eigenstates exhibit multifractality.
2- An exhaustive study of observables for large systems amenable to exact diagonalization

1- The formulation of the problem is relatively less original.
2- Lack of a physical picture for the effects of multifractality on MBL and ETH.

The authors present a detailed study of MBL in a quasiperiodic Fibonacci chain. A quasiperiodic potential is free from Griffiths regions realised in systems with quenched disorder, where local fluctuations produce parts of the chain where the disorder is too weak to be localized. In the absence of such regions, the signatures of MBL are more robust and there are indications that the nature of the MBL-ergodic transition itself changes. The eigenstates of the non-interacting Fibonacci chain exhibit multifractality for all strengths of the potential. These are states which are intermediate between being exponentially localized and extended. Therefore, this system can possibly not be many-body localized or host a form of intermediate behaviour. This question is an interesting one.

In this article, the authors have investigated the spectrum and the eigenstates of the model using exact diagonalization and conclude the existence of a phase transition into an MBL phase. They have studied the transition using a number of static probes such as level statistics, local fermionic density, entanglement entropy and single particle density matrix. They have also studied the dynamics of imbalance and entanglement entropy. The evidence for an MBL phase in this model is very convincing. Although for some quantities like the level statistics (Fig. 4), the finite size scaling near the critical is not very 'clean', largely because the number of different realizations of the quasiperiodic potential is small for the largest system sizes that can be probed using ED. The evidence for the MBL phase is the main result of this work. They also show that the correlated nature of the potential produces certain magic angle states in the MBL phase which leave their imprints in the distribution of the local fermion density and half entanglement entropy. Although this work does not give a physical picture of the nature of MBL in multifractal systems and how it is distinct from conventional randomness, it is nonetheless a useful result and fit for publication.

In section 5.3, the authors present the scaling of the orbital participation entropy with system size for non-interacting and interacting models at low disorder. In the non-interacting case the coefficient of the ln L scaling gives the fractal dimension, which is less than one. On the other hand the authors show the scaling with L in the interacting case and claim this to be a 'multifractal ETH' phase, without any justification. It would be useful if there is a theoretical picture for using this term. Also, it would be desirable for the authors to comment on the relationship of this phase to the non-ergodic, delocalized phase seen in Bethe lattices and random regular graphs (for e.g. Phys. Rev. Lett. 113, 046806 (2014)). In this section, the authors should also include plots of this quantity in the localized phase and comment on any imprints of multifractality that does or does not exist.

In section 6.2 the authors study the growth of entanglement with time for product states. The free Fibonacci chain shows a power law growth, where the power depends on the strength of the quasiperiodicity. At strong quasiperiodicity the 'log-periodic' oscillations are associated with the discrete scale invariance. Based on the discussion in section 3 the discrete scale invariance presumably exists at all strengths of quasiperiodicity, yet the oscillations are absent low 'disorder'. Is this a finite size effect? In non-interacting Anderson localization the saturation value of the entanglement entropy is O(1) and independent of system size. What is the scaling with L of the saturation value of entanglement entropy in the free Fibonacci chain? How does it depend on the fractal dimension?

In the ETH phase, the entanglement entropy grows as a power law in time but the exponent does go to zero at the MBL transition, possibly due to finite-size effects. In the ETH phase does the entanglement entropy saturate to the infinite temperature EE? A plot showing the scaling with L of the saturated value of ETH and MBL phases would provide clarity on this question. In the MBL phase the growth is consistent with being logarithmic in time. The discussion of the fitting procedure is a bit unclear. The distinction between the coefficient of the logarithm between the ETH and MBL phases is not substantial. Could this be because the entanglement growth is much slower in the 'multifractal ETH' phase compared to quenched disorder? A comparison with the binary disorder model studied earlier would be very helpful in clarifying this point. In the plot of A vs h in Fig. 12 the L=18 data deviates quite significantly from the rest of the points. What is causing this?

Overall the manuscript is clearly written, well-structured and motivates the importance of the problem quite well. I recommend it for publishing once the authors address the points raised earlier.

Please refer to the report

1 - An interesting study of a model with some unique properties that distinguish it from others typically studied in the context of MBL.
2 - Comprehensive numerical study that covers many different quantities, with convincing arguments backing up the data.
3 - Likely to stimulate further interest in Fibonacci potentials in the MBL community.

1 - As a result of the large number of different quantities computed by the authors, many are introduced briefly and dispensed with very quickly, making it occasionally hard to follow.
2 - Though the results are solid and interesting, the manuscript would benefit from a stronger punchline as to the significance of the results.

This work takes a novel approach to the study of many-body localization by examining a model which is delocalized (and critical) in the non-interacting limit, and becomes localized only when interactions are switched on.

The authors consider a model of interacting fermions in one dimension subject to a quasiperiodic on-site potential of the Fibonacci form. This potential is balanced on a delicate knife-edge, being neither disordered enough to localize the non-interacting system nor smooth enough to entirely delocalize it: this results in the various interesting features of the non-interacting model that are well summarised by the authors.

Using a simple yet persuasive argument, the authors demonstrate that the addition of interactions acts to ‘smooth’ the Fibonacci potential, disrupting this delicate balance. With numerical simulations, the authors go on to demonstrate that the interacting system appears to follow typical MBL phenomenology, i.e. the system is delocalized at small disorder but there is a many-body localization transition at some critical disorder strength. The authors conduct a comprehensive numerical study and present various different quantities which corroborate this claim. As a consequence of the impressive variety of different quantities that the authors consider, some are presented in a rather brief way that relies on readers being already familiar with MBL phenomenology: some small improvements to the discussion would make the manuscript more accessible to a broader audience.

In summary, the authors have given convincing evidence for MBL phenomenology in this model and reviewed the novel features of the Fibonacci potential in a way that serves as an adept introduction of this model to an audience already familiar with MBL. This work likely to stimulate further interest in models with Fibonacci potentials in the MBL community.

I have a variety of minor comments, questions and suggestions which are included in this report. Based on a satisfactory response to these points, I’d be happy to recommend publication in SciPost Physics.

---

Remarks:

- The perturbative argument given for the ‘smoothing’ of the Fibonacci potential by interactions is plausible and intuitive. This suggests that, in the language of renormalization group, the non-interacting system is controlled by a fixed point unstable to the addition of interactions. The authors cite a variety of RG works which have examined the low-energy properties of the model, but as a future work, it would be interesting to see if RG methods could examine this model in the context of many-body localization.

- In Fig. 10, the bar over the I(t) suggests some form of averaging was used, but this isn’t mentioned anywhere in the text - was this data averaged in some way?

- In Fig. 12, why is the maximum time different for the two largest disorder strengths?

- In future directions, the authors remark that the fate of the transition in the non-interacting limit presents a technical challenge – what precisely is the challenge that this poses?

1 - The current statement in the abstract that ‘the MBL phase presents specific features’ is a little vague. Could the authors say something more definite here, perhaps giving an example of one?

2 - Can the authors clarify what they mean by the sentence “Finally, discrete disorder does not seem to affect the transition as compared to continuous random potential distributions” (p2)?

3 - Can the authors provide some further argument to clarify why there are L+1 Fibonacci words of length L? Currently this is stated without proof (p4).

4 - The authors comment on ‘large statistical errors’ (p4) as a result of the small number of samples available, due to the intrinsic nature of the model. Can the authors provide any further comment on which features of their study may be most strongly affected by these statistical errors?

5 - The existence of secondary peaks in the density (Fig. 5) is interesting and the explanation is ultimately convincing, though this section is a little convoluted to read. The left panel of Fig. 6, however, clearly demonstrates that the secondary peaks are associated with configurations where neighbouring sites have the same sign: I’d suggest leading with this observation before introducing the magic angle states.

6 - Section 5.3 covers a lot of concepts not previously introduced in the paper, and does so very quickly. A slightly expanded discussion of how multifractality can be quantified using the participation entropy would make this section easier to follow, particularly for those unfamiliar with these quantities.

7 - The use of a logarithmic fit to the growth of entanglement entropy in the ETH phase is explained in a confusing way. Although it makes sense after a few reads, can the authors find a clearer way to explain this?

A few typos I spotted:

1) p2: ‘random potentials distributions’ -> ‘random potential distributions’
2) p2: ‘recasted’ -> ‘recast’
3) p3: ‘eigenstates properties’ -> ‘eigenstate properties’
4) p5: ‘free fermions Fibonacci chain’ -> ‘free fermion Fibonacci chain’
5) p10: ‘were’ -> ‘where’
6) p17: ‘quench disordered systems’ -> ‘systems with quenched disorder’ would read better.
7) p7, caption of Fig. 3: the authors say the colour coding is ‘similar’ in both panels, but unless there is some difference, do they mean ‘the same’?
8) p18: ‘disordered strength’ -> ‘disorder strength’