Generalized unitarity method for unstable particles

In theories with unstable particles, unitarity is satisfied by the inclusion of only stable states in unitarity sums. Hence unitarity cuts are not to be taken through unstable particles. This raises a challenge to the generalized unitarity method, whose aim is to reconstruct amplitudes by analyzing sets of unitarity cuts. Nevertheless, under some general physical conditions, and perhaps some methodological modifications, we prove that the method is still reliable for one-loop amplitudes containing resonances. We discuss some simple examples which illustrate these features.


INTRODUCTION
Measurements of parameters of the Standard Model to very high precision is a landmark for the physics program of the Large Hadron Collider [1]. In recent years these have triggered an ongoing stream of research activities dedicated to assess precise predictions within perturbative quantum field theories. The immediate consequence is the development of modern tools to address the calculation of loop integrals. In this regard there is a pressing need to obtain a better knowledge of the analytic structure of such integrals to uncover more streamlined methods for their computation. This implies taking into account the investigation of the so-called cuts of internal propagators associated with intermediate particles in a given scattering amplitude. The major importance of cuts in this context is that it allows one to efficiently probe the analytic structure of loop integrals [2,3].
From unitarity constraints, we know that Feynman integrals should be multi-valued functions, whose discontinuities are precisely described by cuts. Actually, scattering amplitudes can be structured in terms of their singularities, so in principle the investigation of branch cuts and other singularities enables one to calculate loop amplitudes. Fourdimensional amplitudes that are uniquely specified by the nature of their branch cuts are said to be cut-constructible. Modern unitarity methods build on Landau conditions [4] in order to use cuts to set up projectors onto a basis of master integrals [5][6][7][8][9][10][11][12]. For recent applications, see Refs. [13][14][15][16][17][18][19][20][21]. A precise definition of cuts is available for given classes of cuts. Among these, we can quote the so-called unitarity cuts which focus on a particular external channel [22][23][24][25]. Historically, the unitarity method [26][27][28][29][30] was established as a systematic framework for one-loop evaluations, and is applicable to both supersymmetric and non-supersymmetric theories.
The standard practice of the unitarity method requires the replacement of two internal propagators by Dirac delta functions which projects the loop momenta they carry onto their on-shell values. On the other hand, in generalized unitarity [5,6,[31][32][33][34][35][36][37], to be briefly reviewed below, one considers additional cut conditions to constrain other momenta to their associated on-shell values. As a consequence, if the momenta carried by more than two massless propagators take their on-shell values, the solutions to the cut conditions are now complex, which implies that the associated delta functions must give zero [6,8]. This observation has led to the concept that cuts should be computed via contour integration so that the associated contours should be suitably deformed in such a way as to encircle the poles of the cut propagators [8,38,39].
The obvious requirement here is that unitarity must be satisfied to all orders in perturbation theory. However, we know that many of our theories possess unstable particles which do not appear as asymptotic states. Should such unstable particles be incorporated nevertheless in unitarity relations? This issue was addressed by Veltman [40][41][42][43][44].
The conclusion is that one should not take cuts through unstable propagators 1 , hence unstable particles are not enclosed by unitarity sums. This can be generalized to encompass also unstable ghostlike resonances emerging in higher-derivative theories such as quadratic gravity and Lee-Wick theories [45]. As the astute reader might have noticed, this might potentially put some obstructions to unitarity methods. Notwithstanding the foregoing remark, in this paper we will discuss how such methods are still solid for the calculation of loop amplitudes that comprise resonances of any type. Here we will use units such that = c = 1. We take the Minkowski metric as η µν = diag(1, −1, −1, −1).
The knowledge of tree amplitudes can be used to seek information about loop integrands. The action of taking loop propagators on-shell is known as a unitarity cut. It comes from the unitarity constraint of the S-matrix, which is a statement on the generalized optical theorem. That is, for an arbitrary process a → b one has that where dΠ f is the Lorentz-invariant phase space measure [51] and the sum runs over all possible sets f of intermediate states and there is an overall delta function associated with energy-momentum conservation. In the above expression, the A's are (invariant) scattering matrix elements. In a perturbative quantum field theory, when an expansion in powers of a small coupling constant exists, this constraint instructs us that the imaginary part of scattering amplitudes at a given order is obtained from the product of lower-order amplitudes. For instance, in the case of one-loop processes, one finds a product of two tree amplitudes on the right-hand side of Eq. (1). Usually this product presupposes the sum over all possible on-shell states that can cross the cut. One fundamental requirement is that only states from the physical spectrum of the theory are allowed to be included in this sum [40,45]. In a unitarity cut, we restrict the loop-momenta to be on-shell and only physical modes are enclosed in the two on-shell amplitudes on the right-hand side of Eq. (1). The cutting rules also consider integrals of any remaining freedom in the loop momentum after prescribing the so-called cut constraints (and, of course, momentum conservation). Unitarity cuts are efficient tools that enable one to relate the pole structure of the integrand with the branch-cut structure of the associated loop integral. Unitarity cuts can also involve more than two cut lines, which implies that several internal lines are taken on-shell. Here we say that we are able to reconstruct amplitudes from sets of generalized unitarity cuts. It turns out that such a set is overcomplete, which means that we can recourse to different strategies for extracting the relevant information. For example, one interesting approach is to use the method of maximal cuts [50,52]. In this case we consider the maximum possible number of cut lines so that each cut furnishes a small piece of information. That is, we begin with generalized cuts possessing the maximum number of cut propagators. We use the information from such cuts to lay out an initial ansatz for the amplitude. Further cuts with reduced number of cut propagators are then considered and their information is systematically gathered in order to improve such an ansatz. The aim is to find an integrand that reproduces all the unitarity cuts. In principle this helps along the construction of the amplitude [50,[53][54][55].
Here we will study generalized unitarity cuts on the level of the integrand, which can be written as a product of on-shell (tree-level or lower-loop) amplitudes. In particular, we are interested in considering maximal cuts consisting of only three-point tree amplitudes, namely: The information from unitarity cuts can be used most efficiently if a complete basis of integrals is known. Indeed, all one-loop amplitudes in D dimensions can be written as a sum of one-loop scalar integrals I m , m = 1, 2, 3, . . . , D [52]: where R denotes rational terms (contributions that do not have branch cuts), C (i) D are coefficients associated with tree-level amplitudes and I (i) m are m-gon scalar integrals. In D = 4, one-loop integrals reduce to a combination of box, triangle, bubble and tadpole scalar integrals [49,[56][57][58][59][60][61][62]. The latter are related to the coefficients C (l) 1 ; such integrals vanish in dimensional regularization when only massless particles circulate in the loop. In D = 4 power counting demonstrates that the scalar box and triangle integrals do not display UV divergences, but IR divergences due to possible massless corners. The bubble integrals have UV divergences but no IR divergences as both corners are massive. Integrated results of such integrals can be found in several places in the literature, see for instance Refs. [52,58,59,61,[63][64][65][66].
Therefore in four dimensions the following expansion is generically valid to any one-loop amplitude: where K i are sums of external momenta and I n are the associated scalar integrals. The coefficients c n are calculated using generalized cuts. For instance, consider a generic one-loop point amplitude written in the basis above. In this section we are working with only stable particles circulating in the loop. If we cut four propagators then the four dimensional integral becomes trivial: where A tree j are tree-level amplitudes and, using a spectral representation [45], are the cut propagators (or positive-frequency Wightman functions) associated with stable particles. For brevity we have absorbed 2π factors into the definition of the loop integral in Eq. (5). Since in the case of stable particles the spectral function σ(s) has a pole at one-particle states, we can also write that (assuming we are not above a given multi-particle threshold) In other words, the "cut of a propagator" means removing its principal part while preserving the delta function imposing the on-shell condition. For simplicity we have taken all internal propagators to have the same mass m, which can be zero. On the other hand, when applied to the master integrals, the quadruple cut selects the contribution from the box integral with momenta K 1 , K 2 , K 3 , K 4 at the corners. Therefore where I 4 (K 1 , K 2 , K 3 , K 4 ) is the associated 4-point box scalar integral: In particular, the quadruple cut of the scalar box integral is a Jacobian factor. This factor appears on both sides of the equation. Hence, comparing both expressions for ∆ 4 A 1−loop , we see that the coefficient c 4 can be expressed as a product of tree-level amplitudes [49]: where the factor of 1/2 emerges as there are exactly two solutions for the set S of cut conditions determined by the four delta functions of the cut propagators. Hence in principle the quadruple cut of the scalar box integral would suffice to calculate the box coefficient. Furthermore, this implies that the maximal cut in this case reads which is a direct proof of the one-loop form of Eq. (2) for stable particles.
In conclusion, equipped with the integral reduction (4) valid to all one-loop amplitudes, and benefiting from the factorization property of the amplitude, by using unitarity methods one is in a position to reconstruct one-loop amplitudes from tree-level information without the often burdensome Feynman diagram expansion. Moreover, we see that the application of the generalized unitarity method requires exploring further discontinuities which implies that a different number of propagators ought to be put on-shell in comparison with textbook unitarity cuts. And this can only be achieved if there is a contribution of an isolated simple pole at p 2 = m 2 (or p 2 = 0 for massless particles) coming from one-particle states -in other words, if the cut propagators have the expected cut structure, as given by Eq. (7).
Actually, one must be more careful when resorting to generalized unitarity, since the solutions to the cut conditions are generally complex, leading to delta functions that trivially yield zero. The solution is to use contour integration. That is, instead of replacing the propagators by delta functions, one must replace the original contour of integration [8]. In summary, the idea is that, as the support of the delta functions is outside the physical region, the integration procedure is implemented in terms of contour integrals in C 4 , the loop momentum being regarded as a complex vector. Such contours are such that their product encircles the poles in the four-dimensional components of the loop momentum. By performing the four-dimensional loop-momentum integral over each contour, the residue at the corresponding encircled pole is attained. In fact, one defines the product of delta functions to generate exactly this contour integral [8].
Of course, the aforementioned operation does not leave expression (4) intact, as there are terms that integrate to zero in the original contour which no longer necessarily vanish if we integrate over general contours in the complex plane. In order to do away with such spurious terms, one evaluates the integral over a suitable linear combination of new contours in such a way that such additional contributions are always projected out. This produces the coefficients of the box integrals as given in Eq (10). For a careful survey of all subtleties associated with this discussion, see Ref. [8]. See also Ref. [12].

Possible issues with unstable particles
In order to understand what could be the general issues involving unstable particles, let us imagine a given reaction consisting of the scattering of a particles a 1 + a 2 + · · · + a n → A + c 1 + c 2 + · · · + c n (12) producing final products given by a collection of c particles and a particle labeled A. If all final products are described by stable particles, then in principle there is no concern in evaluating an on-shell amplitude such as A(a → A + c) in any given order in perturbation theory. However, if A is heavy enough, its coupling to lighter states in the theory, say labeled by b, makes it decay: and now the scattering process we have to consider is a → b + c with the associated on-shell amplitude A(a → b + c).
In this case A enters the calculation as a virtual particle, not as an external state, which implies, by the Feynman rules, the presence of its propagator 1/(p 2 − m 2 A ) in internal lines of given Feynman diagrams. The point is, if the phase space contains the resonance region ( p b ) 2 = m 2 A , then the results calculated from perturbation theory cannot be trusted close to this region.
In other words, basically a diagram with a single internal A propagator must met one of the following conditions [42]: (i)m A < m b , which always happen if A is stable or if A is an unstable particle with threshold E < m b such that it cannot decay into b. In this case perturbation theory can still be trustworthy 3 ; (i)m A > m b , which implies that A is unstable and can decay to b-particles. In this situation, the phase space contains the resonance, and perturbation theory can no longer be generically trusted.
The most direct approach to solve this problem is to consider a resummed form for the propagator of the unstable particle. However, for gauge theories, the resummation procedure must be done carefully, otherwise one might expect to be confronted with issues associated with gauge invariance and gauge-fixing parameter dependence [42][43][44][67][68][69][70][71]. Another possible approach is provided by the so-called complex-mass scheme (CMS) [72][73][74]. In few words, it corresponds to a suitable generalization of the on-shell renormalization scheme. In the latter, the renormalized mass m is specified by demanding p 2 = m 2 to be the pole position associated with the resummed propagator. This is fine for stable particles -for unstable particles, the self-energy acquires an imaginary part, and as a consequence the renormalized mass does not correspond to the pole position. The modification proposed by the complex-mass scheme is the following: Define a complex renormalized massm by requiring that p 2 =m 2 matches the pole position of the resummed propagator for unstable particles. The fact thatm is complex, and therefore cannot be associated with a physical entity, should be of no concern as renormalized parameters in the Lagrangian do not carry any physical meaning [42]. For a recent discussion on the definition of the mass and width of a normal unstable particle, see Ref. [75].
One can show that this modification put forward by the complex-mass scheme avoids the aforementioned issues appearing in gauge theories as it renders unnecessary the resummation of internal propagators [42][43][44]. Indeed, the bare propagator of the unstable particle A (or its scalar part) within this method acquires the form (in momentum space) By where m A andΓ are real, one can prove that the above propagator can be envisaged as the resummed form of the following propagator in a scheme in which the renormalized mass is given by m A [42]: The self-energy obeys Σ(m 2 A ) = im AΓ . That is, the bare propagator within the complex-mass scheme is intrinsically resummed. In addition, notice that, as m AΓ is evaluated throughout the renormalization procedure, we must envisage this quantity as a function of the coupling constant λ describing the interaction between the unstable particle and the lighter states, m AΓ ∼ O(λ). For more technical details concerning the complex-mass scheme, we refer the reader Refs. [42][43][44].
Given the result (14), one can be tempted to think that the complex-mass scheme allows for an adequate spectral representation for the propagator of the unstable particle such that an unambiguous one-particle state contribution can be identified. That this is not straightforward can be seen as follows. Within the complex-mass scheme, the positive-frequency Wightman function associated with the unstable particle in momentum space reads [42][43][44] Observe that the CMS cut propagator does not quite have the correct cut structure as given by Eq. (7) and, as a result, in principle we cannot connect D +,CMS with physical particles carrying positive energy forward in time. However, at leading orderΓ/m A ∼ λ this is possible; when taking the limitΓ → 0, D +,CMS (p 2 ) turns into a nascent delta function, thereby recovering the standard cut structure for the cut propagator, which allows us to associate with the propagation of positive-energy physical particles. In general, we can write [42][43][44] where the first equality is valid to leading order, and So we see that F (p 2 ,Γ/m A ) corresponds to higher-order contributions. For the phase space as a whole, the function D +,CMS (p 2 ) is suppressed as the imaginary part ofp 0 is small. However, in the region of resonance the small imaginary part yields a non-negligible contribution given by the nascent delta function above. This means that, outside the resonance region, where the CMS cut propagator does not have the correct cut structure (i.e., far from the poles of D +,CMS ), the cut of the CMS propagator of the unstable particle will produce a contribution of higher order in perturbation theory, which can thus be neglected. Only when one is close to the resonance region -which can take place depending on the external momentum configuration of a diagram -that the CMS cut propagator is nonnegligible. These features persist when including corrections to the leading-order result [42][43][44]. In any case, to leading order the cut of the unstable particle propagator is simply the cut through its one-loop correction; in other words, through stable particle propagators. This is a consequence of the fact that at one-loop order the widths in the Complex-Mass Scheme and the traditional on-shell scheme coincide.

When the unitarity method works for unstable particles
As shown in Refs. [40,45], in a theory with unstable particles (of any kind), unitarity is satisfied by the sole inclusion of asymptotically stable states. This means that cuts should not be taken through the unstable particles. Unitaritybased methods represent a kind of generalization of the optical theorem in that they investigate discontinuities of an amplitude in several kinematical channels in order to fully reconstruct loop amplitudes. But if discontinuities of a given loop amplitude are given by the cutting rules, how can one make sense out of the method when one is to cut a propagator associated with an unstable particle?
Let us suppose that all internal propagators of a given one-loop amplitude describe unstable particles. Naively one would conclude that one cannot apply directly unitarity methods to unstable particles. Fortunately, this is not the end of the story. That the method can still be applied to these cases can be observed by recalling the aforementioned discussion of unstable particles within the complex-mass scheme. Indeed, as we have mentioned, at leading order the cut propagator reproduces the nascent delta function typical of stable particles when one is close to the resonance region. Therefore in the complex-mass scheme we write (21) so that, at leading order and close to the resonance region, we find that This implies that the coefficient c 4 of the box integral is still given by Eq. (10) at leading order. In particular, the maximal cut of the one-loop amplitude is also given by the one-loop form of Eq. (2), which represents a proof of this result at leading order to the case of unstable particles running in the loop. Moreover, it is also clear when this procedure cannot be trusted -this is when is off resonance, so that we are not able to put the internal momenta on-shell, see Eq. (19). So when cutting an internal line corresponding to an unstable particle off resonance, the result we obtain is not a contribution to the imaginary part of the scattering amplitude, and, as a consequence, not a valid contribution to the coefficients of the scalar integrals in the above expansion given by Eq. (4). So it is not at all clear whether generalized cuts of propagators associated with unstable particles produce sensible results in this case. We will get back to this off-resonance topic shortly. At higher orders, as the cut of a propagator associated with unstable particles must correspond to the cut through a loop of stable particles when one is close to the resonance region, again we find the correct cut structure. As a result, we believe that Eq. (2) must still be valid to all orders in perturbation theory, but a general proof of this result is beyond the scope of the present work. This is an interesting exploration, and we hope to return to this calculation in the near future.
This discussion shows us that the unitarity method still makes sense in the case of unstable particles running in the loops; in order to implement the technique in a straightforward way, one must ensure that external momentum configurations of an amplitude allows the unstable particle propagator to become resonant. In this case the unstable particle cut propagator will have the correct cut structure to guarantee that unitarity is satisfied. In turn, from previous discussions, we know that the correct strategy for the cut of several propagators is to interpret the corresponding loop integral as a contour integral in C 4 . Moreover, the determination of the contours is such that it must encircle one-particle poles so that we can define the result of integrating over the product of delta functions as given by this contour integral. This is a necessary requirement on the grounds of unitarity -cuts are applied only on the stable particles of the theory so that the sum in Eq. (10) is guaranteed to be over only asymptotically stable states. It is only in this case that one can assert that the result of the integration over the contour |p 2 −m 2 A | = ε that encircle p 2 −m 2 A = 0 will represent an on-shell particle carrying positive energy forward in time. However, it is also clear that when one is off resonance, then one is also away from the pole of the propagator; the associated contour integration must have a vanishing residue in this case. To get a finite result, one must be close to the resonance region; in this case, the operation described in the previous section will yield a well-defined residue.
On the other hand, there is also other situation that the method can be applied without further issues -this is the so-called narrow-width approximation (NWA). In this situation, the coupling to the decay products is taken to be sufficiently small so that only resonance production is significant. In this limit, we can take Here γ = Γm, where Γ is the width of the resonance. In the narrow-width approximation, Γ m and hence that is, near the resonance, we can treat the resonant particle as being on-shell. This means that in this limit the cut taken through the unstable particle with Γ → 0 recovers the result from the cut through the decay products 4 . So effectively the NWA allows us to regard a long-lived resonance as being approximately a stable particle. Moreover, for gauge theories, the NWA does not suffer either from the gauge invariance problem alluded to above [42]. Therefore in this situation the usual reasoning that lies behind the generalized unitarity method can be fully applied.
In other words, for unstable particles the present practice of the unitarity method is valid if the assumption of a resonant unstable propagator is warrant. This can happen depending on external momentum configurations (and this can be proved at least in the Complex-Mass Scheme, as discussed above) or else one should verify whether the narrow-width approximation holds in the particular case under studied. In the complex-mass scheme, for a resonant unstable propagator, one can show that the cut of this propagator follows through the cut of only stable particles, preserving unitarity in Veltman's sense (i.e., by using the Largest Time Equation and employing suitably defined cut propagators). At higher-orders life will not be so simple, but in any case one can still prove that unitarity is satisfied.

Lee-Wick theories
Now let us discuss Lee-Wick-type theories [76][77][78][79][80]. As well known, these class of theories have a ghost mode which can be directly seen from the propagator, which has the form 1 The overall negative sign in the second term signals that this pole is ghost-like. However, the coupling to light particles of the theory makes the heavy ghost state unstable. As discussed above, this implies that generically the associated propagator must be resummed to ensure the validity of perturbation theory. A spectral representation of the corresponding cut propagator can be written as [45] Following the same reasoning employed for normal unstable particles, we find that the cut of internal Lee-Wick propagators of a given one-loop amplitude cannot produce in general a contribution to its discontinuity. However, recall that the structure of a normal resonance propagator is given by with Im[Σ(q)] > 0. Now if the Lagrangian gets modified with a 2 term, the propagator accordingly is changed to Setting Λ → ∞, we obtain the normal resonance. However for large finite Λ we find a heavy-mass resonance, when q 2 ∼ Λ 2 . Near this resonance, the propagator behaves as The residue at this pole is always negative. Furthermore, the sign of the width is always opposite from normal. That is, for both finite m and Λ, we verify the appearance of resonances of both types in the same propagator. In both cases, the imaginary part of the self-energy arising from the coupling to stable states is the same; notwithstanding it manifests itself in distinct ways near the resonances. This means that ghost resonances also obey a unitarity relation as a consequence of the fact that normal resonances satisfy this constraint [45]. The above discussion shows us how to implement unitarity-based methods to one-amplitudes involving unstable ghost modes. In the complex-mass scheme we write the Lee-Wick propagator as On the other hand, as the aforementioned discussion indicates, normal resonances and ghost-like resonances have a similar structure [45] iD(q) ∼ Zi with Z = +1 for a normal resonance and Z = −1 for the ghost resonance. The imaginary part is Z-independent, This implies that the CMS Lee-Wick cut propagator will have the same structure as the propagator associated with a normal unstable particle within the complex-mass scheme. In particular, close to the resonance region, at leading order it will have the form given by Eq. (18), producing the correct cut structure. Hence Eq. (2) is also valid for Lee-Wick theories when one is close to the resonance region. We finally remark that again we can also resort to the NWA in order to apply unitarity techniques to one-loop amplitudes with unstable ghost modes running in the loop. Nevertheless, we emphasize that one must be very careful when dealing with ghost modes in the NWA; in order to reproduce correctly the cuts one must resort to a modification of the contour in performing the loop momentum integration, as originally discussed by Lee and Wick [45].

When the unitarity method does not seem to work for unstable particles
Suppose we wish to study a particular process a + b → c + d which takes place exclusively through loops of unstable particles (of any type) and let us assume that we are off resonance. As asserted above, when one is off resonance the cut through the propagator of an unstable particle will always violate the cut structure. That is, the cut of an unstable particle propagator off resonance yields a contribution of higher order in perturbation theory -such cuts can surely be disregarded. Hence if we use the reasoning above, the cuts of internal unstable propagators will produce a vanishing contribution, resulting in a vanishing amplitude by employing the current practice of the unitarity method to reconstruct it. This is obviously an unsatisfactory answer since we know that amplitudes can be built using standard Feynman rules and Feynman diagrams. So can we take this as an indication that the unitarity method cannot be trusted in this case, as it seems that the associated amplitude (or some of its contributions) could not be determined from the knowledge of its cuts?
There are ways to circumvent this issue. For instance, recall that, in a theory containing unstable particles, unitarity is satisfied by the inclusion of only stable states in unitarity sums. This suggests that, in order to generically implement the generalized unitarity method to a theory containing unstable particles, we must consider the inclusion of only cuts from stable states in unitarity sums. This means that in general one must be able to reformulate the theory in terms of the stable particles only, eliminating from the outset any unstable fields in the Lagrangian. But this will actually introduce non-local vertices in our description. There is one constraint that we should impose in this situation. In order to preserve unitarity, the only acceptable poles in tree-level amplitudes are the ones that come from propagators. Since non-local vertices may generate unphysical poles that would not be be consonant with an exchange of a physical particle, we must impose that such poles have zero residue. Or we must claim that the residues of all such spurious poles must cancel to give zero. Of course other constraints can also be imposed on the non-local vertices, such as proper infrared behavior, valid Ward identities, etc. For a recent interesting discussion of tree-level scattering amplitudes of a particular category of non-local field theories, see Ref. [81]. One-loop unitarity for a class of perturbative scalar quantum field theories with non-local operators of fractional order was established in Ref. [82].

EXAMPLES OF THE USE OF THE UNITARITY METHOD FOR UNSTABLE PARTICLES
We now proceed to discuss with some detail three examples which can be relevant for particle physics in order to see how one can implement the unitarity method when unstable particles run inside loops in scattering amplitudes.

Normal unstable particles
Let us begin our discussions with normal unstable particles. Here we wish to investigate the one-loop helicity amplitude A 1−loop (+ + ++) associated with γ − γ scattering via W loops in the Standard Model [83][84][85][86][87]. It is known that the reaction γγ → γγ via W -boson at one-loop is finite [88]. The standard Feynman-diagram formulation proceeds via box, triangle and bubble diagrams in which one must allow also for unphysical Higgs bosons (when working in a suitable non-linear R ξ gauge) and Faddeev-Popov ghosts in the loops (besides, of course, the W particles). For our study we do not need to consider the unphysical particles; the finiteness of the amplitude will be easily established as we will see. We will use the method of maximal cuts in order to evaluate the amplitude in the expansion in terms of one-loop master integrals. The diagram is depicted in Fig. 1. Since the W boson is heavy, it is unstable and hence in principle we may not be allowed to cut the W internal lines. However, the W mass M is about 80 GeV and its decay width Γ is about 2 GeV [89], hence Γ M and in principle one is justified in resorting to the narrow-width approximation, at least in a primary analysis. In this case, the production and decay of the resonance can be treated approximately in a separate way. As discussed above, the propagator in the NWA has the correct cut structure and hence one can safely use the unitarity method in order to reconstruct the aforementioned box amplitude. On the other hand, as mentioned all one-loop integrals can be written in terms of a basis of scalar one-loop integrals as in Eq. (3), so we are in safe ground here -we can trust the results obtained here. In any case, we can also contemplate our results as an independent check of the helicity amplitude calculated in Ref. [86]; see also Ref. [90].
Using the method of maximal cuts, in the narrow-width approximation the coefficient of the scalar box function may be calculated from the formula (33) The above on-shell tree amplitudes involve one photon and two W bosons. The latter carry explicit SU (2) little-group indices associated with massive spinors. A review of the formalism designed to deal with massive particles can be found in the Appendix. In addition, S is the solution set for the four delta functions of the cut propagators: where we took all external momenta incoming. Notice that the cut conditions imply that Let us calculate the 3-particle amplitude that appears above. Feynman rules will tell us that Let us first choose s = +. Resorting to a bold notation for massive spinors, one finds where we used the Schouten identities and the following relations (which follow from the Schouten identities) [91] r|k 1 |1 k 1 k 2 = M rk 1 1k 2 + rk 2 1k 1 With an almost identical calculation for the s = − case, one finds that For concreteness, let us choose a specific set of helicities for the external photons. Now the cut reads (sum over repeated SU (2) little-group indices is implicit) where we used the on-shell conditions, symmetrization of SU (2) little-group indices and also that |λ I α λ I | β = −M δ α β and λ I λ J = M δ I J . In addition, we have chosen r 1 = p 2 , r 2 = p 1 , r 3 = p 4 and r 4 = p 3 and, as usual, s ij = (p i + p j ) 2 . The convention we use here is the following: We have so far calculated the coefficient associated with the scalar box integral. In order to calculate the coefficients of triangles, bubbles and tadpoles, one must resort to lower-order cuts. For instance, a triple-cut reads The other possible three-particle cut diagrams are obtained from this one by cyclic relabeling of the external particles. The cut conditions are given by Notice that these imply that · p 1 = 0. One finds that where we have chosen r 1 = 2 and r 2 = 1 and we used the cut conditions. As for the two-particle cut, we find As above, the other possible two-particle cut diagrams are obtained from this one by cyclic relabeling of the external particles. The cut conditions are given by One finds that where we used momentum conservation and the cut conditions. Now let us discuss our results. Concerning the triple cut, there are two possible integrals that can contribute, namely the box integral and the triangle integral. However, our result shows the presence of one uncut propagator. So this would exclude triangle integrals from the expansion. To confirm this, let us analyze the 2-particle cut. Again box and triangle integrals contribute, and now also bubble integrals can contribute. Nevertheless, our result shows the presence of two uncut propagators. This confirms the exclusion of triangle integrals from the expansion, and also states the absence of bubble integrals. We can perform a single cut to confirm that there will remain three uncut propagators in the result. Hence, the final answer is that only the box integral is present. So finally we can write where Perm. indicates permutations of external particles, R comprise rational terms and with, as already quoted, Γ M . We can think of the presence of Γ in the above equation as a consequence of the fact that, for unstable particles, we should use a resummed form for its propagator. In any case, one should bear in mind that, as we are considering the decay width to be very small, one must envisage I 4 (p 1 , p 2 , p 3 , p 4 ) in the limit Γ → 0, which ought to be taken at the end of the calculations. Otherwise, one can prove that the coefficient of the box will also display a finite Γ-dependence which is not captured by the unitarity method. As well known, the inclusion of the appropriate dependence on Γ both in the propagators as well as in the corresponding coefficients of the integral is mandatory to ensure the correct gauge cancellations. Nevertheless, as we are considering in this calculation a situation which is dominated by production of on-shell unstable particles with a vanishingly small decay width, finite-width effects are negligible as long as the required precision is taken to be small in comparison with Γ/M 5 .
Rational terms are not detected by unitarity cuts. Hence the above result have potential ambiguities in rational functions. In order to remove such ambiguities one may consider dimensionally regularized representations for the tree amplitudes. This means considering d-dimensional cuts, with d = (4 − 2 ). Photons live in four dimensions, whereas the loop momentum is d-dimensional. In this case one has to be careful when dealing with the summation over states.
A crucial component of generalized unitarity cuts is the sum over physical states. One must be careful in the sum over the physical states of gauge bosons in d dimensions [93]. It is given by the so-called physical state projector: where the ellipsis stand for terms depending on an arbitrary null reference momentum (for massless particles) or on the mass of the particle (for massive particles). In the present case, we will only be concerned with the maximal cut since we already know that only the box integral is present. Here we simply adopt the four-dimensional helicity scheme [94,95] in which all internal and external states (and also polarization vectors) are four-dimensional and loop momentum and phase-space integrals are in d = 4 − 2 dimensions. There are no remaining ambiguities to be considered in our case as the amplitude under consideration vanishes at tree level and there are no ultraviolet divergences.
In general, in the evaluation of the quadruple-cut, we have to discriminate between the dimension of loop momenta and the dimension of the space of physical states; in other words, we should envisage any factor of D emerging from contracting Lorentz indices (δ µ µ = D) as a different quantity in comparison with the dimension d of the loop momenta for which we take d = 4 − 2 . In the limit that D → d one should obtain the same result as before, except that the mass has undergone the shift M 2 → M 2 + µ 2 , where µ α is a vector associated with the (−2 )−dimensional part of the loop momentum. This means that we should use the following modified cut conditions The final result is given by where [48,95,96] and hence We recall that for a complete removal of the ambiguity associated with the rational terms, additional procedures should be carried out [95]; however, such procedures are trivial in the present case since the associated tree-level amplitude vanishes and there are no ultraviolet divergences. Scalar box integrals were explicitly calculated in Refs. [61,65,95]. Furthermore, as promised the helicity amplitude is free from UV divergences. Finally, by exploring the fact that −is 12 s 23 12 23 34 41 2 = 1 and taking into account the different permutations over the external photons, one can easily see that our result agrees perfectly with the ones given in the literature [86,90], apart from an overall phase factor (which is unimportant); one simply needs to be careful with the different conventions on external momenta.

Lee-Wick QED
Now we will discuss a simple example coming from higher-derivative QED. The Lagrangian for the gauge sector reads [79] As well known, Lee-Wick Lagrangians can be rewritten by introducing auxiliary gauge bosons with a very large mass M , very much larger than any other particle masses in our problem. As extensively discussed elsewhere, the coupling of these auxiliary massive gauge bosons to light fields makes them decay, and positive energy is required to excite this resonance [45,97]. Furthermore, this resonance has a "backwards in time" feature in that the propagator has the approximate form (close to the resonance) where we have suppressed Lorentz indices. Notice that there are two minus sign differences from a normal resonance, the −i in the numerator and the −iγ in the denominator. These combined sign differences lead to the distinguishing property of a time-reversed version of a usual unstable particle propagator. This unusual resonance was dubbed a Merlin mode in Refs. [45,98]. There are evidences that point to the stability of theories containing Merlin modes [45,97].
Here we are interested in using the unitarity method to calculate a scattering amplitude involving Merlin particles circulating in the loop. For simplicity, we will consider the narrow-width approximation, Γ M , where Γ is the width of the Merlin particle. The process we have in mind is the muon-electron scattering e + e − → µ + µ − at next-to-leading order, which is one of the simplest in QED processes, but a crucial one in the comprehension of all reactions in e + e − colliders [99]. The calculations with photons running in the loop were carried out in a number of places, see for instance Refs. [100][101][102][103][104][105][106][107][108]. Here we wish to consider solely the one-loop box diagram depicted in Fig. 2. We show how to calculate the coefficient associated with scalar boxes when internal gauge lines are associated with Merlin propagators.
For simplicity, consider the high-energy scattering limit in which the fermions are massless. The 3-particle amplitude involving a fermion, an anti-fermion and a Merlin particle reads (all momenta incoming) where we used that η µν σ αα µ σ ββ ν = 2ε αβ εαβ. Notice that both fermions need to have opposite helicity to give a non-vanishing result. Since p|γ µ |q = q|γ µ |p , one also obtains The 3-point amplitudes involving fermions and a Merlin particle can be computed similarly.
In order to calculate the associated contribution to the coefficient of the scalar box integral, we will resort to the maximal-cut technique, which in the present case means evaluating a quadruple cut. We choose the associated helicities to be h e − = h µ − = −1/2 and h e + = h µ + = +1/2. We have that with the following cut conditions: where p 1 , p 1 (p 2 , p 2 ) are the momenta associated with external µ − and µ + (e − and e + ), respectively. Using the 3-particle amplitudes derived above, we find that As a simple check, one can easily prove that we obtain the same result by resorting to the standard evaluation in terms of polarization sums: where we used the Fierz identity and momentum conservation at each vertex. So we can write that for the contribution coming from Merlin particles running inside the loop, and now Observe the change in the overall sign of the Merlin propagators -as well as in their imaginary parts -in comparison with the W -boson propagators discussed in the previous subsection. Moreover, the presence of Γ can be understood along the same lines as in the normal case -it is important for defining the contour associated with the loop integration but it cannot appear in the final answer obtained after performing the loop integral. Indeed, integrals associated with the Merlin propagators have to be evaluated using the Lee-Wick prescription for integration in the complex 0 plane so that the Wick rotation remains well defined [76][77][78][79][80].

Non-local theories
Previously we have claimed that we can work only with stable particles at the expense of locality. That is, when one is off resonance, a way to deal with the problem of unstable particles is to eliminate them altogether and as a consequence we introduce a non-local description of the problem, albeit one containing only stable modes. Let us briefly discuss this method for the case of light-light scattering. A similar reasoning can also be used for the case of Lee-Wick theories. Since what one has is a non-local interaction, we will discuss the process γγ → γγ within a non-local effective theory. We consider a variant of the theory described in Ref. [109], i.e., we will consider a nonlocal scalar QED. This effective context should account for the issues that one will face when dealing with unstable particles off resonance. The non-local interaction between complex scalars and photons is described by the following Lagrangian density where the non-local coefficient Σ(x − y), which plays the role of a scalar self-energy term, is assumed to be a function of the scalar invariant (x − y) 2 . In addition, the path-ordered exponential U (x, y) is defined as where dω µ is the element of integration along a path connecting points x and y. The path ordering in the definition of U (x, y) is necessary to maintain the gauge-transformation property of U (x, y) for the non-Abelian case. In any case, the path ordering is not required for the photon case [110]. There are also further conditions that should be imposed on the path [109]. The non-local gauge-boson-scalar-scalar vertex can be derived in the standard way and the result is (assuming that the non-local coefficient has analyticity properties resembling standard self-energy functions) where Σ(p 2 ) being the Fourier transform of the self-energy Σ. It is easy to verify that, for the full vertex (containing also the local part, which is not of interest to us here), the Ward identity for dressed scalar propagators is respected [109].
As discussed above, any poles coming from non-local vertices should have zero residues. It is easy to see that this is the case. Furthermore, notice also that This implies the following three-particle amplitudes involving two complex scalars and one photon: q|p|ξ qξ (69) where m is the mass of the scalars. Observe that this is a gauge-invariant amplitude. Since the other gauge interactions are determined from the condition of off-shell gauge invariance, they should not comprise any new on-shell information.
Hence calculation of higher-point tree-level amplitudes may proceed via the usual BCFW recursion relations 6 . The one-loop contribution to the process γγ → γγ proceeding through scalar loops with the above non-local interaction can be calculated in the same way as in the previous case with the W boson. For instance, for the all-plus helicity amplitude one finds that where an extra factor of two was taken into account due to the fact that there is a complex scalar propagating in the loop.

SUMMARY
Here we have discussed the use of unitarity methods in field theories containing resonances. We have shown, through the detailed assessment of three physical situations, how the technique can still be put in practice to such theories at one-loop. Our purpose was to provide an one-loop proof of the validity of unitarity-based methods for unstable particles as a consequence of unitarity itself. That is, generically speaking, if unitarity is satisfied by the inclusion of only stable states in unitarity sums, this implies that, for the unitarity method, one must sum over only the asymptotic states of the theory in Eq. (2). In the complex-mass scheme, this means that, to leading order and close to the resonance region, the cut of the unstable propagator proceeds through the cut of the loop of stable particles. On the other hand, in the NWA, the cut taken through the unstable particle (setting its width to zero) recovers the same result as a cut through the stable decay products. This is the basic requirement for a four-dimensional amplitude involving unstable particles to have cut-constructible parts. Without it, unitarity cuts could not enable one to establish a relationship between the pole structure of the integrand and the branch-cut structure of the loop integral.
Our aim here was not to devise a complete account of all the aspects of the method for unstable particles. Indeed, even though the proof of applicability of unitarity-based methods might be extended to higher orders in perturbation theory, in the present study we have limited ourselves to one-loop order. This is because unitarity cuts provide useful information that can be used in the most efficient way when a complete basis of integrals is known, and this is somewhat straightforward at one loop; the set of master integrals necessary to perform the reduction of generic one-loop tensor integrals is well known, as discussed above. On the other hand, at leading order and close to the resonance region, the CMS cut propagator of the unstable particle turns into a nascent delta function, reproducing the stable particle result, and hence we are allowed to associate the outcome with physical particles carrying positive energy. Such observations allows one to prove Eq. (2) in a straightforward way. In turn, through the examination of simple one-loop examples, we demonstrated explicitly how powerful the method still is when constructing oneloop amplitudes with unstable particles. Further work is recommended to better understand the application of such methods to these theories. One should establish the validity of the method to higher loops: In principle, the proof to higher orders proceeds in much the same way -however, in this case other methods (such as integration-by-parts techniques) should also be employed as integral reduction of the loop integral in terms of a set of master integrals is no longer straightforward and the cuts of internal propagators become more involved.
We have tried to fill this aforementioned gap in the literature of unitarity methods with this primary exploration, and we believe that our study can be useful in the investigations of the Standard Model Effective Field theory (SMEFT) [111], or in the Higgs Effective Field Theory (HEFT) [112]. Indeed, the calculation of loop amplitudes of massive on-shell SMEFT amplitudes focusing on the electroweak sector will obviously involve internal resonances, and an understanding of the unitarity method as a framework to tackle this computation would be most welcome. On the other hand, loop amplitudes of higher-derivative theories also contains resonances, the Merlin modes, and now a careful treatment of those within unitarity methods is available. We believe this will have an important impact on the evaluation of amplitudes of quadratic gravity, a promising conservative ultraviolet completion of quantum gravity. This would indeed be interesting to investigate, and we hope to explore this calculation in subsequent works [113,114].
The gamma matrices in the Weyl basis read where σ µ = (1 σ) andσ µ = (1 − σ). Let us briefly discuss the formalism for massive particles that we used in this work [91,[115][116][117][118]. This is obtained by noting that det p αα = m 2 in the massive case and now p αα has rank 2. So it can be written as the sum of two rank-one matrices: The index I = 1, 2 indicates a doublet of the SU (2) little group. Since det p αα = det λ detλ = m 2 , we simply take det λ = detλ = m. Just like spinor indices, the little group indices are raised and lowered by the SU (2)-invariant tensor ε IJ , ε IJ . It implies that p αα = λ α Iλα I = −λ α Iλα I = |p I p I | pα α = −λ Iα λ α I =λα I λ I α = −|p I p I |.
Comparing this with the usual Dirac equations of motion one is led to the natural identifications for the Dirac spinors: and similarly for the conjugate spinorsū I (p) = −λ α Iλα I v I (p) = λ α Iλα I .
Since the massive spinor bilinears satisfy λ I λ J = mδ I J , λ I λ J = −mε IJ , λ I λ J = mε IJ λ Iλ J = −mδ I J , λ IλJ = mε IJ , λ IλJ = −mε IJ λ I α λ I β = |λ I λ I | = −mδ α β ,λ Iαλβ I = |λ I λ I | = mδβ α it is easy to see that the Dirac spinors obey the usual spin sums. Let us introduce a bold notation to indicate symmetric compositions of the SU (2) little-group indices of massive spinors. One has that p αα = |p p| pα α = −|p p| (86) and the Dirac equation can be rewritten as In addition: We can alternatively write the massive momentum as where k 2 = q 2 = 0 and p 2 = 2k · q = kq qk = m 2 . In terms of bispinors: One has the following identifications The polarization vector of a massive vector boson of momentum p and mass m is given by or, in terms of bispinors: There is an implicit symmetrization on SU (2) indices. These polarizations correspond to transverse and longitudinal modes: