{"id":"3e6c2d71-b84c-4c47-bfd2-3c80f27d865a","arxiv_id":"2411.17090","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Non-extensive HTL corrections increase the Debye mass, weaken the heavy quark potential, and lower the predicted melting temperatures of J/Ψ and Υ, with a magnetic field opposing the effect.","lead":"This paper extends hard thermal loop resummation to Tsallis non-extensive statistics and applies it to heavy quarkonia in a quark-gluon plasma. It finds that non-extensivity increases color screening and lowers quarkonium melting temperatures, while a magnetic field pushes them back up.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"First-order (q−1) expansion is applied at q=1.2 where hard-momentum corrections are O(0.3–1.5) of the local distribution; exact-Tsallis integrals are needed before Table 1 can be trusted.","rationale":"The reader's weakest-assumption diagnosis is exactly the load-bearing issue: a first-order expansion in (q−1) is used to produce quantitative predictions at q=1.2, where the expansion parameter is not small for the hard momenta controlling HTL integrals. I checked the expansion against the exact Tsallis form and confirmed that eqs. (2.5)–(2.6) are the correct first-order terms, but the size of the neglected terms is uncontrolled. This matters because the central numbers — the Debye-mass shifts (3.21)/(3.40) and the Table 1 melting temperatures — are all first-order-in-(q−1) quantities used at (q−1)=0.2. The direction of the effect (increased screening, lower T_melt) is plausible and likely robust, since the exact Tsallis distribution also enhances high-momentum occupancy for q>1, but the magnitude is not. The paper's q=1 results matching known HTL/lattice benchmarks is genuine independent support for the machinery at q=1, but it does not constrain the convergence of the q-expansion. I therefore do not see grounds to reject the paper; the appropriate disposition remains conditional, with the quantitative claims needing either an estimate of higher-order terms or a demonstration that exact-Tsallis results are close to the linearized ones.","tokens_in":39964,"tokens_out":6579,"duration_ms":65495,"concrete_test":"Recompute the zero-chemical-potential retarded and symmetric Debye masses at q=1.1 and q=1.2 using the exact Tsallis distributions, e.g. f_q,B(k) = {[1+(q−1)βk]^{1/(q−1)}}^{-q} − 1? (properly with the q-power in the denominator) and the analogous Fermi form, evaluated numerically in the integrals of eqs. (3.1)–(3.2) and (3.22)–(3.23) for T fixed. Compare the exact ratios m̃²_D,R/m²_D and m̃²_D,F/m²_D with eqs. (3.21) and (3.40); if the difference exceeds ~10% of the claimed (q−1) shift, recompute the Table 1 melting temperatures with the exact Debye masses to see whether the reported T_melt values change by more than their claimed effect.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative claim rests on eqs. (2.4)–(2.6), where the Tsallis distributions are expanded to first order in (q−1), and on the resulting Debye-mass shifts (3.21) and (3.40). The paper itself notes (after eq. 2.5) that the linear expansion holds when k/T is not too large, and asserts that the HTL regime satisfies this. But the HTL integrals receive their dominant contributions from k ≈ 1.5T–3T for the leading term, and from somewhat larger k for the non-extensive correction, where the factor (q−1)(k/T)²/2 is not small: for q=1.2 it is 0.3 at k=√3 T and 0.9 at k=3T, and the relative correction δf_q/f_0 reaches tens of percent at k≈3T and exceeds unity at k≈5T. The numerical scans then use q=1.2, i.e. (q−1)=0.2, while Table 1 reports melting temperatures as sharp numbers with no estimate of the neglected O((q−1)²) terms. Because the claimed 10–15% lowering of T_melt is generated by a first-order correction of the same size, the truncation is not demonstrably controlled. The q=1 benchmark agreement with lattice results shows the machinery is internally consistent, but it does not validate the linearization at q=1.2.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript incorporates Tsallis non-extensive statistics into hard thermal loop (HTL) resummation in the real-time formalism. Starting from the non-extensive quark and gluon distributions (2.1), expanded to first order in (q−1), the authors compute the retarded/advanced and symmetric HTL gluon self-energies and the resummed propagators, both at zero magnetic field and in a finite magnetic field. From the resummed propagators they derive the dielectric permittivity and a complex heavy-quark potential, then solve a Schrödinger equation for the J/Ψ and Υ binding energies and compute decay widths, estimating melting temperatures from the criterion Γ(T_melt)=E_bin(T_melt). The central quantitative results are the Debye-mass shifts (3.21) and (3.40) and Table 1, which show that q>1 lowers the melting temperatures of heavy quarkonia while a magnetic field raises them.","tokens_in":40232,"tokens_out":5702,"duration_ms":54476,"significance":"If the quantitative predictions are reliable, the paper would provide a useful phenomenological extension of HTL resummation that connects non-extensive statistics to quarkonium observables in a QGP. The analytic derivations in Sections 3 and Appendices A–B are carefully laid out, the q→1 limit correctly recovers the standard Debye masses and resummed propagators, and the q=1 benchmark melting temperatures are consistent with the quoted lattice results. The distinction between retarded/advanced and symmetric Debye masses in the non-extensive setting is an interesting and nontrivial result, as is the anisotropic imaginary potential in a magnetic field. However, the quantitative claims for q>1 are not yet controlled because the first-order expansion in (q−1) is used at values as large as q=1.2, where the correction terms are not small; the reported melting temperatures are therefore conditional on a truncation that remains to be justified.","major_comments":[{"comment":"The first-order expansion in (q−1) is used for q=1.1 and q=1.2, but the expansion parameter at the hard momenta that dominate the HTL integrals is not small. The correction to the distribution functions is proportional to (q−1)(k∓μ)^2/(2T^2), which for q=1.2 equals 0.3 at k=√3 T and 0.9 at k=3T; the HTL loop integrals receive their main contributions from precisely this momentum range, so the statement after Eq. (2.5) that the HTL approximation 'satisfies this condition' is not substantiated. Moreover, the symmetric Debye-mass shift in Eq. (3.40) is 1 + (42ζ(3)/π^2 − 3)(q−1) ≈ 1 + 2.115(q−1), which is a 42% increase at q=1.2, while the retarded shift in Eq. (3.21) is 1 + (21ζ(3)/π^2 − 2)(q−1) ≈ 1 + 0.558(q−1). A 42% correction is not a small perturbation, and the neglected O((q−1)^2) terms are not estimated anywhere in the manuscript. Since Table 1 and the central claim that non-extensivity lowers the melting temperatures are produced from these linearized Debye masses, the central quantitative claim is not yet controlled. The authors should either evaluate the exact Tsallis integrals numerically or provide a rigorous truncation-error estimate and restrict the phenomenological conclusions to the range of q where the linearization is demonstrably valid.","section":null},{"comment":"The melting temperatures are quoted to three significant figures with no estimate of uncertainties. The criterion Γ(T_melt)=E_bin(T_melt) is implemented using a Coulomb wave function and the simplified asymptotic form of the real potential described in Section 5; both approximations carry systematic uncertainties that are not quantified, and the q-dependence of those uncertainties is not assessed. The q=1 comparison with lattice QCD is encouraging, but without error estimates the reported differences between q=1, 1.1, and 1.2 (for example, J/Ψ at eB=0: 0.254 → 0.219 GeV) cannot be judged as statistically or systematically significant.","section":null},{"comment":"The non-extensive corrections to the imaginary part of the heavy-quark potential and hence to the decay widths inherit the same uncontrolled linearization. In particular, the symmetric Debye-mass combination appearing in these expressions carries the large coefficient 2.115 in (q−1), so at q=1.2 the imaginary potential is modified by O(40%) corrections while only linear-order terms are retained. Before the conclusion that non-extensivity broadens the decay widths and lowers T_melt is accepted, the authors should demonstrate that higher-order terms in the Tsallis expansion do not change the sign or magnitude of these corrections.","section":null}],"minor_comments":[{"comment":"The caption states that all plots are performed 'at a fixed temperature of T = 0.3 GeV', but the horizontal axes of the same figures are temperature T; please clarify whether 0.3 GeV is a reference scale or remove the phrase.","section":null},{"comment":"There are typographical errors: 'Braatten' should be 'Braaten', and the table label 'T able 1' contains an erroneous space.","section":null},{"comment":"The phrase 'the differece between em2 D,R,B and emD,R reduces' contains a typo ('differece') and also an inconsistent notation: the comparison appears to be between em2 D,R,B and em2 D,R, not between a squared and an unsquared quantity; please make the notation uniform.","section":null},{"comment":"The running coupling is written as α_s(Λ^2, eB) with a logarithm of Λ^2/(Λ^2+eB); since eB has mass dimension two, please specify the units used for eB and clarify the scale-setting prescription for the magnetic-field argument.","section":null}],"recommendation":"major_revision","confidential_remarks":"The paper is within the scope of JHEP and the q=1 benchmark shows that the machinery is internally consistent. The main obstacle is the uncontrolled first-order expansion in (q−1) at q=1.2, which directly affects the quantitative claims in Table 1; this is fixable in principle by computing the exact Tsallis integrals numerically or by deriving a rigorous error bound. I would not recommend rejection, but the resubmission should include such a control calculation or a correspondingly restricted phenomenological range."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Let me give you my read on this paper. The genuinely new thing is a systematic derivation of non-extensive (Tsallis) corrections to the HTL gluon self-energies in the Keldysh real-time formalism, including the splitting between retarded and symmetric Debye masses. The q=1 limit reproduces the standard results, and the appendices show real work. The extension to a magnetic field is also a legitimate addition. That part is solid and worth engaging with.\n\nThe soft spot is exactly where the stress-test note lands. The whole numerical program rests on a first-order expansion in (q-1), and the paper itself admits the expansion holds only when k/T is not too large. But the HTL integrals are dominated by k ~ 1.5T to 3T, and for q=1.2 the correction term (q-1)(k/T)^2/2 is 0.3 to 0.9 there. That is not a small perturbation. So the Debye mass shifts and the melting temperatures in Table 1 are first-order numbers with no control over the neglected terms. The qualitative claim that non-extensivity lowers melting temperatures is plausible, but the quantitative 10-15% lowering is not established.\n\nThere are also standard phenomenological simplifications: Cornell potential, Coulomb wave function, and the Gamma = E_bin criterion for melting. These are common in the literature, but they add model dependence on top of the uncontrolled expansion. No uncertainties are reported anywhere.\n\nThe paper is not a throwaway. The derivation is careful, the references are relevant, and the q=1 comparison with lattice gives a sanity check. But the central quantitative claim needs more than a first-order expansion at q=1.2. I would send it to a serious referee, but the referee should push for exact Tsallis integrals or a controlled estimate of higher-order terms before the melting temperatures can be trusted.","headline":"A careful first-order Tsallis deformation of HTL self-energies, but the uncontrolled expansion at q=1.2 makes the melting temperatures unreliable.","tokens_in":40779,"tokens_out":2458,"would_cite":false,"duration_ms":23552,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Mh","25.75.Nq"],"model":"deepseek-v4-flash","headline":"Non-extensive q-statistics shifts hot-QCD Debye masses and lowers heavy-quarkonium melting temperatures, with a magnetic field pushing them back up.","keywords":["non-extensive statistics","hard thermal loop resummation","quark-gluon plasma","Debye screening mass","heavy quark potential","heavy quarkonium dissociation","magnetic field effects","real-time finite-temperature field theory"],"falsifier":"Recompute the non-extensive HTL self-energies and heavy quark potential to second order in $q-1$ (or with the exact q-deformed distributions) at $q=1.2$ and check whether the Debye-mass shifts and melting temperatures change by more than a perturbatively expected few percent; if the first-order values are not within that tolerance, the predicted lowering of melting temperatures is not established.","tokens_in":39697,"feed_emoji":"🔥","tokens_out":12071,"duration_ms":104718,"temperature":0.7,"pith_summary":"The paper claims that non-extensive statistics changes how a quark-gluon plasma screens color charge, and that this should show up as earlier dissociation of heavy quarkonia. It extends hard thermal loop (HTL) resummation—the standard method for collecting infrared-sensitive thermal corrections into effective propagators—to q-deformed Bose-Einstein and Fermi-Dirac distributions, working to first order in $q-1$ in the real-time formalism. At zero chemical potential the central result is a q-dependent shift of the Debye masses, with the retarded/advanced mass obeying $\\tilde m^2_{D,R}/m^2_D = 1+(21\\zeta(3)/\\pi^2-2)(q-1)$ and the symmetric mass obeying $\\tilde m^2_{D,F}/m^2_D = 1+(42\\zeta(3)/\\pi^2-3)(q-1)$. Feeding these deformed self-energies through the resummed gluon propagator yields a heavy quark potential whose real part is more screened and whose imaginary part is larger as $q$ grows, so binding energies drop, decay widths broaden, and melting temperatures fall; a magnetic field acts in the opposite direction. If the paper is right, quarkonium melting temperatures extracted from ordinary thermal distributions may need to be revised downward in systems that are better described by non-extensive statistics.","feed_headline":"Quarkonium melts earlier in a non-extensive quark-gluon plasma","feed_subtitle":"Non-extensive statistics strengthens screening and decay, while a magnetic field pushes quarkonium melting back up.","key_machinery":"The central object is the q-deformed distribution function, formed by replacing the ordinary exponential with $\\exp_q x=[1+(q-1)x]^{1/(q-1)}$ in the Bose-Einstein and Fermi-Dirac factors (eq. 2.1), and then expanded to first order in $q-1$ (eqs. 2.4–2.6). These distributions enter the real-time bare propagators, and the one-loop HTL integrals over hard momenta turn them into deformed self-energies. The load-bearing identities are the Debye-mass ratios $\\tilde m^2_{D,R}/m^2_D = 1+(21\\zeta(3)/\\pi^2-2)(q-1)$ and $\\tilde m^2_{D,F}/m^2_D = 1+(42\\zeta(3)/\\pi^2-3)(q-1)$; they carry all of the non-extensive physics at leading order. The propagator side of the machinery is the self-consistent resummation equation for the resummed retarded/advanced propagator, whose first-order piece is $\\Pi_{R,(1)}/(G_R^{-1}-\\Pi_{R,(0)})^2$, together with the equation for the symmetric propagator, whose non-extensive correction contains a combination encoding the departure from the equilibrium fluctuation-dissipation relation. This chain—deformed distributions, deformed self-energies, deformed propagators, then dielectric permittivity and potential—is what converts a statistical-mechanical parameter $q$ into an experimentally visible change in quarkonium survival.","core_discovery":"On its own terms, the paper's discovery is that non-extensivity enters the HTL effective theory through a single parameter-dependent shift of the Debye masses, and that the shift is different for the retarded/advanced and symmetric components of the gluon self-energy. With $\\mu=0$, $N_f=3$, $N_c=3$, the paper obtains $\\tilde m^2_{D,R}/m^2_D = 1+(21\\zeta(3)/\\pi^2-2)(q-1)$ and $\\tilde m^2_{D,F}/m^2_D = 1+(42\\zeta(3)/\\pi^2-3)(q-1)$. Since $\\zeta(3)\\approx1.202$, both coefficients are positive, so for $q>1$ both masses increase; the symmetric mass increases more, making $\\tilde m^2_{D,F}>\\tilde m^2_{D,R}$ and breaking the equilibrium equality that connects fluctuation and dissipation. The same deformed self-energies are then used to construct the dielectric permittivity and, through a Fourier convolution of the Coulomb-plus-linear vacuum potential, the in-medium complex heavy quark potential. Solving the quantum mechanical bound-state equation with the real part and folding the imaginary part into the wavefunction, the paper finds $T_{\\rm melt}$ decreases with $q$ for both J/$\\Psi$ and $\\Upsilon$ in zero field, and increases when $eB=15\\,m_\\pi^2$ is switched on: for J/$\\Psi$ it is 0.254 GeV at $q=1$ and $eB=0$, 0.219 GeV at $q=1.2$, and 0.243 GeV at $q=1.2$ with $eB=15\\,m_\\pi^2$.","pith_inferences":["Editorial inference: because only first-order terms in $q-1$ are kept while the numerics use $q=1.2$, the quantitative melting temperatures should be read as indicative; a second-order or exact-$q$ evaluation could shift them by an $O(0.2)$ amount, and this is testable by repeating the calculation.","The paper does not compute radiative quantities, but the split between retarded and symmetric Debye masses implies that the photon and dilepton emission rate, which is controlled by the symmetric propagator, should also carry a $q$-dependent enhancement; measuring the dilepton spectrum could give an independent handle on $q$ in the plasma.","Since $q$ lowers $T_{\\rm melt}$ and $eB$ raises it, the two effects could partially cancel; mapping the dissociation boundary in the $(q,eB)$ plane for J/$\\Psi$ and $\\Upsilon$ would separate non-extensive effects from magnetic-field effects in heavy-ion phenomenology."],"forward_implications":["For any $q>1$, the retarded Debye mass $\\tilde m^2_{D,R}$ exceeds the standard $m^2_D$, so the medium screens the color Coulomb interaction more strongly; the real part of the heavy quark potential flattens and binding energies drop.","The symmetric Debye mass $\\tilde m^2_{D,F}$ grows even faster than the retarded one, so the fluctuation-dissipation relation between the symmetric self-energy and the retarded/advanced ones is violated at first order in $q-1$.","Larger $q$ increases the magnitude of the imaginary part of the potential, broadening quarkonium decay widths; with both smaller binding and larger width, J/$\\Psi$ and $\\Upsilon$ melt at lower temperature.","A magnetic field $eB=15\\,m_\\pi^2$ raises the melting temperatures of both J/$\\Psi$ and $\\Upsilon$ at fixed $q$; at $q=1$ it shifts J/$\\Psi$ from 0.254 to 0.270 GeV, and at $q=1.2$ from 0.219 to 0.243 GeV."],"supporting_citations":[{"why":"Supplies the HTL resummation technique, the perturbative method this paper extends to non-extensive distributions.","marker":"[9, 16]"},{"why":"Supplies the non-extensive distribution functions and the real-time bare propagators that form the starting point of the calculation.","marker":"[64, 65]"},{"why":"Supplies the dielectric-permittivity and Fourier-convolution method that converts resummed gluon propagators into an in-medium heavy quark potential.","marker":"[78–80]"},{"why":"Provides the magnetic-field one-loop quark contribution to the gluon self-energy and the magnetic-field quantized-level structure used in the finite-field case.","marker":"[23]"},{"why":"Supplies the melting criterion $\\Gamma(T_{\\rm melt})=E_{\\rm bin}(T_{\\rm melt})$ used to define the melting temperature.","marker":"[86]"},{"why":"Provides lattice-QCD melting temperature values against which the paper benchmarks its $q=1$ results.","marker":"[87]"}],"fun_headline_variants":["Non-extensive QGP lowers quarkonium melt temperature","Magnetic field offsets non-extensive quarkonium melting","Debye mass split breaks dissipation-equilibrium link in QGP","Non-extensive corrections weaken quarkonium binding energy","Stronger screening from non-extensive HTL melts quarkonia at lower T"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The calculation is linearized in $q-1$, and the numerical scans at $q=1.2$ assume terms of order $(q-1)^2$ are negligible even though $(q-1)=0.2$ makes them comparable in size to the retained first-order corrections.","fun_headline_variants_meta":{"raw":{"variants":["Non-extensive QGP lowers quarkonium melt temperature","Magnetic field offsets non-extensive quarkonium melting","Debye mass split breaks dissipation-equilibrium link in QGP","Non-extensive corrections weaken quarkonium binding energy","Stronger screening from non-extensive HTL melts quarkonia at lower T"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000835,"raw_usage":{"total_tokens":3770,"prompt_tokens":1198,"completion_tokens":2572,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":814,"completion_tokens_details":{"reasoning_tokens":2490}},"tokens_in":814,"tokens_out":2572,"duration_ms":18322,"temperature":1.0,"reasoning_tokens":2490,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T12:31:06.335438+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the non-extensive HTL self-energies and heavy quark potential to second order in $q-1$ (or with the exact q-deformed distributions) at $q=1.2$ and check whether the Debye-mass shifts and melting temperatures change by more than a perturbatively expected few percent; if the first-order values are not within that tolerance, the predicted lowering of melting temperatures is not established.","supporting_citations":[{"cited_title":"Responses of quark-antiquark interaction and heavy quark dynamics to magnetic field","cited_arxiv_id":"2301.09110","evidence_quote":"Provides the magnetic-field one-loop quark contribution to the gluon self-energy and the magnetic-field quantized-level structure used in the finite-field case."}],"review_version":1}