Two-loop doubly massive four-point amplitude involving a half-BPS and Konishi operator

The two-loop four-point amplitude of two massless SU(N) colored scalars and two color singlet operators with different virtuality described by a half-BPS and Konishi operators is calculated analytically in maximally supersymmetric Yang-Mills theory. We verify the ultraviolet behaviour of the unprotected composite operator and exponentiation of the infrared divergences with correct universal values of the anomalous dimensions in the modified dimensional reduction scheme. The amplitude is found to contain lower transcendental weight terms in addition to the highest ones and the latter has no similarity with similar amplitudes in QCD.


Introduction
Scattering amplitudes and correlation functions are the most fundamental objects in any quantum field theory (QFT). Gauge theory is the language of nature which is so far well tested through the standard model of particle physics. Even after many decades of the formal formulation of the Yang-Mills gauge theory [1], it remains a formidable task to go beyond a certain order in perturbation theory. In particular, very little is known in the nonperturbative regime. However, a special class of theories which admits a dual description in the strongly coupled sector are largely explored. The N = 4 supersymmetric Yang-Mills (SYM) belongs to this category. Through the AdS/CFT correspondence [2], certain quantities in strongly coupled sectors of N = 4 SYM is related to the weakly coupled sectors of gravity in Anti-de Sitter space. Beyond the AdS/CFT correspondence, in particular, in the weakly coupled sector, the N = 4 SYM is seen to be a time tested sandbox to explore new ideas and computational techniques. It offers perhaps the best chance to solve an interacting four dimensional QFT. Due to the underlying superconformal symmetries, often it renders many computations much simpler than their non-supersymmetric counterparts. In many cases the resulting conceptual understanding and computational developments eventually help to solve problems in a generic QFT.
Besides on-shell amplitudes, the study of form factors (FFs) has generated a surge of interest in the community. The FFs are a set of quantities which are constructed out of scattering amplitudes of on-shell states consisting of elementary particles of the theory and JHEP05(2019)066 off-shell states described by composite operators. These are calculated by evaluating the quantity of the form p σ 1 1 , . . . , p σn n |O|0 which represents the transition matrix element from vacuum |0 to an on-shell state |p σ 1 1 , . . . , p σn n through an interaction caused by the gauge invariant operator O. p n and σ n represent the corresponding momentum and quantum numbers of the n-th particle. Studying these quantities is of paramount importance. For example, any n-point planar amplitude is factorised [3] into an infrared divergent part described by a product of FF and a finite part, often called hard function. The divergent part exponentiates and the exponential term is described by some universal quantities like the light-like cusp anomalous dimension.
In recent times, several calculations [4][5][6][7][8][9][10][11][12] have been carried out on FF in the context of N = 4 SYM. A very first calculation was done long back in [13] where the two-loop contribution to Sudakov FF of a half-BPS operator belonging to the stress-energy supermultiplet was performed. This was later extended to three loops in [8] where a remarkable connection to the highest transcendental terms of the corresponding FF in quantum chromodynamics (QCD) was found. There is another operator, called Konishi [14], which is the primary operator of the Konishi supermultiplet and belongs to the non-BPS category, has drawn a lot of attention due to many interesting properties it exhibits. It is the simplest gauge invariant operator in N = 4 SYM that is not protected by supersymmetry and consequently, receives non-zero anomalous dimensions to all orders in perturbation theory. The results of the anomalous dimensions up to five loops are known in the literature [15][16][17][18][19][20][21][22][23][24][25]. The two-point FF to two loops and three-point to one-loop were computed in [26] where the former was later extended by us in [11] to three loops and the latter one to two loops by one of us in [27]. In this article, for the first time, we focus on a four-point amplitude of two different composite operators: half-BPS and Konishi. More specifically, we consider the four-point amplitude of two massless scalars (colored) in N = 4 SYM and two color singlet states with different masses represented through half-BPS and Konishi operators, see figure 1. In QCD, two loops four-point amplitudes involving two massive vector bosons V 1 V 2 were computed in [28,29]. The case for two identical operators (half-BPS/Konishi) is being investigated by us in [30]. For similar calculations in QCD involving di-Higgs boson, see [31,32].
The multiloop corrections to the processes involving massive particles are known to be very complicated. Nevertheless, the present multiloop techniques seem to be able to stand the challenges of evaluation of the loop amplitudes. For the first time in N = 4 SYM, by employing state-of-the-art techniques we compute the four-point amplitude involving two massive particles with different virtuality to two loops level analytically. Unlike the most popular and relatively modern method of unitarity, we compute the amplitude by applying the Feynman diagrammatic approach. This approach is particularly useful in the context of regularisation prescription which is required in order to regulate the infrared (IR) divergences present in the theory due to the presence of massless particles. In addition, though the N = 4 SYM is ultraviolet (UV) finite in 4-dimensions, there can be UV divergences in the FF beyond leading order because of the composite operators. In [26], it was shown that the FF of unprotected operators like Konishi calculated in four dimensional helicity (FDH) scheme [33,34] fail to produce the correct anomalous dimensions, instead in mod- ified dimensional reduction (DR) [35,36] it indeed gives the correct results. Due to the similarity between the latter scheme with dimensional regularisation [37] which is mostly used for the radiative corrections following Feynman diagrammatic approach, it is much more convenient to employ the DR scheme.
The degree of transcendentality, τ , of a function f is defined as the number of iterated integrals required to define the function f , e.g. τ (log) = 1 , τ (Li n ) = n , τ (ζ n ) = n and also we define τ (f 1 f 2 ) = τ (f 1 ) + τ (f 2 ). Algebraic factors are assigned degree zero. It is an observed [5,8,9,[38][39][40][41][42][43][44], albeit unproven fact that the results of scattering amplitudes in N = 4 SYM exhibit uniform transcendentality (UT) i.e. those can be expressed in terms of polylogarithmic functions of uniform degree 2L, where L denotes the loop order, with constant coefficients. For the planar amplitudes, this is even true for individual integrals when these are expressed in an appropriate basis of dual conformal integrals [45,46]. However, for non-planar integrals the dual conformal symmetry does not hold true. For four points, the non-planar double ladder integral does not exhibit UT, however, if it is defined with an appropriate loop dependent numerator, it does obey the UT property [46,47]. So, by making an appropriate choice of basis, one can understand the UT property of fourpoint amplitudes in N = 4 SYM [48] and N = 8 supergravity amplitudes [49,50]. Inspired by these observations, a long-standing question is floating around: is the UT property a generic feature of N = 4 SYM amplitudes? In [11,26], it has been shown that UT property breaks down for the Sudakov FF of non-protected operators like Konishi. In this article, we address this question in the context of four-point amplitude involving two different color singlet states described by a half-BPS and Konishi. We see that the UT property does not hold true.
The paper is organised as follows. In section 2, we introduce the Lagrangian, define the half-BPS and Konishi operators and describe the kinematics. In section 3, the FF is defined in terms of the renormalised matrix elements. The regularisation prescription is discussed in section 4. The UV operator renormalisation and infrared factorisation in terms of universal quantities are described in section 4.1 and 4.2, respectively. The methodology of performing two-loop computation is described in section 5. We present our results in section 6 and discuss about the symmetries which the amplitudes exhibit. We introduce the finite remainder function in section 7. In the next section 8, we elaborate on the violation of the uniform transcendental property. Finally we make concluding remarks in section 9.

JHEP05(2019)066
The 2-loop results are presented in the supplementary material, Finite ppbk.m, with the arXiv as well as JHEP submission.

Theoretical framework
The Lagrangian describing the dynamics of N = 4 SYM reads [51][52][53][54] and We consider two composite operators, a half-BPS operator [13,55], O hBPS , belonging to the stress-energy supermultiplet containing the conserved currents of N = 4 SYM and a non-BPS operator, namely, the Konishi, O nBPS , which are given by More specifically, we are interested in the four-point amplitude of two different off-shell states (J) described by these two operators which are produced from two on-shell scalar particles i.e.
where p i and q i are the corresponding 4-momentum with p 2 i = 0 and q 2 i = m 2 i and the corresponding Madelstam variables are defined as The physical region of the phase space is bounded by tu = q 2 1 q 2 2 such that it satisfies where κ is the Källén function defined as The underlying Lagrangian encapsulating the interaction of these off-shell states to the fields of N = 4 SYM is given by In order to represent the four-point scattering amplitude, we define the form factors F = 1 + ∞ n=1 a n F (n) where the component at O(a n ) is connected to the matrix elements as In the above expression, |M (n) is the n-th loop transition matrix element of the production of two off-shell particles having different masses from the on-shell states composed of two colored scalar particles and a is the 't Hooft coupling [3] given by where the γ E ≈ 0.5772 is the Euler constant and N is the Casimir of SU(N) group in adjoint representation. The goal of this article is to compute the FF at one and two loops i.e. F (1) and F (2) . In the subsequent sections, we describe the methodology of computing these quantities.

Regularisation prescription
The N = 4 SYM is UV finite in 4-dimensions i.e. the β-function vanishes to all orders in perturbation theory. This ensures the UV finiteness of the on-shell amplitudes and the FFs of the protected operators, like O hBPS . However, the unprotected operators like O nBPS , do develop UV divergences arising from short distance effects. In addition to that, this theory is not free from soft and collinear divergences (IR) due to the presence of massless fields. The on-shell amplitudes, FF of protected as well as unprotected operators give rise to these IR divergences. In order to identify these divergences, we need to regularise the theory. Among many other regularisation schemes, we employ the DR scheme for our computation. This is very much similar to the dimensional regularisation by 't Hooft and Veltman [37] and preserves SUSY. In order to do that the number of fermionic and bosonic JHEP05(2019)066 degrees of freedom (DOF) should be maintained equal throughout the calculation. This is achieved by changing the number of generations of the scalar and pseudo-scalar from n g = 3 to n g, = 3 + in d = (4 − 2 ) ≡ d space-time dimensions. Along with the (2 − 2 ) DOF of the gauge fields, the total number of bosonic DOF becomes 8 which is equal to the fermionic ones, same as that of 4-dimensions and consequently, preserves SUSY. In this scheme, all the traces of the α, β and Dirac matrices are performed in d dimensions.
In addition to the usual Lie algebra obeyed by the six anti-symmetric matrices α and β, they fulfil in d dimensions. Unlike other available schemes, like FDH scheme, the DR scheme is universally applicable to the amplitudes for the protected as well as unprotected operators.
To maintain the traceless property of the O hBPS , defined in eq. (2.4), in d dimensions, it has to be modified as

UV divergences and operator renormalisation
As mentioned in the previous section 4, the amplitudes of the unprotected composite operators like O nBPS are not UV finite which can be seen by the presence of non-zero anomalous dimensions, γ K , [15][16][17][18][19][20][21][22][23][24][25] in the FF. So, any amplitude involving unprotected operator needs to go through UV renormalisation which is performed by multiplying an overall operator renormalisation constant, This is totally different than the coupling constant or wave-function renormalisation and is solely arising due to the nature of the composite operator. The scale µ appeared in the above renormalisation group (RG) equation is introduced througĥ where,â is the coupling constant appearing in the regularised Lagrangian in d dimensions and µ 0 is the scale introduced to makeâ dimensionless. Because of the vanishing βfunctions in N = 4 SYM, the coupling constant a satisfies very simple RG equation Employing this, the RG equation for Z K , eq. (4.3), can be solved exactly to all orders which is obtained as

JHEP05(2019)066
with γ K 1 = −6 and γ K 2 = 24. The Z K up to two loops reads Z K (a(µ 2 ), ) = 1 + a(µ 2 ) 6 + a 2 (µ 2 ) 18 2 − 12 . (4.7) The UV renormalised matrix element can be written in terms of the bare quantities as where, |M (n) is the n-th loop bare amplitude. On the other hand, we can equally express |M in powers of renormalised quantities |M (n) through By comparing the eq. (4.9) with (4.8), we get the UV renormalised matrix elements in terms of bare ones as Substituting the above relations in eq. (3.1), we obtain the UV renormalised form factors at 1-and 2-loop level which are presented in this article.

Universality of IR divergences
Beyond leading order, the resulting UV renormalised FF contains IR divergences arising from the soft and collinear configurations which appear as poles in the dimensional regularisation parameter . Due to the universal nature of these singularities, these are process independent and depend only on the nature of the external particles containing SU(N) color index. In a seminal paper [56], Catani predicted these poles for n-point two loops scattering amplitudes which are related to universal anomalous dimensions. Later, a formal derivation was presented in [57] exploiting the factorisation and resummation properties of QCD amplitudes, and were subsequently generalised to all loop order in [58,59]. Following [56], we get fin . only the highest transcendental terms which turn out to be with [58]  fin . For the convenice of the readers, we write these below: fin . (4.14) The goal of this article is to compute the quantities F fin and F fin . In the next section, we describe the methodology for the computation of the four-point amplitudes.

Calculation of the amplitudes
In contrast to the most popular method of unitarity for computing the scattering amplitudes in the context of N = 4 SYM, we employ the Feynman diagrammatic approach that carries advantages in light of the regularisation scheme. The Feynman diagrams for the process under consideration are generated using QGRAF [60]. Special care is taken in order to incorporate the Majorana fermions which are its own anti-particle and consequently, destroy the flow of fermionic current in the QGRAF output. This is rectified by an in-house code based on Python. There are 2, 18 and 489 Feynman diagrams at tree, one and two loop, respectively. The QGRAF output is passed through several in-house codes written in symbolic manipulating program FORM [61] in order to apply the Feynman rules, perform Dirac, Lorentz and SU(N) color algebras. We employ the Feynman-'t Hooft gauge (ξ = 1) for the internal gluons. After evaluation of traces of Dirac's γ-matrices and contraction of Lorentz indices every Feynman diagram is expressed as linear combination of a larger number of scalar Feynman integrals which belong to the family of the massless four-point functions with two off-shell legs of different virtualities. Using the liberty of transforming the loop momenta, all the scalar integrals are categorised into 3 and 6 different integral families at 1-and 2-loop, respectively, with the help of REDUZE2 [62,63]. These scalar integrals are reduced to a smaller set of master integrals (MI) employing the integration-by-parts (IBP) [64,65] and Lorentz invariant [66] identities. Being a process involving two massive external legs of different masses, the reduction is resonably complicated. In order to achieve Upon performing the IBP reductions, there are 134 MIs at 2-loop level. Among these, 49 are related to each other by crossing and 2-integrals turn out to be same which are related to each other by simple transformation of the loop momentum. At the end, we have 84 independent master integrals. These are matched with the set of MIs already present in the literature. The integrals for the non-equal masses were first derived in [73][74][75]. A subset of these MIs were computed in [76,77]. An independent derivation of these MIs was performed in [28] where the solutions are optimised for the numerical evaluation. For our current calculation, we use the optimised solutions of the MIs presented in [28] which are available in HepForge [78] in computer readable format. To get the optimised solutions of the MIs [28], a convenient choice of variables x, y, z and m 2 is made: which rationalise the root of κ. The physical region is constrained through The symbol alphabet of the MIs involve 19 letters l i ∈ {x, 1 + x, y, 1 − y, z, 1 − z, −y + z, 1 + y − z, 1 + xy, 1 + xz, xy + z, 2xy + x 2 y + x 2 y 2 + z − x 2 yz, 1 + x + y + xy + xy 2 − z − xz − xyz, 1 + y + xy + y 2 + xy 2 − z − yz − xyz, −xy + z + xz + xyz, The results are presented as Laurent series expansion in to weight 4 which enable us to obtain the 2-loop results to O( 0 ). In the next section 6, we discuss and present the results.

Results of the one & two loops amplitudes
We calculate the 1-and 2-loop form factors, eq. (3.1), to O( 2 ) and O( 0 ) where the results contain maximum weight 4 terms. The results are expressed in terms of logarithms log, classical polylogarithms Li n (n = 2, 3, 4) and Li 2,2 where the latter one is defined in [28]. The square of the Born amplitude is obtained as The above expression is symmetric under the exchange of z ↔ (1 + y − z) which translates to t ↔ u or equivalently p 1 ↔ p 2 . This is in accordance with the expectation. In the next subsection 6.1, we elaborate on this and present the results for 1-loop and 2-loop partially. Complete results can be found from the supplementary material supplied with the arXiv as well as JHEP submission in Mathematica format.

Permutation symmetry
The process under consideration is symmetric under the exchange of initial state particles that is under p 1 ↔ p 2 which translates to the Mandelstam variables as This symmetry should be reflected in the amplitude. Since the results are expressed in terms of functions containing variables defined through eq. (5.1), it is not straightforward to check the symmetry analytically. In order to check this, we go to a new set of variables, Z and Y introduced through The permutation symmetry implies the amplitude should remain invariant under the exchange of Z ↔ Y. In terms of these new variables, the born amplitude eq. (6.1) becomes where the Z ↔ Y symmetry is explicit. Note that the leading order amplitude does depend on . In above eq. (6.4), only the O( 0 ) term is presented. To check this behaviour for 1and 2-loop amplitudes, we express the letters in eq. (5.3) in terms of these variables which turn out to be We see under the exchange of Z ↔ Y, some letters transform among themselves: l 5 ↔ l 8 , l 6 ↔ l 7 , l 10 ↔ l 13 , l 11 ↔ l 12 , l 14 ↔ l 15 , l 16 ↔ l 18 , l 17 ↔ l 19 , (6.6) whereas the remaining ones are unaffected. Using these, we see the 1-loop form factor indeed remains invariant under the exchange of p 1 ↔ p 2 . By making this behaviour explicit, we present the O( 0 ) term of the 1-loop finite part (related to the hard function in QCD), eq. (4.14), below F (1) fin = − 7 − 3iπ + 6 log l 1 + 3 log l 3 − 1 3 π 2 − 2iπ log (l 2 l 9 ) + log l 3 log(l 2 l 9 ) + 2u t + u log 2 l 10 + log 2 l 11 − log 2 l 5 − 2 log l 11 log (l 1 l 3 ) + 2 log l 5 log l 1 l 3 l 2 l 9 + 2Li 2 1 l 10 + 2Li 2 l 1 l 3 l 11 + 2iπ log l 10 l 11 l 5 + t ↔ u . (6.7)

F
(2),τ (2) fin = − 5n 1 3d 1 π 2 +iπ 28 log (l 2 l 9 )+ u t+u (56 log l 5 −32 log l 11 )−48 log l 1 − 4un 5 d 2 ×log l 10 − 6n 4 d 1 log l 3 + 24n 3 d 1 log l 6 +42 log 2 l 1 + u t+u 28 log 2 l 5 −16 log 2 l 11 − 2un 5 d 2 log 2 l 10 + 3n 2 d 1 log 2 l 3 + 6n 3 d 1 log 2 l 6 + 32u t+u log l 11 log(l 1 l 3 )+48 log l 1 log l 3 −14 log l 3 log(l 2 l 9 )− 56u t+u log l 5 log where the {n i } and {d i } are given by n 1 = 11m 4 (t + u) + 5s 2 (t + u) + 11tu(t + u) + s(5t 2 + 22tu + 5u 2 ) − m 2 (t + u)(16s + 11(t + u)) , Im[F (2) fin ] x cosθ=-0.2 cosθ=-0.4 cosθ=-0.6 Figure 2. Behaviour of the real and imaginary parts of F fin as a function of variable x for different values of cos θ. n 2 = 6m 4 (t + u) + 5s 2 (t + u) + 6tu(t + u) + s(5t 2 + 12tu + 5u 2 ) − m 2 (t + u)(11s + 6(t + u)) , − m 2 (t + u)(12s + 7(t + u)) , (6.10) The remaining, namely, the highest transcendental term F (2),τ (4) fin can be found in the supplementary material supplied with the arXiv as well as JHEP submission. The renormalisation scale µ is set equal to m throughout the article which can be restored by using renormalisation group evolution. The 1-loop form factor is calculated to O( 2 ) and 2-loop to O( 0 ) or to weight 4 terms. The t ↔ u symmetry is checked for all these terms. It is checked analytically for the 1-loop terms to O( ) and for the lower transcendental terms of the 2-loop to O( 0 ). On the other hand, for the O( 2 ) term at 1-loop and the highest transcendental part of the 2-loop form factor, numerical checks are performed which is necessitated by the long expression. In order to demonstrate the t ↔ u symmetry at 2-loop, we plot in figure 2, the real and imaginary parts of F (2) fin as a function of x for different choices of cos θ, where θ is the angle between one of the composite operators under consideration and one of the initial scalars in their center of mass frame. Without loss of generality, we set q 2 1 = 100 GeV and q 2 2 = 150 GeV for numerical purpose. The behaviour of the amplitude is clearly seen to be invariant under cos θ → − cos θ, which reflects the t ↔ u symmetry in the final 2-loop expression. In addition to exhibiting the universal infrared structure, this symmetry serves as a very strong check on our calculation, in particular on the finite parts of the FF.

Finite remainders
The iterative structure of form factors in terms of Catani's IR subtraction operators up to two loops level is a result of factorisation at long distances. The observed universality JHEP05(2019)066 of these IR poles lead to exponentiation of these divergences whose coefficients are controlled by universal anomalous dimensions. In the context of the BDS ansatz [3], both the IR divergences and finite parts can be exponentiated and the exponents are expressed in terms of their one-loop counter parts multiplied by a function that depends on the anomalous dimensions such as cusp and collinear ones. For maximally helicity violating (MHV) amplitudes, this iterative structure of the finite terms breaks down starting from two-loop six-point [79,80] and the deviation is captured through a quantity called finite remainder function.
In our case of study, the IR divergences exponentiate, thanks to their iterative structure but not the finite terms. Hence the corresponding finite remainder functions are always non-zero. For the two-point half-BPS operator, the authors of [4] have obtained the finite reminder function up to two loops. Later it was extended to cases with more than two external states [9,27]. Following BDS conjecture, the finite reminder function at two loops is defined as . Note that f (2) ( ) and C (2) are independent of operators as well as external states which is a consequence of the universality of IR divergences. In this article, we compute the finite remainder function for the amplitude under consideration at two loops. The O( 0 ) term of the R (2) , eq. (7.1), can be written as where τ (k) represents the transcendentality of degree k, defined in the introduction, section 1 and R (2),τ (k) represents the terms having transcendentality k. Due to the large size of the terms, we present only the lower transcendental weight terms. The full result can be found from the supplementary material supplied with the arXiv and JHEP submission as well. The lower transcendental terms read as R (2),τ (3) = −2iπ 3 +π 2 (4 log(l 1 )+2 log(l 3 )) , t+u log(l 11 ) +6 log 2 (l 1 ) +12 log(l 1 ) log(l 3 )− 12n 3 d 1 log(l 3 ) log(l 6 )+ 12n 3 d 1 log(l 5 ) log(l 6 )+ 24u t+u Li 2 l 1 l 3 l 11 where the n i and d i are defined in eq. (6.10). Finite remainders are also checked to obey the t ↔ u symmetry. Owing to a large number of cancellation among all the terms, the τ (3) terms of the finite reminder turns out to be surprisingly small compared to other terms.
In the next section, we discuss the behaviour of highest transcendentality terms.

Behaviour of leading transcendental terms
It is an observed [5, 8, 9, 18, 38-44, 81, 82], albeit unproven fact that certain kinds of scattering amplitudes in N = 4 SYM exhibit uniform transcendentality (UT) i.e. those can be expressed in terms of polylogarithmic functions of uniform degree 2L with constant coefficients, where L denotes the loop order. However, this property is no longer true for Sudakov [11,26] or three point [27] form factor of non-protected operators. Naturally, it is expected that the four point amplitude involving a half-BPS and Konishi operator would involve highest (2L) as well as lower (<2L) transcendental terms. Our explicit calculation to two loops shows, eq. (6.7) and (6.8), it is indeed true. Moreover, looking at the one loop finite FF in eq. (6.7) reveals that the coefficients of the highest transcendental terms are indeed simple numerical constant apart from an overall kinematic factor. This is a reflection of the very simple structure of the amplitude in N = 4 SYM.
In [8,11], it was observed that the highest transcendental part of the QCD form factors match exactly with those of half-BPS and Konishi in N = 4 SYM. In spirit of that we intend to examine if a similar thing happens in our case. So, we compare the finite parts of the form factors, F (1) fin and F (2) fin in the limit of same virtuality of the singlet operators with the corresponding quantities of di-Higgs boson production through gluon fusion [31]. A complete mismatch is found, neither highest nor the lower transcendental terms match to each other.
In the limiting case, when masses of the two massive particles become equal, i.e. q 2 1 = q 2 2 = M , the resulting expression is checked against two different four-point amplitude of double O hBPS and O nBPS [30] operators. Comparison with one loop FF of double O hBPS reveals almost none of the terms match whereas for double O nBPS , almost all the highest transcendental terms match exactly but Li 2 (−X) , log 2 (X) and π 2 : fin, KK is the corresponding finite part of the 1-loop FF for two SU(N) scalars to two massive particles of equal masses described by Konishi operator. The subscript τ (2) implies the transcendental 2 terms. The variable X and Y are defined through s = M 2 (1 + X) 2 /X and t = −M 2 Y . At this point it is not fully clear why for certain cases the highest JHEP05(2019)066 transcendental terms match, we need to explore more cases to reveal the underlying reasons behind the coincidence of this kind.

Conclusions and outlook
In this article, we report the very first calculation of two-loop four-point amplitude involving two SU(N) colored massless scalars and two massive particles with different masses represented through a half-BPS and Konishi operators in N = 4 SYM. Following the Feynman diagrammatic approach and using the state-of-the-art techniques the results of the form factors are presented up to weight 4 terms expressed in terms of log , Li n and Li 2,2 functions. Supersymmetry (SUSY) preserving regularisation scheme, the modified dimensional reduction, is employed in order to regulate the ultraviolet divergences arising from the compositeness of the unprotected operator and infrared singularities resulting from massless particles in the theory. Upon verifying the universal behaviour of the infrared divergences, we show that the results contain not only the highest transcendental terms but also all the lower ones. A detailed comparison, particularly, of the highest transcendental terms with several other four-point amplitudes in N = 4 SYM and QCD is performed which does not provide any decisive pattern. More studies [30] of similar kind would be required to shed light on the structures of massive amplitudes.