Investigating R D and R D ∗ anomalies in a left–right model with an inverse seesaw

We investigate R D and R D ∗ anomalies in a low scale left–right symmetric model based on SU ( 3 ) C × SU ( 2 ) L × SU ( 2 ) R × U ( 1 ) B − L with a simpliﬁed Higgs sector consisting of only one bidoublet and one SU ( 2 ) R doublet. The Wilson coefﬁcients relevant to the transition b → c τν are derived by integrating out the charged Higgs H ± boson, which gives the dominant contributions. We emphasize that the charged Higgs effects, with the complex right-handed quark mixing matrix, can account for both R D and R D ∗ anomalies simultaneously, while adhering to a set of significant constraints including, for instance, BR ( B − c → τ − ¯ ν τ ) and B s ( d ) − ¯ B s ( d ) mixing. In relation to this, we show that the predicted values of the D ∗ , τ longitudinal polarizations and P τ ( D ) can be affected for the set of the parameters of the model resolving the R D , D ∗ anomalies.

Recently, with the announcement of the LHCb collaboration new results of R D and R D * in autumn 2022 using the Run 1 data sample of a luminosity of 2 f b −1 , the HFLAV obtained a new averages of R D and R D * for 2023 as [11], R D * = 0.284 ± 0.013.(5) The corresponding averages of SM predictions as reported by HFLAV [10] based on the predictions obtained in Refs.[12][13][14] are given by R SM D * = 0.254 ± 0.005.(7) Clearly, the listed results of R D and R D * in Eqs.(2,3) exceed the SM predictions given above, by 1.4σ and 2.9σ respectively.Moreover, as reported by HFLAV [10], the difference with the SM predictions reported above, corresponds to about 3.4σ [R D − R D * ] combination.In Ref. [15], the authors used all the available experimental data and took the higher derivative QCDSR constraints on HQET parameters into account to fit B q → {D * , D} form factors, and gave the accurate SM prediction of R D * and R D .The result of their investigation indicated a 4σ discrepancy in R D * and R D .All this may serve as a hint for a violation of lepton-flavor universality.
In the literature, many theoretical models have been investigated aiming to explain the R D and R D * anomalies.Based on the remark that the measured ratio R D * is higher than its value in the SM, attempts in many models have focused on enhancing the rate of b → cτ ν transition through extra contributions arising from the new particles mediating the transition.This turns to be much easier than reducing the rate of b → c(e, μ)ν transitions, given the much more severe restraints on the couplings of new physics to muons and electrons [16].
The new particles mediating the b → cτ ν can be spin-0 or spin-1.In addition, they can either carry baryon and lepton number (leptoquarks) or be B/L neutral (charged Higgs and W ). As reported in Ref. [16], the existing models can be classified into three main categories: • Models with charged Higgs [17][18][19][20][21][22].In these class of models, integrating out the charged Higgs, that mediates b → cτ ν transition at tree-level results in scalarscalar operators contributing to decay rates.The contributions of these additional operators can be constrained using the measured B c lifetime.This leads to the upper limits BR(B − c → τ − ντ ) ≤ 40% [23], BR(B − c → τ − ντ ) ≤ 30% [24] and a much stronger bound BR(B − c → τ − ντ ) ≤ 10% was obtained from the LEP data taken at the Z peak [25].Later on, the bounds using the measured B c lifetime had been critically investigated and relaxed upper limit of ≤ 39% [26] and ≤ 60% [27,28] were obtained.An analysis of the LHC sensitivity of the charged Higgs that can explain R D and R D * anomalies showed that the bound 400 GeV < m H ± is more stringent than B c lifetime bound [29].Based on the finding of Ref. [29] and other investigations carried out later in Refs.[30,31], the available charged Higgs mass range for the explanation of the R D and R D * anomalies within 1 σ range is bounded from the above as m H ± ≤ 400 GeV.It should be remarked that m H ± > 400 GeV is ruled out by the τ ν resonance search at the LHC [29,32], the low-mass bottom flavored di-jet search [33,34] and a conventional search for tau sleptons [35] set constraint on the available parameter region.Recently, the analysis has been extended to cover the mass range 180 GeV < m H ± < 400 GeV [30,36].Experimentally, searches for charged Higgs bosons at high-energy colliders such as the Large Hadron Collider (LHC) have set model-dependent limits on their masses.ATLAS and CMS collaborators have searched for charged Higgs bosons that have large couplings only to the third generation, looking for multijet events with one lepton and at least 2 b-jets.They have carried out dedicated searches only for charged Higgses that dominantly couple to tb quarks.In addition, there is no dedicated searches looking for a charged Higgs with sizable couplings to both q b and t b quarks [37].
The production of charged Higgs bosons is influenced by the particle's mass and its couplings to other particles, which are determined by the specific beyond the Standard Model (BSM) scenario being considered.When the charged Higgs boson is heavier than the top quark, the experimental constraints on its mass become weaker.Certain regions of the charged Higgs boson mass parameter space have been excluded by experimental searches, even for masses higher than the top quark mass.However, these constraints depend on the couplings to other particles, particularly top and bottom quarks.We have ensured that our coupling parameter, g H − t b, is below the limits obtained from previous model-independent analyses in Ref. [37].Therefore, the mass of the charged Higgs boson can fall within the range of 200 GeV to 1 TeV.
• Models with heavy charged vector bosons [38,39].In such class of models, integrating out W 's yields new contributions to the vector-vector operators.Simultaneous explanation of both R D and R D * with left-handed neutrinos requires non zero values for the associated new contributions.It should be noted that, these class of models are subjected to constraints originating from the presence of an accompanying Z mediator.Thus, the vertex Z b L s L can be inferred from the W b L c L vertex through the SU (2) invariance.The Z b L s L term in the Lagrangian can lead to tree-level flavor-changing neutral currents (FCNCs) and hence there is a need to some mechanism to subdue the contributions of this term.One way to do that is to assume minimal flavor violation (MFV) as considered in Refs.[40,41].However, adopting this choice in general can not suppress Z bb and Z τ τ vertices.Consequently, LHC direct searches for Z → τ τ resonances can result in stringent restraints on the related parameters of the models.However, as shown in Refs.[41,42], selecting unnaturally high Z widths helps to avoid these severe constraints.• Models with leptoquarks [18,43].Either scalar and vector leptoquarks can be potentially favored as an explanations for the R D ( * ) anomaly [42,44].These models are subjected also to many important constraints including those from b → sνν restraints [42], searches at the LHC for τ τ resonances [41,42], and the measured B c life-time [24,25].However, these constraints turn to be weaker compared to the corresponding ones imposed on the alternative models described above [16].
In this study, we aim to derive the effective Hamiltonian resulting from the four fermion transition b → cτ ν in the minimal left-right model with an inverse-seesaw (LRIS).At tree-level, the transition is mediated by the exchanging of the charged Higgs boson in the model.The large neutrino Yukawa, which measures the strength of the charged Higgslepton interaction, and the charged Higgs-quark interaction, which is proportional to right-handed mixing matrix V C K M R , are notable features of this class of left-right model with an inverse seesaw mechanism that motivate it to solve R D and R D * .
Having the effective Hamiltonian, we will proceed to show the dominant Wilson coefficients contributing to the ratios R D and R D * .Consequently, we will display the dependency of the ratios on the parameters of the model to determine the most relevant ones.Moreover, we reexamine the relation between the observable R D * and BR(B c → τ ντ ), reported in Refs.[24,25], in light of Belle combination [6] and LHCb [8] experimental results, and the projection at Belle II experiment with an integrated luminosity of 50 ab −1 [45].This turns to be very important because the constraint on BR(B c → τ ντ ) affects substantially the contributions from scalar operators [26][27][28].The B s(d) − Bs(d) mixing can lead to possible constraints on the parameter space related to the the ratios R D and R D * .These constraints will be investigated in this study also.Finally, we will also investigate the possibility of resolving the anomalies after including the new contributions to the Wilson coefficients of the total effective Hamiltonian in the presence of the LRIS model and give predictions of the D * and τ longitudinal polarizations.
This paper is organized as follows.In Sect.2, we give a brief review of the left-right model with inverse seesaw.In the review, we discuss the gauge structure and the particle content of the model.We also discuss the fermion interactions with charged Higgs and with the charged W, W gauge bosons.These interactions are required to derive the effective Hamiltonian governs the processes contributing to R D and R D * .Then, in Sect.3, we present the effective Hamiltonian describing b → cτ ν transition in the presence of a general New Physics (NP) beyond the SM.Particularly, we derive the analytic expressions of the Wilson coefficients up to one loop level originating from the charged Higgs mediation in LRIS model under study in this paper.Based on the Hamiltonian, we show the total expressions of R D and R D * .In Sect.4, we give our estimation of R D , R D * , the D * and τ longitudinal polarizations and their dependency on the parameter space.Finally, in Sect.5, we give our conclusion.

Left-right model with inverse seesaw (LRIS)
The gauge sector of the minimal left-right model with an inverse-seesaw (LRIS) is based on the symmetry The particle content of the model includes the usual particles in the SM in addition to extra fermions, scalars and gauge bosons.Regarding the fermion content of the model, it is the same as its counterpart in the conventional left-right models [46][47][48][49][50][51][52][53].The scalar sector consists of only one Bidoublet φ and one doublet χ R , which is more minimal and able to circumvent the stringent FCNC constraints placed on the conventional leftright model.Additionally, we abandon using the left-handed doublet in order to prevent the severe fine tuning in the VEV of the neutral component of the left-handed doublet scalar, v L , caused by light neutrino masses.As a result, our model is SU (2) L × SU (2) R gauge invariant but not symmetric in left-right parity.
The implementation of the IS mechanism for neutrino masses can be carried out via introducing three SM singlet fermions S 1 with B − L charge = −2 and three singlet fermions S 2 with B − L charge = +2.It should be remarked that, this particular choice of the B − L charges of the pair S 1,2 assures that the U (1) B−L anomaly is free.
The scalar sector of the LRIS model contains SU (2) R scalar doublet χ R with B − L charge equals -1, and a scalar bidoublet φ with zero B − L charge.The extra doublet and bidoublet are essential for breaking the symmetries of the model upon having nonvanishing VEVs.We define χ R = v R / √ 2 and assume that v R of an order TeV to break the righthanded electroweak sector together with B − L. Regarding the VEV of the scalar bidoublet φ, we use the parametrization where with v of an order O(100) GeV to break the electroweak symmetry of the SM.Here and after we use the definitions s x = sin x, c x = cos x, t x = tan x and ct x = cot x.
The scalar-fermion interactions in the LRIS model can be inferred from the left-right symmetric Yukawa Lagrangian which can be generally expressed as here i and j are family indices that run from 1 . . .3, y Q and ỹ Q represent the quark Yukawa couplings while y L and ỹL stand for the lepton Yukawa couplings.It should be noted that a Z 4 discrete symmetry is implemented in order to forbid a mixing mass term M Sc 1 S 2 and thus preserving the IS mechanism in a way similar to imposing Z 2 symmetry discussed in [54].Under this discrete Z 2 symmetry, S 1 has a charge −1 and all other fields have a charge +1.The fields in Eq. ( 8) are defined as The dual bidoublet φ and doublet χ R are defined as φ = τ 2 φ * τ 2 and χ R = iτ 2 χ * R respectively.Clearly, from the scalar sector of the model and before symmetry breaking, there are 12 scalar degrees of freedom: 4 of χ R and 8 of φ.After symmetry breaking, six scalars out of these 12 degrees of freedom remain as physical Higgs bosons while the other degrees of freedom are eaten by the neutral gauge bosons: Z μ and Z μ and the charged gauge bosons: W ± μ and W ± μ to acquire their masses.The physical Higgs bosons are two charged Higgs bosons, one pseudoscalar Higgs boson, and the remaining three are C P-even neutral Higgs bosons.For a detailed discussion about the Higgs sector of the LRIS model we refer to Ref. [55].
We turn now to the neutrino sector of LRIS model.After breaking the B − L symmetry, one finds that the neutrino Yukawa interaction terms lead to the following mass terms [54]: where M R and M D are 3 × 3 matrices related to v R and Dirac neutrino masses respectively via M R = y s v R / √ 2 and In addition one may generate very small Majorana masses for S 1,2 fermion through possible non-renormalizable terms.This tiny mass is required in the standard inverse seesaw mechanism for generating light neutrino masses [56][57][58][59].Thus, the Lagrangian of neutrino masses can be expressed as where μ s 10 −6 GeV, as it is suppressed by high nonrenormalizable scale.The inverse seesaw mechanism actually depends on small mass scale μ 2 , which violates the residual Lepton number symmetry after breaking the LR symmetry.The light neutrino masses vanish identically in the limit of μ 2 → 0, and the lepton number is restored.According to 't Hooft criteria, such a small scale is natural.We demonstrated that this mass parameter can be generated using a non-renormalizable dimension 7 operator [54], and despite its small size, it plays an important role in generating very small neutrino masses in this mechanism.In this context, the neutrino mass matrix can be written as ) and M ν is given by Note that in order to avoid a possible large mass term mS 1 S 2 in the Lagrangian (8), that would spoil the above inverse seesaw structure, one assumes that L R , χ R , and S 2 are even under a Z 2 -symmetry, while S 1 is an odd particle.Furthermore, similar to S 2 , the fermion S 1 can acquire a mass via another renormalizable term.As discussed in Ref. [60], S 1 can be a natural candidate for warm dark matter.It is a kind of sterile neutrino that has no mixing with active neutrinos and hence can only interact with Z gauge boson.As a consequence, S 1 is not subjected to all constraints imposed on sterile neutrinos due to their mixing with the active neutrinos.On the the hand, the mass of S 1 has to be of order O(10 keV) to satisfy the combined constraints from Lymanα forest data [61] and phase space arguments for fermionic dark matter [62].This can be accommodated by choosing μ s 1 = O(10 keV) without any conflict with the mass parameter μ s 2 contributing to the neutrino mass matrix given below which is required in our analysis.Moreover, for S 1 being hot dark matter, any related expected constraints will be set on its mass μ s 1 and its couplings to the Z gauge boson which have no effects on the tau couplings to neutrinos or neutrino masses or neutrino mixings needed in this study.
Following the standard procedure, the diagonalization of M ν results in the light and heavy neutrino mass eigen states ν i , ν h j with the mass eigenvalues given by: The light neutrino mass matrix in Eq. ( 13) must be diagonalized by the physical neutrino mixing matrix U MNS [63], i.e., Thus, one can easily show that the Dirac neutrino mass matrix can be defined as: where R is an arbitrary orthogonal matrix.Clearly, for M R at TeV scale and μ s 2 M R the light neutrino masses can be of order eV.The 9 × 9 neutrino mass matrix M ν , can be diagonalized with the help of the matrix V satisfying As a good approximation, the matrix V 3×3 can be given a where Generally F, as one can see from the expression of V 3×3 , is not unitary matrix.This unitarity violation, i.e., the deviation from the standard U MNS matrix, depends on the size of 1  2 F F T [65].The matrix V 3×6 is given by In addition, V 6×3 = (V 3×6 ) † .Finally, the matrix V 6×6 is the one that diagonalize the {ν R , S 2 } mass matrix.
The neutral scalar fields and their masses can be obtained if one expands the neutral components of the bidoublet φ and the doublet χ R around their vacua as follows where In this case, the symmetric mass matrix of the C P-odd Higgs bosons in the basis is given as (21) which can be diagonlized by a unitary matrix in the form with the physical mass of the pseudoscalar boson A, m 2  A , given by Similar to the C P-odd Higgs bosons, the elements of the (3×3) symmetric mass matrix of the C P-even Higgs bosons are given by where the potential parameter and λ 23 = 2λ 2 + λ 3 .This matrix can be diagonalized by a unitary transformation matrix For more details, we refer to Ref. [55].The process under study can receive dominant contributions from the charged Higgs mediation at tree level as we will show in the following.Thus, below, we list the relevant charged Higgs couplings related to our study following Ref.[55].In the flavor basis (φ ± 1 , φ ± 2 , χ ± R ), the charged Higgs bosons symmetric mass matrix takes the form The above matrix can be diagonalized by the unitary matrix, The mass eigenstates basis can be obtained using the rotation Here G ± 1 and G ± 2 represent the massless charged Goldstone bosons and H ± is a massive physical charged Higgs boson.The massless Goldstone bosons are eaten by the charged gauge bosons W μ and W μ to acquire their masses via the familiar Higgs mechanism.On the other hand, the mass of the charged Higgs boson H ± is given by: where α 32 = α 3 − α 2 is a potential parameter Ref. [55].
Clearly, the charged Higgs boson mass can be of the order of hundreds GeV if we pick out the values v R ∼ O(TeV) and α 32 ∼ O(10 −2 ).It is clear to see that, the physical charged Higgs boson is a linear combination of the flavor basis fields φ ± 1 , φ ± 2 , χ ± R , namely given as The effective Lagrangian, in the LRIS model, describing the charged Higgs and the charged gauge bosons W and W couplings to quarks and leptons can be expressed as The Lagrangians can be obtained from expanding L Y , given in Eq. ( 8), and rotating the fields to their corresponding ones in the mass eigenstates basis.It is direct to obtain where In the LRIS model and after electroweak symmetry breaking, quarks and charged leptons acquire their masses via Higgs mechanism.Consequently, we can express the quark Yukawa couplings in terms of the quark masses and CKM matrices in the left and right sectors as ) is the diagonal up (down) quark mass matrix, V L CKM and V R CKM are the CKM matrices in the left and right sectors respectively.The mixing matrices for left and right quarks result in the CKM matrices in the left and right sectors We can choose the bases where V u L = I and as a result, in this basis, V L CKM = V d L .Turning now to the right sector, we follow Ref. [66], where the matrix V R CKM can be written as a unitary matrix similar to V L C K M , but with new mixing angles θ R i j and a new Dirac phase δ R , as follows: here the diagonal matrices K and K contain five of six nonremovable phases in V R CKM .It was emphasized that CP violation and FCNC in the right-handed sector can be under control if the V R CKM is of the form where c θ R

13
= sin(θ R 13 ) and as we will see that, in the following, the phase α in the third column plays an important role in increasing the values of R D * .It is also possible to have a non-vanishing phase in the second column, which turns out to be irrelevant and has no effect on the R D or R D * results.Therefore, we set it to zero.We work in the bases where V d R = I and thus V R CKM = V u R , similar to the left-handed quark sector.
We proceed now to the lepton-neutrinio charged Higgs interactions in the model under study.These interactions contribute to the effective Lagrangian L ν H ± e f f and generally, can be written as with The lepton Yukawa couplings y L and ỹL can be expressed as where M is charged lepton diagonal mass matrix.Finally, the leptons and quarks gauge interactions related to our processes at tree-level can be deduced from the effective Lagrangian given by (44) where φ W is the mixing angle between W L and W R , which is of order 10 −3 .We have now all the ingredient required for deriving the effective Hamiltonian contributing to the transition b → cτ ν.In next section, we derive this Hamiltonian and list the expressions of the Wilson coefficients corresponding to this Hamiltonian.

The effective Hamiltonian relevant to the processes in the LRIS
In the presence of NP beyond SM, the effective Hamiltonian governs | c| = 1 B decays transition relevant to our processes up to one-loop level derived from the diagrams in Fig. 1, can be expressed as where ] and as mentioned in Ref. [43] the tensor operator with chirality ( cσ μν P R b)( ¯ j σ μν P L ν i ) vanishes.In case of = τ we find that the Wilson coefficients at the high energy scale μ = m H ± can be expressed as where i refers to the neutrino flavor, and we set the energy scale μ = 1 TeV.In our numerical analysis, where all results are evaluated at the bottom scale m b = 4.2 GeV, the evolution down of the Wilson coefficients from the scale μ = 1 TeV to the bottom scale can be inferred from the renormalization group evolution (RGE).It should be noted that in Eq. ( 46), the tensor contributions to the Wilson coefficients g L L T (μ) and g R R T (μ) are generated from the one-loop diagrams in the figure.Moreover, in the same equation, we kept only the dominant contributions to the scalar Wilson coefficients g AB S , with A, B run over L , R, which originate from treelevel diagrams in Fig. 1.As we will show below, only these tree-level scalar contributions will have sizable effects on the ratios R D,D * and other observables under concern in this work.
As can be seen from Eq. ( 46), the coefficients g L L S (μ) and g RL S (μ) have the same lepton vertex ) * .Therefore, they are expected to have different values due to receiving different contributions of the quark vertices H ± RL eff u i d j and H ± L R eff u i d j expressed in Eq. (36) in terms of y Q and ỹ Q .Furthermore, Eq. (37) shows that the contributions from the terms proportional to V L C K M to both y Q and ỹ Q are sup-pressed by the small down quark masses appear in the matrix M diag d .In the processes under consideration, the corresponding Wilson coefficient g L L S (μ) (g RL S (μ)) receives contributions only from the elements in the second column (row) of the matrix CKM , which are present in both of y Q and ỹ Q .With the texture of V R CKM given previously, we find that In this regard, one can show that R D , R D * and BR(B − c → τ − ντ ) receive contributions from the terms proportional to m t s θ R 13 e −iα which originate from g L L S (μ) only and not from g RL S (μ).The other contributions generated from the Wilson coefficient g RL S (μ) are proportional to the charm quark mass which is small compared to the top quark mass.
Finally in Eq. ( 46), the quantities g L L(R R) Z ,Z T refer to the suppressed contributions, comparing to photon ones, originating from Z and Z mediating the diagrams.Upon neglecting the small contributions g L L,R R Z,Z T we find the following relations For a charged Higgs of a mass 300 GeV we find that Clearly the tensor contributions to the processes under study can be safely neglected.This is the case also regarding the vector Wilson coefficients g AB V (μ), where AB can be any combination of the L and R chiralities, as they are suppressed by either In the given expressions below, the quantities g L L S , g RL S , g L R S and g R R S refer to the Wilson coefficients at the bottom scale μ = m b .In terms of these quantities, the ratios R M (M = D, D * ) are given as [16,67,68]: It should be noted that in the above expressions of R D,D * we assumed that NP effects are only present in the third generation of leptons (τ, ν τ ).This assumption is motivated by the absence of deviations from the SM for light lepton modes = e or μ.
The D * and τ longitudinal polarizations depend on the same Wilson coefficients affecting the R D,D * ratios.Thus, it is relevant to our investigation to show their predicted values for the set of the parameters of the model resolving the R D,D * anomalies.The two observables have been measured at Belle experiment.Their expressions can be written as [16,67,68] . In our analysis we use the measured values of the D * and τ longitudinal polarizations reported by Belle collaborations namely, F E x pt L (D * ) = 0.60 ± 0.08 ± 0.035 [69] and P E x pt τ (D * ) = −0.38 ± 0.51 +0.21 −0.16 [4,5,70].On the other hand their SM predictions are estimated as F SM L (D * ) = 0.464 ± 0.010, P SM τ (D) = 0.321 ± 0.003 and P SM τ (D * ) = −0.496± 0.015 [14].Although the tau polarization observable P τ (D) is known to be a good discriminator of scalar contributions originated in many beyond SM physics, (for instances the leptoquark (LQ) scenarios [71] ) there is no available experimental measurements of this observable so far.However, it is important to show its prediction in the model under study.
The tree-level charged Higgs boson exchange also modifies the branching ratio of the tauonic decay B − c → τ − ντ as follows [16,67,68] where where V cb stands for the CKM matrix element, τ B c and f B c denote the B − c meson lifetime and decay constant, respectively.The SM prediction of BR(B − c → τ − ντ ) SM = (2.25 ± 0.21) × 10 −2 [72].Unfortunately, no direct constraints from upper bounds on the leptonic B c branching ratios are available from the LHC.In view of this, an estimate of a bound on BR(B − c → τ − ντ ) has been derived from LEP data at the Z peak in Ref. [25].The bound turns to be strong BR(B − c → τ − ντ ) ≤ 10%.Later on, the bounds using the measured B c lifetime had been critically investigated and relaxed upper limit of ≤ 39% [26] and ≤ 60% [27,28] were obtained.Thus we will follow Ref. [72] and take in our analysis the bound: BR(B − c → τ − ντ ) < 60%.The processes B s → + A − A , where − A denotes a charged lepton, can be used to derive constraints on the parameter space on the model under concern in this work as we discuss in the following.These processes are mediated at tree-level by the neutral Higgs (H 0 k = H 0 , h 0 , A 0 ) exchange.Their branching ratios, including tree-level neutral Higgs contributions, can be written as where M Bs , τ B s and f B s stand for the B − s meson lifetime and decay constant, respectively.In the above equation, the neutral Higgs non-vanishing Wilson coefficients only include C bs P,S C bs P,S .Their expressions are listed in Eq. ( 69) in Appendix A.1.On the other hand and within SM, up to one loop-level, we have C bs P = C bs P = C bs A = 0 [73] and The function Y = η Y Y 0 is defined in a way that the NLO QCD effects are included in η Y = 1.0113 [74].The expression of the one loop Inami-Lim function Y 0 is given as [75] Y It should be remarked that, the SM Wilson coefficient C bs A is scale independent as its corresponding effective operator corresponds to conserved vector current with vanishing anomalous dimensions [73].
In our analysis we use the numerical values of the CKM matrix elements reported in Ref. [76].Moreover, the numerical values of the mass and life time of B s are taken from Ref. [77] and we take the value f B s = 0.230 GeV [78].Setting the neutral Higgs Wilson coefficients to zero in Eq. ( 56), we find that We note that the SM prediction for the process B 0 s → μ + μ − obtained here is larger than the well known one by Misiak  .This can be attributed to the fact that the authors of Ref. [79] performed extensive study and included O(α em ) and O(α 2 s ) corrections to the amplitude of the process.Experimentally, from the non observation of the decay B 0 s → τ + τ − , we have the upper limit B(B 0 s → τ + τ − ) < 6.8 × 10 −3 [77].This result allows new physics to have large contributions to B(B 0 s → τ + τ − ) and hence one obtains very loose constraints.This is not the case regarding the process B 0 s → μ + μ − for which the experimental measurements B B s → μ + μ − LHCb = (3.09+0.46+0.15−0.43−0.11 ) × 10 −9 [80].We refer to Ref. [10] for the results reported by the CMS, ATLAS and CDF collaborations.Using the 2σ range of the HFLAV average B B s → μ + μ − HFLAV = (3.45± 0.29) × 10 −9 [11] and with the help of Eq. ( 56), we can derive the required constraints on our parameter space In the model under consideration in this work, the B s(d) − Bs(d) neutral meson mixing receives new contributions from tree-level diagrams mediated by the exchange of neutral Higgs bosons, box diagrams mediated by the charged Higgs only, the W bosons only, both charged Higgs and W or W ± together and finally both W and W ± together.In Appendix A.1.1,we list the set of the operators contributing to the effective Hamiltonian governing the B s(d) − Bs(d) mixing and their corresponding Wilson coefficients.Following the analysis in Ref. [81], we can use the reported experimental and SM values of M B s,d and M SM B s,d respectively to impose the constraints 0.85 03 up to 2σ level.The b → sγ can lead to constraints on the parameter space of the model under consideration.In our investigation of these constraints, we work in leading logarithmic (LL) precision.Charged Higgs can mediate a loop diagram similar to the one in the SM but with replacing the charged W bosons with the charged Higgs.On the other hand, the contributions of the neutral Higgs boson to b → sγ are suppressed and thus can be neglected.This can be explained as the flavor off-diagonal elements in the down sector can be stringently constrained from the tree-level decays.Thus, we are left with the contributions originating from charged Higgs mediating the loop diagram.
In the small t β and α 32 scheme adopted in this study, we find that the effective couplings H ± L R eff cb and H ± L R eff cs appearing in Eq. ( 36) are so tiny, regardless the values of the parameters α and θ 13 .This is not the case if one considers and H ± RL eff cb which can be of order O(10 −2 ) or larger than that.Consequently, in our analysis, we keep only the dominant contributions to b → sγ in the following.Upon, neglecting the operators with mass dimension higher than six, one obtains the same effective Hamiltonian as in the case of the SM [82] In the approximation we adopted above, charged Higgs, propagating in the loop, contributes only to the Wilson coefficients C H ± , where y j = m 2 u j /m  −8y 3 j + 3y 2 j + 12y j − 7 + (18y 2 j − 12y j ) ln y j (y j − 1) 4  , From [83][84][85] we have leading to [83] Here, we used C H ± 7,8 at a matching scale of 1 TeV as input.Again, these constraints are so stringent that the effect of C H ± 7,8 on the flavour anomalies can be mostly neglected.Finally, It is worth to recall that, due to the requirement of having light neutrino masses we found that the contributions of the right neutrino sector to the corresponding Wilson coefficients g R R S (μ) and g L R S (μ) are very small and thus can be safely ignored, leaving us with only g L L S (μ) and g RL S (μ).

Numerical results and analysis
Having discussed all relevant constraints related to the ratios R D,D * , we are ready now to estimate their predictions in the model under concern in this investigation.To proceed, we need first to show, numerically, the RGE running effect of g RL,L L S (μ) from μ = 1 TeV to m b = 4.2 GeV scale.According to the estimation carried out in Refs.[86,87], we have 1.71 0 0 1.71 3 R D and R D * as function of the charged Higgs mass where the parameters as stated in the text are chosen as follows: Furthermore, we have checked that g RL S is about one order of magnitude smaller than g L L S ; thus for real g L L S (i.e., Re(g L L S ) = g L L S ), one finds This expression clearly shows that enhancing the values of R D * to be in the range of the given experimental results, while keeping the limit −1 1 in mind is possible for a range of negative values of g L L S as can be seen from the left plot in Fig. 2.However, as can be remarked from the right plot in Fig. 2, these negative values reduce R D below their allowed 2σ region of the experimental results shown by the shaded green region in the plot.Clearly, we deduce that the phase α of the mixing matrix V R C K M is crucial to solve the concerned anomalies for the processes under consideration.
The attainable values of the Wilson coefficients g L L S and g RL S are affected by the charged Higgs mass.As a result, to enhance these coefficients and thus the values of R D and R D * while adhering to direct search constraints, the charged Higgs masses should not be too heavy, namely of the order of hundreds GeV.According to Eq. ( 32), this can be accomplished by considering v R ∼ O(TeV), α 32 ∼ O(10 −2 ) and t β less than one.Remarkably, from the pre factor 1/(t β −ct β ) 4 R D and R D * as function of the sin of mixing angle θ R 13 left and the phase α of the V R C K M matrix right and other parameters are fixed as in the previous figure in Eq. ( 37), it is direct to see that small t β values can enhance also the quark Yukawa couplings and hence together with the angle θ R 13 and the complex phase α of the right-quark mixing V R C K M , defined in Eq. ( 39), play crucial roles in increasing the values of R D and R D * , and allow them to take values that are compatible with the limit of the experiments at the same time.In our scan of the parameter space we take α 32 ∈ [0.007, 0.016], With this in hand, we show below our results corresponding to the scanned points in the parameter space respecting all the bounds discussed in the previous section.It should be noted that, in these results we only select the points that lead to values of F L (D * ) and P τ (D * ) within their 2σ range of the corresponding experimental results.Moreover, we have checked that, for these parameters in the chosen ranges and values, the contribution of the g RL S are irrelevant and can be safely neglected.This confirms our previous conclusion that only the Wilson coefficient g L L S plays the major role through the terms proportional to m t s θ R 13 e −iα .In the following, we present a set of elucidative plots of R D and R D * versus some selected relevant parameters of the model.It should be noted that these plots were generated through a random scan and therefore do not reflect any specific correlations with the chosen parameter.
In Fig. 3, we display the variation of R D and R D * with the charged Higgs mass.As can be seen from the figure, it is possible to account for the experimental results of R D and R D * within 2σ range, while respecting the the mentioned constraints in the previous section and the 2σ range of F L (D * ) and P τ (D * ) mentioned in the previous section, with charged Higgs masses m H ± can be chosen of order 300 GeV.On the other hand, the dependence of R D and R D * on the sin of the mixing angle θ R 13 and the phase α of the V R C K M matrix is depicted in Fig. 4. It is clear from left plot in the figure that, moderate and large values of s θ R 13 are preferable to satisfy the experimental results of R D and R D * within 2σ range while respecting the bounds and the requirements considered in the Fig. 5 Allowed region in the (θ R 13 , α) plane by all considered constraints in green color.The red region satisfies the 2σ experimental results of R D * and R D together for t β 0.2074, α 32 = 0.00416 which result in m H ± = 305 GeV and the other parameters are fixed as before scan.On the other hand, regarding the phase α, large phases are favored as can be seen from the right plot in the figure.The explicit dependence of R D and R D * on the parameters α 23 , which determines the mass of the charged Higgs M H ± , θ 13 , and the phase α, can be understood through the following approximate analytical expressions: Table 1 A set of selective benchmarks of the parameters and intended predictions resulting from the scan after taking into account the constraints discussed in the previous section and requirements stated in the beginning of this section We can obtain the regions in the (θ R 13 , α) parameter space in which the anomalies are satisfied through varying θ R 13 and α while assigning fixed values of the other parameters.As an example, we take the fixed values t β 0.2074 and α 32 = 0.00416 which result in m H ± = 305 GeV and the other parameters are fixed as before.In Fig. 5, we show the experimentally allowed 2σ region in the (θ R 13 , α) plane satisfying R D * and R D together in red color.In the same plot, the regions in green color are the allowed regions, in the (θ R 13 , α) plane for the set of the input parameters we use, by all constraints discussed above in the previous section.Clearly, the imposed constraints have a sensible effect on the parameter space as large part of this space is excluded by the constraints.Moreover, as shown in the plot, their is a small region in the parameter space which is allowed by all the aforementioned constraints in which R D * and R D anomalies are satisfied together.Specifically, this region is the intersection region of the two colored regions in the plot.It should be noted that this conclusion corresponds to our particular example of the input parameters stated above.Taking other values of the input parameters t β , α 32 and v R may lead to another regions in (θ R 13 , α) plane that respect all constraints and satisfy R D * and R D anomalies together.For instances, in Table 1, we list several benchmarks obtained upon variation of α 32 while keeping other input parameters as before.
In Fig. 6, we present the correlation between R D and R D * for the same set of the parameter space considered in the scan over the values and ranges mentioned in the beginning of this section.Only points that satisfy all the bounds discussed in Sect. 3 and the 2σ range of F L (D * ) and P τ (D * ) are included here.It is remarkable from this figure that both R D and R D * are satisfied for a lot of points in the parameter space, thanks to the complex mixing of right-handed quarks.
Fig. 6 The correlation between R D and R D * for the same set of parameter space considered in Fig. 3.The dashed red line represents the χ 2 = 1.0 contour, which is consistent with the 2023 latest results from HFLAV [11] It should be noted that, if the tbH + Yukawa coupling is non-negligible (c R θ13 is non-zero), there will be nonnegligible ttφ couplings thanks to the SU(2) invariance where φ stands for the neutral extra scalars: H and A. As a consequence, gg → φ → τ τ would also impose constraints on the parameter space.To consider this possibility, we have analyzed the processes gg → φ → τ τ .We have found that the corresponding cross-sections are less than 10 −8 Pb, indicating that they do not impose any significant constraints on our analysis.
We turn now to the discussion of our predictions of the D * and τ longitudinal polarizations and tau polarization observable P τ (D).For this purpose we list in Table 1 a selective set of 6 benchmarks representing part of our results.Recall that the corresponding measurements, reported by Belle collaborations, are F E x pt L (D * ) = 0.60 ± 0.08 ± 0.035 [69] and P E x pt τ (D * ) = −0.38 ± 0.51 +0.21  −0.16 [4,5,70].Clearly, the uncertainty in the measurement of P E x pt τ (D * ) is so large and thus we only restrict ourselves to F L (D * ).The SM prediction is estimated to be F SM L (D * ) = 0.464 ± 0.010 [14].Clearly from Table 1, the new contributions of the model under concern can increase the central value of the SM by about 11% for the listed set of the benchmarks.We also noted there are points in our scan that imply a larger enhancement which makes our predictions closer to the experimentally measured central value.On the other hand, the enhancement in the tau polarization P τ (D) can reach 6 times its SM predicted value as can be noted from the same table.This can serve as test of the model under consideration once this observable is measured.
Searching for new physics via B meson exclusive decays originating at the quark level from the b → s + − transitions has gained a lot of attention in the last decades.In the framework of the SM, these transitions take place only at loop-level.Consequently, their amplitudes are suppressed, for instances, by factors accounting for the integration over the momenta running in the loop.This in turn leads to deviations of the measured branching ratios of B → K μ + μ − , B → K * μ + μ − , and B s → φμ + μ − decays.from their SM predictions.Also, the angular observable P 5 in the B 0 → K * 0 μ + μ − decay, [88,89], has shown tension with the SM values.For instance, ATLAS [90], and LHCb [91,92], measured the value of P 5 in the kinematical region 4.0 < q 2 < 6.0 GeV 2 and found departure from the SM value to be more than 3σ [93].Furthermore Belle [94,95] and CMS [96] measured the value of P 5 for the same decay mode in q 2 bin 4.0 < q 2 < 8.0 GeV 2 and 6.0 < q 2 < 8.68 GeV 2 respectively.Belle measurement shows the deviation of 2.6σ from the SM prediction and CMS measurement shows a discrimination of 1σ from the SM value.Their is a possibility that the tensions between the measurements and the SM predictions of all these observables can be relaxed if one considers NP models affecting the b → s + − transitions.In the literature, various performed model independent global fit analyses [97][98][99][100][101][102][103][104][105][106][107][108][109][110][111][112][113][114] based on the assumption that NP present only in the muon sector revealed that two simple one-dimensional (1D) NP scenarios (S1) C NP 9μ or (S2) C NP 9μ = −C NP 10μ , that give better fit to all the data, with preferences reaching ≈ 5 − 6σ compared to the SM [115].
The observable P 5 is sensitive to the real parts of C NP 9μ and C NP 10μ .On the other hand, the ratios R K ( * ) , defined as

Conclusion
In this work we have explored the possibility of resolving the tension between the SM prediction and the experimental results of the R D and R D * ratios using a low scale left-right symmetric model based on SU The scalar sector of the model contains charged Higgs boson with masses that can be chosen in the order of hundreds GeV without any conflict with direct search constraints.We have shown that integrating out the charged Higgs mediating the tree-level diagrams generates a set of non vanishing scalar Wilson coefficients contributing to the effective Hamiltonian governing the transition b → cτ ν and hence to the ratios R D * and R D and the D * and τ polarizations.
The dependency of the scalar Wilson coefficients on the matrix elements of the quark mixing angle in the right sector turns to be important.We emphasized that the mixing element (V R C K M ) 23 should be complex in order to satisfy both R D and R D * .We have also shown the complex phase associated with this mixing element is essential to accommodate the experimental results of the ratios for charged Higgs masses of order 300 GeV while respecting the constraints from BR(B − c → τ − ντ ), B s(d) − Bs(d) mixing and other relevant constraints discussed above.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited.[Authors' comment: Data sharing is not applicable to this article as no datasets were generated or analysed during the current study].
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.

Appendix
A.1 Wilaon coefficients relevant to the process B s → + A − A Following Ref. [73], the neutral Higgs non-vanishing Wilson coefficients can be expressed as The expressions of Z A,H i j can be found in Ref. [55].The Wilson coefficients at low energy scale μ low , C ( )bs S,P (μ low ), can be obtaining from their corresponding ones listed in Eq. (69) using the relation [73] C ( )bs here m q represents the running quark mass with the appropriate number of active flavors.Finally, it should be noted that since the Wilson coefficients are given at the matching scale, m b and m t must be evaluated at this scale [73].

A.1.1 B q − Bq mixing
The contributions of the NP to the effective hamiltonian generating B = 2 transitions, q = d, s, can be written as The four-quark operators are given as The operators Q1,2,3 can be obtained from the operators Q 1,2,3 by the replacement L ↔ R. At tree-level, we have contributions to H Q=2 eff from only neutral Higgs bosons mediation.Denoting these contributions by C H 0 k i , we find that their expressions can be written as The quantities AB H 0 k qb are defined as before.We remark from the above expressions that, the neutral Higgs bosons contributes only to the scalar color singlet four-quark operators in Eq. ( 80).In Ref. [81], it was pointed out that the contributions of the operators that contains scalar and tensor Dirac structures are highly disfavored by the fits to b → s data and hence their contributions can be neglected.The charged Higgs contributions to H Q=2 eff originate at one loop-level and can be expressed as where the loop functions G i for i = 1, 2, 3 are given as  [118] (see also [119]) where η = α s (μ N P ))/α s (m t ), a j , b (r,i) j and c (r,i) j are magic numbers given in [118] and B i B s are the B parameters that can be found in Table 9 in Ref. [120].It should be noted that the magic numbers for the evolution of the Wilson coefficients C1−3 are the same as the ones for the evolution of C 1−3 [119].Moreover, in the basis Qi , the B B q i parameters and the hadronic matrix elements B q | Qr | Bq are equal to their corresponding ones in the basis Q i [121].The matrix elements are given by with N r = (−5/24, 1/24, 1/4, 1/12) for r = (2,3,4,5).With all this in hand, it is direct to calculate the quantity M q [122,123] and thus derive the bounds on the parameter space using the measured value of M q .
The K − K mixing can be studied in a similar way used above in the B q − Bq mixing.In fact, it is possible to have a non-vanishing phase in the second column of the V R CKM matrix, which turns out to be irrelevant and has no effect on the R D and R * D result.However, this imaginary phase can affect the bound from other flavor observables e.g.K .In particular, through the contributions from the neutral Higgs mediating the tree-level diagrams.In order to provide the ingredient required for the estimation of K , we define where H S=2 are "magic numbers" listed in [124] and B i K are the B parameters collected in Table 8.The matrix elements are given by with N r = (−5/24, 1/24, 1/4, 1/12) for r = (2, 3, 4, 5).In the next step, the Wilson coefficients are evolved down from the mass scale of the Higgs scalars to the scale μ = 2 GeV at which the hadronic matrix elements are evaluated using the RG equations in [124].Doing so, it is quite forward and direct to compute the constraint on the non-vanishing phase in the second column of the V R CKM matrix using the bound from K .

Fig. 1
Fig. 1 Diagrams contributing to H ef f , in Eq. (45) up to one loop-level due to charged Higgs mediation with V can be Z or γ or both of them depending on the fermion lines connecting them.In the figure when A and B represent b and c quarks the other two external lines will repre-

2 W 6 . 4 ×
10 −3 for m W O(1 TeV) or by sin φ W φ W O(10 −3 ) or both of them.Consequently, we are left only with contributions of the scalar Wilson coefficients.

7 and C H ± 8 corresponding to the operators O 7
= e 16π 2 m b sσ μν P R bF μν ; O 8 = g s 16π 2 m b sσ μν T a P R bG a μν

Fig. 2
Fig. 2 Left (right) R D * (R D ) variation with g L L S where the shaded green regions represent the allowed 2σ region of their experimental values

G 2 (= m 2 W /m 2 H
x, y, z) = − 1 (x − y)(1 − z) × x log x 1 − x − y log y 1 − y − x log x z z − x + y log y z z − y , + and x k = m 2 u k /m 2 H + .The next step is to calculate the matrix elements of the operators Q i at the scale μ = μ b and to run the Wilson coefficients from the electroweak scale to the scale μ = μ b .The contribution to the B q − Bq mixing amplitudes induced by a given NP scale coefficient C i (μ = μ N P ), denoted by B q |H B=2 eff | Bq i , as a function of α s (μ N P ) and the scale μ = μ b is given as NP .The calculation of K 0 |H S=2 full | K 0 i can be done in a similar manner to B q |H B=2 eff | Bq i with the replacement of bs ↔ sd in the Wilson coefficients listed above.The K − K mixing can be evaluated using[118]

Table 1 .
B→K ( * ) e + e − ) , are also sensitive to C NP 9μ and C NP 10μ .Recently the updated measurements of R K ( * ) ,[116,117]have put stringent constraints on the NP couplings and the NP models.This in turn lead to a strong constraints on the Wilson coefficients C NP 9μ and C NP 10μ which affect the aforementioned observables.In the model under concern we found that it is not possible to accommodate P 5 and other observables related to the b → s + − transitions through Re(C N P 9μ ) and Re(C N P 10μ ).To make it more clear, we show in Table2As can be seen from the table there is no sizable enhancement due to the new contributions of the model under consideration compared to the corresponding SM Wilson coefficients.