Investigation of Effects of New Physics in $\Lambda_b\to\Lambda_c \tau\bar\nu_\tau$ Decay

Recent experimental results of ${\cal R}(D^{(*)})$ deviate from the standard model (SM) by $3.1\sigma$, suggesting a new physics (NP) that affects the $b\to c \tau \bar\nu_\tau$ transition. Motivated by this, we investigate the possible NP effects in the $\Lambda_b\to\Lambda_c \tau\bar\nu_\tau$ decay. For this purpose, assuming the neutrinos are left-handed, we calculate in detail the helicity amplitudes of $\Lambda_b\to\Lambda_c \ell\bar\nu_\ell$ decays with all possible four-fermion operators. Within the latest results of $\Lambda_b\to\Lambda_c$ form factors from lattice QCD calculations, we study these decays in a model-independent manner. The differential and total branching fractions and other observables are calculated. In SM, we obtain the ratio ${\cal R}(\Lambda_c)=0.33\pm0.01$. Supposing that NP only affects the third generation fermions, we present the correlations among ${\cal R}(D)$, ${\cal R}(D^*)$ and ${\cal R}(\Lambda_c)$. We perform a minimum $\chi^2$ fit of the wilson coefficient of each operator to the latest experimental data of different observables. It is found that the left-handed scalar operator ${\cal O}_{SL}$ affects the branching fraction remarkably, and the ratio ${\cal R}(\Lambda_c)$ can be enhanced by $30\%$. For other operators, the ratio amounts to $0.38\pm0.02$, which is larger than prediction of SM by $20\%$. Using the fitted values of the wilson coefficients of the single NP operators, we also give a prognosis for the physical observables of $\Lambda_b\to\Lambda_c \tau\bar\nu_\tau$. Furthermore, we also study the effects of three typical NP models on $\Lambda_b\to\Lambda_c \tau\bar\nu_\tau$. We hope our results can be tested in the current LHCb experiment and the future high energy experiments.


Introduction
After the discovery of Higgs boson, one of the most important tasks is searching for the possible new physics (NP) beyond the standard model (SM). So far, at the LHC people has not directly observed any signals of new particles beyond SM. However, there are currently a number of anomalies in semi-leptonic b-flavor hadrons decays, which have attracted a great attention in the scientific community. The unexpected deviations seem to appear in both b → c and b → s semi-leptonic decay transitions when different generations of leptons are involved, see refs. [1][2][3] for recent reviews. Unlike the flavor-changing neutral-current b → s transition that has been thought of as an ideal plate for probing the effect of NP, b → c transition is a tree level process by the exchange of a W boson in SM. Primordially, in order to study the charged Higgs contribution, the ratios are defined as [4,5]: For these two ratios, the uncertainties from the Cabibbo-Kobayashi-Maskawa (CKM) matrix element V cb can be removed, and ones from the hadronic transition form factors can also be reduced largely. In SM, based on the heavy-quark effective theory or the lattice QCD approach, they have been studied explicitly [6][7][8][9][10][11], and the theoretical averaged results are given by [12] R(D) SM = 0.299 +0.003 −0.003 , and R(D * ) SM = 0.258 +0.005 −0.005 .
However, from 2012, measurements from BaBar [13,14], Belle [15][16][17] and LHCb [18,19] showed that there are deviations between experimental data and predictions of SM, and the combination results performed by the Heavy Flavour Averaging Group (HFLAV) [12] are R(D) avg = 0.407 ± 0.039 ± 0.024 and R(D * ) avg = 0.306 ± 0.013 ± 0.007 , with a correlation of −53%, and the combined results compatible with SM at the 1.2σ level. Including these new measurements in the global average leads to [12] R(D) avg,new = 0.340 ± 0.027 ± 0.013 and R(D * ) avg,new = 0.299 ± 0.013 ± 0.008 , this work, we shall study the Λ b → Λ c τν τ comprehensively in model-dependence and independence manners, as well as Λ b → Λ c ℓν ℓ (ℓ = e, µ). In this work, we will calculate the total rates, the differential decay distributions, the longitudinal polarizations of the final states Λ c and τ , and the forward-backward asymmetries in the leptonside, with and without contributions of NP. Of course, the correlations between R(D ( * ) ) and R(Λ b ) will be given. Compared with previous studies, there are some improvements in the current work: (i) we calculate the complete amplitude including the contributions of scalar and tensor operators, and the latter was usually neglected in many previous studies; (ii) we present the relations between different conventions of from factors, which are nontrivial in the calculations ; (iii) the Λ b → Λ c form factors we adopted are the latest results from lattice QCD calculations; (iv) in discussing the effects of NP, we refit the each wilson coefficient individually within the latest experimental averaged results, and the form factors of B → D ( * ) in fitting are from the heavy quark effective theory including next-to-leading power O(1/m b ) and next-to-leading order O(α s ) corrections. Our paper is organized as follows: In Sec. 2, the analytical formulae in this work are presented, such as the effective hamiltonian, the form factors and the helicity amplitudes. All definitions of physical observables are also given in this section. In Sec. 3, we will present the numerical results and discussions. At last, we will draw conclusions in Sec. 4.

The effective hamiltonian
Within the presence of NP, the effective Lagrangian can be extended by incorporating new operator basis that includes all possible four-fermion interactions. Due to the absence of experimental evidence of deviations from the SM in tree-level transitions involving light leptons such as the precision measurements of the ratio B(τ → µν τνµ )/B(τ → eν τνe ) = 0.9762±0.0028 [94], the new physics is generally supposed to appear in the third generation fermions. Therefore, without considering the effects from the right-handed neutrinos, the effective Lagrangian of b → cτν τ can be written as where G F is the Fermi constant and V cb is the CKM matrix element. The four-fermion operators can be defined as where P L,R = (1 ∓ γ 5 )/2 and C i (i = V L, V R, SL, SR, and T ) are the corresponding wilson coefficients at the scale µ = m b , with C i = 0 in SM. It should be stressed that recent studies [95][96][97] based on the effective field theory of SM (SMEFT) [98,99] showed that the operator O V R does not contribute to LFU violation at leading order, however we here include this operator for completeness.

Form factors
In the calculation, the most important inputs are the hadronic transition form factors. The hadronic matrix elements of the vector and axial-vector currents between the two spin-half baryons Λ b and Λ c can be parameterized in terms of three form factors, respectively, as [83] Λ c , λ 2 |cγ µ b|Λ b , Of course, these two definitions are physically equivalent, and the relations between two groups of form factors are given as For the various helicity form factors we have used the formulae from the latest lattice QCD calculations [85][86][87], and each form factor can be written as where m f pole is pole mass, and f represents f +,⊥, 0 , g +,⊥, 0 , h +,⊥ and h +,⊥ , respectively. The explicit numerical values of m f pole , a f 0 , and a f 1 can be found in refs. [85][86][87], and we will not list them in the current work. Besides, the definition of the expansion parameter z is

The helicity amplitude
In SM, the semi-leptonic decay Λ b → Λ c ℓ −ν ℓ is considered to be Λ b → Λ c W − off−shell , and the off-shell W − off−shell decays to ℓ −ν l , subsequently. It is known to us that the W − off−shell has four helicities, namely λ W = ±1, 0 (J = 1) and λ W = 0 (J = 0), and only the off-shell W − off−shell has a timelike polarization, with J = 1, 0 denoting the two angular momenta of the rest frame W − off−shell . In order to distinguish the two λ W = 0 states we adopt the notation λ W = 0 for J = 1 and λ W = t for J = 0. In the Λ b -baryon rest frame, choosing the z-axis to be along the W − off−shell (see Fig. 1), we can write the polarization of the W − off−shell as where q µ is the four-momentum of the W − off−shell . In this case, the polarization vectors of the W − off−shell satisfy the orthonormality and completeness relations: ǫ * µ (m)ǫ µ (n) = g mn , m,n ǫ * µ (m)ǫ ν (n)g mn = g µν (25) with g mn = diag(+, −, −, −) for m, n = t, ±, 0. Because the current is not conserved in SM, so it consists of a superposition of a spin-1 and a spin-0 component where the J P content of the vector current J V µ and the axial vector current J A µ are (0 + , 1 − ) and (0 − , 1 + ), respectively. Therefore, we define helicity amplitudes as Due to the conservation of angular momentum, λ 1 = −λ 2 + λ W is satisfied. Wth above definitions, we preform the helicity amplitudes and have (see e.g. Refs. [101][102][103]) From parity or from an explicit calculation, we have If the right-handed current exists, we also have For the (S ∓ P )-type current, the corresponding helicity amplitudes are given by [93] H SP L ± 1 Also, one could define the hadronic helicity amplitudes of the (pseudo-)tensor-type current as and their explicit expressions are given by

Observable
With the above hadronic helicity amplitudes, we then write down the two-fold differential angular decay distribution as In above equation, | p 2 | = Q + Q − /(2M Λ b ) denotes the Λ c momentum in the Λ b rest frame. θ is the polar angle of the lepton, the definition of which can be seen in Fig.1. The auxiliary functions A i are given as When integrating out cos θ in eq. (34), we then obtain the differential decay rate dΓ(Λ b → Λ c ℓν ℓ )/dq 2 and the total branching fraction can be given as τ Λ b denoting the lifetime of Λ b baryon. Similar to R(D ( * ) ), one can define the differential and integrated ratios as Also, the forward-backward asymmetry in the lepton-side is defined as Besides, we can calculate the helicity-dependent differential decay rates, which are given as (normalized by the prefactor N/2) : Based on above helicities, we have two observables, the q 2 -dependent longitudinal polarizations of Λ c baryon and τ lepton, the definitions of which are given respectively as P Λc L (q 2 ) = dΓ λ2=1/2 /dq 2 − dΓ λ2=−1/2 /dq 2 dΓ λ2=1/2 /dq 2 + dΓ λ2=−1/2 /dq 2 ; In eqs. (36)(37)(38)(39)(40)(41)(42)(43)(44)(45), one finds that there are many quadratic terms of cos θ. In order to show up these terms, one defines a convexity parameter as 3 Results and Discussions

Input Parameters
In this section, we shall present the numerical results and discussions only for two cases ℓ − = µ − and τ − . Because the results for the electron are almost same as those of muon mode, we will not discuss the electron case any more. In our calculations, the most important non-perturbative parameters are the hadronic transition form factors that have been specified already in previous section. The other parameters we used are list here as

SM Predictions
Let us present the numerical results of SM firstly. Using the parameters, we calculate the branching fractions of Λ b → Λ c ℓν ℓ as and the ratio R(Λ b ) is given by In Figure.2, the q 2 -dependence of the differential ratios dΓ/dq 2 , the forward-backward asymmetries of the leptonic side A F B (q 2 ), the longitudinal polarization components of the Λ c and leptons (P h L (q 2 ) and P l L (q 2 )), and the convexity parameter C l F (q 2 ) are presented, respectively. As seen from Figure. 2, for the muon mode, the differential ratio shows a step-like behavior when q 2 = m 2 µ . At zero recoil q 2 = (M Λ b − M Λc ) 2 , A ℓ F B approaches to zero for both µ and τ modes. However, at large recoil, the A ℓ F B are −0.12 and −0.35, for both µ and τ modes, respectively. Furthermore, A µ F B is positive in most range, and the negative value arises from the small mass of muon. Due to the large mass of τ lepton, the A τ F B becomes negative when moving away from zero recoil, and it goes through zero point at q 2 = 7.93 GeV 2 . For the parameter C ℓ F , at zero recoil limit, both C µ F and C τ F are zero. At large recoil range, C µ F = −1.4 when q 2 = 0.4 GeV 2 , while C µ F changes to zero quickly when q 2 = m 2 µ due to the effect of mass. This behavior indicates that the cos θ distribution in q 2 ∈ [0.4, 11]GeV 2 is strongly parabolic. On the contrary, the C τ F is quite small in the whole range, which implies a straight-line behavior of the cos θ distribution. For the longitudinal polarization of the Λ c , it is found that at zero recoil the longitudinal polarizations are zero for two modes, while they approach to −1 in the large recoil regions. In the last panel, the q 2 dependence of the longitudinal polarization of the charged lepton are displayed. For the muon mode, the mass of muon lepton can be neglected in the Λ b decays, thus the muon lepton can be viewed as purely left-handed, and P µ L ≃ −1 in whole region. While for the tau mode, the behavior is quite different, P τ L = −0.52 at the zero recoil, and P τ L = 0 at the maximal recoil point.

Model-independence
We now present our discussions in the model-independence manner. Motivated by tensions stated in previous section, many works have been carried out including the effects of NP, by modifying the wilson coefficients or introducing new operators. As we known, in the quark level, Λ b → Λ c τ −ν τ , B → D ( * ) τ −ν τ , B c → J/ψ(η c )τ −ν τ and B − c → τ −ν τ are all induced by the b → cτ −ν τ transition. So, once we assume that R(D ( * ) ) were induced by NP, it should affect the observables of the Λ b → Λ c τ −ν τ decays. In Figure.3, we at first depict the correlations among the R(D), R(D * ) and R(Λ c ) in the presence of each single NP operator given in eq. (8). In the calculations, the form factors of B → D ( * ) we employed are from the heavy quark effective theory including higher power and higher order corrections [27]. The horizontal and vertical bands represent the experimental constraints at 1σ confidence level (C.L.) [12]. It should be emphasized that in the above discussions, we have not considered the constraints from the recent measurement of R(J/ψ) and the upper limit of branching fraction of pure leptonic decay B − c → τ −ν τ . In fact, Figure 2: The q 2 -dependence of differential ratios dBr/dq 2 , the forward-backward asymmetries of the leptonic side A F B (q 2 ), the convexity parameters C l F (q 2 ), and the longitudinal polarization components of the Λ c and leptons in SM. The red (solid) and blue (dashed) lines indicate muon mode and tau mode, respectively. once these constraints were added, all single-operator cannot accommodate the 1σ experimental constraints any more, which have been concluded in ref. [27]. In order to make a general model-independent analysis, the authors in ref. [27] had performed a minimum χ 2 fit of the wilson coefficients to the experimental data of all related observables, namely R(D), R(D * ), R(J/ψ), the longitudinal polarizations of τ lepton and final meson (D * ), with single-operator assumption. As aforementioned, the Belle collaboration released their new results of R(D ( * ) ) [20] at Moriond EW 2019. So, within the latest data, we shall refit the wilson coefficients including 2σ uncertainties, and the updated results are listed in Table.1. The technological details are refereed to ref. [27]. It is noted that in contrast to ref. [27] where all experimental results from Babar, Belle and LHCb were input, we here only input the five averaged results. With the new obtained wilson coefficients, we recalculate R(D), R(D * ) and R(Λ c ), and they are also presented in Table.1, from which one finds that all operators can enhance the R(Λ c ). Typically, for the operator O SL , the R(Λ c ) can be as large as 0.44.  With the new wilson coefficients, we also plot the q 2 -dependence of the differential branching fraction dB(Λ b → Λ c τν τ )/dq 2 and the ratio R(Λ c ) in Fig. 4. It is can be seen that these two observables are sensitive to  the NP operators. Moreover, the effects of the O SL operator are remarkable, and ones of other four operators are very similar, although all of them can enhance the branching fraction. Furthermore, we also show in Fig. 5 the q 2 -dependence of the Λ c and τ longitudinal polarizations as well as the lepton-side forward-backward asymmetry and the convexity parameter C l F (q 2 ). For the O V L operator, since the resulting wilson coefficients (1 + C V L ) appears in both the numerator and the denominator of P Λc,τ L and A τ F B , the contributions of NP effects are cancelled out by each other. As for other operators, due to the interference between SM and NP, the effects of NP can show up. The effects of O SL are significant. To be honest, there are certain uncertainties in the hadronic form factors. However, because the uncertainties from the form factor can be removed largely, these four diagrams can be used to be the probes of the new physics, if there observables can be measured precisely in future, such as in LHCb or the proposing circular electron-positron collider.

Models
In this section, we shall discuss some typical NP models that could affect b → cτν τ transition. As we stated previously, the intention of measurements of R(D ( * ) ) is to search for the charged Higgs-like scalars, such as H ± in 2HDMs. As known well, the charged Higgs boson contributes b → cτν τ at tree level and its effects is enhanced in some 2HDMs by introducing O SL and O SR operators. Unfortunately, these models are disfavored with the combined experimental results, see ref. [13]. Besides that, models with charged Higgses lead to unexpected large effects in the total B c lifetime, and can also be disfavored by the q 2 distribution. In these respects, we will not discuss the Higg-like models in this current work. Since the charged Higgs effects on b → cτν τ in the minimal supersymmetric standard model are the same as those in the 2HDM of type II at the tree level, we will not study it, either. As for the R-parity violating minimal supersymmetric standard model, it may affect b → cτν τ via mediating the slepton and down-squark. However, because its behavior is very similar to the leptoquarks that we will discuss later, we will not discuss this model any more. In the following, we mainly focus on the W ′ model, the scalar leptoquark (LQ) model and the vector LQ model as well.

W ′ Model
One simple way to obtain a new physics contribution to b → cτν τ is to use W ′ gauge boson which couples to the second-and third-generation fermions [104]. Without the right-handed neutrinos, the general Lagrangian that describes the couplings of an extra W ′ boson to quarks and leptons can be written as [105] with u i ∈ (u, c, t), d j ∈ (d, s, b) and ℓ i , ℓ j ∈ (e, µ, τ ). ǫ L,R uidj and ǫ L ℓiνj are the dimensionless coupling parameters that can be constrained by the current data of direct or indirect searches. In addition, for simplicity, all the effective couplings are set to be real. Because we hope the new physics only appear in the τ sector, so only the matrix element ǫ L τ ντ is assumed to be nonzero. If we transfer the eq.(58) to the effective Lagrangian as eq. (7), the wilson coefficients including the contribution of W ′ can be expressed as where M W ′ is the mass of the W ′ boson, and ǫ L,R cb and ǫ L τ ντ are the flavor-dependent couplings of eq. (58). As aforementioned, the O V R cannot contribute to LFU violation in the first order within SMEFT, so we only discuss the effect of the left-handed W ′ currently. In fact, when all current experimental constraints are taken into account, the permitted region of |ǫ R cb ǫ L τ ντ |/M 2 W ′ is barely reduced, though the authors of ref. [106] stated that the right-handed M W ′ could induce large CP asymmetry in B → D * µν decay mode.
To constrain the parameters in eq.(59), many works have been preformed by analysing the R(D ( * ) ) anomalies, such as in ref. [105,[107][108][109]. In ref. [107], within the current data of ATLAS [110] and CMS [111] the authors studied R(D ( * ) ) anomalies and the mono-tau signature at the LHC in the left-handed W ′ model and obtained ǫ L cb ǫ L τ ντ = (0.14 ± 0.03) with W ′ mass in the range [0.5, 3.5] TeV, which was consistent with the value ǫ L cb ǫ L τ ντ = 0.107 obtained in [108]. Very Recently, in ref. [105], the authors performed a χ 2 analysis within the latest experimental data and got the best fitted values ǫ L cb ǫ L τ ντ = 0.11 with M W ′ = 1TeV, which also agrees with previous studies. So, based on the above discussions, in order to maximize the effects of NP, we adopt the range as With above parameters, R(Λ c ) is predicted to be 0.38 ± 0.03. Furthermore, we show R(Λ c ) as a function of ǫ L cb ǫ L τ ντ in Figure.6, where the q 2 dependence of the differential branching fraction is also presented. It is found that the curves with an extra W ′ boson are on the above of ones from SM. Since we here only introduce a left-handed W ′ boson, it only affects the wilson coefficient C V L and cannot produce other operators, so that other observables, such as A F B , P h,l L and C l F , are unchanged, compared with results of SM. In addition, because we here suppose W ′ only couples to the third family fermions, it cannot be produced in the current protonantiproton collider, so it can escape from the current constrain of direct searches from ATLAS and CMS at LHC.

Scalar Leptoquark Model
For the b → cℓν ℓ , another potential tree-level mediators are leptoquaks. In the literatures, there are about ten possible leptoquark models that respect the SM symmetry, assuming the leptoquark be a scalar, a vector or a tensor particle. For the scalar particle, we here discuss a model [52], where the LQ φ transforms under the SM gauge group as (3, 1, − 1 3 ). It is in this model that not only the scalar and vector operators but the tensor operator can be generated. The part of the Lagrangian relevant to the couplings between LQ and SM fermions is given by [52] L where Q L , L denote the left-handed quark and lepton doublet, and u R , ℓ R denote the right-handed up-type quark and lepton singlet, respectively. λ L,R are the Yukawa coupling matrices in flavour space. More detailed discussions can be found in ref. [52]. Such a scalar LQ φ could induce the b → cτν τ transition with operators O V L , O SL and O T , and the corresponding wilson coefficients at scale µ = M φ are given by Because the decay we concerned is at scale µ = m b , all wilson coefficients should be run down to µ = m b scale. Due to the conservation of vector current, C V L is not renormalized. In the cases of C SL and C T , the explicit evolution equations and discussions could be found in refs. [68,112]. It should be emphasized that the model with such a scalar LQ can also accommodate the R K and (g − 2) µ anomalies. Taking M φ = 1TeV as a benchmark, in refs. [51,52,54,68], the authors had fitted the parameters within the experimental data, one point they favored is λ L bντ λ L * cτ = 0.35 and λ L bντ λ R * cτ = −0.03. Therefore, to maximize the effect of LQ, we set the parameter space as Using above parameters, we get R(Λ c ) = 0.41 ± 0.02, which is larger than SM prediction by 25%.By setting M φ = 1 TeV, we plot the changes of R(Λ c ) with λ L bντ λ L * cτ in Figure.7, and the differential branching fraction as a function of q 2 is also displayed. Figure 7: The prediction of R(Λ c ) as a function of λ L bντ λ L * cτ (left panel) and the differential branching fraction as a function of q 2 (right panel) with a scalar leptoquark. The black and red lines indicate the results of SM and the scalar leptoquark model, respectively. The yellow band is from the uncertainties of NP model.

Vector Leptoquark Model
Now, we move to discuss the effect of the vector leptoquark. In this article, we shall explore a model [53] where the introduced vector SU (2) L triplet LQ U µ 3 transforms as (3, 3, 2 3 ) under the SM gauge group. In the mass basis, the part of Lagrangian is given by where i, j = 1, 2, 3 denote the i, j-th generation fermion. The parameter g ij are the couplings of the Q = 2/3 component of the triplet, U (2/3) µ , tod Li and ℓ Lj , and V represents the CKM matrix. From the above Lagrangian, it is found that the b → cτν τ transition proceeds via exchange of LQ with the operator O V L . The corresponding wilson coefficient C V L is given correspondingly by It should be noted that in this model, due to the existence of the second term the R(D ( * ) ), R K ( * ) as well as the angular observable P ′ 5 in B → K * µ + µ − decay can be explained simultaneously with suitable parameter space, as shown in ref. [53]. From the results of χ 2 fits to the measured ratios R(D ( * ) ) and acceptable q 2 spectra done in ref. [54], we learn that at 1σ we have the following two best-fit solution (S1 and S2) It can be seen that the vector LQ only affects the functions A 1 and A V L 2 in eq.(34), and both two terms are proportional to |1 + C V L | 2 . Therefore, although the fit results S 1 and S 2 are quite different, the coefficients (1 + C V L ) they induced have nearly same absolute values with different signs. Currently, it is very hard for us to differentiate these two solutions. With above results, we obtain R(Λ c ) = 0.43 ± 0.03. Also, we plot changes of the R(Λ c ) with λ L bντ λ L * cτ , as well as the differential branching fraction as a function of q 2 in Figure.8. From two figures, we find that both the R(Λ c ) and the branching fractions are enhanced by about 30%. From above analyses, one finds that although the left-handed W ′ and the vector leptoquark particle origin form different models, however for the concerned decay mode, their behaviors are very similar, since both of them only contribute to the operator O V L .

Summary
Recent measurements of ratios R(D ( * ) ) imply that a new physics may appear in the b → cτν τ transition. So, in this work we explored the effect of NP in Λ b → Λ c τν τ decay mode that is also induced by b → cτν τ . To show up the NP effect, we firstly calculated the branching fraction of Λ b → Λ c τν τ and its ratio to Λ b → Λ c µν µ in SM, and R(Λ c ) = 0.33 ± 0.01, within the latest form factors from lattice QCD calculations. When studying the effect of NP, we explore Λ b → Λ c τν τ in a model-independent manner. The differential and total branching fraction, the longitudinal polarizations of final lepton, the forward-backward asymmetries in the lepton-side, the convexity parameters, and the ratio R(Λ c ) have been calculated. Supposing that NP only affects the third generation fermion, we also presented the correlations among R(D), R(D * ) and R(Λ c ). Considering the latest experimental data, we perform a minimum χ 2 fit of the wilson coefficient of each operator, and found that the left-handed scalar operator O SL could enhance the branching fraction by 30%. Using the fitted values of the wilson coefficients, we also give a prognosis for the physical observables of Λ b → Λ c τν τ , such as the ratio R(Λ c ), forward-backward asymmetry and other polarized observables as well as the differential branching fraction. At last, we also study the contributions of three typical NP models, the W ′ model, the scalar leptoquark (LQ) model and the vector LQ model, to the ratio and the differential branching fraction of Λ b → Λ c τν τ . These results can be tested in the current LHCb experiment and the future high energy experiments.