• The Bakamjian-Thomas relativistic quark model for hadron current matrix elements, while non-covariant at finite mass, is successful in the heavy quark limit : form factors are covariant and satisfy Isgur-Wise scaling and Bjorken-Uraltsev sum rules. Motivated by the so-called "1/2 vs. 3/2 puzzle" in B decays to positive parity D**, we examine the implications of the model at finite mass. In the elastic case 1/2^- -> 1/2^-, the HQET constraints for the O(1/m_Q) corrections are analytically fulfilled. A number of satisfying regularities is also found for inelastic transitions. We compute the form factors using the wave functions given by the Godfrey-Isgur potential. For 1/2^- \to 3/2^+ the departures from the heavy quark limit are small, but we find a strong enhancement in 1/2^- -> 1/2^+ (for 0^- -> 0^+). This enhancement is linked to a serious difficulty of the model at finite mass for the inelastic transitions, namely a violation of the HQET constraints at zero recoil formulated by Leibovich et al. These are nevertheless satisfied in the non-relativistic limit for the light quark. We conclude that these HQET rigorous constraints are crucial in the construction of a sensible relativistic quark model of inelastic form factors.
  • We propose a model independent extraction of the hadronic information needed to determine the photon polarization of the b-> s gamma process by the method utilizing the B -> K1 gamma -> K pi pi gamma angular distribution. We show that exactly the same hadronic information can be obtained by using the B -> J/psi K1 -> J/psi K pi pi channel, which leads to a much higher precision.
  • The Bakamjian-Thomas relativistic quark model, describing hadrons with a fixed number of constituents, yields in the heavy quark limit of QCD covariant Isgur-Wise functions and satisfies the whole tower of lowest moment sum rules (Bjorken-Uraltsev type sum rules). We first recall, as well as earlier results, the new formalism presented in our recent papers on Lorentz representations, which provide an elegant framework for the analysis of this model in the heavy quark limit and stress the results which have been already obtained in this direction. Then, we give some very explicit demonstrations of the fact that the Bakamjian-Thomas framework satifies the sum rules by considering simple cases of Isgur-Wise functions. In addition to the specific Bjorken and Uraltsev sum rules, an important sum rule that involves only heavy mesons with light cloud $j^P = {1 \over 2}^-$ and their radial excitations is demonstrated. This latter sum rule is phenomenologically interesting because it constrains the derivatives of the radially excited Isgur-Wise functions at zero recoil. On the other hand, we recall the limitations of the Bakamjian-Thomas scheme. At finite mass, current matrix elements with the current coupled to the heavy quark are no longer covariant, and higher moment sum rules that hold in the heavy quark limit of QCD are not satisfied.
  • It is recalled why, as already stated in a previous paper, there seems to be an inconsistency in identifying the non perturbative ZA,MOM as the renormalisation of the axial current, or equivalently, in setting as normalisation condition that the renormalised vertex=1 at p^2 = mu^2 at some renormalisation scale mu, where p is the momentum in the legs. Indeed, unlike the vector case, the Ward-Takahashi (WT) identity for the axial current is shown to imply both the renormalisation scale independence of ZA and a mu2 dependence of ZA,MOM. This mu^2 dependence is simply related to certain invariants in the pseudoscalar vertex and can persist in the chiral limit due to the spontaneous breaking of chiral symmetry (pion pole). It is seen clearly in the mu^2 dependence of some lattice calculations of ZA,MOM/ZV,MOM near the chiral limit.
  • We pursue the group theoretical method to study Isgur-Wise functions. We apply the general formalism, formerly applied to the baryon case j^P = 0^+ (for \Lambda_b -> \Lambda_c \ell \nu), to mesons with j^P = 1/2^-, i.e. $\overline{B} -> D(D^{(*)})\ell\nu. In this case, more involved from the angular momentum point of view, only the principal series of unitary representations of the Lorentz group contribute. We obtain an integral representation for the IW function xi(w) with a positive measure, recover the bounds for the slope and the curvature of xi(w) obtained from the Bjorken-Uraltsev sum rule method, and get new bounds for higher derivatives. We demonstrate also that if the lower bound for the slope is saturated, the measure is a delta-function, and xi(w) is given by an explicit elementary function. Inverting the integral formula, we obtain the measure in terms of the IW function, allowing to formulate criteria to decide if a given ansatz for the Isgur-Wise function is compatible or not with the sum rule constraints. Moreover, we have obtained an upper bound on the IW function valid for any value of w. We compare these theoretical constraints to a number of forms for \xi(w) proposed in the literature. The "dipole" function \xi(w) = (2/(w+1))^(2c) satisfies all constraints for c \geq 3/4, while the QCD Sum Rule result including condensates does not satisfy them. Special care is devoted to the Bakamjian-Thomas relativistic quark model in the heavy quark limit and to the description of the Lorentz group representation that underlies this model. Consistently, the IW function satisfies all Lorentz group criteria for any explicit form of the meson Hamiltonian at rest.
  • We discuss the possibility to measure in present experiments, especially LHCb, the non leptonic decay branching ratio $B \to D' \pi$, and emphasize phenomenological implications on $B \to D' l \nu$ semileptonic decay. We have estimated by lattice QCD the $D'$ decay constant $f_{D'}$ that parameterizes the $D'$ emission contribution to the Class-III non leptonic decay $B^- \to D^0 \pi^-$. In addition, we provide a new estimate of the decay constants $f_{D_{s,q}}$ which read $f_{D_{s}}=252(3)$ MeV and $f_{D_{s}}/f_{D}=1.23(1)(1)$.
  • We review the infrared properties of the pure Yang-Mills correlators and discuss recent results concerning the two classes of low-momentum solutions for them reported in literature; i.e. decoupling and scaling solutions. We will mainly focuss on the Landau gauge and pay special attention to the results inferred from the analysis of the Dyson-Schwinger equations of the theory and from "{\it quenched}" lattice QCD. The results obtained from properly interplaying both approaches are strongly emphasized.
  • We argue from LQCD that there is a non vanishing v.e.v of $A_a^\mu A^a_\mu$ in QCD in the Landau gauge. We use operator product expansion to provide a clear definition of $A_a^\mu A^a_\mu$ and extract a number both in the quenched and unquenched case.
  • We study the low-momentum ghost propagator Dyson-Schwinger equation (DSE) in Landau gauge, assuming for the truncation a constant ghost-gluon vertex, as it is extensively done, and a simple model for a massive gluon propagator. Then, regular DSE solutions (the zero-momentum ghost dressing function not diverging) appear to emerge and we show the ghost propagator to be described by an asymptotic expression reliable up to the order ${\cal O}(q^2)$. That expression, depending on the gluon mass and the zero-momentum Taylor-scheme effective charge, is proven to fit pretty well the low-momentum ghost propagator obtained through big-volume lattice simulations.
  • We study the quark mass function on hypercubic lattices, in a large range of physical volumes and cutoffs. To avoid the very large Wilson term artefact, we exploit the relation between the quark mass function and the pseudoscalar vertex in the continuum. We extrapolate to the chiral limit. In function of the physical volume, we observe a striking discontinuity in the properties of chiral extrapolation around a physical volume $L_c 6 (GeV}^{-1}=1.2 fm$. It is present in the quark mass function, which collapses to zero, as well as in the pion mass and the quark condensate as directly calculated from the pseudoscalar correlator. It is strongly reminiscent of the phenomenon of chiral symmetry restoration observed by Neuberger and Narayanan at $N_C=\infty$ around the same physical length. In the case of spontaneous symmetry breaking, we confirm that the OPE of the quark mass function, involving the quark condensate, is not operative at the available momenta, even taking into account the unusually large high order corrections to the Wilson coefficient calculated by Chetyrkin and Maier ; the gap remains large, around a factor 2, even at the largest momenta available to us (p \simeg GeV)
  • We discuss the different methods to obtain reliable informations about the deep infra-red behaviour of the gluon and ghost Green functions in QCD. We argue that a clever combination of analytical inputs and numerical ones is necessary. We illustrate this statement about the distinction between two classes of solutions of the ghost propagator Dyson-Schwinger equation (GPDSE). We conclude that the solution II ("decoupling") with a finite renormalised ghost dressing function at zero momentum is strongly favored by lattice QCD, We derive a method to solve numerically the GPDSE using lattice inputs concerning the gluon propagator. We derive an analytical small momentum expansion of the Ghost dressing function. We prove from the large cut-off behaviour of the ghost propagator renormalisation constant, $\widetilde Z_3$, that the bare ghost dressing function is infinite at the infinite cut-off limit.
  • Beyond the deep IR, the analysis of ghost and gluon propagators still keeps very interesting non-perturbative information. The Taylor-scheme coupling can be computed and applied to obtain the $\Lambda_{\rm QCD}$ parameter from Landau gauge lattice simulations. Furthermore, a dimension-two gluon condensate, that can be understood in the instanton liquid model, plays an important role in the game.
  • We study a relation recently derived by K. Kondo at zero momentum between the Zwanziger's horizon function, the ghost dressing function and Kugo's functions $u$ and $w$. We agree with this result as far as bare quantities are considered. However, assuming the validity of the horizon gap equation, we argue that the solution $w(0)=0$ is not acceptable since it would lead to a vanishing renormalised ghost dressing function. On the contrary, when the cut-off goes to infinity, $u(0) \to \infty$, $w(0) \to -\infty$ such that $u(0)+w(0) \to -1$. Furthermore $w$ and $u$ are not multiplicatively renormalisable. Relaxing the gap equation allows $w(0)=0$ with $u(0) \to -1$. In both cases the bare ghost dressing function, $F(0,\Lambda)$, goes logarithmically to infinity at infinite cut-off. We show that, although the lattice results provide bare results not so different from the $F(0,\Lambda)=3$ solution, this is an accident due to the fact that the lattice cut-offs lie in the range 1-3 GeV$^{-1}$. We show that the renormalised ghost dressing function should be finite and non-zero at zero momentum and can be reliably estimated on the lattice up to powers of the lattice spacing ; from published data on a $80^4$ lattice at $\beta=5.7$ we obtain $F_R(0,\mu=1.5$ GeV)$\simeq 2.2$.
  • We propose a group theoretical method to study Isgur-Wise functions. A current matrix element splits into a heavy quark matrix element and an overlap of the initial and final clouds, related to the IW functions, that contain the long distance physics. The light cloud belongs to the Hilbert space of a unitary representation of the Lorentz group. Decomposing into irreducible representations one obtains the IW function as an integral formula, superposition of irreducible IW functions with positive measures, providing positivity bounds on its derivatives. Our method is equivalent to the sum rule approach, but sheds another light on the physics and summarizes and gives all its possible constraints. We expose the general formalism, thoroughly applying it to the case j = 0 for the light cloud, relevant to the semileptonic decay Lambda_b -> Lambda_c + l + nu. In this case, the principal series of the representations contribute, and also the supplementary series. We recover the bound for the curvature of the j = 0 IW function xi_Lambda (w) that we did obtain from the sum rule method, and we get new bounds for higher derivatives. We demonstrate also that if the lower bound for the curvature is saturated, then xi_Lambda (w) is completely determined, given by an explicit elementary function. We give criteria to decide if any ansatz for the Isgur-Wise function is compatible or not with the sum rules. We apply the method to some simple model forms proposed in the literature. Dealing with a Hilbert space, the sum rules are convergent, but this feature does not survive hard gluon radiative corrections.
  • In the heavy quark limit of QCD, using the Operator Product Expansion, the formalism of Falk for hadrons or arbitrary spin, and the non-forward amplitude, as proposed by Uraltsev, we formulate sum rules involving the Isgur-Wise function $\xi_{\Lambda} (w)$ of the baryon transition $\Lambda_b \to \Lambda_c \ell \overline{\nu}_{\ell}$, where the light cloud has $j^P=0^+$ for both initial and final baryons. We recover the lower bound for the slope $\rho_\Lambda^2 = - \xi '_\Lambda (1) \geq 0$ obtained by Isgur et al., and we generalize it by demonstrating that the IW function $\xi_{\Lambda} (w)$ is an alternate series in powers of $(w-1)$, i.e. $(-1)^n \xi_{\Lambda}^{(n)} (1) \geq 0$. Moreover, exploiting systematically the sum rules, we get an improved lower bound for the curvature in terms of the slope, $\sigma_\Lambda^2 = \xi "_\Lambda (1) \geq {3 \over 5} [\rho_\Lambda^2 + (\rho_\Lambda^2)^2]$. This bound constrains the shape of the Isgur-Wise function and it will be compelling in the analysis of future precise data on the differential rate of the baryon semileptonic decay $\Lambda_b \to \Lambda_c \ell \overline{\nu}_{\ell}$, that has a large measured branching ratio, of about 5%.
  • We compute a formula including OPE power corrections to describe the running of a QCD coupling non-perturbatively defined through the ghost and gluon dressing functions. This turns out to be rather accurate. We propose the ``{\it plateau}''-procedure to compute $\Lambda_{\bar{\rm MS}}$ from the lattice computation of the running coupling constant. We show a good agreement between the different methods which have been used to estimate $\Lambda_{\bar{\rm MS}}^{N_f=0}$. We argue that $\Lambda_{\bar{\rm MS}}$ or the strong coupling constant computed with different lattice spacings may be used to estimate the lattice spacing ratio.
  • We solve numerically the Schwinger-Dyson (SD hereafter) ghost equation in the Landau gauge for a given gluon propagator finite at k=0 (alpha_gluon=1) and with the usual assumption of constancy of the ghost-gluon vertex ; we show that there exist two possible types of ghost dressing function solutions, as we have previously inferred from analytical considerations : one singular at zero momentum, satisfying the familiar relation alpha_gluon+2 alpha_ghost=0 between the infrared exponents of the gluon and ghost dressing functions(in short, respectively alpha_G and alpha_F) and having therefore alpha_ghost=-1/2, and another which is finite at the origin (alpha_ghost=0), which violates the relation. It is most important that the type of solution which is realized depends on the value of the coupling constant. There are regular ones for any coupling below some value, while there is only one singular solution, obtained only at a critical value of the coupling. For all momenta k<1.5 GeV where they can be trusted, our lattice data exclude neatly the singular one, and agree very well with the regular solution we obtain at a coupling constant compatible with the bare lattice value.
  • We examine analytically the ghost propagator Dyson-Schwinger Equation (DSE) in the deep IR regime and prove that a finite ghost dressing function at vanishing momentum is an alternative solution (solution II) to the usually assumed divergent one (solution I). We furthermore find that the Slavnov-Taylor identities discriminate between these two classes of solutions and strongly support the solution II. The latter turns out to be also preferred by lattice simulations within numerical uncertainties.
  • The Bakamjian-Thomas relativistic quark model provides a Poincar\'e representation of bound states with a fixed number of constituents and, in the heavy quark limit, form factors of currents satisfy covariance and Isgur-Wise scaling. We compute the Light Cone Distribution Amplitudes of $B$ mesons $\phi_{\pm}^B(\omega)$ as well as the Shape Function $S(\omega)$, that enters in the decay $B \to X_s \gamma$, that are also covariant in this class of models. The LCDA and the SF are related through the quark model wave function. The former satisfy, in the limit of vanishing constituent light quark mass, the integral relation given by QCD in the valence sector of Fock space. Using a gaussian wave function, the obtained $S(\omega)$ is identical to the so-called Roman Shape Function. From the parameters for the latter that fit the $B \to X_s\gamma$ spectrum we predict the behaviour of $\phi_{\pm}^B(\omega)$. We discuss the important role played by the constituent light quark mass. In particular, although $\phi_-^B(0) \not= 0$ for vanishing light quark mass, a non-vanishing mass implies the unfamiliar result $\phi_-^B (0) = 0$. Moreover, we incorporate the short distance behaviour of QCD to $\phi_+^B (\omega)$, which has sizeable effects at large $\omega$. We obtain the values for the parameters $\bar{\Lambda} \cong 0.35$ GeV and $\lambda_B^{-1} \cong 1.43$ GeV$^{-1}$. We compare with other theoretical approaches and illustrate the great variety of models found in the literature for the functions $\phi_{\pm}^B (\omega)$; hence the necessity of imposing further constraints as in the present paper. We briefly review also the different phenomena that are sensitive to the LCDA.
  • The OPE treatment that has been so successful in describing inclusive $\bar B \to l \bar \nu X_c$ decays yields sum rules (in particular the Uraltsev sum rule and its higher moments) implying the dominance of the $P$ wave $j_q = 3/2$ charm states in $X_c$ over their $j_q=1/2$ counterparts. This prediction is supported by other general arguments as well as quark model calculations, which illustrate the OPE results, and by preliminary lattice findings. Its failure would indicate a significant limitation in our theoretical understanding of $\bar B \to l \bar \nu X_c$. Some experimental issues have been clarified since a preliminary version of this note had appeared, yet the verdict on the composition of the final states {\em beyond} $D$, $D^*$ and the two narrow $j_q = 3/2$ resonances remains unsettled. Establishing which hadronic configurations -- $D/D^* + \pi, D/D^* + 2 \pi, ...$ -- contribute, what their quantum numbers are and their mass distributions will require considerable experimental effort. We explain the theoretical issues involved and why a better understanding of them will be of significant value. Having significant contributions from a mass continuum distribution below 2.5 GeV raises serious theoretical questions for which we have no good answer. Two lists are given, one with measurements that need to be done and one with items of theoretical homework. Some of the latter can be done by employing existing theoretical tools, whereas others need new ideas.
  • We revisit the problem of interplay between the strong and the Coulomb interaction in the charged-to-neutral yield ratio for $B {\bar B}$ and $D {\bar D}$ pairs near their respective thresholds in $e^+e^-$ annihilation. We consider here a realistic situation with a resonant interaction in the isospin I=0 channel and a nonresonant strong scattering amplitude in the I=1 state. We find that the yield ratio has a smooth behavior depending on the scattering phase in the I=1 channel. The same approach is also applicable to the $K {\bar K}$ production at the $\phi(1020)$ resonance, where the Coulomb effect in the charged-to-neutral yield ratio is generally sensitive to the scattering phases in both the isoscalar and the isovector channels. Furthermore, we apply the same approach to the treatment of the effect of the isotopic mass difference between the charged and neutral mesons and argue that the strong-scattering effects generally result in a modification to the pure kinematical effect of this mass difference.
  • We consider the constraints of the Slavnov-Taylor identity of the IR behaviour of gluon and ghost propagators and their compatibility with solutions of the ghost Dyson-Schwinger equation and with the lattice picture.
  • From sum rules in the heavy quark limit of QCD, using the non-forward amplitude, we demonstrate that if the slope rho^2 = -xsi'(1) of the Isgur-Wise function xsi(w) attains its lower bound 3/4 (as happens in the BPS limit proposed by Uraltsev), the IW function is completely determined, given by the function xsi(w) = [2/(w+1)]^(3/2).
  • Using previously formulated sum rules in the heavy quark limit of QCD, we demonstrate that if the slope rho^2 = -xi'(1) of the Isgur-Wise function xi(w) attains its lower bound 3/4, then all the derivatives (-1)^L xi^(L)(1) attain their lower bounds (2L+1)!!/2^(2L), obtained by Le Yaouanc et al. This implies that the IW function is completely determined, given by the function xi(w) = [2/(w+1)]^(3/2). Since the so-called BPS condition proposed by Uraltsev implies rho^2 = 3/4, it implies also that the IW function is given by the preceding expression.
  • We extend the formalism of Leibovich, Ligeti, Stewart and Wise in the 1/m_Q expansion of Heavy Quark Effective Theory for the B semileptonic decays into excited D[(3/2)^+] mesons to the opposite parity states D[(3/2)^-]. For D[(3/2)^+] the 1/m_Q current perturbation dominates over the leading term at zero recoil, while for D[(3/2)^-] the 1/m_Q perturbation due to L_mag dominates also at zero recoil. We show that the corresponding 1/m_Q magnetic coupling is proportional to the mixing between the states D[(3/2)^-] and D[(1/2)^-] induced by the tensor force. We point out some subtleties that appear in this respect in HQET.