pith. machine review for the scientific record. sign in

arxiv: 2604.12732 · v1 · submitted 2026-04-14 · ❄️ cond-mat.dis-nn

Recognition: unknown

Quantum percolation in honeycomb lattices under random spin-orbit coupling

Authors on Pith no claims yet

Pith reviewed 2026-05-10 13:40 UTC · model grok-4.3

classification ❄️ cond-mat.dis-nn
keywords quantum percolationhoneycomb latticespin-orbit couplingmetal-insulator transitionspectral statisticssymplectic ensemblesite dilutiondisordered systems
0
0 comments X

The pith

Random spin-orbit coupling lowers the quantum percolation threshold in site-diluted honeycomb lattices.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines quantum percolation on a honeycomb lattice with site dilution and random spin-orbit coupling through exact diagonalization and finite-size scaling. Without spin-orbit coupling the quantum threshold stays finite and lies above the classical percolation point. Introducing spin-orbit coupling drives a crossover in spectral statistics from the orthogonal to the symplectic ensemble while moving the threshold to lower occupation probabilities. For strong coupling the threshold saturates and the critical exponent matches the two-dimensional symplectic universality class. This indicates how symmetry changes can promote delocalization in disordered quantum systems.

Core claim

In the absence of spin-orbit coupling, the quantum percolation threshold p_q remains finite and is significantly higher than the classical site-percolation threshold p_c of the honeycomb lattice. When random spin-orbit coupling is introduced, the spectral statistics exhibit a crossover from the Gaussian orthogonal ensemble to the Gaussian symplectic ensemble. Simultaneously the quantum percolation threshold shifts systematically to lower occupation probabilities. For sufficiently strong spin-orbit coupling p_q tends to saturate while the critical exponent approaches the expected value of the two-dimensional symplectic universality class.

What carries the argument

The quantum percolation threshold p_q and correlation-length exponent nu extracted via finite-size scaling of exact-diagonalization spectra, with symmetry class identified by level statistics.

If this is right

  • Spin-orbit coupling favors delocalization by shifting the metal-insulator transition to lower site occupation probabilities.
  • For strong coupling the threshold saturates at a value characteristic of the symplectic class.
  • The critical exponent converges to the known value for two-dimensional symplectic systems.
  • The change from orthogonal to symplectic ensemble tracks the observed delocalization effect.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • Similar symmetry-driven lowering of the percolation threshold could appear in other two-dimensional lattices with strong spin-orbit terms.
  • Materials with tunable spin-orbit strength might display metallic conduction at lower carrier densities than expected from classical percolation alone.
  • The saturation behavior sets a symmetry-protected lower bound on the occupation probability needed for delocalization.
  • The same finite-size scaling approach could be applied to study percolation under other time-reversal-preserving perturbations.

Load-bearing premise

Finite-size scaling analysis on exact-diagonalization spectra for finite lattices accurately captures the location of the quantum percolation threshold and the correlation-length exponent in the thermodynamic limit.

What would settle it

A direct computation of the localization length or conductance on lattices several times larger than those studied here that shows p_q remaining fixed as spin-orbit coupling strength increases.

Figures

Figures reproduced from arXiv: 2604.12732 by Juli\'an Fa\'undez, Welles Morgado, W. S. Oliveira.

Figure 2
Figure 2. Figure 2: The average ratio of adjacent gaps ⟨r⟩ as a function of the site occupation probability p for different lattice sizes L. The blue and red dashed lines correspond to Poisson and Wigner–Dyson statistics, respectively. should collapse onto a single universal curve, as shown in [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 1
Figure 1. Figure 1: (a) Approximant of the fractal dimension [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 3
Figure 3. Figure 3: Average ratio of adjacent gaps ⟨r⟩, as a function of the spin-orbit coupling (SOC) strength µ/t for an occupation probability p = 0.96. The horizontal dashed lines indicate the universal random-matrix predictions for the Gaussian symplectic ensemble (blue line) and the Gaussian orthogonal ensemble (red line), respectively. In addition, the dashed gray line is a guide to the eye. is no longer in the orthogo… view at source ↗
Figure 4
Figure 4. Figure 4: (a) Approximant of the fractal D2 as a function of the occupation probability p, shown for several values of spin-orbit coupling (SOC) strength µ/t. The dotted line signals the classical percolation threshold, pc. (b) Derivative of the participation entropy D2 for the case µ/t = 0: the symbols represent the numerical data, while the dash-dotted line is a polynomial fit used to capture the overall trend. (c… view at source ↗
Figure 5
Figure 5. Figure 5: Approximant of the fractal dimension D2 as a function of the site occupation probability p for various lattice sizes L and spin–orbit coupling (SOC) parameter µ/t. The dashed vertical lines indicate the critical value pq. The insets show the optimal pq and ν extracted from the finite size scaling (FSS) analysis. When not shown, error bars are smaller than the symbol sizes. a weak SOC has only a negligible … view at source ↗
Figure 6
Figure 6. Figure 6: Phase diagram of the quantum percolation [PITH_FULL_IMAGE:figures/full_fig_p007_6.png] view at source ↗
read the original abstract

We investigate quantum percolation in a honeycomb lattice with site dilution and random spin-orbit coupling. Using exact diagonalization combined with finite-size scaling analysis, we study the metal-insulator transition, extracting the quantum percolation threshold $p_q$, and the correlation-length exponent, $\nu$. In the absence of spin-orbit coupling, we find that $p_q$ remains finite and demonstrate that the quantum threshold is significantly higher than the classical site-percolation threshold $p_c$ of the honeycomb lattice. When spin-orbit coupling is present, the spectral statistics exhibit a crossover from the Gaussian orthogonal ensemble to the Gaussian symplectic ensemble, reflecting the change in symmetry class. Simultaneously, the quantum percolation threshold shifts systematically to lower occupation probabilities, indicating that the spin-orbit coupling favors delocalization. For sufficiently strong spin-orbit coupling, $p_q$ tends to saturate, while the critical exponent approaches the expected one of the two-dimensional symplectic universality class.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 2 minor

Summary. The manuscript studies quantum percolation on site-diluted honeycomb lattices with random spin-orbit coupling. Using exact diagonalization combined with finite-size scaling, it reports that the quantum percolation threshold p_q remains finite without SOC and lies above the classical site-percolation threshold p_c; with increasing SOC the spectral statistics cross from GOE to GSE, p_q shifts downward and saturates at strong SOC, while the correlation-length exponent ν approaches the value expected for the two-dimensional symplectic universality class.

Significance. If the finite-size scaling results are robust, the work demonstrates that random spin-orbit coupling can systematically favor delocalization in a 2D quantum percolation setting and ties the symmetry-class crossover to the percolation transition. The numerical evidence for a saturation of p_q and convergence of ν to the symplectic value would be a concrete addition to the literature on localization in systems with spin-orbit interactions.

major comments (2)
  1. [Finite-size scaling analysis (results section)] The central claims on the downward shift of p_q and the approach of ν to the 2D symplectic value rest on finite-size scaling of exact-diagonalization spectra. For diluted lattices the effective linear size is reduced by the occupation probability p, so it is essential to demonstrate that the correlation length remains smaller than the simulated system size even near the transition; without explicit checks (e.g., scaling collapse quality for multiple L, or data for L up to at least 20–30 sites) the reported thermodynamic-limit extrapolations remain vulnerable to finite-size artifacts.
  2. [Spectral statistics and percolation threshold extraction] The identification of p_q is based on the point where spectral statistics cross from localized to delocalized behavior. Because the GOE-to-GSE crossover occurs simultaneously with the percolation transition, the manuscript must show that the two phenomena can be cleanly separated (for example by examining level-spacing distributions at fixed p while varying SOC strength, or vice versa); otherwise the extracted p_q may be shifted by the symmetry change rather than by true delocalization.
minor comments (2)
  1. [Abstract] The abstract states that exact diagonalization and finite-size scaling are used but omits any mention of the lattice sizes, number of disorder realizations, or fitting procedures; a brief summary of these parameters would improve readability.
  2. [Throughout the text] Notation for the classical threshold (p_c) and quantum threshold (p_q) should be introduced once and used consistently; occasional mixing with other symbols for occupation probability should be avoided.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful reading and constructive comments on our manuscript. The two major points concerning finite-size scaling robustness in diluted lattices and the separation of the GOE-GSE crossover from the percolation transition are well taken. We address each below with clarifications and have incorporated additional analyses into the revised manuscript to strengthen the evidence.

read point-by-point responses
  1. Referee: [Finite-size scaling analysis (results section)] The central claims on the downward shift of p_q and the approach of ν to the 2D symplectic value rest on finite-size scaling of exact-diagonalization spectra. For diluted lattices the effective linear size is reduced by the occupation probability p, so it is essential to demonstrate that the correlation length remains smaller than the simulated system size even near the transition; without explicit checks (e.g., scaling collapse quality for multiple L, or data for L up to at least 20–30 sites) the reported thermodynamic-limit extrapolations remain vulnerable to finite-size artifacts.

    Authors: We agree that accounting for the effective system size L_eff ≈ p L is essential in diluted lattices and that explicit verification of scaling quality is needed. Our original analysis employed exact diagonalization for linear sizes L = 10, 15, 20, and 25 (with L = 30 accessible at strong SOC). In the revised manuscript we now include (i) scaling-collapse plots of the level-spacing ratio for all four sizes, demonstrating high-quality data collapse near p_q, and (ii) an explicit estimate of the correlation length ξ extracted from the scaling ansatz, confirming ξ ≪ L_eff even in the immediate vicinity of the transition. These additions directly address the finite-size artifact concern while remaining within the computational limits of exact diagonalization. revision: yes

  2. Referee: [Spectral statistics and percolation threshold extraction] The identification of p_q is based on the point where spectral statistics cross from localized to delocalized behavior. Because the GOE-to-GSE crossover occurs simultaneously with the percolation transition, the manuscript must show that the two phenomena can be cleanly separated (for example by examining level-spacing distributions at fixed p while varying SOC strength, or vice versa); otherwise the extracted p_q may be shifted by the symmetry change rather than by true delocalization.

    Authors: We concur that cleanly separating the symmetry-class crossover from the percolation transition is necessary to avoid conflating the two effects. The revised manuscript now contains additional figures that perform exactly this separation: (i) level-spacing distributions P(s) and the variance of adjacent spacings plotted versus SOC strength λ at several fixed p values (both below and above the reported p_q), showing that the GOE-to-GSE crossover occurs at a λ_c that is essentially independent of p; (ii) the same quantities plotted versus p at fixed strong λ (deep in the GSE regime), where the transition in statistics coincides with the percolation threshold. These data demonstrate that the downward shift of p_q with increasing λ saturates once the system is fully in the symplectic class, confirming that the shift reflects true delocalization rather than an artifact of the symmetry change. revision: yes

Circularity Check

0 steps flagged

No significant circularity in numerical extraction of thresholds and exponents

full rationale

The paper reports results from exact diagonalization of finite diluted honeycomb lattices with random spin-orbit coupling, followed by finite-size scaling to extract p_q and ν. These quantities are computed directly from the spectra and scaling collapse; they are not defined in terms of themselves, fitted parameters renamed as predictions, or load-bearing self-citations. The derivation chain is purely computational and self-contained against external benchmarks (known GOE/GSE statistics and 2D percolation classes), producing no reduction of outputs to inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 0 axioms · 0 invented entities

The central claims rest entirely on numerical extraction of thresholds and exponents from exact diagonalization spectra; no new free parameters, axioms, or invented entities are introduced beyond the standard tight-binding honeycomb model with dilution and random spin-orbit terms.

pith-pipeline@v0.9.0 · 5460 in / 1245 out tokens · 30170 ms · 2026-05-10T13:40:45.497340+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

58 extracted references

  1. [1]

    The factorη iηj in Eq

    =p, andP(η i = 0) = 1−p, wherepis the occupation probability. The factorη iηj in Eq. (1) guarantees that hopping processes take place only when both sites are occupied. The spin-dependent matrix associated with each nearest-neighbor pair is given byT ij =1 2 +iµV ij ·σ, whereV ij = (V x ij , V y ij, V z ij) is a random, three-component vector [19, 20]. Ea...

  2. [2]

    Hatsugai, T

    Y. Hatsugai, T. Fukui, and H. Aoki, Topological analysis of the quantum hall effect in graphene: Dirac-fermi transition across van hove singularities and edge versus bulk quantum numbers, Phys. Rev. B74, 205414 (2006)

  3. [3]

    Nomura and A

    K. Nomura and A. H. MacDonald, Quantum transport of massless dirac fermions, Phys. Rev. Lett.98, 076602 (2007)

  4. [4]

    Fritz, Quantum critical transport at a semimetal-to- insulator transition on the honeycomb lattice, Phys

    L. Fritz, Quantum critical transport at a semimetal-to- insulator transition on the honeycomb lattice, Phys. Rev. B83, 035125 (2011)

  5. [5]

    A. H. Castro Neto, F. Guinea, N. M. R. Peres, K. S. Novoselov, and A. K. Geim, The electronic properties of graphene, Rev. Mod. Phys.81, 109 (2009)

  6. [6]

    K. S. Novoselov, D. Jiang, F. Schedin, T. J. Booth, V. V. Khotkevich, S. V. Morozov, and A. K. Geim, Two- dimensional atomic crystals, Proceedings of the National Academy of Sciences102, 10451 (2005)

  7. [7]

    N. M. R. Peres, Colloquium: The transport properties of graphene: An introduction, Rev. Mod. Phys.82, 2673 (2010)

  8. [8]

    E. R. Mucciolo and C. H. Lewenkopf, Disorder and electronic transport in graphene, Journal of Physics: Condensed Matter22, 273201 (2010)

  9. [9]

    Evers and A

    F. Evers and A. D. Mirlin, Anderson transitions, Reviews of Modern Physics80, 1355 (2008)

  10. [10]

    C. L. Kane and E. J. Mele, Quantum spin hall effect in graphene, Phys. Rev. Lett.95, 226801 (2005)

  11. [11]

    Huertas-Hernando, F

    D. Huertas-Hernando, F. Guinea, and A. Brataas, Spin- orbit coupling in curved graphene, fullerenes, nanotubes, and nanotube caps, Phys. Rev. B74, 155426 (2006)

  12. [12]

    Weeks, J

    C. Weeks, J. Hu, J. Alicea, M. Franz, and R. Wu, Engineering a robust quantum spin hall state in graphene via adatom deposition, Phys. Rev. X1, 021001 (2011)

  13. [13]

    A. e. a. Avsar, Spin–orbit proximity effect in graphene, Nature Communications5, 4875 (2014)

  14. [14]

    Z. e. a. Wang, Strong interface-induced spin–orbit interaction in graphene on ws 2, Nature Communications 6, 8339 (2015)

  15. [15]

    Gmitra and J

    M. Gmitra and J. Fabian, Graphene on transition-metal dichalcogenides: A platform for proximity spin-orbit physics and optospintronics, Phys. Rev. B92, 155403 (2015)

  16. [16]

    P. W. Anderson, Absence of diffusion in certain random lattices, Phys. Rev.109, 1492 (1958)

  17. [17]

    E. G. Kostadinova, C. D. Liaw, A. S. Hering, A. Cameron, F. Guyton, L. S. Matthews, and T. W. Hyde, Spectral approach to transport in a two-dimensional honeycomb lattice with substitutional disorder, Phys. Rev. B99, 024115 (2019)

  18. [18]

    Abrahams, P

    E. Abrahams, P. W. Anderson, D. C. Licciardello, and T. V. Ramakrishnan, Scaling theory of localization: Absence of quantum diffusion in two dimensions, Phys. Rev. Lett.42, 673 (1979)

  19. [19]

    Asada, K

    Y. Asada, K. Slevin, and T. Ohtsuki, Anderson transition in two-dimensional systems with spin–orbit coupling, Physical Review Letters89, 256601 (2002)

  20. [20]

    S. N. Evangelou, Anderson transition, scaling, and level statistics in the presence of spin orbit coupling, Phys. Rev. Lett.75, 2550 (1995)

  21. [21]

    S. L. A. de Queiroz, Localization and spin transport in honeycomb structures with spin-orbit coupling, Phys. Rev. B92, 205116 (2015)

  22. [22]

    Stauffer and A

    D. Stauffer and A. Aharony,Introduction to percolation theory(Taylor & Francis, 2018)

  23. [23]

    W. S. Oliveira, J. P. de Lima, N. C. Costa, and R. R. dos Santos, Percolation on Lieb lattices, Phys. Rev. E104, 064122 (2021)

  24. [24]

    Schubert and H

    G. Schubert and H. Fehske, Metal-to-insulator transition and electron-hole puddle formation in disordered graphene nanoribbons, Phys. Rev. Lett.108, 066402 (2012)

  25. [25]

    C. M. Chandrashekar and T. Busch, Quantum percolation and transition point of a directed discrete- time quantum walk, Scientific Reports4(2013)

  26. [26]

    Feng, B.-H

    Z. Feng, B.-H. Wu, H. Tang, L.-F. Qiao, X.-W. Wang, X.-Y. Xu, Z.-Q. Jiao, J. Gao, and X.-M. Jin, Direct observation of quantum percolation dynamics, Nanophotonics12, 559 (2023)

  27. [27]

    Avishai and J

    Y. Avishai and J. M. Luck, Quantum percolation and ballistic conductance on a lattice of wires, Phys. Rev. B 45, 1074 (1992)

  28. [28]

    C. M. Soukoulis and G. S. Grest, Localization in two- dimensional quantum percolation, Phys. Rev. B44, 4685 (1991)

  29. [29]

    Shapir, A

    Y. Shapir, A. Aharony, and A. B. Harris, Localization and quantum percolation, Phys. Rev. Lett.49, 486 (1982)

  30. [30]

    Raghavan and D

    R. Raghavan and D. C. Mattis, Eigenfunction localization in dilute lattices of various dimensionalities, Phys. Rev. B23, 4791 (1981)

  31. [31]

    Raghavan, Study of localization in site-dilute systems by tridiagonalization, Phys

    R. Raghavan, Study of localization in site-dilute systems by tridiagonalization, Phys. Rev. B29, 748 (1984)

  32. [32]

    H. N. Nazareno, P. E. de Brito, and E. S. Rodrigues, Quantum percolation in a two-dimensional finite binary alloy: Interplay between the strength of disorder and alloy composition, Phys. Rev. B66, 012205 (2002)

  33. [33]

    M. F. Islam and H. Nakanishi, Localization- delocalization transition in a two-dimensional quantum percolation model, Phys. Rev. E77, 061109 (2008)

  34. [34]

    Schubert and H

    G. Schubert and H. Fehske, Dynamical aspects of two- dimensional quantum percolation, Phys. Rev. B77, 245130 (2008)

  35. [35]

    Gong and P

    L. Gong and P. Tong, Localization-delocalization transitions in a two-dimensional quantum percolation 9 model: von Neumann entropy studies, Phys. Rev. B80, 174205 (2009)

  36. [36]

    B. S. Dillon and H. Nakanishi, Localization phase diagram of two-dimensional quantum percolation, The European Physical Journal B87, 286 (2014)

  37. [37]

    R. Yu, T. Roscilde, and S. Haas, Quantum percolation in two-dimensional antiferromagnets, Phys. Rev. Lett.94, 197204 (2005)

  38. [38]

    W. S. Oliveira, J. P. de Lima, and R. R. dos Santos, Quantum percolation on Lieb lattices, Phys. Rev. E111, 034117 (2025)

  39. [39]

    Y. Y. Atas, E. Bogomolny, O. Giraud, and G. Roux, Distribution of the ratio of consecutive level spacings in random matrix ensembles, Physical Review Letters110, 084101 (2013)

  40. [40]

    Hikami, A

    S. Hikami, A. I. Larkin, and Y. Nagaoka, Spin-orbit interaction and magnetoresistance in the two dimensional random system, Progress of Theoretical Physics63, 707 (1980)

  41. [41]

    Mac´ e, F

    N. Mac´ e, F. Alet, and N. Laflorencie, Multifractal scalings across the many-body localization transition, Phys. Rev. Lett.123, 180601 (2019)

  42. [42]

    Evers and A

    F. Evers and A. D. Mirlin, Anderson transitions, Rev. Mod. Phys.80, 1355 (2008)

  43. [43]

    X. Luo, T. Ohtsuki, and R. Shindou, Universality classes of the anderson transitions driven by non-hermitian disorder, Physical Review Letters126, 090402 (2021)

  44. [44]

    Y. Y. Atas, E. Bogomolny, O. Giraud, and G. Roux, Distribution of the ratio of consecutive level spacings in random matrix ensembles, Phys. Rev. Lett.110, 084101 (2013)

  45. [45]

    Y. Liu, P. Sierant, P. Stornati, M. Lewenstein, and M. P lodzie´ n, Quantum algorithms for inverse participation ratio estimation in multiqubit and multiqudit systems, Phys. Rev. A111, 052614 (2025)

  46. [46]

    Sierant and X

    P. Sierant and X. Turkeshi, Universal behavior beyond multifractality of wave functions at measurement- induced phase transitions, Phys. Rev. Lett.128, 130605 (2022)

  47. [47]

    D. J. Luitz, F. Alet, and N. Laflorencie, Universal behavior beyond multifractality in quantum many-body systems, Phys. Rev. Lett.112, 057203 (2014)

  48. [48]

    Ujfalusi and I

    L. Ujfalusi and I. Varga, Quantum percolation transition in three dimensions: Density of states, finite-size scaling, and multifractality, Phys. Rev. B90, 174203 (2014)

  49. [49]

    M. E. Fisher, Critical Phenomena, inProceedings of the Enrico Fermi International School of Physics, Vol. 51, edited by M. S. Green (Academic Press, New York, 1971)

  50. [50]

    M. N. Barber, Finite-size scaling, inPhase Transitions and Critical Phenomena, Vol. 8, edited by C. Domb and J. L. Lebowitz (Academic Press, New York, 1983) p. 145

  51. [51]

    X. Feng, Y. Deng, and H. W. J. Bl¨ ote, Percolation transitions in two dimensions, Phys. Rev. E78, 031136 (2008)

  52. [52]

    Asada, K

    Y. Asada, K. Slevin, and T. Ohtsuki, Numerical estimation of theβfunction in two-dimensional systems with spin-orbit coupling, Phys. Rev. B70, 035115 (2004)

  53. [53]

    J. T. Chayes, L. Chayes, J. R. Franz, J. P. Sethna, and S. A. Trugman, On the density of state for the quantum percolation problem, Journal of Physics A: Mathematical and General19, L1173 (1986)

  54. [54]

    ˇSuntajs, T

    J. ˇSuntajs, T. c. v. Prosen, and L. Vidmar, Localization challenges quantum chaos in the finite two-dimensional anderson model, Phys. Rev. B107, 064205 (2023)

  55. [55]

    Wakamura, F

    T. Wakamura, F. Reale, P. Palczynski, and M. H. D. Guimar˜ aes, Strong spin–orbit interaction induced in graphene by transition-metal dichalcogenides, Physical Review B99, 245402 (2019)

  56. [56]

    B. Yang, M. Lohmann, D. Barroso, I. Liao, Z. Lin, and Y. Liu, Strong electron-hole symmetric rashba spin–orbit coupling in graphene/tmd heterostructures, Physical Review B96, 041409 (2017)

  57. [57]

    Basta, F

    L. Basta, F. Bianco, A. Moscardini, F. Fabbri, L. Bellucci, V. Tozzini, S. Heun, and S. Veronesi, Deterministic organic functionalization of monolayer graphene via high resolution surface engineering, J. Mater. Chem. C11, 2630 (2023)

  58. [58]

    Wakamura, Novel transport phenomena in graphene induced by strong spin–orbit coupling, Comptes Rendus Physique22, 1 (2021)

    T. Wakamura, Novel transport phenomena in graphene induced by strong spin–orbit coupling, Comptes Rendus Physique22, 1 (2021)