{"id":"ad696cc9-ee3d-4aa7-9256-b369a5db11a0","arxiv_id":"2412.15414","paper_version":2,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"Using the screened massive expansion at one loop, the authors derive a qualitative QCD phase boundary from the maximum of the longitudinal gluon propagator, with the transition temperature vanishing at the lightest infrared quark mass.","lead":"The authors compute the gluon propagator in full QCD at finite temperature and baryon density using a massive perturbation expansion, then use the temperature maximum of its longitudinal component as a deconfinement marker. They obtain a qualitative QCD phase diagram in which the transition temperature falls with chemical potential and vanishes at the lightest effective quark mass.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The load-bearing premise that the longitudinal gluon propagator's maximum marks deconfinement in full QCD is extrapolated from pure Yang-Mills and remains untested; the phase boundary inherits this premise.","rationale":"The reader's conditional verdict already identifies the Tmax-to-deconfinement correspondence as a load-bearing premise, and I agree that this is the most consequential assumption. The screened massive expansion calculation itself is explicit, internally consistent, and the stability analysis is thorough; credit is due for the honest discussion of the absolute scale and of the mass-drop ambiguity. However, the central claim—that Tmax(μ) traces the QCD phase boundary—cannot be established by the one-loop calculation alone, because the criterion is imported from pure Yang-Mills without an independent full-QCD test. The quark-mass modeling is a real limitation, but it is secondary: even with perfect quark masses, the result would only be meaningful if the gluon-propagator maximum actually marks deconfinement in full QCD. The proposed lattice check at μ=0 is feasible and would directly test that premise. Since the paper itself frames the result as conditional on this parallelism, and the reader's verdict already reflects that conditionality, no change to the verdict is needed.","tokens_in":36756,"tokens_out":4736,"duration_ms":45071,"concrete_test":"Compute the Landau-gauge, zero-Matsubara longitudinal gluon propagator Δ_L(ω=0, p→0; T) in 2+1-flavor lattice QCD with physical quark masses at μ=0 and locate its maximum as a function of T. If the maximum occurs near the established crossover region (T ≈ 155–175 MeV), the Tmax criterion survives in full QCD; if it occurs near 77 MeV or is monotonic, the central claim fails. A finite-density variant would repeat the exercise at imaginary chemical potential and analytically continue to real μ, comparing the resulting Tmax(μ) with the lattice chiral crossover curvature.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central move is made in Sec. IVA: 'It is then natural to interpret an analogous change in full QCD as marking either a true phase transition, or a smooth crossover between phases.' The only empirical anchor is quenched lattice QCD (Ref. [82]), where the maximum of the longitudinal propagator coincides with the first-order Yang-Mills transition. Full QCD at physical quark masses is a crossover, with pseudo-critical temperatures from chiral susceptibilities and Polyakov-loop observables in the 155–175 MeV range. The paper's Tmax(0) = 0.117 m0 ≈ 77 MeV is about half those values, and the authors explicitly disclaim quantitative reliability of the absolute scale. The ratio argument—T_YM/T_QCD ≈ 1.57 against the expected 1.54–1.74—is suggestive but does not establish the criterion: a scheme whose one-loop YM transition is also roughly half the physical value will plausibly reproduce a similar ratio even if the gluon-propagator peak in full QCD has no direct relation to deconfinement. Furthermore, the mass-drop procedure in Sec. IVA uses the Tmax curve itself to define the confined region and then imposes lighter quark masses outside it, so the endpoint μ_c = M1 and the disappearance of the humps are consequences of the assumed criterion plus an input mass. If the Tmax-to-deconfinement correspondence fails in full QCD, the entire phase boundary—including Tc(0) and the endpoint—is unsubstantiated. This is an extrapolation across theories rather than an internal inconsistency, but it is the least secure link in the central claim.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript extends the authors' screened massive expansion of the Landau-gauge gluon propagator from pure Yang-Mills theory to full QCD with 2+1 dynamical quarks at finite temperature T and baryon chemical potential μ. Quarks are modeled as massive fields with constant effective infrared masses M1=350 MeV and M2=450 MeV, and the one-loop polarization is evaluated analytically up to a one-dimensional integral. At zero Matsubara frequency and vanishing spatial momentum, the longitudinal propagator is found to be non-monotonic in T for μ<M1, with Tmax(0)≈0.117 m0≈77 MeV; Tmax decreases with μ and extrapolates to zero at μ→M1. Interpreting Tmax as the deconfinement temperature yields a phase diagram in the (T,μ) plane whose normalized shape is stable under parameter changes. The authors clearly state that the absolute scale is not quantitatively reliable.","tokens_in":37089,"tokens_out":6750,"duration_ms":62128,"significance":"The paper's technical core is sound: the finite-density extension of the quark-loop polarization is standard, the one-dimensional integral representation is useful, and the stability analysis in Sec. IIID covers the free parameters of the expansion. The main conceptual contribution is a falsifiable prediction for the shape Tc(μ)/Tc(0) for μ<M1, including a change of concavity that could hint at a critical endpoint, and a simple explanation of why the endpoint sits at the lightest infrared quark mass. The significance is conditional on the assumed correspondence between the longitudinal gluon propagator's maximum and deconfinement in full QCD; this correspondence is presently established only in pure Yang-Mills lattice simulations. The authors are honest about this limitation, and the analytic expressions and explicit caveats are strengths.","major_comments":[{"comment":"The phase boundary rests on the unvalidated assumption that the maximum of the zero-frequency longitudinal gluon propagator marks deconfinement in full QCD. In pure Yang-Mills the criterion is supported by lattice data [82], but in the presence of dynamical quarks the transition is a crossover and the gluon-propagator maximum has no independent lattice or functional check. The ratio T_YM/T_QCD ≈ 1.57 quoted in Sec. IVA is not a discriminating test: the one-loop calculation with temperature-independent parameters underestimates the pure-Yang-Mills Tc by about a factor of two (Sec. IIB), so a common systematic offset can produce the same ratio even if the criterion fails in full QCD. A concrete validation would be to compute Tmax from unquenched gluon-propagator data at μ=0 and compare it with the chiral-susceptibility or Polyakov-loop crossover temperature; until then, the claim should be presented as a model conjecture rather than a QCD phase-diagram prediction.","section":"Sec. IVA, Eq. (38), Fig. 13"},{"comment":"The endpoint of the phase boundary is an input, not a derived result. The effective quark masses M1=350 MeV and M2=450 MeV are free parameters, and Tmax(μ)→0 as μ→M1 follows from the Fermi-distribution threshold in Eq. (43) once M1 is chosen. Similarly, the deconfined-phase masses (M1',M2')=(125,225) MeV are imposed by hand in Sec. IVA, and the disappearance of the humps for μ>M1 in Fig. 13 is a consequence of that mass drop. The model therefore does not independently predict the absolute location μc; it predicts the normalized curve Tc(μ)/Tc(0) and its robustness. The authors should either implement the self-consistent determination of Mf(T,μ) they mention in Sec. IVA, or explicitly restrict the parametric predictions to μ<M1 and to statements about shape.","section":"Secs. IIIA and IVA, Eq. (43), Fig. 12"},{"comment":"The absolute temperature scale Tc(0)=0.117 m0≈77 MeV is set by m0=656 MeV, which is taken from a fit to zero-temperature pure-Yang-Mills lattice data [13]. Because the calculation is one-loop and uses temperature-independent parameters, the analogous pure-Yang-Mills result is Tc≈121 MeV, about half of the physical 270 MeV (Sec. IIB). The authors acknowledge this in Sec. IVA, but the abstract and title still present the result as 'the QCD phase diagram'; the wording should make clear that only the normalized shape, and not the absolute position of the crossover line, is claimed.","section":"Sec. IIIA, Sec. IVA"}],"minor_comments":[{"comment":"The caption says the starred value 10 MeV for the longitudinal m(T) at T=260 MeV is the lowest value that could be reached by the numerical routines, but the text says the routines could not reach lower values; clarify whether 10 MeV is a converged fit value or a numerical bound.","section":"Sec. IIB, Table I"},{"comment":"The panel labels use 'Tc(µ)' for the value Tmax(µ) before the notation is defined; define Tc(µ) in the captions or use Tmax consistently.","section":"Figs. 8 and 9"},{"comment":"The prefactor p²/p² in Eq. (33) is redundant and should be simplified to avoid confusion.","section":"Eq. (33)"},{"comment":"The symbol μ is used both for the renormalization scale in Eq. (31) and for the chemical potential throughout the paper; use a distinct symbol such as μ_R for the renormalization scale.","section":"Eq. (31)"},{"comment":"The phrase 'around twice as small as' should read 'about half of', and several garbled mathematical symbols appear in the text (e.g., '/greaterorapproxeql'); a careful proofreading pass is needed.","section":"Sec. IVA, Sec. V"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is transparent and well organized, and the technical calculation appears sound. My main concern is fit-to-scope rather than correctness: the title and abstract promise a phase diagram, while the load-bearing criterion is an unvalidated extrapolation and the endpoint is an input. I would support publication after a major revision in which the claims are reframed as a model-dependent study and, if possible, a concrete test of the Tmax criterion in full QCD is added or explicitly identified as an open problem."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The new thing here is the finite-density extension of the screened massive expansion to full QCD with quark loops, and the phase diagram is read off from the maximum of the longitudinal gluon propagator at zero Matsubara frequency. The calculation is explicit, the analytic expressions are reduced to a one-dimensional integral, and the stability analysis across the free parameters is genuinely thorough. The authors also correct an error in their earlier code, which is a good sign of care.\n\nWhat the paper does well: the normalized curve Tc(mu)/Tc(0) is robust across mass configurations, the Yang-Mills to full-QCD transition-temperature ratio comes out close to the lattice result, and the qualitative shape—monotonic decrease with mu and endpoint at the lightest quark mass—matches expectations. The authors are candid about the absolute scale being roughly half the known crossover temperature and about the fragility of the humps in Tmax(mu).\n\nThe soft spots are real but mostly acknowledged. The load-bearing premise is that the longitudinal propagator maximum marks deconfinement in full QCD; that is extrapolated from quenched lattice Yang-Mills, and the paper explicitly says it is an assumption. In full QCD at physical quark masses the transition is a crossover, so the correspondence is not automatic. The endpoint mu_c = M1 is literally an input, and the mass-drop procedure uses the Tmax curve itself to define the confined region before imposing lighter masses outside it. So the phase boundary is contingent on the modeling, not a first-principles output. The authors say so themselves, which keeps this from being misleading.\n\nIs the central argument holding up? As a derived quantitative prediction, no; as a qualitative analytic estimate of the phase boundary shape and a ratio, yes. The calculation is sound within its one-loop framework, and the weakness is the interpretive link, not the algebra.\n\nWho benefits: people working on massive expansion or Curci-Ferrari-type approaches, and anyone wanting an analytic route to the qualitative QCD phase boundary with a clear statement of its assumptions. It deserves a serious referee, not a desk reject; a referee should press on the Tmax criterion in full QCD and on the dynamical treatment of quark masses, but the paper already frames those as open problems.","headline":"A careful one-loop extension of the screened massive expansion to finite density QCD with quark loops, whose phase boundary rests on an imported criterion that the paper itself flags as an assumption.","tokens_in":37616,"tokens_out":1698,"would_cite":true,"duration_ms":17833,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["81V05","81T28"],"pacs":["12.38.-t","11.10.Wx"],"model":"deepseek-v4-flash","headline":"The temperature at which the longitudinal gluon propagator peaks traces the QCD deconfinement boundary from about 77 MeV at zero density to zero at the light quark mass.","keywords":["QCD phase diagram","gluon propagator","screened massive expansion","deconfinement","finite temperature field theory","finite chemical potential","Landau gauge","chiral symmetry restoration"],"falsifier":"A lattice calculation of the unquenched zero-frequency longitudinal gluon propagator in 2+1-flavor QCD at vanishing momentum that finds no maximum near T≈77 MeV at μ=0, or a gap-equation computation in which Tmax(μ) does not vanish at μ=M1, would disprove the claim.","tokens_in":36505,"feed_emoji":"⚛️","tokens_out":5691,"duration_ms":47908,"temperature":0.7,"pith_summary":"Using a one-loop version of the screened massive expansion, with light quarks modeled as fields of fixed infrared mass, the paper argues that the QCD deconfinement boundary can be read off the temperature at which the zero-frequency longitudinal gluon propagator, at vanishing spatial momentum, reaches its maximum. At zero chemical potential that maximum sits at Tmax(0)≈77 MeV, roughly half the usually quoted crossover temperature, and the ratio of pure-glue to full-QCD critical temperatures comes out close to the Polyakov-loop ratio. As chemical potential grows below the lightest quark mass M1=350 MeV, Tmax falls monotonically to zero at μ=M1, making the infrared quark mass the critical chemical potential. Above μ=M1 the propagator ceases to be non-monotonic once quark masses drop as expected with chiral symmetry restoration. The shape of the resulting Tc(μ) curve reproduces the expected topology of the QCD phase diagram and is qualitatively insensitive to the expansion parameters.","feed_headline":"Gluon peak traces QCD phase boundary down to zero temperature","feed_subtitle":"One-loop calculation puts deconfinement at 77 MeV at zero density and links the critical chemical potential to the light quark mass.","key_machinery":"The machinery is the screened massive expansion: ordinary QCD perturbation theory reorganized by giving transverse gluons a tree-level mass and quarks fixed infrared masses, with a compensating two-gluon mass counterterm. To one loop the gluon polarization receives gluon, ghost, and quark-loop contributions, and the quark loop at finite temperature and density is evaluated by shifting the fermionic Matsubara frequencies by iμ, leaving closed expressions up to a one-dimensional momentum integral. The order parameter used here is the inverse of the longitudinal propagator at zero Matsubara frequency and vanishing spatial momentum, essentially a Debye mass; its first decrease and then growth with T defines Tmax.","core_discovery":"The central claim is that the deconfinement transition in full QCD is encoded in the non-monotonic temperature dependence of the longitudinal Landau-gauge gluon propagator at zero Matsubara frequency: the temperature Tmax at which ΔL(0,|p|→0) is maximal is the (pseudo)critical temperature Tc(μ). At one loop with effective quark masses M1=350 MeV and M2=450 MeV and gluon mass m0=656 MeV, the paper finds Tmax(0)≈0.117 m0≈77 MeV, a steady decrease with μ for μ<M1, and Tmax→0 as μ→M1. It also finds that the normalized Tmax(μ)/Tmax(0) curve is practically independent of the chosen quark and gluon masses when plotted against μ/M1, and that the two humps appearing for μ between M1 and M2 disappear if the effective quark masses drop to about half their values above the Tmax curve, as chiral symmetry restoration would suggest.","pith_inferences":["If the effective quark masses fall continuously rather than abruptly above Tc(μ), the phase boundary would bend smoothly and the endpoint would sit above μ=M1; this is testable by coupling the expansion to a gap equation for the quark mass.","Because the crossover temperature is gauge- and observable-dependent, a gauge-independent check would be to test whether the longitudinal-propagator maximum in a first-order region of the phase diagram coincides across covariant gauges; the paper leaves this as future work.","The curvature-change point near (μB,T)≈(0.48 GeV, 0.62Tc) could be compared with lattice Taylor-expansion and imaginary-chemical-potential determinations of the crossover line curvature.","The criterion could be extended to gluon spectral functions or transverse-sector observables, which would connect the propagator-maximum method to transport coefficients of the quark-gluon plasma."],"forward_implications":["At μ=0 the full-QCD deconfinement temperature is predicted at about 77 MeV, roughly half the phenomenological 155–175 MeV, and the ratio T_c^YM/T_c^QCD≈1.57 is close to the Polyakov-loop based ratio 1.54.","For μ<M1 the critical temperature decreases monotonically with baryon chemical potential, with a change in concavity near (μ,T)=(0.46M1,0.62Tc) that may signal a change in the nature of the transition.","The critical chemical potential at T=0 equals the lightest infrared quark mass M1≈350 MeV; beyond it the longitudinal propagator is a decreasing function of temperature once quark masses are reduced above the boundary.","Adding a charm-like quark of mass ~1.2 GeV leaves Tc(μ) essentially unchanged for μ<M1, while dropping the strange quark raises Tc(μ) by 4–9% in the same region."],"supporting_citations":[{"why":"Provides the thermal extension of the screened massive expansion in pure Yang-Mills, including the finite-temperature propagator expressions and the observed Tmax behavior.","marker":"[80]"},{"why":"Supplies quenched lattice data showing that the longitudinal gluon propagator peaks at Tc≈270 MeV, the empirical anchor for the Tmax criterion.","marker":"[82]"},{"why":"Contains the suggestion that the longitudinal propagator at fixed momentum may serve as an order parameter for the deconfinement transition.","marker":"[81]"},{"why":"Gives the zero-temperature fit values m0=656 MeV and π0=−0.876 used for the full-QCD parameter choices.","marker":"[43]"},{"why":"Motivates the treatment of quarks with effective infrared masses in the screened massive expansion framework.","marker":"[61]"},{"why":"Lattice quark propagator results showing infrared masses around 350–450 MeV, providing the values for M1 and M2.","marker":"[157]"},{"why":"Supplies the finite-temperature field theory frequency-sum formalism used to evaluate the quark loop at finite chemical potential.","marker":"[158]"},{"why":"Establishes the screened massive expansion of pure Yang-Mills, including the one-loop diagrams and polarization structure used here.","marker":"[40]"}],"fun_headline_variants":["Gluon peak marks deconfinement at 77 MeV","Gluon propagator peak traces QCD phase boundary","Mass-independent gluon peak maps QCD phase diagram","Zero-T limit of gluon peak sets critical chemical potential","Gluon peak predicts deconfinement down to zero temperature"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The result hinges on treating the light quark effective masses as constant inputs, M1=350 MeV and M2=450 MeV, that abruptly drop above the Tmax curve, and on assuming that the longitudinal gluon propagator's maximum marks deconfinement in full QCD exactly as it does in quenched lattice Yang-Mills.","fun_headline_variants_meta":{"raw":{"variants":["Gluon peak marks deconfinement at 77 MeV","Gluon propagator peak traces QCD phase boundary","Mass-independent gluon peak maps QCD phase diagram","Zero-T limit of gluon peak sets critical chemical potential","Gluon peak predicts deconfinement down to zero temperature"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000861,"raw_usage":{"total_tokens":3694,"prompt_tokens":861,"completion_tokens":2833,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":477,"completion_tokens_details":{"reasoning_tokens":2749}},"tokens_in":477,"tokens_out":2833,"duration_ms":20051,"temperature":1.0,"reasoning_tokens":2749,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T11:26:55.347289+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A lattice calculation of the unquenched zero-frequency longitudinal gluon propagator in 2+1-flavor QCD at vanishing momentum that finds no maximum near T≈77 MeV at μ=0, or a gap-equation computation in which Tmax(μ) does not vanish at μ=M1, would disprove the claim.","supporting_citations":[{"cited_title":"Braun, W.-J","cited_arxiv_id":null,"evidence_quote":"Lattice quark propagator results showing infrared masses around 350–450 MeV, providing the values for M1 and M2."},{"cited_title":"Gao and J","cited_arxiv_id":null,"evidence_quote":"Supplies the finite-temperature field theory frequency-sum formalism used to evaluate the quark loop at finite chemical potential."}],"review_version":1}