{"id":"37057f28-d0c6-437a-b2ae-5a1a78b2ac86","arxiv_id":"1908.04639","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":7,"one_line_summary":"An analytical self-similarity parameter Π(s, m_t, y) is derived for small rapidity and used to reproduce pion rapidity spectra in pp and AA collisions with claimed agreement within 10%.","lead":"This paper derives a simple formula for how the number of pions produced in proton-proton and nucleus-nucleus collisions varies with rapidity near the central region, using a self-similarity variable built from four-momentum conservation. It then compares the formula with published data at |y| ≤ 0.3 across a wide range of collision energies and reports agreement within 10%.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq. (14) is internally inconsistent: its y=0 limit contradicts Eqs. (10)-(12) and is imaginary at m1t=mπ for any finite sqrt(s), so the central formula and Fig. 1 fits are unsupported as printed.","rationale":"The Pith reader identified the small-y approximation as the weakest assumption. That concern is legitimate, but I found a more basic internal inconsistency in the paper's central new formula. Eq. (14) does not reduce to the y=0 limit of the derivation in Section 3, even when y is exactly zero, because the factor δ enters inside the square root in the wrong place. For pions (M=0), the printed expression is imaginary at the lower limit of transverse-mass integration for every finite collision energy, since δ=1-4m0²/s<1. This is not a matter of external consensus or of the size of dropped sh(y) terms; it is a direct algebraic contradiction within the manuscript. If Eq. (14) were actually used to produce the curves in Fig. 1, the low-energy AGS and SPS calculations would be ill-defined; if it was not used, then the paper's stated central formula is not the one that was tested. In either case, the central claim as written is unsupported. The paper should be rejected in its current form, but a corrected version that places δ in the numerator inside the square root, regenerates all figures, and quantifies the small-y approximation error could be considered for conditional acceptance. I therefore disagree with the reader's weakest_assumption because the algebraic error is the more load-bearing concern.","tokens_in":5145,"tokens_out":10244,"duration_ms":99687,"concrete_test":"Analytic check: set M=0, y=0 in Eq. (14) and compare with the exact y=0 solution obtained from Eqs. (10)-(12). They differ by δ↔1/δ inside the square root. Then evaluate Eq. (16) at √s=4.31 GeV with the printed Eq. (14); the integrand at m1t=mπ is sqrt(1-1/0.810) <0, so the integral is undefined. Recomputing Fig. 1 with the corrected expression sqrt(1-δ m1²/(m1t² ch²y)) would show whether the claimed <10% agreement survives or disappears.","verdict_should_be":"REJECT","load_bearing_attack":"At y=0, M=0, the authors' Section 3 gives Φ=(m1t chY)/(2m0 sh²Y), ΦM=-m1²/(4m0² sh²Y), so ΦM/Φ² = -m1² sh²Y/(m1t² ch²Y) = -m1² δ/m1t², where δ=sh²Y/ch²Y=1-4m0²/s. Hence Π = Φ chY [1+sqrt(1+ΦM/Φ²)] = [m1t/(2m0δ)] [1+sqrt(1-δ m1²/m1t²)]. Eq. (14) at y=0, M=0 instead gives [m1t/(2m0δ)] [1+sqrt(1-m1²/(m1t²δ))]: the δ inside the square root is inverted. This is not a small-y approximation effect; it is present at y=0. Since δ<1 for any finite s, the printed square root becomes negative at the lower integration boundary m1t=mπ, e.g. √s=5 GeV yields 1-1/0.86<0. Thus the integrand in Eq. (16) is undefined as a real function and Fig. 1 cannot be the output of Eq. (14) as written. The expression derivable from Eqs. (10)-(12) has δ in the numerator inside the root, not the denominator.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek-v4-flash","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)^Δ.","tokens_in":5443,"tokens_out":5302,"duration_ms":54005,"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":[{"comment":"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":"Section 4, Eq. (14)"},{"comment":"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":"Section 3, Eqs. (9)-(12)"},{"comment":"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.","section":"Section 4, after Eq. (16) and Fig. 1"}],"minor_comments":[{"comment":"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":"Abstract"},{"comment":"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.","section":"Section 2, Eq. (4)"},{"comment":"The caption contains 'shirt dash line', which should be 'short dash line'.","section":"Fig. 1 caption"},{"comment":"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.","section":"Section 4, Eqs. (15)-(16)"}],"recommendation":"major_revision","confidential_remarks":"The error in Eq. (14) is likely a transposition of δ, but it is central: the entire figure and the claimed 10% precision depend on this formula. I would ask the authors to correct the expression, redo the calculation and figures, and provide a quantitative comparison measure (residuals, chi-square, or tabulated deviations). The paper also leans heavily on constants fitted in earlier papers by overlapping authors; the revision should make the provenance and numerical values of these constants explicit so the y-dependence can be checked independently."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nQuick 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.\n\nI 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.\n\nNow 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.\n\nBut 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.\n\nNo code or data is provided, so the fits are not independently checkable.\n\nRecommendation: 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.","headline":"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.","tokens_in":6013,"tokens_out":14686,"would_cite":false,"duration_ms":125754,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"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.","keywords":["self-similarity","rapidity distribution","pion production","inclusive spectra","heavy-ion collisions","universal energy scaling","transverse mass","Pomeron intercept"],"falsifier":"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.","tokens_in":4903,"feed_emoji":"⚛️","tokens_out":10085,"duration_ms":90091,"temperature":0.7,"pith_summary":"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.","feed_headline":"Pion yields at mid-rapidity follow one universal formula","feed_subtitle":"Extending the self-similarity parameter to |y|≤0.3 reproduces pp and AA pion data with universal (s/s0)^Δ scaling.","key_machinery":"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}$.","core_discovery":"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.","pith_inferences":["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."],"forward_implications":["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."],"supporting_citations":[{"why":"introduces the self-similarity parameter $\\Pi$ as a minimum over four-velocity combinations, the object this paper extends to nonzero rapidity.","marker":"[7]"},{"why":"provides the analytical $y=0$ solution for $\\Pi$ that the small-$y$ approximation of Eq. (11) reuses.","marker":"[8]"},{"why":"establishes the universal function $F(\\Pi)$ and its relation to Mandelstam variables, which Eq. (5) and the energy scaling inherit.","marker":"[9]"},{"why":"gives the thermal-model rapidity distribution to which the simplified form Eq. (17) is compared.","marker":"[12]"},{"why":"supplies the AGS AuAu pion rapidity data at low energies used to validate the $|y|\\le 0.3$ description.","marker":"[17]"},{"why":"supplies the pp pion rapidity spectra at several beam momenta used in the comparison.","marker":"[18]"},{"why":"supplies the RHIC and SPS AA pion rapidity spectra used in the comparison.","marker":"[19]"},{"why":"fixes the constant $g$ in the quasi-eikonal normalization of $F(\\Pi)$.","marker":"[20]"},{"why":"fix the constants $A_q,C_q,A_g,C_g$ in the universal function $F(\\Pi)$.","marker":"[21, 22]"}],"fun_headline_variants":["One formula fits mid-rapidity pion data across systems","Self-similarity predicts pion yields at small rapidity","Mid-rapidity pion spectra unify under single expression","Universal pion law covers pp and AA at |y|≤0.3","Pion yields at mid-rapidity: one self-similar solution"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"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.","fun_headline_variants_meta":{"raw":{"variants":["One formula fits mid-rapidity pion data across systems","Self-similarity predicts pion yields at small rapidity","Mid-rapidity pion spectra unify under single expression","Universal pion law covers pp and AA at |y|≤0.3","Pion yields at mid-rapidity: one self-similar solution"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000303,"raw_usage":{"total_tokens":1761,"prompt_tokens":979,"completion_tokens":782,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":595,"completion_tokens_details":{"reasoning_tokens":694}},"tokens_in":595,"tokens_out":782,"duration_ms":8250,"temperature":1.0,"reasoning_tokens":694,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T13:36:53.092399+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"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.","supporting_citations":[{"cited_title":"M.Baldin, A","cited_arxiv_id":null,"evidence_quote":"introduces the self-similarity parameter $\\Pi$ as a minimum over four-velocity combinations, the object this paper extends to nonzero rapidity."},{"cited_title":"Baldin, A","cited_arxiv_id":null,"evidence_quote":"provides the analytical $y=0$ solution for $\\Pi$ that the small-$y$ approximation of Eq. (11) reuses."},{"cited_title":"Lykasov, A.I","cited_arxiv_id":null,"evidence_quote":"establishes the universal function $F(\\Pi)$ and its relation to Mandelstam variables, which Eq. (5) and the energy scaling inherit."},{"cited_title":"Schnedermann, J","cited_arxiv_id":null,"evidence_quote":"gives the thermal-model rapidity distribution to which the simplified form Eq. (17) is compared."},{"cited_title":"Kley, et al., E895 Collaboration, Phys.Rev","cited_arxiv_id":null,"evidence_quote":"supplies the AGS AuAu pion rapidity data at low energies used to validate the $|y|\\le 0.3$ description."},{"cited_title":"Abgrall, et al., NA61 /SHINE Collaboration, Eur","cited_arxiv_id":null,"evidence_quote":"supplies the pp pion rapidity spectra at several beam momenta used in the comparison."},{"cited_title":"Cleymans, J","cited_arxiv_id":null,"evidence_quote":"supplies the RHIC and SPS AA pion rapidity spectra used in the comparison."},{"cited_title":"Ter-Martirosyan, Sov.J.Nucl.Phys., 44, (1986) 817","cited_arxiv_id":null,"evidence_quote":"fixes the constant $g$ in the quasi-eikonal normalization of $F(\\Pi)$."}],"review_version":1}