Search

Multiple polylogarithms are equipped with rich algebraic structures including the motivic coaction and the single-valued map which both found fruitful applications in high-energy physics. In recent work arXiv:2312.00697, the current authors presented a conjectural reformulation of the motivic coaction and the single-valued map via zeta generators, certain operations on non-commuting variables in suitable generating series of multiple polylogarithms. In this work, the conjectures of the reference will be proven for multiple polylogarithms that depend on any number of variables on the Riemann sphere., Comment: 39 + 11 pages

Zeta generators are derivations associated with odd Riemann zeta values that act freely on the Lie algebra of the fundamental group of Riemann surfaces with marked points. The genus-zero incarnation of zeta generators are Ihara derivations of certain Lie polynomials in two generators that can be obtained from the Drinfeld associator. We characterize a canonical choice of these polynomials, together with their non-Lie counterparts at even degrees $w\geq 2$, through the action of the dual space of formal and motivic multizeta values. Based on these canonical polynomials, we propose a canonical isomorphism that maps motivic multizeta values into the $f$-alphabet. The canonical Lie polynomials from the genus-zero setup determine canonical zeta generators in genus one that act on the two generators of Enriquez' elliptic associators. Up to a single contribution at fixed degree, the zeta generators in genus one are systematically expanded in terms of Tsunogai's geometric derivations dual to holomorphic Eisenstein series, leading to a wealth of explicit high-order computations. Earlier ambiguities in defining the non-geometric part of genus-one zeta generators are resolved by imposing a new representation-theoretic condition. The tight interplay between zeta generators in genus zero and genus one unravelled in this work connects the construction of single-valued multiple polylogarithms on the sphere with iterated-Eisenstein-integral representations of modular graph forms., Comment: 92 pages. Submission includes ancillary data files. v2: Typos corrected

We study non-holomorphic modular forms built from iterated integrals of holomorphic modular forms for SL$(2,\mathbb Z)$ known as equivariant iterated Eisenstein integrals. A special subclass of them furnishes an equivalent description of the modular graph forms appearing in the low-energy expansion of string amplitudes at genus one. Notably the Fourier expansion of modular graph forms contains single-valued multiple zeta values. We deduce the appearance of products and higher-depth instances of multiple zeta values in equivariant iterated Eisenstein integrals, and ultimately modular graph forms, from the appearance of simpler odd Riemann zeta values. This analysis relies on so-called zeta generators which act on certain non-commutative variables in the generating series of the iterated integrals. From an extension of these non-commutative variables we incorporate iterated integrals involving holomorphic cusp forms into our setup and use them to construct the modular completion of triple Eisenstein integrals. Our work represents a fully explicit realisation of the modular graph forms within Brown's framework of equivariant iterated Eisenstein integrals and reveals structural analogies between single-valued period functions appearing in genus zero and one string amplitudes., Comment: 102 pages plus appendices; submission includes ancillary data files; v2: minor corrections, published version

We introduce a new Lie-algebraic approach to explicitly construct the motivic coaction and single-valued map of multiple polylogarithms in any number of variables. In both cases, the appearance of multiple zeta values is controlled by conjugating generating series of polylogarithms with Lie-algebra generators associated with odd zeta values. Our reformulation of earlier constructions of coactions and single-valued polylogarithms preserves choices of fibration bases, exposes the correlation between multiple zeta values of different depths and paves the way for generalizations beyond genus zero., Comment: 13 pages; v2: minor corrections and clarifications; matches published version

A wealth of information on multiloop string amplitudes is encoded in fermionic two-point functions known as Szeg\"o kernels. In this paper we show that cyclic products of any number of Szeg\"o kernels on a Riemann surface of arbitrary genus may be decomposed into linear combinations of modular tensors on moduli space that carry all the dependence on the spin structure $\delta$. The $\delta$-independent coefficients in these combinations carry all the dependence on the marked points and are composed of the integration kernels of higher-genus polylogarithms. We determine the antiholomorphic moduli derivatives of the $\delta$-dependent modular tensors., Comment: 5.5 + 1.5 pages; v2: version to be published in Physics Review Letters, merged with the supplemental material as appendices

An explicit construction is presented of homotopy-invariant iterated integrals on a Riemann surface of arbitrary genus in terms of a flat connection valued in a freely generated Lie algebra. The integration kernels consist of modular tensors, built from convolutions of the Arakelov Green function and its derivatives with holomorphic Abelian differentials, combined into a flat connection. Our construction thereby produces explicit formulas for polylogarithms as higher-genus modular tensors. This construction generalizes the elliptic polylogarithms of Brown-Levin, and prompts future investigations into the relation with the function spaces of higher-genus polylogarithms in the work of Enriquez-Zerbini., Comment: 54 pages, 2 figures; v2: references added, expanded the discussion of modular properties in sections 3 and 4

The summation over spin structures, which is required to implement the GSO projection in the RNS formulation of superstring theories, often presents a significant impediment to the explicit evaluation of superstring amplitudes. In this paper we discover that, for Riemann surfaces of genus two and even spin structures, a collection of novel identities leads to a dramatic simplification of the spin structure sum. Explicit formulas for an arbitrary number of vertex points are obtained in two steps. First, we show that the spin structure dependence of a cyclic product of Szeg\"o kernels (i.e. Dirac propagators for worldsheet fermions) may be reduced to the spin structure dependence of the four-point function. Of particular importance are certain trilinear relations that we shall define and prove. In a second step, the known expressions for the genus-two even spin structure measure are used to perform the remaining spin structure sums. The dependence of the spin summand on the vertex points is reduced to simple building blocks that can already be identified from the two-point function. The hyper-elliptic formulation of genus-two Riemann surfaces is used to derive these results, and its $SL(2,\mathbb C)$ covariance is employed to organize the calculations and the structure of the final formulas. The translation of these results into the language of Riemann $\vartheta$-functions, and applications to the evaluation of higher-point string amplitudes, are relegated to subsequent companion papers., Comment: 61 + 55 pages, v2: minor corrections, matches published version

The low-energy expansion of closed-string scattering amplitudes at genus one introduces infinite families of non-holomorphic modular forms called modular graph forms. Their differential and number-theoretic properties motivated Brown's alternative construction of non-holomorphic modular forms in the recent mathematics literature from so-called equivariant iterated Eisenstein integrals. In this work, we provide the first validations beyond depth one of Brown's conjecture that equivariant iterated Eisenstein integrals contain modular graph forms. Apart from a variety of examples at depth two and three, we spell out the systematics of the dictionary and make certain elements of Brown's construction fully explicit to all orders., Comment: 45 pages; submission includes ancillary data files; v2: typos corrected / minor improvements, matches published version

Elliptic modular graph forms (eMGFs) are non-holomorphic modular forms depending on a modular parameter $\tau$ of a torus and marked points $z$ thereon. Traditionally, eMGFs are constructed from nested lattice sums over the discrete momenta on the worldsheet torus in closed-string genus-one amplitudes. In this work, we develop methods to translate the lattice-sum realization of eMGFs into iterated integrals over modular parameters $\tau$ of the torus with particular focus on cases with one marked point. Such iterated-integral representations manifest algebraic and differential relations among eMGFs and their degeneration limit $\tau \rightarrow i\infty$. From a mathematical point of view, our results yield concrete realizations of single-valued elliptic polylogarithms at arbitrary depth in terms of meromorphic iterated integrals over modular forms and their complex conjugates. The basis dimensions of eMGFs at fixed modular and transcendental weights are derived from a simple counting of iterated integrals and a generalization of Tsunogai's derivation algebra., Comment: 115 + 35 pages

We present a new method for numerically computing generic multi-loop Feynman integrals. The method relies on an iterative application of Feynman's trick for combining two propagators. Each application of Feynman's trick introduces a simplified Feynman integral topology which depends on a Feynman parameter that should be integrated over. For each integral family, we set up a system of differential equations which we solve in terms of a piecewise collection of generalized series expansions in the Feynman parameter. These generalized series expansions can be efficiently integrated term by term, and segment by segment. This approach leads to a fully algorithmic method for computing Feynman integrals from differential equations, which does not require the manual determination of boundary conditions. Furthermore, the most complicated topology that appears in the method often has less master integrals than the original one. We illustrate the strength of our method with a five-point two-loop integral family.

Feynman diagrams constitute one of the essential ingredients for making precision predictions for collider experiments. Yet, while the simplest Feynman diagrams can be evaluated in terms of multiple polylogarithms -- whose properties as special functions are well understood -- more complex diagrams often involve integrals over complicated algebraic manifolds. Such diagrams already contribute at NNLO to the self-energy of the electron, $t \bar{t}$ production, $\gamma \gamma$ production, and Higgs decay, and appear at two loops in the planar limit of maximally supersymmetric Yang-Mills theory. This makes the study of these more complicated types of integrals of phenomenological as well as conceptual importance. In this white paper contribution to the Snowmass community planning exercise, we provide an overview of the state of research on Feynman diagrams that involve special functions beyond multiple polylogarithms, and highlight a number of research directions that constitute essential avenues for future investigation., Comment: 32+24 pages, 11 figures, contribution to Snowmass 2021

Modern particle physics is increasingly becoming a precision science that relies on advanced theoretical predictions for the analysis and interpretation of experimental results. The planned physics program at the LHC and future colliders will require three-loop electroweak and mixed electroweak-QCD corrections to single-particle production and decay processes and two-loop electroweak corrections to pair production processes, all of which are beyond the reach of existing analytical and numerical techniques in their current form. This article presents a new semi-numerical approach based on differential equations with boundary terms specified at Euclidean kinematic points. These Euclidean boundary terms can be computed numerically with high accuracy using sector decomposition or other numerical methods. They are then mapped to the physical kinematic configuration with a series solution of the differential equation system. The method is able to deliver 8 or more digits precision, and it has a built-in mechanism for checking the accuracy of the obtained results. Its efficacy is illustrated with examples for three-loop self-energy and vertex integrals and two-loop box integrals.

DiffExp is a Mathematica package for integrating families of Feynman integrals order-by-order in the dimensional regulator from their systems of differential equations, in terms of one-dimensional series expansions along lines in phase-space, which are truncated at a given order in the line parameter. DiffExp is based on the series expansion strategies that were explored in recent literature for the computation of families of Feynman integrals relevant for Higgs plus jet production with full heavy quark mass dependence at next-to-leading order. The main contribution of this paper, and its associated package, is to provide a public implementation of these series expansion methods, which works for any family of integrals for which the user provides a set of differential equations and boundary conditions (and for which the program is not computationally constrained.) The main functions of the DiffExp package are discussed, and its use is illustrated by applying it to the three loop equal-mass and unequal-mass banana graph families., Comment: Typos fixed, references added

In this paper we complete the computation of the two-loop master integrals relevant for Higgs plus one jet production initiated in arXiv:1609.06685, arXiv:1907.13156, arXiv:1907.13234. We compute the integrals by defining differential equations along contours in the kinematic space, and by solving them in terms of one-dimensional generalized power series. This method allows for the efficient evaluation of the integrals in all kinematic regions, with high numerical precision. We show the generality of our approach by considering both the top- and the bottom-quark contributions. This work along with arXiv:1609.06685, arXiv:1907.13156, arXiv:1907.13234 provides the full set of master integrals relevant for the NLO corrections to Higgs plus one jet production, and for the real-virtual contributions to the NNLO corrections to inclusive Higgs production in QCD in the full theory., Comment: 32 pages, references added, minor revision

We define linearly reducible elliptic Feynman integrals, and we show that they can be algorithmically solved up to arbitrary order of the dimensional regulator in terms of a 1-dimensional integral over a polylogarithmic integrand, which we call the inner polylogarithmic part (IPP). The solution is obtained by direct integration of the Feynman parametric representation. When the IPP depends on one elliptic curve (and no other algebraic functions), this class of Feynman integrals can be algorithmically solved in terms of elliptic multiple polylogarithms (eMPLs) by using integration by parts identities. We then elaborate on the differential equations method. Specifically, we show that the IPP can be mapped to a generalized integral topology satisfying a set of differential equations in $\epsilon$-form. In the examples we consider the canonical differential equations can be directly solved in terms of eMPLs up to arbitrary order of the dimensional regulator. The remaining 1-dimensional integral may be performed to express such integrals completely in terms of eMPLs. We apply these methods to solve two- and three-points integrals in terms of eMPLs. We analytically continue these integrals to the physical region by using their 1-dimensional integral representation., Comment: The differential equations method is applied to linearly reducible elliptic Feynman integrals, the solutions are in terms of elliptic polylogarithms, JHEP version, 50 pages