Next Article in Journal
Numerical Solutions Caused by DGJIM and ADM Methods for Multi-Term Fractional BVP Involving the Generalized ψ-RL-Operators
Next Article in Special Issue
The International Linear Collider Project—Its Physics and Status
Previous Article in Journal
Towards Fire Prediction Accuracy Enhancements by Leveraging an Improved Naïve Bayes Algorithm
Previous Article in Special Issue
Extending the Fully Bayesian Unfolding with Regularization Using a Combined Sampling Method
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Heavy Quark Expansion of Λb→Λ*(1520) Form Factors beyond Leading Order

1
Dipartimento di Fisica, Università di Torino & INFN, Sezione di Torino, I-10125 Torino, Italy
2
Theoretische Physik 1, Naturwissenschaftlich-Technische Fakultät, Universität Siegen, Walter-Flex-Straße 3, D-57068 Siegen, Germany
Symmetry 2021, 13(4), 531; https://doi.org/10.3390/sym13040531
Submission received: 5 February 2021 / Revised: 1 March 2021 / Accepted: 6 March 2021 / Published: 24 March 2021
(This article belongs to the Special Issue Particle Physics and Symmetry)

Abstract

:
I review the parametrisation of the full set of Λ b Λ * ( 1520 ) form factors in the framework of Heavy Quark Expansion, including next-to-leading-order O ( α s ) and, for the first time, next-to-leading-power O ( 1 / m b ) corrections. The unknown hadronic parameters are obtained by performing a fit to recent lattice QCD calculations. I investigate the compatibility of the Heavy Quark Expansion and the current lattice data, finding tension between these two approaches in the case of tensor and pseudo-tensor form factors, whose origin could come from an underestimation of the current lattice QCD uncertainties and higher order terms in the Heavy Quark Expansion.

1. Introduction

The flavour changing neutral current (FCNC)-mediated b s + transition plays an important role in the search for physics beyond the Standard Model (SM). Its potential has been extensively studied through the B K ( * ) + decays. Interestingly, the LHCb experiment found some discrepancies with respect to the SM predictions in a few observables: R K and R K * , which test universality between the muon and electron final states and the angular coefficient P 5 in the B K * μ + μ angular distribution [1,2,3,4,5,6]. These hints, together with all available b s + data, form a coherent pattern of discrepancies. They can be addressed by introducing New Physics (NP) effects. Low-energy fits point toward a breaking of lepton flavour universality with a combined significance for the NP hypothesis higher than 5 σ [7,8,9,10]. Whether these data show without doubt first signs of NP is not clear yet. Only further measurements with higher statistics or measurements of new processes able to corroborate these data will give a final answer.
A possibility to better understand these data is studying further b s + -mediated decays, among which baryon decays are promising candidates. The decay channel involving ground-state baryons Λ b Λ μ + μ has already received attention both from the experimental and theoretical point of view. The LHCb experiment measured the Λ b Λ μ + μ angular distribution [11,12], finding good agreement between the measured values of the angular observables and thieir SM predictions [13,14,15] for the angular observables. Even though this result might be discouraging for NP searches, Refs. [14,15] showed that this is still consistent with the NP hypothesis. In fact, angular distributions describing baryon decays are very different from those in the meson cases, and NP affects them differently.
Another possibility is studying excited Λ * states. The LHCb experiment used the decay chain Λ b Λ * ( p K ) + to measure R ( p K ) , the universality ratio between muons and electrons, finding results consistent with both the SM expectation and the measured values of R K ( * ) [16]. In this analysis, the various Λ * resonances below a certain mass threshold are not distinguished. However, in Ref. [17] it is shown that the Λ * ( 1520 ) is expected to be the most frequent among the Λ * resonances, with quite narrow mass distribution. Therefore, it is motivated to study in detail the Λ b Λ * ( 1520 ) + decay from both experimental and theoretical point of view, the latter being the main focus of this paper.
The determination of the form factors of the Λ b Λ * ( 1520 ) decays had already been object of study in the literature [18,19,20,21,22]. The Heavy Quark Expansion (HQE) of the form factors up to next-to-leading (NLO) order in α s and leading power in 1 / m b has been employed [20,21], using Quark Models to constrain the unknown hadronic parameters [19]. Recently, a lattice QCD determination of the full base of form factors has become available [22], but constrained to the low-recoil region.
In this work, I investigate the compatibility between the HQE form factors and the recent lattice QCD determination. At this scope, I perform a HQE of form factors including NLO α s and next-to-leading power (NLP) 1 / m b corrections, the latter being not known so far in the literature. The results are then matched onto the lattice QCD calculation. In the HQE the number of independent, hadronic parameters is reduced compared to the lattice QCD case, introducing strict correlations among the form factors. I perform a fit to the lattice QCD results to obtain the central values, uncertainties, and correlations among the HQE parameters. The comparison between lattice QCD results and HQE predictions shows a tension between the two methods for tensor and pseudo-tensor form factors, whose origin is not yet completely determined.
This paper is organised as follows: in Section 2 I present the HQE of the form factors; in Section 3 I discuss the fit to lattice data; in Section 4 I conclude. Appendix A and Appendix B report details on the calculation of the form factors and Appendix C contains the covariance matrix for the fitted values of the HQE parameters.

2. Setup

In the following, I investigate the form factors mediating the transition Λ b ( p , s b ) Λ * ( 1520 ) ( k , η ( λ Λ ) , s Λ ) , where p and k are the momenta of the initial and final states, respectively, s b and s Λ are the rest-frame helicities of the two baryons and η ( λ Λ ) is the polarisation vector of the Λ * for each polarisation state λ Λ . Since in this work I refer to the Λ * ( 1520 ) , only, I denote this state as Λ * in the following. It is worth noticing that the Λ * is considered stable in this discussion. The subsequent Λ * decay has to be taken into account when comparing to experimental data. In Refs. [20,21] the Λ * N K ¯ ( The Λ * N K ¯ is the Λ * decay mode with larger branching fraction and the one that will be employed by the LHCb collaboration to reconstruct the Λ * in future experimental analysis.)decay is discussed, and the four-dimensional differential decay width of the decay chain Λ b Λ * ( N K ¯ ) + is presented.
I define the helicity form factors for Λ b ( p , s b ) Λ * ( k , η ( λ Λ ) , s Λ ) as
Λ * ( k , η ( λ Λ ) , s Λ ) | s ¯ γ μ b | Λ b ( p , s b ) = + u ¯ α ( k , η ( λ Λ ) , s Λ ) i F i ( q 2 ) Γ V , i α μ u ( p , s b ) , Λ * ( k , η ( λ Λ ) , s Λ ) | s ¯ γ μ γ 5 b | Λ b ( p , s b ) = u ¯ α ( k , η ( λ Λ ) , s Λ ) i G i ( q 2 ) γ 5 Γ A , i α μ u ( p , s b ) , Λ * ( k , η ( λ Λ ) , s Λ ) | s ¯ i σ μ ν q ν b | Λ b ( p , s b ) = u ¯ α ( k , η ( λ Λ ) , s Λ ) i T i ( q 2 ) Γ T , i α μ u ( p , s b ) , Λ * ( k , η ( λ Λ ) , s Λ ) | s ¯ i σ μ ν q ν γ 5 b | Λ b ( p , s b ) = u ¯ α ( k , η ( λ Λ ) , s Λ ) i T i 5 ( q 2 ) γ 5 Γ T 5 , i α μ u ( p , s b ) ,
where u ¯ α is the spin 3 / 2 projector of a Rarita-Schwinger object [23]. The Dirac structures Γ L , i α μ , with L = V , A , T , T 5 are given in Appendix A. In the cases L = V , A I adopt the parametrisation in Ref. [24], while for L = T , T 5 I follow and adapt the parametrisation in Refs. [20,22]. In the following I use the convention σ μ ν = i 2 ( γ μ γ ν γ ν γ μ ) .

The Heavy Quark Expansion

In the low-recoil limit, a HQE of the Λ b Λ * form factors can be performed. At leading power in 1 / m b and leading order in α s , the hadronic matrix element for Λ b Λ * transitions reads:
Λ * ( k , η , s Λ ) | s ¯ Γ μ b | Λ b ( p , s b ) = 4 u ¯ α ( k , η , s Λ ) ζ α ( k , η , s Λ ) Γ μ u ( m Λ b v , s b ) ,
where v is the velocity of the initial state and Γ μ denotes a Dirac structure. In the following, I focus on the cases Γ μ = γ μ , γ μ γ 5 , i σ μ ν q ν , i σ μ ν q ν γ 5 . The most general decomposition for the leading-order and leading-power contribution ζ α reads
ζ α = v α [ ζ 1 + ζ 2 v ] ,
where ζ 1 and ζ 2 are the leading Isgur-Wise (IW) functions.
The discussion of 1 / m b and α s corrections closely follows Refs. [25,26,27]. In this spirit, I replace the (axial-)vector and (pseudo-)tensor currents with:
s ¯ γ μ b s ¯ J V ( A ) μ h v = ( 1 + C 0 ( v ) ) s ¯ γ μ ( γ 5 ) h v ± C 1 ( v ) v μ s ¯ ( γ 5 ) h v + 1 2 m b s ¯ Δ J V ( A ) μ h v ,
s ¯ ( γ 5 ) i σ μ ν q ν b s ¯ J T ( 5 ) μ h v = ( 1 + C 0 ( t ) ) s ¯ ( γ 5 ) i σ μ ν q ν h v + 1 2 m b s ¯ Δ J T ( 5 ) μ h v .
The matching coefficients at NLO read [25,26]
C 0 ( v ) ( μ ) = α s C F 4 π 3 log μ m b + 4 + O ( α s 2 ) , C 1 ( v ) ( μ ) = + α s C F 2 π + O ( α s 2 ) , C 0 ( t ) ( μ ) = α s C F 4 π 5 log μ m b + 4 + O ( α s 2 ) .
For numerical purposes, the scale of the Wilson coefficients is set to μ 2 GeV . The NLP 1 / m b corrections due to the expansion of the current are parametrised as
Λ * ( k , η , s Λ ) | s ¯ Δ J V ( A , T , T 5 ) μ b | Λ b ( p , s b ) = 4 i u ¯ α ( k , η , s Λ ) ζ α β [ O i V ( A , T , T 5 ) ] β μ u ( m Λ b v , s b ) ,
where
ζ α β = g α β ζ 1 SL + v ζ 2 SL + v α v β ζ 3 SL + v ζ 4 SL + v α γ β ζ 5 SL + v ζ 6 SL .
The functions ζ 1 6 SL are the subleading Isgur-Wise functions, and they correspond to all the independent Dirac structures that can appear in ζ α β . The possible operators O i μ β are listed in Ref. [25]. Out of the possible six of them, only the operator [ O 1 Γ ] β μ arises at order 1 / m b , while the others are suppressed by O ( α s / m b ) and therefore are beyond the precision here required. Therefore, the only contributions that I consider for this analysis come from:
O 1 V β μ = + γ μ γ β , O 1 A β μ = γ 5 γ μ γ β , O 1 T β μ = + i σ μ ν q ν γ β , O 1 T 5 β μ = + i γ 5 σ μ ν q ν γ β ,
By means of Dirac algebra, and using the properties of Rarita-Schwinger objects in Ref. [23], inserting Equation (9) in Equation (8) yields
Λ * ( k , η , s Λ ) | s ¯ Δ J V μ b | Λ b ( p , s b ) = 2 { 2 u ¯ α ( k , η , s Λ ) u ( m Λ b v , s b ) g α μ ( ζ 1 SL + ζ 2 SL ) + u ¯ α ( k , η , s Λ ) γ μ ( m Λ b v , s b ) v α ( ζ 3 SL ζ 4 SL 2 ζ 2 SL 2 ζ 5 SL ) + 2 u ¯ α ( k , η , s Λ ) u ( m Λ b v , s b ) v α v μ ( ζ 4 SL + 2 ζ 6 SL ) } ,
Λ * ( k , η , s Λ ) | s ¯ Δ J A μ b | Λ b ( p , s b ) = 2 { 2 u ¯ α ( k , η , s Λ ) γ 5 u ( m Λ b v , s b ) g α μ ( ζ 1 SL + ζ 2 SL ) u ¯ α ( k , η , s Λ ) γ 5 γ μ ( m Λ b v , s b ) v α ( ζ 3 SL + ζ 4 SL + 2 ζ 2 SL + 2 ζ 5 SL ) + 2 u ¯ α ( k , η , s Λ ) γ 5 u ( m Λ b v , s b ) v α v μ ( ζ 4 SL 2 ζ 6 SL ) } ,
Λ * ( k , η , s Λ ) | s ¯ Δ J T μ b | Λ b ( p , s b ) = 2 { u ¯ α ( k , η , s Λ ) γ μ u ( m Λ b v , s b ) v α × [ 2 m Λ b ζ 1 SL + 2 m Λ * ζ 2 SL + ( m Λ b + m Λ * ) ( ζ 3 SL + 2 ζ 6 SL ) + s m Λ * 2 + m Λ b m Λ * q 2 m Λ b ζ 4 SL ] + 2 u ¯ μ ( k , η , s Λ ) u ( m Λ b v , s b ) × [ ( m Λ b m Λ * ) ζ 1 SL m Λ * 2 m Λ * m Λ b q 2 m Λ b ζ 2 SL ] + 2 u ¯ α ( k , η , s Λ ) γ μ u ( m Λ b v , s b ) k α ( ζ 1 SL ζ 2 SL ) + u ¯ α ( k , η , s Λ ) u ( m Λ b v , s b ) v α k μ ( ζ 3 SL + 2 ζ 2 SL + ζ 4 SL + 2 ζ 6 SL ) + 4 u ¯ α ( k , η , s Λ ) u ( m Λ b v , s b ) v μ k α ζ 2 SL + u ¯ α ( k , η , s Λ ) u ( m Λ b v , s b ) v α v μ × [ m Λ b ζ 3 SL 2 m Λ b ζ 2 SL + ( 2 m Λ * m Λ b ) ζ 4 SL + 2 m Λ b ζ 6 SL ] } ,
Λ * ( k , η , s Λ ) | s ¯ Δ J T 5 μ b | Λ b ( p , s b ) = 2 { u ¯ α ( k , η , s Λ ) γ 5 γ μ u ( m Λ b v , s b ) v α [ 2 m Λ b ζ 1 SL + 2 m Λ * ζ 2 SL + ( m Λ b m Λ * ) ( ζ 3 SL + 2 ζ 6 SL ) m Λ * 2 m Λ b m Λ * q 2 m Λ b ζ 4 SL ] + 2 u ¯ μ ( k , η , s Λ ) γ 5 u ( m Λ b v , s b ) × ( m Λ b + m Λ * ) ζ 1 SL + m Λ * 2 + m Λ * m Λ b q 2 m Λ b ζ 2 SL + 2 u ¯ α ( k , η , s Λ ) γ 5 γ μ u ( m Λ b v , s b ) k α ( ζ 1 SL + ζ 2 SL ) + u ¯ α ( k , η , s Λ ) γ 5 u ( m Λ b v , s b ) v α k μ ( ζ 3 SL 2 ζ 2 SL ζ 4 SL + 2 ζ 6 SL ) 4 u ¯ α ( k , η , s Λ ) γ 5 u ( m Λ b v , s b ) v μ k α ζ 2 SL + u ¯ α ( k , η , s Λ ) γ 5 u ( m Λ b v , s b ) v α v μ × [ m Λ b ζ 3 SL + 2 m Λ b ζ 2 SL + ( 2 m Λ * + m Λ b ) ζ 4 SL + 2 m Λ b ζ 6 SL ] } .
The number of independent subleading IW functions can be reduced by using equations of motions. In particular, for this decay, the relation v β ζ α β = 0 gives the following conditions:
ζ 1 SL + ζ 3 SL + ζ 6 SL = 0 ,
ζ 2 SL + ζ 4 SL + ζ 5 SL = 0 .
I choose to retain as independent quantities the subleading IW functions ζ 1 SL , ζ 2 SL , ζ 3 SL and ζ 4 SL .
Corrections to the form factors arise also by non-local insertions of the heavy quark Lagrangian at order 1 / m b . Following the discussions in Refs. [25,28,29], non-local insertion of the kinetic operator yields a universal shift proportional to the tree-level matrix elements. Hence I reabsorb such shift in a redefinition of the leading order IW function ζ 1 . Non-local insertion of the chromomagnetic operator can be parametrised as [30]
R μ ν α u ¯ α ( k , η , s Λ ) Γ 1 + v 2 i σ μ ν u ( m Λ b v , s b ) ,
where the object R μ ν α is an antisymmetric tensor in the indices μ and ν and contains the velocity v. Γ symbolises all the possible Dirac structures which mediate Λ b Λ * transitions. Using equations of motion, it can be shown that no possible form of R μ ν α gives a non-zero contribution for the chromomagnetic operator.
The expressions of the form factors in terms of the leading and subleading IW functions are obtained by matching the helicity amplitudes with their HQE expansion. For the vector current, by comparing Equations (A17)–(A20) to Equations (A33)–(A36), I get
F 1 / 2 , 0 = s + 2 m Λ b 3 / 2 m Λ * 3 / 2 { ζ 1 s ( 1 + C 0 ( v ) ) + C 1 ( v ) s + 2 m Λ b ( m Λ b + m Λ * ) + ζ 2 s ( 1 + C 0 ( v ) ) m Λ * 2 + m Λ * m Λ b q 2 m Λ b ( m Λ b + m Λ * ) + C 1 ( v ) s + 2 m Λ b ( m Λ b + m Λ * ) ζ 1 SL λ + m Λ b 2 ( m Λ * 2 m Λ b 2 + q 2 ) m Λ b 2 ( m Λ * + m Λ b ) ζ 2 SL ( m Λ * 2 m Λ b 2 + q 2 ) ( m Λ * + m Λ b ) s ζ 3 SL ( 2 m Λ * 2 + 3 m Λ * m Λ b + m Λ b 2 2 q 2 ) 2 m Λ b 2 ( m Λ * + m Λ b ) ζ 4 SL ( m Λ * 2 + 3 m Λ * m Λ b + 2 m Λ b 2 q 2 ) 2 m Λ b 2 ( m Λ * + m Λ b ) } ,
F 1 / 2 , t = s s + 2 m Λ b 3 / 2 m Λ * 3 / 2 { ζ 1 ( 1 + C 0 ( v ) ) C 1 ( v ) m Λ * 2 m Λ b 2 q 2 2 m Λ b ( m Λ b m Λ * ) ζ 2 1 2 m Λ b ( m Λ b m Λ * ) 2 ( 1 + C 0 ( v ) ) ( m Λ * 2 m Λ * m Λ b q 2 ) + C 1 ( v ) ( m Λ * 2 m Λ b 2 q 2 ) 1 m Λ b m Λ * s m Λ * 2 m Λ b 2 ζ 1 SL ζ 2 SL + 2 m Λ * 2 m Λ * m Λ b m Λ b 2 2 q 2 2 m Λ b 2 ( m Λ b m Λ * ) ζ 3 SL m Λ * 2 + m Λ * m Λ b 2 m Λ b 2 q 2 2 m Λ b 2 ( m Λ b m Λ * ) ζ 4 SL ,
F 1 / 2 , = s + 2 m Λ b 3 / 2 m Λ * 3 / 2 s ( 1 + C 0 ( v ) ) ( ζ 1 ζ 2 ) m Λ * ( ζ 1 SL + ζ 2 SL ) + s 2 m Λ b ( ζ 3 SL + ζ 4 SL ) ,
F 3 / 2 , = s + 2 m Λ b 3 / 2 m Λ * 1 / 2 ζ 1 SL + ζ 2 SL ,
and for the axial vector current, by matching Equations (A21)–(A24) to Euqationa (A37)–(A40) I obtain
G 1 / 2 , 0 = s 2 m Λ b 3 / 2 m Λ * 3 / 2 { ζ 1 s + ( 1 + C 0 ( v ) ) + C 1 ( v ) s 2 m Λ b ( m Λ b m Λ * ) ζ 2 s + ( 1 + C 0 ( v ) ) m Λ * 2 m Λ * m Λ b q 2 m Λ b ( m Λ b m Λ * ) C 1 ( v ) s 2 m Λ b ( m Λ b m Λ * ) ζ 1 SL λ + m Λ b 2 ( m Λ * 2 m Λ b 2 + q 2 ) m Λ b 2 ( m Λ b m Λ * ) + ζ 2 SL ( m Λ * 2 m Λ b 2 + q 2 ) ( m Λ b m Λ * ) s + ζ 3 SL ( 2 m Λ * 2 3 m Λ * m Λ b + m Λ b 2 2 q 2 ) 2 m Λ b 2 ( m Λ b m Λ * ) + ζ 4 SL ( m Λ * 2 3 m Λ * m Λ b + 2 m Λ b 2 q 2 ) 2 m Λ b 2 ( m Λ b m Λ * ) } ,
G 1 / 2 , t = s + s 2 m Λ b 3 / 2 m Λ * 3 / 2 { ζ 1 ( 1 + C 0 ( v ) ) C 1 ( v ) m Λ * 2 m Λ b 2 q 2 2 m Λ b ( m Λ b + m Λ * ) + ζ 2 1 2 m Λ b ( m Λ b + m Λ * ) 2 ( 1 + C 0 ( v ) ) ( m Λ * 2 + m Λ * m Λ b q 2 ) C 1 ( v ) ( m Λ * 2 m Λ b 2 q 2 ) 1 m Λ b + m Λ * s m Λ * 2 m Λ b 2 ζ 1 SL + ζ 2 SL + 2 m Λ * 2 + m Λ * m Λ b m Λ b 2 2 q 2 2 m Λ b 2 ( m Λ b + m Λ * ) ζ 3 SL + m Λ * 2 m Λ * m Λ b 2 m Λ b 2 q 2 2 m Λ b 2 ( m Λ b + m Λ * ) ζ 4 SL ,
G 1 / 2 , = s 2 m Λ b 3 / 2 m Λ * 3 / 2 s + ( 1 + C 0 ( v ) ) ( ζ 1 + ζ 2 ) + m Λ * ( ζ 1 SL ζ 2 SL ) + s + 2 m Λ b ( ζ 3 SL ζ 4 SL ) ,
G 3 / 2 , = s 2 m Λ b 3 / 2 m Λ * 1 / 2 ζ 1 SL ζ 2 SL .
For the tensor current, the comparison between Equations (A25)–(A27) and Equations (A41)–(A43) yields
T 1 / 2 , 0 = s + m Λ b 3 / 2 m Λ * 1 / 2 { ( 1 + C 0 ( t ) ) ( ζ 1 ζ 2 ) s m Λ * 2 + m Λ b 2 q 2 m Λ b ( ζ 1 SL ζ 2 SL ) s 2 m Λ b ( ζ 4 SL + ζ 3 SL ) } , T 1 / 2 , = s + m Λ b 3 / 2 m Λ * 1 / 2 { ( 1 + C 0 ( t ) ) ζ 1 + m Λ * ( m Λ b + m Λ * ) q 2 m Λ b ( m Λ * + m Λ b ) ζ 2 s + m Λ * m Λ * ( m Λ b m Λ * ) + q 2 m Λ b ( m Λ * + m Λ b ) ζ 1 SL + m Λ * ( m Λ b m Λ * ) m Λ * + m Λ b ζ 2 SL
+ s 2 m Λ b ζ 3 SL + m Λ * ( m Λ b + m Λ * ) q 2 m Λ b ( m Λ * + m Λ b ) ζ 4 SL } ,
T 3 / 2 , = + s + m Λ b 3 / 2 m Λ * 1 / 2 ζ 1 SL ( m Λ b m Λ * ) + m Λ * 2 m Λ b m Λ * q 2 m Λ b ζ 2 SL ,
while for the axial-tensor form factors the comparison between Equations (A29)–(A31) and Equations (A44)–(A46) gives
T 1 / 2 , 0 5 = s m Λ b 3 / 2 m Λ * 1 / 2 { ( 1 + C 0 ( t ) ) ( ζ 2 + ζ 1 ) s + m Λ * 2 + m Λ b 2 q 2 m Λ b ( ζ 1 SL ζ 2 SL ) s + 2 m Λ b ( ζ 3 SL ζ 4 SL ) } , T 1 / 2 , 5 = s m Λ b 3 / 2 m Λ * 1 / 2 { ( 1 + C 0 ( t ) ) ζ 1 m Λ * ( m Λ b + m Λ * ) q 2 m Λ b ( m Λ b m Λ * ) ζ 2 s + + m Λ * m Λ * ( + m Λ b + m Λ * ) q 2 m Λ b ( m Λ b m Λ * ) ζ 1 SL + m Λ * ( m Λ b + m Λ * ) m Λ b m Λ * ζ 2 SL
s + 2 m Λ b ζ 3 SL + m Λ * ( m Λ b + m Λ * ) q 2 m Λ b ( m Λ b m Λ * ) ζ 4 SL } ,
T 3 / 2 , 5 = s m Λ b 3 / 2 m Λ * 1 / 2 ζ 1 SL ( m Λ b + m Λ * ) + m Λ * 2 + m Λ b m Λ * q 2 m Λ b ζ 2 SL .
The expressions in Equations (17)–(30) have been checked against the results in Ref. [21], where the HQE for the form factors including NLO α s corrections is presented. I find agreement with the expressions reported in Ref. [21], apart from a sign flip in the term proportional to C 1 ( v ) ζ 2 in G 1 / 2 , t . This misalignment is not expected to invalidate the analysis in Ref. [21] since it appears in a term that is α s suppressed. In the following analysis, I stick to my findings and adopt the signs in Equation (22).

3. Form Factors Relations and Comparison with Lattice QCD Results

The first lattice QCD calculation for the full basis of Λ b Λ * form factors is presented in Ref. [22]. The calculation is performed in the low-recoil region, very close to the zero-recoil point q max 2 = ( m Λ b m Λ * ) 2 . Two lattice points per form factor are obtained, allowing to determine the normalisation and the slope of each form factor. In the kinematical limit where the lattice QCD computation is valid, it is more convenient to substitute the variable q 2 with the adimensional variable w = p · k / ( 2 m Λ b m Λ * ) = ( m Λ b 2 + m Λ * 2 q 2 ) / ( 2 m Λ b m Λ * ) , where the zero-recoil point corresponds to w = 1 . The continuum extrapolation of the results in Ref. [22] is performed using the following functional form for each of the form factor f i :
f i = F i + A i ( w 1 ) .
Values for the coefficients F i and A i and their covariance matrix are given in ancillary files of Ref. [22].
Given that both the results in Ref. [22] and the parametrisation based on HQE in Section 2 are valid in the low-recoil region, the former is suitable to extract the unknown, hadronic parameters for the leading and sub-leading IW functions introduced in Section 2. I want to stress that it is not possible to extrapolate these results to the large-recoil region without any further information on the form factors valid at low q 2 .
The form factor base employed in Ref. [22] differs from the one presented in Section 2 and the matching between the two is given in Appendix B. In the following, I denote with capital letters the HQE base and with lower cases the base of Ref. [22].

3.1. Relations in the Zero-Recoil Point

I study the form factors first at the zero-recoil point. At this particular kinematical configuration, all the axial-vector and pseudo-tensor HQE form factors become zero because they are weighted by the factor s . The remainder is further simplified since the terms associated with ζ 1 , ζ 2 , ζ 3 SL , ζ 4 SL , are always proportional to s , hence vanish, leaving only ζ 1 SL and ζ 2 SL to determine the form factors in the zero-recoil point. Even more interestingly, from Equations (17)–(30) it can be seen that only the combination ζ 1 SL + ζ 2 SL appears, with different weights for each of the form factor. Therefore, HQE provides predictions for ratios of form factors at zero-recoil independent of the IW functions. They read:
F 1 / 2 , 0 F 1 / 2 , = F 1 / 2 , 0 F 3 / 2 , = 2 m Λ b m Λ * m Λ b + m Λ * = 1.15 , F 1 / 2 , F 3 / 2 , = 1 , T 1 / 2 , 0 T 1 / 2 , = 2 m Λ b + m Λ * m Λ b m Λ * = 3.48 , T 1 / 2 , 0 T 3 / 2 , = 2 m Λ * m Λ b m Λ * = 0.74 , T 1 / 2 , T 3 / 2 , = m Λ * m Λ b + m Λ * = 0.21 ,
where errors on the baryon masses have been neglected. These predictions have to be compared with the results in Ref. [22]. Extracting pseudodata for the parameters entering the lattice QCD form factors, according to their central values and covariance matrix, I find (in the HQE basis)
F 1 / 2 , 0 latt F 1 / 2 , latt = 0.63 ± 1.82 , F 1 / 2 , 0 latt F 3 / 2 , latt = 0.94 ± 0.14 , F 1 / 2 , latt F 3 / 2 , latt = 1.48 ± 0.39 , T 1 / 2 , 0 latt T 1 / 2 , latt = 5.62 ± 7.86 , T 1 / 2 , 0 latt T 3 / 2 , latt = + 1.03 ± 0.20 , T 1 / 2 , latt T 3 / 2 , latt = 0.18 ± 0.05 .
The central values showed above are the medians of the distributions and the errors correspond to the 68 % intervals. Concerning the ratios F 1 / 2 , 0 latt / F 1 / 2 , latt and T 1 / 2 , 0 latt / T 1 / 2 , latt , I stress that their uncertainties are very large to reflect the fact that their distributions are highly non gaussian. For the other ratios, the gaussian approximation well describes lattice data. The predictions for the ratios obtained in the HQE framework and using the lattice QCD results agree within 1.1 σ .

3.2. Form Factors Parametrisation and Fit

The HQE of the Λ b Λ * form factors presented in Section 2 depends on 6 unknown IW functions. The HQE does not predict the w dependence of the IW functions, and it has to be inferred from other principles. Since the formalism of HQE is valid mainly in the low-recoil region for b s transitions, the form factors can be expanded around the zero-recoil point. Substituting q 2 with w, the IW functions ζ i can be expanded as
ζ i = n = 0 N ζ i ( n ) n ! ( w 1 ) n .
The truncation order N of the expansion depends on the precision required and how far from w = 1 the form factors are evaluated. The parameters ζ i ( n ) are unknown and have to be fixed using some dynamical information. Notice that since the Λ * is not a ground state baryon, the HQE does not predict any normalisation for the leading IW functions.
For convenience, I perform a fit to the lattice QCD data using the base in Ref. [22], using the chiral and continuum extrapolation there provided to obtain pseudo-points for the form factors. Since the lattice QCD data do not provide information on the curvature of the form factors, it is useful to express the HQE form factors as in the form of Equation (31). At this scope, I use the parametrisation in Equation (34) up to N = 1 and then re-expand the full form factors up to the first order in w 1 . After this procedure, it can be noticed that (i) the parameters ζ 1 ( 1 ) and ζ 2 ( 1 ) appear always in the combination ζ 1 ( 1 ) + ζ 2 ( 1 ) and (ii) the parameters ζ 3 SL , ( 1 ) and ζ 4 SL , ( 1 ) appear always in the combination ζ 4 SL , ( 1 ) ζ 3 SL , ( 1 ) . Therefore these parameters cannot be determined on their own, but only the combinations ζ 1 ( 1 ) + ζ 2 ( 1 ) and ζ 4 SL , ( 1 ) ζ 3 SL , ( 1 ) are determined. This makes the number of independent, unknown HQE parameters 10.
Before discussing the fit results, a couple of technical comments are in order:
  • Given the available information in Ref. [22], it is possible to use two pseudo-points for each form factor. I choose to evaluate them at w = 1.02 and w = 1.04 . Given that the HQE parametrisation depends on fewer parameters than the lattice QCD one, it is possible to perform a fit to only a subset of the lattice QCD data. In particular, I choose to employ the data on the vector and axial-vector form factors and provide predictions for the tensor and pseudo-tensor form factors based on the fit results. I comment on the consequences of this choice in the following. I stress that this is a common procedure and has been already employed in Ref. [31].
  • The HQE based form factors are affected by uncertainties due to the unknown contributions from higher orders of expansion. By naive dimensional arguments these contributions are expected to be roughly O ( few % ) . Hence, I introduce an uncorrelated 1 % uncertainty on all HQE expressions of the form factors to take these effects into account. Comments on this choice can be found later in the text.
The fit is performed with a χ 2 minimisation, yielding at the minimum χ min 2 / d . o . f . 0.25 . This low value is a direct consequence of the poor current knowledge of the exact size of the theory uncertainties and their correlations. If the theory uncertainties were uncorrelated, the fit procedure would indicate that their natural size is smaller than the one inferred from HQE. However, such a low value for the χ min 2 / d . o . f . could also indicate that strong correlations among unknown higher-order terms in the HQE have been neglected. At the current status, it is not possible to obtain a more precise estimation of theory uncertainties and their correlations. Therefore, I retain the conservative choice of an uncorrelated 1 % uncertainty on all the form factors in the analysis.
The best-fit point for the hadronic parameters and their uncertainties is shown in Table 1, and the correlation matrix among the parameters is given in Appendix C. With this result, I compare the HQE form factors to the lattice QCD ones. The comparison is given in Figure 1, showing excellent agreement between the two parametrisations.
I use the results of the fit to obtain predictions for the tensor and pseudo-tensor form factors. The comparison between them and the lattice QCD computation is shown in Figure 2 ( The lattice QCD results employed here for tensor and pseudo-tensor form factors differ by a global sign with respect to the first version of Ref. [22]. This sign inconsistency will be fixed in a forthcoming update of Ref. [22].). The form factors h and h ˜ show a tension between lattice QCD data and HQE predictions. This tension manifests itself also when introducing the tensor and pseudo-tensor form factors in the fit. This procedure increases the χ 2 / d . o . f . , making it much higher than 1. This corroborates the choice of excluding them from the fit procedure.
Sources of these tensions can be looked for in the lattice QCD data and HQE framework. The hypotheses are mainly two: (i) the uncertainties on the lattice QCD parameters describing tensor and pseudo-tensor form factors are underestimated, and (ii) missing corrections in the HQE cause a shift in the hadronic parameters. In the case of (i), I checked explicitly the result of inflating lattice QCD uncertainties by 20 % for h and h ˜ . In Figure 3 the results of this test are shown, proving that the compatibility slightly improves, even if it is still poor in the case of h ˜ . Concerning (ii), the most important corrections in the HQE, beside the ones already discussed in this work, arise at order O ( α s / m b ) , and could produce a O ( few % ) shift in the central value of the form factors parameters. However, assessing the quantitative impact of these corrections requires understanding how they affect the form factors in a correlated way.
From these estimates, it seems that neither (i) nor (ii) can explain the tension on their own. Most likely, a combination of the two effects might help to reconcile the HQE for the tensor and pseudo-tensor form factors with the current lattice QCD data.

4. Conclusions

I revisit the Heavy Quark Expansion of the Λ b Λ * ( 1520 ) form factors including next-to-leading order α s corrections and for the first time next-to-leading power 1 / m b corrections. In this framework, form factors depend on unknown hadronic parameters which are obtained by fitting the form factor parametrisation here discussed to a recent lattice QCD computation [22]. I perform the fit using data for vector and axial-vector form factors, finding good agreement between the lattice QCD calculation and the Heavy Quark Expansion predictions. The fit results are used to predict tensor and pseudo-tensor form factors, showing tensions between the Heavy Quark Expansion-based predictions and the lattice QCD data. I discuss two possible sources of the tensions: an underestimation of the uncertainties on the lattice QCD parameters involved in these form factors and missing higher-order terms in the Heavy Quark Expansion, e.g. at order O ( α s / m b ) . Most likely, only a combination of these two effects could reconcile lattice QCD determination and Heavy Quark Expansion based parametrisation of tensor and pseudo-tensor form factors. Until then, it is not possible to claim high precision in the Heavy Quark Expansion parametrisation of Λ b Λ * ( 1520 ) form factors.
Besides this, I want to point out the need to extend the calculation of the Λ b Λ * ( 1520 ) form factors to the high-recoil region. Quark models [19] are available, although without a consistent treatment of uncertainties. It is therefore needed to perform up-to-date calculations of the Λ b Λ * ( 1520 ) form factors using e.g. sum rules at q 2 0 . Estimates of this type will allow to extrapolate the form factors to the high-recoil region and to assess the magnitude of their curvature. These studies will be crucial for future LHCb analysis of Λ b Λ * ( 1520 ) + decays.

Funding

This research received no external funding.

Data Availability Statement

The lattice QCD data employed in this analysis are available at 2009.09313.

Acknowledgments

I thank Stefan Meinel and Gumaro Rendon for insightful discussions on lattice QCD data. I am also grateful to Bernat Capdevila, Thorsten Feldmann, Paolo Gambino, Martin Jung and Danny van Dyk for useful discussions. I also acknowledge useful inputs on the text from Nico Gubernari. This work is supported by Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under grant 396021762—TRR 257 “Particle Physics Phenomenology after the Higgs Discovery” and by the Italian Ministry of Research (MIUR) under grant PRIN 20172LNEEZ.

Conflicts of Interest

The author declare no conflict of interest.

Appendix A. Details on the Form Factor Parametrisation

Concerning the vector and the axial vector currents, we follow the notation of Ref. [24]. For the vector current we have:
Γ V , ( 1 / 2 , t ) α μ = 4 m Λ b m Λ * s + 2 m Λ * s + s p α m Λ b m Λ * q 2 q μ q 2 ,
Γ V , ( 1 / 2 , 0 ) α μ = 4 m Λ b m Λ * s 2 m Λ * s + s p α m Λ b + m Λ * s + ( p + k ) μ m Λ b 2 m Λ * 2 q 2 q μ ,
Γ V , ( 1 / 2 , ) α μ = 4 m Λ b m Λ * s 2 m Λ * s + s p α γ μ 2 m Λ * s + p μ 2 m Λ b s + k μ ,
Γ V , ( 3 / 2 , ) α μ = 4 m Λ b m Λ * s 4 i ϵ α μ p k s + s γ 5 + Γ V , ( 1 / 2 , ) ,
while for the axial vector current:
Γ A , ( 1 / 2 , t ) α μ = 4 m Λ b m Λ * s 2 m Λ * s + s p α m Λ b + m Λ * q 2 q μ q 2 ,
Γ A , ( 1 / 2 , 0 ) α μ = 4 m Λ b m Λ * s + 2 m Λ * s + s p α m Λ b m Λ * s ( p + k ) μ m Λ b 2 m Λ * 2 q 2 q μ ,
Γ A , ( 1 / 2 , ) α μ = 4 m Λ b m Λ * s + 2 m Λ * s + s p α γ μ + 2 m Λ * s p μ 2 m Λ b s k μ ,
Γ A , ( 3 / 2 , ) α μ = 4 m Λ b m Λ * s + 4 i ϵ α μ p k s + s γ 5 Γ A , ( 1 / 2 , ) .
Concerning the tensor currents, we modify the parametrisation in Ref. [20] by rescaling each structure with suitable factors. We have:
Γ T , ( 1 / 2 , 0 ) α μ = 4 m Λ b m Λ * s + q 2 s + s p α ( p + k ) μ m Λ b 2 m Λ * 2 q 2 q μ ,
Γ T , ( 1 / 2 , ) α μ = 4 m Λ b m Λ * s + m Λ b + m Λ * s p α γ μ 2 m Λ * s + p μ 2 m Λ b s + k μ ,
Γ T , ( 3 / 2 , ) α μ = 4 m Λ b m Λ * s + g α μ + m Λ * s p α γ μ 2 1 m Λ * k μ + 2 m Λ * s + p μ + 2 m Λ b s + k μ ,
and
Γ T 5 , ( 1 / 2 , 0 ) α μ = 4 m Λ b m Λ * s q 2 s + s p α ( p + k ) μ m Λ b 2 m Λ * 2 q 2 q μ ,
Γ T 5 , ( 1 / 2 , ) α μ = 4 m Λ b m Λ * s m Λ b m Λ * s + p α γ μ + 2 m Λ * s p μ 2 m Λ b s k μ ,
Γ T 5 , ( 3 / 2 , ) α μ = 4 m Λ b m Λ * s g α μ m Λ * s + p α γ μ + 2 1 m Λ * k μ 2 m Λ * s p μ + 2 m Λ b s k μ .
I define the helicity amplitudes as
A Γ ( s b , s Λ , λ Λ , λ q ) = Λ * ( s Λ , η ( λ Λ ) ) | s ¯ Γ μ ϵ μ * ( λ q ) b | Λ b ( s b ) ,
where ϵ μ * ( λ q ) are a basis of polarisation vectors for the virtual W exchange with polarisation states λ q = { t , 0 , + 1 , 1 } . For details see Ref. [24]. The physical helicity amplitudes are identified by the following set:
A Γ ( + 1 / 2 , + 3 / 2 , + 1 ) A Γ ( + 1 / 2 , + 1 / 2 , + 1 , + 1 ) , A Γ ( + 1 / 2 , + 1 / 2 , 0 ) 2 3 A Γ ( + 1 / 2 , + 1 / 2 , 0 , 0 ) + 1 3 A Γ ( 3 / 2 ) ( + 1 / 2 , 1 / 2 , + 1 , 0 ) , A Γ ( + 1 / 2 , + 1 / 2 , t ) 2 3 A Γ ( + 1 / 2 , + 1 / 2 , 0 , t ) + 1 3 A Γ ( 3 / 2 ) ( + 1 / 2 , 1 / 2 , + 1 , t ) , A Γ ( + 1 / 2 , 1 / 2 , 1 ) 2 3 A Γ ( + 1 / 2 , 1 / 2 , 0 , 1 ) + 1 3 A Γ ( 3 / 2 ) ( + 1 / 2 , + 1 / 2 , 1 , 1 ) ,
where Γ represents the four possible currents. In the case of the vector and axial vector current, the helicity amplitudes are already calculated in Ref. [24]. For convenience, I adapt them to this case and report them here:
A V ( + 1 / 2 , + 3 / 2 , + 1 ) = 4 m Λ b m Λ * F 3 / 2 , ,
A V ( + 1 / 2 , + 1 / 2 , 0 ) = 2 2 m Λ b m Λ * 3 q 2 ( m Λ b + m Λ * ) F 1 / 2 , 0 ,
A V ( + 1 / 2 , + 1 / 2 , t ) = 2 2 m Λ b m Λ * 3 q 2 ( m Λ b m Λ * ) F 1 / 2 , t ,
A V ( + 1 / 2 , 1 / 2 , 1 ) = 4 3 m Λ b m Λ * F 1 / 2 , ,
and for the axial vector current:
A A ( + 1 / 2 , + 3 / 2 , + 1 ) = 4 m Λ b m Λ * G 3 / 2 , ,
A A ( + 1 / 2 , + 1 / 2 , 0 ) = 2 2 m Λ b m Λ * 3 q 2 ( m Λ b m Λ * ) G 1 / 2 , 0 ,
A A ( + 1 / 2 , + 1 / 2 , t ) = 2 2 m Λ b m Λ * 3 q 2 ( m Λ b + m Λ * ) G 1 / 2 , t ,
A A ( + 1 / 2 , 1 / 2 , 1 ) = 4 3 m Λ b m Λ * G 1 / 2 , .
In the case of tensor currents, with the definitions in Equations (A11)–(A14), I find
A T ( + 1 / 2 , + 3 / 2 , + 1 ) = 2 m Λ b m Λ * T 3 / 2 , ,
A T ( + 1 / 2 , + 1 / 2 , 0 ) = 2 3 m Λ b m m Λ * q 2 T 1 / 2 , 0 ,
A T ( + 1 / 2 , 1 / 2 , 1 ) = + 2 3 m Λ b m Λ * ( m Λ b + m Λ * ) T 1 / 2 , ,
A T ( + 1 / 2 , + 1 / 2 , t ) = 0 ,
and
A T 5 ( + 1 / 2 , + 3 / 2 , + 1 ) = + 2 m Λ b m Λ * T 3 / 2 , 5 ,
A T 5 ( + 1 / 2 , + 1 / 2 , 0 ) = + 2 3 m Λ b Λ * q 2 T 1 / 2 , 0 5 ,
A T 5 ( + 1 / 2 , 1 / 2 , 1 ) = + 2 3 m Λ b m Λ * ( m Λ b m Λ * ) T 1 / 2 , ,
A T 5 ( + 1 / 2 , + 1 / 2 , t ) = 0 .
In the heavy quark expansion, the helicity amplitudes concerning the vector current read:
A V ( + 1 / 2 , + 3 / 2 , + 1 ) = 2 s + m Λ b ( ζ 1 SL + ζ 2 SL ) ,
A V ( + 1 / 2 , + 1 / 2 , 0 ) = 2 s + m Λ b m Λ * 3 q 2 { s m Λ b ( 1 + C 0 ( v ) ) ( m Λ * 2 + m Λ * m Λ b q 2 ) + 1 2 C 1 ( v ) s + ζ 2 + s ( 1 + C 0 ( v ) ) ( m Λ * + m Λ b ) + 1 2 m Λ b C 1 ( v ) s + ζ 1 m Λ b 2 ( m Λ * 2 m Λ b 2 + q 2 ) + s s + ζ 1 SL ( m Λ * 2 m Λ b 2 + q 2 ) ζ 2 SL s 2 m Λ b 2 [ ( 2 m Λ * 2 + 3 m Λ * m Λ b + m Λ b 2 2 q 2 ) ζ 3 SL ( m Λ * 2 + 3 m Λ * m Λ b + 2 m Λ b 2 q 2 ) ζ 4 SL ] } ,
A V ( + 1 / 2 , + 1 / 2 , t ) = 2 s s + m Λ b m Λ * 3 q 2 { 1 m Λ b [ ( 1 + C 0 ( v ) ) ( m Λ * 2 + m Λ * m Λ b + q 2 ) + 1 2 C 1 ( v ) ( m Λ * 2 + m Λ b 2 + q 2 ) ] ζ 2 + [ ( 1 + C 0 ( v ) ) ( m Λ b m Λ * ) + C 1 ( v ) m Λ * 2 + m Λ b 2 + q 2 2 m Λ b ] ζ 1 + m Λ * 2 q 2 m Λ b 2 ζ 1 SL + ζ 2 SL + 1 2 m Λ b 2 [ ( 2 m Λ * 2 m Λ * m Λ b m Λ b 2 2 q 2 ) ζ 3 SL ( m Λ * 2 + m Λ * m Λ b 2 m Λ b 2 q 2 ) ζ 4 SL ] } ,
A V ( + 1 / 2 , 1 / 2 , 1 ) = 2 s + 3 m Λ b m Λ * { s ( 1 + C 0 ( v ) ) ( ζ 2 ζ 1 ) + m Λ * ( ζ 1 SL + ζ 2 SL ) s 2 m Λ b ( ζ 3 SL + ζ 4 SL ) } ,
and for the axial vector current:
A A ( + 1 / 2 , + 3 / 2 , + 1 ) = 2 s m Λ b ( ζ 1 SL ζ 2 SL ) ,
A A ( + 1 / 2 , + 1 / 2 , 0 ) = 2 s 3 q 2 m Λ b m Λ * { s + m Λ b [ ( 1 + C 0 ( v ) ) ( m Λ * 2 + m Λ * m Λ b + q 2 ) + s 2 C 1 ( v ) ] ζ 2 + s + ( m Λ b m Λ * ) ( 1 + C 0 ( v ) ) + s 2 m Λ b C 1 ( v ) ζ 1 + m Λ b 2 ( m Λ b 2 m Λ * 2 q 2 ) s + s m Λ b 2 ζ 1 SL + ( m Λ * 2 m Λ b 2 + q 2 ) ζ 2 SL , s + 2 m Λ b 2 [ ( 2 m Λ * 2 3 m Λ * m Λ b + m Λ b 2 2 q 2 ) ζ 3 SL + ( m Λ * 2 3 m Λ * m Λ b + 2 m Λ b 2 q 2 ) ζ 4 SL ] }
A A ( + 1 / 2 , + 1 / 2 , t ) = 2 s + s 3 q 2 m Λ b m Λ * { 1 m Λ b [ ( m Λ * 2 + m Λ * m Λ b q 2 ) ( 1 + C 0 ( v ) ) + 1 2 ( m Λ b 2 m Λ * 2 + q 2 ) C 1 ( v ) ] ζ 2 + [ ( m Λ b + m Λ * ) ( 1 + C 0 ( v ) ) + m Λ b 2 m Λ * 2 + q 2 2 m Λ b C 1 ( v ) ] ζ 1 + m Λ * 2 q 2 m Λ b 2 [ ζ 1 SL ζ 2 SL ] + 1 2 m Λ b 2 ( 2 m Λ * 2 + m Λ * m Λ b m Λ b 2 2 q 2 ) ζ 3 SL + ( m Λ * 2 m Λ * m Λ b 2 m Λ b 2 q 2 ) ζ 4 SL ]
A A ( + 1 / 2 , 1 / 2 , 1 ) = 2 s 3 m Λ b m Λ * { s + ( 1 + C 0 ( v ) ) ( ζ 1 + ζ 2 ) + m Λ * ( ζ 1 SL ζ 2 SL ) + s + 2 m Λ b ( ζ 3 SL ζ 4 SL ) } .
In the case of the tensor current, the non-zero helicity amplitudes in the heavy quark limit are
A T ( + 1 / 2 , + 3 / 2 , + 1 ) = 2 s + m Λ b ( m Λ b m Λ * ) ζ 1 SL + m Λ * m Λ b m Λ * 2 + q 2 m Λ b ζ 2 SL ,
A T ( + 1 / 2 , + 1 / 2 , 0 ) = 2 q 2 s + 3 m Λ b m Λ * { s ( ζ 2 ζ 1 ) + m Λ * 2 + m Λ b 2 q 2 m Λ b ( ζ 2 SL + ζ 1 SL ) + s 2 m Λ b ( ζ 3 SL + ζ 4 SL ) } ,
A T ( + 1 / 2 , 1 / 2 , 1 ) = 2 s + 3 m Λ b m Λ * { s ( m Λ * + m Λ b ) ζ 1 + m Λ * 2 + m Λ * m Λ b q 2 m Λ b ζ 2 + m Λ * ( m Λ * m Λ b m Λ * 2 + q 2 ) m Λ b ζ 1 SL + m Λ * ( m Λ b m Λ * ) ζ 2 SL + s 2 m Λ b ( m Λ b + m Λ * ) ζ 3 SL + m Λ * 2 + m Λ * m Λ b q 2 m Λ b ζ 4 SL ,
and for the tensor axial current:
A T 5 ( + 1 / 2 , + 3 / 2 , + 1 ) = 2 s m Λ b ( m Λ b + m Λ * ) ζ 1 SL + m Λ * 2 + m Λ * m Λ b q 2 m Λ b ζ 2 SL ,
A T 5 ( + 1 / 2 , + 1 / 2 , 0 ) = 2 q 2 s 3 m Λ b m Λ * { s + ( ζ 1 + ζ 2 ) + m Λ * 2 + m Λ b 2 q 2 m Λ b ( ζ 1 SL + ζ 2 SL ) + s + 2 m Λ b ( ζ 3 SL + ζ 4 SL ) } ,
A T 5 ( + 1 / 2 , 1 / 2 , 1 ) = 2 s 3 m Λ b m Λ * { s + ( m Λ b m Λ * ) ζ 1 + m Λ * m Λ b m Λ * 2 + q 2 m Λ b ζ 2 + m Λ * ( m Λ * 2 + m Λ * m Λ b q 2 ) m Λ b ζ 1 SL + m Λ * ( m Λ * + m Λ b ) ζ 2 SL + s + 2 m Λ b ( m Λ b m Λ * ) ζ 3 SL + m Λ * m Λ b m Λ * 2 + q 2 m Λ b ζ 4 SL .

Appendix B. Relations with Lattice Form Factors

The definitions of the form factors here employed differ from other conventions in the literature. In particular, the translation with the Lattice determination in Ref. [22] is needed. I find
F 1 / 2 , t = 1 2 s 4 m Λ b m Λ * f 0 , F 1 / 2 , 0 = 1 2 s + 4 m Λ b m Λ * f + , F 1 / 2 , = 1 2 s + 4 m Λ b m Λ * f , F 3 / 2 , = 1 2 s + 4 m Λ b m Λ * f , G 1 / 2 , t = 1 2 s + 4 m Λ b m Λ * g 0 , F 1 / 2 , 0 = 1 2 s 4 m Λ b m Λ * g + , G 1 / 2 , = 1 2 s 4 m Λ b m Λ * g , G 3 / 2 , = 1 2 s 4 m Λ b m Λ * g , T 1 / 2 , 0 = s + 1 / 2 m Λ * 4 m Λ b h + , T 1 / 2 , = s + 1 / 2 m Λ * 4 m Λ b h , T 3 / 2 , = s + 4 m Λ b m Λ * ( m Λ b + m Λ * ) h , T 1 / 2 , 0 5 = s 1 / 2 m Λ * 4 m Λ b h ˜ + , T 1 / 2 , 5 = s 1 / 2 m Λ * 4 m Λ b h ˜ , T 3 / 2 , 5 = s 4 m Λ b m Λ * ( m Λ b m Λ * ) h ˜ .

Appendix C. Correlations between the Fit Parameters

The correlation matrix for the HQE parameters is reported in Table A1. The order of the various correlation coefficients is the same as in Table 1.
Table A1. Correlation matrix for the HQE parameters.
Table A1. Correlation matrix for the HQE parameters.
1 0.879 0.440 0.0458 0.0460 0.120 0.0619 0.363 0.337 0.0312
0.879 1 0.160 0.0109 0.0585 0.130 0.0936 0.325 0.343 0.121
0.440 0.160 1 0.00723 0.00712 0.0101 0.0465 0.211 0.139 0.218
0.0458 0.0109 0.00723 1 0.861 0.512 0.382 0.221 0.0500 0.611
0.0460 0.0585 0.00712 0.861 1 0.406 0.435 0.214 0.0603 0.707
0.120 0.130 0.0101 0.512 0.406 1 0.887 0.0941 0.119 0.783
0.0619 0.0936 0.0465 0.382 0.435 0.887 1 0.160 0.0287 0.871
0.363 0.325 0.211 0.221 0.214 0.0941 0.160 1 0.966 0.164
0.337 0.343 0.139 0.0500 0.0603 0.119 0.0287 0.966 1 0.0606
0.0312 0.121 0.218 0.611 0.707 0.783 0.871 0.164 0.0606 1

References

  1. Aaij, R.; Mazeliauskas, A. Search for lepton-universality violation in B+K++ decays. Phys. Rev. Lett. 2019, 122, 191801. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Aaij, R.; Adeva, B.; Adinolfi, M.; Ajaltouni, Z.; Akar, S.; Albrecht, J.; Alessio, F.; Alexander, M.; Ali, S.; Alkhazov, G.; et al. Test of lepton universality with B0K*0+ decays. J. High Energy Phys. 2017, 8, 55. [Google Scholar] [CrossRef] [Green Version]
  3. Aaij, R.; Adeva, B.; Adinolfi, M.; Affolder, A.; Ajaltouni, Z.; Akar, S.; Albrecht, J.; Alessio, F.; Alexander, M.; Ali, S.; et al. Test of lepton universality using B+K++ decays. Phys. Rev. Lett. 2014, 113, 151601. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Aaij, R.; Beteta, C.A.; Ackernley, T.; Adeva, B.; Adinolfi, M.; Afsharnia, H. Measurement of CP-Averaged Observables in the B0K*0μ+μ Decay. Phys. Rev. Lett. 2020, 125, 011802. [Google Scholar] [CrossRef]
  5. The Lhcb Collaboration. Angular analysis of the B0K*0μ+μ decay using 3 fb−1 of integrated luminosity. J. High Energy Phys. 2016, 2, 104. [Google Scholar] [CrossRef]
  6. LHCb Collaboration. Measurement of Form-Factor-Independent Observables in the Decay B0K*0μ+μ. Phys. Rev. Lett. 2013, 111, 191801. [Google Scholar] [CrossRef] [Green Version]
  7. Aebischer, J.; Altmannshofer, W.; Guadagnoli, D.; Reboud, M.; Stangl, P.; Straub, D.M. B-decay discrepancies after Moriond 2019. Eur. Phys. J. C 2020, 80, 252. [Google Scholar] [CrossRef] [Green Version]
  8. Algueró, M.; Capdevila, B.; Crivellin, A.; Descotes-Genon, S.; Masjuan, P.; Matias, J.; Novoa Brunet, M.; Virto, J. Emerging patterns of New Physics with and without Lepton Flavour Universal contributions. Eur. Phys. J. C 2019, 79, 714. [Google Scholar] [CrossRef]
  9. Hurth, T.; Mahmoudi, F.; Neshatpour, S. Implications of the new LHCb angular analysis of BK*μ+μ: Hadronic effects or new physics? Phys. Rev. D 2020, 102, 055001. [Google Scholar] [CrossRef]
  10. Ciuchini, M.; Fedele, M.; Franco, E.; Paul, A.; Silvestrini, L.; Valli, M. Lessons from the B0,+K*0,+μ+μ angular analyses. Phys. Rev. D 2021, 103, 015030. [Google Scholar] [CrossRef]
  11. Aaij, R.; Adeva, B.; Adinolfi, M.; Affolder, A.; Ajaltouni, Z.; Akar, S.; Albrecht, J.; Alessio, F.; Alexander, M.; Ali, S.; et al. Differential branching fraction and angular analysis of Λb0→Λμ+μ decays. J. High Energy Phys. 2018, 6, 115, Erratum in 2018, 9, 145. [Google Scholar] [CrossRef] [Green Version]
  12. Aaij, R.; Adeva, B.; Adinolfi, M. Angular moments of the decay Λb0→Λμ+μ at low hadronic recoil. J. High Energy Phys. 2018, 09, 146. [Google Scholar] [CrossRef] [Green Version]
  13. Detmold, W.; Meinel, S. Λb→Λ+ form factors, differential branching fraction, and angular observables from lattice QCD with relativistic b quarks. Phys. Rev. D 2016, 93, 074501. [Google Scholar] [CrossRef] [Green Version]
  14. Meinel, S.; van Dyk, D. Using Λb→Λμ+μ data within a Bayesian analysis of |ΔB|=|ΔS|=1 decays. Phys. Rev. D 2016, 94, 013007. [Google Scholar] [CrossRef] [Green Version]
  15. Blake, T.; Meinel, S.; van Dyk, D. Bayesian Analysis of bsμ+μ Wilson Coefficients using the Full Angular Distribution of Λb→Λ(→pπ)μ+μ Decays. Phys. Rev. D 2020, 101, 035023. [Google Scholar] [CrossRef] [Green Version]
  16. Bernet, R.; Muller, K.; Owen, P.; Serra, N.; Steinkamp, O.; LHCb Collaboration. Test of lepton universality with Λb0pK+ decays. J. High Energy Phys. 2020, 05, 040. [Google Scholar] [CrossRef]
  17. Aaij, R. Observation of J/ψp Resonances Consistent with Pentaquark States in Λb0J/ψKp Decays. Phys. Rev. Lett. 2015, 115, 072001. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Hiller, G.; Knecht, M.; Legger, F.; Schietinger, T. Photon polarization from helicity suppression in radiative decays of polarized Lambda(b) to spin-3/2 baryons. Phys. Lett. B 2007, 649, 152–158. [Google Scholar] [CrossRef] [Green Version]
  19. Mott, L.; Roberts, W. Rare dileptonic decays of Λb in a quark model. Int. J. Mod. Phys. A 2012, 27, 1250016. [Google Scholar] [CrossRef] [Green Version]
  20. Descotes-Genon, S.; Novoa-Brunet, M. Angular analysis of the rare decay Λb→Λ(1520)(→NK)+. J. High Energy Phys. 2019, 06, 136, Erratum in 2019, 6, 102. [Google Scholar] [CrossRef] [Green Version]
  21. Das, D.; Das, J. The Λb→Λ*(1520)(→NK¯)+ decay at low-recoil in HQET. J. High Energy Phys. 2020, 07, 002. [Google Scholar] [CrossRef]
  22. Meinel, S.; Rendon, G. Λb→Λ*(1520)ℓ+ℓ- form factors from lattice QCD. 2020. [Google Scholar]
  23. Falk, A.F. Hadrons of arbitrary spin in the heavy quark effective theory. Nucl. Phys. B 1992, 378, 79–94. [Google Scholar] [CrossRef]
  24. Böer, P.; Bordone, M.; Graverini, E.; Owen, P.; Rotondo, M.; Van Dyk, D. Testing lepton flavour universality in semileptonic Λb→Λc* decays. J. High Energy Phys. 2018, 6, 155. [Google Scholar] [CrossRef] [Green Version]
  25. Neubert, M. Heavy quark symmetry. Phys. Rept. 1994, 245, 259–396. [Google Scholar] [CrossRef] [Green Version]
  26. Grinstein, B.; Pirjol, D. Exclusive rare BK*+ decays at low recoil: Controlling the long-distance effects. Phys. Rev. D 2004, 70, 114005. [Google Scholar] [CrossRef] [Green Version]
  27. Böer, P.; Feldmann, T.; van Dyk, D. Angular Analysis of the Decay Λb→Λ(→Nπ)+. J. High Energy Phys. 2015, 1, 155. [Google Scholar] [CrossRef] [Green Version]
  28. Falk, A.F.; Neubert, M. Second order power corrections in the heavy quark effective theory. 2. Baryon form-factors. Phys. Rev. D 1993, 47, 2982–2990. [Google Scholar] [CrossRef] [Green Version]
  29. Falk, A.F.; Neubert, M. Second order power corrections in the heavy quark effective theory. 1. Formalism and meson form-factors. Phys. Rev. D 1993, 47, 2965–2981. [Google Scholar] [CrossRef] [Green Version]
  30. Leibovich, A.K.; Stewart, I.W. Semileptonic Lambda(b) decay to excited Lambda(c) baryons at order Lambda(QCD) / m(Q). Phys. Rev. D 1998, 57, 5620–5631. [Google Scholar] [CrossRef] [Green Version]
  31. Bernlochner, F.U.; Ligeti, Z.; Robinson, D.J.; Sutcliffe, W.L. New predictions for Λb→Λc semileptonic decays and tests of heavy quark symmetry. Phys. Rev. Lett. 2018, 121, 202001. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Comparison between the lattice results in Ref. [22] (grey band) and the fit results for the HQE form factors (red band) for the vector and axial-vector form factors. The two bands represent the 68 % interval.
Figure 1. Comparison between the lattice results in Ref. [22] (grey band) and the fit results for the HQE form factors (red band) for the vector and axial-vector form factors. The two bands represent the 68 % interval.
Symmetry 13 00531 g001
Figure 2. Comparison between the lattice results in Ref. [22] (grey band) and the predictions for tensor and pseudo-tensor form factors based on the fir results in Table 1 (red band). The two bands represent the 68 % interval.
Figure 2. Comparison between the lattice results in Ref. [22] (grey band) and the predictions for tensor and pseudo-tensor form factors based on the fir results in Table 1 (red band). The two bands represent the 68 % interval.
Symmetry 13 00531 g002
Figure 3. Comparison between the lattice QCD results in Ref. [22] (grey band) and the predictions for tensor and pseudo-tensor form factors based on the fir results in Table 1 (red band) with lattice QCD uncertainties inflated of 20 % . The two bands represent the 68 % interval.
Figure 3. Comparison between the lattice QCD results in Ref. [22] (grey band) and the predictions for tensor and pseudo-tensor form factors based on the fir results in Table 1 (red band) with lattice QCD uncertainties inflated of 20 % . The two bands represent the 68 % interval.
Symmetry 13 00531 g003
Table 1. Best fit points for the HQE parameters.
Table 1. Best fit points for the HQE parameters.
ParameterBest Fit Point
ζ 1 ( 0 ) 0.454 ± 0.070
ζ 2 ( 0 ) 0.0303 ± 0.0552
ζ 1 ( 1 ) + ζ 2 ( 1 ) 0.113 ± 0.024
ζ 1 SL , ( 0 ) 0.125 ± 0.038
ζ 1 SL , ( 1 ) 0.0487 ± 0.0614
ζ 2 SL , ( 0 ) 0.0110 ± 0.0363
ζ 2 SL , ( 1 ) 0.00362 ± 0.06184
ζ 3 SL , ( 0 ) 0.228 ± 0.190
ζ 4 SL , ( 0 ) 0.0883 ± 0.185
ζ 4 SL , ( 1 ) ζ 3 SL , ( 1 ) 0.0267 ± 0.0773
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bordone, M. Heavy Quark Expansion of Λb→Λ*(1520) Form Factors beyond Leading Order. Symmetry 2021, 13, 531. https://doi.org/10.3390/sym13040531

AMA Style

Bordone M. Heavy Quark Expansion of Λb→Λ*(1520) Form Factors beyond Leading Order. Symmetry. 2021; 13(4):531. https://doi.org/10.3390/sym13040531

Chicago/Turabian Style

Bordone, Marzia. 2021. "Heavy Quark Expansion of Λb→Λ*(1520) Form Factors beyond Leading Order" Symmetry 13, no. 4: 531. https://doi.org/10.3390/sym13040531

APA Style

Bordone, M. (2021). Heavy Quark Expansion of Λb→Λ*(1520) Form Factors beyond Leading Order. Symmetry, 13(4), 531. https://doi.org/10.3390/sym13040531

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop