Gravitational Scattering in the High-Energy Limit

Any gravitational scattering amplitude takes a remarkably simple factorized form at tree level in multi-Regge kinematics (MRK), where the produced particles are strongly ordered in rapidity. Very recently, it was shown that also the scattering equations have a very simple structure in MRK. In this paper we study Einstein gravity amplitudes in MRK in the framework of the scattering equations. We present a new derivation of the multi-Regge factorization of tree-level amplitudes with any number of external gravitons and any helicity configuration.


Introduction
Inspired by Witten's twistor string theory [1], enormous advances have been made in understanding the mathematical structures behind scattering amplitudes. Among them, the tree-level S-matrix in many theories is expressed as a multiple contour integral in a large class of new formulations, e.g. [2][3][4][5][6] (see also [7] for a recent review). Along this line, recent progress was achieved by Cachazo, He and Yuan (CHY) who proposed a new framework where any tree-level scattering amplitudes in massless field theories in arbitrary spacetime dimensions is reformulated as a multiple integral over the moduli space of Riemann spheres with n marked points [8][9][10][11]. The integral is completely localized on the zeroes of the scattering equations that are independent of the theory. Various worldsheet models were constructed based on ambitwistor strings [12][13][14][15] that provide in particular an approach to extend the CHY formalism to loop level [16][17][18][19][20].
As the backbone of the CHY formalism, the scattering equations are universal for all theories and play a fundamental role in this content. It is thus important to reveal mathematical structures behind the scattering equations. In the single-soft limit, expanding the scattering equations as well as other ingredients of CHY formulas around the soft momentum provides a beautiful framework to produce various soft theorems in many theories, for example the soft graviton theorem up to sub-sub-leading order [10,[21][22][23][24]. More

JHEP02(2019)112
interestingly, based on the special behavior of the solutions of the scattering equations in the double-soft limit, many new universal double-soft theorems were obtained [25][26][27]. Similarly, the collinear limits of the amplitudes were also investigated up to sub-leading order in Yang-Mills, gravity and cubic scalar theories [28]. It was found that the solutions of the scattering equations can be interpreted as the zeros of the Jacobi polynomials in a two-parameter family of kinematics [29]. In the so-called "positive region" of the space of kinematic invariants, the scattering equations can be interpreted as the equilibrium equations for a stable system of n−3 particles on a finite real interval [30].
Very recently, an investigation of the high-energy limit of the scattering equations has been initiated [31]. It was observed that in multi-Regge kinematics (MRK) where the final state particles are strongly ordered in rapidity, each solution of the scattering equations displays the same hierarchy as the rapidity ordering. More remarkably, in four dimensions in the spinor-helicity formalism, the scattering equations can be exactly solved at the leading order in the multi-Regge limit for each "helicity configuration" for any number of external particles. As a result, multi-Regge factorization of gluon amplitudes was exactly derived from the scattering equations.
It is natural to expect a similar simplification to appear in gravity in the multi-Regge limit. Indeed, a compact formula for n graviton scattering at tree level in MRK was obtained from t-channel unitarity methods by Lipatov more than three decades ago [32,33]. This paper aims at extending the analysis in [31] from gauge theory amplitudes to Einstein gravity theory with a goal to provide a new alternative understanding of the high-energy limit of gravitational scattering. We first translate the Lipatov formula in the spinor-helicity language, and then derive it using the four-dimensional scattering equations. While the framework presented in [31] is applicable to graviton amplitudes because of the universality of the scattering equations, this is highly non-trivial since graviton amplitudes have a rather complicated structure even in the MHV sector. This paper is organized as follows. In section 2 we briefly review the multi-Regge factorization of tree-level graviton amplitudes and the scattering equation formalism, which are the two most important ingredients of this paper. In section 3 we derive the factorized form of graviton amplitudes in the multi-Regge regime. We start by studying amplitudes in the MHV sector in section 3.1, and then extend the analysis to all helicity sectors in section 3.2. We include two appendices with technical proofs omitted throughout the main text.

Preliminaries
The goal of this paper is to perform an investigation of the multi-Regge behavior of treelevel scattering amplitudes in Einstein gravity in the framework of the scattering equations. Therefore, we provide a brief review of multi-Regge kinematics and the scattering equations in this section before presenting the main result in subsequent sections. We follow the notations of ref. [31] in this paper.

Multi-Regge kinematics
For a 2 → (n−2) scattering, multi-Regge kinematics is defined as the regime where the final-state particles are strongly ordered in rapidity while having comparable transverse momenta, i.e., where k ⊥ a denote the transverse momenta, and in four dimensions we define the complexified transverse momenta as k ⊥ a = k x a + ik y a . Employing lightcone coordinates k a = (k + a , k − a , k ⊥ a ) with k ± a = k 0 a ± k z a , the strong ordering in rapidity is equivalent to a strong ordering in k + -components as follows: It is convenient to work in the center-of-momentum frame where two incoming particles are back-to-back on the z-axis, where s is the square of the center-of-mass energy, and we take a convention of considering all momenta outgoing. In this kinematical regime, the tree-level scattering amplitude in gravity takes a surprisingly simple factorized form: any n graviton amplitude is given by only one Feynman graph with two kinds of effective vertices (see figure 1) [32,33]. More precisely, one has where we define q a = a−1 i=2 k i with 4 ≤ a ≤ n. Here some overall factor including the Gravitational constant κ 2 = 8πG N has been stripped off. The effective graviton-graviton- which is manifestly a double copy of the gluon-gluon-(Reggeized gluon) vertex defined as Similarly, the effective Reggeon-Reggeon-graviton vertex can also be obtained as the double copy of gauge theory vertices: where C µ is the famous Lipatov vertex of (Reggeized gluon)-(Reggeized gluon)-gluon in QCD [34] 1 In this paper, Reggeon denotes 'Reggeized graviton'. with (q ⊥ ) µ ≡ (0, 0; q ⊥ ), while N µ is the so-called QED Bremsstrahlung vertex: The contractions between vertices and the polarization tensors of external gravitons give the gravitational impact factor and gravitational Lipatov vertex appearing in formula (2.4), i.e., (2.11) At this point, we would like to make some comments on the effective vertices (2.5) and (2.8) in gravity in MRK. These effective vertices have also been derived from an effective action (cf. e.g. [35][36][37]). It is extremely remarkable that the double copy relation between gravity and gauge theories was uncovered for the first time in MRK. In general kinematics, Kawai, Lewellen and Tye (KLT) found that a closed string amplitude can be expressed in terms of sums of products of two open string amplitudes [38]. In the field theory limit, the KLT relation naturally implies the double copy relation between amplitudes in gravity and Yang-Mills.
In the following we would like to translate the Lipatov formula (2.4) to modern language, say spinor-helicity variables. In four dimensions, any massless momentum can be written as a product of two spinors with opposite chirality, i.e. k αα i = λ α iλα i . To be explicit,

JHEP02(2019)112
in this paper two-component spinors are defined as (2.12) Similarly, we can also write the polarization vectors in terms of spinor-helicity variables as follows: where it is free to choose the reference spinor λ x orλ x . A polarization tensor can be expressed in terms of the symmetric-traceless tensor product of two polarization vectors. Therefore we can write the graviton polarization tensor in spinor variables as 2 where x and y stand for arbitrary reference spinors.
Using the spinor-helicity variables defined above, it is easy to compute the gravitational impact factors and Lipatov vertices defined in (2.11). A straightforward calculation gives We see from (2.16) that helicity is conserved by the impact factors, like in gauge theory.

Scattering equation formalism
The aim of this paper is to show that the multi-Regge limit of graviton amplitudes summarized in the previous section can be elegantly derived from the scattering equations. Let us continue to review the scattering equation formalism in the following. The cornerstone is of course the scattering equations [8,9,[39][40][41][42][43][44]  In four dimensions, the use of spinor-helicity variables usually leads to significant simplifications. Indeed, it has been shown that in four dimensions the scattering equations can be decomposed into different 'helicity sectors'. One of constructions is [13] where N is a subset of {1, . . . , n} with length k ∈ {2, . . . , n−2} and P is the corresponding complement. We refer to equations (2.19) as the four-dimensional scattering equations of sector k; an important property is that only the equations in the sector k are needed for N k−2 MHV amplitudes. Similarly, the system of the four-dimensional scattering equations owns a GL(2, C) = SL(2, C) × GL(1, C) symmetry.
As shown in the previous section, the multi-Regge limit is most naturally defined in terms of lightcone variables. It is thus natural to write the scattering equations in terms of lightcone coordinates. While one can get equations with lightcone coordinates by simply substituting (2.12) into the spinor-valued scattering equations (2.19), we show that one can obtain a nicer form by fixing the redundancy and rescaling variables and equations. First, let us use the GL(2, C) redundancy to fix four variables as follows: Here we always use the convention where {1, 2} ⊆ N, and we define N ≡ N\{1, 2}. Moreover, we follow the convention that elements of P and N are denoted by small and capital letters respectively, e.g. i ∈ P and I ∈ N. Second, we perform a rescaling for the t a variables as follows: Then let us also perform a rescaling for the scattering equations according to: As a consequence, we obtain a set of equations that contain only the terms linear in k + a .

JHEP02(2019)112
Explicitly, we have (2.23) We would like to emphasize that no limit has been applied to these equations, and they are completely equivalent to the four-dimensional scattering equations in (2.19), up to fixing the GL(2, C) redundancy according to (2.20) and performing the rescaling in (2.21).
In terms of lightcone variables, the formula for tree-level N k−2 MHV graviton amplitudes reads [13] If we assume that gravitons 1 and 2 carry negative helicity and use N and P to collect the babels of negative and positive helicity gravitons, the integrand function in formula (2.24) takes the following simple form: where det ′ denotes the minor with any one column and row deleted, and H is the symmetric k × k matrix defined as 27) and the H is the symmetric (n−k) × (n−k) matrix: In formula (2.24), we have eliminated four of the scattering equations and identify them with the momentum conservation delta-functions, i.e.,

Gravitational scattering in MRK
The goal of this section is to present an alternative derivation of the multi-Regge factorization of graviton amplitudes based on the scattering equations. It was conjectured in [31] that in MRK all solutions of the scattering equations satisfy Expanding the formula for graviton amplitudes (2.24) to leading power in MRK according to (3.1), it gets significantly simplified 2) where the new variables are defined as Here we fix {3, n} ⊆ P as a convention. The leading order approximation of the scattering equations in the multi-Regge limit is given by They have the following unique solution: In the following we show that we can obtain the multi-Regge factorization of graviton amplitudes by localizing integrals in (3.2) to the unique solution (3.5).
In [31], it is shown that for any function F(τ a , ζ a ) of τ a and ζ a we have .
(3.6) For gluon amplitudes, F is just a constant factor without the dependence of ζ a and τ a [31], i.e.,

JHEP02(2019)112
This paper focuses on the gravitational scattering whose F function takes The main task of the rest part of the paper is to calculate this quantity on the support of the unique solution (3.5) of the four-dimensional scattering equations in MRK.

MHV sector
Let us first consider the F function defined in (3.8) in the MHV sector. The experience from the MHV sector will be useful for evaluating the determinants in other N k MHV sectors. In this case, P = {3, . . . , n}, the solution (3.5) becomes ζ i = τ i = −1, and we have where we define (3.12) Let us choose to delete the first column and row corresponding to the particle label '3' from the matrix H. Then the reduced determinant can be written as (3.14) This matrix is nothing but the leading order approximation of Hodges' matrix 3 [45,46] in the multi-Regge limit. We can observe a lot of nice properties. In particular, a conspicuous feature is that the entriesφ ij are equal when j < i for each i-th row. As we now show in the following, this implies further simplification. By performing some elementary row/column transformations of matrix, we have This is almost an upper triangular matrix. We find it remarkable that we can nicely compute its determinant by employing the so-called the matrix determinant lemma (cf. [47,48], see also appendix A of this paper). We first simply decompose the matrix into an upper triangular part and a matrix that has only non-zero element x n in the lower left corner.
To be more precise, we writē where ϕ is nothing but the matrixφ ′ with replacing the first element of the last row x n by zero. Then by making use of the matrix determinant lemma, we have In general, it seems difficult to exactly find the inverse of the matrix ϕ. Fortunately, it is not hard to obtain the entry (ϕ −1 ) 1,n−3 by induction (see appendix B for the details of the derivation), Plugging it into (3.18) immediately gives We derive the correct multi-Regge factorization of any MHV amplitude. In order to extend the analysis in the MHV sector to other helicity sectors, at this point let us summarize some key technical ideas that have been used above. First, we can transform the matrix into a near upper triangular form by some elementary row and column operations since the matrix has a special structure in MRK. Then the matrix determinant lemma can be used to compute its determinant. We show in the next section that this technique is useful for any other helicity configuration.

All helicity configurations
Let us first consider the positive helicity part P. In MRK, the entries of (n−k) × (n−k) matrix H take For diagonal elements, we have where one has used the momentum conservation (2.30) in the third line, and u i is defined as Using the scattering equations (3.4) and their solution (3.5), one can obtain (3.27) Let us choose to delete the first column and row corresponding to particle label '3', then the reduced determinant becomes nicely where Here labels satisfy 3 < i 1 < i 2 < · · · < n. In the case of the MHV sector, since c ij = x j (i < j), this matrix is identical to the matrixφ in (3.14). More remarkably, they have a similar structure and share many properties. As a consequence, we can use the same

JHEP02(2019)112
technique to calculate the determinant as in the MHV sector. Here we summarize the result as follows (see appendix B for a detailed derivation): Let us now discuss the k × k matrix H. It is reasonable to expect that the similar structure appears in this matrix such that we can compute its determinant using the matrix determinant lemma. Let us first compute the entries of the H in MRK By using the scattering equations (3.4) and their unique solution (3.5), it is easy to obtain Then we have where one choose to remove the first column and row corresponding to particle label '1', and H ′ is defined as where I 1 < · · · < I m ∈ N, m = k−2. This matrix again displays a similar structure asφ in (3.14). Hence we can calculate its determinant by the similar technique based on the matrix determinant lemma. The final result is

JHEP02(2019)112
The details of the derivation can be found in appendix B.
Putting everything together, we obtain that any N k MHV amplitude of graviton factorizes as in agreement with Lipatov formula (2.4). Let us conclude this section by making some comments. First, we have assumed in previous sections that the two incoming particles 1 and 2 carry the same helicity. Here we show that all conclusions hold for the case where the gravitons 1 and 2 have opposite helicities. For example, let us consider amplitude M n with helicity configuration (1 + , 2 − , . . . , n − ). In this case, in MRK, we have (3.45) The factor (1 n) −8 combines with the impact factor C(1 + ; n − ) to give This shows that the amplitude

JHEP02(2019)112
Similarly, a combination of the factor (1 n) −1 and the impact factor C(1 + ; n − ) gives (3.50) which exactly agrees with the result in [32,33]. It is also easy to obtain the result for other cases where particles with other spins in the supergravity multiplet are involved in the same way.

Conclusions
We have initiated the investigation of the gravitational scattering in the multi-Regge regime in the framework of the scattering equations. Unlike gauge theory, the evaluation of the determinants of the two matrices is involved in the formula of graviton amplitudes. In general, it seems impossible to obtain the exact compact results of these determinants even in the MHV sector. We have shown that the two matrices get greatly simplified in the multi-Regge limit, and finally we can obtain compact expressions for any N k MHV sector for any number of external particles. As a consequence, we provide an elegant derivation of the tree-level multi-Regge factorization of gravitational scattering amplitudes. It should be emphasized that our analysis in this paper is based on the asymptotic behaviour of the solutions to the scattering equations in the multi-Regge limit. In [31], it is conjectured that all solutions of the scattering equations admit the same hierarchy as the rapidity ordering in MRK, as shown in (3.1) in this paper. While we do currently not have a rigorous mathematical proof, the result in this paper provides very strong support to the validity of the conjecture. As a next step, it would also be interesting to investigate the multi-Regge limit of amplitudes for more theories along this path.
While this paper has been concentrated on the multi-Regge limit, it would be interesting to study graviton amplitudes in various generalizations of MRK, where two or more produced particles have comparable rapidities. We leave this study for future work.
Acknowledgments I would especially like to thank Xiaoran Zhao for stimulating discussions, and Claude Duhr for useful discussions and a careful reading of the manuscript, as well as collaboration on related topics. I am also grateful to Song He and Brenda Penante for discussions. I would like to acknowledge the hospitality of CERN Theory Division in Geneva, Galileo Galilei Institute in Firenze, ITP, CAS in Beijing, and Tianjin University in Tianjin, where part of the work was done. This work was supported by the "Fonds Spécial de Recherche" (FSR) of the UCLouvain.

A Matrix determinant lemma
Let A be an invertible matrix, u and v be two column vectors. Then the matrix determinant lemma states that [47] (cf. also [48])

JHEP02(2019)112
Proof. Let us first see the special case of the identity matrix, i.e. A = 1. Note This ends the proof for the case of A = 1. Then from we can prove the lemme, i.e.,

B Proof of three identities
This appendix provides the details of deriving three identities used in section 3, i.e. eqs. (3.19), (3.30) and (3.41).

B.1 Equation (3.19)
In this section, we focus on the following triangular matrix: Our goal is to find (ϕ −1 ) 1,n−3 . Let us denote the last column of the inverse of this matrix as α = (α 4 , α 5 , . . . , α n ) T , and then it satisfies the following equation In the following we show that its solution is First, we can easily get from x n α n = 1 Then we assume all α j with j > i are given by (B.3), and let us solve α i from the following equation:

JHEP02(2019)112
Plugging the values of α j (j > i) given in (B.3) into this equation gives which agrees with eq. (B.3). Finally, solving the last equation gives eq. (3.19), i.e., In the following we consider the reduced determinant of the H k×k whose indices take values from the set P. First fo all, it is useful to introduce a new notation related to particle labels: I ℓ i ∈ N denotes the smallest number that satisfies I ℓ i > i ∈ P. For example, ℓ 3 = 1 because 3 < I 1 ∈ N. Then, by abuse of multiple subscripts, we can rewrite ζ i and τ i in terms of ζ I and τ I as follows: , if i > I m , (B.9) Let us also rewrite the matrix (3.29): where particle labels in P have been reordered as 3 < i 1 < · · · < i p < n with p = (n−k)−2. By performing some elementary row and column transformations, we have Using the matrix determinant lemma, we have

JHEP02(2019)112
whereΦ is nothing but H ′′ with replacing the first element of the last row x n by zero. Now our task becomes to calculate the entry in the upper right corner of the inverse of the matrixΦ. Let us denote the last column of the inverse ofΦ as: Clearly, α i 1 = (Φ −1 ) 1,p+1 . Then it can be determined by the following linear equations: 14) The solution is Using the definition of c ij and α j given by (B.15), we have j∈P,j>i where we denote Next, we calculate these tow terms for two cases respectively: the label i is bigger than the label of any negative-helicity particle or not.

JHEP02(2019)112
• Let us first consider the case of the label i is less than the largest label carried by negative-helicity particles, i.e. i < I m . In this case, we have • In the other case, i.e. i > I m , it is easy to obtain In both cases, as expected, we obtain the same results for T i and Z i . By inserting them into (B.20), we find j∈P,j>i Open Access. This article is distributed under the terms of the Creative Commons Attribution License (CC-BY 4.0), which permits any use, distribution and reproduction in any medium, provided the original author(s) and source are credited.