Phenomenology of using lattice QCD calculations
In a recent paper we studied the effect of new-physics operators with different Lorentz structures on the semileptonic decay. This decay is of interest in light of the puzzle in the semileptonic decays. In this work we add tensor operators to extend our previous results and consider both model-independent new physics (NP) and specific classes of models proposed to address the puzzle. We show that a measurement of can strongly constrain the NP parameters of models discussed for the puzzle. We use form factors from lattice QCD to calculate all observables. The tensor form factors had not previously been determined in lattice QCD, and we present new lattice results for these form factors here.
A major part of particle physics research is focused on searching for physics beyond the standard model (SM). In the flavor sector a key property of the SM gauge interactions is that they are lepton flavor universal. Evidence for violation of this property would be a clear sign of new physics (NP) beyond the SM. In the search for NP, the second and third generation quarks and leptons are quite special because they are comparatively heavier and are expected to be relatively more sensitive to NP. As an example, in certain versions of the two Higgs doublet models (2HDM) the couplings of the new Higgs bosons are proportional to the masses and so NP effects are more pronounced for the heavier generations. Moreover, the constraints on new physics, especially involving the third generation leptons and quarks, are somewhat weaker allowing for larger new physics effects.
The charged-current decays have been measured by the BaBar Lees:2013uzd (), Belle Huschle:2015rga (); Abdesselam:2016cgx () and LHCb Aaij:2015yra () Collaborations. It is found that the values of the ratios , where , considerably exceed their SM predictions.
This ratio of branching fractions has certain advantages over the absolute branching fraction measurement of decays, as this is relatively less sensitive to form factor variations and several systematic uncertainties, such as those on the experimental efficiency, as well as the dependence on the value of , cancel in the ratio.
These values are also in good agreement with the phenomenological prediction Sakaki:2013bfa ()
which is based on form factors extracted from experimental data for the differential decay rates using heavy-quark effective theory. See also Ref. Wang:2017jow () for a recent analysis of form factors using light-cone sum rules.
A calculation of is not yet available from lattice QCD. The phenomenological prediction using form factors extracted from experimental data is Fajfer:2012vx ()
The averages of and measurements, evaluated by the Heavy-Flavor Averaging Group, are HFAGWinter2016 ()
where the first uncertainty is statistical and the second is systematic. and exceed the SM predictions by 3.3 and 1.9, respectively. The combined analysis of and , taking into account measurement correlations, finds that the deviation is 4 from the SM prediction HFAGWinter2016 (); Ricciardi:2016pmh ().
Since lattice QCD results are not yet available for the form factors at nonzero recoil and for the tensor form factor, we use the phenomenological form factors from Ref. Sakaki:2013bfa () for both channels in our analysis. For , we have compared the phenomenological results for and to the results obtained from a joint BGL -expansion fit DeTar () to the FNAL/MILC lattice QCD results Bailey:2015rga () and Babar Aubert:2009ac () and Belle experimental data Glattauer:2015teq (), and we found that the differences between both sets of form factors are below 5% across the entire kinematic range. The constraints on the new-physics couplings from the experimental measurement of obtained with both sets of form factors are practically identical.
There have been numerous analyses examining NP explanations of the measurements Fajfer:2012jt (); Crivellin:2012ye (); Datta:2012qk (); Becirevic:2012jf (); Deshpande:2012rr (); Celis:2012dk (); Choudhury:2012hn (); Tanaka:2012nw (); Ko:2012sv (); Fan:2013qz (); Biancofiore:2013ki (); Celis:2013jha (); Duraisamy:2013kcw (); Dorsner:2013tla (); Sakaki:2013bfa (); Sakaki:2014sea (); Bhattacharya:2014wla (). The new physics involves new charged-current interactions. In the neutral-current sector, data from decays also hint at lepton flavor non-universality – the so called puzzle: the LHCb Collaboration has found a 2.6 deviation from the SM prediction for the ratio in the dilepton invariant mass-squared range 1 GeV GeV Aaij:2014ora (). There are also other, not necessarily lepton-flavor non-universal anomalies in decays, most significantly in the angular observable Aaij:2013qta (); Aaij:2015oid (). Global fits of the experimental data prefer a negative shift in one of the Wilson coefficients, Blake:2016olu (). Common explanations of the and anomalies have been proposed in Refs. Bhattacharya:2014wla (); Calibbi:2015kma (); Greljo:2015mma (); Boucenna:2016wpr (); Bhattacharya:2016mcc (); Barbieri:2016las ().
The underlying quark level transition in the puzzle can be probed in both and decays. Recently, the decay was discussed in the standard model and with new physics in Ref. Gutsche:2015mxa (); Woloshyn:2014hka (); Shivashankara:2015cta (); Dutta:2015ueb (); Faustov:2016pal (); Li:2016pdv (); Celis:2016azn (). decays could be useful to confirm possible new physics in the puzzle and to point to the correct model of new physics.
In Ref. Shivashankara:2015cta () the following quantities were calculated within the SM and with various new physics operators:
where represents or . In this paper we work with the ratio , defined as
We also consider the forward-backward asymmetry
where is the angle between the momenta of the lepton and baryon in the dilepton rest frame.
This paper improves upon the earlier work Shivashankara:2015cta () in several ways:
We add tensor interactions in the effective Lagrangian.
In addition to and , we also calculate the forward-backward asymmetry (12) in the SM and with new physics.
In addition to analyzing the effects of individual new physics-couplings, we study specific models that introduce multiple new-physics couplings simultaneously. We consider a 2-Higgs doublet model, models with new vector bosons, and several leptoquark models.
The paper is organized in the following manner: In Sec. 2 we introduce the effective Lagrangian to parametrize the NP operators and give the expressions for the decay distribution in terms of helicity amplitudes. In Sec. 3, we present the new lattice QCD results for the tensor form factors. The model-independent phenomenological analysis of individual new-physics couplings is discussed in Sec. 4, while explicit models are considered in Sec. 5. We conclude in Sec. 6.
2.1 Effective Hamiltonian
where is the Fermi constant, is the Cabibbo-Kobayashi-Maskawa (CKM) matrix element, and we use . We consider that the above Hamiltonian is written at the energy scale.
If the effective interaction is written at the cut-ff scale then running down to the scale will generate new operators and new contributions, which have been discussed in Refs. Feruglio:2016gvd (); Feruglio:2017rjo (). These new contributions can strongly constrain models but to really calculate their true impacts we have to consider specific models where there might be cancellations between various terms.
The SM effective Hamiltonian corresponds to . In Eq. (LABEL:eq:Heff), we have assumed the neutrinos to be always left chiral. In general, with NP the neutrino associated with the lepton does not have to carry the same flavor. In the model-independent analysis of individual couplings (Sec. 4) we will not consider this possibility. Specific models will be discussed in Sec. 5.
2.2 Decay process
The process under consideration is
The differential decay rate for this process can be represented as Tanaka:2012nw ()
and the helicity amplitude is written as
Here, (, ) indicate the helicity of the virtual vector boson (see Appendix A), and are the helicities of the baryon and lepton, respectively, and for and for .
The scalar-type, vector/axial-vector-type, and tensor-type hadronic helicity amplitudes are defined as
The leptonic amplitudes are defined as
Above, are the polarization vectors of the virtual vector boson (see Appendix A). The explicit expressions for the hadronic and leptonic helicity amplitudes are presented in the following.
Hadronic helicity amplitudes
In this paper, we use the helicity-based definition of the form factors, which was introduced in Feldmann:2011xf (). The matrix elements of the vector and axial vector currents can be written in terms of six helicity form factors , , , , , and as follows:
The matrix elements of the scalar and pseudoscalar currents can be obtained from the vector and axial vector matrix elements using the equations of motion:
In our numerical analysis, we use GeV, GeV Olive:2016xmw (). The matrix elements of the tensor currents can be written in terms of four form factors , , , ,
The matrix elements of the current can be obtained from the above equation by using the identity
In the following, only the non-vanishing helicity amplitudes are given. The scalar and pseudo-scalar helicity amplitudes associated with the new physics scalar and pseudo-scalar interactions are
The parity-related amplitudes are
For the vector and axial-vector helicity amplitudes, we find
We also have the relations
The tensor helicity amplitudes are
The other non-vanishing helicity amplitudes of tensor type are related to the above by
Leptonic helicity amplitudes
In the following, we define
The scalar and pseudoscalar leptonic helicity amplitudes are
the vector and axial-vector amplitudes are
and the tensor amplitudes are
Here we have the relation
The angle is defined as the angle between the momenta of the lepton and baryon in the dilepton rest frame.
2.3 Differential decay rate and forward-backward asymmetry
From the twofold decay distribution (14), we obtain the following expression for the differential decay rate by integrating over :