Pith. sign in

REVIEW 3 major objections 4 minor 3 cited by

Efficient Computation of One-Loop Feynman Integrals and Fixed-Branch Integrals to High Orders in $\epsilon$

T0 review · 3 major / 4 minor · reviewed 2026-08-10 · deepseek-v4-flash

Pith's one-line read The dimension-changing transformation rewrites one-loop Feynman integrals and fixed-branch integrals as a single contour integral over an auxiliary mass, so that solving differential equations once in one dimension yields the full…

desk verdict The DCT trick is a real efficiency win for high-order one-loop epsilon expansions, but the unproved contour condition needs scrutiny before the method is used outside the paper's safe zone. read the letter →

arxiv 2412.21054 v1 pith:C2JEK2Y3 submitted 2024-12-30 hep-ph

classification hep-ph MSC 81T1881T1581-08 PACS 11.15.Bt
keywords dimension-changingtransformationone-loopFeynmanintegralsfixed-branchepsilonexpansionauxiliarymassflowmastercontourdeformationdimensionalregularization
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper introduces a dimension-changing transformation (DCT) that relates a one-loop Feynman integral or a fixed-branch integral (a Feynman-parameter integral representing a branch of a multi-loop diagram) in any spacetime dimension to the same family of integrals carrying an auxiliary mass in a different dimension, through a single contour integral over that mass. Because the dimensional regulator $\epsilon$ enters the contour integral only through the factor $\eta^{\delta-1}/\Gamma(\delta)$, the $\epsilon$-expansion can be extracted by expanding that factor and integrating the auxiliary-mass series term by term, so the differential equations in the auxiliary mass need to be solved only once, at one nonsingular dimension. The transformation is implemented in an open-source C++ package and demonstrated on a five-point and a six-point one-loop family and a two-loop six-point fixed-branch family, producing Laurent expansions to $O(\epsilon^{10})$ with 14-15 digit agreement against existing results. If correct, the method removes the need to solve differential equations at many $\epsilon$ values for one-loop and fixed-branch integrals, which is the main computational bottleneck for high-precision predictions.

What carries the argument

The central object is the dimension-changing transformation, Eq. (10), a contour-integral identity expressing an integral in $D$ dimensions as an integral over the auxiliary mass $\eta$ of the same integral in $D_0=D-2\delta$ dimensions, with the contour taken along the negative imaginary axis. The work it does is to separate the $\epsilon$-dependence from the hard part of the calculation: after the auxiliary-mass integrals are expanded in power series around $\eta=0$, $\eta=\infty$, and regular points, each term is integrated against $\eta^{\delta-1}$ analytically, and $\Gamma(\delta)$ and $\eta^{\delta-1}$ are expanded in $\epsilon$ before integration. The identity is obtained by splitting the loop momentum into $D_0$ parallel and $D-D_0$ orthogonal components, integrating out the orthogonal angular variables, and then deforming the radial $\eta$ integration path from the positive real axis to the negative imaginary axis. The paper also gives the auxiliary-mass differential equation used to generate the series expansions, with boundary conditions at $\eta\to\infty$.

What would settle it

Set up a one-loop integral with a Breit-Wigner propagator mass $m^2 = M^2 - iM\Gamma$ and physical scattering kinematics, compute the auxiliary-mass integral along the negative imaginary axis, and scan the fourth quadrant numerically for poles; if a pole is found, the DCT contour integral in Eq. (10) will not reproduce the original integral, and comparing DCT with a direct numerical integration at that point would reveal the discrepancy.

Watch

Extended reading notes

Core claim

The central claim is the identity in Eq. (10): $$I^\Delta_{\vec\nu}(X) = \frac{1}{\Gamma(\delta)} \int_{-i\infty}^{-i0^+} d\eta\, \$eta^{{\delta-1}}$\, I^{\$\Delta$-\delta}_{\vec\nu}(X,\eta),$$ where $I^{\Delta-\delta}_{\vec\nu}(X,\eta)$ is the fixed-branch integral with auxiliary mass $\eta$ in a dimension $D_0=2(\Delta-\delta)$, and $I^\Delta_{\vec\nu}(X)$ is the same integral in a different dimension $D=2\Delta$ after removing the auxiliary mass. The same relation reduces to a one-loop Feynman integral when there is a single branch. The authors argue that, because the $\epsilon$-dependence sits in the prefactor and in $\eta^{\delta-1}$, one can solve the auxiliary-mass differential equations once at a conveniently chosen non-integer dimension and then produce the full Laurent series in the dimensional regulator at very low cost. Their numerical examples confirm that this one-solution-then-transform route reproduces known results to 14-15 digits, including cases with vanishing Gram determinants.

Load-bearing premise

The whole derivation rests on the claim that the auxiliary-mass integral $I^{D_0}(\eta)$ has no poles in the fourth quadrant of the $\eta$-plane, so the contour may be rotated from the positive real axis to the negative imaginary axis; the paper states this condition but never proves it for arbitrary kinematics, complex masses, or general fixed-branch integrals.

Editorial extensions

If this is right

  • High-order $\epsilon$ expansions of one-loop master integrals can be obtained by solving the auxiliary-mass differential equations at a single value of $\epsilon$, so the wall time grows only with the cost of the DCT integration, not with repeated differential-equation solves.
  • The same efficiency extends to fixed-branch integrals, which appear in the representation of multi-loop integrals; this opens a route to high-precision multi-loop predictions once the fixed-branch decomposition is known.
  • The method handles real and complex kinematic configurations, including Breit-Wigner propagators, provided the auxiliary-mass integral has no fourth-quadrant poles, so it applies to physical scattering regions rather than only Euclidean kinematics.
  • Exceptional cases such as vanishing Gram determinants are handled by the construction of the differential-equation system, and arbitrary precision is available through the package's support for high-precision complex arithmetic.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Editorial inference: the same dimension-changing trick could be applied to any master-integral family whose auxiliary-mass differential equations are already solved, converting a stock of existing solutions at one dimension into $\epsilon$-expansions at many dimensions without new differential-equation work.
  • Editorial inference: if the fourth-quadrant analyticity condition can be proven for broad classes of kinematics, DCT may become a general-purpose tool for evaluating Feynman integrals directly at complex or non-integer dimensions, not just at $D=4-2\epsilon$.
  • Editorial inference: the Laplace-transform variant sketched in the outlook suggests a family of dimension-changing identities parameterized by different integration weights; comparing their numerical convergence might yield even faster contour choices for specific integrals.
Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

3 major / 4 minor

Summary. This paper introduces a dimension-changing transformation (DCT) that relates one-loop Feynman integrals and fixed-branch integrals (FBIs) with an auxiliary mass in D0 dimensions to the corresponding integrals without the auxiliary mass in D dimensions, via an integral over the auxiliary parameter η (Eqs. (4)-(5) and (10)). The authors propose to solve auxiliary-mass differential equations at a single non-singular dimension and then evaluate the DCT integral term by term, which yields the full ϵ expansion at low cost. The method is implemented in an open-source C++ package and demonstrated on a five-point one-loop topology, a massless six-point one-loop topology, and a two-loop three-branch FBI example, with comparisons to AMFlow and LoopTools.

Significance. The DCT identity is derived from first principles with no fitted parameters, and the manuscript provides an open-source implementation and numerical evidence of efficiency gains over AMFlow for high ϵ orders. If the contour-deformation assumption is justified, the method would be a useful tool for high-precision one-loop and FBI computations. However, the paper currently relies on an unpublished companion paper for the FBI differential equations, and the central analytic-continuation step is asserted rather than proved, so the significance is conditional on resolving these points.

major comments (3)
  1. [II.A (Eq. (5)) and II.B (Eq. (10))] The deformation of the integration contour from the positive real axis to the negative imaginary axis is stated to be valid provided the auxiliary-mass integral has no poles in the fourth quadrant, but no proof or sufficient condition is given. For one-loop integrals, the poles of I^{D0}_ν(η) in the η-plane are threshold singularities; for complex masses or general kinematic regions their location is not obvious. For FBIs, the integrand in Eq. (6) has singularities at η = -F(X,y), and nothing in Sections II.B or III shows that these singularities cannot enter the fourth quadrant for allowed complex branch parameters X or complex masses. If a fourth-quadrant pole is crossed, Eq. (5)/(10) differs from the original integral by the residue of that pole, and the package would return an incorrect result without warning. The numerical examples in Section IV all use real X and Feynman-prescription kinematics and therefore cannot expose this issue. Please provide a proof or a precise, checkable condition for the absence of fourth-quadrant poles, and add a numerical test with complex masses and/or complex branch parameters.
  2. [III, Eq. (11)] The differential equations for FBIs are quoted from Ref. [40], which is unpublished and cited with a placeholder arXiv number (arXiv:2412.xxxxx). Because the FBI differential equations are essential for the method, the manuscript is not self-contained and cannot be independently checked. The authors should either include a derivation of Eq. (11) and the construction of the matrix S and the constants C_b and z_α, or cite a publicly available version of [40], and should clearly state which results are inherited from [40] and which are new in this work.
  3. [Table II] Table II lists the precision of the hexagon calculation using 1000-digit boost::mpc as '-3 -5 -4 ...' for the ϵ orders. This is inconsistent with the definition precision = -log10(ε) with ε the relative error, since negative values would imply relative errors larger than one. If these are typographical errors (e.g., missing decimal points or exponents), please correct them; as written, the table does not support the claimed high-precision validation of the hexagon example.
minor comments (4)
  1. [III, after Eq. (10)] The sentence beginning 'since ∆ − δ is fixed' should have 'Since' capitalized at the start of the sentence.
  2. [IV, first example] The word 'diferential' should be spelled 'differential' in the description of solving the η differential equations.
  3. [References] Reference [40] is listed with a placeholder arXiv number 'arXiv:2412.xxxxx'; this should be updated once the companion paper is available.
  4. [Fig. 3 caption] The caption states 'Time spent for computing a pentagon family...' but the y-axis is labeled 't(ms)'; please define the symbol t explicitly in the caption.

Circularity Check

0 steps flagged · score 1.0 of 10

DCT identity is a self-contained change of variables; no circular step found.

full rationale

The central claim is the dimension-changing identity, Eq. (4) for one-loop Feynman integrals and Eq. (10) for fixed-branch integrals. Both are derived directly in the text rather than assumed: Eq. (4) follows from splitting the loop momentum into parallel and transverse components, integrating over the transverse angular variables, and substituting l_perp^2 = eta; Eq. (10) follows from the Feynman-parameter representation by explicitly evaluating the integral over eta, giving Gamma(delta) I^{Delta+delta}_nu(X). The epsilon-dependence enters only through eta^{delta-1} and Gamma(delta), which are expanded analytically, while the auxiliary-mass input I^{Delta-delta}_nu(eta) is solved once at a fixed non-integer dimension D0 = 7/13. No parameter appearing in the examples is fitted to the target integrals, and the reported high-epsilon coefficients are not imposed by the input data by construction. The contour deformation from Eq. (4) to Eq. (5) is justified by the stated condition that I^{D0}_nu(eta) has no fourth-quadrant poles; this condition is asserted rather than proven for arbitrary complex kinematics and fixed-branch integrals, which is a rigor gap but not a circularity. The numerical inputs are obtained by solving differential equations cited from the authors' companion work [40] and following the AMFlow method [36,39]; these are dependencies with overlapping authorship, but they supply the integrand, not the claimed transformation. The DCT identity is independently cross-checked against directly computed AMFlow expansions in the examples. Under the review rules, self-citation and same-author validation are not by themselves circularity. Therefore no step in the derivation chain reduces to its own input by construction; the paper is assigned a low score reflecting only minor self-citation dependencies, not actual circular reasoning.

Assumptions & free parameters 0 free parameters · 5 assumptions · 0 invented entities

No free parameters are fitted. The computation relies on the analytic-continuation and contour assumptions summarized above, on the companion paper [40] for fixed-branch integrals, and on AMFlow as solver and validator. No new particles, forces, dimensions, or other invented entities are introduced; the auxiliary mass eta is the standard AMFlow device.

assumptions (5)
  • domain assumption Dimensional regularization admits analytic continuation so that the DCT relation holds for any values of D0 and D.
    Invoked immediately after Eq. (4) in Section II.A without proof. This is the basis for using a non-integer value D0 = 7/13 in the examples.
  • domain assumption The auxiliary-mass integral I^{D0}(eta) has no poles in the fourth quadrant for the kinematics of interest.
    Needed to rotate the DCT contour from the positive real axis to the negative imaginary axis in Eq. (5) and Eq. (10). It is stated as a condition but not proved for general configurations.
  • domain assumption The differential equation system for fixed-branch integrals, Eq. (11), and the associated Gram-matrix construction are correct as taken from Ref. [40].
    The companion paper [40] is unpublished and its arXiv identifier is a placeholder. For B = 1 it reduces to the AMFlow one-loop equation, but for general FBIs it is imported rather than derived here.
  • domain assumption Master integrals of fixed-branch integrals can always be chosen as corner integrals, so the differential-equation system closes with z0 and C not both vanishing.
    Stated in Section III without proof. It is used to guarantee a solvable closed system of differential equations for general fixed-branch integrals.
  • domain assumption AMFlow correctly solves the boundary-value problems for the auxiliary-mass integrals at the chosen dimension D0.
    AMFlow is the numerical engine used to produce the input integrals and also the validation baseline. It is a separate, established method, but it is from the same group and is taken as correct.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Efficient Computation of One-Loop Feynman Integrals and Fixed-Branch Integrals to High Orders in $\epsilon$." pith.science (2026). https://pith.science/paper/C2JEK2Y3

@misc{pith2026241221054,
  author       = {Pith},
  title        = {Pith review of: Efficient Computation of One-Loop Feynman Integrals and Fixed-Branch Integrals to High Orders in $\epsilon$},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/C2JEK2Y3}},
  note         = {Machine review of arXiv:2412.21054}
}
abstract

We propose a novel method, called the dimension-changing transformation (DCT), to compute one-loop Feynman integrals and recently introduced fixed-branch integrals to arbitrary orders in $\epsilon$. The DCT relates one-loop Feynman integrals or fixed-branch integrals in one spacetime dimension to their corresponding quantities with auxiliary mass in any other dimension, making the expansion to high orders in $\epsilon$ highly efficient. We applied this method to several examples to demonstrate its validity and efficiency. The approach introduced in this work has been implemented in an open-source C++ package, available at \href{https://gitlab.com/multiloop-pku/dct}{https://gitlab.com/multiloop-pku/dct}.

Figures

Figures reproduced from arXiv: 2412.21054 by the authors.

Figure 1
Figure 1. FIG. 1. Different choices of integration paths for recover [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. A pentagon diagram in [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Time spent for computing a pentagon family with [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗
Figures from the paper (1 more)
Figure 5
Figure 5. Figure 5: FIG. 5. A two-loop diagram with three branches, red lines [PITH_FULL_IMAGE:figures/full_fig_p005_5.png]

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 3 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Recursive construction of scalar one-loop integrals in dimensional regularisation

    hep-th 2026-07 conditional novelty 8.0 of 10

    A recursion based on hyperbolic simplex volumes expresses every epsilon-expansion coefficient of scalar one-loop Feynman integrals in terms of multiple polylogarithms.

  2. Fast evaluation of Feynman integrals for Monte Carlo generators

    hep-ph 2025-07 conditional novelty 6.0 of 10

    A new numerical integrator evaluates one- and two-loop five-point Feynman integrals with complex masses in milliseconds, using variable-by-variable integration and a new branch-cut prescription.

  3. AMFlow 2.0: significant algorithmic and software improvements for Feynman integral evaluation

    hep-ph 2026-07 accept novelty 5.0 of 10

    AMFlow 2.0 cuts symbolic and numerical cost of multi-loop Feynman integral evaluation via an FT recursion mode, a C++ DE solver, and modern IBP reducers, demonstrated on a three-loop five-point family.

Reference graph

Works this paper leans on

44 extracted references · 10 canonical work pages · cited by 3 Pith papers

  1. [40]

    Huang, R.-J

    L.-H. Huang, R.-J. Huang, and Y.-Q. Ma, Tame multi-leg Feynman integrals beyond one loop , [arXiv:2412.xxxxx]

  2. [1]

    F. V. Tkachov, A theorem on analytical calculability of 4-loop renormalization group functions , Phys. Lett. B 100 (1981) 65–68 [ InSPIRE]

  3. [2]

    K. G. Chetyrkin and F. V. Tkachov, Integration by parts: The algorithm to calculate β-functions in 4 loops , Nucl. Phys. B 192 (1981) 159–204 [ InSPIRE]

  4. [3]

    Laporta, High-precision calculation of multiloop Feynman integrals by difference equations, Int

    S. Laporta, High-precision calculation of multiloop Feynman integrals by difference equations, Int. J. Mod. Phys. A 15 (2000) 5087–5159 [ hep-ph/0102033] [InSPIRE]

  5. [4]

    Gluza, K

    J. Gluza, K. Kajda, and D. A. Kosower, Towards a Basis for Planar Two-Loop Integrals , Phys. Rev. D 83 (2011) 045012 [ arXiv:1009.0472] [InSPIRE]

  6. [5]

    R. M. Schabinger, A New Algorithm For The Generation Of Unitarity-Compatible Integration By Parts Relations, JHEP 01 (2012) 077 [arXiv:1111.4220] [InSPIRE]

  7. [6]

    von Manteuffel and C

    A. von Manteuffel and C. Studerus, Reduze 2 - Distributed Feynman Integral Reduction, [arXiv:1201.4330] [InSPIRE]

  8. [7]

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

Show all 44 references
  1. [8]

    von Manteuffel and R

    A. von Manteuffel and R. M. Schabinger, A novel approach to integration by parts reduction , Phys. Lett. B 744 (2015) 101–104 [ arXiv:1406.4513] [InSPIRE]

  2. [9]

    K. J. Larsen and Y. Zhang, Integration-by-parts reductions from unitarity cuts and algebraic geometry , Phys. Rev. D 93 (2016) 041701 [ arXiv:1511.01071] [InSPIRE]

  3. [10]

    Peraro, Scattering amplitudes over finite fields and multivariate functional reconstruction, JHEP 12 (2016) 6 030 [arXiv:1608.01902] [InSPIRE]

    T. Peraro, Scattering amplitudes over finite fields and multivariate functional reconstruction, JHEP 12 (2016) 6 030 [arXiv:1608.01902] [InSPIRE]

  4. [11]

    Mastrolia and S

    P. Mastrolia and S. Mizera, Feynman Integrals and Intersection Theory, JHEP 02 (2019) 139 [arXiv:1810.03818] [InSPIRE]

  5. [12]

    Liu and Y.-Q

    X. Liu and Y.-Q. Ma, Determining arbitrary Feynman integrals by vacuum integrals , Phys. Rev. D 99 (2019) 071501 [arXiv:1801.10523] [InSPIRE]

  6. [13]

    X. Guan, X. Liu, and Y.-Q. Ma, Complete reduction of integrals in two-loop five-light-parton scattering amplitudes, Chin. Phys. C 44 (2020) 093106 [arXiv:1912.09294] [InSPIRE]

  7. [14]

    Klappert and F

    J. Klappert and F. Lange, Reconstructing rational functions with FireFly , Comput. Phys. Commun. 247 (2020) 106951 [ arXiv:1904.00009] [InSPIRE]

  8. [15]

    Peraro, FiniteFlow: multivariate functional reconstruction using finite fields and dataflow graphs , JHEP 07 (2019) 031 [ arXiv:1905.08019] [InSPIRE]

    T. Peraro, FiniteFlow: multivariate functional reconstruction using finite fields and dataflow graphs , JHEP 07 (2019) 031 [ arXiv:1905.08019] [InSPIRE]

  9. [16]

    Frellesvig, F

    H. Frellesvig, F. Gasparotto, S. Laporta, M. K. Mandal, P. Mastrolia, L. Mattiazzi, and S. Mizera, Decomposition of Feynman Integrals on the Maximal Cut by Intersection Numbers , JHEP 05 (2019) 153 [arXiv:1901.11510] [InSPIRE]

  10. [17]

    Y. Wang, Z. Li, and N. Ul Basat, Direct reduction of multiloop multiscale scattering amplitudes , Phys. Rev. D 101 (2020) 076023 [ arXiv:1901.09390] [InSPIRE]

  11. [18]

    A. V. Smirnov and F. S. Chukharev, FIRE6: Feynman Integral REduction with modular arithmetic , Comput. Phys. Commun. 247 (2020) 106877 [arXiv:1901.07808] [InSPIRE]

  12. [19]

    Klappert, F

    J. Klappert, F. Lange, P. Maierh¨ ofer, and J. Usovitsch, Integral reduction with Kira 2.0 and finite field methods , Comput. Phys. Commun. 266 (2021) 108024 [arXiv:2008.06494] [InSPIRE]

  13. [20]

    Boehm, M

    J. Boehm, M. Wittmann, Z. Wu, Y. Xu, and Y. Zhang, IBP reduction coefficients made simple , JHEP 12 (2020) 054 [ arXiv:2008.13194] [InSPIRE]

  14. [21]

    Heller and A

    M. Heller and A. von Manteuffel, MultivariateApart: Generalized partial fractions, Comput. Phys. Commun. 271 (2022) 108174 [ arXiv:2101.08283] [InSPIRE]

  15. [22]

    Bendle, J

    D. Bendle, J. Boehm, M. Heymann, R. Ma, M. Rahn, L. Ristau, M. Wittmann, Z. Wu, H. Xu, and Y. Zhang, pfd-parallel, a Singular/GPI-Space package for massively parallel multivariate partial fractioning , Comput. Phys. Commun. 294 (2024) 108942 [arXiv:2104.06866] [InSPIRE]

  16. [23]

    ’t Hooft and M

    G. ’t Hooft and M. J. G. Veltman, Scalar One Loop Integrals, Nucl. Phys. B 153 (1979) 365–401 [ InSPIRE]

  17. [24]

    Passarino and M

    G. Passarino and M. J. G. Veltman, One Loop Corrections for e+ e- Annihilation Into mu+ mu- in the Weinberg Model, Nucl. Phys. B 160 (1979) 151–207 [InSPIRE]

  18. [25]

    Z. Bern, L. J. Dixon, and D. A. Kosower, Dimensionally regulated pentagon integrals, Nucl. Phys. B 412 (1994) 751–816 [hep-ph/9306240] [InSPIRE]

  19. [26]

    G. J. van Oldenborgh and J. A. M. Vermaseren, New Algorithms for One Loop Integrals , Z. Phys. C 46 (1990) 425–438 [ InSPIRE]

  20. [27]

    Hahn and M

    T. Hahn and M. Perez-Victoria, Automatized one loop calculations in four-dimensions and D-dimensions , Comput. Phys. Commun. 118 (1999) 153–165 [hep-ph/9807565] [InSPIRE]

  21. [28]

    R. K. Ellis and G. Zanderighi, Scalar one-loop integrals for QCD , JHEP 02 (2008) 002 [ arXiv:0712.1851] [InSPIRE]

  22. [29]

    van Hameren, OneLOop: For the evaluation of one-loop scalar functions , Comput

    A. van Hameren, OneLOop: For the evaluation of one-loop scalar functions , Comput. Phys. Commun. 182 (2011) 2427–2438 [ arXiv:1007.4716] [InSPIRE]

  23. [30]

    Denner, S

    A. Denner, S. Dittmaier, and L. Hofer, Collier: a fortran-based Complex One-Loop LIbrary in Extended Regularizations, Comput. Phys. Commun. 212 (2017) 220–238 [arXiv:1604.06792] [InSPIRE]

  24. [31]

    Cullen, J

    G. Cullen, J. P. Guillet, G. Heinrich, T. Kleinschmidt, E. Pilon, T. Reiter, and M. Rodgers, Golem95C: A library for one-loop integrals with complex masses , Comput. Phys. Commun. 182 (2011) 2276–2284 [arXiv:1101.5595] [InSPIRE]

  25. [32]

    H. H. Patel, Package-X: A Mathematica package for the analytic calculation of one-loop integrals , Comput. Phys. Commun. 197 (2015) 276–290 [ arXiv:1503.01469] [InSPIRE]

  26. [33]

    Badger, M

    S. Badger, M. Becchetti, E. Chaubey, R. Marzucca, and F. Sarandrea, One-loop QCD helicity amplitudes for pp → ttj to O( ε2), JHEP 06 (2022) 066 [arXiv:2201.12188] [InSPIRE]

  27. [34]

    J. M. Henn, A. Matijaˇ si´ c, and J. Miczajka,One-loop hexagon integral to higher orders in the dimensional regulator, JHEP 01 (2023) 096 [ arXiv:2210.13505] [InSPIRE]

  28. [35]

    Buccioni, P

    F. Buccioni, P. A. Kreer, X. Liu, and L. Tancredi, One loop QCD corrections to gg → ttH at O(ϵ2), [arXiv:2312.10015] [InSPIRE]

  29. [36]

    Liu, Y.-Q

    X. Liu, Y.-Q. Ma, and C.-Y. Wang, A Systematic and Efficient Method to Compute Multi-loop Master Integrals, Phys. Lett. B 779 (2018) 353–357 [arXiv:1711.09572] [InSPIRE]

  30. [37]

    Liu and Y.-Q

    X. Liu and Y.-Q. Ma, Multiloop corrections for collider processes using auxiliary mass flow , Phys. Rev. D 105 (2022) L051503 [ arXiv:2107.01864] [InSPIRE]

  31. [38]

    Liu and Y.-Q

    Z.-F. Liu and Y.-Q. Ma, Determining Feynman Integrals with Only Input from Linear Algebra , Phys. Rev. Lett. 129 (2022) 222001 [ arXiv:2201.11637] [InSPIRE]

  32. [39]

    Liu and Y.-Q

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

  33. [41]

    https://www.boost.org/doc/libs/1_76_0/libs/ multiprecision/doc/html/index.html

    The Boost Development Team, Boost.Multiprecision Library, 2023 [ InSPIRE]. https://www.boost.org/doc/libs/1_76_0/libs/ multiprecision/doc/html/index.html. Accessed on 2024-06-13

  34. [42]

    https://gmplib.org/

    The GMP Development Team, GMP: The GNU Multiple Precision Arithmetic Library , 2023 [ InSPIRE]. https://gmplib.org/. Accessed on 2024-06-13

  35. [43]

    https://www.mpfr.org/

    The MPFR Development Team, MPFR: The GNU Multiple Precision Floating-Point Reliable Library , 2023 [InSPIRE]. https://www.mpfr.org/. Accessed on 2024-06-13

  36. [44]

    https://www.multiprecision.org/

    The MPC Development Team, MPC: The GNU Multiple Precision Complex Arithmetic Library , 2023 [InSPIRE]. https://www.multiprecision.org/. Accessed on 2024-06-13

Pith tools

Reviewed August 10, 2026 · model on record in the stance chip above.