Lepton-pair production in di-pion τ lepton decays
LLepton-pair production in di-pion τ lepton decays J. L. Gutiérrez Santiago ∗ a , G. López Castro † a , P. Roig ‡ a a Departamento de Física, Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional,Apdo. Postal 14-740, 07000 Ciudad de México, México.
Abstract:
We study the τ − → ν τ π − π (cid:96) + (cid:96) − ( (cid:96) = e, µ ) decays, which are O ( α ) -suppressed with respect to the dominant di-pion tau decay channel. Both the inner-bremsstrahlung and the structure- (and model-)dependent contributions are consid-ered. In the (cid:96) = e case, structure-dependent effects are O (1%) in the decay rate,yielding a clean prediction of its branching ratio, . × − , measurable with BaBaror Belle(-II) data. For (cid:96) = µ , both contributions have similar magnitude and we geta branching fraction of (1 . ± . × − , reachable by the end of Belle-II operation.These decays allow to study the dynamics of strong interactions with simultaneousweak and electromagnetic probes; their knowledge will contribute to reducing back-grounds in lepton flavor/number violating searches.
1. Introduction
Among the different charged lepton flavours, the τ lepton has the largest mass, enough to produce a varietyof hadronic states which provide an ideal environment to study the dynamics of the strong interactions as far asphase space allows. In this sense, τ lepton decays involve energy domains where the resonance degrees of freedombecome relevant. By means of semileptonic τ decays, we can study the hadronization of the weak chargedcurrents and use the resulting hadronic vertices either to test the fundamental parameters of the StandardModel (SM) or to understand the properties of Quantum Chromodynamics and the electroweak sectors in aclean way [1, 2].This article studies the semileptonic five-body decay τ − → ν τ π − π (cid:96) + (cid:96) − with the lepton pair ( (cid:96) = e or µ )produced via a virtual photon. The corresponding radiative case τ − → ν τ π − π γ has been analyzed using theResonance Chiral Theory (RChT) [3, 4, 5] and vector meson dominance [6, 7] approaches to describe vector andaxial-vector form factors involved in the W − → π − π γ vertex (see also ref. [8], which includes isospin-breakingand electromagnetic radiative corrections according to the Hidden Local Symmetry model, last updated in ref.[9]). Interestingly, this weak vertex involves also the interplay with strong and electromagnetic interactions.From a more practical point of view, the study of this decay is useful because it may pollute searches forprocesses involving lepton flavor violation in the charged sector or lepton number violation. In addition, itcould serve to verify the radiative corrections used in the contribution of the hadronic vacuum polarization(HVP) entering the anomalous magnetic moment of the muon ( a µ ) obtained using hadronic τ decays data[10, 11, 12, 7, 13, 14, 8, 15, 16, 5].Specifically, the decay τ − → ν τ π − π (cid:96) + (cid:96) − , with π − misidentified as muon and undetected π , may be animportant background that can mimic the signal in searches for lepton flavour violation (LFV) processes ofthe form τ − → (cid:96) ( (cid:48) ) − (cid:96) + (cid:96) − . Currently the branching ratios of these LFV decay modes have upper bounds of O (10 − ) [17], while their SM predictions are unmeasurably small [18, 19]. These decays (for (cid:96) = µ ) can alsobe misidentified as lepton number violating processes of the type τ − → ν τ µ − µ − π + [20]. To avoid these decayspolluting new physics searches, it will be most useful to include them in the Monte Carlo Generator TAUOLA[21], where other tau decay modes including (cid:96) + (cid:96) − pairs [22, 23] have recently been incorporated [24].Ref. [10] first took advantage of the clean LEP tau data samples to evaluate a HV P,LOµ using tau data.At the level of precision attained in the last twenty years (see ref. [25] and references therein), one requirescomputing the (model-dependent) electromagnetic and isospin breaking corrections relating σ ( e + e − → hadrons) in the isovector channel to the hadronic tau decay distributions, which were taken into account in subsequentevaluations [11, 12, 7, 13, 14, 8, 15, 16, 5]. Particularly, refs. [12, 7, 5] highlighted that different observablesof τ − → ν τ π − π γ decays for a not-so-low cut on photon energies (so that the inner bremsstrahlung (IB) part ∗ jlgutierrez@fis.cinvestav.mx † glopez@fis.cinvestav.mx ‡ proig@fis.cinvestav.mx a r X i v : . [ h e p - ph ] D ec oes not saturate the observables) can reduce substantially the model-dependent error in τ -based evaluations of a HV P,LOµ . We hope that the analysis presented in this work of the τ − → ν τ π − π (cid:96) + (cid:96) − decays can be helpful inreducing the error of a HV P,LOµ using tau data. We will find that the (cid:96) = µ case is promising in this respect, as themodel-dependent contributions are of the same size of the IB part. Conversely, its low branching ratio, ∼ − will challenge the Belle-II analysis [26]. For the (cid:96) = e case the situation will be the opposite, with a ∼ − branching ratio (measurable already with BaBar or Belle data), that has little ( O (%) ) model-dependence.The IB part is model-independent [27], while the structure-dependent part is not. The vector and axial-vectorform factors entering the latter can be computed using the RChT framework. Including operators contributing-upon integrating resonances out- to the O ( p ) Chiral Perturbation Theory (ChPT) [28] low-energy couplings–as in ref. [12]– all free parameters are related to the pion decay constant (after applying short-distance QCDconstraints), which results in controlled ( O (20%) ) errors for our prediction in the (cid:96) = µ mode. In the (cid:96) = e case, structure-dependent effects are as small as (uncomputed) one-loop QED corrections, which set the size ofour corresponding uncertainty.The paper is organized as follows: we start with a very short review of the radiative process τ − → ν τ π − π γ in order to introduce our conventions and recall the main features of the process under study. In section 3 wedescribe the amplitude for lepton pair production τ − → ν τ π − π (cid:96) + (cid:96) − . Section 4 deals with both, vector andaxial-vector, structure-dependent amplitudes. We derive the corresponding basis for the relevant (vector) caseat this order in the chiral expansion in section 4.1, relegating the axial-vector structure (which only appears atthe next order) to section 4.2. The branching ratio and the (cid:96) + (cid:96) − invariant mass spectrum for both channelsare predicted in section 5. Finally, we provide our conclusions in section 6.
2. The radiative τ − → ν τ π − π γ decay The τ − → ν τ π − π decay is the dominant channel among tau decays. The precise measurement of the di-pionmass spectrum allows to extract the weak pion form factor, which can be compared to the electromagnetic pionform factor measured in electron-positron collisions via the Conserved Vector Current (CVC) hypothesis. Thissame spectrum is useful to extract, in a clean way, information on the tower of vector resonance parametersthat are produced in the hadronization of the isovector current. On the other hand the radiative di-pion taulepton decay, τ − → ν τ π − π γ , provides additional information on the hadronization of the weak current and itis necessary to account for the di-pion observables at the few percent level. Finally, the lepton-pair productioninduced by the virtual photon in the radiative decay, namely τ − → ν τ π − π γ ∗ ( → (cid:96) + (cid:96) − ) , serves to furtherscrutinize the hadronization of the weak current in an extended kinematical domain. As previously pointed out,it allows also to quantify one of the sources of background in searches of three charged lepton flavor violatingdecays of tau leptons.In order to fix our conventions in the study of lepton-pair production, let us first consider the radiativedi-pion decays of the tau lepton. The matrix element for the process τ − ( P ) → ν τ ( q ) π − ( p − ) π ( p ) γ ( k ) , has thegeneral structure [12]: M ( γ ) = eG F V ∗ ud (cid:15) ∗ µ ( k ) (cid:34) ( p − p − ) ν f + ( t )2 P · k u ( q ) γ ν (1 − γ ) (cid:0) /P − /k + M τ (cid:1) γ µ u τ ( P ) . + ( V µν − A µν ) u ( q ) γ ν (1 − γ ) u τ ( P ) (cid:35) . (1)Here, G F is the Fermi coupling constant, V ud the Cabibbo-Kobayashi-Maskawa quark mixing matrix element, e is the magnitude of the electron charge and (cid:15) ∗ µ ( k ) the photon polarization four-vector. The first term in Eq.(1) corresponds to the photon emission off the tau lepton. It is described in terms of the charged pion vectorform factor, defined as (cid:104) π − π | dγ µ u | (cid:105) = √ f + ( t ) ( p − − p ) µ , with t = ( p − + p ) the invariant-mass of the di-pion system. Later, in our numerical analysis we will use f + ( x ) as derived in Refs. [29, 30].The second term in Eq. (1) contains the (structure-dependent) vector V µν and axial-vector A µν components.They describe the hadronization of the weak current involving an additional photon: W − ( P − q ) → π − ( p − ) π ( p ) γ ( k ) . (2)2ccording to Low’s soft-photon theorem [27], the leading terms (LO) of the radiative amplitude in the photonenergy expansion are fixed in terms of the non-radiative amplitude and the gauge invariance requirement of thetotal amplitude (equivalently, k µ V µν = f + ( t )( p − − p ) ν ). This leads to ( t (cid:48) = ( p − + p + k ) = t + 2( p − + p ) · k ) V LOµν = f + ( t (cid:48) ) p − µ p − · k ( p − + k − p ) ν − f + ( t (cid:48) ) g µν + f + ( t (cid:48) ) − f + ( t )( p − + p ) · k ( p − + p ) µ ( p − p − ) ν . In addition, the full vector contribution to the radiative amplitude contains model-dependent gauge-invariantterms of O ( k ) , such that V µν = V LOµν + (cid:98) V µν . Low’s soft-photon theorem [27] is manifestly satisfied since [12]: V LOµν = f + ( t ) p − µ p − · k ( p − − p ) ν + f + ( t ) (cid:18) p − µ k ν p − · k − g µν (cid:19) + 2 df + ( t ) dt (cid:18) p − µ p · kp − · k − p µ (cid:19) ( p − − p ) ν + O ( k ) . Since axial-vector contributions A µν are not present in the non-radiative amplitude, they start at O ( k ) in thephoton energy expansion. Thus, they are model-dependent and must be manifestly gauge-invariant: k µ A µν = 0 .Structure-dependent vector and axial-vector contributions to τ − → ν τ π − π (cid:96) + (cid:96) − decays, which reproduce thecorresponding terms of the radiative decay in the case of the real photon ( k = 0 ), are considered in section 4.
3. The τ − → ν τ π − π (cid:96) + (cid:96) − amplitude: model-independent contribution In this section we consider the leading terms (in the virtual photon momentum expansion) of the amplitudefor lepton-pair production in the di-pion tau lepton decay. As stated before, they depend only upon the formfactors and electromagnetic properties of the particles involved in the non-radiative amplitude and are fixed fromthe gauge invariance requirement. The contributions to this part of the amplitude are given by the diagramsshown in Fig. 1. τ − (cid:96) + (cid:96) − ν τ π − π γ W − (a) τ − ν τ π π − (cid:96) + (cid:96) − W − γ (b) τ − ν τ π π − (cid:96) + (cid:96) − W − γ (c) Figure 1: Feynman diagrams contributing to the τ − → ν τ π − π (cid:96) + (cid:96) − decay: (a) Inner Bremsstrahlung (IB) offthe tau lepton; (b) IB off the π − meson, (c) contact term.The matrix element for the decay τ − ( P ) → ν τ ( q ) π − ( p − ) π ( p ) (cid:96) + ( p (cid:96) + ) (cid:96) − ( p (cid:96) − ) has a similar structure to theradiative amplitude: M ( (cid:96)(cid:96) ) = e k G F V ∗ ud l µ (cid:40) ( p − − p ) ν k − P · k f + ( t ) u ν ( q ) γ ν (1 − γ ) (cid:0) /P − /k + M τ (cid:1) γ µ u τ ( P )+ ( V µν − A µν ) L ν (cid:41) , (3)where we have defined the weak L ν = u ν ( q ) γ ν (1 − γ ) u τ ( P ) and electromagnetic l µ = u ( p (cid:96) − ) γ µ v ( p (cid:96) + ) leptoniccurrents. Owing to Dirac’s equation and k = p (cid:96) + + p (cid:96) − , we have l µ k µ = 0 . The variable t is still defined as theinvariant mass of the two-pion system, however for virtual photons t (cid:48) ≡ ( P − q ) = t + 2( p − + p ) · k + k .3s in the case of the radiative amplitude, we split the vector contribution into two terms: V µν = V LOµν + (cid:98) V µν . (4)The leading order LO (model-independent) contribution is fixed from the diagrams in Fig. 1 and the gauge-invariance requirement. One gets: V LOµν = f + ( t (cid:48) ) 2 p − µ + k µ p − · k + k ( p − + k − p ) ν − f + ( t (cid:48) ) g µν + (cid:18) f + ( t (cid:48) ) − f + ( t )2( p + p − ) · k + k (cid:19) (cid:16) p − + p ) µ + k µ (cid:17) ( p − − p ) ν . (5)The terms proportional to k µ can be dropped ( k µ l µ = 0 ) in the previous expression owing to the conservation ofthe electromagnetic current. Note also the presence of k terms appearing in denominators of the propagatorsof charged particles.This model-independent amplitude, also known as the inner bremsstrahlung term in this paper, is the leadingcontribution at low photon momenta given the off-shell propagators of charged particles, and the enhancementfactor provided by the photon pole propagator. This feature makes model-independent contributions moreimportant for observables in e + e − , which has a lower threshold than for muon-pair production. On the otherhand, model or structure-dependent contributions, described by (cid:98) V µν − A µν terms, start at O ( k ) and may becomeimportant for large photon momenta. The effects of model-dependent contributions become more visible in the µ + µ − pair production, allowing to explore the rich dynamics of strong interactions in the intermediate energyregime. We will consider the later contributions in the following section.
4. Structure-dependent contributions
The structure-dependent piece of vector contributions to the amplitude, (cid:98) V µν , as well as axial-vector contribu-tions A µν are computed in the framework of resonance chiral theory. They are analogous to similar amplitudesin the radiative di-pion decays of tau leptons, but now they involve a virtual photon ( k = ( p (cid:96) + + p (cid:96) − ) ).The hadronic vertex W ∗ ( W ) → π − ( p − ) π ( p ) γ ∗ ( k ) with two virtual gauge bosons (characterized by W = t (cid:48) = ( p − + p + k ) (cid:54) = m W , k = ( p (cid:96) + + p (cid:96) − ) (cid:54) = 0 ) can be parameterized in terms of two different sets of formfactors, vector and axial-vector. In addition to the squared momenta of these virtual bosons ( t (cid:48) , k ), the formfactors can depend upon two independent kinematical variables which can be taken as ( t, s − ) or ( t, s ), where t = ( p − + p ) is the invariant mass of the di-pion system and s − , ≡ ( p − , + k ) , or equivalently p − · k or p · k .Once we chose ( t, t (cid:48) , k ) as relevant variables, either p − · k or p · k can be chosen as the remaining kinematicalscalar to describe the hadronic vertex.In this section we discuss the Lorentz structure of the vector and axial-vector contributions to the hadronicvertex. We first compute the expressions for the vector form factors within RChT and end this subsection bydiscussing the short-distance QCD constraints on the resonance couplings. In the second part we comment onthe axial-vector contributions. Since they play a subleading role numerically (as checked in the real photon case[6, 7, 5]), we will consider only the terms arising from the axial anomaly and the Wess-Zumino contributionsfor the case of real photons as computed in [12]. As stated before, the most general form of the structure-dependent part of the vector contributions to thehadronic vertex can be built out of the metric tensor g µν and the three independent momenta p − , p and k byimposing gauge-invariance. Including the leading order terms (5), we get : V µν = v (cid:0) p − · kg µν − p µ − k ν (cid:1) + v ( p · kg µν − p µ k ν ) + v (cid:0) p · kp µ − − ( p − · kp µ (cid:1) p ν − + v (cid:0) p · kp µ − − p − · kp µ (cid:1) ( p − + p + k ) ν + v (cid:0) k g µν − k µ k ν (cid:1) + v (cid:0) k p µ − − p − · kk µ (cid:1) p ν + v (cid:0) p · kk µ − k p µ (cid:1) p ν − − f + ( t (cid:48) ) g µν + f + ( t (cid:48) ) − f + ( t )2 ( p + p − ) · k + k (2( p + p − ) + k µ ) µ ( p − p − ) ν + f + ( t (cid:48) )2 k · p − + k (cid:0) p µ − + k µ (cid:1) ( p − + k − p ) ν . (6)4he Lorentz-invariant form factors v , ··· , depend in general upon the four invariant variables discussed above.They encode the information about the dynamics of the strong, weak and electromagnetic interactions involvedin the W ∗ π − π γ ∗ vertex. They are the coefficients of (explicitly gauge-invariant) Lorentz vector tensor struc-tures. In the case of a real photon, v , , do not contribute to the amplitude given that k = 0 and (cid:15) · k = 0 and one recovers the results of Ref. [12] for the radiative amplitude, as it should be. Although terms pro-portional to k µ do not contribute to physical results owing to current conservation, we keep them becauseexplicit calculations of the hadronic vertex gives rise to such structures and in order to exhibit explicitly gaugeinvariance.Next, we compute the different Feynman diagrams appearing in Figure 2 within the RChT framework [3, 4],which ensures the low-energy behaviour of ChPT [28] and includes resonances as dynamical degrees of freedomupon their approximate U (3) flavor symmetry.Besides the kinetic terms for the resonances (that we do not quote), we have used the interaction Lagrangiangiven by (see ref. [3] for further details) L V = F V √ T r (cid:16) (cid:101) V µν f µν + (cid:17) + i G V √ T r (cid:16) (cid:101) V µν u µ u ν (cid:17) , L A = F A √ T r (cid:16) (cid:101) A µν f µν − (cid:17) , where T r stands for a trace in flavor space. Resonance fields are represented by the antisymmetric tensors[31, 32] (cid:101) V µν and (cid:101) A µν (not to be confused with the tensor contributions to the decay amplitude defined in Eq.(1)), the coupling to the weak charged V − A current proceeds through the f µν ± tensors, and the u µ tensorscouple the resonances to either the vector part of the W boson or (derivatively) to pion fields. The couplingconstants of resonances F V , G V and F A can be fixed from short-distance constraints in terms of f = F π ∼ MeV. Short-distance QCD constraints on the spin-one correlators [3, 4, 33] predict the former in terms of thelatter as F V = √ f, G V = f / √ , F A = f . The associated uncertainties are discussed below.In figure 2 we show the Feynman diagrams that contribute to the vector tensor amplitudes in Eq. (6). Forconvenience and later comparison, we quote here the results in the real photon case [12] and defer to the end ofthis subsection their generalization for a virtual photon: v = F V G V f m ρ (cid:2) m ρ D − ρ ( t (cid:48) ) + tD − ρ ( t ) + tm ρ D − ρ ( t ) D − ρ ( t (cid:48) ) (cid:3) + F V f m ρ (cid:2) − − m ρ D − ρ ( t (cid:48) ) + t (cid:48) D − ρ ( t (cid:48) ) (cid:3) + F A f m a (cid:20) m a − m π + t (cid:21) D − a (cid:2) ( p − + k ) (cid:3) ,v = F V G V tf m ρ (cid:2) − D − ρ ( t ) − m ρ D − ρ ( t ) D − ρ ( t (cid:48) ) (cid:3) + F V f m ρ (cid:2) − − m ρ D − ρ ( t (cid:48) ) − t (cid:48) D − ρ ( t (cid:48) ) (cid:3) + F A f m a (cid:2) m a − m π − p − · k (cid:3) D − a (cid:2) ( p − + k ) (cid:3) ,v = F A f m a D − a (cid:2) ( p − + k ) (cid:3) ,v = − F V G V f D − ρ ( t ) D − ρ ( t (cid:48) ) + F V f m ρ D − ρ ( t (cid:48) ) , (7)where the propagators are given by D ρ ( s ) = m ρ − s − im ρ Γ ρ ( s ) ,D a ( s ) = m a − s − im a Γ a ( s ) . (8)The off-shell widths of the ρ (770) and a (1260) mesons that appear in the above expressions and are used inthis paper, are obtained within RChT. In the first case it includes the ππ and K ¯ K cuts [34, 35] and in thesecond case the π [36, 37] and K ¯ Kπ [38] cuts.In the case of a virtual photon, we have contributions from the same Feynman diagrams shown in figure2, taking due care of k (cid:54) = 0 . As in the case of a real photon, some of these contributions would appear in5 π − γ ∗ (a) γ ∗ π π − (b) γ ∗ π π − (c) π γ ∗ π − (d) π γ ∗ π − (e) γ ∗ π π − (f) γ ∗ π π − (g) π γ ∗ π − a − (h) Figure 2: Contributions to V µν in RChT of the hadronic vertex W ∗− → π − π γ ∗ vertices. Those involvingresonances are highlighted with a thick dot. Insertion of the weak charged current is represented by the squaredot. Resonances ( ρ unless specified) are represented with a double line and γ ∗ stands for γ ∗ → (cid:96) + (cid:96) − .the leading order term V LOµν given in Eq. (5). As it was discussed above, in this case the structure-dependentcontributions can be described in terms of seven form factors v , ··· , , defined in Eq. (6) as the coefficients ofgauge-invariant structures. An explicit evaluation of them leads to: v = F V G V f (cid:2) D − ρ (cid:0) k (cid:1) + 2 D − ρ ( t (cid:48) ) + tD − ρ ( t ) D − ρ (cid:0) k (cid:1) + tD − ρ ( t ) D − ρ ( t (cid:48) ) (cid:3) + F V f (cid:2) − D − ρ (cid:0) k (cid:1) − D − ρ ( t (cid:48) ) + (cid:0) t (cid:48) − k (cid:1) D − ρ ( t (cid:48) ) D − ρ (cid:0) k (cid:1)(cid:3) + F A f m a (cid:18) m a − m π + t (cid:19) D − a (cid:104) ( p − + k ) (cid:105) ,v = F V G V tf (cid:2) − D − ρ ( t ) D − ρ (cid:0) k (cid:1) − D − ρ ( t ) D − ρ ( t (cid:48) ) (cid:3) + F V f (cid:2) − D − ρ (cid:0) k (cid:1) − D − ρ ( t (cid:48) ) − (cid:0) t (cid:48) − k (cid:1) D − ρ ( t (cid:48) ) D − ρ (cid:0) k (cid:1)(cid:3) + F A f m a (cid:0) m a − m π − k · p − (cid:1) D − a (cid:104) ( p − + k ) (cid:105) ,v = F A f m a D − a (cid:104) ( p − + k ) (cid:105) ,v = F V f D − ρ ( t (cid:48) ) D − ρ (cid:0) k (cid:1) − F V G V f D − ρ ( t ) D − ρ ( t (cid:48) ) , = F V f (cid:2) − D − ρ (cid:0) k (cid:1) − D − ρ ( t (cid:48) ) − k · ( p − p − ) D − ρ ( t (cid:48) ) D − ρ (cid:0) k (cid:1)(cid:3) + F A f m a (cid:18) m a − m π + t (cid:19) D − a (cid:104) ( p − + k ) (cid:105) ,v = 2 F V G V f D − ρ ( t ) D − ρ (cid:0) k (cid:1) + F V f D − ρ ( t (cid:48) ) D − ρ (cid:0) k (cid:1) ,v = 2 F V G V f D − ρ ( t ) D − ρ (cid:0) k (cid:1) + F V f D − ρ ( t (cid:48) ) D − ρ (cid:0) k (cid:1) + F A f m a D − a (cid:104) ( p − + k ) (cid:105) . (9)Note that the above expressions reduce to the corresponding Eqs. (7) for v , ··· , in the case of real photons( k → ). Also, we observe that, owing to the conservation of electromagnetic current, the form factors v , , do not contribute to the vector tensor terms in Eq. (6) in the real photon case.The short-distance constraints for two-point Green functions (which include the set of relations F V = √ f , G V = f / √ , F A = f ) get modified when including three-point Green functions in both intrinsic parity sectors[39, 40, 41, 42, 43, 44, 38, 36, 45, 46, 47, 48, 49, 50, 51, 52, 53, 5]. The consistent set of relations in thismore general case includes F V = √ f [49] that –through the appropriate asymptotic behaviour of the spin-onecorrelators– implies, G V = f / √ and F A = √ f . As studied extensively in ref. [5] for the τ − → ν τ π − π γ decays, shifting from F V = √ f, G V = f √ , F A = f , (10)to F V = √ f, G V = f √ , F A = √ f , (11)gives a rough estimate of the uncertainty in the calculation with the interaction Lagrangian (7) due to missinghigher-order terms in the chiral expansion. We will take relations (10) as the reference ones but evaluatealternatively with (11) to assess our model-dependent error. The most general form of the axial-vector weak current contribution in τ − ( P ) → ν τ ( q ) π − ( p − ) π ( p ) (cid:96) + ( p (cid:96) + ) (cid:96) − ( p (cid:96) − ) can be built out of the rank-four antisymmetric Levi-Civita tensor and the three independent momenta in the W ∗ ( W ) → π − ( p − ) π ( p ) γ ∗ ( k ) vertex. Making use of Schouten’s identity and the gauge invariance condition k µ A µν = 0 , one gets [54, 12, 51] the same result as in the real photon case: A µν = ia ε µνρσ ( p − p − ) ρ k σ + ia W ν (cid:15) µλρσ p − λ p ρ k σ + ia ε µνρσ k ρ W σ + ia ( p + k ) ν ε µλρσ p − λ p ρ k σ , (12)where W = P − q = p − + p + k . As in the vector tensor case, the axial-vector form factors a ··· are Lorentzinvariant functions that depend upon two kinematical Lorentz invariants (in addition to t (cid:48) = W and k ).At O ( p ) only a and a , from the Wess-Zumino-Witten functional ([55, 56]), contribute and are obviouslythe same as in the real photon case at this order (see Fig. 3). They are given by [12] a = (cid:2) π f (cid:3) − ,a = − (cid:2) π f (cid:0) t (cid:48) − m π (cid:1)(cid:3) − . (13)As in the case of radiative τ − → ν τ π − π γ decays, we expect the corresponding contributions to the decayobservables in lepton pair production to be negligibly small [6, 7, 5]. Therefore, we do not consider necessaryto compute all remaining axial-vector contributions which will introduce in addition further (although small)uncertainties in our computation. 7 π − π γπ − π Figure 3: Contributions from the Wess-Zumino-Witten to the W ∗− → π − π γ ∗ vertex [55] [56]. The weakcharged current is represented by the square dot.
5. Branching ratio and lepton-pair spectrum
As it is well known, the unpolarized squared amplitude of a five-body decay as τ − → ν τ π − π (cid:96) + (cid:96) − , dependson eight independent kinematical variables. Depending upon the specific observable we are interested in, it willbecome necessary to integrate some or all of these kinematical variables. In this paper we find convenient to usethe set of invariant variables described in Ref. [57] (see Appendix A in which we have defined these variablesand have calculated one of the non-trivial scalar products, a subtlety when there are more than four particlesin the final state for decay processes). More specifically, we will compute the invariant mass distribution of thelepton pair ( k -distribution) and the corresponding branching fraction for the tau decays under consideration.The kinematical domain of the lepton pair distribution is the interval (cid:104) m (cid:96) , ( M τ − m π ) (cid:105) , being (cid:96) = e, µ .In order to distinguish among the different contributions, we split the total decay observables into three terms:1) the IB piece, i.e. the inner bremsstrahlung or model-independent contributions obtained with vanishing v , ··· , and a , ··· , form factors; 2) the VV (AA) model-dependent part, corresponding to the terms with non-vanishingvector (axial-vector) form factors and, 3) the IB-V, IB-A and V-A pieces, which correspond to the interferencesof IB, V and A contributions.Table 1 shows the results of different contributions to the branching ratios of e + e − and µ + µ − pair production(the numerical errors in the integration are shown within parentheses). These calculations were obtained usingthe vector form factors given in Eq. (9), the axial-vector form factors of Eq. (13) and the short-distancecontraints exhibited in Eq. (10); for comparison, we also show within square brackets the results obtained usingthe vector form factors of Eq. (7), corresponding to the real photon case. In the third (fifth) column of Table1 we also show the results for various contributions to the branching fraction obtained using the vector (9)and axial-vector (13) form factors, but subject to the short-distance constraints on the couplings constants ofresonances shown in Eq. (11). As it was explained at the end of Section 4.1, shifting the values of couplingconstants according to the prescriptions on short-distance constraints allows us to assess an important part oftheoretical uncertainties.The results shown in Table 1 exhibit the suppression expected since lepton-pair production is O ( α ) withrespect to the dominant τ − → π − π ν τ decay . Also, the µ + µ − pair production is further suppressed withrespect to e + e − production given that the later is largely dominated by model-independent contributions,which are enhanced and peaked at lower invariant mass values of the lepton-pair invariant mass due to thevirtual photon propagator. As pointed out before, the axial-vector contributions are suppressed in all cases.Our final predictions for the branching fractions are: BR( τ − → ν τ π − π e + e − ) = (2 . ± . × − , (14) BR( τ − → ν τ π − π µ + µ − ) = (1 . ± . × − . (15)The associated errors cover the results shown in the different columns of Table 1 .The normalized (to the total τ decay width) lepton-pair invariant mass distributions are shown in Figure4 for electron-positron and in Figure 5 for µ + µ − production. Both distributions are peaked very close tothe corresponding threshold ( k = 4 m (cid:96) ) for lepton pair production, with an enhanced peaking for e + e − production, due to the /k dependence of the squared amplitude. The second peak in the plots corresponds Noteworthy, the large inner bremsstrahlung contributions coming from photon emission off the τ − lepton and π − meson, isalmost cancelled by their interference, which yields physical results (e. g., branching ratios of order α ). This type of cancellationagrees with that observed in τ − → ν τ π − (cid:96) + (cid:96) − decays [22]. We note that these errors are larger than /N C , typical of a large- N C expansion, for the structure-dependent contributions. + (cid:96) − = e + e − (cid:96) + (cid:96) − = µ + µ − Contribution (cid:96) + (cid:96) − = e + e − using (11) for (cid:96) + (cid:96) − = µ + µ − using (11) for F V , F A and G V F V , F A and G V IB 2.213(11) × − × − (cid:2) × − (cid:3) (cid:2) × − (cid:3) VV 6.745(36) × − × − × − × − (cid:2) × − (cid:3) (cid:2) × − (cid:3) AA 1.91(1) × − × − (cid:2) × − (cid:3) (cid:2) × − (cid:3) IB-V − × − -1.02(18) × − × − × − (cid:2) − × − (cid:3) (cid:2) × − (cid:3) IB-A 9.1(4.5) × − × − (cid:2) × − (cid:3) (cid:2) × − (cid:3) V-A 5.2(2.1) × − × − − × − -1.65(5) × − (cid:2) × − (cid:3) (cid:2) − × − (cid:3) Total 2.245(13) × − × − × − × − (cid:2) × − (cid:3) (cid:2) × − (cid:3) Table 1:Contributions to the branching ratio of τ − → ν τ π − π (cid:96) + (cid:96) − decays. IB, VV, and AA stand for the InnerBremsstrahlung, Vector, and Axial-Vector contributions, respectively, while IB-V, IB-A, and V-A correspond totheir interferences. Columns three and five display the branching ratios obtained using Eq. (11) for the relationsof resonance couplings with the pion decay constant entering the structure-dependent vector form factors, whilethe second and four columns correspond to the use of relations (10).9igure 4: Contributions to the normalized invariant mass distribution for e + e − pair production (interferencesare not displayed). A double logarithmic scale was used. The second peak is due to the ρ (770) dominance ofthe virtual photon propagator.to the ρ (770) − γ ∗ couplings dominance in the vector form factors. It is also clear that the model-dependentcontributions are more visible in the µ + µ − than in e + e − production, which is also related to the suppressionof inner bremsstrahlung for large photon virtualities.
6. Conclusions
We have calculated for the first time the branching ratios and lepton-pair mass distributions of the five-body decays τ − → ν τ π − π (cid:96) + (cid:96) − ( (cid:96) = e, µ ). As expected, these observables are of O ( α ) with respect to thecorresponding dominant di-pion τ lepton decay. For the (cid:96) = e case, a clear inner bremsstrahlung (IB) dominanceis observed due to the small (cid:96) + (cid:96) − invariant mass ( k ) threshold values. On the other hand, for (cid:96) = µ , bothcontributions, structure-dependent and IB, are of the same order.The structure-dependent contributions corresponding to the W − → π − π γ ∗ effective vertex, were calcu-lated using the Resonance Chiral Theory framework. Such an approach considers the lightest resonances asactive degrees of freedom giving the low-energy chiral limit of QCD and ensuring an appropriate short-distancebehaviour. The structure-dependent vector form factors coincide (in the limit k → ) with their counterpartscomputed in the case of the radiative τ − → ν τ π − π γ decays [12]. We expect axial-vector structure-dependentcontributions to be negligible and we stick to their values provided by the Wess-Zumino-Witten anomalousterms.Within this framework, we get BR (cid:0) τ − → ν τ π − π e + e − (cid:1) = 2 . × − (which is essentially free ofhadronic uncertainties) and BR (cid:0) τ − → ν τ π − π µ + µ − (cid:1) = 1 . × − . The estimated theoretical uncertain-ties are associated to different relations between resonances couplings and the pion decay constant, obtainedfrom the short-distance behavior of two- and three-point Green functions. While the branching fraction for e + e − channel allows to conclude that it could be discovered already with BaBar or Belle data, the µ + µ − casewill challenge the capabilities of Belle-II. On the other hand, the measurement of the µ + µ − spectrum, whichis more sensitive to structure-dependent contributions, can be useful to test previous calculations of radiative10igure 5: Contributions to the normalized invariant mass distribution for µ + µ − pair production (inteferencesare not displayed). The second peak is due to the ρ (770) dominance of the virtual photon propagator.corrections to di-pion tau lepton decays. Therefore, it has the potential of reducing the uncertainties on thedominant piece of the hadronic vacuum polarization part of a µ using tau data.Finally, the addition of the matrix elements derived in this work to the Monte Carlo generator TAUOLA[58] will be useful in improving background rejection for searches of three-prong lepton flavor or lepton numberviolating tau decays. Acknowledgements
J. L. G. S. thanks Conacyt for his Ph. D. scholarship. G. L. C. received support from Ciencia de FronteraConacyt project No. 428218. P. R. is indebted for the funding received through Fondo SEP-Cinvestav 2018(project No. 142) and Cátedra Marcos Moshinsky (2020).
Appendix A. Five-body kinematics
The kinematics of the five-body decay process τ − ( P ) → ν τ ( q ) π − ( p − ) π ( p ) (cid:96) + ( p (cid:96) + ) (cid:96) − ( p (cid:96) − ) where the leptonpair (cid:96) + (cid:96) − is either e + e − or µ + µ − and k = (cid:96) + + (cid:96) − is the momentum of the virtual photon, is described interms of eight independent variables. All the scalar products of two-momenta can be written in terms of theseindependent kinematical variables. Following reference [57] we choose these variables as follows: s = ( P − q ) , s = ( P − q − p − ) , s = ( P − q − p − − p ) ,u = ( P − p − ) , u = ( P − p ) , u = ( P − p (cid:96) − ) ,t = ( P − p − − p ) , t = ( P − p − − p − p (cid:96) − ) , and the auxiliary variables s = M τ , s = m (cid:96) + , u = s and t = u .In general, for decay processes with n particles in the final state, it can be shown that we will have n − independent invariants. In the case n ≥ , it is well known that some scalar products cannot be written directly11n terms of the s i , t i and u i variables. This is the case for the p − · p (cid:96) + and p · p (cid:96) + scalar products. Followingreference [57] and making use of symmetry considerations, p − · p (cid:96) + reads p − · p (cid:96) + = a (cid:0) P · p − − q · p − − p − − p − · p (cid:1) + b (cid:0) P · p − − p − − p − · p (cid:1) + c ( P · p − ) . (A.1)with a , b and c given as follows, a = 1 Z (cid:2) − ADF + AEs − BDE + BF t + C (cid:0) D − M τ t (cid:1)(cid:3) ,b = 1 Z (cid:2) − B ( AF + CD ) + M τ AC + B E + DF s − M τ Es (cid:3) ,c = 1 Z (cid:2) A F − A ( BE + CD ) + BCt + DEs − F s t (cid:3) , where the capital letters are defined in terms of the already known scalar products in the following way: A = ( P − q − p − − p ) · ( P − p − − p ) ,B = P · ( P − q − p − − p ) ,C = p (cid:96) + · ( P − q − p − − p ) ,D = P · ( P − p − − p ) ,E = p (cid:96) + · ( P − p − − p ) ,F = P · p (cid:96) + ,Z = M τ A − ABD + t (cid:0) B − M τ s (cid:1) + D s . Then, once we have calculated p − · p (cid:96) + it is straighforward to obtain (a very lenghty expression for) p · p (cid:96) + ,which we do not quote. References [1] M. Davier, A. Hocker, and Z. Zhang, “The Physics of Hadronic Tau Decays,”
Rev. Mod. Phys. , vol. 78, pp. 1043–1109, 2006.[2] A. Pich, “Precision Tau Physics,”
Prog. Part. Nucl. Phys. , vol. 75, pp. 41–85, 2014.[3] G. Ecker, J. Gasser, A. Pich, and E. de Rafael, “The Role of Resonances in Chiral Perturbation Theory,”
Nucl. Phys. B ,vol. 321, pp. 311–342, 1989.[4] G. Ecker, J. Gasser, H. Leutwyler, A. Pich, and E. de Rafael, “Chiral Lagrangians for Massive Spin 1 Fields,”
Phys. Lett. B ,vol. 223, pp. 425–432, 1989.[5] J. Miranda and P. Roig, “New τ -based evaluation of the hadronic contribution to the vacuum polarization piece of the muonanomalous magnetic moment,” .[6] A. Flores-Tlalpa, G. Lopez Castro, and G. Sanchez Toledo, “Radiative two-pion decay of the tau lepton,” Phys. Rev. D , vol. 72,p. 113003, 2005.[7] F. Flores-Baez, A. Flores-Tlalpa, G. Lopez Castro, and G. Toledo Sanchez, “Long-distance radiative corrections to the di-piontau lepton decay,”
Phys. Rev. D , vol. 74, p. 071301, 2006.[8] F. Jegerlehner and R. Szafron, “ ρ − γ mixing in the neutral channel pion form factor F eπ and its role in comparing e + e − with τ spectral functions,” Eur. Phys. J. C , vol. 71, p. 1632, 2011.[9] M. Benayoun, L. Delbuono, and F. Jegerlehner, “BHLS , a New Breaking of the HLS Model and its Phenomenology,” Eur.Phys. J. C , vol. 80, no. 2, p. 81, 2020. [Erratum: Eur.Phys.J.C 80, 244 (2020)].[10] R. Alemany, M. Davier, and A. Hocker, “Improved determination of the hadronic contribution to the muon (g-2) and to alpha(M(z)) using new data from hadronic tau decays,”
Eur. Phys. J. C , vol. 2, pp. 123–135, 1998.[11] V. Cirigliano, G. Ecker, and H. Neufeld, “Isospin violation and the magnetic moment of the muon,”
Phys. Lett. B , vol. 513,pp. 361–370, 2001.[12] V. Cirigliano, G. Ecker, and H. Neufeld, “Radiative tau decay and the magnetic moment of the muon,”
JHEP , vol. 08, p. 002,2002.[13] A. Flores-Tlalpa, F. Flores-Baez, G. Lopez Castro, and G. Toledo Sanchez, “Model-dependent radiative corrections to τ − → π − π ν revisited,” Nucl. Phys. B Proc. Suppl. , vol. 169, pp. 250–254, 2007.[14] M. Davier, A. Hoecker, G. Lopez Castro, B. Malaescu, X. Mo, G. Toledo Sanchez, P. Wang, C. Yuan, and Z. Zhang, “TheDiscrepancy Between tau and e+e- Spectral Functions Revisited and the Consequences for the Muon Magnetic Anomaly,”
Eur. Phys. J. C , vol. 66, pp. 127–136, 2010.[15] M. Davier, A. Hoecker, B. Malaescu, and Z. Zhang, “Reevaluation of the Hadronic Contributions to the Muon g-2 and toalpha(MZ),”
Eur. Phys. J. C , vol. 71, p. 1515, 2011. [Erratum: Eur.Phys.J.C 72, 1874 (2012)].[16] M. Davier, A. Höcker, B. Malaescu, C.-Z. Yuan, and Z. Zhang, “Update of the ALEPH non-strange spectral functions fromhadronic τ decays,” Eur. Phys. J. C , vol. 74, no. 3, p. 2803, 2014.[17] P. Zyla et al. , “Review of Particle Physics,”
PTEP , vol. 2020, no. 8, p. 083C01, 2020.[18] G. Hernández-Tomé, G. López Castro, and P. Roig, “Flavor violating leptonic decays of τ and µ leptons in the Standard Modelwith massive neutrinos,” Eur. Phys. J. C , vol. 79, no. 1, p. 84, 2019. [Erratum: Eur.Phys.J.C 80, 438 (2020)].
19] P. Blackstone, M. Fael, and E. Passemar, “ τ → µµµ at a rate of one out of tau decays?,” Eur. Phys. J. C , vol. 80, no. 6,p. 506, 2020.[20] G. Lopez Castro and N. Quintero, “Lepton number violating four-body tau lepton decays,”
Phys. Rev. D , vol. 85, p. 076006,2012. [Erratum: Phys.Rev.D 86, 079904 (2012)].[21] O. Shekhovtsova, T. Przedzinski, P. Roig, and Z. Was, “Resonance chiral Lagrangian currents and τ decay Monte Carlo,” Phys. Rev. D , vol. 86, p. 113008, 2012.[22] A. Guevara, G. López Castro, and P. Roig, “Weak radiative pion vertex in τ − → π − ν τ (cid:96) + (cid:96) − decays,” Phys. Rev. , vol. D88,no. 3, p. 033007, 2013.[23] A. Flores-Tlalpa, G. López Castro, and P. Roig, “Five-body leptonic decays of muon and tau leptons,”
JHEP , vol. 04, p. 185,2016.[24] S. Antropov, S. Banerjee, Z. Wąs, and J. Zaremba, “TAUOLA Update for Decay Channels with e + e − Pairs in the FinalState,” .[25] T. Aoyama et al. , “The anomalous magnetic moment of the muon in the Standard Model,”
Phys. Rept. , vol. 887, pp. 1–166,2020.[26] W. Altmannshofer et al. , “The Belle II Physics Book,”
PTEP , vol. 2019, no. 12, p. 123C01, 2019. [Erratum: PTEP 2020,029201 (2020)].[27] F. Low, “Bremsstrahlung of very low-energy quanta in elementary particle collisions,”
Phys. Rev. , vol. 110, pp. 974–977, 1958.[28] J. Gasser and H. Leutwyler, “Chiral Perturbation Theory to One Loop,”
Annals Phys. , vol. 158, p. 142, 1984.[29] D. Gómez Dumm and P. Roig, “Dispersive representation of the pion vector form factor in τ → ππν τ decays,” Eur. Phys. J.C , vol. 73, no. 8, p. 2528, 2013.[30] S. Gonzàlez-Solís and P. Roig, “A dispersive analysis of the pion vector form factor and τ − → K − K S ν τ decay,” Eur. Phys.J. C , vol. 79, no. 5, p. 436, 2019.[31] E. Kyriakopoulos, “Vector-meson interaction hamiltonian,”
Phys. Rev. , vol. 183, pp. 1318–1323, Jul 1969.[32] Y. Takahashi and R. Palmer, “Gauge-independent formulation of a massive field with spin one,”
Phys. Rev. D , vol. 1, pp. 2974–2976, May 1970.[33] S. Weinberg, “Precise relations between the spectra of vector and axial vector mesons,”
Phys. Rev. Lett. , vol. 18, pp. 507–509,1967.[34] F. Guerrero and A. Pich, “Effective field theory description of the pion form-factor,”
Phys. Lett. B , vol. 412, pp. 382–388,1997.[35] D. Gomez Dumm, A. Pich, and J. Portoles, “The Hadronic off-shell width of meson resonances,”
Phys. Rev. D , vol. 62,p. 054014, 2000.[36] D. Dumm, P. Roig, A. Pich, and J. Portoles, “ τ → πππν τ decays and the a (1260) off-shell width revisited,” Phys. Lett. B ,vol. 685, pp. 158–164, 2010.[37] I. Nugent, T. Przedzinski, P. Roig, O. Shekhovtsova, and Z. Was, “Resonance chiral Lagrangian currents and experimentaldata for τ − → π − π − π + ν τ ,” Phys. Rev. D , vol. 88, p. 093012, 2013.[38] D. Dumm, P. Roig, A. Pich, and J. Portoles, “Hadron structure in τ → KKπν τ decays,” Phys. Rev. D , vol. 81, p. 034031,2010.[39] P. Ruiz-Femenia, A. Pich, and J. Portoles, “Odd intrinsic parity processes within the resonance effective theory of QCD,”
JHEP , vol. 07, p. 003, 2003.[40] V. Cirigliano, G. Ecker, M. Eidemuller, A. Pich, and J. Portoles, “The < VAP > Green function in the resonance region,”
Phys. Lett. B , vol. 596, pp. 96–106, 2004.[41] V. Cirigliano, G. Ecker, M. Eidemuller, R. Kaiser, A. Pich, and J. Portoles, “The < SPP > Green function and SU(3) breakingin K(l3) decays,”
JHEP , vol. 04, p. 006, 2005.[42] V. Cirigliano, G. Ecker, M. Eidemuller, R. Kaiser, A. Pich, and J. Portoles, “Towards a consistent estimate of the chirallow-energy constants,”
Nucl. Phys. B , vol. 753, pp. 139–177, 2006.[43] V. Mateu and J. Portoles, “Form-factors in radiative pion decay,”
Eur. Phys. J. C , vol. 52, pp. 325–338, 2007.[44] Z.-H. Guo, “Study of τ − → V P − ν τ in the framework of resonance chiral theory,” Phys. Rev. D , vol. 78, p. 033004, 2008.[45] Z.-H. Guo and P. Roig, “One meson radiative tau decays,”
Phys. Rev. D , vol. 82, p. 113016, 2010.[46] K. Kampf and J. Novotny, “Resonance saturation in the odd-intrinsic parity sector of low-energy QCD,”
Phys. Rev. D , vol. 84,p. 014036, 2011.[47] D. Gomez Dumm and P. Roig, “Resonance Chiral Lagrangian analysis of τ − → η ( (cid:48) ) π − π ν τ decays,” Phys. Rev. D , vol. 86,p. 076009, 2012.[48] Y.-H. Chen, Z.-H. Guo, and H.-Q. Zheng, “Study of η − η (cid:48) mixing from radiative decay processes,” Phys. Rev. D , vol. 85,p. 054018, 2012.[49] P. Roig and J. J. Sanz Cillero, “Consistent high-energy constraints in the anomalous QCD sector,”
Phys. Lett. B , vol. 733,pp. 158–163, 2014.[50] P. Roig, A. Guevara, and G. López Castro, “
V V (cid:48) P form factors in resonance chiral theory and the π − η − η (cid:48) light-by-lightcontribution to the muon g − ,” Phys. Rev. D , vol. 89, no. 7, p. 073016, 2014.[51] A. Guevara, G. López-Castro, and P. Roig, “ τ − → η ( (cid:48) ) π − ν τ γ decays as backgrounds in the search for second class currents,” Phys. Rev. D , vol. 95, no. 5, p. 054015, 2017.[52] A. Guevara, P. Roig, and J. Sanz-Cillero, “Pseudoscalar pole light-by-light contributions to the muon ( g − in ResonanceChiral Theory,” JHEP , vol. 06, p. 160, 2018.[53] L.-Y. Dai, J. Fuentes-Martín, and J. Portolés, “Scalar-involved three-point Green functions and their phenomenology,”
Phys.Rev. D , vol. 99, no. 11, p. 114015, 2019.[54] J. Bijnens, G. Ecker, and J. Gasser, “Radiative semileptonic kaon decays,”
Nucl. Phys. B , vol. 396, pp. 81–118, 1993.[55] J. Wess and B. Zumino, “Consequences of anomalous Ward identities,”
Phys. Lett. B , vol. 37, pp. 95–97, 1971.[56] E. Witten, “Global Aspects of Current Algebra,”
Nucl. Phys. B , vol. 223, pp. 422–432, 1983.[57] R. Kumar, “Covariant phase-space calculations of n-body decay and production processes,”
Phys. Rev. , vol. 185, pp. 1865–1875,1969.[58] S. Jadach, Z. Was, R. Decker, and J. H. Kuhn, “The tau decay library TAUOLA: Version 2.4,”
Comput. Phys. Commun. ,vol. 76, pp. 361–380, 1993.,vol. 76, pp. 361–380, 1993.