{"id":"03d477dd-f41d-433c-b66f-3e0bceb29aed","arxiv_id":"2412.07522","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Presents epsilon-factorised master integrals, boundary values, and numerical routines for the planar and non-planar electroweak double-box families relevant to NNLO Moller scattering.","lead":"This paper computes ten families of two-loop Feynman integrals needed for electroweak corrections to electron-electron (Moller) scattering, including the hardest elliptic cases. A generalist should care because these integrals are a bottleneck for matching the precision of the MOLLER weak-mixing-angle experiment at next-to-next-to-leading order.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Boundary constants at (t=0, m^2=∞) are fixed partly by PSLQ over a restricted constant space and validated only at one undocumented AMFlow point; a wrong boundary value would propagate through all solutions.","rationale":"I read the paper as a serious computational achievement: the transformation matrices U_X, the ε-factorised matrices A_X, boundary values, and C++ evaluation routines are all provided as supplementary files, so the ε-factorisation claim is directly checkable and the central construction is plausible. The reader's weakest assumption correctly identifies the boundary values as the least secure input: they are partly PSLQ-fitted under an explicit finite-dimensional constant-space ansatz, and the numerical validation is described only as 'perfect agreement' without quantitative details. This is not an accusation of error; rather, the omission of validation details prevents a reviewer from confirming that the boundary ansatz is complete for all topologies, especially the elliptic ones where the m²→∞ degeneration could in principle leave extra constants. The proposed multi-point, high-precision AMFlow comparison and path-limit PSLQ test would settle the question. I therefore agree with the reader's conditional verdict and see no reason to change it.","tokens_in":37732,"tokens_out":10713,"duration_ms":104265,"concrete_test":"Use AMFlow to evaluate the original I-basis at three to five points spanning the Møller region (e.g., (s,t) = (0.0112, -0.002), (0.010, -0.005), (0.020, -0.001) in GeV², with m=m_Z) with 30-50 digit working precision, convert to the J-basis via the supplied U_X, and compare ε^0...ε^4 coefficients against topo_X_numeric.cc. Separately, take the limit to (xt, x_{1/m²})=(0,0) along two different paths (xt→0 then m²→∞, and m²→∞ then xt→0) and verify that the boundary constants in topo_*_symbolic.mpl are recovered, and confirm each is an exact rational combination of the constants in eq. (54). A discrepancy or a failed PSLQ relation would invalidate the boundary ansatz.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim is the ε-factorised master-integral basis and its evaluation for ten topologies. The most load-bearing unverified input is the boundary vector at (t=0, m^2=∞) in Section 4.2. It is obtained from power counting, vanishing conditions, t-independence, and PSLQ assuming every boundary constant lies in the Q-span of eq. (54): {1, iπ, ζ2, ζ3, iπζ2, ζ4, iπζ3}. The paper does not state which integrals were fixed by PSLQ, the rational relations found, or the precision and number of kinematic points used in the AMFlow check of Section 5, which is described only as 'perfect agreement'. If any boundary constant has a component outside this finite-dimensional ansatz (for example ln2 or a period surviving the m^2→∞ degeneration), the solution of the ε-factorised differential equation is shifted by a homogeneous term, and every reported master integral would be wrong. A single-point, low-precision comparison could still appear consistent. This is load-bearing because the boundary values are an input to the integration, and the central claim depends on their correctness.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper computes the ten planar and non-planar two-loop double-box integral families that arise in NNLO electroweak corrections to Møller scattering, for configurations where three gauge bosons are exchanged between the fermion lines and at least one of them is a photon. For each topology the authors derive a differential system via IBP reduction (Kira), construct (heuristically but verifiably) an epsilon-factorised master-integral basis, fix boundary values at (t=0, m^2=infinity), and integrate the system order-by-order in epsilon. The simpler topologies are expressed in multiple polylogarithms after rationalisation of the square roots; topologies A, A-tilde and B-tilde are given as iterated integrals, with A-tilde and B-tilde involving elliptic periods. C++ numerical routines and Maple files are supplied, and the results are reported to agree perfectly with AMFlow at the benchmark point of eq. (62). Section 6 extracts the leading large logarithms ln(s/m^2), and the paper explicitly states that the existence of an epsilon-factorised basis is in general open and that its construction here is a verified guess.","tokens_in":37913,"tokens_out":20779,"duration_ms":171682,"significance":"The result, if correct, is a genuine step forward for NNLO electroweak corrections: it provides the first complete evaluation of the planar and non-planar double-box integrals for Møller scattering in all mass configurations with up to two massive gauge bosons, including the elliptic non-planar topologies A-tilde and B-tilde. The paper ships a complete, machine-checkable data package - transformation matrices, epsilon-factorised connection matrices, alphabets, boundary values, and C++ routines - and the epsilon-factorised property is verifiable by direct substitution into eq. (43). There are no free parameters: boundary constants are fixed by power-counting and vanishing constraints plus PSLQ over the explicit constant space of eq. (54). The reported cross-check against AMFlow, an independent implementation by a different group, anchors the calculation externally, and the authors honestly flag that a canonical basis is not guaranteed to exist and that their basis resulted from a checked heuristic. If the AMFlow comparison is quantified, I would regard the computational claims as verified rather than merely plausible.","major_comments":[{"comment":"The determination of the boundary constants at (t=0, m^2=infinity) is load-bearing, because these constants are the input to the integration of the epsilon-factorised system and a wrong boundary vector shifts every master integral at all orders. The manuscript states that PSLQ was used for 'some master integrals (of intermediate complexity)' with 'about 50 digits', but it does not specify which integrals were fixed by PSLQ, which by power counting or t-independence, what rational relations were found, and - most importantly - how the high-precision numerical values feeding PSLQ were produced (AMFlow, sector decomposition, or a series expansion at the boundary point). This provenance matters: if the boundary constants are determined from AMFlow values and the final integrals are then validated against AMFlow, the validation is partly circular; naming a different input method would break that loop. Please also address the ansatz of eq. (54): state how constants outside the Q-span of {1, i*pi, zeta2, zeta3, i*pi*zeta2, zeta4, i*pi*zeta3} (for example ln 2 from alternating MZVs in a degenerate massless or t=0 limit) were excluded, or report that PSLQ simply found no relation for such cases.","section":"Section 4.2"},{"comment":"The sentence 'we compared our results to the results of the program AMFlow and found perfect agreement' is not quantified. For a calculation of this size the comparison should state the number of master integrals and epsilon-orders compared, the number of agreeing digits, and the kinematic point or points used. Quantification matters for the boundary-value argument: the path-ordered exponential of an epsilon-factorised flat connection is invertible at each order in epsilon, so the epsilon^j coefficient of any master integral at a generic point equals the corresponding j-th-order boundary constant plus explicit iterated-integral terms; a high-precision check at the benchmark point of eq. (62) would therefore expose an error in any boundary constant at any order. The current one-line statement makes the strength of this check unverifiable. Please also report whether the comparison covers all master integrals or only the selection in Table 3, and consider adding a second kinematic point to guard against accidental near-degeneracies.","section":"Section 5"},{"comment":"Many master integrals of the elliptic topologies A-tilde and B-tilde are defined through functions F that are 'determined by a triangular system of first-order differential equations', with explicit expressions deferred to the supplementary files. The paper does not state the initial conditions that fix these F's at the boundary point (t=0, m^2=infinity), nor whether those initial values are part of the PSLQ/regularity fit, nor whether the modular-transformation freedom of eq. (44) was used to ensure at most simple poles at the chosen boundary point. Since the integration path starts at the boundary point, the initial values of the F's are needed for the master integrals to be uniquely defined and reproducible. Please state these initial conditions explicitly (or state unambiguously where they are imposed in the electronic files) and confirm that the boundary values of the elliptic master integrals and of the F's are mutually consistent.","section":"Appendix B.2.1 and B.2.2"}],"minor_comments":[{"comment":"Eq. (66) defines the leading-log one-form as omega_tilde_1 = d ln L with L = ln(s/m^2), but eq. (67) gives J_LL = sum_j (1/j!)(epsilon L)^j M_tilde_1^j J^(0)_boundary. Solving dJ = epsilon M_tilde_1 d ln L J along the m^2-flow produces factors (epsilon ln L)^j, not (epsilon L)^j; eq. (67) is the standard solution for omega_tilde_1 = d ln(m^2) (or d ln(m^2/s)). This is likely a typo, but as printed the derivation of eq. (67) does not follow from eq. (66). Please correct the definition and ensure the electronic files define M_tilde_1 consistently with it.","section":"Section 6, eqs. (66)-(67)"},{"comment":"The statement 'numerical evaluations with about 50 digits are sufficient' would be more informative if the paper gave the achieved PSLQ residual and at least one representative rational relation, which would allow a reader to gauge the identification confidence.","section":"Section 4.2"},{"comment":"The text says boundary values follow from power counting and t-independence, which would presumably set many constants to zero; yet Table 3 shows J_A-tilde_45 and J_B-tilde_61 with vanishing terms through epsilon^3 and a non-zero epsilon^4 term. Please explain how these higher-order boundary constants are fixed, since the listed constraints alone do not obviously determine them.","section":"Section 4.2 and Table 3"},{"comment":"The definitions J_A-tilde_41, J_A-tilde_44, J_A-tilde_48 and J_B-tilde_57 involve derivatives d/dm^2 of elliptic master integrals and the Wronskians W^(X)_m2 of eq. (38); a sentence stating how these derivatives are evaluated inside the numerical routines would help readers of the C++ code.","section":"Appendix B.2.1 and B.2.2"},{"comment":"The elliptic one-form sets H_A-tilde and H_B-tilde appear only in the supplementary files; please state in the text at least the number of independent one-forms in each set and how they are expressed in terms of the periods and their m^2-derivatives.","section":"Section 4.1"},{"comment":"The claim that the evaluation routines are 'significantly faster than AMFlow' would benefit from a concrete timing example (e.g., seconds per phase-space point at a given precision).","section":"Section 5"},{"comment":"Minor typographical and wording points: in Section 2.5, 'coordinates in a plane 1' should read 'coordinates in the plane' or similar; in Section 1, 'This paper is an example, how techniques...' should be 'This paper is an example of how techniques...'; and the caption of Fig. 1 could say 'photon lines' for clarity. The kinematic region is also repeated verbatim in eq. (5) and eq. (50).","section":"General presentation"}],"recommendation":"major_revision","confidential_remarks":"To the editor: the manuscript is a technical computation paper of a kind this journal handles well; the methods are standard in the elliptic-Feynman-integral community and the results will be of interest to the precision-electroweak and multi-loop communities. The self-citation pattern (refs. [47,48,60,81,89]) is heavy but topically justified, since the computation builds directly on those techniques; I do not see a problem there. The single most important revision request is to quantify the AMFlow comparison in Section 5: if the agreement is at high precision (say 20 or more digits) for all master integrals, the single benchmark point is in principle sufficient to validate the boundary values, because each order of the epsilon-expansion at a generic point contains the corresponding boundary constant with coefficient one. Please require the authors to state the precision, the list of compared integrals, the kinematic coverage, and the provenance of the PSLQ inputs."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Short version: this is a genuinely useful and correct-looking computation. The authors construct epsilon-factorised master-integral bases for ten double-box topologies relevant to Møller scattering, including the non-planar two-mass cases that require elliptic curves. That goes beyond what was available; the three-photon planar case was known, but the massive and non-planar families, and the elliptic treatment for tilde-A and tilde-B, are new. The paper delivers numerical evaluation routines, and the results are checked against AMFlow, an independent program from a different group. That is real evidence, and it anchors the whole calculation.\n\nThe integration setup is the strong part. The IBP reduction, differential equations, and the epsilon-factorised bases are shown explicitly, and Maple/C++ files are attached. The bases are found heuristically, as the authors admit, but that is normal in this field and not a flaw: the result is checkable, and the appendix lists all master integrals in full.\n\nThe soft spots are in the documentation rather than the substance. Boundary values at (t=0, m^2=∞) are fixed partly by power counting and partly by PSLQ over a small set of MZVs. The paper says 50 digits suffice and that AMFlow agrees 'perfectly', but it does not state how many digits the AMFlow comparison reached, how many kinematic points were checked, or which integrals were fixed by PSLQ. A referee will want that spelled out. Even so, the worry that a wrong boundary constant could escape the AMFlow check is mostly theoretical: the benchmark point (62) is far from the boundary, so any constant outside the PSLQ ansatz would shift the solutions and show up unless the AMFlow comparison were extremely loose. The paper should say how tight it is.\n\nAnother gap: the elliptic one-forms for tilde-A and tilde-B live only in the supplementary files, which is inconvenient but acceptable. The paper also clearly states what it does not do: three-massive-boson exchange (genus two) is left out.\n\nBottom line: this is a paper for people working on two-loop electroweak corrections and Feynman integrals. It deserves a serious referee. The only required change is to document the AMFlow validation and the PSLQ boundary constants in detail. I would accept it after that.","headline":"Genuinely useful two-loop integral computation; central claims hold up, weak spot is thin documentation of boundary-value and AMFlow checks.","tokens_in":38482,"tokens_out":4383,"would_cite":true,"duration_ms":40951,"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 constructs an $\\varepsilon$-factorised basis for all ten planar and non-planar double-box integral families relevant to Møller scattering at NNLO, evaluating the elliptic ones as iterated integrals of elliptic periods.","keywords":["electroweak corrections","Møller scattering","double-box integrals","NNLO","epsilon-factorised differential equations","elliptic Feynman integrals","iterated integrals","multiple polylogarithms"],"falsifier":"Evaluate any elliptic master integral, for example $J^{\\tilde{A}}_{45}$ or $J^{\\tilde{B}}_{61}$, at a second point well inside the region $-t \\lesssim s \\ll m^2$ using direct sector-wise numerical integration of the original Feynman parameter integral, and compare the $\\varepsilon^0$ through $\\varepsilon^4$ coefficients with the published routines. A wrong boundary constant or a region-dependent analytic continuation would show up as a mismatch at low order in $\\varepsilon$, while the single benchmark point would not expose it.","tokens_in":37468,"feed_emoji":"⚛️","tokens_out":7861,"duration_ms":78810,"temperature":0.7,"pith_summary":"This paper computes the two-loop double-box Feynman integrals that enter the next-to-next-to-leading-order electroweak corrections to electron–electron (Møller) scattering. Its claim is that all ten planar and non-planar integral families — five topologies each, with zero, one, or two massive gauge bosons and at least one photon — can be brought into an $\\varepsilon$-factorised form, where the dimensional regulator factors out of the differential equation. If the claim is right, the hardest two-loop objects in this process are no longer an obstacle: the simpler families are expressed as multiple polylogarithms, while the two elliptic non-planar families are evaluated as iterated integrals of elliptic periods, with numerical routines that are fast and stable in the heavy-boson limit. That matters because low-energy Møller scattering is a proposed way to measure the weak mixing angle, and missing two-loop corrections are a main theory bottleneck.","feed_headline":"All ten Møller double-box topologies get epsilon-factorised forms","feed_subtitle":"Planar families become polylogarithms; the two elliptic non-planar cases reduce to iterated elliptic-period integrals.","key_machinery":"The load-bearing mechanism is the $\\varepsilon$-factorised (canonical) differential equation. Starting from a pre-canonical basis produced by integration-by-parts reduction, the authors seek a transformation $U$ such that $dJ=\\varepsilon A J$ with $A$ independent of $\\varepsilon$; the check is a direct computation of $U\\tilde{A}U^{-1}-U\\,dU^{-1}$. The construction is guided by maximal cuts in the Baikov representation and uses the square roots $r_1,\\dots,r_8$, rationalisation variable changes, and three elliptic curves ($E(a)$ for $\\tilde{B}$, $E(b)$ and $E(c)$ for $\\tilde{A}$) whose periods satisfy rational differential equations in $s,t,m^2$, so the elliptic one-forms can be written without explicit root expressions. Once the basis is found, iterated integration gives the $\\varepsilon$-expansion; leading large logarithms $\\ln(s/m^2)$ are extracted from a single matrix and the boundary vector.","core_discovery":"The central discovery is an explicit master-integral basis $J = U I$ for every topology, with the property $dJ = \\varepsilon A J$ where $A$ depends on the kinematics but not on $\\varepsilon$. The authors construct it by an educated guess informed by the maximal cut in the loop-by-loop Baikov representation and by elliptic-period calculus, then verify the factorization directly. The payoff is that the integrals can be integrated order by order in $\\varepsilon$: topologies $E$ and $\\tilde{E}$ reduce to harmonic polylogarithms, $B,C,D,\\tilde{C},\\tilde{D}$ to multiple polylogarithms after rationalising the relevant square roots, while $A$ requires iterated integrals, and the non-planar topologies $\\tilde{A}$ and $\\tilde{B}$ involve elliptic curves and are expressed as iterated integrals of elliptic periods. Boundary values are fixed at $t=0$, $m^2=\\infty$ using power counting, vanishing conditions, and a PSLQ reconstruction over a small set of constants, after which a benchmark comparison with an independent numerical solver is reported.","pith_inferences":["The same transformation-building pipeline could be applied to the three-massive-boson double boxes: the maximal-cut guide plus the ten-family alphabets give a concrete ansatz, but the boundary-value procedure would need an enlarged constant set beyond the seven constants used here.","The residual boundary-value risk could be closed by deriving boundary constants from the small-$s/m^2$ expansion of the elliptic periods themselves, and by documenting the precision of the independent numerical comparison at more than one kinematic point.","If the numerical routines are as fast as reported, the limiting step toward a complete NNLO electroweak Møller prediction becomes the remaining amplitude assembly and real-emission contributions rather than the double-box integrals."],"forward_implications":["NNLO electroweak corrections to Møller scattering can now include all double-box diagrams with three exchanged gauge bosons where at least one is a photon.","The two elliptic non-planar topologies are no longer a separate obstruction; they are evaluated in the $\\varepsilon$-expansion as iterated integrals of elliptic periods with supplied numerical routines.","Large logarithms $\\ln(s/m^2)$ at each order in $\\varepsilon$ are read off from a matrix and the boundary values, giving the leading-logarithm behaviour relevant in the low-energy Møller region.","Because the constructed basis is crossing-invariant, the same integrals cover Bhabha scattering, Drell-Yan production, and quark-pair production in $e^+e^-$ annihilation.","The explicit basis and alphabets provide a concrete starting point for the remaining three-massive-boson double-box integrals, whose non-planar case is known to involve a genus-two curve."],"supporting_citations":[{"why":"introduces the $\\varepsilon$-factorised form of differential equations that the new basis must achieve; the paper's central construction is a new basis attaining this form.","marker":"[88]"},{"why":"supplies the kinematic conventions, notation, and integration setup for Møller box integrals that this work extends to the double boxes.","marker":"[81]"},{"why":"provides the elliptic-curve technique for $\\varepsilon$-factorised differential equations involving elliptic periods, used for topologies $\\tilde{A}$ and $\\tilde{B}$.","marker":"[47]"},{"why":"gives the rational-derivative formalism for elliptic periods, including second derivatives in $m^2$, used to avoid explicit expressions for curve roots.","marker":"[48]"},{"why":"shows that the non-planar double box with three massive bosons is related to a genus-two curve, explaining why that case is deferred and marking the boundary of the method.","marker":"[80]"},{"why":"is the independent numerical solver used in the benchmark comparison of Section 5 to validate the boundary values and $\\varepsilon$-expansions.","marker":"[110]"}],"fun_headline_variants":["All Møller double-box topologies factorise in ε","Double-box integrals yield to ε-factorised basis","Møller two-loop integrals: explicit ε-form for all topologies","Elliptic and polylog: complete ε-form for Møller double-box"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The central fragile point is the assignment of boundary constants at $t=0$, $m^2=\\infty$: they are fixed partly by power counting and partially reconstructed numerically with PSLQ, and the paper validates the full set against an independent numerical solver at only one kinematic point.","fun_headline_variants_meta":{"raw":{"variants":["All Møller double-box topologies factorise in ε","Double-box integrals yield to ε-factorised basis","Møller two-loop integrals: explicit ε-form for all topologies","Elliptic and polylog: complete ε-form for Møller double-box"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000439,"raw_usage":{"total_tokens":2154,"prompt_tokens":797,"completion_tokens":1357,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":413,"completion_tokens_details":{"reasoning_tokens":1280}},"tokens_in":413,"tokens_out":1357,"duration_ms":11763,"temperature":1.0,"reasoning_tokens":1280,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T18:45:37.906915+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Evaluate any elliptic master integral, for example $J^{\\tilde{A}}_{45}$ or $J^{\\tilde{B}}_{61}$, at a second point well inside the region $-t \\lesssim s \\ll m^2$ using direct sector-wise numerical integration of the original Feynman parameter integral, and compare the $\\varepsilon^0$ through $\\varepsilon^4$ coefficients with the published routines. A wrong boundary constant or a region-dependent analytic continuation would show up as a mismatch at low order in $\\varepsilon$, while the single benchmark point would not expose it.","supporting_citations":[{"cited_title":"Box integrals with fermion bubbles for low-energy measurements of the weak mixing angle","cited_arxiv_id":"2312.06773","evidence_quote":"supplies the kinematic conventions, notation, and integration setup for Møller box integrals that this work extends to the double boxes."}],"review_version":1}