REVIEW
The effect of Coulomb interaction on Levy source parameters
High Energy Physics - Phenomenology · Submitted 29 Oct 2026 · 28 references
Pith's verdict: An honest, clearly-scoped numerical scan of Coulomb-induced shifts in the fitted Levy index that delivers qualitative trends, but the central attribution to Coulomb is not isolated because there is no no-Coulomb baseline — a gap the authors acknowledge.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
A simulation of charged pions shows the fitted Lévy stability index drops, making the source tail heavier.
The reading
What carries the argument
The load-bearing machinery is the pair-source distribution $D(r)=\int d^3\sigma\, S(\sigma+r/2)S(\sigma-r/2)$, computed from the final phase-space coordinates and fitted by the spherically symmetric Lévy distribution $L(r;\alpha,R)$, whose index $\alpha$ controls the tail. Initial positions are sampled with the Kanter method so that the Cartesian components follow a symmetric $\alpha_0$-stable law. Particles then evolve under the Yukawa force with a time-dependent Debye length $\lambda_D(t)=\sqrt{T(t)/(4\pi\alpha_{\rm em} n(t))}$, propagated by leapfrog integration with an adaptive time step; the two observables, $\alpha_{\rm final}/\alpha_0$ and the fitted slope $m$, convert the shape change into a single number.
Extended reading notes
Core claim
On its own terms, the paper establishes a directional claim: after Coulomb evolution, fitting the pair-source distribution $D(r)$ with a one-dimensional Lévy distribution returns $\alpha_{\rm final} < \alpha_0$ in the studied setups. The suppression is quantified by the ratio $\alpha_{\rm final}/\alpha_0$ and by the slope $m$ of the linear fit to the $\alpha_{\rm final}$ versus $\alpha_0$ scatter. The paper reports that both observables decrease as the particle number grows, decrease more strongly as the initial source size shrinks, and drop substantially as $\alpha_0$ moves from 1 toward 1.5, while remaining stable against the upper limit of the fit range. The authors interpret this as the Coulomb interaction making the source tail heavier, and they caution that late-time positions contain free streaming and that only the power-law tail is fitted, so the result is presented as qualitative.
Load-bearing premise
The load-bearing premise is that the observed decrease in the fitted $\alpha$ is caused by the Coulomb interaction; this is not isolated, because no Coulomb-free simulation is reported and the authors acknowledge that free streaming plus a tail-only one-dimensional fit could influence the result.
Editorial extensions
If this is right
- Femtoscopic models that omit or approximate Coulomb will systematically obtain smaller $\alpha$ values than the true source, so part of the model-data discrepancy may reflect missing Coulomb evolution.
- The correction is largest for small, dense sources and high particle numbers, meaning the systematic bias depends on centrality and system size rather than being a constant offset.
- Since the distortion grows as $\alpha_0$ moves from 1 to 1.5, Gaussian-like sources receive a larger Coulomb correction than Cauchy-like sources.
- The near-flat dependence on $r_{\rm max}$ means the conclusion is not an artifact of the specific upper fit cutoff, although the tail-only one-dimensional fit remains a caveat.
- The paper's own proposed refinement—back-extrapolated emission points, a short-distance regulator, and a three-dimensional Lévy fit—would turn the qualitative direction into a quantitative correction table.
Reading between the lines
- Editorial inference: a Coulomb-off control run would show how much of the decrease is interaction physics versus free-streaming plus tail-only fitting; if the control also drops, the magnitude would shrink even if the direction stands.
- Editorial inference: if the effect holds in refined simulations, a shift of roughly 8-10% at small source sizes is of the same order as some discrepancies between measured and modeled $\alpha$, making screened Coulomb a candidate explanation in small, dense systems.
- Editorial extension: the same simulation could report the fitted Lévy scale parameter $R$ alongside $\alpha$, testing whether Coulomb distorts the source width as well as its tail.
Original abstract
Femtoscopy plays an important role in the study of heavy-ion collisions and particle sources generated within them. In both experimental measurements and theoretical models, a L\'evy distribution is used to describe the shape of the source, as the Gaussian distribution has proven to be inaccurate based on experimental observations. Despite this, there are still discrepancies between experimental data and the results of theoretical models. One reason may be that these theoretical models either approximate the Coulomb interaction poorly or omit it entirely. In our research, we focused on how the Coulomb effect influences the stability index $\alpha$ value of the L\'evy distribution. For this, we used a custom-developed simulation framework written in C++ and Python. The C++ module integrates the equations of motion for relativistic pions using the Leapfrog algorithm, while the Python post-processing script generates the two-particle spatial source function from the phase-space coordinates. We investigate the ratio and slope of the generated and pair-source fitted $\alpha$ as a function of the number of particles, fitting range, and the initial source size and density.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
Femtoscopy plays an important role in the study of heavy-ion collisions and particle sources generated within them. In both experimental measurements and theoretical models, a L\'evy distribution is used to describe the shape of the source, as the Gaussian distribution has proven to be inaccurate based on experimental observations. Despite this, there are still discrepancies between experimental data and the results of theoretical models. One reason may be that these theoretical models either approximate the Coulomb interaction poorly or omit it entirely. In our research, we focused on how the Coulomb effect influences the stability index $\alpha$ value of the L\'evy distribution. For this, we used a custom-developed simulation framework written in C++ and Python. The C++ module integrates the equations of motion for relativistic pions using the Leapfrog algorithm, while the Python post-processing script generates the two-particle spatial source function from the phase-space coordinates. We investigate the ratio and slope of the generated and pair-source fitted $\alpha$ as a function of the number of particles, fitting range, and the initial source size and density.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | Femtoscopy plays an important role in the study of heavy-ion collisions and particle sources generated within them. | background 1.00 | 0.04 | 0.02 | asserted 1.00 |
| s1 | In both experimental measurements and theoretical models, a L\'evy distribution is used to describe the shape of the source, as the Gaussian distribution has proven to be inaccurate based on experimental observations. | background 1.00 | 0.06 | 0.03 | asserted 1.00 |
| s2 | Despite this, there are still discrepancies between experimental data and the results of theoretical models. | gap 0.63 | 0.08 | 0.02 | asserted 0.99 |
| s3 | One reason may be that these theoretical models either approximate the Coulomb interaction poorly or omit it entirely. | gap 0.80 | 0.05 | 0.03 | hedged possibility 1.00 |
| s4 | In our research, we focused on how the Coulomb effect influences the stability index $\alpha$ value of the L\'evy distribution. | contribution 1.00 | 0.88 | 0.08 | asserted 0.77 |
| s5 | For this, we used a custom-developed simulation framework written in C++ and Python. | method 1.00 | 0.89 | 0.02 | asserted 0.99 |
| s6 | The C++ module integrates the equations of motion for relativistic pions using the Leapfrog algorithm, while the Python post-processing script generates the two-particle spatial source function from the phase-space coordinates. | method 1.00 | 0.87 | 0.02 | asserted 1.00 |
| s7 | We investigate the ratio and slope of the generated and pair-source fitted $\alpha$ as a function of the number of particles, fitting range, and the initial source size and density. | method 0.55 | 0.94 | 0.07 | asserted 0.84 |
What was examined, and what was not
- arXiv listing abstract v1
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
Figures
Figures reproduced from arXiv:2610.39626 by the authors.
References
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.