{"id":"a446bc3e-27dd-45f2-92b9-b318d765aba7","arxiv_id":"2506.15363","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"The authors derive analytic expressions for all 303 three-loop master integrals with one massive propagator appearing in the mixed QCD-electroweak corrections to the quark form factor, expressed through generalized polylogarithms with interdependent arguments.","lead":"This paper computes all 303 three-loop master integrals needed for the O(alpha alpha_s^2) corrections to the quark form factor, covering Feynman diagrams with a single massive Z or W boson. The results are expressed as generalized polylogarithms and verified numerically, providing a key ingredient for next-to-next-to-next-to-leading order mixed QCD-electroweak predictions for Drell-Yan production at the LHC.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Boundary constants for w-letter integrals are fixed by PSLQ against an unproven constant basis; the analytic claim is therefore conditional until that basis is established.","rationale":"Section 3.3.1 explicitly flags the missing constant basis for the w-letters, making it the only self-admitted gap in an otherwise standard differential-equations computation. The AMFlow agreement at 50-digit precision is strong evidence that the final numerical values are correct; however, the headline claim is about analytic GPL expressions, and PSLQ identification of constants is only as reliable as the assumed basis. Because the paper does not identify which MIs are affected or provide an independent derivation of those boundary constants, the analytic claim is not fully proven. This matches the reader's weakest assumption and supports a CONDITIONAL verdict rather than unconditional acceptance. A dedicated constant-basis computation or independent boundary derivation would close the gap.","tokens_in":37548,"tokens_out":10310,"duration_ms":104849,"concrete_test":"Establish the constant basis for GPLs with letters {w3,w4} up to the maximal weight used (weight 6 or 7) by extending the sixth-root-of-unity evaluation tables of Ref. [117] to these quadratic letters, and check every PSLQ-reconstructed boundary constant in the ancillary files against this basis. If every constant is exactly reproduced, the analytic claim is fully supported; if any constant is missing or mismatched, the affected master integrals are incorrect.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is that all 303 master integrals admit analytic GPL expressions. The weakest link is the determination of boundary constants for integrals whose GPL alphabet contains the quadratic letters w3,w4=(-3±√5)/2 (Section 3.3.1). These constants are reconstructed from AMFlow numerics using PSLQ under the assumption that they lie in an 'anticipated' basis of MZVs, ln(2), Li_n(1/2) and cyclotomic constants. The paper states that establishing the full basis of constant GPLs for {w3,w4} is deferred to future investigation. If that basis is incomplete, PSLQ can return a spurious rational relation at finite precision, and the resulting boundary constants—and hence the affected analytic master integrals—would be incorrect. The AMFlow checks validate the numerical values of the MIs, but they do not prove that the identified constants are the correct analytic ones; a wrong constant would typically be detected at other kinematic points, so the checks make gross errors unlikely, yet the analytic identification itself remains unproven. The paper also does not state which of the 303 MIs rely on the w-letter reconstruction, so the scope of the conditional part is unclear.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper presents the computation of 303 three-loop master integrals that appear in the O(alpha alpha_s^2) corrections to the quark form factor in Feynman diagrams containing a single massive vector boson. The authors perform an IBP reduction to obtain a basis of master integrals, set up first-order differential equations in x = -s/m_V^2, decouple the subsystems into factorizable higher-order equations, and solve them by the method of variation of constants. Because the square roots appearing in the problem cannot be rationalized by a single change of variables, they employ several transformations (x, x_l, x_n, x_i) and express the results in terms of generalized polylogarithms with a simple alphabet but multiple interdependent arguments. Boundary conditions are fixed partly by Feynman-parameter evaluations and regularity conditions and partly by AMFlow numerics combined with PSLQ reconstruction. The analytic expressions are provided in ancillary files, a numerical table at x = 1/11 is given, and the results are checked against AMFlow to 50-digit precision at several kinematic points.","tokens_in":37753,"tokens_out":6738,"duration_ms":67737,"significance":"If correct, this is a substantial technical contribution: the 303 master integrals are a necessary ingredient for the three-loop mixed QCD-electroweak corrections to the quark form factor, with direct relevance to Drell-Yan phenomenology at the LHC. The treatment of simultaneously non-rationalizable square roots via concurrent transformations, leading to GPLs with several interdependent arguments, is methodologically interesting and likely to be useful beyond this specific calculation. The main strengths are the complete reduction to 303 master integrals, the explicit GPL alphabets, and the high-precision numerical validation against AMFlow. The principal weakness is that boundary constants for integrals involving the quadratic letters w3 and w4 are reconstructed from AMFlow numerics via PSLQ under an anticipated but unproven constant basis; the paper itself states that establishing this basis is left to future work. The numerical checks validate the numerical values of the integrals but do not by themselves establish the exact analytic form of those boundary constants.","major_comments":[{"comment":"The boundary conditions for master integrals whose GPL alphabet contains the quadratic letters w3 and w4 are fixed by PSLQ under an 'anticipated' set of constants, and the paper states that establishing the full basis of constant GPLs involving the alphabet {w1, ..., w4} is planned for future investigation. Because these constants enter the analytic expressions of the affected master integrals, the central claim that all 303 master integrals admit analytic GPL representations is conditional on that constant basis being complete. If the basis is incomplete, PSLQ can produce a spurious rational relation at finite precision, and the resulting boundary constants could be incorrect even if the numerical checks at a few kinematic points pass. Please either prove the constant relations needed for the reconstruction, or determine the affected boundary constants through an independent non-PSLQ method (for example, by imposing regularity or evaluating the integrals at additional points in different kinematic regions), or clearly identify which of the 303 master integrals rely on the w-letter reconstruction and state the analytic results with that caveat. The AMFlow checks validate the numerical values of the integrals; they do not by themselves prove that the PSLQ-identified constants are the correct analytic ones.","section":"3.3.1"},{"comment":"The numerical checks are described only as 'perfect agreement with AMFlow output across several kinematic points', without specifying the points, the number of points, or which master integrals were checked. Given that boundary constants for a subset of the integrals are obtained by PSLQ reconstruction, the validation should report the actual kinematic points, the regions of x they cover (for example x > 4, 0 < x < 4, and negative x), and the achieved precision for a representative sample of the affected integrals. This information is necessary to assess the strength of the statement that all 303 master integrals are correct and to ensure that the branch choices and analytic continuations of the GPLs have been exercised.","section":"4.1"}],"minor_comments":[{"comment":"The definitions of the integral families I13 and I14 are missing their closing braces; in the current text both lines end with a comma and no closing brace.","section":"3.1"},{"comment":"The label B5,5 is used twice, once for the limit q^2 -> 0 and once for the limit q^2 -> -1, with two different gamma-function expressions; please disambiguate these labels.","section":"3.3.1"},{"comment":"The boundary-condition block refers to diagrams B4,1, B5,1, T6,1, and so on with the statement that the thick line represents the massive propagator, but no figure showing these diagrams is included in the text; without a diagram-to-label correspondence the boundary conditions cannot be matched to the master integrals.","section":"3.3.1"},{"comment":"The notation for GPLs with multiple arguments should be defined more precisely; the paper should state explicitly how GPLs with argument x, x_l, x_n, and x_i are combined in the final expressions and how the alphabets containing r3, r4, w3, w4, i, and -i are handled in the GPL definition.","section":"2"},{"comment":"The sentence describing integrand sizes says 'tens of megabits', which is a data-size unit; it should read 'tens of megabytes' or 'millions of characters'.","section":"3.3"},{"comment":"The table of 303 numerical entries in the main text is very large; a smaller representative subset in the text with the full table in the ancillary files would improve readability.","section":"4"}],"recommendation":"major_revision","confidential_remarks":"The computational framework appears sound, and the 50-digit AMFlow checks provide strong evidence that the numerical values of the master integrals are correct. The main obstacle to acceptance is the conditional status of the PSLQ-reconstructed boundary constants for the w-letter alphabets. If the authors can either establish the needed constant basis or clearly delimit the claim and identify the affected integrals, I would be supportive of publication. A more detailed validation section reporting the actual kinematic points and checked integrals would also strengthen the paper."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Bottom line: this is a real computational advance, not a repackaging. The authors compute 303 three-loop master integrals for single-massive-propagator topologies that feed O(alpha alpha_s^2) corrections to the quark form factor. Those topologies have not appeared in the literature, and the GPL alphabet with multiple interdependent arguments, handled by concurrent transformations, is a genuine technical development. The AMFlow checks at 50-digit precision across several kinematic points are exactly the right kind of evidence, and the boundary conditions are fixed by external numerics and independent Feynman-parameter integrals, so nothing is circular.\n\nThe soft spot is the one the authors flag in Section 3.3.1: boundary constants for integrals whose GPL alphabet includes w3,w4 = (-3 +/- sqrt(5))/2 are reconstructed by PSLQ under an \"anticipated\" basis of MZVs, ln(2), Li_n(1/2), and cyclotomic constants. They state that establishing the full constant basis for these letters is planned for future investigation. That makes the analytic expressions for those MIs conditional: if the basis is incomplete, PSLQ can lock onto a wrong rational relation at finite precision, and the constants—and whatever MIs depend on them—would be incorrect. The multiple AMFlow checks make gross errors unlikely, but they validate numerical values, not the identification of the analytic constants. The paper also does not say which of the 303 MIs rely on the w-letter reconstruction, so the reader cannot gauge the scope of the conditional part. That is worth fixing in revision, even if the fix is a table listing the affected MIs.\n\nMinor point: the full analytic expressions live only in ancillary files. That is normal in this field and acceptable; the tables in the text give a numerical reference point at x=1/11, which helps.\n\nThe citation pattern looks fine; the self-citations point to closely related Drell-Yan work, which is relevant rather than padding.\n\nWho this is for: anyone computing mixed QCD-EW corrections at three loops, and people working on multi-loop integrals with non-rationalizable roots. It is a technical paper and should go to peer review. I would recommend acceptance conditional on either establishing the w-letter constant basis or clearly separating which results depend on it and adding independent support for those results.","headline":"Novel 303-master-integral computation with strong numeric checks, but the analytic claim is conditional while the w-letter boundary constants rest on an unproven PSLQ basis.","tokens_in":38259,"tokens_out":1963,"would_cite":true,"duration_ms":21023,"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":"All 303 master integrals for the three-loop quark form factor with a single massive boson are computed analytically in generalized polylogarithms, despite square roots that resist a single rationalizing transformation.","keywords":["quark form factor","master integrals","three-loop","mixed QCD-electroweak corrections","generalized polylogarithms","differential equations","Drell-Yan process","boundary constants"],"falsifier":"Take one of the master integrals whose boundary constants were fixed through the $w$-letter alphabet and evaluate its closed form at a kinematic point that was not used in the PSLQ reconstruction (for example, a negative-$x$ threshold point), comparing against a direct high-precision numerical evaluation of the original Feynman integral by an independent method such as sector decomposition; disagreement at the claimed precision would show the assumed constant basis is incomplete, while agreement would corroborate it.","tokens_in":37352,"feed_emoji":"⚛️","tokens_out":15214,"duration_ms":124817,"temperature":0.7,"pith_summary":"The paper aims to supply the analytic master integrals needed for the three-loop mixed strong-electroweak ($O(\\alpha\\alpha_s^2)$) corrections to the quark form factor, specifically the 303 integrals coming from Feynman diagrams with a single massive vector boson in the loop. It claims that every one of these integrals can be written as a combination of generalized polylogarithms, even though the intermediate calculation runs into several square roots that cannot all be removed by one change of variables. The resolution is to apply several rationalizing transformations side by side, obtaining polylogarithms with a small set of allowed integration kernels (letters) but with several interdependent arguments. If correct, these integrals, together with known integrals for the other topologies, complete the master-integral data needed for N3LO mixed QCD-electroweak Drell-Yan predictions at the LHC.","feed_headline":"All 303 quark form factor master integrals solved analytically","feed_subtitle":"They clear the way for N3LO mixed QCD-electroweak Drell-Yan corrections at the LHC.","key_machinery":"The load-bearing mechanism is the first-order differential-equation system for the 303 master integrals, organized into an upper block-triangular form and solved bottom-up by decoupling each block into higher-order equations whose operators factorize into first-order pieces. Rather than searching for a canonical epsilon-form, the paper keeps the system in this factorized form and uses variation of constants. The rationalizing variables $x_l$, $x_n$, and $x_i$ are the central objects that make the final answer simple: each removes one of the square roots $\\sqrt{(4-x)x}$, $\\sqrt{1-4x}$, and the arctangent structure $x^{-3/2}\\tan^{-1}(\\sqrt{x})$, and the final polylogarithms have the small alphabets listed above. Boundary constants that cannot be obtained by Feynman-parameter evaluation or regularity are reconstructed with PSLQ from high-precision values supplied by an auxiliary mass flow method, using known tables for sixth-root-of-unity constants up to weight six.","core_discovery":"The central discovery is that the square-root obstructions in these three-loop integrals are individually rationalizable, just not simultaneously. The paper defines the variables $x_l$, $x_n$, and $x_i$ through $x=(1+x_l)^2/x_l=x_n/(1+x_n)^2=-x_i^2$, uses different transformations in different sectors, and then integrates the mixed-argument polylogarithmic inhomogeneities piecewise. The upshot is that all 303 master integrals are expressible in generalized polylogarithms (iterated integrals built from logarithmic kernels) with alphabets $\\{-1,0,1\\}$ for $x$, $\\{-1,0,1,r_3,r_4\\}$ for $x_l$, $\\{-1,0,1,r_3,r_4,w_3,w_4\\}$ for $x_n$, and $\\{-1,0,1,i,-i\\}$ for $x_i$, where $r_3,r_4$ are the primitive sixth roots of unity and $w_3,w_4=(-3\\pm\\sqrt{5})/2$. Boundary conditions are fixed by Feynman-parameter values, regularity, and PSLQ reconstruction from high-precision numerical values; the resulting expressions agree with an independent numerical method to at least 50 digits.","pith_inferences":["A natural extension left implicit in the paper is that once the constant basis for the $w_3,w_4$ alphabet is established and proven complete, the same PSLQ-based boundary reconstruction could be applied to other multi-scale integral families with multiple square roots, avoiding canonical bases altogether.","If the full set of master integrals for the $O(\\alpha\\alpha_s^2)$ quark form factor is expressible in generalized polylogarithms, the complete form factor amplitude is likely also polylogarithmic; the paper computes only the single-massive-boson sector, so this remains to be checked for the contributions built on [94].","A testable extension would be to derive explicit analytic-continuation rules for the $w_3,w_4$ polylogarithms across the thresholds $s=0$ and $s=-4m_V^2$ and verify them numerically, since the mixed-argument representation makes continuation less transparent than in a single-variable representation."],"forward_implications":["Together with the integrals presented in [94] for the other topologies, these 303 master integrals provide the complete set needed to obtain the three-loop $O(\\alpha\\alpha_s^2)$ quark form factor and, from it, N3LO mixed QCD-electroweak Drell-Yan predictions.","Because the expressions are analytic generalized polylogarithms with small alphabets, they can be evaluated numerically to high precision quickly with standard libraries, making them directly usable in cross-section codes.","The computation demonstrates a workable alternative to canonical-form reductions: when square roots cannot all be rationalized by one transformation, solving the system bottom-up with several rationalizing variables and piecewise integration still yields compact analytic results.","The reported agreement with an independent numerical method to at least 50 digits across several kinematic points indicates the analytic results can be used reliably in the physical region."],"supporting_citations":[{"why":"Introduces the differential-equation method for Feynman integrals that is the backbone of the computation.","marker":"[82]"},{"why":"Extends the differential-equation method to multiloop integrals and underlies the decoupling strategy used here.","marker":"[84]"},{"why":"Defines the generalized polylogarithms in which all final results are expressed.","marker":"[91]"},{"why":"Provides the previously computed master integrals for the other three-loop mixed QCD-EW topologies that these 303 integrals complement.","marker":"[94]"},{"why":"The integration-by-parts reduction code that produces the 303 master integrals from the 25 integral families.","marker":"[96]"},{"why":"Provides the alternative reduction implementation used to improve the reduction and to obtain the integral basis.","marker":"[98]"},{"why":"The integer-relation (PSLQ) algorithm used to reconstruct analytic boundary constants from high-precision values.","marker":"[112]"},{"why":"Introduces the auxiliary mass flow method that supplies the high-precision numerical values used for boundary reconstruction and checks.","marker":"[113]"},{"why":"Provides the tables of constants at sixth roots of unity up to weight six used for the r_3,r_4 alphabets.","marker":"[117]"}],"fun_headline_variants":["303 mixed QCD-EW three-loop integrals now analytic","All quark form factor three-loop master integrals solved","Three-loop form factor integrals expressed via polylogs","Square roots handled: 303 three-loop integrals solved","Quark form factor: three-loop integrals in closed form"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The boundary constants for integrals whose polylogarithms use the $w_3,w_4$ letters are reconstructed numerically under the assumption that they lie in a specific, not-yet-proven-complete basis of multiple zeta values and related constants; if that basis misses a constant, the affected analytic integrals would be wrong even though they currently match high-precision numerics.","fun_headline_variants_meta":{"raw":{"variants":["303 mixed QCD-EW three-loop integrals now analytic","All quark form factor three-loop master integrals solved","Three-loop form factor integrals expressed via polylogs","Square roots handled: 303 three-loop integrals solved","Quark form factor: three-loop integrals in closed form"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000384,"raw_usage":{"total_tokens":2021,"prompt_tokens":925,"completion_tokens":1096,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":541,"completion_tokens_details":{"reasoning_tokens":1019}},"tokens_in":541,"tokens_out":1096,"duration_ms":10884,"temperature":1.0,"reasoning_tokens":1019,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T19:35:20.886204+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Take one of the master integrals whose boundary constants were fixed through the $w$-letter alphabet and evaluate its closed form at a kinematic point that was not used in the PSLQ reconstruction (for example, a negative-$x$ threshold point), comparing against a direct high-precision numerical evaluation of the original Feynman integral by an independent method such as sector decomposition; disagreement at the claimed precision would show the assumed constant basis is incomplete, while agreement would corroborate it.","supporting_citations":[{"cited_title":"Evaluating Multiple Polylogarithm Values at Sixth Roots of Unity up to Weight Six","cited_arxiv_id":"1512.08389","evidence_quote":"Provides the tables of constants at sixth roots of unity up to weight six used for the r_3,r_4 alphabets."}],"review_version":2}