Heavy + heavy and heavy + light pseudoscalar to vector semileptonic transitions

Using a symmetry-preserving regularisation of a vector$\times$vector contact interaction (SCI), we complete a systematic treatment of twelve semileptonic transitions with vector meson final states: $D\to \rho$, $D_{(s)}\to K^\ast$, $D_s\to \phi$, $B\to \rho$, $B_s\to K^\ast$, $B_{(s)}\to D_{(s)}^\ast$, $B_c \to B_{(s)}^\ast, J/\psi, D^\ast$; and thereby finalise a unified analysis of semileptonic decays of heavy+heavy and heavy+light pseudoscalar mesons to both pseudoscalar and vector meson final states. The analysis is marked by algebraic simplicity, few parameters, and the ability to consistently describe systems from Nambu-Goldstone modes to heavy+heavy mesons. Regarding the behaviour of the transition form factors, the SCI results compare well wherever sound experimental or independent theory analyses are available; hence, the SCI branching fraction predictions should be a reasonable guide. Considering the ratios $R(D_{(s)}^{(\ast)})$, $R(J/\psi)$, $R(\eta_c)$, whose values are key tests of lepton universality in weak interactions, the SCI values agree with Standard Model predictions. The $B_{(s)}\to D_{(s)}^\ast$ transitions are used to predict the precursor functions that evolve into the universal Isgur-Wise function in the heavy-quark limit, with results that conform with those from other sources where such are available. The study also exposes effects on the transition form factors that flow from interference between emergent hadron mass from the strong interaction and Higgs boson couplings via current-quark masses, including flavour symmetry violation.

( * ) (s) ), R(J/ψ), R(η c ), whose values are key tests of lepton universality in weak interactions, the SCI values agree with Standard Model predictions. The B (s) → D * (s) transitions are used to predict the precursor functions that evolve into the universal Isgur-Wise function in the heavy-quark limit, with results that conform with those from other sources where such are available. The study also exposes effects on the transition form factors that flow from interference between emergent hadron mass from the strong interaction and Higgs boson couplings via current-quark masses, including flavour symmetry violation.

Introduction
Nature has two known mechanisms for mass generation. In connection with quantum chromodynamics (QCD), that associated with the Higgs boson (HB) is responsible for the current-quark masses, which range from a few MeV for the lightest quarks to nearly 200 GeV for the top quark. Within the Standard Model of particle physics (SM), each current mass is generated by a distinct HB coupling; so, one has a parametric representation but not a satisfying explanation for this hierarchy of scales [1]. The other source of mass appears to be a dynamical feature of QCD; namely, emergent hadron mass (EHM), which is thought to be responsible for an array of phenomena that include [2][3][4][5] the generation of nuclear-size masses for baryons constituted from light quarks and nearly-massless pseudoscalar Nambu-Goldstone bosons whose existence is crucial to the stability of known nuclei.
Pseudoscalar mesons are special because they are massless in the absence of HB couplings. Consequently, weak-interaction mixing between quark flavours, parametrised using the Cabibbo-Kobayashi-Maskawa (CKM) matrix, links both Nature's sources of mass via the study of semileptonic decays of heavy+heavy and heavy+light pseudoscalar mesons. In such transitions, the mass of the initial state owes to HB couplings and constructive EHM plus HB interference, with the following HB:EHM+HB mass budgets [6]: B c (87 : 13), B (s) (80 : 20), D (s) (70 : 30). When the final states are also pseudoscalar mesons, these mass budgets can be dramatically reversed, e.g. [7]: K(20 : 80), π(5 : 95). In contrast, vector meson mass budgets are more like those of baryons, being nonzero even in the absence of HB couplings, with EHM alone delivering the bulk of their masses in lighter-quark cases, e.g., ρ(97 : 3) [8].
Thus, whilst semileptonic pseudoscalar-to-pseudoscalar decays have received most attention, there are also good reasons to consider pseudoscalar-to-vector transitions that go beyond their use as independent constraints on the CKM matrix elements.
The typically large disparity in masses between initial and final states is a significant challenge in the study of heavy pseudoscalar meson semileptonic transitions. No framework with a traceable link to QCD can directly surmount this difficulty today. Nevertheless, using a variety of methods, attempts are being made for pseudoscalar  and vector [9-13, 15, 16, 18, 20-28, 30, 32-35, 49, 50, 52-55] final states. Given the challenges involved, unified treatments of pseudoscalar and vector final states are of particular value.
Using a symmetry-preserving formulation of a vector × vector contact interaction (SCI) [56,57], a study of ten pseudoscalar-to-pseudoscalar semileptonic transitions was described in Ref. [6]: D → π, D (s) → K, B → π, B s → K, B (s) → D (s) , B c → B (s) , η c . Combined, these transitions provide information on four elements of the CKM matrix: |V cd |, |V cs |, |V ub |, |V cb |. Herein, we extend that analysis to a large array of analogous transitions with vector meson final states, all sensitive to the same set of CKM matrix elements, viz. the following twelve decays: D → ρ, D (s) → K * , D s → φ, B → ρ, B s → K * , B (s) → D * (s) , B c → B * (s) , J/ψ, D * . As Ref. [6] provided a benchmark for Ref. [48], which used continuum Schwinger function methods (CSMs) to deliver a coherent, parameter-free treatment of pseudoscalar-to-pseudoscalar transitions, so will the current study serve for the extension of Ref. [48] to vector final states. This is because contemporary implementations of the SCI preserve the essence of more sophisticated treatments of the continuum bound-state problem yet introduce an algebraic simplicity. Widespread use has revealed that, when interpreted judiciously, SCI predictions provide a valuable quantitative guide, see e.g. Refs. [58][59][60][61][62][63][64][65][66][67][68][69]. Thus, SCI results provide a means by which one may check the validity of algorithms employed in calculations that rely (heavily) upon high performance computing.
Section 2 presents some general material on the structure of currents and definitions of form factors in semileptonic pseudoscalar-to-vector meson transitions. Section 3, augmented by an extensive appendix, explains our implementation of the SCI, including constraint of the ultraviolet cutoff and determination of values for the interaction-dependent current-quark masses. Section 4 reports our analysis of s) , J/ψ, D * , including a discussion of the Isgur-Wise function. Section 6 considers the evolution of semileptonic transition form factors as the Higgs mechanism of current-quark mass generation becomes a more significant component of the final-state meson's mass, viz. it discusses the issues of environment sensitivity and flavour symmetry breaking in electroweak transitions. Section 7 provides a summary and perspective.

Transition Form Factors
Semileptonic decays of pseudoscalar mesons (P) to vector meson (V) final states are described by the following transition matrix element: where (t − is the largest accessible value of t in the identified physical decay process.) The scalar functions in Eq. (1) are the semileptonic transition form factors, which express all effects of hadron structure on the transitions. Ensuring the absence of kinematic singularities in Eqs. (1), symmetries guarantee Once results for the transition form factors are available, one can calculate the associated decay branching fractions from the differential decay width for P → V + ν : where: m 2 ≤ t ≤ t − , m is the lepton mass; and It is plain from Eq. (4) that the contribution from A 0 (t) to any cross-section is kinematically suppressed. After integrating Eq. (3) to obtain the required partial widths, one quotes the branching fractions, B P →V ν , with respect to the total width determined from the pseudoscalar meson's lifetime [70].

Matrix Elements
When employing CSMs, the leading-order approximation to a matrix element like that in Eq. (1) is provided by the rainbow-ladder (RL) truncation [71], which we choose to illustrate using the D 0 → ρ − transition: where N c = 3 and the trace is over spinor indices. Three distinct types of matrix-valued function feature in Eq. (6): propagators for the dressed-quarks involved in the transition process, here S f (t), f = u, d, c; Bethe-Salpeter amplitudes for the initial-and final-state mesons, Γ M , M = D 0 , ρ − ; and the dressed cd weak transition vertex, W cd µ . The P → V weak transition vertex has two pieces, viz. vector and axial-vector: Physically, the vector part, V cd µ , must exhibit poles at Q 2 +m 2 D * ,D * 0 = 0 and A cd µ has poles at Q 2 +m 2 D,D1 = 0. The presence of such poles is a prerequisite for any valid analysis of P → V semileptonic transitions; and as will become apparent, they are manifest in our treatment.
The general structure of Eq. (6) is the same for any quark+antiquark interaction treated in RL truncation. However, in writing Eq. (6), we have implicitly assumed a SCI for all integral equations relevant to the problem. This is signalled by the simplicity of the arguments of the Bethe-Salpeter amplitudes and weak transition vertex, as elucidated in Appendix A.

Weak D (s) Semileptonic Transitions
In the isospin-symmetry limit, there are four essentially distinct such processes: D 0 → ρ − , D + s → K 0 * , D + →K * 0 , D s → φ. The first two measure c → d and the last two, c → s. Each provides information on the environmental sensitivity of these transitions.  On the physical domain, all D (s) transition form factors are monotonically increasing functions of t that can reliably be interpolated using where m P is the calculated mass of the initial-state pseudoscalar meson. Our predictions for each of the transitions considered in this subsection are described by Eq. (8) and the appropriate interpolation coefficients in Table 1A.
The SCI D → ρ transition form factors are depicted in Fig. 1A, wherein they are compared with monopole forms inferred from data [73,CLEO]. Given the simplicity of those fits, the agreement is good. Using the calculated form factors, Eq. (3) yields the fractions in Table 1B -Rows 1-2. Contrasted with model results in Ref. [33], the individual SCI branching fractions are 30% smaller. On the other hand, those models produce R µ/e = 0.95 (2), in agreement with the SCI prediction. Combin-ing SCI and model results [20,33,54], one obtains (1) ; (10) outcomes which suggest that the Ref. [70,PDG] value for B D + ρ 0 µ + νµ is too large.
Turning to D + s →K * , data on the e + ν e mode are available in Refs. [76,77,CLEO] and [78,BESIII]. SCI predictions for the form factors are depicted in Fig. 1B. They are reliably interpolated using Eq. (8) with the coefficients in Table 1A -Row 2. Little empirical information is available on the form factors; but Ref. [78,BESIII] reports r 2 = 0.77 (29), r V = 1.67 (38), values which agree with the SCI results in Table 1A -Row 2.
Our predictions for the branching fractions to the kinematically allowed semileptonic final states are listed in Table 1B -Row 2. No experimental results are available for the µ + ν µ mode; but the SCI values are commensurate with model estimates collected in Ref. [33]. Combining those results [20,33,54] with the SCI predictions, one obtains  (21), r V = 1.49(05). Whilst the value of r V agrees with the SCI result, that for r 2 is significantly larger. In this connection, the SCI predictions are compatible with the model results collected in Ref. [33]; and combined with those values, one obtains r 2 = 0.64 (24), r V = 1.44 (14). Apparently, with r 2 , there is some tension between experiment and theory. Table 1B -Rows 4 -5 list our predictions for the D 0 → K * − , D + →K * 0 branching fractions to the allowed semileptonic final states. They compare favourably with the averages reported in Ref. [70,PDG]. Combined with commensurate model estimates [20,54], one obtains Our D s → φ transition form factors are drawn in Fig. 2B. Reliable interpolations are obtained by using Eq. (8) with the coefficients in Table 1A -Row 4. Lattice QCD (lQCD) results for the form factors [53] are also drawn in Fig. 2B. They agree with our predictions. This is further highlighted by a comparison between values computed at the maximum recoil point: 1.08 lQCD [53] 0.74(12) 1.72 (21) 1.14 (6) .
SCI predictions for D s → φ transition form factors, with legend as in Panel A. Thinner, like-texture and -colour curves within shaded bands are lQCD results from Ref. [53].
states. They agree with the experimental results tabulated in Ref. [70,PDG]. However, there is only one measurement of the µ + ν µ mode and it has a large error. Given the agreement between SCI and lQCD determinations of the semileptonic form factors, there is merit in improving the precision of the D s → φ measurements.

Weak B (s,c) Semileptonic Transitions
In the isospin-symmetry limit, there are eight distinct processes, which we choose to be: ; and the last, c → s (iv ); plus, naturally, their environmental sensitivity.
All B (s,c) transition form factors are monotonically increasing functions of t on their respective physical domains. However, for decays involving a ρ or K * in the final state, those domains can be large; consequently, we found that a reliable interpolation of A 1 (t) needs a t 3 term in the denominator: with m P the calculated mass of the initial state pseudoscalar meson. Thus, whilst our predictions for almost all transition form factors considered in this subsection are well described by Eq. (8) and the appropriate interpolation coefficients in Table 1A, there are two cases where Eq. (14) and c = 0 are necessary and we note this explicitly in the associated discussion. Fig. 3A. Accurate interpolations are obtained using the coefficients in Table 1A - (14).
Results from a simulation of quenched lQCD using unphysically large light-quark masses are available [81]. Fig. 3A in order to illustrate both their precision and t-domain coverage. Nothing more recent is available from lQCD.
The SCI form factors in Fig. 3A yield the branching fractions to allowed semileptonic final states listed in Table 1B -Row 7. Where comparison is possible, the SCI results are roughly 50% larger than those quoted in Ref. [70,PDG]. Notably, SCI predictions for B 0 → π − exceed experiment by a factor of ∼ 2.8 [6, Table 3b] because the maximum recoil value of the vector transition form factor is too large [6, Table 3a]. We therefore compared the form factor maximum recoil values in Table 1A -Row 5 with the results in Refs. [10,15,18,23,52], finding that the SCI values are 39(9)% larger. This is sufficient to explain the overestimated branching fractions. Given the huge disparity in mass-scales between the initial and final states in the B 0 → ρ − transition, it is not too surprising that the SCI description is imperfect in this case. Nevertheless, it is markedly better than the analogous treatment of B 0 → π − ; hence, the form factors should be a semiquantitatively sound guide, given that those for B 0 → π − are broadly compatible with experiment [6, Fig. 5].
Results for B s → K * transition form factors are available from a simulation of unquenched lQCD, with four points on t/GeV 2 ∈ [14,19] at each of three pion masses, m π /GeV = 0.31, 0.34, 0.52, using lattice nonrelativistic QCD to describe the b-quark [82]. Those points were then fitted and extrapolated therein so as Panel B. SCI B s →K * transition form factor V (t) -dashed green curve. The thinner, like-textured and -coloured curve within the shaded band is a lQCD result from Ref. [82].
though the lQCD curves again fall below on t 8 GeV 2 .
Here it is worth recalling that all the lower-t lQCD results involve a long reaching extrapolation from a small number of points on t/GeV 2 ∈ [14,19].
Working with the SCI form factors in Figs. 3B, 3C and the appropriate form of Eq. (3), one finds the branching fractions listed in Table 1B SCI predictions for the B + c → D * 0 semileptonic transition form factors are drawn in Fig. 4. In this case, interpolations are obtained using Eq. (8) and the coefficients in Table 1A  Refs. [9, 11-13, 16, 50] yields B Bc→D * e + νe = 0.14(5) , viz. a central value that matches the SCI prediction. Including the SCI result in the average, one finds B Bc→D * e + νe = 0.14(4) .
Further, R τ /µ = 0.58 can be read from Ref. [22], a value which complements the SCI result. Notwithstanding these remarks, regarding B + c → D * 0 semileptonic transitions, it is plain that theory precision can improve, calculations with closer links to QCD are necessary, and the want of data is sorely felt.

Cabibbo inhibited b → c transitions: singly heavy
In this sector, we first consider B 0 → D * − semileptonic transitions. The SCI results, drawn in Fig. 5A, are accurately interpolated using Eq. (8) and the coefficients in Table 1A -Row 8. They yield the branching fractions in Table 1A -Row 10. Notably, whilst the τ ν τ fraction is commensurate with the current data average [70,PDG], the ν , = e, µ, fractions are significantly larger. Consequently, the value of R τ /µ -often denoted R(D * ) -is of special interest. The SCI result matches fairly with the value considered to be the SM prediction [83]: 0.252 (3). It therefore confirms a 2.9σ tension between theory and experiment on this ratio, which is a key test of lepton universality in Nature's weak interactions.
On the other hand, there are many model studies of B 0 (s) → D * − (s) semileptonic transitions and the results are widely scattered. For instance, considering an analysis that uses two related methods in a simultaneous treatment of both sets [34], one finds The SCI results fit well within these bands. Moreover, profiting from the kindred study of P → P transitions [6, Table 3], we arrive at the SCI prediction  Lattice results for the B 0 (s) → D * − (s) transition form factors are now available [87]. Curves developed from those results are drawn alongside the SCI predictions in Figs. 5B, 5C: within their uncertainties, the curves agree with the SCI predictions. There is, perhaps, some tension between the results for the vector form factor: the t = 0 value of the lQCD curve is unexpectedly large when compared with the many observableconstrained model calculations, e.g., Refs. [27,34]. Nevertheless, even in this case, the absolute relative difference between the SCI and lQCD curves is just 23 (14)

Isgur-Wise Function
It is common to discuss B 0 (s) → D * − (s) transitions in terms of form factors related to the Isgur-Wise function [88], defined as follows: where , with w = 1 corresponding to the zero recoil point. This is because, in the true heavy-quark limit, In the present cases, one cannot a priori be certain that the conditions defining the heavy-quark symmetry limit are realised because the c-quark mass is not necessarily very much greater than the corresponding u, d, squark masses, e.g., referring to Appendix A - Table 2, the s/c ratio of dressed-quark masses is ≈ 1/3. Consequently, the degree to which the heavy-quark limit is accessible in Nature is an empirical question.
Regarding the transitions discussed in Sec. 5.2, one finds the following zero recoil values of the functions in Eq. (20):  (13). They are commensurate with the SCI values. Evidently, the SCI predicts that with the c-quark involved in vector-meson final-states, there is a mean absolute relative deviation of 14(8)% from the heavy-quark limit result ξ(1) = 1.
In Fig. 6 we plot the precursor Isgur-Wise functions defined in Eq. (20), each divided by its zero recoil value. All drawn curves can reliably be described using   Here, "reliably" means the mean absolute relative difference between a given function and the fit is < 0.2%.

Cabibbo inhibited b → c transitions: doubly heavy
SCI results for the B c → J/ψ semileptonic transition form factors are drawn in Fig. 7, wherein they are compared with predictions obtained using a systematic, symmetry-preserving continuum approach to strong-interaction bound-state problems [49]. Notwithstanding the SCI's algebraic simplicity, its results agree well with  the sophisticated calculations: the mean error-weighted χ 2 values are 0.14, 0.27, 1.25, 1.15 for V , A 0 , A 2 , A 1 , respectively. Notably, too, the predictions in Ref. [49] and the lQCD results in Ref. [91] agree within mutual uncertainties; hence, the SCI curves also agree with those obtained using lQCD.
The SCI form factors in Fig. 7 are accurately interpolated using Eq. (8) with the coefficients in Table 1A -Row 10 and produce the branching fractions listed in Table 1B -Row 12. Unsurprisingly, the SCI predictions agree with those in Ref. [49]; hence, confirm the 2σ discrepancy between the SM prediction for R(J/ψ) and the measurement in Ref. [72,LHCb].
The kindred SCI analysis of B c → η c transitions [6] produces R(η c ) = 0.25, a value 7% larger than the SCI result for R(J/ψ) in Table 1B -Row 12. This size increase agrees semiquantitatively with that found in Ref. [49]. ing Eq. (8) with the coefficients in Table 1A -Row 11. In this instance, we have plotted V (t)/2 because the vector form factor is markedly enhanced compared with all axial form factors. In fact, on the physically accessible domain, V BcB * (t) is significantly larger than all other form factors considered above. This result is consistent with that found in Ref. [30].
The form factors in Fig. 8A produce the branching fractions listed in Table 1B -Row 13. No data are available. On the other hand, the SCI value for the eν e final state is consistent with the results in Refs. [16,25,30], the average of which is 1.32 (12) . Including the SCI value, the average becomes 1.32 (9) . Concerning the ratio of branching fractions, Ref. [30] reports 0.95, matching the SCI prediction.

Cabibbo favoured c → s transition
Our predictions for the B + c → B * 0 s transition form factors are drawn in Fig. 8B. The functions are accurately interpolated using Eq. (8) with the coefficients in Table 1A -Row 12. Here, for the same reasons as with B + c → B * 0 , we have drawn V (t)/2, whose magnitude and behaviour is similar to that found elsewhere [30].
Using the form factors in Fig. 8B, one obtains the branching fractions listed in Table 1B -Row 14. In the absence of data, we note that the SCI result for the eν e final state is consistent with the values in Refs. [16,25,30], which average to 19.8(2.9) . Including the SCI prediction in the average, one obtains 19.3(2.6) . Regarding the ratio of branching fractions, Ref. [30] reports 0.93, matching the SCI prediction.

Environment Sensitivity
Considering Nature's two mass generating mechanisms, it is worth considering the evolution of the form factors with increasing mass of the valence-quark that is a spectator to the transition; namely, as the currentmass contributed by the Higgs boson becomes a more significant part of the final-state meson's mass when compared with the EHM component. Restricted to the {u, d, s} quark sector, such effects express SU(3) flavour symmetry violation. Four classes can be defined for the semileptonic decays considered herein.
with spectator quark, respectively, u, s, b. These are analogues of the pseudoscalar-to-pseudoscalar tran-  Fig. 9B, the general trend is for the zero-recoil value of each form factor to decrease with increasing spectator mass. In the pseudoscalar-topseudoscalar case, the pattern in Fig. 9B can be understood by recalling that f + (t) is kindred to an elastic form factor when the current-quark masses in the initial-and final-states are similar. No link to an elastic form factor can be made in any other case drawn in Fig. 9, so the behaviour expresses novel features of EHM-HB interference in electroweak transitions. spectator flavour Restricting our considerations to cases with u/d, s spectator quarks, the median absolute relative difference between maximum recoil values is ard = 7.2% and the associated mean is ard = 6.5(4.4)%. Working with available information on leptonic decay constants - Table 2, the analogous results are ard = 19% and ard = 17(6)%, viz. 2.7-times larger.
Notably, the Ward-Green-Takahashi identities satisfied by the two pieces of the weak transition vertex -Ref. [6,Eq. (29)] and Eq. (A.21) -involve differences and sums of quark current-masses multiplied by vertices whose structure can be greatly affected by EHM. Consequently, one should anticipate complex interference effects. The SCI predictions can be checked, e.g., by expanding the array of transitions studied in Ref. [49] and extending that analysis to vector-meson final-states.

Summary and Perspective
Working with a symmetry-preserving regularisation of a vector×vector contact interaction (SCI), we presented a unified treatment of twelve independent pseudoscalarto-vector meson (P → V ) semileptonic transitions: , J/ψ, D * , each with four form factors; and the masses and leptonic decay constants of twentythree mesons that are either involved in the transitions or contribute in a material way to the weak transition vertex. In completing this analysis, we simultaneously unified P → V semileptonic transitions with ten kindred pseudoscalar-to-pseudoscalar (P → P ) transitions, which were calculated elsewhere [6] using precisely the same framework; and thereby finalised a comprehensive, coherent analysis of the semileptonic decays of heavy+heavy and heavy+light pseudoscalar mesons.
Our implementation of the SCI has four parameters, which are values of a mass-dependent quark+antiquark coupling strength chosen at the current-masses of the u/d, s, c, b quarks. The merits of the approach are its algebraic simplicity; paucity of parameters; and simultaneous applicability to a wide variety of systems and processes, sometimes involving large disparities in mass-scales between initial and final states.
Concerning meson masses, which are long-wavelength properties of the systems, the agreement between SCI predictions and experiment is good [ Table 2]. The comparison yields a median absolute relative difference ard = 2% and mean ard = 4%. Including the leptonic decay constants, ard = 3% and ard = 10%. Since leptonic decay constants are dominated by ultraviolet momenta, they are more of a challenge for the SCI.
Concerning the t dependence of the P → V transition form factors, wherever experiment or solid theory results are available for comparison, the SCI results compare well [Secs. 4,5]. The stiffness found with P → P transitions is largely avoided in P → V form factors because the SCI does not support a tensor term in the vector-meson Bethe-Salpeter amplitude when the quark+antiquark scattering kernel is treated in rainbowladder truncation. The least satisfactory comparisons are found with B → ρ and B s → K * ; but owing to the huge disparity between mass-scales of the initial and final states, these transitions present difficulties for all available methods. In every other case, our analysis suggests that the SCI branching fraction predictions should be a sound guide. Of special interest are the SCI results for the branching fraction ratios R(D ( * ) (s) ), R(J/ψ), R(η c ), whose values are key tests of lepton universality in Nature's weak interactions. In all cases, the SCI values agree with accepted Standard Model predictions.
Working with B (s) → D * (s) transitions, we provided predictions for each of the functions that evolves into the universal Isgur-Wise function in the heavy-quark limit [Sec. 5.3] and combined those results with the anal-ogous functions in B (s) → D (s) transitions [ Fig. 6]. The SCI predictions for zero-recoil values are commensurate with lattice-QCD calculations and the functional form matches the curve inferred from experiment. Recapitulating an observation from Ref. [6], when compared with B → D data from the Belle Collaboration, the SCI curve produces χ 2 /datum= 1.9.
Since the SCI delivers a qualitatively and semiquantitatively reliable expression of both Nature's mass generating mechanisms, i.e., emergent hadron mass and Higgs boson effects, it enables a comparison between the impacts of these effects on the entire array of electroweak transitions [ Fig. 9]. Exploiting this and focusing on the evolution of the maximum-recoil value of each form factor as the mass of the spectator quark, m spec , is varied, we found that some patterns emerged. For instance, V (0) tends to increase as m spec increases, whereas A 1,0 (0), f + (0) decrease; and A 2 (0) increases for c → d, s transitions, but decreases for b → u, c. These responses are expressions of the information encoded in the Ward-Green-Takahashi identities satisfied by the two pieces of each weak transition vertex.
Our analysis has highlighted the need for an extension of studies using symmetry-preserving formulations of realistic interactions in the analysis of semileptonic transitions involving heavy+heavy and heavy+light pseudoscalar mesons in the initial state. There is a scarcity of such theory in this area, especially for transitions with a vector meson in the final state. Work is therefore underway to adapt the framework exploited in Refs. [49] for use in treating B c → B (s) , B c → D, B (s) → D (s) , and all the transitions considered herein, with the exception of the J/ψ final state, which was studied in Ref. [48].
The defining quantity isG . After roughly twenty years of study, it is now known that, owing to the emergence of a gluon mass-scale in QCD [7,[92][93][94],G is nonzero and finite at infrared momenta; hence, can be written as follows: In QCD [94]: m G ≈ 0.5 GeV, α IR ≈ π. Following Ref. [66], we retain this value of m G and, profiting from the fact that a SCI cannot support relative momentum between bound-state constituents, simplify the tensor structure in Eqs. (A.1): A rudimentary form of confinement is implemented in the SCI by including an infrared regularisation scale, Λ ir , when defining all integral equations relevant to bound-state problems [95]. This artifice excises momenta below Λ ir , thereby eliminating quark+antiquark production thresholds [96]. The usual choice is Λ ir = 0.24 GeV = 1/[0.82 fm] [57], i.e., which relates to a length scale that is roughly the same as the proton radii [97].
All integrals appearing in SCI bound-state equations require ultraviolet regularisation in a step that breaks the link between infrared and ultraviolet scales that is characteristic of QCD. The associated ultraviolet mass-scales, Λ uv , thereby become physical parameters, which may be interpreted as upper bounds on the domains whereupon distributions within the associated systems are effectively momentum-independent. For instance, the ρ-meson is larger in size than the B *meson; hence, one should expect 1/Λ ρ uv > 1/Λ B * uv . As explained elsewhere [6,66,98] and sketched below, this observation leads to a completion of the SCI through introduction of a scale-dependent coupling.
For a quark of flavour f , the SCI gap equation is where m f is the associated quark current-mass. Using a Poincaré-invariant regularisation, the solution is with M f , the dynamically generated dressed-quark mass, obtained by solving where (τ 2 uv = 1/Λ 2 uv , τ 2 ir = 1/Λ 2 ir ) Here, the "iu" superscript stresses that the function depends on both the infrared and ultraviolet cutoffs and Γ (α, y) is the incomplete gamma-function. In general, functions of the following type arise in solving SCI bound-state equations: iu n (σ), n ∈ Z ≥ . Pseudoscalar mesons are generated as quark+antiquark bound-states, fḡ. They are described by a Bethe-Salpeter amplitude, whose SCI form is [6]: , Q is the bound-state's total momentum, Q 2 = −m 2 P , m P is the meson's mass. Additional details concerning these states as described by the SCI are provided elsewhere [6,Sec. 2]. It is nevertheless useful to specify the kernel: In arriving at the kernel in Eq. (A.10), we followed Ref. [6] in using the SCI identity which is a consequence of requiring that there are no quadratic or logarithmic divergences in the treatment of integral equations; namely, that shifting integration variables is permitted [99]. This condition is kindred to those implemented via dimensional-regularisation.
The SCI Bethe-Salpeter amplitude for a vector meson with polarisation λ is [57]: It is normalised canonically by rescaling such that where, with the trace over spinor indices and t + = t+Q: 15) In terms of the appropriate canonically normalised Bethe-Salpeter amplitudes, pseudoscalar and vector meson leptonic decay constants are given by the following formulae: With our normalisation, the empirical value of the pion's leptonic decay constant is f π = 0.092 GeV [70]. Improving upon the SCI introduced in Ref. [57], Ref. [6] kept all light-quark parameter values therein but determined the s-quark current mass, m s , and Kmeson ultraviolet cutoff, Λ K uv , through a least-squares fit to measured values of m K , f K whilst imposing the relation: This procedure eliminated one parameter by imposing the physical constraint that any increase in the momentum-space extent of a hadron wave function should be matched by a reduction in the effective coupling between the constituents. One useful consequence is that critical over-binding is avoided. The procedure was repeated for the c-quark/D-meson andb-quark/B-meson. The complete set of results is collected in Table 2.
The evolution of Λ uv with m P reported in Table 2 Table 3. Calculated axial-vector meson masses are also reported in Table 3, but the leptonic decay constants are omitted because empirical comparisons are unavailable. Kindred SCI predictions for these quantities are reported elsewhere [66, Table 1].
Considering the absolute value of the relative differences between the SCI results in Table 2 and listed comparison values, the median is 3.3% and the mean is 7.6% with a standard deviation of 10%. A more detailed discussion is presented in connection with Fig. 1 in Ref. [6]. Evidently, as noted therein, the SCI is not a precision tool for electroweak physics; but when employed judiciously, it is qualitatively reliable and usually a good quantitative guide.
Turning now to the electroweak vertex in Eq. (7), a detailed discussion of the vector part is contained in Ref. [6,Sec. 3.2]. Herein, therefore, we describe our SCI treatment of the axial-vector part. In paralleling that of the vector vertex, it effectively serves to recapitulate that analysis. We continue to use the cd transition as our exemplar.
Considering Eq. (A.22b), one readily obtains the following algebraic equation for the two terms in the pseudoscalar vertex: These identities are guaranteed so long as the regularisation ensures there are no logarithmic or quadratic divergences. As noted above, this is typically achieved by using Eq. (A.12). It is worth remarking that when refocused on the axial-vector vertex relevant for neutron β-decay and working in the chiral limit, the analogue of Eq. (A.28b) entails that the du transition vertex exhibits a pion pole at Q 2 = 0 whose residue is intimately connected with the dressed-quark mass. This has long been a known feature of QCD [105]