{"id":"fe43e45c-d79d-4cbc-8931-964875898092","arxiv_id":"1908.08363","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"A RG-optimized perturbation theory calculation of the cold, dense QCD pressure at two loops reduces renormalization-scale uncertainty and matches higher-order pQCD results better than standard two-loop pQCD.","lead":"This paper applies a renormalization-group-based resummation method to compute the pressure of cold, dense quark matter at two-loop order. It finds that the resummed result is less sensitive to the arbitrary renormalization scale than standard perturbation theory and closely tracks higher-order perturbative predictions.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"NLO RGOPT pressure in this application is tied to an ad hoc renormalization-scheme parameter B2; the claimed agreement with N3LO pQCD is not shown to be scheme-independent.","rationale":"I read the central claim as: at NLO (order g) the RGOPT pressure for cold, dense, three-flavor QCD is closer to the state-of-the-art N3LO pQCD result and has smaller renormalization-scale dependence than standard pQCD at the same order. For that claim to be well founded, the NLO RGOPT pressure must be well-defined and not effectively tuned. The obstruction is real: the MOP and reduced RG equations have no real solutions at NLO, and the paper introduces B2 through Eq. (4.6), fixing it by the contact condition Eq. (4.7). That condition is a reasonable but ad hoc scheme-choice prescription. The numerical pressure depends on B2, and Fig. 3 shows that B2 g^2 is not negligible in the low-µ regime. Moreover, the alternative NLO prescription in Sec. IV.C gives a different central pressure, so the method has at least two inequivalent NLO outputs; choosing one of them is part of what produces the claimed agreement. This matches the reader's identified weakest assumption. I do not see a flaw in the derivation of the RGOPT subtraction coefficients or in the LO RG-invariance property, and the comparison to Ref. [16] is appropriate as a benchmark. The concern is about robustness of the NLO numerical result, not about internal consistency of the formalism, so the appropriate verdict remains conditional pending a scheme-independence check. I agree with the reader's assessment and recommend no change to the verdict.","tokens_in":23,"tokens_out":7147,"duration_ms":134604,"concrete_test":"Recompute the NLO RGOPT normalized pressure at M = 2µ for µ = 0.6, 1.0, 1.5, 2.0 GeV using the same Eqs. (3.13), (3.14), (4.1), and (4.7), but replace the tangency condition (4.7) with an independent RSC selection rule, such as minimizing |B2 g^2| subject to real m. If any of the four P/Pfg values shifts by more than the NLO-RGOPT-versus-NLO-pQCD difference read from Fig. 5 at the same µ, then the claimed agreement is not robust under scheme choice.","verdict_should_be":"UNCHANGED","load_bearing_attack":"At NLO the optimization equations (4.1) and (2.12) have no real solutions, so the calculation is defined only after introducing the RSC mass redefinition (4.6) and fixing B2 by the tangency condition (4.7). That condition is a mathematical 'closest to MS-scheme' criterion, not a consequence of RG invariance; nothing guarantees it selects the physical scheme. The paper's own Fig. 3 shows B2 g^2 is a nontrivial function of µ and M, growing towards small µ and becoming large at M = µ, so the NLO pressure depends on this choice. Internal evidence of prescription sensitivity is the alternative IV.C variant: Eq. (4.13) uses a perturbative solution of the RG equation without B2 and produces a different central pressure (the paper says the 'exact' B2-based pressure is 'sensibly lower' than the approximate one). Therefore the headline agreement in Figs. 4 and 5 could be an artifact of the tangency RSC rather than a robust property of RGOPT.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript applies renormalization group optimized perturbation theory (RGOPT) to the quark contribution to the QCD pressure at zero temperature and finite chemical potential, for three massless flavors, working through two-loop (order g) level. At leading order, order g^0, the method produces a nontrivial, exactly renormalization-group-invariant pressure and optimized mass. At next-to-leading order the optimization equations have no real solutions, so the authors introduce a renormalization-scheme-change parameter B2, fixed by the tangency condition in Eq. (4.7). The resulting NLO pressure is compared with perturbative QCD results of Refs. [13] and [16], including the partial N3LO alpha_s^3 ln^2 alpha_s contribution, and is claimed to be in much better agreement with those higher-order results than standard pQCD at the same order, while also showing reduced renormalization-scale sensitivity. An alternative, simpler NLO prescription based on a perturbative expansion of the optimized mass is also explored and is reported to agree with Eq. (4.11) to better than 1.5%. The paper is the first RGOPT application to in-medium QCD at nonzero chemical potential and is aimed at the cold-dense regime relevant to neutron-star equations of state, where lattice QCD encounters the sign problem.","tokens_in":18648,"tokens_out":27266,"duration_ms":243499,"significance":"If the central claim is robust, the paper provides a useful resummation framework for cold and dense QCD, a regime of direct relevance to neutron-star physics and currently inaccessible to lattice simulations. The manuscript has real strengths: the LO result is exactly RG invariant, the calculation is presented in explicit analytic form, the comparison with known perturbative results and with the large-N limit of the GN model is informative, and the limitations of the approach (residual scale dependence, loss of reliability at low mu) are discussed candidly. The central comparison is not circular in the narrow sense: the variational mass and B2 are not fixed by matching the N3LO pressure. However, the NLO claim rests on an extra scheme-change parameter whose physical uniqueness and scheme independence are not established, and the paper presents two NLO variants whose numerical outputs differ. These features make the central quantitative claim currently conditional rather than fully demonstrated.","major_comments":[{"comment":"The NLO RGOPT pressure, which is the central result of the paper, is defined only after introducing the renormalization-scheme-change parameter B2, fixed by the tangency condition (4.7), because the optimization equations otherwise have no real solutions. The manuscript does not establish that this prescription is unique, nor does it test whether the quantitative agreement with Eq. (4.11) survives within a perturbatively reasonable range of B2 or under a different RSC choice. Since Fig. 3 shows that B2 g^2 is a non-negligible function of mu and M, the close agreement with the higher-order pQCD result could in principle be an artifact of the contact condition. A sensitivity study is needed before the headline claim can be considered robust.","section":"Section IV.B, Eqs. (4.6)-(4.7), Fig. 3"},{"comment":"There is an unresolved tension between the two NLO RGOPT variants. The text states that the exact NLO pressure is \"sensibly lower\" than the other approximations, while the alternative prescription based on the perturbative mass of Eq. (4.13) agrees with Eq. (4.11) to better than 1.5%. The abstract and conclusions present \"the NLO RGOPT pressure\" as being in much better agreement with the higher-order perturbative result without distinguishing these versions. The authors should specify which curve is the claimed RGOPT prediction and quantify the difference between the two NLO pressures; otherwise the central claim is ambiguous.","section":"Section IV.C, Figs. 4 and 8"}],"minor_comments":[{"comment":"Please clarify the factor of Nc in the stated replacement below Eq. (4.4): with s1 as defined in Eq. (3.11) and Nc=3, the expression -1/2 - 8 pi^2 s1 is not equal to 11/84; the value 11/84 corresponds to -1/2 - (8 pi^2/Nc) s1. The displayed formula or the definition of s1 should be adjusted accordingly.","section":"Section IV.A, Eq. (4.4)"},{"comment":"The comparison is with the partial N3LO result of Ref. [16], which contains only the alpha_s^3 ln^2 alpha_s contribution, not the complete N3LO pressure. This qualification should be made explicit in the abstract and conclusions, where \"higher-order perturbative results\" could be read as implying a complete next-order calculation.","section":"Abstract and Section V"},{"comment":"The labels \"pQCD O(g^3)\" should indicate that this is the leading-logarithm contribution at N3LO from Ref. [16], in order to avoid suggesting that the full N3LO pressure has been computed.","section":"Captions of Figs. 4 and 8"},{"comment":"There is a typo in the Introduction: \"Beam Energy Scam\" should be \"Beam Energy Scan.\"","section":"Section I"},{"comment":"Equation (4.3) is first derived as the LO formal solution, but Section IV.B refers to solving \"the MOP equation (4.3) at two-loop order\". Please clarify that in the NLO case one solves the analogous numerical derivative condition, not the literal LO closed form.","section":"Section IV.B"}],"recommendation":"major_revision","confidential_remarks":"The main concern for the editor is the prescription dependence of the NLO result. If the authors can provide a robustness test for B2 and clarify the relationship between the exact and approximate NLO variants, the paper would be considerably stronger. The factor-of-Nc issue in Eq. (4.4) should also be corrected."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"First, the thing to know: this is the first RGOPT application to QCD at T=0 with finite chemical potential, and the main numerical claim is credible. At NLO (two-loop) the optimized pressure sits much closer to the N3LO pQCD result of Gorda et al. (the alpha_s^3 ln^2 alpha_s term) than ordinary pQCD at the same order, and the scale band M=mu..4mu is moderately narrower. The calculation is not circular: the N3LO pressure is used only for comparison, not to fix any parameter; the input alpha_s(M0=1.5 GeV)=0.326 is taken from Ref. [16]. I think the authors also deserve credit for being explicit that the LO result is a poor approximation and that the NLO scale improvement is only moderate.\n\nWhat is genuinely new: RGOPT has been applied before to QCD at mu=0 and to other models at finite temperature, but not to cold dense QCD. The LO piece is RG invariant by construction and gives a non-trivial pressure at order alpha_s^0, which is a real feature of the method. The NLO part uses the known two-loop in-medium pressure and the optimization is done carefully.\n\nSoft spots. The biggest one is the B2 parameter. At NLO the MOP and RG equations have no real solutions, so the authors introduce a scheme change m -> m(1+B2 g^2) and fix B2 by the tangency condition (4.7). That is a well-defined mathematical prescription, but it is not independently physical, and Fig. 3 shows B2 g^2 is not tiny at low mu. The paper's own alternative variant, Eq. (4.13), which avoids B2 by using a perturbative solution of the RG equation, gives a central pressure that agrees even more closely with N3LO pQCD but lies at a different value from the B2-based curve. So the headline agreement is not purely an artifact of B2, but the numerical NLO result is prescription-dependent at the few-percent level. A referee should ask for a scan over B2 or a cleaner argument that the tangency condition is the right one.\n\nThere is an apparent typo in Eq. (4.13): the final coefficient 1 - sqrt(43/8) is negative, while the first expression evaluated at pF ~ mu is positive (1 - sqrt(43)/8). The authors clearly used the positive version in their numerics, but the printed formula needs fixing. Minor point: including s1 at LO is a reasonable but slightly ad hoc way to improve the LO curve; it does not affect the NLO conclusion.\n\nWho should read it: people who build neutron-star EoS from pQCD and people who care about resummation methods at finite density. It is not a final answer but a first step. I would send it to peer review with a referee who knows the RGOPT literature and ask for the B2 clarification and the typo fix. If those come back clean, I would cite it.","headline":"First RGOPT calculation for cold, dense QCD gives a plausible NLO pressure that tracks N3LO pQCD and shrinks scale dependence, but the NLO curve leans on an extra renormalization-scheme parameter fixed by a tangency condition, so treat the headline agreement as promising rather than proven.","tokens_in":19203,"tokens_out":5538,"would_cite":true,"duration_ms":51121,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper claims that a variational resummation, RGOPT, already at two-loop order gives a cold dense QCD pressure that matches higher-order perturbative results and reduces renormalization-scale sensitivity.","keywords":["renormalization group optimized perturbation theory","QCD pressure","cold dense QCD","finite chemical potential","resummation","equation of state","neutron stars","renormalization scale dependence"],"falsifier":"Compute the RGOPT pressure at the next perturbative order, $O(g^2)$, at $T=0$ and finite $\\mu$ with the same $B_2$ contact prescription; if the optimized pressure no longer tracks the N3LO pQCD curve and the scale-dependence band widens instead of shrinking, the NLO agreement was an artifact of the scheme-choice prescription rather than a genuine resummation property.","tokens_in":18201,"feed_emoji":"⚛️","tokens_out":9499,"duration_ms":81376,"temperature":0.7,"pith_summary":"The paper tries to establish that renormalization group optimized perturbation theory (RGOPT) gives a usable first-principles approximation to the pressure of cold, dense QCD with three massless quark flavors, a regime where lattice QCD cannot currently operate. At leading order the method already produces a non-perturbative pressure that is exactly invariant under changes of the renormalization scale, whereas ordinary perturbation theory at that order is just the free-gas result. At next-to-leading order the RGOPT pressure is claimed to lie much closer to the known higher-order perturbative result, which includes an $\\alpha_s^3\\ln^2\\alpha_s$ term, than standard perturbative QCD at the same two-loop order does. The same calculation also narrows the renormalization-scale uncertainty compared with same-order pQCD. If correct, RGOPT supplies an alternative equation of state for neutron-star matter at high baryonic density.","feed_headline":"Resummed QCD pressure tracks higher-order pQCD at cold density","feed_subtitle":"A variational resummation cuts scale uncertainty in dense quark matter equations of state.","key_machinery":"The central machinery is the RGOPT interpolation: deform the QCD Lagrangian by rescaling the coupling $g\\to\\delta g$ and adding a variational mass term $m(1-\\delta)^a\\bar\\psi_f\\psi_f$, with the interpolation exponent fixed by the reduced renormalization-group equation to $a=\\gamma_0/(2b_0)$. After subtracting zero-point terms $m^4\\sum_k s_k g^{k-1}$ to restore perturbative RG invariance, the arbitrary mass $m$ is fixed by the principle of minimal sensitivity, i.e. stationarity of the pressure with respect to $m$. Because the NLO stationarity and RG equations have no real solutions, the paper adds a renormalization-scheme-change parameter $B_2$ through $m\\to m'(1+B_2 g^2)$ and fixes it by requiring tangency of the two optimization curves. The resulting dressed mass acts as an in-medium variational parameter, not a physical mass, and it carries the resummation into the pressure.","core_discovery":"The central claim is that the next-to-leading order ($O(g)$) RGOPT pressure, evaluated at $T=0$ and finite quark chemical potential $\\mu$, is in much better agreement with the higher-order perturbative predictions at order $g^2$ and $g^3\\ln^2g$ than standard pQCD at the same order. The leading-order ($O(g^0)$) pressure, built from the one-loop term plus a renormalization-scheme subtraction, is already non-trivial and exactly renormalization-group invariant. At NLO the optimization equations admit real solutions only after introducing a renormalization-scheme-change parameter $B_2$ fixed by a contact condition; with that prescription the optimized pressure lies closer to the higher-order pQCD curves than NLO pQCD does at the central scale $M=2\\mu$, and its residual scale variation over $M=\\mu$ to $M=4\\mu$ is moderately smaller, especially above $\\mu\\simeq1$ GeV. The paper also gives a simpler perturbative variant whose pressure agrees with the N3LO expression to better than 1.5% for $\\mu > 0.6$ GeV at the central scale.","pith_inferences":["The growth of $|B_2 g^2|$ toward small $\\mu$ could serve as an internal diagnostic for where the resummed expansion breaks down, an implication the paper does not develop.","If an $O(g^2)$ RGOPT evaluation keeps the pressure close to the N3LO pQCD curve and further narrows the scale band, that would support the idea that the variational mass is the natural expansion variable for cold dense QCD.","Extending the calculation to beta-equilibrated massive quarks and comparing the resulting equation of state against astrophysical constraints on neutron stars would test whether this resummation is reliable where pQCD and lattice QCD are weakest."],"forward_implications":["Already at one loop, the RGOPT pressure is non-trivial and exactly renormalization-scale invariant, where ordinary pQCD at the same order is just the free-gas pressure.","At two loops, the RGOPT pressure and quark number density show a moderately reduced scale uncertainty compared with same-order pQCD in the perturbative region above roughly 1 GeV, with about a 25% improvement near $\\mu\\simeq2$ GeV.","The NLO RGOPT pressure lies closer to the $O(g^2)$ and $O(g^3\\ln^2g)$ pQCD results than NLO pQCD does, which the paper reads as evidence that the resummation captures part of the higher-order physics.","The method needs no lattice input, so it can be extended to massive quarks and used to build equations of state for neutron-star matter in the density range blocked for lattice QCD by the sign problem.","Residual scale dependence is expected to shrink further at NNLO, since RGOPT preserves perturbative RG invariance up to the next order."],"supporting_citations":[{"why":"Supplies the state-of-the-art N3LO perturbative result, including the $\\alpha_s^3\\ln^2\\alpha_s$ term, that the RGOPT pressure is compared with.","marker":"[16]"},{"why":"Provides the NNLO pQCD in-medium pressure and the original renormalization-scale dependence analysis that the paper uses as a baseline.","marker":"[13]"},{"why":"Gives the finite-density quark contributions combined with vacuum terms to build the perturbative pressure before optimization.","marker":"[12]"},{"why":"Establishes the RGOPT construction, including the renormalization-scheme-change parameter and the contact condition used to fix $B_2$.","marker":"[32]"},{"why":"Supplies the vacuum contributions and subtraction coefficients $s_k$ used at one- and two-loop order in the QCD pressure.","marker":"[34]"},{"why":"Provides the world-average value of $\\alpha_s$ used to set the running coupling input at the reference scale.","marker":"[33]"}],"fun_headline_variants":["Optimized QCD pressure cuts scale uncertainty at high density","Renormalization-group improved pressure tracks higher-order pQCD","RGOPT pressure rivals lattice at high baryon density","New resummation yields accurate dense QCD pressure"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that at next-to-leading order the missing real solutions can be repaired by a single renormalization-scheme-change parameter $B_2$ fixed by the contact condition; if that scheme choice is not the physically correct one, the close agreement with higher-order pQCD could be partly an artifact of the prescription.","fun_headline_variants_meta":{"raw":{"variants":["Optimized QCD pressure cuts scale uncertainty at high density","Renormalization-group improved pressure tracks higher-order pQCD","RGOPT pressure rivals lattice at high baryon density","New resummation yields accurate dense QCD pressure"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000661,"raw_usage":{"total_tokens":3041,"prompt_tokens":985,"completion_tokens":2056,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":601,"completion_tokens_details":{"reasoning_tokens":1990}},"tokens_in":601,"tokens_out":2056,"duration_ms":15553,"temperature":1.0,"reasoning_tokens":1990,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T11:40:39.061484+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the RGOPT pressure at the next perturbative order, $O(g^2)$, at $T=0$ and finite $\\mu$ with the same $B_2$ contact prescription; if the optimized pressure no longer tracks the N3LO pQCD curve and the scale-dependence band widens instead of shrinking, the NLO agreement was an artifact of the scheme-choice prescription rather than a genuine resummation property.","supporting_citations":[{"cited_title":"Gorda, A","cited_arxiv_id":null,"evidence_quote":"Supplies the state-of-the-art N3LO perturbative result, including the $\\alpha_s^3\\ln^2\\alpha_s$ term, that the RGOPT pressure is compared with."},{"cited_title":"Kneur and A","cited_arxiv_id":null,"evidence_quote":"Establishes the RGOPT construction, including the renormalization-scheme-change parameter and the contact condition used to fix $B_2$."},{"cited_title":"Kneur and A","cited_arxiv_id":null,"evidence_quote":"Supplies the vacuum contributions and subtraction coefficients $s_k$ used at one- and two-loop order in the QCD pressure."}],"review_version":1}