{"id":"9b00fc50-a7cd-4f61-ac5b-f89802af4509","arxiv_id":"1909.00796","paper_version":2,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":3,"one_line_summary":"A fit of the b-quark pole and running masses from HERA beauty data is claimed to match the NNLO perturbative QCD relation, but the comparison uses the wrong strong-coupling scale and overstates the agreement.","lead":"This paper fits the b-quark mass from HERA data in two different renormalization schemes, pole mass and running mass, at next-to-next-to-leading order. The authors report a 99.98 percent compatibility between the two schemes and the perturbative QCD prediction, a result that would matter for global fits of the proton's quark and gluon content.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Equation (20) evaluates the NNLO mass relation at alpha_s(M_Z), but the relation is defined at alpha_s(M_b); with alpha_s(M_b) ~ 0.21 the right-hand side becomes ~4.0 GeV, not 4.38 GeV, so the claimed 99.98% compatibility is a scale artifact.","rationale":"I read the central claim in good faith: the paper performs NNLO PDF fits to HERA inclusive and beauty data, extracts M_b=4.66+/-0.14 GeV and m_b=4.40+/-0.10 GeV, and claims these are compatible at the 0.02 GeV level with the standard NNLO perturbative relation. The derivation of Eq. (17) is a standard result and is not the problem. The load-bearing step is the numerical evaluation in Eq. (20). The paper's own Sec. IV states that the alpha_s in this relation is renormalized at mu=M_b, yet Sec. VI uses alpha_s(M_Z)=0.118. Running the coupling to M_b changes alpha_s/pi from 0.0376 to ~0.067, which moves the bracket from ~0.934 to ~0.86 and changes the predicted m_b(M_b) from ~4.38 to ~4.0 GeV. The 0.02 GeV agreement is therefore an artifact of the scale choice. The reader's weakest_assumption identifies exactly this issue, and I agree with the REJECT verdict. The fitted masses and PDF uncertainty study may be of interest, but they do not rescue the headline compatibility claim.","tokens_in":11273,"tokens_out":7167,"duration_ms":72353,"concrete_test":"Recompute the right-hand side of Eq. (20) with alpha_s^{(5)}(M_b) obtained by running alpha_s(M_Z)=0.118 down to M_b=4.66 GeV in the five-flavour theory (alpha_s ~ 0.21), and with N_L=4 in the coefficient of Eq. (17). If the result is ~4.0 GeV instead of 4.38 GeV, the claimed 'excellent compatibility' is not supported by the quoted numbers.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The compatibility check in Sec. VI, Eq. (20), inserts alpha_s(M_Z)=0.118 into the NNLO pole-to-MSbar relation of Eq. (17). But Sec. IV defines the coupling in exactly this relation as alpha_s^{(N_L+N_H)}(M_b), i.e. renormalized at the pole-mass scale, and the left-hand side is m_b(M_b), the running mass at that same scale. At mu=M_b with five active flavours alpha_s is ~0.21, not 0.118. Replacing 0.118/pi in Eq. (20) by alpha_s(M_b)/pi changes the bracket from ~0.934 to ~0.86, giving a predicted m_b(M_b) of ~4.0 GeV against the fitted 4.40 GeV. The discrepancy is then ~0.4 GeV, far larger than the quoted uncertainty of +/-0.10 GeV. The value N_L=3 used in Eq. (20) is also questionable at this scale, where four massless flavours plus the b quark is the appropriate counting, although the dominant problem is the scale of alpha_s. The central claim therefore rests on evaluating a mu=M_b relation at mu=M_Z, and it fails when the relation is evaluated at its proper scale.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript determines the b-quark pole mass and MSbar running mass from HERA combined inclusive and beauty production data using xFitter at NNLO (the HBPoleMass and HBRunMass fits), and it derives a standard NNLO relation between the pole and MSbar masses (Eqs. (15)-(18)). The fits give Mb = 4.66 ± 0.14 GeV and mb = 4.40 ± 0.10 GeV. The paper then claims that these two fitted values satisfy the NNLO relation at the 0.02 GeV level, quoting a \"compatibility up to approximately 99.98%\", and it studies the impact of the b-quark mass on the gluon and valence PDFs.","tokens_in":11539,"tokens_out":8215,"duration_ms":82443,"significance":"If the central compatibility claim were correct, this would be a useful cross-check of the NNLO pole-to-MSbar mass relation using HERA beauty data. The paper has strengths: it uses standard open-source tools (xFitter, QCDNUM, APFEL), it reports fit-quality numbers and uncertainties, and the PDF-sensitivity study concerning the b-quark mass is a potentially useful byproduct. The fitted mass values are plausible and broadly consistent with the PDG range. However, the central quantitative claim is invalid: Eq. (20) evaluates a relation defined at the b-quark mass scale with alpha_s(M_Z), and a correct evaluation gives a discrepancy of roughly 0.4 GeV rather than 0.02 GeV. The advertised compatibility is therefore not established.","major_comments":[{"comment":"Equation (20) evaluates the NNLO mass relation of Eq. (17) with alpha_s(M_Z) = 0.118, but Eqs. (15)-(17) define alpha_s as alpha_s^{(N_L+N_H)}(M_b), i.e. the strong coupling renormalized at the b-quark pole-mass scale. At this scale alpha_s is approximately 0.21, not 0.118. Replacing 0.118/pi in Eq. (20) by alpha_s(M_b)/pi changes the bracket from roughly 0.93 to roughly 0.86, so the theoretical prediction for m_b(M_b) from M_b = 4.66 GeV becomes approximately 4.0 GeV. The difference from the fitted m_b = 4.40 ± 0.10 GeV is then about 0.4 GeV, far outside the quoted uncertainty. The claimed compatibility therefore rests on evaluating a mu = M_b relation at mu = M_Z and fails when the relation is evaluated at its proper scale.","section":"Sec. IV and Sec. VI, Eq. (20)"},{"comment":"The text states a \"precision of 1 part in 10^2\" and a \"compatibility up to approximately 99.98%\". These statements are arithmetically inconsistent: |4.40 - 4.38| = 0.02 GeV relative to 4.40 GeV is 0.45%, not 0.02%, so the compatibility is about 99.55% even if the value 4.38 were correct. Moreover, evaluating Eq. (17) with alpha_s(M_Z) = 0.118 and N_L = 3 gives 4.35 GeV rather than 4.38 GeV, making the numerical agreement still weaker. These arithmetic errors compound the scale error identified above.","section":"Sec. VI, Eq. (20) and following text"},{"comment":"The compatibility check compares two mass values extracted from the same HERA data set with the same NNLO FONLL implementation. The HBRunMass variant (FONLL-C RUNMON) presumably uses a running-mass scheme internally, and the manuscript does not document whether that scheme conversion is independent of the perturbative relation that Eq. (20) is testing. Without such documentation, the agreement between the two fitted masses is a weaker cross-check than the text implies, even setting aside the scale error.","section":"Secs. V-VI and Table II"}],"minor_comments":[{"comment":"The sentence beginning \"Measurements of open b-quark production in DIS...\" is repeated nearly verbatim in the first two paragraphs of the Introduction.","section":"Sec. I"},{"comment":"In the bullet defining the heavy-quark scale, \"mu_r = mu_f = mu_r = sqrt(Q^2 + 4m_b^2)\" contains a typo; the last scale should presumably be a heavy-quark scale such as mu_b.","section":"Sec. V"},{"comment":"The notation for beauty contributions is inconsistent: sigma_red^{b bar b} is used in Eqs. (3)-(4), while F^{bb} appears in Eq. (5) and the surrounding text.","section":"Sec. II, Eqs. (3)-(5)"},{"comment":"The figure captions refer to the \"PPDs analysis\" where \"PPDFs analysis\" is meant.","section":"Sec. VI, Figs. 2-5"},{"comment":"The quantities delta(x d_v)/(x d_v) and similar ratios shown in the figures are not defined in the text; a definition should be given in the captions or in Sec. VI.","section":"Sec. VI, Figs. 2-5"},{"comment":"Equation (18) is presented as the inverse of Eq. (17), but the scale and flavor-counting conventions used there are not stated; the text should specify that alpha_s is evaluated at mu = m_b in that formula.","section":"Sec. IV, Eq. (18)"}],"recommendation":"reject","confidential_remarks":"For the editor: the central compatibility claim is not supported. The scale error in Eq. (20) is unambiguous and changes the conclusion from agreement to a discrepancy of roughly 0.4 GeV. In addition, the two fitted masses come from the same dataset and the manuscript does not document the internal scheme conversion of FONLL-C RUNMON, so even a corrected numerical test would be less informative than claimed. I also note that the reference list contains many self-citations; most are conference proceedings and may be defensible, but the imbalance is worth checking. My recommendation is driven by the technical errors rather than by the citation pattern."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper's central compatibility claim does not survive contact with its own equations. Eq. (17) defines alpha_s as alpha_s^(N_L+N_H)(M_b), the strong coupling at the pole mass scale. Eq. (20) plugs in alpha_s(M_Z)=0.118. That is not a negligible distinction: with alpha_s(M_b) ~ 0.21, the bracket in Eq. (20) drops from about 0.934 to 0.86, giving a predicted m_b(M_b) near 4.0 GeV against the fitted 4.40 GeV. The discrepancy is then roughly 0.4 GeV, about four times the quoted +/-0.10 uncertainty. The stress-test note is right, and the reader's strongest claim holds up.\n\nWhat the paper does well: the fits are a routine but legitimate NNLO xFitter/QCDNUM/APFEL exercise using HERA inclusive and beauty data. The mass extractions are reported with parameter tables and uncertainties, and the derivations in Sec. IV are clearly presented, even though the pole-MSbar relation is a standard result rather than new. The PDF sensitivity curves may be of casual interest.\n\nThe soft spots beyond the scale error: the 99.98% figure is arithmetically wrong regardless of scale, since |4.40-4.38|/4.38 is about 0.46%, not 0.02%. The predicted value carries no uncertainty. The choice N_L=3 is also questionable at mu=M_b, where four light flavors plus the b is the natural counting. And the 'for the first time' claim is not substantiated; pole/running compatibility checks are common in the heavy-quark literature.\n\nI would desk-reject this. The main advertised result is an artifact, and the ancillary PDF study is not enough to carry the paper. If the authors corrected the scale, propagated uncertainties, and honestly reported the resulting tension, that would be a legitimate — if far less exciting — contribution. As it stands, the paper's own equations contradict its headline claim.","headline":"The advertised NNLO compatibility check evaluates a relation defined at the b-quark scale with alpha_s(M_Z), so the claimed 99.98% agreement is a scale artifact; the underlying PDF fit is competent but the paper's central claim fails.","tokens_in":12117,"tokens_out":11570,"would_cite":false,"duration_ms":111400,"reading_group":"no","serious_thinker":"no","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["12.38.Bx","13.60.Hb","14.65.Fy"],"model":"deepseek-v4-flash","headline":"The paper claims that separately fitted b-quark pole and running masses agree with the NNLO perturbative QCD relation to within 0.02 GeV.","keywords":["b-quark pole mass","MSbar running mass","next-to-next-to-leading order","HERA beauty production","parton distribution functions","gluon distribution","QCD mass renormalization","strong coupling constant"],"falsifier":"Recompute Eq. (20) with $\\alpha_s$ evaluated at $\\mu = M_b \\approx 4.66$ GeV, where $\\alpha_s \\approx 0.21$, instead of $\\alpha_s(M_Z) = 0.118$; the right-hand side then gives about 4.0 GeV rather than 4.38 GeV, which differs from the fitted 4.40 GeV by roughly 0.4 GeV and removes the claimed 99.98% compatibility.","tokens_in":11041,"feed_emoji":"⚛️","tokens_out":11689,"duration_ms":93028,"temperature":0.7,"pith_summary":"This paper attempts to show that two standard ways of quoting the b-quark mass, the pole mass (the propagator-pole mass) and the MSbar running mass (a scale-dependent short-distance mass), satisfy the next-to-next-to-leading-order conversion relation of perturbative QCD. The authors derive the NNLO formula from the mass renormalization constants and then fit both masses from HERA inclusive and beauty-production data. The pole mass comes out $4.66 \\pm 0.14$ GeV and the running mass $4.40 \\pm 0.10$ GeV, and inserting these into the pQCD relation predicts $4.38$ GeV, an agreement to $0.02$ GeV. They also report that freeing the b-quark mass in the PDF fit improves the uncertainty band and shifts the gluon and valence distributions. A sympathetic reader would care because a genuine 0.02 GeV agreement would validate the NNLO conversion and show that beauty data can pin down a Standard Model parameter directly.","feed_headline":"Two b-quark mass fits match NNLO QCD at 0.02 GeV","feed_subtitle":"HERA beauty data yield pole and running masses that NNLO pQCD reconciles within 0.02 GeV.","key_machinery":"The carrying identity is the NNLO conversion formula derived from the ratio of the on-shell and MSbar mass renormalization factors, $Z^{\\mathrm{OS}}_m / Z^{\\overline{\\mathrm{MS}}}_m$. With the standard color factors $C_F = 4/3$, $C_A = 3$, $T_R = 1/2$ and one heavy flavor, the relation becomes $\\overline{m}_b(M_b) = M_b\\left[1 - \\frac{4}{3}\\left(\\frac{\\alpha_s}{\\pi}\\right) + \\left(\\frac{\\alpha_s}{\\pi}\\right)^2(1.0414\\,N_L - 14.3323)\\right]$. The paper inserts the phenomenological masses into this identity as a consistency check. The companion machinery is a QCD fit over HERA inclusive and beauty data with the b mass treated as an extra free parameter.","core_discovery":"On its own terms, the paper's central claim is that the separately fitted pole mass and MSbar running mass satisfy the NNLO pQCD relation: $4.40 \\simeq 4.66\\left[1 - \\frac{4}{3}\\left(\\frac{\\alpha_s}{\\pi}\\right) + \\left(\\frac{\\alpha_s}{\\pi}\\right)^2 (1.0414\\,N_L - 14.3323)\\right] = 4.38$. The paper calls this compatibility up to approximately $99.98\\%$, with a difference of $|4.40 - 4.38| = 0.02$ GeV. It further claims that the running-mass fit is preferred: the $\\chi^2$ improves by about $0.8\\%$ and the mass uncertainty shrinks by about $4\\%$ relative to the pole-mass fit, and that the gluon distribution is sensitive to the mass choice.","pith_inferences":["The authors do not state this, but if the NNLO relation is evaluated at the b-quark scale rather than at $M_Z$, the extracted pole and running masses are not compatible; the 0.02 GeV agreement depends on inserting $\\alpha_s(M_Z) = 0.118$.","A decisive extension would be a simultaneous fit of the pole mass and the strong coupling so that the scale dependence becomes an explicit fit parameter rather than an external input.","The modest $\\chi^2$ preference for the running mass may reflect PDF-mass correlations rather than a genuinely better description; separating the two requires fits at different mass scales and comparing the pulls of the beauty data."],"forward_implications":["If the compatibility is genuine, HERA beauty-production data alone can fix the b-quark mass in either scheme with an NNLO uncertainty below 0.15 GeV.","The small $\\chi^2$ improvement and reduced uncertainty for the running mass imply that future fits should treat the beauty mass as scale dependent rather than as a fixed pole mass.","The gluon distribution and valence ratios change noticeably when the beauty mass is freed, so combined inclusive-plus-beauty fits must report the mass parameter alongside the PDFs.","Agreement with the world average supports using this approach as a cross-check on the b mass."],"supporting_citations":[{"why":"Supplies the HERA I+II combined inclusive DIS data that form the main cross-section input for all three QCD fits.","marker":"[1]"},{"why":"Supplies the combined beauty vertex production data that constrain the b-quark mass.","marker":"[2]"},{"why":"Provides the heavy-flavour scheme implementation used to include beauty contributions to the structure functions at NNLO.","marker":"[11]"},{"why":"Defines the heavy-flavour matching method on which the treatment is based.","marker":"[12]"},{"why":"Provides the open-source QCD fitting program used to perform the parameter fits.","marker":"[13]"},{"why":"Provides the world average masses used as the external comparison for the extracted values.","marker":"[22]"}],"fun_headline_variants":["b-quark pole and running masses agree at NNLO within 0.02 GeV","NNLO pQCD matches b-quark mass fits to 0.02 GeV","b-quark mass: NNLO theory and data reconcile at 0.02 GeV","0.02 GeV difference: b-quark mass fits compatible at NNLO"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the strong coupling in the NNLO conversion formula can be evaluated at the Z-boson scale, $\\alpha_s(M_Z)=0.118$, even though the formula itself is defined at the b-quark mass scale; replace it with $\\alpha_s(M_b)\\approx 0.21$ and the claimed agreement disappears.","fun_headline_variants_meta":{"raw":{"variants":["b-quark pole and running masses agree at NNLO within 0.02 GeV","NNLO pQCD matches b-quark mass fits to 0.02 GeV","b-quark mass: NNLO theory and data reconcile at 0.02 GeV","0.02 GeV difference: b-quark mass fits compatible at NNLO"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000237,"raw_usage":{"total_tokens":1520,"prompt_tokens":974,"completion_tokens":546,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":590,"completion_tokens_details":{"reasoning_tokens":457}},"tokens_in":590,"tokens_out":546,"duration_ms":242942,"temperature":1.0,"reasoning_tokens":457,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T05:37:12.020416+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute Eq. (20) with $\\alpha_s$ evaluated at $\\mu = M_b \\approx 4.66$ GeV, where $\\alpha_s \\approx 0.21$, instead of $\\alpha_s(M_Z) = 0.118$; the right-hand side then gives about 4.0 GeV rather than 4.38 GeV, which differs from the fitted 4.40 GeV by roughly 0.4 GeV and removes the claimed 99.98% compatibility.","supporting_citations":[{"cited_title":"8 % in the quality of the fit for determination of the b-quark MS running mass mb relative to the b-quark pole mass Mb","cited_arxiv_id":null,"evidence_quote":"Supplies the HERA I+II combined inclusive DIS data that form the main cross-section input for all three QCD fits."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the combined beauty vertex production data that constrain the b-quark mass."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Defines the heavy-flavour matching method on which the treatment is based."},{"cited_title":"Measurement of Charm and Beauty Dijet Cross Sections in Photoproduction at HERA using the H1 Vertex Detector","cited_arxiv_id":"hep-ex/0605016","evidence_quote":"Provides the open-source QCD fitting program used to perform the parameter fits."},{"cited_title":"The pQCD analysis to extract PDFs and $\\alpha_s^{{\\rm NLO}}(M^2_Z)$ from inclusive jet-hadron production data","cited_arxiv_id":"1904.04285","evidence_quote":"Provides the world average masses used as the external comparison for the extracted values."}],"review_version":1}