The Λb→Λ(→pπ-)μ+μ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\varvec{\Lambda _b\rightarrow \Lambda (\rightarrow p\pi ^-)\mu ^+\mu ^-}$$\end{document} decay in the aligned two-Higgs-doublet model

The rare baryonic decay Λb→Λ(→pπ-)μ+μ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _b\rightarrow \Lambda (\rightarrow p\pi ^-)\mu ^+\mu ^-$$\end{document} provides valuable complementary information compared to the corresponding mesonic b→sμ+μ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b\rightarrow s\mu ^+\mu ^-$$\end{document} transition. In this paper, using the latest high-precision lattice QCD calculation of the Λb→Λ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda _b\rightarrow \Lambda $$\end{document} transition form factors, we study this interesting decay within the aligned two-Higgs-doublet model, paying particularly attention to the effects of the chirality-flipped operators generated by the charged scalars. In order to extract the full set of angular coefficients in this decay, we consider the following ten angular observables, which can be derived from the analysis of the subsequent parity-violating Λ→pπ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\Lambda \rightarrow p\pi ^-$$\end{document} decay: the differential branching fraction dB/dq2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{d}\mathcal{B}/\mathrm{d}q^2$$\end{document}, the longitudinal polarization fraction FL\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$F_\mathrm{L}$$\end{document}, the lepton-, hadron- and combined lepton–hadron-side forward–backward asymmetries AFBℓ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_\mathrm{FB}^\ell $$\end{document}, AFBΛ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_\mathrm{FB}^\Lambda $$\end{document} and AFBℓΛ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_\mathrm{FB}^{\ell \Lambda }$$\end{document}, as well as the other five asymmetry observables Yi\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$Y_i$$\end{document} (i=2,3s,3sc,4s,4sc\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$i=\mathrm{2,\,3s,\,3sc,\,4s,\,4sc}$$\end{document}). Detailed numerical comparisons are made between the SM and NP values for these angular observables. It is found that, under the constraints from the inclusive B→Xsγ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B\rightarrow X_s\gamma $$\end{document} branching fraction and the latest global fit results of b→sℓℓ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$b\rightarrow s\ell \ell $$\end{document} data, the contributions of right-handed semileptonic operators O9,10′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$O^{\prime }_{9,10}$$\end{document}, besides reconciling the P5′\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$P_5^{\prime }$$\end{document} anomaly observed in B0→K∗0μ+μ-\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$B^0\rightarrow K^{*0}\mu ^+\mu ^-$$\end{document} decay, could also enhance the values of dB/dq2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\mathrm{d}\mathcal{B}/\mathrm{d}q^2$$\end{document} and AFBℓ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$A_\mathrm{FB}^\ell $$\end{document} in the bin [15,20]GeV2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$[15,20]~\mathrm{GeV}^2$$\end{document}, leading to results consistent with the current LHCb measurements.


Introduction
The rare semileptonic b-hadron decays induced by the flavour-changing neutral current (FCNC) transition b → s + − do not arise at tree level and, due to the Glashowa e-mail: qyhu@mails.ccnu.edu.cn b e-mail: xqli@mail.ccnu.edu.cn c e-mail: yangyd@mail.ccnu.edu.cn Iliopoulos-Maiani (GIM) mechanism [1], are also highly suppressed at higher orders within the Standard Model (SM). In many extensions of the SM, on the other hand, new TeVscale particles can participate in the SM loop diagrams and lead to measurable effects in these rare decays. As a consequence, they play an important role in testing the SM and probing New Physics (NP) beyond it [2,3].
While no solid evidence of NP at all has been found in direct searches at high-energy colliders, it is interesting to note that several persistent deviations from the SM predictions have been observed in rare B-meson decays [3]. Specific to the b → s + − mesonic decays, these include the angular observable P 5 in the kinematical distribution of B 0 → K * 0 μ + μ − [4][5][6][7][8], the lepton-flavour-universalityviolation ratio R K of the decay widths for B → K μ + μ − and B → K e + e − [9][10][11], as well as the differential decay rates for B → K ( * ) μ + μ − [12][13][14] and B s → φμ + μ − [15][16][17]. Motivated by these anomalies and using the other available data on such rare mesonic decays, several global analyses have been made [18][19][20][21][22][23][24][25][26][27], finding that a negative shift in the Wilson coefficient C 9 improves the agreement with the data. However, due to the large hadronic uncertainties involved in exclusive modes, it remains quite unclear whether these anomalies indicate the smoking gun of NP, or are caused merely by underestimated hadronic power corrections [27][28][29][30][31][32][33][34] or even just by statistical fluctuations. In order to further understand the origin of the observed anomalies, it is very necessary to study other processes mediated by the same quark-level b → s + − transition.
In this respect, the rare baryonic b → μ + μ − decay is of particular interest for the following two reasons. Firstly, due to the spin-half nature of b and baryons, there is the potential to improve the currently limited understanding of the helicity structure of the underlying effective weak Hamiltonian [35][36][37]. Secondly, exploiting the full angular distribution of the four-body b → (→ pπ − )μ + μ − decay, one can obtain information on the underlying shortdistance Wilson coefficients of effective four-fermion operators, which is complementary to that obtained from the corresponding mesonic decays [38][39][40]. Experimentally, this decay was observed firstly by the CDF collaboration with 24 signal events and a statistical significance of 5.8 Gaussian standard deviations [41]. Later, the LHCb collaboration published the first measurements of the differential branching fractions as well as three angular observables of this decay [42]. As the b baryons account for around 20% of the b-hadrons produced at the LHC [43], refined measurements of this decay will be available in the near future. On the theoretical side, this decay is challenged by the hadronic uncertainties due to the b → transition form factors and the non-factorizable spectator dynamics [38,[44][45][46]. As the theory of QCD factorization at low q 2 [47,48] is not yet fully developed for the baryonic decay, we neglect all the non-factorizable spectator-scattering effects. For the factorizable nonlocal hadronic matrix elements of the operators O 1 -O 6 , O 8 , we absorb them into the effective Wilson coefficients C eff 7 (q 2 ) and C eff 9 (q 2 ) [47][48][49][50][51][52]. For previous studies of this decay, the reader is referred to Refs. .
Interestingly, it has been observed by Meinel and Dyk [81] that the b → (→ pπ − )μ + μ − decay prefers a positive shift to the Wilson coefficient C 9 , which is opposite in sign compared to that found in the latest global fits of only mesonic decays [22,26,27]. This suggests that a simple shift in C 9 alone could not explain all the current data and needs more thorough analyses. In our previous paper [82], we have studied the B 0 → K * 0 μ + μ − decay in the aligned two-Higgs-doublet model (A2HDM) [83], and found that the angular observable P 5 could be increased significantly to be consistent with the experimental data in the case when the charged-scalar contributions to C H ± 7 and C H ± 9,10 are sizeable, but C H ± 9,10 0. In order to further understand the anomalies observed in the b → s + − mesonic decays, in this paper, we shall study the b → (→ pπ − )μ + μ − decay in the A2HDM. As the b polarization in the LHCb setup has been measured to be small and compatible with zero [84], and the polarization effect will be averaged out for the symmetric ATLAS and CMS detectors, we consider only the case of unpolarized b decay. In order to reduce as much as possible the uncertainties arising from input parameters and transition form factors, we shall calculate all of the angular observables in some appropriate combinations [38][39][40]. For the b → transition form factors, we use the latest highprecision lattice QCD calculation [85], which is extrapolated to the whole q 2 region using the Bourrely-Caprini-Lellouch parametrization [86]. These results are also consistent with those of the recent QCD light-cone sum rule calculation [46], but with much smaller uncertainties in most of the kinematic range.
Our paper is organized as follows. In Sect. 2, we give a brief overview of the A2HDM. In Sect. 3, we present the theoretical framework for b → (→ pπ − )μ + μ − decay, including the effective weak Hamiltonian, the b → transition form factors, and the observables of this decay. In Sect. 4, we give our numerical results and discussions. Our conclusions are made in Sect. 5. Some relevant formulae for the Wilson coefficients are collected in the appendix.

The aligned two-Higgs-doublet model
We consider the minimal version of 2HDM, which is invariant under the SM gauge group and includes, besides the SM matter and gauge fields, two complex scalar SU (2) L doublets, with hypercharge Y = 1/2 [83,87]. In the Higgs basis, the two doublets can be parametrized as GeV is the nonzero vacuum expectation value, and G ± , G 0 are the massless Goldstone fields. The remaining five physical degrees of freedom are given by the two charged fields H ± (x) and the three neutral ones ϕ 0 with the orthogonal transformation R fixed by the scalar potential [83,87,88].
The most general Yukawa Lagrangian of the 2HDM is given by [83] a (x) are the charge-conjugated scalar doublets with hypercharge Y = − 1 2 ;Q L andL L are the lefthanded quark and lepton doublets, and u R , d R and R the corresponding right-handed singlets, in the weak-interaction basis. All fermionic fields are written as 3-dimensional vectors and the Yukawa couplings M f and Y f ( f = u, d, ) are therefore 3 × 3 matrices in flavour space. Generally, the couplings M f and Y f cannot be diagonalized simultaneously and the non-diagonal elements will give rise to unwanted tree-level FCNC interactions. In the fermion mass-eigenstate basis, with diagonal mass matrices M f , the tree-level FCNCs can be eliminated by requiring the alignment in flavour space of the Yukawa matrices [83]: where ς f ( f = u, d, ) are arbitrary complex parameters and could introduce new sources of CP violation beyond the SM.

Effective weak Hamiltonian
The effective weak Hamiltonian for b → s + − transition is given by [52] where G F is the Fermi coupling constant, and we have neglected the doubly Cabibbo-suppressed contributions to the decay amplitude. The operators O i≤6 are identical to P i given in Ref. [91], and the remaining ones read where e (g s ) is the electromagnetic (strong) coupling constant, T a the generator of SU (3) C in the fundamental representation, andm b denotes the b-quark running mass in the MS scheme. Within the SM, O 7,9,10 play the leading role in b → s + − transition, while the factorizable contributions from O 1−6,8 can be absorbed into the effective Wilson coefficients C eff 7 (q 2 ) and C eff 9 (q 2 ) [25]: 1,c (q 2 ) + C 8 F where the basic fermion loop function is given by [47] h(m q , and the functions F 2) of Ref. [47], while F (7,9) 1,c (q 2 ) and F (7,9) 2,c (q 2 ) are provided in Ref. [92] for low q 2 and in Ref. [93] for high q 2 . 1 The quark masses appearing in these functions are defined in the pole scheme. The contribution from O 7 is suppressed bym s /m b and those from O 9,10 are zero within the SM.
In the A2HDM, the charged-scalar exchanges lead to additional contributions to C 7,9,10 and make the contributions of chirality-flipped operators O 7,9,10 to be significant, through the Z 0 -and photon-penguin diagrams shown in Fig. 1. Since we have neglected the light lepton mass, there is no contribution from the SM W -box diagrams with the W ± bosons replaced by the charged scalars H ± . The new contributions to the Wilson coefficients read [82] 1 Here we incorporate only the leading contributions from an operator product expansion (OPE) of the nonlocal product of O 1−6,8 with the quark electromagnetic current, because the first and second-order corrections in /m b from the OPE are already well suppressed in the high-q 2 region [49,50]. Although non-factorizable spectator-scattering effects (i.e., corrections that are not described using hadronic form factors) are expected to play a sizeable role in the low-q 2 region [47,48], we shall neglect their contributions because there is presently no systematic framework in which they can be calculated for the baryonic decay [46]. As a consequence, our predictions in the low-q 2 region are affected by a hitherto unquantified systematic uncertainty.
.3) Fig. 1 Z -and photon-penguin diagrams involving the charged-scalar exchanges in the A2HDM with the functions C ( ) i,XY (i = 7, 9, 10; X, Y = u, d) given by Eqs. (A.1)-(A.10). Assuming ς u,d to be real, one has C H ± 7 =m s m b C H ± 7 , and we shall therefore neglect C H ± 7 in the following discussion.

Transition form factors
In order to obtain compact forms of the helicity amplitudes [38], we adopt the helicity-based definitions of the b → transition form factors, which are given by [38,44] for the vector and axial-vector currents, respectively, and for the tensor and pseudo-tensor currents, respectively. Here q = p − p and s ± = (m b ± m ) 2 − q 2 . The helicity form factors satisfy the endpoint relations f All these ten form factors have been recently calculated using (2 + 1)-flavour lattice QCD [85].
• The other five asymmetry observables which, along with the previous observables, determine all the ten angular coefficients K nλ . Here Y 2 also has a zero-crossing point, which lies in the low q 2 region.
In order to compare with the experimental data [97], we also consider the binned differential branching fraction defined by (3.22) and the binned normalized angular coefficients defined by where the numerator and denominator should be binned separately. As the theoretical calculations are thought to break down close to the narrow charmonium resonances, we make no predictions for these observables in this region. Finally, it should be noted that, unlike the strong decay K * → K π in the mesonic counterpart B → K * + − , the subsequent weak decay → pπ − is parity violating, with the asymmetry parameter α being known from experiment [98]. This fact distinguishes the signal with an intermediate baryon from the direct b → pπ − μ + μ − decay, and facilitates the full angular analysis of b → (→ pπ − )μ + μ − decay [38,39].

Input parameters
Firstly we collect in Table 1 the theoretical input parameters entering our numerical analysis throughout this paper. These include the SM parameters such as the electromagnetic and strong coupling constants, gauge boson, quark and hadron masses, 2 as well as the CKM matrix elements. The Weinberg mixing angle θ W is given by sin 2 θ W = 1 − M 2 W /M 2 Z . For the b → transition form factors, we use the latest high-precision lattice QCD calculation with 2 + 1 dynamical flavours [85]. The q 2 dependence of these form factors are parametrized in a simplified z expansion [86], modified to account for pion-mass and lattice-spacing dependences. All relevant formulae and input parameters can be found in Eqs. (38) and (49) and Tables III-V and IX-XII of Ref. [85]. To compute the central value, statistical uncertainty, and total systematic uncertainty of any observable depending on the form factors, such as the differential branching fraction and angular observables given in Eqs. (3.18)-(3.21), as well as the corresponding binned observables and the zerocrossing points, we follow the same procedure as specified in Eqs. (50)-(55) of Ref. [85].  Table 2 The Wilson coefficients at the scale μ b = 4.2 GeV, to NNLL accuracy in the SM

Results within the SM
For the short-distance Wilson coefficients at the low scale μ b = 4.2 GeV, we use the numerical values collected in Table 2, which are obtained at the next-to-next-to-leading logarithmic (NNLL) accuracy within the SM [3,91,100-103]. We show in Fig. 2 the SM predictions for the differential branching fraction and angular observables as a function of the dimuon invariant mass squared q 2 , where the central values are plotted as red solid curves and the theoretical uncertainties, which are caused mainly by the b → transform form factors, are labelled by the red bands. The latest experimental data from LHCb [97], where available, are also included in the figure for comparison. 3 The SM predictions for the corresponding binned observables are presented in Table 3.
As can be seen from Fig. 2 and Table 3, in the bin [15,20] GeV 2 where both the experimental data and the lattice QCD predictions for the b → transition form factors are most precise, the measured differential branching fraction [97], (1.20 ± 0.27) × 10 −7 GeV −2 , exceeds the SM prediction, (0.766 ± 0.069) × 10 −7 GeV −2 , by about 1.6 σ . Although being not yet statistically significant, it is interesting to note that the deviation is in the opposite direction to what has been observed in the B → K ( * ) μ + μ − [12][13][14] and B s → φμ + μ − [15][16][17] decays, where the measured differ- 3 For the differential branching fraction, the error bars are shown both including and excluding the uncertainty from the normalization mode b → J/ψ [98]. ential branching fractions favor, on the other hand, smaller values than their respective SM predictions. Also in this bin, the lepton-side forward-backward asymmetry measured by LHCb [97], −0.05 ± 0.09, is found to be about 3.3 σ higher than the SM value, −0.349 ± 0.013. As detailed in Ref. [81], combining the current data for b → (→ pπ − )μ + μ − decay with that for the branching ratios of B s → μ + μ − and inclusive b → s + − decays, Meinel and Dyk found that their fits prefer a positive shift to the Wilson coefficient C 9 , which is opposite in sign compared to that found in the latest global fits of only mesonic decays [22,26,27]. This means that a simple shift in C 9 alone could not explain all the current data. Especially, a negative shift in C 9 , as found in global fits of only mesonic observables, would further lower the predicted b → μ + μ − differential branching fraction.
Our SM predictions for the zero-crossing points of angular observables A FB , A FB and Y 2 read, respectively, (4.1) The zero-crossing points of the other observables Y i (i = 3s, 3sc, 4s, 4sc), which correspond to the case when the relative angular momentum between the pπ − system and the dimuon system is (l, m) = (1, ±1), are plagued by large theoretical uncertainties. The observables Y 3s and Y 3sc are predicted to be very small within SM and are, therefore, potentially good probes of NP beyond the SM [40].  Fig. 2 The b → (→ pπ − )μ + μ − observables as a function of the dimuon invariant mass squared q 2 , predicted both within the SM (central values: red solid curves, theoretical uncertainties: red bands) and in the A2HDM (case A: blue bands and case B: green bands). The corresponding experimental data from LHCb [97], where available, are represented by the error bars

Results in the A2HDM
In this subsection, we shall investigate the impact of A2HDM on the b → (→ pπ − )μ + μ − observables. For simplicity, the alignment parameters ς u,d are assumed to be real. As in our previous paper [82], we use the inclusive B → X s γ branching fraction [104,105] and the latest global fit results of b → s data [26,81] to restrict the model parameters ς u,d . Under these constraints, numerically, the chargedscalar contributions to the Wilson coefficients can be divided into the following two cases [82]: Case A : C H ± 7,9,10 are sizeable, but C H ± 9,10 0 ; Case B : C H ± 7 and C H ± 9,10 are sizeable, but C H ± 9,10 0.
They are associated to the (large |ς u |, small |ς d |) and (small |ς u |, large |ς d |) regions, respectively; see Ref. [82] for more details. This means that the charged-scalar exchanges contribute mainly to left-and right-handed semileptonic operators in case A and case B, respectively. The influences of these two cases on the b → (→ pπ − )μ + μ − observables are shown in Fig. 2, where the blue (in case A) and red (in case B) bands are obtained by varying randomly the model parameters within the ranges allowed by the global fits [26,81,82], with all the other input parameters taken at their respective central values.
In case A, the impact of A2HDM is found to be negligibly small on the hadron-side forward-backward asymmetry A FB and the observables Y i (i = 3s, 3sc, 4s, 4sc). For the differential branching fraction, on the other hand, visible enhancements are observed relative to the SM prediction, especially in the high q 2 region. For the remaining observables, the A2HDM only affects them in the low q 2 region, but the effect is diluted by the SM uncertainty. In order to see clearly the A2HDM effect in case A, we give in Table 4 the values of the binned observables in the bin [15,20] GeV 2 , including also the SM predictions, the A2HDM effect in case B, as well as the LHCb data (where available) for comparison. Although being improved a little bit, the deviations between the LHCb data and the theoretical values for the differential branching fraction and the lepton-side forward-backward asymmetry are still at 1.3 σ and 3.2 σ , respectively. Including the A2HDM in case A, there are only small changes on the zerocrossing points: In case B, however, the A2HDM has a significant influence on almost all the observables, as shown in Fig. 2. The most prominent observation is that it can enhance both the differential branching fraction and the lepton-side forward-backward asymmetry in the bin [15,20] GeV 2 , being now compatible with the experimental measurements at 0.2 σ and 1.3 σ , respectively (see also Table 4). The magnitude of the hadronside forward-backward asymmetry tends to become smaller in the whole q 2 region in this case, but is still in agreement with the LHCb data, with the large experimental and theoretical uncertainties taken into account. In the high (whole) q 2 region, a large effect is also observed on the asymmetry observable Y 3s (Y 4s ). Adding up the A2HDM effect in case B, the zero-crossing points are now changed to which are all significantly enhanced compared to the SM predictions (see Eq. (4.1)) and the results in case A (see Eq. (4.2)). It should be noticed that our predictions for the zero-crossing points given by Eqs. (4.1)-(4.3) are most severely affected by the hitherto unquantified systematic uncertainty coming from the non-factorizable spectatorscattering contributions at large hadronic recoil, a caveat emphasized already in Sect. 3.1.
Combining the above observations with our previous studies -the angular observable P 5 in B 0 → K * 0 μ + μ − decay could be increased significantly to be consistent with the experimental data in case B [82], we could, therefore, conclude that the A2HDM in case B is a promising alternative to the observed anomalies in b-hadron decays.

Conclusions
In this paper, we have investigated the A2HDM effect on the rare baryonic b → (→ pπ − )μ + μ − decay, which is mediated by the same quark-level b → sμ + μ − transition as in the mesonic B → K ( * ) μ + μ − decays. In order to extract all the ten angular coefficients, we have considered the differential branching fraction dB/dq 2 , the longitudinal polarization fraction F L , the lepton-, hadron-and combined lepton-hadron-side forward-backward asymmetries A FB , A FB and A FB , as well as the other five asym-metry observables Y i (i = 2, 3s, 3sc, 4s, 4sc). For the b → transition form factors, we used the most recent high-precision lattice QCD calculations with 2+1 dynamical flavours.
Taking into account constraints on the model parameters ς u,d from the inclusive B → X s γ branching fraction and the latest global fit results of b → s data, we found numerically that the charged-scalar exchanges contribute either mainly to the left-or to the right-handed semileptonic operators, labelled case A and case B, respectively. The influences of these two cases on the b → (→ pπ − )μ + μ − observables are then investigated in detail. While there are no significant differences between the SM predictions and the results in case A, the A2HDM in case B is much favored by the current data. Especially in the bin [15,20] GeV 2 where both the experimental data and the lattice QCD predictions are most precise, the deviations between the SM predictions and the experimental data for the differential branching fraction and the lepton-side forward-backward asymmetry could be reconciled to a large extend. Also in our previous paper [82], we have found that the angular observable P 5 in B 0 → K * 0 μ + μ − decay could be increased significantly to be consistent with the experimental data in case B. Therefore, we conclude that the A2HDM in case B is a very promising solution to the currently observed anomalies in b-hadron decays.
Finally, it should be pointed out that more precise experimental measurements of the full angular observables, especially with a finer binning, as well as a systematic analysis of non-factorizable spectator-scattering effects in b → (→ pπ − )μ + μ − decay, would be very helpful to further deepen our understanding of the quark-level b → sμ + μ − transition. and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .