Semileptonic b-baryon decays within a new physics approach

The excess in the measurement of the branching fractions of B¯→D(D∗)τ-ν¯τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\bar{B}\rightarrow D(D^*){\tau }^- \bar{\nu }_{\tau }$$\end{document} and Bc→J/ψτ-ν¯τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B_{c}\rightarrow J/\psi {\tau }^{-} \bar{\nu }_{\tau }$$\end{document} from the standard model expectation hints the existence of new physics beyond the standard model. Motivated by these indications, we study similar modes mediated by b→cτν¯τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b\rightarrow c \tau {\bar{\nu }_{\tau }}$$\end{document} transitions, in particular, the semileptonic b-baryon decay modes Ξb→Ξcτ-ν¯τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varXi _b\rightarrow \varXi _{c} \tau ^- \bar{\nu }_{\tau }$$\end{document} and Σb→Σc(∗)τ-ν¯τ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varSigma _b \rightarrow \varSigma _c^{(*)} {\tau }^- \bar{\nu }_{\tau }$$\end{document}. We consider a general low energy effective Hamiltonian approach, which includes both standard model and new physics contributions. Within different new physics scenarios, we investigate the impact of the new contributions on these modes and present predictions for various semileptonic q2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$q^2$$\end{document}-spectra.


Introduction
The standard model (SM) of particle physics provides a unified framework of the fundamental particles and their interactions. Although the SM has been successful in describing a wide range of experimental measurements, it is not a complete theory. The limitations of the SM have led to various new physics (NP) beyond the SM searches via both direct and indirect means. NP has been hinted by precision measurements such as that of the muon anomalous magnetic moment (g − 2) μ . Flavor anomalies observed in b-hadron decays represent one of the indirect indications of NP beyond the SM. A number of measurements of b-decay observables have been found to disagree with SM predictions. These discrepancies have been seen particularly in decays mediated by b → s + − and b → cτ −ν τ transitions. One of the key properties of the SM is lepton flavor universality (LFU) which states that the couplings of leptons to gauge bosons are flavor independent. In the b sector, some of the observables that have been found to violate LFU are the ratio of branching fractions R D [8] also deviates from the SM prediction R SM J/ψ = 0.289 ± 0.01 [9] by about 2σ . These observations suggest the presence of NP beyond the Standard Model (SM) and motivate the study of other decays mediated by the same b → cτ −ν τ transitions. The τ polarization P D * τ = −0.38 ± 0.51 +0. 21 −0.16 [10,11] and the D * − polarization F D * L = 0.60 ± 0.08 ± 0.04 [12] as measured by the Belle collaboration also provide additional aspects to analyse NP in these transitions.
To explain the anomalies in the B meson modes, various NP models have been proposed, some of which can be found in [13][14][15][16][17][18][19][20][21]. In our work, we analyse the semileptonic b → cτ −ν τ transitions involving heavy b-baryons, c , in a model-independent approach. The semileptonic decays of b-baryons have not been studied as extensively as the b-meson ones though. Insights on weak interaction properties of the heavy baryons, the underlying dynamics and sensitivity to NP can be obtained from such decays. The semileptonic heavy b-baryon decays can also complement the sensitivity search of NP to that of the meson modes. The b-baryon decay modes which are of half-integer spin provide an auxiliary environment to test the observed anomalies. In addition, the CKM matrix parameter V cb can be determined from the semileptonic decays of heavy baryons mediated by b → cτ − ν τ transitions complimentary to heavy meson semileptonic decays. Decays involving heavy baryons are also good grounds for obtaining information on heavy quark physics. There is also less contamination from non-perturbative QCD effects with regards to semileptonic decays as leptons are involved. From the experimental side, experiments such as the Tevatron and LHC have made significant progress in the study of heavy baryons containing a b quark with a considerable accumulation of data. For the baryons of interest in this work, the Ξ 0 b and Ξ ± b baryons were first observed by CDF and D0 [22,23]. LHCb [24][25][26][27][28][29][30][31][32] and CDF [33] have also measured their masses, lifetimes and branching ratios. Clear signals of four strongly-decaying baryon states, have been obtained by CDF [34,35] and LHCb [36]. Since the Σ b baryon predominantly decays via the strong interaction, it may be difficult to measure its weak branching fraction. However, analysing the semileptonic decay modes may be insightful as they can be more sensitive to new particles and any deviation from the SM prediction will be a clear indication of NP [37].
The experimental progress in the study of heavy b-baryons which contain a single heavy quark necessitates theoretical progress. Important theoretical studies on b-baryons and their decays have indeed been carried out [38][39][40]. The semileptonic mode involving the spin-1/2 baryon Λ b has been considerably studied. The form factors of Λ b → Λ c mode have been computed using lattice QCD (LQCD) [41]. Exclusive semileptonic decays of spin- have been examined in the spectator quark model [42], the relativistic three-quark model [39] and a non-relativistic quark model [43]. The decay modes Λ b → Λ c −ν and Ξ b → Ξ c −ν were analysed in [44] within a light-front constituent quark model and similarly in [45] within a combined non-relativistic constituent quark model and heavy quark effective theory study. The authors in [46] studied the weak transition of Σ b → Σ c and Ω b → Ω c in the lightfront quark model considering the quark-diquark picture for heavy baryons. Similarly, in [37], the authors studied the Λ b → Λ c and Σ b → Σ c weak decays in the light-front quark model by considering the two spectator light quarks as individual ones rather than a diquark. In the work of [47], the semileptonic decays of heavy baryons with scalar and axial-vector diquarks in the relativistic quark model framework were analysed. The semileptonic decays framework of the relativistic quark-diquark model based on the quasipotential approach in [48]. In these works, SM estimates of decay rates, branching fractions, longitudinal and transverse asymmetries, the longitudinal to transverse decay ratio, the CKM parameter |V cb | and various other asymmetry parameters such as forward-backward asymmetry, convexity parameter, polarizations of the daughter baryon and the charged lepton, and ratio of branching fractions were determined.
There have also been some model-independent studies to probe NP in the heavy baryon decays. The Λ b → Λ c τ −ν τ mode has already been analysed with several NP signatures such as those that can be found in [49][50][51][52][53][54][55][56][57][58][59][60][61][62][63]. In [64], the mode Ξ b → Ξ c τ −ν τ decay was analyzed in the SM and in various NP scenarios with vector and scalar type of interactions. The observables R D and R D * were used to constrain the NP parameter space. The decay modes [65] within a model-independent effective field theory formalism. Real vector and scalar NP couplings were considered and the allowed NP parameter space was obtained from the experimental values of the observables R D ( * ) . In [66], new physics was probed in the baryon decays The new couplings assumed to be complex were constrained using the experimental measurements of B(B + c → τ + ν τ ), R l π , R D * , R J/ψ and F D * L . The authors in [67] studied the decay modes −ν in a model-independent effective field theory formalism. The NP couplings were assumed to be complex and the allowed NP parameter space was constrained from the experimental values of the observables B(B → D ( * ) −ν ), R D ( * ) , R J/ψ . In these analyses, predictions for various q 2 -dependent observables such as differential decay rate, branching fraction, forward-backward asymmetry of the charged lepton, convexity parameter, ratio of branching fractions, polarizations of the daughter baryon and the charged lepton were presented.
In our analysis of the heavy b-baryon semileptonic decay we consider a model-independent effective field theory framework, with a low energy effective Hamiltonian that includes both SM and NP contributions. For the amplitudes of the semileptonic decays of heavy b-hadrons, the hadronic matrix elements of weak currents are generally parametrized in terms of form factors that embody nonperturbative QCD effects. For b → c transitions, the heavy quark effective theory (HQET) [68][69][70][71][72][73][74][75] which is based on the 1/m Q expansion of the QCD Lagrangian and heavy quark symmetry, is the right approach for hadrons containing a single heavy quark. In HQET, the form factors can be simplified and described by the universal Isgur-Wise functions in the heavy quark expansion. Here, we analyze the decay modes without taking into account effects of radiative corrections to the HQET form factors. In our work, we use these form factors which have been obtained in the relativistic quark model (RQM) framework of [47], with the quark-diquark picture where a heavy baryon is seen as a bound state of a heavy quark and a light diquark system. We also use helicity amplitude formalism [76][77][78] to analyze the semileptonic transitions For the new couplings, we obtain constraints using the experimental measurements of R D ( * ) , R J/ψ , F D * L , P D * τ and the upper bound of B(B + c → τ + ν τ ). We then investigate the impact of these new contributions and present the predictions for various q 2 -dependent observables. It is to be noted that new physics sensitivity of the Σ b → Σ * c τ −ν τ mode has not been investigated so far. The helicity amplitudes in the presence of NP for this 1/2 + → 3/2 + transition has been worked out thoroughly and presented here.
The paper is organized as follows. In Sect. 2, we briefly discuss the theoretical framework used for the analysis, including form factors and helicity amplitudes. We also define various q 2 -dependent observables in this section. To assess the validity of the RQM framework, we first compare the theoretical predictions obtained using RQM with those obtained using lattice QCD results in Sect. 3. In Sect. 4, we present and discuss the NP sensitivity of the q 2 -dependent observables for the decay modes, considering one new coupling at a time and also using the global fits of [54]. We summarize our findings and conclude in Sect. 5.

Theoretical framework
We consider the most general low energy effective Hamiltonian relevant for b → c −ν transitions at the b-quark mass scale, with only left-handed neutrinos, given by [54] where G F is the Fermi constant and V cb is the CKM matrix element. The fermionic operators are given by and their corresponding Wilson coefficients denoted by C V L ,R , C S L ,R , C T are the vector, scalar and tensor type couplings that encode the NP contributions. The new contributions are assumed to be present only in the = τ mode as indicated by the LFU ratios. In our analysis, we consider only vector and scalar type of interactions and the NP couplings are assumed to be real. It has been shown that within the context of the SM effective field theory (SMEFT) [79][80][81], the vector operator with a right-handed quark current does not contribute to LFU violation [54,55,[82][83][84]. Hence, we do not include the effects of C V R in our work.

Hadronic matrix elements and form factors
The hadronic matrix elements of the vector and axial vector currents for the decays B b → B ( * ) c −ν can be parametrized in terms of various form factors expressed as functions of velocities of baryons [47]. For 1/2 → 1/2 transition, the parametrization is given by where where v and v are the four-velocities of the baryons B b and B c , respectively.
Using equations of motion, the hadronic matrix elements of the scalar and pseudoscalar currents can be obtained from those of the vector and axial vector currents and are given by where m b and m c denote the mass of b and c quarks respectively evaluated at the renormalization scale, μ = m b .
For 1/2 → 3/2 transition, the parametrization is given by where u B * c ,λ is the Rarita-Schwinger spinor for the B * c baryon. Again using equations of motion, the scalar and pseudoscalar matrix elements can be obtained as before and are given by In the heavy quark limit, the form factors for the semileptonic decay of Ξ b → Ξ c (scalar diquark picture) can be expressed as where the parameter Λ = (m B b − m b ) and ζ(ω) is the leading order Isgur-Wise (IW) function. The additional function χ(ω) arises from the 1/m b correction to the HQET Lagrangian. Near the zero recoil point of the daughter baryon, the functions ζ(ω) and χ(ω) can be approximated by where ρ 2 gives the slope and c gives the curvature of the IW functions. The values of these parameters are c (axial-vector diquark picture), the form factors in the heavy quark limit can be expressed in terms of the Isgur-Wise function ζ 1 (ω) as Again, near the zero recoil point of the daughter baryon, the ζ 1 (ω) function can be approximated by where ρ 2 ζ 1 gives the slope and c ζ 1 gives the curvature of the IW functions. The values of these parameters are ρ 2 ζ 1 = 2.17 and c ζ 1 = 3.62. Further details of the Isgur-Wise functions used in this work can be found in [47].

Helicity amplitudes
The helicity amplitudes are defined as [49,56,76,85,86] where λ 2 and λ W are the helicities of the daughter baryon and the virtual vector boson, respectively; * μ is the polarization vector of the virtual vector boson.
In the parent baryon (B b ) rest frame, the vector (V ) and axial vector ( A) helicity amplitudes can be obtained in terms of form factors and NP couplings. The scalar (S) and pseudoscalar (P) helicity amplitudes corresponding to new scalar and pseudoscalar interactions can be obtained using equations of motion.
For 1/2 → 1/2 transition, the helicity amplitudes are given by For 1/2 → 3/2 transition, the helicity amplitudes are given by The remaining helicity amplitudes can be obtained using parity relations The total helicity amplitudes are given by

q 2 -dependent observables
The twofold angular distribution for the decay including NP contributions can be written as [49] where θ is the angle of the lepton with respect to the W momentum in the rest frame of the W boson, |p B ( * ) , and    On integrating out cos θ from Eq. (16), the differential decay rate for the B b → B ( * ) c −ν decay including NP contributions can be obtained as where The differential branching fraction is given by Apart from these, other interesting q 2 -dependent observables are defined, such as:

Ratio of branching fractions
Forward-backward asymmetry of the charged lepton Convexity parameter Longitudinal polarization of the charged lepton where dΓ λ τ =±1/2 /dq 2 are helicity-dependent differential decay rates.

Comparison of RQM and LQCD predictions
In this section, we first confront the theoretical predictions obtained in the relativistic quark model framework [47] with those obtained using lattice QCD [41]. Since lattice results are available for the well-studied mode Λ b → Λ c τ −ν τ , we therefore present the comparison for form factors calculated in RQM and in LQCD for the semileptonic Λ b to Λ c decay at q 2 = 0 and q 2 = q 2 max in Table 1. Such a comparison can help to comprehend the different approaches used in the determination of form factors for baryons. From Table 1, it can be seen that the RQM form factor predictions and those from lattice results are in reasonable agreement at the zero recoil point of the daughter baryon. This is observed in the  [87] also. The relation between the RQM form factors and those used in lattice calculations are given in Appendix A.
In the RQM framework, the uncertainties coming from the form factors were estimated to be not more than 5% [48,88]. These uncertainties were shown to have their origin in the uncertainties of the baryon wave functions used in the calculation of the form factors. Using RQM form factors and considering a 5% uncertainty coming from the form factor inputs, we find that the predicted total decay rate divided by the square of the CKM matrix element is 91 ± 0.29) ps −1 and the relative error is 2.9%. With lattice form factors, the predic- Here, the relative error is about 2%. When the branching ratio is compared, we find that the prediction using lattice form factors within an effective Lagrangian approach (ELA) [50] is smaller compared to that obtained using RQM. Predictions for the ratio in RQM are somewhat smaller compared to lattice predictions. However, it should be noted that the sensitivity of these ratios to form factor uncertainties is less as these cancel out partially in these ratios. The predictions compared here are presented in Table 2.
In Fig. 1, we also display the differential branching fraction prediction for Λ b → Λ c τ −ν τ decay with RQM and LQCD calculations within the SM. It can be observed that there is a sizeable difference in the differential distribution and this may be attributed to the variation in the q 2dependence of the form factors. However, close to the point of zero recoil, the two predictions are in agreement with each other.

New physics sensitivity
We now proceed to study the NP sensitivity of the various q 2 -dependent observables defined in the earlier section for the modes The input parameters used for our numerical analysis are listed in Table 3.
For the SM calculation of decay amplitudes, the theoretical uncertainties that arise from the input parameters such as the CKM matrix element V cb and the form factors are taken into consideration. We consider the V cb uncertainty as given Fig. 1 Comparison of the differential branching fraction for Λ b → Λ c τ −ν τ decay in the SM obtained using RQM framework (red) and lattice calculations (green) in Table 3 and a 5% uncertainty coming from the form factor inputs as determined in [48,88]. We constrain the new couplings C V L , C S L , C S R using the experimental measurements of R D ( * ) , R J/ψ , F D * L and P D * τ . A 30% constraint is also imposed from the upper bound of B(B + c → τ + ν τ ) [90]. The relevant relations pertaining to these observables are given below.
The differential decay rate ofB → Dτ −ν τ is given by [19] dΓ where the hadronic helicity amplitudes H V,0 , H V,t and H S are expressed in terms of the Caprini et al. parametrized HQET form factors [91]. These form factors are evaluated using parameters obtained from lattice QCD calculations [92].  The differential decay rate ofB → D * τ −ν τ is given by [19,93] where H V,0 , H V,± , H V,t and H S are the hadronic helicity amplitudes. Here λ D ( * ) ( We useB → D * HQET form factors parametrized by Caprini et al. [91] where the fitted parameters are determined by HFLAV [1]. The differential decay rate expression for B c → J/ψτ −ν τ is similar to that ofB → D * τ −ν τ with appropriate substitutions for masses and form factors. For this mode, we use the form factors obtained in [93,94] by employing perturbative QCD approach. The longitudinal polarization of τ (P D * τ ) and of D * (F D * L ) in B → D * τ −ν τ are respectively given by where λ τ and λ D * denote the helicity of τ and D * , respectively.
The branching fraction of B + c → τ + ν τ is given by We use f B c = 489 ± 4 ± 3MeV from [95] and the values of other constants are taken from [89]. In our work, we consider one new coupling at a time and the best-fit values are obtained by performing a χ 2 fitting. The χ 2 function is defined as where O ex p i are the experimental values of the observables, O th i (C k ) are the theoretical predictions for the observables with new couplings C k , and C is the covariance matrix which takes into account the correlation of R D and R D * . The χ 2 function is minimized to obtain the best-fit values for each NP coupling. The NP effects are obtained by imposing a 1σ constraint from the measured values of R D ( * ) , R J/ψ , F D * L and P D * τ . The experimental and SM theoretical values of the observables used for obtaining the constraints are listed in Table 4. The obtained best-fit values are listed in Table 5.
The q 2 -spectra of various observables for the c τ −ν τ decay modes in the SM case and in the presence of the NP couplings C V L , C S L , C S R Fig. 2 The q 2 -dependence of various observables for the Ξ b → Ξ c τ −ν τ decay mode in the presence of vector (C V L ) and scalar (C S L , C S R ) NP couplings are presented in Figs. 2, 3 and 4. The predictions are distinguished by blue (SM), orange (C V L ), purple (C S L ) and red (C S R ) colors, respectively. For the mode Ξ b → Ξ c τ −ν τ , the differential branching ratio D B R deviates more noticeably from the SM in the presence of the C V L coupling than in the case with the scalar couplings. The differential decay rate dΓ /dq 2 is enhanced over the whole q 2 region in the presence of both the vector and scalar couplings for decay mode, dΓ /dq 2 is enhanced with respect to the vector coupling but is largely unaffected by the scalar couplings. In the case of and all the other observables behave SM-like for all three decay modes, as the NPdependent factor, (1 + C V L ) 2 cancels out in these quantities. R Ξ c (Σ c ) displays a distinctive deviation from the SM prediction in the higher q 2 region for C S L and C S R couplings. Measuring R Ξ c (Σ c ) in the higher q 2 region may thus further substantiate the observed anomalies in b-decays. No sizable deviation is observed for R Σ * c in the presence of scalar couplings. The forward-backward asymmetry A τ F B has a SM zero-crossing point at q 2 ≈ 8 GeV 2 for Ξ b → Ξ c τ −ν τ . This shifts to a slightly higher q 2 value in the presence of scalar couplings. There is no significant deviation from the SM prediction in the presence of scalar couplings. For the Σ b → Σ * c τ −ν τ decay mode, A τ F B has a SM zero-crossing point at q 2 ≈ 6.6 GeV 2 . The zero-crossing shifts slightly to a lower q 2 value with C S L coupling and to a higher q 2 value with C S R coupling. The convexity parameter C τ F displays a distinctive deviation from the SM prediction in case of C S L and C S R couplings for the decay modes considered, with a more pronounced deviation in the case of the Σ b → Σ c τ −ν τ mode. The measurement of this observable may further endorse existence of NP beyond the SM. The longitudinal polarization P τ L (q 2 ) has a zero-crossing at q 2 ≈ 4.6 GeV 2 in the SM decay modes using Fits 1 and 2 of [54] for Ξ b → Ξ c τ − ν τ decay. It shifts to a higher q 2 value in case of both the scalar couplings. The deviation from the SM prediction is more prominent at large q 2 . For Σ b → Σ c τ − ν τ decay, P τ L (q 2 ) shows a distinct deviation from the SM prediction in the case of the scalar couplings. Thus, measuring P τ L (q 2 ) may provide insight into the nature of NP and distinguish between the vector and scalar type contributions. For Σ b → Σ * c τ − ν τ , a SM zero-crossing point is observed at q 2 ≈ 3.5 GeV 2 . This shifts slightly to a lower q 2 value for C S L coupling and to a higher q 2 value for C S R coupling.
In addition to our results, we also present the NP predictions obtained using the Wilson coefficients from the global fits (Fits 1 and 2) of Murgui et al. [54] for all three decay modes in Fig. 5. Performing a global χ 2 fit to the data available in b → cτ −ν τ transitions, the authors of [54] obtained three fits, Fit 1, 2 and 3, where the coefficient C V R was assumed to be lepton-flavour universal also. The experimental measurements of the ratios R D , R D * and F D * L were used to constrain the new couplings. An upper bound constraint from B(B + c → τ + ν τ ) was also imposed. Here, we show the results using the numerical values of C V L , C S L and C S R from Fits 1 and 2 tabulated in Table 3 of [54] In this paper, we do not consider Fit 3 as it is found to be already disregarded by the measured differential distributions. For all the three decay modes, NP scenarios using Fits 1 and 2 are clearly distinguishable. It is observed that predictions with Fit 1 are much closer to the SM than those with Fit 2, which is similar to the observations of Murgui et al. for b-meson modes.

Summary and conclusion
We have analyzed the semileptonic b-baryon decays Ξ b → Ξ c τ −ν τ and Σ b → Σ ( * ) c τ −ν τ within a model-independent effective theory framework, consisting of both SM and NP contributions. The helicity amplitudes method was used to describe the hadronic transitions. Form factors expressed in terms of Isgur-Wise functions in the heavy quark limit have been adopted to parametrize the decay amplitudes. These form factors were obtained within the relativistic quark model. The SM predictions in the relativistic quark model and those obtained using lattice QCD results were first compared for an assessment of the validity of the different approaches. It was found that form factors calculated in RQM and on the lattice are compatible at the zero recoil point of the final baryon. The total decay rate, branching ratios and R Λ c differed slightly in the two cases. On analyzing NP effects, we considered one new coupling at a time and using a χ 2 fitting, best-fit values of these couplings were obtained. Predictions for the q 2 -dependence of various observables such as the differential branching ratio, ratio of branching fractions, forward-backward asymmetry of the charged lepton, convexity parameter and lepton polarization have been presented in both the SM and NP scenarios. In order to determine the allowed NP parameter space, we imposed a 1σ constraint from the experimental measurements of the LFU ratios R D ( * ) , R J/ψ , the D * − polarization F D * L and the τ polarization P D * τ . We observed NP sensitivity for most of the observables of interest, with the deviations from the SM prediction being more discernible in the scenario with new scalar couplings than that with new vector couplings. In particular, observables such as the ratio of branching fractions, convexity parameter and lepton polarization were found to be more sensitive to the new contributions for some of the decay modes considered. Using previously obtained significant fits of the new couplings, predictions for the different q 2dependent observables were also presented. It was observed that NP with these fits mostly displayed distinct variations from the SM. Hence, the study of semileptonic b → c −ν transitions of half-integer spin b-baryons such as those considered in this work can be invaluable in revealing the presence and ascertaining the exact nature of NP. Also, the study of such baryon modes augments that of the b-meson modes as they furnish a complementary environment in the search for NP.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: Data used in the manuscript has been cited and can be found in the references section. All results are given within the main text itself. There is no other associated data.] Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecomm ons.org/licenses/by/4.0/. Funded by SCOAP 3 . SCOAP 3 supports the goals of the International Year of Basic Sciences for Sustainable Development.

A Form factors for the semileptonic Λ b → Λ c decay
The form factors used in Eq. (3) are related to those used in [41] by the following relations.