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 →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
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.
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 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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)
- [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.
- [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.
- [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)
- [III, after Eq. (10)] The sentence beginning 'since ∆ − δ is fixed' should have 'Since' capitalized at the start of the sentence.
- [IV, first example] The word 'diferential' should be spelled 'differential' in the description of solving the η differential equations.
- [References] Reference [40] is listed with a placeholder arXiv number 'arXiv:2412.xxxxx'; this should be updated once the companion paper is available.
- [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
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
assumptions (5)
- domain assumption Dimensional regularization admits analytic continuation so that the DCT relation holds for any values of D0 and D.
- domain assumption The auxiliary-mass integral I^{D0}(eta) has no poles in the fourth quadrant for the kinematics of interest.
- 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].
- 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.
- domain assumption AMFlow correctly solves the boundary-value problems for the auxiliary-mass integrals at the chosen dimension D0.
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
Forward citations
Cited by 3 Pith papers
-
Recursive construction of scalar one-loop integrals in dimensional regularisation
A recursion based on hyperbolic simplex volumes expresses every epsilon-expansion coefficient of scalar one-loop Feynman integrals in terms of multiple polylogarithms.
-
Fast evaluation of Feynman integrals for Monte Carlo generators
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.
-
AMFlow 2.0: significant algorithmic and software improvements for Feynman integral evaluation
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
-
[40]
L.-H. Huang, R.-J. Huang, and Y.-Q. Ma, Tame multi-leg Feynman integrals beyond one loop , [arXiv:2412.xxxxx]
-
[1]
F. V. Tkachov, A theorem on analytical calculability of 4-loop renormalization group functions , Phys. Lett. B 100 (1981) 65–68 [ InSPIRE]
work page 1981
-
[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]
work page 1981
-
[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]
arXiv 2000
- [4]
-
[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]
arXiv 2012
-
[6]
A. von Manteuffel and C. Studerus, Reduze 2 - Distributed Feynman Integral Reduction, [arXiv:1201.4330] [InSPIRE]
-
[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]
arXiv 2014
Show all 44 references
-
[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]
2015 arXiv
-
[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]
2016 arXiv
-
[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]
2016 arXiv
-
[11]
Mastrolia and S
P. Mastrolia and S. Mizera, Feynman Integrals and Intersection Theory, JHEP 02 (2019) 139 [arXiv:1810.03818] [InSPIRE]
2019 arXiv
-
[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]
2019 arXiv
-
[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]
2020 arXiv
-
[14]
Klappert and F
J. Klappert and F. Lange, Reconstructing rational functions with FireFly , Comput. Phys. Commun. 247 (2020) 106951 [ arXiv:1904.00009] [InSPIRE]
2020 arXiv
-
[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]
2019 arXiv
-
[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]
2019 arXiv
-
[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]
2020 arXiv
-
[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]
2020 arXiv
-
[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]
2021 arXiv
-
[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]
2020 arXiv
-
[21]
Heller and A
M. Heller and A. von Manteuffel, MultivariateApart: Generalized partial fractions, Comput. Phys. Commun. 271 (2022) 108174 [ arXiv:2101.08283] [InSPIRE]
2022 arXiv
-
[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]
2024 arXiv
-
[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]
1979
-
[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]
1979
-
[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]
1994 arXiv
-
[26]
G. J. van Oldenborgh and J. A. M. Vermaseren, New Algorithms for One Loop Integrals , Z. Phys. C 46 (1990) 425–438 [ InSPIRE]
1990
-
[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]
1999 arXiv
-
[28]
R. K. Ellis and G. Zanderighi, Scalar one-loop integrals for QCD , JHEP 02 (2008) 002 [ arXiv:0712.1851] [InSPIRE]
2008 arXiv
-
[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]
2011 arXiv
-
[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]
2017 arXiv
-
[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]
2011 arXiv
-
[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]
2015 arXiv
-
[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]
2022 arXiv
-
[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]
2023 arXiv
-
[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]
-
[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]
2018 arXiv
-
[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]
2022 arXiv
-
[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]
2022 arXiv
-
[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]
2023 arXiv
-
[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
2023
-
[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
2023
-
[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
2023
-
[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
2023
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.