{"id":"8c7fdda9-e73b-4698-b21b-93b6a802c590","arxiv_id":"1908.08799","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":4,"one_line_summary":"The Tsallis entropic index q is derived from QCD parameters via fractal self-similarity, giving q = 1 + 1/7, approximately 1.14, consistent with high-energy data.","lead":"The paper argues that scale invariance in Yang-Mills theories creates fractal self-similar structures, and that this leads to Tsallis non-extensive statistics with the index q fixed by QCD parameters. It reports q = 1.14 for QCD, matching values fitted to high-energy collision data.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The value q=1.14 is fixed by equating the QCD beta function with Eq. (49), but Eq. (49) is not derived from the preceding formalism; without it the matching—and hence the prediction—does not follow.","rationale":"The paper promises a first derivation of q from QCD parameters. The derivation has several conjectural steps, but the crucial one is Eq. (49), because it is the only equation that converts the undetermined exponent α=1/(q-1) into the specific number 7. Inspecting the surrounding text, Eq. (48) is the only displayed scale derivative of the effective coupling; for fixed ε_j it is linear in g and lacks the required g^2 and 16π^2, while for ε_j=λ/2 (the assignment used in Fig. 3) the q-exponential factors become constants and the derivative of Eq. (37) vanishes. Either way Eq. (49) does not follow. The reader's weakest assumption, Eq. (30), is also a legitimate concern—it is a self-similarity postulate rather than a derived result—but even if one accepts it, the q-value remains undetermined without Eq. (49). I therefore see the missing derivation of Eq. (49) as the more load-bearing problem. This does not change the REJECT verdict; if anything it sharpens the reason: at the exact point where the quantitative prediction is made, the calculation is absent. If a correct derivation of Eq. (49) were supplied, the claim would become testable and the paper could be reconsidered.","tokens_in":14106,"tokens_out":14365,"duration_ms":121580,"concrete_test":"Compute μ dg/dμ explicitly from Eq. (37) with λ=λ_o/μ and the assignments ε5=ε6=λ/2 stated in Sec. 4 and Fig. 3, and check whether the result can be reduced to -g^3/[16π^2(q-1)]. If not—and the displayed Eq. (48) gives no g^2 or 16π^2—then Eq. (49) is an additional unproven assumption and the q=1.14 matching does not follow.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (49) is the sole point where the numerical value of q is fixed: matching to Eq. (50) gives 1/(q-1)=7 and q=1.14. But Eq. (49) is asserted, not derived. The preceding formula Eq. (48) is μ∂g/∂μ = -Σ_{j=5,6} g μ ε_j [1+(q-1)μ ε_j/λ_o]^{-1}. For fixed ε_j this derivative is linear in g, with no g^2 to produce g^{N~+1}=g^3 and no 16π^2; if instead ε_j=λ/2 as stated for Fig. 3, the q-exponential factors in Eq. (37) become constants and the derivative of the effective coupling vanishes. Neither reading gives β_{\\bar g}=-g^3/[16π^2(q-1)]. The sentence 'The beta function can be easily calculated' (Sec. 4) covers the whole gap, and Eq. (47) does not fill it because ∂Γ/∂g is never computed. The self-similarity condition Eq. (30) fixes the q-exponential form but leaves α=1/(q-1) undetermined; the number 1.14 comes entirely from Eq. (49). The central claim therefore rests on an unproven beta-function identity.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper claims that the scaling properties of Yang-Mills theory produce a self-similar fractal hierarchy of truncated n-point functions, which at high perturbative order can be described by Tsallis non-extensive statistics. The authors derive a q-exponential energy distribution for effective partons, introduce an effective coupling, and compute its beta function. Matching this beta function to the one-loop QCD result gives q=1.14, in agreement with the experimental value 1.14±0.01. The same formalism is then used to obtain a fractal dimension D=0.69 and a multiplicity growth M∼E^0.31, which are compared with intermittency and multiplicity data.","tokens_in":14411,"tokens_out":13510,"duration_ms":116858,"significance":"If the derivation were valid, this would be an important result: it would turn the Tsallis parameter q from a phenomenological fit parameter into a prediction determined by the number of colors and flavors, and it would provide a field-theoretic rationale for power-law tails and intermittency in high-energy collisions. The paper has real strengths: the phase-space density calculation in the appendix leading to Eq. (82) is explicit and correct, and the predicted value q=1.14 is a sharp, falsifiable number. However, the central step connecting the fractal model to the QCD beta function is not derived, and several of the applications rely on algebraic inconsistencies; as it stands the significance is conditional on a major revision of the derivation.","major_comments":[{"comment":"Eq. (49), which fixes the numerical value of q through the matching with Eq. (50), is asserted rather than derived. Eq. (48) has a right-hand side linear in g and contains no 16π² normalization; if ε5=ε6=λ/2 as used in Fig. 3, the q-exponential factors in Eq. (48) are constants, so the g^{N~+1} term in Eq. (49) cannot be obtained from the preceding formulas. The sentence 'The beta function can be easily calculated' is the only transition between these equations. Since this is the sole point where q=1.14 is obtained, the central claim of the paper is unsupported unless a complete derivation of Eq. (49) is supplied.","section":"Sec. 4, Eqs. (48)-(49)"},{"comment":"The self-similarity condition (4N−5)+αν=α is imposed by hand; it is not derived from the renormalization-group equations. This condition is load-bearing because it converts the ideal-gas power law into the q-exponential form and identifies α with 1/(q−1). The parameters N, ν, and α are never connected to QCD parameters, and the final q is fixed only by the later beta-function matching in Sec. 4; thus the claimed first-principles determination of q from field-theory parameters is not established.","section":"Sec. 3, Eq. (30)"},{"comment":"The average energy quoted in Eq. (55) does not follow from the distribution in Eq. (34). For a q-exponential with exponent −1/(q−1), the normalized mean is λ/(3−2q) (for q<3/2), not λ/(2q−1). In addition, the derivation of R=(q−1)/(2q−1) from Eqs. (55)-(57) sets E=λr/(q−1) together with λ/E=r, which are mutually inconsistent except for special values; with the corrected mean, the resulting fractal dimension is not 0.69. The fractal-dimension and multiplicity claims are therefore not supported by the model as presented.","section":"Sec. 6, Eqs. (55)-(62)"},{"comment":"The approximation (1−ε²/M²)^{4N−5}≈1 in Eq. (24) is only valid for ε/M≪1, but the resulting power-law distribution is used for all ε up to M and is integrated in Sec. 6. No error estimate is given, and this approximation is one of the steps that turns a polynomial factor into a q-exponential; it needs a controlled justification.","section":"Sec. 3, Eq. (24)"}],"minor_comments":[{"comment":"There are numerous typographical errors (e.g., 'cathegory', 'Satandard Model', 'autovectors', 'chech'), and Ref. [11] is left incomplete as 'Phys. Reports 14 (1974) incomplete!!!'. The manuscript needs a careful proofread.","section":"Throughout"},{"comment":"The sentence 'Using the value q−1 = 1.14' should be 'q = 1.14' (or 'q−1 = 0.14'); as written it is inconsistent with the q=1.14 obtained in Sec. 4.","section":"Sec. 6"},{"comment":"The comparison in Fig. 3 uses G=3.67 and G=752 with no justification for these values; since G and ε0 are free parameters, the agreement shown in the plots is not a predictive test.","section":"Fig. 3"}],"recommendation":"reject","confidential_remarks":"The central result is appealing, but the gap between Eq. (48) and Eq. (49) is too large for a standard revision: the asserted beta function is structurally inconsistent with the preceding formula, and the numerical q prediction rests entirely on this step. The algebra in Sec. 6 also needs repair. I would need to see an actual derivation of the beta function from the effective coupling and a corrected calculation of the fractal dimension before reconsidering the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"One line: the new thing here is the claim that QCD's Tsallis q is fixed by the one-loop beta function to q=1.14; the supporting calculation is not shown and I don't see how to reproduce it from the preceding equations.\n\nThe paper's opening move, taking the scaling properties of Yang-Mills vertices as a fractal self-similarity and deriving a q-exponential distribution, is mostly a restatement of the thermofractal program the authors have published before. That part is fine as a review, but the genuinely new step is the beta-function matching that turns q into a function of Nc and Nf. That is a nice idea, and if it held up it would be an important result.\n\nIt doesn't yet hold up. The leap from Eq. (48) to Eq. (49) is the whole ballgame, and the paper just says 'easily calculated.' Eq. (48) is a derivative with a q-exponential factor; it is not transparently a g^{N+1} term with a 16π^2 normalization. The stress-test note is right: neither reading of Eq. (48) gives β = -g^3/[16π^2(q-1)] without extra assumptions. The self-similarity condition Eq. (30) is also imposed, not derived from the RG; it forces the exponent but nothing fixes α from the fractal parameters. So the final q=1.14 comes entirely from an identity that is asserted. That is a load-bearing flaw, not a cosmetic one.\n\nThere are also minor slips: Eq. (46) uses M∂Γ/∂λ with mixed symbols, and the text says 'using the value q−1 = 1.14' where it must mean q=1.14. These are annoying but secondary.\n\nThe paper is readable and the ambition is reasonable. It deserves a referee's time in the sense that a full derivation of Eq. (49) would settle whether the result is real. I would send it out rather than desk reject, with an explicit request to show the derivation. In its current form, I would not cite it.\n\nFor a reading group: maybe, as an example of how a plausible-sounding unification can hide a gap.","headline":"The q=1.14 prediction rests on an unshown beta-function identity; the fractal-to-Tsallis part is mostly prior work, but the paper deserves refereeing to force the derivation.","tokens_in":14974,"tokens_out":3493,"would_cite":false,"duration_ms":32824,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81T13","81T17","82B30"],"pacs":["12.38.-t","05.45.Df"],"model":"deepseek-v4-flash","headline":"This paper claims that the Tsallis index of QCD is not a free parameter: fractal self-similarity of Yang-Mills vertices fixes 1/(q-1) = (11/3)c1 - (4/3)c2 = 7, so q = 1.14, matching experiment.","keywords":["Tsallis statistics","non-extensive statistics","quantum chromodynamics","Yang-Mills theory","fractals","self-similarity","renormalization group","particle multiplicity"],"falsifier":"Measure the Tsallis index from transverse-momentum spectra in a process where the one-loop $\\beta$-function coefficient differs from 7, such as a QCD process at very high scale with an effectively different number of flavors or a quenched calculation with $N_f = 0$; if the fitted $q$ remains $1.14$ rather than following $1 + 1/[(11/3)N_c - (2/3)N_f]$, the central identification is wrong. Alternatively, compute truncated $n$-point functions at successive scales on the lattice and test whether $(4N-5) + \\alpha\\nu = \\alpha$ holds.","tokens_in":13881,"feed_emoji":"⚛️","tokens_out":9963,"duration_ms":84505,"temperature":0.7,"pith_summary":"This paper tries to prove that the fractal-like self-similarity of Yang-Mills interactions under scale changes is the physical origin of Tsallis non-extensive statistics in high-energy collisions. Starting from the ordinary phase space of an ideal gas of N partons and imposing that the same distribution repeats at every level of a self-similar hierarchy, it derives the q-exponential distribution $P(\\varepsilon/\\lambda) = [1+(q-1)\\varepsilon/\\lambda]^{-1/(q-1)}$. The paper then identifies the Tsallis index with the coefficient of the one-loop $\\beta$ function, so for QCD with three colors and three flavors $1/(q-1) = 7$, i.e. $q = 1.14$, consistent with the measured value $q = 1.14 \\pm 0.01$. If this is right, the widely used Tsallis parameter is not a phenomenological fit constant but a prediction of QCD's color and flavor content, with further observable consequences for fractal dimension and multiplicity growth.","feed_headline":"One formula fixes QCD's Tsallis index at q = 1.14","feed_subtitle":"Scale-invariant Yang-Mills hierarchies turn a free fit parameter into a prediction from colors and flavors.","key_machinery":"The load-bearing object is the self-similarity condition $(4N-5) + \\alpha\\nu = \\alpha$, which says a parton's energy distribution, written in the scale-free variable $\\chi = \\varepsilon/\\Lambda$, is the same function at every level of the fractal structure. Combined with the approximation $(1-\\varepsilon/M)^{4N-5} \\approx (1+\\varepsilon/M)^{-(4N-5)}$, this converts the ideal-gas power law into the q-exponential form with $1/(q-1) = \\alpha$. The final step compares the resulting $\\beta$ function $\\beta = -\\frac{1}{16\\pi^2}\\frac{1}{q-1}g^{\\tilde N+1}$, with $\\tilde N = 2$, to the one-loop QCD $\\beta$ function $\\beta_{\\rm QCD} = -\\frac{g^3}{16\\pi^2}\\left(\\frac{11}{3}c_1 - \\frac{4}{3}c_2\\right)$, giving $1/(q-1) = 7$.","core_discovery":"The central claim is that renormalizable Yang-Mills field theory generates a fractal-like hierarchy of self-similar vertices, and that this hierarchy is exactly what Tsallis statistics describes. Concretely, the paper obtains $1/(q-1) = (11/3)c_1 - (4/3)c_2$, which for QCD is $7$, so $q = 1 + 1/7 = 1.14$, consistent with $q = 1.14 \\pm 0.01$ from experiment. The same derivation yields a fractal (Hausdorff) dimension $D \\approx 0.69$, matching intermittency measurements, and a multiplicity law $M \\propto E^{1-D} = E^{0.31}$, matching the power-law exponent roughly $0.302$ extracted from proton-proton collisions. In the authors' framing, the q-exponential is not an ad hoc fit function but the effective coupling of the fractal gauge theory.","pith_inferences":["A direct test the paper does not report: fitting Tsallis distributions in processes with a different effective number of flavors should shift $q$ according to $1 + 1/[(11/3)N_c - (2/3)N_f]$; for quenched QCD this would give roughly $q \\approx 1.09$, not $1.14$.","Since the one-loop beta function is used, the derivation suggests $q$ may run with resolution scale as the effective degrees of freedom change, an energy-dependent $q$ not worked out in the paper.","The self-similarity condition $(4N-5) + \\alpha\\nu = \\alpha$ could be checked in a non-perturbative numerical calculation by comparing truncated $n$-point functions at successive scales; if the exponent relation is not scale-independent, the link between $q$ and the beta function would break."],"forward_implications":["The Tsallis index $q$ in high-energy fits becomes a derived quantity: for QCD, $q = 1.14$, so long-tail transverse-momentum spectra need no independent non-extensivity parameter.","The q-exponential acts as an effective coupling in the vertex recursion, suppressing parton energies far above the scale $\\lambda$ and favoring configurations with $\\varepsilon_5 \\approx \\varepsilon_6 \\approx \\lambda/2$.","The fractal dimension of the Yang-Mills hierarchy is fixed, $D \\approx 0.69$, matching the value obtained from intermittency analyses of high-energy distributions.","Particle multiplicity grows with energy as $M \\propto E^{1-D} \\approx E^{0.31}$, consistent with the power-law behavior observed in proton-proton collisions.","For any Yang-Mills theory, the same reasoning determines $q$ from its color and flavor content through the one-loop beta-function coefficient."],"supporting_citations":[{"why":"Supplies the Callan-Symanzik equation whose scale invariance underlies the self-similar/fractal hierarchy of vertices.","marker":"[13–15]"},{"why":"One-loop QCD beta function, the source of the coefficient $(11/3)c_1 - (4/3)c_2$ that fixes $1/(q-1)$.","marker":"[31, 32]"},{"why":"Defines the Tsallis $q$-entropy and q-exponential distribution that the paper derives for parton energies.","marker":"[7]"},{"why":"Earlier thermofractal construction showing that a fractal structure yields non-extensive statistics, used as the interpretive basis for the QCD fractal.","marker":"[17]"},{"why":"Experimental fits giving $q = 1.14 \\pm 0.01$, the value the derivation reproduces.","marker":"[33–35]"},{"why":"Intermittency analyses whose measured fractal dimension $D \\approx 0.69$ is compared with the prediction.","marker":"[39, 40]"},{"why":"Proton-proton multiplicity data with power-law exponent around $0.302$, compared with the predicted $1-D = 0.31$.","marker":"[49]"}],"fun_headline_variants":["Fractals in gauge theory predict q = 1.14 exactly","QCD's q emerges from fractal Yang-Mills hierarchy","No free fit: Tsallis index derived from color and flavor","Fractal dimension D=0.69 appears from gauge self-similarity","q=1.14 emerges from field theory's fractal skeleton"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The argument depends on the assumption that the same scale-free function governs a parton's energy distribution at every level of the fractal hierarchy, together with the approximation that turns $(1-\\varepsilon/M)^{4N-5}$ into $(1+\\varepsilon/M)^{-(4N-5)}$; if that exponent-matching condition is not exact, the q-exponential form and the $\\beta$-function identification do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Fractals in gauge theory predict q = 1.14 exactly","QCD's q emerges from fractal Yang-Mills hierarchy","No free fit: Tsallis index derived from color and flavor","Fractal dimension D=0.69 appears from gauge self-similarity","q=1.14 emerges from field theory's fractal skeleton"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000611,"raw_usage":{"total_tokens":2825,"prompt_tokens":906,"completion_tokens":1919,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":522,"completion_tokens_details":{"reasoning_tokens":1714}},"tokens_in":522,"tokens_out":1919,"duration_ms":12036,"temperature":1.0,"reasoning_tokens":1714,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:49:39.424815+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure the Tsallis index from transverse-momentum spectra in a process where the one-loop $\\beta$-function coefficient differs from 7, such as a QCD process at very high scale with an effectively different number of flavors or a quenched calculation with $N_f = 0$; if the fitted $q$ remains $1.14$ rather than following $1 + 1/[(11/3)N_c - (2/3)N_f]$, the central identification is wrong. Alternatively, compute truncated $n$-point functions at successive scales on the lattice and test whether $(4N-5) + \\alpha\\nu = \\alpha$ holds.","supporting_citations":[{"cited_title":"Deppman, Phys","cited_arxiv_id":null,"evidence_quote":"Earlier thermofractal construction showing that a fractal structure yields non-extensive statistics, used as the interpretive basis for the QCD fractal."},{"cited_title":"Sarkisyan, Aditya Nath Mishra, Raghunath Sahoo and Alexander S","cited_arxiv_id":null,"evidence_quote":"Proton-proton multiplicity data with power-law exponent around $0.302$, compared with the predicted $1-D = 0.31$."}],"review_version":1}