Explore the vast range of titles available.

A bstract We explore inequality constraints as a new tool for numerically evaluating Feynman integrals. A convergent Feynman integral is non-negative if the integrand is non-negative in either loop momentum space or Feynman parameter space. Applying various identities, all such integrals can be reduced to linear sums of a small set of master integrals, leading to infinitely many linear constraints on the values of the master integrals. The constraints can be solved as a semidefinite programming problem in mathematical optimization, producing rigorous two-sided bounds for the integrals which are observed to converge rapidly as more constraints are included, enabling high-precision determination of the integrals. Positivity constraints can also be formulated for the ϵ expansion terms in dimensional regularization and reveal hidden consistency relations between terms at different orders in ϵ . We introduce the main methods using one-loop bubble integrals, then present a nontrivial example of three-loop banana integrals with unequal masses, where 11 top-level master integrals are evaluated to high precision.

A bstract We reformulate differential equations (DEs) for Feynman integrals to avoid doubled propagators in intermediate steps. External momentum derivatives are dressed with loop momentum derivatives to form tangent vectors to unitarity cut surfaces, in a way inspired by unitarity-compatible IBP reduction. For the one-loop box, our method directly produces the final DEs without any integration-by-parts reduction. We further illustrate the method by deriving maximal-cut level differential equations for two-loop nonplanar five-point integrals, whose exact expressions are yet unknown. We speed up the computation using finite field techniques and rational function reconstruction.

A bstract We compute classical gravitational observables for the scattering of two spinless black holes in general relativity and N =8 supergravity in the formalism of Kosower, Maybee, and O’Connell (KMOC). We focus on the gravitational impulse with radiation reaction and the radiated momentum in black hole scattering at O ( G 3 ) to all orders in the velocity. These classical observables require the construction and evaluation of certain loop-level quantities which are greatly simplified by harnessing recent advances from scattering amplitudes and collider physics. In particular, we make use of generalized unitarity to construct the relevant loop integrands, employ reverse unitarity, the method of regions, integration-by-parts (IBP), and (canonical) differential equations to simplify and evaluate all loop and phase-space integrals to obtain the classical gravitational observables of interest to two-loop order. The KMOC formalism naturally incorporates radiation effects which enables us to explore these classical quantities beyond the conservative two-body dynamics. From the impulse and the radiated momentum, we extract the scattering angle and the radiated energy. Finally, we discuss universality of the impulse in the high-energy limit and the relation to the eikonal phase.

A bstract We use N = 8 supergravity as a toy model for understanding the dynamics of black hole binary systems via the scattering amplitudes approach. We compute the conservative part of the classical scattering angle of two extremal (half-BPS) black holes with minimal charge misalignment at O ( G 3 ) using the eikonal approximation and effective field theory, finding agreement between both methods. We construct the massive loop integrands by Kaluza-Klein reduction of the known D -dimensional massless integrands. To carry out integration we formulate a novel method for calculating the post-Minkowskian expansion with exact velocity dependence, by solving velocity differential equations for the Feynman integrals subject to modified boundary conditions that isolate conservative contributions from the potential region. Motivated by a recent result for universality in massless scattering, we compare the scattering angle to the result found by Bern et. al. in Einstein gravity and find that they coincide in the high-energy limit, suggesting graviton dominance at this order.

A bstract We present the computation of a full set of planar five-point two-loop master integrals with one external mass. These integrals are an important ingredient for two-loop scattering amplitudes for two-jet-associated W-boson production at leading color in QCD. We provide a set of pure integrals together with differential equations in canonical form. We obtain analytic differential equations efficiently from numerical samples over finite fields, fitting an ansatz built from symbol letters. The symbol alphabet itself is constructed from cut differential equations and we find that it can be written in a remarkably compact form. We comment on the analytic properties of the integrals and confirm the extended Steinmann relations, which govern the double discontinuities of Feynman integrals, to all orders in ϵ . We solve the differential equations in terms of generalized power series on single-parameter contours in the space of Mandelstam invariants. This form of the solution trivializes the analytic continuation and the integrals can be evaluated in all kinematic regions with arbitrary numerical precision.

A bstract We describe a systematic framework for computing the conservative potential of a compact binary system using modern tools from scattering amplitudes and effective field theory. Our approach combines methods for integration and matching adapted from effective field theory, generalized unitarity, and the double-copy construction, which relates gravity integrands to simpler gauge-theory expressions. With these methods we derive the third post-Minkowskian correction to the conservative two-body Hamiltonian for spinless black holes. We describe in some detail various checks of our integration methods and the resulting Hamiltonian.

A bstract We compute ϵ-factorized differential equations for all dimensionally-regularized integrals of the nonplanar hexa-box topology, which contribute for instance to 2-loop 5-point QCD amplitudes. A full set of pure integrals is presented. For 5-point planar topologies, Gram determinants which vanish in 4 dimensions are used to build compact expressions for pure integrals. Using unitarity cuts and computational algebraic geometry, we obtain a compact IBP system which can be solved in 8 hours on a single CPU core, overcoming a major bottleneck for deriving the differential equations. Alternatively, assuming prior knowledge of the alphabet of the nonplanar hexa-box, we reconstruct analytic differential equations from 30 numerical phase-space points, making the computation almost trivial with current techniques. We solve the differential equations to obtain the values of the master integrals at the symbol level. Full results for the differential equations and solutions are included as supplementary material.

A bstract We study the singularity structure of two-loop QED amplitudes for the production of multiple off-shell photons in massless electron-positron annihilation and develop counterterms that remove their infrared and ultraviolet divergences point by point in the loop integrand. The remainders of the subtraction are integrable in four dimensions and can be computed in the future with numerical integration. The counterterms capture the divergences of the amplitudes and factorize in terms of the Born amplitude and the finite remainder of the one-loop amplitude. They consist of simple one- and two-loop integrals with at most three external momenta and can be integrated analytically in a simple manner with established methods. We uncover novel aspects of fully local IR factorization, where vertex and self energy subdiagrams must be modified by new symmetrizations over loop momenta, in order to expose their tree-like tensor structures and hence factorization of IR singularities prior to loop integration. This work is a first step towards isolating locally the hard contributions of generic gauge theory amplitudes and rendering them integrable in exactly four dimensions with numerical methods.

A bstract Quantum Electrodynamics (QED) serves as a useful toy model for classical observables in gravitational two-body systems with reduced complexity due to the linearity of QED. We investigate scattering observables in scalar QED at the sixth order in the charges (two-loop order) in a classical regime analogous to the post-Minkowskian expansion in General Relativity. We employ modern scattering amplitude tools and extract classical observables by both eikonal methods and the formalism of Kosower, Maybee, and O’Connell (KMOC). In addition, we provide a simplified approach to extracting the radial action beyond the conservative sector.

A bstract We introduce a constructive method for defining a global loop-integrand basis for scattering amplitudes, encompassing both planar and nonplanar contributions. Our approach utilizes a graph-based framework to establish a well-defined, non-redundant basis of integrands. This basis, constructed from a chosen set of non-redundant graphs together with a selection of irreducible scalar products, provides clear insights into various physical properties of scattering amplitudes and proves useful in multiple contexts, such as on-shell Ward identities and manifesting gauge-choice independence. A key advantage of our integrand basis is its ability to streamline the generalized unitarity method. Specifically, we can directly read off the coefficients of basis elements without resorting to ansätze or solving linear equations. This novel approach allows us to lift generalized unitarity cuts — expressed as products of tree amplitudes — to loop-level integrands, facilitating the use of the tree-level double copy to generate complete gravitational integrands at any loop order. This method circumvents the difficulties in identifying complete higher-loop-order gauge-theory integrands that adhere to the color-kinematics duality. Additionally, our cut-based organization is well-suited for expansion in hard or soft limits, aiding in the exploration of ultraviolet or classical limits of scattering amplitudes.