Berends-Giele Recursion

• Berends-Giele Recursion is a recursive algorithm that constructs off-shell currents to efficiently compute tree-level scattering amplitudes in gauge and string theories.
• It systematically assembles currents via cubic and quartic vertices and extends to frameworks like bi-adjoint scalars and string disk integrals with higher-derivative corrections.
• Its efficient implementation with dynamic programming reduces computational overhead and underpins essential constructs such as the inverse KLT kernel and BCJ numerators.

The Berends-Giele (BG) recursion relation is a fundamental recursive algorithm for constructing off-shell gauge-theory and scalar currents, forming the basis of efficient computation of tree-level scattering amplitudes in a wide range of quantum field theories and string theories. Originally devised for gluon amplitudes in Yang-Mills theory, the BG recursion generalizes to bi-adjoint scalars, string disk integrals, theories with higher-derivative deformations, and applications in gravitational, celestial, and loop-level contexts.

1. Foundations and General Structure

The BG recursion was originally developed to construct color-ordered tree amplitudes by recursively assembling off-shell currents from smaller ones, systematically summing over all partitions of external legs. For a generic color-ordered theory, one defines an off-shell current for an ordered set $P = \{1,2,\dots,n\}$ as $J(P)$, built by sewing together lower-point currents via cubic (and where relevant, quartic) vertices.

The canonical BG recursion for pure Yang-Mills reads: $$J^\mu(P) = \frac{1}{s_P} \sum_{P=X\cup Y} \left[ J_\nu(X) J_\rho(Y) V^{\nu\rho\mu}(k_X, k_Y) \right] + \text{quartic terms}$$ with $s_P = (k_{p_1} + \cdots + k_{p_n})^2$ and $V^{\nu\rho\mu}$ the color-ordered three-gluon vertex. For scalar bi-adjoint $\varphi^3$ theory, the basic recursion simplifies to

$J(P) = \frac{i}{P^2} \sum_{P=A\cup B} J(A) J(B)$

where all subcurrents are constructed with the same rule recursively down to the single-leg base case $J(i) = 1$ (Du et al., 2022). The recursion generates all tree diagrams with a known ordering and tracks off-shell kinematics (Mafra et al., 2015, Macrelli et al., 2019, Garozzo et al., 2018).

2. Double-Color and String Theory Extensions

The BG recursion is adapted for double-color-ordered objects in bi-adjoint scalar theory, and for the $\alpha'$-expansion of open superstring disk integrals ("$Z$-theory"). The key building block in this context is the double current $J(P)$0, constructed recursively as

$J(P)$1

valid only if $J(P)$2 are permutations of the same label set (Mafra, 2016). For $J(P)$3-theory, the all-orders $J(P)$4 expansion of the full string disk integral $J(P)$5 is given as: $J(P)$6 The string-corrected recursion incorporates higher-derivative corrections,

$J(P)$7

with commutator-like sums and Riemann zeta values $J(P)$8 as coefficients of higher orders (Mafra et al., 2016). This recursion solves the non-abelian $J(P)$9-theory equation of motion and compactly generates all orders in $$J^\mu(P) = \frac{1}{s_P} \sum_{P=X\cup Y} \left[ J_\nu(X) J_\rho(Y) V^{\nu\rho\mu}(k_X, k_Y) \right] + \text{quartic terms}$$0.

3. Key Recursion Properties and Computational Realization

The recursion is constructed in terms of “words” (ordered multiplets of external labels) and relies fundamentally on operations of deconcatenation. The general scheme is implemented efficiently in symbolic and numerical codes. For practical computation, a dynamic programming approach caches previously computed subcurrents and incrementally constructs larger currents, reducing computational overhead from exponential to quartic scaling in leg number for QCD-like theories (Badger et al., 2012, Cruz-Martinez et al., 10 Feb 2025). For $$J^\mu(P) = \frac{1}{s_P} \sum_{P=X\cup Y} \left[ J_\nu(X) J_\rho(Y) V^{\nu\rho\mu}(k_X, k_Y) \right] + \text{quartic terms}$$1-theory, specific pseudocode prescribes an algorithm of precomputing kinematic invariants, initializing base cases, and recursively building higher-order contributions, employing shuffle and reflection symmetry constraints in the double-current basis (Mafra et al., 2016).

In bi-adjoint or $$J^\mu(P) = \frac{1}{s_P} \sum_{P=X\cup Y} \left[ J_\nu(X) J_\rho(Y) V^{\nu\rho\mu}(k_X, k_Y) \right] + \text{quartic terms}$$2-theory, the $$J^\mu(P) = \frac{1}{s_P} \sum_{P=X\cup Y} \left[ J_\nu(X) J_\rho(Y) V^{\nu\rho\mu}(k_X, k_Y) \right] + \text{quartic terms}$$3 recursion efficiently produces the inverse KLT kernel, required for the double-copy construction of gravity amplitudes (Mafra, 2016). The effective current expansion underlies the construction of BCJ numerators, crucial for manifesting color-kinematics duality (Wu et al., 2021, Garozzo et al., 2018).

4. Connection to Gravitational, Massive, and Deformed Theories

The BG recursion admits gravity analogues, both in gauge-theoretic and manifestly double-copy formulations. For pure gravity, the recursion for the all-plus or all-minus graviton off-shell current involves sums over spanning trees of the complete graph, aligning with the matrix-tree theorem and determinant representations: $$J^\mu(P) = \frac{1}{s_P} \sum_{P=X\cup Y} \left[ J_\nu(X) J_\rho(Y) V^{\nu\rho\mu}(k_X, k_Y) \right] + \text{quartic terms}$$4 Specializations of weights and cocycles recover the “half-soft” function and MHV amplitudes (Delfino et al., 2014, Krasnov et al., 2013, Hasuwannakit et al., 18 Jul 2025). For double-copy constructions, the classical KLT relations at the current level hold exactly in the “MHV” configurations (Cho et al., 2021).

In higher-derivative (e.g., $$J^\mu(P) = \frac{1}{s_P} \sum_{P=X\cup Y} \left[ J_\nu(X) J_\rho(Y) V^{\nu\rho\mu}(k_X, k_Y) \right] + \text{quartic terms}$$5 and $$J^\mu(P) = \frac{1}{s_P} \sum_{P=X\cup Y} \left[ J_\nu(X) J_\rho(Y) V^{\nu\rho\mu}(k_X, k_Y) \right] + \text{quartic terms}$$6) deformed Yang-Mills, or string-inspired effective actions, the currents are recursively deformed by insertions of extra tensor structures tied to string corrections, and pre-factors involving $$J^\mu(P) = \frac{1}{s_P} \sum_{P=X\cup Y} \left[ J_\nu(X) J_\rho(Y) V^{\nu\rho\mu}(k_X, k_Y) \right] + \text{quartic terms}$$7 and $$J^\mu(P) = \frac{1}{s_P} \sum_{P=X\cup Y} \left[ J_\nu(X) J_\rho(Y) V^{\nu\rho\mu}(k_X, k_Y) \right] + \text{quartic terms}$$8 (Garozzo et al., 2018, Mafra et al., 2016). The recursion remains of the same schematic type, but with more elaborate kinematic kernels and auxiliary currents. For massive higher-spin states (e.g., massive spin-2 states in open type-I superstring), the BG recursion assembles symmetric tensor currents bilinearly in lower-rank gluon field strengths, with cyclic sum structures for color-ordered amplitudes (Huang, 2024).

5. Loop-Level Generalizations and Algebraic Frameworks

The classical BG recursion extends to quantum and loop-level recursion relations via homological and algebraic frameworks. The homological perturbation lemma for $$J^\mu(P) = \frac{1}{s_P} \sum_{P=X\cup Y} \left[ J_\nu(X) J_\rho(Y) V^{\nu\rho\mu}(k_X, k_Y) \right] + \text{quartic terms}$$9- or $s_P = (k_{p_1} + \cdots + k_{p_n})^2$0-algebras organizes both the tree- and loop-level recursion as a minimal-model computation, with loop order governed by insertions of BV Laplacian terms. The resulting “quantum BG recursion” generates all quantum corrections recursively, with single-loop terms arising as quantum descendants grafted onto tree-level diagrams (Garg et al., 2024, Jurco et al., 2019).

More abstractly, the BG recursion is identified as the explicit realization of the minimal $s_P = (k_{p_1} + \cdots + k_{p_n})^2$1 or $s_P = (k_{p_1} + \cdots + k_{p_n})^2$2 model for the BV complex of a field theory, with each off-shell current corresponding to a sequence of homotopy brackets summing over cubic tree structures (Macrelli et al., 2019).

6. Graphical, Algebraic, and Celestial Extensions

Advanced formulations interpret the BG recursion in terms of graph-theoretic identities. Graph-based BCJ relations encode the recursion as combinatorial sums over chain and star graphs; the BG recursion naturally organizes the proof of these fundamental identities (Du et al., 2022, Mafra, 2016). The recursion admits a celestial generalization, where Mellin-transformed off-shell currents on the celestial sphere are constructed recursively, with the sewing operation adapted to the celestial context for loop integrands (Tao, 2023). The celestial BG recursion preserves features such as $s_P = (k_{p_1} + \cdots + k_{p_n})^2$3 OPE singularities and correctly captures the conformal structure in collinear limits.

7. Algorithmic and Computational Considerations

The BG recursion is a workhorse in amplitude computations for QCD and string theory. Efficient implementations employ dynamic programming, memoization of subcurrents, and bottom-up approaches for faster scaling, especially at high multiplicity (Badger et al., 2012, Cruz-Martinez et al., 10 Feb 2025). The methodology supports computations in arbitrary spacetime dimension, over finite fields and for massive and deformed kinematics, and is now realized in high-performance computing environments including GPU acceleration (Cruz-Martinez et al., 10 Feb 2025).

Analytical and numerical stability in the presence of soft and collinear singularities is preserved, and in QCD, the recursion organizes color and helicity sums efficiently, allowing automated amplitudes as building blocks for higher-loop integrands and cross-section calculations (Garg et al., 2024).

---

Table: BG Recursion Variants (Key Theories and Features)

Theory	BG Current Structure	Notable Features
Yang-Mills	$s_P = (k_{p_1} + \cdots + k_{p_n})^2$4; cubic and quartic	Fundamental cubic/ quartic recursions
Bi-adjoint scalar	$s_P = (k_{p_1} + \cdots + k_{p_n})^2$5 (single-color), $s_P = (k_{p_1} + \cdots + k_{p_n})^2$6 (double)	Double-color recursion, inverse KLT
Gravity	$s_P = (k_{p_1} + \cdots + k_{p_n})^2$7, $s_P = (k_{p_1} + \cdots + k_{p_n})^2$8	Tree-sum/determinant structure, half-soft function
$s_P = (k_{p_1} + \cdots + k_{p_n})^2$9-theory (string)	$V^{\nu\rho\mu}$0, string corrections	$V^{\nu\rho\mu}$1, $V^{\nu\rho\mu}$2-expanded recursion
Light-front QCD	Fragmentation functions $V^{\nu\rho\mu}$3	Adapted recursion in light-front variables
Quantum/loop	$V^{\nu\rho\mu}$4, $V^{\nu\rho\mu}$5	Loop-level recursion via homotopy algebra

• "Non-abelian $V^{\nu\rho\mu}$6-theory: Berends-Giele recursion for the $V^{\nu\rho\mu}$7-expansion of disk integrals" (Mafra et al., 2016)
• "Berends-Giele recursion for double-color-ordered amplitudes" (Mafra, 2016)
• "Note on graph-based BCJ relation for Berends-Giele currents" (Du et al., 2022)
• "Recursion for Differential Cross-Section from the Optical Theorem" (Garg et al., 2024)
• "Gravity MHV amplitudes via Berends-Giele currents" (Hasuwannakit et al., 18 Jul 2025)
• "Berends-Giele currents in Bern-Carrasco-Johansson gauge for $V^{\nu\rho\mu}$8- and $V^{\nu\rho\mu}$9-deformed Yang-Mills amplitudes" (Garozzo et al., 2018)
• "Scattering Amplitude Recursion Relations in BV Quantisable Theories" (Macrelli et al., 2019)
• "Accelerating Berends-Giele recursion for gluons in arbitrary dimensions over finite fields" (Cruz-Martinez et al., 10 Feb 2025)
• "Higher-Point Gauge-Theory Couplings of Massive Spin-2 States in 4-Dimensional String Theories" (Huang, 2024)
• "Pure connection formalism for gravity: Recursion relations" (Delfino et al., 2014)
• "Weighted Laplacians, cocycles and recursion relations" (Krasnov et al., 2013)
• "Celestial Berends-Giele current" (Tao, 2023)
• "Loop Amplitudes and Quantum Homotopy Algebras" (Jurco et al., 2019)
• "Berends-Giele recursions and the BCJ duality in superspace and components" (Mafra et al., 2015)