REVIEW 2 major objections 5 minor 8 references
L\'evy $\alpha$-stable generalization of the ReBB model of elastic proton-proton and proton-antiproton scattering
T0 review · 2 major / 5 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read The Lév y α-stable generalization of the Bialas-Bzdak (LBB) model replaces the Gaussian distributions in the ReBB model with Lévy α-stable shapes, recovering ReBB at α_L = 2 and aiming to describe low- and high-|t| elastic pp and p\bar{p}…
desk verdict Clean model extension, but the Lévy parameter α_L<2 violates the unitarity bound on constituent probabilities unless A_ab is re-constrained. 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 Lévy α-stable distribution $L(\vec{x}|\alpha_L, R_L) = \frac{1}{(2\pi)^2} \int d^2\vec{q}\, e^{i\vec{q}\cdot\vec{x}} e^{-|q^2 R_L^2|^{\alpha_L/2}}$, which replaces the Gaussian both in the constituent-constituent inelastic probability (Eq. 2) and in the quark-diquark distribution (Eq. 3). Its defining convolution stability keeps the Glauber expansion algebraically closed: the inelastic probability for two constituents is again a Lévy α-stable function with the combined scale $S_{ab}$. The elastic amplitude is then obtained by averaging over constituent positions and applying the unitary real-part formula $\tilde{T}_{el}(s,b) = i\left(1 - e^{i\alpha_R \tilde{\sigma}_{in}(s,b)}\sqrt{1 - \tilde{\sigma}_{in}(s,b)}\right)$, followed by a Fourier transform to momentum space.
What would settle it
Fit the LBB model to the combined TOTEM low-|t| and high-|t| elastic pp data at $\sqrt{s} = 8$ TeV: if the best-fit Lév y index $\alpha_L$ is statistically compatible with 2, or if the fit does not reach CL $\geq 0.1\%$, or if any extracted $A_{ab}$ makes the constituent-constituent inelastic probability $\sigma_{ab}^{in}(0)$ exceed 1, the paper's central expectation is falsified.
Extended reading notes
Core claim
The central claim is that the Lévy α-stable generalization of the Bialas-Bzdak model is made by changing both (i) the inelastic scattering probabilities of two constituents and (ii) the quark-diquark distribution inside the proton from Gaussian shapes to Lévy α-stable shapes. The constituent-constituent inelastic probability is defined as a convolution of two Lévy α-stable distributions, which is again Lévy α-stable with scale parameter $S_{ab} = (R_a^{\alpha_L} + R_b^{\alpha_L})^{1/\alpha_L}$, and the proton's internal distribution is a Lévy α-stable separation distribution between the quark and the diquark. The model reduces to the ReBB model when $\alpha_L = 2$, and the paper argues that the new free parameter $\alpha_L$ is expected to come out below 2, motivated by the strong non-exponential low-|t| behavior seen by TOTEM and by the good performance of a simpler Lévy α-stable model with $\alpha \approx 1.959$.
Load-bearing premise
The load-bearing premise is that the inelastic constituent-constituent scattering probability defined by the Lév y α-stable convolution remains between 0 and 1 for all parameter values with $\alpha_L < 2$, so that the Glauber expansion and the unitarized amplitude describe a physical probability.
Editorial extensions
If this is right
- The LBB model should describe simultaneously the low-|t| and high-|t| domains of elastic pp and p\bar{p} differential cross sections with a Lévy index $\alpha_L < 2$, resolving the statistical failure of the ReBB model at √s = 8 TeV.
- Because $\alpha_L < 2$ produces heavier tails in the impact-parameter distribution, the model can generate the strong non-exponential low-|t| behavior observed by TOTEM, which the Gaussian ReBB model cannot reproduce.
- After fits reach statistical acceptance, the model can be used to study the discrepancy between ATLAS and TOTEM total cross-section measurements.
- With Coulomb-nuclear interference effects included, the model can extract the Odderon contribution to the parameter $\rho_0$ at √s = 13 TeV.
- If fits yield $\alpha_L$ near the value around 1.959 seen in the simpler model, that would support the Lévy α-stable framework across SPS, Tevatron, and LHC energies.
Reading between the lines
- The same Lév y α-stable replacement of Gaussian distributions could be applied to other hadronic scattering models built on Glauber expansions, potentially inheriting the same unification of low- and high-|t| behavior.
- A physically meaningful $\alpha_L < 2$ may indicate that the proton's parton cloud has power-law tails, connecting the model's parameters to generalized central-limit-theorem arguments rather than being purely phenomenological.
- A necessary check not performed in this paper is that $\sigma_{ab}^{in}(\vec{s})$ stays within [0,1] for all allowed parameters with $\alpha_L < 2$; for Lév y distributions the peak at $\vec{s}=0$ grows with $\Gamma(2/\alpha_L)/\alpha_L$, so the normalization $A_{ab}$ must be small enough to preserve unitarity.
- The model can be tested at other LHC energies (e.g., 2.76 and 13 TeV) to see whether a single $\alpha_L$ describes the energy evolution or whether the Lév y index itself runs with $\sqrt{s}$.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This short proceedings paper introduces a Lévy α-stable generalization of the Bialas–Bzdak (ReBB) model of elastic proton-proton and proton-antiproton scattering. The constituent-constituent inelastic scattering probabilities of Eq. (2) and the quark-diquark distribution of Eq. (3) are replaced by symmetric Lévy α-stable distributions, with the same stability index α_L for all ingredients. The paper notes that α_L = 2 recovers the Gaussian ReBB model, and it argues that, based on earlier low-|t| analyses, α_L < 2 should allow a simultaneous description of low- and high-|t| data. No data fitting is performed; the paper is a model proposal together with motivation from the authors' previous work.
Significance. If the LBB model indeed describes both the low-|t| non-exponential behavior and the high-|t| dip region of elastic pp and p̄p scattering with a single value of α_L < 2, it would be a useful tool for studying the ATLAS–TOTEM tension and the Odderon contribution. The mathematical generalization is cleanly specified: the convolution identity in Eq. (2) is valid for symmetric stable distributions, the normalization of Eq. (3) is correct, and the α_L = 2 limit reproduces the Gaussian ReBB forms. The paper is honest about being a proposal and explicitly says that applying the model to data is the next step. Its significance is prospective rather than demonstrated; no new experimental insight is obtained yet. The main weakness is the absence of any check that the constituent-level probabilities remain physical for α_L < 2.
major comments (2)
- [Section 3, Eq. (2)] The paper does not enforce the probabilistic bound 0 ≤ σ_ab^in ≤ 1. For the symmetric stable density normalized as L(0|α,R) = Γ(2/α)/(2π α R^2), the zero-impact value is σ_ab^in(0) = A_ab f(α_L) with f(α) = 2Γ(2/α)/α. Since f(2)=1 and f(α)>1 for every α<2, using the ReBB value A_qq=1 (as in the fits displayed in Figs. 1 and 2) gives σ_qq^in > 1 at b=0 for the very regime α_L<2 that motivates the paper. Because the stable density is continuous and maximal at 0, the violation occurs on a set of positive measure, so (1−σ_ab^in) in Eq. (1) can become negative and the square root in Eq. (5) imaginary. The manuscript should state and impose the constraint A_ab ≤ α_L/[2Γ(2/α_L)] and discuss whether the expected α_L ≈ 1.959 is compatible with the ReBB calibration A_qq=1.
- [Section 4 (Summary)] The Summary asserts that the LBB model 'is expected to describe simultaneously the low-|t| and high-|t| domains of elastic pp and p̄p dσ/dt' with α_L < 2. This is presented as a central outcome of the paper, but no data comparison, no fit, and no model calculation of a differential cross section is performed anywhere in the manuscript. The expectation is a hypothesis motivated by the ReBB results and by Ref. [8], not a demonstrated property of the LBB model. Please rephrase the claim as a program to be carried out in future work, and specify the intended kinematic domains and the fitting procedure already at the proposal stage.
minor comments (5)
- [Section 2, second paragraph] The text says '√s is the squared center of mass energy'; this should read '√s is the center-of-mass energy'.
- [Section 3, Eq. (2)] The Fourier convention for L(⃗x|α_L,R_L) should be stated explicitly, e.g., with the measure d²q/(2π)², so that the normalization check leading to L(0|α,R)=Γ(2/α)/(2π α R^2) is transparent.
- [Figures 1 and 2] The captions do not identify which symbols correspond to the TOTEM and ATLAS data and which curves are the ReBB fit; please add explicit legends or state this in the captions.
- [Section 3, final paragraph] The sentence 'The power of a simple Lévy α-stable model for elastic scattering was demonstrated in Ref. [8]' should explicitly repeat that the demonstration concerns low-|t| data only, as done earlier in the same paragraph, to avoid giving the impression that high-|t| data were already described by that model.
- [References] In the Introduction, 'studies published in 2021 and 2022 [3, 4]' is slightly misleading because Ref. [4] is dated 2021; please correct the year description or the citation grouping.
Circularity Check
No significant circularity: the LBB construction is a self-contained model generalization, with no fit-derived predictions.
full rationale
The paper does not perform any fits and makes no empirical prediction that could reduce to an input. Its central content is a model construction: Eqs. (2) and (3) define the Lévy α-stable generalization directly in terms of the stable distribution L, and the α_L = 2 limit is a mathematical identity because the Gaussian is the α = 2 stable distribution. The cited references [7,8] are authored by the present authors, but they are used for standard convolution properties and for prior empirical fits to SPS, Tevatron, and LHC data, not as a substitute for the construction. The construction itself is displayed in the paper and does not depend on the cited papers for its validity. The noted unitarity-bound issue for α_L < 2 is a physical correctness concern, not a circularity, because it does not involve the derivation re-entering its own assumptions. Therefore no circular step can be exhibited, and the appropriate score is 0.
Assumptions & free parameters
free parameters (4)
- alpha_L
- R_q, R_d, R_qd (scale parameters)
- A_ab (normalization of constituent inelastic probabilities)
- lambda = m_q/m_d
assumptions (3)
- standard math Convolution of symmetric Levy alpha-stable distributions is a stable distribution with the same index alpha_L.
- domain assumption The proton is a bound state of a constituent quark and a constituent diquark, the p=(q,d) picture.
- domain assumption The inelastic constituent probability sigma_ab^in remains in [0,1] for all parameters and alpha_L, ensuring unitarity of the Glauber expansion and the amplitude in Eq. (5).
Cite this review
Pith. "Pith review of L\'evy $\alpha$-stable generalization of the ReBB model of elastic proton-proton and proton-antiproton scattering." pith.science (2026). https://pith.science/paper/WAWL2MKU
@misc{pith2026250201911,
author = {Pith},
title = {Pith review of: L\'evy $\alpha$-stable generalization of the ReBB model of elastic proton-proton and proton-antiproton scattering},
year = {2026},
howpublished = {\url{https://pith.science/paper/WAWL2MKU}},
note = {Machine review of arXiv:2502.01911}
}
abstract
The L\'evy $\alpha$-stable generalization of the ReBB model of elastic proton-proton and proton-antiproton scattering is presented. The motivation for the future use of this model in describing experimental data is discussed.
Figures
Reference graph
Works this paper leans on
- [8]
- [1]
-
[2]
R. J. Glauber and G. Matthiae. Nucl. Phys. B, 21:135–157, 1970
work page 1970
- [3]
- [4]
- [5]
- [6]
- [7]
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.