REVIEW 3 major objections 3 minor 42 references
Spectroscopy of the Fractal Hofstadter Energy Spectrum
T0 review · 3 major / 3 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read Hofstadter's fractal spectrum resolved in twisted bilayer graphene
desk verdict First direct STS of Hofstadter subbands in TBG: raw data are strong and the model comparison is convincing, but the unverified gate-to-density calibration underpins every quantitative label. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the Hofstadter subband structure that emerges when a flat moiré band is placed in a perpendicular magnetic field with rational flux per unit cell Φ/Φ0 = 1/q. Each flat band folds into q subbands in a q-times-larger magnetic unit cell, and the gaps between subbands carry the topological labels (t, s), with the filling given by ν = t Φ/Φ0 + s; the integer t is the Chern number read off from the density-versus-field slope via the Streda formula. The argument for self-similarity rests on a discrete scaling transformation—multiplying energy and density by ¼(Φ0/Φ) at integer Φ0/Φ—that maps spectra at different rational fluxes onto each other near ν = ±4. A companion continuum-model LDOS calculation, including heterostrain and a filling-dependent Hartree potential, is used to match the measured dI/dV maps and to explain why certain Hofstadter gaps vanish.
What would settle it
A Hall-effect or capacitance measurement on the same device that shows the Vg-to-ν conversion is nonlinear or offset would shift every (t, s) label and break the self-similarity collapse, even if the raw spectral peaks are genuine. Alternatively, counting Hofstadter subbands at a flux where 1/q is not one of the reported values, say q = 7, and finding the number of subbands differs from q would falsify the fractionalization claim.
Extended reading notes
Core claim
The paper's central claim is that the flat bands of twisted bilayer graphene near the second magic angle (θ ≈ 0.6°) are an ideal arena in which Hofstadter's butterfly can be resolved directly by tunneling spectroscopy. At rational flux Φ/Φ0 = 1/q, the magnetic unit cell is enlarged q-fold, and each flat moiré band fractionalizes into q Hofstadter subbands; the resulting gaps appear at fillings ν = t Φ/Φ0 + s, where (t, s) are integer topological invariants. The paper further claims that STS maps at Φ/Φ0 = 1/6, 1/5, and 1/4 collapse onto one another when energy and density are rescaled by the linear factor ¼(Φ0/Φ), demonstrating discrete self-similarity near ν = ±4, and that non-integer flux ratios break this scaling. The measurements also show density-dependent changes in the spectrum, such as Hofstadter gaps that appear and disappear with filling, which the authors attribute to electron-electron interactions and model using a Hartree-corrected continuum calculation.
Load-bearing premise
The entire gap-labeling scheme assumes the back-gate voltage converts linearly and knownly to electron density, with the zero-field integer fillings ν = 4N resetting the unbound charge density; the paper gives no independent capacitance or density calibration.
Editorial extensions
If this is right
- At any rational flux 1/q, each flat moiré band should split into q resolvable subbands, giving a spectroscopic fingerprint for identifying flux ratios in moiré materials.
- The gap labels (t, s) assign Chern numbers to each spectroscopic gap, so STS maps can serve as a local, density-resolved probe of topology.
- The discrete self-similarity near ν = ±4 means that a single measured spectrum at one rational flux can be rescaled to predict the spectrum at another rational flux, within the single-particle picture.
- Correlation effects, captured by a density-dependent Hartree potential, change the Hofstadter spectrum qualitatively, so the fractal structure cannot be fully described by non-interacting models.
Reading between the lines
- If the density-axis calibration (Vg to ν) is accurate, the same (t, s) assignment should predict the quantized Hall conductance of each gap; a transport measurement on the same device would directly test this.
- The observed discrete self-similarity suggests that STS energy-density maps could be used as a quantitative local probe of flux density, potentially imaging spatial variations of magnetic field or strain.
- Because the Hartree-corrected model reproduces the vanishing gaps, a similar analysis could predict which correlated insulating states survive at finite field in other moiré flat-band systems.
- A direct test of the fractionalization claim would be to count Hofstadter subbands at a rational flux not reported here, such as q = 7 at lower field, and verify that the number of subbands equals q.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript reports scanning tunneling spectroscopy of twisted bilayer graphene near the second magic angle in a magnetic field, claiming the first direct spectroscopic observation of the fractal Hofstadter spectrum in a moiré material. At rational flux Φ/Φ0 = 1/q, the authors observe the fractionalization of flat moiré bands into q Hofstadter subbands, label spectroscopic gaps by topological invariants (t, s) via Streda's formula, and show that energy-versus-density maps near ν = ±4 at different rational fluxes can be collapsed onto each other by a discrete linear scaling transformation. The paper also reports interaction-driven modifications of the Hofstadter spectrum, including vanishing gaps and density-dependent subband splitting. The central experimental evidence is a set of dI/dV(Vs, Vg) maps, supported by LDOS calculations from a Bistritzer–MacDonald-based model with heterostrain and a filling-dependent Hartree potential.
Significance. If the claims hold, this would be a landmark result: the first direct spectroscopic access to Hofstadter subbands and their fractal self-similarity in a moiré system, with the gap labels providing Chern numbers through Streda's formula. The raw data show clear field-dependent subband splitting and nontrivial gap dynamics, and the paper should be credited for making the code available through a Code Ocean capsule, for comparing data to an explicit model, and for reporting an anomalous vanishing gap that the model qualitatively reproduces. However, the absolute density calibration underpinning every quantitative label is not established, and the model comparison depends on an in-preparation companion paper and unspecified interaction parameters. The significance is therefore conditional: the experimental phenomenology appears rich, but the headline claims about topological labels and discrete self-similarity are not yet supported to the standard required.
major comments (3)
- [Spectroscopy of Hofstadter Gaps; Fig. 3] The mapping from back-gate voltage Vg to filling ν is asserted but never calibrated. The spectral labels and all quantitative claims rely on the formula ν = t Φ/Φ0 + s, on Streda's formula for Chern numbers, and on the self-similarity transformation ν → ¼(Φ0/Φ)(ν − νref) + νref in Fig. 4. The Methods section states only that a p-doped Si back-gate was used; no dielectric thickness, capacitance, Hall density measurement, or compressibility check is given, and the zero-field suppressions at ν = 4N are used as density anchors. An offset or a few-percent nonlinearity in n(Vg) would shift every inferred filling; for example, δν ≈ 0.2 would turn the reported (4,4) gap label into a non-integer pair and invalidate the topological assignment. The paper must provide a capacitance calibration or an independent density measurement, and an uncertainty analysis showing that plausible gate-oxide or offset variations do not change any assigned (t, s).
- [Model comparison, Figs. 3b,d] The LDOS model used to reproduce the anomalous vanishing of the (4,4) gap is not self-contained. The main text specifies only that the model includes heterostrain and a filling-dependent Hartree potential and cites Ref. 33, an in-preparation paper, for the formalism; the functional form of the Hartree term, the screening parameter, and the precise parameter values for Devices A1 and A2 are not given. The statement that calculations without interactions show 'weaker connections' is not quantified. Because this model comparison is the primary evidence that the vanishing gap is a real interaction-enhanced effect, the paper should include the model equations, the parameter values used, and a sensitivity study over the Hartree interaction strength and strain parameters; otherwise the agreement cannot be independently assessed.
- [Self-Similar Fractal Hofstadter Spectrum; Fig. 4] The discrete self-similarity claim is supported only by visual inspection of the scaled data. No quantitative similarity metric, correlation coefficient, or residual analysis is provided for the collapse in the bottom panels of Figs. 4a,c, and the factor ¼ in the density transformation is not derived in the main text. The violations at non-integer Φ0/Φ (Figs. 4b,d) are likewise identified by eye. Since the same density transformation inherits the uncalibrated ν axis from the previous comment, the collapse could in part reflect a trivial linear rescaling of coordinates. A quantitative measure of feature alignment, with and without the calibrated density axis, is needed to support the self-similarity claim.
minor comments (3)
- [General] The text uses lowercase 'v' in several places where the filling ν is meant (e.g., 'near v = -4' and 'v = -4 V' in the self-similarity section); these should be corrected to ν.
- [Fig. 5 caption and related text] The reported gate voltages for fixed fillings, such as Vg = 16.8 V for ν = +4 and Vg = 3 V for ν = 0, imply a specific zero-bias offset that is never justified; a consistency check between these implied fillings and the uncalibrated Vg axis would help the reader evaluate the density scale.
- [Data and Code Availability] The Code Ocean capsule is mentioned but no capsule identifier or DOI is given, which makes the code availability statement hard to verify; please provide a persistent identifier.
Circularity Check
No load-bearing circularity found; the main spectroscopic claims are compared against, not derived from, the data. Minor self-citation to an in-preparation same-group theory paper is a reproducibility concern rather than a circular step.
full rationale
The paper's central claims - flat moiré bands fractionalizing into q Hofstadter subbands, the (t, s) gap labels, Chern numbers via Streda's formula, and the discrete self-similarity scaling near ν = ±4 - are checked against the measured dI/dV data using standard external theory (Hofstadter's 1976 spectrum, Streda's formula, the Bistritzer-MacDonald model). The gap-labeling relation ν = t Φ/Φ0 + s is a textbook/Streda relation and is not defined in terms of the measured spectra; the data determine the filling at which a gap appears, and the relation assigns integer invariants. The self-similarity scaling in Fig. 4 uses the theoretically expected factor ¼(Φ0/Φ)(ν − νref) + νref; this is not a fitted collapse, and the raw data must still show matching spectral features for the comparison to be meaningful. The principal weakness is the unverified linear back-gate calibration Vg → ν, on which all filling labels and the scaled density axes depend. A calibration offset or nonlinearity would shift every inferred ν, invalidating the topological labels and the self-similarity collapse. However, that is an experimental-systematics/correctness risk, not a circular reduction: no equation in the paper is equivalent to its own input by construction. The only self-citation of note is Ref. 33, an in-preparation paper by the same group, used for the LDOS simulations. Those simulations include device-specific heterostrain and a filling-dependent Hartree potential whose parameters are not fully specified in the main text, making the model comparison hard to verify independently. But the main text does not claim that these parameters were fitted to the very dI/dV data being 'predicted,' so a circularity finding cannot be sustained under the quoted-evidence rule. Overall, the derivation chain is not circular; the in-preparation same-group citation is a minor reproducibility concern.
Assumptions & free parameters
free parameters (3)
- Per-device heterostrain epsilon =
Device A1: 0.04%; Device A2: 0.07%
- Hartree interaction strength (screening parameter) =
not stated in main text
- Gate-voltage to filling conversion factor =
not stated in main text
assumptions (4)
- domain assumption The Bistritzer-MacDonald continuum model accurately describes the flat bands at the second magic angle.
- domain assumption At rational flux Φ/Φ0=1/q, each moiré flat band folds into q Hofstadter subbands with gap fillings ν = t Φ/Φ0 + s.
- domain assumption dI/dV measured at an AA site is proportional to the local density of states of the intrinsic TBG spectrum.
- domain assumption The moiré lattice constant is uniform and the flux relation Φ/Φ0 = AmB/(h/e) holds with Am from the STM topograph.
Cite this review
Pith. "Pith review of Spectroscopy of the Fractal Hofstadter Energy Spectrum." pith.science (2026). https://pith.science/paper/QQAB4UDY
@misc{pith2026250104777,
author = {Pith},
title = {Pith review of: Spectroscopy of the Fractal Hofstadter Energy Spectrum},
year = {2026},
howpublished = {\url{https://pith.science/paper/QQAB4UDY}},
note = {Machine review of arXiv:2501.04777}
}
read the original abstract
Hofstadter's butterfly, the predicted energy spectrum for non-interacting electrons confined to a two-dimensional lattice in a magnetic field, is one of the most remarkable fractal structures in nature. At rational ratios of magnetic flux quanta per lattice unit cell, this spectrum shows self-similar distributions of energy levels that reflect its recursive construction. For most materials, Hofstadter's butterfly is predicted under experimental conditions that are unachievable using laboratory-scale magnetic fields. More recently, electrical transport studies have provided evidence for Hofstadter's butterfly in materials engineered to have artificially large lattice constants, such as those with moir\'e superlattices. Yet to-date, direct spectroscopy of the fractal energy spectrum predicted by Hofstadter nearly 50 years ago has remained out of reach. Here we use high-resolution scanning tunneling microscopy / spectroscopy (STM / STS) to probe the flat electronic bands in twisted bilayer graphene near the predicted second magic angle, an ideal setting for spectroscopic studies of Hofstadter's spectrum. Our study shows the fractionalization of flat moir\'e bands into discrete Hofstadter subbands and discerns experimental signatures of self-similarity of this spectrum. Moreover, our measurements uncover a spectrum that evolves dynamically with electron density, displaying phenomena beyond that of Hofstadter's original model due to the combined effects of strong correlations, Coulomb interactions, and the quantum degeneracy of electrons in twisted bilayer graphene.
Reference graph
Works this paper leans on
-
[1]
Hofstadter, D. R. Energy levels and wave functions of Bloch electrons in rational and irrational magnetic fields. Phys. Rev. B 14, 2239–2249 (1976)
1976
-
[2]
Wannier, G. H. A Result Not Dependent on Rationality for Bloch Electrons in a Magnetic Field. Phys. Status Solidi B 88, 757–765 (1978)
work page 1978
-
[3]
Quantised Hall effect in a two-dimensional periodic potential
Streda, P. Quantised Hall effect in a two-dimensional periodic potential. J. Phys. C Solid State Phys. 15, L1299 (1982)
work page 1982
-
[4]
Gerhardts, R. R., Weiss, D. & Wulf, U. Magnetoresistance oscillations in a grid potential: Indication of a Hofstadter-type energy spectrum. Phys. Rev. B 43, 5192–5195 (1991)
work page 1991
-
[5]
Albrecht, C. et al. Evidence of Hofstadter’s Fractal Energy Spectrum in the Quantized Hall Conductance. Phys. Rev. Lett. 86, 147–150 (2001)
work page 2001
-
[6]
Geisler, M. C. et al. Detection of a Landau Band-Coupling-Induced Rearrangement of the Hofstadter Butterfly. Phys. Rev. Lett. 92, 256801 (2004)
work page 2004
-
[7]
Yankowitz, M. et al. Emergence of superlattice Dirac points in graphene on hexagonal boron nitride. Nat. Phys. 8, 382–386 (2012)
work page 2012
-
[8]
Dean, C. R. et al. Hofstadter’s butterfly and the fractal quantum Hall effect in moiré superlattices. Nature 497, 598–602 (2013)
2013
Show all 42 references
-
[9]
Hunt, B. et al. Massive Dirac Fermions and Hofstadter Butterfly in a van der Waals Heterostructure. Science 340, 1427–1430 (2013)
2013
-
[10]
Ponomarenko, L. A. et al. Cloning of Dirac fermions in graphene superlattices. Nature 497, 594–597 (2013)
2013
-
[11]
& MacDonald, A
Bistritzer, R. & MacDonald, A. H. Moiré bands in twisted double-layer graphene. Proc. Natl. Acad. Sci. 108, 12233–12237 (2011)
2011
-
[12]
Lu, X. et al. Multiple flat bands and topological Hofstadter butterfly in twisted bilayer graphene close to the second magic angle. Proc. Natl. Acad. Sci. 118, e2100006118 (2021)
2021
-
[13]
& MacDonald, A
Bistritzer, R. & MacDonald, A. H. Moir\’e butterflies in twisted bilayer graphene. Phys. Rev. B 84, 035440 (2011)
2011
-
[14]
Yu, G. L. et al. Hierarchy of Hofstadter states and replica quantum Hall ferromagnetism in graphene superlattices. Nat. Phys. 10, 525–529 (2014)
2014
-
[15]
Krishna Kumar, R. et al. High-temperature quantum oscillations caused by recurring Bloch states in graphene superlattices. Science 357, 181–184 (2017)
2017
-
[16]
Barrier, J. et al. Long-range ballistic transport of Brown-Zak fermions in graphene superlattices. Nat. Commun. 11, 5756 (2020)
2020
-
[17]
& LeRoy, B
Yankowitz, M., Ma, Q., Jarillo-Herrero, P. & LeRoy, B. J. van der Waals heterostructures combining graphene and hexagonal boron nitride. Nat. Rev. Phys. 1, 112–125 (2019)
2019
-
[18]
Spanton, E. M. et al. Observation of fractional Chern insulators in a van der Waals heterostructure. Science 360, 62–66 (2018)
2018
- [19]
-
[20]
Yoo, H. et al. Atomic and electronic reconstruction at the van der Waals interface in twisted bilayer graphene. Nat. Mater. 18, 448–453 (2019)
2019
-
[21]
Kazmierczak, N. P. et al. Strain fields in twisted bilayer graphene. Nat. Mater. 20, 956–963 (2021)
2021
-
[22]
Wong, D. et al. Cascade of electronic transitions in magic-angle twisted bilayer graphene. Nature 582, 198–202 (2020)
2020
-
[23]
Cao, Y. et al. Correlated insulator behaviour at half-filling in magic-angle graphene superlattices. Nature 556, 80–84 (2018)
2018
-
[24]
Cao, Y. et al. Unconventional superconductivity in magic-angle graphene superlattices. Nature 556, 43–50 (2018)
2018
-
[25]
Yankowitz, M. et al. Tuning superconductivity in twisted bilayer graphene. Science 363, 1059–1064 (2019)
2019
-
[26]
Lu, X. et al. Superconductors, orbital magnets and correlated states in magic-angle bilayer graphene. Nature 574, 653–657 (2019)
2019
-
[27]
B., Glicofridis, P
Chan, H. B., Glicofridis, P. I., Ashoori, R. C. & Melloch, M. R. Universal Linear Density of States for Tunneling into the Two-Dimensional Electron Gas in a Magnetic Field. Phys. Rev. Lett. 79, 2867–2870 (1997)
1997
-
[28]
P., Pfeiffer, L
Eisenstein, J. P., Pfeiffer, L. N. & West, K. W. Coulomb barrier to tunneling between parallel two-dimensional electron systems. Phys. Rev. Lett. 69, 3804–3807 (1992)
1992
-
[29]
Liu, X. et al. Visualizing broken symmetry and topological defects in a quantum Hall ferromagnet. Science 375, 321–326 (2022)
2022
-
[30]
Hu, Y. et al. High-Resolution Tunneling Spectroscopy of Fractional Quantum Hall States. Preprint at https://doi.org/10.48550/arXiv.2308.05789 (2023)
2023 doi
-
[31]
MacDonald, A. H. Landau-level subband structure of electrons on a square lattice. Phys. Rev. B 28, 6713–6717 (1983)
1983
-
[32]
Nuckolls, K. P. et al. Strongly correlated Chern insulators in magic-angle twisted bilayer graphene. Nature 588, 610–615 (2020)
2020
-
[33]
Scheer, Jonah Herzog-Arbeitman, Kevin P
Michael G. Scheer, Jonah Herzog-Arbeitman, Kevin P. Nuckolls, Ali Yazdani, & Biao Lian. Hofstadter Band Theory for Continuum Models. in preparation (2024)
2024
-
[34]
Bi, Z., Yuan, N. F. Q. & Fu, L. Designing flat bands by strain. Phys. Rev. B 100, 035448 (2019)
2019
-
[35]
& Bernevig, B
Lian, B., Xie, F. & Bernevig, B. A. Landau level of fragile topology. Phys. Rev. B 102, 041402 (2020)
2020
-
[36]
How Long Is the Coast of Britain? Statistical Self-Similarity and Fractional Dimension
Mandelbrot, B. How Long Is the Coast of Britain? Statistical Self-Similarity and Fractional Dimension. Science 156, 636–638 (1967)
1967
-
[37]
Fractal Geometry: Mathematical Foundations and Applications
Falconer, K. Fractal Geometry: Mathematical Foundations and Applications. (John Wiley & Sons, 2014)
2014
-
[38]
E., Ashoori, R
Dial, O. E., Ashoori, R. C., Pfeiffer, L. N. & West, K. W. High-resolution spectroscopy of two- dimensional electron systems. Nature 448, 176–179 (2007)
2007
-
[39]
& MacDonald, A
Nomura, K. & MacDonald, A. H. Quantum Hall Ferromagnetism in Graphene. Phys. Rev. Lett. 96, 256602 (2006)
2006
-
[40]
Young, A. F. et al. Spin and valley quantum Hall ferromagnetism in graphene. Nat. Phys. 8, 550–556 (2012)
2012
-
[41]
& Bernevig, B
Herzog-Arbeitman, J., Song, Z.-D., Regnault, N. & Bernevig, B. A. Hofstadter Topology: Noncrystalline Topological Materials at High Flux. Phys. Rev. Lett. 125, 236804 (2020)
2020
-
[42]
Wong, D. et al. A modular ultra-high vacuum millikelvin scanning tunneling microscope. Rev. Sci. Instrum. 91, 023703 (2020). Acknowledgements We thank B. Andrei Bernevig for useful discussions. This experimental work reported here was primarily supported by the DOE-BES grant D...
2020
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.