Radiative corrections of O ( α ) to B − → V 0 − ν̄ decays

The O(α) electromagnetic radiative corrections to the B− → V 0 −ν̄ (V is a vector meson and a charged lepton) decay rates are evaluated using the cutoff method to regularize virtual corrections and incorporating intermediate resonance states in the real-photon amplitude to extend the region of validity of the soft-photon approximation. The electromagnetic and weak form factors of hadrons are assumed to vary smoothly over the energies of virtual and real photons under consideration. The cutoff dependence of radiative corrections upon the scale that separates the longand short-distance regimes is found to be mild and is considered as an uncertainty of the calculation. Owing to partial cancellations of electromagnetic corrections evaluated over the threeand four-body regions of phase space, the photoninclusive corrected rates are found to be dominated by the short-distance contribution. These corrections will be relevant for a precise determination of the b quark mixing angles by testing isospin symmetry when measurements of semileptonic rates of charged and neutral B mesons at the few percent level become available. For completeness, we also provide numerical values of radiative corrections in the three-body region of the Dalitz plot distributions of these decays.


Introduction
A precise determination of the |V ub | quark mixing at a few percent level is crucial for future tests of the Standard Model (SM) picture of CP violation [1,2] as it fixes one of the sides of the db unitarity triangle.A combination of inclusive and exclusive channels of charmless leptonic and semileptonic decays of b hadrons can serve this purpose in two ways.On the one hand, they can provide a consistency check of the models employed to extract |V ub |; this in turn should yield reduced uncertainties by using the average of independent a e-mail: stostado@fis.cinvestav.mxb e-mail: glopez@fis.cinvestav.mx and consistent determinations.As an example, let us recall that the consistency check of the set of data for superallowed Fermi transitions (kaon semileptonic decays), which require a detailed calculation of all sources of isospin breaking corrections, allows a determination of |V ud | (|V us |) at the per mille (one percent) level [3].Given that different charmless semileptonic decays of charged and neutral B → (P, V ) transitions are related by isospin symmetry, similar consistency checks can be attempted with future measurements in order to reach a better accuracy in the extraction of |V ub |.
Several leptonic and semileptonic exclusive decays induced by the b → u ν transition can be used to extract |V ub | in an independent way.The cleanest channel owing to the better control on theoretical and experimental input is B → π ν , yielding |V ub | excl = (3.28 ± 0.29) × 10 −3 [3].This relies on the most precise measurements of the total and differential rates [4][5][6][7][8][9][10][11][12][13], as well as on the lattice [14,15] and the light-cone QCD sum rules (LCSR) [16][17][18] calculations 1 of the form factor.When compared to the most precise determination from inclusive b → u ν transitions one is led to discrepant results |V ub | = |V ub | incl − |V ub | excl = (1.13 ± 0.36) × 10 −3 [3], with errors added in quadrature.Improvement in both, experiment and theory, will play a crucial role in understanding and solving this discrepancy. 2On the other hand, as long as new measurements of B + → τ + ν τ [21][22][23][24][25] are getting in better agreement with the SM prediction [26], it will become a useful and independent test in the determination of |V ub |.Along this goal, a recent measurement of the ratio ( b → )/( b → c ) of partially integrated semileptonic rates by the LHCb collaboration [27], combined with the exclusive value of |V cb | [3] and Lattice QCD calculations of the relevant hadronic form 1 New LCSR calculations of the B → π form factors may actually further decrease the extracted value of |V ub | by up to 10 % [19]. 2 An interesting suggestion that duality violations may occur close to thresholds and would lead to smaller values of |V ub | incl , has been put forward recently [20].
In this work we are concerned with the calculation of electromagnetic radiative corrections to B ± → V 0 semileptonic transitions, a useful input in testing isospin symmetry in those B decays.In the isospin symmetry limit, the following relations hold 3 : which seem to be supported by experimental data within current experimental errors (see Table 1).Here we will assume the narrow width approximation for vector mesons; in practice, vector mesons are unstable particles that are reconstructed from the invariant-mass distribution of suitable decay channels [29][30][31].For instance, B ± → ρ 0 ± ν must be extracted by choosing a narrow window in the 2π invariant-mass distribution of the full B ± → π + π − ± ν process.Radiative corrections to this decay requires the consideration of additional electromagnetic interactions involving final state particles, which are not present in decays with neutral vector mesons.This and the modeling of the hadronic current increases considerably the difficulty of calculations.Therefore, considering vector mesons as asymptotic states is an intrinsic limitation of our present calculation.Departures from the isospin symmetry relations (1) are expected at the few percent level owing to effects of electromagnetism and u-d quark mass difference.The observable effects of isospin breaking manifest in hadron masses, hadronic form factors [29,40,41] and long-distance electromagnetic radiative corrections.Future improved measurements of semileptonic B → (π, ρ, ω) transitions will require a consideration of isospin breaking effects.
In this paper we contribute to this goal with the calculation of long-distance radiative corrections to B − → (ρ 0 , ω) − ν decays, which to the best of our knowledge have not been considered before.Our results can be applied to B − → D * 0 − ν decays as well, within the validity of the underlying assumptions.We compute both, the long-distance corrections to the 3 Similar relations hold by replacing ρ → π mesons.
photon-inclusive decay rates and the three-body region of the Dalitz plot distribution of these decays.These electromagnetic corrections should affect the equality between the former two decay rates in Eq. ( 1), but not the ratio between the latter two.Instead, ρ 0 -ω mixing, which originates mainly from the u-d quark mass difference, affects only the decay rates with (ρ 0 , ω) mesons in the final state [40,41].
Electromagnetic corrections for B meson decays are more complicated than in light hadrons due to the presence of hard photons.Here we deal with this problem by cutting photon momenta in virtual corrections at some relatively small cutoff scale where photon-meson interactions can be loosely approximated by scalar QED.For real-photon emission, we include the effects of resonance contributions in order to extend the validity of the soft-photon approximation.This sets a limitation of our calculation, which may weakens the validity of the approach for charmless semileptonic transitions but can work better for the charmfull case.

Radiative corrections to B − → V 0 semileptonic transitions
There are only a few works related to the calculation of radiative corrections to semileptonic decays of B mesons [42][43][44][45][46].Among the main reasons one finds the experimental difficulty to reach the few percent level accuracy in the measurements of the branching fractions and the theoretical limitation imposed by our knowledge of the meson-photon interactions in B meson decays.A first calculation in the case of B → P ν transitions was made by incorporating the resummation of real soft-photon emission and virtual corrections [42].Other work [43,44] follows an approach similar to the ones adopted in analogous kaon semileptonic decays [47][48][49][50][51][52][53][54][55].An important limitation in the attempts to compute those radiative corrections has to do with the poor knowledge of the photon-meson interactions at wavelength in the transition between the long-and short-distance regimes.Thus, most of these calculations, as well as the present one, suffer from limitations related to the used approximations, which may be well justified for light mesons and baryons but not necessarily for B mesons.
The O(α) electromagnetic radiative corrections involve the emission and/or re-absorption of virtual and real photons, with virtual-photon energies ranging from zero to infinity.By long-distance (LD) corrections we mean, emission/absorption of real and virtual photons with small momenta k (< 1 ) such that they cannot resolve the structure of hadrons; at these low photon momenta, the pointlike approximation for hadrons and the use of scalar QED for describing photon-meson interactions should be a good approximation.For photon momenta larger than a scale 2 , photons resolve the charged components of mesons and the description in terms of the electroweak (EW) theory to describe short-distance (SD) interactions is justified.
The transition region 1 (afew − hundreds MeV s) < k < 2 (a few GeV) between the short-and long-distance approximations is more difficult to evaluate and a description based on exchange of meson resonances may be an appropriate model.To the best of our knowledge, an extension of the long-distance approximation to cover photon energies around 1 GeV has been done only in Ref. [56] by including resonance degrees of freedom in τ → Mν(γ ) and M → ν (γ ) decays.Reference [56] has found that the long-distance corrections to the ratio R τ/π ≡ (τ → πν)/ (π → μν) computed in the scalar QED approximation using a scale μ cut to cut the integration over virtual photons differs from the improved calculation that include resonances by less than 0.4 %.Also, it was found [56] that the dependence upon the scale μ cut (assumed to separate short-and long-distance regions) largely cancel in the ratio R τ/π and, actually, it becomes insignificant for μ cut > 1.5 GeV's.However, the radiative corrections to the individual decay rates keeps a substantial dependence upon μ cut [56].
In practice a single scale (which plays the same role as μ cut ) can be used to separate the long-and short-distance regimes, and further assumptions as regards a small variation of hadronic structure have to be adopted.Thus, in the absence of a widely accepted prescription that perfectly matches the corrections in these two regions, an unavoidable dependence of radiative corrections δ T upon the cutoff scale naturally arises. 4As long as the radiative correction δ T = δ LD ( ) + δ SD ( ) does not strongly depend on the cutoff scale, we may be more confident upon the model assumptions used to compute the long-distance corrections.In the present paper we will adhere to this procedure in the calculation of the integrated rates and assign an uncertainty due to the dependence owing to un-matching of long-and shortdistance radiative corrections.Since our results will be suitable for photon-inclusive rates, in order to extend the validity of real-photon corrections to higher photon momenta, we also consider some of their model-dependent terms.Certainly, this and the difficulty to assess the uncertainty of other possible structure-dependent contributions is one of the limitations of our approach.
The integrated rate of B − → V 0 − ν decays (V = ρ, ω, D * ) including corrections of O(α) will be written as In the first line of (2), the superscripts in the rates denotes the order in α and the subscripts the region of phase space corresponding to three-and four-body kinematics (the four-body region is accessible only when real-photons are emitted).Up to O(α), the LD and SD corrections (superscripts denote the order in α) can be factorized according to the second line in Eq. ( 2).
The SD corrections to semileptonic decays are finite in the electroweak theory [57][58][59].In the dominant logarithmic approximation the O(α) corrections depend upon the separation scale as follows (m Z is the Z boson mass): It is customary to choose the lower cutoff of SD corrections as the mass of the decaying particle; however, we will keep it explicitly to study the cutoff scale dependence of the full radiative corrections.The large EW logarithms can be summed to all orders using the renormalization group [60]; the dominant part of O(αα s ) corrections have also been calculated [57,58] and resummed to all orders [61,62].For consistency, here we keep only the leading approximation (3) in our evaluation.
In order to compute the long-distance part of the radiative corrections, let us consider for definiteness the decay , where V is a vector meson; the first character within parentheses denote the four-momentum of each particle and the second its corresponding mass.At the tree level, the decay amplitude is where G F is the Fermi constant, V qb (q = u, c) the relevant CKM matrix element, L ν = ū γ ν (1 − γ 5 )v ν represents the leptonic current, and c V = 1/ √ 2 is the Clebsch-Gordan coefficient for V = (ρ 0 , ω) mesons.The hadronic matrix element can be parametrized in terms of four, q 2 = (P − P V ) 2dependent form factors, 5 which we choose as The four-vector ϕ denotes the vector-meson polarization, with the orthogonality condition P V • ϕ = 0.The form factors are specific to the B → V transition, although hereof Correspond to the self-energies of charged particles and c, d to vertex contributions we will disregard their superscripts.At higher orders, the factorization of the amplitude is not valid and the effective form factors may depend upon another invariant variable, for example u = (P − p) 2 .

Virtual-photon corrections
The virtual QED corrections of O(α) are shown in Fig. 1.
We will consider that the momenta of virtual photons are not very large (typically < 1.5 GeV) so that scalar QED can be used for the electromagnetic vertices of the B − meson and the weak hadronic vertex does not vary strongly.The charged meson (δ Z B ) and lepton (δ Z ) self-energies (Fig. 1a, b) contribute to the wave function renormalization of the charged particles.In terms of M 0 their contributions to the decay amplitude are where Here, are Passarino-Veltman functions corresponding to the scalar two-point integral and its derivative [63,64].Because of the infrared (IR) singularity in B (M) 0 , we have provided a fictitious mass λ to the photon.The ultraviolet (UV) singularity contained in B 0 can be regulated by using the cutoff as the maximum scale at which the long-distance approximation is expected to be valid; the expressions for the Passarino-Veltman functions required in this work are shown in Appendix A.
In order to treat the remaining corrections shown in Fig. 1c, d we follow Sirlin's prescription [65] to separate their modelindependent (basically scalar QED) and model-dependent contributions.This leaves us with the amplitude where D μν (k), (P + k) and S F ( p + k) correspond to the propagators of the photon (chosen in this work in the Feynman gauge), the pseudoscalar meson and the lepton, respectively.Thus, the model-independent part of the virtual correction corresponding to the photon exchange between the initial meson and the charged lepton (Fig. 1c) gives rise to the first term in square brackets of Eq. ( 7).The modeldependent part is encoded in T μλ , where gauge invariance requires k μ T μλ = 0, which can be verified by applying the generalized Ward identity.This model-dependent contribution can be included in a re-definition of the hadronic form factors as done for instance in [52][53][54][55]65] or evaluated explicitly as in Ref. [56] in a given model.For the model-independent part of the virtual corrections we have (remember u = (P − p) 2 ) where C i [m 2 , M 2 , u, m 2 , λ 2 , M 2 ] (i = 0, 1, 2) denote the three-point functions.We have used again λ to isolate the IR divergence contained in C 0 , while C 1,2 can be evaluated in the limit λ → 0, and be written in terms of -dependent scalar two-point functions as is shown in Appendix A. The last term in Eq. ( 8), where we have defined mV ub W μ (P V , P) ū Pγ μ (1 − γ 5 )v ν , corresponds to a non-factorizable (NF) amplitude, it is IR safe and gives a negligible contribution for light charged leptons in the final state.
The regulated infrared divergences can be isolated and written explicitly, while the UV divergences are implicit in the Passarino-Veltman functions.By collecting the results of the virtual corrections, Eqs. ( 6)-( 8), we can write the radiatively corrected rate as follows: where β(E) is the velocity (energy) of the charged lepton and E V the energy of the vector meson in the rest frame of the decaying particle.The expression for the function F 2 (E) can be found in Appendix A. In the above result, d 0 /dEdE V corresponds to the differential decay rate of B − → V 0 − ν at the tree level, while d 1 N F /dEdE V arises from the interference of the tree-level amplitude and the last term in Eq. ( 8).

Real-photon corrections
In this subsection we consider the real-photon corrections.As is well known, the cancellation of IR divergences occurs when we add incoherently the decay rates of a process with virtual and real-photon corrections.We will first consider the soft-photon approximation (SPA) for the decay amplitude and assume scalar QED for the electromagnetic vertex of B − meson.This model-independent approximation may be questionable and in principle it would be more appropriate only for the end-point region of the lepton spectrum [42].Then we will also include model-dependent terms of order zero and one in the photon momentum k.The general form of the decay amplitude can be written as M γ = M Low +M MD , where M Low is the Low soft-photon amplitude [66,67] which contains terms of order k −1 and k 0 , only; the second term M M D is explicitly model dependent and of O(k).In the case of K 3 decays the model-dependent terms of O(k 0 ) were found to be negligible 6 ; however, they deserve a proper study in B decays.

Model-independent corrections
By taking into account the contributions in Fig. 2 one gets the decay amplitude in the Low approximation [66,68] 6 The expression of those terms depends on the specific model for the weak form factors; in the case of K 3 decays their effects in the radiative corrections are at the one per mille level [51].

Fig. 2 Feynman diagrams for real-photon emission
where ε denotes the polarization four-vector of the real photon, with k•ε = 0. We have defined D λ = (ε• P/k• P)k λ −ε λ , L and W are the same as in Eq. ( 4), and As required by Low's theorem [66,67], Eq. ( 10) contains terms up to O(k 0 ) which are fixed by gauge invariance and depend only upon the form factors of the non-radiative decay amplitude.The partial derivative in the last term of Eq. ( 10) must be taken over explicit q 2 -dependent form factors in W .For the purposes of numerical calculations, we will consider the evaluation of the last two terms in Eq. ( 10) as a part of the model-dependent contribution (see the next subsection).
In the remainder of this subsection, we will focus only on the first term of Eq. ( 10) which we call the model-independent piece of real-photon corrections.
After integrating over the momentum of the massless particles, we write the differential decay width corresponding to the model-independent term (MI) as where |M 0 | 2 is the unpolarized squared amplitude at lowest order, and (E, E V ) are the lepton and vector-meson energies in the rest frame of the decaying particle ( Note, however, that, in radiative decays, the phase-space region accessible to (E, E V ) is larger than in the three-body decays and will depend upon the energies carried out by the massless particles through the invariant mass variable x ≡ ( p + k) 2 [47][48][49].
In a similar way as in Refs.[47][48][49], the first term in Eq. ( 12) contains the contribution proportional to the zeroth order decay width and carries all the infrared divergences of realphoton corrections, regulated as in the previous case by λ, 2 and the region of (E, E V ) accessible in the four-body decay can be found in [47][48][49].The integrals (P f ≡ P (14) required for K 3 decays, were first evaluated in [47][48][49] and verified in [50].The coefficients C m,n of the additional terms that appear in Eq. ( 12) can be found in Appendix B. The second term in Eq. ( 12) is infrared finite and can be evaluated by setting the photon mass to zero.
The region of integration of the variable x is identical to the one considered in [47][48][49], after tanking into account the corresponding changes in the hadron masses.The domain of the Dalitz plot (E, E ρ ) is shown in Fig. 4 for the B − → ρ 0 μ − ν μ channel (hereafter we will refer to this channel; under the corresponding approximations, the results can be applied to channels with any neutral vector meson in the final state, as well).The region denoted R III is accessible to the muon and the ρ meson of the three-and four-body decays, while R IV−III is allowed only for the corresponding radiative decay.The evaluation of the integrals in Eq. ( 12) is done by considering appropriately the region of phase space corresponding to three-and four-body regions of the Dalitz plot.

Model-dependent corrections
In addition to the amplitude M Low for real-photon emission obtained by attaching the photon in all possible ways, one can have also model-dependent terms which are of O(k) and higher and are contained in M M D .For relatively low photon energies, it has been pointed out [69] that the diagrams shown in Fig. 3 must be taken into account in B → D ν γ decays since on-shell D * (→ Dγ ) intermediate states can have an important contribution as non-radiative events.For the decay channels considered in this paper 3 Model-dependent real-photon corrections; R 0 denotes a pseudoscalar resonance there are no nearby narrow intermediate states decaying into V 0 γ states that can produce a similar important effect.
For hard photons, the structure of electromagnetic vertices in terms of elementary components becomes important and an approach based on QCD degrees of freedom like the Soft Collinear Effective Theory [45,46] looks more appropriate.Our evaluation of the model-dependent contributions to M M D for moderately energetic photons follows closely Ref. [69].
The general form of the MD decay amplitude becomes From Fig. 3 we get [69] where P 1 = P −k and P 2 = P V +k denote the four momenta of the corresponding resonances, and m B * (m R ) the mass of B * − (R 0 ) (if the R resonance is produced on-shell, we have to use → m R R ).In the above expression J μ denotes the electromagnetic current and (V ν − A ν ) denotes the usual V-A weak currents [69].We define the transition matrix elements of the electromagnetic current as with g PV γ denoting the V Pγ coupling, and η the polarization four vector of B * .As expected, the MD amplitudes in Eq. ( 15) are of O(k); as noticed in Ref. [69], when the intermediate resonances become on-shell the corresponding decay amplitude in ( 16) become of O(k 0 ).In our case, we do not expect MD contributions with on-shell resonances, although it is interesting to evaluate the effect given that in the softphoton limit some contributions behave as (m 2 V − m 2 P ) −1 .The B → R weak matrix element is defined as usual with F 1,0 being the vector and scalar form factors, respectively.We have neglected the k dependence in evaluating the matrix element, since this would induce O(k 2 ) terms in the MD amplitude, which are not considered in our approximation.
Since the main decay of the B * meson is electromagnetic (B * → Bγ ), there are not available calculations of its weak B * → V matrix element.Thus, we have to rely on the Heavy Quark Symmetry [70,71] to estimate the B * − → V 0 weak transition.This symmetry allows to relate this matrix element to the B − → V 0 matrix element, in the limit of an infinitely heavy b quark.Using the algebra of heavy quarks operators in the HQET [72], one can derive: where v = p B /M is the four-velocity of the heavy meson and η the polarization four-vector of the V 0 vector meson.Since we are interested in getting an estimate of this modeldependent contribution to the real-photon radiative correction of O(α), this approximation should be good enough.
In the following subsection, we define the different contributions to the radiative corrections.As noticed before, we call the model-dependent (MD) part of the radiative amplitude to the sum of the last two terms in Eq. ( 10) and the amplitude defined in Eqs. ( 15)-( 16).

Order α radiative corrections
By including the effects of LD virtual-and real-photon corrections that were obtained in Eq. ( 9) and section 2.2, we can write the following expression for the photon-inclusive Dalitz plot distribution: The first term in the above expression, δ1 L D (E, E ρ , 2 ), contains the effects of virtual photons and bremsstrahlung contributions of (Eq.12) that should be integrated in the phasespace region R III of (E, E V ) compatible with the kinematics of non-radiative (three-body) decays.This piece of radiative corrections is infrared finite, model independent [65], but depend upon the UV cutoff scale that determines the range of validity of the LD approximation, as mentioned at the beginning of this section.
The second term in Eq. ( 20) contains the contribution of the MI terms of the real-photon corrections of Eq. ( 12) eval- The region R III is allowed to three-and four-body (radiative) decays, while R IV−III is accessible only in the radiative mode uated in the R IV−III region of phase space (see Fig. 4) [47][48][49]73,74].The last two terms in Eq. ( 20) arise from the interference of model-independent and model-dependent (MI-MD) contributions, and the square of the model-dependent (M D) contributions, respectively, integrated over the full kinematical range.Let us note that we have included Eq. ( 15) and the last two terms of Eq. ( 10) as the components of the MD amplitude.
On the other hand, the short-distance corrections (photon momenta above ) can be calculated using the electroweak theory in terms of elementary field components.At O(α), the dominant logarithmic term of this correction to the decay rate is given by Eq. ( 3) [57,58], where the mass of the Z boson serves as a natural UV cutoff in the EW theory.As is well known, the cutoff dependence does not cancel completely in the sum of δ1 L D and δ 1 S D corrections.Actually, the -dependent term remaining in the sum of LD and SD corrections, δ 1 ( ) = α ln(m Z / )/(2π) changes from 0.63 to 0.33 % when varies from 0.4 GeV up to the B meson mass scale.Note that the dependence can be eliminated, as done for instance, in the calculation of radiative corrections to K 3 [50] and semileptonic tau decays [75,76] within the framework of chiral perturbation theory; there, a local counterterm can be chosen and added to the LD contributions in order to achieve a complete matching between LD and SD corrections.Within the approach used in the present paper, we will consider the effect of this cutoff dependence upon the observable decay rate as an uncertainty of the calculation.

Results
In this section we present the results of radiative corrections to the rates of B − → V 0 − ν decays, where V 0 is a neutral vector meson and a charged lepton.The integrals involved in virtual-and real-photon contributions were evaluated numerically.

Corrections to B
In this case we use the form factors of the B → V transitions as provided by lattice [77] and light-cone QCD sum rules [16][17][18] calculations, in the cases of V = ρ 0 , ω mesons.In the full range of q 2 values, the LD radiative corrections turned out to be rather insensitive to these two different models for the form factors; we will assign an uncertainty due to the numerical evaluation and to account for the different yields obtained by using the above parametrizations for the form factors.For the input required in model-dependent terms, we use the form factors F B→π 1,0 (q 2 ) given in [78,79], and the following values of coupling constants for V Pγ vertices: Table 2 shows the results for the short-distance (SD), longdistance (LD), and the total radiative corrections δ 1 T ( ), at three different values of the scale , for the B − → ρ 0 − ν decay rates.From this table we quote the following radiative corrections to the total decay: We have taken as central value for δ 1 T ( ) the result corresponding to the cutoff scale = m ρ .The first (larger) error is obtained by varying in the range shown in Table 2, while the second one is associated to the use of two different parametrizations of the form factors [16][17][18]77] and the numerical evaluation.In the absence of a rigorous prescription to compute radiative corrections in the intermediate energy regime that gives a perfect matching between scalar QED and the Standard Electroweak theory, we have to deal with this rather small uncertainty.
A few comments are in order: 1.As can be observed from Table 2, the LD corrections turn out to be very small, at the few per mille level.This is the case because we are including the effects of radiative corrections coming from radiative events falling in the four-body region, which almost cancel the LD contributions coming from the three-body region.Actually, events from the R IV−III region, see Fig. 4, contribute δ 1 L D (R IV−III ) = +0.0096(+0.0352) to the LD corrections in Table 2 of the rates into a muon (electron) channels.This trend is in agreement with similar conclusions gotten in Ref. [74] in the case of K − 3 decay.2. Inclusion of model-dependent contributions gives a correction ∼ −0.02 %, an order of magnitude smaller than the model-independent ones, for both leptons.The smallness of such a contribution is due to a partial cancellation between interference (M I − M D) and model-dependent (M D) corrections.On the one hand, the correction associated to M I − M D, ∼ −0.07 %, is dominated by the interference of the last two terms of the Low amplitude with the first line of Eq. ( 10).On the other hand, the MD correction, ∼ +0.05 %, is dominated by the MD amplitude in Eqs. ( 15)-( 16).Each of these corrections are found to contribute at the per mille level, as was noticed in K 3 decays [51,[53][54][55], but the relative sign between them makes the final result even smaller.3. Recently, bremsstrahlung simulations for specific channels in K 3 decays were implemented in PHOTOS [82,83].This novelty incorporates terms of O(k 0 ) in the decay amplitude which were found to be negligible compared to the contribution of the universal emission kernel [82,83].The approximation considered in Refs.[82,83] would correspond to neglecting the last term in Eq. ( 10) in the amplitude of B + → V 0 transitions considered in this paper.When we compare the results of LD corrections obtained using our full calculation with the one in PHO-TOS' approximation, we get a difference of +0.05 % for both leptonic channels.4. Given the equality of the matrix elements in Eq. ( 1) and the near equality of ρ 0 and ω meson masses, we expect the radiative corrections to B − → ω(782) − ν to be the same as in Eq. ( 21).An explicit numerical evaluation, based on the previous formulas, confirms this expectation.Thus, we conclude that the dominant source of isospin breaking to the second of the equalities in Eq.
(1) should be given by the effect of ρ-ω mixing [40,41], and not by radiative corrections.Note that the ratio of (B − → ρ 0 )/(B → ω) semileptonic rates is rather independent of the cutoff scale and form factor models of the B − → V 0 transition.For completeness, we provide also the evaluations of the LD radiative corrections to the Dalitz plot distribution.These corrections can be useful for experimental studies aiming a precise extraction of the weak form factors [47][48][49][52][53][54][55] which focus in the three-body region of the phase space.For this purpose, we re-write Eq. ( 20) as follows: with an obvious definition of the LD radiative correction function 1 L D (E, E ρ ) to the Dalitz plot.In Table 3 we show the values of the 1  L D (E, E ρ ) function for different energies of final states particles in the three-body region of the phase space.The large and negative values of radiative corrections for more relativistic charged leptons tend to dominate the contribution of the three-body phase-space region of radiative corrections.Despite the similar values of the total radiative corrections for both leptons, it is clear from Table 3 that the corrections can vary up to an order of magnitude from one leptonic channel to the other at the same point of the Dalitz plot.

Corrections to the
We can apply the same formalism to compute the LD QED radiative corrections to the decay B − →D * 0 − ν , by replac-ing the corresponding masses and hadronic form factors.For the numerical evaluation of virtual and model-independent real-photon corrections we use the form factor parametrization given in [84] and the results of fits obtained in Refs.[85,86].The relation between those parametrizations and the form factors in Eq. ( 5) are given in [87].In addition, for model-dependent real-photon contributions, we use the value g D 0 D * 0 γ = (1.0±0.9)×10−3 MeV −1 for the coupling constant extracted from D * 0 and D * ± decays [3].The form factors F B→D 1,0 are taken from Ref. [88].Since photons are not very hard in this case (E max γ ≈ 1.28 GeV), our approximations for the real-photon amplitude should be more reliable.
Table 4 summarizes our results for the radiative corrections to B − → D * 0 − ν , from which we can extract As in the previous case, the first uncertainty corresponds to variations of the (un)matching scale and the second one to the use of different form factor models and the numerical evaluation.Similar to the previous case, both corrections turn out to be equal for both leptonic channels.In this case, we chose (arbitrarily) the central value for the matching scale = m D * 0 , in analogy to Refs.[43,44] in the case of (pseudo)scalar charmed meson, but we have included the range of within m D * 0 /2 and 2m D * 0 .As in the case of charmless B − → V 0 transitions, a partial cancellation occurs between LD corrections evaluated in the three-and four-body regions of phase space.In this case, the LD effects from the R IV−III region of the Dalitz plot contribute +0.0051 and +0.0223 to the radiative corrections into muon and electron channels in Table 4, respectively.It is worth to mention that such cancellation of LD corrections has been noticed previously in radiative corrections to K − → π 0 − ν decays [74].The interference term M I − M D gives a correction ∼ −0.10 %, while the M D term ∼ +0.06 %, with the final result ∼ −0.04 %, for the two leptonic channels; the origin of such values is explained in the previous subsection.As in the case of B + → (ρ 0 , ω) transitions, the difference between our full calculation of the LD corrections and the one obtained in the PHOTOS' approximation [82,83] is +0.05 % in the integrated rate.
In Table 5 we show the contribution of LD radiative correction function 1  L D (E, E D * ), defined in Eq. ( 22), to the Dalitz plot distributions of B − → D * 0 − ν decays.Similarly, the predominance of radiative corrections of negative sign in most of the three-body region of phase space turns out to give radiative correction of negative sign for the LD corrections contributions to the decay rate coming from the R III region of the Dalitz plot.

Conclusions
We have calculated the long-distance (LD) corrections to the charmless and charmfull semileptonic B − →V 0 − ν decays involving the ground state vector mesons V 0 and = μ and e flavored leptons.We have used a formalism where the LD corrections are approximated by scalar QED, and assumed valid up to a scale , which serves as a separation cutoff of LD and short-distance (SD) corrections.Under these assumptions, the cutoff dependent corrections of O(α) depend mildly upon the scale .We have associated an uncertainty to radiative corrections in the decay rate due to this cutoff dependence by allowing it to vary between m V /2 and approximately 2m V .Also, in agreement with a similar analysis for K 3 decays [51], the effects of using different models for the weak form factors are estimated to play a subleading role compared to the dependence upon the separation scale .
The use of scalar QED in the hadronic vertices may be justified for relatively low values of the cutoff in the photon momenta.In the case of real-photon corrections we have included, in addition to the Low soft-photon amplitude, some model-dependent contributions that originate in the exchange of meson resonances.This allows one to cover harder real photons that are suitable to describe photon-inclusive rates.These model-dependent contributions, are found to be ten times smaller than the one obtained from the soft-photon approximation.
We have found that there exist large cancellations between the contributions of LD corrections to the decay rates coming from the three-and four-body regions of the Dalitz plot.These cancellations occurs for the integrated rates of photoninclusive B − → V 0 transitions.This happens for both charmless and charmfull transitions, and for both flavors of leptons.This is in agreement with previous calculations of radiative corrections to K ± 3 decay rates [74].Our calculation of virtual corrections assumes a single separation scale of LD and SD corrections long-distance, which avoids the consideration of structure-dependent effects of photon-hadron interactions in the region populated by resonances.Although the calculation of intermediate distance effects deserves further studies (beyond the scope of this paper), currently we miss a precise prescription that matches the calculations based on the fundamental and the effective low-energy theories used to compute, respectively, the SD and LD radiative corrections.At present, any calculation can be based at most on a model-dependent approach [56].
Our results can be useful for future/improved measurements of the different charged decay channels of B → V − ν branching ratios in order to test consistency of data with isospin symmetry.To complete this test, the calculation of radiative corrections to B0 → V + − ν (and remaining isospin breaking corrections in form factors) will be necessary, which we plan as a future project since the electromagnetic vertex of the charged vector meson requires a careful treatment.Actually, the determination of the V ub and V cb matrix elements at the few percent level, requires the consideration of the full radiative corrections of O(α).

Fig. 4
Fig.4 Phase-space region accessible to B − → ρ 0 μ − ν (γ ) decays.The region R III is allowed to three-and four-body (radiative) decays, while R IV−III is accessible only in the radiative mode

Table 2
Contribution of the short-distance (SD) and long-distance (LD) radiative corrections to B − → ρ 0 − ν ( = μ or e) at three different matching scales.The LD corrections include contributions from three-and four-body regions, which partially cancel each other

Table 4
Same description as in Table2for B − → D * 0 − ν decay rates