## Linear power corrections to single top production and decay at the LHC in the narrow width approximation

## Sergei Makarov, a Kirill Melnikov, a Paolo Nason, b Melih A. Ozcelik c

- a Institute for Theoretical Particle Physics, KIT, 76128 Karlsruhe, Germany
- b INFN, Sezione di Milano-Bicocca, and Universit` a di Milano-Bicocca, Piazza della Scienza 3, 20126 Milano, Italy
- c Universit´ Paris-Saclay, CNRS, IJCLab, 91405 Orsay, France e

E-mail: sergei.makarov@kit.edu , kirill.melnikov@kit.edu ,

paolo.nason@mib.infn.it , melih.ozcelik@ijclab.in2p3.fr

Abstract: We consider top quark production and decay in the narrow width approximation and study if the polarisation effects, that manifest themselves in correlations of angular distributions of particles from top quark decays and final state jets in the production subprocess, are affected by linear power corrections. We find that, in general, the answer to this question is affirmative. We also discuss how these non-perturbative corrections affect polarisation observables used to study single top production at the LHC. Finally, we point out that generic kinematic distributions of leptons from top quark decays are affected by linear power corrections, which may have implications for proposals to extract the top quark mass from such leptonic observables. On the other hand, we demonstrate that the distribution of the 'out-of-collision-plane' component of the positron momentum is free from linear power corrections, making it an interesting candidate for the top quark mass measurement.

## Contents

| 1   | Introduction                                                          |   1 |
|-----|-----------------------------------------------------------------------|-----|
| 2   | The narrow width approximation                                        |   2 |
| 3   | QCD corrections to the production sub-process                         |   5 |
| 3.1 | Real emission contributions                                           |   6 |
| 3.2 | Virtual corrections and renormalisation in the production sub-process |  10 |
| 4   | Corrections to the decay sub-process                                  |  12 |
| 5   | Redefinition of the top quark mass parameter                          |  14 |
| 6   | Results and corrections to observables                                |  17 |
| 7   | Conclusions                                                           |  24 |
| A   | Alternative derivation of power corrections                           |  24 |

## 1 Introduction

In a recent series of papers [1, 2] we have studied O (Λ QCD ) power corrections to top quark production in hadron collisions using the approach based on infra-red renormalons. 1 In this paper we extend these analyses by accounting for top quark decays. Such an extension is non-trivial. Indeed, since the top quark width Γ t serves as an infra-red regulator, its interplay with the non-perturbative QCD parameter Λ QCD and its proxy in the context of renormalon calculus - the gluon mass λ - is important. Since Γ t ≫ Λ QCD ∼ λ , top quarks decay before hadronisation. A renormalon-based analysis of power corrections in such a case is technically very challenging, because the produced top quarks are off-shell, and diagrams where real or virtual gluons connect production and decay stages of off-shell top quarks have to be considered. Although this problem represents an interesting challenge for future work, we believe that it makes sense to start by considering the opposite case Γ t ≪ Λ QCD ∼ λ , which can be studied in the narrow width approximation. Even if the phenomenological relevance of such an analysis is limited, it provides, for the very first time, an estimate of the non-perturbative corrections to a full physical process with unstable particles at a hadron collider and, as such, might be quite valuable for modelling the non-perturbative effects.

1 An in-depth discussion of infra-red renormalons in QCD can be found in review [3]. For a detailed description of how these methods can be used in the context of hadron collider applications, see ref. [4].

The narrow width approximation leads to important technical simplifications since, in this case, QCD radiative corrections cannot connect top quark production and top quark decay sub-processes [5, 6]. In fact, the QCD corrections to top quark production and decay process are known through next-to-next-to-leading order in this approximation [7, 8]. However, even when the narrow width approximation is used, the production and decay sub-processes are not independent; the communication between them occurs because of momentum conservation and polarisation effects . 2 For the case of single top production, that we study in this paper, this implies that e.g. the direction of the outgoing light-quark jet and the direction of the positron in the top quark decay are not independent in spite of the fact that they originate from widely separated sub-processes. As we will show below, linear power corrections affect angular distributions in single top production. In fact, we find that they vanish when positrons from top decays and light jets from the production process are collinear to each other, but are non-trivial functions of momenta of final state particles otherwise. We will also show that the O (Λ QCD ) effects that we discuss in this paper impact the various observables designed to study top quark polarisation [14-16] differently, so that for each of them a dedicated study is required.

The paper is organised as follows. In the next section we describe the narrow width approximation and present the result for the differential cross section of single top production followed by top decay in a way that is useful for the subsequent analysis. In sections 3 and 4 we discuss the calculation of O ( λ ) corrections to the production and decay sub-processes. In section 5, O ( λ ) corrections due to mass-parameter redefinition are studied. In section 6, corrections to observables are discussed and the final formulas are derived. We conclude in section 7. In the appendix we illustrate an alternative approach to the calculation of power corrections where, at variance with the method used in the previous publications [1, 2], we deal directly with the amplitude of the process rather than with its square.

## 2 The narrow width approximation

In this paper we consider the following partonic process

$$u ( p _ { u } ) + b ( p _ { i } ) \to d ( p _ { d } ) + t ( p _ { t } ) \, & \quad \ > b ( p _ { f } ) + \nu ( p _ { 1 } ) + e ^ { + } ( p _ { 2 } ).$$

It is shown in Fig. 1. The amplitude for this process can be written as

$$\mathcal { M } = A _ { \text{dec} } ^ { i } \frac { i \delta ^ { i j } ( \vec { p } _ { t } + m _ { t } ) } { p _ { t } ^ { 2 } - m _ { t } ^ { 2 } + i m _ { t } \Gamma _ { t } } A _ { \text{prod} } ^ { j },$$

where i, j are colour indices. In the on-shell p 2 t → m 2 t limit, the quantities A dec and A prod correspond to on-shell amplitudes for the respective sub-processes from which the top quark spinors are removed.

2 Theoretical studies of spin correlations in top quark pair production and polarisation effects in single top production have a long history, see e.g. refs [9-13] for original and more recent work.

Figure 1 : The Born amplitude for single top production and decay.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00632v2-figure-44876f83-5a39-43a2-b86f-37d7cfda5202.jpg)

The expression for the cross section becomes

$$\mathrm d \sigma _ { P D } = \frac { 1 } { \mathcal { N } } \frac { \mathrm d \Phi ( p _ { u }, p _ { i } ; p _ { d }, \{ p _ { \text{dec} } \} ) } { ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) ^ { 2 } + m _ { t } ^ { 2 } \Gamma _ { t } ^ { 2 } } \sum _ { \text{spins} } A _ { \text{dec} } ^ { i } ( \mathfrak { p } _ { t } + m _ { t } ) A _ { \text{prod} } ^ { i } \bar { A } _ { \text{prod} } ^ { j } ( \mathfrak { p } _ { t } + m _ { t } ) \bar { A } _ { \text{dec} } ^ { j }, \quad ( 2. 3 )$$

where the suffix PD indicates production and decay, dΦ is the phase space of the process in eq. (2.1), the sum over spins includes external particles only and { p dec } describes the momenta p f, 1 2 , that refer to particles originating from the 'decay' of the virtual top quark. The normalisation factor N includes all the averaging factors needed to compute the cross section of the process in eq. (2.1) as well as the relevant flux factor.

To simplify eq. (2.3), we note that since in the narrow width approximation no colour transfer between production and decay amplitudes occurs, the following equation holds

$$A _ { \text{dec} } ^ { i } \dots \bar { A } _ { \text{dec} } ^ { j } = \frac { \delta ^ { i j } } { N _ { c } } A _ { \text{dec} } ^ { k } \dots \bar { A } _ { \text{dec} } ^ { k },$$

where N c = 3 is the number of colours. Hence,

$$\mathrm d \sigma _ { P D } = \frac { 1 } { \mathcal { N } \, N _ { c } } \frac { \mathrm d \Phi ( p _ { u }, p _ { i } ; p _ { d }, \{ p _ { \text{dec} } \} ) } { ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) ^ { 2 } + m _ { t } ^ { 2 } \Gamma _ { t } ^ { 2 } } \, \sum _ { \text{spins} } A _ { \text{dec} } ^ { i } ( \varphi _ { t } + m _ { t } ) A _ { \text{prod} } ^ { j } \bar { A } _ { \text{prod} } ^ { j } ( \varphi _ { t } + m _ { t } ) \bar { A } _ { \text{dec} } ^ { i }. \quad ( 2. 5 )$$

To proceed further, we factorise the phase space

$$\mathrm d \Phi ( p _ { u }, p _ { i } ; p _ { d }, \{ p _ { \text{dec} } \} ) = \frac { \mathrm d p _ { t } ^ { 2 } } { 2 \pi } \, \mathrm d \Phi ( p _ { u }, p _ { i } ; p _ { d }, p _ { t } ) \, \mathrm d \Phi ( p _ { t } ; \{ p _ { \text{dec} } \} ),$$

and make use of the fact that we work in the narrow width approximation which implies that the following equation holds

$$\frac { 1 } { ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) ^ { 2 } + m _ { t } ^ { 2 } \Gamma _ { t } ^ { 2 } } \Big | _ { \Gamma _ { t } / m _ { t } \to 0 } = \frac { 2 \pi } { 2 m _ { t } \Gamma _ { t } } \delta ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ).$$

We obtain

$$\mathrm d \sigma _ { P D } = \frac { \mathrm d \Phi ( p _ { u }, p _ { i } ; p _ { d }, p _ { t } ) \, \mathrm d \Phi ( p _ { t } ; \{ p _ { \text{dec} } \} ) } { \mathcal { N } \, N _ { c } \, 2 m _ { t } \Gamma _ { t } } \, \sum _ { \text{spins} } A _ { \text{dec} } ^ { i } ( \mathfrak { p } _ { t } + m _ { t } ) A _ { \text{prod} } ^ { j } \bar { A } _ { \text{prod} } ^ { j } ( \mathfrak { p } _ { t } + m _ { t } ) \bar { A } _ { \text{dec} } ^ { i }. \ ( 2. 8 )$$

It is easy to see that the following (matrix) equation holds 3

$$\sum _ { \text{spins} } ( \vec { p } _ { t } + m _ { t } ) \bar { A } _ { \text{dec} } ^ { i } A _ { \text{dec} } ^ { i } ( \vec { p } _ { t } + m _ { t } ) = X \, \vec { \phi } _ { t, D },$$

where the sum extends over the polarisations of the top quark decay products and where we have defined the top quark spin density matrix

$$\phi _ { t, D } = ( \phi _ { t } + m _ { t } ) \frac { \left ( 1 + \gamma _ { 5 } \phi _ { D } \right ) } { 2 }.$$

In eq. (2.9), X is a function of scalar products constructed out of the top quark momentum p t and the momenta of its decay products, and s µ D is a space-like unit vector that also depends on p t and the momenta of the particles in top decay. To find X , we compute the trace of both sides of eq. (2.9) and obtain

$$\sum _ { \text{spins} } 2 m _ { t } \text{Tr} \left [ ( \phi _ { t } + m _ { t } ) \bar { A } _ { \text{dec} } ^ { i } A _ { \text{dec} } ^ { i } \right ] = 2 m _ { t } \text{ } X.$$

We can write this equation in a better way by using the formula for the differential decay width of the unpolarised top quark

$$d \Gamma _ { t } = \frac { d \Phi ( p _ { t } ; \{ p _ { \text{dec} } \} ) } { 4 m _ { t } N _ { c } } \, \sum _ { \text{spins} } \text{Tr} \left [ ( \vec { p } _ { t } + m _ { t } ) \bar { A } _ { \text{dec} } ^ { i } A _ { \text{dec} } ^ { i } \right ].$$

It follows that

$$X = \sum _ { \text{spins} } \text{Tr} \left [ ( \vec { p } _ { t } + m _ { t } ) \bar { A } _ { \text{dec} } ^ { i } A _ { \text{dec} } ^ { i } \right ] = 4 m _ { t } N _ { c } \, \frac { \text{d} \Gamma _ { t } } { \text{d} \Phi ( p _ { t } ; \{ p _ { \text{dec} } \} ) }.$$

Finally, using this result as well as eq. (2.9), we re-write eq. (2.8) as

$$d \sigma _ { P D } = 2 \frac { d \Gamma _ { t } } { \Gamma _ { t } } \times d \sigma _ { t } ( s _ { D } ).$$

In the above equation, d σ t ( s D ) is the cross section for producing a single top quark whose spin is aligned with the axis s D ; it is defined as

$$\mathrm d \sigma _ { t } ( s _ { D } ) = \frac { \mathrm d \Phi ( p _ { u }, p _ { i } ; p _ { d }, p _ { t } ) } { \mathcal { N } } \, \sum _ { \text{spins} } \text{Tr} \left [ \phi _ { t, D } A _ { \text{prod} } ^ { j } \bar { A } _ { \text{prod} } ^ { j } \right ].$$

We note that we will refer to the spin axis s D as the top quark spin vector below. It is important to emphasise that for the process in eq. (2.1) this vector depends on the momenta of the top quark decay products. It is computed in the appendix and given in eq. (A.15). The calculation of the top quark polarised cross section in eq. (2.14) proceeds in the standard way. The only difference with the unpolarised case is that instead of the density matrix ( / p t + m t ) one has to use / ρ t,D .

3 This is due to the fact that helicities of all other external particles are fixed, so that the top quark must be in a pure spin state. See appendix for more details.

The formula for the cross section shown in eq. (2.14) is suitable for computing QCD corrections to the production process followed by tree-level decay. When tree-level production is followed by the QCD-corrected decay, it is more convenient to deal with the decay of the polarised top quark. Following the discussion above, we find

$$d \sigma _ { P D } = 2 \frac { d \Gamma _ { t } ( s _ { P } ) } { \Gamma _ { t } } \times d \sigma _ { t },$$

where this time d σ t is the unpolarised top quark production cross section and dΓ ( t s P ) is the decay rate of a polarised top quark

$$d \Gamma _ { t } ( s _ { P } ) = \frac { d \Phi ( p _ { t } ; \{ p _ { \text{dec} } \} ) } { 4 m _ { t } N _ { c } } \, \sum _ { \text{spins} } \text{Tr} \left [ \phi _ { t, P } \bar { A } _ { \text{dec} } ^ { i } A _ { \text{dec} } ^ { i } \right ].$$

The corresponding spin density matrix reads

$$\phi _ { t, P } = ( \phi _ { t } + m _ { t } ) \frac { ( 1 + \gamma _ { 5 } \oint _ { P } ) } { 2 }.$$

The polarisation vector of the top quark s P is given in eq. (A.16). We note that at leading order the two representations of the full differential cross section, given in eqs (2.14) and (2.16), are equivalent, so that the following equation holds

$$d \sigma _ { t } \, d \Gamma _ { t } ( s _ { P } ) = d \sigma _ { t } ( s _ { D } ) \, d \Gamma _ { t }.$$

This equation allows us to use one representation to study radiative corrections to the production process and another one to study radiative corrections to top decay. Furthermore, eq. (2.19) will be useful for studying effects related to the redefinition of the top quark mass. Such a redefinition is needed to remove O ( λ ) corrections caused by the fact that the pole mass of the top quark is used in perturbative computations and that this mass parameter itself receives O ( λ ) corrections when expressed through a short-distance mass [17, 18].

Eq. (2.19) will also be useful for simplifying the final result. In particular, for single top production, the above equation assumes a particularly simple form 4

$$d \sigma _ { t } \, \text{d} \Gamma _ { t } ( s _ { P } ) = \text{d} \sigma _ { t } ( s _ { D } ) \, \text{d} \Gamma _ { t } = \frac { 1 } { 2 } \text{d} \sigma _ { t } \text{d} \Gamma _ { t } \left ( 1 - s _ { D } \cdot s _ { P } \right ).$$

Hence, by choosing the explicit form of one of the two spin vectors and pretending that the other one is general, one can put the polarisation information either to the production or to the decay sub-process of the full process in eq. (2.1).

## 3 QCD corrections to the production sub-process

We consider the production sub-process

$$u ( p _ { u } ) + b ( p _ { i } ) \to d ( p _ { d } ) + t ( p _ { t } ),$$

4 For the derivation of this formula, see eq. (A.10) and the discussion before it.

with the assumption that the top quark is polarised. As explained in ref. [19], corrections to the light-quark line do not produce linear O ( λ ) contributions. For this reason, we only need to consider the QCD corrections to the heavy-quark line. To describe this process, we can use eq. (2.14) and compute the standard perturbative contributions to the polarised cross section. In our discussion, we will assume that the reader is familiar with the computation for the unpolarised case reported in ref. [1], and we will mostly emphasise the differences between the two cases in what follows.

## 3.1 Real emission contributions

We begin with the real emission corrections to the process in eq. (3.1)

$$u ( p _ { u } ) + b ( p _ { i } ) \to d ( q _ { d } ) + t ( q _ { t } ) + g ( k ).$$

The gluon g k ( ) is massive, k 2 = λ 2 . We remark that we have used slightly different notations for the four-momenta of the top quark and the down quark in eqs (3.1) and (3.2). This is done for future convenience since, as we will see, momenta q t and q d will absorb the recoil due to the emitted soft gluon. We will also consider the top quark to be polarised, and we will denote its spin vector with s µ D .

The calculation proceeds along the lines described in ref. [1] where the case of stable top quark was discussed. The main element of that discussion was the Low-Burnett-Kroll theorem [20, 21] that can be used to describe soft radiation in QCD with next-to-leadingpower accuracy. 5 This theorem was derived from the transversality of the amplitude with respect to the gluon momentum which allowed us to relate soft gluon emission by the external particles and the structure-dependent radiation. The fact that the top quark spinor represents a state with a particular polarisation plays no role in this argument. Hence, upon writing

$$\mathcal { A } _ { \text{prod} } & = g _ { s } T _ { i j } ^ { a } \epsilon _ { \mu } \mathcal { M } ^ { \mu },$$

and repeating all the steps described in ref. [1], we arrive at the following result

$$\cdots \cdots _ { r } \cdots \cdots _ { r - 1 } \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ \mathcal { M } ^ { \mu } & = J _ { t } ^ { \mu } \bar { u } _ { t } \mathfrak { N } ( q _ { t } + k, p _ { i }, q _ { d },.. ) u _ { i } + J _ { i } ^ { \mu } \bar { u } _ { t } \mathfrak { N } ( q _ { t }, p _ { i } - k, q _ { d },.. ) u _ { i } \\ & \quad + \bar { u } _ { t } \left [ \mathfrak { S } _ { t } ^ { \mu } \mathfrak { N } ( q _ { t }, p _ { i }, q _ { d },.. ) + \mathfrak { N } ( q _ { t }, p _ { i }, q _ { d },.. ) \mathfrak { S } _ { i } ^ { \mu } \right ] u _ { i } \\ & \quad - \bar { u } _ { t } \left [ \frac { \partial \mathfrak { N } ( q _ { t }, p _ { i }, q _ { d },.. ) } { \partial q _ { t, \mu } } + \frac { \partial \mathfrak { N } ( q _ { t }, p _ { i }, q _ { d },.. ) } { \partial p _ { i, \mu } } \right ] u _ { i }. \\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \end{$$

In the above expression, N is the tree-level amplitude for single top production from which top and bottom spinors have been removed, and

$$\text{mno move to the movie}, \, \text{and} \\ J _ { t } ^ { \mu } & = \frac { 2 q _ { t } ^ { \mu } + k ^ { \mu } } { d _ { t } }, \quad \text{S} _ { t } ^ { \mu } = \frac { \sigma ^ { \mu \nu } k _ { \nu } } { d _ { t } }, \\ J _ { i } ^ { \mu } & = \frac { 2 p _ { i } ^ { \mu } - k ^ { \mu } } { d _ { i } }, \quad \text{S} _ { i } ^ { \mu } = \frac { \sigma ^ { \mu \nu } k _ { \nu } } { d _ { i } }, \\ \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots$$

are top (bottom) currents and spin operators, respectively. The two quantities in the above equation, d t = ( q t + ) k 2 -m 2 t and d i = ( p i -k ) 2 , are the denominators of top and bottom propagators. Furthermore, σ µν = [ γ µ , γ ν ] / 2.

5 An extension of this theorem to processes with polarised particles can be found in ref. [22].

We can further simplify the expression in eq. (3.4) by combining the first two terms, expanded to first subleading power in k , with the last two terms. We obtain

$$\mathcal { M } ^ { \mu } & = J ^ { \mu } \bar { u } _ { t } \mathbf N ( q _ { t }, p _ { i }, q _ { d },.. ) u _ { i } + \bar { u } _ { t } ( L ^ { \mu } \mathbf N ( q _ { t }, p _ { i }, q _ { d },.. ) ) u _ { i } \\ & + \bar { u } _ { t } \left [ \mathbf S _ { t } ^ { \mu } \mathbf N ( q _ { t }, p _ { i }, q _ { d },.. ) + \mathbf N ( q _ { t }, p _ { i }, q _ { d },.. ) \mathbf S _ { i } ^ { \mu } \right ] u _ { i },$$

where we have introduced the notation

$$J ^ { \mu } = J ^ { \mu } _ { t } + J ^ { \mu } _ { i }, \ \ L ^ { \mu } = L ^ { \mu } _ { t } - L ^ { \mu } _ { i },$$

with

$$L ^ { \mu } _ { t } = \frac { 2 } { d _ { t } } \left ( p ^ { \mu } _ { t } k ^ { \nu } \frac { \partial } { \partial q ^ { \nu } _ { t } } - ( p _ { t } k ) \frac { \partial } { \partial q _ { t, \mu } } \right ), \quad L ^ { \mu } _ { i } = \frac { 2 } { d _ { i } } \left ( p ^ { \mu } _ { i } k ^ { \nu } \frac { \partial } { \partial p ^ { \nu } _ { i } } - ( p _ { i } k ) \frac { \partial } { \partial p _ { i, \mu } } \right ).$$

Eq. (3.6) expresses the amplitude with the emission of a soft gluon through the elastic amplitude and its derivatives. However, we observe further simplifications if we square the amplitude and sum over the polarisations of the external particles, or consider external states with definite helicities. We find

$$& | \mathcal { M } | ^ { 2 } = - g _ { \mu \nu } \mathcal { M } ^ { \mu } \mathcal { M } ^ { \nu, + } = - J ^ { \mu } J _ { \mu } \text{Tr} \left [ \phi _ { t, D } \mathbb { N } \phi _ { i } \bar { \mathbb { N } } \right ] \\ & - J ^ { \mu } \text{Tr} \left [ \phi _ { t, D } \mathbb { N } \phi _ { i } L _ { \mu } \bar { \mathbb { N } } \right ] - J _ { \mu } \text{Tr} \left [ \phi _ { t, D } ( L ^ { \mu } \mathbb { N } ) \phi _ { i } \bar { \mathbb { N } } \right ] \\ & + J _ { \mu } \text{Tr} \left [ [ \mathbb { S } _ { t } ^ { \mu }, \phi _ { t, D } ] \mathbb { N } \phi _ { i } \bar { \mathbb { N } } \right ] + J _ { \mu } \text{Tr} \left [ \phi _ { t, D } \mathbb { N } [ \phi _ { i }, \mathbb { S } _ { t } ^ { \mu } ] \bar { \mathbb { N } } \right ].$$

Note that the only difference between eq. (3.9) and a similar expression for |M| 2 in ref. [1] is that the density matrix / ρ t,D defined in eq. (2.10) appears in eq. (3.9) instead of ( / p t + m t ). Since,

$$[ \mathbf S _ { t } ^ { \mu }, \phi _ { t, D } ] = ( - L _ { t } ^ { \mu } - S _ { t } ^ { \mu } ) \, \phi _ { t, D } = - ( L ^ { \mu } + S _ { t } ^ { \mu } ) \phi _ { t, D }, \quad [ \mathbf S _ { i } ^ { \mu }, \phi _ { i } ] = - L _ { i } ^ { \mu } \phi _ { i } = ( L ^ { \mu } + S _ { t } ^ { \mu } ) \phi _ { i }, \ ( 3. 1 0 )$$

where

$$S _ { t } ^ { \mu } = \frac { 2 } { d _ { t } } \left ( s _ { D } ^ { \mu } k ^ { \nu } \frac { \partial } { \partial s _ { D } ^ { \nu } } - ( s _ { D } k ) \frac { \partial } { \partial s _ { D, \mu } } \right ),$$

we obtain

$$& | \mathcal { M } | ^ { 2 } = - J ^ { \mu } J _ { \mu } \text{Tr} \left [ \phi _ { t, D } \mathbf N \phi _ { i } \bar { \mathbf N } \right ] - J ^ { \mu } \text{Tr} \left [ \phi _ { t, D } \mathbf N \phi _ { i } L _ { \mu } \bar { \mathbf N } \right ] - J ^ { \mu } \text{Tr} \left [ \phi _ { t, D } ( L _ { \mu } \mathbf N ) \phi _ { i } \bar { \mathbf N } \right ] \\ & - J ^ { \mu } \text{Tr} \left [ ( ( L _ { \mu } + S _ { t, \mu } ) \phi _ { t, D } ) \mathbf N \phi _ { i } \bar { \mathbf N } \right ] - J ^ { \mu } \text{Tr} \left [ \phi _ { t, D } \mathbf N ( ( L _ { \mu } + S _ { t, \mu } ) \phi _ { i } ) \bar { \mathbf N } \right ].$$

We emphasise that in the above formulas, starting from eq. (3.10), derivatives with respect to the four-momenta of partons that appear in the operators L t/i do not act on the polarisation vector s µ D .

Making use of the fact that L µ is a linear differential operator, and that the only dependence on s µ D is in the density matrix / ρ t,D , we combine the last four terms to obtain a derivative of the leading order polarised amplitude. The final result reads

$$| \mathcal { M } | ^ { 2 } = - \left ( J ^ { \mu } J _ { \mu } + J _ { \mu } ( L ^ { \mu } + S _ { t } ^ { \mu } ) \right ) F _ { p } ( p _ { u }, p _ { i }, q _ { t }, q _ { d }, s _ { D } ),$$

where

$$F _ { p } ( p _ { u }, p _ { i }, q _ { t }, q _ { d }, s _ { D } ) = \text{Tr} \left [ \phi _ { t, D } \mathbf N \phi _ { i } \bar { \mathbf N } \right ],$$

is the matrix element squared for the production process with polarised top quark and where, by construction, derivatives with respect to momenta do not act on / s D . We also note that one can obtain the unpolarised result by simply setting s D → 0 in eq. (3.14) and multiplying the result by a factor 2. By a slight abuse of notation, we will write

$$2 F _ { p } ( p _ { u },.., q _ { d }, s _ { D } \to 0 ) \equiv F _ { p } ( p _ { u },.., q _ { d } ),$$

in what follows.

Similarly to the stable-top case, the first term on the right hand side in eq. (3.13) contributes to the top quark production cross section starting at order O ( λ 0 ), whereas the second and the third ones contribute at order O ( λ ). To extract O ( λ ) contributions from the first term in eq. (3.13), we redefine the momenta of various particles to remove the momentum k from the energy-momentum conserving delta-function. However, an important difference between stable and unstable top cases occurs at this point because the top quark momentum appears both in the production and in the decay phase spaces. Hence, redefinition of the top quark momentum in the production process leads to the redefinition of the top quark momentum in the decay and one needs to understand how to deal with it.

To adhere as much as possible to the discussion of single top production and top decay that were studied separately in ref. [1], we employ the momenta redefinitions used there. We write

$$q _ { t } = p _ { t } - k + \frac { ( p _ { t } k ) } { ( p _ { t } p _ { d } ) } p _ { d }, \quad q _ { d } = p _ { d } - \frac { ( p _ { t } k ) } { ( p _ { t } p _ { d } ) } p _ { d }.$$

Repeating the steps described in ref. [1], we find

$$\mathrm d \Phi _ { P } ( p _ { u }, p _ { i } ; q _ { d }, q _ { t }, k ) = \mathrm d \Phi _ { P } ( p _ { u }, p _ { i } ; p _ { d }, p _ { t } ) \ [ \mathrm d k ] _ { \lambda } \ \left ( 1 + \frac { ( k p _ { d } ) } { ( p _ { t } p _ { d } ) } - \frac { ( k p _ { t } ) } { ( p _ { t } p _ { d } ) } + \mathcal { O } ( \lambda ^ { 2 } ) \right ), \ \ ( 3. 1 7 )$$

where

$$[ d k ] _ { \lambda } = \frac { d ^ { 4 } k } { ( 2 \pi ) ^ { 3 } } \delta _ { + } ( k ^ { 2 } - \lambda ^ { 2 } ).$$

In the stable-top case, once the above transformation is performed and relevant matrix elements squared are expanded in k , integration over k becomes possible. The same happens when decay is considered except that we need to account for the change in the differential decay width and the decay phase space introduced by momenta redefinitions in eq. (3.16).

Momenta redefinitions are only relevant for the leading O ( λ -2 ) term in eq. (3.13). Its contribution to the cross section is proportional to

$$\mathrm d \Phi _ { P } ( \{ q _ { \ln } \} ; p _ { d }, p _ { t } ) \ \mathrm d \Gamma _ { t } ( q _ { t }, \{ q _ { d e c } \} ) \ J _ { \mu } ^ { ( 0 ) } J ^ { ( 0 ), \mu } \ F _ { p } (.., q _ { t }, q _ { d }, s _ { D } ( q _ { t }, q _ { 2 } ) ) | _ { q _ { t } \rightarrow p _ { t } +.., q _ { d } \rightarrow p _ { d } +.. },$$

where dΦ P is the phase space of the production subprocess after momenta redefinition and expansion in k , and J (0) µ is the leading power contribution to the eikonal current given

in eq. (3.7). The required momenta shifts are shown in eq. (3.16). We note that the differential width that appears in the above expression depends on the original top quark momentum q t . Furthermore, momenta redefinitions also affect the spin vector of the top quark s D as it depends on the top quark momentum q t .

To understand how to expand eq. (3.19) to first sub-leading order in the gluon momentum, we note that the top quark momentum redefinition in eq. (3.16) can be interpreted as a Lorentz transformation. Indeed, we can write eq. (3.16) as follows

$$q _ { t } ^ { \mu } & = \Lambda ^ { \mu \nu } p _ { t, \nu },$$

where

$$\Lambda ^ { \mu \nu } = g ^ { \mu \nu } - \frac { k ^ { \mu } p _ { d } ^ { \nu } - p _ { d } ^ { \mu } k ^ { \nu } } { ( p _ { t } p _ { d } ) } = g ^ { \mu \nu } + \delta \Lambda ^ { \mu \nu }.$$

Since δ Λ µν is an anti-symmetric traceless matrix, it can be interpreted as an infinitesimal Lorentz transformation.

This observation has important implications for the calculation of the differential decay width dΓ ( t q , t { q dec } ) that appears in eq. (3.19). Indeed, after momentum transformation, it becomes

$$\mathrm d \Gamma _ { t } ( q _ { t }, \{ q _ { \text{dec} } \} ) = \mathrm d \Gamma _ { t } ( \Lambda p _ { t }, \{ q _ { \text{dec} } \} ) = \mathrm d \Gamma _ { t } ( \Lambda p _ { t }, \{ \Lambda p _ { \text{dec} } \} ) = \mathrm d \Gamma _ { t } ( p _ { t }, \{ p _ { \text{dec} } \} ), \quad ( 3. 2 2 )$$

where we also transformed momenta of the final state particles in the decay and made use of the fact that the differential width is invariant under Lorentz transformations.

Although the above result implies that momenta redefinitions in the production do not generate O ( k ) corrections in the unpolarised decay width that appears in eq. (3.19), the need to redefine momenta of final-state particles in the decay has implications for the spin vector s D . Indeed, s D depends on q t and q 2 and, therefore, changes when the above momenta transformations are performed. It is easy to see that this change is described by a Lorentz boost

$$s _ { D } ^ { \mu } ( q _ { t }, q _ { 2 } ) = \Lambda ^ { \mu \nu } s _ { D, \nu } ( p _ { t }, p _ { 2 } ),$$

where Λ µν is given in eq. (3.21).

We are now in a position to write the result for the O ( λ ) contribution to the differential cross section of the process in eq. (2.1) due to an emission of a soft gluon in the production sub-process. It arises as a sum of the correction to the production matrix element described in eq. (3.13), correction to the production phase-space shown in eq. (3.17) and corrections to the leading order term shown in eq. (3.19) and discussed afterwards. Many of these contributions do not involve modifications of the spin vector s D and are identical to contributions studied in ref. [1]. New spin-dependent contributions arise because of the spin-operator S t in eq. (3.13) and because of the Lorentz boost of the spin-vector in eq. (3.23) that has to be inserted into the function F p in eq. (3.19) and expanded in k .

Hence, we write

$$\mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { P D } ^ { R, \text{prod} } \right ] = 2 \frac { \mathrm d \Gamma _ { t } } { \Gamma _ { t } } \, \mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { t } ^ { R } ( s _ { D } ) \right ] _ { k } + \mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { P D, \text{ new} } ^ { R, \text{prod} } \right ] _ { k },$$

where subscript k indicates that the integration over gluon momentum is still to be performed. Furthermore, the first term on the right hand side in eq. (3.24) is the same as in the no-decay case [1] except that one employs the polarised cross section for single top production, and the second term is new. It reads

$$\mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { P D, \text{new} } ^ { R, \text{prod} } \right ] _ { k } = - \frac { \mathrm d \Phi _ { P } ( \dots ; p _ { t }, p _ { d } ) } { \mathcal { N } } [ \mathrm d k ] _ { \lambda } \, 2 \, \frac { \mathrm d \Gamma _ { t } } { \Gamma _ { t } } \, J _ { \alpha } ^ { ( 0 ) } \, \mathcal { O } ^ { ( s ), \alpha } \, F _ { p } ( \dots, p _ { t }, p _ { d }, s _ { D } ( p _ { t }, p _ { 2 } ) ), \ \ ( 3. 2 5 )$$

where

$$\mathcal { O } ^ { ( s ), \alpha } = S _ { t } ^ { \alpha } - J ^ { ( 0 ), \alpha } s _ { D, \mu } \, \delta \Lambda ^ { \mu \nu } \frac { \partial } { \partial s _ { D } ^ { \nu } }.$$

To complete the computation of the real-emission contributions we need to integrate over the gluon momentum. This is straightforward since in both old and new contributions in eq. (3.24) the dependence on the gluon momentum k is exposed and one can use the formulas in appendix A of ref. [1] to calculate the relevant integrals over k . Hence, we find

$$\mathcal { T } _ { \lambda } \left [ d \sigma _ { P D } ^ { R, \text{prod} } \right ] = 2 \frac { \text{d} \Gamma _ { t } } { \Gamma _ { t } } \ \mathcal { T } _ { \lambda } \left [ d \sigma _ { t } ^ { R } ( s _ { D } ) \right ] + \mathcal { T } _ { \lambda } \left [ d \sigma _ { P D, \text{new} } ^ { R, \text{prod} } \right ].$$

The first term on the right hand side of eq. (3.27) can be found in eq. (2.31) of ref. [1], where F LO should be replaced with F p ( ..., s D ), and the second (new) term reads

$$\mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { P D, \text{ new} } ^ { R, \text{prod} } \right ] = \frac { \alpha _ { s } C _ { F } } { 2 \pi } \frac { \pi \lambda } { m _ { t } } \, 2 \, \frac { \mathrm d \Gamma _ { t } } { \Gamma _ { t } } \, s _ { D } ^ { \mu } \, W _ { \mu \nu } \, \frac { \partial } { \partial s _ { D, \nu } } \, \mathrm d \sigma _ { t } ( s _ { D } ).$$

The rank-two tensor in eq. (3.28) can be written as

$$W ^ { \mu \nu } = \omega _ { d t } ^ { \mu \nu } - \omega _ { i t } ^ { \mu \nu } + \frac { 2 m _ { t } ^ { 2 } ( p _ { i } p _ { d } ) } { ( p _ { t } p _ { i } ) ( p _ { t } p _ { d } ) } \omega _ { i d } ^ { \mu \nu },$$

where the quantity

$$\omega _ { x y } ^ { \mu \nu } = \frac { p _ { x } ^ { \mu } p _ { y } ^ { \nu } - p _ { y } ^ { \mu } p _ { x } ^ { \nu } } { ( p _ { x } p _ { y } ) },$$

describes an infinitesimal boost in the ( p x -p y ) plane. Eq. (3.28) provides an additional real emission contribution to the O ( λ ) terms when single top production process is combined with top quark decay in the narrow width approximation.

## 3.2 Virtual corrections and renormalisation in the production sub-process

As explained in ref. [1], the O ( λ ) contributions from the virtual corrections can be treated in the same way as the real emission ones. In this section we study the virtual corrections to the production sub-process of the process in eq. (2.1) and use the polarised top-quark spinor to describe the influence of the top quark decay.

Similarly to the real emission case, the O ( λ ) contributions to the virtual corrections can only arise from the region of soft k ∼ λ loop momenta. Our goal, therefore, is to establish a similar soft expansion of one-loop virtual corrections to the single top production process u p ( u ) + b p ( i ) → d p ( d ) + t ( p t ). We again consider corrections to the heavy-quark line since corrections to the light-quark line do not produce O ( λ ) contributions [19]. We write

$$\mathcal { A } _ { \text{virt} } = g _ { s } ^ { 2 } C _ { F } \delta _ { i j } \mathcal { M } _ { \text{virt} },$$

where i, j are the colour indices of the top quark and the bottom quark. Proceeding as in ref. [1], we find

$$\cdots _ { \, \iota ^ { - j } }, \, \cdots \, \cdots \, \cdots \\ & \mathcal { M } _ { \text{virt} } = \int \frac { \mathrm d ^ { 4 } k } { ( 2 \pi ) ^ { 4 } } \frac { - i } { k ^ { 2 } - \lambda ^ { 2 } } \left [ J _ { t } ^ { \alpha } J _ { i, \alpha } \, \bar { u } _ { t } \left ( \mathbf N ( p _ { t }, p _ { i },.. ) + k ^ { \mu } D _ { p, \mu } \mathbf N ( p _ { t }, p _ { i },.. ) \right ) u _ { i } \\ & - J _ { t } ^ { \alpha } \bar { u } _ { t } \mathbf N ( p _ { t }, p _ { i },.. ) \mathbf S _ { i, \alpha } u _ { i } + J _ { i } ^ { \alpha } \bar { u } _ { t } \mathbf S _ { t, \alpha } \mathbf N ( p _ { t }, p _ { i },.. ) u _ { i } - ( J _ { t } ^ { \alpha } + J _ { i } ^ { \alpha } ) \bar { u } _ { t } D _ { p, \alpha } \mathbf N u _ { i } \right ],$$

$$( 3. 3 2 )$$

where d t = ( p t + ) k 2 -m 2 t and d i = ( p i + ) , k 2

$$D _ { p } ^ { \mu } = \frac { \partial } { \partial p _ { t, \mu } } + \frac { \partial } { \partial p _ { i, \mu } },$$

and

$$J _ { t } ^ { \alpha } & = \frac { 2 p _ { t } ^ { \alpha } + k ^ { \alpha } } { d _ { t } }, \ \mathbf S _ { t } ^ { \alpha } = \frac { \sigma ^ { \alpha \beta } k _ { \beta } } { d _ { t } }, \\ J _ { i } ^ { \alpha } & = \frac { 2 p _ { i } ^ { \alpha } + k ^ { \alpha } } { d _ { i } }, \ \mathbf S _ { i } ^ { \alpha } = \frac { \sigma ^ { \alpha \beta } k _ { \beta } } { d _ { i } }. \\ \text{real emission case the randomone on the non momentum has been}$$

Similar to the real emission case, the dependence on the loop momentum has been made explicit so that the integration over k becomes straightforward. However, it is better to square the matrix element before integrating over the loop momentum k . We do this following ref. [1] and accounting for the fact that the top quark is polarised. We find

$$\text{wing ref. } [ 1 ] & \text{ and accounting for the fact that the top quark is polarised.  We find} \\ & \delta _ { \text{virt} } [ \mathcal { M } \mathcal { M } ^ { + } ] = \int \frac { d ^ { 4 } k } { ( 2 \pi ) ^ { 4 } } \frac { - i } { k ^ { 2 } - \lambda ^ { 2 } } \left [ 2 J _ { t } ^ { \alpha } J _ { i, \alpha } \text{Tr} \left [ \phi _ { t, D } \mathbb { N } \phi _ { i } \bar { N } \right ] \\ & + J _ { t } ^ { \alpha } J _ { i, \alpha } k ^ { \mu } \text{ Tr} \left [ \phi _ { t, D } ( D _ { p, \mu } \mathbb { N } ) \phi _ { i } \bar { N } + \phi _ { t, D } \mathbb { N } \phi _ { i } ( D _ { p, \mu } \bar { N } ) \right ] \\ & - ( J _ { t } ^ { \alpha } + J _ { i } ^ { \alpha } ) \text{Tr} \left [ \phi _ { t, D } ( D _ { p, \alpha } \mathbb { N } ) \phi _ { i } \bar { N } + \phi _ { t, D } \mathbb { N } \phi _ { i } ( D _ { p, \alpha } \bar { N } ) \right ] \\ & + J _ { i } ^ { \alpha } \text{Tr} \left [ [ \phi _ { t, D }, \mathbb { N } _ { t, \alpha } ] \mathbb { N } \phi _ { i } \bar { N } \right ] - J _ { t } ^ { \alpha } \text{Tr} \left [ \phi _ { t, D } \mathbb { N } [ \mathbb { N } _ { i, \alpha }, \phi _ { i } ] \bar { N } \right ] \right ], \\ \mathbb { \sim } \mathcal { M } - \mathcal { M } \rightharpoonup \mathcal { M } \dots \text{ the above equation content in all } \mathcal { M } / \mathcal { M } \text{ overwritten to } \mathcal { M } \mathcal { M } +$$

where M = M 0 + M virt . The above equation contains all O ( λ ) corrections to MM + .

We can further simplify eq. (3.35) following steps already discussed in the previous section where the real emission contribution was considered. Indeed, using

$$& [ \phi _ { t, D }, \mathbf S _ { t } ^ { \alpha } ] = ( L _ { t } ^ { \alpha } + S _ { t } ^ { \alpha } ) \phi _ { t, D }, \\ & [ \phi _ { i }, \mathbf S _ { i } ^ { \alpha } ] = L _ { i } ^ { \alpha } \phi _ { i },$$

we arrive at

$$\text{we arrive at} \\ & \delta [ \mathcal { M } \mathcal { M } ^ { + } ] _ { \text{virt} } = \int \frac { d ^ { 4 } k } { ( 2 \pi ) ^ { 4 } } \frac { - i } { k ^ { 2 } - \lambda ^ { 2 } } \left [ 2 J _ { t } ^ { \alpha } J _ { i, \alpha } \text{Tr} \left [ \phi _ { t, D } \mathbf N \phi _ { i } \bar { N } \right ] \\ & + J _ { t } ^ { \alpha } J _ { i, \alpha } k ^ { \mu } \text{ Tr} \left [ \phi _ { t, D } ( D _ { p, \mu } \mathbf N ) \phi _ { i } \bar { N } + \phi _ { t, D } \mathbf N \phi _ { i } ( D _ { p, \mu } \bar { N } ) \right ] \\ & - ( J _ { t } ^ { \alpha } + J _ { i } ^ { \alpha } ) \text{Tr} \left [ \phi _ { t, D } ( D _ { p, \alpha } \mathbf N ) \phi _ { i } \bar { N } + \phi _ { t, D } \mathbf N \phi _ { i } ( D _ { p, \alpha } \bar { N } ) \right ] \\ & + J _ { i } ^ { \alpha } \text{Tr} \left [ ( ( L _ { t, \alpha } + S _ { t, \alpha } ) \phi _ { t, D } ) \mathbf N \phi _ { i } \bar { N } \right ] + J _ { t } ^ { \alpha } \text{Tr} \left [ \phi _ { t, D } \mathbf N ( L _ { i, \alpha } \phi _ { i } ) \bar { N } \right ] \right ].$$

We stress that, similar to the real emission case, derivatives with respect to momenta do not act on the spin vector s D that appears in the density matrix / ρ t,D .

We note that the difference between the result in eq. (3.37) and a similar result computed for the unpolarised case in ref. [1], is the appearance of the spin-dependent density matrix / ρ t,D everywhere in eq. (3.37) and the presence of an additional term that contains the operator S t,α . This term evaluates to

$$\int \frac { \mathrm d ^ { 4 } k } { ( 2 \pi ) ^ { 4 } } \frac { - i } { k ^ { 2 } - \lambda ^ { 2 } } \ J _ { i } ^ { \alpha } \ S _ { t, \alpha } \text{Tr} \left [ \phi _ { t, D } \mathbf N \phi _ { i } \bar { \mathbf N } \right ] = - \frac { 1 } { 8 \pi ^ { 2 } } \frac { \pi \lambda } { m _ { t } } \, \frac { ( p _ { i } s D ) } { ( p _ { i } p _ { t } ) } p _ { t } ^ { \nu } \frac { \partial } { \partial s _ { D } ^ { \nu } } \, \text{Tr} \left [ \phi _ { t, D } \mathbf N \phi _ { i } \bar { \mathbf N } \right ], \, ( 3. 3 8 )$$

as follows from the integrals collected in appendix A of ref. [1]. Hence, we can write the virtual contribution as

$$\mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { P D } ^ { V, \text{prod} } \right ] = 2 \frac { \mathrm d \Gamma _ { t } } { \Gamma _ { t } } \ \mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { t } ^ { V } ( s _ { D } ) \right ] + \mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { P D, \text{new} } ^ { V, \text{prod} } \right ],$$

where the first term on the right hand side can be found in eq. (3.16) of ref. [1] and the second term is new. We note that in the first term we again need to replace F LO with F p ( ..., s D ) and ( / p t + m t ) with / ρ t,D to account for the polarisation effects. The second term in eq. (3.39) is a new contribution. It reads

$$\mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { P D, \text{new} } ^ { V, \text{prod} } \right ] = - \frac { \alpha _ { s } C _ { F } } { 2 \pi } \frac { \pi \lambda } { m _ { t } } \, 2 \, \frac { \mathrm d \Gamma _ { t } } { \Gamma _ { t } } \, s _ { D, \mu } \, \omega _ { i t } ^ { \mu \nu } \, \frac { \partial } { \partial s _ { D } ^ { \nu } } \, \mathrm d \sigma _ { t } ( s _ { D } ).$$

Finally, contributions due to wave function and (explicit) mass renormalisation are not affected by the fact that the top quark is polarised. Hence, we conclude that these contributions can be borrowed from ref. [1] without any modification

$$\mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { P D } ^ { \text{Ren,prod} } \right ] = 2 \frac { \mathrm d \Gamma _ { t } } { \Gamma _ { t } } \ \mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { t } ^ { \text{Ren} } ( s _ { D } ) \right ].$$

To summarise, O ( λ ) contributions to the cross section of the process in eq. (2.1) caused by the radiation of real and virtual gluons and the renormalisation in the production subprocess are obtained as the sum of the contributions given in eqs (3.27), (3.39) and (3.41). Each of these contributions is written as the sum of two terms: the 'old one' that are identical to the stable-top production case discussed in ref. [1], except for the fact that one has to employ there the polarised leading order cross section, and the 'new one' which is entirely due to the fact that there are spin correlations between production and decay processes. When single top production was considered in isolation, 'old corrections' were cancelling against the redefinition of the top quark mass parameter; a similar cancellation also exists in the current case. However, before discussing this point, we need to compute the O ( λ ) power correction to the decay sub-process. We do this in the next section.

## 4 Corrections to the decay sub-process

In this section we explain how the O ( λ ) power correction to the top quark decay subprocess is computed. Following the discussion in section 2, the top quark is polarised and

its polarisation vector s P is determined by the kinematics of the production sub-process. We note that corrections to the decay of an unpolarised top quark considered in isolation can be found in appendix B in ref. [1].

Our starting point is eq. (2.16). The momenta assignments differ from the ones in appendix B of ref. [1]; for this reason we emphasise that we consider the decay process

$$t ( p _ { t } ) \to \nu ( p _ { 1 } ) + e ^ { + } ( p _ { 2 } ) + b ( p _ { f } ).$$

The calculation of the real-emission contributions proceeds similarly to the case of the single top production and along the lines of appendix B of ref. [1]. We use the following momenta assignment to describe the real-emission process

$$t ( p _ { t } ) \to \nu ( p _ { 1 } ) + e ^ { + } ( q _ { 2 } ) + b ( q _ { f } ) + g ( k ),$$

with k 2 = λ 2 . We again use the Low-Burnett-Kroll theorem [20-22] shown in eq. (3.13) where for the process in eq. (4.2) we have

$$J ^ { \mu } = J _ { t } ^ { \mu } + J _ { f } ^ { \mu }, \ \ J _ { t } ^ { \mu } = \frac { 2 p _ { t } ^ { \mu } - k ^ { \mu } } { d _ { t } }, \ \ J _ { f } ^ { \mu } = \frac { 2 q _ { f } ^ { \mu } + k ^ { \mu } } { d _ { f } },$$

with d t = ( p t -k ) 2 -m 2 t and d f = ( q f + ) . k 2

In order to factorise the phase space of the gluon momentum from the rest of the decay phase space, we employ a momentum mapping. In variance with the case of the production sub-process, this mapping does not need to involve the top quark momentum and, hence, the production process remains unaffected. Following ref. [1], we write

$$q _ { f } = p _ { f } - k + \frac { ( p _ { f } k ) } { ( p _ { 2 } p _ { f } ) } \ p _ { 2 }, \quad q _ { 2 } = \left ( 1 - \frac { ( p _ { f } k ) } { ( p _ { 2 } p _ { f } ) } \right ) p _ { 2 }.$$

Using this transformation, the phase space changes as follows [1]

$$\mathrm d \Phi _ { D } ( p _ { t } ; p _ { 1 }, q _ { 2 }, q _ { f }, k ) = \mathrm d \Phi _ { D } ( p _ { t } ; p _ { 1 }, p _ { 2 }, p _ { f } ) \ [ \mathrm d k ] _ { \lambda } \left ( 1 + \frac { ( p _ { 2 } k ) } { ( p _ { f } p _ { 2 } ) } - \frac { ( p _ { f } k ) } { ( p _ { f } p _ { 2 } ) } + \mathcal { O } ( \lambda ^ { 2 } ) \right ).$$

Since these momenta transformations do not impact p t and, therefore, the 'spin' of the top quark as defined by the production process, the only addition to the unpolarised case for the width arises because of the analog of the S µ t term in eq. (3.13) which is already O ( λ ) and, hence, can be easily integrated over k . We therefore find

$$\mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { P D } ^ { R, \mathrm d e c } \right ] = 2 \frac { \mathcal { T } _ { \lambda } \left [ \mathrm d \Gamma _ { t } ^ { R } ( s _ { P } ) \right ] } { \Gamma _ { t } } \, \mathrm d \sigma _ { t } + \mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { P D, \ \text{new} } ^ { R, \mathrm d e c } \right ],$$

where in the first term formulas from unpolarised case can be employed except that the leading order matrix element squared should be replaced with the polarised one. The second term is new; after integration over the momentum of the soft gluon it evaluates to

$$\mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { P D, \text{new} } ^ { R, \text{dec} } \right ] = - \frac { \alpha _ { s } C _ { F } } { 2 \pi } \frac { \pi \lambda } { m _ { t } } \, 2 \, \frac { \mathrm d \sigma _ { t } } { \Gamma _ { t } } \, s _ { P, \mu } \, \omega _ { t f } ^ { \mu \nu } \frac { \partial } { \partial s _ { P } ^ { \nu } } \, \mathrm d \Gamma _ { t } ( s _ { P } ).$$

We also need to compute virtual corrections and perform mass and wave function renormalisation for the decay sub-process. The virtual corrections are computed in the same way as what was described for the production sub-process and in appendix B of ref. [1]. We find

$$\mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { P D } ^ { V, \text{dec} } \right ] = 2 \frac { \mathcal { T } _ { \lambda } \left [ \mathrm d \Gamma _ { t } ^ { V } ( s _ { P } ) \right ] } { \Gamma _ { t } } \, \mathrm d \sigma _ { t } + \mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { P D, \text{new} } ^ { V, \text{dec} } \right ],$$

where

$$\mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { P D, \text{new} } ^ { V, \text{dec} } \right ] = \frac { \alpha _ { s } C _ { F } } { 2 \pi } \frac { \pi \lambda } { m _ { t } } \, 2 \, \frac { \mathrm d \sigma _ { t } } { \Gamma _ { t } } \, s _ { P, \mu } \, \omega _ { t f } ^ { \mu \nu } \frac { \partial } { \partial s _ { P } ^ { \nu } } \, \mathrm d \Gamma _ { t } ( s _ { P } ).$$

The renormalisation contributions are not affected by the polarisation of the top quark and, therefore, can be directly borrowed from the results in appendix B in ref. [1] except that the differential decay width has to be computed for the polarised top quark. Hence, we write

$$\mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { P D } ^ { \text{Ren}, \text{dec} } \right ] = 2 \frac { \mathcal { T } _ { \lambda } \left [ \mathrm d \Gamma _ { t } ^ { \text{Ren} } ( s _ { P } ) \right ] } { \Gamma _ { t } } \, \mathrm d \sigma _ { t }.$$

## 5 Redefinition of the top quark mass parameter

In cases when the top quark production and the top quark decay are considered separately, it is known [1] that the cancellation of O ( λ ) contributions is only possible if the production cross section and the decay rate are expressed through a short-distance top quark mass and not through the pole mass . For reasons of technical convenience, we performed the renormalisation in the on-shell scheme, similar to what was done in ref. [1]. 6 To derive the final results, we need to switch to a short-distance mass parameter. We explain below how to do this in case when top quark production and decay are considered simultaneously.

As explained in ref. [1], we can switch to a short-distance mass parameter by redefining momenta of final-state particles. Our goal will be to do this in such a way that, when spin correlations are neglected, we obtain formulas which are identical to the ones in ref. [1], where production and decay are considered separately.

We begin with the momenta transformations for particles that originate from top decay and write

$$q _ { f } ^ { \mu } = \tilde { p } _ { f } ^ { \mu } - \kappa q _ { t } + \kappa \, \frac { ( \tilde { p } _ { f } q _ { t } ) } { ( \tilde { p } _ { f } \tilde { p } _ { 2 } ) } \, \tilde { p } _ { 2 } ^ { \mu }, \quad q _ { 2 } ^ { \mu } = \tilde { p } _ { 2 } ^ { \mu } \left ( 1 - \kappa \, \frac { ( \tilde { p } _ { f } q _ { t } ) } { ( \tilde { p } _ { f } \tilde { p } _ { 2 } ) } \right ).$$

This momentum transformation leads to the following change in the decay phase space [1]

$$\mathrm d \Phi _ { D } ( q _ { t } ; q _ { f }, q _ { 2 }, q _ { 1 } ) = \mathrm d \Phi _ { D } ( ( 1 + \kappa ) q _ { t } ; \tilde { p } _ { f }, \tilde { p } _ { 2 }, q _ { 1 } ) \left ( 1 + \kappa \, \frac { ( \tilde { p } _ { 2 } q _ { t } ) } { ( \tilde { p } _ { f } \tilde { p } _ { 2 } ) } - \kappa \, \frac { ( \tilde { p } _ { f } q _ { t } ) } { ( \tilde { p } _ { f } \tilde { p } _ { 2 } ) } + \mathcal { O } ( \lambda ^ { 2 } ) \right ). \ \ ( 5. 2 )$$

It follows from the above equation that the mass of the decaying top quark becomes

$$\tilde { m } _ { t } = \sqrt { ( 1 + \kappa ) ^ { 2 } q _ { t } ^ { 2 } } = ( 1 + \kappa ) m _ { t }.$$

Hence,

$$\tilde { m } _ { t } - m _ { t } = \kappa m _ { t },$$

6 Acalculation that directly uses the mass parameter defined in a short-distance scheme, see the appendix.

which implies that κm t is the shift in the mass parameter and where κ is defined as in ref. [1],

$$\kappa = \frac { \alpha _ { s } C _ { F } } { 2 \pi } \frac { \pi \lambda } { m _ { t } }.$$

To proceed further, it is convenient to define the top quark momentum that appears in the decay phase space

$$\tilde { p } _ { t } & = ( 1 + \kappa ) q _ { t }.$$

When the production and decay processes are considered together, the top quark momentum appears in the phase space of the production sub-process; hence, the above redefinition will modify the production phase space and the matrix element. We begin with the analysis of the production phase space and write 7

$$d \Phi _ { P } ( \dots ; q _ { t }, q _ { d } ) & = \frac { d ^ { 4 } q _ { t } } { ( 2 \pi ) ^ { 3 } } \delta ( q _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) \ [ d q _ { d } ] \ ( 2 \pi ) ^ { 4 } \delta ( p _ { u } + p _ { i } - q _ { t } - q _ { d } ) \\ & = ( 1 - 2 \kappa ) \frac { d ^ { 4 } \tilde { p } _ { t } } { ( 2 \pi ) ^ { 4 } } \delta ( \tilde { p } _ { t } ^ { 2 } - \tilde { m } _ { t } ^ { 2 } ) [ d q _ { d } ] ( 2 \pi ) ^ { 4 } \delta ( p _ { u } + p _ { i } - \tilde { p } _ { t } + \kappa \tilde { p } _ { t } - q _ { d } ).$$

We then perform one more momentum redefinition but this time we change the momenta in such a way that the top quark remains on a (new) mass shell. We write

$$\tilde { p } _ { t } = p _ { t } + \kappa p _ { t } - \frac { \kappa m _ { t } ^ { 2 } } { ( p _ { t } p _ { d } ) } p _ { d }, \quad q _ { d } = \left ( 1 + \frac { \kappa m _ { t } ^ { 2 } } { ( p _ { t } p _ { d } ) } \right ) p _ { d }.$$

This gives

$$\mathrm d \Phi _ { P } ( p _ { u }, p _ { i } ; q _ { t }, q _ { d } ) = \mathrm d \Phi _ { P } ( p _ { u }, p _ { i } ; p _ { t }, p _ { d } ) \left ( 1 + \frac { \kappa m _ { t } ^ { 2 } } { ( p _ { t } p _ { d } ) } + \mathcal { O } ( \lambda ^ { 2 } ) \right ),$$

where p 2 t = ˜ m 2 t .

We note that the change ˜ p t → p t impacts the decay phase space again . However, it is easy to solve this problem because this momentum change can be written as a Lorentz transformation

$$\tilde { p } _ { t } ^ { \mu } & = \Lambda _ { m } ^ { \mu \nu } p _ { t, \nu },$$

where

$$\Lambda _ { m } ^ { \mu \nu } = g ^ { \mu \nu } + \kappa \, \omega _ { t d } ^ { \mu \nu }.$$

It follows that the decay phase transforms as follows

$$d \Phi _ { D } \left ( \tilde { p } _ { t } ; \{ \tilde { p } _ { \text{dec} } \} \right ) = \text{d} \Phi _ { D } ( \Lambda _ { m } p _ { t }, \{ \Lambda _ { m } p _ { \text{dec} } \} ) = \text{d} \Phi _ { D } ( p _ { t }, \{ p _ { \text{dec} } \} ).$$

We have worked out the momenta transformations required to modify the mass of a top quark in a process where it is produced and then decays. We now need to combine the several transformations and write down formulas that elucidate phase-space and matrixelement transformations of the full process. We use the representation shown in eq. (2.14)

̸

7 To derive this formula, one needs to account for the fact that q 2 t = m 2 t a priori. Hence, δ q ( 2 t -m 2 t ) = δ ((1 -2 κ )(˜ p 2 t -˜ m 2 t )) = (1 + 2 κ δ p ) (˜ 2 t -˜ m 2 t ).

where unpolarised decay width and polarised production cross section are combined. We use the invariance of the decay matrix element squared F d under Lorentz transformations and find

$$& d \Phi _ { P } ( p _ { u }, p _ { i } ; q _ { t }, q _ { d } ) F _ { p } ( p _ { u }, p _ { i } ; q _ { t }, q _ { d }, s _ { D } ( q _ { t }, q _ { 2 } ) ) \, d \Phi _ { D } ( q _ { t } ; q _ { f }, q _ { 2 }, q _ { 1 } ) F _ { d } ( q _ { t } ; q _ { f }, q _ { 2 }, q _ { 1 } ) \\ & = d \Phi _ { P } \, j _ { p } \, F _ { p } \Big ( \dots ; \ p _ { t } - \kappa \xi _ { d } \, p _ { d }, \, ( 1 + \kappa \xi _ { d } ) \, p _ { d }, \, \Lambda _ { m } s _ { D } \Big ( \, ( 1 - \kappa ) \, p _ { t }, ( 1 - \kappa \xi _ { 2 } ) \, p _ { 2 } ) \Big ) \\ & \quad \times d \Phi _ { D } \, j _ { d } \, F _ { d } \Big ( \, ( 1 - \kappa ) \, p _ { t } ; \ p _ { f } + \kappa \delta p _ { f }, \, ( 1 - \kappa \xi _ { 2 } ) \, p _ { 2 }, \ p _ { 1 } \Big ) + \mathcal { O } ( \kappa ^ { 2 } ),$$

where the phase spaces dΦ p,d depend on the transformed momenta { p } , Λ m is the boost defined in eq. (5.11) and

$$\begin{array} { c c c } \cdot \cdot & \cdot \\ j _ { p } = 1 + \frac { \kappa m _ { t } ^ { 2 } } { ( p _ { t } p _ { d } ) }, & j _ { d } = 1 + \kappa \, \frac { ( p _ { 2 } p _ { t } ) } { ( p _ { f } \tilde { p } _ { 2 } ) } - \kappa \, \frac { ( p _ { f } p _ { t } ) } { ( p _ { f } p _ { 2 } ) }, \\ & \xi _ { d } = \frac { m _ { t } ^ { 2 } } { ( p _ { t } p _ { d } ) }, & \xi _ { 2 } = \frac { ( p _ { f } p _ { t } ) } { ( p _ { f } p _ { 2 } ) }, & \delta p _ { f } ^ { \mu } = - p _ { t } ^ { \mu } + \frac { ( p _ { f } p _ { t } ) } { ( p _ { f } p _ { 2 } ) } \, p _ { 2 } ^ { \mu }. \\ & \cdot \cdot \cdot \cdot \cdot \end{array} \cdot \end{array} \cdot \end{array} \cdot \end{array}$$

To obtain O ( λ ) correction to the cross section of the process in eq. (2.1) related to mass redefinition, we expand eq. (5.13) in κ and keep linear terms. These terms can be combined into three groups:

- 1. the term that originates from the expansion of F p caused by the O ( κ ) contribution to the matrix Λ m acting on s D ;
- 2. all O ( κ ) terms that appear from the expansion of the second line in eq. (5.13) but without a term discussed in the previous item and without a correction to the argument of the spin vector s D ;
- 3. terms that originate from the expansion of the third line in eq. (5.13) in κ and terms that originate from the expansion of the argument of spin vector s D in function F p .

We now discuss these three groups of terms separately. The term in the first item is new. Terms in the second item provide the required contribution to cancel all 'old' O ( λ ) corrections to the production sub-process discussed in section 3. Note that this cancellation occurs for the polarised matrix element squared F p ( ..., s D ) since this is what appears in eq. (5.13).

The contribution of the third group of terms should, in principle, cancel all 'old' O ( λ ) terms to the decay sub-process, described in section 4. However, it follows from eq. (5.13) that this contribution lacks the polarisation vector s P , which is present in the similar contributions in section 4. Hence, to claim this cancellation, we need to put it back into the decay matrix element squared. This is possible because the following identity holds

$$& F _ { p } \Big ( \dots ; p _ { t }, p _ { d }, s _ { D } \Big ( \, ( 1 - \kappa ) \, p _ { t }, \ ( 1 - \kappa \xi _ { 2 } ) \, p _ { 2 } \Big ) \Big ) F _ { d } \Big ( \, ( 1 - \kappa ) \, p _ { t } ; \ p _ { f } + \kappa \delta p _ { f }, \ ( 1 - \kappa \xi _ { 2 } ) p _ { 2 }, \ p _ { 1 } \Big ) \\ & = F _ { p } \Big ( \dots ; p _ { t }, p _ { d } \Big ) F _ { d } \Big ( \, ( 1 - \kappa ) \, p _ { t } ; \ p _ { f } + \kappa \delta p _ { f }, \ ( 1 - \kappa \xi _ { 2 } ) p _ { 2 }, \ p _ { 1 }, \ s _ { P } ( p _ { t }, p _ { d } ) \Big ),$$

thanks to the relation between polarised production and decay cross sections shown in eq. (2.20).

We conclude that the only new term that we need to consider is the term in the first item that arises from the boost of the spin vector. It evaluates to

$$\mathcal { T } _ { \lambda } \left [ \mathrm d \sigma _ { \max } ^ { \text{new} } \right ] = - \frac { \alpha _ { s } C _ { F } } { 2 \pi } \frac { \pi \lambda } { m _ { t } } \, 2 \, \frac { \mathrm d \Gamma _ { t } } { \Gamma _ { t } } s _ { D, \mu } \, \omega _ { t d } ^ { \mu \nu } \, \frac { \partial } { \partial s _ { D } ^ { \nu } } \, \mathrm d \sigma _ { t } ( s _ { D } ).$$

Other terms that arise from the mass redefinition combine with 'old' contributions to the production and decay sub-processes and cancel in the same way as discussed in ref. [1].

## 6 Results and corrections to observables

The final result is obtained by combining the O ( λ ) contributions to single top production and decay process derived in sections 3, 4 and 5. As we argued extensively during the calculation many O ( λ ) contributions cancel in the sum; the only ones that survive involve polarisation effects which is an important new feature of a process with a long-lived particle that is first produced and then decays. They are obtained by adding eqs (3.28, 3.40, 4.7, 4.9, 5.16). We find

$$\mathcal { T } _ { \lambda } [ d \sigma _ { P D } ] = \frac { \alpha _ { s } C _ { F } } { 2 \pi } \ \frac { \pi \lambda } { m _ { t } } \ 2 \ \frac { \text{d} \Gamma _ { t } } { \Gamma _ { t } } \ s _ { D, \mu } \left ( 2 \omega _ { t i } ^ { \mu \nu } + 2 \omega _ { d t } ^ { \mu \nu } + \frac { 2 m _ { t } ^ { 2 } ( p _ { i } p _ { d } ) } { ( p _ { t } p _ { i } ) ( p _ { t } p _ { d } ) } \omega _ { i d } ^ { \mu \nu } \right ) \frac { \partial } { \partial s _ { D } ^ { \nu } } \ \text{d} \sigma _ { t } ( s _ { D } ). \ \ ( 6. 1 )$$

We can use eq. (2.20) as well as the relations s P · p t = s D · p t = 0, to find

$$\mathcal { T } _ { \lambda } [ \mathrm d \sigma _ { P D } ] = - \frac { \alpha _ { s } C _ { F } } { 2 \pi } \ \frac { \pi \lambda } { m _ { t } } \ \frac { \mathrm d \Gamma _ { t } \mathrm d \sigma _ { t } } { \Gamma _ { t } } \ \frac { 2 m _ { t } ^ { 2 } ( p _ { i } p _ { d } ) } { ( p _ { t } p _ { i } ) ( p _ { t } p _ { d } ) } \ s _ { D, \mu } \omega _ { i d } ^ { \mu \nu } s _ { P, \nu }.$$

The expression in eq. (6.2) assumes a particularly simple form in the top quark rest frame . Indeed, in this case

$$s _ { D } ^ { \mu } = ( 0, \vec { n } _ { 2 } ) \,, \ \ s _ { P } ^ { \mu } = ( 0, \vec { n } _ { d } ),$$

where ⃗ n 2 and ⃗ n d are unit vectors aligned with directions of the positron and d -quark in this frame, respectively. This implies that

$$\frac { 2 m _ { t } ^ { 2 } ( p _ { i } p _ { d } ) } { ( p _ { t } p _ { i } ) ( p _ { t } p _ { d } ) } \ s _ { D, \mu } \omega _ { i d } ^ { \mu \nu } s _ { P, \nu } = 2 \left ( ( \vec { n } _ { 2 } \cdot \vec { n } _ { i } ) - ( \vec { n } _ { 2 } \cdot \vec { n } _ { d } ) ( \vec { n } _ { i } \cdot \vec { n } _ { d } ) \right ) = 2 \left [ \vec { n } _ { 2 } \times \vec { n } _ { d } \right ] \cdot [ \vec { n } _ { i } \times \vec { n } _ { d } ], \ \ ( 6. 4 )$$

where ⃗ n i is the direction of the incoming b quark in the top quark rest frame. We conclude that, in the top quark rest frame, eq. (6.2) takes a remarkably simple form

$$\mathcal { T } _ { \lambda } [ \mathrm d \sigma _ { P D } ] = - \frac { \alpha _ { s } C _ { F } } { 2 \pi } \ \frac { \pi \lambda } { m _ { t } } \ \frac { \mathrm d \Gamma _ { t } \mathrm d \sigma _ { t } } { \Gamma _ { t } } \ 2 \ [ \vec { n } _ { 2 } \times \vec { n } _ { d } ] \cdot [ \vec { n } _ { i } \times \vec { n } _ { d } ].$$

For the case, when the observable does not depend on decay momenta, we recover our previous result presented in ref. [1]. Indeed, integrating over the total phase space of the decay products and using the fact that

$$\int d \Gamma _ { t } \, \vec { n } _ { 2 } = 0,$$

in the top quark rest frame, we recover the stable-top-quark result [1]

$$\mathcal { T } _ { \lambda } [ \sigma _ { P D } ] = 0.$$

Nevertheless, eq. (6.5) shows that, in general, there is an O ( λ ) contribution to the differential cross section related to polarisation effects. However, because we have used momenta redefinitions to derive this result, it is important to account for them also in the observables since they are defined using the original momenta.

To this end, we consider an observable X and study the following integral

$$O _ { X } = \int \mathrm d \sigma _ { P D } \, X.$$

In principle, the observable X is generic; however, we would like to focus upon observables that are used in practice to study polarisation effects in single top production [14-16]. For this reason, we assume that the observable X depends on the top quark momentum, the d -quark momentum, the incoming b -quark momentum (the collision axis) and the positron momentum

$$X = X ( q _ { t }, q _ { d }, q _ { i }, q _ { 2 } ).$$

When different contributions to eq. (6.8) are studied and different mappings are performed, there will be shifts in the arguments of the function X that are proportional to the gluon momentum k or to the mass-redefinition parameter κ . We are interested in terms that originate in the expansion of the function X in these small parameters.

There are three contributions that affect the arguments of X : real radiation in production, real radiation in decay and mass redefinition. As the first step, we summarise the momenta redefinitions for each of these contributions. Since none of these momenta redefinitions changes the collision axis, we will not show q i among the arguments of X in what follows. We find:

- · Radiation in the production subprocess:

$$X \rightarrow X \left ( \Lambda p _ { t }, \left ( 1 - \frac { ( p _ { t } k ) } { ( p _ { t } p _ { d } ) } \right ) p _ { d }, \Lambda p _ { 2 } \right ),$$

where

$$\Lambda ^ { \mu \nu } = g ^ { \mu \nu } + \frac { p _ { d } ^ { \mu } k ^ { \nu } - k ^ { \mu } p _ { d } ^ { \nu } } { ( p _ { t } p _ { d } ) }.$$

- · Radiation in the decay subprocess:

$$X \rightarrow X \left ( p _ { t }, p _ { d }, \left ( 1 - \frac { ( p _ { f } k ) } { ( p _ { f } p _ { 2 } ) } \right ) p _ { 2 } \right ),$$

- · Mass redefinition:

$$X \rightarrow X \left ( ( 1 - \kappa ) \, \Lambda _ { m } p _ { t }, \left ( 1 + \kappa \frac { m _ { t } ^ { 2 } } { ( p _ { t } p _ { d } ) } \right ) p _ { d }, \left ( 1 - \kappa \frac { ( p _ { f } p _ { t } ) } { ( p _ { f } p _ { 2 } ) } \right ) \Lambda _ { m } p _ { 2 } \right ), \quad \quad ( 6. 1 3 )$$

where

$$\Lambda _ { m } ^ { \mu \nu } = g ^ { \mu \nu } + \kappa \frac { p _ { t } ^ { \mu } p _ { d } ^ { \nu } - p _ { d } ^ { \mu } p _ { t } ^ { \nu } } { ( p _ { t } p _ { d } ) }.$$

We then expand X in series for each of the three contributions and integrate over the gluon momentum k where appropriate. This is straightforward, and the matrix element squared is only needed in the eikonal approximation. The result reads

$$\delta X = \frac { \alpha _ { s } C _ { F } } { 2 \pi } \, \frac { \pi \lambda } { m _ { t } } \left [ \left ( p _ { t } ^ { \mu } - \frac { 2 m _ { t } ^ { 2 } } { ( p _ { t } p _ { i } ) } p _ { i } ^ { \mu } \right ) \frac { \partial X } { \partial p _ { t } ^ { \mu } } - 2 p _ { 2, \nu } \left ( \omega _ { t d } ^ { \nu \mu } + \frac { m _ { t } ^ { 2 } ( p _ { d } p _ { i } ) } { ( p _ { t } p _ { d } ) ( p _ { t } p _ { i } ) } \omega _ { d i } ^ { \nu \mu } \right ) \frac { \partial X } { \partial p _ { 2 } ^ { \mu } } \right ].$$

We note that the above result assumes that the mass parameter does not appear in the definition of the observable; if this is not the case, the mass parameter needs to be replaced with √ p 2 t .

Eq. (6.15) is applicable to any observable; the only constraint is that it can only depend on the momenta of final-state particles shown in eq. (6.9). Hence, we conclude that the complete linear correction to the expectation value of such an observable reads

$$& \mathcal { T } _ { \lambda } [ \mathcal { O } _ { X } ] = \frac { \alpha _ { s } C _ { F } } { 2 \pi } \, \frac { \pi \lambda } { m _ { t } } \, \int \frac { \mathrm d \Gamma _ { t } \mathrm d \sigma _ { t } } { \Gamma _ { t } } \, \left [ \, - \, \frac { 2 m _ { t } ^ { 2 } ( p _ { i } p _ { d } ) } { ( p _ { t } p _ { i } ) ( p _ { t } p _ { d } ) } \, s _ { D, \mu } \omega _ { i d } ^ { \mu \nu } s _ { P, \nu } \, X \\ & + ( 1 - s _ { D } \cdot s _ { P } ) \left [ \left ( p _ { t } ^ { \mu } - \frac { 2 m _ { t } ^ { 2 } } { ( p _ { t } p _ { i } ) } p _ { i } ^ { \mu } \right ) \frac { \partial X } { \partial p _ { t } ^ { \mu } } - 2 p _ { 2, \nu } \left ( \omega _ { t d } ^ { \nu \mu } + \frac { m _ { t } ^ { 2 } ( p _ { d } p _ { i } ) } { ( p _ { t } p _ { d } ) ( p _ { t } p _ { i } ) } \omega _ { i d } ^ { \nu \mu } \right ) \frac { \partial X } { \partial p _ { 2 } ^ { \mu } } \right ] \, \right ].$$

We will now analyse this general formula. First we note that one can consider observables that depend on the top quark momentum, but are inclusive with respect to the momenta of its decay products. Then X is a function of p t only. For such observables, we can integrate over the momenta of the top quark decay products. Then, considering the integrand in the top quark rest frame and using eq. (6.6), we conclude that the first term on the right hand side in eq. (6.16) vanishes. The second term then coincides with the correction to observable discussed in ref. [1] and the last term vanishes if X is a function of p t only.

There are also observables that are designed to study polarisation effects in single top production. Perhaps the simplest observable that belongs to this class is the one used by the CMS collaboration where the asymmetry between the direction of the outgoing light jet in single top production ( d -jet in our case) and the direction of positron in top decay is studied in the top quark rest frame [15]. We can construct such an observable by simply multiplying the production and decay spin polarisation vectors, s P and s D . Since in the top rest frame

$$s _ { D } \cdot s _ { P } = - \vec { n } _ { 2 } \cdot \vec { n } _ { d } = - \cos \theta _ { d 2 },$$

any function of this variable will provide a probe of polarisation effects; the observable used by the CMS collaboration corresponds to

$$X _ { \text{CMS} } = \theta ( - s _ { D } s _ { P } ) - \theta ( s _ { D } s _ { P } ).$$

Using eq. (6.15), it is easy to show that

$$\delta X _ { \text{CMS} } = 0,$$

which then implies that the only relevant term in eq. (6.16) that contributes for such observables is the first term in the integrand of eq. (6.16).

It is interesting to note that one can arrive at the same result without any computation. In fact, there is a simple argument that can be used to argue that for any observable X that depends upon the directions of p t , p d and p 2 only 8 there cannot be any change in X after the remapping described in this paper. To illustrate this argument, consider radiation in the production. According to eq. (6.10) the momenta redefinitions lead to

$$X ( p _ { t }, p _ { d }, p _ { 2 } ) \to X \left ( \Lambda p _ { t }, \left ( 1 - \frac { ( p _ { t } k ) } { ( p _ { t } p _ { d } ) } \right ) p _ { d }, \Lambda p _ { 2 } \right ) = X \left ( \Lambda p _ { t }, \left ( 1 + \frac { ( p _ { d } k ) } { ( p _ { t } p _ { d } ) } \right ) p _ { d }, \Lambda p _ { 2 } \right ), \\ \cdot \quad \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \end{array}$$

where in the last step we used the fact that the observable X depends on the direction of p d . This implies that the exact form of rescaling is irrelevant, and we can change it at will. Since

$$\left ( 1 + \frac { ( p _ { d } k ) } { ( p _ { t } p _ { d } ) } \right ) p _ { d } = \Lambda p _ { d },$$

we find

$$X \left ( \Lambda p _ { t }, \left ( 1 + \frac { ( p _ { d } k ) } { ( p _ { t } p _ { d } ) } \right ) p _ { d }, \Lambda p _ { 2 } \right ) = X \left ( \Lambda p _ { t }, \Lambda p _ { d }, \Lambda p _ { 2 } \right ) = X \left ( p _ { t }, p _ { d }, p _ { 2 } \right ),$$

where in the last step Lorentz invariance of the observable was used. Hence, we conclude that the momenta redefinitions employed in the description of the real emission in production do not change an observable which depends on directions of final-state particles. The same reason also applies to the momenta transformations employed to describe radiation in decay and the mass redefinition.

To complete the analysis of the CMS asymmetry, we need to understand the fate of the first term in eq. (6.16). Considering this term in the top rest frame, we find that it involves the following integral

$$\int d \Gamma _ { t } \left [ \vec { n } _ { 2 } \times \vec { n } _ { d } \right ] X ( \vec { n } _ { 2 } \cdot \vec { n } _ { d } ).$$

Since for any function X

$$\cdots \int d \Gamma _ { t } \, \vec { n } _ { 2 } \, X ( \vec { n } _ { 2 } \cdot \vec { n } _ { d } ) \sim \vec { n } _ { d },$$

the integral in eq. (6.23) vanishes. We conclude that the asymmetries in single top production studied by the CMS collaboration [14, 15] are not affected by the non-perturbative effects that can be modelled with renormalons.

A more complex polarisation observable was studied by the ATLAS collaboration [16]. To define it, a reference system in the top rest frame is introduced, where the three axes

8 We note that s D and s P belong to this category.

are 9

$$\vec { e } _ { z } = \vec { n } _ { d }, \quad \vec { e } _ { y } = \frac { \vec { n } _ { i } \times \vec { n } _ { d } } { | \vec { n } _ { i } \times \vec { n } _ { d } | }, \quad \vec { e } _ { x } = \vec { e } _ { y } \times \vec { e } _ { z } = \frac { \vec { n } _ { i } \times \vec { n } _ { d } } { | \vec { n } _ { i } \times \vec { n } _ { d } | } \times \vec { n } _ { d }.$$

The observable Q is defined as follows

$$Q ( \vec { n } _ { 2 }, \{ \vec { e } \} ) = 4 \theta ( \vec { n } _ { 2 } \cdot \vec { e } _ { z } ) + 2 \theta ( \vec { n } _ { 2 } \cdot \vec { e } _ { x } ) + \theta ( \vec { n } _ { 2 } \cdot \vec { e } _ { y } ).$$

We now determine the expectation value of Q at leading order and the non-perturbative correction to it. First, writing the leading order cross section using the reference frame described above, we obtain

$$d \sigma _ { P D } = d \sigma _ { t } \frac { d \Gamma _ { t } } { \Gamma _ { t } } \left ( 1 + \vec { e } _ { z } \cdot \vec { n } _ { 2 } \right ).$$

If we integrate over the top quark decay products without imposing any cuts on final-state particles, the following equations hold

$$\int \mathrm d \Gamma _ { t } \, \theta ( \vec { n } _ { 2 } \cdot \vec { a } ) = \frac { 1 } { 2 } \Gamma _ { t }, \quad \int \mathrm d \Gamma _ { t } \, \vec { n } _ { 2 } \, \theta ( \vec { n } _ { 2 } \cdot \vec { a } ) = \frac { 1 } { 4 } \Gamma _ { t } \, \vec { a },$$

where ⃗ a is an arbitrary unit vector. We use eq. (6.28) together with the leading order cross section in eq. (6.27) to find

$$\langle Q ( \vec { n } _ { 2 }, \{ \vec { e } \} ) \rangle = \frac { \int d \sigma _ { P D } \, Q ( \vec { n } _ { 2 }, \{ \vec { e } \} ) } { \int d \sigma _ { P D } } = \frac { 9 } { 2 }.$$

To compute the power corrections to this result, we need to combine the corrections to the cross section and to the observable. We begin with the latter. The correction to the observable is computed using eq. (6.15). To apply this equation to the observable Q , we should write it in a Lorentz-covariant form. To this end, we write

$$Q ( \vec { n } _ { 2 }, \{ \vec { e } \} ) & = 4 \theta ( \hat { Q } _ { z } ) + 2 \theta ( \hat { Q } _ { x } ) + \theta ( \hat { Q } _ { y } ),$$

with

$$\hat { Q } _ { x } & = \vec { n } _ { 2 } \cdot \vec { e } _ { x } = \frac { 1 } { | \vec { n } _ { i } \times \vec { n } _ { d } | } \, \left [ \, ( \vec { n } _ { 2 } \cdot \vec { n } _ { d } ) \, ( \vec { n } _ { i } \cdot \vec { n } _ { d } ) - ( \vec { n } _ { i } \cdot \vec { n } _ { 2 } ) \, \right ], \\ \hat { Q } _ { y } & = \vec { n } _ { 2 } \cdot \vec { e } _ { y } = \frac { 1 } { | \vec { n } _ { i } \times \vec { n } _ { d } | } \frac { p _ { t } ^ { 2 } } { ( p _ { t } p _ { 2 } ) \, ( p _ { t } p _ { i } ) \, ( p _ { t } p _ { d } ) } \, \epsilon _ { \mu \nu \rho \sigma } \, p _ { t } ^ { \mu } p _ { 2 } ^ { \nu } p _ { i } ^ { \rho } p _ { d } ^ { \sigma }, \\ \hat { Q } _ { z } & = \vec { n } _ { 2 } \cdot \vec { e } _ { z } = \vec { n } _ { 2 } \cdot \vec { n } _ { d }, \\. & \quad. & \quad. & \quad. & \quad. & \quad. & \quad.$$

and note that the covariant generalisation of the scalar product of two vectors in the top rest frame is given by

$$\vec { n } _ { i } \cdot \vec { n } _ { j } = 1 - \frac { p _ { t } ^ { 2 } \left ( p _ { i } p _ { j } \right ) } { \left ( p _ { t } p _ { i } \right ) \left ( p _ { t } p _ { j } \right ) }.$$

9 We note that in our calculation ⃗ n i denotes the direction of the incoming b -quark in the top quark rest frame, whereas in the ATLAS paper [16] the direction of the incoming light quark is chosen to define the reference system. These vectors are not back-to-back in the top quark rest frame but, thanks to momentum conservation, in this reference frame their vector products with ⃗ n d are the same up to a sign.

After that, the calculation becomes straightforward. We obtain

$$\delta \hat { Q } _ { x } & = \, 0, & \delta \hat { Q } _ { y } & = \, 0, & \delta \hat { Q } _ { z } & = \, 0,$$

$$( 6. 3 3 )$$

so that also in this case there is no change in the observable

$$\delta Q & = 0.$$

Written in the reference system defined in eq. (6.25), the non-perturbative shift in the cross section shown in eq. (6.5) becomes

$$\mathcal { T } _ { \lambda } [ \mathrm d \sigma _ { P D } ] = \frac { \alpha _ { s } C _ { F } } { 2 \pi } \, \frac { \pi \lambda } { m _ { t } } \, \frac { \mathrm d \Gamma _ { t } \mathrm d \sigma _ { t } } { \Gamma _ { t } } \, 2 \, | \vec { n } _ { i } \times \vec { n } _ { d } | \, ( \vec { e } _ { x } \cdot \vec { n } _ { 2 } ).$$

We then integrate the product of this quantity with the observable Q over the top quark decay products, and find

$$& \int \mathcal { T } _ { \lambda } [ \text{d} \sigma _ { P D } ] \ Q ( \vec { n } _ { 2 }, \{ \vec { e } \} ) \ = \frac { \alpha _ { s } C _ { F } } { 2 \pi } \ \frac { \pi \lambda } { m _ { t } } \, \text{d} \sigma _ { t } \frac { 2 | \vec { n } _ { i } \times \vec { n } _ { d } | } { \Gamma _ { t } } \ \vec { e } _ { x } \cdot \int \text{d} \Gamma _ { t } \ Q ( \vec { n } _ { 2 }, \{ \vec { e } \} ) \, \vec { n } _ { 2 } \\ & = \frac { \alpha _ { s } C _ { F } } { 2 \pi } \ \frac { \pi \lambda } { m _ { t } } \, \text{d} \sigma _ { t } \frac { 4 | \vec { n } _ { i } \times \vec { n } _ { d } | } { \Gamma _ { t } } \ \vec { e } _ { x } \cdot \int \text{d} \Gamma _ { t } \ \theta ( \vec { n } _ { 2 } \cdot \vec { e } _ { x } ) \, \vec { n } _ { 2 } = \frac { \alpha _ { s } C _ { F } } { 2 \pi } \ \frac { \pi \lambda } { m _ { t } } \, \text{d} \sigma _ { t } \, | \vec { n } _ { i } \times \vec { n } _ { d } |.$$

We use this equation to determine the non-perturbative correction to the expectation value of the observable Q

$$\langle Q \rangle = \frac { 1 } { \sigma _ { t } } \int \mathrm d \sigma _ { t } \left ( \frac { 9 } { 2 } + \frac { \alpha _ { s } C _ { F } } { 2 \pi } \ \frac { \pi \lambda } { m _ { t } } \ | \vec { n } _ { i } \times \vec { n } _ { d } | \right ),$$

at fixed center-of-mass collision energy √ s . We note that in the center-of-mass frame of partonic collision, the absolute value of the vector product of ⃗ n i and ⃗ n d reads

$$| \vec { n } _ { i } \times \vec { n } _ { d } | = \sqrt { \frac { 4 m _ { t } ^ { 2 } \ s \, t \, u } { \left ( s - m _ { t } ^ { 2 } \right ) ^ { 2 } \left ( m _ { t } ^ { 2 } - t \right ) ^ { 2 } } } = \frac { 2 m _ { t } \, s \, p _ { d \perp } } { \left ( m _ { t } ^ { 2 } - t \right ) \left ( s - m _ { t } ^ { 2 } \right ) },$$

where p d ⊥ is the transverse momentum of the d -jet relative to the collision axis. We note that in the last step we used the fact that t u = s p d 2 ⊥ .

Integrating over the scattering angle, we find

$$\frac { 1 } { \sigma _ { t } } \int \text{d} \sigma _ { t } \, | \vec { n } _ { i } \times \vec { n } _ { d } | = f _ { Q } ( s, m _ { t }, m _ { W } ),$$

where the function f Q reads

$$f _ { Q } ( s, m _ { t }, m _ { W } ) = \frac { \dot { \ } \cdot } { } \frac { \pi m _ { t } m _ { W } \sqrt { s } \sqrt { \bar { s } } \left ( - m _ { t } ^ { 4 } + m _ { t } ^ { 2 } \left ( m _ { W } ^ { 2 } + s \right ) - 2 m _ { t } m _ { W } \sqrt { s } \sqrt { \bar { s } } + m _ { W } ^ { 2 } s \right ) } { \left ( m _ { t } ^ { 2 } - m _ { W } ^ { 2 } \right ) ^ { 2 } \left ( s - m _ { t } ^ { 2 } \right ) ^ { 2 } }, \ \ ( 6. 4 0 )$$

and we have defined the quantity

$$\bar { s } = s - m _ { t } ^ { 2 } + m _ { W } ^ { 2 }.$$

The non-perturbative correction to the expectation value of the variable Q in proton collisions is obtained by convoluting the above result with parton distribution functions. In principle, since the function f Q ( s, m , m t W ) depends on the center-of-mass energy, parton distribution functions do not decouple. However, in practice, f ( s, m , m t W ) is a slowly changing function of s . Indeed, it changes from the value

$$\lim _ { s \to m _ { t } ^ { 2 } } f ( s, m _ { t }, m _ { W } ) = \frac { \pi } { 4 } \approx 0. 7 8 5,$$

at the threshold, to

$$\i o l, \, \text{to} \quad \lim _ { s \to \infty } f ( s, m _ { t }, m _ { W } ) = \frac { \pi \, m _ { t } \, m _ { W } } { ( m _ { t } + m _ { W } ) ^ { 2 } } \approx 0. 6 8,$$

at s = ∞ for physical values of m t and m W . Hence, we find the following estimate for the ATLAS variable Q in proton-proton collisions

$$\langle Q \rangle \approx \frac { 9 } { 2 } + \frac { \alpha _ { s } C _ { F } } { 2 \pi } \ \frac { \pi \lambda } { m _ { t } } \ \frac { \pi } { 4 }.$$

The above result does not account for realistic event selection criteria which in many ways introduce additional directions into the integration over top quark decay products. However, it does illustrate the point that non-perturbative effects that we discuss in this paper have a small but direct impact on the measured values of the top quark polarisation observables at hadron colliders.

An outstanding problem in collider physics is the measurement of the top quark mass with an ultrahigh precision in a credible way [23]. The tricky issue is the control (or lack of it) of non-perturbative corrections, which is very hard to do for exclusive observables that are used currently for the highest-precision measurements. In this regard, suggestions were made to study lepton observables from top quark decay because they are considered to be less prone to contaminations by non-perturbative effects. Interestingly, our analysis allows us to make exact statements to this effect, albeit in the narrow width approximation.

To this end, consider the following quantity

$$L _ { \perp } & = | \vec { p } _ { 2 } \cdot \vec { e } _ { y } |,$$

where ⃗ p 2 is the positron three-momentum in the top quark rest frame. The vector ⃗ e y is defined in eq. (6.25); it is orthogonal to the collision plane of single top production process. Hence, L ⊥ measures the component of the lepton momentum that points outside of the collision plane. We are interested in computing the average value of L ⊥ . Since

$$\int d \Gamma _ { t } \, \vec { n } _ { 2 } \, \theta ( \vec { p } _ { 2 } \cdot \vec { e } _ { y } ) \sim \vec { e } _ { y },$$

it follows that neither the term s P · s D in the leading order cross section, nor the power correction T λ [d σ PD ] receive contribution from the above integral. Since it is also easy to check that L ⊥ does not receive any corrections from momenta redefinitions, δL ⊥ = 0, it follows that

$$\langle L _ { \perp } \rangle = \frac { 1 } { \Gamma _ { t } } \int \text{d} \Gamma _ { t } \ L _ { \perp } + \mathcal { O } ( \lambda ^ { 2 } ) = \frac { 1 } { 2 \Gamma _ { t } } \int \text{d} \Gamma _ { t } \, \frac { ( p _ { 2 } p _ { t } ) } { m _ { t } } = \frac { m _ { t } ^ { 2 } + m _ { W } ^ { 2 } } { 8 m _ { t } } + \mathcal { O } ( \lambda ^ { 2 } ), \quad \text{ (6.47)}$$

where in the last step we employed the narrow width approximation for the W -boson to integrate over the positron momentum.

In summary, the point of the above calculation is to demonstrate that there are no linear power corrections to L ⊥ , so that the (short-distance) top quark mass can be determined from this observable with very high precision. An obvious reservation is that the above calculation is valid in the case when no fiducial cuts are imposed on final-state particles but it is also obvious that to perform such measurements in practice, cutting-edge simulations are required that account for perturbative and parton shower effects.

## 7 Conclusions

In this paper we have studied linear power corrections to the process of single top production followed by the top quark decay. Our primary interest is the impact of top quark instability on these corrections. Working in the narrow width approximation, we have found that linear power corrections do affect the top quark production cross section if the top quark is allowed to decay, at variance with the case of a stable top quark that was studied earlier in ref. [1].

The non-perturbative corrections that we have found in this article do affect measurements of the top quark polarisation in such processes, and also influence the kinematic distributions of leptons in top quark decays that were suggested as 'clean' observables for measuring the top quark mass. However, the particular form of power corrections, that we derived in this paper, allows us to show that 'out of the collision plane' component of the positron momentum from top quark decays does not receive linear non-perturbative corrections. Since the average value of this observable depends on m t , it is an interesting candidate for measuring the top quark mass.

Finally, the results discussed in this paper are obtained in the narrow width approximation for the top quark which corresponds to an unphysical limit Γ t ≪ Λ QCD . The next important step is to extend these result to the physical case Γ t ≫ Λ QCD . Then, the analysis becomes significantly more complicated because top quark production process and top quark decay do not factorise any more. Nevertheless, we hope that our understanding of non-perturbative power corrections to top quark production processes achieved in this paper as well as in refs [1, 2] will allow us to successfully analyse this challenging problem.

## Acknowledgments

The research of K.M. and S.M. was supported by the German Research Foundation (DFG, Deutsche Forschungsgemeinschaft) under grant 396021762-TRR 257. P.N acknowledges support by the Humboldt Foundation. M.A.O. was supported by the ANR PIA funding under grant ANR-20-IDEES-0002.

## A Alternative derivation of power corrections

The goal of this appendix is to discuss an alternative derivation of power corrections to the single top production and decay process. Here we deal directly with the amplitudes

as opposed to amplitudes squared, as was done in earlier papers [1, 2]. Below, we first discuss the Born amplitude and cross section and then continue with the real emission and virtual corrections to the top quark production and decay process. We use directly a short-distance mass scheme for the top quark mass in the calculation and we explain how to do this by considering the self-energy insertion in the (nearly) on-shell top quark line.

## The Born cross section

The Born diagram for single top production and decay is shown in Fig 1. The colour structure of this process is quite simple. We will ignore it for now and reconstruct it at the end. In the narrow width approximation, we write the Born amplitude as

$$\mathcal { B } _ { P D } = \frac { i } { p _ { t } ^ { 2 } - m _ { t } ^ { 2 } + i m _ { t } \Gamma _ { t } } \bar { u } ( p _ { f } ) N _ { D } \left [ \sum _ { \lambda _ { s } = \pm 1 } u ( p _ { t }, s, \lambda _ { s } ) \bar { u } ( p _ { t }, s, \lambda _ { s } ) \right ] N _ { P } u ( p _ { i } ), \quad \text{(A.1)}$$

where the subscripts P and D denote production and decay, and we have displayed explicitly the b i -t -b f fermion line. The functions N P ( N D ) contain all remaining structures pertinent to the production and decay processes. We assume that a quantisation axis s , satisfying the conditions s 2 = -1 and p t · s = 0, has been chosen for the top quark spin. We denote the signs of the top quark spin along the quantisation axis s with λ s . We then define

$$\mathcal { B } _ { P } ( s, \lambda _ { s } ) & = \bar { u } ( p _ { t }, s, \lambda _ { s } ) N _ { P } u ( p _ { i } ), \\ \mathcal { B } _ { D } ( s, \lambda _ { s } ) & = \bar { u } ( p _ { f } ) N _ { D } u ( p _ { t }, s, \lambda _ { s } ),$$

and write the Born cross section for single top production and decay as follows

$$\lambda \cdot \cdot \cdot \cdot \cdot \\ | A | ^ { 2 } = & \frac { 1 } { 2 m _ { t } \Gamma _ { t } } 2 \pi \delta ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) | \mathcal { B } _ { P D } | ^ { 2 }, \\ | \mathcal { B } _ { P D } | ^ { 2 } = & \sum _ { \lambda _ { s }, \lambda _ { s } ^ { \prime } } [ \mathcal { B } _ { P } ( s, \lambda _ { s } ) \mathcal { B } _ { P } ^ { * } ( s, \lambda _ { s } ^ { \prime } ) ] [ \mathcal { B } _ { D } ( s, \lambda _ { s } ) \mathcal { B } _ { D } ^ { * } ( s, \lambda _ { s } ^ { \prime } ) ],$$

where we have used the narrow width approximation, see eq. (2.7). The form of the amplitude in eq. (A.3) is the expected product of spin correlation matrices. In single top production process, a further simplification occurs since there are choices of the top quark spin quantisation axes such as B D ( s, -1) = 0 or B P ( s, -1) = 0. The fact that such quantisation axes must exist is a consequence of the fact that the helicities of all massless particles in single top production and decay are fixed by the charged-current interactions, so that also the top quark must be in a pure spin state. 10 If we call such a quantisation axis for the production process s P , we can write

$$\mathcal { B } _ { P } ( s, \lambda _ { s } ) = \bar { u } ( p _ { t }, s, \lambda _ { s } ) B = & \bar { u } ( p _ { t }, s, \lambda _ { s } ) \left [ \frac { 1 + \gamma _ { 5 } \not \, p _ { P } } { 2 } + \frac { 1 - \gamma _ { 5 } \not \, p _ { P } } { 2 } \right ] B \\ = & \bar { u } ( p _ { t }, s, \lambda _ { s } ) \frac { 1 + \gamma _ { 5 } \not \, p _ { P } } { 2 } B. \\ \frac { 1 0 \mu \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot.$$

10 This suggests that this property should also be valid in a class of single top production processes with the addition of colour-neutral particles with definite spin.

In eq. (A.4) B denotes whatever is left of B P ( s, λ s ) when the ¯ spinor is removed, and the u last step follows from the fact that

$$\bar { u } ( p _ { t }, s, \lambda _ { s } ) \left [ \frac { 1 - \gamma _ { 5 } \not \, p _ { P } } { 2 } \right ],$$

is an eigenstate of the projection of the top quark spin operator on the axis s P with an eigenstate -1 2, which by assumption does not contribute to the production process. /

Upon squaring the amplitude, we find

$$| \mathcal { B } _ { P } ( s, \lambda _ { s } ) | ^ { 2 } = & \bar { B } \frac { 1 + \gamma _ { 5 } \not \, P } { 2 } u ( p _ { t }, s, \lambda _ { s } ) \bar { u } ( p _ { t }, s, \lambda _ { s } ) \frac { 1 + \gamma _ { 5 } \not \, P } { 2 } B \\ = & \bar { B } \frac { 1 + \gamma _ { 5 } \not \, P } { 2 } \left [ ( \psi _ { t } + m _ { t } ) \frac { 1 + \lambda _ { s } \gamma _ { 5 } \not \, P } { 2 } \right ] \frac { 1 + \gamma _ { 5 } \not \, P } { 2 } B.$$

Working out the simple Dirac algebra we find

$$\frac { 1 + \gamma _ { 5 } \not \, _ { P } } { 2 } \left [ ( \phi _ { t } + m _ { t } ) \frac { 1 + \lambda _ { s } \gamma _ { 5 } \not \, _ { P } } { 2 } \right ] \frac { 1 + \gamma _ { 5 } \not \, _ { P } } { 2 } = \frac { 1 - \lambda _ { s } s \cdot s _ { P } } { 2 } ( \phi _ { t } + m _ { t } ) \frac { 1 + \gamma _ { 5 } \not \, _ { P } } { 2 }. \quad ( \text{A.7} )$$

We then insert this result into eq. (A.6) and obtain

$$| \mathcal { B } _ { P } ( s, \lambda _ { s } ) | ^ { 2 } = \frac { 1 - \lambda _ { s } s \cdot s _ { P } } { 2 } \sum _ { \pm \lambda _ { s } ^ { \prime } } \mathcal { B } _ { P } ^ { * } ( s _ { P }, \lambda _ { s } ^ { \prime } ) \mathcal { B } _ { P } ( s _ { P }, \lambda _ { s } ^ { \prime } ) = \frac { 1 - \lambda _ { s } s \cdot s _ { P } } { 2 } \, | \mathcal { B } _ { P } | ^ { 2 } \,, \quad \text{(A.8)}$$

where we have introduced the notation |B P | 2 = ∑ λ s |B P ( s, λ s ) | 2 .

A similar formula can be derived for the decay amplitude,

$$| \mathcal { B } _ { D } ( s, \lambda _ { s } ) | ^ { 2 } = \frac { 1 - \lambda _ { s } s \cdot s _ { D } } { 2 } \left | \mathcal { B } _ { D } \right | ^ { 2 },$$

where the proper quantisation axis s D differs from the one in the production. Combining the results for the production and decay amplitudes, we obtain

$$| \mathcal { B } _ { P D } | ^ { 2 } = \frac { 1 - s _ { P } \cdot s _ { D } } { 2 } \, | \mathcal { B } _ { P } | ^ { 2 } \, | \mathcal { B } _ { D } | ^ { 2 } \,.$$

This result can be derived from eq. (A.3) by choosing the quantisation axis to be either s P or s D in which case only a single term λ s = 1 contributes to sums over spin projections, and using either eq. (A.8) or (A.9). We also note that the differential decay width takes the form

$$d \Gamma _ { t } = \frac { 1 } { 4 m _ { t } } d \Phi _ { D } | \mathcal { B } _ { D } | ^ { 2 },$$

where we had to divide by two for the spin average. On the other hand our expression for the differential cross section (ignoring spin and colour averages for the initial fermions) is given by

$$\ell \lambda _ { J } \\ \ell \lambda _ { P D } = \frac { 1 } { 2 m _ { t } \Gamma _ { t } } \mathrm d \Phi _ { P } \mathrm d \Phi _ { D } | \mathcal { B } _ { D } | ^ { 2 } | \mathcal { B } _ { P } ( s _ { D } ) | ^ { 2 } = 2 \frac { \mathrm d \Gamma _ { t } } { \Gamma _ { t } } \mathrm d \sigma ( s _ { D } ),$$

which agrees with eq. (2.14). For the case of single top production depicted in Fig. 1, we can easily identify the quantisation axes s P and s D . We begin by computing s D . The decay amplitude is proportional to

$$\mathcal { B } _ { D } & \sim \bar { u } _ { f } \gamma ^ { \mu } ( 1 - \gamma _ { 5 } ) u _ { t } \, \bar { u } _ { 1 } \gamma _ { \mu } ( 1 - \gamma _ { 5 } ) v _ { 2 } = - \bar { u } _ { f } \gamma ^ { \mu } ( 1 - \gamma _ { 5 } ) u _ { t } \, \bar { u } _ { 2 R } \gamma _ { \mu } ( 1 + \gamma _ { 5 } ) v _ { 1 R } \\ & = - \left [ \bar { u } _ { 2 R } ( 1 - \gamma _ { 5 } ) u _ { t } \right ] \, \left [ \bar { u } _ { f } ( 1 + \gamma _ { 5 } ) v _ { 1 R } \right ],$$

where we have introduced the charge conjugate spinors u 2 R and v 1 R for the positron and the neutrino, and the last step uses a Fierz identity. The subscript R on the conjugate spinors is to remind that they are right-handed, i.e. ¯ u 2 ,R γ 5 = -¯ u 2 ,R . We write

$$\Pr \max \infty \colon \max \min \min \max \max \max \max \max \\ \bar { u } _ { 2, R } ( p _ { 2 } ) u ( p _ { t }, s ) = & \frac { 1 } { 2 } \bar { u } _ { 2, R } ( p _ { 2 } ) \left ( 2 + \gamma _ { 5 } \frac { m _ { t } } { ( p _ { 2 } p _ { t } ) } \left ( \vec { p } _ { 2 } - \frac { ( p _ { 2 } p _ { t } ) } { m _ { t } ^ { 2 } } ( \vec { p } _ { t } - m _ { t } ) \right ) \right ) u ( p _ { t }, s ) \\ = & \frac { 1 } { 2 } \bar { u } _ { 2, R } ( p _ { 2 } ) \left ( 1 + \gamma _ { 5 } \frac { m _ { t } } { ( p _ { 2 } p _ { t } ) } \left ( \vec { p } _ { 2 } - \frac { ( p _ { 2 } p _ { t } ) } { m _ { t } ^ { 2 } } \vec { p } _ { t } \right ) \right ) u ( p _ { t }, s ) \\ = & \bar { u } _ { 2, R } ( p _ { 2 } ) \frac { 1 + \gamma _ { 5 } \vec { p } _ { D } } { 2 } u ( p _ { t }, s ), \\ \text{where } s _ { R } \text{ reads}$$

where s D reads

$$s _ { D } ^ { \mu } & = \frac { m _ { t } } { ( p _ { 2 } p _ { t } ) } \, p _ { 2 } ^ { \mu } - \frac { 1 } { m _ { t } } \, p _ { t } ^ { \mu },$$

and satisfies the conditions s 2 D = -1 and s D · p t = 0. We note that in deriving eq. (A.14), we used the fact that the spinors ¯ u 2 ,R and u p , s ( t ) satisfy the respective Dirac equations, and that ¯ u 2 ,R γ 5 = -¯ u 2 ,R as follows from its definition. It follows from eq. (A.14) that the top quark in the decay is polarised along the axis s D , and we will refer to this quantity as the top quark spin vector in the decay.

Repeating the same calculation for the production amplitude, we easily find that top quarks are produced polarised along the quantisation axis which is given by the following equation

$$s _ { P } ^ { \mu } = \frac { m _ { t } } { ( p _ { d } p _ { t } ) } \, p _ { d } ^ { \mu } - \frac { 1 } { m _ { t } } \, p _ { t } ^ { \mu }.$$

Again, we will refer to this vector as the top quark spin vector in the production. Furthermore, in the following we will use a simplified notation, where omitting the λ s argument implies that it is taken equal to one.

## Real corrections in production

We will use the letter q rather than p to indicate momenta of particles that are affected by recoil when a soft gluon is emitted. We will also denote the top spin vector as s q , since it must be orthogonal to q t . Momentum conservation is given by

$$p _ { i } + p _ { u } = q _ { t } + q _ { d } + k.$$

The difference between q 's and p 's (and between s and s q ) are of order k , so we can change q 's into p 's and s q into s when dealing with subleading terms.

It was mentioned several times that to compute linear power corrections, we only need to consider gluon radiation off the heavy quark line. We split this contribution into two

diagrams, one that describes radiation off the final-state top quark, and the other one that describes gluon radiation off the b quark in the initial state.

We begin by computing the contribution to the amplitude of the gluon emission from the final state top quark. It is given by

$$\mathcal { R } ^ { \mu } _ { P, t } = & \bar { u } ( q _ { t }, s _ { q } ) \gamma ^ { \mu } \frac { \tilde { q } _ { t } + k + m _ { t } } { ( q _ { t } + k ) ^ { 2 } - m _ { t } ^ { 2 } } N _ { P } ( q _ { t } + k, q _ { i } ) u ( q _ { i } ) \\ = & \bar { u } ( q _ { t }, s _ { q } ) \ \frac { 2 q _ { t } ^ { \mu } + k ^ { \mu } + \sigma ^ { \mu \nu } k _ { \nu } } { ( q _ { t } + k ) ^ { 2 } - m _ { t } ^ { 2 } } \ N _ { P } ( q _ { t } + k, q _ { i } ) u ( q _ { i } ), \\. \quad. \quad. \quad. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \. \.$$

where the gluon polarisation vector has been omitted. For simplicity, we have omitted the arguments of R P, f . We should remind the reader, however, that it depends upon all the q and p momenta, and upon the spin vector s q . The arguments in N P show that this function depends on the q momenta. N P is similar to the Born diagram case, except that the incoming top quark momentum is off-shell. Nevertheless, it is a well-defined function of the external momenta. We are interested in the leading O ( k -1 ) and next-to-leading O ( k 0 ) terms in the limit of small gluon momentum k . When performing the manipulations below, we will always discard terms that vanish in the k → 0 limit.

We focus on the term σ µν k ν in eq. (A.18) acting on the ¯ spinor. u Consider the following equation

$$\bar { u } ( q _ { t }, s _ { q } ) ( 1 + a _ { \mu } \sigma ^ { \mu \nu } k _ { \nu } ) = \bar { u } _ { a },$$

where u a is defined as

$$u _ { a } \equiv ( 1 + k _ { \nu } \sigma ^ { \nu \mu } a _ { \mu } ) u ( q _ { t }, s _ { q } ),$$

and a µ is an arbitrary four-vector. We note that, up to an irrelevant phase, a general Lorentz transformation of a spinor is given by the following expression

$$\hat { S } ( \Lambda ) u ( p, s ) & = u ( \Lambda p, \Lambda s ), & ( \text{A.21} )$$

where for an infinitesimal transformation

$$\Lambda _ { \alpha \beta } = g _ { \alpha \beta } + \omega _ { \alpha \beta }, \ \omega _ { \alpha \beta } = - \omega _ { \beta \alpha },$$

the spinor transformation matrix ˆ (Λ) reads S

$$\hat { S } ( \Lambda ) = e ^ { \frac { 1 } { 4 } \omega _ { \alpha \beta } \sigma ^ { \alpha \beta } } \approx 1 + \frac { 1 } { 4 } \omega _ { \alpha \beta } \sigma ^ { \alpha \beta }.$$

If we choose

$$\omega _ { \mu \nu } = 2 ( a _ { \mu } k _ { \nu } - a _ { \nu } k _ { \mu } ),$$

we find the following equation for the spinor u a

$$u _ { a } & = u ( \Lambda q _ { t }, \Lambda s _ { q } ). & ( \text{A.25} )$$

Then, writing

$$\sigma ^ { \nu \mu } k _ { \nu } a _ { \mu } u ( q _ { t }, s _ { q } ) & = u ( \Lambda q _ { t }, \Lambda s _ { q } ) - u ( q _ { t }, s _ { q } ), & ( \Lambda. 2 6 )$$

and expanding the right-hand side in powers of k through linear terms, we find

$$\text{expanding one thing} \text{ use in power as of any way} \text{ } \text{ } \text{ } k _ { \nu } \sigma ^ { \nu \mu } a _ { \mu } u = & \frac { \partial u } { \partial p _ { t } ^ { \rho } } \Lambda ^ { \sigma \rho } p _ { t, \rho } + \frac { \partial u } { \partial s ^ { \sigma } } \Lambda ^ { \sigma \rho } s _ { \rho }, \\ & \frac { \partial u } { \partial p _ { t } ^ { \sigma } } \Lambda ^ { \sigma \rho } p _ { t, \rho } = & 2 a _ { \mu } \left ( p _ { t } ^ { \mu } k ^ { \sigma } \frac { \partial } { \partial p _ { t } ^ { \sigma } } - ( p _ { t } k ) \frac { \partial } { \partial p _ { t, \mu } } \right ) u = d _ { t } \, a _ { \mu } L _ { t } ^ { \mu } \, u, \\ & \frac { \partial u } { \partial s ^ { \sigma } } \Lambda ^ { \sigma \rho } s _ { \rho } = & 2 a _ { \mu } \left ( s ^ { \mu } k ^ { \sigma } \frac { \partial } { \partial s ^ { \sigma } } - ( s k ) \frac { \partial } { \partial s _ { \mu } } \right ) u = d _ { t } \, a _ { \mu } S _ { t } ^ { \mu } \, u, \\ \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text { } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text
 \end{cases}$$

where d t = 2( p k t ) is the denominator of the top propagator without the subleading k 2 = λ 2 term. Eqs (A.27) implicitly define the operators L t and S t as given by

$$L ^ { \mu } _ { t } = \frac { 2 } { d _ { t } } \left ( p ^ { \mu } _ { t } \, k ^ { \nu } \frac { \partial } { \partial p ^ { \nu } _ { t } } - ( p _ { t } k ) \, \frac { \partial } { \partial p _ { t, \mu } } \right ), \quad S ^ { \mu } _ { t } = \frac { 2 } { d _ { t } } \left ( s ^ { \mu } \, k ^ { \nu } \frac { \partial } { \partial s ^ { \nu } } - ( s k ) \, \frac { \partial } { \partial s _ { \mu } } \right ).$$

Since the vector a is arbitrary, from eq. (A.27) we infer the following result

$$k _ { \nu } \sigma ^ { \nu \mu } u = d _ { t } ( L _ { t } ^ { \mu } + S _ { t } ^ { \mu } ) u.$$

Next, using the definition of the current J t in eq. (3.5), and discarding terms of order k , we write the amplitude as

$$\mathcal { R } ^ { \mu } _ { P, t } = J ^ { \mu } _ { t } \bar { u } ( q _ { t }, s _ { q } ) N ( q _ { t } + k, q _ { i } ) u ( q _ { i } ) + [ ( L ^ { \mu } _ { t } + S ^ { \mu } _ { t } ) \bar { u } ( p _ { t }, s ) ] N u ( p _ { i } ). \quad \text{(A.30)}$$

We note that if the arguments of the function N are not written explicitly, it is to be understood as N p , p ( t i ). To simplify the leading term, we write

$$J _ { t } ^ { \mu } \bar { u } ( q _ { t }, s _ { q } ) N ( q _ { t } + k, q _ { i } ) u ( q _ { i } ) = J _ { t } ^ { \mu } \bar { u } ( q _ { t }, s _ { q } ) N ( q _ { t }, q _ { i } ) u ( q _ { i } ) + \frac { 2 p _ { t } ^ { \mu } } { d _ { t } } \bar { u } ( p _ { t }, s ) \ k ^ { \alpha } \frac { \partial N } { \partial p _ { t } ^ { \alpha } } u ( p _ { i } ). \ ( A. 3 1 )$$

As stated earlier, the function N q ( t + k, q i ) is a well-defined function of its arguments and can be constructed from Feynman graphs. It is not uniquely defined, however, if momentum conservation is violated, and this is exactly what happens in eq. (A.31) both in the leading term N q , q ( t i ) and when derivative with respect to p t in the last term is taken. To interpret this equation, we need to assume that N is extended in some way to account for the momentum non-conservation. The ambiguity introduced by such arbitrariness must cancel in the end, since it was not present in the initial formula. We will see later that this, in fact, is the case. We finally write

$$\mathcal { R } _ { P, f } ^ { \mu } = J _ { t } ^ { \mu } \mathcal { B } _ { P } ( s _ { q }, q ) + \bar { u } ( p _ { t }, s ) \left [ \left ( L _ { t } ^ { \mu } + \frac { \partial } { \partial p _ { t, \mu } } \right ) N \right ] u ( p _ { i } ) + [ ( L _ { t } ^ { \mu } + S _ { t } ^ { \mu } ) \bar { u } ( p _ { t }, s ) ] N u ( p _ { i } ), \ ( A. 3 2 )$$

where

$$\mathcal { B } _ { P } ( s _ { q }, q ) = \bar { u } ( q _ { t }, s _ { q } ) N ( q _ { t }, q _ { i } ) u ( q _ { i } ). & & ( \text{A.33} )$$

A similar calculation can be performed for the radiation off the b -quark in the initial state. We obtain

$$\mathcal { R } ^ { \mu } _ { P, i } = J ^ { \mu } _ { i } \mathcal { B } _ { P } ( s _ { q }, q ) + \bar { u } ( p _ { t }, s ) \left [ \left ( - L ^ { \mu } _ { i } + \frac { \partial } { \partial p _ { i, \mu } } \right ) N \right ] u ( p _ { i } ) - \bar { u } ( p _ { t }, s ) [ L ^ { \mu } _ { i } u ( p _ { i } ) ], \quad ( A. 3 4 )$$

where J i and L i are defined in eqs (3.5) and (3.8). Unlike the case of radiation off the top quark, no term analogous to the S t operator arises here, since the u p ( i ) spinor is a helicity eigenstate, and helicity is Lorentz invariant. We thus find

$$\mathcal { R } _ { P, t } + \mathcal { R } _ { P, i } = & J \, \mathcal { B } _ { P } ( s _ { q }, q ) + ( L _ { t } + S _ { t } - L _ { i } ) \left [ \bar { u } ( p _ { t }, s ) N u ( p _ { i } ) \right ] \\ + & \bar { u } ( p _ { t }, s ) \left [ \left ( \frac { \partial } { \partial p _ { t } } + \frac { \partial } { \partial p _ { i } } \right ) N \right ] u ( p _ { i } ).$$

where, as usual, J = J t + J i .

We note that all terms that appear in the first line in eq. (A.35) vanish if we multiply the equation by k µ , while the term on the second line does not. 11 On the other hand, since this term is non-singular in the soft limit, its lack of transversality must be compensated by non-singular contributions caused by the radiation from internal lines, that must have the form

$$\mathcal { R } ^ { \mu } _ { P, \text{int} } & = \bar { u } ( p _ { t }, s ) N ^ { \mu } _ { \text{reg} } u ( p _ { i } ).$$

Current conservation implies

$$[ \mathcal { R } _ { P, t } + \mathcal { R } _ { P, i } + \mathcal { R } _ { P, \text{int} } ] \cdot k = \bar { u } ( p _ { t }, s ) \, k \cdot \left [ \frac { \partial N } { \partial p _ { t } } + \frac { \partial N } { \partial p _ { i } } + N _ { \text{reg} } \right ] u ( p _ { i } ) = 0. \quad \quad \text{(A.37)}$$

In order for this equation to hold for any value of k we therefore must have

$$N _ { \text{reg} } = - \left ( \frac { \partial N } { \partial p _ { t } } + \frac { \partial N } { \partial p _ { i } } \right ).$$

Thus, the full result for gluon emission in production reads

$$\mathcal { R } _ { P } = \mathcal { R } _ { P, f } + \mathcal { R } _ { P, i } + \mathcal { R } _ { P, \text{int} } = J \, \mathcal { B } _ { P } ( s _ { q }, q ) + ( L _ { t } + S _ { t } - L _ { i } ) \, \mathcal { B } _ { P } ( s ). \quad \text{ (A.39)}$$

We now introduce the mapping from q - to p -momenta. We employ the mapping already used in ref. [1], and discussed at length near eq. (3.16). Since the mass of the top is not changed by the mapping, it must be possible to write it as a Lorentz transformation Λ which is given in eq. (3.21). In the present context, we should remember that we also need a transformation for s q , that can be conveniently chosen to be given by the same Lorentz transformation, so that the identities s 2 q = s 2 and q t · s q = p t · s hold. For convenience we report here the complete mapping transformation:

$$\cdot \quad \cdot \quad \cdot \quad & \cdot \quad \cdot \\ q _ { t } & = \Lambda p _ { t } = p _ { t } - k + \frac { ( p _ { t } k ) } { ( p _ { t } p _ { d } ) } p _ { d }, \quad q _ { d } = p _ { d } - \frac { ( p _ { t } k ) } { ( p _ { t } p _ { d } ) } p _ { d }, \\ s _ { q } & = \Lambda s = s + \frac { ( s k ) } { ( p _ { t } p _ { d } ) } p _ { d } - \frac { ( s p _ { d } ) } { ( p _ { t } p _ { d } ) } k. \\ \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot \quad \cdot$$

$$( A. 4 0 )$$

We recall that also the decay momenta must change, since the top quark momentum has changed. However, since this change is the Lorentz transformation Λ, the decay amplitude does not change. Our final result is then

$$\mathcal { R } _ { P } ^ { \mu } & = J ^ { \mu } \mathcal { B } _ { P } ( s, p ) + D _ { P, r } ^ { \mu } \mathcal { B } _ { P } ( s, p ), & \text{(A.41)}$$

11 In fact it does vanish in the single top production case. It does not necessarily vanish if we consider some associated production process, and we prefer to keep the discussion general.

where

$$D _ { P, r } & = J D _ { \text{rec} } + ( L _ { t } + S _ { t } - L _ { i } ), \\ D _ { \text{rec} } & = \left ( - k + \frac { ( p _ { t } k ) } { ( p _ { t } p _ { d } ) } p _ { d } \right ) \cdot \frac { \partial } { \partial p _ { t } } - \frac { ( p _ { t } k ) } { ( p _ { t } p _ { d } ) } p _ { d } \cdot \frac { \partial } { \partial p _ { d } } + \left ( \frac { ( k s ) } { ( p _ { t } p _ { d } ) } p _ { d } - \frac { ( p _ { d } s ) } { ( p _ { t } p _ { d } ) } k \right ) \cdot \frac { \partial } { \partial s }.$$

D rec is the differential operator associated with the momenta and spin mappings, and it can be immediately read out of eq. (A.40). It is straightforward to verify that D P,r preserves physical conditions, such as the momentum conservation, the on-shell conditions and the spin transversality condition,

$$D ^ { \mu } _ { P, r } \left ( p _ { i } + p _ { u } - p _ { t } - p _ { d } \right ) ^ { \nu } & = 0, \quad \ \ \ D ^ { \mu } _ { P, r } \, p _ { i } ^ { 2 } & = 0, \quad \ \ \ D ^ { \mu } _ { P, r } \, p _ { d } ^ { 2 } & = 0, \\ D ^ { \mu } _ { P, r } \, p _ { t } ^ { 2 } & = 0, \quad \ \ \ D ^ { \mu } _ { P, r } \, s \cdot p _ { t } & = 0, \quad \ \ \ D ^ { \mu } _ { P, r } \, s ^ { 2 } & = 0. \quad \ \ \ \ ( \Lambda. 4 3 )$$

Thus, eq. (A.41) depends upon B P ( s, p ) evaluated with momenta and spin satisfying the physical conditions, since the derivative acts in a direction tangent to the manifold where the Born amplitude is unambiguously defined.

To obtain the full amplitude for the top production and decay process, we should multiply eq. (A.41) with the decay amplitude. Since, as discussed earlier, we can assume that the momentum mapping satisfies the equation

$$\mathcal { B } _ { D } ( s _ { q }, q ) = \mathcal { B } _ { D } ( s, p ),$$

the full amplitude for the production and decay reads

$$\mathcal { R } ^ { \mu } _ { P D } = \sum _ { \lambda _ { s } = \pm 1 } \left [ J ^ { \mu } \mathcal { B } _ { P } ( s, \lambda _ { s } ) \mathcal { B } _ { D } ( s, \lambda _ { s } ) + \mathcal { B } _ { D } ( s, \lambda _ { s } ) D ^ { \mu } _ { P, r } \mathcal { B } _ { P } ( s, \lambda _ { s } ) \right ]. \quad \text{(A.45)}$$

As we explained earlier, if we choose s = s D only contribution with λ s = +1 survives in the sum. Thus, upon squaring the above formula we arrive at

$$- g _ { \mu \nu } \mathcal { R } _ { P D } ^ { \mu } \mathcal { R } _ { P D } ^ { \nu ; + } = | \mathcal { R } _ { P D } ^ { \mu } | ^ { 2 } = | \mathcal { B } _ { D } | ^ { 2 } \left \{ - J ^ { 2 } | \mathcal { B } _ { P } ( s ) | ^ { 2 } - J \cdot D _ { P ; r } | \mathcal { B } _ { P } ( s ) | ^ { 2 } \right \} _ { s = s _ { D } }. \quad \text{(A.46)}$$

The k -dependence is exposed in the above formula and, after momenta redefinitions, integration over gluon momentum factorises from the rest of phase space. The needed integrals in k can be found in ref. [1]. The result for the linear term in λ arising from the integration is given by

$$\mathcal { T } _ { \lambda } \left [ | R _ { P D } ^ { \mu } | ^ { 2 } \right ] = - | \mathcal { B } _ { P D } | ^ { 2 } \mathcal { T } _ { \lambda } \left [ \int \frac { d ^ { 3 } k } { 2 k ^ { 0 } ( 2 \pi ) ^ { 3 } } \ J ^ { 2 } \right ] + | \mathcal { B } _ { D } | ^ { 2 } \left [ \tilde { D } _ { P, r } | \mathcal { B } _ { P } ( s ) | ^ { 2 } \right ] _ { s = s _ { D } }, \quad \text{(A.47)}$$

where

$$\tilde { D } _ { P, r } = - \mathcal { T } _ { \lambda } \left [ \int \frac { \mathrm d ^ { 3 } k } { 2 k ^ { 0 } ( 2 \pi ) ^ { 3 } } \ J \cdot D _ { P, r } \right ].$$

We will not show the result of the integration of the first term in the above equation because, as we will see later, it can be combined with other contributions and argued to

cancel in a way similar to what was found in ref. [1]. Computing the required integral explicitly for the second term in eq. (A.47), we find

$$\tilde { D } _ { P, r } & = - \beta \left [ \frac { ( p _ { i } s ) ( 2 m _ { t } ^ { 2 } p _ { d } - ( p _ { t } p _ { d } ) p _ { t } ) - ( p _ { d } s ) ( 2 m _ { t } ^ { 2 } p _ { i } - ( p _ { t } p _ { i } ) p _ { t } ) } { ( p _ { t } p _ { d } ) ( p _ { t } p _ { i } ) } \cdot \frac { \partial } { \partial s } } \\ & + \left ( \frac { m _ { t } ^ { 2 } } { ( p _ { t } p _ { d } ) } p _ { d } - \frac { m _ { t } ^ { 2 } } { ( p _ { t } p _ { i } ) } p _ { i } \right ) \cdot \frac { \partial } { \partial p _ { t } } - \frac { m _ { t } ^ { 2 } } { ( p _ { t } p _ { i } ) } p _ { i } \cdot \frac { \partial } { \partial p _ { i } } - \frac { m _ { t } ^ { 2 } } { ( p _ { t } p _ { d } ) } p _ { d } \cdot \frac { \partial } { \partial p _ { d } } \right ],$$

where

$$\beta = \frac { 1 } { 2 ( 2 \pi ) ^ { 2 } } \frac { \lambda \pi } { m _ { t } }.$$

## Virtual corrections in production

We continue with the discussion of the virtual correction to the heavy line in the production subprocess. We write it as follows

$$\mathcal { V } _ { P } = \int \frac { d ^ { 4 } k } { ( 2 \pi ) ^ { 4 } } \frac { - i } { k ^ { 2 } - \lambda ^ { 2 } + i \epsilon } \ F _ { V P } ( k, \dots ),$$

where

$$F _ { V P } & = \left [ \bar { u } ( p _ { t }, s ) \gamma ^ { \mu } \frac { \vec { p } _ { t } + \vec { k } + m _ { t } } { ( p _ { t } + k ) ^ { 2 } - m _ { t } ^ { 2 } + i \epsilon } N ( p _ { t } + k, p _ { i } + k ) \frac { \vec { p } _ { i } + \vec { k } } { ( p _ { i } + k ) ^ { 2 } + i \epsilon } \gamma _ { \mu } u ( p _ { i } ) \\ & + \bar { u } ( p _ { t }, s ) \gamma ^ { \mu } \frac { \vec { p } _ { t } + \vec { k } + m _ { t } } { ( p _ { t } + k ) ^ { 2 } - m _ { t } ^ { 2 } + i \epsilon } N _ { \mu } u ( p _ { i } ) + \bar { u } ( p _ { t }, s ) N _ { \mu } \frac { \vec { p } _ { i } + \vec { k } } { ( p _ { i } + k ) ^ { 2 } + i \epsilon } \gamma ^ { \mu } u ( p _ { i } ) \right ].$$

The first line provides a contribution where a virtual gluon is emitted by an incoming bottom and absorbed by the outgoing top quark, and the terms in the second line describe contributions where virtual gluons are emitted by either bottom or top quarks and are absorbed by the internal lines of the diagrams. Potential contributions where gluons are emitted and absorbed by internal lines are not shown as they cannot produce O ( λ ) corrections [1]. Using the Dirac equations, and neglecting contributions that cannot produce O ( λ ) corrections, we rewrite the above expression as follows

$$F _ { V P } & = \bar { u } ( p _ { t }, s ) \frac { 2 p _ { t } ^ { \mu } + k ^ { \mu } + \sigma ^ { \mu \nu } k _ { \nu } } { ( p _ { t } + k ) ^ { 2 } - m _ { t } ^ { 2 } + i \epsilon } N ( p _ { t } + k, p _ { i } + k ) \frac { 2 p _ { t, \mu } + k _ { \mu } - \sigma _ { \mu \rho } k ^ { \rho } } { ( p _ { i } + k ) ^ { 2 } + i \epsilon } u ( p _ { i } ) \\ & + \bar { u } ( p _ { t }, s ) \frac { 2 p _ { t } ^ { \mu } } { ( p _ { t } + k ) ^ { 2 } - m _ { t } ^ { 2 } + i \epsilon } N _ { \mu } u ( p _ { i } ) + \bar { u } ( p _ { t }, s ) N _ { \mu } \frac { 2 p _ { i } ^ { \mu } } { ( p _ { i } + k ) ^ { 2 } + i \epsilon } u ( p _ { i } ) \,. \\ \quad \,. \quad \,. \quad \,.$$

Using again eq. (A.38) and eq. (A.29) we arrive at

$$\mathcal { V } _ { P } = & \int \frac { d ^ { 4 } k } { ( 2 \pi ) ^ { 4 } } \frac { - i } { k ^ { 2 } - \lambda ^ { 2 } + i \epsilon } \left \{ J _ { t } \cdot J _ { i } \, \mathcal { B } _ { P } ( s ) + \left [ ( J _ { i } \cdot ( L _ { t } + S _ { t } ) + J _ { t } \cdot L _ { i } ) \mathcal { B } _ { P } ( s ) \right ] \\ - & \bar { u } ( p _ { t }, s ) \left ( J _ { t } \cdot \frac { \partial N } { \partial p _ { t } } + J _ { i } \cdot \frac { \partial N } { \partial p _ { i } } \right ) u ( p _ { i } ) \right \},$$

where the currents J t and J i have now changed appropriately for eq. (A.52), and the definitions of L t , S t and L i are given in eqs (A.28) and (3.8) except that the newly defined denominators d , d t i should be used there.

The second term on the right hand side of eq. (A.54) is not a total derivative. To remedy this, we assume that spinors can be written as functions of their momentum alone. 12 To do this, in the expression for u p , s, λ ( t s ) we systematically replace the mass m t with √ p 2 t . This implies the following modification in the density matrix

$$\sum _ { \lambda _ { s } = \pm 1 } u ( p _ { t }, s, \lambda _ { s } ) \bar { u } ( p _ { t }, s, \lambda _ { s } ) = ( \phi _ { t } + \sqrt { p _ { t } ^ { 2 } } ) \,.$$

The replacement m t → √ p 2 t is important for simplifying eq. (A.54), because it implies

$$p _ { t } \cdot \frac { \partial } { \partial p _ { t } } u ( p _ { t }. s ) = \frac { 1 } { 2 } u ( p _ { t }. s ).$$

This result easily follows from the fact that the mass dimension of a spinor is 1 / 2 and that once the mass is eliminated in favour of √ p 2 t , p t becomes the only mass scale that appears in the formula for the spinor. Hence, we find

$$& \bar { u } ( p _ { t }, s ) \left ( J _ { t } \cdot \frac { \partial N } { \partial p _ { t } } + J _ { i } \cdot \frac { \partial N } { \partial p _ { i } } \right ) u ( p _ { i } ) \\ & = \left ( J _ { i } \cdot \frac { \partial } { \partial p _ { i } } + J _ { t } \cdot \frac { \partial } { \partial p _ { t } } \right ) [ \bar { u } ( p _ { t }, s ) N u ( p _ { i } ) ] - \left ( \frac { 1 } { d _ { i } } + \frac { 1 } { d _ { t } } \right ) \bar { u } ( p _ { t }, s ) N u ( p _ { i } ).$$

Defining

$$D _ { P, v } = J _ { i } \cdot ( L _ { t } + S _ { t } ) + J _ { t } \cdot L _ { i } - J _ { i } \cdot \frac { \partial } { \partial p _ { i } } - J _ { t } \cdot \frac { \partial } { \partial p _ { t } },$$

we write eq. (A.54) as

$$\mathcal { V } _ { P } = \int \frac { \mathrm d ^ { 4 } k } { ( 2 \pi ) ^ { 4 } } \frac { - i } { k ^ { 2 } - \lambda ^ { 2 } + i \epsilon } \left [ \left ( J _ { t } \cdot J _ { i } + \frac { 1 } { d _ { i } } + \frac { 1 } { d _ { t } } \right ) \mathcal { B } _ { P } ( s ) + D _ { P, v } \mathcal { B } _ { P } ( s ) \right ].$$

We note that the derivative D P,v violates the physicality constraint related to momentum conservation; we will see that it is restored once the mass renormalisation is accounted for.

Since the dependence on the gluon momentum k is exposed in eq. (A.59), we can integrate over it. Similar to the discussion of the real-emission contribution, we leave terms proportional to B P ( s ) as they are, since we will argue later that their cancellation is already demonstrated in ref. [1]. We obtain

$$\mathcal { T } _ { \Lambda } [ \mathcal { V } _ { P } ] = \mathcal { T } _ { \lambda } \left [ \int \frac { d ^ { 4 } k } { ( 2 \pi ) ^ { 4 } } \frac { - i } { k ^ { 2 } - \lambda ^ { 2 } + i \epsilon } \left ( J _ { t } \cdot J _ { i } + \frac { 1 } { d _ { i } } + \frac { 1 } { d _ { t } } \right ) \right ] \mathcal { B } _ { P } ( s ) + \tilde { D } _ { P, v } \mathcal { B } _ { P } ( s ), \\ \frac { \tilde { D } _ { P, v } } { ( p _ { t } p _ { i } ) } = \frac { \beta } { ( p _ { i } s ) } \left ( - ( p _ { t } s ) \, p _ { t } \cdot \frac { \partial } { \partial s } + m _ { t } ^ { 2 } \, p _ { i } \cdot \frac { \partial } { \partial p _ { i } } + m _ { t } ^ { 2 } \, p _ { i } \cdot \frac { \partial } { \partial p _ { t } } \right ).$$

12 This choice does not affect the L t + S t and the L i derivatives, since they act as Lorentz transformations on the argument of the spinor, leaving the mass and the constraints on the spin parameter s unchanged.

We then multiply eq. (A.60) by the decay amplitude, set s = s D to get rid of the spin summation, compute the interference with the Born amplitude, and finally obtain

$$2 | \mathcal { B } _ { P D } | ^ { 2 } \mathcal { T } _ { \lambda } & \left [ \int \frac { d ^ { 4 } k } { ( 2 \pi ) ^ { 4 } } \frac { - i } { k ^ { 2 } - \lambda ^ { 2 } + i \epsilon } \left ( J _ { t } \cdot J _ { i } + \frac { 1 } { d _ { i } } + \frac { 1 } { d _ { t } } \right ) \right ] \\ & + | \mathcal { B } _ { D } | ^ { 2 } \left [ \tilde { D } _ { P, v } | \mathcal { B } _ { P } ( s ) | ^ { 2 } \right ] _ { s = s _ { D } }.$$

## Real corrections in decay

We can describe radiation in the top quark decay following the approach used in the discussion of the production process. The only essential difference is in the momentum mapping, that we choose to coincide with the one discussed in ref. [1]. Details are given in eq. (4.4). Both the initial top quark momentum and its spin vector are not affected by the mapping. The decay amplitude expanded through linear terms in the gluon momentum reads

$$\mathcal { R } _ { D } ^ { \mu } ( s ) = J ^ { \mu } \mathcal { B } _ { D } ( s ) + D _ { D, r } ^ { \mu } \mathcal { B } _ { D } ( s ), \quad D _ { D, r } ^ { \mu } = J ^ { \mu } D _ { \text{rec} } - ( L _ { t } ^ { \mu } + S _ { t } ^ { \mu } - L _ { f } ^ { \mu } ), \quad ( A. 6 2 )$$

where the currents are given in eq. (4.3). The differential operator associated with the mapping can be immediately read out of eq. (4.4). It reads

$$D _ { r e c } = \left ( - k + \frac { ( p _ { f } k ) } { ( p _ { 2 } p _ { f } ) } p _ { 2 } \right ) \cdot \frac { \partial } { \partial p _ { f } } - \frac { ( p _ { f } k ) } { ( p _ { 2 } p _ { f } ) } p _ { 2 } \cdot \frac { \partial } { \partial p _ { 2 } }.$$

The operators L t and S t are defined in eq. (A.27) but now d t = -2( p k t ) has to be used there. The definition of L f is the same as the one for L i after the replacement i → f is performed. The operator D D,r is easily seen to preserve the physicality conditions for the decay. Proceeding as for eq. (A.46), we get

$$| \mathcal { R } _ { P D } ^ { \mu } | ^ { 2 } = | \mathcal { B } _ { P } | ^ { 2 } \left \{ - J ^ { 2 } | \mathcal { B } _ { D } ( s ) | ^ { 2 } - J \cdot D _ { D, r } | \mathcal { B } _ { D } ( s ) | ^ { 2 } \right \} _ { s = s _ { P } }.$$

After the momentum mapping, the integration over k factorises and can be performed. Defining

$$\tilde { D } _ { D, r } = & \mathcal { T } _ { \lambda } \left [ \int \frac { \mathrm d ^ { 3 } k } { 2 k ^ { 9 } ( 2 \pi ) ^ { 4 } } ( - J \cdot D _ { D, r } ) \right ] = \beta \left [ \frac { ( p _ { f } s ) } { ( p _ { t } p _ { f } ) } p _ { t } \cdot \frac { \partial } { \partial s } \\ & + \frac { ( p _ { t } p _ { f } ) } { ( p _ { 2 } p _ { f } ) } p _ { 2 } \cdot \left ( \frac { \partial } { \partial p _ { 2 } } - \frac { \partial } { \partial p _ { f } } \right ) + \left ( p _ { t } - \frac { m _ { t } ^ { 2 } } { ( p _ { t } p _ { f } ) } p _ { f } \right ) \cdot \left ( \frac { \partial } { \partial p _ { f } } + \frac { \partial } { \partial p _ { t } } \right ) \right ],$$

we find for the total amplitude squared

$$\mathcal { T } _ { \lambda } \left [ | \mathcal { R } _ { P D } | ^ { 2 } \right ] = | \mathcal { B } _ { P D } | ^ { 2 } \mathcal { T } _ { \lambda } \left [ \int \frac { \mathfrak { d } ^ { 3 } k } { 2 k ^ { 0 } ( 2 \pi ) ^ { 4 } } ( - J ^ { 2 } ) \right ] + | \mathcal { B } _ { P } | ^ { 2 } \left [ \tilde { D } _ { D, r } | \mathcal { B } _ { D } ( s ) | ^ { 2 } \right ] _ { s = s _ { P } }. \quad ( \text{A.66} )$$

Figure 2 : The graph with the top quark self-energy insertion. The complex conjugate diagram should also be added.

![image](https://ac.marcodsn.me/arxiv-markdown-images/2408.00632v2-figure-c4cb2f60-94d1-4ee6-8a67-95181a49fdad.jpg)

## Virtual corrections in decay

The calculation of the virtual corrections in the decay process closely follows the production case, and the result can be obtained from eq. (A.60) after the substitution i → f . We obtain

$$\mathcal { T } _ { \lambda } [ \mathcal { V } _ { D } ] & = \mathcal { T } _ { \lambda } \left [ \int \frac { d ^ { 4 } k } { ( 2 \pi ) ^ { 4 } } \frac { - i } { k ^ { 2 } - \lambda ^ { 2 } + i \epsilon } \left ( J _ { t } \cdot J _ { f } + \frac { 1 } { d _ { f } } + \frac { 1 } { d _ { t } } \right ) \right ] \mathcal { B } _ { D } ( s ) + \tilde { D } _ { D, v } \mathcal { B } _ { D } ( s ), \\ \tilde { D } _ { D, v } & = \frac { \beta } { ( p _ { t } p _ { f } ) } \left ( - ( p _ { f } s ) \, p _ { t } \cdot \frac { \partial } { \partial s } + m _ { t } ^ { 2 } \, p _ { f } \cdot \frac { \partial } { \partial p _ { f } } + m _ { t } ^ { 2 } \, p _ { f } \cdot \frac { \partial } { \partial p _ { t } } \right ). \\ \dots \quad \lambda \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

The final contribution is

$$2 | \mathcal { B } _ { P D } | ^ { 2 } \mathcal { T } _ { \lambda } & \left [ \int \frac { \mathrm d ^ { 4 } k } { ( 2 \pi ) ^ { 4 } } \frac { - i } { k ^ { 2 } - \lambda ^ { 2 } + i \epsilon } \left ( J _ { t } \cdot J _ { f } + \frac { 1 } { d _ { f } } + \frac { 1 } { d _ { t } } \right ) \right ] \\ & \quad + | \mathcal { B } _ { P } | ^ { 2 } \left [ \tilde { D } _ { D, v } | \mathcal { B } _ { D } ( s ) | ^ { 2 } \right ] _ { s = s _ { P } }. \\ \cdot \quad \cdot \quad \cdot \quad \cdot$$

## Top quark self-energy contribution

Finally, we need to account for the top quark self-energy correction shown in Fig. 2. Since we perform the calculation in a short-distance scheme (e.g. the MS scheme), no O ( λ ) correction can arise from the mass counter-terms, and we do not need to account for them. Following the discussion in section 7 in ref. [1], we find that the term of order λ from the self-energy insertion to the left of the cut line in Fig. 2 reads

$$\Sigma = \kappa \left [ - \frac { 1 } { 4 m _ { t } } ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) - ( \psi _ { t } - m _ { t } ) + m _ { t } \right ].$$

where

$$\kappa = 4 \pi \alpha _ { s } C _ { F } \beta = \frac { \alpha _ { s } C _ { F } } { 2 \pi } \, \frac { \pi \lambda } { m _ { t } }. \quad \quad \quad \text{(A.70)}$$

The product of the denominators of the top propagators in Fig. 2 yields

$$\frac { 1 } { ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) ^ { 2 } + ( m _ { t } \Gamma _ { t } ) ^ { 2 } } \frac { 1 } { ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) + i m _ { t } \Gamma _ { t } }.$$

When this expression is combined with the complex conjugate of Fig. 2, it becomes

$$\frac { 2 ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) } { ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) ^ { 2 } + ( m _ { t } \Gamma _ { t } ) ^ { 2 } ) ^ { 2 } } = - \frac { \partial } { \partial p _ { t } ^ { 2 } } \frac { 1 } { ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) ^ { 2 } + ( m _ { t } \Gamma _ { t } ) ^ { 2 } } \approx - \frac { \pi } { m _ { t } \Gamma _ { t } } \delta ^ { \prime } ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ), \quad ( A. 7 2 )$$

where δ ′ ( p 2 t -m 2 t ) is the derivative of a delta function with respect to p 2 t . Thus, the net effect of the self-energy correction amounts to the following replacement in the Born cross section

$$i ( \tilde { \varphi } _ { t } + m _ { t } ) \delta ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) \to - i ( \tilde { \varphi } _ { t } + m _ { t } ) i \Sigma i ( \tilde { \varphi } _ { t } + m _ { t } ) \delta ^ { \prime } ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) \\ = i \kappa \left [ \frac { 3 } { 2 } ( \tilde { \varphi } _ { t } + m _ { t } ) - m _ { t } \right ] \delta ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) + i \kappa ( \tilde { \varphi } _ { t } + m _ { t } ) 2 m _ { t } ^ { 2 } \delta ^ { \prime } ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ),$$

where the relation δ ′ ( x x ) = -δ x ( ) was used.

In the self-energy computation, one needs to consider slightly off-shell momenta of the top quark. This means that we cannot set p 2 t = m 2 t in the factor involving the derivative of the delta function. We thus rewrite this term as follows

$$\psi _ { t } + m _ { t } = ( \psi _ { t } + \sqrt { p _ { t } ^ { 2 } } ) + \frac { m _ { t } ^ { 2 } - p _ { t } ^ { 2 } } { m _ { t } + \sqrt { p _ { t } ^ { 2 } } }.$$

Inserting this identity in eq. (A.73) we obtain

$$( \varphi _ { t } + m _ { t } ) \delta ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) \to \kappa _ { \frac { 3 } { 2 } } ( \varphi _ { t } + m _ { t } ) \delta ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) + \kappa ( \varphi _ { t } + \sqrt { p _ { t } ^ { 2 } } ) 2 m _ { t } ^ { 2 } \delta ^ { \prime } ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ). \quad ( A. 7 5 )$$

The first term in the above equation leads to a term proportional to |B PD | 2 , and it can be set aside with all other terms of this form. Using for consistency eq. (A.55), we obtain for the second term

$$| \mathcal { B } _ { P D } | ^ { 2 } ( 2 m _ { t } ^ { 2 } \kappa \delta ^ { \prime } ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) ) = | \mathcal { B } _ { P D } | ^ { 2 } \left [ \delta \left ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } \left ( 1 - 2 \kappa \right ) \right ) - \delta ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } ) \right ]. \quad \text{(A.76)}$$

In order to compute the result we need to perform a change of variables in the first term. We first relabel all the p 's into q 's and s into s q . In particular

$$\delta \left ( p _ { t } ^ { 2 } - m _ { t } ^ { 2 } \left ( 1 - 2 \kappa \right ) \right ) \rightarrow \delta \left ( q _ { t } ^ { 2 } - m _ { t } ^ { 2 } \left ( 1 - 2 \kappa \right ) \right ).$$

Then we recall that in the unpolarised case [1] two different mappings were used for singletop production process and for top decay. The momenta and spin transformations for the decay are given in eq. (5.1), and we report them here for convenience

$$q _ { t } = p _ { t } ( 1 - \kappa ), \ \ q _ { f } = p _ { f } - \kappa p _ { t } + \kappa \frac { ( p _ { f } p ) } { ( p _ { f } p _ { 2 } ) } p _ { 2 }, \ \ q _ { 2 } = p _ { 2 } \left ( 1 - \kappa \frac { ( p _ { f } p _ { t } ) } { ( p _ { f } p _ { 2 } ) } \right ), \ \ s _ { q } = s.$$

Notice that the spin is unchanged, since the top momentum is only rescaled. For the production we use

$$q _ { t } = p _ { t } - \kappa \frac { p _ { t } ^ { 2 } } { ( p _ { d } p _ { t } ) } p _ { d }, \quad q _ { d } = p _ { d } \left ( 1 + \kappa \frac { p _ { t } ^ { 2 } } { ( p _ { d } p _ { t } ) } \right ), \quad s _ { q } = \Lambda _ { m } s,$$

where Λ m = 1 + κω td (see eq. (5.11)) and ω td is given in eq. (3.30). The transformation for s q can be seen to satisfy the condition s q · q t = s · p t and s 2 q = s 2 , as will become clear in the following. We should modify these transformations in such a way that the top quark momentum transforms in the same way in both production and decay processes. This is achieved by using the fact that Λ m is a Lorentz transformation, and that the top momentum mapping can be written as the product of a rescaling times Λ m . It is sufficient then to apply Λ m also to all decay products. In summary, we write

$$q _ { t } & = ( 1 - \kappa ) \Lambda _ { m } p _ { t }, \quad q _ { d } = p _ { d } \left ( 1 + \kappa \frac { p _ { t } ^ { 2 } } { ( p _ { d } p _ { t } ) } \right ), \quad s _ { q } = \Lambda _ { m } s, \\ q _ { f } & = \Lambda _ { m } \left ( p _ { f } - \kappa p _ { t } + \kappa \frac { ( p _ { f } p _ { t } ) } { ( p _ { f } p _ { 2 } ) } p _ { 2 } \right ), \quad q _ { 2 } = \left ( 1 - \kappa \frac { ( p _ { f } p _ { t } ) } { ( p _ { f } p _ { 2 } ) } \right ) \Lambda _ { m } p _ { 2 }, \quad q _ { 1 } = \Lambda _ { m } p _ { 1 },$$

where only O ( κ ) terms need to be retained on the right hand sides of the above equations. Notice that the first line is just a rewriting of eq. (A.79), while the second line is the Lorentz transformation applied to the eq. (A.78). Momenta modifications induce changes in the squared amplitudes. Depending on whether momenta in the production or decay undergo these transformations, we group such terms into modification of production and decay amplitudes and write

$$D ^ { m } \left | \sum _ { \pm 1 } \mathcal { B } _ { D } ( \lambda _ { s } ) \mathcal { B } _ { P } ( \lambda _ { s } ) \right | ^ { 2 } = | \mathcal { B } _ { P } | ^ { 2 } [ D ^ { m } | \mathcal { B } _ { D } ( s ) | ^ { 2 } ] _ { s = s _ { P } } + | \mathcal { B } _ { D } | ^ { 2 } [ D ^ { m } | \mathcal { B } _ { P } ( s ) | ^ { 2 } ] _ { s = s _ { D } }. \quad ( A. 8 1 )$$

The differential operator D m is associated with the transformations shown in eq. (A.80). It can be written as

$$D ^ { m } | \mathcal { B } _ { D } ( s ) | ^ { 2 } & = \kappa \left [ - p _ { t } \cdot \frac { \partial } { \partial p _ { t } } + \left ( - p _ { t } + \frac { ( p _ { f } p _ { t } ) } { ( p _ { f } p _ { 2 } ) } p _ { 2 } \right ) \cdot \frac { \partial } { \partial p _ { f } } - \frac { ( p _ { f } p _ { t } ) } { ( p _ { f } p _ { 2 } ) } p _ { 2 } \cdot \frac { \partial } { \partial p _ { 2 } } \right ] | \mathcal { B } _ { D } ( s ) | ^ { 2 }, \ \left ( \text{A.8} 2 \right ) \\ D ^ { m } | \mathcal { B } _ { P } ( s ) | ^ { 2 } & = \kappa \left [ ( - p _ { t } + \omega _ { t d } p _ { t } ) \cdot \frac { \partial } { \partial p _ { t } } + \frac { p _ { t } ^ { 2 } } { ( p _ { d } p _ { t } ) } p _ { d } \cdot \frac { \partial } { \partial p _ { d } } + ( \omega _ { t d } s ) \cdot \frac { \partial } { \partial s } \right ] | \mathcal { B } _ { P } ( s ) | ^ { 2 },$$

where we have simply dropped from D m the derivatives with respect to variables not contained in the corresponding amplitude and, in the case of the derivative of the decay amplitude, we have removed the Lorentz transformation, since it affects all the decay momenta and the spin vector s , and the decay amplitude is Lorentz invariant.

## Assembling everything

We begin by considering the non-derivative terms that appear in eqs (A.47, A.61, A.66, A.68, A.82, A.83). They arise from the dominant terms in the cross sections, from Jacobians due to momenta transformations and in the calculation of the virtual contributions in production and decay. Their calculation is straightforward, and can be carried out along the lines of ref. [1], where it has been shown that they cancel.

We continue with terms that contain derivatives with respect to momenta and spins of external particles. Combining such contributions in the virtual corrections for production, eq. (A.61), the real emission contributions in production, eq. (A.47), and the corresponding

part of the self-energy correction, eq. (A.83), we find

$$| \mathcal { B } _ { D } | ^ { 2 } \left [ \left ( D ^ { m } + \frac { \kappa } { \beta } [ \tilde { D } _ { P, v } + \tilde { D } _ { P, r } ] \right ) | \mathcal { B } _ { P } ( s ) | ^ { 2 } \right ] _ { s = s _ { D } }.$$

It is a matter of simple algebra to verify that the derivative operator D m + κ/βD ˜ P,v , with D m restricted to the production variables, and, independently, ˜ D P,r preserve all physicality conditions, c.f. eq. (A.43). The result shown in eq. (A.84) becomes

$$| \mathcal { B } _ { D } | ^ { 2 } \kappa \left [ \left \{ s \cdot \left ( 2 \omega _ { t i } + 2 \omega _ { d t } + \frac { 2 m _ { t } ^ { 2 } ( p _ { i } p _ { d } ) } { ( p _ { t } p _ { i } ) ( p _ { t } p _ { d } ) } \omega _ { i d } \right ) \cdot \frac { \partial } { \partial s } \right \} | \mathcal { B } _ { P } ( s ) | ^ { 2 } \right ] _ { s = s _ { D } }.$$

For the decay contribution, assembling eqs (A.66), (A.68) and (A.82) we get

$$| \mathcal { B } _ { P } | ^ { 2 } \left [ \left ( D ^ { m } + \frac { \kappa } { \beta } [ \tilde { D } _ { D, v } + \tilde { D } _ { D, r } ] \right ) | \mathcal { B } _ { D } ( s ) | ^ { 2 } \right ] _ { s = s _ { P } } = 0.$$

The sum of eqs (A.85) and (A.86) agrees with eq. (6.1).

## References

- [1] S. Makarov, K. Melnikov, P. Nason and M.A. Ozcelik, Linear power corrections to single top production processes at the LHC , JHEP 05 (2023) 153 [ 2302.02729 ].
- [2] S. Makarov, K. Melnikov, P. Nason and M.A. Ozcelik, Linear power corrections to top quark pair production in hadron collisions , JHEP 01 (2024) 074 [ 2308.05526 ].
- [3] M. Beneke, Renormalons , Phys. Rept. 317 (1999) 1 [ hep-ph/9807443 ].
- [4] S. Ferrario Ravasio, P. Nason and C. Oleari, All-orders behaviour and renormalons in top-mass observables , JHEP 01 (2019) 203 [ 1810.10931 ].
- [5] V.S. Fadin, V.A. Khoze and A.D. Martin, How suppressed are the radiative interference effects in heavy instable particle production? , Phys. Lett. B 320 (1994) 141 [ hep-ph/9309234 ].
- [6] K. Melnikov and O.I. Yakovlev, Top near threshold: All alpha-S corrections are trivial , Phys. Lett. B 324 (1994) 217 [ hep-ph/9302311 ].
- [7] E.L. Berger, J. Gao, C.P. Yuan and H.X. Zhu, NNLO QCD Corrections to t-channel Single Top-Quark Production and Decay , Phys. Rev. D 94 (2016) 071501 [ 1606.08463 ].
- [8] J. Campbell, T. Neumann and Z. Sullivan, Single-top-quark production in the t -channel at NNLO JHEP , 02 (2021) 040 [ 2012.01574 ].
- [9] G.L. Kane, G.A. Ladinsky and C.P. Yuan, Using the Top Quark for Testing Standard Model Polarization and CP Predictions , Phys. Rev. D 45 (1992) 124.
- [10] J.A. Aguilar-Saavedra and S. Amor Dos Santos, New directions for top quark polarization in the t -channel process , Phys. Rev. D 89 (2014) 114009 [ 1404.1585 ].
- [11] G. Mahlon and S.J. Parke, Improved spin basis for angular correlation studies in single top quark production at the Tevatron , Phys. Rev. D 55 (1997) 7249 [ hep-ph/9611367 ].
- [12] G. Mahlon and S.J. Parke, Single top quark production at the LHC: Understanding spin , Phys. Lett. B 476 (2000) 323 [ hep-ph/9912458 ].

- [13] M. Baumgart and B. Tweedie, A New Twist on Top Quark Spin Correlations , JHEP 03 (2013) 117 [ 1212.4888 ].
- [14] CMS collaboration, Measurement of top quark polarization in t-channel single-top production , CMS-PAS-TOP-13-001 (2013) .
- [15] CMS collaboration, Measurement of Top-Quark Polarization in t-channel Single-Top Production , in 7th International Workshop on Top Quark Physics , 11, 2014 [ 1412.0106 ].
- [16] ATLAS collaboration, Measurement of the polarisation of single top quarks and antiquarks produced in the t-channel at √ s = 13 TeV and bounds on the tWb dipole operator from the ATLAS experiment , JHEP 11 (2022) 040 [ 2202.11382 ].
- [17] I.I.Y. Bigi, M.A. Shifman, N.G. Uraltsev and A.I. Vainshtein, The Pole mass of the heavy quark. Perturbation theory and beyond , Phys. Rev. D 50 (1994) 2234 [ hep-ph/9402360 ].
- [18] M. Beneke and V.M. Braun, Heavy quark effective theory beyond perturbation theory: Renormalons, the pole mass and the residual mass term , Nucl. Phys. B 426 (1994) 301 [ hep-ph/9402364 ].
- [19] F. Caola, S. Ferrario Ravasio, G. Limatola, K. Melnikov and P. Nason, On linear power corrections in certain collider observables , JHEP 01 (2022) 093 [ 2108.08897 ].
- [20] F.E. Low, Bremsstrahlung of very low-energy quanta in elementary particle collisions , Phys. Rev. 110 (1958) 974.
- [21] T.H. Burnett and N.M. Kroll, Extension of the low soft photon theorem , Phys. Rev. Lett. 20 (1968) 86.
- [22] H.W. Fearing, Note on the generalization of the burnett-kroll soft-photon theorem to polarized cross-sections , Phys. Rev. D 7 (1973) 243.
- [23] P. Azzi et al., Report from Working Group 1: Standard Model Physics at the HL-LHC and HE-LHC CERN Yellow Rep. Monogr. , 7 (2019) 1 [ 1902.04070 ].
Sergei Makarov Kirill Melnikov Paolo Nason Melih A. Ozcelik
Linear power corrections to single top production and decay at the LHC in the narrow width approximation
We consider top quark production and decay in the narrow width approximation
and study if the polarisation effects, that manifest themselves in correlations
of angular distributions of particles from top quark decays and final state
jets in the production sub-process, are affected by linear power corrections.
We find that, in general, the answer to this question is affirmative. We also
discuss how these non-perturbative corrections affect polarisation observables
used to study single top production at the LHC. Finally, we point out that
generic kinematic distributions of leptons from top quark decays are affected
by linear power corrections, which may have implications for proposals to
extract the top quark mass from such leptonic observables. On the other hand,
we demonstrate that the distribution of the ``out-of-collision-plane''
component of the positron momentum is free from linear power corrections,
making it an interesting candidate for the top quark mass measurement.