{"id":"d805d9a6-f713-4ae9-996d-d5d956c81e7c","arxiv_id":"2501.12777","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Excited bottomonia in a quark-gluon plasma acquire temperature-dependent widths with no significant mass shifts, according to a Gaussian model fit to lattice NRQCD correlators.","lead":"Using a fine lattice and non-relativistic QCD, the authors study excited bottomonium states in a hot quark-gluon plasma up to T=250 MeV. They find the states keep their vacuum masses while their thermal width grows, a key observable for heavy-ion collision experiments.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The reported in-medium widths are fit parameters of an unmotivated Gaussian spectral ansatz (Eq. 6); a Lorentzian or asymmetric peak can reproduce the same continuum-subtracted correlators, so the 'confirm nonzero thermal widths' claim and the screening conclusion are not yet established.","rationale":"The reader's weakest assumption identifies the load-bearing concern correctly: the in-medium width is not directly observable, and Eq. 6 imposes a specific line shape that is acknowledged by the authors to be unmotivated. The paper is transparent about this limitation and defers to the companion paper for Breit-Wigner comparisons, but the proceedings itself goes beyond the evidence when it says the results 'confirm nonzero thermal widths'. A Lorentzian refit of the same correlators is the natural arbiter: if the width and mass survive that change of ansatz, the qualitative conclusion is robust; if not, the screening conclusion cannot be sustained. Since the reader's CONDITIONAL verdict already captures this state of affairs, I do not recommend changing the verdict. The abstract could be softened to 'indicate' instead of 'confirm' until the companion analysis is public, but that is a wording adjustment rather than a change in verdict.","tokens_in":11360,"tokens_out":6060,"duration_ms":68857,"concrete_test":"Re-fit the same continuum-subtracted correlators used for Fig. 4 with a Breit-Wigner spectral function rho_med(omega) = A Gamma / [(omega-M)^2 + Gamma^2] + A_cut delta(omega-omega_cut), using identical fit ranges and the same number of parameters, for every channel and temperature. If the extracted Lorentzian M and Gamma agree with M_med and sqrt(2 ln 2) Gamma_med within quoted uncertainties, the Gaussian choice is not load-bearing; if they differ by more than the errors, the 'nonzero thermal widths' and the screening conclusion are artifacts of the Gaussian shape.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central quantitative result — 'nonzero thermal widths' — is not a measured quantity but a parameter of the Gaussian ansatz in Eq. 6, rho_med = A_med exp[-(omega-M_med)^2/(2 Gamma_med^2)] + A_cut delta(omega-omega_cut). Under this ansatz, the continuum-subtracted effective mass obeys M_sub_eff(tau) approximately M_med - Gamma_med^2 tau plus tail terms, so the reported width is essentially the slope of M_sub_eff in the middle-tau region. Any broadened spectral shape — Lorentzian, skewed peak, or a peak plus a temperature-dependent low-energy component — can produce a similar near-linear slope. The authors themselves state that the Gaussian form is 'not physically motivated' and that the interpretation of Gaussian widths as bottomonium widths is 'not well-defined', yet the abstract says the results 'confirm nonzero thermal widths'. The mass-shift conclusion is less fragile because it tracks the peak position, but M_med is extracted from the same five-parameter fit and could be biased by the ansatz. The screening conclusion in Section 4 therefore rests on an unvalidated shape assumption: it is plausible but not established by this proceedings.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This proceedings contribution reports a lattice NRQCD study of bottomonium correlators and their in-medium properties at temperatures from T ≈ 133 to 250 MeV, using a fine lattice spacing a = 0.0493 fm with (2+1)-flavor HISQ gauge configurations near the physical point. Extended operators (Gaussian-smeared and wave-function-optimized) are used to improve overlap with the ground and excited states up to 3S and 3P. The authors extract in-medium masses and widths by fitting a continuum-subtracted correlator with a Gaussian-plus-delta spectral ansatz (Eq. 6). They report that the in-medium masses are consistent with vacuum masses across the temperature range, while the extracted widths increase with temperature and are larger for excited states. From this they conclude that screening may not be the likely source of bottomonium dissociation in the medium. The paper also checks that the results are insensitive to the choice of extended operators.","tokens_in":11522,"tokens_out":7327,"duration_ms":71903,"significance":"If the central claims are reliable, the results are a valuable addition to the quarkonium-in-medium literature: they extend lattice NRQCD determinations of in-medium properties to the 3S and 3P excited states, which are rarely accessed, and they use a fine lattice spacing with (near-)physical light quarks. The use of two independent extended-operator constructions and the explicit check of operator dependence are strengths. The zero mass shift and the increasing width hierarchy are consistent with the general expectation of sequential thermal broadening and provide a useful constraint on model calculations. However, because the width is defined through a specific Gaussian ansatz that the paper itself calls 'not physically motivated,' the quantitative claims and the screening conclusion are only as strong as that model assumption. The manuscript acknowledges this limitation in the body but the abstract and conclusion state the results more categorically. As a proceedings paper, the contribution is worth publishing after the framing of the claims is brought into line with the method's uncertainties.","major_comments":[{"comment":"The abstract states 'Our results confirm nonzero thermal widths' and Sec. 4 repeats that 'nonzero widths are observed,' but the width is a fit parameter of the Gaussian-plus-delta ansatz in Eq. (6), which the text itself describes as 'not physically motivated' with an interpretation that 'is not well-defined.' Since alternative spectral shapes (e.g., Lorentzian or asymmetric peaks) can reproduce the same continuum-subtracted correlator, the data as presented do not confirm a nonzero physical width in a model-independent sense. Please rephrase the headline claims to explicitly incorporate the ansatz dependence, or add a quantitative test showing that the extracted width is stable under plausible alternative parameterizations; the screening conclusion in Sec. 4 rests on this point.","section":"Abstract and Sec. 4"},{"comment":"The continuum-subtraction procedure assumes a single temperature-independent exponential C_cont(τ) obtained from a vacuum fit. If the in-medium spectral function develops a low-frequency tail or an asymmetric broadening, this subtraction can distort M_sub_eff and generate an artificial linear slope that is then interpreted as a width. The statement that small-τ effective masses are temperature independent only constrains the high-frequency part of the spectrum and does not establish temperature independence of the continuum in the frequency region relevant to the subtraction. The authors should discuss this systematic uncertainty or explicitly refer to tests in Ref. [27], since it directly affects both the extracted width and the mass shift.","section":"Sec. 3.2, Eqs. (3)-(4)"},{"comment":"The five-parameter ansatz (A_med, M_med, Γ_med, A_cut, ω_cut) is fitted to a correlator with at most sixteen time slices at the highest temperature, leaving few degrees of freedom. The paper does not report the number of fitted points, χ²/dof, or the sensitivity of Γ_med to the fit-range boundaries and to the omission of the largest-τ points. The special treatment at T = 167 MeV (setting A_cut = 0 and omitting 2–4 points) is described, but no analogous stability information is given for the other temperatures. To support the quantitative widths in Fig. 4, a fit-stability analysis (or a clear reference to such an analysis in the companion paper) is needed.","section":"Sec. 3.2, Eq. (6) and Fig. 4"},{"comment":"The conclusion that 'all bottomonium states below the open-bottom threshold can exist as well-defined quasi-states above the crossover temperature, including 3P states' is stronger than the evidence shown in this proceedings. For the 3S and 3P states at high temperature, the effective masses in Fig. 2 show a steep drop without a discernible plateau, and a Gaussian-fitted M_med alone does not establish a well-defined quasi-particle peak. A quantitative criterion (e.g., M_med > Γ_med, or a clear separation between the peak and the continuum contribution) should be stated, or the conclusion should be softened to say that the data are consistent with such quasi-states within the assumed Gaussian parameterization.","section":"Sec. 4"}],"minor_comments":[{"comment":"Subject-verb agreement: 'The temperature dependence ... are presented' should be 'is presented.'","section":"Abstract and Sec. 1"},{"comment":"The definition of τ_vac_min is internally inconsistent: the text says δM(τ) is 'less than the statistical uncertainty' but then gives 'δM(τ) < 25% × σ_Meff(τ).' Please clarify which criterion was used.","section":"Sec. 3.1"},{"comment":"The gray-shaded area is referenced in the text but not explained in the Fig. 4 caption; the reader cannot tell which temperature regions correspond to the single-exponential fits. Please describe it in the caption or in the text.","section":"Sec. 3.2 and Fig. 4"},{"comment":"The 'width at half maximum height' is defined as √(2 ln 2) Γ_med, which is actually the half-width at half maximum (HWHM), not the full width at half maximum (FWHM). Please clarify the convention to avoid confusion with the usual physical width.","section":"Sec. 3.2, after Eq. (6)"},{"comment":"The phrase 'the rotation matrix Ω_α,ij is computed at zero temperature and uniformly applied across all temperatures' is slightly awkward; consider 'computed at zero temperature and then used at all temperatures.'","section":"Sec. 2"}],"recommendation":"major_revision","confidential_remarks":"This manuscript is a proceedings contribution based on the companion paper Ref. [27] (arXiv:2501.11257), which presumably contains the full systematic analysis. The data and the operator construction are commendable, but the abstract and conclusion overstate the model-independence of the extracted widths. The requested revisions are within the scope of a proceedings: rephrasing the central claims, making the fit-stability information explicit (or referencing it precisely), and softening the 'well-defined quasi-states' statement. If the companion paper indeed provides the missing systematics, a properly framed proceedings version would be acceptable."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First thing to know: this is a proceedings write-up of the same group's companion paper (Ref. [27]), and the central in-medium numbers already appear there. What this contribution actually adds is the fine lattice spacing (0.0493 fm) at every temperature, wave-function-optimized operators covering up to 3P states, and a consistency check between two operator choices. The vacuum spectrum is benchmarked against PDG and the correlator data look solid. That is real work and worth a look.\n\nThe soft spot is exactly where the stress-test points. The abstract says the results 'confirm nonzero thermal widths.' That overstates what a fit can deliver. The widths are parameters of the Gaussian ansatz in Eq. (6), which the authors themselves describe as 'not physically motivated' and whose width interpretation they call 'not well-defined.' The continuum-subtracted effective mass has a near-linear slope in the middle time range, and any broadened peak shape — Lorentzian, skewed, peak plus a low-energy component — can reproduce that slope. So the width numbers are not a measured property of the states yet. The mass-shift result is more robust, since it follows the peak position, but it still comes from the same five-parameter fit with unquantified systematic uncertainties. The conclusion that screening is not the main dissociation mechanism is plausible, but it rests on the shape assumption and on checks in the companion paper, not on anything this proceedings establishes by itself.\n\nTo be fair, the body is transparent about the ansatz problem. The gap is between the careful caveat in Sec. 3.2 and the confident abstract/conclusion wording. For a proceedings, this is an acceptable progress report, and the underlying calculation deserves serious referee time. If I were the editor, I would send it out with the request that the abstract be reworded and the fit-model systematics either shown or explicitly referenced with numbers. The companion paper is where the real claim will stand or fall.","headline":"Fine-lattice NRQCD progress report; the width claim is a Gaussian fit parametrization, not a measured quantity, so the abstract overstates the data.","tokens_in":12229,"tokens_out":3039,"would_cite":false,"duration_ms":29337,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Gc","12.38.Mh","25.75.Nq"],"model":"deepseek-v4-flash","headline":"Lattice QCD shows excited bottomonia keep their vacuum masses up to 250 MeV while gaining thermal widths, indicating broadening rather than screening drives dissociation.","keywords":["bottomonium","quark-gluon plasma","lattice NRQCD","thermal width","in-medium mass","excited states","spectral function","quarkonium suppression"],"falsifier":"Fit the same continuum-subtracted correlators with a Breit-Wigner (Lorentzian) in-medium peak, the parameterization the authors describe as more natural and say they address in their companion paper. If the Lorentzian fit yields nonzero mass shifts, or widths that differ from the Gaussian widths by more than the statistical uncertainties, the paper's central claims would be contradicted; an even more direct test is a model-independent reconstruction of $\\rho(\\omega,T)$ from the same correlators.","tokens_in":11025,"feed_emoji":"⚛️","tokens_out":12711,"duration_ms":101749,"temperature":0.7,"pith_summary":"Excited bottomonia are heavy quark-antiquark states whose survival in the quark-gluon plasma is a classic test of color screening. This paper uses lattice QCD with non-relativistic bottom quarks to track the 1S, 2S, 3S and 1P, 2P, 3P bottomonium states at temperatures from 133 to 250 MeV. It finds that, within uncertainties, all of these states keep their vacuum masses while acquiring thermal widths that grow with temperature and with the size of the state. If correct, this means screening of the quark-antiquark force is not what dissolves these states; instead, the dominant effect is in-medium broadening that damps the bound state. That would change how quarkonium suppression in heavy-ion collisions is interpreted.","feed_headline":"Hot plasma broadens bottomonia but leaves their masses intact","feed_subtitle":"Up to 3S and 3P states persist with vacuum masses, pointing to broadening, not screening.","key_machinery":"The machinery is the NRQCD Euclidean correlator $C(\\tau,T) = \\int_{-\\infty}^{+\\infty} d\\omega\\, \\rho(\\omega,T) e^{-\\tau\\omega}$, which gives a temporal range of $1/T$ rather than $1/(2T)$ and thus better sensitivity to thermal effects. Extended operators—Gaussian-smeared sources and wave-function-optimized sources obtained from solving a discretized three-dimensional bound-state wave equation—provide clean overlap with the targeted excited states. The in-medium spectral function is modeled as $\\rho_{\\mathrm{med}}(\\omega,T) = A_{\\mathrm{med}}(T)\\exp[-(\\omega - M_{\\mathrm{med}}(T))^2/(2\\Gamma_{\\mathrm{med}}^2)] + A_{\\mathrm{cut}}(T)\\delta(\\omega - \\omega_{\\mathrm{cut}}(T))$, with the vacuum continuum removed by subtracting a single exponential $A e^{-M\\tau}$ from the zero-temperature correlator. Fitting this ansatz to the continuum-subtracted correlators yields the in-medium mass $M_{\\mathrm{med}}(T)$ and width $\\Gamma_{\\mathrm{med}}(T)$ that carry all the paper's conclusions.","core_discovery":"The paper's central claim is that the in-medium masses of $\\Upsilon(nS)$ and $\\chi_{b0}(nP)$ for $n = 1, 2, 3$ are equal to their vacuum values within uncertainties across $T \\simeq 133$–$250$ MeV, while the extracted widths increase with temperature and obey the hierarchy $\\Gamma(1S) < \\Gamma(2S) < \\Gamma(3S)$ and $\\Gamma(1P) < \\Gamma(2P) < \\Gamma(3P)$. The authors take the near-zero mass shifts together with nonzero widths as evidence that screening may not be the likely source of bottomonium dissociation in the medium. They further show that these in-medium properties do not depend on which of the two extended-operator constructions is used, indicating that the extraction is not sensitive to the operator choice.","pith_inferences":["If the quasi-state picture holds, the suppression hierarchy observed in heavy-ion experiments could be read as a dynamical rate effect—larger widths for higher excited states mean faster dissociation—rather than as static screening by a shorter screening length.","Because the paper concedes the Gaussian ansatz is not physically motivated, the reported widths are plausibly qualitative upper bounds; a Breit-Wigner or fully reconstructed spectral function could change the width values without disturbing the near-zero mass shifts.","A natural next step would be to convert the extracted widths into dissociation rates and compare with open-quantum-system models of quarkonium in the plasma, making the broadening-versus-screening distinction testable against experimental yields."],"forward_implications":["Excited bottomonia up to the 3S and 3P states persist as well-defined quasi-states, with masses equal to their vacuum values, across the whole temperature range $T \\simeq 133$–$250$ MeV.","Thermal widths grow with temperature and follow the size hierarchy $\\Gamma(1S) < \\Gamma(2S) < \\Gamma(3S)$ and $\\Gamma(1P) < \\Gamma(2P) < \\Gamma(3P)$, giving a quantitative sequential-broadening pattern.","The combination of unchanged masses and growing widths implies that screening of the real part of the quark-antiquark potential is not the dominant cause of bottomonium dissociation; in-medium damping is.","The extracted in-medium properties are insensitive to the choice between Gaussian-smeared and wave-function-optimized operators, so the results are not an artifact of the operator construction."],"supporting_citations":[{"why":"Derives the NRQCD correlator relation used in Eq. 2, doubling the usable temporal range.","marker":"[21]"},{"why":"Introduces the Gaussian-smeared extended operators and the Gaussian spectral-function ansatz used for the in-medium extraction.","marker":"[25]"},{"why":"Introduces the wave-function-optimized operators and the continuum-subtraction procedure for excited bottomonia.","marker":"[26]"},{"why":"The companion paper where the Breit-Wigner parameterization and fuller analysis are presented.","marker":"[27]"},{"why":"Shows the finite-temperature quark-antiquark potential has a large imaginary part, the dissipative effect the authors contrast with screening.","marker":"[11]"},{"why":"Lattice evidence that the real part of the static potential shows no screening, supporting the paper's conclusion.","marker":"[13]"},{"why":"Reports unscreened forces in the quark-gluon plasma, reinforcing the same interpretation.","marker":"[14]"},{"why":"Experimental observation of excited Upsilon suppression in heavy-ion collisions that motivates the study.","marker":"[5]"}],"fun_headline_variants":["Hot QCD broadens bottomonia but holds masses flat","Excited bottomonia widen in plasma, masses stay put","Thermal widths grow for bottomonia up to 3S and 3P","No mass shift for bottomonia in hot matter, only width","Operator choice doesn't alter bottomonium thermal widths"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The extraction assumes the in-medium part of the spectral function is a Gaussian peak plus a delta tail, and that the vacuum continuum is a single exponential with fitted amplitude and mass; if the true spectral shape differs from this ansatz, the reported widths and the conclusion that screening is not the dissociation mechanism do not follow.","fun_headline_variants_meta":{"raw":{"variants":["Hot QCD broadens bottomonia but holds masses flat","Excited bottomonia widen in plasma, masses stay put","Thermal widths grow for bottomonia up to 3S and 3P","No mass shift for bottomonia in hot matter, only width","Operator choice doesn't alter bottomonium thermal widths"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000153,"raw_usage":{"total_tokens":1187,"prompt_tokens":907,"completion_tokens":280,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":523,"completion_tokens_details":{"reasoning_tokens":195}},"tokens_in":523,"tokens_out":280,"duration_ms":3480,"temperature":1.0,"reasoning_tokens":195,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-10T16:48:45.747222+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fit the same continuum-subtracted correlators with a Breit-Wigner (Lorentzian) in-medium peak, the parameterization the authors describe as more natural and say they address in their companion paper. If the Lorentzian fit yields nonzero mass shifts, or widths that differ from the Gaussian widths by more than the statistical uncertainties, the paper's central claims would be contradicted; an even more direct test is a model-independent reconstruction of $\\rho(\\omega,T)$ from the same correlators.","supporting_citations":[],"review_version":1}