Mid-rapidity dependence of hadron production in $p-p$ and $A-A$ collisions

A.I. Malakhov; G.I. Lykasov

REVIEW 3 major objections 4 minor 22 references

Mid-rapidity dependence of hadron production in $p-p$ and $A-A$ collisions

T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read The paper derives an analytic self-similarity parameter for small nonzero rapidity that reproduces pion rapidity spectra in pp and AA collisions to better than 10% at |y|≤0.3 with universal energy scaling.

desk verdict Extends a self-similarity parameter to nonzero rapidity, but the printed central formula has an inverted δ that makes the integrand imaginary at the lower limit, so the data fits are unsupported as written. read the letter →

arxiv 1908.04639 v1 pith:X3TGXU6V submitted 2019-08-13 hep-ph

A.I. Malakhov , G.I. Lykasov This is my paper · ORCID

classification hep-ph

keywords self-similarityrapiditydistributionpionproductioninclusivespectraheavy-ioncollisionsuniversalenergyscalingtransversemassPomeronintercept

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper tries to show that the self-similarity description of inclusive hadron production, which was previously solved only at midrapidity $y=0$, can be extended into the small-nonzero-rapidity region without introducing new free parameters. The key result is an analytic expression for the self-similarity parameter $\Pi(s,m_{1t},y)$ at $|y|\ll 1$, obtained by symmetrizing the two scalar products between initial and produced four-velocities. Plugging this $\Pi$ into the already fitted universal function $F(\Pi)$ reproduces measured pion rapidity spectra in $pp$ and $AA$ collisions at $|y|\le 0.3$ to better than 10%, across beam energies from a few GeV to 200 GeV per nucleon. The reason this would matter is that it turns the $y$-dependence of hadron production at mid-rapidity into a universal, energy-scaling statement, $(s/s_0)^{\Delta}$, rather than a collection of model fits.

What carries the argument

The central object is the self-similarity parameter $\Pi$, defined as the minimum, over the fractions $N_A,N_B$ of four-momenta transmitted by the colliding nuclei, of half the norm of the weighted sum of the initial four-velocities $u_A,u_B$ (Eq. 3). The mechanism that carries the argument is the small-$y$ symmetrisation: replacing $(u_A\cdot u_1)$ and $(u_B\cdot u_1)$ by their common form $(m_{1t}/m_1)\cosh(y)\cosh(Y)$ turns $\Phi_A$ and $\Phi_B$ into a single $\Phi$, so the known $y=0$ algebraic solution for $N$ can be reused with only $\cosh(y)$ modifications. This yields Eq. (14), a parameter-free analytic $\Pi(s,m_{1t},y)$. All subsequent results follow from inserting this $\Pi$ into the previously fitted universal spectrum $F(\Pi)$, whose gluon and quark terms carry the constants and the energy scaling $(s/s_0)^{\Delta}$.

What would settle it

Compute the exact minimisation of Eq. (3) using the full scalar products $(u_A\cdot u_1)=(m_{1t}/m_1)\cosh(Y+y)$ and $(u_B\cdot u_1)=(m_{1t}/m_1)\cosh(Y-y)$, and compare the resulting $\Pi$ with Eq. (14) at $y=0.3$ over the transverse-mass range that dominates the integral in Eq. (16). A difference exceeding about 10% would show that the symmetrised analytic form, not the exact kinematics, is responsible for the reported data agreement.

Watch

Extended reading notes

Core claim

The paper's central claim is that at small but nonzero rapidity the self-similarity parameter takes the closed form $\Pi(s,m_{1t},y)\simeq \frac{m_{1t}\cosh(y)}{2m_0\delta}\left[1+\sqrt{1+\frac{M^2-m_1^2}{m_{1t}^2\cosh^2(y)\delta}}\right]$ with $\delta=1-4m_0^2/s$, and that this single expression, inserted into the universal function $F(\Pi)$ of Eq. (5), accounts for the pion rapidity spectra in $pp$ and $AA$ collisions at $|y|\le 0.3$ with a claimed precision better than 10%. The derivation reuses the zero-rapidity solution for the minimising fraction $N$, because at $y\ll 1$ the two scalar products $(u_A\cdot u_1)$ and $(u_B\cdot u_1)$ are both replaced by $(m_{1t}/m_1)\cosh(y)\cosh(Y)$, making the invariant functions $\Phi_A$ and $\Phi_B$ equal. The paper further claims that the integrated rapidity distribution inherits a universal energy dependence of the form $(s/s_0)^{\Delta}$ with $\Delta\simeq 0.12$, the excess of the subcritical Pomeron intercept over unity. A simplified low-$y$ closed form for $d\sigma/dy$ is also given, Eq. (17), and is compared with the thermal-model form containing longitudinal and transverse flow.

Load-bearing premise

The derivation's load-bearing premise is that for $|y|\le 0.3$ the two scalar products $(u_A\cdot u_1)$ and $(u_B\cdot u_1)$ may be replaced by their symmetric average $(m_{1t}/m_1)\cosh(y)\cosh(Y)$; if the omitted $\sinh(y)$ terms shift $\Pi$ by more than the claimed 10%, the data comparison and the universal energy scaling lose their support.

Editorial extensions

If this is right

At $|y|\le 0.3$, pion rapidity spectra can be computed with no additional parameters; the $y$-dependence enters only through $\cosh(y)$ in $\Pi(s,m_{1t},y)$.
The integrated rapidity distribution inherits a universal energy scaling $(s/s_0)^{\Delta}$ with $\Delta\simeq 0.12$, so spectra at different collision energies should collapse onto one curve.
The simplified closed form Eq. (17) gives an analytic approximation for $d\sigma/dy$ at small $y$, structurally similar to thermal-model results with longitudinal and transverse flow.
The same $\Pi$ can describe projectile mesons and any hadronic species, because the scalar-product identity Eq. (13) does not depend on the projectile mass.
Deviations from this description are expected to appear mainly at $|y|>0.3$, where the asymmetric terms and nuclear thermal effects enter.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

Because the symmetrisation error grows like $\sinh(y)$, one can compute the exact $\Pi$ numerically to map where the 10% band fails; this would replace the empirical $|y|\le 0.3$ window with a quantitative boundary.
Eq. (17) is structurally identical to a thermal-model rapidity distribution with longitudinal and transverse flow, so the fitted constant $C_q$ can be converted into an effective inverse slope in rapidity; comparing it with independently extracted flow velocities would separate geometric from collective-flow contributions to the mid-rapidity plateau.
The derivation is species-blind apart from $m_1$ and $M$, so the same $\Pi$ should describe kaon and antiproton spectra at $|y|\le 0.3$; those data would test whether the universal function $F(\Pi)$ is truly universal or tuned to pions.
If the $(s/s_0)^{\Delta}$ scaling holds across energies, it provides a parameter-free baseline that other mid-rapidity signals, such as broadening or enhancement, would have to be measured against.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

Colleague,

Quick take: this is a modest extension of the self-similarity parameter to small nonzero rapidity, and the idea is fine, but the printed central formula is wrong. Eq. (14) has δ = 1 - 4m0²/s in the wrong place in the square root, so it becomes imaginary at the lower integration limit for any finite √s. The data fits in Fig. 1 cannot be reproduced from the equations as written.

I checked against Section 3. At y=0, M=0, the paper's own definitions give Φ = m1t chY/(2m0 sh²Y) and ΦM/Φ² = -δ m1²/m1t², so 1 + ΦM/Φ² = 1 - δ m1²/m1t². Eq. (14) instead has sqrt(1 - m1²/(m1t² δ)). That's an inverted δ. The likely source is the N formula just above: N = 1 + ... Φ doesn't satisfy Eq. (6); the correct root is N = Φ(1 + sqrt(1 + ΦM/Φ²)). Fix that and the Π expression follows, but Eq. (14) still has the δ error.

Now the fair parts. The step from y=0 to finite y is genuinely new, and it is done without extra fitted parameters. The authors are transparent about which constants come from earlier papers. That's good practice.

But even after a typo fix, the paper oversells itself. The claim of better than 10% agreement is not backed by error bars, residuals, or any quantitative measure. The small-y approximation drops terms of order th(y) ≈ 0.29 at |y|=0.3, so the claimed validity range is at the edge of the expansion. The conclusion's universal energy scaling is approximate; Eq. (17) is a simplified limit, not a proof.

No code or data is provided, so the fits are not independently checkable.

Recommendation: this deserves a referee, not a desk reject, because the error is concrete and fixable, and the underlying extension is plausible. But the revision must correct Eq. (14), redo the fits, and add actual agreement measures. As printed, I would not cite it.

Referee Report

3 major / 4 minor

Summary. The manuscript extends the self-similarity description of inclusive hadron production from y=0 to small nonzero rapidity. From four-momentum conservation (Eq. 2) and the minimization condition, the authors derive an analytic expression for the self-similarity parameter Π(s, m_1t, y), Eq. (14), claimed to be valid for |y| ≲ 0.3. Inserting this Π into the previously constructed universal function F(Π) of Eq. (5) with constants taken from earlier papers, they compute pion rapidity spectra via Eq. (16) for pp, AuAu, and PbPb collisions and compare with AGS, RHIC, SPS, and NA61/SHINE data in Fig. 1. They claim a satisfactory description with precision better than 10% and a universal energy dependence of the form (s/s0)^Δ.

Significance. If the corrected version of Eq. (14) holds, the paper would provide a parameter-free analytic dependence of the self-similarity parameter on rapidity, obtained from conservation laws rather than from fitting the y-dependent data. That would be a useful extension of the previous y=0 formalism. The constants in F(Π) are fixed from earlier analyses, so the y-dependence is in principle predictive. However, the central formula as printed contains an algebraic error that makes it inconsistent with the derivation and yields imaginary values at the lower integration boundary, and the claimed 10% precision is not quantitatively demonstrated. The significance of the paper therefore depends on a substantial revision and re-analysis, not merely on presentation changes.

major comments (3)

[Section 4, Eq. (14)] At y=0 and M=0 (pions), Eq. (14) reduces to Π = [m_1t/(2m_0 δ)] [1 + sqrt(1 - m_1²/(m_1t² δ))], with δ = 1 - 4m_0²/s. However, substituting Eqs. (10)-(12) into Π = N chY gives Π = [m_1t/(2m_0 δ)] [1 + sqrt(1 - δ m_1²/m_1t²)]. The factor δ appears inverted inside the square root in the printed equation. Since δ < 1 for any finite s, the printed expression becomes imaginary at m_1t = m_π, e.g., at √s = 5 GeV, so the integrand in Eq. (16) is not a real function and Fig. 1 cannot have been produced from Eq. (14) as written. This is a load-bearing algebraic error that must be corrected before the data comparison can be assessed.
[Section 3, Eqs. (9)-(12)] The derivation replaces (u_A·u_1) and (u_B·u_1) with the symmetric expression (m_1t/m_1) ch(y) ch(Y), dropping terms proportional to sh(Y) sh(y). At the claimed boundary |y| = 0.3, sh(y)/ch(y) ≈ 0.29, and for the AGS energies chY is of order 2, so thY ≈ 0.9; the dropped terms are not numerically small. The manuscript provides no estimate of the resulting error in Π, so the stated validity range |y| ≤ 0.3 is not justified by the derivation as presented.
[Section 4, after Eq. (16) and Fig. 1] The claim of 'a precision less than 10%' is not supported by any residual plot, uncertainty band, chi-square value, or numerical comparison table. Fig. 1 displays curves and data points but no quantitative agreement measure. Similarly, the universal energy dependence (s/s0)^Δ rests on the unquantified assertion that the first term of Eq. (5) dominates at low y; without quantitative evidence, the central claims of the paper are not verifiable from the material presented.

minor comments (4)

[Abstract] The phrase 'at not large rapiditiesy' contains a typographical error and should read 'at not large rapidities y'; the wording 'is illustrated within a good agreement' is awkward and should be rephrased.
[Section 2, Eq. (4)] The notation Aα(NA_A · Aα(NB_B is typeset incorrectly; it should presumably read A^{α(N_A)} A^{α(N_B)} F(Π), with the exponents and parentheses clearly displayed.
[Fig. 1 caption] The caption contains 'shirt dash line', which should be 'short dash line'.
[Section 4, Eqs. (15)-(16)] Eq. (16) defines dσ/dy as a differential cross section, while Fig. 1 is labeled dN/dy; the conversion between cross section and yield, and the overall normalization used for the curves, are not specified and should be stated for reproducibility.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the new y-dependent form of Pi is derived from four-momentum conservation and compared with external data, not fitted to it.

full rationale

The new element in the paper is Eq. (14), an analytical expression for Pi(s,m_1t,y). The derivation in Sec. 3 starts from the conservation law Eq. (2), introduces the scalar products (u_A.u_1) and (u_B.u_1), and at small y replaces them by (m_1t/m_1) cosh(y) cosh(Y), Eq. (9). Equations (10)-(12) are then algebraic developments written out in the paper rather than imported as a black box. The y-dependence therefore does not reduce to a fit parameter: no constant is adjusted to the |y|<=0.3 data, and the comparisons with AGS, RHIC/SPS, and NA61/SHINE data are external checks. Self-citations to [7-9] supply the starting model and the previously fixed constants in F(Pi), Eq. (5), but those are inputs whose values are not re-derived from the y data; hence the central y-dependent prediction is not equivalent to its inputs by construction. The energy factor (s/s0)^Delta is contained in F(Pi) and is explicitly presented as a known input; the paper says it is used and confirmed, not newly derived. No circular step satisfying the hard-evidence requirement was found. A separate observation that the printed Eq. (14) may be algebraically inconsistent with the y=0 limit of Eqs. (10)-(12) would be a correctness or consistency concern, not a circularity. Because the result is not forced by a self-citation chain or by identifying the input with the output, the circularity score is 0.

Assumptions & free parameters 7 free parameters · 5 assumptions · 0 invented entities

The central comparison rests on the self-similarity conservation law and minimization principle from earlier papers, on the universal function F(Π) with constants fitted by overlapping authors, and on a small-rapidity approximation that is not error-controlled. No new physical entities are introduced. The genuinely new content is the analytical y-dependence of Π.

free parameters (7)

Aq = 3.68 (GeV/c)^-2
This is the coefficient of the quark term in F(Π), Eq. (5); it was fitted to inclusive spectra in ref. [21] by overlapping authors and used unchanged here.
Cq = 0.147
This is the exponent in the quark exponential term of F(Π); it has the same fitted origin as Aq.
Ag = 1.7249 (GeV/c)^-2
This is the coefficient of the gluon term in F(Π), taken from refs. [21,22] by overlapping authors.
Cg = 0.289
This is the exponent in the gluon exponential term of F(Π), taken from refs. [21,22].
g = 21 mb
This is the normalization constant in Eq. (5), computed in the quasi-eikonal approximation of ref. [20]; it is not fitted in this paper.
Delta = 0.12
This is the excess of the subcritical Pomeron intercept over 1; it is a phenomenological input to Eq. (5) determined from pp data.
s0 = not specified in text
This is the reference energy scale in the factor (s/s0)^Delta; it is inherited from ref. [9] and not defined in this paper.

assumptions (5)

domain assumption Four-momentum conservation takes the form (N_A P_A + N_B P_B - p_1)^2 = (N_A m0 + N_B m0 + M)^2, with fractions N_A, N_B and a bookkeeping mass M for quantum number conservation (Eq. 2).
This is the foundational ansatz of the self-similarity approach; the paper uses it without derivation.
domain assumption The self-similarity parameter is Π = min over N_A, N_B of (1/2)[(u_A N_A + u_B N_B)^2]^{1/2}, following the minimization principle of ref. [7] (Eq. 3).
This defines the variable whose universality the paper asserts.
domain assumption The inclusive spectrum is a universal function F(Π) of the form in Eq. (5), combining quark and gluon exponentials and a Pomeron factor.
This functional form and its constants are imported from refs. [9,21,22].
domain assumption For y<<1, the scalar products obey (u_A·u_1) ≈ (u_B·u_1) ≈ (m_1t/m_1) ch(y) ch(Y), so terms with sh(y) are dropped (Eqs. 9-10).
This small-rapidity approximation is what allows Φ_A=Φ_B and the reuse of the y=0 solution; its error grows with y and is not quantified in the paper.
domain assumption The Mandelstam-variable representation of Π from ref. [9] remains valid at small nonzero y and is used to write Eq. (14).
The paper relies on this prior result to express Π in terms of s and m_1t.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Mid-rapidity dependence of hadron production in $p-p$ and $A-A$ collisions." pith.science (2026). https://pith.science/paper/X3TGXU6V

@misc{pith2026190804639,
  author       = {Pith},
  title        = {Pith review of: Mid-rapidity dependence of hadron production in $p-p$ and $A-A$ collisions},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/X3TGXU6V}},
  note         = {Machine review of arXiv:1908.04639}
}

abstract

The calculation of inclusive spectra of pions produced in $pp$ and $AA$ collisions as a function of rapidity $y$ is presented within the self-similarity approach. It is shown that at not large rapidities $y$ one can obtain the analytical form of the self-similarity function $\Pi(y,p_t)$ dependent of $y$ and hadron transverse momentum $p_t$. A satisfactory description of data on the rapidity spectra at $|y|\leq$ 0.3 is illustrated within a good agreement. The universal energy dependence of these spectra is also shown.

Figures

Figures reproduced from arXiv: 1908.04639 by the authors.

**Figure 1.** Left: pion rapidity y-spectra in AuAu collision at √ s = 4.31 A GeV (solid line), 3.84 A GeV (long dash line), 3.32 A GeV (shirt dash line), 2.7 A GeV (dashed-dotted line) or the initial kinetic energies per nucleon about Ekin = 8, 6, 4, 2 GeV, respectively. They are compared to the AGS data [17]. Middle: pion y-spectra in AuAu collision (solid line, RHIC data) at √ s = 200 A GeV and PbPb collision (SPS data) at √ s… view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

22 extracted references · 22 canonical work pages

[1]

Fermi, Phys

E. Fermi, Phys. Rev. 92, 452 (1953)

work page 1953
[2]

I. Ya. Pomeranchuk, Izv. Dokl. Akad. Nauk Ser.Fiz. 78, 889 (1951)

work page 1951
[3]

Landau, Izv

L.D. Landau, Izv. Akad. Nauk Ser. Fiz. 17, 51 (1953)

work page 1953
[4]

Hagedorn, Supplemento al Nuovo Cimento 3, 147 (1965)

R. Hagedorn, Supplemento al Nuovo Cimento 3, 147 (1965)

work page 1965
[5]

Baldin, L.A

A.M. Baldin, L.A. Didenko, Fortsch.Phys. 38, 261 (1990)

work page 1990
[6]

Baldin, A.I

A.M. Baldin, A.I. Malakhov, and A. N. Sissakian, Phys. Part. Nucl.29 (Suppl. 1), 4 (2001)

work page 2001
[7]

M.Baldin, A

A. M.Baldin, A. A. Baldin. Phys. Particles and Nuclei, 29 No3, 232 (1998)

work page 1998
[8]

Baldin, A

A.M. Baldin, A.. Malakhov. JINR Rapid Communications, No.1(87)-98, pp.5-12 (1998)

work page 1998

Show all 22 references

[9]

Lykasov, A.I

G.I. Lykasov, A.I. Malakhov, Eur. Phys. J. A 54, 187 (2018)

2018
[10]

Lykasov, A

G. Lykasov, A. Malakhov, Eur. Phys. J. (Web of Conf.) 204, 01022 (2019)

2019
[11]

Malakhov, G

A. Malakhov, G. Lykasov, Eur. Phys. J. (Web of Conf.) 204, 01021 (2019)

2019
[12]

Schnedermann, J

E. Schnedermann, J. Sollfrank, U. Heinz, Phys.Rev.C48,2462 (1993)

1993
[13]

G. Wilk, Z. Wlodarczyk, Phys.Lett. 84, 2770 (2000)

2000
[14]

Bugaev, J.Phys.G:Nucl.Phys., 28, 1981 (2002)

K.A. Bugaev, J.Phys.G:Nucl.Phys., 28, 1981 (2002)

2002
[15]

Bugaev, M

K.A. Bugaev, M. Gadzicki, M.I. Gorenstein, Phys.Lett. B 544, 127 (2002)

2002
[16]

Cleymans, G.I

J. Cleymans, G.I. Lykasov, A.S. Parvan, et al., Phys.Lett. B 723, 351 (2013)

2013
[17]

Kley, et al., E895 Collaboration, Phys.Rev

J.L. Kley, et al., E895 Collaboration, Phys.Rev. C 68, 054905 (2003)

2003
[18]

Abgrall, et al., NA61 /SHINE Collaboration, Eur

N. Abgrall, et al., NA61 /SHINE Collaboration, Eur. Phys. J. C 74, 2794 (2014)

2014
[19]

Cleymans, J

J. Cleymans, J. Struempfer, L. Tirko, Phys.Rev. C 78, 017901 (2008)

2008
[20]

Ter-Martirosyan, Sov.J.Nucl.Phys., 44, (1986) 817

K.A. Ter-Martirosyan, Sov.J.Nucl.Phys., 44, (1986) 817

1986
[21]

Grinyuk, G.I

A.A. Grinyuk, G.I. Lykasov, A.V . Lipatov, N.P. Zotov, Phys.Rev.D87, (2013) 074017

2013
[22]

Artemenkov, G.I

D.A. Artemenkov, G.I. Lykasov, A.I. Malakhov, Int.J.Mod.Phys. A30 (2015) 1550127. 6

2015

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [1]

Fermi, Phys

E. Fermi, Phys. Rev. 92, 452 (1953)

work page 1953

[2] [2]

I. Ya. Pomeranchuk, Izv. Dokl. Akad. Nauk Ser.Fiz. 78, 889 (1951)

work page 1951

[3] [3]

Landau, Izv

L.D. Landau, Izv. Akad. Nauk Ser. Fiz. 17, 51 (1953)

work page 1953

[4] [4]

Hagedorn, Supplemento al Nuovo Cimento 3, 147 (1965)

R. Hagedorn, Supplemento al Nuovo Cimento 3, 147 (1965)

work page 1965

[5] [5]

Baldin, L.A

A.M. Baldin, L.A. Didenko, Fortsch.Phys. 38, 261 (1990)

work page 1990

[6] [6]

Baldin, A.I

A.M. Baldin, A.I. Malakhov, and A. N. Sissakian, Phys. Part. Nucl.29 (Suppl. 1), 4 (2001)

work page 2001

[7] [7]

M.Baldin, A

A. M.Baldin, A. A. Baldin. Phys. Particles and Nuclei, 29 No3, 232 (1998)

work page 1998

[8] [8]

Baldin, A

A.M. Baldin, A.. Malakhov. JINR Rapid Communications, No.1(87)-98, pp.5-12 (1998)

work page 1998