REVIEW 2 major objections 4 minor 40 references
Re-entrant localization induced by short-range hopping in the fractal Rosenzweig-Porter Model
T0 review · 2 major / 4 minor · reviewed 2026-08-12 · deepseek-v4-flash
Pith's one-line read In the fractal Rosenzweig-Porter model, the fractal dimension D2 of bulk eigenstates is a non-monotonic function of nearest-neighbor hopping: increasing κ first suppresses D2 and later re-enters the ergodic phase.
desk verdict A solid extension of the fractal RP model with a genuinely new non-monotonic re-entrant effect, held back mainly by an unproved homogenization step that should be checked directly. 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 mechanism is the counting argument that compares the number of fractal levels $L^{{1-f(b)}}$ in a spectral interval of width $L^{{-b}}$ with the number M_b ∼ $L^{{1+b+2k}}$ of levels that the nearest-neighbor hopping shifts by that same energy. The crossover scale b* = |2k|/(1+d) separates intervals where the fractal disorder's counting function f(b)=d·b remains intact from intervals where hybridization homogenizes the levels and changes f(b) to b − (1−d)b*. Plugging this modified f(b) into the self-consistency equation 1 + 2a − f(a) = γ for the Lorentzian broadening Γ_d yields the phase diagram of Eqs. (4)–(6). In the strong-hopping regime the short-range term creates blocks of size ξ ∼ $κ^{2}$ and the argument becomes a block-RP comparison between Γ_1 and the block level spacing.
What would settle it
For a fixed disorder strength γ inside the predicted re-entrant window (e.g., d=0.6, γ≈1.5), compute the fractal dimension D2 from IPR scaling at system sizes L=$2^{7}$ to $2^{14}$. The theory predicts D2(κ) first decreases below its κ=0 value and then increases; observing D2 to be monotonically increasing in κ, or the localized phase to disappear, would falsify the central claim. A second, sharper falsifier is the level-spacing ratio: the theory predicts a return to Poisson statistics in the fractal IV region (k>0, γ>2−k), whereas a Wigner-Dyson ratio there would contradict the phase diagram.
Extended reading notes
Core claim
On its own terms, the paper establishes that adding a short-range kinetic term to the fractal Rosenzweig-Porter model produces a re-entrant localization phase diagram. For a fixed disorder strength γ, starting from κ=0, the fractal dimension D2 of bulk eigenstates first decreases with increasing κ, signifying a tendency toward localization, and then increases again beyond a critical κ, re-entering the ergodic (or a higher-fractal) phase. The analytical treatment splits the (κ,γ) plane into three regions: for strong hopping (κ ∼ L^k, k>0) the short-range term dominates and produces a block-like fractal phase with D2 = 2k at large γ; for intermediate weak hopping (−(1+d)/(2d) < k < 0) competition between the two delocalizing mechanisms produces four phases including a new 'fractal II' regime; and for very weak hopping (k < −(1+d)/(2d)) the fractal disorder alone controls the diagram. The mechanism is the interplay between the local-in-energy Fermi-golden-rule broadening Γ_d, which is sensitive to the fine level-spacing structure of the fractal disorder, and the local-in-space nearest-neighbor hopping, which is insensitive to that structure. When the two mechanisms compete, the effective level statistics are reshuffled and the system can transiently localize.
Load-bearing premise
The counting argument in Eqs. (14)–(15) assumes that once the number of levels shifted by the nearest-neighbor hybridization M_b in a spectral window exceeds the number of fractal levels $L^{{1-f(b)}}$, every level in that window is fully hybridized and its level statistics become locally homogeneous; if real fractal disorder does not mix completely at this crossover, the boundaries of the fractal II region and the re-entrant loop would shift.
Editorial extensions
If this is right
- For strong hopping (k > 0), the model has no Anderson localization transition; instead, beyond γP = 2 − k it enters an extended phase with Poisson level statistics and fractal dimension D2 = 2k, and for k ≥ 1/2 it is ergodic for all γ.
- For intermediate hopping (−(1+d)/(2d) < k < 0), the phase diagram contains four distinct phases — ergodic, fractal I, fractal II, and localized — separated by transitions at γET = 1, γFT = 1 + ((2−d)/(1+d))|2k|, and γAT = 2 + ((2−2d)/(1+d))|2k|.
- For very weak hopping (k < −(1+d)/(2d)), the phase diagram coincides with that of the fractal RP model alone, with γAT = 2/d.
- The re-entrant behavior becomes weaker as the fractal dimension d of the disorder approaches 1, where the on-site disorder effectively becomes uncorrelated.
Reading between the lines
- The same competition between an energy-local broadening and a space-local hopping should appear in any single-particle model with a multifractal (Cantor-like) on-site potential, making the re-entrant effect a general feature rather than a special property of the RP ensemble.
- The counting criterion δh_max ∼ L^{2k+1} > δtyp ∼ L^{-1/d} (i.e., 2k > −(1+1/d)) predicts where re-entrant behavior appears; this criterion could be used to design other random-matrix models that exhibit a non-monotonic phase diagram without solving for the full D2.
- If this mechanism transfers to many-body Fock space, a local perturbation that is insensitive to the fractal structure of the many-body spectrum could be used to tune a system from ergodic to many-body localized and back, offering a potential control knob for quantum information storage.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies the fractal Rosenzweig-Porter model with an added nearest-neighbor hopping term of amplitude κ∼L^k. Its central claim is that, for a fixed disorder strength γ, the bulk eigenstate fractal dimension D2 is a non-monotonic function of κ: increasing κ from zero first reduces D2, driving the system from a fractal toward a localized state, and then, beyond a critical κ, the trend reverses and D2 grows back toward ergodic behavior. The authors propose an analytical phase diagram, Eqs. (4)-(6), built from a Lorentzian eigenfunction ansatz and a self-consistency equation for the broadening Γ, together with a counting argument for how nearest-neighbor hopping modifies the level-counting function f(b) of the fractal diagonal disorder (Eqs. (14)-(15)). They support the phase diagram with exact diagonalization data for D2(γ) at several fixed k, reporting good agreement.
Significance. If correct, this is a valuable counterexample to the standard expectation that increasing kinetic coupling delocalizes: it shows that a spatially local hopping term can first enhance localization by destroying the fine fractal level structure before eventually dominating and restoring ergodicity. The manuscript is careful with limiting cases: the piecewise D2 formulas match at all phase boundaries and reduce to the known RP and fractal-RP results in the appropriate limits. The numerics in Figs. 2 and 3 agree with the proposed formulas over the system sizes and parameter values shown. The main weakness is that the derivation of the modified f(b) in Region 2 is heuristic: it relies on a strong mixing assumption that is not proven and, as detailed in the major comments, uses a scaling input for the pair-difference distribution that appears inconsistent with the model's own fractal disorder construction. These issues are load-bearing for the analytical phase diagram but appear fixable within the manuscript's scope.
major comments (2)
- [Weak nearest-neighbour hopping; Eqs. (14)-(15)] The step from the counting inequality L^{1-f(b)} ≪ M_b to the conclusion that 'all the levels in the interval L^{-b} will be hybridized ... and redistributed homogeneously' is an unproven strong-mixing assumption. M_b counts shifted levels over the entire sample, whereas L^{1-f(b)} counts levels in a single energy window; the inference that every level in every such window participates in a shift event requires an ergodicity/uniformity property of the near-degenerate pairs that is not established. Because the entire fractal II region and the re-entrant boundaries γ_FT and γ_AT in Eq. (5) are computed from the f(b) of Eq. (15), this assumption is load-bearing. I would ask the authors either to derive it from the disorder statistics or to verify f(b) directly by numerical level counting in the dressed spectrum; the D2 agreement in Fig. 3 alone does not test this intermediate step.
- [Weak nearest-neighbour hopping; after Eq. (12)] The statement that 'the marginal distribution P(h_m - h_{m±1}) ... is regular close to the origin' is inconsistent with the fractal disorder construction used in the paper. For L points whose sorted spacings are Pareto with a hard cutoff δ_typ∼L^{-1/d}, the pairwise difference distribution of the shuffled sequence has cumulative probability C(x)∼x^d for x≫δ_typ and vanishes for x<δ_typ; the density diverges as x^{d-1} rather than being regular. Consequently p_b in Eq. (14) should scale as L^{d(b+2k)}, not L^{b+2k}. Repeating the counting then gives b*=|k| instead of |2k|/(1+d), which shifts γ_FT and γ_AT in Eq. (5) and changes the extent of the fractal II and re-entrant regions. The authors should either correct the counting to use the actual correlation exponent of the disorder (and check whether the final phase diagram survives with re-fit parameters), or state and use a different disorder construction for which the regular-difference assumption holds.
minor comments (4)
- [Fig. 1] The caption of Fig. 1(b) states that the horizontal dotted lines denote the γ values chosen for panel (b), but the values of γ are not given in the text or figure; please include them.
- [Eq. (1)] The notation 'h_n = R_mn = 0' should presumably denote ensemble averages, e.g. with an overline; as written it is ambiguous.
- [Fig. 1(a) caption] The caption describes the result as an 'exact analytical computation'; given the heuristic counting argument used for Region 2, 'analytical' would be more accurate than 'exact'.
- [Before Eq. (11)] The sentence preceding Eq. (11) ends with 'The corresponding energy shift is' followed directly by the equation; the text should explicitly state that Eq. (11) is the second-order perturbative shift.
Circularity Check
No significant circularity: the re-entrant phase diagram follows from prior fractal-RP results plus counting arguments; the target D2(κ) is not an input to the derivation.
full rationale
The paper's central claim, the non-monotonic behavior of D2 as κ increases, is not contained in, or forced by, the cited prior results. Region 3, Eq. (6), is explicitly the known fractal-RP result of Ref. [29] (same authors), but it is used as the k→−∞ limit, not as the re-entrant prediction. Region 2, Eqs. (14)-(16), is derived by comparing the number of fractal levels in an energy window, L^{1−f(b)}, with the number of nearest-neighbor-induced shift events, M_b∼L^{1+b+2k}, and then solving the self-consistency equation (8) with the modified f(b). The b* and b** parameters are fixed by d and k through Eq. (14), with no parameter fitted to the IPR data; the numerical D2 values in Figs. 2-3 are obtained independently by exact diagonalization. The assumption that all levels in the window are fully hybridized once M_b exceeds L^{1−f(b)} is a physical ansatz rather than a circular step; if wrong it would shift the phase boundaries, but that is a correctness risk, not an equivalence to inputs. The reductions of Eqs. (4) and (6) to the standard RP and 1D Anderson limits are consistency checks. Self-citations [17,29,36] are used as published building blocks, not as unverified uniqueness claims or fitted inputs, so they do not make the derivation circular.
Assumptions & free parameters
assumptions (5)
- domain assumption Unperturbed RP eigenfunctions have Lorentzian form and Gamma_d satisfies the self-consistency equation 1 + 2a - f(a) = gamma (Eqs. 7-8).
- domain assumption Fractal disorder is realized by Pareto-distributed level spacings randomly reshuffled, so nearest-neighbor differences h_m - h_{m+1} have a regular distribution near zero.
- domain assumption For kappa >> 1 the eigenstates form hybridized blocks of size xi ~ kappa^2 = L^{2k}, using the 1D Anderson localization length.
- ad hoc to paper On energy scales finer than L^{-b*}, hybridization homogenizes the fractal level set, giving f(b) = b - (1-d)b* for b > b* (Eq. 15).
- standard math Only power-law L-dependencies matter; slower dependencies are neglected in saddle-point counting (footnote 33).
Cite this review
Pith. "Pith review of Re-entrant localization induced by short-range hopping in the fractal Rosenzweig-Porter Model." pith.science (2026). https://pith.science/paper/XKMVY3FY
@misc{pith2026241116851,
author = {Pith},
title = {Pith review of: Re-entrant localization induced by short-range hopping in the fractal Rosenzweig-Porter Model},
year = {2026},
howpublished = {\url{https://pith.science/paper/XKMVY3FY}},
note = {Machine review of arXiv:2411.16851}
}
abstract
Typically, metallic systems localized under strong disorder exhibit a transition to \imk{delocalization} %finite conduction as kinetic terms increase. In this work, we reveal the opposite effect~--~increasing kinetic terms leads to an unexpected \imk{reduction of mobility, }%suppression of conductivity, enhancing localization of the system, and even lead to re-entrant delocalization transitions. Specifically, we add a nearest-neighbor hopping with amplitude \(\kappa\) to the Rosenzweig-Porter (RP) model with fractal on-site disorder and surprisingly see that, as \(\kappa\) grows, the system initially tends to localization from the fractal phase, but then re-enters the ergodic phase. We build an analytical framework to explain this re-entrant behavior, supported by exact diagonalization results. The interplay between the spatially local $\kappa$ term, insensitive to fractal disorder, and the energy-local RP coupling, sensitive to fine-level spacing structure, drives the observed re-entrant behavior. This mechanism offers a novel pathway to re-entrant localization phenomena in many-body quantum systems.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[1]
J. M. Deutsch, Quantum statistical mechanics in a closed system, Phys. Rev. A 43, 2046 (1991)
1991
-
[2]
Srednicki, Chaos and quantum thermalization, Phys
M. Srednicki, Chaos and quantum thermalization, Phys. Rev. E 50, 888 (1994)
1994
- [3]
-
[4]
I. V. Gornyi, A. D. Mirlin, and D. G. Polyakov, Interact- ing electrons in disordered wires: Anderson localization and low-t transport, Phys. Rev. Lett. 95, 206603 (2005)
2005
-
[5]
D. A. Abanin, E. Altman, I. Bloch, and M. Serbyn, Col- loquium: Many-body localization, thermalization, and entanglement, Rev. Mod. Phys. 91, 021001 (2019)
2019
-
[6]
D. J. Luitz, N. Laflorencie, and F. Alet, Many-body local- ization edge in the random-field Heisenberg chain, Phys. Rev. B 91, 081103 (2015)
2015
-
[7]
Sels and A
D. Sels and A. Polkovnikov, Dynamical obstruction to localization in a disordered spin chain, Phys. Rev. E 104, 054105 (2021)
2021
-
[8]
ˇSuntajs, J
J. ˇSuntajs, J. Bonˇ ca, T. c. v. Prosen, and L. Vidmar, Quantum chaos challenges many-body localization, Phys. Rev. E 102, 062144 (2020)
2020
Show all 40 references
-
[9]
De Roeck, F
W. De Roeck, F. Huveneers, M. M¨ uller, and M. Schiulaz, Absence of many-body mobility edges, Phys. Rev. B 93, 014203 (2016)
2016
-
[10]
J. Z. Imbrie, Multi-Scale Jacobi Method for Anderson Localization, Communications in Mathematical Physics 341, 491 (2016)
2016
-
[11]
W. D. Roeck, L. Giacomin, F. Huveneers, and O. Prosniak, Absence of normal heat conduction in strongly disordered interacting quantum chains (2024), arXiv:2408.04338 [math-ph]
2024 arXiv
-
[12]
Sierant, M
P. Sierant, M. Lewenstein, A. Scardicchio, L. Vidmar, and J. Zakrzewski, Many-body localization in the age of classical computing (2024), arXiv:2403.07111 [cond- mat.dis-nn]
2024 arXiv
-
[13]
V. N. Smelyanskiy, K. Kechedzhi, S. Boixo, S. V. Isakov, H. Neven, and B. Altshuler, Nonergodic delocalized states for efficient population transfer within a narrow band of the energy landscape, Phys. Rev. X 10, 011017 (2020)
2020
-
[14]
Kechedzhi, V
K. Kechedzhi, V. N. Smelyanskiy, J. R. McClean, V. S. Denchev, M. Mohseni, S. V. Isakov, S. Boixo, B. L. Altshuler, and H. Neven, Efficient population transfer via non-ergodic extended states in quantum spin glass, arXiv:1807.04792 (2018)
2018 arXiv
-
[15]
P. W. Anderson, Absence of diffusion in certain random lattices, Phys. Rev. 109, 1492 (1958)
1958
-
[16]
Evers and A
F. Evers and A. D. Mirlin, Anderson transitions, Rev. Mod. Phys 80, 1355 (2008)
2008
-
[17]
V. E. Kravtsov, I. M. Khaymovich, E. Cuevas, and M. Amini, A random matrix model with localization and ergodic transitions, New J. Phys. 17, 122002 (2015)
2015
-
[18]
Facoetti, P
D. Facoetti, P. Vivo, and G. Biroli, From non-ergodic eigenvectors to local resolvent statistics and back: A random matrix perspective, Europhys. Lett. 115, 47003 (2016)
2016
-
[19]
Truong and A
K. Truong and A. Ossipov, Eigenvectors under a generic perturbation: Non-perturbative results from the random matrix approach, Europhys. Lett. 116, 37002 (2016)
2016
-
[20]
Monthus, Statistical properties of the Green function in finite size for Anderson localization models with multi- fractal eigenvectors, J
C. Monthus, Statistical properties of the Green function in finite size for Anderson localization models with multi- fractal eigenvectors, J. Phys. A: Math. Theor.50, 295101 (2017)
2017
-
[21]
Bogomolny and M
E. Bogomolny and M. Sieber, Eigenfunction distribu- tion for the Rosenzweig-Porter model, Phys. Rev. E 98, 032139 (2018)
2018
-
[22]
von Soosten and S
P. von Soosten and S. Warzel, Non-ergodic delocalization in the Rosenzweig–Porter model, Letters in Mathemati- cal Physics , 1 (2018)
2018
-
[23]
Venturelli, L
D. Venturelli, L. F. Cugliandolo, G. Schehr, and M. Tarzia, Replica approach to the generalized Rosenzweig-Porter model, SciPost Phys. 14, 110 (2023)
2023
-
[24]
Hopjan and L
M. Hopjan and L. Vidmar, Scale-invariant critical dy- namics at eigenstate transitions, Phys. Rev. Res. 5, 043301 (2023)
2023
-
[25]
De Tomasi, I
G. De Tomasi, I. M. Khaymovich, F. Pollmann, and S. Warzel, Rare thermal bubbles at the many-body lo- calization transition from the Fock space point of view, Phys. Rev. B 104, 024202 (2021)
2021
-
[26]
Tarzia, Many-body localization transition in Hilbert space, Phys
M. Tarzia, Many-body localization transition in Hilbert space, Phys. Rev. B 102, 014208 (2020)
2020
-
[27]
Hopjan and L
M. Hopjan and L. Vidmar, Scale-invariant survival prob- ability at eigenstate transitions, Phys. Rev. Lett. 131, 060404 (2023)
2023
-
[28]
Altshuler and V
B. Altshuler and V. Kravtsov, Random Cantor sets and mini-bands in local spectrum of quantum systems, An- nals of Physics , 169300 (2023), in press
2023
-
[29]
Sarkar, R
M. Sarkar, R. Ghosh, and I. M. Khaymovich, Tuning the phase diagram of a rosenzweig-porter model with fractal disorder, Phys. Rev. B 108, L060203 (2023)
2023
-
[30]
Here we don’t consider the models with deterministic quasiperiodic potentials as they are crucially sensitive to random perturbations and immediately undergo the localization transition for any random perturbation
-
[31]
Rosenzweig and C
N. Rosenzweig and C. E. Porter, ”Repulsion of energy levels” in complex atomic spectra, Phys. Rev. B 120, 1698 (1960)
1960
-
[32]
Kutlin and C
A. Kutlin and C. Vanoni, Investigating finite-size ef- fects in random matrices by counting resonances, SciPost Phys. 18, 090 (2025)
2025
-
[33]
Here and further, we neglect all the L-dependencies slower than any power and look at the model in terms of the saddle-point approximation
-
[34]
See Supplemental Material at [URL will be inserted by publisher] for technical details, which includes Refs. []
-
[35]
De Tomasi and I
G. De Tomasi and I. M. Khaymovich, Non-Hermitian Rosenzweig-Porter random-matrix ensemble: Obstruc- tion to the fractal phase, Phys. Rev. B 106, 094204 (2022)
2022
-
[36]
A. K. Das, A. Ghosh, and I. M. Khaymovich, Absence of mobility edge in short-range uncorrelated disordered model: Coexistence of localized and extended states, Phys. Rev. Lett. 131, 166401 (2023)
2023
-
[37]
Barney, M
R. Barney, M. Winer, C. L. Baldwin, B. Swingle, and V. Galitski, Spectral statistics of a minimal quantum glass model, SciPost Phys. 15, 084 (2023)
2023
-
[38]
Izrailev, S
F. Izrailev, S. Ruffo, and L. Tessieri, Classical represen- tation of the one-dimensional Anderson model, Journal of physics A: Mathematical and general 31, 5263 (1998). 7
1998
-
[39]
Tang and I
W. Tang and I. M. Khaymovich, Non-ergodic delocal- ized phase with Poisson level statistics, Quantum 6, 733 (2022)
2022
-
[40]
Appendix A: Details of F ractal RP model Let us provide a brief description of the fractal RP model
Note that b >−2k >0 is not achievable, because |hm − hm±1| ≤O(1). Appendix A: Details of F ractal RP model Let us provide a brief description of the fractal RP model. In this special RP model, we consider the diago- nal elements (hn) to be chosen from certain (multi)fractal di...
Reviewed August 12, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.