{"id":"660bb505-fa11-4764-9776-00dc36b7f2c2","arxiv_id":"2505.10602","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":8,"one_line_summary":"At one loop in the Refined Gribov-Zwanziger framework, the quark mass function agrees roughly with lattice QCD after fitting the gluon propagator, but the dressing function disagrees in the infrared.","lead":"Quarks are coupled minimally to the Refined Gribov-Zwanziger gluon field, and the one-loop quark propagator is computed. The quark mass function roughly matches lattice data after fitting the gluon propagator, while the quark dressing function fails in the deep infrared, as in the Curci-Ferrari model.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The 'prediction' of the quark mass function is weakened by scanning g0 and mu and selecting the pair that minimizes chi2 to the quark data; the parameters are not fully fixed before the quark comparison.","rationale":"I agree with the reader's weakest assumption. The paper is technically careful, presents analytic expressions, checks the lambda -> 0 limit against the Curci-Ferrari result, and openly reports the Z(p) failure. The main issue is epistemic rather than computational: the abstract and conclusions claim a prediction using parameters fixed by the gluon sector, but Sec. IV describes a scan over g0 and mu followed by selection of the pair that minimizes the quark mass function discrepancy. This selection, plus the use of the lattice quark mass function to set m_psi, means the agreement in Fig. 4 is not an independent test of minimal coupling. A fixed, pre-registered choice of (g0, mu) is needed; the concrete test above would settle whether the agreement survives. The reader's CONDITIONAL verdict is appropriate; I would not move to REJECT because the computation is valid, the selection is disclosed, and the framework's comparison with the Curci-Ferrari model is informative. No change to the reader's verdict is needed.","tokens_in":22509,"tokens_out":5564,"duration_ms":55242,"concrete_test":"Re-run the analysis with (g0, mu) chosen before any use of quark data: take mu = 1 GeV (the renormalization scale used for the quenched gluon fit in Ref. [87]) and determine g0 together with {lambda^2, m^2, M^2} by fitting only D(p) to the unquenched lattice data [91]. Then compute M(p) and chi2_M against [5] once, with m_psi = M_Lt(p = mu). If this chi2_M is substantially larger than the value at the selected (g0 = 7, mu = 4 GeV), the headline agreement is largely an artifact of selection; if it is comparable, the minimal-coupling prediction survives.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim (Sec. IV and conclusions) is that, after fitting the unquenched gluon propagator, the quark mass function M(p) is predicted with the same fixed parameters. The procedure, however, selects (g0, mu) using the quark data themselves: for each g0 in 0-9 (step 0.05) and each mu = 1, 2, 3, 4 GeV, one fits the RGZ masses {lambda^2, m^2, M^2} to the gluon lattice data [91], sets m_psi = M_Lt(p = mu) from the quark lattice [5], computes M(p), and then retains the pair minimizing chi2_M (Sec. IV). Thus {lambda, m, M} are fixed per pair, but the choice of the pair is a fit to the quark mass function. Moreover, the quark mass input m_psi comes from the same lattice data being 'predicted'. The selected point g0 = 7, mu = 4 GeV also gives N_c g0^2/(4 pi)^2 ~ 0.93 at p = 0, so the one-loop expansion is not manifestly under control in the deep IR. The Z(p) failure is honestly disclosed and plausibly two-loop, but the M(p) agreement cannot be called a parameter-free confirmation until the selection effect is removed.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes the one-loop quark propagator in the Refined Gribov-Zwanziger (RGZ) framework with quarks minimally coupled to the gauge sector, extending the authors' earlier one-loop gluon computation. After deriving the analytic one-loop expressions (Appendix A), the authors fit the unquenched RGZ gluon propagator to lattice data [91], treating the RGZ masses {lambda^2, m^2, M^2} as free parameters for each choice of g0 and the renormalization scale mu, with a toy running coupling defined in Eq. (24). Using the fitted parameters, they compare the quark mass function M(p) and dressing function Z(p) with lattice data [5], reporting good agreement for M(p) at the selected point g0=7, mu=4 GeV, while Z(p) fails qualitatively in the infrared. They also compare their results with the Curci-Ferrari model limit and with a non-minimal matter coupling proposal. The central claim is that the RGZ non-perturbative gluon effects, transmitted through minimal coupling at one loop, suffice to predict the quark mass function once the gluon propagator is fitted.","tokens_in":22855,"tokens_out":2481,"duration_ms":25764,"significance":"If the central claim holds, the paper provides a nontrivial demonstration that minimal coupling transmits the RGZ non-perturbative gluon dynamics to the quark sector at one loop, without invoking additional matter-sector horizon terms. The analytic computation is explicit and self-contained, with a clear and verified limit lambda -> 0 reproducing the Curci-Ferrari expressions of Ref. [92]; the appendix contains the complete one-loop integrals. The paper also honestly discloses the failure of the one-loop dressing function and attributes it to missing two-loop effects, consistent with the known Curci-Ferrari experience. However, the significance of the 'prediction' of M(p) is substantially tempered by the fitting methodology in Sec. IV: the pair (g0, mu) is selected by minimizing chi^2 to the very quark data being compared, and the quark input mass m_psi is taken from the same lattice data, so the agreement is not parameter-free. The toy running coupling of Eq. (24) is also ad hoc and, at the selected point, gives an expansion parameter approaching 0.93 at p=0, making the one-loop truncation in the deep IR questionable.","major_comments":[{"comment":"The claim that M(p) is 'predicted' after fixing parameters to the gluon data is weakened by the selection procedure. For each g0 in steps of 0.05 and each mu in {1,2,3,4} GeV, the RGZ masses are refitted to the gluon lattice data, and the pair (g0, mu) is then chosen by minimizing chi^2_M against the quark mass function data. Consequently, the final comparison is a two-parameter scan selected on the quark observable, and the mass parameters are not fully fixed before the quark sector is used. The authors should either present the full chi^2_M landscape to show that the agreement is not a sharp minimum, or determine g0 and mu from the gluon fit alone (or from another criterion independent of the quark data).","section":"Sec. IV"},{"comment":"The toy running coupling in Eq. (24) is an ad hoc prescription with the IR regulator Lambda placed at the would-be Landau pole. At the selected parameters g0=7, mu=4 GeV, the value of N_c g^2/(4 pi)^2 is 0.257 at the renormalization point but grows to roughly 0.93 as p -> 0, so the one-loop expansion is not manifestly under control in the deep infrared, where the mass function comparison is performed. The paper should quantify the sensitivity of the M(p) agreement to this running prescription, for instance by comparing with the fixed-coupling limit Lambda -> infinity and with other choices of Lambda.","section":"Sec. IV; Eq. (24)"},{"comment":"The quark mass input m_psi is set to the lattice mass function at the renormalization point, m_psi = M_Lt(p=mu), using the same lattice data set [5] that is later used to assess agreement. This means the analytical M(p) is forced to match the lattice data at one point by construction, and the comparison is not fully independent. The authors should clarify how much of the agreement in Fig. 4 is driven by this matching condition and, ideally, test the sensitivity to varying m_psi around the lattice value.","section":"Sec. III; Sec. IV"}],"minor_comments":[{"comment":"There is a typo in the sentence preceding Eq. (13): it reads 'is given by e Eq. (13)' and should be 'is given by Eq. (13)'.","section":"Sec. II, below Eq. (13)"},{"comment":"The phrase 'This is agreement with the analogue computation' is grammatically awkward; it should read 'This is in agreement with the analogous computation in the Curci-Ferrari model.'","section":"Abstract and Sec. V"},{"comment":"The chi^2_M definition normalizes by the lattice value M_Lt(p_i) rather than by the lattice statistical error; the authors should state whether the quoted agreement accounts for the uncertainties of the quark lattice data, and if not, should acknowledge that the chi^2 value is not a statistical measure.","section":"Footnote 1, Sec. IV"},{"comment":"The captions of Figs. 3-5 would benefit from stating whether the curves include the global multiplicative factors A_D and A_Z, and from specifying units for the horizontal axis (GeV) to match the text.","section":"Sec. IV, Figs. 3-5"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is a straightforward one-loop computation in an established framework and is within the scope of the journal. The main concern is not the calculation itself but the strength of the 'prediction' claim in view of the parameter selection; I would encourage the editor to ask for the additional analysis described in the major comments without treating this as a fundamental flaw. The treatment of the literature is fair and complete."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: the advertised prediction of the quark mass function is not parameter-free, but the one-loop calculation itself is honest, nontrivial, and technically believable. The paper is the first one-loop quark propagator in the RGZ framework with minimal coupling, and it includes the first unquenched one-loop gluon fit in that framework.\n\nCredit where due: the integrals in App. A are plausible, and the lambda->0 limit exactly reproduces the quark-propagator expression of the Curci-Ferrari model [92], which is a strong cross-check. The Z(p) failure is disclosed up front, and the comparison with the Curci-Ferrari two-loop story is reasonable. The authors also compare their minimal coupling with the non-minimal horizon-like matter coupling, arguing with some justice that minimal coupling may suffice.\n\nThe soft spot is exactly what the stress test says: Sec. IV scans g0 from 0 to 9 and mu = 1,2,3,4 GeV, fits {lambda,m,M} to the unquenched gluon data, then picks the pair that minimizes chi2 against the quark mass function. At that point the quark data has been used twice: once to set m_psi = M_lattice(p=mu), and once to select g0 and mu. Calling the mass function a prediction overstates it. The chi2 also neglects lattice errors. And the toy running in Eq. (24), with Lambda placed at the would-be Landau pole, is ad hoc, though the authors acknowledge this. One more check: the chosen point (g0=7, mu=4 GeV) gives N_c g0^2/(4 pi)^2 ~ 0.26 at the renormalization point, but in the deep IR, where the agreement is strongest, the expansion parameter is ~0.9. So the one-loop expansion is not manifestly under control there.\n\nThose are limitations of a defensible first step, not hidden fatal flaws. The calculation stands as a genuine one-loop result, the CF limit is a good check, and the authors are transparent about what has been fitted versus predicted.\n\nWho is this for? People working in the Gribov-Zwanziger program, the Curci-Ferrari model, or infrared QCD correlation functions. A referee should check the appendix integrals and the fitting protocol, but the paper deserves referee time. I would send it to peer review.","headline":"A genuine one-loop RGZ quark calculation with an honest Curci-Ferrari limit check, but the advertised 'prediction' of the quark mass function is softened by scanning g0 and mu against the quark data.","tokens_in":23379,"tokens_out":2363,"would_cite":false,"duration_ms":20721,"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":"Minimal quark coupling transmits Refined Gribov-Zwanziger effects at one loop, yielding a lattice-matching quark mass function.","keywords":["quark propagator","refined Gribov-Zwanziger","Landau gauge","one-loop","quark mass function","unquenched gluon propagator","minimal coupling","Gribov copies"],"falsifier":"Fix the renormalization scale and coupling before looking at quark data (for example, $\\mu = 2$ GeV and $g_0 = 4$), fit only the unquenched gluon-lattice data, and then compare the resulting one-loop $M(p)$ against lattice points without further adjustment; if agreement disappears or the mass function develops the wrong infrared sign, the claimed prediction is an artifact of the selection. A second falsifier is to compute the two-loop quark dressing function in the RGZ framework: if $Z(p)$ still fails to reproduce the infrared shape, the analogy with the CF model breaks down.","tokens_in":22295,"feed_emoji":"⚛️","tokens_out":3813,"duration_ms":36202,"temperature":0.7,"pith_summary":"This paper claims that the non-perturbative effects encoded in the Refined Gribov-Zwanziger (RGZ) gluon propagator reach the quark sector at one loop through ordinary minimal coupling: no extra horizon-like matter term is needed. The authors compute the one-loop quark propagator in Landau gauge, fit the unquenched gluon propagator to lattice data to fix the RGZ mass parameters, and then use those parameters unchanged to predict the quark mass function $M(p)$, finding good agreement with lattice results. The quark dressing function $Z(p)$ fails to reproduce lattice data in the infrared, a failure the authors expect a two-loop computation to cure by analogy with the CF model.","feed_headline":"One loop carries Gribov effects into the quark mass function","feed_subtitle":"With RGZ parameters fixed by gluon data, the predicted quark mass function matches lattice data.","key_machinery":"The load-bearing object is the tree-level RGZ gluon propagator form factor $D_0(p)=\\frac{p^2+M^2}{(p^2+M^2)(p^2+m^2)+\\lambda^4}$, whose complex poles and finite infrared value carry the information about Gribov-copy elimination and condensate formation. Feeding this propagator into the standard one-loop quark self-energy diagram, together with a renormalization scheme that pins propagators to their tree-level form at a scale $\\mu$, produces the one-loop quark mass and dressing functions; the mass parameters $\\{\\lambda, m, M\\}$ are fixed by the gluon fit, not by the quark data.","core_discovery":"In the Refined Gribov-Zwanziger framework, the one-loop quark propagator inherits non-perturbative information from the gauge sector through a single gluonic loop containing the RGZ-modified gluon propagator. With the RGZ parameters fixed by fitting the unquenched gluon propagator to lattice data, the quark mass function $M(p)$ emerges as a parameter-free prediction and agrees reasonably with the lattice, both qualitatively and quantitatively. The same computation predicts a quark dressing function $Z(p)$ that is qualitatively wrong in the deep infrared, mirroring the one-loop CF model result; the paper argues that two-loop corrections, again following the CF development, are likely to repair this.","pith_inferences":["Editorial inference: if a future two-loop computation indeed restores the infrared shape of $Z(p)$, the minimal coupling prescription would become a parameter-free bridge from gluon physics to quark physics, strengthening the case that no non-minimal matter coupling is needed.","Editorial inference: the reported agreement for $M(p)$ rests partly on a scan over the coupling and renormalization scale that selects the pair minimizing the discrepancy with quark lattice data; a stricter test would fix these quantities before inspecting the quark sector.","Editorial inference: the same one-loop machinery could be applied to different quark masses or to finite temperature to test whether the RGZ parameters remain universal across flavors and environments.","Editorial inference: a direct one-loop comparison between the minimal and non-minimal couplings would quantify what the extra mass parameters in the non-minimal proposal actually buy, and whether the horizon-like matter term is phenomenologically distinguishable."],"forward_implications":["If the central claim is correct, minimal coupling alone is sufficient to transmit infrared gluonic effects to the quark mass function at one loop, so no additional matter horizon term is required for this observable.","The same fixed RGZ parameters can be reused to predict other matter-sector correlators, such as quark-gluon vertices, without introducing new free parameters.","The one-loop failure of the quark dressing function singles out two-loop corrections as the next decisive test of the minimal coupling prescription.","The unquenched RGZ gluon propagator remains finite at zero momentum and fits the unquenched lattice data, extending the established quenched success to the case with dynamical quarks.","In the limit where the Gribov parameter vanishes, the one-loop results reduce exactly to the CF model expressions, providing a consistency check across the two frameworks."],"supporting_citations":[{"why":"Provides the one-loop RGZ gluon propagator and the renormalization scheme that the present work extends to the unquenched case.","marker":"[87]"},{"why":"Supplies the unquenched lattice gluon-propagator data used to fit the RGZ mass parameters.","marker":"[91]"},{"why":"Supplies the lattice quark-propagator data against which the predicted mass function and dressing function are compared.","marker":"[5]"},{"why":"Gives the one-loop CF model quark results that the RGZ computation must reproduce in the zero-Gribov-parameter limit and that exhibit the same $Z(p)$ shortcoming.","marker":"[92]"},{"why":"Shows in the CF model that two-loop corrections restore the quark dressing function, motivating the paper's expectation for the RGZ case.","marker":"[28]"},{"why":"Provides the analogous one-loop minimal-coupling computation for scalar matter in the RGZ framework, which the present fermionic calculation extends.","marker":"[88]"}],"fun_headline_variants":["One-loop RGZ quark propagator matches lattice mass function","Predicting quark mass function from RGZ gluon data alone","One-loop quark dressing function fails in infrared, like Curci-Ferrari","RGZ one-loop quark propagator: mass matches, dressing misses","Quark mass from RGZ gluon loop predicts lattice; dressing does not"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The claim that the quark mass function is predicted rests on transferring RGZ parameters fixed from gluon data to the quark sector without re-tuning, but in practice the paper scans over coupling and renormalization scale and selects the pair that minimizes the quark-data discrepancy.","fun_headline_variants_meta":{"raw":{"variants":["One-loop RGZ quark propagator matches lattice mass function","Predicting quark mass function from RGZ gluon data alone","One-loop quark dressing function fails in infrared, like Curci-Ferrari","RGZ one-loop quark propagator: mass matches, dressing misses","Quark mass from RGZ gluon loop predicts lattice; dressing does not"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000637,"raw_usage":{"total_tokens":2951,"prompt_tokens":979,"completion_tokens":1972,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":595,"completion_tokens_details":{"reasoning_tokens":1880}},"tokens_in":595,"tokens_out":1972,"duration_ms":13268,"temperature":1.0,"reasoning_tokens":1880,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T21:10:21.023110+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Fix the renormalization scale and coupling before looking at quark data (for example, $\\mu = 2$ GeV and $g_0 = 4$), fit only the unquenched gluon-lattice data, and then compare the resulting one-loop $M(p)$ against lattice points without further adjustment; if agreement disappears or the mass function develops the wrong infrared sign, the claimed prediction is an artifact of the selection. A second falsifier is to compute the two-loop quark dressing function in the RGZ framework: if $Z(p)$ still fails to reproduce the infrared shape, the analogy with the CF model breaks down.","supporting_citations":[{"cited_title":"Confinement and dynamical chiral symmetry breaking in a non-perturbative renormalizable quark model","cited_arxiv_id":"1303.7134","evidence_quote":"Supplies the unquenched lattice gluon-propagator data used to fit the RGZ mass parameters."}],"review_version":1}