{"id":"99521faa-0f95-43cd-a051-cdcbf358491a","arxiv_id":"2411.08802","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"NLO QCD corrections reduce the predicted widths of Z to two jets plus one photon by 6.03% and to two jets plus two photons by 12.39%.","lead":"This paper computes next-to-leading-order QCD corrections to two rare Z boson decays, where the Z turns into a quark-antiquark pair plus one or two photons. The corrections change the predicted decay rates by about 6% and 12%, which matters for precision Standard Model tests at future colliders.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The quoted NLO shifts rest on an undocumented local K=5 cut; the claimed 4<K<6 plateau is not shown and no Monte Carlo uncertainties are reported, so the 6.03% and 12.39% values could be biased.","rationale":"The paper does a genuine new NLO calculation and has several internal checks (pole cancellation, dipole subtraction, alpha-stability, gauge invariance, and agreement with MadGraph at LO and with Ref. [13] for Z to q qbar at NLO). Those checks support the soundness of the infrastructure. The weakest point is the K-factor cut, exactly as the reader identified. The paper says there is a plateau but supplies no evidence for it; with no numerical uncertainties anywhere in Tables I and II, the quoted precision of 6.03% and 12.39% cannot be assessed. This does not invalidate the calculation, but it justifies the CONDITIONAL verdict: the central numbers should not be used at face value until the K dependence is demonstrated to be negligible. I see no reason to move the verdict to reject or accept; the condition is precisely that the stability evidence be shown.","tokens_in":17215,"tokens_out":3684,"duration_ms":35986,"concrete_test":"Rerun the Z to q qbar gamma gamma (and Z to q qbar gamma) NLO integrals with the same phase-space generator and statistics for K = 4, 4.5, 5, 5.5, 6, and 10, reporting Gamma_NLO and RI with Monte Carlo errors at each K. If the RI values move by more than about 0.2% (or more than one combined sigma) across 4<K<6, the quoted central values are K-dependent. As a secondary check, re-evaluate a sample of the discarded points (K>5) in higher precision to confirm they are numerical artifacts rather than physical large virtual corrections.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the size of the NLO QCD correction. The paper's only guard against numerical instabilities in the pentagon amplitudes is Eq. (1), a local K-factor cut that discards phase-space points where |2 Re(M_virt M_LO*)|/|M_LO|^2 exceeds K=5. The text states a plateau for 4<K<6 but gives no table, plot, or numerical spread across that interval. This matters because the virtual contribution is signed; removing points with large |M_virt| preferentially removes negative corrections and can shift RI systematically toward more negative values. The q qbar gamma gamma channel has a 12.39% correction, so even a 1-2% bias is a large fraction of the claimed effect. Additionally, the paper reports no Monte Carlo statistical uncertainty for Gamma_LO, Gamma_NLO, or RI, so the reader cannot tell whether the plateau is flat within errors or whether the two quoted decimals are meaningful. The K-cut is thus the most load-bearing unverified assumption in the otherwise internally consistent calculation.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript reports a calculation of the NLO QCD corrections to the Z-boson decay channels Z -> q qbar, Z -> q qbar gamma, and Z -> q qbar gamma gamma, with the aim of providing predictions for the corresponding two-jet-plus-photon(s) observables at the HL-LHC and at future e+e- colliders. The amplitudes are generated with FeynArts/FormCalc and processed with in-house FORM routines, tensor integrals are reduced with OVReduce, scalar integrals are evaluated with OneLoop, and phase-space integration is performed with AMCI/VEGAS. Infrared singularities are handled with Catani-Seymour dipole subtraction, including a finite alpha parameter. The authors list several validation checks: pole cancellation, dipole cancellation, alpha variation, Ward identities, agreement with MadGraph at tree level, and an exact match of the NLO Z -> q qbar width to Ref. [13]. Their central numerical results are the relative increments RI = +3.77% for Z -> q qbar, -6.03% for Z -> q qbar gamma, and -12.39% for Z -> q qbar gamma gamma, together with jet-level widths and a set of differential distributions.","tokens_in":17395,"tokens_out":6435,"duration_ms":60947,"significance":"If the quoted results are robust, this is a useful calculation: the two negative NLO corrections are the first complete one-loop QCD results for these rare Z-decay channels, and the paper provides concrete numbers for observables that could be tested at future Tera-Z factories and possibly at the HL-LHC. The calculation is a direct Standard Model computation with no free parameters fitted to the target observables, and the listed internal checks (pole cancellation, Ward identities, MadGraph comparison, exact NLO Z -> q qbar match) are appropriate and give reasonable confidence in the technical implementation. The main weakness is numerical documentation: the local K-factor stability is asserted but not demonstrated, and no Monte Carlo statistical uncertainties are quoted, so the precision implied by the two-decimal percentages is not currently supported. These are fixable presentation and validation gaps rather than apparent errors in the derivation.","major_comments":[{"comment":"The local K-factor cut defined in Eq. (1) is the only stabilization described for the pentagon amplitudes, and the text states that there is a plateau region for the final results between 4 < K < 6 without showing the actual K-dependence of Gamma_NLO or RI. Because the cut removes phase-space points where the pointwise ratio |2 Re(M_virt M_LO*)|/|M_LO|^2 exceeds K, and because the virtual-Born interference is signed, the removed region can have a nonzero integral and can therefore bias the total correction if the plateau is not exactly flat. Please include a table or plot of Gamma_NLO and RI as a function of K, with statistical errors at each K value, or alternatively demonstrate that the total discarded weight is numerically negligible. This is necessary to verify the quoted 12.39% (and, if the cut is used there, the 6.03%) correction.","section":"§III, Eq. (1); §IV, Table I"},{"comment":"No Monte Carlo statistical uncertainties are reported for Gamma_LO, Gamma_NLO, RI, or any of the differential distributions. The headline numbers are quoted to two decimals (6.03% and 12.39%), but without error bars the reader cannot assess whether the K=5 value is stable within integration errors, whether the alpha-variation check in §III.C.3 is meaningfully passed, or whether the final digits are significant. Please quote the statistical uncertainty for each width and for the derived RI. This documentation is essential for the central numerical claim of the paper.","section":"§IV, Tables I and II; Figs. 8–16"}],"minor_comments":[{"comment":"The sentence 'following Ref[]' contains an empty citation; please insert the appropriate reference for the alpha-modified dipole kernel.","section":"§III.B"},{"comment":"There are several typographical errors: 'prosesses' in §IV, 'diagramns' in the Fig. 7 caption, 'to to mass' in the Fig. 15(d) caption, and 'spitting matrix' in §III.B should be corrected.","section":"Throughout"},{"comment":"The quantity RI is called the 'relative increment'; for negative corrections this is more naturally called a relative change, and the caption could be adjusted for clarity.","section":"§IV, Table I"},{"comment":"The paper does not state the number of phase-space points or the VEGAS iterations used to obtain the quoted widths; adding this information would improve reproducibility.","section":"§IV.A"},{"comment":"The claims of agreement with MadGraph and of an exact match to Ref. [13] are not quantified; please state the numerical comparison and the achieved agreement.","section":"§III.C"}],"recommendation":"major_revision","confidential_remarks":"I concur with the stress-test concern: the K-cut plateau is the most load-bearing unverified element of the numerical claim, and the absence of Monte Carlo error bars makes the two-decimal percentages unsupported as presented. However, because the authors state that they scanned the plateau and because the rest of the calculation is standard and well cross-checked, this is a documentation gap rather than a demonstrated error. A major revision that adds the K-scan data and statistical uncertainties would address the concern directly. I do not see a circularity problem: no parameter is fitted to the target observables, and the calculation is an honest direct SM computation."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a workmanlike first NLO QCD calculation for Z -> q qbar gamma and Z -> q qbar gamma gamma. The claimed corrections, -6.03% and -12.39%, are credible but rest on a stability claim the paper doesn't actually document.\n\nWhat's genuinely new: no previous NLO QCD result for these two decay modes, as far as I can tell. The authors do the full one-loop calculation with spinor helicity, Catani-Seymour dipole subtraction, and a set of checks that are worth taking seriously: pole cancellation, dipole cancellation, alpha-parameter variation, Ward identities, tree-level agreement with MadGraph, and an exact match to the known NLO width for Z->q qbar. That last cross-check is strong evidence that the setup is right.\n\nThe soft spot is the local K-factor cut in Eq. (1). They throw away points where |2 Re(M_virt M_LO*)|/|M_LO|^2 > K, set K=5, and say there's a plateau for 4<K<6. But they don't show the plateau, and there are no Monte Carlo statistical errors anywhere in the paper. The stress-test note is right that this is the load-bearing unverified assumption. The virtual contribution is signed; cutting on large |M_virt| can bias the correction. For the 12.39% effect, even a 1-2% bias is a sizable fraction. That said, this is a standard numerical stabilization technique, and the internal consistency checks make a gross error unlikely. The fix is easy: show the K-scan or at least report the integration error.\n\nMinor things: the b-quark treated as massless is standard but worth stating, the HL-LHC cross-section estimates are rough, and the paper's organization is a bit repetitive. None of these bother me.\n\nWho it's for: phenomenologists working on Tera-Z physics and rare Z decays. It deserves a serious referee, but the referee should ask for the plateau scan and Monte Carlo errors before the numbers are taken at face value.","headline":"First NLO QCD calculation for Z -> q qbar gamma and Z -> q qbar gamma gamma gives credible percent-level corrections; the main caveat is an unshown K=5 cut plateau and no Monte Carlo errors.","tokens_in":17940,"tokens_out":2125,"would_cite":true,"duration_ms":134859,"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 reports the first complete next-to-leading-order QCD calculation of the rare $Z$-boson decays $Z\\to q\\bar{q}\\gamma$ and $Z\\to q\\bar{q}\\gamma\\gamma$, and finds that the corrections reduce the partial widths by $6.03\\%$ and…","keywords":["Z boson rare decays","NLO QCD corrections","two jets plus photon","two jets plus two photons","decay width","dipole subtraction","future e+e- colliders"],"falsifier":"Recompute the NLO widths with the local K-factor cut removed or moved far outside the reported $4<K<6$ plateau, or with a different infrared-subtraction scheme, and check whether the $-6.03\\%$ and $-12.39\\%$ shifts survive within Monte Carlo errors; an independent NLO calculation of $Z\\to q\\bar{q}\\gamma\\gamma$ would settle the claim directly.","tokens_in":1857,"feed_emoji":"⚛️","tokens_out":2298,"duration_ms":75595,"temperature":0.7,"pith_summary":"This paper computes the first complete next-to-leading-order (NLO) QCD corrections to the rare $Z$-boson decay channels $Z\\to q\\bar{q}\\gamma$ and $Z\\to q\\bar{q}\\gamma\\gamma$, treating the final quarks as massless and applying photon isolation cuts suitable for hadron-collider analyses. It finds that the NLO corrections reduce the inclusive partial width of $Z\\to q\\bar{q}\\gamma$ by about $6.03\\%$ and of $Z\\to q\\bar{q}\\gamma\\gamma$ by about $12.39\\%$ relative to leading order. At the jet level the reductions are larger, reaching about $20\\%$ for the two-jet-plus-one-photon mode and more than $30\\%$ for the one-jet-plus-two-photon mode. These NLO-improved widths are the numbers that future high-statistics electron-positron Z-pole measurements will need to compare against.","feed_headline":"Rare Z decays shrink by 6% to 12% at next-to-leading order","feed_subtitle":"One-loop QCD widths for Z to two jets plus one or two photons give the numbers future Z-pole runs will compare against.","key_machinery":"The calculation is carried by a one-loop QCD amplitude machinery: helicity amplitudes built from spinor products and vector currents, dimensional regularization in the 't Hooft-Veltman scheme, on-shell renormalization for the massless quarks, and a dipole subtraction scheme that cancels infrared singularities between virtual and real-emission contributions before numerical Monte Carlo integration. For the pentagon diagrams that appear in $Z\\to q\\bar{q}\\gamma\\gamma$, the paper applies a local K-factor cut, discarding rare phase-space points where the ratio of the virtual amplitude to the Born amplitude exceeds a fixed threshold, set to $K=5$ after a reported plateau for $4<K<6$. The calculation is checked by pole cancellation, gauge invariance through Ward identities, agreement with an independent tree-level generator, and reproduction of the known NLO width of $Z\\to q\\bar{q}$.","core_discovery":"On the paper's own terms, the central result is quantitative: the one-loop QCD corrections to $Z\\to q\\bar{q}\\gamma$ and $Z\\to q\\bar{q}\\gamma\\gamma$ are negative and sizable, shifting the partial widths from $6.30$ MeV to $5.92$ MeV and from $11.3$ keV to $9.90$ keV respectively, while the same correction to the benchmark $Z\\to q\\bar{q}$ width is positive, $3.77\\%$. The paper argues that these shifts are large enough that any high-precision test of $Z$-boson decays must include them, and that the corrected channels have observable rates at a high-luminosity hadron collider or at a future electron-positron Z-pole machine.","pith_inferences":["Inference: the larger negative shift in the two-photon channel suggests that each additional photon roughly doubles the size of the NLO QCD correction, a pattern that could be tested by computing the NLO correction to $Z\\to q\\bar{q}\\gamma\\gamma\\gamma$.","Inference: because the local K-factor cut only removes extreme virtual-to-Born ratios, publishing the integrated width as a function of $K$ over a wider range would make the reported plateau directly auditable by other groups.","Inference: the same dipole-subtraction setup can be applied to other rare $Z$-boson decays with photons in the final state, such as $Z$ to $q\\bar{q}$ plus a Higgs boson, where NLO QCD corrections have not yet been computed.","Inference: experimental analyses that count jets should use the jet-level NLO widths rather than the inclusive widths, since real gluon emission moves events among one-, two-, and three-jet topologies and steepens the effective corrections."],"forward_implications":["The corrected inclusive partial widths are $5.92$ MeV for $Z\\to q\\bar{q}\\gamma$ and $9.90$ keV for $Z\\to q\\bar{q}\\gamma\\gamma$, corresponding to NLO QCD shifts of $-6.03\\%$ and $-12.39\\%$.","At jet level the NLO corrections are larger: $-13.75\\%$ and $-20.26\\%$ for the one-jet and two-jet modes with one photon, and $-31.80\\%$ and $-21.47\\%$ for the one-jet and two-jet modes with two photons.","The NLO processes produce events with one, two, or three jets accompanied by one or two photons, so jet-clustering and photon-isolation cuts become part of the predicted observables.","At a high-luminosity hadron collider the estimated production cross sections, roughly $0.13$ nb for $q\\bar{q}\\gamma$ and $0.22$ pb for $q\\bar{q}\\gamma\\gamma$, make these decay channels experimentally testable.","For future electron-positron Z-pole runs, the $Z\\to q\\bar{q}\\gamma$ partial width alone is larger than the current theoretical uncertainty assigned to unaccounted higher-order corrections, so the NLO shift is material at the projected precision."],"supporting_citations":[{"why":"Supplies the dipole subtraction scheme used to cancel infrared singularities between real and virtual contributions.","marker":"[62]"},{"why":"Supplies the spinor-helicity formalism used for tree-level and one-loop helicity amplitudes.","marker":"[48]"},{"why":"Provides the known NLO width of $Z\\to q\\bar{q}$ that the paper reproduces as a validation check.","marker":"[13]"},{"why":"Supplies the scalar one-loop integral library used in the numerical evaluation of loop amplitudes.","marker":"[56]"},{"why":"Supplies the phase-space generation routine used in the Monte Carlo integration.","marker":"[60]"},{"why":"Provides the prescription to discard phase-space points with unnaturally large virtual-to-Born ratios, the basis of the local K-factor cut.","marker":"[61]"},{"why":"Supplies the tensor reduction routine used to reduce one-loop tensor integrals to scalar integrals.","marker":"[53]"}],"fun_headline_variants":["NLO QCD trims Z to jets plus γ by 6%, plus 2γ by 12%","Rare Z decay widths drop 6-12% under one-loop QCD","Z to two jets with photons: NLO cuts widths by up to 12%","NLO QCD reduces Z->qqγ width 6%, Z->qqγγ width 12%"],"cache_read_input_tokens":20096,"weakest_assumption_plain":"The load-bearing premise is that the local K-factor cut at $K=5$ does not bias the integrated widths; the paper reports a plateau for $4<K<6$ but does not show how the results behave outside that range.","fun_headline_variants_meta":{"raw":{"variants":["NLO QCD trims Z to jets plus γ by 6%, plus 2γ by 12%","Rare Z decay widths drop 6-12% under one-loop QCD","Z to two jets with photons: NLO cuts widths by up to 12%","NLO QCD reduces Z->qqγ width 6%, Z->qqγγ width 12%"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001101,"raw_usage":{"total_tokens":4554,"prompt_tokens":870,"completion_tokens":3684,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":486,"completion_tokens_details":{"reasoning_tokens":3584}},"tokens_in":486,"tokens_out":3684,"duration_ms":27002,"temperature":1.0,"reasoning_tokens":3584,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T21:18:53.244055+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the NLO widths with the local K-factor cut removed or moved far outside the reported $4<K<6$ plateau, or with a different infrared-subtraction scheme, and check whether the $-6.03\\%$ and $-12.39\\%$ shifts survive within Monte Carlo errors; an independent NLO calculation of $Z\\to q\\bar{q}\\gamma\\gamma$ would settle the claim directly.","supporting_citations":[{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Supplies the dipole subtraction scheme used to cancel infrared singularities between real and virtual contributions."},{"cited_title":"Kleiss and W","cited_arxiv_id":null,"evidence_quote":"Supplies the tensor reduction routine used to reduce one-loop tensor integrals to scalar integrals."}],"review_version":1}