A Integrals in terms of modified Bessel functions

Finite-volume effects due to spatially non-local operators


Spatially non-local matrix elements are useful lattice-QCD observables in a variety of contexts, for example in determining hadron structure. To quote credible estimates of the systematic uncertainties in these calculations, one must understand, among other things, the size of the finite-volume effects when such matrix elements are extracted from numerical lattice calculations. In this work, we estimate finite-volume effects for matrix elements of non-local operators, composed of two currents displaced in a spatial direction by a distance . We find that the finite-volume corrections depend on the details of the matrix element. If the external state is the lightest degree of freedom in the theory, e.g. the pion in QCD, then the volume corrections scale as , where is the mass of the light state. For heavier external states the usual form is recovered, but with a polynomial prefactor of the form that can lead to enhanced volume effects. These observations are potentially relevant to a wide variety of observables being studied using lattice QCD, including parton distribution functions, double-beta-decay and Compton-scattering matrix elements, and long-range weak matrix elements.

parton distribution functions, lattice QCD, finite volume

I Introduction

One of the fundamental goals in theoretical nuclear physics is the prediction of hadron structure from first-principles calculations based on the underlying gauge theory of the strong nuclear force, quantum chromodynamics (QCD). Quarks and gluons, the degrees of freedom of QCD, are confined into color-singlet bound states that are observed experimentally. The internal structure of these hadrons, encoded in the spatial, momentum and angular momentum distributions of the constituent quarks and gluons, is inherently nonperturbative and poorly understood. Forthcoming data from the 12 GeV upgrade at JLab Dudek et al. (2012) and a future electron-ion-collider Accardi et al. (2016) will provide experimental insight into hadron structure, but a complete understanding of the experimental data requires a concomitant improvement in our theoretical understanding.

Observables related to hadronic structure are most naturally formulated using light-cone quantization, and this provides a serious challenge for lattice QCD, which is necessarily restricted to Euclidean-signature correlation functions. As a specific example, parton distribution functions (PDFs), which capture the distribution of the longitudinal momentum of a hadron among its constituent quarks and gluons, are defined via light-like separated fields and thus cannot be directly accessed in a Euclidean spacetime, where defines a single point rather than a cone. In the past two decades, various ideas have been proposed to overcome this challenge, and thereby calculate PDFs and similar observables from lattice QCD Aglietti et al. (1998); Abada et al. (2001); Detmold and Lin (2006); Braun and Mueller (2008); Ji (2013); Radyushkin (2017a, b); Chambers et al. (2017); Hansen et al. (2017).

Although the details differ, these methods generally require the evaluation of matrix elements of non-local operators, frequently using hadronic states with high momentum. Preliminary results for several of these ideas have now appeared Lin et al. (2015); Alexandrou et al. (2015); Chen et al. (2016); Alexandrou et al. (2017a); Lin et al. (2017a); Chen et al. (2018); Orginos et al. (2017); Chambers et al. (2017); Ma and Qiu (2018); Bali et al. (2018); Alexandrou et al. (2018a) (see Ref. Lin et al. (2017b) for a recent review), but an understanding of all systematic uncertainties is not yet feasible. In general, the systematic uncertainties associated with such calculations include discretization effects, which may be significant for high-momentum states, uncertainties associated with the momentum of the hadron Chen et al. (2017a); Bali et al. (2016), truncation errors arising from perturbative renormalization or matching Constantinou and Panagopoulos (2017); Alexandrou et al. (2018b); Radyushkin (2017c), and finite-volume effects.

In addition to systematically reducing computational uncertainties, it is crucial to shore up the theoretical foundations of the approach. Here, significant progress has been made in understanding issues associated with the renormalization of the Wilson line operator Ishikawa et al. (2016); Monahan and Orginos (2017); Ji et al. (2017); Ishikawa et al. (2017); Green et al. (2017); Alexandrou et al. (2017b); Chen et al. (2017b, 2018), the existence of factorization theorems Zhang et al. (2018); Izubuchi et al. (2018), and the role of the Euclidean signature in lattice calculations Briceño et al. (2017).

In this work we study the finite-volume artifacts associated with spatially-displaced currents. We do so by perturbatively studying a matrix element of a product of two currents in a toy theory with one light degree of freedom (corresponding to the pion in QCD) and one heavy (corresponding to a nucleon or heavy meson). Proposals to extract hadronic structure observables from lattice-QCD calculations using these types of operators, in place of those defined with a Wilson line, appeared in Ma and Qiu (2018); Braun and Mueller (2008), with the first numerical results for pion distribution amplitudes presented in Bali et al. (2018).

To perform this type of finite-volume analysis for Wilson-line based operators would require a low-energy representation for these, more complicated, non-local objects Chen and Ji (2001); Kivel and Polyakov (2002, 2008); Diehl et al. (2005). One possible avenue is to build an effective field theory based on the method of Ref. Green et al. (2017), in which an auxiliary heavy-quark field enables one to write gauge-invariant extended operators in terms of products of quark bilinears. However, this goes beyond the scope of the present work and we focus our attention on composite bilinear currents, whose hadronic representation is more straightforward.

For lattice calculations of hadronic masses, and other properties defined through local operators, finite-volume effects lead to corrections of the form  Lüscher (1986); Colangelo et al. (2003); Colangelo and Durr (2004); Beane (2004); Colangelo and Haefeli (2004); Beane and Savage (2004); Bijnens et al. (2014), where is the linear extent of the cubic spatial volume.4 Numerical data is expected to be described by this leading exponential form, with a power-law prefactor, provided one performs the calculation with asymptotically large volumes, . (In practice is generally sufficient.) However, in the presence of a second infrared (IR) length scale, such as the current separation in a spatially-extended operator, one naturally expects the finite-volume effects to be modified.

For matrix elements of composite currents, we show that finite-volume effects take the form


where the left-hand side represents the difference between the finite-volume matrix element (obtained via lattice QCD) and its infinite-volume limit. The external states here are zero-momentum, single-particle states, labeled by their physical mass, ; is the displacement vector within the composite current and is its magnitude. To derive this result, we assume .

The right-hand side of Eq. (1) gives the leading finite-volume effects. We focus on two terms, one scaling with the mass of the external state and the other with the mass of the lightest degree of freedom. In the case where these two coincide the first term scales as and is expected to dominate the volume effects once becomes a non-negligible fraction of . By contrast if , as in the case of a nucleon or heavy meson, then the second term dominates. Both terms have polynomial prefactors, denoted and , with terms scaling as . Such factors can also have a significant impact on the size of finite-volume corrections if is non-negligible compared to the box size. Finally, the ellipsis in Eq. (1) represents subleading exponentials.

Figure 1: Finite- versus infinite-volume behavior of non-local matrix elements. These plots were made using the tree-level result derived in the next section with pions as the external states. Subfigure shows how, as is varied, the finite-volume matrix element develops large deviations from its infinite-volume counterpart. For , indicated by the shaded region, high energy scales are sample so that the effective field theory is expected to break down. Subfigure shows the fractional difference between finite- and infinite-volume matrix elements, . Finally, (c) shows the finite-volume matrix element, , as a function of , together with its infinite-volume limit. For fixed , finite-volume effects for decay with as . These are enhanced by a prefactor relative to the typical, , finite-volume effects for local operators.

To better understand these volume effects we note that, while the infinite-volume matrix element generally decays as a function of , its finite-volume counterpart must be periodic, with periodicity . Thus, as we illustrate in Fig. 1(a), the differences between the finite- and infinite-volume objects become arbitrarily large as approaches . We are interested in the onset of this effect for . As we show in Fig. 1(b), in the case where the external state is light, the finite-volume effects grow exponentially with , exhibiting deviations for when . For this same volume, Fig. 1(c) shows that these volume effects can be removed by fitting to a decaying exponential in at fixed . We stress that the details of these features hold only for matrix elements built from products of local currents.

The origin of periodicity for matrix elements built from products of local operators is straightforward: Given that the quark fields and the gauge links are periodic in all spatial directions, the same must be true for any local current constructed from these fields and links. The periodicity property , with a unit vector in the , or direction, is then directly inherited by matrix elements constructed from products of such currents at different locations.

However, this argument does not hold for non-local quark bilinears connected by Wilson lines, i.e. the type of non-local operator used to extract quasi- and pseudo-PDFs. Defining as the straight Wilson line connecting the points and , given by


one can construct a gauge-invariant non-local operator by contracting this with quark and antiquark fields at and , respectively. The quark fields and the gauge fields are periodic in the spatial directions, but for fixed x there is no periodicity in the coordinate . In particular, wrapping around the torus times gives


where the factor in parenthesis on the right-hand side breaks the naive periodicity relation, .

On the one hand, this additional factor may lead to matrix-elements of this operator being closer to their infinite-volume counterparts than in the case of products of currents that satisfy -periodicity. On the other hand, the fact that the boundary conditions are felt by both the quark fields and individual links leads us to expect that large values of may generate enhanced volume effects in this case as well.

As we describe in detail in Sec. II below, we derive Eq. (1) and Fig. 1 using a toy theory with two relativistic scalar particles. A rigorous demonstration that the scaling also holds in QCD would require first defining a specific matrix element, then developing a low-energy effective-field-theory description (based in chiral perturbation theory) and finally calculating finite-volume corrections. However, since our result only relies on the appearance of scalar propagators with the light particle mass, together with the scale that characterizes the operator non-locality, we expect that a more realistic description would change only the detailed form of and , and not the overall exponential behavior.

To close the introduction, we comment on a number of other examples in which finite-volume effects on non-local matrix elements have already been discussed in the literature.

The authors of Ref. Christ et al. (2015) describe formalism for removing finite-volume effects in neutral kaon mixing. The starting point is a matrix element reminiscent of that considered here, defined with external kaon states and two insertions of the weak hamiltonian. In contrast to the matrix elements in this study, however, the currents are also separated in Euclidean time. By summing over time slices, the authors demonstrate how to identify a finite-volume version of . In a second step, the leading finite-volume effects are removed using a generalization of the Lellouch-Lüscher formalism Lellouch and Lüscher (2001). The step of identifying the finite-volume version of relies on picking out a single term in the temporally-summed correlator. This term in isolation has power-law volume effects associated with on-shell intermediate states, i.e. effects parametrically larger than those identified in the present study. It is these volume artifacts that are corrected via the extended Lellouch-Lüscher formalism.

In a different application, in Ref. Hansen et al. (2017), one of us considered an approach for extracting total decay and transition rates from temporally displaced currents with single-particle external states. The method requires estimating a smeared-out inversion of the Laplace transform, for example by using the Backus-Gilbert method Backus and Gilbert (1968, 1970); Press et al. (2007). As discussed in detail in Ref. Hansen et al. (2017), this smearing suppresses finite-volume effects in the target observable. The infinite-volume observable must then be extracted by identifying an optimal trajectory in the two-coordinate plane of box-size, , and smearing width . The enhanced volume effects identified here will likely influence this optimal trajectory, but the detailed consequences are not clear and are the subject of future work.

The remainder of this article is organized as follows. In Sec. II we explain the set-up of our calculation, including the detailed definition of the toy theory and the external currents. We then summarize the general framework for calculating finite-volume effects in non-local matrix elements and apply this to the tree-level diagram of Fig. 2(b). In Sec. III we provide a detailed study of the finite-volume effects in one-loop diagrams, focusing on the case where the external state in the matrix element is a heavy particle, for example a nucleon. We then extract the large-volume scaling of these diagrams in Sec. IV, and deduce the result summarized by Eq. (1). In Sec. V we briefly conclude and outline possible future work. Technical details of certain functions used in the analysis are discussed in the appendix.

Ii Set-up and simple example

Figure 2: (a) Feynman rules for the EFT described in Sec. II. The dashed lines denote the lighter particle, , and the solid lines denote the heavier particle, . (b) The leading order contribution to the matrix element with external states. (c) Contact interactions that may arise in EFTs.

To determine finite-volume corrections to field-theoretic observables analytically, it is often useful to work with a low-energy effective theory (EFT), e.g. chiral-perturbation theory (PT) for QCD. This is natural because quarks are confined, and only the hadrons (the degrees of freedom in PT) can propagate long distances and feel the finite-volume boundary conditions.

As a first step, in this work we study a toy theory that is expected to capture the basic scaling of the finite-volume corrections. We present a formal expression [Eq. (14) below] that gives the finite-volume effects of spatially non-local currents from an arbitrary Feynman diagram. Using this result, we determine the finite-volume corrections from the tree-level diagram shown in Fig. 2(b) and, in the following two sections, consider the role of one-loop corrections in volume effects.

We consider a theory with two scalar particles, and , with physical pole masses satisfying . Here plays the role of the pion in QCD and that of the nucleon or a heavy meson. Using baryon and heavy meson PT Weinberg (1991, 1990); Wise (1992) as an inspiration, these states couple via a vertex whose Feynman rule is given in Fig. 2(a). In QCD, pions are pseudo-Goldstone bosons, and thus have derivative couplings to heavier particles. However, the exact form of these couplings does not change the leading exponential behavior of the theory, so here we only consider a momentum-independent coupling, labeled . Similarly, we restrict attention to scalars, as details associated with spin and isospin are not expected to change the overall scaling of the finite-volume corrections.

The and states couple to a renormalized external current, given by


that generates the additional Feynman rules shown in Fig. 2(a).

The renormalization factors and , are inherited from the mass terms in the Lagrangian (with the scheme fixed by ). The three-point renormalization is inherited in a similar way from , with the scheme that the amputated three-point function equals (its tree-level value) when all . A similar scheme can be used for the term, although the coupling does not appear in the Lagrangian.5 Finally, the ellipsis in Eq. (4) stands for terms with higher orders in and .

The final step is to define a power-counting scheme for the theory. We take . As we consider matrix elements with two insertions of the local current, leading-order (LO) contributions scale as and next-to-leading order (NLO) as .

We are now ready to set-up our general approach for determining finite-volume effects in Feynman diagrams contributing to matrix elements of spatially non-local operators. We define the infinite-volume matrix element as


where is a single-particle state with momentum p, either a or a to be specified below. Now note that any diagram, , contributing to this quantity can be written as


where we have introduced the shorthand


In Eq. (6), is the standard integrand, constructed according to the usual Feynman rules, and containing all couplings and symmetry factors. The separation of the factor simplifies the relation to Euclidean-signature quantities. In particular, from the Wick rotation we find


where is the usual integrand that one would construct with Euclidean Feynman rules.

As an example, for the leading order diagram show in Fig. 2(b), the Minkowski integrand is


and Wick rotation gives


consistent with the usual Feynman rules. Each loop introduces a factor of to the Minkowski integrand, but in our convention this is factored out to preserve as defined in the two signatures.6

We now give our general expression for the finite-volume effects from spatially non-local currents. From the Poisson summation formula follows that the finite-volume residue for any given diagram can be written as


where and the notation under the sum indicates that the only point omitted is when all three vectors vanish. Introducing we reach a very compact form for the residue


Heuristically, M parametrizes the space of images that enforce the finite-volume boundary conditions and the smallest nonzero values (the nearest neighbors) give the dominant finite-volume effects. These can be in the n direction, corresponding to effects on the Fourier transform from q to , as well as the directions, corresponding to finite-volume effects within the diagram.

Returning again to the leading-order diagram, Fig. 2(b), and using the Euclidean form of Eq. (14), we reach


In Appendix A.1 we review standard tools for writing these integrals in terms of modified Bessel functions. In particular we find it convenient to define




where in the second step we used that and therefore that . We comment that, up to a proportionality constant, this is just the finite-volume meson propagator in position space. [See also Ref. Cherman et al. (2017).]

Keeping only the term, we find that the dominant finite-volume effect scales as


where the arrow indicates the asymptotic limit.7 The key scaling is given by stripping off the coupling and other prefactors


This is the main result of this section and corresponds to the first term on the right-hand side of Eq. (1).

Note that the infinite-volume prediction of this diagram can also be read from this expression by replacing with . This implies, in particular, that the diagram diverges in the limit , as illustrated in Fig. 1. However, given that we consider a toy EFT that is necessarily written in terms of hadrons, one cannot expect to accurately describe behavior of physical amplitudes for short distances of the scale . We thus require , to ensure that the finite- and infinite-volume matrix elements are well described by the EFT.

We close this section by commenting on the diagram shown in Fig. 2(c). As this is only a contact interaction it introduces no finite-volume effects to the matrix element. To understand this in detail requires including a renormalization factor for the product of currents, to accommodate divergences when the two overlap. We have studied this to ensure that no unexpected issues arise.

Iii Beyond leading order

We now turn our attention to the case that the heavy particle, denoted , appears in the external state. The leading-order contribution to this matrix element is given by Diagram 2(b), with the dotted propagator replaced by a propagator, and the result is Eq. (19) with the substitutions and ,


We take and , implying that effects of can be safely ignored. In particular, as we show in this section, the leading finite-volume effects the for the matrix element with a external state are generated by the next-to-leading-order corrections to this result.

In Fig. 3 we show the one-loop corrections to the matrix element. Here we omit diagrams that give finite-volume corrections to the external states. These give corrections of and, since we are interested in volume effects enhanced by the non-locality scale , can be safely dropped. We first derive general integral expressions for the diagrams in Fig. 3, restricting attention to the case where the external particle is at rest in the finite volume and highlighting Diagram 3(a) as a specific example. Generally, the integrals that arise in evaluating these diagrams cannot be carried out analytically. To study their asymptotic behavior, in Sec. IV we separate the expressions into analytic parts that dominate the volume scaling together with numerically determined functions that are slowly-varying and give only subleading corrections to the scaling.

Figure 3: Next-to-leading-order contributions to the matrix elements when the external state is the heavy particle. Corrections to the wavefunction renormalization of the external states are not shown.

Diagram 3(a):

We illustrate the calculation of the one-loop diagrams in Fig. 3 by outlining the derivation of finite-volume effects for Diagram 3(a). The calculation of Diagrams 3(b) to 3(h) proceeds in a similar fashion but Diagrams 3(i) and 3(j) require special treatment, as we discuss below and in Appendix A.3. From Eq. (14) we identify the finite-volume residue for 3(a) as


where the notation below the summation indicates that only the term is omitted from the sum. We separate the and dependence by shifting to reach


Next we use a Feynman parameter to reduce the second integral. Starting with the identity


we shift to reach




In the second step we have set the Euclidean external momentum on-shell, .

At this stage we have written the loop in terms of products of two integrals of the kind given in Eq. (16). Substituting the definition of then gives


where we have shifted the summed integer vectors. Taking the external state to be at rest in the finite-volume, i.e. setting , then gives


Diagrams 3(b) to 3(h):

This set of diagrams is amenable to the same approach as Diagram 3(a). In Appendix A.2, we present a simple generalization of the technique presented above for Diagram 3(a) that allows for a rapid derivation of the finite-volume effects for these diagrams. The results for are


The key feature for these diagrams is that one can factorize the dependence on the current momentum, , from that on the internal loop momentum, . In all cases this results in two modified Bessel functions, corresponding to the two momenta after an appropriate shift has been performed. Note that the sum of the indices on the two Bessel functions always corresponds to the number of internal propagators.

Diagrams 3(i) and (j)

These two diagrams cannot be factorized into two separate momentum integrals and must be studied using a different approach. In Appendix A.3, we evaluate these diagrams and place upper bounds on their values. We demonstrate that the finite-volume artifacts associated with these are smaller than those for 3(a) to (h). As we are only interested in the dominant finite-volume effects, we ignore the contributions from Figs. 3(i) and (j) from here on.

Iv Asymptotic behavior

In this section we study the asymptotic behavior of Eqs. (28) through (35). As mentioned above, we assume that and ignore corrections that decrease with the volume as or more rapidly. As mentioned above, the matrix element must be periodic with periodicity . Thus, as approaches , the finite-volume effects become arbitrarily large [see also Fig. 1 above]. Here we are not directly interested in this regime of extreme volume effects but rather in the region of with . This motivates us to take the asymptotic forms of the Bessel functions, i.e. to take the arguments as large.

Combining the asymptotic form of the Bessel functions with the definition of Eq. (16), we find


Given this exponential suppression, terms with n chosen to minimize will dominate the sum. In addition, terms scaling as , i.e. with the mass of the heavier particle, will be highly suppressed and we drop such contributions throughout.

In Eqs. (28) through (35), only appear. We thus give their explicit forms for convenience


where we use to indicate that the two sides agree up to terms suppressed by additional powers of .

The asymptotic forms of the one loop diagrams can be determined using a similar approach to that for the leading order diagram, Eq. (21). We identify the dominant terms in the sums over n and m assuming with . The only additional subtlety is that the integrals over Feynman parameters are found to be numerically dominated by . Factoring out this dependence we reach the following



Figure 4: Plot of the function vs for set to the nucleon-pion mass ratio. The figure shows that the functions are slowly varying compared to the corresponding exponentials and powers of appearing in Eqs. (40) to (47). Thus the leading scaling is given by treating these functions as constant.

As we show in Fig. 4, is a slowly varying function of its argument. Thus the leading scaling can be read from the given expressions. We deduce that, in all cases, the finite-volume effects are suppressed by a factor of . In particular, the dominant finite-volume contributions come from Diagrams (a), (b) and (e) with the leading effect for small driven by the factor appearing in Diagram (a).

V Summary

We have presented the first steps towards understanding the finite-volume artifacts that arise in matrix elements of spatially non-local operators. These operators are relevant for a wide variety of observables being studied using lattice QCD, including parton distribution functions, double-beta-decay and Compton-scattering matrix elements, and long-range weak matrix elements. In particular, matrix elements of products of spatially-separated currents represent one approach to determining hadron structure directly from lattice QCD Ma and Qiu (2018); Braun and Mueller (2008); Bali et al. (2018).

We considered a toy model involving two scalar particles, one analogous to the pion in QCD and one analogous to the nucleon or a heavy meson, and determined the finite-volume matrix elements of these states with two spatially-separated scalar currents at one-loop in perturbation theory. As expected, we find that these matrix elements are contaminated by larger finite-volume artifacts than is the case for matrix elements of local operators. The results are summarized in Eq. (1). There are two terms that potentially dominate the finite-volume artifacts, one scaling with the mass of the external state and the other with the mass of the lightest degree of freedom. When these two coincide, the dominant finite-volume correction scales as , provided is a non-negligible fraction of . By contrast, if the external state is significantly heavier than the lightest particle, as in the case of a nucleon or heavy meson, then the leading finite volume artifacts scale as . In both cases these exponential factors have polynomial prefactors, including terms scaling as that could have a significant impact on the size of finite-volume corrections.

Future extensions of this work include implementing the approach for specific QCD matrix elements using PT, including flavor-changing currents, and more carefully studying the dependence on the external state momenta (especially at next-to-leading order). A more involved extension would be to apply the approach to operators involving Wilson lines, such as those relevant for determining quasi- and pseudo-PDFs. This requires a PT-based representation of such operators Chen and Ji (2001); Kivel and Polyakov (2002, 2008); Diehl et al. (2005), potentially built around the methods of Ref. Green et al. (2017). Finally, here we have only considered periodic boundary conditions. Previous works Briceno et al. (2014); Cherman et al. (2017); Korber and Luu (2016) have shown that particular choices of twisted boundary conditions Bedaque (2004); Byers and Yang (1961) can be used to reduce the size of exponentially suppressed finite-volume artifacts. This may also prove useful in reducing volume corrections for matrix elements of spatially non-local operators.


R. A. B. and J. V. G. acknowledge support from U.S. Department of Energy contract DE-AC05-06OR23177, under which Jefferson Science Associates, LLC, manages and operates Jefferson Lab. J. V. G. is supported in part by the U.S. Department of Energy through DOE contract No. DE-SC0008791 and also through the JSA Graduate Fellowship Program. C. J. M. is supported in part by the U.S. Department of Energy through Grant Number DE-FG02-00ER41132. The authors would like to thank Robert Edwards for persistently asking questions that inspired this project. The authors would also like to thank Alberto Accardi, Martin Savage, Joseph Karpie, and David Richards for helpful conversations during the course of this work.

Appendix A Integrals in terms of modified Bessel functions

In Sec. III we have shown that complicated diagrams can be written as integrals of products of modified Bessel functions. Although these integrals are well documented and derived in the literature (see, for example, Ref. Bijnens et al. (2014)), we review the derivation of the necessary functions in this appendix.

a.1 Derivation of Eq. (16)

We begin by deriving Eq. (16), the identity relating the function


to the modified Bessel function, .

Beginning with the definition of the function


we observe


Next we complete the square in the four-vector, , to write


where . Performing the integral over then gives


where we set from here on.

Finally we perform the variable substitution to reach


Here we have used the fact the and are symmetric and antisymmetric respectively, and have introduced the modified Bessel function, . Note that .

a.2 Loops in terms of

Figures 3(a)-(h) can be written as integrals over products of defined in Eq. (49). Here we show our general method for doing this for all integrals of the form


First one inserts a Feynman parameter integral to combine the denominators


where we performed the variable transformation , used the on-shell condition for the external states, , and also substituted .

Using the functions defined in Eq. (56), we arrive at


From this, it is straightforward to arrive at the expressions given for Figures 3(b)-(h) in Eqs. (29)-(35).

a.3 Detailed calculation of diagrams (i) and (j)

Diagram 3(j):

The contribution of Diagram 3(j) is given by