Form factors of two-hadron states from a covariant finite-volume formalism

In this work we develop a Lorentz-covariant version of the previously derived formalism for relating finite-volume matrix elements to $\textbf 2 + \mathcal J \to \textbf 2$ transition amplitudes. We also give various details relevant for the implementation of this formalism in a realistic numerical lattice QCD calculation. Particular focus is given to the role of single-particle form factors in disentangling finite-volume effects from the triangle diagram that arise when $\mathcal J$ couples to one of the two hadrons. This also leads to a new finite-volume function, denoted $G$, the numerical evaluation of which is described in detail. As an example we discuss the determination of the $\pi \pi + \mathcal J \to \pi \pi$ amplitude in the $\rho$ channel, for which the single-pion form factor, $F_\pi(Q^2)$, as well as the scattering phase, $\delta_{\pi\pi}$, are required to remove all power-law finite-volume effects. The formalism presented here holds for local currents with arbitrary Lorentz structure, and we give specific examples of insertions with up to two Lorentz indices.


I. INTRODUCTION
In recent years, interest in hadron spectroscopy has increased significantly, primarily due to various experimental discoveries of unconventional excitations. 1 This has led to an abundance of theoretical proposals as to the underlying nature of the unexpected states.Possible explanations range from multi-hadron molecules to compact multi-quark configurations, to kinematic singularities arising from specific Feynman-diagram topologies [1][2][3][4][5].In many cases, experimental data alone is not sufficient to distinguish between available explanations, and thus many questions remain unresolved.
In some cases, theoretical calculations can provide access to experimentally-unavailable quantities that may shed light onto the structure of the QCD spectrum.With this goal in mind, in this work we present a detailed framework that will allow for the rigorous lattice-QCD calculation of transition amplitudes, mediated by electroweak or other external currents, involving two hadrons each in the initial and final states.We abbreviate our process of interest by 2 + J → 2, where each bold-faced 2 counts the QCD-stable hadrons in the state and J represents a generic, local external current.
The approach discussed here is based on prior formalism developed by two of us in Ref. [6]. 2 In Sec.II we present a slightly modified version of this formalism in which all infinite-volume quantities are Lorentz covariant and the singleparticle matrix elements that enter, abbreviated 1 + J → 1, are expressed in terms of standard Lorentz-invariant form factors.After extracting the 2 + J → 2 transition amplitudes, one can proceed to determine form factors as well as distribution functions 3 of bound states or resonances that couple to the asymptotic states.From the form factors and distribution functions, in turn, one can obtain various structural parameters, e.g. the charge or even the gluonic [12,13] radius of a given state.
The primary reason why a non-trivial formalism is required to extract multi-hadron observables from lattice QCD calculations is that the latter are performed in a finite spatial volume, usually a cube defined with periodic boundary conditions on the quark and gluon fields.This complicates the determination of scattering and transition amplitudes, x t 5 n j R O K q 7 h N 6 f l 6 k V e R x H t o X 1 0 h F x 0 h q r o C t V Q H V H 0 g J 7 Q C 3 q 1 H q 1 n 6 8 1 6 n 4 4 W r H x n F / 2 B 9 f E N v H q a E Q = = < / l a t e x i t > < l a t e x i t s h a x t 5 n j R O K q 7 h N 6 f l 6 k V e R x H t o X 1 0 h F x 0 h q r o C t V Q H V H 0 g J 7 Q C 3 q 1 H q 1 n 6 8 1 6 n 4 4 W r H x n F / 2 B 9 f E N v H q a E Q = = < / l a t e x i t > < l a t e x i t s h a x t 5 n j R O K q 7 h N 6 f l 6 k V e R x H t o X 1 0 h F x 0 h q r o C t V Q H V H 0 g J 7 Q C 3 q 1 H q 1 n 6 8 1 6 n 4 4 W r H x n F / 2 B 9 f E N v H q a E Q = = < / l a t e x i t > < l a t e x i t s h a x t 5 n j R O K q 7 h N 6 f l 6 k V e R x H t o X 1 0 h F x 0 h q r o C t V Q H V H 0 g J 7 Q C 3 q 1 H q 1 n 6 8 1 6 n 4 4 W r H x n F / 2 B 9 f E N v H q a E Q = = < / l a t e x i t > FIG. 1. Road map of the formal approach outlined in this work.See also Fig. 2 of Ref. [6].The four red arrows merging together represent how the present approach combines various finite-and infinite-volume information to extract the 2 + J → 2 amplitudes.Analytically continuing these to the resonance-pole location gives a robust, model-independent definition of the resonance form factor.
because there is no simple relation between the finite-volume QCD eigenstates and the asymptotic multi-particle states that arise in the infinite-volume limit of the theory.However, in certain cases, it is possible to derive relations between finite-and infinite-volume observables.These have been implemented with great success to access a wide variety of scattering quantities directly from numerical LQCD calculations.See Ref. [14] for a recent review.The most well-established such relation is that derived by Lüscher in Refs.[15,16] over three decades ago.In these seminal papers he showed how elastic two-particle scattering amplitudes can be extracted from the finite-volume energy spectrum below the lowest lying three-or four-particle threshold.Since then, the idea has been generalized to all possible two-body systems, in particular to multiple two-particle channels built form any number of particle species, including particles with any intrinsic spin [7,[17][18][19][20][21][22][23][24][25][26][27][28][29][30][31][32].These formal ideas, together with significant algorithmic developments, have resulted in a proliferation of scattering amplitudes determined directly from lattice QCD [33][34][35][36][37][38][39][40][41][42][43][44][45][46].A key limitation to the methods currently being implemented is the restriction to two-particle states, but the formal extension to three-particle systems has received significant attention recently and is progressing [47][48][49][50][51][52].
In this work we turn our attention to the prospect of determining 2+J → 2 transition amplitudes from finite-volume matrix elements.This was previously considered in Ref. [6].In contrast to that work, here we restrict attention to kinematics such that only one two-particle channel is open, and take the two particles in the channel to be scalars.In addition, we only consider flavor-conserving external currents, so that the initial and final two-particle states contain the same particles.Just as in Ref. [6], the two-particle states are composed of QCD-stable (pseudo-)scalars.Relaxing these restrictions to provide a Lorentz-covariant formalism for any number of two-particle channels, including those with intrinsic spin, is expected to be straightforward.
As in Ref. [6] in this work we derive a mapping between finite-volume matrix elements of two-particle states and the infinite-volume 2 + J → 2 amplitude.The result is summarized by the flow-chart shown in Fig. 1.We find that, given the following quantities: • the two-particle finite-volume spectrum, • the 1 + J → 1 form factors, • the finite-volume two-particle matrix elements of J , one can systematically constrain the 2 + J → 2 amplitude in the kinematic window in which only the accommodated channel contributes.Our relation requires the generalized Lellouch-Lüscher factors [6,54,60], that enter multiplicatively in the conversion, as well as a new finite-volume function, denoted G, that appears in an additive correction, together with the single-particle form factor as well as the two-to-two scattering amplitude.
A simple limiting case of our result is the one in which the single-particle form factors vanish.In this limit the finite-volume artifacts associated with the G function also vanish and one recovers a Lellouch-Lüscher-like relation in which the correction factor appears twice, once each for the initial and final two-particle states.However, when the single-particle form factors are nonzero, the term containing G is expected to give the dominant finite-volume effects.In particular, the analysis of Ref. [53] showed that, in the case of weak interactions, the finite-volume effects on the ground-state 2 + J → 2 matrix element are given by an expansion in powers of 1/L. 4 In these expressions, the diagrams that appear as our G give 1/L 2 corrections while all other terms contribute with additional powers of 1/L.
The purpose of this work is to improve certain technical details of the formalism and to provide more concrete information on the implementation procedure.We begin by providing a covariant version of the formalism in Sec.II, where we also discuss three key examples involving ππ states.Then, in Sec.III, we explain in detail our approach for evaluating G and outline why this is more challenging than the more-standard finite-volume functions relevant for two-to-two scattering.In addition to the standard threshold singularities, in this section we discuss and illustrate the presence of triangle singularities in the G function.In Sec.IV we conclude and provide an outlook for future studies.Finally we include two appendices: Appendix A gives various details relevant for the derivation of the improved formalism used here.Appendix B includes various technical aspects regarding the evaluation of the finite-volume functions discussed in the main text.

II. COVARIANT REPRESENTATION OF THE FORMALISM
In this section we revisit the formalism derived in Ref. [6] and present a modified form in which all infinite-volume quantities are Lorentz covariant.We focus here only on the final result, and in Appendix A we explain the (minimal) modifications to the original derivation that lead to this new form.
This section is divided into three subsections: In Sec.II A we review the required notation and give the quantization condition, as well as the generalized Lellouch-Lüscher matrix, for two-particle states in a finite volume.Then, in Sec.II B, we provide a full description of our covariant 2 + J → 2 formalism.Finally, in Sec.II C, we consider a handful of specific examples to show how our general expressions reduce for a particular system with specified quantum numbers.

A. Kinematics and the quantization condition
We denote the 4-momentum of the incoming state in the finite-volume frame by P µ i ≡ (E i , P i ) and that of the outgoing state by P µ f ≡ (E f , P f ).The center-of-momentum (c.m.) energies corresponding to these are then given by This also defines the metric used for the Minkowski-signature 4-vector dot products throughout.Following the notation of Ref. [22] we use to indicate quantities defined in either the incoming or the outgoing c.m. frame.As explained below, we often use an i or an f index in addition to the , in order to completely specify the frame.
In this work we accommodate all values of 3-momenta allowed by the periodic boundary conditions, i.e.P i = 2πd i /L and P f = 2πd f /L where d i and d f are 3-vectors of integers.The energies and 3-momenta can differ between the initial and final states due to the momentum carried by the external current, P µ f − P µ i [see Fig. 2(a)].The physical quantities discussed below depend on Lorentz scalars.For most systems we will primarily be sensitive to spacelike values of the momentum transfer, motivating us to introduce which is positive for spacelike P µ f − P µ i .As mentioned above, we restrict our attention here to values of E i and E f such that only a single two-particle channel can propagate.Within the single channel considered, we accommodate both identical and non-identical scalars and allow these to be non-degenerate in the latter case, with physical masses m 1 and m 2 .We assume, however, that the current, J , is flavor conserving so that the same two particles appear in the initial and final state. 5e now turn to the kinematic variables describing individual particles within the two-particle channel.For the remainder of this subsection, take P µ ≡ (E, P) to simultaneously represent the initial and final state 4-momenta.Denoting the 3-momentum of particle 2 (with mass m 2 ) by k, the corresponding on-shell 4-vector is k µ = (ω k2 , k), where with k = |k|.
Next note that, in order to satisfy the specified total 4-momentum (P µ ), particle 1 must carry P µ − k µ ≡ (E − ω k2 , P − k).Thus, for general k, one cannot simultaneously require that the particle momenta sum to P µ and that particle 1 is on shell.The latter holds only when the temporal component, E − ω k2 , is equal to To better understand when the on-shell condition ) as a boost matrix with boost velocity β β β ≡ P/E.We then define and observe where k ≡ |k |. 6 By contrast, the 4-momentum of particle 1 boosts to where we have used [Λ β β β ] µ ν P ν = (E , 0).We deduce that the c.m. frame on-shell condition is E = ω k1 + ω k2 where The advantage of working in this frame is that the on-shell condition reduces to a constraint on the value of k .In particular the particle pair is on shell if and only if k = q with the latter defined by Finally, the initial or final state-indices must be applied to all of these quantities once the total 4-momentum is associated with a particular state.In particular if we take P µ → P µ i , then the corresponding quantities above become The only quantities that do not inherit a frame index are the finite-volume frame momentum, k, and the corresponding energy, ω k2 .
With this notation in hand we now give the quantization condition for two scalar particles in a finite volume.This is written as a determinant condition involving the on-shell two-particle scattering amplitude, M(P 2 ), represented as a diagonal matrix in angular-momentum space.For any fixed values of P and L, the finite-volume energy spectrum is given by solutions to [22,28,32] det[M(P 2 ) −1 + F (P, L)] = 0 . ( This holds in the region 0 < E < E th , where E th is the energy of the lowest-lying multi-particle threshold that we ignore, which could be a two-, three-, or four-particle threshold.The relation holds up to neglected corrections of the form e −mL , where m is the physical mass of the lightest degree of freedom in the theory.
+ F q 0 l p 5 g 5 h j 9 y P n 8 A D K a N k g = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " 7 e B S i K + F q 0 l p 5 g 5 h j 9 y P n 8 A D K a N k g = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " 7 e B S i K + F q 0 l p 5 g 5 h j 9 y P n 8 A D K a N k g = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " 7 e B S i K + F q 0 l p 5 g 5 h j 9 y P n 8 A D K a N k g = = < / l a t e x i t > r z 5 3 5 h u P e j m g 8 D L e 7 8 f y X t h y q j S j v N t L S 2 v r K 6 t r z 5 3 5 h u P e j m g 8 D L e 7 8 f y X t h y q j S j v N t L S 2 v r K 6 t r z 5 3 5 h u P e j m g 8 D L e 7 8 f y X t h y q j S j v N t L S 2 v r K 6 t r z 5 3 5 h u P e j m g 8 D L e 7 8 f y X t h y q j S j v N t L S 2 v r K 6 t x + z 0 i W r 6 j k A f 2 B 9 / g B O D p q g < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " O + x + z 0 i W r 6 j k A f 2 B 9 / g B O D p q g < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " O + x + z 0 i W r 6 j k A f 2 B 9 / g B O D p q g < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " O + The black circles depict the → 2 scattering amplitude in the absence of the current.The crossed circles represent various couplings of the external current.Those with two hadronic external legs correspond to the standard 1 + J → 1 matrix element, while those with four external hadrons represent diagrams that are two-particle irreducible in the channel carrying the total momentum.The solid lines denote fully-dressed propagators of the low-energy degrees of freedom (the hadrons).In the second line we separate this into a contribution with on-shell singularities together with the divergence-free amplitude.(d) The diagrammatic representation of the new finite-volume function Gµ 1 ...µn , defined in Eq. ( 27).
The precise definition of M(P 2 ) is given by where δ ( ) (q ) is the scattering phase shift, and ξ = 1/2 if the particles are identical and ξ = 1 otherwise.The remaining ingredient is F (P, L), a matrix of finite-volume geometric functions defined by where In Appendix B 2 we review an efficient method to evaluate this, based on analytic expressions for the integrals defined using the cutoff functions introduced in Ref. [22].
To close this subsection, we introduce one additional finite-volume matrix This object, introduced in this form in Ref. [60], is the generalization of the Lellouch-Lüscher factor [54] that relates finite-volume matrix elements of two-particle states to the corresponding infinite-volume decay and transition amplitudes.
B. Relating finite-volume matrix elements with 2 + J → 2 transition amplitudes We are now ready to present our improved finite-volume formalism.The approach that we advocate here differs from that of Ref. [6] in two key ways.
First, the separation of singularities, required to disentangle finite-volume effects in the 2 + J → 2 amplitude, is done here using Lorentz invariant poles of the form 1/(k 2 −m 2 ).In the previous work we instead used As long as we consistently modify the pole form everywhere, it turns out that either choice is valid.The advantage of the present approach is that it ensures all infinite-volume quantities are Lorentz covariant and also simplifies the form of the new finite-volume function, G, arising due to the triangle diagram shown in Fig. 2(b) and defined explicitly in Fig. 2(d) and in Eq. ( 27) below.
Second, we treat the single-particle matrix elements in a simpler way here than we did in Ref. [6].Our approach requires isolating the 1 + J → 1 matrix element in order to express the finite-volume effects of the triangle diagram [Fig.2(b)].In our previous publication, a complicated scheme was presented to project the matrix element on shell.Though correct, we have come to realize that this procedure is unnecessary.The reason, as we explain in more detail below, is that one can decompose the matrix elements into kinematically determined tensor structures and form factors. Projecting the kinematic factors on shell is unnecessary, and removing this step gives a simpler approach that leads to the same infinite-volume observables.
We begin by introducing notation for the physical 2 + J → 2 matrix element that we are after [see also Fig. 2(c)] Here the initial and final states are standard two-particle asymptotic states with the usual relativistic normalization convention and J µ1•••µn (0) is a generic local current insertion.As is shown in Fig. 2(a), the initial state is built from particle 1 [with mass m 1 and on-shell momentum . The final state is built from the same pair, now carrying momenta P f − k and k .Following the discussion of the previous sub-section, the c.m. frame 3-momenta are denoted by k i and k f and have magnitudes equal to q i and q f respectively.The label "conn."emphasizes that only fully connected diagrams, those shown in Fig. 2(c), contribute to the definition of the 2 + J → 2 matrix element.
A consequence of the on-shell constraints is that, once total energy and momenta are fixed, the two-particle states only have directional degrees of freedom, k i and k f .However, in contrast to the scattering amplitude, M, for W a decomposition in spherical harmonics is not useful.The directional degrees of freedom sweep across pole singularities due to the diagrams in the second line of Fig. 2(c), implying that the decomposition is only defined in the sense of distributions.More importantly, these long-distance singularities guarantee that higher partial waves will not be suppressed.
The issue is easily resolved by removing the singular terms before decomposing in harmonics.This was already discussed in detail in Ref. [6] where the quantity W df was first introduced.In this work we define an alternative, Lorentz-covariant version of the divergence-free amplitude with the same symbol [see again the second line of Fig. 2 where w µ1•••µn is the single-particle matrix element of the external current, defined in detail below, and Here P (cos θ) are the standard Legendre polynomials, satisfying m=− 4πY * m (x)Y m (ŷ) = (2 + 1)P (x • ŷ).Unlike W, W df admits a uniformly convergent decomposition in spherical harmonics where , and the repeated harmonic indices on the right-hand side are summed.The subscript "df", short for divergence-free, refers only to the lack of kinematic singularities arising from a long lived state between the 2 → 2 and 1 + J → 1 transitions.In particular, as we show in Sec.III B 2 below, W df does have so-called triangle singularities as a function of E i and E f associated with diagram of the kind depicted in Fig. 2(b).These pose no problem to the harmonic decomposition but must, of course, be understood in order to successfully extract and interpret both W df and W.
The scattering amplitude, M, is defined with powers of (k /q ), referred to as barrier factors.These must be included due to the manner in which the factors of w µ1•••µn and M arise in the triangle diagram of Fig. 2(b).In particular, the loop is summed over all finite-volume momenta in such a way that the current insertion, w µ1•••µn , as well as the external factors of M, are sampled at off-shell values of momenta (i.e.values for which k 2 = m 2 i ).Nonetheless, the power-law finite-volume effects are governed by the on-shell values of w µ1•••µn and M only, and the off-shell contributions can be absorbed into other infinite-volume quantities.The catch here is that the on-shell projection, effected via k → q , amounts to replacing k Y m ( k ) with q Y m ( k ).The latter has spurious singularities near k = 0 and thus more care is needed.The inclusion of barrier factors resolves the issue.
The single-particle matrix element, w µ1•••µn , appearing in Eq. ( 16), is a function of (P f − k, P i − k) in the first term and (P f − k , P i − k ) in the second.Using the first term as a reference, the explicit definition is given in a three step processes.First, the fully on-shell version is defined via a single-particle matrix element Second, this is formally continued to an off-shell quantity in the context of some generic effective field-theory.The latter object is then decomposed into kinematic tensors and form factors where the sum runs over all possible tensor structures for the given current.Third, and finally, a partial on-shell projection is performed to define the version of w µ1•••µn that appears in Eq. (16).In this step we set the virtualities within the form factors to be on shell ] and also set k 0 = ω k2 everywhere.We reach This definition is not completely on shell because, within K (j) µ1•••µn only, the virtualities (P i − k) 2 and (P f − k) 2 may differ from m 2  1 .In what follows we will consider sums and integrals over the spectator momentum, k.With this in mind, it is convenient to rewrite K (j)  µ1•••µn as a sum of terms that isolate the dependence on this quantity where C (j) µ n +1 ...µn (P f , P i ) has no indices for n = n and Here the superscript ω indicates that the 4-momenta in K ω are on shell.This will be important in Sec.III below, where we introduce various off-shell versions of K in our formulas for evaluating G.
Having defined all of the infinite-volume quantities that enter our formalism, we now turn our attention to the finite volume.As mentioned in the introduction, we restrict attention to a finite cubic spatial volume, with periodicity L applied to the fields in each of the three spatial directions.In this set-up, we consider a matrix element in which the local current J is sandwiched between two finite-volume states, each of which has the quantum numbers of the two-particle system.As we demonstrate in Appendix A, following the derivation of Ref. [6], this LQCD observable is related to the infinite-volume 2 + J → 2 transition amplitude via where W L,df directly determines W df via In the case of distinct particles that both couple to the current one must subtract two terms.In the second the particle labels 1 and 2 are swapped everywhere in the definition of G.In addition, the f (j) (Q 2 ) will take on different values if the particle species differs.These should be given a species label in the case that the current couples to both.See Eq. (A22) as well as Ref. [6] for explicit expressions.
In these expressions we have suppressed angular-momentum indices on R, Each object carries the set m ; m and these are contracted between adjacent factors in the usual matrix multiplication.The trace in Eq. ( 25) is also over this angular-momentum space.Note, by contrast, that the index set µ 1 • • • µ n is not summed but rather fixed to common values for all objects appearing in these equations.
Finally, the matrix ) is defined diagrammatically in Fig. 2(d) and has the explicit form where This differs from the form presented in Ref. [6] due to the aforementioned modifications: The poles are Lorentz invariant and the 1 + J → 1 matrix element is expressed in terms of tensor structures leading to K ω .Note that the modifications to G are directly connected to those in the definition of W df , Eq. ( 16).We have altered these two intermediate objects in such a way that W is unchanged.

C. Examples
In this final subsection we show how the construction outlined above may be applied to specific, phenomenologically well-motivated examples.
We begin with the electro-magnetic form factors of a charged ρ meson.The ρ decays predominantly to the vectorisovector ππ state.Indeed for heavier than physical light-quark masses (such that 4M π > M ρ ) and in the iso-spin symmetric theory, this is the only possible QCD decay channel.If the light-quark mass is further increased, the two-pion threshold eventually exceeds the ρ mass (2M π > M ρ ) and the latter becomes a stable particle.In this case, one can extract the form factors of the ρ directly from finite-volume matrix elements.See, for example, Ref. [69].
To determine the analogous observable at quark masses for which the ρ → ππ decay occurs, it is necessary to first extract the (π + π 0 ) J=1 + j µ → (π + π 0 ) J=1 amplitude for a wide range of kinematic points.As depicted in Fig. 1, by fitting the dependence of the initial-and final-state energies to a functional ansatz, one may analytically continue these amplitudes to the complex-valued pole to obtain the ρ form factors. Detailing the steps of this continuation will be the focus of future work.Here we focus on the extraction of the (π + π 0 ) J=1 + j µ → (π + π 0 ) J=1 amplitude for real ππ energies.
By interpolating isospin-one initial and final states (I = 1, m I = ±1), we project the system to a sector where all even angular momenta vanish.Thus, regardless of the values of P f and P i , we will always have a finite-volume irrep that contains J P = 1 − , with the next contamination coming from J ≥ 3. Taking the latter to be negligible, we can approximate total angular momentum as a good quantum number.However, even in this simple case, the azimuthal component, m J , is not a good quantum number in the finite volume.In general the positive and negative helicity states mix, but one can readily construct linear combinations of these that are invariant under transformations of the cubic group or its little groups [70].
Considering first the incoming state, we restrict attention to a specified set of finite-volume quantum numbers: Λ i and µ i , labeling the irrep and row respectively, of the little group defined by P i .In addition we assume that, within this irrep, the interactions are dominated by i = 1 and a particular m i value so that the matrices can be truncated to single entries.Doing the same for the final-state, Eq. ( 25) becomes8 where δ I=1 ππ is the elastic scattering phase shift for ππ → ππ in the ρ channel and with φ d Λµ (E , L) defined for the irreps that couple to = 1 ππ in, for example, Ref. [60].We stress that all instances of E f and E i in these expressions are to be evaluated at any pair of finite-volume energies, E f,n (L) and E i,n (L) respectively, satisfying Eq. ( 10).In Eq. ( 29), W µ L,df;Λ f µ f ;Λiµi is the subduced version of W µ L,df .As discussed in Ref. [6], this can be obtained from W µ L,df; f m f ; imi by rotating these into the helicity basis and then using the subduction matrices [70].This subduction procedure requires no approximation.However each irrep couples to an infinite tower of partial waves, and only by neglecting these above a certain maximum value does one reach useful expressions.
Compared to Eq. ( 25), in Eq. ( 29) we have dropped the trace, since we are ignoring all but one partial wave, and have solved for W µ L,df .For the latter step there is a potential sign ambiguity that one must address.Note that R I=1 Λµ (P, L) = |R I=1 Λµ (P, L)|e −2iδ I=1 ππ (q ) as is shown, for example, in Eqs. ( 132)-(134) of Ref. [59].The phases in R I=1 Λµ (P f , L) and R I=1 Λµ (P i , L) precisely generate the Watson phases within W µ L,df (P f , P i , L) as they must, since the finite-volume matrix element is real.The remaining ± ambiguity is constrained by the known value in the Q 2 → 0 limit, but in certain cases a remaining ambiguity may survive.
The final step is to convert this to the infinite-volume quantity W µ df via The general form of Eq. ( 26) is overly complicated for this application but still applies with corresponding to the usual (spacelike) pion form factor.Here we have also used the standard relation In these expressions we are neglecting the electro-magnetic form factor of the neutral pion which is expected to be small but non-zero for Q 2 = 0. Finally we remark that a factor of i may appear in this relation depending on the conventions for Euclidean or Minkowski gamma matrices.As the same gamma factors appear in all terms of our formalism changing conventions just amounts to multiplying both sides of µ-indexed equations by a common factor.

(π
The electromagnetic current has both I = 0 and I = 1 components, but G parity guarantees that all matrix elements of the iso-scalar part between two ππ states must vanish.If we take the angular momentum to be unconstrained then the incoming ππ state may, in general carry isospin I i = 0, 1, 2. The iso-vector current then couples the fixed incoming isospin as follows: 0 → 1, 1 → 0, 1, 2 and 2 → 1, 2. If we further restrict attention to |m I | = 1 states (i.e.π + π 0 states) then this reduces to 1 → 1, 2 and 2 → 1, 2. An alternative way to distinguish these possibilities is by fixing orbital angular momentum: J = 0 ⇔ I = 2 and J = 1 ⇔ I = 1.
In this way we identify four possible transitions involving π + π 0 states, the p-wave to p-wave matrix element considered above as well as s to s, s to p and p to s.It turns out that all four of these transitions are described by Eqs. ( 29) to (31) provided that we can neglect the effects of finite-volume mixing with J ≥ 2 states.The only modifications are that one must generally project to different irreps Λ to access the J = 0 components and that the scattering amplitudes must correspond to the isospin of the state with 3. Gluonic structure Thirty years ago, Jaffe and Manohar identified a structure function that provides a measure of the gluon distribution within hadrons [71].This has since lead to lattice calculations of the leading moments of these distributions, for example within the φ-meson.Thus far, the calculations are restricted to heavy quark masses where the φ is stable within QCD [13,72].Similarly, calculations of gluonic moments for light nuclei are already underway, again for values of the light quark masses that lead to the nuclei being deeply bound [12].The formalism presented here will allow for future calculations closer to the physical point by accommodating the finite-volume effects of loosely-bound as well as resonant states.
In order to extract the leading moment of the gluon structure function, one must evaluate the traceless part of the product of two gluon energy-momentum tensors [O µν ∼ G µβ G βν ].As one might expect, this is more complicated than the case considered above in part because it is a rank-two tensor.A starting point in extracting gluonic moments of resonances from LQCD would likely be to consider the ρ, discussed in Sec.II C 1 above.In this case, the relation between the finite-volume matrix elements and the transition amplitude is very similar to Eq. ( 29).The only distinction arises in relating W µν df and W µν L,df .To do so one must determine the scalar (G), vector (G µ ), and tensor (G µν ) contributions to the finite-volume G-function and combine these with the relevant gluonic form-factors of the single-pion state.This concludes our discussion of the covariant formalism.The aim of the section was to provide a procedure by which three inputs: (i) single-particle form factors, (ii) 2 → 2 scattering amplitudes and (iii) finite-volume kinematic functions, can be combined with finite-volume two-particle matrix elements to extract the infinite-volume 2 + J → 2 transition amplitudes.In this recipe the ingredient that remains most obscure is the new finite-volume function G µ1•••µn (P f , P i , L), defined in Eq. ( 27).Thus, in the next section, we give a detailed description of how this can be efficiently evaluated numerically.
Our aim is to evaluate where we have introduced the collective index σ We also take the convention that if σ is written as a low (high) index, then all of the Lorentz indices it contains are also understood to be low (high).The sum is straightforward to calculate numerically once a cutoff function has been included.We comment here that the ultraviolet divergences match between the sum and the integral, meaning that the difference has an unambiguous limit as the cutoff is removed and thus that G σ is a universal quantity with no scheme dependence.
Evaluating the integral part of G σ turns out to be significantly more challenging.The integrand contains singularities associated with on-shell intermediate states and, although these are perfectly integrable given the i pole prescription, numerical evaluations converge very slowly for standard numerical techniques.Thus it is highly advantageous to find analytical representations to the extent possible.
For the case of P i = P f , it turns out that one can provide exact analytical expressions for the integral, as discussed in the following subsection.For the generic case, with P i = P f , we have not managed to obtain fully analytic results.Instead, we write the integral as the sum of two terms.The first includes all singularities and can be evaluated analytically to the level of Feynman parameters.The second is defined with a smooth integrand such that standard numerical integration is effective.Our approach for evaluating G σ in the case of P i = P f case is detailed in Sec.III B below.

A. Pi = P f
The function G σ simplifies considerably when the initial and final momenta coincide, P i = P f ≡ P , i.e. when the momentum transfer vanishes, P µ f − P µ i = 0.A particularly helpful feature of these kinematics is that a natural preferred frame emerges: the simultaneous rest frame of the incoming and outgoing particles, in which the spatial part of P vanishes.
Another consequence of P i = P f is that the product of poles within D(m, k) becomes a double pole of the form Focusing on the factor within parenthesis, note that this can be rewritten as where k and q are defined with respect to the P rest frame.If D(m, k) were to only contain a single pole, then, after acting with the sum-integral difference, only the leadingorder, singular term would be relevant.This is because the sum-integral difference of smooth functions leads only to exponentially suppressed volume dependence that we neglect.However, in this case the first two terms in the expansion are important as they generate double and single poles upon squaring.This leads to where At this stage it is useful to decompose the angular dependence within the tensors into a single set of spherical harmonics As we explain in Appendix B 1, C σ,JM (β β β, k ) can be efficiently calculated by writing the factors within Such factors can then be written as spherical harmonics and, using Clebsch-Gordon coefficients, these can be combined with the external factors of Y f m f (k ) and Y * imi (k ) to identify the a final harmonic basis.As a final note, we stress that it is possible to unambiguously separate the dependence on β β β = P/E and k within C σ , i.e. one can vary k while holding β β β constant.This will be important for the manipulations performed below.
The construction of C σ,JM (β β β, k ) is discussed in detail in Appendix B 1. As a specific simple example, here we consider the case σ = [µ; 10; 10].The numerator within G σ is then The current insertion, k µ , can be written as a combination of = 0 and = 1 spherical harmonics.Combining this with the two = 1 harmonics from the external states, one finds that C σ,JM is zero for J > 3.
The JM = 00 component only has a non-zero contribution in the k 0 component.Isolating this contribution and substituting the definition of the boost matrix, we reach The remaining nonzero coefficients, arising for J ≤ 3, take a more complicated form in general, but simplify for P ∝ ẑ.
If we additionally focus on the µ = z component, then only three coefficients survive We revisit this case below to show how it enters our final construction for G µ 10;10 .Returning to the main line of argument, we now substitute Eq. ( 40) into the expression for G σ (P, P, L), Eq. ( 38), to reach The next step is to expand the k dependence in the first factor within the summand about the pole location In the second term we have rewritten the derivative with respect to k 2 (to be evaluated afterwards at k 2 = q 2 ) directly as a derivative with respect to q 2 .This is possible because C σ,JM (β β β, k ) has no implicit q 2 dependence and it is formally possible to vary q , and thus E , while holding β β β constant.
Combining this with the definition of D(m, k , q ) one can show Remarkably, the operator in parenthesis vanishes when acting on 1/ω q2 so that this reduces to It follows that Eq. ( 44) can be rewritten as To reach our final form we make two additional modifications.First we introduce a cutoff function on the sumintegral difference to enable effective numerical evaluation.Second, we re-express our results in terms of dimensionless quantities r = k L/(2π) and x = q L/(2π).Then, shuffling around terms and introducing a new geometric function, we conclude where These two equations give the main result of this sub-section.In Appendix B 4 we give some details about the evaluation of Z We give two final comments concerning the new kinematic function, Z (n) .First we note an advantage of the decomposition over a single spherical harmonic performed in Eq. (40).It is now straightforward to use the symmetries Im Z (1) 00 JM (top two panels) and Z (2) JM (bottom two).All curves show the function plotted versus E for fixed spatial momentum, P = (2π/L)[001], for which only the M = 0 components are non-zero (up to M = 4).The real parts are shown as solid curves and, for J = 0, the imaginary part is indicated with the dashed curve.As discussed in detail in Appendix B 5, for n = 1 the odd J are exponentially suppressed and indistinguishable from zero on the plotted scale.  of the finite-volume system to identify for which values of JM , Z JM will be nonzero.This is discussed in detail in Appendix B 5.
Second, we remark that the cutoff function used here is designed with the property that the O(α) correction cancels the pole, and thus generates a smooth quantity within the sum-integral difference.If, for example, one were to instead use e −α(r 2 −x 2 ) for all n values, this would still be formally valid, but would lead to corrections of the form α/L k , making it more difficult to estimate the α → 0 limit.In fact these can be systematically subtracted and, as we shown in Appendix B 3, this proves to be an efficient alternative approach for evaluating lim α→0 Z (n) JM .
Suppressing the arguments of Z (n) , the final result for G σ can be written as Note that, in the case of mass-degenerate particles, Z JM = 0 for all odd J.If the particles are at rest in the finite-volume frame then all odd-J functions vanish as does J = 2.This holds for both n = 1 and n = 2 for both degenerate and non-degenerate particles.[See again Appendix B 5.] In Fig. 4 we plot the real and imaginary parts of G µ 10;10 (P, P, L) for P = (2π/L)[001].

B. Pi = P f
We now turn our attention to the more challenging general case of P i = P f .Note that this is realized if any of the four components of the 4-vectors differ, in particular also for P i = P f but E i = E f .As with P i = P f , here the evaluation of the sum is straightforward, while the integral is significantly more challenging.One of the major complications is that the two poles do not coincide in general as one varies the integration variable, k, but may overlap on a two-dimensional subspace for certain choices of external momenta.The contribution of this double-pole submanifold to the integral must be treated with care.
Though we have not found a fully analytic determination of the integral entering G σ , we do have a recipe that gives the desired quantity accurately and with high efficiency.The approach is to rewrite the three-dimensional integral in terms of a d 4 k integral plus a second d 3 k integral with a smooth, singularity-free integrand.The four-dimensional integral, which carries all of the singularity structure, can then be reduced to a one-dimensional integral over a Feynman parameter.The second, smooth term can directly be evaluated using standard numerical integrators.
To give the relevant expressions, we first introduce an extension of the cutoff function entering the definition of where ᾱ = L 4 /(2π) 4 α and r 2 i (2π/L) 2 = k 2 i , x 2 i (2π/L) 2 = q 2 i , etc.We then write where The sum can be evaluated directly as written, and thus we make no further modifications to S σ .The remainder of this section is dedicated to I σ .
1. Separation into IA;σ(P f , Pi) and IN ;σ ( ᾱ, P f , Pi) As summarized above, our approach is to split the integral into a singular part that can be evaluated semianalytically, denoted I A;σ (P f , P i ), and a smooth remainder that is well-suited to numerical evaluation, denoted I N ;σ (ᾱ, P f , P i ).To proceed we define and then introduce D r (m, k) via In Appendix B 6 we give an explicit expression for D r .This term mops up the contributions from the two P -dependent poles in D c .Here the subscripts c and r stand for covariant and remainder, respectively.The idea is to perform the integral of D c semi-analytically and that of D r numerically, and then to take the difference.
In these relations we have neglected the possible factors of k µ and the spherical harmonics multiplying D(m, k).To include these, we first introduce a tensor, M , defined such that where N = n + f + i .A more explicit definition, together with various examples, is given in Appendix B 1.This simply amounts to recasting the factors of k within the harmonics as boosted factors of k µ .
To incorporate K ω ν1•••ν N , note that D r receives two contributions, one each from the poles at k 0 = E f + ω P f k1 − i and k 0 = E i + ω Pik1 − i .(The third pole, at k 0 = ω k2 − i , generates the term we are after, D(m, k).)We thus define where the integral here is a closed clockwise contour encircling the pole indicated.The definitions with an i subscript are given by making the replacement f → i everywhere, and explicit expressions for D ri (m, k) and D rf (m, k) are given in App.B 6.With these quantities in hand, Eq. ( 57) generalizes to Here the integral on the left-hand side is a closed clockwise contour encircling the three poles below the real axis and we have also introduced The next step is to address the issue of ultraviolet convergence for these integrals.Equation ( 56) is manifestly convergent, due to the inclusion of the cutoff function H(ᾱ, k).But to reach an integral that can be evaluated analytically it is convenient to introduce a second form of ultraviolet regularization.We explain the approach first for the special case of two indices, N = 2.Here the integral has a logarithmic divergence that can be removed by taking Eq. ( 62) and subtracting from this the same equation defined with m 1,2 → Λ The regularization scale, Λ, is chosen so that the integrands that depend on it are smooth functions of k with no need for an i prescription.This holds for any Λ satisfying 2Λ > max[E f , E i ] though in practice it is useful to take the cutoff well above this minimum value.On the left-hand side of Eq. ( 64) we have used that, for N = 2, the k 0 integral can be extended to the entire real axis, with a vanishing arc at negative complex infinity.We additionally note that, as a result of the subtraction, the left-hand side and also both square bracketed terms on the right-hand side vanish as 1/|k| 5 in the limit |k| → ∞.These thus give convergent integrals with respect to d 3 k.This approach can be extended to any number of k µ factors, simply by forming more complicated linear combinations to cancel all divergent powers where we have introduced c 0 = 1 and Λ 0 = {m 1 , m 2 }.As above, for j > 0 we take Λ j such that the corresponding integrands are smooth functions of k (Λ j>0 > max[E f , E i ]).In all cases, the linear combinations are constructed such that (i) the k 0 integral extends to the entire real axis with a vanishing contribution from the arc at negative complex infinity, and (ii) the left-hand side and each of the two sums on the right-hand side give convergent integrals with respect to d 3 k.In the following we sometimes refer to this as a Pauli-Villars-like regulator.We give a general algorithm for forming these linear combinations in Appendix B 8.
The final step is to multiply both sides of Eq. ( 65) by the cutoff function, and solve for the desired integral, defined in Eq. (56).We deduce where Equations ( 66)-( 68) are the main results of this subsection.
We emphasize here that the integrands in the definition of I N are smooth for all real k, and the integrals are ultraviolet convergent.For the second and third terms this follows from the fact that H(ᾱ, k) decays exponentially, together with the observation that K r is a smooth function, as is D(Λ j , k) for j > 0. For the first term, smoothness is guaranteed because [H( ᾱ, k)−1] vanishes at the pole and ultraviolet convergence follows from the careful construction of the linear combination.As the integrands are smooth and the integrals are convergent, I N ;σ (ᾱ, P f , P i ) is well suited to numerical evaluation.
It is instructive to consider a few specific examples of this construction, beginning with n = 0, i = f = 0, in which no factors of k µ appear in the numerator of I σ .In this case the integrals are directly convergent, without any need for additional subtraction terms (i.e. the sum over j in Eq. ( 65) reduces to the j = 0 term).Equations ( 67) and ( 68) then reduce to This no-index version of I N is plotted in the left panel of Fig. 5.We emphasize here that I N , and thus also the original integral, diverges in the ᾱ → 0 limit.This is because in the original integral the covariant propagators are evaluated at on-shell k (k 2 = m 2 ), so the propagators scale as 1/|k| and the integrand as d 3 k/|k| 3 .In other words the convergence of I A is always better than that of the original integral by two powers of k, resulting from the off-shell integration of k µ .For fixing the subtraction scheme in Eqs. ( 67) and (68), it is only necessary that I A be rendered finite.[See also Appendix B 8.] We close with one final example: σ = [µ = 0; 10; 10], corresponding to a factor of Y 10 k 0 Y * 10 in the numerator.This leads to an I A;σ integral with an integrand scaling as d 4 kk 2 /k 6 , i.e. diverging as log Λ. Performing a single subtraction of the same integral with m → Λ is therefore sufficient to render the result finite.In fact, to improve the numerical evaluation of I N ;σ , and to test our general method, here we choose to make two subtractions.As explained in Appendix B 8, one possible choice is to add an integral evaluated at Λ = 3m (with coefficient −35/27) and a second at Λ = 6m (with coefficient 8/27).Implementing this in Eqs. ( 67) and (68) leads to convergent forms of I A;σ and I N ;σ respectively, with integrands scaling as dk/k 5 .I N ;σ in this scheme is plotted in the right panel of Fig. 5.
As we include additional factors of k µ the expressions complicate, first because we need additional terms in the sum over j (to maintain convergent integrals) and second because the numerical integrals depend on multiple vector components.However, we find that no conceptual issues arise and the task amounts to coding Eq. ( 68) with an efficient numerical integrator.We give some details on our approach in Appendix B 6, but consider I N ;σ (ᾱ, P f , P i ) as a numerically known function for the remainder of the main text.
Thus it remains only to evaluate I A;ν1•••ν N (P f , P i ), to which we now turn.(q ?i q ?f /m) x c w h / 4 H z + A P 2 I k P w = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " B U g F X x 2 9 I t T 3 C P n e + H g 0 1 m h p x c w h / 4 H z + A P 2 I k P w = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " B U g F X x 2 9 I t T 3 C P n e + H g 0 1 m h p x c w h / 4 H z + A P 2 I k P w = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " B U g F X x 2 9 I t T 3 C P n e + H g 0 1 m h p In the left panel we consider the case of f , m f = i, mi = 00 with no factors of kµ in the numerator.For these kinematics the function requires no subtractions and we directly evaluate IN at ᾱ = 1/3 4 , for fixed external 3-momenta and evenly spaced values of E i , as indicated.In the right panel we plot the function with a numerator factor of 3(k z ) 2 k0.In this case one requires Pauli-Villars-like subtractions, as described in the main text and summarized in the caption.

Evaluating IA;σ(P f , Pi)
We first use the tensor , introduced in Eq. ( 59) above, to define a version of I A with no spherical-harmonic indices The integrals on the right-hand side can then be written as where Here we have used the fact that the sum over j gives a convergent integral, and is thus equal to the δ → 0 limit of the integral in 4 − δ dimensions.Then, at fixed delta, one can exchange the orders of summation and integration, leading to Eq. ( 72).
To evaluate the integral in Eq. ( 73), we first perform the standard Wick rotation on the k 0 integration contour (counter-clockwise to the imaginary axis), and similarly rotate the time components P i,0 and P f,0 .We then make the variable redefinitions to reach where . Note here that the indices are not contracted between ξ ν1•••ν N and the momenta but rather are held fixed on both sides of the equation.
In the second step we have introduced the generating functional As we show in Appendix B 7, this reduces to where In these results we have analytically continued back to real P i,0 and P f,0 and expressed all quantities in terms of the 4-vectors P f amd P i as well as the Lorentz invariants Q 2 = −(P i − P f ) 2 , s i = P 2 i , and s f = P 2 f .The corresponding analytic continuation of Eq. ( 76) is given by Taken together, Eqs. ( 78)- (80) give the main result of this subsection.
As above, it is instructive to consider I A for N = 0, i.e. with no factors of k µ in the numerator.In this case Eq. ( 80) is evaluated with no χ-derivatives and gives As noted at the end of the previous subsection, this integral is convergent in the δ → 0 limit, and simplifies to where we have substituted with The final analytic step is to evaluate the integral with respect to y.We do so via the identity where in the second equality we have introduced a short-hand for the result of the integral.Note that, as long as Imf (x) is nonzero, we can safely set = 0 in these expressions.In the case that the imaginary part does vanish, then we use the relation lim →0 arctan(a/ ) = sign[a]π/2.Applying this to Eq. (82) we deduce where F (1) (x) [also given in Eq. (B98)] is  FIG.6.The real and imaginary parts of IA, generated using the single-parameter integral of Eq. ( 88), with m1 = m2 and all other kinematics as indicated in the legend.The fixed external coordinates are chosen such that a singularity arises at E f,c ≈ 2.07.As explained in the text, the discontinuity is induced when a pole in the Feynman parameter, x, crosses into the integrated region.We vary the value of (in the i pole prescription) to illustrate how the singularity arises as → 0.

In this case Imy
Thus the substitution lim →0 arctan(a/ ) = sign[a]π/2 is only needed when the argument of the square root is positive.As we explain in detail in Appendix B 9, evaluating the remaining integral over x reveals that I A (P f , P i ) has triangle singularities that arise whenever P f and P i take on values for which all three particles in the triangle of Fig. 2(b) can go on shell.
More precisely, we show in Appendix B 9 that the singularity locations are governed by the discriminant of the polynomial A 2 + B ≡ ax 2 + bx + c, given by 2 , and in addition, x c and y c ≡ y + (x c ) fall in the integrated range.It can be shown that these conditions are equivalent to the ones found using Landau's singularity classification [73].
At values of P f and P i satisfying these conditions, the real part of I A has a step-function discontinuity of height and the imaginary part shows a logarithmic divergence.In Fig. 6 we illustrate how these singularities form in the → 0 limit of the i pole prescription.In Fig. 7 we show the singularity structure as a function of E f for various fixed values of E i .In particular one sees that, for sub-thresold E i , I A is a smooth function away from the two-particle production threshold.As E i approaches 2m a step forms in ReI A and ImI A develops a log divergence.Then, as E i is further increased, the location of the singularity in E f moves towards and eventually collides with the two-particle threshold.
To complete this subsection we would like to comment on the behavior of this singularity for some special set of kinematics.First, in the case of identical initial and final 3-momenta, i.e.P i = P f , I A does not have any other singularities apart from those arising at threshold, and therefore the G-function will not exhibit a triangle singularity.This is consistent with our analysis of the P f = P i case, and with the numerical example shown in Fig. 8.In other words, given that all the external momenta in the triangle diagram are time-like, the condition of all three internal propagators to be on shell cannot be realized.A second example is the special case of m 1 = m 2 and s i = s f = s.Solving X(s, s, Q 2 ) = 0 in this simple case leads to a singular manifold given by 1.9  7. The real and imaginary parts of IA, generated using the single-parameter integral of Eq. ( 88), evaluated piecewise in order to work exactly at = 0 as described in the paragraph containing Eq. (87).As indicated, the various curves correspond to fixed values of E i , chosen to illustrate the behavior of the triangle singularity.For E i < 2m, IA is a smooth function of E f away from threshold.As E i approaches threshold from below, ReIA forms a step-function singularity and ImIA a logarithmic divergence.When E i is then further increased, this singularity moves to lower values of E f , eventually colliding with the threshold cusp.
In this case it is easier to visualize that this condition is equivalent to that of all three intermediate particles in the triangle diagram, Fig 2(b), going on-shell.(See also Refs.[74,75].) This concludes our discussion of I A within the main text.In Appendix B 10 we extend the results here by explicitly evaluating I A;ν1 , I A;ν1ν2 and I A;ν1ν3ν3 .For these integrals we find that the above-threshold discontinuities persist, but are milder when factors of k µ appear in the numerator.

Examples of G
Having discussed the integral entering G σ (P f , P i , L) in great detail we are now ready to put everything together and evaluate the complete function.We do so for two different examples of external kinematics.First, in Fig. 8, we consider the case of P i = P f = (2π/L)[001] with σ = [µ = 0; 10; 10].Then, in Fig. 9 we take P i = [000] and P f = (2π/L)[001] with σ = [00; 00], i.e. with no numerator factors.In both cases we fix mL = 6 and plot Re G σ for all values of E i and E f in the region of interest.
As explained in the figure captions, each example illustrates important issues and features that arise.In Fig. 8 we consider various diagonal slices of the E i , E f plane.We find that the result for E i = E f is in perfect agreement with the P i = P f result determined by combining the various Z JM functions.This provides a strong check on the two different methods.Figure 8 also illustrates that double poles arise along the E i = E f line, but these split to single poles as the slice is rotated away from this singular choice.
In addition, Fig. 8 illustrates the results of using the Pauli-Villars-like regulator to separate the integral into I A and I N .As discussed towards the end of Sec.III B 1, the original integral contains a factor scaling as Y 10 k µ Y 10 ∼ |k| 3 in the numerator leading to a log Λ divergence in I A .Following Appendix B 8 we handle this by evaluating Eqs.(67) and (68) with {Λ 1 , Λ 2 } = {3m, 6m} and {c 1 , c 2 } = {−35/27, 8/27}.This removes not only the divergence but also a dk/k 3 term in the integrand to further optimize the numerical convergence of I N .We already gave the result for I N in this prescription in the right panel of Fig. 5.Here the result is combined with the sum and I A to reach G σ .With all building blocks summed together, the dependence on the Pauli-Villars parameters, Λ i , cancels (between I A and I N ) as does the dependence on the smooth cutoff parameter, ᾱ, (between I N and the sum).
Turning now to Fig. 9, this result displays two additional features of G σ .First we see that, when P i and P f differ, R X T M T F p g M m r a k J w F 1 9 e J p 2 z h u s 0 3 L v z e v O q j K O C j t A x O k U u u k B N d I t a q I 0 o y t A z e k V v 1 p P 1 Y r 1 b H / P W F a u c q a E / s D 5 / A G I G l T w = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 R X T M T F p g M m r a k J w F 1 9 e J p 2 z h u s 0 3 L v z e v O q j K O C j t A x O k U u u k B N d I t a q I 0 o y t A z e k V v 1 p P 1 Y r 1 b H / P W F a u c q a E / s D 5 / A G I G l T w = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 R X T M T F p g M m r a k J w F 1 9 e J p 2 z h u s 0 3 L v z e v O q j K O C j t A x O k U u u k B N d I t a q I 0 o y t A z e k V v 1 p P 1 Y r 1 b H / P W F a u c q a E / s D 5 / A G I G l T w = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 A T P 4 M V 6 s J 6 s V + t t 1 p q z 5 j P 7 4 E d Z 7 1 / + l 5 A + < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " A T P 4 M V 6 s J 6 s V + t t 1 p q z 5 j P 7 4 E d Z 7 1 / + l 5 A + < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " A T P 4 M V 6 s J 6 s V + t t 1 p q z 5 j P 7 4 E d Z 7 1 / + l 5 A + < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " A T P 4 M V 6 s J 6 s V + t t 1 p q z 5 j P 7 4 E d Z 7 1 / + l 5 A + < / l a t e x i t > [+0.10, +0.15] < l a t e x i t s h a 1 _ b a s e 6 4 = " j H s U L T 4 5 + Y X V 8 X t k y l e i p M o y T c 8 < l a t e x i t s h a 1 _ b a s e 6 4 = " 6 L      < l a t e x i t s h a 1 _ b a s e 6 4 = " u E v J J 5 K 7 d F R S u 7 r l 7 S t y Z B P + A 8 e w T N 4 s R 6 s J + v V e p u O 5 q z Z z i 7 4 A e v 9 C 1 4 J l G M = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " u E v J J 5 K 7 d F R S u 7 r l 7 S t y Z B P + A 8 e w T N 4 s R 6 s J + v V e p u O 5 q z Z z i 7 4 A e v 9 C 1 4 J l G M = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " u E v J J 5 K 7 d F R S u 7 r l 7 S t y Z B P + A 8 e w T N 4 s R 6 s J + v V e p u O 5 q z Z z i 7 4 A e v 9 C 1 4 J l G M = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " u E v J J 5 K 7 d F R S u 7 r l 7 S t y Z B P + < l a t e x i t s h a 1 _ b a s e 6 4 = " 6 6 Z X k 1 3 P w 0 s 8 F z 0 < l a t e x i t s h a 1 _ b a s e 6 4 = " 7 t 0 B y X g 6 s 9 a S U c z s o z 8 w P r 4 B 6 o + a F A = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " 7 t 0 B y X g 6 s 9 a S U c z s o z 8 w P r 4 B 6 o + a F A = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " 7 t 0 B y X g 6 s 9 a S U c z s o z 8 w P r 4 B 6 o + a F A = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " 7 t 0 B y X g 6 8. Contour plot representation of Re G µ=0 10;10 (P f , Pi, L) for P f = Pi = [001] and mL = 6.The grey diagonal dashed lines in the left panel indicate slices defined by E f = s(E i − 2m) + 2m, for s = 1.0, 0.8, 0.6.We plot Re G σ along these slices in the right panel as indicated.The top right panel corresponds to P f = P f and the plot matches Fig. 4 well within the expected e −mL discrepancies.Here we also separately show the contributions from IA (dashed orange) and the remaining contribution to Re G σ (dashed red).
the non-interacting two-particle poles appear in different locations for E i and E f .Interactions shift the finite-volume energies away from these singularities so that G σ , like the Lellouch-Lüscher factors, will generally be evaluated away from the divergent locations.However, as with all finite-volume kinematic functions, this implicit knowledge of the non-interacting spectrum is a key ingredient in the all-orders correction of the scattering-state volume effects.
Second, we see the appearance of triangle singularities inherited through I A .Such features are simply part of the correct definition of G σ .Indeed, because the singularity structure is directly induced by the infinite-volume diagram of Fig. 2(b), it also appears within the infinite-volume 2 + J → 2 transition amplitude itself.The steps, cusps and log divergences of G σ are present in both W and W df for exactly the same kinematics.Thus, understanding the features is crucial to extracting and interpreting the infinite-volume observables that we are after.We discuss the role of discontinuities within the transition amplitude in more detail in future work.

IV. CONCLUSION
In this work we have presented a modified version of the finite-volume formalism for studying 2 + J → 2 transition amplitudes.This is closely related to the approach of Ref. [6], but differs in that all infinite-volume quantities are Lorentz covariant and the 1 + J → 1 matrix elements have been reformulated in terms of standard form factors.As explained in Sec.II B and Appendix A, the new result is reached by making minor adjustments to the derivation presented in Ref. [6].For example, in that work finite-volume effects are expressed as sums over poles of the form 1/[2ω k (k 0 − ω k )], and here the same effects are expressed via invariant poles, 1/(k 2 − m 2 ).
These changes lead to a modified form of the finite-volume function, denoted G, with an added benefit that the new form is easier to evaluate numerically.As described in Sec.III B 2, the Lorentz covariant structure allow us to write the integral appearing in G in terms of Feynman parameters.This reduction is also crucial to revealing the analytic structure of G, including the triangle singularities described in great detail Sec.III B 2 and Appendix B 9, and illustrated in Figs. 6, 7 and 9.We also recall that the form of G presented in Ref. [6] carries four sets of spherical harmonic indices, resulting from a cumbersome description of the 1 + J → 1 matrix elements.By contrast, our improved expression carries only the angular momentum indices of the external states, together with Lorentz indices to describe the current insertion.
To avoid proliferation of flavor and channel indices, in this work we restricted attention to kinematics for which a single channel of two scalar particles is open.Accommodating multiple channels is straightforward given the results 9. Contour plot representation of Re G00;00(P f , Pi, L) for Pi = [000], P f = [001] and mL = 6.Allowing the spatial momenta to differ means that the finite-volume spectrum must be different for the incoming and outgoing states, as is apparent from the different positions of the poles, corresponding to non-interacting levels (vertical and horizontal dashed lines in the left panel).Another feature of differing spatial momenta is the appearance of triangular singularities, emphasized in the middle and right panels.On the far right we plot three slices, running over the step discontinuity in Re G σ .
of Ref. [6].Incorporating particles with spin has yet to be considered for these types of observables and is the subject of future work.This final generalization is of great importance given the phenomenological interest in two-nucleon matrix elements [76,77].From our previous experience with spinning particles [6,32], we expect that the extension will be relatively straightforward.
Looking to less trivial extensions, it would be of great interest to extended these ideas in order to develop an approach for extracting non-local matrix elements of two-particle systems.This would make it possible to extract distribution functions of resonant states, following the methods of Refs.[10], and would open the possibility for lattice QCD calculations of two-body contributions to double-beta decays [78,79].Matrix elements of non-local operators suffer from different types of finite-volume artifacts.These depend crucially on whether the operators are displaced in Euclidean time, as in Ref. [80], or in a spatial direction, as considered for example in Ref. [11].Finally, extensions of this work to energies for which three or more particles can go on-shell should be feasible in the future, especially given the recent progress in understanding the finite-volume spectrum of three-particle states [47][48][49][50][51][52].
Returning to the present formalism, several open questions persist that we plan to address in future work.For example, W df is defined using the partial-wave basis, while the form factors of resonances or bound states are more naturally described using Lorentz decomposition.It is always possible to relate the partial-wave and Lorentz bases.However, due to the reduction of rotational symmetry in the finite-volume, we do not expect a one-to-one correspondence between the finite-volume matrix elements and the different Lorentz components of W df .For example, in the ρ channel with nonzero spatial momenta in the finite-volume frame, the different helicity components mix to different finite-volume irreps.This means that the components are sampled by different finite-volume quantization conditions and thus at different energies.As a result, just as is done in the analysis of coupled-channel scattering [40][41][42][43][44][45], it will be necessary to perform global fits of the matrix elements using Eq. ( 25).This requires a detailed understanding of the analytic structure of these amplitudes in which triangle singularities play a crucial role.
P A E L / B q P B r P x p v x P h s t G P O d f f g F 4 + M b 1 + e U n w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " i M + h p o y w K e D l P S P A E L / B q P B r P x p v x P h s t G P O d f f g F 4 + M b 1 + e U n w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " i M + h p o y w K e D l P S P A E L / B q P B r P x p v x P h s t G P O d f f g F 4 + M b 1 + e U n w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " i M + h p o y w K e D l P S In this appendix we describe how the derivation of Ref. [6] must be modified to give the covariant version of the formalism, presented in Sec.II B above.As in our previous derivation, we restrict attention to a finite cubic spatial volume, implemented by requiring all fields to have periodicity L in the three spatial directions.
Within this set-up we introduce a finite-volume three-point function, C 2→2 L (P f , P i ), defined as the sum over all possible diagrams connecting the initial and final states to the inserted current.See also Fig. 10(a).For E i , E f below the next multi-hadron threshold (such that only a single two-particle channel can propagate), all volume effects scaling as a power of 1/L are captured by the skeleton expansion shown in the figure .Here the label V within the loops stands for volume and denotes that the diagrams are defined with the spatial momenta summed over the discrete set allowed by the periodic boundary conditions, k = 2πn/L with n a 3-vector of integers.The corresponding diagrams in an infinite volume are given by replacing these sums with integrals and are represented by a loop with no label.
The power-like volume effects of C 2→2 L (P f , P i ) are encoded in the skeleton expansion of Fig. 10(a), built from fully dressed hadron propagators (indicated by the simple black lines) and two-particle Bethe-Salpeter kernels (indicated by open circles).The vertices with a current insertion are given by the same diagrammatic set defining propagators and kernels, but with the current attached at all possible locations.In the kinematic window of interest, the difference between the finite-and infinite-volume definitions of propagators and kernels as e −mL .
The set-up here is identical to that of Ref. [6].Indeed the only modifications required are in the evaluation of the two-particle loop, in which the current couples to one of the two particles.The relevant diagram is shown in Fig. 10(b).To explain how the analysis is altered, we begin by recalling the finite-volume reside of this diagram, given in Eq. ( 27) of Ref. [6] G Here L(P f , k) and R † (P f , k) are generic endcap functions to be replaced with Bethe-Salpeter kernels or the overlap to the interpolators in the final derivation.In contrast to Ref. [6], we define G L,µ1•••µn with these endcaps accompanied by factors of i.The difference arises because we formulate the derivation here with Minkowski momenta (in contrast to the Euclidean conventions of the previous publication).In this set-up the i factors give a more natural extension to the Bethe-Salpeter kernels, multiplied by this factor due to the weight, e iS , in the quantum path integral.Our second notational modification is to represent the 1 + J → 1 insertion with a set of Lorentz indices, and to make explicit that the quantity is off-shell, i.e. that (P [See also Eq. ( 21) above.]Our third and final alteration is to restrict attention to a single channel, thus removing the a and b indices from Eq. ( 27) or Ref. [6].
We now express L(P f , k), R † (P f , k) and w off µ1•••µn (P f − k, P i − k) in terms of their on-shell counterparts, plus corrections.For the endcap functions this is done exactly as in Ref. [6], by first defining and then recombining the components with the Y m harmonics defined in Eq. ( 13) It now remains only to separate the propagators by defining where is the non-interacting covariant propagator.Our conventions are such that ∆(k 2 ) has unit residue at the single particle pole, implying S(k 2 ) is smooth and finite near k 2 = m 2 1 .This form of separation, in which D(k 2 ) remains Lorentz invariant, is the second key distinction relative to our earlier formalism.
At this stage we have separated the endcaps [Eqs.(A4), (A5)] the single-particle form factor [Eq. (A11)] and the covariant propagator [Eq.(A15)] into on-shell terms plus corrections.Substituting these four identities into Eq.(A1) then gives the analog of Eq. ( 43) of Ref. [6].The final step is to rearrange terms according to singularity structure.Of the 64 terms (reach by multiplying 4 binomials as well as the form-term separation of f (j) ) all but 17 are smooth at both poles, and thus give only exponentially suppressed contributions to the sum-integral difference.Those that are singular break into three classes, 8 with only the 1/(k 2 f + m 2 − i ) singularity, 8 with 1/(k 2 i + m 2 − i ), and a single maximally singular term that leads to the appearance of G σ (P f , P i , L).
To give an explicit form for the single-pole terms we need to introduce one final piece of notation.We define the operator δ df via and similarly for δ df acting on the left.We can make use of this to rewrite Eq. (A1) as where D(m, k) is defined in Eq. ( 28).On the second line the square bracketed quantities have been decomposed in spherical harmonics and carry implicit indices that are contracted with F (P, L), defined in Eq. ( 12).We have included the inverted symmetry factor, ξ −1 = 2 for identical particles, to compensate the factor in the definition of F .Note that this arises naturally from the fact that the current couples to each of the two particles when these are identical.This implies that W df is defined with four subtraction terms, given by coupling the current to each of the four external propagators.When W df is then projected to definite angular momenta, these terms become pairwise redundant leading to the factors of 2. In the case of non-identical particles that both couple one must sum over the two choices of species mass within ∆(k 2 f ) as well as the alternative mass assignments within D(m, k), i.e. {m 1 , m 2 } → {m 2 , m 1 }.
Using Eq. ( 23), we can write the first term on the right-hand side in terms of G µ1•••µn; f m f ; imi (P f , P i , L), defined in Eq. (27).Following the steps outlined in Ref. [6], one finally arrives at Eq. ( 25) -the relationship between the finitevolume matrix element and iW df;µ1•••µn , defined in Eq. ( 16).The latter emerges through the identity δ df Wδ df ≡ W df .
We close by giving explicit expressions for the case that the current has a non-negligible coupling to both of the two particles.When two distinct species, 1 = 2, each admit a 1 + J → 1 transition, then the definition of W df , Eq. ( 16), is replaced by an expression with four subtractions As a specific example, we return to the case of σ = [µ = z; 10; 10] and P = [00d z ], already discussed in the main text.For these indices, Eq. (B3) reduces to Here we have used the fact that, with µ fixed to z, we only need to sum over one row of the boost matrix, [Λ −β β β ] z ν = (β z γ, 0, 0, γ).The integrals are now trivial to evaluate.For example, the JM = 00 component reduces to where we have used β z = P z /E and γ = E/E .This matches Eq. ( 42) of the main text.An alternative approach for determining C σ,JM (β β β, k ) is to note that the product of two spherical harmonics can be written in terms of Clebsch-Gordon coefficients and that each factor of k ν can be rewritten using the identity where the last equality defines T ν m (k ).Substituting this into the left-hand side of Eq. (B2) leads to sums over products of spherical harmonics on that side of the equation.These can be pairwise combined using Clebsch-Gordon coefficients, until the left-hand side is reduced to a sum over a single harmonic.Then one can use the orthogonality of spherical harmonics to match this, term by term, to the right-hand side and thereby determine the values of C σ,JM (β β β, k ).
We next consider the case of P i = P f .As explained in the main text, here we find it more useful to convert indices in the other direction, i.e. to trade all dependence on spherical harmonics for additional momentum 4-vectors.The key distinction between this case and that discussed above is that we no longer have a natural c.m. frame.The rest frames of P i and P f differ, and expressing the integrand in either frame leads to ugly expressions.This, together with the need to reach covariant expressions that we can evaluate semi-analytically, led us to introduce M ν1•••ν N σ , in Eq. ( 71) above.The definition can be re-expressed as The second, and final aim of this subsection, is to derive a useful expression for M ν1•••ν N σ .First we introduce a new set of tensors, denoted T , that allow us to express the right-hand side in terms of 4-vectors The exact definition of the can be inferred by comparing Eqs.(B8) and (B9).Note that they contain the √ 4π/(q ) prefactors and also encode the combinations of k components needed to form the various spherical harmonics.For example , is closely related to T ν m , introduced above.The final step is to boost all 4-vectors to the finite-volume frame.We deduce Here we give some further details on the evaluation of I N ;σ (α, P f , P i ), defined in Eq. ( 68) of the main text.Recall that this is a 3-dimensional integral with a smooth integrand, to be evaluated numerically.When added to the semi-analytic expression, I A;σ , it gives the full integral entering the definition of G σ .
The complication we wish to address here is that the definition in the main text requires evaluating a large number of integrals defined with factors of k µ (but with no spherical harmonics) in the integrand.These are then contracted with the tensor M ν1•••ν N σ to reach the final expression.While highly advantageous for the analytic integral, the Lorentz-index-based expressions lead to a costly determination of I N ;σ .For example, for f = i = 1 and a single current index, the covariant form contains factors of k ν1 k ν2 k ν3 that naively require the evaluation of 64 terms.
To improve this situation, it is preferable to move M ν1•••ν N σ back inside the integrals defining I N ;σ , and thereby rewrite the integrands to only carry the final σ index.This procedure is a bit subtle, because the time components of the k µ are evaluated at various values.To proceed we first recall that I N ,σ generally involves terms evaluated at the physical masses, m 1 , m 2 , together with regulating integrals evaluated at some higher scale Λ j .For the physical-mass terms, time components are evaluated at k 0 ∈ ω k2 , E f + ω P f k1 , E i + ω Pik1 , and for the regulating integrals at . Some of these four vectors are recombined into the harmonics, Y m (k ), and, because the time and space components mix upon boosting k µ = [Λ β β β ] µ ν k ν , we end up with strange spatial components in some of the harmonics.
To give concrete expressions it is convenient to define Here Ω k represents any of the possible choices made for the temporal component of the 4-vectors.In each of these three cases, an implicit mass dependence enters and the second argument, Λ, refers to this mass dependence.In the following we will use m or Λ 0 to indicate that the ωs are evaluated at their physical masses and Λ j>0 to indicate evaluation at an unphysical value m 1 = m 2 = Λ j>0 .In short, the first two entries in N σ simply serve to indicate the value at which all k 0 are evaluated.We stress that the tensor does not depend on these parameters but only on the c.m. frame energies E i and E f as well as the boost velocities β β β i and β β β f .Thus, the only modification to the spherical harmonics is that they now depend on k ,Ω , defined via If we set Ω k = ω k2 and Λ = m, then we exactly recover the spherical harmonic definitions used everywhere else in this work.
With our new numerator function in hand, we are ready to give our final form for I N ;σ where ) This is identically equal to the quantity defined in Eq. ( 68) of the main text.The only difference is that we have absorbed M ν1•••ν N σ inside the integrands, via the new function N σ .To complete the specification we require explicit expressions for D rf (m, k) and D ri (m, k) The motivation for these quantities is discussed in the text around Eq. ( 60), where an implicit definition is also given.We close this subsection with one final simplification to I N .We show that it is always possible to simplify the numerical evaluation from a three-down to a two-dimensional integral.To show this we must prove the following key identity Begin the proof by considering the generic integral assumed to be convergent.Here we have separated the integrand into a Lorentz-scalar function, F, together with an arbitrary function of the c.m. frame momenta, G.The latter is also Lorentz invariant in the vacuous sense, i.e. because its arguments carry a frame label.Now act on both sides with Λ β β β , on each index To simplify the right-hand side we perform a change of integration variable k µ ≡ [Λ β β β ] µ ν k ν and also define P µ f ≡ [Λ β β β ] µ ν P ν f and Here we have used the Lorentz-invariance of the various building blocks, including the fact that k i and k i must be unchanged if we replace k, P f and P i with their boosted counterparts.This result, which can be rewritten as is just a statement of Lorentz-covariance for X .To conclude our demonstration of Eq. (B40) we note that I N ≡ I − I A and that the two terms on the right-hand side each satisfy the functional form of X , shown in Eq. (B41).In the case of I, the original integral defining G, one takes F(P f , P i , k) = Θ(k 0 )(2π)δ(k 2 − m 2 2 ) This function is only invariant under orthochronous transformations, as is standard when one discards the anti-particle pole, but this is sufficient for the present argument.For both I and I A the spherical harmonics, as well as the cutoff function H(α, k) can be absorbed into the definition of G. Again the key point is that these objects are frameindependent because they carry a frame label, k (k, P ) = k (k , P ).We deduce that I N must satisfy Eq. (B44).
To see the power of this identity we take β β β = β β β i , implying and thus that only one external direction enters the integral.In this case the integration coordinate is simply transformed to k i .We are then left with where the sum over k runs over all possible choices for the temporal component, as discussed above.The key point here is that, once the spherical harmonics and the factors of k µ i are factored out, the remaining integrand can only depend on k i through its magnitude and a single angle.This is because no other direction is defined in the system once we have expressed all coordinates in the rest frame of the incoming state.Defining cos θ ≡ k i • P f and picking an aribtrary orientation for the azimuthal angle φ, we reach Here the φ integral can be evaluated analytically using the rotation properties of the spherical harmonics.Thus only the k i and θ integrals need to be performed numerically.
7. Evaluating I χ (P f , Pi, m, δ) In this subsection we show how the the generating functional, I χ (P f , P i , m, δ), defined in Eq. ( 77) of the main text, reduces to the results given in Eqs. ( 78)- (79).We begin by inserting Feynman-parameter integrals into the definition to reach where in the second equality we have rewritten the integral with a dummy parameter α, to be set to zero after differentiation.Next we multiply both sides by exp[−k 2 E ] and integrate with respect to dk E k 4−δ E to write Solving for A n and evaluating the remaining integral using Schwinger parameters, we deduce A n (χ 2 ) n (2n)!J n (P f , P i , m, δ) , (B60) where we have returned to the Minkowski signature for χ, P f and P i and have also defined (B61) To conclude we simplify J n (P f , P i , m, δ) by evaluating the momentum integral This directly gives Eq. ( 78) in the main text.

Determining cj and Λj
As discussed in Secs.II C 1 and II C 2, an immediate application of the proposed formalism is electromagnetic reactions coupling two-pion states: (π + π 0 ) Ii + j µ → (π + π 0 ) I f .For this case we require the function G σ for σ = [ f m f ; i m i ] (no current index) as well as σ = [µ; f m f ; i m i ] (one index) for i , f ≤ 1.This requires evaluating I A (with no indices) through I A;ν1ν2ν3 .This set depends on only two scalar integrals, J 0 and J 1 .The integral defining J 0 is convergent so that we only need the c 0 = 1, Λ 0 = m term, i.e. no subtraction is required.The integral defining J 1 , by contrast, has a logarithmic divergence (arising from d 4 k E k 2 E /k 6 E ).This is removed using the subtraction given in Eq. ( 64), corresponding to c 1 = −1 with Λ 1 = Λ equal to any value exceeding 2m.
Though the integrals of direct interest are rendered convergent by (up to) one simple subtraction, we think it useful here to give the recipe for general I A,ν1•••ν N .In a nut shell the approach requires identifying the divergent part of J n and, by substituting this into the relation to I A,ν1•••ν N , to identify an expression of the form I div A;ν1•••ν N (P f , P i , δ) ≡ N/2 j=0 c j I div A;ν1•••ν N (P f , P i , Λ j , δ) . (B66) The coefficients c j and the scales Λ j are then tuned such that this quantity vanishes.Note that N/2 terms must be tuned to vanish.This is because the integral I A;ν1•••ν N depends on the scalar integrals up to n = N/2 .The integral then scales as d 4 k E k 2n E /k 6 E and generates a series of divergences of the form Λ 2n−2 , Λ 2n−4 , • • • , log(Λ), so a total of n = N/2 terms.
It turns out that one does not need to tune both the regularization scales Λ j and the coefficients c j .We thus choose the recipe of setting Λ j = 2 j−1 Λ [Λ 1 = Λ , Λ 2 = 2Λ , Λ 3 = 4Λ , • • • ] and tuning only the n different c j terms.
The latter is achieved by studying the 1/δ terms.In particular one can show that J n ∼ M (m, x, y) 2n−2 /δ.Thus, for n = 1 the divergence is M (m, x, y)-independent and is removed by setting c 1 = −1 as explained above.For n = 2 the divergence scales as M (m, x, y) 2 leading to the linear combination M (m, x, y) 2 + c 1 M (Λ, x, y) 2 − (1 + c 1 )M (2Λ, x, y) 2 where we have already enforced 1 + c 1 + c 2 = 0.This leads to a cancellation in all terms except for the explicit mass and Λ dependence M (m, x, y) 2 + c 1 M (Λ, x, y) 2 − (1 + c 1 )M (2Λ, x, y) 2 = m 2 + c 1 Λ 2 − 4(1 + c 1 )Λ 2 , (B67) and requiring this to vanish gives Here we restrict our attention to the degenerate case, m 1 = m 2 = m.While the subtraction is always mass-independent for integrals up to I A;ν1ν2ν3 , beyond this the scheme in the non-degenerate case becomes more complicated.Note that the choices of coefficients that remove J 2 divergences also automatically remove those in J 1 .This is because j c j M (Λ j , x, y) 2 contains both Λ j -dependent and independent pieces.Thus the vanishing of j c j M (Λ j , x, y) 2 guarantees that the same holds for j c j , and the latter is the condition for removing divergences in J 0 .
As we now show, this pattern continues to all orders, so that it is always sufficient to determine c j by tuning away the highest J n divergences.For the general-n, degenerate case, the system that we need to solve is It is straightforward to invert this matrix and read-off the values of c i to regulate an integral with any number of indices.In Table II we give the values up to n = 5, assuming Λ j = 2 j−1 (3m).
Returning to the case of different masses, here one must instead solve The key distinction is that x and y dependence now appears in both M (0, x, y) 2 and the mass terms on the right-hand side.Thus it is not possible to express the right-hand side in terms of a matrix product in which M (0, x, y) 2k is factored off.The upshot is that, in this case, the c j depend on the kinematics s i , s f , q 2 , the masses m 1 , m 2 and the generating parameter χ.Again we stress that this is not of immediate concern as it is only relevant for form factors with many indices.
9. Triangle singularities within IA(P f , Pi) In this appendix we give a more detailed discussion of the singularities that arise away from threshold in I A (P f , P i ), as summarized around Eqs. ( 90) and (91) of the main text and in Figs.(B73) These arise when the three particles in the triangle loop of Fig. 2(b) can all go on shell.As Landau described in Ref. [73], the on-shell condition is realized at critical points (x c , y c ) of M (m, x, y) 2 defined by three conditions: ∂M (m, x c , y c ) 2 /∂x = 0, ∂M (m, x c , y c ) 2 /∂y = 0, and M (m, x c , y c ) 2 = 0. Since M (m, x, y) 2 is, at most, a second degree polynomial in x and y, solutions to the three conditions can be found analytically.
To see the role of these critical points in practice, we integrate Eq. (B73) with respect to y, to reach where with y ± (x) ≡ (1/2)(A ± √ A 2 + B + i ).Here A and B are known functions of the kinematic variables and the Feynman parameter x, defined in Eqs.(85) and (86).
We next note that the three Landau conditions given above are satisfied whenever A 2 + B = 0 and, in addition d[A 2 + B]/dx = 0. Noting that A 2 + B is a quadratic polynomial in x we see that the conditions are equivalent to A 2 + B ∝ (x − x c ) 2 .Before considering this special case we take our general form and substitute A 2 + B = 2[(P i • P f ) 2 /s 2 i − s f /s i ](x − x 1 )(x − x 2 ) For concreteness, here we have assumed (P i • P f ) 2 /s 2 i − s f /s i > 0. In addition we have split the integral in x into a sum over regions where Im √ A 2 + B = 0 (x ∈ [0, x 1 ] ∪ [x 2 , 1]) as well as the region where it is non-zero (x ∈ [x 1 , x 2 ]).This separation assumes 0 < x 1 < x 2 < 1.If these are instead complex valued, or outside the range of integration, then one can directly evaluate the integral over the entire range.For general values of x 1 and x 2 , the integral over each region is well defined even with = 0 and can be directly evaluated.
The final case to consider is when the external kinematics are tuned to critical values P f = P f,c and P i = P i,c for which x 1 = x 2 = x c ∈ [0, 1].This is equivalent to the Landau conditions mentioned above, and corresponds to the apex of the (x − x 1 )(x − x 2 ) parabola sitting on the x-axis.Perturbing the kinematics in one direction shifts the parabola down, opening a finite region of [x 1 , x 2 ] that must be integrated in isolation.Perturbing the kinematics in the other direction shifts the parabola upwards, causing the roots to become complex such that we can directly integrate x from 0 to 1.We now demonstrate that, as one approaches P f,c and P i,c from the side of real x 1 , x 2 , the integral ´x2 x1 dx has a nonzero limit due to a singularity in the integration range.As a consequence, the limit has a different value when approached from opposite sides.This manifests as a step singularity in the real part of I A .In addition, the imaginary part diverges as log|P f (i) − P f (i),c |.
The magnitude of the discontinuity is given by evaluating the integral between x 1 and x 2 , for kinematics such that 0 < x 1 < x 2 < 1, and then taking the limit x 1 , x 2 → x c .This can easily be done by noting that, in this region, y + (x) and y − (x) are complex conjugates of each other so that the integrand simplifies to Disc(I A ) = ˆx2 x1 dxF (1)  We see that only the first term will contribute in the limit x 1 → x 2 , and that it only contributes when y c = y + (x c ) is in the y integration region, i.e. 0 < y c < 1 − x c .Evaluating the remaining integral, we conclude Here we provide compact expressions here for the I A;σ (P f , P i ) integrals with up to three Lorentz indices.Starting with Eqs.(78) and (80) one can show I A;ν1 (P f , P i , m) = P f,ν1 I (1,1) (P f , P i , m) + P i,ν1 I (1,2) (P f , P i , m) , (B82) I A;ν1ν2 (P f , P i , m) = P f,ν1 P f,ν2 I (2,1) (P f , P i , m) + P i,ν1 P i,ν2 I (2,2) (P f , P i , m) + P [i,ν1 P f,ν2] I (2,3) (P f , P i , m) − g ν1ν2 4 I (2,4) (P f , P i , m) , (B83) I A;ν1ν2ν3 (P f , P i , m) = P f,ν1 P f,ν2 P f,ν3 I (3,1) (P f , P i , m) + P i,ν1 P i,ν2 P i,ν3 I (3,2) (P f , P i , m) where the brackets in the indices denote a sum over permutations, even when the indices are identical.The definition is such that, for n indices within a pair of square brackets, the sum runs over n! terms (some of which may vanish).Some examples include P [f,0 P i,0] = 2P f,0 P i,0 , P [f,0 P f,0 P i,0] = 6P 2 f,0 P i,0 , P [f,1 P i,2] = P f,1 P i,2 + P f,2 P i,1 , g [00 P i,1] = 2g 00 P i,1 . (B85)

2 + J ! 2 <and form factors 1 + J ! 1 <
l a t e x i t s h a 1 _ b a s e 6 4 = " b o 0 j 1 L j g t q + 2 H x e R n h h o l o b 7 S 8 E = " > A A A C D X i c b V D L S s N A F J 3 4 r P U V d e l m s A q C U J I i 6 L L o R l x V s A 9 o Q p l M J + 3 Q y S T M 3 I i l 9 A f c + C t u X C j i 1 r 0 7 / 8 Z J G 0 R b D w w c z r m X u e c E i e A a H O f L W l h c W l 5 Z L a w V 1 z c 2 t 7 b t n d 2 G j l N F W Z 3 G I l a t g G g m u G R1 4 C B Y K 1 G M R I F g z W B w m f n N O 6 Y 0 j + U t D B P m R 6 Q n e c g p A S N 1 7 E M P 2 D 0 E I a 7 g E + x F B P q U C H y N P Y j x j 9 W x S 0 7 Z m Q D P E z c n J Z S j 1 r E / v W 5 M 0 4 h J o I J o 3 X a d B P w R U c C p Y O O i l 2 q W E D o g P d Y 2 V J K I a X 8 0 S T P G R 0 b p 4 j B W 5 k n A E / X 3 x o h E W g + j w E x m 9 + p Z L x P / 8 9 o p h O f + i M s k B S b p 9 K M w F d h E z a r B X a 4 Y B T E 0 h F D F z a 2 Y 9 o k i F E y B R V O C O x t 5 n j Q q Z d f w m 9 N S 9 S K v o 4 D 2 0 Q E 6 R i 4 6 Q 1 V 0 h W q o j i h 6 Q E / o Bb 1 a j 9 a z 9 W a 9 T 0 c X r H x n D / 2 B 9 f E N v 5 2 a E w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " b o 0 j 1 L j g t q+ 2 H x e R n h h o l o b 7 S 8 E = " > A A A C D X i c b V D L S s N A F J 3 4 r P U V d e l m s A q C U J I i 6 L L o R l x V s A 9 o Q p l M J + 3 Q y S T M 3 I i l 9 A f c + C t u X C j i 1 r 0 7 / 8 Z J G 0 R b D w w c z r m X u e c E i e A a H O f L W l h c W l 5 Z L a w V 1 z c 2 t 7 b t n d 2 G j l N F W Z 3 G I l a t g G g m u G R 1 4 C B Y K 1 G M R I F g z W B w m f n N O 6 Y 0 j + U t D B P m R 6 Q n e c g p A S N 1 7 E M P 2 D 0 E I a 7 g E + x F B P q U C H y N P Y j x j 9 W x S 0 7 Z m Q D P E z c n J Z S j 1 r E / v W 5 M 0 4 h J o I J o 3 X a d B P w R U c C p Y O O i l 2 q W E D o g P d Y 2 V J K I a X 8 0 S T P G R 0 b p 4 j B W 5 k n A E / X 3 x o h E W g + j w E x m 9 + p Z L x P / 8 9 o p h O f + i M s k B S b p 9 K M w F d h E z a r B X a 4 Y B T E 0 h F D F z a 2 Y 9 o k i F E y B R V O C O x t 5 n j Q q Z d f w m 9 N S 9 S K v o 4 D 2 0 Q E 6 R i 4 6 Q 1 V 0 h W q o j i h 6 Q E / o B b 1 a j 9 a z 9 W a 9 T 0 c X r H x n D / 2 B 9 f E N v 5 2 a E w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " b o 0 j 1 L j g t q + 2 H x e R n h h o l o b 7 S 8 E = " > A A A C D X i c b V D L S s N A F J 3 4 r P U V d e l m s A q C U J I i 6 L L o R l x V s A 9 o Q p l M J + 3 Q y S T M 3 I i l 9 A f c + C t u X C j i 1 r 0 7 / 8 Z J G 0 R b D w w c z r m X u e c E i e A a H O f L W l h c W l 5 Z L a w V 1 z c 2 t 7 b t n d 2 G j l N F W Z 3 G I l a t g G g m u G R 1 4 C B Y K 1 G M R I F g z W B w m f n N O 6 Y 0 j + U t D B P m R 6 Q n e c g p A S N 1 7 E M P 2 D 0 E I a 7 g E + x F B P q U C H y N P Y j x j 9 W x S 0 7 Z m Q D P E z c n J Z S j 1 r E / v W 5 M 0 4 h J o I J o 3 X a d B P w R U c C p Y O O i l 2 q W E D o g P d Y 2 V J K I a X 8 0 S T P G R 0 b p 4 j B W 5 k n A E / X 3 x o h E W g + j w E x m 9 + p Z L x P / 8 9 o p h O f + i M s k B S b p 9 K M w F d h E z a r B X a 4 Y B T E 0 h F D F z a 2 Y 9 o k i F E y B R V O C O x t 5 n j Q q Z d f w m 9 N S 9 S K v o 4 D 2 0 Q E 6 R i 4 6 Q 1 V 0 h W q o j i h 6 Q E / o Bb 1 a j 9 a z 9 W a 9 T 0 c X r H x n D / 2 B 9 f E N v 5 2 a E w = = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " b o 0 j 1 L j g t q+ 2 H x e R n h h o l o b 7 S 8 E = " > A A A C D X i c b V D L S s N A F J 3 4 r P U V d e l m s A q C U J I i 6 L L o R l x V s A 9 o Q p l M J + 3 Q y S T M 3 I i l 9 A f c + C t u X C j i 1 r 0 7 / 8 Z J G 0 R b D w w c z r m X u e c E i e A a H O f L W l h c W l 5 Z L a w V 1 z c 2 t 7 b t n d 2 G j l N F W Z 3 G I l a t g G g m u G R 1 4 C B Y K 1 G M R I F g z W B w m f n N O 6 Y 0 j + U t D B P m R 6 Q n e c g p A S N 1 7 E M P 2 D 0 E I a 7 g E + x F B P q U C H y N P Y j x j 9 W x S 0 7 Z m Q D P E z c n J Z S j 1 r E / v W 5 M 0 4 h J o I J o 3 X a d B P w R U c C p Y O O i l 2 q W E D o g P d Y 2 V J K I a X 8 0 S T P G R 0 b p 4 j B W 5 k n A E / X 3 x o h E W g + j w E x m 9 + p Z L x P / 8 9 o p h O f + i M s k B S b p 9 K M w F d h E z a r B X a 4 Y B T E 0 h F D F z a 2 Y 9 o k i F E y B R V O C O x t 5 n j Q q Z d f w m 9 N S 9 S K v o 4 D 2 0 Q E 6 R i 4 6 Q 1 V 0 h W q o j i h 6 Q E / oB b 1 a j 9 a z 9 W a 9 T 0 c X r H x n D / 2 B 9 f E N v 5 2 a E w = = < / l a t e x i t > amplitudes l a t e x i t s h a 1 _ b a s e 6 4 = " M e 3 I U 5 e s O 6 j Y 3 W X c 7 7 0 pE O r + Q F Y = " > A A A C D X i c b V D L S s N A F J 3 U V 6 2 v q E s 3 g 1 U Q h J K I o M u i G 3 F V w T 6 g C W U y n b R D J 5 M w c y O W 0 h 9 w 4 6 + 4 c a G I W / f u / B s n b R B t P T B w O O d e 5 p 4 T J I J r c J w v q 7 C w u L S 8 U l w t r a 1 v b G 7 Z 2 z s N H a e K s j q N R a x a A d F M c M n q w E G w V q I Y i Q L B m s H g M v O b d 0 x pH s t b G C b M j 0 h P 8 p B T A k b q 2 A c e s H s I Q u z i Y + x F B P q U C H y N P Y j x j 9 W x y 0 7 F m Q D P E z c n Z Z S j 1 r E / v W 5 M 0 4 h J o I J o 3 X a d B P w R U c C p Y O O S l 2 q W E D o g P d Y 2 V J K I a X 8 0 S T P G h 0 b p 4 j B W 5 k n A E / X 3 x o h E W g + j w E x m 9 + p Z L x P / 8 9 o p h O f + i M s k B S b p 9 K M w F d h E z a r B X a 4 Y B T E 0 h F D F z a 2 Y 9 o k i F E y B J V O C O 1 _ b a s e 6 4 = " M e 3 I U 5 e s O 6 j Y 3 W X c 7 7 0 p E O r +Q F Y = " > A A A C D X i c b V D L S s N A F J 3 U V 6 2 v q E s 3 g 1 U Q h J K I o M u i G 3 F V w T 6 g C W U y n b R D J 5 M w c y O W 0 h 9 w 4 6 + 4 c a G I W / f u / B s n b R B t P T B w O O d e 5 p 4 T J I J r c J w v q 7 C w u L S 8 U l w t r a 1 v b G 7 Z 2 z s N H a e K s j q N R a x a A d F M c M n q w E G w V q I Y i Q L B m s H g M v O b d 0 x p H s t b G C b M j 0 h P 8 p B T A k b q 2 A c e s H s I Q u z i Y + x F B P q U C H y N P Y j x j 9 W x y 0 7 F m Q D P E z c n Z Z S j 1 r E / v W 5 M 0 4 h J o I J o 3 X a d B P w R U c C p Y O O S l 2 q W E D o g P d Y 2 V J K I a X8 0 S T P G h 0 b p 4 j B W 5 k n A E / X 3 x o h E W g + j w E x m 9 + p Z L x P / 8 9 o p h O f + i M s k B S b p 9 K M w F d h E z a r B X a 4 Y B T E 0 h F D F z a 2 Y 9 o k i F E y B J V O C O 1 _ b a s e 6 4 = " M e 3 I U 5 e s O 6 j Y 3 W X c 7 7 0 p E O r +Q F Y = " > A A A C D X i c b V D L S s N A F J 3 U V 6 2 v q E s 3 g 1 U Q h J K I o M u i G 3 F V w T 6 g C W U y n b R D J 5 M w c y O W 0 h 9 w 4 6 + 4 c a G I W / f u / B s n b R B t P T B w O O d e 5 p 4 T J I J r c J w v q 7 C w u L S 8 U l w t r a 1 v b G 7 Z 2 z s N H a e K s j q N R a x a A d F M c M n q w E G w V q I Y i Q L B m s H g M v O b d 0 x p H s t b G C b M j 0 h P 8 p B T A k b q 2 A c e s H s I Q u z i Y + x F B P q U C H y N P Y j x j 9 W x y 0 7 F m Q D P E z c n Z Z S j 1 r E / v W 5 M 0 4 h J o I J o 3 X a d B P w R U c C p Y O O S l 2 q W E D o g P d Y 2 V J K I a X8 0 S T P G h 0 b p 4 j B W 5 k n A E / X 3 x o h E W g + j w E x m 9 + p Z L x P / 8 9 o p h O f + i M s k B S b p 9 K M w F d h E z a r B X a 4 Y B T E 0 h F D F z a 2 Y 9 o k i F E y B J V O C O 1 _ b a s e 6 4 = " M e 3 I U 5 e s O 6 j Y 3 W X c 7 7 0 p E O r +Q F Y = " > A A A C D X i c b V D L S s N A F J 3 U V 6 2 v q E s 3 g 1 U Q h J K I o M u i G 3 F V w T 6 g C W U y n b R D J 5 M w c y O W 0 h 9 w 4 6 + 4 c a G I W / f u / B s n b R B t P T B w O O d e 5 p 4 T J I J r c J w v q 7 C w u L S 8 U l w t r a 1 v b G 7 Z 2 z s N H a e K s j q N R a x a A d F M c M n q w E G w V q I Y i Q L B m s H g M v O b d 0 x p H s t b G C b M j 0 h P 8 p B T A k b q 2 A c e s H s I Q u z i Y + x F B P q U C H y N P Y j x j 9 W x y 0 7 F m Q D P E z c n Z Z S j 1 r E / v W 5 M 0 4 h J o I J o 3 X a d B P w R U c C p Y O O S l 2 q W E D o g P d Y 2 V J K I a X8 0 S T P G h 0 b p 4 j B W 5 k n A E / X 3 x o h E W g + j w E x m 9 + p Z L x P / 8 9 o p h O f + i M s k B S b p 9 K M w F d h E z a r B X a 4 Y B T E 0 h F D F z a 2 Y 9 o k i F E y B J V O C O up to O(e m⇡L ) corrections < l a t e x i t s h a 1 _ b a s e 6 4 = " V B c y s l B U w 2 G S u Q 3 A y A F r o H W l G m k = " > A A A C K n i c b V D L T g I x F O 3 g C / E 1 6 t J N I 5 H g Q j J j T H S J u n F h I i b y S B g k n X K B h s 4 j b c d I J s P v u P F X 3 L D Q E L d + i B 1 g o e h J m p y c e 2 7 b c 9 y Q M 6 k s a 2 J k l p Z X V t e y 6 7 m N z a 3 t H X N 3 r y a D S F C o 0o A H o u E S C Z z 5 U F V M c W i E A o j n c q i 7 g + t 0 X n 8 C I V n g P 6 h h C C 2 P 9 H z W Z Z Q o L b X N y 4 K j 4 F n F U Y h V g B P H I 6 p P C Y / v k i I 8 x i d e 2 w k Z v k 2 O H S d X G I 1 G M z M N h A C a X i B x 0 j b z V s m a A v 8 l 9 p z k 0 R y V t j l 2 O g G N P P A V 5 U T K p m 2 F q h U T o R j l k O S c S E J I 6 I D 0 o K m p T z y Q r X g a N c F H W u n g b i D 0 8 R W e q j 8 3 Y u J J O f R c 7 U y j y M V Z K v 4 3 a 0 a q e 9 G K m R 9 G C n w 6 e 6 g b 8 b S V t D f c Y W l k P t S E U M H 0 X z H t E 0 G o 0 u 3 m d A n 2 Y u S / p H Z a s j W / P 8 u X r + Z 1 Z N E B O k R F Z K N z V E Y 3 q I K q i K I X 9 I b e 0 Y f x a o y N i f E 5 s 2 a M + c 4 + + g X j 6 x t w 5 K f l < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " V B c y s l B U w 2 G S u Q 3 A y A F r o H W l G m k = " > A A A C K n i c b V D L T g I x F O 3 g C / E 1 6 t J N I 5 H g Q j J j T H S J u n F h I i b y S B g k n X K B h s 4 j b c d I J s P v u P F X 3 L D Q E L d + i B 1 g o e h J m p y c e 2 7 b c 9 y Q M 6 k s a 2 J k l p Z X V t e y 6 7 m N z a 3 t H X N 3 r y a D S F C o 0 o A H o u E S C Z z 5 U F V M c W i E A o j n c q i 7 g + t 0 X n 8 C I V n g P 6 h h C C 2 P 9 H z W Z Z Q o L b X N y 4 K j 4 F n F U Y h V g B P H I 6 p P C Y / v k i I 8 x i d e 2 w k Z v k 2 O H S d X G I 1 G M z M N h A C a X i B x 0 j b z V s m a A v 8 l 9 p z k 0 R y V t j l 2 O g G N P P A V 5 U T K p m 2 F q h U T o R j l k O S c S E J I 6 I D 0 o K m p T z y Q r X g a N c F H W u n g b i D 0 8 R W e q j 8 3 Y u J J O f R c 7 U y j y M V Z K v 4 3 a 0 a q e 9 G K m R 9 G C n w 6 e 6 g b 8 b S V t D f c Y W l k P t S E U M H 0 X z H t E 0 G o 0 u 3 m d A n 2 Y u S / p H Z a s j W / P 8 u X r + Z 1 Z N E B O k R F Z K N z V E Y 3 q I K q i K I X 9 I b e 0 Y f x a o y N i f E 5 s 2 a M + c 4 + + g X j 6 x t w 5 K f l < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " V B c y s l B U w 2 G S u Q 3 A y A F r o H W l G m k = " > A A A C K n i c b V D L T g I x F O 3 g C / E 1 6 t J N I 5 H g Q j J j T H S J u n F h I i b y S B g k n X K B h s 4 j b c d I J s P v u P F X 3 L D Q E L d + i B 1 g o e h J m p y c e 2 7 b c 9 y Q M 6 k s a 2 J k l p Z X V t e y 6 7 m N z a 3 t H X N 3 r y a D S F C o 0 o A H o u E S C Z z 5 U F V M c W i E A o j n c q i 7 g + t 0 X n 8 C I V n g P 6 h h C C 2 P 9 H z W Z Z Q o L b X N y 4 K j 4 F n F U Y h V g B P H I 6 p P C Y / v k i I 8 x i d e 2 w k Z v k 2 O H S d X G I 1 G M z M N h A C a X i B x 0 j b z V s m a A v 8 l 9 p z k 0 R y V t j l 2 O g G N P P A V 5 U T K p m 2 F q h U T o R j l k O S c S E J I 6 I D 0 o K m p T z y Q r X g a N c F H W u n g b i D 0 8 R W e q j 8 3 Y u J J O f R c 7 U y j y M V Z K v 4 3 a 0 a q e 9 G K m R 9 G C n w 6 e 6 g b 8 b S V t D f c Y W l k P t S E U M H 0 X z H t E 0 G o 0 u 3 m d A n 2 Y u S / p H Z a s j W / P 8 u X r + Z 1 Z N E B O k R F Z K N z V E Y 3 q I K q i K I X 9 I b e 0 Y f x a o y N i f E 5 s 2 a M + c 4 + + g X j 6 x t w 5 K f l < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " V B c y s l B U w 2 G S u Q 3 A y A F r o H W l G m k = " > A A A C K n i c b V D L T g I x F O 3 g C / E 1 6 t J N I 5 H g Q j J j T H S J u n F h I i b y S B g k n X K B h s 4 j b c d I J s P v u P F X 3 L D Q E L d + i B 1 g o e h J m p y c e 2 7 b c 9 y Q M 6 k s a 2 J k l p Z X V t e y 6 7 m N z a 3 t H X N 3 r y a D S F C o 0 o A H o u E S C Z z 5 U F V M c W i E A o j n c q i 7 g + t 0 X n 8 C I V n g P 6 h h C C 2 P 9 H z W Z Z Q o L b X N y 4 K j 4 F n F U Y h V g B P H I 6 p P C Y / v k i I 8 x i d e 2 w k Z v k 2 O H S d X G I 1 G M z M N h A C a X i B x 0 j b z V s m a A v 8 l 9 p z k 0 R y V t j l 2 O g G N P P A V 5 U T K p m 2 F q h U T o R j l k O S c S E J I 6 I D 0 o K m p T z y Q r X g a N c F H W un g b i D 0 8 R W e q j 8 3 Y u J J O f R c 7 U y j y M V Z K v 4 3 a 0 a q e 9 G K m R 9 G C n w 6 e 6 g b 8 b S V t D f c Y W l k P t S E U M H 0 X z H t E 0 G o 0 u 3 m d A n 2 Y u S / p H Z a s j W / P 8 u X r + Z 1 Z N E B O k R F Z K N z V E Y 3 q I K q i K I X 9 I b e 0 Y f x a o y N i f E 5 s 2 a M + c 4 + + g X j 6 x t w 5 K f l < / l a t e x i t > Lüscher formalism m 2 b g 7 a + s P D w z g w 7 8 w a p 4 N q 4 7 r d T W l l d W 9 8 o b 1 a 2 t n d 2 9 6 r 7 B 4 8 6 y R T g g m d 4 h T c n d V 6 c d + d j 3 l p y i p l D + C P n 8 w c w q Y 7 N < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " A m z L Q u H 8 i d z p z 8 D U H 8 b p 4 8 f 9 v P I = " > A A A B 7 n i c b Z B N S 8 N A E I Y n 9 a v W r 6 p H L 4 t F 8 K A l E U G P R S 8 e K 1 h b a E P Z b C f t k s 0 m 7 G 6 E E v o j v H h Q x K u / x 5 v / x m 2 b g 7 a + s P D w z g w 7 8 w a p 4 N q 4 7 r d T W l l d W 9 8 o b 1 a 2 t n d 2 9 6 r 7 B 4 8 6 y R T g g m d 4 h T c n d V 6 c d + d j 3 l p y i p l D + C P n 8 w c w q Y 7 N < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " A m z L Q u H 8 i d z p z 8 D U H 8 b p 4 8 f 9 v P I = " > A A A B 7 n i c b Z B N S 8 N A E I Y n 9 a v W r 6 p H L 4 t F 8 K A l E U G P R S 8 e K 1 h b a E P Z b C f t k s 0 m 7 G 6 E E v o j v H h Q x K u / x 5 v / x m 2 b g 7 a + s P D w z g w 7 8 w a p 4 N q 4 7 r d T W l l d W 9 8 o b 1 a 2 t n d 2 9 6 r 7 B 4 8 6 y R T g g m d 4 h T c n d V 6 c d + d j 3 l p y i p l D + C P n 8 w c w q Y 7 N < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " A m z L Q u H 8 i d z p z 8 D U H 8 b p 4 8 f 9 v P I = " > A A A B 7 n i c b Z B N S 8 N A E I Y n 9 a v W r 6 p H L 4 t F 8 K A l E U G P R S 8 e K 1 h b a E P Z b C f t k s 0 m 7 G 6 E E v o j v H h Q x K u / x 5 v / x m 2 b g 7 a + s P D w z g w 7 8 w a p 4 N q 4 7 r d T W l l d W 9 8 o b 1 a 2 t n d 2 9 6 r 7 B 4 8 6 y R T y l C g 1 S U J T W S 6 u 5 r 1 S / M 8 b Z D q + 9 H P K 0 0 w T j m e D 4 o x B L W C Z C 4 y o J F i z i S E I S 2 p 2 h X i E J M L a p F c 3 I b j z J y + S 7 l n T N f z u v N G 6 r u K o g U N w D E 6 B C y 5 A C 9 y C N u g A D B 7 B M 3 g F b 9 a T 9 W K 9 W y l C g 1 S U J T W S 6 u 5 r 1 S / M 8 b Z D q + 9 H P K 0 0 w T j m e D 4 o x B L W C Z C 4 y o J F i z i S E I S 2 p 2 h X i E J M L a p F c 3 I b j z J y + S 7 l n T N f z u v N G 6 r u K o g U N w D E 6 B C y 5 A C 9 y C N u g A D B 7 B M 3 g F b 9 a T 9 W K 9 W y l C g 1 S U J T W S 6 u 5 r 1 S / M 8 b Z D q + 9 H P K 0 0 w T j m e D 4 o x B L W C Z C 4 y o J F i z i S E I S 2 p 2 h X i E J M L a p F c 3 I b j z J y + S 7 l n T N f z u v N G 6 r u K o g U N w D E 6 B C y 5 A C 9 y C N u g A D B 7 B M 3 g F b 9 a T 9 W K 9 W FIG. 2. (a) The kinematics of the process considered, as described in the text.(b) The triangle diagram that appears due to the single insertion of the external current.(c) The diagrammatic representation of the 2 + J → 2 transition amplitude.The black circles depict the → 2 scattering amplitude in the absence of the current.The crossed circles represent various couplings of the external current.Those with two hadronic external legs correspond to the standard 1 + J → 1 matrix element, while those with four external hadrons represent diagrams that are two-particle irreducible in the channel carrying the total momentum.The solid lines denote fully-dressed propagators of the low-energy degrees of freedom (the hadrons).In the second line we separate this into a contribution with on-shell singularities together with the divergence-free amplitude.(d) The diagrammatic representation of the new finite-volume function Gµ 1 ...µn , defined in Eq. (27).
FIG.4.Example plots for G σ (P f , Pi, L) for the case of Pi = P f = (E, P) with P = (2π/L)[001].The real parts are shown as solid curves and the imaginary parts (multiplied by a factor to make the functional form visible) are dashed.At all non-interacting energy levels the function diverges as a double pole with a positive coefficient.

FIG. 5 .
FIG. 5. Examples of the numerical integral I σN ( ᾱ, P f , Pi), plotted as a function of E f with all other arguments fixed.In the left panel we consider the case of f , m f = i, mi = 00 with no factors of kµ in the numerator.For these kinematics the function requires no subtractions and we directly evaluate IN at ᾱ = 1/3 4 , for fixed external 3-momenta and evenly spaced values of E i , as indicated.In the right panel we plot the function with a numerator factor of 3(k z ) 2 k0.In this case one requires Pauli-Villars-like subtractions, as described in the main text and summarized in the caption.
1 p P 1 Y r 1 b H / P W F a u c q a E / s D 5 / A G I G l T w = < / l a t e x i t > [+0.20, +0.25] < l a t e x i t s h a 1 _ b a s e 6 4 = " 1 l k c k / U c H F o M + o v a l Y x 3 L 8 9 h 9 / g = " > A A A B 9 X i c d Z D L S g M x F I Y z 9 V b r r e r S T b A I g j J k S t V 2 V 3 T j s o K 9 w H Q s m T T T h m Y y Q 5 J R y t D 3 c O N C E b e + i z v f x k x b i o o e C P n 5 / n P I y e / H n C m N 0 K e V W 1 p e W V 3 L r x c 2 N r e 2 d 4 q 7 e y 0 / l a t e x i t > [ 0.25, 0.15] < l a t e x i t s h a 1 _ b a s e 6 4 = " V 4 B K K 0 T Y / 8 v 1 F 0 s 2 L h / o y J T C a N E = " > A A A B + H i c d Z D L S g M x F I Y z 9 V b r p a M u 3 Q S L 4 M I O m d J q u y u 6 c V n B X m A 6 l E y a a U M z m S H J C L X 0 S d y 4 U M S t j + L O t z G 9 U F T 0 h 8 D h + 8 / h n P x B w p n S C H 1 a m b X 1 j c 2 t 7 H Z u Z 3 d v P 2 8 f H L Z U n E p C m y T m s e w E W F H

1 [
+0.20, +0.25] < l a t e x i t s h a 1 _ b a s e 6 4 = " d R r C n C 0 a S y R L 2 A p 7 Z g 8 U i Y o 9 e n 8 = " 5 8 0 S 4 6 L H P e 6 X K x d z O r I g 3 1 w A I 6 A C 8 5 B D V y B O m g A A u 7 B I 3 g G L 9 a D 9 W S 9 W m / T 0 Z w 1 2 9 k F P 2 C 9 f w F h H Z R l < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " d R r C n C 0 a S y R L 2 A p 7 Z g 8 U i Y o 9 e n 8 = " 5 8 0 S 4 6 L H P e 6 X K x d z O r I g 3 1 w A I 6 A C 8 5 B D V y B O m g A A u 7 B I 3 g G L 9 a D 9 W S 9 W m / T 0 Z w 1 2 9 k F P 2 C 9 f w F h H Z R l < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " d R r C n C 0 a S y R L 2 A p 7 Z g 8 U i Y o 9 e n 8 = " 5 8 0 S 4 6 L H P e 6 X K x d z O r I g 3 1 w A I 6 A C 8 5 B D V y B O m g A A u 7 B I 3 g G L 9 a D 9 W S 9 W m / T 0 Z w 1 2 9 k F P 2 C 9 f w F h H Z R l < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " d R r C n C 0 a S y R L 2 A p 7 Z g 8 U i Y o 9 e n 8 = " 5 8 0 S 4 6 L H P e 6 X K x d z O r I g 3 1 w A I 6 A C 8 5 B D V y B O m g A A u 7 B I 3 g G L 9 a D 9 W S 9 W m / T 0 Z w 1 2 9 k F P 2 C 9 f w F h H Z R l < / l a t e x i t > [+0.10, +0.15] 8 e w T N 4 s R 6 s J + v V e p u O 5 q z Z z i 7 4 A e v 9 C 1 4 J l G M = < / l a t e x i t > [+✏, +0.05] < l a t e x i t s h a 1 _ b a s e 6 4 = " F b t 5 5 W e T G 0 b w y J A u D G N p n t B w Q r 9 P Z D h S a h g F p n O 8 s f r t j e F f n p v q s O J n T C S p p o J M P w p T D n U M x 4 H A L p O U a D 4 0 B S a S m V 0 h 6 W O J i T a x 5 U 0 I 8 9 v / L 5 q n t o N s 5 6 Z c r F 3 O 4 s i B A 3 A I j o E D L k A N X I M 6 a A A C 7 s E j e A Y v 1 o P 1 Z L 1 a b 9 P W B W s 2 s w d + y H r / A g f a l y U = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " 6 6 Z X k 1 3 P w 0 s 8 F z 0 y 5 W e T G 0 b w y J A u D G N p n t B w Q r 9 P Z D h S a h g F p n O 8 s f r t j e F f n p v q s O J n T C S p p o J M P w p T D n U M x 4 H A L p O U a D 4 0 B S a S m V 0 h 6 W O J i T a x 5 U 0 I 8 9 v / L 5 q n t o N s 5 6 Z c r F 3 O 4 s i B A 3 A I j o E D L k A N X I M 6 a A A C 7 s E j e A Y v 1 o P 1 Z L 1 a b 9 P W B W s 2 s w d + y H r / A g f a l y U = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " 6 6 Z X k 1 3 P w 0 s 8 F z 0 y 5 W e T G 0 b w y J A u D G N p n t B w Q r 9 P Z D h S a h g F p n O 8 s f r t j e F f n p v q s O J n T C S p p o J M P w p T D n U M x 4 H A L p O U a D 4 0 B S a S m V 0 h 6 W O J i T a x 5 U 0 I 8 9 v / L 5 q n t o N s 5 6 Z c r F 3 O 4 s i B A 3 A I j o E D L k A N X I M 6 a A A C 7 s E j e A Y v 1 o P 1 Z L 1 a b 9 P W B W s 2 s w d + y H r / A g f a l y U = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " 6 6 Z X k 1 3 P w 0 s 8 F z 0 y 5 W e T G 0 b w y J A u D G N p n t B w Q r 9 P Z D h S a h g F p n O 8 s f r t j e F f n p v q s O J n T C S p p o J M P w p T D n U M x 4 H A L p O U a D 4 0 B S a S m V 0 h 6 W O J i T a x 5 U 0 I 8 9 v / L 5 q n t o N s 5 6 Z c r F 3 O 4 s i B A 3 A I j o E D L k A N X I M 6 a A A C 7 s E j e A Y v 1 o P 1 Z L 1 a b 9 P W B W s 2 s w d + y H r / A g f a l y U = < / l a t e x i t > [ 0.15, 0.10] < l a t e x i t s h a 1 _ b a s e 6 4 = " B n z v U 4 C e x x + Y 6 8 l 6 t d 6 m o z l r t r M L f s B 6 / w J n b 5 R p < / l a t e x i t > mL = 6 = [10; 10; µ = 0] l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " Y HX F U h 1 t Y U P W 0 u h h y Q L 2 a X R 8 u m U = " > A A A B 8 n i c b V D L S g N B E J y N r x h f U Y9 e B o P g K e 6 K o M e g C B 4 j m A c k a 5 i d z C Z D 5 r H M 9 A p h y W d 4 8 a C I V 7 / G m 3 / j J N m D J h Y 0 F F X d d H d F i e A W f P / b K 6 y s r q 1 Y x I a 8 c y c p 2 S w N A u e l P x P 6+ T Q n w V Z l w l K T B F 5 4 v i V G D Q e P o / 7 n P D K I i x I 4 Q a 7 m 7 F d E g M o e B S K r k Q g s W X l 0 n z v B r 4 1 e D + o l K 7 z u M o o i N 0 j E 5 R g C 5 R D d 2 h O m o g i j R 6 R q / o z Q P v x X v 3 P u a t B S + f O U R / 4 H 3 + A P J I k Q U = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " Y H X F U h 1 t Y U P W 0 u h h y Q L 2 a X R 8 u m U = " > A A A B 8 n i c b V D L S g N B E J y N r x h f U Y9 e B o P g K e 6 K o M e g C B 4 j m A c k a 5 i d z C Z D 5 r H M 9 A p h y W d 4 8 a C I V 7 / G m 3 / j J N m D J h Y 0 F F X d d H d F i e A W f P / b K 6 y s r q 1 Y x I a 8 c y c p 2 S w N A u e l P x P 6 + T Q n w V Z l w l K T B F 5 4 v i V G D Q e P o / 7 n P D K I i x I 4 Q a 7 m 7 F d E g M o e B S K r k Q g s W X l 0 n z v B r 4 1 e D + o l K 7 z u M o o i N 0 j E 5 R g C 5 R D d 2 h O m o g i j R 6 R q / o z Q P v x X v 3 P u a t B S + f O U R / 4 H 3 + A P J I k Q U = < / l a t e x i t > E ?f m < l a t e x i t s h a 1 _ b a s e 6 4 = " k Z Y b L g yI f W v Y S s n q j S 2 o u W K Z M I M = " > A A A B / H i c b V B N S 8 N A E N 3 4 W e t X t E c v i 0 X w V B I R 9 F g U w W M F + w F t D J v t p l 2 6 m 4 T d i R B C / C t e P C j i 1 R / i z X / j t s 1 B W x 8 M P N 6 b Y W Z e k A i u w X G + r Z X V t f W N z c p W d X t n d 2 / f P j j s 6 D h V l L V p L G L V C 4 h m g k e s D R w E 6 y W K E R k I 1 g 0 m 1 1 O / + 8 i U 5 n F 0 D 1 n C P E l G E Q 8 5 J W A k 3 6 4 N Q k V o f u O H D w M N R B W 5 L H y 7 7 j S c G f A y c U t S R y V a v v 0 1 G M Y 0 l S w C K o j W f d d J w M u J A k 4 F K 6 q D V L O E 0 A k Z s b 6 h E Z F M e / n s + A K f G G W I w 1 i Z i g D P 1 N 8 T O Z F a Z z I w n Z L A W C 9 6 U / E / r 5 9 C e O n l P E p S Y B G d L w p T g S H G 0 y T w k C t G Q W S G E K q 4 u R X T M T F p g M m r a k J w F 1 9 e J p 2 z h u s 0 3 L v z e v O q j K O C j t A x O k U u u k B N d I t a q I 0 o y t A z e k V v 1 p P 1 Y r 1 b H / P W F a u c q a E / s D 5 / A G I G l T w = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " k Z Y b L g y I f W v Y S s n q j S 2 o u W K Z M I M = " > A A A B / H i c b V B N S 8 N A E N 3 4 W e t X t E c v i 0 X w V B I R 9 F g U w W M F + w F t D J v t p l 2 6 m 4 T d i R B C / C t e P C j i 1 R / i z X / j t s 1 B W x 8 M P N 6 b Y W Z e k A i u w X G + r Z X V t f W N z c p W d X t n d 2 / f P j j s 6 D h V l L V p L G L V C 4 h m g k e s D R w E 6 y W K E R k I 1 g 0 m 1 1 O / + 8 i U 5 n F 0 D 1 n C P E l G E Q 8 5 J W A k 3 6 4 N Q k V o f u O H D w M N R B W 5 L H y 7 7 j S c G f A y c U t S R y V a v v 0 1 G M Y 0 l S w C K o j W f d d J w M u J A k 4 F K 6 q D V L O E 0 A k Z s b 6 h E Z F M e / n s + A K f G G W I w 1 i Z i g D P 1 N 8 T O Z F a Z z I w n Z L A W C 9 6 U / E / r 5 9 C e O n l P E p S Y B G d L w p T g S H G 0 y T w k C t G Q W S G E K q 4 uR X T M T F p g M m r a k J w F 1 9 e J p 2 z h u s 0 3 L v z e v O q j K O C j t A x O k U u u k B N d I t a q I 0 o y t A z e k V v 1 p P 1 Y r 1 b H / P W F a u c q a E / s D 5 / A G I G l T w = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 64 = " k Z Y b L g y I f W v Y S s n q j S 2 o u W K Z M I M = " > A A A B / H i c b V B N S 8 N A E N 3 4 W e t X t E c v i 0 X w V B I R 9 F g U w W M F + w F t D J v t p l 2 6 m 4 T d i R B C / C t e P C j i 1 R / i z X / j t s 1 B W x 8 M P N 6 b Y W Z e k A i u w X G + r Z X V t f W N z c p W d X t n d 2 / f P j j s 6 D h V l L V p L G L V C 4 h m g k e s D R w E 6 y W K E R k I 1 g 0 m 1 1 O / + 8 i U 5 n F 0 D 1 n C P E l G E Q 8 5 J W A k 3 6 4 N Q k V o f u O H D w M N R B W 5 L H y 7 7 j S c G f A y c U t S R y V a v v 0 1 G M Y 0 l S w C K o j W f d d J w M u J A k 4 F K 6 q D V L O E 0 A k Z s b 6 h E Z F M e / n s + A K f G G W I w 1 i Z i g D P 1 N 8 T O Z F a Z z I w n Z L A W C 9 6 U / E / r 5 9 C e O n l P E p S Y B G d L w p T g S H G 0 y T w k C t G Q W S G E K q 4 uR X T M T F p g M m r a k J w F 1 9 e J p 2 z h u s 0 3 L v z e v O q j K O C j t A x O k U u u k B N d I t a q I 0 o y t A z e k V v 1 p P 1 Y r 1 b H / P W F a u c q a E / s D 5 / A G I G l T w = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 64 = " k Z Y b L g y I f W v Y S s n q j S 2 o u W K Z M I M = " > A A A B / H i c b V B N S 8 N A E N 3 4 W e t X t E c v i 0 X w V B I R 9 F g U w W M F + w F t D J v t p l 2 6 m 4 T d i R B C / C t e P C j i 1 R / i z X / j t s 1 B W x 8 M P N 6 b Y W Z e k A i u w X G + r Z X V t f W N z c p W d X t n d 2 / f P j j s 6 D h V l L V p L G L V C 4 h m g k e s D R w E 6 y W K E R k I 1 g 0 m 1 1 O / + 8 i U 5 n F 0 D 1 n C P E l G E Q 8 5 J W A k 3 6 4 N Q k V o f u O H D w M N R B W 5 L H y 7 7 j S c G f A y c U t S R y V a v v 0 1 G M Y 0 l S w C K o j W f d d J w M u J A k 4 F K 6 q D V L O E 0 A k Zs b 6 h E Z F M e / n s + A K f G G W I w 1 i Z i g D P 1 N 8 T O Z F a Z z I w n Z L A W C 9 6 U / E / r 5 9 C e O n l P E p S Y B G d L w p T g S H G 0 y T w k C t G Q W S G E K q 4 u R X T M T F p g M m r a k J w F 1 9 e J p 2 z h u s 0 3 L v z e v O q j K O C j t A x O k U u u k B N d I t a q I 0 o y t A z e k V v 1 p P 1 Y r 1 b H / P W F a u c q a E / s D 5 / A G I G l T w = < / l a t t e x i t s h a 1 _ b a s e 6 4 = " 1 q W6 3 d W D d f r s 1 9 0 O 0 O 4 F + v u D y S o = " > A A A B 9 X i c b V B N S w M x E J 3 1 s 9 a v q k c v w S J 4 K l k R F E Q o e v F Yw X 5 A u 5 Z s m m 1 D k + y S Z J W y 9 H 9 4 8 a C I V / + L N / + N a b s H b X 0 w 8 H h v h p l 5 Y S K 4 s R h / e 0 v L K 6 t r 6 4 W N 4 u b W 9 s 5 u a W + / Y e J U U 1 a n 3 t y I 6 I J p Q 6 4 I q u h D 8 + Z c X S e O 0 4 u O K f 3 d W r l 7 n c R T g E I 7 g B H w 4 h y r c Q g 3 q Q E H D M 7 z C m / f k v X j v 3 s e s d c n L Z w 7 g D 7 z P H 8 4 p k W M = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " 1 q W 6 3d W D d f r s 1 9 0 O 0 O 4 F + v u D y S o = " > A A A B 9 X i c b V B N S w M x E J 3 1 s 9 a v q k c v w S J 4 K l k R F E Q o e v F Yw X 5 A u 5 Z s m m 1 D k + y S Z J W y 9 H 9 4 8 a C I V / + L N / + N a b s H b X 0 w 8 H h v h p l 5 Y S K 4 s R h / e 0 v L K 6 t r 6 4 W N 4 u b W 9 s 5 u a W + / Y e J U U 1 a n 3 t y I 6 I J p Q 6 4 I q u h D 8 + Z c X S e O 0 4 u O K f 3 d W r l 7 n c R T g E I 7 g B H w 4 h y r c Q g 3 q Q E H D M 7 z C m / f k v X j v 3 s e s d c n L Z w 7 g D 7 z P H 8 4 p k W M = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " 1 q W 6 3d W D d f r s 1 9 0 O 0 O 4 F + v u D y S o = " > A A A B 9 X i c b V B N S w M x E J 3 1 s 9 a v q k c v w S J 4 K l k R F E Q o e v F Yw X 5 A u 5 Z s m m 1 D k + y S Z J W y 9 H 9 4 8 a C I V / + L N / + N a b s H b X 0 w 8 H h v h p l 5 Y S K 4 s R h / e 0 v L K 6 t r 6 4 W N 4 u b W 9 s 5 u a W + / Y e J U U 1 a n 3 t y I 6 I J p Q 6 4 I q u h D 8 + Z c X S e O 0 4 u O K f 3 d W r l 7 n c R T g E I 7 g B H w 4 h y r c Q g 3 q Q E H D M 7 z C m / f k v X j v 3 s e s d c n L Z w 7 g D 7 z P H 8 4 p k W M = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " 1 q W 6 3d W D d f r s 1 9 0 O 0 O 4 F + v u D y S o = " > A A A B 9 X i c b V B N S w M x E J 3 1 s 9 a v q k c v w S J 4 K l k R F E Q o e v F Yw X 5 A u 5 Z s m m 1 D k + y S Z J W y 9 H 9 4 8 a C I V / + L N / + N a b s H b X 0 w 8 H h v h p l 5 Y S K 4 s R h / e 0 v L K 6 t r 6 4 W N 4 u b W 9 s 5 u a W + / Y e J U U 1 a n 3 t y I 6 I J p Q 6 4 I q u h D 8 + Z c X S e O 0 4 u O K f 3 d W r l 7 n c R T g E I 7 g B H w 4 h y r c Q g 3 q Q E H D M 7 z C m / f k v X j v 3 s e s d c n L Z w 7 g D 7 z P H 8 4 p k W M = < / l a t e x i t >P i = [000] < l a t e x i t s h a 1 _ b a s e 6 4 = " Q X m G O o k r 6 n 5 D 3 S r E J V B k e Q J + l R w = " > A A A B / H i c b V D L S s N A F J 3 4 r P U V 7 d L N Y B F c l Y k I u h G K b l x W s A 9 I Q 5 h M J + 3 Q y Y O Z G z G E + i t u X C j i1 g 9 x 5 9 8 4 b b P Q 1 g M X D u f c y 7 3 3 B K k U G g j 5 t l Z W 1 9 Y 3 N i t b 1 e 2 d 3 b 1 9 + + C w o 5 N M M d 5 m 5 m X T O G g 5 p O H f n 9 e Z 1 G U c F H a F j d I o c d I G a 6 B a 1 U B s x l K N n 9 I r e r C f r x X q 3 P u a t K 1 Y 5 U 0 N / Y H 3 + A G 9 H k 0 8 = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " Q X m G O o k r 6 n 5 D 3 S r E J V B k e Q J + l R w = " > A A A B / H i c b V D L S s N A F J 3 4 r P U V 7 d L N Y B F c l Y k I u h G K b l x W s A 9 I Q 5 h M J + 3 Q y Y O Z G z G E + i t u X C j i1 g 9 x 5 9 8 4 b b P Q 1 g M X D u f c y 7 3 3 B K k U G g j 5 t l Z W 1 9 Y 3 N i t b 1 e 2 d 3 b 1 9 + + C w o 5 N M M d 5 m 5 m X T O G g 5 p O H f n 9 e Z 1 G U c F H a F j d I o c d I G a 6 B a 1 U B s x l K N n 9 I r e r C f r x X q 3 P u a t K 1 Y 5 U 0 N / Y H 3 + A G 9 H k 0 8 = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " Q X m G O o k r 6 n 5 D 3 S r E J V B k e Q J + l R w = " > A A A B / H i c b V D L S s N A F J 3 4 r P U V 7 d L N Y B F c l Y k I u h G K b l x W s A 9 I Q 5 h M J + 3 Q y Y O Z G z G E + i t u X C j i1 g 9 x 5 9 8 4 b b P Q 1 g M X D u f c y 7 3 3 B K k U G g j 5 t l Z W 1 9 Y 3 N i t b 1 e 2 d 3 b 1 9 + + C w o 5 N M M d 5 m 5 m X T O G g 5 p O H f n 9 e Z 1 G U c F H a F j d I o c d I G a 6 B a 1 U B s x l K N n 9 I r e r C f r x X q 3 P u a t K 1 Y 5 U 0 N / Y H 3 + A G 9 H k 0 8 = < / l a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " Q X m G O o k r 6 n 5 D 3 S r E J V B k e Q J + l R w = " > A A A B / H i c b V D L S s N A F J 3 4 r P U V 7 d L N Y B F c l Y k I u h G K b l x W s A 9 I Q 5 h M J + 3 Q y Y O Z G z G E + i t u X C j i1 g 9 x 5 9 8 4 b b P Q 1 g M X D u f c y 7 3 3 B K k U G g j 5 t l Z W 1 9 Y 3 N i t b 1 e 2 d 3 b 1 9 + + C w o 5 N M M d 5 m 5 m X T O G g 5 p O H f n 9 e Z 1 G U c F H a F j d I o c d I G a 6 B a 1 U B s x l K N n 9 I r e r C f r x X q 3 P u a t K 1 Y 5 U 0 N / Y H 3 + A G 9 H k 0 8 = < / l a t e x i t >P f = [001] < l a t e x i t s h a 1 _ b a s e 6 4 = " m + u K H 1 Y H c 8 K C K p R N 4 Y D p e Y h + r j 0 = " > A A A C A X i c b V B N S 8 N A E N 3 4 W e t X 1 I v g Z b U I n k o i g l 6 E o h e P F e w H N C F s t p t 2 6 W Y T d i d i C f X i X / H i Q R G v / g t v / h u 3 b Q 7 a + m D g 8 d 4 M M / P C V H A N j v N t L S w u L a + s l t b K 6 x u b W 9 v 2 z m 5 T J 5 m i r E E T k a h 2 S D Q T X L I G c B C s n S p G 4 l C w V j i 4 H v u t e 6 Y 0 T + Q d D F P m x 6 Q n e c Q p A S M F 9 j 7 2 g D 1 A G O F 6 k H u H 0 Q h f 4 o 7 j u H 5 g V 5 y q M w G e J 2 5 B K q h A P b C / v G 5 C s 5 h J o I J o 3 X G d F P y c K O B U s F H Z y z R L C R 2 Q H u s Y K k n M t J 9 P P h j h Y 6 N 0 c Z Q o U x L w R P 0 9 k Z N Y 6 2 E cm s 6 Y Q F / P e m P x P 6 + T Q X T h 5 1 y m G T B J p 4 u i T G B I 8 D g O 3 O W a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " m +u K H 1 Y H c 8 K C K p R N 4 Y D p e Y h + r j 0 = " > A A A C A X i c b V B N S 8 N A E N 3 4 W e t X 1 I v g Z b U I n k o i g l 6 E o h e P F e w H N C F s t p t 2 6 W Y T d i d i C f X i X / H i Q R G v / g t v / h u 3 b Q 7 a + m D g 8 d 4 M M / P C V H A N j v N t L S w u L a + s l t b K 6 x u b W 9 v 2 z m 5 T J 5 m i r E E T k a h 2 S D Q T X L I G c B C s n S p G 4 l C w V j i 4 H v u t e 6 Y 0 T + Q d D F P m x 6 Q n e c Q p A S M F 9 j 7 2 g D 1 A G O F 6 k H u H 0 Q h f 4 o 7 j u H 5 g V 5 y q M w G e J 2 5 B K q h A P b C / v G 5 C s 5 h J o I J o 3 X G d F P y c K O B U s F H Z y z R L C R 2 Q H u s Y K k n M t J 9 P P h j h Y 6 N 0 c Z Q o U x L w R P 0 9 k Z N Y 6 2 E cm s 6 Y Q F / P e m P x P 6 + T Q X T h 5 1 y m G T B J p 4 u i T G B I 8 D g O 3 O W a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " m +u K H 1 Y H c 8 K C K p R N 4 Y D p e Y h + r j 0 = " > A A A C A X i c b V B N S 8 N A E N 3 4 W e t X 1 I v g Z b U I n k o i g l 6 E o h e P F e w H N C F s t p t 2 6 W Y T d i d i C f X i X / H i Q R G v / g t v / h u 3 b Q 7 a + m D g 8 d 4 M M / P C V H A N j v N t L S w u L a + s l t b K 6 x u b W 9 v 2 z m 5 T J 5 m i r E E T k a h 2 S D Q T X L I G c B C s n S p G 4 l C w V j i 4 H v u t e 6 Y 0 T + Q d D F P m x 6 Q n e c Q p A S M F 9 j 7 2 g D 1 A G O F 6 k H u H 0 Q h f 4 o 7 j u H 5 g V 5 y q M w G e J 2 5 B K q h A P b C / v G 5 C s 5 h J o I J o 3 X G d F P y c K O B U s F H Z y z R L C R 2 Q H u s Y K k n M t J 9 P P h j h Y 6 N 0 c Z Q o U x L w R P 0 9 k Z N Y 6 2 E cm s 6 Y Q F / P e m P x P 6 + T Q X T h 5 1 y m G T B J p 4 u i T G B I 8 D g O 3 O W a t e x i t > < l a t e x i t s h a 1 _ b a s e 6 4 = " m +u K H 1 Y H c 8 K C K p R N 4 Y D p e Y h + r j 0 = " > A A A C A X i c b V B N S 8 N A E N 3 4 W e t X 1 I v g Z b U I n k o i g l 6 E o h e P F e w H N C F s t p t 2 6 W Y T d i d i C f X i X / H i Q R G v / g t v / h u 3 b Q 7 a + m D g 8 d 4 M M / P C V H A N j v N t L S w u L a + s l t b K 6 x u b W 9 v 2 z m 5 T J 5 m i r E E T k a h 2 S D Q T X L I G c B C s n S p G 4 l C w V j i 4 H v u t e 6 Y 0 T + Q d D F P m x 6 Q n e c Q p A S M F 9 j 7 2 g D 1 A G O F 6 k H u H 0 Q h f 4 o 7 j u H 5 g V 5 y q M w G e J 2 5 B K q h A P b C / v G 5 C s 5 h J o I J o 3 X G d F P y c K O B U s F H Z y z R L C R 2 Q H u s Y K k n M t J 9 P P h j h Y 6 N 0 c Z Q o U x L w R P 0 9 k Z N Y 6 2 E cm s 6 Y Q F / P e m P x P 6 + T Q X T h 5 1 y m G T B J p 4 u i T G B I 8 D g O 3 O W FIG. 10. (a)The finite-volume three-point correlator used to derive the 2 + J → 2 formalism.All symbols common to Fig.2have the same definition.The open circles, new to this figure, denote Bethe-Salpeter kernels, defined to include all diagrams besides the s-channel two-particle reducible set, shown explicitly.The two-particle loops shown explicitly are evaluated in a finite volume as indicated with the V label.(b) The finite-volume loop within the correlator that leads to the appearance of the G-function.

6 and 7 . 1 M
The task is to study discontinuities of the integral I A (P f , P i , m, (m, x, y) 2 .