Amplitude relations in heterotic string theory and Einstein-Yang-Mills

We present all-multiplicity evidence that the tree-level S-matrix of gluons and gravitons in heterotic string theory can be reduced to color-ordered single-trace amplitudes of the gauge multiplet. Explicit amplitude relations are derived for up to three gravitons, up to two color traces and an arbitrary number of gluons in each case. The results are valid to all orders in the inverse string tension α′ and generalize to the ten-dimensional superamplitudes which preserve 16 supercharges. Their field-theory limit results in an alternative proof of the recently discovered relations between Einstein-Yang-Mills amplitudes and those of pure Yang-Mills theory. Similarities and differences between the integrands of the Cachazo-He-Yuan formulae and the heterotic string are investigated.


Introduction
Modern results on scattering amplitudes have revealed a fascinating interplay between perturbative gauge interactions and gravity. A steadily growing list of examples suggests that gravitational S-matrices can be assembled from squares of suitably-arranged building blocks of amplitudes in Yang-Mills (YM) gauge theory. This double-copy structure is invisible to the field equations or a Lagrangian description but emerges naturally from string theory. Given that the massless vibration modes of open and closed strings include gauge bosons and gravitons, respectively, the double-copy construction of gravity ties in with the geometric intuition that closed strings can be formed by joining the endpoints of two open strings. For tree-level amplitudes, the Kawai-Lewellen-Tye (KLT) relations [1] provided a precise prescription for squaring open-string amplitudes to those of the closed string. In the limit of vanishing string length, this immediately propagates to tree amplitudes in perturbative gauge theory and gravity.
While a higher-genus analogue of the KLT relations is currently out of reach in string theory, a conjectural loop-level incarnation is available in field theory [2][3][4]. The backbone of this double-copy construction is the duality between color and kinematics due to Bern, Carrasco and Johansson (BCJ) which organizes gauge-theory amplitudes in terms of cubic graphs and treats color and kinematic degrees of freedom on equal footing. At tree level, this duality implies the BCJ relations between color-ordered gauge-theory amplitudes [2] which were for instance elegantly proven from string theory [5,6]. At loop level, the BCJ JHEP11(2016)074 duality has set new standards for multiloop calculations in various supergravities, see for instance [7][8][9][10]. Explicit realizations of the kinematic constituents which obey the BCJ duality are challenging to find in general, and a variety of examples at tree level and loop level has been generated via string-theory techniques [11][12][13][14][15][16][17].
Less is known about mixed scattering amplitudes involving both gravitons and gauge bosons, starting with the results of [18][19][20] from the late 1990's. More recently, a doublecopy description of gluon-graviton amplitudes in supergravities has been established in [21,22], see [23] for a review. Also here, string-theory methods provide a successful line of attack: elegant amplitude relations for mixed open-and closed-string scattering in the type-I superstring were extracted from monodromy properties of the disk worldsheet [6,24,25], valid to all orders in the inverse string tension α . In the field-theory limit α → 0 of this setup, tree amplitudes of Einstein gravity minimally coupled to YM theory (EYM) involving n gluons and a single graviton are reduced to pure gauge-theory amplitudes [25], 1 A EYM (1, 2, . . . , n; p) = n−1 l=1 ( p · x l ) A YM (1, 2, . . . , l, p, l + 1, . . . , n) . (1.1) The cyclically ordered gluon legs 1, 2, . . . , n refer to a single trace of nonabelian gauge-group generators t a 1 t a 2 . . . t an , and the graviton leg labelled by p is associated with a polarization vector p . The latter contracts a region variable which encompasses several lightlike gluon momenta k i , Another approach to manifest double-copy structures is the Cachazo-He-Yuan (CHY) formalism [26][27][28]. Similar to string theory, tree-level amplitudes in a variety of theories and spacetime dimensions are represented as integrals over the moduli space of punctured spheres. The CHY formulae for EYM amplitudes [29,30] were recently exploited in [31,32] to give an alternative proof of (1.1), along with different sorts of generalizations to multiple graviton insertions and additional relations for double-trace amplitudes in [31].
In absence of gluon-graviton couplings, the CHY integrands for gluon and graviton scattering are essentially identical to those in the heterotic and type-II string theory, respectively. The integrands for the gauge supermultiplet in the pure-spinor incarnation of the CHY formalism [33] were shown in [34] to reproduce the string-theory integrands [11,35,36] derived from the open pure-spinor superstring [37]. The CHY formulae for EYM amplitudes [29,30] involving both gluons and gravitons, however, do not resemble any stringtheory integrands in an obvious manner. This work contributes to clarifying the relation between the CHY integrands for EYM in [29,30] and string theory.
In this paper, we will study gluon-graviton scattering in the heterotic string, where also the gauge multiplet arises from closed strings, and the gravity sector descends from a chiral half of the bosonic string. The heterotic integrand for tree-level amplitudes of gluons and gravitons shares the overall structure of the CHY description in [29,30] but differs in JHEP11(2016)074 the contributions from the bosonic string. As we will show, the field-theory relation (1.1) of [25] with a single graviton insertion literally carries over to single-trace amplitudes of the heterotic string, including the full-fledged α -dependence on both sides. For several gravitons, however, we will see that the expressions in the heterotic theory involve additional tensor structures in their α -corrections as compared to the field-theory relations from the CHY formalism [31]. Two-graviton amplitudes, for instance, will be reduced to purely gluonic single-trace amplitudes via 2 where the contractions ( p · q)( q · p) of the graviton momenta p, q in the second line are subleading in α and therefore absent in the corresponding field-theory relation [31]. In the field-theory limit α → 0, the relations (1.3) and (1.4) of the heterotic string are easily seen to reproduce the EYM relations (1.1) due to [25] and their two-graviton counterparts given in [31] and (3.16). Analogous statements apply to the amplitude relation (3.26) involving three gravitons and n gluons. The heterotic-string approach to gluon-graviton scattering is complementary to the type-I results of [6,24,25] and governed by worldsheets of sphere topology instead of disks. Accordingly, our methods to derive heterotic amplitude relations (1.3) and (1.4) are based on integrand manipulations similar to those in the CHY formalism [31] and do not involve any monodromy considerations. While multitrace amplitudes in the type-I theory appear at higher orders of the genus expansion, the heterotic string also incorporates multitrace contributions at tree level. This will be exploited to reduce color-ordered doubletrace amplitudes of the heterotic string involving at most one graviton to their single-trace counterparts. The main results for the double-trace sector can be found in (4.2) and (4.18), and their derivation closely follows the integrand manipulations of the CHY formalism [31].
The amplitude relations of this work also connect gluon-graviton scattering in the heterotic string with pure gluon amplitudes of the type-I string, thereby hinting an extension of the heterotic-type-I duality [38] along the lines of [39]: based on recent results on the 2 Note that the terms with i = 1 in the second and third line of (1.4) yield the color ordering A het (q, 1, 2, . . . , n) with q adjacent to n. Moreover, color-ordered amplitudes with adjacent gravitons can be found in each line of (1.4) including the contributions n−1 i=1 ( p · xi) ( q · xi)[A het (1, . . . , i, p, q, i + 1, . . . , n) + A het (1, . . . , i, q, p, i + 1, . . . , n)] from the i = j terms in the first line.

JHEP11(2016)074
α -expansion of tree-level amplitudes [40,41], single-trace gluon amplitudes in the heterotic string were identified in [39] to be single-valued images of their type-I counterparts. The single-valued map truncates the multiple zeta values (MZVs) in their α -expansion to the subclass which can be obtained from single-valued polylogarithms [42,43]. Hence, any instance of A het (. . .) on the right hand side of (1.3), (1.4) and later such relations may be read as sv A I (. . .), where sv denotes the single-valued map [42,43] and A I refers to color-ordered amplitudes (2.3) of the type-I superstring [35,36]. The availability of such connections with single-valued open-string amplitudes in presence of gravitons and multiple traces was pointed out in [39], and we provide an explicit realization for a variety of cases.
From the pure-spinor methods in the supersymmetric side of the heterotic-string calculation [35,36], it is clear that our results including (1.3) and (1.4) apply to the entire superamplitudes involving multiplets of ten-dimensional SYM and half-maximal supergravity. Moreover, they are valid for any dimensional reduction of the ten-dimensional theories including any helicity sector upon specialization to four dimensions. Finally, with appropriate constraints on the external states, the results hold in all string compactifications which admit a CFT description, independently on the number of supercharges they preserve.

Review
In this section, we review selected aspects of tree-level amplitudes in the type-I superstring and the heterotic string. As a starting point, we recall that the vertex operators of the gauge-and gravity multiplet in the heterotic string are given by 3

Type-I superstring
The n-point correlation functions of the supersymmetric vertex operators V SUSY (z) have been evaluated in [35,36] by means of the pure-spinor formalism 5 [37]. These correlators yield the following massless amplitudes of the type-I superstring, where S n−3 denotes the set of permutations σ of n−3 legs. All the polarization dependence in (2.3) is captured by the color-ordered tree amplitudes A YM (. . .) in ten-dimensional SYM [17,46], whereas α enters through the iterated integrals F σ (α ) whose definition will become clear from the subsequent discussion. The structure of (2.3) can be related to the KLT formula [1] by introducing holomorphic and antiholomorphic Parke-Taylor factors In this notation, the supersymmetric correlator of [35,36] takes the form [47] K with the following expressions for the Koba-Nielsen factor I n and the Mandelstam invariants s ij , The (n−3)!×(n−3)! matrix S 0 [σ|τ ] 1 in (2.5) has entries of order s n−3 ij and can also be found in the n-point supergravity tree upon replacing Cz(. . .) by another copy of A YM (. . .) [1,48], see [49][50][51] for its description in the momentum-kernel formalism. Note that, up to the Koba-Nielsen factor (2.6), the bosonic components of (2.5) yield the reduced Pfaffian of the CHY formalism [28] while the fermionic components 6 yield its supersymmetrization along the lines of [33,34]. Note that total derivatives inz have been dropped in the derivation of (2.5) [35,36], they do not contribute to the worldsheet integrals in presence of the Koba-Nielsen factor (2.6).
Color-ordered amplitudes of the type-I superstring are obtained as an iterated integral

JHEP11(2016)074
where the order of the gauge-group generators Tr(t a 1 t a 2 . . . t an ) in the accompanying trace is reflected by the integration domain D(12 . . . n) ≡ {(z 1 , z 2 , . . . , z n ) ∈ R n , −∞ < z 1 < z 2 < . . . < z n < ∞} . (2.8) The formal division by vol(SL(2, R)) instructs to remove the redundancy due to Moebius transformations z → az+b cz+d with a b c d ∈ SL(2, R). This amounts to dropping the integrations w.r.t. any three punctures z i , z j , z k and to multiply by the Jacobian z ij z ik z jk . Analogous statements hold for the inverse of vol(SL(2, C)) seen in the following section. Note that the functions F σ (α ) in (2.3) can be assembled by inserting the correlator (2.5) into (2.7).

Heterotic string
Color-dressed tree-level amplitudes in heterotic string theory involving n gluons and r gravitons are computed from the prescription [44] M het (1, 2, . . . , n; n+1, . . . , n+r) In view of the supersymmetric constituents V SUSY (z) in the vertex operators (2.1) and (2.2), it is clear that the type-I correlator K n+r ({z j }) given by (2.5) can always be factored out from the integrand. In the following, we will discuss a convenient treatment of the additional ingredients J a and µ ∂X µ of the heterotic string which enter through the second line of M het (1, 2, . . . , n; n+1, . . . , n+r) = C n+r Note that the definition (2.5) of the supersymmetric correlator K n+r ({z j }) already incorporates the Koba-Nielsen factor (2.6), that is why we will be interested in the reduced correlation function

Current correlators and color ordering
Tree-level correlation functions among Kac-Moody currents J a boil down to products of color traces Tr(t a 1 t a 2 . . . t an ) accompanied by Parke-Taylor factors (z 12 z 23 . . . z n,1 ) −1 in the corresponding cyclic order. Expressions for the full correlators with a systematic discussion of their multi-trace structure can be found in [53,54]. We will mostly be interested in the

JHEP11(2016)074
single-and double-trace contributions with an obvious generalization to arbitrary products of traces. Upon insertion into (2.10), this leads to the following expressions for the color-ordered single-trace amplitudes involve an additional minus sign relative to the single-trace sector which can be traced back to the current correlator.

Single-trace gauge amplitudes as the central building block
In absence of gravitons, (2.14) yields the following single-trace amplitudes for the gaugemultiplet which will serve as the central building block for the remainder of this work. It will be shown to capture the supersymmetric polarizations in a variety of single-and double-trace cases (2.14) and (2.15) involving both gluons and gravitons. For this purpose, the second line of their integrand needs to be expanded in (n+r)-particle Parke-Taylor factors (2.4) such as of the gravitons' half-polarizations j and momenta k j as well as α . The rearrangement in (2.17) will be explicitly performed in section 3 for single-trace amplitudes (2.14) with up to three gravitons, leading to amplitude relations A het (1, 2, . . . , n; n+1, . . . , n+r) = . . , n+r)) . (2.18)
(2.20) Note that the rearrangements in (2.17) and (2.19) usually include integration by parts, where the total derivatives in z do not contribute to the integrals of (2.10) and will be dropped.

Type-I versus heterotic string and BCJ relations
For the single-trace amplitudes of the gauge multiplet, an intriguing connection between the type-I expression (2.7) and its heterotic counterpart in (2.16) has been found in [39]. They are related by adjusting the integration measures via and the α -expansion of the heterotic-string amplitude follows from the single-valued image [42,43] of the accompanying MZVs 7 [39] A het (1, 2, . . . , n) = sv A I (1, 2, . . . , n) .

(2.22)
As mentioned in the introduction, the image of the single-valued map sv only comprises those MZVs which arise from single-valued polylogarithms at unity, see [42,43] or [39,41] for further details. The identity (2.22) can be inserted into all the amplitude relations (2.18) and (2.20) to be derived in the following, establishing an explicit connection between a variety of heterotic-string amplitudes and the type-I superstring, in lines with the general argument of [39]. As a consequence of (2.22), heterotic single-trace amplitudes (2.16) were concluded in [39] to obey the BCJ relations known from field-theory amplitudes [2], whose repeated application yields the expression (2.6) for the Koba-Nielsen factor as well as the following correlators (2.11) of interest for amplitudes with up to three gravitons: Note that the summation range j =p is understood to also include the graviton legs different from p, for instance j ∈ {1, 2, . . . , n, q, r} in the three-graviton correlator (2.27). Any contraction of the form ( p · q ) is accompanied by an inverse power of α , and the obvious r-graviton generalization will have a maximum of r/2 such factors. Note that the right hand side of (2.25) coincides with the diagonal elements of the C-matrix which constitutes the off-diagonal block of the Pfaffian in the CHY formalism [27,28]. Since these diagonal elements have played a crucial role for the derivation of EYM relations in [31,32], it is not surprising that (2.25) to (2.27) will form the backbone of very similar derivations in heterotic string theory.

Single-trace amplitudes with graviton insertions
In this section, we will reduce single-trace tree-level amplitudes (2.14) of the heterotic string involving n gluons and up to three gravitons to purely gluonic single-trace trees (2.16). We derive the amplitude relations (1.3) and (1.4) for one and two gravitons as well as the analogous identity (3.26) for three gravitons. The computations closely resemble the manipulations in the CHY formalism which have been performed in [31] to obtain the analogous amplitude relations for the EYM field-theory limits. Among the α -corrections to two-and three-graviton amplitudes, we identify a few additional tensor structures at subleading orders in α which are thereby absent in the field-theory relations of [31].

One graviton
Single-trace amplitudes (2.14) involving n gluons and a single graviton at leg n+1 → p with momentum k p → p can be simplified by rewriting the reduced correlator (2.25) in terms of region momenta x j defined in (1.2), We have used transversality ( p · p) = 0, momentum conservation k n = −p − n−1 j=1 k j and rearranged the dependence on z j , see [31,32] for analogous manipulations in the CHY formalism. The z-dependence on the right hand side has the suitable form to convert an n-particle Parke-Taylor factor (2.4) into its (n+1)-particle counterparts,  [25] and later from the CHY formalism [31,32]. Note that gauge invariance of (1.3) under p → p follows from the BCJ relations (2.23) among single-trace amplitudes in the heterotic string.

Two gravitons
As a first step towards the single-trace case with two gravitons and n gluons, we express the relevant reduced correlator in a form similar to (3.1): with the definition adjusted to n color-ordered gluons, the correlator (2.26) takes the form after eliminating ( p · k n ) and ( q · k n ) via transversality and momentum conservation. We refer to the graviton momenta via k p → p and k q → q, respectively. It remains to cast the product of (3.4) with C z (1, 2, . . . , n) into the desired form (2.17), and we will now discuss two kinds of techniques for this purpose.

Algebraic rearrangements
For all the terms in (3.4) involving at least one factor of (3.3), one can literally repeat the manipulations in section 4.1 and section 4.2 of [31]: the z-dependence of C pp and C qq can be merged with C z (1, 2, . . . , n) to yield an (n+2)-particle Parke-Taylor factor by iterating elementary manipulations such as 1 Intermediate steps are facilitated by a choice of SL(2, C)-frame with z n → ∞, and the frame independent form of the results is given by [31] C z (1, 2, . . . , n) . . , i, p, i+1, . . . , j, q, j+1, . . . , n) The summation range for the ordered set σ on the right hand side of (3.7) involves the shuffle product ¡ of α ≡ {a 1 , a 2 , . . . , a i } and β ≡ {b 1 , b 2 , . . . , b j }. It is defined to collect all permutations of the union α ∪ β that preserve the individual orderings among the a l and b l , which is equivalent to the recursion (3.8)

Integration by parts
The double pole z −2 p,q in the last two terms of (3.4) requires an integration by parts to combine with C z (1, 2, . . . , n) and to yield (n+2)-particle Parke-Taylor factors. The idea is to discard a total derivative of the Koba-Nielsen factor (2.6) which simplifies in an SL(2, C)-frame where z n → ∞, and integrates to zero in the open-and closed-string measures in (2.21). The right hand side makes contact with the scattering equations except for the offset s pq ↔ s pq − 1 2α . By virtue of (3.9), the last two contributions in (3.4) along with ( p · q ) − 2α ( p · q)( q · p) can be rewritten as

JHEP11(2016)074
The factor of z i,n zq,n at the end of the first step has been inserted to restore covariance under SL(2, C), and it drops out as z i,n zq,n → 1 in the frame z n → ∞ of interest to (3.9). Up to the offset s pq ↔ s pq − 1 2α , the manipulations in (3.10) have also been performed in section 4.3 of [31]. In view of its repeated appearance in the three-graviton counterpart, we introduce the shorthand for the ubiquitous combination arising from (3.10). In passing to the second line, we have expanded the denominator as a geometric series, and all the α -corrections are seen to involve the gauge invariant linearized field strength s ij ( i · j ) − ( i · k j )( j · k i ). By (3.10) and (3.11), the net contribution from the double-pole part of (3.4) is given by (3.12) Note that the same rearrangements could have been performed with p and q exchanged, and the hidden symmetry of the right hand side of (3.12) under p ↔ q follows from further integrations by parts or from BCJ relations (2.23) among the resulting A het (. . .). We will explicitly symmetrize the final form for the amplitude relation in p and q.

Tachyon suppression and transcendentality properties
Note that α enters (3.13) and (3.14) in two different places, firstly through the geometric series (3.11) represented by E pq and secondly through the expansion of A het (. . .) on the right hand side. By (2.22) [39], the latter descends from the type-I results of [35,36,40,47,58], where the powers of α tie in with the transcendental weights of the accompanying MZVs. This property of gluonic single-trace amplitudes of the heterotic string is known as uniform transcendentality.
The series representation of E pq in (3.11) obviously violates uniform transcendentality, and so does any single-trace amplitude with r ≥ 2 graviton insertions, starting with (3.13). Only the one-graviton case (1.3) preserves uniform transcendentality which may be attributed to the particularly simple tensor structures. This ties in with the breakdown of uniform transcendentality in the bosonic string and its implications for purely gravitational amplitudes in the heterotic string [60].
One might be worried that the pole of E pq ∼ (1 − 2α s pq ) −1 at (k p + k q ) 2 = α signals tachyon propagation. However, it is clear from level matching that such poles cannot occur, and their cancellation rests on the compensating zeros of the single-trace amplitudes as (k p + k q ) 2 → α . In the above four-point example (3.14), the pole of E 34 ∼ (1 − 2α s 34 ) −1 coincides with a zero of In the field-theory limit, the two α -dependent quantities on the right hand side of (3.13) reduce to A het → A YM and E pq → ( p · q ), respectively, and we are left with the EYM amplitude relation A EYM (1, 2, . . . , n; p, q) = n−1 1=i≤j ( p · x i ) ( q · x j ) A YM (1, . . . , i, p, i+1, . . . , j, q, j+1, . . . , n) known from [31].

Three gravitons
For three gravitons and n gluons, the reduced correlator (2.27) takes the form p,q C rr + ( r · q) z q,n z q,r z r,n + ( r · p) z p,n z p,r z r,n + cyc(p, q, r) + C rr + ( r ·q) z q,n z q,r z r,n + ( r ·p) z p,n z p,r z r,n C qq + ( q ·p) z p,n z p,q z q,n + ( q ·r) z r,n z r,q z q,n × C pp + ( p ·q) z q,n z q,p z p,n + ( p ·r) z r,n z r,p z p,n , (3.17) which serves as a convenient starting point towards amplitude relations. One can identify nine tensor structures which are inequivalent under permutations of the graviton polarizations p , q , r and momenta p, q, r: The three classes of tensor structures (vi), (viii) and (ix) were observed to drop out from the field-theory relation found in [31], and we will identify them to be subleading in α similar to the terms ( p · q)( q · p) in the two-graviton case. In the following subsections, the reduction (2.17) to (n+3)-particle Parke-Taylor factors will be performed for all of (i) to (ix).

Tensor structures (v) and (vi)
It turns out that the representatives ( p · q )C rr and ( p · q)( q · p)C rr of the two tensor structures (v) and (vi) in (3.18) multiply the same worldsheet functions (including a double pole in z p,q ) and conspire to E pq defined in (3.11). We can apply the integration by parts identity (3.9) with an additional leg j inserted into the Parke-Taylor factor C z (1, 2, . . . , n): (1, 2, . . . , j, r, j+1, . . . , n) + cyc(p, q, r) Note that the result closely resembles section 5.2 of [31] up to the α -corrections in ( p · q ) → E pq .

Tensor structures (vii) and (viii)
The representatives ( p · q )( r · q) and ( p · q)( q · p)( r · q) of the tensor structures (vii) and (viii) in (3.18) turn out to involve the same dependence on the punctures and again combine to E pq : p,q ( r · q) z q,n z q,r z r,n + perm(p, q, r)

JHEP11(2016)074
which is valid in the frame where z n → ∞ and slightly generalizes (3.9). Note that the analogous scattering equation has been exploited in section 5.1 of [31], where the factor of (z p,q z q,r z r,p ) −1 in the last line is cancelled by other parts of the CHY integrands which do not have any heterotic-string counterpart.

Tensor structures (ix)
It is convenient to combine the last tensor structures (ix) in (3.18) with the terms ∼ (z p,q z q,r z r,p ) −1 in the last line of (3.21): s pr z p,q z q,r z r,p + perm(p, q, r) In passing to the second line, the tensor structures (ix) and E pq ( r ·q)s pr have been combined to yield the totally antisymmetric and gauge invariant generalization of (3.11), (3.24) which reproduces the Lorentz trace of three linearized field strengths at the leading order in α . In passing to the third line of (3.23), we are anticipating the integration by parts identity (4.15), see the later section 4.1 for its derivation. Moreover, the following definition of multiparticle Mandelstam invariants has been used: (3.25) Note that the entire contribution in (3.23) is subleading in α and does not have any analogue in EYM and the CHY computation of [31]. Moreover, permutation invariance has been obscured in the last step of (3.23), but it can be checked via integration by parts or BCJ relations and will be reinstated by manual averaging in the final result. reduces three-graviton single-trace trees to purely gluonic ones,

JHEP11(2016)074 4 Towards multitrace amplitudes
In this section, we reduce multitrace amplitudes (2.15) involving n gluons and at most one graviton to their single-trace counterparts in the gauge sector. The main challenge is posed by the conversion (2.19) of products of Parke-Taylor factors and graviton contributions to a single Parke-Taylor factor. This rearrangement will be addressed via integration by parts which resembles the application of scattering equations in the analogous CHY analysis [31] but, as we shall see, organizes the low-energy regime in a different manner.

No graviton
In absence of gravitons, the general formula (2.15) for double-trace amplitudes reduces to . . , n) . (4.1) As we will demonstrate below, suitable rearrangements of the Parke-Taylor factors on the right hand side via integration by parts lead to the following reduction to single-trace amplitudes:     (1,2,3,4,5) .
The right hand sides of the above relations do not manifest cyclicity within the traces or symmetry under their exchange. These properties can be checked to hold by means of BCJ relations (2.23) of the single-trace amplitudes. In the field-theory limit, double-and single-trace amplitudes of the heterotic string reduce to those of EYM and YM, respectively, i.e. A het The additional factor of 2α in the conversion between A het (2) and 2α A EYM stems from the gravitational constant in the coupling between graviton and gluon which we suppress in the normalization conventions for A EYM chosen in this work and [31]. Then, the EYM relations of [31] are recovered in the low-energy limit of (4. (4.8) The following subsections are devoted to the derivation of the all-multiplicity identity (4.2).

Total derivatives
In this section, we generalize the total derivative (3.9) which was used in (3.10) to remove the Parke-Taylor factor C z (p, q) from the two-graviton integrand. This will be done in an SL(2, C)-frame where z r → ∞ and Koba-Nielsen derivatives simplify to ∂ ∂z i I n = 2α I n j =i,r s ij z ij excluding j = r. We recapitulate the two-particle case (3.9) for legs 1, 2, . . . , r and p, q as and find analogous identities relevant to (4.4) and (4.5) when C z (1, 2, . . . , r) is augmented by a second Parke-Taylor factor with legs p, q, t and p, q, t, u,

Parke-Taylor rearrangements
Once we drop the total derivatives in the above identities and cast the contributions from see (3.23) for an earlier application of (4.15). These integrations by parts can be viewed as the string-theory counterparts of the cross-ratio relations of [61], also see [62,63] for related work on CHY integrals involving multiple Parke-Taylor factors. Note that (4.13) is instrumental for the reduction of disk integrals to an (n−3)! basis in the open superstring [35,36] and the open bosonic string [60]. However, as exemplified by the following section, (4.13) is not yet sufficient to address the most general integrand with suitable SL(2, C) weights.

Comparison with CHY and breakdown of uniform transcendentality
The above integration by parts identities (4.13) which reduce products of two Parke-Taylor factors to a single one are seen to introduce factors of (1 − 2α s 12...r ) −1 into the amplitude relation (4.2). It is instructive to compare this with the EYM integrands for gluonic double-trace amplitudes in the CHY formalism [29,30]. These CHY integrands augment the two Parke-Taylor factors by an explicit numerator factor of s 12...r which can be thought of as compensating the α -dependent part of the string-theory specific denominator (1 − 2α s 12...r ) −1 : a sequence of scattering equations [31,61] only yields s −1 12...r in the place of (s 12...r − 1 2α ) −1 when converting a product of two Parke-Taylor factors to a single one. In the field-theory limit, the net effect of (1 − 2α s 12...r ) −1 → 1 is the same as the conspiration between the outcome ∼ s −1 12...r of scattering equations and the additional numerator s 12...r in the CHY integrand. That is why the heterotic string and the CHY formalism converge on the EYM amplitude relations (4.8), although the opening line for their integrands differ by s 12...r .
Similar to the discussion in section 3.2.4, the tentative tachyon poles of (1−2α s 12...r ) −1 at (k 1 + . . . +k r ) 2 = α are cancelled by the zeros of the accompanying worldsheet integrals in the single-trace amplitudes, in lines with level matching. A geometric-series expansion of the denominator factors shows explicitly that double-trace amplitudes of gluons do not exhibit uniform transcendentality. Since the complexity of the integrations by parts and the associated denominators grows with the number of traces and graviton insertions (see the next section) we expect that uniform transcendentality is violated in any multitrace amplitude of the heterotic string.
It will be elaborated in the following paragraphs that suitable integrations by parts combine the Parke-Taylor factors in the second line with the graviton contribution R n ( p ; {k j , z j }) given in (2.25)

A new class of integrations by parts
The leftover task is to express the first line of (4.22) ∼ ( p · x r ) in Parke-Taylor form, where contributions with one or two factors of C z (. . .) are equally acceptable in view of (4.2). It turns out that a sequence of integrations by parts along the lines of section 4.1.1 yields (1, 2, . . . , l, p, l+1, . . . , r) C z (r+1, . . . , n) (4.23) when integrated against the (n+1)-particle Koba-Nielsen factor in (4.17). The antisymmetry of the left hand side under exchange of (r, 1, 2, . . . , r−1) ↔ (r+1, . . . , n) is not reflected in the presentation of the right hand side. Accordingly, the contributions of (4.23) to the amplitude relation (4.18) -specifically its third and fourth line -do not manifest the symmetry under exchange of the gluon traces. Still, the hidden exchange symmetry can be verified by expanding the double-trace amplitudes in a single-trace basis via (4.2) and applying BCJ relations (2.23).

Comparison with CHY
The CHY integrand for multitrace EYM amplitudes with any number of graviton insertions has been given in [30]. The string-theory correlator R n ( p ; {k j , z j }) for a single graviton was observed earlier on to match with the diagonal elements C pp of the C-matrix in the reduced Pfaffian [27]. Indeed, C pp appears in the double-trace CHY integrand for one graviton [30], but it is accompanied by a factor of s 12...p as well as extra summands spelt out in [31]. It turns out that the additional ingredients of the CHY integrand besides R n ( p ; {k j , z j }) compensate for the different organizations of the integral over C z (1, 2, . . . , r) z r,r+1 zr,pz p,r+1 C z (r+1, . . . , n) in the two setups. In string theory, the products of Parke-Taylor factors in the middle line of the integration by parts identity (4.23) are subleading in α and do not contribute to the EYM relation (4.20). The counterpart of this middle line in the analogous scattering equations, however, does contribute to the field-theory order and cancels some of the extra terms of the CHY integrand which do not appear in the string-theory correlator. This example gives rise to suspect that a systematic handle on string integrals with multiple Parke-Taylor factors and additional cross ratios as in (4.23) might allow for a streamlined approach to EYM relations involving multiple traces, without the need to track spurious sectors of the CHY integrands.

Conclusions
We have presented novel relations between tree-level amplitudes in heterotic string theory involving gluons and gravitons, including their supersymmetry multiplets and αcorrections. Purely gluonic single-trace amplitudes were identified as the fundamental building blocks to capture both single-and double-trace amplitudes with additional graviton insertions and an arbitrary number of gluons in each trace. The derivation of such amplitude relations is based on algebraic rearrangements of the worldsheet integrands and integration by parts. These techniques have led us to explicit results for single-trace amplitudes with up three gravitons in (1.3), (3.13) and (3.26) as well as double-trace amplitudes with at most one graviton, see (4.2) and (4.18).
The algorithmic derivation of the results in this work gives rise to expect that any multitrace amplitude with an arbitrary number of graviton insertions can be reduced to single-trace amplitudes among gluons by iterating the above techniques. 8 Instead of performing the tedious but by now straightforward analysis for additional multigraviton and multitrace cases and displaying the resulting formulae of increasing length, we hope to address the general problem through a recursion in the number of gravitons and traces in the near future.
Our results extend the connections between tree-level amplitudes in different string vacua found in [39]. In this reference, gluonic single-trace amplitudes in heterotic string theory were related to single-valued [42,43] gluon trees of the type-I as well as graviton trees of the type-II superstring. On these grounds, our findings pave the way to reducing the entire tree-level S-matrix for gluons and gravitons in heterotic string theory to purely gluonic type-I input, augmented by the additional polarization vectors from the gravitons. This goes beyond the heterotic-type-I duality [38] and might carry some general lessons on the web of string dualities where each theory is interpreted as a separate vacuum of an overarching M-theory [64].
Finally, it would be interesting to directly compare the α -expansions of single-trace amplitudes involving the same gluon and graviton states between the heterotic string and the type-I superstring. However, tentative matchings between the heterotic-string results of this work and type-I expressions in [6,24,25] can only occur in specific subsectors of their α -expansion: the single-valued image of the type-I amplitudes and the uniformtranscendentality parts 9 on the heterotic side. Uniform transcendentality refers to the correspondence between powers of α and the weights of the accompanying MZVs, and we observe it to be violated by multitrace amplitudes and by single-trace amplitudes with two or more graviton states in heterotic string theory.
In summary, this work contributes new structural insights into the tree-level S-matrix of string theories and is expected to open the door to numerous further research directions.