REVIEW 3 major objections 5 minor 21 references
Bitstring probabilities for Rydberg ladder vacua collapse onto a Fermi-like curve, and the shots needed to tame low-probability states grow exponentially with atom number.
Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →
T0 review · grok-4.5
2026-07-30 14:01 UTC pith:3EPXLQZX
load-bearing objection Useful empirical shot-budget diagnostic for Rydberg-ladder vacua; Fermi collapse is a fit, not a law, and the device-scale extrapolations ride on binning-sensitive slopes. the 3 major comments →
Finite size scaling of bitstring probability distributions for Rydberg arrays
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
For the vacuum of Rydberg ladders the cumulative probability Σ(p_Λ, N_q) of all bitstrings with probability at most p_Λ can be approximately collapsed, for moderate p_Λ, onto a single Fermi-Dirac-like function of −ln(p_Λ). Both the largest single probability and the cutoff at which Σ reaches one-half decay exponentially with atom number N_q, so the number of shots needed to keep the unresolved tail below a chosen accuracy likewise grows exponentially with N_q.
What carries the argument
The cumulative probability distribution Σ(p_Λ, N_q) = sum of all p_{n} ≤ p_Λ. It converts the raw bitstring histogram into a single curve whose mid-region collapse and inflection-point scaling supply the finite-size shot-cost law.
Load-bearing premise
The exponential scalings and Fermi-like collapse fitted on ladders of at most a few dozen atoms with a billion DMRG samples continue to hold when extrapolated to the sixty-to-one-hundred-twenty qubit regime.
What would settle it
Compute or measure Σ(p_Λ) for a Rydberg ladder larger than those already studied (for example 60 atoms) at the same blockade ratio; if the mid-region no longer collapses onto the same Fermi form or if the half-probability cutoff fails to follow the reported exponential, the scaling claim is false.
If this is right
- Shot budgets required to keep vacuum observables accurate can be read off from the exponential fit of p* versus N_q rather than guessed.
- The same collapse supplies a practical error bar: choose the desired tail weight Σ = z and read the minimal shots as 1/p_Λ(z).
- Region-I maximal probabilities already limit useful system size to roughly 120 qubits at a billion shots; Region-II half-filling tightens that to roughly 60 qubits.
- Because the exponential grows slower than the Hilbert-space dimension 2^{N_q}, larger devices remain usable if shots are allocated intelligently.
Where Pith is reading between the lines
- If the Fermi collapse is not an artifact of the Rydberg blockade, the same mid-region scaling may appear in other gapped lattice models whose ground-state bitstring spectra are similarly sparse.
- The documented sensitivity of the slope β to binning suggests that an exact, unbinned definition of Σ would tighten the extrapolated qubit limits.
- Observables that are dominated by the highest-probability bitstrings may still be accessible beyond the 60-qubit half-filling limit, while observables sensitive to the tail will hit the wall earlier.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper studies bitstring probability distributions of the ground state (vacuum) of Rydberg atom ladders obtained from large DMRG/MPS sampling (10^9 shots). It introduces the cumulative distribution Σ(p_Λ,N_q), the total probability mass on bitstrings with p≤p_Λ, and reports that for intermediate p_Λ the curves for successive system sizes N_q can be approximately collapsed onto a single Fermi–Dirac-like function of −ln(p_Λ). From exponential fits of p_max, p*_1/2 and related anchors versus N_q (up to ~44 atoms at fixed ρ=2, R_b/a=2.00), the authors conclude that the number of shots needed to drive Σ below a fixed small threshold grows exponentially with N_q, with rough extrapolations to ~60–120 qubits, and discuss implications for vacuum observables.
Significance. If the mid-region collapse and the exponential shot-cost scaling hold beyond the fitted window, the work would give a practical, system-size-aware diagnostic for how many measurements are required to control the low-probability tail when estimating vacuum observables on Rydberg simulators—an issue of direct interest for analog quantum simulation of many-body and lattice-gauge models. Strengths include the concrete construction of large bitstring dictionaries from DMRG, the clear empirical collapse of Σ in the mid-probability window (Figs. 1–2), and the explicit connection of cumulative mass to shot budgets. The Fermi identification and the device-scale extrapolations are, however, phenomenological fits rather than derived scaling laws, so the significance is that of a useful empirical organizing principle whose range of validity still needs to be delimited.
major comments (3)
- [Region II, Eqs. (3)–(4); Appendix; Fig. 4] Abstract, Region II (Eqs. 3–4) and Conclusion: the central claim that Σ collapses onto a Fermi function of −ln(p_Λ) and that the shot number needed to reduce Σ to a fixed small value grows exponentially with N_q rests on two-parameter phenomenological fits. The Appendix (“Dependence of Parameter Fitting”) and Figs. 8–9 document that the slope β (and the Region-III exponent B) vary substantially with histogram binning and fit-window choice. Because the same β and the fitted anchors p*, p*_1/2, p*_0.01 are used both to collapse the curves and to produce the exponential extrapolations in Fig. 4 (~60–120 qubits), the load-bearing shot-cost statement inherits that variance. The manuscript should either (i) adopt an exact (unbinned) cumulative construction and report error bands on β and on the extrapolated N_q, or (ii) clearly demote the Fermi form and the 60–120-qubit numbers to illustrative
- [Region III–IV; Fig. 4; Conclusion] All quantitative results are obtained at a single Hamiltonian point (ρ=2, R_b/a=2.00) for ladders up to ~44 atoms, with Region III/IV explicitly truncated by the 10^9-shot floor (text around Figs. 5–7). The exponential persistence of p_max, p*_1/2 and p*_0.01 outside this window, and the claim that the same collapse governs the low-p tail that dominates shot cost, are therefore uncontrolled extrapolations. At minimum the paper should test one additional detuning/blockade point or geometry, and/or push the single-occurrence floor with a substantially larger shot budget on a subset of sizes, so that the device-scale statement is not read off a single finite-shot fit.
- [Region III, Eqs. (5)–(6)] Region III (Eq. 5) and the comparison (Eq. 6) treat the power-law form and the Fermi form as interchangeable via B≃β and A≃(p*)^{-β}. Given the documented binning dependence of both exponents, this identification is not robust enough to underwrite the error-estimate procedure advocated in the text (Σ=z as a shot-budget proxy). Either the equivalence should be dropped or it should be replaced by a binning-independent statement (e.g., direct interpolation of the empirical Σ without parametric collapse).
minor comments (5)
- [Eq. (2); Abstract] Notation for the cumulative is inconsistent: Σ(p_Λ) in Eq. (2) versus Σ(p_Λ,N_q) in the abstract and later text; fix one form throughout.
- [Figure 1] Figure 1 caption and main text refer to “4 rungs to 22 rungs” while elsewhere sizes are given in atoms (20–44); state N_q explicitly in every figure caption.
- [Introduction / Hamiltonian] The Hamiltonian parameters Ω and Δ used for the DMRG ground states are never stated numerically; a short table or sentence would make the dictionaries reproducible.
- [Region I–II] Typos and wording: “these large are due to” (Region I); “Taking −ln(p_Λ)→ϵ and →−ln(p*)→μ” is missing the left-hand variable; “ap ∗” spacing; “finitely discrete states” is unclear.
- [Introduction] Self-citations [8,21] supply related diagnostics; a brief sentence distinguishing the new Σ collapse from those earlier filtered-probability tools would help the reader.
Circularity Check
Empirical DMRG fits and constructed collapse; mild reuse of fitted p*/β as scaling predictors, not definitional circularity.
specific steps
-
fitted input called prediction
[Region II, Eqs. (3)–(4), Fig. 4; Conclusion]
"Σ(p_Λ)≃ (p_Λ/p*)^β/(1+(p_Λ/p*)^β). Taking −ln(p_Λ)→ϵ and −ln(p*)→μ the cumulative probability distributions can be well-approximated by Fermi-Dirac... We show that the number of shots necessary to reduce Σ(p_Λ,N_q) to some low enough value grows exponentially with N_q."
p* and β are fitted per N_q to the mid-region of Σ; the collapse and the Fermi resemblance are then presented as structural findings, and the exponential shot-cost claim is the same exponential fit of those p* values (Fig. 4) restated as N_sh∼1/p*. The prediction is statistically forced by the fit rather than independently derived. Content remains empirical (new DMRG dictionaries), so this is mild, not definitional.
-
self definitional
[COLLAPSE OF Σ; Fig. 2 construction]
"The distributions of different system sizes will collapse by taking p_Λ to (p_Λ/p*)^B where B is the slope of the line fit to each individual distribution and then plotting the distributions against this adjusted independent variable."
The collapse coordinate is defined from the same per-curve p* and fitted slope B that characterize each Σ. Any family of similar sigmoids will overlie under that rescaling by construction; the residual empirical claim is only that shapes (β) are similar enough across N_q. Partial self-definition of the collapse, not of the underlying Σ data.
full rationale
The paper’s load-bearing content is new numerical sampling: bitstring dictionaries from DMRG/MPS (10^9 shots) on Rydberg ladders, from which Σ(p_Λ,N_q) is computed directly. The mid-region Fermi-like form (Eqs. 3–4), the collapse via per-curve p* and β, and the exponential fits of p_max and p*_1/2 vs N_q (Fig. 4) are empirical observations and extrapolations, not algebraic rewrites of prior equations into themselves. Self-citations ([8], [21]) supply related diagnostics and motivation; they do not force the scaling claims. Mild circularity appears only in the usual fit-then-read-off sense: the same fitted p* and slope used to construct the collapse are then cited as evidence of universal Fermi behavior and of exponential shot cost (N_sh ∼ 1/p*). That is proportionate to score 2, not a self-definitional or self-citation chain. No uniqueness theorem, smuggled ansatz, or renaming of a known closed-form result is load-bearing.
Axiom & Free-Parameter Ledger
free parameters (6)
- Fermi slope β (and collapse exponent B) =
O(1), system- and protocol-dependent (Figs. 8–9)
- Anchor probabilities p*, p*_1/2, p*_0.01 =
read off per N_q from Σ curves
- Power-law prefactor A and exponent B in Region III
- Histogram binning and fit-window size =
e.g. 300 bins, ±20 points (appendix)
- Sampling shot budget N_sh = 10^9 =
10^9
- Rydberg ladder parameters (ρ=2, R_b/a=2.00, Ω, Δ implied) =
ρ=2, R_b/a=2.00
axioms (5)
- domain assumption DMRG/MPS ground states and 10^9 samples faithfully represent the true bitstring probability distribution down to p~10^{-9} for the studied sizes.
- domain assumption The Rydberg Hamiltonian (Eq. 1) with van der Waals 1/r^6 interactions on a two-leg ladder is the correct microscopic model for the vacuum under study.
- ad hoc to paper Mid-region Σ(p_Λ) is well approximated by a two-parameter Fermi–Dirac function of −ln(p_Λ).
- ad hoc to paper Exponential fits of p_max, p*_1/2, and p*_0.01 versus N_q continue outside the fitted window to ~60–120 qubits.
- domain assumption Cumulative probability Σ is the right figure of merit for vacuum observable error from unsampled bitstrings.
invented entities (2)
-
Four phenomenological regions (I–IV) of Σ
no independent evidence
-
Identification of collapsed Σ with a Fermi–Dirac distribution in ε=−ln(p)
no independent evidence
read the original abstract
We calculate the probabilities $p_{\{n\}}$ of the measured bitstrings $\{n\}$ for the vacuum of Rydberg ladders with $N_q$ atoms. As $N_q$ increases, the $p_{\{n\}}$ decrease but become more dense in the low $p$ region raising the possibility that their smallness could be compensated by their large number. The importance of the low probability states can be estimated from the cumulative probability distribution $\Sigma(p_{\Lambda},N_q)$, which is the probability to observe any state having a probability $p\leq p_{\Lambda}$. For not too large values of $p_{\Lambda}$, it is possible to approximately collapse the $\Sigma(p_{\Lambda},N_q)$ for successive $N_q$ into a function resembling the Fermi function when plotted as a function of $-\ln(p_{\Lambda})$. We show that the number of shots necessary to reduce $\Sigma(p_{\Lambda},N_q)$ to some low enough value grows exponentially with $N_q$. We discuss the implications for calculating observables associated with the vacuum.
Figures
Reference graph
Works this paper leans on
-
[1]
E. Altmanet al., Quantum Simulators: Architectures and Opportunities, PRX Quantum2, 017003 (2021), arXiv:1912.06938 [quant-ph]
Pith/arXiv arXiv 2021
-
[2]
C. W. Baueret al., Quantum Simulation for High- Energy Physics, PRX Quantum4, 027001 (2023), arXiv:2204.03381 [quant-ph]
Pith/arXiv arXiv 2023
-
[3]
S. Boixo, S. V. Isakov, V. N. Smelyanskiy, R. Bab- bush, N. Ding, Z. Jiang, M. J. Bremner, J. M. Marti- nis, and H. Neven, Characterizing quantum supremacy in near-term devices, Nature Phys.14, 595 (2018), arXiv:1608.00263 [quant-ph]
Pith/arXiv arXiv 2018
-
[4]
A. L. Shaw, D. K. Mark, J. Choi, R. Finkelstein, P. Scholl, S. Choi, and M. Endres, Experimental Sig- natures of Hilbert-Space Ergodicity: Universal Bitstring Distributions and Applications in Noise Learning, Phys. Rev. X15, 031001 (2025), arXiv:2403.11971 [quant-ph]
Pith/arXiv arXiv 2025
-
[5]
T. Guhr, A. Muller-Groeling, and H. A. Weidenmuller, Random matrix theories in quantum physics: Com- mon concepts, Phys. Rept.299, 189 (1998), arXiv:cond- mat/9707301
arXiv 1998
-
[6]
D. K. Mark, F. Surace, A. Elben, A. L. Shaw, J. Choi, G. Refael, M. Endres, and S. Choi, Maximum En- tropy Principle in Deep Thermalization and in Hilbert- Space Ergodicity, Phys. Rev. X14, 041051 (2024), arXiv:2403.11970 [quant-ph]
Pith/arXiv arXiv 2024
-
[7]
W.-K. Mok, T. Haug, A. L. Shaw, M. Endres, and J. Preskill, Optimal Conversion from Classical to Quan- tum Randomness via Quantum Chaos, Phys. Rev. Lett. 134, 180403 (2025), arXiv:2410.05181 [quant-ph]
Pith/arXiv arXiv 2025
-
[8]
A. Kaufman, M. Asaduzzaman, Z. Ozzello, B. Senseman, J. Corona, and Y. Meurice, Diagnosing Device Perfor- mance in Rydberg-Ladder Gauge Simulators with Cu- mulative Probabilities and Filtered Mutual Information, (2025), arXiv:2507.14128 [quant-ph]
arXiv 2025
-
[9]
F. M. Surace, P. P. Mazza, G. Giudici, A. Lerose, A. Gambassi, and M. Dalmonte, Lattice gauge theories and string dynamics in Rydberg atom quantum simula- tors, Phys. Rev. X10, 021041 (2020), arXiv:1902.09551 [cond-mat.quant-gas]
Pith/arXiv arXiv 2020
-
[10]
S. Notarnicola, M. Collura, and S. Montangero, Real- time-dynamics quantum simulation of (1+1)-dimensional lattice QED with Rydberg atoms, Phys. Rev. Res.2, 013288 (2020), arXiv:1907.12579 [cond-mat.quant-gas]
Pith/arXiv arXiv 2020
-
[11]
A. Celi, B. Vermersch, O. Viyuela, H. Pichler, M. D. Lukin, and P. Zoller, Emerging Two-Dimensional Gauge Theories in Rydberg Configurable Arrays, Phys. Rev. X 10, 021057 (2020), arXiv:1907.03311 [quant-ph]
Pith/arXiv arXiv 2020
-
[12]
Y. Meurice, Theoretical methods to design and test quan- tum simulators for the compact Abelian Higgs model, Phys. Rev. D104, 094513 (2021), arXiv:2107.11366 [quant-ph]
Pith/arXiv arXiv 2021
-
[13]
B. Senseman, Z. Ozzello, K. Heitritter, Y. Meurice, and S. Mrenna, Prolegomena to a hybrid classi- cal/Rydberg simulator for hadronization (QuPyth), (2022), arXiv:2212.02476 [quant-ph]
arXiv 2022
-
[14]
P. Fromholz, M. Tsitsishvili, M. Votto, M. Dalmonte, A. Nersesyan, and T. Chanda, Phase diagram of Rydberg-dressed atoms on two-leg triangular ladders, Phys. Rev. B106, 155411 (2022), arXiv:2207.00385 [cond-mat.stat-mech]
Pith/arXiv arXiv 2022
-
[15]
D. Gonzalez-Cuadraet al., Observation of string breaking on a (2 + 1)D Rydberg quantum simulator, Nature642, 321 (2025), arXiv:2410.16558 [quant-ph]
Pith/arXiv arXiv 2025
-
[16]
J. C. Halimeh, M. Aidelsburger, F. Grusdt, P. Hauke, and B. Yang, Cold-atom quantum simulators of gauge theories, Nature Phys.21, 25 (2025), arXiv:2310.12201 [cond-mat.quant-gas]
Pith/arXiv arXiv 2025
-
[17]
D. Gonz´ alez-Cuadra, T. V. Zache, J. Carrasco, B. Kraus, and P. Zoller, Hardware Efficient Quantum Simulation of Non-Abelian Gauge Theories with Qudits on Ryd- berg Platforms, Phys. Rev. Lett.129, 160501 (2022), arXiv:2203.15541 [quant-ph]
Pith/arXiv arXiv 2022
-
[18]
J. Choiet al., Preparing random states and benchmark- ing with many-body quantum chaos, Nature613, 468 (2023), arXiv:2103.03535 [quant-ph]
Pith/arXiv arXiv 2023
-
[19]
D. K. Mark, F. M. Surace, T. Schuster, A. L. Shaw, W. Gong, S. Choi, and M. Endres, Observation of bal- listic plasma and memory in high-energy gauge theory dynamics, (2025), arXiv:2510.11679 [quant-ph]
arXiv 2025
-
[20]
M. D’Arcangelo, L.-P. Henry, L. Henriet, D. Loco, N. Gouraud, S. Angebault, J. Sueiro, J. Forˆ et, P. Mon- march´ e, and J.-P. Piquemal (Fresnel team), Leveraging analog quantum computing with neutral atoms for sol- vent configuration prediction in drug discovery, Phys. Rev. Res.6, 043020 (2024), arXiv:2309.12129 [quant-ph]
Pith/arXiv arXiv 2024
-
[21]
A. Kaufman, J. Corona, Z. Ozzello, B. Senseman, M. Asaduzzaman, and Y. Meurice, Improved entangle- ment entropy estimates from filtered bitstring probabili- ties, Phys. Rev. A112, 032430 (2025), arXiv:2411.07092 [quant-ph]. 7 Dependence of Parameter Fitting When calculating the slopes for the fits in Regions II and III, there are computational choices mad...
Pith/arXiv arXiv 2025
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.