{"id":"dd753ce9-d268-4053-98d7-4d8f226ad7e6","arxiv_id":"2610.23675","paper_version":1,"verdict":"ACCEPT","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"The previously unknown hard-scale four-loop coefficient β̃_E1 = 6236 ± 40 is determined, completing the full O(g^6) pressure of hot SU(3) Yang-Mills theory.","lead":"This paper computes the complete order g^6 term in the weak-coupling expansion of the pressure of hot SU(3) Yang-Mills theory, a problem open for over two decades. The result agrees with lattice simulations and provides the first perturbative description that captures all three momentum scales of the quark-gluon plasma.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek/deepseek-v4-flash","headline":"Cancellation of unknown O(epsilon) coefficients from three-loop integrals B and S is asserted but not demonstrated; if it fails, the finite four-loop result is incomplete.","rationale":"The reader's weakest assumption was the numerical integration and subtraction procedure, but the paper provides consistency checks (pole cancellation, gauge independence, lattice comparison) that support that. A more subtle load-bearing concern is the unverified cancellation of the unknown O(epsilon) coefficients b and s. If that cancellation is not exact, the central numerical result could be unreliable. The paper acknowledges that the full cancellation is not shown here but states that it occurs; given the novelty and precision claimed, this should be demonstrated before acceptance. No other critical flaws are found: the algorithmic description is credible, the dependence on previous results is properly cited, and the lattice agreement is a strong external test. The recommended verdict is CONDITIONAL because the missing verification can be supplied (likely in the forthcoming paper), and the result is probably correct, but the current submission does not fully substantiate it.","tokens_in":55063,"tokens_out":3162,"duration_ms":27590,"concrete_test":"Independently derive the full O(g^6) pressure expression p_E by combining the four-loop hard result (including b and s symbolically) with the renormalization of the three-loop pressure and the EQCD/MQCD contributions. Verify analytically that the coefficients of b and s sum to zero. This can be done using the known renormalization constants and the three-loop pressure from Ref. [27] or [13]. If the cancellation is confirmed, the result stands; if not, the finite coefficient must be corrected.","verdict_should_be":"CONDITIONAL","load_bearing_attack":"The finite coefficients C_n in Table I explicitly depend on the unknown O(epsilon) coefficients b and s of the three-loop basketball B and spectacle S sum-integrals (Eq. (14)). The paper states that 'it turns out to be sufficient to keep b and s in symbolic form, as they cancel in the full pressure after accounting for the respective terms generated by the renormalization of the 3-loop hard pressure p_3' (End Matter, paragraph after Eq. (14)). However, no explicit verification of this cancellation is provided in the manuscript, and the referenced forthcoming paper is not yet available. If the cancellation fails or relies on an unstated assumption about the renormalization of p_3, the central value beta_E1 = 6236 +/- 40 would be contaminated by the unknown b and s. This is a load-bearing step because the final result's finite part is obtained only after this cancellation, and the paper does not supply the algebraic check needed to confirm it.","agreement_with_reader":"disagree"},"referee_report":{"model":"deepseek/deepseek-v4-flash","summary":"This paper presents the complete O(g^6) term in the weak-coupling expansion of the pressure of hot SU(3) Yang-Mills theory. The authors evaluate nine previously unknown nonfactorizing four-loop sum-integrals using the thermal loop-tree duality method with local UV and IR subtractions, thereby completing the hard-scale contribution (p_E) for general SU(N_c). Combining this with known chromoelectric (EQCD) and chromomagnetic (MQCD) contributions, they obtain the full O(g^6) pressure for SU(3) and demonstrate the analytic cancellation of infrared divergences between the different scales. The main result is the four-loop hard-scale coefficient β̃_E1 = 6236 ± 40. The resulting pressure shows quantitative agreement with 4D lattice data over a broad temperature range after EQCD resummation, providing a stringent test of the dimensional reduction paradigm.","tokens_in":55393,"tokens_out":2493,"duration_ms":24830,"significance":"If the result holds, this constitutes a landmark achievement in thermal field theory, resolving a problem that had been open for over two decades. It is the first complete O(g^6) calculation in a gauge theory and provides a direct, parameter-free test of the dimensional-reduction effective-field-theory framework at the order where all three momentum scales (hard, soft, ultrasoft) contribute. The agreement with lattice data is striking and supports the expectation that the weak-coupling expansion becomes quantitatively reliable once contributions from all scales are included. The work also demonstrates the power of the thermal loop-tree duality method for multiloop sum-integrals.","major_comments":[{"comment":"The finite coefficients C_n for several of the ten nonfactorizing sum-integrals depend on the unknown O(ϵ) coefficients b and s of the three-loop basketball B and spectacle S sum-integrals (Eq. (14)). The paper asserts that 'it turns out to be sufficient to keep b and s in symbolic form, as they cancel in the full pressure after accounting for the respective terms generated by the renormalization of the 3-loop hard pressure p_3'. However, no explicit verification of this cancellation is provided in the manuscript, and the referenced forthcoming paper is not available. This step is load-bearing because the finite part of the hard-scale contribution (and hence β̃_E1 = 6236 ± 40) is obtained only after this cancellation. The authors should include a demonstration of the cancellation (or a reference to a published source that contains it) before the central result can be considered fully established.","section":"End Matter, Table I and surrounding text","authors_acknowledge":null,"authors_acknowledge_p":null}],"minor_comments":[{"comment":"The manuscript appears to contain extensive duplicated text, with the same paragraphs (including Figure 3 and its caption) repeated many times. This is clearly a compilation error and should be corrected to make the paper readable.","section":"Pages 5-6","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The abstract states 'nine previously unknown nonfactorizing four-loop sum-integrals', while the main text notes that ten such integrals arise but I_6 is already known. This is consistent, but the phrasing could be clarified to avoid confusion.","section":"Abstract and Introduction","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The uncertainty in the magnetic constant B_G (250) dominates the total uncertainty in q_c. The authors might comment on whether improved lattice determinations of B_G are expected in the near future, as this would sharpen the comparison with lattice data.","section":"Equation (10) and Figure 3","authors_acknowledge":null,"authors_acknowledge_p":null},{"comment":"The definitions of Z_n ≡ ζ'(-n)/ζ(-n) should be spelled out explicitly, as the notation may be unfamiliar to some readers.","section":"End Matter, Table I","authors_acknowledge":null,"authors_acknowledge_p":null}],"recommendation":"major_revision","confidential_remarks":"The calculation appears to be technically sound and the results are highly significant. The one load-bearing issue is the undocumented cancellation of the unknown O(ϵ) coefficients b and s. I consider this a major but fixable gap: the authors should include the algebraic cancellation check in the revision. The text duplication on pages 5-6 is clearly a production error and should be corrected. If the cancellation is verified, the paper deserves rapid publication as a major result in thermal field theory."},"author_rebuttal":null,"desk_editor":{"model":"deepseek/deepseek-v4-flash","letter":"Hi [Colleague],\n\nThe short version: this paper delivers the complete O(g^6) pressure of hot SU(3) Yang-Mills theory, solving a problem that has been open for over twenty years. The result is genuinely new and passes multiple nontrivial checks. I think it is solid and deserves a serious referee, though there is one detail a referee should ask to see.\n\nWhat is new: the hard-scale coefficient β̃_E1 = 6236 ± 40 is obtained by evaluating nine previously unknown four-loop sum-integrals. The methodology—thermal loop-tree duality with local UV and IR subtractions—is sensible and builds on recent advances. The paper includes several strong consistency checks: the gauge parameter ξ cancels in the hard-scale contribution, the pole coefficients agree with earlier effective-theory predictions, and the I6 sum-integral is independently verified against a prior scalar-field result. Combining the hard-scale result with known EQCD and MQCD contributions yields a pressure that agrees quantitatively with four-dimensional lattice data over a broad temperature range when EQCD resummation is applied. That agreement is real and supports the dimensional-reduction paradigm.\n\nThe soft spots: the numerical evaluation uses Monte Carlo integration with proprietary software (Symbolica), and the code and raw data are not provided. This limits immediate reproducibility, though the methodological description is detailed enough for experts to follow. The stress-test concern about the cancellation of unknown O(epsilon) coefficients b and s from the three-loop basketball and spectacle integrals is real but, in my reading, less worrying than it first sounds. The paper states that these cancel after including renormalization terms from the three-loop hard pressure, and it cites published references [28, 69] for the structure of that renormalization. The claim is plausible and consistent with the overall renormalization scheme. Still, the explicit algebraic check is not shown, and a referee should ask for it or for a clear statement that it will appear in the forthcoming longer paper. I do not consider this a load-bearing flaw; the paper would still be a major result even if the numerical part is correct and the cancellation holds by construction.\n\nBottom line: this is an important paper for thermal field theory and lattice QCD. The central result—the full O(g^6) pressure—is a milestone. I would recommend sending it to peer review, not a desk reject. A referee should press for the explicit cancellation check and for more details on the numerical uncertainties, but neither issue is fatal.\n\nBest,\n[Your Name]","headline":"The paper delivers the long-awaited complete O(g^6) pressure of hot Yang-Mills theory, a genuinely new result that passes multiple consistency checks; a concern about cancellation of unknown O(epsilon) coefficients is real but minor.","tokens_in":55848,"tokens_out":2447,"would_cite":false,"duration_ms":25626,"reading_group":"yes","serious_thinker":"yes","would_accept_peer_review":true},"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek/deepseek-v4-flash","headline":"The full order $g^6$ term in the weak-coupling expansion of the pressure of hot SU(3) Yang-Mills theory has been determined, resolving a two-decade-old problem.","keywords":["thermal field theory","Yang-Mills theory","pressure","perturbation theory","four-loop","sum-integrals","loop-tree duality","dimensional reduction"],"falsifier":"A future recalculation of the nine new sum-integrals using an independent method (e.g., coordinate-space or sector-decomposition techniques) that yields a value for $\\tilde{\\beta}_{\\rm E1}$ outside the quoted uncertainty range $6236 \\pm 40$ would challenge the result.","tokens_in":56611,"feed_emoji":"🔥","tokens_out":3126,"duration_ms":31103,"temperature":0.7,"pith_summary":"The paper claims to have determined the previously missing hard-scale contribution to the order $g^6$ term in the weak-coupling expansion of the pressure of hot SU(3) Yang-Mills theory, completing a calculation that has been open for more than two decades. The authors evaluate nine unknown four-loop thermal sum-integrals using the thermal loop-tree duality method with novel local ultraviolet and infrared subtractions. Combining their result with known chromoelectric and chromomagnetic contributions, they obtain the complete $O(g^6)$ pressure. The result shows quantitative agreement with four-dimensional lattice simulations over a broad range of temperatures, providing a stringent test of the dimensional reduction effective-field-theory paradigm.","feed_headline":"Hot Yang-Mills pressure complete at order g^6","feed_subtitle":"New four-loop coefficient β̃_E1 = 6236±40 agrees with lattice data and tests dimensional reduction.","key_machinery":"The key machinery is the thermal loop-tree duality method extended to four-loop thermal sum-integrals, combined with a local subtraction procedure for ultraviolet and infrared divergences. The infrared counterterms are constructed by identifying soft subgraphs associated with vanishing Matsubara frequencies. After local subtractions, the finite remainders of ten nonfactorizing four-loop sum-integrals are evaluated numerically using Monte Carlo integration. The counterterms are evaluated analytically using dimensional regularization and integration-by-parts reduction.","core_discovery":"The core discovery is that the full order $g^6$ pressure of hot SU(3) Yang-Mills theory can now be expressed analytically, with the previously unknown hard-scale coefficient $\\tilde{\\beta}_{\\rm E1}$ determined to be $6236 \\pm 40$. This coefficient completes the weak-coupling expansion up to $O(g^6)$ and, when combined with EQCD resummation, yields a perturbative prediction that agrees with lattice data down to temperatures only a few times the critical temperature $T_c \\approx 260$ MeV. The work demonstrates that the infrared divergences between the hard, soft, and ultrasoft sectors cancel as predicted by dimensional reduction, confirming the effective-field-theory paradigm at the first order where all three momentum scales contribute.","pith_inferences":["If the thermal loop-tree duality method is as robust as demonstrated, it could be used to compute other four-loop thermodynamic quantities such as the entropy density or trace anomaly, providing further constraints on the equation of state.","The agreement with lattice data may suggest that the $O(g^7)$ contribution is small and the expansion is converging; a five-loop calculation could test this directly.","The result may help refine effective models for the quark-gluon plasma used in heavy-ion collision phenomenology by providing a precise perturbative baseline."],"forward_implications":["The complete $O(g^6)$ pressure provides a benchmark for lattice QCD thermodynamics and can be used to extract the nonperturbative magnetic coefficient with improved precision.","The successful application of thermal loop-tree duality at four loops opens the door to computing higher-order perturbative corrections, such as the $O(g^7)$ term.","The agreement with lattice data supports the conjecture that the weak-coupling expansion becomes quantitatively reliable once contributions from all three momentum scales are included.","The method can be extended to full QCD with dynamical quarks and nonzero chemical potentials, relevant for neutron-star merger phenomenology.","The independent confirmation of the lattice-based estimate $\\tilde{\\beta}_{\\rm E1} \\approx 6000 \\pm 600$ demonstrates the maturity of both perturbative and lattice approaches to finite-temperature gauge theories."],"supporting_citations":[{"why":"Reduces the 65 four-loop vacuum diagrams to 21 canonical sum-integrals, providing the starting point for the present evaluation.","marker":"[28]"},{"why":"Introduces the thermal loop-tree duality method extended to nonzero temperature, which is used to evaluate the nonfactorizing sum-integrals.","marker":"[40]"},{"why":"Provides the four-dimensional lattice data for the SU(3) Yang-Mills pressure against which the perturbative result is compared.","marker":"[8]"},{"why":"Summarizes the weak-coupling expansion of the hot QCD pressure and gives the earlier lattice-based estimate of $\\tilde{\\beta}_{\\rm E1} \\approx 6000 \\pm 600$ that the present work confirms.","marker":"[27]"},{"why":"Extracts the $g^6 \\ln(1/g)$ term in the pressure within the dimensional-reduction framework and defines the notation used for the hard-scale coefficient.","marker":"[23]"},{"why":"Determine the nonperturbative magnetic coefficients from three-dimensional lattice simulations that enter the complete $O(g^6)$ pressure.","marker":"[24, 25]"},{"why":"Provides the four-loop pressure of massless $O(N)$ scalar field theory, including the known result for one of the ten sum-integrals ($I_6$) that is validated in the present work.","marker":"[47]"}],"fun_headline_variants":["Full g^6 pressure of hot Yang-Mills now known","Two-decade gap in hot Yang-Mills pressure filled at g^6","Lattice-verified g^6 pressure for SU(3) Yang-Mills","Complete g^6 term resolves hot Yang-Mills pressure","g^6 pressure of hot Yang-Mills: completed after 20 years"],"cache_read_input_tokens":0,"weakest_assumption_plain":"The numerical evaluation of the ten nonfactorizing four-loop sum-integrals relies on a local infrared subtraction procedure whose validity for thermal soft regions is assumed but not rigorously proven in the text.","fun_headline_variants_meta":{"raw":{"variants":["Full g^6 pressure of hot Yang-Mills now known","Two-decade gap in hot Yang-Mills pressure filled at g^6","Lattice-verified g^6 pressure for SU(3) Yang-Mills","Complete g^6 term resolves hot Yang-Mills pressure","g^6 pressure of hot Yang-Mills: completed after 20 years"]},"model":"deepseek/deepseek-v4-flash","effort":"low","cost_usd":7.2e-05,"raw_usage":{"cost":7.2115176e-05,"is_byok":false,"cost_details":{"upstream_inference_cost":7.2115176e-05,"upstream_inference_prompt_cost":2.0822256e-05,"upstream_inference_completions_cost":5.129292e-05},"total_tokens":1667,"prompt_tokens":860,"completion_tokens":807,"prompt_tokens_details":{"audio_tokens":0,"video_tokens":0,"cached_tokens":256,"cache_write_tokens":0},"completion_tokens_details":{"audio_tokens":0,"image_tokens":0,"reasoning_tokens":711}},"tokens_in":860,"tokens_out":807,"cost_source":"openrouter_usage","duration_ms":8627,"image_pages":[],"temperature":1.0,"billed_attempts":1,"reasoning_tokens":711,"transport_provider":"openrouter","cache_read_input_tokens":0,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-10-19T13:33:02.672620+00:00","model_set":{"reader":"deepseek/deepseek-v4-flash"},"falsifier":"A future recalculation of the nine new sum-integrals using an independent method (e.g., coordinate-space or sector-decomposition techniques) that yields a value for $\\tilde{\\beta}_{\\rm E1}$ outside the quoted uncertainty range $6236 \\pm 40$ would challenge the result.","supporting_citations":[{"cited_title":"The $g^6$ pressure of hot Yang-Mills theory: Canonical form of the integrand","cited_arxiv_id":"2408.15830","evidence_quote":"Reduces the 65 four-loop vacuum diagrams to 21 canonical sum-integrals, providing the starting point for the present evaluation."},{"cited_title":"Cosmological phase transitions without high-temperature expansions","cited_arxiv_id":"2507.07014","evidence_quote":"Introduces the thermal loop-tree duality method extended to nonzero temperature, which is used to evaluate the nonfactorizing sum-integrals."},{"cited_title":"Precision SU(3) lattice thermodynamics for a large temperature range","cited_arxiv_id":"1204.6184","evidence_quote":"Provides the four-dimensional lattice data for the SU(3) Yang-Mills pressure against which the perturbative result is compared."},{"cited_title":"Schr¨ oder, Weak-coupling expansion of the hot QCD pressure, PoSJHW2005, 029 (2006), arXiv:hep- ph/0605057","cited_arxiv_id":null,"evidence_quote":"Summarizes the weak-coupling expansion of the hot QCD pressure and gives the earlier lattice-based estimate of $\\tilde{\\beta}_{\\rm E1} \\approx 6000 \\pm 600$ that the present work confirms."},{"cited_title":"The pressure of hot QCD up to g^6 ln(1/g)","cited_arxiv_id":"hep-ph/0211321","evidence_quote":"Extracts the $g^6 \\ln(1/g)$ term in the pressure within the dimensional-reduction framework and defines the notation used for the hard-scale coefficient."},{"cited_title":"Four-loop pressure of massless O(N) scalar field theory","cited_arxiv_id":"hep-ph/0703307","evidence_quote":"Provides the four-loop pressure of massless $O(N)$ scalar field theory, including the known result for one of the ten sum-integrals ($I_6$) that is validated in the present work."}],"review_version":1}