Angular distribution of polarised Λb baryons decaying to Λℓ+ℓ−

Rare b → sℓ+ℓ− flavour-changing-neutral-current processes provide important tests of the Standard Model of particle physics. Angular observables in exclusive b → sℓ+ℓ− processes can be particularly powerful as they allow hadronic uncertainties to be controlled. Amongst the exclusive processes that have been studied by experiments, the decay Λb → Λℓ+ℓ− is unique in that the Λb baryon can be produced polarised. In this paper, we derive an expression for the angular distribution of the Λb → Λℓ+ℓ− decay for the case where the Λb baryon is produced polarised. This extends the number of angular observables in this decay from 10 to 34. Standard Model expectations for the new observables are provided and the sensitivity of the observables is explored under a variety of new physics models. At low-hadronic recoil, four of the new observables have a new short distance dependence that is absent in the unpolarised case. The remaining observables depend on the same short distance contributions as the unpolarised observables, but with different dependence on hadronic form-factors. These relations provide possibilities for novel tests of the SM that could be carried out with the data that will become available at the LHC or a future e+e− collider.


Introduction
Rare b → s + − have been studied extensively by experiments at the B-factories as well as experiments at the Tevatron and Large Hadron Collider (LHC). Amongst the b → s + − processes that have been studied, the decay Λ b → Λµ + µ − is unique for two reasons: it is the only baryonic decay that has been studied; and the Λ baryon decays weakly leading to new hadron-side observables. The angular distribution of Λ b → Λµ + µ − decays has been studied in refs. [1,2] for the case of unpolarised Λ b baryons. The resulting angular distribution is described by 10 angular observables. The decay rate and lepton side angular distribution has also been studied in the SM and in several extensions of the SM (NP models) in refs. [3][4][5][6][7][8][9][10][11][12]. If the Λ b is produced polarised, a much larger number of observables are measurable.
In e + e − collisions, Λ b baryons can be produced with large longitudinal polarisations. The longitudinal polarisation of Λ b baryons and b-quarks produced via e + e − → Z 0 (→ bb) decays has been studied by the LEP experiments in refs. [17][18][19]. The production of Λ b baryons with longitudinal polarisation is forbidden in strong interactions, due to parity conservation. The Λ b can, however, be produced with transverse polarisation in pp collisions. In this paper, we focus on the transverse polarisation of the Λ b baryon. The transverse polarisation of Λ b baryons produced in pp collisions at √ s = 7 and 8 TeV has been studied by the LHCb and CMS experiments in refs. [20] and [21], respectively. The LHCb experiment measures P Λ b = 0.06 ± 0.07 ± 0.02 at √ s = 7 TeV. The CMS experiment measures P Λ b = 0.00 ± 0.06 ± 0.02 combining data from √ s = 7 and 8 TeV. In both cases, the production polarisation is determined from the observed angular distribution of Λ b → J/ψ Λ decays. Whilst the measured transverse production polarisation is small, polarisations of O(10%) cannot be excluded. Polarised Λ b baryons can also be obtained from decays of heavier b-baryons, for example in decays of the Σ ( * ) b [22]. The only existing measurements of the angular distribution of the Λ b → Λ + − decay come from the LHCb experiment [23]. Due to the limited size of their dataset, LHCb only studied a subset of the angular distribution that could be accessed from single angle projections on the lepton-and hadron-side. With the much larger data sets that will be available at the LHC experiments after run 2 of the LHC, the experiments will be able to probe the full angular distribution. However, the sheer number of observables involved will most likely require an analysis of the moments of the angular distribution (see for example ref. [24]) rather than the conventional approach of fitting for the angular observables. This approach is discussed in section 6, where we provide the weighting functions needed to extract the observables.

Angular distribution
The angular distribution of the Λ b → Λ + − decay has been previously studied in refs. [1,2]. In this paper we extend those studies to include the case where the Λ b baryon is produced with a transverse polarisation. We start by expanding the differential decay rate for the Λ b → Λ + − decay in terms of generalised helicity amplitudes

JHEP11(2017)138
which depends on five angles, Ω = (θ l , φ l , θ b , φ b , θ), and the dilepton invariant mass squared, q 2 . The angular basis is illustrated in figure 1. The helicity basis is defined starting from the normal vector between the direction of the Λ b baryon in the lab-frame and the beamaxis of the experiment (n =p Λ b ×p beam ). This is an appropriate choice when considering transverse production polarisation of the Λ b baryon. Equation 2.1 involves three sets of helicity amplitudes: H m,J λ Λ ,λ (q 2 ) describing the decay of the Λ b baryon into a Λ baryon with helicity λ Λ and a dilepton pair with helicity λ ; h m,J λ 1 ,λ 2 describing the decay of the dilepton system to leptons with helicities λ 1 and λ 2 ; and h Λ λp,0 describing the decay Λ → pπ to a proton with helicity λ p . The index J refers to the spin of the dilepton system, which can either be zero or one. When J = 0, λ = 0, and when J = 1, λ can take the values −1, 0, +1. The helicity labels λ p , λ Λ , λ 1 and λ 2 can take the values ±1/2. Angular momentum conservation in the Λ b decay requires |λ Λ −λ | = 1/2. The factor (−1) J+J originates from the structure of the Minkowski metric tensor, see ref. [25] for details. The remaining index, m = V, A, denotes the decay of the dilepton system by either a vector or an axial-vector current. The term ρ λ Λ −λ ,λ Λ −λ is the polarisation density matrix for the transverse polarisation of the Λ b . The matrix is a two-by-two matrix (with Tr(ρ) = 1) given by Finally, the D j m,m (φ, θ, −φ) are Wigner-D functions. An explicit form of the Wigner-D functions is given in appendix A.

Lepton system amplitudes
There are two sets of amplitudes for the dilepon system, with either a vector or an axialvector current, where γ µ is a Dirac γ-matrix and ε µ is a polarisation vector. These amplitudes evaluate to [1] Figure 1. The Λ b → Λ + − decay is described by five angles: the angle, θ, between the direction of the Λ baryon and the normal vectorn in the Λ b rest-frame; and two sets of helicity angles, describing the decays of the Λ baryon (θ b , φ b ) and the dilepton system (θ l , φ l ). For transverse production polarisationn is chosen to bep Λ b ×p beam . The helicity angles are then defined with respect to this normal vector through the coordinate systems (x Λ ,ŷ Λ ,ẑ Λ ) and (x ¯ ,ŷ ¯ ,ẑ ¯ ). Thê z axis points in the direction of the Λ/dilepton system in the Λ b rest-frame. The angle between the two decay planes in the Λ b rest frame is χ = φ l + φ b . The angles θ l , θ b and χ are sufficient to parameterise the angular distribution of the decay in the case of zero production polarisation where m l is the lepton mass and β l is the lepton velocity in the dilepton rest frame (| p l |/E l ), i.e.
The amplitudes with J = 0 vanish in the case that the lepton mass is zero (when β l = 1). Under the Parity transformation (2.6)

Hadron system amplitudes
On the hadron side, the Λ decay amplitudes can be expressed in terms of the well known Λ asymmetry parameter [26]

Helicity and transversity amplitudes
After replacing the lepton and hadron-side amplitudes with the expressions given in sections 2.2 and 2.1, the angular distribution can be expanded in terms of 10 helicity amplitudes, , α Λ and a set of kinematic factors that come from the lepton-side amplitudes. For the remainder of this paper it is convenient to absorb a common factor of q 2 from the lepton-side amplitudes into these helicity amplitudes, i.e.
By absorbing this factor, the only kinematic dependence outside of H m,J λ Λ ,λ (q 2 ) comes from factors of β l .
The helicity amplitudes can be replaced by a corresponding set of transversity amplitudes for the decay that separate the vector and axial-vector contributions on the hadronside: the amplitudes A R, L 1 and A R, L 0 depend only on the vector contribution to H λ Λ ,λ (i.e. on Λ|sγ µ b|Λ b ); and the amplitudes A R, L ⊥1 and A R, L ⊥0 depend only on the axial-vector contribution to H λ Λ ,λ (i.e. on Λ|sγ µ γ 5 b|Λ b ). To do this, we start by re-writing the original helicity amplitudes as where the indices L and R refer to left-and right-handed chiralities of the dilepton system, respectively. This is followed by the replacements Here, the subscript t refers to the time-like polarisation vector of the dilepton system.

JHEP11(2017)138
3 Observables Expanding out the sum in eq. (2.1), gives 34 different angular terms Integrating this expression over Ω yields the differential decay rate as a function of q 2 , This can be used to define a set of normalised angular observables

Angular terms
The first ten angular terms are These terms are accessible even if the Λ b baryon is unpolarised and have been previously studied in refs. [2,27]. There is a straightforward relationship between our observables and those of ref. [2], with The remaining 24 terms are only non-vanishing if P Λ b is non-zero. Terms K 11 through K 16 have a similar dependence to K 1 through K 6 . These are The observables K 13 and K 16 are trivially related to K 6 and K 3 through K 13 = −P Λ b K 6 and K 16 = −P Λ b K 3 and can therefore be used as an experimental consistency check or to determine P Λ b . The observables K 11 , K 12 , K 14 and K 15 have a similar structure to K 1 , K 2 , K 4 and K 5 but, unlike in those observables, the amplitudes with λ = 0 enter with a different relative sign to those with λ = ±1. The observables K 17 through K 34 also involve new combinations of amplitudes that are not accessible if the Λ b baryon is unpolarised. They are

JHEP11(2017)138
3) The angular terms K 29 and K 31 are zero in the massless lepton limit.

Angular distribution of Λ b → J/ψ Λ
The angular distribution of the Λ b → J/ψ Λ decay is a limiting case of eq. (2.1), with a pure vector current in the dilepton system. In this limit, the expression collapses to the one JHEP11(2017)138 given in refs. [28,29] with β l ∼ 1. The amplitudes a ± and b ± in refs. [28,29] are related to the ones in this paper by (5.1)

(6.5)
More compact expressions can also be found for many of the other observables. Note, the different sets of weighting functions can lead to different experimental precision on the normalised moments. In general, the longer form of the weighting functions provides the best precision.

Standard model predictions
In order to describe the SM contribution to the decay amplitudes, an effective field theory approach is used. The Hamiltonian for the decay is factorised into local four-fermion operators and Wilson coefficients (see for example ref. [31]). The Wilson coefficients describe the short-distance contributions from the heavy SM particles. Numerical values for the SM predictions, in the case that P Λ b = 1, are provided in appendix B in two q 2 ranges: at large hadronic recoil, in the range 1 < q 2 < 6 GeV 2 /c 4 , and at low hadronic recoil, in the range 15 < q 2 < 20 GeV 2 /c 4 . To evaluate SM predictions for the different angular observables we use the EOS flavour tool [32]. At low hadronic recoil, the SM calculations employ an operator product expansion of the four-quark contributions to the matrix element in powers of Λ QCD / q 2 [33]. At large recoil, EOS uses some of the known α s corrections to charm loop processes. However, potentially large contributions from hard spectator scattering [34] and soft gluon emission [35] are neglected. The formfactors for the Λ b → Λ transition are taken from a recent Lattice QCD calculation in ref. [27]. These form-factors enable the observables to be computed with high-precision. The form-factors at large hadronic recoil have also been calculated in the framework of light-cone-sum-rules, see for example Refs. [37] and [38]. The SM Wilson coefficients are JHEP11(2017)138 computed in EOS to NNLO in QCD. The Λ b lifetime and CKM matrix elements are taken from the latest experimental values [26]. The quark masses are taken in the MS scheme. Tables 2 and 3 in appendix B also provide 68% confidence level intervals for the SM predictions. To evaluate these intervals: the form-factors from ref. [27] have been varied within their full covariance matrix; the Λ b lifetime, the Λ asymmetry parameter and CKM matrix elements are varied within their experimental precision [26,36]; the scale dependence of Wilson coefficients C i (µ) is explored by varying the scale, µ, in the range m b /2 < µ < 2m b ; and in keeping with ref. [39] a 3% correction to the amplitudes from hadronic matrix elements is considered (see also ref. [40]).  [41][42][43]. In the low-recoil range the sensitivity to C NP 9 is reduced. In ref. [2], the authors point out that the observables at low hadronic recoil place constraints on six combinations of Wilson coefficients

Low-hadronic recoil
where C V contains contributions from C 7 and C 9 . The primed coefficients correspond to right-handed currents whose contribution is vanishingly small in the SM. The shortdistance dependence of K 1 -K 34 on ρ ± 1 , ρ ± 3 , ρ 2 and ρ 4 is provided for completeness in appendix D.
If the Λ b is unpolarised, the decay rate is insensitive to the short-distance contribution Im(ρ 2 ) but provides sensitivity to ρ ± 1 , Re(ρ 2 ), ρ ± 3 , Re(ρ 4 ) and Im(ρ 4 ). The polarised observables also depend on these short-distance contributions but have different formfactor dependencies. This permits a new set of checks of the OPE and the form-factors. The short-distance combination Im(ρ 2 ) can also be determined from M 19 , M 25 , M 30 and M 34 . Furthermore, in K 1 -K 10 the short-distance contributions ρ + 1 and ρ − 1 always appear together as a sum. Using the polarised observables , ρ + 1 and ρ − 1 can be separated, e.g. by using

JHEP11(2017)138
where f V ⊥ and f A ⊥ are helicity form-factors (see for example ref. [44]). A similar trick can be used to separate ρ + 3 and ρ − 3 using K 24 and K 8 . It is also possible to form new short-distance relationships, in which the form-factors cancel by taking ratios of the K i , 3) The short-distance combinations ρ 2 and ρ 4 can then be determined up-to their overall normalisation, independent of the hadronic form-factors, using eq. (7.3) and the relationship from ref. [2]. Similarly, one can form short-distance relationships that depend only on ρ ± 1 and ρ ± 3 Alternatively, it is possible to form ratios that depend only on the form-factors and not on the short-distance physics. For example, allow the ratios f V 0 /f V ⊥ and f A 0 /f A ⊥ to be determined independent of the ρ i .

Photon-polarisation at large hadronic-recoil
At very large hadronic recoil (q 2 1 GeV 2 /c 4 ), the angular distribution of the Λ b → Λµ + µ − decay is sensitive primarily to the Wilson coefficients C 7 and C 7 due to a pole-like enhancement of the amplitudes. The observable K 33 is proportional to Re(C 7 C 7 ) and can therefore provide a null test of the size of C 7 (in the same way as the S 3 observable in the B 0 → K * 0 µ + µ − decay). In this case, however, the observable is suppressed by the size of P Λ b . Table 1 indicates the typical precision on the angular moments that could be achieved at the LHCb experiment. The experimental precision has been estimated using pseudoexperiments corresponding approximately to the expected signal yield in the current and in a future LHCb dataset. Experimental backgrounds and non-uniform angular acceptance have been neglected in this estimate. However, these are expected to have only a small impact on the experiments sensitivity. The sensitivity that can be achieved with the large datasets that will be available at an upgraded LHCb experiment is interesting event for modest values of P Λ b .

JHEP11(2017)138
Obs. Run Table 1. Expected experimental precision on the angular moments of the Λ b → Λµ + µ − decay at the LHCb experiment. The four columns correspond to: the observed yield of 300 Λ b → Λµ + µ − candidates with 15 < q 2 < 20 GeV 2 /c 4 in the LHC run 1 dataset [23]; an expected yield of ∼1000 candidates at the end of run 2 of the LHC; an expected yield of ∼8 000 candidates in 50 fb −1 of integrated luminosity with an upgraded LHCb experiment; and an expected yield of ∼50 000 candidates in 300 fb −1 with the proposed LHCb phase II upgrade.

Conclusion
In this paper we have derived an expression for the angular distribution of the Λ b → Λµ + µ − in the case of non-zero production polarisation. This extends the number of observables in the decay from 10 to 34. These observables can be determined from moments of the Λ b → Λµ + µ − angular distribution. Explicit expressions have been provided for the observables in terms of the angular moments to enable an experiment to determine the new observables from their dataset. A phenomenological analysis has also been performed to illustrate how these observables might vary in extensions of the Standard Model. The analysis shows that there is interesting new sensitivity that can be gained if the Λ b baryon is produced polarised.

JHEP11(2017)138
A Wigner D-functions The Wigner D-functions are where the α, β and γ correspond to the Euler rotation angles needed to rotate between the reference frame of the mother particle and the helicity frame of its daughters. The relevant small d-functions are

B Numerical results
Standard Model predictions for the angular observables with P Λ b = 1 are provided in tables 2 and 3. Predictions are provided in two q 2 ranges: at large hadronic recoil, in the range 1 < q 2 < 6 GeV 2 /c 4 , and at low hadronic recoil, in the range 15 < q 2 < 20 GeV 2 /c 4 . The SM predictions are evaluated using the EOS flavour-tool. For any other choice of P Λ b , predictions for M 11 -M 34 can be achieved by multiplying the values in tables 2 and 3 by the new value of P Λ b .   Table 2. Predictions from EOS for the angular observables of the Λ b → Λµ + µ − decay with P Λ b = 1 in the range 1 < q 2 < 6 GeV 2 /c 4 . The SM calculation is described in the text. The observables M 31 and M 34 vanish due to the small size of the muon mass. Observables that depend on the imaginary part of the product of two transversity amplitudes also tend to be vanishingly small, due to the small strong phase difference between pairs of amplitudes in the SM.  Table 3. Predictions from EOS for the angular observables of the Λ b → Λµ + µ − decay with P Λ b = 1 in the range 15 < q 2 < 20 GeV 2 /c 4 . The SM calculation is described in the text. The observables M 31 and M 34 vanish due to the small size of the muon mass. Observables that depend on the imaginary part of the product of two transversity amplitudes also tend to be vanishingly small, due to the small strong phase difference between pairs of amplitudes in the SM.      . Variation of the polarisation dependent angular observables of the Λ b → Λµ + µ − decay from their SM central values in the low-recoil region (15 < q 2 < 20 GeV 2 /c 4 ) with a NP contribution to Re(C 9 ) or Re(C 9 ). The SM point is at (0, 0). To illustrate the size of the effects, P Λ b = 1 is used.