Pith. sign in

REVIEW 3 major objections 4 minor 54 references

Multi-loop Feynman integrals can be computed as solutions of a second-order PDE in variational form, and the finite element method solves them over an entire phase-space region in one computation.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-01 15:29 UTC pith:X5REPSYB

load-bearing objection Genuinely new idea — FEM on a second-order PDE for Feynman integrals — but the benchmark uses exact analytic boundary/load data, so the practical claim is not yet tested. the 3 major comments →

arxiv 2607.18406 v1 pith:X5REPSYB submitted 2026-07-20 hep-ph hep-th

Feynman Integrals Meet Second-Order Partial Differential Equations

Xiang Chen , Hantian Zhang This is my paper
classification hep-ph hep-th MSC 65N3035J1581Q30
keywords Feynman integralssecond-order PDEsfinite element methodvariational formulationmaster integralsintegration-by-parts reductiontwo-loop four-point integralsphase-space evaluation
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper aims to change how Feynman integrals are evaluated: instead of tracing each kinematic point along a one-dimensional path from known boundary values, it proposes to treat the whole region of phase space as a single equilibrium problem. Each master integral obeys a second-order PDE, and a finite-element discretization of the corresponding variational form turns the calculation into one sparse linear solve that returns the integrals everywhere in the region at once. The proof of concept is a two-loop non-planar four-point integral family; the top-sector integral is reproduced to average relative error 1.4×10⁻⁴ in 14 seconds on a single CPU core. If it holds up, this would give collider phenomenologists a fast, global alternative to point-by-point transport for precision predictions at future e⁺e⁻ machines. The open question the paper leaves is the cost of the boundary data, which currently come from analytic or path-based methods.

Core claim

The paper establishes that multi-loop Feynman integrals can be treated as the equilibrium solution of a second-order PDE over the whole (s,t) phase-space region, rather than as a one-dimensional transport problem. The evidence is a proof-of-concept: the massless two-loop non-planar four-point family is reduced to 12 master integrals, 55 PDEs after ε-expansion, and solved over a bounded region in one sparse linear solve. On a graded mesh with 75,008 interior nodes, the O(ε⁰) top-sector integral matches the known analytic result to average relative error 1.4×10⁻⁴ and maximal error 8.1×10⁻³, in about 14 seconds, with algebraic convergence rate around 0.8 in the average. The authors stress the e

What carries the argument

The central object is the second-order PDE system ∇²I = M(s,t,ε)I for the vector of master integrals, rewritten in variational form as −∫∇Iₙ·∇v − ∫V Iₙ v = ∫f v, with prescribed boundary values γₙ on the boundary of the (s,t) domain. Here V is a potential and f is a load built from lower-sector master integrals; the same integration-by-parts identities used for first-order differential equations produce the system, so no new reduction cost is incurred. The workhorse is a piecewise-linear nodal-basis finite-element discretization: it assembles sparse matrices element by element and solves the resulting linear system once. The variational form only requires the solution and test functions to h

Load-bearing premise

The load-bearing premise is that accurate values of every master integral on the boundary of the phase-space region can be supplied at acceptable cost; the paper treats those boundary values as given input and does not measure their price or error.

What would settle it

Run the identical 12-master-integral finite-element solve with boundary values produced by a path-based transport code instead of analytic formulas, then compare interior values against an independent high-precision evaluation. If the interior error is much larger than the reported 1.4×10⁻⁴ while the analytic-boundary solve stays accurate, the method's practical value collapses.

Watch this falsifier — get emailed when new claim-graph text bears on it.

If this is right

  • If the convergence rate holds at higher loops, the method supplies the full (s,t)-dependence of master integrals from one linear solve, replacing per-point evaluations.
  • The second-order system is built from the same integration-by-parts reduction as first-order DEs, so no additional reduction work is needed to apply the method.
  • Higher-order finite elements and spectral bases are natural upgrades that the paper says could raise the algebraic convergence rate, potentially to exponential in smooth regions.
  • Massive thresholds, where integrals have kinks but remain continuous, stay in the square-integrable-derivative space and can be handled by partitioning the phase space or using complex masses.
  • Coupled master integrals are covered by off-diagonal potential blocks, as demonstrated for the two coupled top-sector integrals solved simultaneously.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The practical bottleneck the paper leaves implicit is boundary data: if accurate boundary values can only come from the same expensive analytic machinery the method is meant to bypass, the one-shot interior solve gains less than the headline numbers suggest. A quantitative cost/accuracy comparison of boundary input versus interior error would settle this.
  • An immediate testable extension is to feed the boundary data from the path-based transport methods cited in the paper rather than from analytic solutions, and compare interior errors on the same 75k-node mesh.
  • Since the benchmark uses a known analytic result, applying the same FEM machinery to a case where no analytic expression exists—for instance, a two-loop four-point integral with internal masses—would show whether the boundary-data requirement can be met in practice.
  • The same variational formulation might be used not just for master integrals but for the amplitudes themselves, if integration-by-parts reduction can be performed pointwise on the mesh; that would skip the separate assembly of master integrals.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

3 major / 4 minor

Summary. The paper proposes to compute multi-loop Feynman integrals by reformulating the system of first-order differential equations as a second-order PDE, deriving a variational form, and solving it with the finite element method over a bounded two-dimensional phase-space region. The method is demonstrated on a massless two-loop non-planar four-point integral family with 12 master integrals. Using analytic results from AsyInt for boundary conditions and lower-sector master integrals, the authors report algebraic convergence of the O(ϵ^0) top-sector integral with average relative errors near 1e-4 and runtimes of up to 14.2 s for the sparse linear solve. The paper claims this is a new 'equilibrium' approach that solves the integral over a broad phase-space region 'once and for all', in contrast to one-dimensional path-based transport methods.

Significance. If the second-order PDE construction is correct and the required boundary/load data can be obtained cheaply, the method could complement existing path-based differential-equations approaches by providing a genuinely global solution over a region. The paper is a proof of concept with concrete numerical benchmarks, and the absence of fitted parameters is a strength: the Galerkin system is assembled directly from IBP-derived coefficients. However, the practical value hinges on two premises that the manuscript does not establish: the PDE system in Eq. (1) is asserted rather than derived, and the numerical experiments rely on exact analytic inputs for all boundary conditions and lower-sector integrals. The reported runtimes exclude the cost of generating those inputs, so the end-to-end computational advantage over direct evaluation is not demonstrated.

major comments (3)
  1. [Methodology, Eq. (1)] The construction of M(s,t,ϵ) is asserted without a derivation or an explicit example. The claim that 'the IBP relations required for constructing the second-order PDEs are identical to the ones for first-order DEs' needs justification: differentiating a first-order system introduces products of coefficient matrices and their derivatives, and the resulting M may be singular or lead to a non-elliptic operator. No ellipticity or well-posedness check is given. This is load-bearing because the Galerkin convergence observed in the benchmarks cannot be expected for indefinite or degenerate operators. Please provide the explicit derivation of Eq. (1) for the top-sector 2×2 system (or a small worked example) and discuss the properties of the operator -Δ - V over the chosen domain.
  2. [Numerical experiments, Tables I–II and captions] The central claim of solving the integrals 'once and for all' is not tested end-to-end. The runtimes in Tables I and II are only LinearSolve times, while the paper states that 'the analytic solutions of these 12 MIs are taken from AsyInt and serve as the boundary conditions and as benchmarks.' Thus both γ_n in Eq. (4) and f in Eq. (2) are exact analytic inputs. The cost of obtaining these inputs via path-based transport or asymptotic methods is never quantified. Without that cost, the 14 s solve cannot be compared with direct evaluation, and the method may amount to interpolating someone else's computation. Please report the total wall-clock time including IBP reduction, mesh generation, assembly, and the numerical determination of boundary conditions and lower-sector integrals, or clearly state that the proof of concept only measures the discretization error of a manufactured problem.
  3. [Numerical experiments, error reporting] Only the O(ϵ^0) term of the top-sector integral I_NPL is benchmarked (with the statement that I_NPL^(num) behaves similarly). However, the method solves 55 PDEs over five ϵ-orders, and the load functions for higher sectors depend on the FEM solutions of lower sectors. Since the lower-sector solutions are analytic inputs in this experiment, the error propagation through the sector hierarchy is not tested. Please report per-sector and per-ϵ-order errors for the full family, at least for the largest mesh, to support the claim that the whole system is solved reliably.
minor comments (4)
  1. [Eq. (6)] The decomposition In = u + u0 is confusingly written. 'u0|Ω/∂Ω = 0' is not a standard notation; please clarify that u0 is a finite-extension of the boundary data into the interior, e.g. by harmonic extension or by nodal interpolation on the boundary layer.
  2. [Page 3, error definition] The phrase 'the absolute value is equivalent to the L2-norm for complex numbers' is imprecise; for a complex number the absolute value is the modulus, not an L2 norm.
  3. [Summary] The statement 'requiring only a limited number of boundary conditions as input' understates the data requirement: every sector needs boundary values on ∂Ω and the load function f requires all lower-sector master integrals over the entire domain. Please qualify this sentence.
  4. [Methodology, Eq. (9)–(11)] The assembly of the Galerkin matrices is described only in words. A short pseudo-code or a reference to a standard FEM assembly routine would improve reproducibility.

Circularity Check

0 steps flagged

No circularity: boundary data and lower-sector integrals are external inputs; benchmark is an independent manufactured-solution test.

full rationale

The derivation chain is: master integrals satisfy first-order differential equations (standard, cited to Refs. [1-4]); differentiating these and using IBP gives the second-order PDEs of Eqs. (1)-(2), with M, V, and f constructed from reduction data, not from the target solution. The variational formulation (Eqs. (3)-(4)) and Galerkin/FEM discretization (Eqs. (8)-(12)) are standard numerical machinery with no fitted parameters. Boundary conditions gamma_n are explicitly declared input data in Eq. (4), and the paper states: 'The analytic solutions of these 12 MIs are taken from AsyInt and serve as the boundary conditions and as benchmarks for comparison with the numerical results on the interior nodes.' This is a manufactured-solution test: the analytic solution supplies Dirichlet data on the boundary and is then compared with the FEM interior solution. The interior result is not used to construct the PDE, the mesh, or the boundary data, so the agreement is a genuine check of the PDE formulation and discretization. The likely practical dependency — that accurate boundary and lower-sector data must be obtained by other path-based or analytic methods — is acknowledged ('For various cases where analytic boundary conditions are not available, they can be computed numerically to high precision using path-based transport methods'). That is an input-cost or end-to-end limitation, not a circular reduction. The self-citation to AsyInt [49] supplies independent test data and is not load-bearing for the derivation of the PDE method itself. No step in the paper reduces its claimed prediction to its own inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 4 axioms · 0 invented entities

The method introduces no fitted free parameters and no new physical entities. It relies on domain assumptions: existence of the second-order PDE from IBP, well-posedness of the variational problem, accurate external boundary/lower-sector data, and empirical FEM convergence. These are the honest inputs the paper pulls from outside its own derivation.

axioms (4)
  • domain assumption Feynman master integrals satisfy a linear second-order PDE ∇²I = M(s,t,ϵ)I obtainable by iterated IBP reduction.
    Invoked in Eq. (1)-(2); the paper asserts this follows from taking derivatives and applying IBP iteratively, but the construction of M is not shown for the benchmark family.
  • domain assumption The Dirichlet problem for the variational form is well-posed on the chosen bounded phase-space domain.
    Required for the Galerkin method (Eq. (7)); the paper assumes H^1 regularity and does not prove coercivity or well-posedness for arbitrary V and f.
  • domain assumption Boundary conditions γ_n and lower-sector master integrals can be precomputed accurately by existing path-based/analytic methods.
    Used throughout the numerical experiments: 'The analytic solutions of these 12 MIs are taken from AsyInt and serve as the boundary conditions.' No cost or error analysis of this input is given.
  • domain assumption Linear FEM converges algebraically for these solutions on the chosen domains.
    Convergence rates in Eq. (15)-(17) are empirical; no a priori error estimate for this nonstandard PDE is established.

pith-pipeline@v1.3.0-alltime-deepseek · 9645 in / 14840 out tokens · 130314 ms · 2026-08-01T15:29:33.809392+00:00 · methodology

0 comments
read the original abstract

We propose a second-order partial differential equation method to solve multi-loop Feynman integrals as an equilibrium problem. As a proof-of-concept demonstration, we perform a Galerkin discretization of the corresponding variational form and employ the finite element method to compute two-loop four-point Feynman integrals. This method can solve the integral over a broad region of multi-dimensional phase space once and for all. This work establishes a new connection between perturbative quantum field theory and modern partial differential equation methods, with the potential to streamline a wide range of phenomenological applications.

Figures

Figures reproduced from arXiv: 2607.18406 by Hantian Zhang, Xiang Chen.

Figure 1
Figure 1. Figure 1: FIG. 1. Feynman diagram of massless non-planar integral. [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Convergence in terms of average and maximal errors [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Runtime of [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

54 extracted references · 42 linked inside Pith

  1. [1]

    A. V. Kotikov, Differential equations method: New tech- nique for massive Feynman diagrams calculation, Phys. Lett. B254, 158 (1991)

    1991

  2. [2]

    Remiddi, Differential equations for Feynman graph amplitudes, Nuovo Cim

    E. Remiddi, Differential equations for Feynman graph amplitudes, Nuovo Cim. A110, 1435 (1997), arXiv:hep- th/9711188

    arXiv 1997

  3. [3]

    Gehrmann and E

    T. Gehrmann and E. Remiddi, Differential equations for two-loop four-point functions, Nucl. Phys. B580, 485 (2000), arXiv:hep-ph/9912329

    Pith/arXiv arXiv 2000

  4. [4]

    J. M. Henn, Multiloop integrals in dimensional regu- larization made simple, Phys. Rev. Lett.110, 251601 (2013), arXiv:1304.1806 [hep-th]

    Pith/arXiv arXiv 2013

  5. [5]

    Caola, A

    F. Caola, A. Von Manteuffel, and L. Tancredi, Diphoton Amplitudes in Three-Loop Quantum Chromodynamics, Phys. Rev. Lett.126, 112004 (2021), arXiv:2011.13946 [hep-ph]

    Pith/arXiv arXiv 2021

  6. [6]

    Gehrmann, J

    T. Gehrmann, J. Henn, P. Jakubˇ c ´ ık, J. Lim, C. C. Mella, N. Syrrakos, L. Tancredi, and W. J. Torres Bobadilla, Graded transcendental functions: an application to four- point amplitudes with one off-shell leg, JHEP12, 215, arXiv:2410.19088 [hep-th]

    Pith/arXiv arXiv

  7. [7]

    Chicherin, T

    D. Chicherin, T. Gehrmann, J. M. Henn, P. Wasser, Y. Zhang, and S. Zoia, All Master Integrals for Three- Jet Production at Next-to-Next-to-Leading Order, Phys. Rev. Lett.123, 041603 (2019), arXiv:1812.11160 [hep- ph]

    Pith/arXiv arXiv 2019

  8. [8]

    Abreu, J

    S. Abreu, J. Dormans, F. Febres Cordero, H. Ita, and B. Page, Analytic Form of Planar Two-Loop Five-Gluon Scattering Amplitudes in QCD, Phys. Rev. Lett.122, 082002 (2019), arXiv:1812.04586 [hep-ph]

    Pith/arXiv arXiv 2019

  9. [9]

    Badger, D

    S. Badger, D. Chicherin, T. Gehrmann, G. Heinrich, J. M. Henn, T. Peraro, P. Wasser, Y. Zhang, and S. Zoia, Analytic form of the full two-loop five-gluon all-plus he- licity amplitude, Phys. Rev. Lett.123, 071601 (2019), arXiv:1905.03733 [hep-ph]

    Pith/arXiv arXiv 2019

  10. [10]

    Abreu, D

    S. Abreu, D. Chicherin, H. Ita, B. Page, V. Sotnikov, W. Tschernow, and S. Zoia, All Two-Loop Feynman In- tegrals for Five-Point One-Mass Scattering, Phys. Rev. Lett.132, 141601 (2024), arXiv:2306.15431 [hep-ph]

    Pith/arXiv arXiv 2024

  11. [11]

    Y. Liu, A. Matijaˇ si´ c, J. Miczajka, Y. Xu, Y. Xu, and Y. Zhang, Analytic computation of three-loop five-point Feynman integrals, Phys. Rev. D112, 016021 (2025), arXiv:2411.18697 [hep-ph]

    Pith/arXiv arXiv 2025

  12. [12]

    Chicherin, Y

    D. Chicherin, Y. Wu, Z. Wu, Y. Xu, S.-Q. Zhang, and Y. Zhang, Complete computation of all three-loop five- point massless planar integrals, (2025), arXiv:2512.17330 [hep-ph]

    arXiv 2025

  13. [13]

    J. M. Henn, A. Matijaˇ si´ c, J. Miczajka, T. Peraro, Y. Xu, and Y. Zhang, A computation of two-loop six-point Feyn- man integrals in dimensional regularization, JHEP08, 027, arXiv:2403.19742 [hep-ph]

    Pith/arXiv arXiv

  14. [14]

    Abreu, P

    S. Abreu, P. F. Monni, B. Page, and J. Usovitsch, Planar six-point Feynman integrals for four-dimensional gauge theories, JHEP06, 112, arXiv:2412.19884 [hep-ph]

    Pith/arXiv arXiv

  15. [15]

    J. Henn, A. Matijaˇ si´ c, J. Miczajka, T. Peraro, Y. Xu, and Y. Zhang, Complete Function Space for Planar Two- Loop Six-Particle Scattering Amplitudes, Phys. Rev. Lett.135, 031601 (2025), arXiv:2501.01847 [hep-ph]

    Pith/arXiv arXiv 2025

  16. [16]

    R. N. Lee, A. V. Smirnov, and V. A. Smirnov, Solving dif- ferential equations for Feynman integrals by expansions near singular points, JHEP03, 008, arXiv:1709.07525 [hep-ph]. 6

    Pith/arXiv arXiv

  17. [17]

    Moriello, Generalised power series expansions for the elliptic planar families of Higgs + jet production at two loops, JHEP01, 150, arXiv:1907.13234 [hep-ph]

    F. Moriello, Generalised power series expansions for the elliptic planar families of Higgs + jet production at two loops, JHEP01, 150, arXiv:1907.13234 [hep-ph]

    Pith/arXiv arXiv 1907

  18. [18]

    Liu, Y.-Q

    X. Liu, Y.-Q. Ma, and C.-Y. Wang, A Systematic and Ef- ficient Method to Compute Multi-loop Master Integrals, Phys. Lett. B779, 353 (2018), arXiv:1711.09572 [hep- ph]

    Pith/arXiv arXiv 2018

  19. [19]

    M. Fael, F. Lange, K. Sch¨ onwald, and M. Steinhauser, A semi-analytic method to compute Feynman integrals applied to four-loop corrections to the MS-pole quark mass relation, JHEP09, 152, arXiv:2106.05296 [hep-ph]

    Pith/arXiv arXiv

  20. [20]

    M. Fael, F. Lange, K. Sch¨ onwald, and M. Steinhauser, Massive Vector Form Factors to Three Loops, Phys. Rev. Lett.128, 172003 (2022), arXiv:2202.05276 [hep-ph]

    Pith/arXiv arXiv 2022

  21. [21]

    Armadillo, R

    T. Armadillo, R. Bonciani, S. Devoto, N. Rana, and A. Vicini, Evaluation of Feynman integrals with arbi- trary complex masses via series expansions, Comput. Phys. Commun.282, 108545 (2023), arXiv:2205.03345 [hep-ph]

    Pith/arXiv arXiv 2023

  22. [22]

    Davies, G

    J. Davies, G. Mishima, M. Steinhauser, and D. Well- mann, Double Higgs boson production at NLO in the high-energy limit: complete analytic results, JHEP01, 176, arXiv:1811.05489 [hep-ph]

    Pith/arXiv arXiv

  23. [23]

    Davies, G

    J. Davies, G. Mishima, K. Sch¨ onwald, M. Steinhauser, and H. Zhang, Higgs boson contribution to the leading two-loop Yukawa corrections to gg→HH, JHEP08, 259, arXiv:2207.02587 [hep-ph]

    Pith/arXiv arXiv

  24. [24]

    X. Chen, X. Guan, C.-Q. He, Z. Li, X. Liu, and Y.-Q. Ma, Complete two-loop unrenormalized electroweak cor- rections to e+e-→HZ, Phys. Rev. D114, 013003 (2026), arXiv:2209.14953 [hep-ph]

    Pith/arXiv arXiv 2026

  25. [25]

    Armadillo, R

    T. Armadillo, R. Bonciani, S. Devoto, N. Rana, and A. Vicini, Two-loop mixed QCD-EW corrections to neu- tral current Drell-Yan, JHEP05, 072, arXiv:2201.01754 [hep-ph]

    Pith/arXiv arXiv

  26. [26]

    Bi, L.-H

    H.-Y. Bi, L.-H. Huang, R.-J. Huang, Y.-Q. Ma, and H.- M. Yu, Electroweak Corrections to Double Higgs Produc- tion at the LHC, Phys. Rev. Lett.132, 231802 (2024), arXiv:2311.16963 [hep-ph]

    Pith/arXiv arXiv 2024

  27. [27]

    Davies, K

    J. Davies, K. Sch¨ onwald, M. Steinhauser, and H. Zhang, Analytic next-to-leading order electroweak corrections to Higgs boson pair production at high energies, (2026), arXiv:2603.08789 [hep-ph]

    arXiv 2026

  28. [28]

    Abadaet al.(FCC), FCC Physics Opportunities: Fu- ture Circular Collider Conceptual Design Report Volume 1, Eur

    A. Abadaet al.(FCC), FCC Physics Opportunities: Fu- ture Circular Collider Conceptual Design Report Volume 1, Eur. Phys. J. C79, 474 (2019)

    2019

  29. [29]

    Abadaet al.(FCC), FCC-ee: The Lepton Collider: Future Circular Collider Conceptual Design Report Vol- ume 2, Eur

    A. Abadaet al.(FCC), FCC-ee: The Lepton Collider: Future Circular Collider Conceptual Design Report Vol- ume 2, Eur. Phys. J. ST228, 261 (2019)

    2019

  30. [30]

    Benediktet al.(FCC), Future Circular Collider Feasibility Study Report: Volume 1, Physics, Exper- iments, Detectors, Eur

    M. Benediktet al.(FCC), Future Circular Collider Feasibility Study Report: Volume 1, Physics, Exper- iments, Detectors, Eur. Phys. J. C85, 1468 (2025), arXiv:2505.00272 [hep-ex]

    Pith/arXiv arXiv 2025

  31. [31]

    CEPC Conceptual Design Report: Volume 1 - Accelera- tor, (2018), arXiv:1809.00285 [physics.acc-ph]

    Pith/arXiv arXiv 2018

  32. [32]

    Donget al.(CEPC Study Group), CEPC Conceptual Design Report: Volume 2 - Physics & Detector, (2018), arXiv:1811.10545 [hep-ex]

    M. Donget al.(CEPC Study Group), CEPC Conceptual Design Report: Volume 2 - Physics & Detector, (2018), arXiv:1811.10545 [hep-ex]

    Pith/arXiv arXiv 2018

  33. [33]

    de Blaset al., Physics Briefing Book: In- put for the 2026 update of the European Strategy for Particle Physics 10.23731/CYRM-2025-008 (2025), arXiv:2511.03883 [hep-ex]

    J. de Blaset al., Physics Briefing Book: In- put for the 2026 update of the European Strategy for Particle Physics 10.23731/CYRM-2025-008 (2025), arXiv:2511.03883 [hep-ex]

    arXiv 2026

  34. [34]

    Denner, S

    A. Denner, S. Dittmaier, M. Roth, and L. H. Wieders, Electroweak corrections to charged-current e+ e- —>4 fermion processes: Technical details and further results, Nucl. Phys. B724, 247 (2005), [Erratum: Nucl.Phys.B 854, 504–507 (2012)], arXiv:hep-ph/0505042

    Pith/arXiv arXiv 2005

  35. [35]

    X. Chen, D. Chicherin, E. Fox, N. Glover, M. Marcoli, V. Sotnikov, H. Sun, H. Zhang, and S. Zoia, Four-Jet Rate in Electron-Positron Annihilation at Orderαs4, Phys. Rev. Lett.136, 251901 (2026), arXiv:2602.18185 [hep-ph]

    arXiv 2026

  36. [36]

    Bredt, T

    P. Bredt, T. Banno, M. H¨ ofer, S. Iguro, W. Kilian, Y. Ma, J. Reuter, and H. Zhang, Chasing the Two- Higgs-Doublet Model via Electroweak Corrections at e+e- Colliders, Phys. Rev. Lett.136, 081801 (2026), arXiv:2509.05421 [hep-ph]

    arXiv 2026

  37. [37]

    Liu and Y

    Y. Liu and Y. Zhang, CHESS: CHEbyshev pSeudo- Spectral transport for Feynman integral differential equa- tions, (2026), arXiv:2606.26691 [hep-ph]

    Pith/arXiv arXiv 2026

  38. [38]

    Abreu, A

    S. Abreu, A. Guerreiro, and B. Page, Chebyshev Ap- proximations of Feynman Integrals for Collider Physics, (2026), arXiv:2607.02411 [hep-ph]

    Pith/arXiv arXiv 2026

  39. [39]

    Hidding, DiffExp, a Mathematica package for com- puting Feynman integrals in terms of one-dimensional series expansions, Comput

    M. Hidding, DiffExp, a Mathematica package for com- puting Feynman integrals in terms of one-dimensional series expansions, Comput. Phys. Commun.269, 108125 (2021), arXiv:2006.05510 [hep-ph]

    Pith/arXiv arXiv 2021

  40. [40]

    Liu and Y.-Q

    X. Liu and Y.-Q. Ma, AMFlow: A Mathematica pack- age for Feynman integrals computation via auxiliary mass flow, Comput. Phys. Commun.283, 108565 (2023), arXiv:2201.11669 [hep-ph]

    Pith/arXiv arXiv 2023

  41. [41]

    R. M. Prisco, J. Ronca, and F. Tramontano, LINE: Loop Integrals Numerical Evaluation, JHEP07, 219, arXiv:2501.01943 [hep-ph]

    Pith/arXiv arXiv

  42. [42]

    Petit Ros` as and W

    P. Petit Ros` as and W. J. Torres Bobadilla, Fast evalu- ation of Feynman integrals for Monte Carlo generators, JHEP09, 210, arXiv:2507.12548 [hep-ph]

    Pith/arXiv arXiv

  43. [43]

    Baur and C

    G. Baur and C. Duhr, IterInt: Evaluating iterated inte- grals via differential equations, (2026), arXiv:2606.02744 [hep-ph]

    Pith/arXiv arXiv 2026

  44. [44]

    Czakon and L

    M. Czakon and L. Tancredi, Solution of Canonical Differ- ential Equations for Integrals on Arbitrary Geometries, (2026), arXiv:2606.30354 [hep-ph]

    Pith/arXiv arXiv 2026

  45. [45]

    Huang, X

    R.-J. Huang, X. Liu, and Y.-Q. Ma, AMFlow 2.0: signifi- cant algorithmic and software improvements for Feynman integral evaluation, (2026), arXiv:2607.08477 [hep-ph]

    Pith/arXiv arXiv 2026

  46. [46]

    P. G. Ciarlet,The Finite Element Method for Elliptic Problems(Society for Industrial and Applied Mathemat- ics, 2002)

    2002

  47. [47]

    J. B. Tausk, Non-planar massless two-loop Feynman di- agrams with four on-shell legs, Phys. Lett. B469, 225 (1999), arXiv:hep-ph/9909506

    Pith/arXiv arXiv 1999

  48. [48]

    R. N. Lee, LiteRed 1.4: a powerful tool for reduction of multiloop integrals, J. Phys. Conf. Ser.523, 012059 (2014), arXiv:1310.1145 [hep-ph]

    Pith/arXiv arXiv 2014

  49. [49]

    Zhang, Massive two-loop four-point Feynman in- tegrals at high energies with AsyInt, JHEP09, 069, arXiv:2407.12107 [hep-ph]

    H. Zhang, Massive two-loop four-point Feynman in- tegrals at high energies with AsyInt, JHEP09, 069, arXiv:2407.12107 [hep-ph]

    Pith/arXiv arXiv

  50. [50]

    B. Q. Guo and I. Babuˇ ska, Theh-pversion of the fi- nite element method. part 1: The basic approximation results, Computational Mechanics1, 21 (1986)

    1986

  51. [51]

    A. T. Patera, A spectral element method for fluid dy- namics: Laminar flow in a channel expansion, Journal of Computational Physics54, 468 (1984)

    1984

  52. [52]

    Actis, A

    S. Actis, A. Denner, L. Hofer, J.-N. Lang, A. Scharf, and S. Uccirati, RECOLA: REcursive Computation of One-Loop Amplitudes, Comput. Phys. Commun.214, 140 (2017), arXiv:1605.01090 [hep-ph]. 7

    Pith/arXiv arXiv 2017

  53. [53]

    Buccioni, J.-N

    F. Buccioni, J.-N. Lang, J. M. Lindert, P. Maierh¨ ofer, S. Pozzorini, H. Zhang, and M. F. Zoller, OpenLoops 2, Eur. Phys. J. C79, 866 (2019), arXiv:1907.13071 [hep- ph]

    Pith/arXiv arXiv 2019

  54. [54]

    Armadillo, S

    T. Armadillo, S. Devoto, M. Dradi, and A. Vicini, Towards the two-loop electroweak corrections to the Drell-Yan process: the complete fermionic contributions, (2026), arXiv:2605.28667 [hep-ph]

    Pith/arXiv arXiv 2026