Loop Order Research Articles

All-loop group-theory constraints for four-point amplitudes of SU(N), SO(N), and Sp(N) gauge theories

In the decomposition of gauge-theory amplitudes into kinematic and color factors, the color factors (at a given loop order L) span a proper subspace of the extended trace space (which consists of single and multiple traces of generators of the gauge group, graded by powers of N). Using an iterative process, we systematically construct the L-loop color space of four-point amplitudes of fields in the adjoint representation of SU(N), SO(N), or Sp(N). We define the null space as the orthogonal complement of the color space. For SU(N), we confirm the existence of four independent null vectors (for L ≥ 2) and for SO(N) and Sp(N), we establish the existence of seventeen independent null vectors (for L ≥ 5). Each null vector corresponds to a group-theory constraint on the color-ordered amplitudes of the gauge theory.

Hidden zeros for particle/string amplitudes and the unity of colored scalars, pions and gluons

Recent years have seen the emergence of a new understanding of scattering amplitudes in the simplest theory of colored scalar particles — the Tr(ϕ3) theory — based on combinatorial and geometric ideas in the kinematic space of scattering data. In this paper we report a surprise: far from the toy model it appears to be, the “stringy” Tr(ϕ3) amplitudes secretly contains the scattering amplitudes for pions, as well as non-supersymmetric gluons, in any number of dimensions. The amplitudes for the different theories are given by one and the same function, related by a simple shift of the kinematics. This discovery was spurred by another fundamental observation: the tree-level Tr(ϕ3) field theory amplitudes have a hidden pattern of zeros when a special set of non-planar Mandelstam invariants is set to zero. These zeros are not manifest in Feynman diagrams but are made obvious by the connection of these amplitudes to the new understanding of associahedra arising from “causal diamonds” in kinematic space. Furthermore, near these zeros, the amplitudes simplify, by factoring into a non-trivial product of smaller amplitudes. Remarkably the amplitudes for pions and gluons are observed to also vanish in the same kinematical locus. These properties for Tr(ϕ3) amplitudes hold and further generalize to the “stringy” Tr(ϕ3) amplitudes. The “kinematic causal diamond” picture suggests a unique shift of the kinematic data that preserves the zeros, and this shift is precisely the one that unifies colored scalars, pions, and gluons into a single object. We will focus in this paper on explaining the hidden zeros and factorization properties and the connection between all the colored theories, working for simplicity at tree level. Subsequent works will describe this new formulation for the Non-linear Sigma Model and non-supersymmetric Yang-Mills theory, at all loop orders.

Nonlinear Sigma model amplitudes to all loop orders are contained in the Tr(Φ3) theory

Scattering amplitudes for the simplest theory of colored scalar particles—the Tr(Φ3) theory—have recently been the subject of active investigations. In this work we describe an unanticipated wider implication of this work: the Tr(Φ3) theory secretly contains nonlinear sigma model (NLSM) amplitudes to all loop orders. The NLSM amplitudes are obtained from Tr(Φ3) amplitudes by a unique shift of kinematic variables. We show that this shifted kinematics produces amplitudes for a cubic theory with a linear term in the potential, with extrema spontaneously breaking U(N)→U(N−k)×U(k). The Goldstone amplitudes for this theory coincide with those of pions in the U(N)×U(N)→U(N) chiral Lagrangian to all orders in the planar limit. We also give a purely on-shell understanding of this correspondence, showing that integrands defined by the kinematic shifts have the correct residues on poles and appropriately produce the Adler zero. Finally, we discuss how similar kinematic shifts produce certain infinite classes of mixed amplitudes of pions and Tr(Φ3) scalars, most of which are not interpretable from the Lagrangian description. Published by the American Physical Society 2024

KeldyshQFT: A C++ codebase for real-frequency multiloop functional renormalization group and parquet computations of the single-impurity Anderson model.

We provide a detailed exposition of our computational framework designed for the accurate calculation of real-frequency dynamical correlation functions of the single-impurity Anderson model in the regime of weak to intermediate coupling. Using quantum field theory within the Keldysh formalism to directly access the self-energy and dynamical susceptibilities in real frequencies, as detailed in our recent publication [Ge et al., Phys. Rev. B 109, 115128 (2024)], the primary computational challenge is the full three-dimensional real-frequency dependence of the four-point vertex. Our codebase provides a fully MPI+OpenMP parallelized implementation of the functional renormalization group (fRG) and the self-consistent parquet equations within the parquet approximation. It leverages vectorization to handle the additional complexity imposed by the Keldysh formalism, using optimized data structures and highly performant integration routines. Going beyond the results shown in the previous publication, the code includes functionality to perform fRG calculations in the multiloop framework, up to arbitrary loop order, including self-consistent self-energy iterations. Moreover, implementations of various regulators, such as hybridization, interaction, frequency, and temperature, are supplied.

Two-loop anomalous dimensions for small-R jet versus hadronic fragmentation functions

We study the collinear fragmentation of highly energetic jets defined with a small jet radius. In particular, we investigate how the corresponding fragmentation functions differ from their hadronic counterpart defined in the common MS¯\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\overline{\ extrm{MS}} $$\\end{document} scheme. We find that the anomalous dimensions governing the perturbative evolution of the two fragmentation functions differ starting at the two loop order. We compute for the first time the new anomalous dimensions at two loops and confirm our predictions by comparing the inclusive small-R jet spectrum against a fixed order perturbative calculation at O\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{O} $$\\end{document}(αs2\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ {\\alpha}_s^2 $$\\end{document}). To investigate the dependence of the anomalous dimension on the kinematic cutoff variable, we study the fragmentation functions of Cambridge jets defined with a transverse momentum cutoff as opposed to an angular cutoff R. We further study the evolution of the small-R fragmentation function with an alternative cutoff scale, proportional to zR, representing the maximum possible transverse momentum of emissions within a jet. In these cases we find that the two-loop anomalous dimensions coincide with the MS¯\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\overline{\ extrm{MS}} $$\\end{document} DGLAP ones, highlighting a correspondence between the MS¯\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\overline{\ extrm{MS}} $$\\end{document} scheme and a transverse-momentum cutoff.

Rational terms of UV origin to all loop orders

Numerical approaches to computations typically reconstruct the numerators of Feynman diagrams in four dimensions. In doing so, certain rational terms arising from the (D − 4)-dimensional part of the numerator multiplying ultraviolet (UV) poles in dimensional regularisation are not captured and need to be obtained by other means. At one-loop these rational terms of UV origin can be computed from a set of process-independent Feynman rules. Recently, it was shown that this approach can be extended to two loops. In this paper, we show that to all loop orders it is possible to compute rational terms of UV origin through process-independent vertices that are polynomial in masses and momenta.

Nonplanar integrated correlator in N=4 SYM

The integrated correlator of four superconformal stress-tensor primaries in SU(N) N=4 super Yang-Mills (SYM) theory in the perturbative limit takes a remarkably simple form, where the L-loop coefficient is given by a rational multiple of ζ(2L+1). In this paper, we extend the previous analysis of expressing the perturbative integrated correlator as a linear combination of periods of f-graphs, graphical representations for loop integrands, to the nonplanar sector at four loops. At this loop order, multiple zeta values make their first appearance when evaluating periods of nonplanar f-graphs, but cancel nontrivially in the weighted sum. The remaining single zeta value, along with the rational number prefactor, makes a perfect agreement with the prediction from supersymmmetric localization. Published by the American Physical Society 2024

Geometry from integrability: multi-leg fishnet integrals in two dimensions

We generalise the geometric analysis of square fishnet integrals in two dimensions to the case of hexagonal fishnets with three-point vertices. Our results support the conjecture that fishnet Feynman integrals in two dimensions, together with their associated geometry, are completely fixed by their Yangian and permutation symmetries. As a new feature for the hexagonal fishnets, the star-triangle identity introduces an ambiguity in the graph representation of a given Feynman integral. This translates into a map between different geometric interpretations attached to a graph. We demonstrate explicitly how these fishnet integrals can be understood as Calabi-Yau varieties, whose Picard-Fuchs ideals are generated by the Yangian over the conformal algebra. In analogy to elliptic curves, which represent the simplest examples of fishnet integrals with four-point vertices, we find that the simplest examples of three-point fishnets correspond to Picard curves with natural generalisations at higher loop orders.

On amplitudes and field redefinitions

We derive an off-shell recursion relation for correlators that holds at all loop orders. This allows us to prove how generalized amplitudes transform under generic field redefinitions, starting from an assumed behavior of the one-particle-irreducible effective action. The form of the recursion relation resembles the operation of raising the rank of a tensor by acting with a covariant derivative. This inspires a geometric interpretation, whose features and flaws we investigate.

Calabi-Yau periods for black hole scattering in classical general relativity

The high-precision description of black hole scattering in classical general relativity using the post-Minkowskian (PM) expansion requires the evaluation of single-scale Feynman integrals at increasing loop orders. Up to 4PM, the scattering angle and the impulse are expressible in terms of polylogarithmic functions and Calabi-Yau (CY) twofold periods. As in QFT, periods of higher dimensional CY n-folds are expected at higher PM order. We find at 5PM in the dissipative leading order self-force sector (5PM-1SF) that the only nonpolylogarithmic functions are the K3 periods encountered before and the ones of a new hypergeometric CY threefold. In the 5PM-2SF sector further CY twofold and threefold periods appear. Griffiths transversality of the CY period motives allows one to transform the differential equations for the master integrals into ε-factorized form and to solve them in terms of a well-controlled function space, as we demonstrate in the 5PM-1SF sector. Published by the American Physical Society 2024

Loop amplitudes in the Coulomb branch of N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 4 super-Yang-Mills theory

We study four point planar loop amplitudes at an arbitrary point in the Coulomb branch of N\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\mathcal{N} $$\\end{document} = 4 super-Yang-Mills theory. We study two particle unitary cuts up to four loop order. We explicitly verify that bubble and triangle graphs do not contribute at one loop level and show that the results hold at higher loop level as well. We also write down an all loop recursion relation for two particle reducible graphs for four point amplitudes.

Integration-by-parts identities and differential equations for parametrised Feynman integrals

Integration-by-parts (IBP) identities and differential equations are the primary modern tools for the evaluation of high-order Feynman integrals. They are commonly derived and implemented in the momentum-space representation. We provide a different viewpoint on these important tools by working in Feynman-parameter space, and using its projective geometry. Our work is based upon little-known results pre-dating the modern era of loop calculations [16–19, 30, 31]: we adapt and generalise these results, deriving a very general expression for sets of IBP identities in parameter space, associated with a generic Feynman diagram, and valid to any loop order, relying on the characterisation of Feynman-parameter integrands as projective forms. We validate our method by deriving and solving systems of differential equations for several simple diagrams at one and two loops, providing a unified perspective on a number of existing results.

EFT matching from analyticity and unitarity

We present a new on-shell method for the matching of ultraviolet models featuring massive states onto their massless effective field theory. We employ a dispersion relation in the space of complex momentum dilations to capture, in a single variable, the relevant analytic structure of scattering amplitudes at any multiplicity. Multivariate complex analysis and crossing considerations are therefore avoided. Remarkably, no knowledge about the infrared effective field theory is required in dimensional regularisation. All matching information is extracted from the residues and iscontinuities of the ultraviolet scattering amplitudes, which unitarity expresses in terms of lower-point and lower-loop results, respectively. This decomposition into simpler building blocks could deliver new insight in the structure of the effective field theories obtained from classes of ultraviolet scenarios and facilitate computations at higher loop orders.

Double excitations in the AdS(5)/CFT(4) integrable system and the Lagrange operator

It is argued that the integrable model for the planar spectrum of the AdS/CFT correspondence can accommodate for the full spectrum of excitations Dαα˙,ϕ[IJ],ψI,ψ¯I,Fαβ,F˜α˙β˙ (with I,J∈1…4) if double excitations are allowed for all three lowering operators of the internal SU(4) group as well as the supersymmetries. We present a tree-level analysis of related creation amplitudes in the nested Bethe ansatz as well as in the original level-1 picture in which excitations of various flavours scatter by a true S-matrix. In the latter case, the creation amplitudes for all double excitations we encounter take a perfectly universal form.Building on these ideas we work out Bethe solutions and states relevant in the mixing problem concerning the on-shell Lagrangian of N=4 super Yang-Mills theory. Owing to the very existence of double excitations, the chiral Yang-Mills field strength tensor can be represented by the four fermions {ψ31,ψ32,ψ41,ψ42} moving on a spin chain of length two. Our analysis remains restricted to leading order in the coupling, where the conformal eigenstate corresponding to the on-shell Lagrangian only comprises the pure Yang-Mills action. It should eventually be possible to augment our analysis to higher loop orders by incorporating coupling corrections in the relevant ingredients from the Bethe ansatz.Finally, it was recently realised how structure constants for operators containing the hitherto hidden half of the excitations can be computed by the hexagon formalism. We use this for a first test of our conjecture for the on-shell Lagrangian, namely that its three-point function with two half-BPS operators of equal length ought to vanish.

Towards the full Heisenberg-Euler effective action at large N

We study the Heisenberg-Euler effective action in constant electromagnetic fields F¯\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ \\overline{F} $$\\end{document} for QED with N charged particle flavors of the same mass and charge e in the large N limit characterized by sending N → ∞ while keeping Ne2 ∼ eF¯\\documentclass[12pt]{minimal} \\usepackage{amsmath} \\usepackage{wasysym} \\usepackage{amsfonts} \\usepackage{amssymb} \\usepackage{amsbsy} \\usepackage{mathrsfs} \\usepackage{upgreek} \\setlength{\\oddsidemargin}{-69pt} \\begin{document}$$ e\\overline{F} $$\\end{document} ∼ N0 fixed. This immediately implies that contributions that scale with inverse powers of N can be neglected and the resulting effective action scales linearly with N . Interestingly, due to the presence of one-particle reducible diagrams, even in this limit the Heisenberg-Euler effective action receives contributions of arbitrary loop order. In particular for the special cases of electric- and magnetic-like field configurations we construct an explicit expression for the associated effective Lagrangian that, upon extremization for two constant scalar coefficients, allows to evaluate its full, all-order result at arbitrarily large field strengths. We demonstrate that our manifestly nonperturbative expression correctly reproduces the known results for the Heisenberg-Euler effective action at large N , namely its all-loop strong field limit and its low-order perturbative expansion in powers of the fine-structure constant.

Full three-loop renormalisation of an abelian chiral gauge theory with non-anticommuting γ5 in the BMHV scheme

In this work we present a complete three-loop renormalisation of an abelian chiral gauge theory within the Breitenlohner-Maison/’t Hooft-Veltman (BMHV) scheme of dimensional regularisation (DReg). In this scheme the γ5-matrix appearing in gauge interactions is a non-anticommuting object, leading to a breaking of gauge and BRST invariance. Employing an efficient method based on the quantum action principle, we obtain the complete three-loop counterterm action which serves both to render the theory finite and to restore gauge and BRST invariance. The UV singular counterterms involve not only higher order ϵ-poles but also new counterterm structures emerging at the three-loop level for the first time; the finite symmetry-restoring counterterms are restricted to the same structures as at lower loop orders, just with different coefficients, aligning with our expectations. Both the singular and the finite counterterms include structures which cannot be obtained by the standard multiplicative renormalisation. Our results demonstrate that a rigorous treatment of chiral gauge theories with γ5 defined in the BMHV scheme at the multi-loop level is possible and that the obtained counterterm action is suitable for computer implementations, allowing automated calculations without ambiguities caused by γ5.

No- π schemes for multicoupling theories

We show that even ζ-functions may be removed from the β-functions of general multicoupling theories up to high loop order by means of coupling redefinitions. For theories whose β-function is determined by the anomalous dimensions of the fields, such as supersymmetric theories, this corresponds to a renormalization scheme change to a momentum subtraction scheme. Published by the American Physical Society 2024

Antisymmetric tensor field and Cheshire Cat smile of the local conformal symmetry

The conformal version of the antisymmetric second-order tensor field in four spacetime dimensions does not have gauge invariance extensively discussed in the literature for more than half a century. Our first observation is that, when coupled to fermions, only the conformal version provides renormalizability of the theory at the one-loop level. General considerations are supported by the derivation of one-loop divergences in the fermionic sector, indicating good chances for asymptotic freedom. The arguments concerning one-loop renormalizability remain valid in the presence of self-interactions and the masses for both fermion and antisymmetric tensor fields. In the flat spacetime limit, even regardless the conformal symmetry has gone, there is an expectation to meet renormalizability in all loop orders.

Nonanalyticity and on-shell factorization of inflation correlators at all loop orders

The dynamics of quantum fields during cosmic inflation can be probed via their late-time boundary correlators. The analytic structure of these boundary correlators contains rich physical information of bulk dynamics, and is also closely related to cosmological collider observables. In this work, we study a particular type of nonanalytic behavior, called nonlocal signals, for inflation correlators with massive exchanges at arbitrary loop orders. We propose a signal-detection algorithm to identify all possible sources of nonlocal signals in an arbitrary loop graph, and prove that the algorithm is exhaustive. We then present several versions of the on-shell factorization theorem for the leading nonlocal signal in graphs with arbitrary number of loops, and provide the explicit analytical expression for the leading nonlocal signal. We also generalize the nonlocal-signal cutting rule to arbitrary loop graphs. Finally, we provide many explicit examples to demonstrate the use of our results, including an n-loop melon graph and a variety of 2-loop graphs.

Nicolai maps with four-fermion interactions

Nicolai maps offer an alternative description of supersymmetric theories via nonlinear and nonlocal transformations characterized by the so-called ‘free-action’ and ‘determinant-matching’ conditions. The latter expresses the equality of the Jacobian determinant of the transformation with the one obtained by integrating out the fermions, which so far have been considered only to quadratic terms. We argue that such a restriction is not substantial, as Nicolai maps can be constructed for arbitrary nonlinear sigma models, which feature four-fermion interactions. The fermionic effective one-loop action then gets generalized to higher loops and the perturbative tree expansion of such Nicolai maps receives quantum corrections in the form of fermion loop decorations. The ‘free-action condition’ continues to hold for the classical map, but the ‘determinant-matching condition’ is extended to an infinite hierarchy in fermion loop order. After general considerations for sigma models in four dimensions, we specialize to the case of ℂPN symmetric spaces and construct the associated Nicolai map. These sigma models admit a formulation with only quadratic fermions via an auxiliary vector field, which does not simplify our analysis.