This content will become publicly available on January 2, 2027

Title: Analytic regression of Feynman integrals from high-precision numerical sampling
A<sc>bstract</sc> In mathematics or theoretical physics one is often interested in obtaining an exact analytic description of some data which can be produced, in principle, to arbitrary accuracy. For example, one might like to know the exact analytical form of a definite integral. Such problems are not well-suited to numerical symbolic regression, since typical numerical methods lead only to approximations. However, if one has some sense of the function space in which the analytic result should lie, it is possible to deduce the exact answer by judiciously sampling the data at a sufficient number of points with sufficient precision. We demonstrate how this can be done for the computation of Feynman integrals. We show that by combining high-precision numerical integration with analytic knowledge of the function space one can often deduce the exact answer using lattice reduction. A number of examples are given as well as an exploration of the trade-offs between number of datapoints, number of functional predicates, precision of the data, and compute. This method provides a bottom-up approach that neatly complements the top-down Landau-bootstrap approach of trying to constrain the exact answer using the analytic structure alone. Although we focus on the application to Feynman integrals, the techniques presented here are more general and could apply to a wide range of problems where an exact answer is needed and the function space is sufficiently well understood.  more » « less
Award ID(s):
2019786
PAR ID:
10664242
Author(s) / Creator(s):
; ; ; ;
Publisher / Repository:
Springer Nature
Date Published:
Journal Name:
Journal of High Energy Physics
Volume:
2026
Issue:
1
ISSN:
1029-8479
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. A<sc>bstract</sc> Two-loop multi-leg form factors in off-shell kinematics require knowledge of planar and nonplanar double box Feynman diagrams with massless internal propagators. These are complicated functions of Mandelstam variables and external particle virtualities. The latter serve as regulators of infrared divergences, thus making these observables finite in four space-time dimensions. In this paper, we use the method of canonical differential equations for the calculation of (non)planar double box integrals in the near mass-shell kinematical regime, i.e., where virtualities of external particles are much smaller than the Mandelstam variables involved. We deduce a basis of master integrals with uniform transcendental weight based on the analysis of leading singularities employing the Baikov representation as well as an array of complementary techniques. We dub the former asymptotically canonical since it is valid in the near mass-shell limit of interest. We iteratively solve resulting differential equations up to weight four in terms of multiple polylogarithms. 
    more » « less
  2. A<sc>bstract</sc> We use the embedding formalism to construct conformal fields inDdimensions, by restricting Lorentz-invariant ensembles of homogeneous neural networks in (D+ 2) dimensions to the projective null cone. Conformal correlators may be computed using the parameter space description of the neural network. Exact four-point correlators are computed in a number of examples, and we perform a 4D conformal block decomposition that elucidates the spectrum. In a non-unitary example the decomposition precisely matches OPE coefficients for the self-correlator, but not for the mixed correlator. In others, the analysis is facilitated by recent approaches to Feynman integrals. Generalized free CFTs are constructed using the infinite-width Gaussian process limit of the neural network, enabling a realization of the free boson. The extension to deep networks constructs conformal fields at each subsequent layer, with recursion relations relating their conformal dimensions and four-point functions. Numerical approaches are discussed. 
    more » « less
  3. null (Ed.)
    A bstract Precision studies of the Higgs boson at future e + e − colliders can help to shed light on fundamental questions related to electroweak symmetry breaking, baryogenesis, the hierarchy problem, and dark matter. The main production process, e + e − → HZ , will need to be controlled with sub-percent precision, which requires the inclusion of next-to-next-to-leading order (NNLO) electroweak corrections. The most challenging class of diagrams are planar and non-planar double-box topologies with multiple massive propagators in the loops. This article proposes a technique for computing these diagrams numerically, by transforming one of the sub-loops through the use of Feynman parameters and a dispersion relation, while standard one-loop formulae can be used for the other sub-loop. This approach can be extended to deal with tensor integrals. The resulting numerical integrals can be evaluated in minutes on a single CPU core, to achieve about 0.1% relative precision. 
    more » « less
  4. Martinsson, Gunnar (Ed.)
    Abstract Gauss–Legendre quadrature, Clenshaw–Curtis quadrature and the trapezoid rule are powerful tools for numerical integration of analytic functions. For nearly singular problems, however, these standard methods become unacceptably slow. We discuss and generalize some existing methods for improving on these schemes when the location of the nearby singularity is known. We conclude with an application to some nearly singular surface integrals that arise in three-dimensional viscous fluid flow. 
    more » « less
  5. A<sc>bstract</sc> We compute the tree-level connected six-point function of identical scalar fluctuations of the AdS2string worldsheet dual to the half-BPS Wilson line in planar$$ \mathcal{N} $$ N = 4 Super Yang-Mills. The calculation can be carried out analytically in the conformal gauge approach, where the boundary reparametrization mode of the string plays a crucial role. We also study the analytic continuation of the six-point function to an out-of-time-order configuration, which is related to a 3-to-3 scattering amplitude in flat space. As a check of our results, we also numerically compute the six-point function using the Nambu-Goto action in static gauge, finding agreement with the conformal gauge answer. 
    more » « less