Radiative corrections of $O(\alpha)$ to $B^- \rightarrow V^0 \ell^- \bar{\nu}_{\ell}$ decays

The $O(\alpha)$ electromagnetic radiative corrections to the $B^- \rightarrow V^0 \ell^- \bar{\nu}_{\ell}$ ($V$ is a vector meson and $\ell$ 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 $\Lambda$ that separates the long- and 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 three- and four-body regions of phase space, the photon-inclusive 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.


I. INTRODUCTION
A precise determination of the |V ub | quark mixing at a few of percent level, is crucial for future tests of the Standard Model (SM) picture of CP violation [1] 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 and consistent determinations. As an example, let us remind 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 [2]. 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 inputs is B → πℓν ℓ , yielding |V ub | excl = (3.28 ± 0.29) × 10 −3 [2]. This relies on the most precise measurements of the total and differential rates [3][4][5], as well as on the lattice [6] and the Light-cone QCD sum rules (LCSR) [7] calculations 1 of the form factor. * stostado@fis.cinvestav.mx † glopez@fis.cinvestav.mx 1 New LCSR calculations of the B → π form factors may actually 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 [2], with errors added in quadrature. Improvement in both, experiment and theory, will play a crucial role in understanding and solving this discrepancy 2 . On the other hand, as long as new measurements of B + → τ + ν τ [10] are getting in better agreement with the SM prediction [11], 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 [12], combined with the exclusive value of |V cb | [2] and Lattice QCD calculations of the relevant hadronic form factors [13], have provided the first competitive determination of |V ub | from baryon semileptonic decays: |V ub | Λ b = (3.27 ± 0.15 exp ± 0.16 Latt ± 0.06 V cb ) × 10 −3 [12].
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 I). Here we will assume the narrow width approximation for vector mesons; in practice, vector mesons are unstable particles that further decrease the extracted value of |V ub | by up to 10% [8]. 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 [9]. 3 Similar relations hold by replacing ρ → π mesons. are reconstructed from the invariant-mass distribution of suitable decay channels [14,15]. 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 modelling 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.  [16], BELLE [17] and CLEO [18], and the lifetimes [2] of B mesons as inputs.

Channel
PDG 2014 [2] 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 [14,19] and long-distance electromagnetic radiative corrections. Future improved measurements of semileptonic B → (π, ρ, ω) transitions will require 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 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 first two decay rates in Eq. (1), but not the ratio between the last 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 [19].
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 softphoton 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.
There are only a few works related to the calculation of radiative corrections to semileptonic decays of B mesons [20][21][22]. 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 [20]. Other works [21] follow an approach similar to the ones adopted in analogous kaon semileptonic decays [23][24][25][26][27]. 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 point-like 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 (a few − 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 [28] by including resonance degrees of freedom in τ → M ν(γ) and M → ℓν ℓ (γ) decays. Ref [28] 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 [28] 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 [28].
In practice a single scale Λ (which plays the same role as µ cut ) can be used to separate the long-and shortdistance regimes, and further assumptions about 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 4 . As 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 longdistance 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 short-distance 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 realphotons 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 [29,30]. 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 [31]; the dominant part of O(αα s ) corrections have also been calculated [29] and resummed to all orders [32]. 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 B − (P, M ) → V 0 (P V , m V )ℓ − (p, m)ν ℓ (p ′ , 0), 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 Clebsh-Gordan coefficient for V = (ρ 0 , ω) mesons. The hadronic matrix element can be parametrized in terms of four, 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 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 .

A. 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 ℓ ) selfenergies ( Fig. 1 a and b) contribute to the wave function renormalization of the charged particles. In terms of M 0 their contributions to the decay amplitude are where [· · ·] are Passarino-Veltman functions corresponding to the scalar two-point integral and its derivative [33]. 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 Figure 1c,d we follow Sirlin's prescription [34] to separate their model-independent (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. 1 c) gives rise to the first term in square brackets of Eq. (7). The model-dependent 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 [26,27,34] or evaluated explicitly as in Ref. [28] 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 the Appendix A. The last term in Eq.
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 the 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 treelevel amplitude and the last term in Eq. (8).

B. 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 [20]. Then, we will also include model-dependent terms of order zero and

FIG. 2. Feynman diagrams for real photon emission.
one in the photon momentum k. The general form of the decay amplitude can be written as where M Low is the Low's soft-photon amplitude [35] which contains terms of order k −1 and k 0 , only; the second term M MD 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's approximation [35,36] 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 [35], Eq. (10) contains terms up to O(k 0 ) which are fixed by gauge-invariance and depend only upon the form factors of the nonradiative decay amplitude. The partial derivative in the last term of Eq. (10) must be taken over explicit q 2dependent 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 modeldependent contribution (see 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 vectormeson energies in the rest frame of the decaying par- . Note however that, in radiative decays, the phase space region accessible to (E, E V ) is larger than in the threebody decays and will depend upon the energies carried out by the massless particles through the invariant mass variable x ≡ (q + k) 2 [23].
In a similar way as in Ref. [23], the first term in Eq. (12) contains the contribution proportional to the zeroth order decay width and carries all the infrared divergences of real-photon corrections, regulated as in the previous case by λ, Here, x max = (M − E − E V ) 2 − |p V − p| 2 and the region of (E, E V ) accessible in the four-body decay can be found in [23]. The integrals (P f ≡ P V + p + p ′ + k) required for K ℓ3 decays, were first evaluated in [23] and verified in [24]. The coefficients C m,n of the additional terms that appear in Eq. (12) can be found in the 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 [23], after tanking into account the corresponding changes in the hadron masses. The domain of the Dalitz plot (E, E ρ ) is shown in Figure 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 real photon corrections; R 0 denotes a pseudoscalar resonance.

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 MD . For relatively low photon energies, it has been pointed out [37] that the diagrams shown in Figure 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 B ± → (ρ 0 , ω, D * 0 )ℓ ± ν ℓ , there are not 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 [22] looks more appropriate. Our evaluation of the model-dependent contributions to M MD for moderately energetic photons follows closely Ref. [37].
The general form of the MD decay amplitude becomes From Figure 3 we get [37] V 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 onshell, 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 [37].
We define the transition matrix elements of the electromagnetic current as 2 k η , with g P V γ 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 it was noticed in Ref. [37], when the intermediate resonances become onshell 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 soft-photon 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 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 [38] 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 [39], one can derive: where v = P B /m B 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 model-dependent 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).

C. Order α radiative corrections
By including the effects of LD virtual and real photon corrections that were obtained in Eq. (9) and section 2, we can write the following expression for the photoninclusive Dalitz plot distribution The first term in the above expression,δ 1 LD (E, E ρ , Λ 2 ), contains the effects of virtual photons and bremsstrahlung contributions of (Eq. 12) that should be integrated in the phase-space region R III of (E, E V ) compatible with the kinematics of non-radiative (threebody) decays. This piece of radiative corrections is infrared finite, model-independent [34], but depend upon the UV cutoff scale Λ that determines the range of validity of the LD approximation, as it was mentioned at the beginning of this section. 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) [29], where the mass of the Z boson serves as a natural UV cutoff in the EW theory. As is well known, the cutoffdependence does not cancel completely in the sum ofδ 1 LD and δ 1 SD 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 [24] and semileptonic tau decays [42] 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.

III. 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. The region R III is allowed to three-and four-body (radiative) decays, while R IV −III is accesible only in the radiative mode.
A. Corrections to B − → (ρ 0 , ω)ℓ −ν ℓ decays In this case we use the form factors of the B → V transitions as provided by lattice [43] and light-cone QCD sum rules [7] 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 [44], and the following values of coupling constants for V P γ vertices: g B − B * − γ = (1.0 ± 0.1) × 10 −3 MeV −1 [45] and g πργ = (7.3 ± 0.8) × 10 −4 MeV −1 [46]. Table II shows the results for the short-distance (SD), long-distance (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: (1.62 ± 0.10 ± 0.04)% , for ℓ = µ (1.63 ± 0.11 ± 0.04)% , for ℓ = e .
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 II, while the second one is associated to the use of two different parametrizations of the form factors [7,43] 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 it can be observed from Table II, 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 II. Contribution of the short-distance (SD) and longdistance (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 II of the rates into a muon (electron) channels. This trend is in agreement with similar conclusions gotten in Ref. [41] in the case of K − ℓ3 decay.  (16). Each of these corrections are found to contribute at the per mille level, as was noticed in K ℓ3 decays [25,27], 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 PHO-TOS [47,48]. 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 [47,48]. The approximation considered in Refs. [47,48] 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 PHOTOS' 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 formulae, 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 [19], 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 [23,26,27] 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 LD (E, E ρ ) to the Dalitz plot. In Table III we show the values of the ∆ 1 LD (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 III that  We can apply the same formalism to compute the LD QED radiative corrections to the decay B − → D * 0 ℓ −ν ℓ , by replacing 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 [49] and the results of fits obtained in refs. [50,51]. The relation between those parametrizations and the form factors in Eq. (5) are given in [52]. In addition, for modeldependent 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 [2]. The form factors F B→D 1,0 are taken from Ref. [53]. 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.
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   Table II for 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 ref. [21] 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 IV, respectively. It is worth to mention that such cancellation of LD corrections has been noticed previously in radiative corrections to K − → π 0 ℓ −ν ℓ decays [41]. 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 [47,48] is +0.05% in the integrated rate.
In Table V we show the contribution of LD radiative correction function ∆ 1 LD (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.

IV. 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 [25], the effects of using different models for the weak form factors are estimated to play a sublead-  Table III  ing 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's softphoton amplitude, some model-dependent contributions that originate in the exchange of meson resonances. This allows 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 photon-inclusive 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 [41]. Our calculation of virtual corrections assumes a single separation scale Λ of LD and SD corrections longdistance, which avoids the consideration of structuredependent effects of photon-hadron interactions in the region populated by resonances. Althought 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 in a model-dependent approach [28].
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 toB 0 → 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(α).
where λ represents a ficticious photon mass, introduced to regulate the IR-divergence. Similar expressions to Eqs. (A1),(A3) can be obtained for the scalar integrals associated to the loops involving the charged meson by proper replacements of momenta and masses. The function C 0 can be split into two terms where β = |p|/E is the charged lepton velocity in the rest frame of the decaying particle. The first term is IR-divergent while the second is IR-finite and given by with χ(y) = uy 2 + 2(p · P − p 2 )y + p 2 . The three-point functions C 1,2 can be reduced to scalar two point integrals as In this Appendix we provide the expressions for the non-vanishing coefficients that enter the definition of the model-independent real-photon corrections in Eq. (12). They are different from the analogous coefficientes semileptonic kaon decays [23].
In the above expressions we have defined P 2 and have introduced the following reduced form factors: v = 2V /(M +m ρ ), a 1 = (M +m ρ )A 1 , a 2 = A 2 /(M +m ρ ) and a = 2m ρ A/q 2 .