
The BakamjianThomas relativistic quark model for hadron current matrix
elements, while noncovariant at finite mass, is successful in the heavy quark
limit : form factors are covariant and satisfy IsgurWise scaling and
BjorkenUraltsev sum rules. Motivated by the socalled "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 GodfreyIsgur 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
nonrelativistic 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 BakamjianThomas relativistic quark model, describing hadrons with a
fixed number of constituents, yields in the heavy quark limit of QCD covariant
IsgurWise functions and satisfies the whole tower of lowest moment sum rules
(BjorkenUraltsev 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
BakamjianThomas framework satifies the sum rules by considering simple cases
of IsgurWise 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 IsgurWise functions at zero recoil. On the
other hand, we recall the limitations of the BakamjianThomas 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 WardTakahashi (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 IsgurWise 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
BjorkenUraltsev 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 deltafunction, 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
IsgurWise 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 BakamjianThomas 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 ClassIII 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 YangMills correlators and
discuss recent results concerning the two classes of lowmomentum 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 DysonSchwinger 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 lowmomentum ghost propagator DysonSchwinger equation (DSE) in
Landau gauge, assuming for the truncation a constant ghostgluon vertex, as it
is extensively done, and a simple model for a massive gluon propagator. Then,
regular DSE solutions (the zeromomentum 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 zeromomentum Taylorscheme
effective charge, is proven to fit pretty well the lowmomentum ghost
propagator obtained through bigvolume 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 infrared 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 DysonSchwinger 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 cutoff behaviour of the ghost
propagator renormalisation constant, $\widetilde Z_3$, that the bare ghost
dressing function is infinite at the infinite cutoff limit.

Beyond the deep IR, the analysis of ghost and gluon propagators still keeps
very interesting nonperturbative information. The Taylorscheme coupling can
be computed and applied to obtain the $\Lambda_{\rm QCD}$ parameter from Landau
gauge lattice simulations. Furthermore, a dimensiontwo 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
cutoff 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 cutoff. 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 cutoffs lie in
the range 13 GeV$^{1}$. We show that the renormalised ghost dressing function
should be finite and nonzero 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 IsgurWise 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 IsgurWise 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 nonforward amplitude,
as proposed by Uraltsev, we formulate sum rules involving the IsgurWise
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 $(w1)$, 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 IsgurWise 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 nonperturbatively 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 SchwingerDyson (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 ghostgluon 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 DysonSchwinger 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 SlavnovTaylor
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 BakamjianThomas 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 IsgurWise
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 socalled
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
nonvanishing 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 chargedtoneutral 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
chargedtoneutral 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 strongscattering
effects generally result in a modification to the pure kinematical effect of
this mass difference.

We consider the constraints of the SlavnovTaylor identity of the IR
behaviour of gluon and ghost propagators and their compatibility with solutions
of the ghost DysonSchwinger equation and with the lattice picture.

From sum rules in the heavy quark limit of QCD, using the nonforward
amplitude, we demonstrate that if the slope rho^2 = xsi'(1) of the IsgurWise
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 IsgurWise 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 socalled 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.