9512.net

甜梦文库

甜梦文库

当前位置：首页 >> >> # Heavy Quark Masses, Mixing Angles, and Spin-Flavour Symmetry

.0.0 CERN-TH.7225/94

arXiv:hep-ph/9404296v1 19 Apr 1994

HEAVY QUARK MASSES, MIXING ANGLES, AND SPIN-FLAVOUR SYMMETRY

MATTHIAS NEUBERT? Theory Division, CERN CH-1211 Geneva 23, Switzerland

ABSTRACT In a series of three lectures, I review the theory and phenomenology of heavy quark masses and mixing angles and the status of their determination. In addition, I give an introduction to heavy quark symmetry, with the main emphasis on the development of heavy quark e?ective ?eld theory and its application to obtain a model-independent determination of | Vcb |.

(Lectures presented at TASI-93, Boulder, Colorado, 1993)

CERN-TH.7225/94 April 1994

?

Address before Oct. 1993: Stanford Linear Accelerator Center, Stanford University, Stanford, California 94309

Introduction The rich phenomenology of weak decays has always been a source of information about the structure of elementary particle interactions. A long time ago, β - and ?decay experiments revealed the nature of the e?ective ?avour-changing interactions at low momentum transfer. Today, we are in a similar situation: Weak decays of hadrons containing heavy quarks are employed for tests of the standard model and measurements of its parameters. In particular, they o?er the most direct way to determine the weak mixing angles and to test the unitarity of the Kobayashi–Maskawa matrix. On the other hand, hadronic weak decays also serve as a probe of that part of strong-interaction phenomenology which is least understood: the con?nement of quarks and gluons into hadrons. In fact, it is this intricate interplay between weak and strong interactions that makes this ?eld challenging and attractive to many theorists. Over the last decade or so, a lot of information has been collected about heavy quark decays from experiments on the Υ(4s) resonance, and more recently at e+ e? and hadron colliders. This has led to a rather detailed knowledge of the ?avour sector of the standard model and many of the parameters associated with it. There ?d mixing1,2 , b → u have been several great discoveries in this ?eld, such as Bd –B transitions3 ? 6 , and rare decays induced by penguin operators7 . Yet there is much more to come. Hopefully, the approval of the ?rst B -meson factory at SLAC has opened the way for a bright future for B -physics. At the same time, upgrades of the existing facilities at Cornell, Fermilab, and LEP will provide a wealth of data within the coming years. However, my main purpose here is to talk about theory, and fortunately there has been a lot of progress and enthusiasm in this ?eld in recent years. This is related to the discovery of heavy quark symmetry8 ? 13 and the development of the heavy quark e?ective theory14 ? 24 , which is a low-energy e?ective theory that describes the strong interactions of a heavy quark with light quarks and gluons. The excitement about these developments is caused by the fact that they allow (some) model-independent predictions in an area in which “progress” in theory often meant nothing more than the construction of a new model, which could be used to estimate some stronginteraction hadronic matrix elements. Therefore, I hope that you do not mind it if the part about the heavy quark e?ective theory constitutes the main body of these notes; this is where the main progress has been achieved from the theoretical point of view. Also, I have recently ?nished a long review article on heavy quark symmetry25 , which clearly simpli?es the task of writing up lecture notes. In the ?rst lecture, I will discuss the theory of heavy quark masses, their de?nition in perturbation theory, and their determination from QCD sum rules. The second lecture provides an introduction to the standard model description of quark mixing, di?erent parametrizations of the Kobayashi–Maskawa matrix, the status of the determination of the mixing angles, and the physics of the unitarity triangle. The third lecture, which covers more than half of these notes, is devoted to a review of the ideas on heavy quark symmetry and the formalism of the heavy quark 1

e?ective theory. In particular, I will discuss the theory of the model-independent determination of | Vcb| from exclusive semileptonic B → D ? ? ν ? decays. At the end of each lecture you will ?nd some suggestions for little exercises, which are typically related to the derivation of important equations given in the notes. I have tried to select problems that are fun to solve and contain some interesting pieces of physics. You are invited to see if I am right.

2

1. Heavy Quark Masses Because of con?nement at large distance scales, quarks and gluons do not appear among the physical states of QCD. There is thus no natural, physical de?nition of quark masses. Rather, several de?nitions are possible and have been proposed, and it is often a matter of convenience which one to use. Most of these de?nitions are tied to the framework of perturbation theory. In this lecture, I will discuss some of the most common mass de?nitions and their interrelation. Special emphasis is put on the discussion of the running quark mass in the context of the renormalization group. I will then discuss how values for the bottom and charm quark masses are obtained from QCD sum rules. Very brie?y, I will touch upon the subtle issue of an infrared renomalon in the pole mass, which implies an intrinsic uncertainty in these determinations. 1.1. Quark Mass De?nitions For heavy quarks, the nonrelativistic bound-state picture suggests the notion of a pole mass mpole de?ned, order by order in perturbation theory, as the position of Q the singularity in the renormalized quark propagator. One can show that the pole mass is gauge-invariant, infrared ?nite, and renormalization-scheme independent26 . It is thus a meaningful “physical” parameter as long as the heavy quark is not exactly on-shell. For instance, the pole mass appears in the formula for the energy levels in quarkonium systems. An alternative gauge-invariant de?nition is provided by the running quark mass in some subtraction scheme. One usually works with the modi?ed minimal subtraction (MS) scheme27 and denotes the running mass by mQ (?). Running quark masses are used when there is a large momentum scale ? ? mQ in the problem, so as to absorb large logarithms, which would otherwise render perturbation theory invalid. Of course, one can use other subtraction schemes such as minimal subtraction28 (MS). Sometimes, it may even be convenient to use a gauge-dependent de?nition of a heavy quark mass such as the so-called Euclidean mass mEucl Q , which is de?ned 29 2 2 by a subtraction at the Euclidean point p = ?mQ . What is important is that these perturbative de?nitions are related to each other in a calculable way, order by order in perturbation theory. One may also think of obvious nonperturbative de?nitions of a heavy quark mass. For instance, one could de?ne mQ to be one half of the mass of the ground state in ? ) quarkonium system, or as the mass of the lightest hadron the corresponding (QQ containing the heavy quark, or as that mass minus some ?xed binding energy, etc. The problem is, of course, that these de?nitions are ad hoc, and relating them to each other or to the perturbative de?nitions given above would require to solve QCD, a task that is presently beyond our calculational skills. There is, however, a bridge between the two classes of de?nitions, which is provided by the heavy quark e?ective theory, which is an e?ective ?eld theory appropriate to describe the soft interactions of an almost on-shell heavy quark with light quarks and gluons. There, 3

? , which corresponds to the e?ective mass of the light one can de?ne a parameter Λ degrees of freedom in a hadron HQ containing a single heavy quark Q, in terms of a gauge-invariant, infrared ?nite, and renormalization-scheme independent hadronic 31 ? matrix element30 . One can then de?ne a heavy quark mass m? m? Q as Q = mHQ ? Λ. ? In perturbation theory, and up to corrections of order 1/mQ , the mass mQ coincides with the pole mass. However, the above de?nition is more general, as it does not rely on a perturbative expansion. 1.2. Quark Masses in Perturbation Theory Let me now discuss the concept of heavy quark masses to one-loop order in perturbation theory. The bare quark mass mbare appearing in the QCD Lagrangian Q is related to the renormalized mass mren by a counter term, Q mbare = mren Q Q ? δmQ , (1)

where mbare and δmQ are divergent quantities. The counter term is chosen such Q that the renormalized mass is ?nite. Its value depends, however, on the subtraction prescription. In the MS scheme, the renormalized mass will depend on some subtraction scale ?, i.e. mren Q = mQ (?). To calculate δmQ in perturbation theory, one has to evaluate the heavy quark self-energy Σ(p) shown in Fig. 1. The relation is δmQ = divergent part of Σ(p / = mQ ) + scheme-dependent ?nite terms. (2)

To calculate the self-energy at one-loop order, it is convenient to use dimensional regularization32,33 , i.e. to work in D = 4 ? 2? space-time dimensions. Then the renormalized coupling constant is related to the bare one by34

bare 2 αs = Zα αs (?) ?2? = αs ?2? + O (αs ).

(3)

From a straightforward calculation, one obtains at one-loop order

?? αs m2 3 ? 2? Q Σ(p / = mQ ) = mQ , Γ(?) 2 3π 4π? 1 ? 2? αs 1 4 ?2 = mQ + O (?) , + ln 2 + π ? ? mQ 3

(4)

where 1/? ? = 1/? ? γ + ln 4π . Notice that this result is still ?-independent, since the ? dependence of αs /? ? cancels the ?-dependent logarithm. A scale dependence appears only when one subtracts the ultraviolet divergence. In the MS scheme, the mass counter term is de?ned to subtract the 1/? ?-pole in the self-energy, i.e. δmMS Q = mQ αs , π? ?

mQ (?) = mQ 1 + 4

αs . π? ?

(5)

The pole mass, on the other hand, is de?ned so as to absorb the entire self-energy on-shell. Hence δmpole = Σ(p / = mQ ) , Q 4 ?2 αs 1 pole + ln 2 + mQ = mQ 1 + π ? ? mQ 3 Comparing the two results, one obtains 1? mQ (mQ ) = mpole Q 4αs (mQ ) , 3π ?2 αs . ln mQ (?) = mQ (mQ ) 1 ? π m2 Q

.

(6)

(7)

The ?rst equation relates di?erent perturbative de?nitions of the heavy quark mass, namely the pole mass to the running mass in the MS scheme evaluated at the scale ? = mQ . I have written this relation as a perturbative expansion in powers of αs (mQ ). This expansion is also known to two-loop order from a calculation by Gray et al.35 Similar relations exist in other subtraction schemes, e.g.

pole MS 1? mMS Q (mQ ) = mQ pole Eucl mEucl Q (mQ ) = mQ

αs (mMS Q ) 4 ? γ + ln 4π π 3 αs (mEucl Q ) 2 ln 2 , 1? π

, (8)

where the gauge-dependent Euclidean mass is de?ned in the Landau gauge29 .

i

(p) =

p

p

=

+

:::

Figure 1: The quark self-energy to one-loop order in QCD. The second relation in (7) determines the running of the heavy quark mass. In QCD, quarks become “lighter” at high-energy scales. The simple one-loop calculation presented above is, however, not su?cient to scale mQ (?) up to very large scales. For instance, using (7) to extrapolate the mass of the bottom quark up to a typical grand uni?cation scale MGUT ? 1016 GeV, one would obtain a meaningless result: mb (MGUT ) αs αs M2 ? 1 ? 71 =1? ln GUT =? (9) 2 mb (mb ) π mb π Of course, the problem of large logarithms in perturbation theory is not speci?c to this case, but is encountered frequently. Large logarithms can be controlled in a 5

systematic way using the beautiful machinery of the renormalization group. It is worth going through the solution in great detail. 1.3. Renormalization-Group Improvement For the purpose of this section, I will switch to a slightly more general notation and rewrite (7) as m(?) = m(?0 ) 1 ? ?2 αs 2 ln 2 + O (αs ) . π ?0 (10)

If ? ? ?0 , the logarithm can be so large that αs ln ?2 /?2 0 becomes of order unity, and ordinary perturbation theory breaks down. In fact, one can show that at any given order in perturbation theory there will be large logarithms of the type αs ln ?2 ?2 0

n

?

ln ?2 /?2 0 ln ?2 /Λ2 QCD

n

,

(11)

where I have used the fact that the running coupling constant scales like αs (?) ? 1/ ln(?2 /Λ2 QCD ). It is necessary to resum these “leading logarithms” to all orders in perturbation theory. This is achieved by solving the renormalization-group equation (RGE) for the running quark mass, ? d m(?) = γ (αs ) m(?) . d? (12)

The anomalous dimension γ has a perturbative expansion in the renormalized coupling constant: αs αs 2 + γ1 γ (αs ) = γ0 + ... . (13) 4π 4π The coe?cients γi are known to three-loop order36,37 . For our purposes, however, it is su?cient to note that26 γ0 = ?8 , γ1 = ? 404 40 + nf , 3 9 (14)

where nf denotes the number of quark ?avours with mass below ?. Throughout these notes I will evaluate the QCD coe?cients for Nc = 3 colours, and only display the dependence on the number of ?avours explicitly. The value of γ0 follows from the one-loop result (10). The next step is to rewrite the RGE in the form of a partial di?erential equation, making explicit the scale dependence of the renormalized coupling constant. This gives ? ? ? + β (αs ) ? γ (αs ) m(?) = 0 . (15) ?? ?αs (?) 6

The β -function β (αs ) = ? αs αs ?αs (?) = ?2αs β0 + β1 ?? 4π 4π

2

+ ...

(16)

describes the running of the coupling constant. The one- and two-loop coe?cients are38 ? 42 2 38 β0 = 11 ? nf , β1 = 102 ? nf . (17) 3 3 The exact solution of the RGE can be written in the form m(?) = U (?, ?0 ) m(?0 ) , with the evolution operator27,43,44 αs (?) U (?, ?0 ) = exp αs (?0 ) dα γ (α ) . β (α ) (18)

(19)

The trick is to obtain a perturbative expansion in the exponent of this expression by inserting the expansions for the anomalous dimension and β -function. After a simple calculation, one ?nds U (?, ?0 ) = αs (?0 ) αs (?) γ0 /2β0 1+ αs (?0 ) ? αs (?) γ1 β0 ? β1 γ0 2 + O (αs ) . 2 4π 2β0 (20)

In this result, the running coupling constant has two-loop accuracy. The corresponding expression is β1 ln ln(?2 /Λ2 4π QCD ) 1 ? , αs (?) = 2 β0 ln(?2 /Λ2 β0 ln(?2 /Λ2 QCD ) QCD ) (21)

where ΛQCD is a scheme-dependent scale parameter? . The factor containing the ratio of the running coupling constants in (20) sums the leading logarithms (11) to all orders in perturbation theory. Keeping just this factor corresponds to the so-called leading logarithmic approximation (LLA), which is often used to attack problems containing widely separated scales. Note that to calculate this factor, one only needs to compute the one-loop coe?cient of the anomalous dimension. This is a rather trivial task, as it su?ces to calculate the 1/?-pole in the quark self-energy. Once the leading scaling behaviour has been factored out, perturbation theory becomes well-behaved, i.e. the next-to-leading corrections in the parenthesis in (20) obey a perturbative expansion free of large logarithms. This is why the approach is

?

The value of αs (?0 ) at some reference scale ?0 can be used to eliminate the dependence on ΛQCD . Nowadays, it is convenient to choose ?0 = mZ , since αs (mZ ) is known with high accuracy.

7

called “renormalization-group improved perturbation theory”. The terms of order αs , which I have shown explicitly, correspond to the next-to-leading logarithmic approximation (NLLA). Including them, one sums logarithms of the type αs ?2 αs ln 2 ?0

n

(22)

to all orders. To achieve such an accuracy, it is necessary to calculate the two-loop coe?cients of the anomalous dimension and β -function. Combining the above results and setting ?0 = mQ , one obtains the running quark mass to next-to-leading order in renormalization-group improved perturbation theory. The result is αs (?) mQ (?) = mQ (mQ ) αs (mQ ) where

4/β0

1+S

αs (mQ ) ? αs (?) 2 + O (αs ) , π

(23)

5 34 107 S=? + ? 2 . 6 3β0 β0

(24)

It is to be understood that the number of ?avours changes as ? crosses a quark threshold45 . For instance, when one uses (23) to scale the running bottom quark mass up to large scales, nf changes from 5 to 6 when ? crosses the top quark mass. At the same time, the QCD scale parameter ΛQCD in expression (21) for the running coupling constant changes, so that αs (?) is a continuous function. From the point of view of convergence of perturbation theory, eq. (23) can be used to evaluate the running mass at an arbitrarily large scale ?. However, I have derived this result assuming that there are only QCD interactions. But at high-energy scales, other gauge and Yukawa interactions become important as well. They are most signi?cant for the running of the top quark mass, as I will discuss in detail in the next section. In the case of the bottom quark, the result that renormalization e?ects lower the running mass at high-energy scales remains true. In fact, in many extensions of the standard model it is possible to obtain a uni?cation of the bottom quark and the tau lepton masses at a typical grand uni?cation scale: mb (MGUT ) = mτ (MGUT ) for MGUT ? 1016 GeV . (25)

For more details on this, I refer to the lectures by L. Hall in this volume. 1.4. Running Top Quark Mass The interplay of gauge and Yukawa interactions leads to very interesting e?ects on the evolution of the top quark mass46,47 . In my discussion below, I will focus on the standard model. The analysis for extensions of the standard model such as supersymmetry proceeds in a similar way48 ? 50 . 8

Let me write the running top quark mass in terms of the Yukawa coupling λt (?): v (? ) mt (?) = λt (?) √ , 2 (26)

where v denotes the vacuum expectation value of the Higgs ?eld, normalized so that v (mW ) ? 246 GeV. It is convenient to de?ne a running coupling constant αt (?) as αt (?) = λ2 t (? ) . 4π (27)

The most important contributions to the evolution of the running mass mt (?) come from the QCD interactions as well as from the large top Yukawa coupling. The relevant one-loop diagrams are depicted in Fig. 2. They lead to the RGE ? d αs (?) 3 αt (?) ln mt (?) = ?8 + + ... , d? 4π 2 4π (28)

where the ellipses represent other contributions, which I will neglect. Combining this with the RGE for v (?), ? one obtains ? d αt (?) ln v (?) = ?3 + ... , d? 4π (29)

d αs (?) αt (?) ln αt (?) = ?16 +9 + ... . d? 4π 4π

(30)

To get an idea of the structure of this equation, suppose ?rst that the QCD coupling constant is not running, i.e. αs = const. Then the RGE has an infrared-stable ?xed αs , meaning that irrespective of the value of αt (?) at large scales, the point at αt = 16 9 Yukawa coupling is attracted into the ?xed point as ? decreases. This is illustrated √ in Fig. 3. The resulting ?xed-point value of the top quark mass, mt = 4 2παs v , 3 is entirely determined by the low-energy group structure of the theory, which is responsible for the one-loop coe?cients in (30).

g

H

+

t t t t

+

:::

Figure 2: Dominant one-loop contributions to the self-energy of the top quark in the standard model.

9

4 3 2 1 0 -32

9/16 αt(?)/αs

-24

-16 ln(?/M)

-8

0

Figure 3: The ?xed-point structure of the RGE (30) in the case αs = const., with 4αs /π = 0.1. For running αs (?), one has to ?nd the simultaneous solution of (30) and the evolution equation for the strong coupling constant, ? d αs (?) ln αs (?) = ?14 . d? 4π (31)

Depending on the initial conditions, one ?nds that there are still “quasi-?xed point” solutions, where the running of αt (?) is connected to the running of αs (?), but independent of the value of the Yukawa coupling at large scales. In leading logarithmic approximation, the exact solution is not too hard to obtain. It reads αt (?) = αt (M ) αs (?) αs (M )

8/7

9 αt (M ) 1+ 2 αs (M )

αs (?) αs (M )

1/7

?1

?1

,

(32)

where M ? ? denotes some large mass scale, at which the initial conditions are imposed. In order to illuminate the structure of this equation, let me distinguish three cases: (i) αt (M ) ? αs (M ): In this case, the top Yukawa coupling is weak, and (32) can be approximated by αs (?) αt (?) ? αt (M ) αs (M )

8/7

,

(33)

which is nothing but the standard QCD evolution [cf. (23)].

10

α (M ): (ii) αt (M ) = 2 9 s This case corresponds to the quasi-?xed point obtained by Pendleton and Ross46 , where (32) simpli?es to 2 αt (?) = αs (?) . (34) 9 However, to obtain this solution requires a ?ne-tuning at the large scale M .

0.3

?=1016 GeV

0.2 αt(?)

mt=210 GeV

0.1

?=mt

0 0.02

0.04

0.06

0.08 αs(?)

0.1

0.12

Figure 4: Evolution of the coupling constants for di?erent initial values of αt (M ). I assume that M = 1016 GeV. (iii) αt (M ) ? αs (M ): This is the most interesting case of a large top Yukawa coupling. Then the approximate solution of (32) reads αs (M ) 2 αt (?) ? αs (?) 1 ? 9 αs (?)

1/7 ?1

.

(35)

This is the quasi-?xed point discovered by Hill47. At low scales ? ? M , αt (?) follows αs (?) and becomes independent of the initial value αt (M ). The dynamics that generates the top Yukawa coupling at some large scale M cannot be tested. An example of this mechanism is provided by the top condensation model of Bardeen, Hill, and Lindner51 . Using (35), one obtains for the ?xed-point value of the top quark mass 2v (mt ) mt (mt ) ? 3 αs (M ) π αs (mt ) 1 ? αs (mt )

1/7 ?1/2

? 210 GeV ,

(36)

where I have assumed M ? 1016 GeV. Again, the result is almost entirely determined by the group structure of the low-energy theory. The dependence on M is very weak. 11

In Fig. 4, I show the evolution of the coupling constants according to (32) for a set of initial values of αt (M ). The quasi-?xed point behaviour is clearly visible once αt (M ) is large enough. A more careful analysis including all standard model contributions in the RGE gives47 mt ? 225 GeV (for M ? 1016 GeV). Somewhat smaller values mt < 200 GeV are obtained in the minimal supersymmetric standard model48 ? 50 , depending however upon the value of tan β . 1.5. Determination of mb and mc from QCD Sum Rules In this section, I will discuss the extraction of the masses of the bottom and charm quarks from QCD sum rules, and in particular from the analysis of quarkonium spectra. The main idea of sum rules was developed in the pioneering papers of Shifman, Vainshtein, and Zakharov29 . Their idea was to use quark–hadron duality to obtain a prediction of hadron properties from calculations in the quark–gluon theory of QCD. Consider, as an example, the vacuum polarization induced by the electromagnetic current of a bottom quark: i d4 x eiq·x 0 | T {j ? (x), j ν (0)} | 0 = (q ? q ν ? g ?ν q 2 ) Π(Q2 ) , (37)

where j ? = ? b γ ? b, and Q2 = ?q 2 . The Lorentz structure of the correlator follows from current conservation. The function Π(Q2 ) satis?es a once-subtracted dispersion relation d 1 1 ? ds Π(Q2 ) = ImΠ(s) , (38) 2 dQ π (s + Q2 )2 where by the optical theorem ImΠ(s) = 3s 1 σs (e+ e? → b ? b) 2 12πeb 4πα2 (39)

is related to a measurable cross section. Here, eb = ?1/3 is the electric charge of the b-quark, and “b ? b” is a short-hand notation for “hadrons containing b ? b”. In the sum rule analysis one considers the moments of the correlation function, which are given by 1 d n 1 Mn = ? ds s?n?1 ImΠ(s) . (40) Π(Q2 ) 2 = 2 Q =0 n! dQ π In principle, these moments can be extracted from experiment. As long as Q2 is not close to the resonance region, the function Π(Q2 ) can be calculated in perturbative QCD. Since Q2 = 0 is far away from the physical threshold for bottomonium production, this assumption is justi?ed in the present case. Hence, perturbative QCD should provide a good approximation to a calculation of the moments Mn . The leading and next-to-leading perturbative contributions to the correlator are shown in Fig. 5. They give rise to the spectral density29 ImΠpert (s) = 1 v (3 ? v 2 ) 4αs π 2 3 + v Θ(v 2 ) 1 + ? 4π 2 3π 2v 4 12 π2 3 ? 2 4 , (41)

where v = 1 ? 4m2 b /s is the relative velocity of the heavy quarks. In this expression one uses the Euclidean mass for mb in order to minimize the e?ect of radiative corrections [see (8)]. For s ? 4m2 b , the perturbative spectral density leads to the well-known cross section σs (e+ e? → b ? b)pert

s?4m2 b

→

4πα2 αs . Nc e2 b 1+ 3s π

(42)

However, the above form of the spectral density includes threshold e?ects, which become important when s comes closer to the threshold region.

q

Figure 5: Perturbative contributions to the correlator Π(Q2 ). The current operators are represented by circles. In QCD sum rules, one usually adds to the perturbative contributions nonperturbative corrections, which are power suppressed and appear at higher orders in the operator product expansion (OPE) of correlation functions such as Π(Q2 ). In the case at hand, the leading nonperturbative corrections are proportional to the gluon condensate αs G2 ? 0.04 GeV4 , which in sum rules is treated as a phenomenological parameter29 . In the case of the bottomonium system, one is in the fortunate situation that the nonperturbative contributions to the spectral function ?4 are very small, as they are suppressed by the ratio αs G2 /m4 b ? 10 . Their e?ect on the mass of the bottom quark is below 1% and can safely be neglected. The idea is now to construct two equivalent representations for the moments Mn , a “theoretical” one based on expression (41) for the spectral function, and a “phenomenological” one based on a measurement of the e+ e? → b ? b cross section:

th Mn

1 ? π 1 π

∞

ds s?n?1 ImΠpert (s) ,

exp Mn =

4m2 b ∞

2 MΥ(1 s)

ds s?n?1 ImΠexp (s) .

(43)

th The moments Mn are very sensitive functions of mb . By comparing them to a phenomenological expression that uses detailed experimental information in the b ? b vector channel, one can extract an accurate value for the bottom quark mass. Experimentally, six resonances have been identi?ed in this channel. Their masses MR and

13

electronic widths ΓR (e+ e? ) are known rather precisely52 . To evaluate the integrals over the experimental cross section, one then makes the following approximation: ImΠexp (s) ? 3 4α2 e2 b

2 MR ΓR (e+ e? ) δ (s ? MR ) + ImΠpert (s) Θ(s ? s0 ) ,

(44)

R

where the sum is over the narrow resonances Υ(1s) to Υ(6s), and the continuum above some threshold value s0 is approximated by perturbation theory. This can be justi?ed by quark–hadron duality. Inserting this ansatz into (43) and equating the two representations for the moments, one obtains the sum rules 1 π s0 4m2 b ds s?n?1 ImΠpert (s) ? 3 4πα2e2 b

?2n?1 MR ΓR (e+ e? ) . R

(45)

th exp The goal is to ?nd an optimal set of parameters (mb , s0 ) such that Mn ? Mn for many values of n. In practice, the perturbative calculation breaks down for large values of n, since the Coulomb corrections proportional to αs /v in (41) then become too large. Typically, one is limited to values n < 10. As mentioned above, the nonperturbative corrections are very small (below 1%) in the case of the bottomonium system. They become sizeable, however, in the charmonium system, where a similar analysis can be performed. Therefore, mb can be extracted with larger accuracy than mc . Let me then quote the numerical results obtained by various authors following the strategy outlined above. Values for the Euclidean mass of the bottom quark have been obtained by Shifman et al.29 (mEucl = 4.23 ± 0.05 GeV), Guberina et al.53 b (mEucl = 4.19 ± 0.06 GeV), and Reinders54 (mEucl = 4.17 ± 0.02 GeV). The di?erences b b in these values are mainly due to changes in the experimental input numbers. The result of Reinders is the most up-to-date one. Values for the Euclidean mass of the charm quark have been obtained many years ago by Novikov et al.55 (mEucl ? c 29 Eucl 1.25 GeV) and Shifman et al. (mc = 1.26 ± 0.02 GeV). They have not changed since then. Converting these results into pole masses using (8), one ?nds54

mpole = 4.55 ± 0.05 GeV , b pole mc = 1.45 ± 0.05 GeV .

(46)

The increase in the error is due to an additional uncertainty in the value of αs (mQ ). Let me ?nally mention some alternative approaches to extract pole masses from QCD sum rules. Voloshin has proposed to resum the large Coulomb corrections proportional to (αs /v )n to all orders in perturbation theory, using a nonrelativistic approach56 . This allows one to go to higher moments, thereby suppressing the resonance contributions to the sum rule. The disadvantage is that some relativistic corrections are not taken into account. The result is a rather large value of the pole mass of the bottom quark57 , mpole = 4.79 ± 0.03 GeV . b 14 (47)

Another alternative is to obtain heavy quark masses from the study of sum rules for heavy–light bound states such as the B - and B ? -mesons. From such an analysis, Narison obtains58 mpole = 4.56 ± 0.05 GeV . (48) b A more recent analysis using the heavy quark expansion to order 1/mQ gives25,59,60 mb = 4.71 ± 0.07 GeV , mc = 1.37 ± 0.12 GeV .

(49)

It should be noticed, however, that these sum rules are more sensitive to nonperturbative corrections than are the sum rules for the quarkonium systems. 1.6. Infrared Renomalons As emphasized at the beginning of this lecture, the fact that at low energies quarks are always con?ned into hadrons prohibits a physical, on-shell de?nition of quark masses. Although for heavy quarks the notion of a pole mass is widely used, such a concept becomes meaningless beyond perturbation theory. No precise de?nition of the pole mass can be given once nonperturbative e?ects are taken into account. It is interesting that indications for an intrinsic ambiguity in the pole mass can already be found within the context of perturbation theory, when one studies the asymptotic behaviour of the perturbative series for the quark self-energy61,62 . It can be shown that the existence of so-called infrared renomalons63 ? 67 generates a factorial divergence in the expansion coe?cients. Roughly speaking, the reason is that self-energy corrections to the gluon propagator in the√ diagram depicted in Fig. 1 e?ectively introduce the running coupling constant αs ( k 2 ), where k 2 is the √ square of the virtual gluon momentum in Euclidean space. Since αs ( k 2 ) increases as k 2 decreases, small loop momenta become more important. One may reorganize the perturbative expansion as an expansion in the small coupling constant αs (mQ ) using ∞ n √ m2 β0 Q αs ( k 2 ) ? αs (mQ ) αs (mQ ) ln 2 . (50) k n=0 4π When one then performs the loop integral over k 2 , the extra logarithms give rise to a growth of the expansion coe?cients proportional to n!. In such a situation, it is in principle not possible to improve the accuracy of perturbation theory by including more and more terms in the perturbative series. Starting from some order n, the size of the corrections increases. The best that can be achieved is to truncate the series at an optimal value of n. This introduces an intrinsic error, which depends exponentially on 1/αs (mQ ). In the case of the pole mass, one can show that the irreducible uncertainty in the value of mpole is of order61,62 Q = ?mpole Q 8 mQ exp 3β0 ? 15 2π 8ΛQCD ? . β0 αs (mQ ) 3β0 (51)

pole It is thus of the order of the con?nement scale ΛQCD . The fact that ?mpole Q /mQ vanishes in the mQ → ∞ limit justi?es the concept of a pole mass for heavy quarks, at least in an approximate sense that is limited to the context of perturbation theory. Although one has to be careful when using (51) to obtain an estimate of the intrinsic uncertainty in the pole mass, I will insert ΛQCD ? 150 MeV to ?nd ?mpole ? Q 50 MeV. This uncertainty should be kept in mind when considering numerical results for heavy quark masses obtained, for instance, using QCD sum rules.

1.7. Exercises ? Derive the one-loop expression (4) for the on-shell quark self-energy in QCD. The necessary loop integrals in D space-time dimensions can be found in any reasonable textbook on quantum ?eld theory. ? Show that (19) is the exact solution of the RGE (15), and derive the next-toleading logarithmic approximation (20) for the evolution operator. ? By calculating the 1/?-poles of the diagrams shown in Fig. 2, derive the oneloop RGE (28) for the running top quark mass. ? Solve the RGE (30) for the case αs = const., and show that there is an infrared-stable ?xed point at αt = 16 αs . 9 ? Derive eqs. (32)–(35).

16

2. Quark Mixing in the Standard Model In this lecture, I give an introduction to ?avour-changing decays and quark mixing in the standard model. I will discuss quark mixing in the cases of two and three generations, introduce some useful parametrizations of the Kobayashi–Maskawa matrix, and brie?y review the status of the direct experimental determination of the entries in this matrix from tree-level processes. I will then turn to the geometrical interpretation of the mixing matrix provided by the unitarity triangle, and ?nish with some remarks on rare decays. Since you probably have heard many lectures on these subjects, my presentation will be rather short. Let me refer those of you who want more details to two excellent review articles that I like very much: A Top Quark Story by Buras and Harlander68 , and CP Violation by Nir69 . 2.1. Cabibbo–Kobayashi–Maskawa Matrix Let me brie?y remind you of some facts about the ?avour sector of the standard model. Below mass scales of order mW ? 80 GeV, the standard model gauge group SUC (3) × SUL (2) × UY (1) is spontaneously broken to SUC (3) × Uem (1), since the scalar Higgs doublet φ acquires a vacuum expectation value v 0 =√ ; v ? 246 GeV . (52) 2 1 This gives masses to the W - and Z -bosons, as well as to the quarks and leptons. The quark masses arise from the quark Yukawa couplings to the Higgs doublet, which in the unbroken theory are assumed to be of the most general form that is invariant under local gauge transformations. The Yukawa interactions are written in terms of the weak eigenstates q ′ of the quark ?elds, which have de?nite transformation properties under SUL (2) × UY (1). After the symmetry breaking, one rede?nes the quark ?elds so as to obtain the mass terms in the canonical form. This has an interesting e?ect on the form of the ?avour-changing charged-current interactions. In the weak basis, these interactions have the form φ = φ+ φ0 Lint d′L g ? ′ ′ ? ′ ?? ′ ? = ? √ (? + h.c. uL , c ?L , tL ) γ ? sL ? W? 2 ′ bL

? ? ? ?

(53)

In terms of the mass eigenstates q , however, this becomes Lint

The Kobayashi–Maskawa mixing matrix VKM

?

dL g ? ? ?L ) γ ? VKM ? = ? √ (? uL , c ?L , t ? sL ? W? + h.c. 2 b

L

(54)

Vud Vus Vub ? ? ? ? Vcd Vcs Vcb ? Vtd Vts Vtb 17

?

(55)

is a unitary matrix in ?avour space. In the general case of n quark generations, VKM would be an n × n matrix. I will now discuss the structure of this matrix for the cases of two and three generations. 2.1.1. Mixing Matrix for Two Generations In this case, V is a 2 × 2 unitary matrix and can be parametrized by one angle and three phases: V = = cos θC eiα sin θC eiβ ? sin θC eiγ cos θC ei(β +γ ?α) eiα 0 0 eiγ cos θC sin θC ? sin θC cos θC 1 0 0 ei(β ?α) . (56)

The three phases are not observable, however, as they can be absorbed into a rede?nition of the phases of the quark ?elds uL, cL , and sL relative to dL . After this rede?nition, the matrix takes the standard form due to Cabibbo70 : VC = cos θC sin θC ? sin θC cos θC . (57)

The Cabibbo angle θC can be extracted from K → π e+ νe decay. Experimentally, one ?nds52 sin θC ? 0.22. 2.1.2. Mixing Matrix for Three Generations A 3 × 3 unitary matrix can be parametrized by three Euler angles and six phases, ?ve of which can be removed by adjusting the relative phases of the left-handed quark ?elds. Hence, three angles θij and one observable phase δ remain in the quark mixing matrix, as was ?rst pointed out by Kobayashi and Maskawa71 . For completeness, I note that in the general case of n generations, it is easy to show that there are 1 n(n ? 1) angles and 1 (n ? 1)(n ? 2) observable phases72 . Whereas therefore 2 2 the original Cabibbo matrix was real and had only one parameter, the Kobayashi– Maskawa matrix of the standard model is complex and can be parametrized by four parameters. The imaginary part of the mixing matrix is necessary to describe CP violation in the standard model. In general, CP is violated in ?avour-changing decays if there is no degeneracy of any two quark masses, and if the quantity JCP = 0, where

? ? JCP = | Im(Vij Vkl Vil Vkj ) | ;

i = k, j = l

(58)

is invariant under phase rede?nitions of the quark ?elds. One can show73 that all CP-violating amplitudes in the standard model are proportional to JCP . 18

Here, one uses the short-hand notation cij = cos θij and sij = sin θij . Some advantages of this parametrization are the following:

I will now discuss two of the most convenient parametrizations of the mixing matrix. The “standard parametrization” recommended by the Particle Data Group52 is74 ? ? c12 c13 s12 c13 s13 e?iδ ? ? (59) VKM = ? ?s12 c23 ? c12 s23 s13 eiδ c12 c23 ? s12 s23 s13 eiδ s23 c13 ? . s12 s23 ? c12 c23 s13 eiδ ?c12 s23 ? s12 c23 s13 eiδ c23 c13

? | Vub| = s13 is given by a single angle, which experimentally turns out to be very small. ? Because of this, several other entries are given by single angles to an accuracy of better than four digits. They are: Vud ? c12 , Vus ? s12 , Vcb ? s23 , and Vtb ? c23 . ? The CP-violating phase δ appears together with the small parameter s13 , making explicit the fact that CP violation in the standard model is a small e?ect. Indeed, one ?nds JCP = | s13 s23 s12 sδ c2 13 c23 c12 | . (60)

For many purposes and applications, it is more convenient to use an approximate parametrization of the Kobayashi–Maskawa matrix, which makes explicit the strong hierarchy that is observed experimentally. Setting c13 = 1 (experimentally, it is known that c13 > 0.99998) and neglecting s13 compared to terms of order unity, one ?nds ? ? c12 s12 s13 e?iδ ? ?s12 c23 c12 c23 s23 ? (61) VKM ? ? ? . iδ s12 s23 ? c12 c23 s13 e ? c12 s23 c23 Now denote s12 = λ ? 0.22. Experiments indicate that s23 = O (λ2 ) and s13 = O (λ3). Hence, it is natural to de?ne s23 = A λ2 and s13 e?iδ = A λ3 (ρ ? iη ), with A, ρ and η of order unity. An expansion in powers of λ then leads to the Wolfenstein parametrization75

? ? ? ?

VKM ? ?

1?

λ2 2

λ 1?

λ2 2 2

A λ3 (1 ? ρ ? iη ) ? A λ

?λ

A λ3 (ρ ? iη ) A λ2 1

?

? ? ? + O (λ 4 ) . ?

(62)

It nicely exhibits the hierarchy of the mixing matrix: The entries in the diagonal are close to unity, Vus and Vcd are of order 20%, Vcb and Vts are of order 4%, and Vub and Vtd are of order 1% and thus the smallest entries in the matrix. Some care has to be taken when one wants to calculate the quantity JCP in the Wolfenstein 19

which shows that JCP is generically of order 10?4 for λ ? 0.22. As a consequence, CP violation in the standard model is a small e?ect. In principle, the elements in the ?rst two rows of the mixing matrix are accessible in so-called direct (tree-level) processes, i.e. in weak decays of hadrons containing the corresponding quarks. In practice, as I will discuss below, the entries | Vud| and | Vus | are known to an accuracy of better than 1% from such decays, whereas | Vcd |, | Vcs|, and | Vcb| are known to about 10–20%. The element | Vub| has an uncertainty of about a factor of 2. The same is true for | Vtd |, which is obtained from a rare process, ?d mixing. There is no direct information on | Vts | and | Vtb |. The large namely Bd –B uncertainty in | Vub| and | Vtd | translates into an uncertainty of the Wolfenstein parameters ρ and η . A more precise determination of these parameters will be the challenge to experiments and theory over the next decade. 2.2. Experimental Information on VKM from Tree-Level Processes

parametrization, since the result is of order λ6 and thus beyond the accuracy of the approximation. However, taking i = u, j = d, k = t, and l = b in (58), one obtains the correct answer JCP ? A2 η λ6 ? 1.1 × 10?4 A2 η , (63)

In this section, I will brie?y review what is known about entries of the Kobayashi– Maskawa matrix from tree-level processes. Rare decays, which are induced by loopdiagrams involving heavy particles, will be discussed later. 2.2.1. Determination of | Vud | From a comparison of the very accurate measurements of super-allowed nuclear β -decay with ?-decay, including a detailed analysis of radiative corrections76,77 , one obtains52 | Vud| = 0.9744 ± 0.0010 . (64) 2.2.2. Determination of | Vus |

0 An analysis of Ke3 decays, i.e. K + → π 0 e+ νe and KL → π ? e+ νe , including isospin- and SU(3)-breaking e?ects78 , yields | Vus | = 0.2196 ± 0.0023. Alternatively, from an analysis of hyperon decays, including SU(3)-breaking corrections79 , one obtains | Vus | = 0.222 ± 0.003. This leads to the combined value52

| Vus | = 0.2205 ± 0.0018 .

(65)

Let me add a note on Ke3 decays in this context. The fact that the theoretical uncertainty in the description of these exclusive hadronic decays is only about 1% is at ?rst sight surprising. The reasons are the following: The approximate SU(3) ?avour symmetry of QCD implies that, at zero momentum transfer, the hadronic 20

K →π form factor f+ that parametrizes K → π transitions induced by a vector current equals unity, up to symmetry-breaking corrections. The Ademollo–Gatto theorem80 states that these corrections are of second order in the symmetry-breaking parameter (ms ? mq ), where q = u or d. Chiral perturbation theory can be employed to calculate the leading symmetry-breaking corrections (chiral logarithms) in a modelindependent way78 . In my third lecture, I will discuss that there are analogues to all these ingredients in the case of semileptonic B -decays.

2.2.3. Determination of | Vcd| From the combination of data on single charm production in deep-inelastic neutrino–nucleon scattering with the semileptonic branching ratios of charmed mesons, one deduces52 | Vcd| = 0.204 ± 0.017 . (66) 2.2.4. Determination of | Vcs | Here the usefulness of deep-inelastic scattering data is limited, since the extraction of | Vcs | would depend upon an assumption about the strange quark content of the nucleon. It is better to proceed in analogy to the determination of | Vus |, ? e+ νe . The experimental data are i.e. to use the semileptonic De3 decays D → K compared to various model calculations of the decay width81 ? 83 . This procedure leads, however, to a signi?cant amount of model dependence. The result is52 | Vcs | = 1.00 ± 0.20 . (67)

You may wonder why the theoretical description of the decay width is so uncertain. The reason is that, in general, our capabilities to calculate hadronic form factors in QCD are rather limited. An understanding of the nonperturbative con?ning forces would be necessary to make the connection between quark and hadron properties, which is a prerequisite for any calculation of hadronic matrix elements. In the case of Ke3 decays, a symmetry of QCD helps us eliminate most of the hadronic uncertainties. For De3 decays, there is no such ?avour symmetry. Hence, one has to rely on phenomenological models. Let me brie?y mention some of them: ? Isgur, Scora, Grinstein, and Wise (ISGW) have proposed a nonrelativistic constituent quark model83 . They solve the Schr¨ odinger equation for a Coulomb plus linear potential to obtain meson wave functions. Weak decay form factors are calculated at maximum momentum transfer, where both mesons have a common rest frame, by computing overlap integrals of the wave functions of the initial and ?nal meson. The nonrelativistic approximation breaks down if 2 q 2 ? qmax . ? Bauer, Stech, and Wirbel (BSW) have constructed a relativistic model81 , which uses light-cone wave functions to calculate weak decay form factors 21

at zero momentum transfer. In order to extrapolate to q 2 > 0, they assume nearest-pole dominance. ? K¨ orner and Schuler84 (KS) use a variation of the BSW approach, in which the 2 q dependence of the form factors is adjusted according to asymptotic QCD power-counting rules85 . ? There are a variety of other models, for instance by Suzuki86 , Altomari and Wolfenstein87 , and Grinstein et al.82 , which are more or less closely related to one of the above. Disadvantages of these models are that their relation to the underlying theory of QCD is not obvious, that it is hard to obtain an estimate of their intrinsic uncertainty or range of applicability, that they rely on ad hoc assumptions (for instance, concerning the q 2 dependence of form factors), and that it is di?cult to improve them in a systematic way. Another stream of research focuses on analytical or numerical approaches that bear a closer relation to ?eld theory. QCD sum rules29 provide a fully relativistic approach based on QCD. They rely, however, on the assumption of quark–hadron duality and need some phenomenological input. Nonperturbative e?ects are modelled in a simple way by introducing few universal numbers, the so-called vacuum condensates. There have been extensive applications of sum rules to the study of meson weak decay form factors88 ? 96 . These analyses can be done for all physical values of q 2 . However, in spite of the many successes of QCD sum rules it must be said that the only known approach to low-energy QCD that truly starts from ?rst principles is lattice gauge theory. For more details on this, I refer to the lectures by P. Lepage in this volume. Let me just mention that, as far as weak decay form factors are concerned, these computations are extremely complex and CPU-consuming. They are not yet competitive with (less rigorous) analytical approaches. 2.2.5. Determination of | Vub| Since b → u transitions are strongly suppressed in nature, their discovery in inclusive B → Xu ? ν ? decays was one of the breakthroughs in B -physics in recent 3?6 years . In order to account for CP violation in the standard model, it is necessary that Vub be di?erent from zero; otherwise JCP = 0, and the observed CP violation in the kaon system97 could not be explained by the Kobayashi–Maskawa mechanism. The ?rst direct observation of a b → u transition was provided by two fully reconstructed events, one B 0 → π + ?? ν ? decay and one B + → ω ?+ ν decay, reported by the ARGUS collaboration5. However, the number of exclusive charmless B -decays that have been observed so far is too small to obtain a reasonable determination of | Vub|. A high luminosity B -factory would improve this situation. Given enough statistics, the idea is to obtain a model-independent measurement of | Vub| by using heavy quark symmetry to relate the decay form factors in B → Xu ? ν ? 22

and D → Xd ? ν ? transitions98 , at the same value of the recoil energy of the light hadron Xq . In the limit of in?nite heavy quark masses (mb , mc → ∞), the ratio of these form factors is ?xed in a model-independent way. In the real world, however, there are nonperturbative power corrections to this limit. As a consequence, even in an ideal measurement the ratio | Vub/Vcd | can only be determined up to hadronic corrections of order ΛQCD ΛQCD + ... , (68) ? 1+c mc mb with a coe?cient c of order unity. Model calculations99,100 show that these corrections are of order 15% for X = π or ρ. I believe that it is fair to say that even at a high-luminosity B -factory, the prospects for getting a determination of | Vub| from exclusive decays, which is more reliable than that from inclusive decays, are rather limited. Nevertheless, it is certainly worth while to pursue this strategy as an alternative.

0.5 dΓ/dPl [arb. units] 0.4 0.3 0.2 0.1 0 0 0.5 1

B -> Xu l ν

B -> Xc l ν

1.5 Pl [GeV]

2

2.5

Figure 6: Sketch of the lepton spectrum in inclusive semileptonic B -decays. I assume | Vub/Vcb | ? 0.08. The b → u signal has been multiplied by a factor of 5 to be visible on this plot. The present determinations of | Vub| are based on measurements of the lepton momentum spectrum in inclusive B → Xq ? ν ? decays, where Xq is any hadronic state containing a q -quark, with q = c or u. The expected signal is shown in Fig. 6. Over most of the kinematic region, the spectrum is dominated by b → c transitions, which are strongly enhanced with respect to b → u decays. The only place to observe charmless transitions is in a small window above the kinematic limit for B → D?ν ?, but below the kinematic limit for B → π ? ν ? decays: mB m2 1? D 2 m2 B ≤ p? ≤ 23 mB m2 π , 1? 2 2 mB (69)

i.e. 2.31 GeV ≤ p? ≤ 2.64 GeV. Indeed, the ARGUS and CLEO collaborations have reported signals in this region3 ? 6 . An extraction of | Vub| from these data is di?cult, however, as the shape of the spectrum close to the kinematic endpoint is dominated by nonperturbative e?ects. This can be seen as follows. The conventional description of inclusive decays of hadrons containing a heavy quark starts from the free-quark decay model, in which the hadronic decay B → Xq ? ν ? is modelled by the quark decay b → q ? ν ?. One usually calculates the decay distributions to order αs in perturbation theory by including the e?ects of real and virtual gluon emission101 ? 103 . Computing then the lepton spectrum, one ?nds that the kinematic limit for b → u decays is given by pmax = mb /2 ? 2.35 GeV, where I neglect the u-quark mass ? and take mb ? 4.7 GeV for the sake of argument. The point is that the region in which the b → u signal is observed experimentally is almost not populated in the free-quark decay model. Nonperturbative bound-state e?ects, such as the motion of the b-quark inside the B -meson, are responsible for the population of the spectrum beyond the parton model endpoint. The way in which such e?ects are incorporated into phenomenological approaches104,105 is to a large extent model-dependent. Only very recently has there been progress towards a model-independent description of the distributions near the kinematic endpoint in the context of QCD106 ? 110 . It has been shown that the leading nonperturbative e?ects can be parametrized in terms of a universal structure function, which describes the momentum distribution of the b-quark inside the B -meson. The hope is that it will be possible to get an accurate prediction for this function using various theoretical approaches. I believe that, within a year or two, it should be possible to reduce the theoretical uncertainty in the analysis of the inclusive decay spectrum by a factor of 2 or so. Alternatively, one can describe the lepton spectrum close to the endpoint as the sum of contributions from semileptonic B -decays into exclusive ?nal states. These exclusive decays are then described using one of the various bound-state models81 ? 84 discussed above. Of course, there is again a strong model dependence associated with this approach. Because of these theoretical uncertainties, the current value of | Vub | has a model dependence of almost a factor of 2. In 1991, the ARGUS collaboration obtained5 | Vub/Vcb | = 0.11 ± 0.012, using the approach of Altarelli et al.104 , which is based on the parton model. Using instead bound-state models to describe the spectrum by a sum over several exclusive decay channels, values for | Vub/Vcb | between 0.10 and 0.22 are obtained. However, the most recent data reported by the CLEO collaboration lead to signi?cantly smaller values6 : | Vub/Vcb | = 0.076 ± 0.008 using the approach of Altarelli et al.104 , and 0.05 < | Vub/Vcb | < 0.11 using bound-state models. Thus, I think a reasonable value to quote is Vub = 0.08 ± 0.03 . Vcb (70)

24

2.2.6. Determination of | Vcb| As for | Vcb|, it can be extracted from both inclusive and exclusive B -decays. I start with a discussion of inclusive decays. The idea is to compare the total semileptonic branching ratio to the parton model formula Br(B → Xq ? ν ?) =

5 G2 m2 F mb c 2 + ηu | Vub|2 , τ η | V | f B c cb 2 3 192π mb

(71)

where ηc and ηu contain the short-distance QCD corrections, and f (x) = 1 ? 8x + 8x3 ? x4 ? 12x2 ln x (72)

is a phase-space correction due to the mass of the charm quark. The QCD correction factor for b → c decays is101,102 ηc ? 0.87. The contribution from b → u transitions is very small and, at the present level of accuracy, can be neglected. Let me note that it has recently become clear how to compute in a systematic way the nonperturbative bound-state corrections to the parton model result (71). For the total inclusive decay rates, these corrections turn out to be of order111 ? 119 (ΛQCD /mb )2 . Numerically, they are very small and play only a minor role in the analysis120 . The disadvantage of using inclusive decays to extract | Vcb| is that there are substantial theoretical uncertainties due to the fact that the total rate scales as the ?fth power of the bottom quark mass. For instance, an uncertainty of only ±150 MeV in the value of mb translates into an uncertainty of ±8% in the value of | Vcb|. For this reason, the last time the particle data group121 used inclusive decays to obtain a value for | Vcb| was in 1988. A compilation of more recent analyses of inclusive decay spectra has been given by Stone122 . He ?nds | Vcb| τB 1.5 ps

1/2

= 0.040 ± 0.005 ,

(73)

where I have normalized the result to a value of the B -meson lifetime that is close to the world average τB0 = (1.48 ± 0.10) ps reported last year by Danilov123 , as well as to the most recent value τB0 = (1.52 ± 0.13) ps obtained from an average of LEP results124 . The error quoted in (73) does not take into account the theoretical uncertainty in the value of mb . Ultimately, a more precise determination must come from the analysis of exclusive decay modes. With the development of the heavy quark expansion, it has become clear that certain symmetries of QCD in the heavy quark limit can help to remove much of the model dependence from the theoretical description of exclusive semileptonic B -decays. This will be explained in detail in my third lecture. In particular, for the decay mode B → D ? ? ν ? the theoretical uncertainties in the momentum distribution of the D ? -meson close to the zero recoil limit can be controlled

25

at the level of a few per cent125 . The most recent analysis of this spectrum by the CLEO collaboration leads to the value126 | Vcb| τB 1.5 ps

1/2

= 0.039 ± 0.006 ,

(74)

which I will use in the analysis below. In contrast to (73), here the uncertainty is dominated by the experimental systematic errors and, to a lesser extent, the statistical ones. The theoretical uncertainty is well below 10%. There is thus room for a signi?cant improvement of the accuracy within the next few years, as more data become available. Finally, let me mention that by combining the above result with the value of | Vus | given in (65), one obtains A= for the Wolfenstein parameter A. 2.3. Unitarity Constraints The fact that the quark mixing matrix must be unitary can be employed to put limits on the elements that are not easily accessible to experiments. One has to assume, however, that there are only three quark generations. Unitarity then implies that | Vtj |2 = 1 ? | Vuj |2 ? | Vcj |2 . (76) 0.003 < | Vtd | < 0.016 , 0.032 < | Vts | < 0.048 , 0.9991 < | Vtb | < 0.9994 . Vcb = 0.80 ± 0.12 2 Vus (75)

Inserting here the known values of | Vuj | and | Vcj |, one obtains

(77)

Note, in particular, that | Vtb | is constrained to an accuracy of better than 10?3 . Unitarity also imposes a rather tight constraint on | Vcs|. Using that | Vcs |2 = 1 ? | Vcd|2 ? | Vcb|2 , 0.976 < | Vcs | < 0.981 . (78) (79) one ?nds This should be compared to the value (67) extracted from a direct measurement. 2.4. The Unitarity Triangle A simple but beautiful way to visualize the implications of unitarity is provided by the so-called unitarity triangle127 , which uses the fact that the unitarity equation

? Vij Vik =0

(j = k )

(80)

26

can be represented as the equation of a closed triangle in the complex plane. There are six such triangles, all of which have the same area128 |A? | = 1 JCP . 2 (81)

Under phase reparametrizations of the quark ?elds, the triangles change their orientation in the complex plane, but their shape remains una?ected. Most useful from the phenomenological point of view is the triangle relation

? ? ? Vud Vub + Vcd Vcb + Vtd Vtb = 0,

(82)

since it contains the most poorly known entries in the Kobayashi–Maskawa matrix. It has been widely discussed in the literature127 ? 134 . In the standard parametriza? tion, Vcd Vcb is real, and the unitarity triangle has the form shown in Fig. 7a. It is ? useful to rescale the triangle by dividing all sides by Vcd Vcb . The rescaled triangle is shown in Fig. 7b. It has the coordinates (0, 0), (1, 0), and (ρ, η ), where ρ and η appear in the Wolfenstein parametrization (62). Unitarity amounts to the statement that the triangle is closed, and CP is violated when the area of the triangle does not vanish, i.e. when all the angles are not zero.

( ; )

Vud Vub

Vtd Vtb

Rb

Rt

Vcd Vcb

(a)

(0; 0) (b)

(1; 0)

Figure 7: (a) The unitarity triangle; (b) its rescaled form in the ρ–η plane. The angle γ coincides with the phase δ of the standard parametrization (59). To determine the shape of the triangle, one can aim for measurements of the two sides Rb and Rt , and of the angles α, β , and γ . So far, experimental information is available only on the sides of the triangle. The current value of | Vub | in (70) implies 0.24 < Rb < 0.53 . (83)

As mentioned above, the uncertainty in this number is mainly a theoretical one, so any signi?cant improvement must originate from new developments in the theory of inclusive B -decays. 27

Let me now discuss the determination of | Vtd | and of the side Rt of the unitarity ?d -mesons, triangle. The process of interest here is the mixing between Bd - and B which in the standard model is a rare process mediated through box diagrams with a virtual top quark exchange. The ARGUS and CLEO collaborations have reported measurements135,136 of the mixing parameter xd de?ned as xd = ?mBd τBd , where ?d system. ?mB denotes the mass di?erence between the mass eigenstates in the Bd –B The weighted average of their results is xd = 0.68 ± 0.08. Combining this with the Bd -meson lifetime123 , τBd = (1.48 ± 0.10) ps, one obtains ?mB = (0.46 ± 0.06) ps?1 . Recently, direct measurements of ?mB have been reported by some of the LEP experiments. The average value presented at the Moriond Conference is137 ?mB = .063 ?1 (0.519+0 ?0.059 ) ps . I combine these results to obtain ?mB = (0.49 ± 0.04) ps?1 . The theoretical expression for ?mB in the standard model is ?mB =

2 G2 F mW 2 ? 2 ηQCD mB (BB fB ) S (mt /mW ) | Vtd Vtb | , 6π 2

(84)

(85)

where ηQCD ? 0.55 contains the next-to-leading order QCD corrections138 , and S (mt /mW ) is a function of the top quark mass139,140 . For 100 GeV < mt < 300 GeV, 2 it can be approximated by141 S (mt /mW ) ? 0.784 (mt/mW )1.52 . The product (BB fB ) parametrizes the hadronic matrix element of a local four-quark operator between B -meson states. Recently, there has been a lot of activity and improvement in theoretical calculations of the B -meson decay constant fB , using both lattice gauge theory and QCD sum rule calculations. For a summary and discussion of the extensive literature on this subject, I refer to the review articles by Buras and Harlander68 and by myself25 . Here I shall use the range BB fB = (200 ± 40) MeV , which covers most theoretical predictions. Solving (85) for | Vtd |, one obtains | Vtd | = 0.0088 200 MeV BB fB

1/2 0.76 1/2 1/2

(86)

170 GeV mt

?mB 0.5 ps?1

.

(87)

Using the numbers given above, as well as mt = (170 ± 30) GeV, I ?nd | Vtd | = 0.0087 ± 0.0023 . The corresponding range of values for Rt is 0.75 < Rt < 1.44 . (89) (88)

The error in this value will hopefully be reduced soon, with a measurement of the top quark mass at Fermilab. Further improvements could come from a measurement 28

?s mixing. In the latter of fB in B + → τ + ντ decays, or from the observation of Bs –B case, one could use the relation Vts xs = xd Vtd

2

× 1 + SU(3)-breaking corrections

(90)

to obtain a rather clean measurement of | Vtd |. However, in the standard model one expects large values xs ? 10–30, which will be very hard to observe in near-future experiments.

1.2 1 0.8 0.6 0.4 0.2 0 -0.6 -0.4 -0.2 0 0.2 ρ 0.4 0.6 0.8 1 η

Figure 8: Experimental constraints on the unitarity triangle in the ρ–η plane. The region between the solid (dashed) circles is allowed by the measurement of Rb (Rt ) discussed above. The dotted curves show the constraint following from the measurement of the ε parameter in the kaon system. The shaded region shows the allowed range for the tip of the unitarity triangle. The base of the triangle has the coordinates (0, 0) and (1, 0). In Fig. 8, I show the constraints imposed by the above results for Rb and Rt in the ρ–η plane. Another constraint can be obtained from the study of CP violation in K -decays. 2.5. CP Violation in K -Decays ? The experimental result on the parameter ε that measures CP violation in K –K mixing implies that the unitarity triangle lies in the upper half plane, provided the so-called BK parameter is positive. Most theoretical predictions142 ? 150 fall in the range BK = 0.65 ± 0.20, supporting this assertion. The arising constraint in the ρ–η 29

plane has the form of a hyperbola approximately given by68 η (1 ? ρ) A2 mt mW

1.52

+ Pc A2 BK = 0.50 ,

(91)

where Pc ? 0.66 corresponds to the contributions from box diagrams containing charm quarks, and A = 0.80 ± 0.12 according to (75). I have included this bound in Fig. 8. In principle, the measurement of the ratio Re(ε′ /ε) in the kaon system151,152 could provide a determination of η independent of ρ. In practice, however, the theoretical calculations153 ? 155 of this ratio are a?ected by large uncertainties, so that there currently is no useful bound to be derived. Further information on the parameters ρ and η could be extracted from rare K -decays, such as KL → π 0 ν ν ?, K + → π+ν ν ?, and KL → ?+ ?? . The corresponding branching ratios vary between 10?11 and 10?10 . The experimental detection of such decays will be very hard. 2.6. Potential Improvements in the Determination of ρ and η The main conclusion to be drawn from Fig. 8 is that, given the present theoretical ? mixing, and experimental uncertainties in the analysis of charmless B -decays, B –B and CP violation in the kaon system, there is still a large region in parameter space that is allowed for the Wolfenstein parameters ρ and η . This has important implications. For instance, the allowed region for the angle β of the unitarity triangle (see Fig. 7) is 6.9? < β < 31.8? , corresponding to 0.24 < sin 2β < 0.90 . (92)

Below I will discuss that the CP asymmetry in the decay B → ψ KS , which is one of the favoured modes to search for CP violation at a future B -factory, is proportional to sin 2β . Obviously, the prospects for discovering CP violation with such a machine crucially depend on whether sin 2β is closer to the upper or lower bound in (92). A more reliable determination of the shape of the unitarity triangle is thus of the utmost importance. Let me brie?y discuss what I think are the potential improvements in this determination in the near future. It is not unrealistic to assume the following: ? The top quark will be discovered at Fermilab. ? The uncertainty in the value of | Vcb| will be reduced to ±0.04, mainly through a better control of the experimental systematic errors. ? Both theoretical and experimental improvements will reduce the uncertainty in the determination of | Vub | by a factor of 2. ? Improved theoretical calculations, in particular using lattice gauge theory, 1/2 will allow the determination of the nonperturbative parameters BB fB and BK with an accuracy of 10%. 30

1.2 1 0.8 0.6 0.4 0.2 0 -0.6 -0.4 -0.2 0 0.2 ρ 0.4 0.6 0.8 1 η

Figure 9: Constraints on the unitarity triangle obtained assuming a more precise determination of the various input parameters, as explained in the text. For the purpose of illustration, I shall assume that mt = 170 GeV, | Vcb| = 0.039 ± 1/2 0.004, | Vub/Vcb | = 0.080 ± 0.015, BB fB = (200 ± 20) MeV, and BK = 0.65 ± 0.07. This leads to Rb = 0.38 ± 0.08 and Rt = 1.09 ± 0.20. Similarly, the constraint derived from the measurement of ε becomes more restrictive. This ?ctitious scenario is illustrated in Fig. 9. Obviously, the parameters of the triangle (the value of η , in particular) would be much more constrained than they are at present. Notice also that once a su?cient accuracy is obtained, there is a potential for tests of the standard model. If the three bands in Fig. 9 did not overlap, this would be an indication of new physics. 2.7. Rare B -Decays Let me brie?y touch upon the interesting subject of rare B -decays, i.e. decays induced by loop diagrams. There are nice review articles on this subject by Bertolini156 and Ali157 . Rare B -decays are dominated by short-distance physics, namely the top quark contribution in penguin and box diagrams. The decay rates are usually rather sensitive to the mass of the top quark141 . Long-distance e?ects are much less important than in rare K -decays. Hence, rare B -decays are “clean” from the theoretical point of view. They are ideal to determine the parameters of the unitarity triangle. In addition, these decays can provide important tests of the standard model, since they are often sensitive probes of new physics. So far, the only rare B -decay that has been observed is7 B → K ? γ . In general, decays of the type B → Xs γ are mediated by penguin diagrams such as the ones shown in Fig. 10. These decays receive very large, but calculable, short-distance 31

Of more interest from the experimental point of view are exclusive rare decay modes. Since B → K γ is forbidden by angular momentum conservation, the lightest ?nal state that can appear is K ? γ . To obtain from (93) a prediction for Br(B → K ? γ ) requires an estimate of the ratio RK ? = Γ(B → K ? γ ) . Γ(B → Xs γ )

QCD corrections158 ? 164 , which strongly enhance the decay rate. In addition, it is important to take into account the e?ect of gluon bremsstrahlung, which a?ects the total decay rate as well as the photon spectrum165 . All such e?ects are reasonably well under control for the inclusive B → Xs γ decay rate. For mt = 170 ± 30 GeV, the prediction is68 Br(B → Xs γ ) = (2.8–4.2) × 10?4 . (93)

(94)

This is where hadronic uncertainties enter the analysis. Estimates of RK ? presented in the literature vary between 4% and 40%. I believe that a reasonable value is165 ? 168 RK ? ? 10–15%. This leads to Br(B → K ? γ ) ? (3–6) × 10?5 . This is in agreement with the result reported by the CLEO collaboration7: Br(B → K ? γ ) = (4.5 ± 1.9 ± 0.9) × 10?5 . (96) (95)

W

b

t

t

s

b

s

W

Figure 10: Penguin diagrams contributing to the rare decay B → Xs γ . Further rare B -decays, for which there are “clean” theoretical predictions (which depend, however, on the top quark mass), include B → Xs ν ν ? with a branching ratio ? with a branching ratio of about 10?5 , and Bs → ? ? ? with of 10?5 –10?4, B → Xs ? ? a branching ratio of about 10?7 . 2.8. Measurements of the Angles of the Unitarity Triangle Let me ?nish this short introduction into the phenomenology of rare decays with a very beautiful application, namely the direct measurement of CP violation 32

in decays of neutral B -mesons into CP eigenstates f . For these processes, there are often no (or only very small) hadronic uncertainties. The CP asymmetries arise from an interference of mixing and decay. When such a decay is dominated by a single weak decay amplitude, the time-integrated asymmetry obeys a very simple relation to one of the angles ? of the unitarity triangle68,69,129 :

∞

0 ∞ 0

? 0 (t) → f ) dt Γ(B 0 (t) → f ) ? Γ(B Γ(B 0 (t) ? 0 (t) → f ) → f ) + Γ(B

dt

= ± sin(2?)

x , 1 + x2

(97)

where x = ?mB τB . Examples of such decays are: Bd Bd Bd Bd Bs Bs → ψ KS → φ KS → D+D? → π+π? → ρ KS → φ KS (? sin 2β ) , (? sin 2β ) , (? sin 2β ) , (sin 2α) , (? sin 2γ ) , (sin 2β ) .

(98)

The decay mode Bd → ψ KS is particularly “clean” and is often referred to as the “gold-plated” mode of a future B -factory. 2.9. Exercises ? Show that the generalization of the Kobayashi–Maskawa matrix to n quark 1 generations can be parametrized by 2 n(n ? 1) angles and 1 (n ? 1)(n ? 2) 2 observable phases. ? Show that the quantity JCP de?ned in (58) is invariant under phase rede?nitions of the quark ?elds. Verify that in the standard parametrization JCP is given by (60). ? Derive that the maximum lepton momentum in the decay B → X ? ν ? is given 1 2 2 by 2 mB (1 ? mX /mB ) [cf. (69)]. ? Convince yourself that all six unitarity triangles have the same area (81). ? Show that the decay B → K γ is forbidden by angular momentum conservation.

33

3. Heavy Quark Symmetry In this lecture, I give an introduction to one of the most active and exciting recent developments in theoretical particle physics: the discovery of heavy quark symmetry and the development of the heavy quark e?ective ?eld theory. For a more detailed description, I refer to my recent review article in Physics Reports 25 . Earlier introductions into this ?eld were given by Georgi169 , Grinstein170 , Isgur and Wise171 , and Mannel172 . 3.1. The Physical Picture 8 ? 13 There are several reasons why the strong interactions of systems containing heavy quarks are easier to understand than those of systems containing only light quarks. The ?rst is asymptotic freedom, the fact that the e?ective coupling constant of QCD becomes weak in processes with large momentum transfer, corresponding to interactions at short-distance scales38,39 . At large distances, on the other hand, the coupling becomes strong, leading to nonperturbative phenomena such as the con?nement of quarks and gluons on a length scale Rhad ? 1/ΛQCD ? 1 fm, which determines the typical size of hadrons. Roughly speaking, ΛQCD ? 0.2 GeV is the energy scale that separates the regions of large and small coupling constant. When the mass of a quark Q is much larger than this scale, mQ ? ΛQCD , Q is called a heavy quark. The quarks of the standard model fall naturally into two classes: up, down and strange are light quarks, whereas charm, bottom and top are heavy quarks? . For heavy quarks, the e?ective coupling constant αs (mQ ) is small, implying that on length scales comparable to the Compton wavelength λQ ? 1/mQ the strong interactions are perturbative and much like the electromagnetic interactions. In ? ), whose size is of order λQ /αs (mQ ) ? Rhad , fact, the quarkonium systems (QQ are very much hydrogen-like. Since the discovery of asymptotic freedom, their properties could be predicted173 before the observation of charmonium, and later of bottomonium states. Things are more complicated for systems composed of a heavy quark and other light constituents. The size of such systems is determined by Rhad , and the typical momenta exchanged between the heavy and light constituents are of order ΛQCD . The heavy quark is surrounded by a most complicated, strongly interacting cloud of light quarks, antiquarks, and gluons. This cloud is sometimes referred to as the “brown muck”, a term invented by Isgur to emphasize that the properties of such systems cannot be calculated from ?rst principles (at least not in a perturbative way). In this case, it is the fact that mQ ? ΛQCD , or better λQ ? Rhad , i.e. the fact that the Compton wavelength of the heavy quark is much smaller than the size of the hadron, which leads to simpli?cations. To resolve the quantum numbers of the heavy quark would require a hard probe. The soft gluons which couple to the

?

Ironically, the top quark is of no relevance to my discussion here, since it is too heavy to form hadronic bound states before it decays.

34

“brown muck” can only resolve distances much larger than λQ . Therefore, the light degrees of freedom are blind to the ?avour (mass) and spin orientation of the heavy quark. They only experience its colour ?eld, which extends over large distances because of con?nement. In the rest frame of the heavy quark, it is in fact only the electric colour ?eld that is important; relativistic e?ects such as colour magnetism vanish as mQ → ∞. Since the heavy quark spin participates in interactions only through such relativistic e?ects, it decouples. That the heavy quark mass becomes irrelevant can be seen as follows: As mQ → ∞, the heavy quark and the hadron that contains it have the same velocity. In the rest frame of the hadron, the heavy quark is at rest, too. The wave function of the “brown muck” follows from a solution of the ?eld equations of QCD subject to the boundary condition of a static triplet source of colour at the location of the heavy quark. This boundary condition is independent of mQ , and so is the solution for the con?guration of the light degrees of freedom. It follows that, in the mQ → ∞ limit, hadronic systems which di?er only in the ?avour or spin quantum numbers of the heavy quark have the same con?guration of their light degrees of freedom. Although this observation still does not allow us to calculate what this con?guration is, it provides relations between the properties of such particles as the heavy mesons B , D , B ? and D ? , or the heavy baryons Λb and Λc (to the extent that corrections to the in?nite quark mass limit are small in these systems). Isgur and Wise realized that these relations result from new symmetries of the e?ective strong interactions of heavy quarks at low energies13 . The con?guration of light degrees of freedom in a hadron containing a single heavy quark Q(v, s) with velocity v and spin s does not change if this quark is replaced by another heavy quark Q′ (v, s′ ) with di?erent ?avour or spin, but with the same velocity. Both heavy quarks lead to the same static colour ?eld. It is not necessary that the heavy quarks Q and Q′ have similar masses. What is important is that their masses are large compared to ΛQCD . For Nh heavy quark ?avours, there is thus an SU(2Nh ) spin-?avour symmetry group. These symmetries are in close correspondence to familiar properties of atoms: The ?avour symmetry is analogous to the fact that di?erent isotopes have the same chemistry, since to good approximation the wave function of the electrons is independent of the mass of the nucleus. The electrons only see the total nuclear charge. The spin symmetry is analogous to the fact that the hyper?ne levels in atoms are nearly degenerate. The nuclear spin decouples in the limit me /mN → 0. The heavy quark symmetry is an approximate symmetry, and corrections arise since the quark masses are not in?nite. These corrections are of order ΛQCD /mQ . The condition mQ ? ΛQCD is necessary and su?cient for a system containing a heavy quark to be close to the symmetry limit. In many respects, heavy quark symmetry is complementary to chiral symmetry, which arises in the opposite limit of small quark masses, mq ? ΛQCD . There is an important distinction, however. Whereas chiral symmetry is a symmetry of the QCD Lagrangian in the limit of vanishing quark masses, heavy quark symmetry is not a symmetry of the Lagrangian (not even an approximate one), but rather a symmetry of an e?ective theory which 35

is a good approximation to QCD in a certain kinematic region. It is realized only in systems in which a heavy quark interacts predominantly by the exchange of soft gluons. In such systems the heavy quark is almost on-shell; its momentum ?uctuates around the mass shell by an amount of order ΛQCD . The corresponding changes in the velocity of the heavy quark vanish as ΛQCD /mQ → 0. The velocity becomes a conserved quantity and is no longer a dynamical degree of freedom21 . Results derived based on heavy quark symmetry are model-independent consequences of QCD in a well-de?ned limit. The symmetry-breaking corrections can, at least in principle, be studied in a systematic way. To this end, it is however necessary to recast the QCD Lagrangian for a heavy quark, ? (i D LQ = Q / ? mQ ) Q , into a form suitable for taking the limit mQ → ∞. 3.2. Heavy Quark E?ective Theory 14 ? 24 In particle physics, it is often the case that the e?ects of a very heavy particle become irrelevant at low energies. It is then useful to construct a low-energy e?ective theory, in which this heavy particle no longer appears. Eventually, this e?ective theory will be easier to deal with than the full theory. A familiar example is Fermi’s theory of the weak interactions. For the description of weak decays of hadrons, one can approximate the weak interactions by point-like four-fermion couplings governed by a dimensionful coupling constant GF . Only at energies much larger than the masses of hadrons can one resolve the structure of the intermediate vector bosons W and Z . Technically, the process of removing the degrees of freedom of a heavy particle involves the following steps174 ? 176 : One ?rst identi?es the heavy particle ?elds and “integrates them out” in the generating functional of the Green functions of the theory. This is possible since at low energies the heavy particle does not appear as an external state. However, although the action of the full theory is usually a local one, what results after this ?rst step is a nonlocal e?ective action. The nonlocality is related to the fact that in the full theory the heavy particle (with mass M ) can appear in virtual processes and propagate over a short but ?nite distance ?x ? 1/M . Thus a second step is required to get to a local e?ective Lagrangian: The nonlocal e?ective action is rewritten as an in?nite series of local terms in an operator product expansion (OPE)177,178 . Roughly speaking, this corresponds to an expansion in powers of 1/M . It is in this step that the short- and long-distance physics is disentangled. The long-distance physics corresponds to interactions at low energies and is the same in the full and the e?ective theory. But short-distance e?ects arising from quantum corrections involving large virtual momenta (of order M ) are not reproduced in the e?ective theory, once the heavy particle has been integrated out. In a third step, they have to be added in a perturbative way using renormalization-group techniques. This procedure is called “matching”. It leads to 36 (99)

renormalizations of the coe?cients of the local operators in the e?ective Lagrangian. An example is the e?ective Lagrangian for nonleptonic weak decays, in which radiative corrections from hard gluons with virtual momenta in the range between mW and some renormalization scale ? ? 1 GeV give rise to Wilson coe?cients, which renormalize the local four-fermion interactions179 ? 181 . The heavy quark e?ective theory (HQET) is constructed to provide a simpli?ed description of processes where a heavy quark interacts with light degrees of freedom by the exchange of soft gluons. Clearly, mQ is the high-energy scale in this case, and ΛQCD is the scale of the hadronic physics one is interested in. However, a subtlety arises since one wants to describe the properties and decays of hadrons which contain a heavy quark. Hence, it is not possible to remove the heavy quark completely from the e?ective theory. What is possible, however, is to integrate out the “small components” in the full heavy quark spinor, which describe ?uctuations around the mass shell. The starting point in the construction of the low-energy e?ective theory is the observation that a very heavy quark bound inside a hadron moves with essentially the hadron’s velocity v and is almost on-shell. Its momentum can be written as

? ? p? Q = mQ v + k ,

(100)

where the components of the so-called residual momentum k ? are much smaller than mQ . Interactions of the heavy quark with light degrees of freedom change the residual momentum by an amount of order ?k ? ? ΛQCD , but the corresponding changes in the heavy quark velocity vanish as ΛQCD /mQ → 0. In this situation, it is appropriate to introduce “large” and “small” component ?elds hv and Hv by hv (x) = eimQ v·x P+ Q(x) , Hv (x) = eimQ v·x P? Q(x) , where P+ and P? are projection operators de?ned as P± = It follows that Q(x) = e?imQ v·x hv (x) + Hv (x) . (103) Because of the projection operators, the new ?elds satisfy v / hv = hv and v / Hv = ?Hv . In the rest frame, hv corresponds to the upper two components of Q, while Hv corresponds to the lower ones. Whereas hv annihilates a heavy quark with velocity v , Hv creates a heavy antiquark with velocity v . If the heavy quark was on-shell, the ?eld Hv would be absent. In terms of the new ?elds, the QCD Lagrangian (99) for a heavy quark takes the following form: ? v iv · D hv ? H ?v i D ? v (iv · D + 2mQ ) Hv + h ?v i D LQ = h / ⊥ Hv + H / ⊥ hv , 37 (104) 1±/ v . 2 (102) (101)

? where D⊥ = D ? ? v ? v·D is orthogonal to the heavy quark velocity: v·D⊥ = 0. In the ? rest frame, D⊥ = (0, D ) contains the spatial components of the covariant derivative. From (104), it is apparent that hv describes massless degrees of freedom, whereas Hv corresponds to ?uctuations with twice the heavy quark mass. These are the heavy degrees of freedom, which will be eliminated in the construction of the e?ective theory. The ?elds are mixed by the presence of the third and fourth terms, which describe pair creation or annihilation of heavy quarks and antiquarks. As shown in the ?rst diagram in Fig. 11, in a virtual process a heavy quark propagating forward in time can turn into a virtual antiquark propagating backward in time, and then ? is turn back into a quark. The energy of the intermediate quantum state hhH larger than the energy of the incoming heavy quark by at least 2mQ . Because of this large energy gap, the virtual quantum ?uctuation can only propagate over a short distance ?x ? 1/mQ . On hadronic scales set by Rhad = 1/ΛQCD , the process essentially looks like a local interaction of the form 1 ?v i D iD / ⊥ hv , (105) h /⊥ 2mQ

where I have simply replaced the propagator for Hv by 1/2mQ . A more correct treatment is to integrate out the small component ?eld Hv , thereby deriving a nonlocal e?ective action for the large component ?eld hv , which can then be expanded in terms of local operators. Before doing this, let me mention a second type of virtual corrections that involve pair creation, namely heavy quark loops. An example is depicted in the second diagram in Fig. 11. Heavy quark loops cannot be described in terms of the e?ective ?elds hv and Hv , since the quark velocities in the loop are not conserved, and are in no way related to hadron velocities. However, such shortdistance processes are proportional to the small coupling constant αs (mQ ) and can be calculated in perturbation theory. They lead to corrections which are added onto the low-energy e?ective theory in the matching procedure.

Figure 11: Virtual ?uctuations involving pair creation of heavy quarks. Time ?ows to the right. On a classical level, the heavy degrees of freedom represented by Hv can be eliminated using the equations of motion of QCD. Substituting (103) into (i D / ? mQ ) Q = 0 gives iD / hv + (i D / ? 2mQ ) Hv = 0 , (106) and multiplying this by P± one derives the two equations ? iv · D hv = i D / ⊥ Hv , (iv · D + 2mQ ) Hv = i D / ⊥ hv . 38 (107)

The second can be solved to give (η → +0) Hv = 1 iD / ⊥ hv , (2mQ + iv · D ? iη ) (108)

which shows that the small component ?eld Hv is indeed of order 1/mQ . One can now insert this solution into the ?rst equation to obtain the equation of motion for hv . It is easy to see that this equation follows from the nonlocal e?ective Lagrangian ? v iv · D hv + h ?v i D Le? = h /⊥ 1 iD / ⊥ hv . (2mQ + iv · D ? iη ) (109)

Clearly, the second term precisely corresponds to the ?rst class of virtual processes depicted in Fig. 11. Mannel, Roberts, and Ryzak have derived this Lagrangian in a more elegant way by manipulating the generating functional for QCD Green’s functions containing heavy quark ?elds182 . They start from the ?eld rede?nition (103) and couple the large component ?elds hv to external sources ρv . Green’s functions with an arbitrary number of hv ?elds can be constructed by taking derivatives with respect to ρv . No sources are needed for the heavy degrees of freedom represented by Hv . The functional integral over these ?elds is Gaussian and can be performed explicitly, leading to the nonlocal e?ective action Se? = d4 x Le? ? i ln ? , (110)

with Le? as given in (109). The appearance of the logarithm of the determinant ? = exp 1 Tr ln [2mQ + iv · D ? iη ] 2 (111)

is a quantum e?ect not present in the classical derivation given above. However, in this case the determinant can be regulated in a gauge-invariant way, and by choosing the axial gauge v · A = 0 one shows that ln ? is just an irrelevant constant182,183 . Because of the phase factor in (103), the x dependence of the e?ective heavy quark ?eld hv is weak. In momentum space, derivatives acting on hv produce powers of the residual momentum k , which is much smaller than mQ . Hence, the nonlocal e?ective Lagrangian (109) allows for a derivative expansion in powers of iD/mQ . Taking into account that hv contains a P+ projection operator, and using the identity gs P+ i D /⊥ iD / ⊥ P+ = P+ (iD⊥ )2 + σαβ Gαβ P+ , (112) 2 where [iD α , iD β ] = igs Gαβ is the gluon ?eld-strength tensor, one ?nds that19,23 ? v iv · D hv + Le? = h gs ? 1 ? hv (iD⊥ )2 hv + hv σαβ Gαβ hv + O (1/m2 Q) . 2mQ 4mQ 39 (113)

Figure 12: Feynman rules of the heavy quark e?ective theory (i, j and a are colour indices). It has become standard to represent a heavy quark by a double line. The residual momentum k is de?ned in (100). In the limit mQ → ∞, only the ?rst terms remains: ? v iv · D hv . L∞ = h (114) This is the e?ective Lagrangian of HQET. From it are derived the Feynman rules depicted in Fig. 12. Let me take a moment to study the symmetries of this Lagrangian21 . Since there appear no Dirac matrices, interactions of the heavy quark with gluons leave its spin unchanged. Associated with this is an SU(2) symmetry group, under which L∞ is invariant. The action of this symmetry on the heavy quark ?elds becomes most transparent in the rest frame, where the generators S i of spin SU(2) can be chosen as Si = 1 2 σi 0 0 σi = 1 γ5 γ 0 γ i , 2 (115)

where σ i are the Pauli matrices. In a general frame, one de?nes a set of three orthonormal vectors ei orthogonal to v , and takes the generators of the spin symmetry as S i = 1 γ / v/ ei . These matrices satisfy the commutation relations of SU(2) 2 5 and commute with v /: [S i , S j ] = i?ijk S k , [/ v, S i] = 0 . (116)

and preserves the on-shell condition v / hv = hv , since

An in?nitesimal SU(2) transformation hv → (1 + i? · S ) hv leaves the Lagrangian invariant, ? v [iv · D, i? · S ] hv = 0 , δ L∞ = h (117) v (1 + i? · S ) hv = (1 + i? · S ) / / v hv = (1 + i? · S ) hv . (118)

There is another symmetry of HQET, which arises since the mass of the heavy quark does not appear in the e?ective Lagrangian. When there are Nh heavy quarks 40

moving at the same velocity, one can simply extend (114) by writing

Nh

L∞ =

i=1

? i iv · D hi . h v v

(119)

This is clearly invariant under rotations in ?avour space. When combined with the spin symmetry, the symmetry group becomes promoted to SU(2Nh ). This is the heavy quark spin-?avour symmetry13,21 . Its physical content is that, in the mQ → ∞ limit, the strong interactions of a heavy quark become independent of its mass and spin. Let me now have a look at the operators appearing at order 1/mQ in (113). They are easiest to identify in the rest frame. The ?rst operator, Okin = 1 ? 1 ? hv (iD⊥ )2 hv → ? hv (iD )2 hv , 2mQ 2mQ (120)

is the gauge-covariant extension of the kinetic energy arising from the o?-shell residual motion of the heavy quark. The second operator is the non-Abelian analogue of the Pauli term, which describes the colour-magnetic interaction of the heavy quark spin with the gluon ?eld: gs ? gs ? hv σαβ Gαβ hv → ? hv S · Bc hv . (121) Omag = 4mQ mQ

1 ijk jk i Here S is the spin operator de?ned in (115), and Bc = ?2 ? G are the components of the colour-magnetic ?eld. This chromo-magnetic interaction is a relativistic e?ect, which scales like 1/mQ . This is the origin of the heavy quark spin symmetry.

3.3. Spectroscopic Implications The spin-?avour symmetry leads to many interesting relations between the properties of hadrons containing a heavy quark. The most direct consequences concern the spectroscopy of such states185 . In the mQ → ∞ limit, the spin of the heavy quark and the total angular momentum j of the light degrees of freedom are separately conserved by the strong interactions. Because of heavy quark symmetry, the dynamics is independent of the spin and mass of the heavy quark. Hadronic states can thus be classi?ed by the quantum numbers (?avour, spin, parity, etc.) of the light degrees of freedom. The spin symmetry predicts that, for ?xed j = 0, there is a doublet of degenerate states with total spin J = j ± 1 . The ?avour symmetry 2 relates the properties of states with di?erent heavy quark ?avour. Consider, as an example, the ground-state mesons containing a heavy quark. In this case the light degrees of freedom have the quantum numbers of an antiquark, and the degenerate states are the pseudoscalar (J = 0) and vector (J = 1) mesons. In the charm and bottom systems, one knows experimentally mB? ? mB ? 46 MeV , mD? ? mD ? 142 MeV . 41 (122)

These mass splittings are in fact reasonably small. To be more speci?c, at order 1/mQ one expects hyper?ne corrections to resolve the degeneracy, for instance mB? ? 2 2 2 mB ∝ 1/mb . This leads to the re?ned prediction m2 B ? ? mB ? mD ? ? mD ? const. The data are compatible with this:

2 2 m2 B ? ? mB ? 0.49 GeV , 2 2 m2 D ? ? mD ? 0.55 GeV .

(123)

One can also study excited meson states, in which the light constituents carry orbital angular momentum. It is tempting to interpret D1 (2420) with J P = 1+ and D2 (2460) with J P = 2+ as the spin doublet corresponding to j = 3 . The small mass 2 ? ? mD ? 35 MeV supports this assertion. One then expects di?erence mD2 1

2 2 2 2 m2 ? ? mB ? mD ? ? mD ? 0.17 GeV B2 1 1 2

(124)

for the corresponding states in the bottom system. The fact that this mass splitting is smaller than for the ground-state mesons is not unexpected. For instance, in the nonrelativistic constituent quark model the light antiquark in these excited mesons is in a p-wave state, and its wave function at the location of the heavy quark vanishes. Hence, in this model hyper?ne corrections are strongly suppressed. A typical prediction of the ?avour symmetry is that the “excitation energies” for states with di?erent quantum numbers of the light degrees of freedom are approximately the same in the charm and bottom systems. For instance, one expects mBS ? mB ? mDS ? mD ? 100 MeV , mB1 ? mB ? mD1 ? mD ? 557 MeV , ? ? mB ? mD ? ? mD ? 593 MeV . mB2 2

(125)

The ?rst relation has been very nicely con?rmed by the discovery of the BS meson at LEP186 . The observed mass124 mBS = 5.368 ± 0.005 GeV corresponds to an excitation energy mBS ? mB = 89 ± 5 MeV. 3.4. Weak Decay Form Factors Of particular interest are the relations between the weak decay form factors of heavy mesons, which parametrize hadronic matrix elements of currents between two meson states containing a heavy quark. These relations have been derived by Isgur and Wise13 , generalizing ideas developed by Nussinov and Wetzel10 and by Voloshin and Shifman11,12 . For the purpose of this discussion, it is convenient to work with a mass-independent normalization of meson states, M (p ′ )| M (p ) = 2 p0 (2π )3 δ 3 (p ? p ′ ) , mM (126)

instead of the more conventional, relativistic normalization M (p′ )|M (p) = 2p0 (2π )3 δ 3 (p ? p ′ ) . 42 (127)

In the ?rst case, p0 /mM = v 0 depends only on the meson velocity. In fact, it is more natural for heavy quark systems to use velocity rather than momentum variables. I will thus write |M (v ) instead of |M (p) . The relation to the conventionally ?1/2 normalized states is |M (v ) = mM M (p) .

Figure 13: The action of an external heavy quark current, as seen by the light degrees of freedom in the initial state. Consider now the elastic scattering of a pseudoscalar meson, P (v ) → P (v ′), induced by an external vector current coupled to the heavy quark contained in P . Before the action of the current, the light degrees of freedom in the initial state orbit around the heavy quark, which acts as a source of colour moving with the meson’s velocity v . On average, this is also the velocity of the “brown muck”. The action of the current is to replace instantaneously (at t = t0 ) the colour source by one moving at velocity v ′ , as indicated in Fig. 13. If v = v ′ , nothing really happens; the light degrees of freedom do not realize that there was a current acting on the heavy quark. If the velocities are di?erent, however, the “brown muck” suddenly ?nds itself interacting with a moving (relative to its rest frame) colour source. Soft gluons have to be exchanged in order to rearrange the light degrees of freedom and build the ?nal state meson moving at velocity v ′ . This rearrangement leads to a form factor suppression, which re?ects the fact that as the velocities become more and more di?erent, the probability for an elastic transition decreases. The important observation is that, in the mQ → ∞ limit, the form factor can only depend on the Lorentz boost γ = v · v ′ that connects the rest frames of the initialand ?nal-state mesons. That the form factor is independent of the heavy quark mass can also be seen from the following intuitive argument: The light constituents in the initial and ?nal state carry momenta that are typically of order ΛQCD v and ′ ΛQCD v ′ , respectively. Thus, their momentum transfer is q 2 ? Λ2 QCD (v · v ? 1), which is in fact independent of mQ . The result of this discussion is that in the e?ective theory, which provides the appropriate framework to consider the limit mQ → ∞ with the quark velocities kept ?xed, the hadronic matrix element describing the scattering process can be written as13 ? v′ γ ? hv |P (v ) = ξ (v · v ′ ) (v + v ′ )? , (128) P (v ′ )| h with a form factor ξ (v · v ′ ) that does not depend on mQ . Since the matrix element 43

is invariant under complex conjugation combined with an interchange of v and v ′ , the function ξ (v · v ′ ) must be real. That there is no term proportional to (v ? v ′ )? can be seen by contracting the matrix element with (v ? v ′ )? , and using v /hv = hv ′ ? ? v = hv′ . and hv′ / One can now use the ?avour symmetry to replace the heavy quark Q in one of the meson states by a heavy quark Q′ of a di?erent ?avour, thereby turning P into another pseudoscalar meson P ′ . At the same time, the current becomes a ?avour-changing vector current. In the in?nite mass limit, this is a symmetry transformation, under which the e?ective Lagrangian is invariant. Hence, the matrix element (129) P ′ (v ′ )| ? h′v′ γ ? hv |P (v ) = ξ (v · v ′ ) (v + v ′ )? is still determined by the same function ξ (v · v ′ ). This universal form factor is called the Isgur–Wise function, after the discoverers of this relation13 . ? ′ v ? hv is conserved in For equal velocities, the vector current J ? = ? h′v γ ? hv = h v the e?ective theory, irrespective of the ?avour of the heavy quarks, since ? ′ (iv · D hv ) + (iv · D h ? ′ ) hv = 0 i?? J ? = h v v

(130)

by the equation of motion, iv·D hv = 0, which follows from the e?ective Lagrangian (114). The conserved charges NQ′ Q = d3 x J 0 (x) = d3 x h′? v (x) hv (x) (131)

are generators of the ?avour symmetry. The diagonal generators simply count the number of heavy quarks, whereas the o?-diagonal ones replace a heavy quark by another: NQQ |P (v ) = |P (v ) and NQ′ Q |P (v ) = |P ′(v ) . It follows that P ′ (v )| NQ′Q |P (v ) = P (v )|P (v ) = 2v 0 (2π )3 δ 3 (0 ) , (132) and from a comparison with (129) one concludes that the Isgur–Wise function is normalized at the point of equal velocities: ξ (1) = 1 . (133)

This can easily be understood in terms of the physical picture discussed above: When there is no velocity change, the light degrees of freedom see the same colour ?eld and are in an identical con?guration before and after the action of the current. There is no form factor suppression. Since Erecoil = mP ′ (v · v ′ ? 1) is the recoil energy of the daughter meson P ′ in the rest frame of the parent meson P , the point v · v ′ = 1 is referred to as the zero recoil limit. The heavy quark spin symmetry leads to additional relations among weak decay form factors. It can be used to relate matrix elements involving vector mesons to those involving pseudoscalar mesons. In the e?ective theory, a vector meson with longitudinal polarization ?3 is related to a pseudoscalar meson by |V (v, ?3 ) = 2 S3 Q | P (v ) , 44 (134)

where S3 Q is a Hermitian operator that acts on the spin of the heavy quark Q. It follows that ?′ V ′ (v ′ , ?3 )| ? h′v′ Γ hv |P (v ) = P ′(v ′ )| 2 [S3 Q′ , hv′ Γ hv ] |P (v ) = P ′ (v ′ )| ? h′ ′ (2S 3 Γ) hv |P (v ) ,

v

(135)

where Γ can be an arbitrary combination of Dirac matrices, and S 3 is a matrix representation of the operator S3 Q′ . It is easiest to evaluate this expression in the rest frame of the ?nal-state meson, where [cf. (115)] v ′? = (1, 0, 0, 0) , ?? 3 = (0, 0, 0, 1) , S3 = 1 γ5 γ 0 γ 3 . 2 (136)

It is then straightforward to obtain the following commutation relations for the ? ′ ′ γ ? (1 ? γ5 ) hv : components of the weak current (V ? A)? = h v

3 3 0 0 2 [S3 Q′ , V ? A ] = A ? V , 0 0 3 3 2 [S3 Q′ , V ? A ] = A ? V ,

2 2 1 1 2 [S3 Q′ , V ? A ] = ?i(A ? V ) .

1 1 2 [S3 Q′ , V ? A ] =

i(A2 ? V 2 ) ,

(137)

Using (135) and (137), one can relate the matrix element of the weak current between a pseudoscalar and a vector meson to the matrix element of the vector current between two pseudoscalar mesons given in (129). The result can be written in the covariant form13 :

′ ′ V ′ (v ′ , ?)| ? h′v′ γ ? (1 ? γ5 ) hv |P (v ) = i??ναβ ?? ν vα vβ ξ (v · v )

? ??? (v · v ′ + 1) ? v ′? ?? · v ξ (v · v ′ ) , (138)

where ?0123 = ??0123 = ?1. Once again, the matrix element is completely described in terms of the universal Isgur–Wise form factor. Equations (129) and (138) summarize the relations imposed by heavy quark symmetry on the weak decay form factors describing the semileptonic decay processes B → D?ν ? and B → D ? ? ν ?. These relations are model-independent consequences of QCD in the limit where mb , mc ? ΛQCD . They play a crucial role in the determination of | Vcb|, as I will discuss later in this lecture. 3.5. The 1/mQ Expansion At tree level, eq. (113) de?nes the Lagrangian of HQET as a series of local, higher-dimension operators multiplied by powers of 1/mQ . The expression for Hv in (108) can be used to derive a similar expansion for the full heavy quark ?eld

45

Q(x): Q(x) = e?imQ v·x 1 + = e?imQ v·x 1 iD / ⊥ hv (x) (2mQ + iv · D ? iη ) iD /⊥ 1+ + . . . hv (x) . 2mQ

(139)

This relation contains the recipe for the construction (at tree level) of any operator in HQET that contains one or more heavy quark ?elds. For instance, the vector current V ? = q ? γ ? Q composed of a heavy and a light quark is represented as V ? (x) = e?imQ v·x q ?(x) γ ? 1 + iD /⊥ + . . . hv (x) . 2mQ (140)

Matrix elements of this current can be parametrized by hadronic form factors, and the purpose of using an e?ective theory is to make the mQ dependence of these form factors explicit. Consider, as an example, the matrix element of V ? (0) between a heavy meson M (v ) and some light ?nal state ?: ? | V ? | M (v ) = ? | q ? γ ? hv |M (v ) + 1 ?|q ? γ? i D / ⊥ hv |M (v ) + . . . . 2mQ (141)

It would be nice if the matrix elements on the right-hand side of this equation were independent of mQ . Then the second term would give the leading power correction to the ?rst one. However, the equation of motion for hv derived from (113) contains 1/mQ corrections, too. This leads to an mQ dependence of any hadronic matrix element of operators containing such ?elds. Another way to say this is that the eigenstates of the e?ective Lagrangian Le? (supplemented by the standard QCD Lagrangian for the light quarks and gluons) depend, at higher order in 1/mQ , on the heavy quark mass. This is no surprise, since the e?ective Lagrangian is equivalent to the Lagrangian of the full theory. It is better to organize things in a slightly di?erent way by working with the eigenstates of only the leading-order e?ective Lagrangian L∞ in (114), treating the higher-dimension operators perturbatively as power corrections24,184,187 . Then the equation of motion satis?ed by hv is simply iv · D hv = 0 , (142)

and the states of the e?ective theory are truly independent of mQ . However, these states are now di?erent from the states of the full theory. For instance, instead of (141) one should write ? | V ? | M (v )

QCD

= ?|q ? γ ? hv |M (v ) ∞ 1 + ?|q ? γ? i D / ⊥ hv |M (v ) 2mQ 46

∞

1 ? | i dy T { q ? γ ? hv (0), ? hv (iD⊥ )2 hv (y ) } |M (v ) ∞ 2mQ gs ? | i dy T { q ? γ ? hv (0), ? hv σαβ Gαβ hv (y ) } |M (v ) ∞ + 4mQ + O (1/m2 (143) Q) . + The matrix elements on the right-hand side are independent of the heavy quark mass. The mass dependence of the states |M (v ) QCD is re?ected in HQET by the appearance of the third and fourth term, which arise from an insertion of the ?rstorder power corrections in the Lagrangian into the leading-order matrix element of the current. This insertion can be thought of as being a correction to the wave function of the heavy meson. From now on I will omit the subscripts on the states. I will instead use the symbol “? =” and write V ? (0) ? ? γ ? hv (0) + =q 1 q ? γ? i D / ⊥ hv (0) 2mQ 1 ? v (iD⊥ )2 hv (y ) } + i dy T { q ? γ ? hv (0), h 2mQ gs ? v σαβ Gαβ hv (y ) } + . . . (144) i dy T { q ? γ ? hv (0), h + 4mQ

to indicate that the operators on the two sides of this equation have to be evaluated between di?erent states. 3.6. Matching and Running In Section 3.2, I have discussed the ?rst two steps in the construction of HQET. Integrating out the small components in the heavy quark spinor ?elds, a nonlocal e?ective action was derived. This was then expanded in a series of local operators of increasing dimension to obtain an e?ective Lagrangian. A similar expansion could be written down for any external current. The e?ective Lagrangian and the e?ective currents derived that way correctly reproduce the long-distance physics of the full theory. They cannot describe the short-distance physics correctly, however. The reason is obvious: The heavy quark participates in strong interactions through its coupling to gluons. These gluons can be soft or hard, i.e. their virtual momenta can be small, of the order of the con?nement scale, or large, of the order of the heavy quark mass. But hard gluons can resolve the nonlocality of the propagator of the small component ?elds Hv . Their e?ects are not taken into account in the na¨ ?ve version of the OPE, which was used in the derivation of the e?ective Lagrangian in (113) and the e?ective vector current in (144). So far, the e?ective theory provides an appropriate description only at scales ? ? mQ . In this section, I will discuss the systematic treatment of short-distance corrections. A new feature of such corrections is that through the running coupling constant they induce a logarithmic dependence on the heavy quark mass11 . The 47

important observation is that αs (mQ ) is small, so that these e?ects can be calculated in perturbation theory. Consider, as an example, matrix elements of the vector current V ? = q ? γ ? Q. In QCD this current is partially conserved and needs no renormalization188 . Its matrix elements are free of ultraviolet divergences. Still, these matrix elements can have logarithmic dependence on mQ from the exchange of hard gluons with virtual momenta of the order of the heavy quark mass. If one goes over to the e?ective theory by taking the limit mQ → ∞, these logarithms diverge. Consequently, the vector current in the e?ective theory does require a renormalization18 . Its matrix elements depend on an arbitrary renormalization scale ?, which separates the regions of short- and long-distance physics. If ? is chosen such that ΛQCD ? ? ? mQ , the e?ective coupling constant in the region between ? and mQ is small, and perturbation theory can be used to compute the short-distance corrections. These corrections have to be added to the matrix elements of the e?ective theory, which only contain the long-distance physics below the scale ?. Schematically, then, the relation between matrix elements in the full and in the e?ective theory is V ? (mQ )

QCD

= C0 (mQ , ?) V0(?)

∞

+

C1 (mQ , ?) V 1 (? ) 2mQ

∞

+ ... ,

(145)

where I have indicated that matrix elements of V ? in the full theory depend on mQ , whereas matrix elements of operators in the e?ective theory are mass-independent, but do depend on the renormalization scale. The Wilson coe?cients Ci (mQ , ?) are de?ned by this relation. Order by order in perturbation theory, they can be computed from a comparison of the matrix elements in both theories. Since the e?ective theory is constructed to reproduce correctly the low-energy behaviour of the full theory, this “matching” procedure is independent of any long-distance physics, such as infrared singularities, nonperturbative e?ects, the nature of the external states used in the matrix elements, physical cuts, etc. Only at high energies do the two theories di?er, and these di?erences are corrected for by the short-distance coe?cients. The calculation of the coe?cient functions in perturbation theory uses the powerful methods of the renormalization group, which I have discussed in my ?rst lecture. It is in principle straightforward, yet in practice rather tedious. Much of the recent work on heavy quark symmetry has been devoted to this subject. In this section, I will discuss as an illustration the wave-function renormalization of the heavy quark ?eld in HQET, and the short-distance expansion of a heavy–light vector current to leading order in 1/mQ . For a more comprehensive presentation of short-distance corrections, I refer to my review article25 . In quantum ?eld theory, the parameters and ?elds of the Lagrangian have no direct physical signi?cance. They have to be renormalized before they can be related to observable quantities. In an intermediate step the theory has to be regularized. The most convenient regularization scheme in QCD is dimensional regularization28,32,33 , in which the dimension of space-time is analytically continued 48

to D = 4 ? 2?, with ? being in?nitesimal. Loop integrals that are logarithmically divergent in four dimensions become ?nite for ? > 0. From the fact that the action S = dD x L(x) is dimensionless, one can derive the mass dimensions of the ?elds and parameters of the theory. For instance, one ?nds that the “bare” coupling conbare bare stant αs is no longer dimensionless if D = 4: dim[ αs ] = 2?. In a renormalizable theory, it is possible to rewrite the Lagrangian in terms of renormalized quantities, in such a way that Green’s functions of the renormalized ?elds remain ?nite as 1/2 ? → 0. For QCD, one introduces renormalized quantities by Qbare = ZQ Qren , 1/2 bare ren Abare = ZA Aren , αs = ?2? Zα αs , etc., where ? is an arbitrary mass scale introduced to render the renormalized coupling constant dimensionless. Similarly, in 1/2 HQET one de?nes a renormalized heavy quark ?eld by hbare = Zh hren v v . From now on, the superscript “ren” will be omitted. As a warm-up, let me sketch the calculation of the wave-function renormalization for a heavy quark in HQET. In the MS scheme, Zh can be computed from the 1/? ?pole in the quark self-energy shown in Fig. 14: ? Σ(v · k ) 1 ?1 . 1 ? Zh = -pole of ? ? ?v · k (146)

As long as v · k < 0, the self-energy is infrared ?nite and real. The result is gaugedependent, however. In Feynman gauge, one obtains at one-loop order Σ(v · k ) = ?

2 4igs 3 ∞

8ig 2 =? s 3

dD t 1 D 2 (2π ) (t + iη )[v · (t + k ) + iη ]

dρ

0 ∞

1 dD t D (2π ) [t2 + 2ρ v · (t + k ) + iη ]2 ρ2 + ρω 4π?2

??

2αs = Γ(?) 3π

dρ

0

,

(147)

where ρ is a dimensionful Feynman parameter, and ω = ?2v · k > 0 acts as an infrared cuto?. A straightforward calculation leads to 4αs ω2 ? Σ(v · k ) = Γ(1 + ?) ?v · k 3π 4π?2 =

?? 0 2 1

dz z ?1+2? (1 ? z )??

??

4αs ω Γ(2?) Γ(1 ? ?) 3π 4π?2

,

(148)

where I have substituted ρ = ω (1 ? z )/z . From an expansion around ? = 0, one obtains 2αs Zh = 1 + . (149) 3π ? ? This result was ?rst derived by Politzer and Wise18 . In the meantime, the calculation was also done at the two-loop order189 ? 191 . 49

Figure 14: Self-energy ?iΣ(v · k ) of a heavy quark in HQET. The velocity v is conserved by the strong interactions. Similar to the ?elds and coupling constants, any composite operator built from quark and gluon ?elds may require a renormalization beyond that of its component ?elds. Such operators can be divided into three classes: gauge-invariant operators that do not vanish by the equations of motion (class-I), gauge-invariant operators that vanish by the equations of motion (class-II), and operators which are not gaugeinvariant (class-III). In general, operators with the same dimension and quantum numbers mix under renormalization. However, things simplify if one works with the background ?eld technique192 ? 195 , which serves for quantizing gauge theories, preserving explicit gauge invariance. This o?ers the advantage that a class-I operator cannot mix with class-III operators, so that only gauge-invariant operators need be considered196 . Furthermore, class-II operators are irrelevant since their matrix elements vanish by the equations of motion. It it thus su?cient to consider class-I operators only. For a set {Oi } of n class-I operators that mix under renormalization, one de?nes an n × n matrix of renormalization factors Zij by (summation over bare j is understood) Oi = Zij Oj , such that the matrix elements of the renormalized operators Oj remain ?nite as ? → 0. In contrast to the bare operators, the renormalized ones depend on the subtraction scale via the ? dependence of Zij : ? where

?1 γik = Zij ?

d d ?1 bare Oj = ?γik Ok , Oi = ? Zij d? d? d Zjk d?

(150)

(151)

are called the anomalous dimensions. It is sometimes convenient to introduce a com? is the pact matrix notation, in which O is the vector of renormalized operators, Z matrix of renormalization factors, and γ ? denotes the anomalous dimension matrix. Then the running of the renormalized operators is controlled by the renormalizationgroup equation (RGE) d ? +γ ? O = 0. (152) d? As an example, consider the short-distance expansion of the heavy–light vector current V ? = q ? γ ? Q to leading order in 1/mQ . In the full theory, the vector current is not renormalized. But the e?ective current operators in HQET do require

50

renormalization. The general form of the short-distance expansion reads V ? (mQ ) ? = Ci (mQ , ?) Oi(?) + O (1/mQ )

?1 bare = Ci (mQ , ?) Zij (?) Oj + O (1/mQ) ,

(153)

where I have indicated the ? dependence of the renormalized operators. This is the correct generalization of (145) in the case of operator mixing. In general, a complete set of operators with the same quantum numbers appears on the right-hand side. In the case at hand, there are two such operators, namely O1 = q ? γ ? hv , O2 = q ? v ? hv . (154)

From (152) and the fact that the product Ci (mQ , ?) Oi (?) must be ?-independent, one can derive the RGE satis?ed by the coe?cient functions. It reads ? d ?γ ? t C (mQ , ?) = 0 , d? (155)

where γ ? t denotes the transposed anomalous dimension matrix, and I have collected the coe?cients into a vector. The solution of the RGE proceeds as I described in my ?rst lecture; however, it is technically more complicated in the case when there is operator mixing. One obtains ? (?, mQ ) C (mQ , mQ ) , C (mQ , ?) = U with the evolution matrix27,43,44 αs (?) ? (?, mQ ) = Tα exp U αs (mQ ) dα γ ? t (α ) . β (α ) (157) (156)

Here “Tα ” means an ordering in the coupling constant such that the couplings increase from right to left (for ? < mQ ). This is necessary since, in general, the anomalous dimension matrices at di?erent values of α do not commute: [? γ (α1 ), γ ? (α2 )] = 0. In the present case, however, these complications are absent, since the anomalous dimension matrix turns out to be proportional to the unit matrix. I will now illustrate the solution of (156) to next-to-leading order in renormalization-group improved perturbation theory. There are two ingredients to this calculation. First, the Wilson coe?cients at the high-energy scale ? = mQ must be calculated to one-loop order, C (mQ , mQ ) = C(0) + C(1) αs (mQ ) + ... , 4π (158)

by matching QCD onto the e?ective theory. Then, in order to control the running of the HQET operators, the anomalous dimension matrix γ ? must be calculated to twoloop accuracy. Let me start with a comparison of the one-loop matrix elements of the 51

vector current in the full and in the e?ective theory. It is legitimate to perform the matching calculation with on-shell quark states. Then the matrix elements can be written as V ? = u ?q Γ? uQ , where the heavy quark spinor satis?es v / uQ = uQ . There will be infrared divergences in the calculation, which I regulate by introducing a ?ctitious gluon mass λ. At tree level, the vertex function in both theories is simply Γ? = γ ? . The one-particle irreducible one-loop diagrams are shown in Fig. 15. They have to be supplemented by an on-shell wave-function renormalization of the external lines. The complete one-loop vertex function in QCD is25 m2 αs 2αs ? 11 Q ln 2 ? γ? + v . (159) 2π λ 6 3π As required by the nonrenormalization theorem for partially conserved currents, this result is gauge-independent and ultraviolet ?nite188 . In the limit of degenerate quark masses, the vector is conserved. The space integral over its time component is the generator of a ?avour symmetry. It cannot be renormalized. When the symmetry is softly broken by the presence of mass splittings, this is only relevant for small loop momenta but does not a?ect the ultraviolet region. Thus, there can only be a ?nite renormalization. Notice, however, that (159) does contain an infrared singularity, because the matrix element was calculated using unphysical states. Γ? QCD = 1 +

Figure 15: One-loop corrections to the matrix elements of the vector current in QCD and in HQET. The wavy line represents the current. The external momenta are on-shell. In the e?ective theory, the one-loop vertex function for any heavy–light current q ? Γ hv is found to be25 αs 1 ?2 5 1+ Γ, (160) + ln 2 + 2π ? ? λ 6 where the Dirac structure of Γ is completely arbitrary. It is generally true that these operators are renormalized multiplicatively and irrespectively of their Dirac structure11 . Since there is no approximate ?avour symmetry relating light and heavy quarks in the e?ective theory, it is not unexpected that the matrix elements of the bare currents are ultraviolet divergent. In the MS scheme, the currents are ? with components renormalized by a diagonal matrix Z αs , Z12 = Z21 = 0 . (161) Z11 = Z22 = 1 + 2π ? ? 52

From (153), it then follows that the renormalized one-loop vertex function is Γ? HQET αs = 1+ 2π ?2 5 ln 2 + λ 6 C1 (mQ , ?) γ ? + C2 (mQ , ?) v ? . (162)

Notice that the λ dependence is the same as in (159). The short-distance coe?cients, which follow from a comparison of the two results, are independent of the infrared regulator: C1 (mQ , ?) = 1 + C2 (mQ , ?) = αs π ln mQ 4 , ? ? 3 (163)

2αs . 3π

In particular, one obtains for the matching conditions at ? = mQ : C1 (mQ , mQ ) = 1 ? 4αs (mQ ) , 3π 2αs (mQ ) . C2 (mQ , mQ ) = 3π

(164)

This completes the ?rst step of the calculation. The next step is to calculate the one- and two-loop coe?cients of the anomalous dimension, as de?ned in (13). From (160), it follows that the anomalous dimension matrix is proportional to the unit matrix, with the one-loop coe?cient γ0 = ?4 . (165)

This is the so-called hybrid anomalous dimension of heavy–light currents derived by Voloshin and Shifman11 . In the context of the e?ective theory, it has been calculated by Politzer and Wise18 . The two-loop coe?cient γ1 has been calculated by Ji and Musolf189 , and by Broadhurst and Grozin190 . They obtain γ1 = ? 254 56 2 20 ? π + nf . 9 27 9 (166)

The ?nal result for the Wilson coe?cients can be derived using (20). It is possible to write the answer in the factorized form Ci (mQ , ?) = Ci(mQ ) Khl (?), where the dependence on the two scales has been completely separated. One ?nds αs (mQ ) 4 S? π 3 ?2/β0 2αs (mQ ) C2 (mQ ) = αs (mQ ) , 3π 2/β0 αs (?) S , Khl (?) = αs (?) 1? π C1 (mQ ) = αs (mQ )

?2/β0

1+

,

(167)

53

where S=

381 + 28π 2 ? 30nf 153 ? 19nf γ1 β0 ? β1 γ0 ? . = 3 2 8β0 (33 ? 2nf )2 36 (33 ? 2nf )

(168)

3.7. Flavour-Changing Heavy Quark Currents The OPE of currents in the e?ective theory becomes considerably more complicated in the case of ?avour-changing currents. The weak current c ? γ ? (1 ? γ5 ) b is of this form, however, so it is necessary to consider this case in detail. The complications arise from the fact that there are two di?erent heavy quark masses, mb and mc . Thus the calculation of the Wilson coe?cients becomes a two-scale problem. In addition, the coe?cients will depend in a nontrivial way on the velocity transfer w = v · v′. There are two obvious ways of performing the transition from QCD to an e?ective low-energy theory in which both the bottom and the charm quarks are treated as heavy quarks (in the sense of HQET): The transition can either be done in a single step, or by considering ?rst an intermediate theory with a static bottom quark, but a dynamical charm quark. The latter becomes heavy in a second step. If one could solve perturbation theory to all orders, both treatments would lead to the same results for the Wilson coe?cients. The calculation di?ers in both cases, however, and the results also di?er if the perturbation series is truncated. To see what the di?erences are, suppose ?rst the two heavy quarks have similar masses, i.e. mb ? mc ? m, with m being some average mass. It is then natural to remove at the same time the dynamical degrees of freedom of both heavy quarks. Let me consider the cases of the vector current V ? = c ? γ ? b and of the axial vector ? ? current A = c ? γ γ5 b explicitly and denote them collectively by J ? . Each of these currents obeys an expansion of the form J (mb , mc ) ? =

? 3

Cj (mb , mc , m, w, ?) Jj (?) + O (1/m) ,

j =1

(169)

? c ′ Γj hb are local current operators in the e?ective theory with matrices where Jj = h v v Γ1 = γ ? , for the vector current, and Γ1 = γ ? γ 5 , Γ2 = v ? γ 5 , Γ3 = v ′? γ5 (171) Γ2 = v ? , Γ3 = v ′? (170)

for the axial vector current. In the above expression, I have indicated that matrix elements of the original current J ? between states of the full theory depend on mb and mc , whereas matrix elements of the operators Jj between states of the e?ective theory are mass-independent, but do depend on the renormalization scale ?. The short-distance coe?cients are functions of the heavy quark masses, the 54

renormalization scale, and the matching scale m. The dependence on m would disappear if one could sum the perturbative series to all orders. The advantage of this ?rst approach is its simplicity. At leading order in the 1/m expansion, only three current operators contribute. Matrix elements of higher-dimension operators are suppressed by powers of ΛQCD /m. The short-distance coe?cients contain the dependence on mb and mc correctly to a given order in αs , via matching at ? = m. The only disadvantage is the residual dependence on the matching scale m, which arises when one calculates to ?nite order in perturbation theory. Although in a next-to-leading order calculation the scale in the leading anomalous scaling factor is determined, one has no control over the scale in the next-to-leading corrections proportional to αs (m). Since mb > m > mc , this introduces an uncertainty of order 2 αs ln(mb /mc ). The alternative approach is to consider ?rst, as an intermediate step for mb > ? > mc , an e?ective theory with only a heavy bottom quark. Denoting the e?ective (k ) current operators of dimension (3 + k ) in this theory by Ji (mc , ?), indicating that their matrix elements depend both on the charm quark mass and the renormalization scale, I write the short-distance expansion as an expansion in 1/mb : J (mb , mc ) ? =

? ∞ k =0

1 mk b

Di (mb , ?) Ji (mc , ?) .

i

(k )

(k )

(172)

Since the velocity of the charm quark is still a dynamical degree of freedom, in this intermediate e?ective theory there are only two dimension-three operators, namely J1 = c ? Γ1 hb v,

(0)

J2 = c ? Γ2 hb v.

(0)

(173)

A major complication arises from the fact that matrix elements of the higher(k ) dimension operators Ji will, in general, scale like mk c . This compensates the k prefactor 1/mb . Consequently, these operators cannot be neglected even at leading order in the heavy quark expansion. One would thus have to deal with an in?nite number of operators in order to keep track of the full dependence on the heavy quark masses. Ignoring this di?culty for the moment, one may use (172) to scale the currents from ? = mb down to ? = mc , where the matching is done onto the ?nal e?ective theory with two heavy quarks. In this step, the operator basis collapses considerably. When terms of order ΛQCD /mQ (here I use mQ generically for mc or mb ) are neglected on the level of matrix elements, only the three operators Jj (?) in (169) remain. Each of the operators of the intermediate theory has an expansion in terms of these operators, with coe?cients Eij :

(k ) Ji (mc , mc ) ? = 3 j =1

mk c Eij (mc , w, ?) Jj (?) + O (1/mc ) .

(k )

(174)

Combining this with (172), one obtains J (mb , mc ) ? =

? 3

Cj (mb , mc , w, ?) Jj (?) + O (1/mQ) ,

j =1

(175)

55

with evolution coe?cients197

∞

Cj (mb , mc , w, ?) =

k =0

mc mb

k i

Di (mb , mc ) Eij (mc , w, ?) .

(k )

(k )

(176)

In this expression, the matching scale m of the ?rst approach does not appear. The coe?cients depend either on αs (mb ) or αs (mc ), i.e. the scaling in the intermediate region mb > ? > mc is properly taken into account. To achieve this, however, it would be necessary to consider an in?nite number of operators in the intermediate e?ective theory. This is, of course, not manageable. It is important to realize, however, that the short-distance coe?cients in (169) and (175) must agree, and that this equality must hold order by order in an expansion in the mass ratio mc /mb . Using this fact, one can combine the two approaches into a consistent next-to-leading order calculation197 . For details of the very tedious calculation, the reader is referred to the literature22,197 ? 201 . Below, I will only give numerical results. The short-distance coe?cients can be written in the factorized form Ci (mb , mc , w, ?) = Ci (mb , mc , w ) Khh(w, ?) .

(5) (5)

(177)

I denote the coe?cients for the vector current by Ci , and those for the axial vector current by Ci5 . The ?-dependent function Khh (w, ?) is universal and normalized to unity at zero recoil: Khh (1, ?) = 1. The renormalization-group invariant coe?cients (5) Ci contain all dependence on the heavy quark masses. In Table 1, I give the numerical values of these coe?cients, which are obtained using mb = 4.80 GeV and mc = 1.45 GeV for the heavy quark masses, as well as ΛMS = 0.25 GeV (for nf = 4) in the two-loop expression for αs (?). The corresponding coupling constants are αs (mb ) ? 0.20 and αs (mc ) ? 0.32. Table 1: Short-distance coe?cients for b → c transitions. w 1.0 1.1 1.2 1.3 1.4 1.5 1.6 1.7 1.8 C1 1.136 1.107 1.081 1.056 1.033 1.011 0.991 0.972 0.953 C2 ?0.085 ?0.080 ?0.077 ?0.073 ?0.070 ?0.067 ?0.064 ?0.062 ?0.059 C3 ?0.021 ?0.021 ?0.020 ?0.019 ?0.018 ?0.018 ?0.017 ?0.017 ?0.016

5 C1 0.985 0.965 0.946 0.927 0.910 0.894 0.878 0.864 0.850 5 C2 ?0.122 ?0.115 ?0.109 ?0.103 ?0.098 ?0.094 ?0.089 ?0.086 ?0.082 5 C3 0.042 0.040 0.038 0.036 0.035 0.033 0.032 0.031 0.030

Of particular interest for the model-independent determination of | Vcb| to be 5 discussed below is the value of the short-distance coe?cient C1 (w ) at zero recoil. 56

One ?nds197

5 ηA ≡ C1 (1) = x6/25 1 + 1.561

αs (mc ) ? αs (mb ) αs (mc ) ? π π 1 ?12/25 25 14 ?9/25 +z ? x + x + 54 27 18

2z 2 8 + ln z 3 1?z 8 ln x , (178) 25

where x = αs (mc )/αs (mb ), and z = mc /mb . For ΛMS = 0.25 ± 0.05 GeV and z = 0.30 ± 0.05, this yields ηA = 0.986 ± 0.006. Let me stop, at this point, the discussion of short-distance e?ects. There have been important calculations that I have no time to mention here, such as the renormalization of the higher-dimension operators in the e?ective Lagrangian19,23,202 , and the short-distance expansion of currents to next-to-leading order in203 1/mQ . A discussion of these calculations can be found in my review article25 . 3.8. Covariant Representation of States The purpose of the OPE discussed in the previous sections was to disentangle the short-distance physics related to length scales set by the Compton wavelengths of the heavy quarks from con?nement e?ects relevant at large distances. This procedure makes explicit the mQ dependence of any Green’s function of the full theory, which contains one or more heavy quark ?elds. Hadronic matrix elements of heavy quark currents have a 1/mQ expansion as shown in (145). The HQET matrix elements on the right-hand side of this equation contain the long-distance physics associated with the interactions of the cloud of light degrees of freedom among themselves and with the background colour-?eld provided by the heavy quarks. These hadronic quantities depend in a most complicated way on the “brown muck” quantum numbers of the external states, the quantum numbers of the current, and on the heavy quark velocities. They are related to matrix elements of e?ective current operators in HQET. Recall that these matrix elements are independent of the heavy quark masses, if the states of the e?ective theory are taken to be the eigenstates of the leading-order Lagrangian L∞ in (114). Matrix elements evaluated using these states have a well-de?ned behaviour under spin-?avour symmetry transformations. When combined with the requirements of Lorentz covariance, restrictive constraints on their structure can be derived. I will now introduce a very elegant formalism due to Bjorken204 and Falk et al.22,205 , which allows one to derive the general form of such matrix elements in a straightforward manner. The clue is to work with a covariant tensor representation of states with de?nite transformation properties under the Lorentz group and the heavy quark spin-?avour symmetry. The eigenstates of HQET can be thought of as the “would-be hadrons” built from an in?nitely heavy quark dressed with light quarks, antiquarks and gluons. In such a state, both the heavy quark and the cloud of light degrees of freedom have well-de?ned transformation properties under the Lorentz group. The heavy quark can be represented by a spinor uh (v, s) satisfying v / uh (v, s) = uh (v, s), where v is the velocity of the hadron. Because of heavy quark symmetry, the wave function 57

of the state (when properly normalized) is independent of the ?avour and spin of the heavy quark, and the states can be characterized by the quantum numbers of the “brown muck”. In particular, for each con?guration of light degrees of freedom with total angular momentum j ≥ 0 and parity P , there is a degenerate doublet of states with spin-parity J P = (j ± 1 )P . Following Falk205 , I discuss the cases of 2 integral and half integral j separately. Hadronic states with integral j have odd fermion number and correspond to baryons; states with half-integral j have even fermion number and correspond to mesons. First consider the heavy baryons. In this case, the “brown muck” is an object with spin-parity j P that can be represented by a totally symmetric, traceless tensor A?1 ...?j subject to the transversality condition v? A?1 ...?j = 0. States are said to have “natural” parity if P = (?1)j , and “unnatural” parity otherwise. The composite heavy baryon can be represented by the tensor wave function ψ ?1 ...?j = uh A?1 ...?j . (179)

Under a connected Lorentz transformation Λ, this object transforms as a spinortensor ?eld ?j ν1 ...νj 1 ψ ?1 ...?j → Λ? , (180) ν1 . . . Λνj D (Λ) ψ

i where D (Λ) = exp(? 4 ω?ν σ ?ν ) is the usual spinor representation of Λ. A heavy quark spin rotation Λ, on the other hand, acts only on uh ; hence

ψ ?1 ...?j → D (Λ) ψ ?1 ...?j .

(181)

Here Λ is restricted to spatial rotations (in the rest frame). The in?nitesimal form of D (Λ) was considered in (117). The simplest but important case j P = 0+ corresponds 1+ . It can be represented to the ground-state ΛQ -baryon with total spin-parity J P = 2 by a spinor uΛ . Since the light degrees of freedom are in a con?guration of total spin zero, the spin of the baryon is carried by the heavy quark, and the spinor uΛ coincides with the heavy quark spinor. Hence ψΛ = uΛ (v, s) = uh (v, s) . (182)

For j ≥ 0, the object ψ ?1 ...?j does not transform irreducibly under the Lorentz group, but is a linear combination of two components with total spin j ± 1 . These 2 correspond to a degenerate doublet of physical states, which only di?er in the orientation of the heavy quark spin relative to the angular momentum of the light degrees of freedom. The nonrelativistic quark model suggests that one should identify the + 1+ P states with j P = 1+ with the ΣQ (J P = 2 ) and Σ? = 3 ) baryons. In the Q (J 2 quark model, these states contain a heavy quark and a light vector diquark with no orbital angular momentum. Unlike the ΛQ -baryons, they have unnatural parity. (?) This implies that decays between ΛQ and ΣQ must be described by parity-odd form factors206,207 . 58

in Next consider the heavy mesons. I shall only discuss the case j P = 1 2 detail. Since quarks and antiquarks have opposite intrinsic parity, the corresponding physical states with “natural” parity are the ground-state pseudoscalar (J P = 0? ) and vector (J P = 1? ) mesons. As before, the heavy quark is represented by a spinor uh (v, s) subject to the condition v / uh (v, s) = uh (v, s). The light degrees of freedom as a whole transform under the Lorentz group as an antiquark moving at velocity v . They are described by an antifermion spinor v ?? (v, s′ ) satisfying v ?? (v, s′ ) / v = ′ ?v ?? (v, s ). The ground-state mesons can be represented by the composite object ψ = uh v ?? , which is a 4 × 4 Dirac matrix with two spinor indices, one for the heavy quark and one for the “brown muck”. Under a connected Lorentz transformation Λ, the meson wave function ψ transforms as ψ → D (Λ) ψ D?1 (Λ) , ? ψ. ψ → D (Λ) (183)

?

? whereas under a heavy quark spin rotation Λ

(184)

The composite ψ represents a linear combination of the physical pseudoscalar and vector meson states. It is easiest to identify these states in the rest frame, where uh has only upper components, whereas v ?? has only lower components. The nonvanishing components of ψ are thus contained in a 2 × 2 matrix, which can be written as a linear combination of the identity I and the Pauli matrices σ i . Let me choose the quantization axis in 3-direction and work with the rest-frame spinor basis 1 ?0? ? ? ? , u h (? ) = ? ?0? 0

? ?

Then a basis of states is:

0 ?1? ? ? ? , u h (? ) = ? ?0? 0 ( ?↓ + ?↑ ) = ? ( ?↓ ? ?↑ ) = ? √

? ?

0 ?0? ? ? ? , v? (↑) = ? ?0? 1 0 I 0 0 0 σ3 0 0 , ,

? ?

v? (↓) =

? ? ?0? ? ? ? ? ?1?

0

.

(185)

0

1 2 ( ?↑ ) = ? √ 2 √ 1 2 ( ?↓ ) = ? √ 2

0 σ 1 + iσ 2 0 0 0 σ 1 ? iσ 2 0 0

, . (186)

Let me furthermore de?ne two transverse polarization vectors ?± and a longitudinal polarization vector ?3 by 1 ?? ± = √ (0, 1, ±i, 0) , 2 59 ?? 3 = (0, 0, 0, 1) . (187)

It is then obvious to identify the pseudoscalar (P ) and vector (V ) meson states as22,204,205 P (v = 0) = ? V (v = 0, ?) = 1 + γ0 γ5 , 2 1 + γ0 ?. / 2

(188)

The second state in (186) has longitudinal polarization, whereas the last two states have transverse polarization. To get familiar with this representation, consider the action of the spin operator Σ on ψ . A matrix representation of the components Σi in the rest frame is Σi = 1 γ γ 0 γ i , and the action of the operator Σi on the meson wave function is Σi ψ = 2 5 [Σi , ψ ]. Using this, one ?nds Σ2 P Σ2 V (?) Σ3 V (?± ) Σ3 V (?3 ) = = = = Σ3 P = 0 , 2 V (?) , ±V (?± ) , 0,

(189)

which shows that P has total spin zero, and V has total spin one. Next consider the action of the heavy quark spin operator S. It has the same matrix representation as Σ, but only acts on the heavy quark spinor in ψ : Si ψ = S i ψ , with S i = Σi . It follows that 1 S3 P = V (?3 ) , 2 1 3 S V (?3 ) = P , 2 1 3 S V (?± ) = ± V (?± ) , (190) 2 in accordance with the spin assignments for the heavy quark in (186). In a general frame, the tensor wave functions in (188) can be readily generalized in a Lorentz-covariant way by replacing γ 0 with v /. The covariant representation of states can be used to determine in a very e?cient way the structure of hadronic matrix elements in the e?ective theory. The goal is to ?nd a minimal form factor decomposition consistent with Lorentz covariance, parity invariance of the strong interactions, and heavy quark symmetry. The ?avour symmetry is manifest when one uses mass-independent wave functions. The correct transformation properties under the spin symmetry are guaranteed when one collects a spin-doublet of states into a single object. In case of the ground-state pseudoscalar and vector mesons, for instance, one introduces a covariant tensor wave function M(v ) that represents both P (v ) and V (v, ?) by M (v ) = 1+/ v 2 ?γ5 ; pseudoscalar meson, /; ? vector meson. 60 (191)

(1 ± / v ). This is It has the important property M(v ) = P+ M(v ) P? , where P± = 1 2 often used to simplify expressions. Consider now a transition between two heavy mesons, M (v ) → M ′ (v ′ ), mediated ? ′ ′ Γ hv , which changes a heavy quark Q by a renormalized e?ective current operator h v ′ into another heavy quark Q . According to the Feynman rules of HQET, the “heavy quark part” of the decay amplitude is simply proportional to u ?′h Γ uh ; interactions of the heavy quarks with gluons do not modify the Dirac structure of Γ. Since the heavy quark spinors are part of the tensor wave functions associated with the ′ hadron states, it follows that the amplitude must be proportional to M (v ′ ) Γ M(v ). This is a Dirac matrix with two indices representing the light degrees of freedom. Since the total matrix element is a scalar, these indices must be contracted with those of a matrix Ξ. Hence one may write M ′ (v ′ )| ? h′v′ Γ hv |M (v ) = Tr Ξ(v, v ′ , ?) M (v ′ ) Γ M(v ) .

′

(192)

The matrix Ξ contains all long-distance dynamics. It is a most complicated object, only constrained by the symmetries of the e?ective theory. Heavy quark symmetry requires that it be independent of the spins and masses of the heavy quarks, as well as of the Dirac structure of the current. Hence, Ξ can only be a function of the meson velocities and of the renormalization scale ?. Lorentz covariance and parity invariance imply that Ξ must transform as a scalar with even parity. This allows the decomposition Ξ(v, v ′, ?) = Ξ1 + Ξ2 / v + Ξ3 / v ′ + Ξ4 / v/ v′ , (193)

with coe?cients Ξi = Ξi (v · v ′ , ?). But using the projection properties of the tensor wave functions, one ?nds that under the trace Ξ(v, v ′, ?) → Ξ1 ? Ξ2 ? Ξ3 + Ξ4 ≡ ?ξ (v · v ′ , ?) . Therefore22 , M ′ (v ′ )| ? h′v′ Γ hv |M (v ) = ?ξ (v · v ′ , ?) Tr M (v ′ ) Γ M(v ) .

′

(194)

(195)

The sign is chosen such that the universal form factor ξ (v · v ′ , ?) coincides with the Isgur–Wise function, which is the single form factor that describes semileptonic weak decay processes of heavy mesons in the in?nite quark mass limit. Equation (195) summarizes in a compact way the results derived in Section 3.4. Using the explicit form of the meson wave functions, one can readily recover the relations (128), (129), and (138). The only new feature is that the Isgur–Wise function depends on the renormalization scale ?. This is necessary to compensate the scale dependence of the Wilson coe?cients, which multiply the renormalized current operators in the short-distance expansion.

61

3.9. Meson Decay Form Factors One of the most important applications of heavy quark symmetry is to derive relations between the form factors parametrizing the exclusive weak decays B → D ? ν ? ? and B → D ? ν ?. A detailed theoretical understanding of these processes is necessary for a reliable determination of the element | Vcb| of the Kobayashi–Maskawa matrix. Let me start by introducing a set of six hadronic form factors hi (w ), which parametrize the relevant meson matrix elements of the ?avour-changing vector and axial vector currents V ? = c ? γ ? b and A? = c ? γ ? γ5 b: D (v ′)| V ? |B (v ) = h+ (w ) (v + v ′ )? + h? (w ) (v ? v ′ )? ,

′ D ? (v ′ , ?)| V ? |B (v ) = i hV (w ) ??ναβ ?? ν vα vβ ,

(196)

D ? (v ′ , ?)| A? |B (v ) = hA1 (w ) (w + 1) ??? ? hA2 (w ) v ? + hA3 (w ) v ′? ?? · v . Here w = v · v ′ is the velocity transfer of the mesons, and ? denotes the polarization vector of the D ? -meson. At leading order in the 1/mQ expansion, one can derive expressions for these form factors by using the short-distance expansion (169) of the ?avour-changing currents, as well as the tensor formalism outlined above. According to (177), the ? dependence of the Wilson coe?cients of any bilinear heavy quark current can be factorized into a universal function Khh (w, ?), which is normalized at zero recoil. The ? dependence of this function has to cancel against that of the Isgur–Wise function. One can use this fact to de?ne a renormalization-group invariant Isgur–Wise form factor as Neglecting terms of order 1/mQ , one then obtains197

3

ξren (w ) ≡ ξ (w, ?) Khh(w, ?) ,

ξren(1) = 1 .

(197)

M ′ (v ′ )| J ? |M (v ) = ?ξren (w )

i=1

Ci (w ) Tr M (v ′ ) Γi M(v ) ,

(5)

′

(198)

where the Dirac matrices Γi have been de?ned in (170) and (171). It is now straightforward to evaluate the traces to ?nd w+1 h+ (w ) = C1 (w ) + C2 (w ) + C3 (w ) ξren (w ) , 2 w+1 C2 (w ) ? C3 (w ) ξren (w ) , h? (w ) = 2 hV (w ) = C1 (w ) ξren(w ) ,

5 hA1 (w ) = C1 (w ) ξren(w ) , 5 hA2 (w ) = C2 (w ) ξren(w ) , 5 5 hA3 (w ) = C1 (w ) ξren (w ) . (w ) + C 3

(199)

62

This is the correct generalization of (129) and (138) in the presence of short-distance corrections. The fact that, to leading order in 1/mQ , the meson form factors are given in terms of a single universal function ξren (w ) was the discovery of Isgur and Wise13 . Short-distance QCD corrections a?ect the form factors in a calculable way. Their e?ects are contained in the various combinations of short-distance coe?cients, which can be evaluated using the numerical results given in Table 1. 3.10. Power Corrections and Luke’s Theorem Using the covariant tensor formalism and the short-distance expansions of the e?ective Lagrangian and currents beyond the leading order in 1/mQ , one can investigate in a systematic way the structure of power corrections to the relations derived in the previous section. I have given a simple (tree-level) example of the structure of power corrections in (143). For the more complicated case of meson weak decay form factors, the analysis at order 1/mQ was performed by Luke24 . Later, shortdistance corrections were included to all orders in perturbation theory203,208 . Falk and myself have also analysed the structure of 1/m2 Q corrections for both meson 187 and baryon weak decay form factors . I shall not discuss these rather technical issues in detail, but nevertheless give you the main results. Luke has shown that, for transitions between two heavy ground-state (pseudoscalar or vector) mesons, the 1/mQ corrections can be parametrized by a set of four additional universal functions of the velocity transfer w . The most important outcome of his analysis concerns the zero recoil limit, where an analogue of the Ademollo–Gatto theorem80 can be proved. This is Luke’s theorem24 , which states that the matrix elements describing the leading 1/mQ corrections to meson decay amplitudes vanish at zero recoil. As a consequence, in the limit v = v ′ there are no terms of order 1/mQ in the hadronic matrix elements in (196). This theorem is independent of the structure of the Wilson coe?cients and thus valid to all orders in perturbation theory187,208,209 . There is considerable confusion in the literature about the implications of this result. It is often claimed that the theorem would protect any meson decay rate, or even all form factors that are normalized in the spin-?avour symmetry limit, from ?rst-order power corrections at zero recoil. However, these claims are erroneous. The reason is simple but somewhat subtle. Luke’s theorem only protects the form factors h+ and hA1 in (196), since all the others are multiplied by kinematic factors which vanish for v = v ′ . In fact, an explicit calculation shows that the 1/mQ corrections to h? , hV , hA2 , and hA3 do not vanish at zero recoil. The fact that these functions are kinematically suppressed does not imply that they could not contribute to physical decay rates. This is often overlooked. Consider, as an example, the process B → D ? ν ? in the limit of vanishing lepton mass. By angular momentum conservation, the two pseudoscalar mesons must be in a relative p-wave in order to match the helicities of the lepton pair. The amplitude is proportional to the velocity |vD | of the D -meson in the B -meson rest frame, which leads to a factor (w 2 ? 1) in 63

the decay rate. In such a situation, form factors that are kinematically suppressed can contribute210 . Indeed, the B → D ? ν ? decay rate is proportional to (w 2 ? 1)3/2 h+ (w ) ?

2 mB ? mD h? (w ) . mB + mD

(200)

Both form factors, h+ and h? , contribute with similar strength to the rate near w = 1, although h? is kinematically suppressed in (196). Consequently, the decay rate at zero recoil does receive corrections of order 1/mQ . The situation is di?erent for B → D ? ? ν ? transitions125 . Because the vector meson has spin one, the decay can proceed in an s-wave, and there is no helicity suppression near zero recoil. One ?nds that close to w = 1 the decay rate is proportional to (w 2 ? 1)1/2 |hA1 (w )|2 . Since the form factor hA1 is protected by Luke’s theorem, these transitions are ideally suited for a precision measurement of | Vcb| from an extrapolation of the momentum spectrum of the D ? -meson to zero recoil12,125 . There, the normalization of the decay rate is known in a model-independent way up to corrections of order 1/m2 Q: hA1 (1) = ηA + δ1/m2 + . . . , (201)

where ηA ? 0.99 contains the short-distance corrections [cf. (178)]. One expects higher-order power corrections to be of order δ1/m2 ? (ΛQCD /mc )2 ? 3%, but of course such a na¨ ?ve estimate could be too optimistic. For a precision measurement of | Vcb|, it is important to know the structure of such corrections in more detail. Although in principle straightforward, the analysis of 1/m2 Q corrections in HQET is a tedious enterprise187 . Three classes of corrections have to be distinguished: matrix elements of local dimension-?ve current operators, “mixed” corrections resulting from the combination of corrections to the current and to the Lagrangian, and corrections from one or two insertions of operators from the e?ective Lagrangian into matrix elements of the leading-order currents. Within these 2 classes, one can distinguish corrections proportional to 1/m2 b , 1/mc , or 1/mb mc . More than thirty new universal functions are necessary to parametrize the secondorder power corrections to meson form factors. When radiative corrections are neglected, eleven combinations of these functions contribute to the hadronic form factors hi (w ). The general results greatly simplify at zero recoil, however. There the equation of motion and the Ward identities of the e?ective theory can be used to prove that matrix elements of local dimension-?ve current operators, as well as matrix elements of time-ordered products containing a dimension-four current and an insertion from the e?ective Lagrangian, can be expressed in terms of only two parameters λ1 and λ2 , which are related to the 1/mQ corrections to the physical meson masses. Moreover, the conservation of the ?avour-diagonal vector current in the full theory forces certain combinations of the universal functions to vanish at zero recoil. The consequence is that whenever a form factor is protected by Luke’s theorem, the structure of second-order power corrections at zero recoil becomes rather simple. The quantity δ1/m2 in (201) is of this type. A careful estimate based on some mild model assumptions gives187,211 ?3% < δ1/m2 < ?1%. One can 64

combine this result with (201) to obtain one of the most important, and certainly most precise predictions of HQET: hA1 (1) = 0.97 ± 0.04 . 3.11. Properties of the Isgur–Wise Function The establishment of heavy quark symmetry as an exact limit of the strong interactions enables one to derive approximate relations between decay amplitudes, and normalization conditions for certain form factors, which are similar to the relations and normalization conditions that can be derived for Goldstone-boson scattering amplitudes from the low-energy theorems of current algebra. HQET provides the framework for a systematic investigation of the corrections to the limit of an exact spin-?avour symmetry. The output of such a model-independent analysis is a short-distance expansion of decay amplitudes, in which the dependence on the heavy quark masses is explicit. At each order in the 1/mQ expansion, the longdistance physics is parametrized by a minimal set of universal form factors, which are independent of the heavy quark masses. As presented so far, the analysis is completely model-independent. Since hadronic decay processes are of a genuine nonperturbative nature, however, it is clear that predictions that can be made based on symmetries only are limited. In particular, very little can be said on general grounds about the properties of the form factors of the e?ective theory. But there is a lot of information contained in these functions. Much like the hadron structure functions probed in deep-inelastic scattering, they are fundamental nonperturbative quantities in QCD, which describe the properties of the light degrees of freedom in the background of the colour ?eld provided by the heavy quarks. Since a static colour source is the most direct way to probe the strong interactions of quarks and gluons at large distances, a theoretical understanding of the universal form factors would not only enlarge the predictive power of HQET, but would also teach us in a very direct way about the nonperturbative nature of the strong interactions. The leading-order Isgur–Wise function ξ (w ) plays a central role in the description of the weak decays of heavy mesons. It contains the long-distance physics associated with the strong interactions of the light degrees of freedom and cannot be calculated from ?rst principles. Nevertheless, some important properties of this function can be derived on general grounds, such as its normalization at zero recoil, which is a consequence of current conservation. According to (128), the Isgur–Wise function is the elastic form factor of a ground-state heavy meson in the limit where power corrections are negligible. As such, ξ (w ) must be a monotonically decreasing function of the velocity transfer w = v · v ′ , which is analytic in the cut w -plane with a branch point at w = ?1, corresponding to the threshold q 2 = 4m2 Q for heavy quark pair production. However, being obtained from a limiting procedure, the Isgur–Wise function can have stronger singularities than the physical elastic form factor. In fact, the short-distance corrections contained in the function Khh (w, ?) 65 (202)

lead to an essential singularity at w = ?1 in the renormalized Isgur–Wise function de?ned in (197). When using a phenomenological parametrization of the universal form factor, one should incorporate the above properties. Some legitimate forms suggested in the literature are: ξBSW (w ) = 2 exp w+1

2?2

? (2?2 ? 1)

w?1 , w+1

ξISGW (w ) = exp ? ?2 (w ? 1) , ξpole (w ) = 2 w+1 . (203)

The ?rst function is the form factor derived by Rieckert and myself210 from an analysis of the BSW model81 , the second one corresponds to the ISGW model83 , and the third one is a pole-type ansatz. Of particular interest is the behaviour of the Isgur–Wise function close to zero recoil, which is determined by the slope parameter ?2 > 0 de?ned by ξ ′ (1) = ??2 , so that ξ (w ) = 1 ? ?2 (w ? 1) + O [(w ? 1)2 ] . (204)

It is important to realize that the kinematic region accessible in semileptonic decays is small (1 < w < 1.6). As long as ?2 is the same, di?erent functional forms of ξ (w ) will give similar results. A precise knowledge of the slope parameter would thus basically determine the Isgur–Wise function in the physical region. Bjorken has shown that ?2 is related to the form factors of transitions of a ground-state heavy meson into excited states, in which the light degrees of freedom + 3+ carry quantum numbers j P = 1 or 2 , by a sum rule which is an expression of 2 quark–hadron duality: In the in?nite mass limit, the inclusive sum of the probabilities for decays into hadronic states is equal to the probability for the free quark transition. If one normalizes the latter probability to unity, the sum rule has the form204,212,213 1= w+1 | ξ (w )| 2 + 2 + (w ? 1) 2

2 m

l

|ξ (l) (w )|2

n

|τ1/2 (w )|2 + (w + 1)2

(m)

|τ3/2 (w )|2 (205)

(n)

+ O [(w ? 1) ] ,

where l, m, n label the radial excitations of states with the same spin-parity quantum numbers. The sums are understood in a generalized sense as sums over discrete states and integrals over continuum states. The terms on the right-hand side of the sum rule in the ?rst line correspond to transitions into states with “brown muck” ? . The ground state gives a contribution proportional to quantum numbers j P = 1 2 66

the Isgur–Wise function, and excited states contribute proportionally to analogous functions ξ (l) (w ). Because at zero recoil these states must be orthogonal to the ground state, it follows that ξ (l) (1) = 0, and one can conclude that the corresponding contributions to (205) are of order (w ? 1)2 . The contributions in the second line + + or 3 . Because of the change in correspond to transitions into states with j P = 1 2 2 parity, these are p-wave transitions. The amplitudes are proportional to the velocity |vf | = (w 2 ?1)1/2 of the ?nal state in the rest frame of the initial state, which explains the suppression factor (w ? 1) in the decay probabilities. The functions τj (w ) are the analogues of the Isgur–Wise function for these transitions. Transitions into excited states with quantum numbers other than the above proceed via higher partial waves and are suppressed by at least a factor (w ? 1)2 . For w = 1, eq. (205) reduces to the normalization condition for the Isgur–Wise function. The Bjorken sum rule is obtained by expanding in powers of (w ? 1) and keeping terms of ?rst order. Taking into account the de?nition of the slope parameter ?2 in (204), one ?nds that204,212 ?2 = 1 + 4

m

|τ1/2 (w )|2 + 2

(m)

n

|τ3/2 (w )|2 >

(n)

1 . 4

(206)

(w +1) of the ?rst term in (205) Notice that the lower bound is due to the prefactor 1 2 and is of purely kinematic origin. In the analogous sum rule for ΛQ -baryons, this factor is absent, and consequently the slope parameter of the baryon Isgur–Wise function is only subject to the trivial constraint210,214 ?2 > 0. Based on various model calculations, there is a general belief that the contributions of excited states in the Bjorken sum rule are sizeable, and that ?2 is substantially larger than 1/4. For instance, Blok and Shifman have estimated the contributions of the lowest-lying excited states to (206) using QCD sum rules and ?nd that215 0.35 < ?2 < 1.15. The experimental observation that semileptonic B decays into excited D ?? -mesons have a large branching ratio of about216 2.5% gives further support to the importance of such contributions. Voloshin has derived another sum rule involving the form factors for transitions into excited states, which is the analogue of the “optical sum rule” for the dipole scattering of light in atomic physics. It reads217 mM ? mQ (m) (m) (n) (n) = E1/2 |τ1/2 (w )|2 + 2 E3/2 |τ3/2 (w )|2 , (207) 2 m n where Ej are the excitation energies relative to the mass mM of the ground-state heavy meson. The important point is that one can combine this relation with the Bjorken sum rule to obtain an upper bound for the slope parameter ?2 : 1 mM ? mQ , (208) ?2 < + 4 2Emin where Emin denotes the minimum excitation energy. In the quark model, one expects§ that Emin ? mM ? mQ , and one may use this as an estimate to obtain ?2 < 0.75.

Strictly speaking, the lowest excited “state” contributing to the sum rule is D + π , which has an

§

67

The above discussion of the sum rules ignores renormalization e?ects. However, both the slope parameter ?2 in (206) and the heavy quark mass mQ in (207) are renormalization-scheme dependent quantities. Although there exist some qualitative ideas of how to account for the ? dependence of the various parameters212,217 , there is currently no known way to include renormalization e?ects quantitatively. One should therefore consider the bounds on ?2 as somewhat uncertain. To account for this, I relax the upper bound derived from the Voloshin sum rule and conclude that 0.25 < ?2 < 1.0 , (209) where it is expected that the actual value is close to the upper bound. Recently, de Rafael and Taron claimed to have derived an upper bound ?2 < 0.48 from general analyticity arguments218 . If true, this had severely constrained the form of the Isgur–Wise function near zero recoil. It took quite some time to become clear what went wrong with their derivation: The e?ects of resonances below the threshold for heavy meson pair production invalidate the argument219 ? 222 . It is possible to derive a new bound, which takes into account the known properties of the Υ-states223 . However, it is too loose to be of any phenomenological relevance, and one is thus left with the sum rule result (209).

Figure 16: QCD sum rule prediction for the Isgur–Wise function in the kinematic region accessible in semileptonic decays93 . The dashed lines indicate the bounds on the slope at v · v ′ = 1 derived by Bjorken and Voloshin [see (209)]. Any sensible calculation of the Isgur–Wise function has to respect these bounds. As an example, I show in Fig. 16 the result of a QCD sum rule analysis of the

excitation energy spectrum with a threshold at mπ . However, one expects that this spectrum is broad, so that this contribution will not invalidate the upper bound for ?2 derived here.

68

renormalized universal form factor93 ? 95 . Similar predictions have been obtained for the functions that describe the 1/mQ corrections to the meson form factors96 . 3.12. Model-Independent Determination of | Vcb| One of the most important results of HQET is the prediction (202) of the normalization of the hadronic form factor hA1 at zero recoil. It can be used to obtain a model-independent measurement of the element | Vcb| of the Kobayashi–Maskawa matrix. The semileptonic decay B → D ? ? ν ? is ideally suited for this purpose125 . Experimentally, this is a particularly clean mode, since the reconstruction of the D ? -meson mass provides a powerful rejection against background. From the theoretical point of view, it is ideal since the decay rate at zero recoil is protected by Luke’s theorem against ?rst-order power corrections in 1/mQ . In terms of the hadronic form factors de?ned in (196), one ?nds

w →1

lim √

At zero recoil, the normalization of hA1 is known in a model-independent way with an accuracy of 4%. Ideally, then, one can extract | Vcb| with a theoretical uncertainty of well below 10% from an extrapolation of the spectrum to w = 1. Presently, the proposal to measure | Vcb| close to zero recoil poses quite a challenge to the experimentalists. First, there is the fact that the decay rate vanishes at zero recoil because of phase space. Therefore the statistics gets worse as one tries to measure close to w = 1. However, I do not believe that this will be an important √ limitation of the method. The phase-space suppression is proportional to w 2 ? 1 and is in fact a rather mild one. When going from the endpoint wmax ? 1.5 down to w = 1.05, the change in the statistical error in | Vcb|2 due to the variation of the phase-space factor is not even a factor of 2. A more serious problem is related to the fact that, for experiments working on the Υ(4s) resonance, the zero-recoil limit corresponds to a situation where both the B - and the D ? -mesons are approximately at rest in the laboratory. Then the pion in the subsequent decay D ? → D π is very soft and can hardly be detected. Thus, the present experiments have to make cuts which disfavour the zero recoil region, leading to large systematic uncertainties for values of w smaller than about 1.15. This second problem would be absent at an asymmetric B -factory, where the rest frame of the parent B -meson is boosted relative to the laboratory frame. In view of these di?culties, one presently has to rely on an extrapolation over a wide range in w to obtain a measurement of | Vcb|. In general, the di?erential decay rate can be written as25,125 √ dΓ(B → D ? ? ν ?) G2 F 2 2 3 ?) m ? η w 2 ? 1 (w + 1)2 = ( m ? m B D A D dw 48π 3 4w 1 ? 2wr + r 2 | Vcb |2 ξ 2 (w ) , (211) × 1+ w + 1 (1 ? r )2 69

2 G2 dΓ(B → D ? ? ν ?) F | Vcb | 2 = (mB ? mD? )2 m3 D ? |hA1 (1)| . 3 2 d w 4 π w ?1

1

(210)

where ηA = 0.99 is the short-distance correction to the form factor hA1 (w ) at zero recoil, and r = mD? /mB . Equation (211) is written in such a way that the deviations from the heavy quark symmetry limit are absorbed into the form factor ξ (w ), which in the absence of symmetry-breaking corrections would be the Isgur–Wise function. Since everything except | Vcb| and ξ (w ) is known, a measurement of the di?erential decay rate is equivalent to a measurement of the product | Vcb | ξ (w ). However, theory predicts the normalization of ξ (w ) at zero recoil:

?1 ξ (1) = ηA hA1 (1) = 1 + δ1/m2 = 0.98 ± 0.04 ,

(212)

where the uncertainty comes from power corrections of order 1/m2 Q . Using this information, | Vcb| and ξ (w ) can be obtained separately from a measurement of the di?erential decay rate. I have applied this strategy for the ?rst time125 to the combined sample224 of the data on B 0 → D ?+ ? ν ? decays collected until 1989 by the ARGUS and CLEO collaborations. For the extrapolation to zero recoil, I used the parametrizations given in (203) for the function ξ (w ), treating its slope at zero recoil as a free parameter. The result obtained for | Vcb| was | Vcb| τB0 1.5 ps

1/2

= 0.039 ± 0.006 .

(213)

Since this original analysis the data have changed. In particular, the branching ratio for D ?+ → D 0 π + has increased from52 55% to225 68%. This lowers the decay rate for B 0 → D ?+ ? ν ?, and correspondingly decreases | Vcb| by 10%. However, the new data recently reported by the ARGUS216 and CLEO126 collaborations give a larger branching ratio than the old data, indicating that further changes in the analysis must have taken place. It is thus not possible to simply rescale the result for | Vcb| given above. In Fig. 17, I show the new ARGUS data216 for the product | Vcb| ξ (w ). From an unconstrained ?t using again the parametrizations in (203), the following value is obtained: τB0 1/2 | Vcb| = 0.049 ± 0.008 . (214) 1.5 ps However, the ?t gives very large values for the slope parameter ?2 , between 1.9 and 2.3. Since the slope of ξ (w ) agrees with the slope of the Isgur–Wise function up to power corrections of order ΛQCD /mc , I believe that such large values cannot be tolerated in view of the Voloshin sum rule. In Fig. 17, I therefore show a ?t to the data which uses the pole ansatz in (203) with the constraint that ?2 ≤ 1. This leads to a signi?cantly smaller value: | Vcb| τB0 1.5 ps

1/2

= 0.037 ± 0.006 . 70

(215)

0.06

| Vcb| ξ(w)

0.04

0.02

0

1

1.1

1.2

1.3 w

1.4

1.5

Figure 17: ARGUS data216 for the product | Vcb| ξ (w ) as a function of the recoil w = v · v ′ , assuming τB0 = 1.5 ps. Using the normalization condition (212), | Vcb| is obtained from an extrapolation to w = 1. The ?t curve is explained in the text. Very recently, the CLEO collaboration has reported new results for the recoil spectrum with higher statistics126 . They have applied tight cuts in order to reduce the systematic errors. These data are shown in Fig. 18. The curves correspond to the results of a ?t using the various functions in (203). Without any constraint on the slope at zero recoil, one obtains | Vcb| τB0 1.5 ps

1/2

= 0.039 ± 0.006 ,

(216)

together with values 1.0 ± 0.4 < ?2 < 1.2 ± 0.7. It is reassuring that the result for the slope parameter is in agreement with the bound derived from the Voloshin sum rule. Because of the careful analysis of systematic uncertainties, I consider these results to be the best numbers for | Vcb| and ?2 that are currently available. The model-independent determination of | Vcb| is one of the many applications of heavy quark symmetry that have been explored over the last years. I hope that the above presentation gives you a ?avour of the potential of the HQET formalism. For a more comprehensive discussion, I refer to my review paper25 . 3.13. Exercises ? Go through the derivation of the e?ective Lagrangian of HQET and prove eqs. (104)–(113). 71

0.06

| Vcb| ξ(w)

0.04

0.02

0

1

1.1

1.2

1.3 w

1.4

1.5

Figure 18: CLEO data126 for the product | Vcb| ξ (w ) as a function of the recoil w . ? Derive the Feynman rules depicted in Fig. 12 from the e?ective Lagrangian (114). ? Convince yourself that (138) is the correct covariant expression for the pseudoscalar-to-vector transition matrix element. ? Derive the one-loop expression given in (149) for the wave-function renormalization constant Zh . ? By performing the appropriate Dirac traces in (198), derive the expressions given in (199) for the meson form factors hi (w ). ? Derive the Bjorken sum rule (206) by taking the limit w → 1 in (205).

72

Concluding Remarks In these lectures, I have presented an introduction to the theory and phenomenology of heavy quark masses and mixing angles, as well as to the exciting new ?eld of heavy quark spin-?avour symmetry. The synthesis of these topics is interesting in itself, in that it combines rather mature subjects with one of the most active ?elds of research in particle physics. As far as quark masses and mixing angles are concerned, the theoretical concepts have been set many years ago. Yet new theoretical tools have to be developed to determine these parameters more and more accurately. Heavy quark e?ective ?eld theory, on the other hand, has received broad attention in recent years only, although some ideas of a spin-?avour symmetry for heavy quarks had been around for much longer. These new developments start to have an impact on the precision with which one can determine some elements of the Kobayashi–Maskawa matrix. I have discussed this in detail for the extraction of | Vcb | from the semileptonic decays B → D ? ? ν ?. Besides reviewing the status of the theoretical developments, my purpose was to convince you that B -physics is a rich and diverse area of research. Presently, this ?eld is characterized by a very fruitful interplay between theory and experiments, which has led to many signi?cant discoveries and developments on both sides. Heavy quark physics has the potential to determine important parameters of the ?avour sector of the electroweak theory, and at the same time provides stringent tests of the standard model. On the other hand, weak decays of hadrons containing a heavy quark are an ideal laboratory to study the nature of nonperturbative phenomena in QCD. When I presented these lectures in the summer of 1993, the future of B -physics was uncertain. Today, the prospects for further signi?cant developments in this ?eld look rather promising. With the approval of the ?rst asymmetric B -factory at SLAC, and with ongoing B -physics programs at the existing facilities at Cornell, Fermilab, and LEP, there are clearly B eautiful times ahead of us!

Acknowledgements: It is a great pleasure to thank the organizers of TASI-93 for the invitation to present these lectures and for providing a stimulating and relaxing atmosphere, which helped to initiate many physics discussions. Coming to TASI-93 also gave me the opportunity to meet, for the ?rst time, my collaborator Zoltan Ligeti, and to complete a piece of work with him. As far as fun and sports are concerned, I will remember an honest attempt to climb Longs Peak with Zoltan, and a most pleasant day of rock climbing in the Boulder area with Tom DeGrand. Writing lecture notes, on the other hand, is a painful and time-consuming experience. It is only because of the continuous support and queries of Stuart Raby that I ?nally ?nished this manuscript. I very much appreciated his patience with me. Many discussions with colleagues helped me to prepare the lectures and had an impact on the form of these notes. I am particularly grateful to my collaborators Adam Falk, Zoltan Ligeti, Yossi Nir, Michael Luke, and Thomas Mannel. During 73

the school, I have enjoyed discussions with many of the students, as well as with Andy Cohen, Chris Hill, and Peter Lepage. Before coming to Boulder, talking to Michael Peskin was an invaluable asset in preparing the lectures. In addition, I would like to thank David Cassel and Persis Drell for keeping me updated about the CLEO analysis of semileptonic B -decays. Finally, it is my great pleasure to thank Suzy Vascotto for her careful reading of the manuscript, and Maria Girone for teaching me PAW.

74

References 1. 2. 3. 4. 5. 6. 7. 8. 9. 10. 11. 12. 13. 14. 15. 16. 17. 18. 19. 20. 21. 22. 23. 24. 25. 26. 27. 28. 29. 30. 31. 32. 33. H. Albrecht et al. (ARGUS collaboration), Phys. Lett. B 192, 245 (1987). C. Albajar et al. (UA1 collaboration), Phys. Lett. B 186, 247 (1987). H. Albrecht et al. (ARGUS collaboration), Phys. Lett. B 234, 409 (1990). R. Fulton et al. (CLEO collaboration), Phys. Rev. Lett. 64, 16 (1990). H. Albrecht et al. (ARGUS collaboration), Phys. Lett. B 255, 297 (1991). J. Bartelt et al. (CLEO collaboration), Phys. Rev. Lett. 71, 4111 (1993). R. Ammar et al. (CLEO collaboration), Phys. Rev. Lett. 71, 674 (1993). E.V. Shuryak, Phys. Lett. B 93, 134 (1980); Nucl. Phys. B 198, 83 (1982). J.E. Paschalis and G.J. Gounaris, Nucl. Phys. B 222, 473 (1983); F.E. Close, G.J. Gounaris, and J.E. Paschalis, Phys. Lett. B 149, 209 (1984). S. Nussinov and W. Wetzel, Phys. Rev. D 36, 130 (1987). M.B. Voloshin and M.A. Shifman, Yad. Fiz. 45, 463 (1987) [Sov. J. Nucl. Phys. 45, 292 (1987)]. M.B. Voloshin and M.A. Shifman, Yad. Fiz. 47, 801 (1988) [Sov. J. Nucl. Phys. 47, 511 (1988)]. N. Isgur and M.B. Wise, Phys. Lett. B 232, 113 (1989); 237, 527 (1990). E. Eichten and F. Feinberg, Phys. Rev. D 23, 2724 (1981). W.E. Caswell and G.P. Lepage, Phys. Lett. B 167, 437 (1986). E. Eichten, in: Field Theory on the Lattice, edited by A. Billoire et al., Nucl. Phys. B (Proc. Suppl.) 4, 170 (1988). G.P. Lepage and B.A. Thacker, in: Field Theory on the Lattice, edited by A. Billoire et al., Nucl. Phys. B (Proc. Suppl.) 4, 199 (1988). H.D. Politzer and M.B. Wise, Phys. Lett. B 206, 681 (1988); 208, 504 (1988). E. Eichten and B. Hill, Phys. Lett. B 234, 511 (1990); 243, 427 (1990). B. Grinstein, Nucl. Phys. B 339, 253 (1990). H. Georgi, Phys. Lett. B 240, 447 (1990). A.F. Falk, H. Georgi, B. Grinstein, and M.B. Wise, Nucl. Phys. B 343, 1 (1990). A.F. Falk, B. Grinstein, and M.E. Luke, Nucl. Phys. B 357, 185 (1991). M.E. Luke, Phys. Lett. B 252, 447 (1990). M. Neubert, SLAC preprint SLAC-PUB-6263 (1993), to appear in Phys. Rep. R. Tarrach, Nucl. Phys. B 183, 384 (1981). W.A. Bardeen, A.J. Buras, D.W. Duke, and T. Muta, Phys. Rev. D 18, 3998 (1978); A.J. Buras, Rev. Mod. Phys. 52, 199 (1980). G. ’t Hooft, Nucl. Phys. B 61, 455 (1973). M.A. Shifman, A.I. Vainshtein, and V.I. Zakharov, Nucl. Phys. B 147, 385 (1979); 147, 448 (1979). A.F. Falk, M. Neubert, and M.E. Luke, Nucl. Phys. B 388, 363 (1992). M. Neubert, Phys. Rev. D 46, 3914 (1992). G. ’t Hooft and M. Veltman, Nucl. Phys. B 44, 189 (1972). C.G. Bollini and J.J. Giambiagi, Phys. Lett. B 40, 566 (1972); Nuovo Cim. 75

34. 35. 36. 37. 38. 39. 40. 41. 42. 43. 44. 45. 46. 47. 48. 49. 50. 51. 52. 53. 54. 55. 56. 57. 58. 59. 60. 61. 62. 63.

B 12, 20 (1972). D.J. Gross, in: Methods in Field Theory, edited by C.R. Balian and J. ZinnJustin (North-Holland, Amsterdam, 1976), p. 141. N. Gray, D.J. Broadhurst, W. Grafe, and K. Schilcher, Z. Phys. C 48, 673 (1990). O.V. Tarasov, JINR preprint P2-82-900 (1982), unpublished. S.G. Gorishny, A.L. Kataev, and S.A. Larin, Phys. Lett. B 135, 457 (1984). D.J. Gross and F. Wilczek, Phys. Rev. Lett. 30, 1343 (1973). H.D. Politzer, Phys. Rev. Lett. 30, 1346 (1973). W.E. Caswell, Phys. Rev. Lett. 33, 244 (1974). D.R.T. Jones, Nucl. Phys. B 75, 531 (1974). A.A. Belavin and A.A. Migdal, JETP Lett. 19, 181 (1974) [ZhETF Pis. Red. 19, 317 (1974)]. E.G. Floratos, D.A. Ross, and C.T. Sachrajda, Nucl. Phys. B 129, 66 (1977) [E: 139, 545 (1978)]. A.J. Buras, M. Jamin, M.E. Lautenbacher, and P.H. Weisz, Nucl. Phys. B 370, 69 (1992). T. Appelquist and J. Carazzone, Phys. Rev. D 11, 2856 (1975). B. Pendleton and G.G. Ross, Phys. Lett. B 98, 291 (1981). C.T. Hill, Phys. Rev. D 24, 691 (1981); C.T. Hill, C.N. Leung, and S. Rao, Nucl. Phys. B 262, 517 (1985). J. Bagger, S. Dimopoulos, and E. Mass? o, Phys. Rev. Lett. 55, 920 (1985). M. Carena et al., Nucl. Phys. B 369, 33 (1992); M. Carena, S. Pokorski, and C.E.M. Wagner, Nucl. Phys. B 406, 59 (1993). V. Barger, M.S. Berger, and P. Ohmann, Phys. Rev. D 47, 1093 (1993). W.A. Bardeen, C.T. Hill, and M. Lindner, Phys. Rev. D 41, 1647 (1990). K. Hikasa et al., Review of Particle Properties, Phys. Rev. D 45, Part II (1992). B. Guberina, R. Meckbach, R.D. Peccei, and R. R¨ uckl, Nucl. Phys. B 184, 476 (1981). L.J. Reinders, Phys. Rev. D 38, 947 (1988); L.J. Reinders, H.R. Rubinstein, and S. Yazaki, Phys. Rep. 127, 1 (1985). V.A. Novikov et al., Phys. Rep. 41, 1 (1978). M.B. Voloshin, Yad. Fiz. 29, 1368 (1979) [Sov. J. Nucl. Phys. 29, 703 (1979)]. M.B. Voloshin and Yu.M. Zaitsev, Usp. Fiz. Nauk 152, 361 (1987) [Sov. Phys. Usp. 30, 553 (1987)]. S. Narison, Phys. Lett. B 210, 238 (1988); 308, 305 (1993). E. Bagan, P. Ball, V.M. Braun, and H.G. Dosch, Phys. Lett. B 278, 457 (1992). M. Neubert, Phys. Rev. D 46, 1076 (1992). I.I. Bigi, M.A. Shifman, N.G. Uraltsev, and A.I. Vainshtein, CERN preprint CERN-TH.7171/94 (1994). M. Beneke and V.M. Braun, Munich preprint MPI-PhT/94-9 (1994). G. ’t Hooft, in: The Whys of Subnuclear Physics, Proceedings of the 15th 76

64. 65. 66. 67. 68. 69.

70. 71. 72. 73. 74. 75. 76. 77. 78. 79. 80. 81. 82. 83. 84. 85. 86. 87. 88.

89. 90. 91.

International School on Subnuclear Physics, Erice, Sicily, 1977, edited by A. Zichichi (Plenum Press, New York, 1979), p. 943. B. Lautrup, Phys. Lett. B 69, 438 (1977). G. Parisi, Phys. Lett. B 76, 65 (1978); Nucl. Phys. B 150, 163 (1979). F. David, Nucl. Phys. B 234, 237 (1984); 263, 637 (1986). A.H. Mueller, Nucl. Phys. B 250, 327 (1985); Phys. Lett. B 308, 355 (1993). A.J. Buras and M.K. Harlander, in: Heavy Flavours, edited by A.J. Buras and M. Lindner (World Scienti?c, Singapore,1992), p. 58. Y. Nir, in: The Third Family and the Physics of Flavour, Proceedings of the 20th SLAC Summer Institute on Particle Physics, Stanford, California, 1992, edited by L. Vassilian (SLAC Report No. 412, Stanford, 1993), p. 81; Y. Nir and H.R. Quinn, Annu. Rev. Nucl. Part. Sci. 42, 221 (1992). N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963). M. Kobayashi and K. Maskawa, Prog. Theor. Phys. 49, 652 (1973). C. Jarlskog, in: CP Violation, edited by C. Jarlskog (World Scienti?c, Singapore, 1989), p. 3. C. Jarlskog, Phys. Rev. Lett. 55, 1039 (1985); Z. Phys. C 29, 491 (1985). L.L. Chau and W.-Y. Keung, Phys. Rev. Lett. 53, 1802 (1984). L. Wolfenstein, Phys. Rev. Lett. 51, 1945 (1983). A. Sirlin and R. Zucchini, Phys. Rev. Lett. 57, 1994 (1986); A. Sirlin, Phys. Rev. D 35, 3423 (1987). W. Jaus and G. Rasche, Phys. Rev. D 35, 3420 (1987). H. Leutwyler and M. Roos, Z. Phys. C 25, 91 (1984). J.F. Donoghue, B.R. Holstein, and S.W. Klimt, Phys. Rev. D 35, 934 (1987). M. Ademollo and R. Gatto, Phys. Rev. Lett. 13, 264 (1964). M. Wirbel, B. Stech and M. Bauer, Z. Phys. C 29, 637 (1985); M. Bauer, B. Stech, and M. Wirbel, Z. Phys. C 34, 103 (1987). B. Grinstein, M.B. Wise, and N. Isgur, Phys. Rev. Lett. 56, 298 (1986). N. Isgur, D. Scora, B. Grinstein, and M.B. Wise, Phys. Rev. D 39, 799 (1989). J.G. K¨ orner and G.A. Schuler, Z. Phys. C 38, 511 (1988) [E: 41, 690 (1989)]; 46, 93 (1990). S.J. Brodsky and G.P. Lepage, Phys. Rev. D 22, 2157 (1980). M. Suzuki, Nucl. Phys. B 258, 553 (1985). T. Altomari and L. Wolfenstein, Phys. Rev. Lett. 58, 1583 (1987). T.M. Aliev and V.L. Eletskij, Yad. Fiz. 38, 1537 (1983) [Sov. J. Nucl. Phys. 38, 936 (1983)]; T.M. Aliev, V.L. Eletskij, and Y.I. Kogan, Yad. Fiz. 40, 823 (1984) [Sov. J. Nucl. Phys. 40, 527 (1984)]. A.A. Ovchinnikov, V.A. Slobodenyuk, Z. Phys. C 44, 433 (1989); A.A. Ovchinnikov, Yad. Fiz. 50, 1433 (1989) [Sov. J. Nucl. Phys. 50, 891 (1989)]. V.N. Baier and A.G. Grozin, Z. Phys. C 47, 669 (1990). P. Ball, V.M. Braun, H.G. Dosch, and M. Neubert, Phys. Lett. B 259, 481 (1991); P. Ball, V.M. Braun, and H.G. Dosch, Phys. Rev. D 44, 3567 (1991); 77

92.

93. 94. 95. 96. 97. 98. 99. 100. 101. 102. 103. 104. 105. 106. 107. 108. 109. 110. 111. 112.

113. 114. 115. 116.

Phys. Lett. B 273, 316 (1991). P. Ball, Phys. Lett. B 281, 133 (1992); in: Perturbative QCD and Hadronic Interactions, Proceedings of the 27th Rencontres de Moriond, edited by J. Tran Thanh Van, Gif-sur-Yvette, France, 1992 (Editions Fronti` eres, Gifsur-Yvette, 1992), p. 343. M. Neubert, Phys. Rev. D 45, 2451 (1992); 47, 4063 (1993). M. Neubert, V. Rieckert, B. Stech, and Q.P. Xu, in: Heavy Flavours, edited by A.J. Buras and M. Lindner (World Scienti?c, Singapore, 1992), p. 286. A.V. Radyushkin, Phys. Lett. B 271, 218 (1991). M. Neubert, Z. Ligeti, and Y. Nir, Phys. Lett. B 301, 101 (1993); Phys. Rev. D 47, 5060 (1993); 49, 1302 (1994). J.H. Christenson, J.W. Cronin, V.L. Fitch, and R. Turlay, Phys. Rev. Lett. 13, 138 (1964); Phys. Rev. B 140, 74 (1965). N. Isgur and M.B. Wise, Phys. Rev. D 42, 2388 (1990). C.O. Dib and F. Vera, Phys. Rev. D 47, 3938 (1993). G. Burdman, Z. Ligeti, M. Neubert, and Y. Nir, Phys. Rev. D 49, 2331 (1994). A. Ali and E. Pietarinen, Nucl. Phys. B 154, 519 (1979). G. Corb` o, Nucl. Phys. B 212, 99 (1983); N. Cabibbo, G. Corb` o, and L. Maiani, Nucl. Phys. B 155, 93 (1979). M. Jezabek and J.H. K¨ uhn, Nucl. Phys. B 320, 20 (1989). G. Altarelli et al., Nucl. Phys. B 208, 365 (1982). A. Bareiss and E.A. Paschos, Nucl. Phys. B 327, 353 (1989); C.H. Jin, W.F. Palmer, and E.A. Paschos, Dortmund preprint DO-TH-93/21 (1993). M. Neubert, Phys. Rev. D 49, 3392 (1994). M. Neubert, CERN preprint CERN-TH.7113/93 (1993), to appear in Phys. Rev. D 49, no. 9 (1994). I.I. Bigi, M.A. Shifman, N.G. Uraltsev, and A.I. Vainshtein, CERN preprint CERN-TH.7129/93 (1993). A.F. Falk, E. Jenkins, A.V. Manohar, and M.B. Wise, San Diego preprint UCSD/PTH-93-38 (1993). T. Mannel and M. Neubert, CERN preprint CERN-TH.7156/94 (1994). J. Chay, H. Georgi, and B. Grinstein, Phys. Lett. B 247, 399 (1990). I.I. Bigi, N.G. Uraltsev, and A.I. Vainshtein, Phys. Lett. B 293, 430 (1992); I.I. Bigi, M.A. Shifman, N.G. Uraltsev, and A.I. Vainshtein, Phys. Rev. Lett. 71, 496 (1993); I.I. Bigi et al., in: Proceedings of the Annual Meeting of the Division of Particles and Fields of the American Physical Society, Batavia, Illinois, 1992, edited by C. Albright et al. (World Scienti?c, Singapore, 1993), p. 610. B. Blok, L. Koyrakh, M.A. Shifman, and A.I. Vainshtein, Phys. Rev. D 49, 3356 (1994). A.V. Manohar and M.B. Wise, Phys. Rev. D 49, 1310 (1994). A.F. Falk, M. Luke, and M.J. Savage, Phys. Rev. D 49, 3367 (1994). T. Mannel, Nucl. Phys. B 413, 396 (1994). 78

117. L. Koyrakh, Minnesota preprint TPI-MINN-93/47-T (1993). 118. A.F. Falk, Z. Ligeti, M. Neubert, and Y. Nir, CERN preprint CERNTH.7124/93 (1993), to appear in Phys. Lett. B. 119. S. Balk, J.G. K¨ orner, D. Pirjol, and K. Schilcher, Mainz preprint MZ-TH/9332 (1993). 120. M. Luke and M.J. Savage, Phys. Lett. B 321, 88 (1994). 121. G.P. Yost et al., Review of Particle Properties, Phys. Lett. B 204 (1988). 122. S. Stone, to appear in the 2nd edition of: B Decays, edited by S. Stone (World Scienti?c, Singapore, 1991). 123. M. Danilov, to appear in: Proceedings of the International Europhysics Conference on High Energy Physics, Marseille, France, July 1993, edited by J. Carr and M. Perrottet (Editions Fronti` eres, Gif-sur-Yvettes). 124. T. Hessing, to appear in: Electroweak Interactions and Uni?ed Theories, Proceedings of the 29th Rencontres de Moriond, M? eribel, France, March 1993. 125. M. Neubert, Phys. Lett. B 264, 455 (1991). 126. G. Crawford et al. (CLEO collaboration), to appear in: Proceedings of the 16th International Symposium on Lepton and Photon Interaction, Ithaca, New York, Aug. 1993. 127. J.D. Bjorken, presented in discussions at the Workshop on Experiments, Detectors and Experimental Areas for the Supercollider, Berkeley, California, July 1987 (Proceedings edited by R. Donaldson and M.G.D. Gilchriese, World Scienti?c, Singapore, 1988). 128. C. Jarlskog and R. Stora, Phys. Lett. B 208, 268 (1988). 129. I.I. Bigi, V.A. Khoze, N.G. Uraltsev, and A.I. Sanda, in: CP Violation, edited by C. Jarlskog (World Scienti?c, Singapore, 1989), p. 175. 130. C.O. Dib, I. Dunietz, F.J. Gilman, and Y. Nir, Phys. Rev. D 41, 1522 (1990); F.J. Gilman and Y. Nir, Annu. Rev. Nucl. Part. Sci. 40, 213 (1990). 131. C.S. Kim, J.L. Rosner, and S.-P. Yuan, Phys. Rev. D 42, 96 (1990). 132. M. Schmidtler and K.R. Schubert, Z. Phys. C 53, 347 (1992). 133. M. Lusignoli, L. Maiani, G. Martinelli, and L. Reina, Nucl. Phys. B 369, 139 (1992). 134. G.R. Harris and J.L. Rosner, Phys. Rev. D 45, 946 (1992). 135. H. Albrecht et al. (ARGUS collaboration), Z. Phys. C 55, 357 (1992). 136. J. Bartelt et al. (CLEO collaboration), Phys. Rev. Lett. 71, 1680 (1993). 137. D. Abbaneo, to appear in: Electroweak Interactions and Uni?ed Theories, Proceedings of the 29th Rencontres de Moriond, M? eribel, France, March 1993. 138. A.J. Buras, M. Jamin, and P.H. Weisz, Nucl. Phys. B 347, 491 (1990). 139. T. Inami and C.S. Lim, Progr. Theor. Phys. 65, 297 (1981); 65, 1772 (1981). 140. A.J. Buras, Phys. Rev. Lett. 46, 1354 (1981). 141. G. Buchalla, A.J. Buras, and M.K. Harlander, Nucl. Phys. B 349, 1 (1991). 142. J.F. Donoghue, E. Golowich, and B.R. Holstein, Phys. Lett. B 119, 412 (1982). 79

143. A. Pich and E. de Rafael, Phys. Lett. B 158, 477 (1985); Nucl. Phys. B 358, 311 (1991). 144. R. Decker, in: Hadronic Matrix Elements and Weak Decays, edited by A.J. Buras, J.-M. G? erard, and W. Huber, Nucl. Phys. B (Proc. Suppl.) 7a, 180 (1989), and references therein. 145. W.A. Bardeen, A.J. Buras, and J.-M. G? erard, Phys. Lett. B 211, 343 (1988). 146. M.B. Gavela et al., Phys. Lett. B 206, 113 (1988); Nucl. Phys. B 306, 677 (1988). 147. G.W. Kilcup, S.R. Sharpe, R. Gupta, and A. Patel, Phys. Rev. Lett. 64, 25 (1990). 148. C. Bernard and A. Soni, in: Lattice 89, edited by N. Cabibbo et al., Nucl. Phys. B (Proc. Suppl.) 17, 495 (1990). 149. C. Sachrajda, in: QCD-20 Years Later, Proceedings of the Workshop on QCD, Aachen, Germany, 1992, edited by P.M. Zerwas and H.A. Kastrup (World Scienti?c, Singapore, 1993), p. 668, and references therein. 150. S.R. Sharpe and A. Patel, Washington preprint UW-PT-93-1 (1993); S.R. Sharpe, to appear in: Proceedings of the 11th Symposium on Lattice Field Theory (Lattice 93), Dallas, Texas, Oct. 1993. 151. H. Burckhardt et al. (NA31 collaboration), Phys. Lett. B 206, 169 (1988). 152. J.R. Patterson et al. (E731 collaboration), Phys. Rev. Lett. 64, 1491 (1990). 153. G. Buchalla, A.J. Buras, and M.K. Harlander, Nucl. Phys. B 337, 313 (1990). 154. A.J. Buras, M. Jamin, and M.E. Lautenbacher, Nucl. Phys. B 408, 209 (1993). 155. M. Ciuchini, E. Franco, G. Martinelli, and L. Reina, Phys. Lett. B 301, 263 (1993); Rome preprint 93/913 (1993). 156. S. Bertolini, in: Higgs Particle(s): Physics Issues and Searches in HighEnergy Collisions, Proceedings of the 8th INFN Eloisatron Project Workshop, Erice, Sicily, 1989, edited by A. Ali (Plenum Press, New York, 1990), p. 243. 157. A. Ali, to appear in the 2nd edition of: B Decays, edited by S. Stone (World Scienti?c, Singapore, 1991). 158. S. Bertolini, F. Borzumati, and A. Masiero, Phys. Rev. Lett. 59, 180 (1987). 159. N.G. Deshpande et al., Phys. Rev. Lett. 59, 183 (1987). 160. B. Grinstein, R. Springer, and M.B. Wise, Phys. Lett. B 202, 138 (1988); Nucl. Phys. B 339, 269 (1990). 161. R. Grigjanis, P.J. O’Donnel, M. Sutherland, and H. Navelet, Phys. Lett. B 213, 355 (1988); 223, 239 (1989); 237, 252 (1990). 162. G. Cella, G. Curci, G. Ricciardi, and A. Vicere, Phys. Lett. B 248, 181 (1990). 163. M. Misiak, Phys. Lett. B 269, 161 (1991); Nucl. Phys. B 393, 23 (1993). 164. M. Ciuchini et al., Phys. Lett. B 316, 127 (1993). 165. A. Ali and C. Greub, Z. Phys. C 49, 431 (1991); Phys. Lett. B 259, 182 (1991); 287, 191 (1992). 166. A. Ali and T. Mannel, Phys. Lett. B 264, 447 (1991) [E: 274, 526 (1992)]. 80

167. P. Ball, Munich preprint TUM-T31-43/93 (1993). 168. S. Narison, CERN preprint CERN-TH.7166/94 (1994). 169. H. Georgi, in: Perspectives in the Standard Model, Proceedings of the Theoretical Advanced Study Institute in Elementary Particle Physics (TASI-91), Boulder, Colorado, 1991, edited by R.K. Ellis, C.T. Hill, and J.D. Lykken (World Scienti?c, Singapore, 1992), p. 589. 170. B. Grinstein, in: High Energy Phenomenology, Proceedings of the Workshop on High Energy Phenomenology, Mexico City, Mexico, 1991, edited by M.A. P? eres and R. Huerta (World Scienti?c, Singapore, 1992), p. 161. 171. N. Isgur and M.B. Wise, in: Heavy Flavours, edited by A.J. Buras and M. Lindner (World Scienti?c, Singapore, 1992), p. 234. 172. T. Mannel, in: QCD-20 Years Later, Proceedings of the Workshop on QCD, Aachen, Germany, 1992, edited by P.M. Zerwas and H.A. Kastrup (World Scienti?c, Singapore, 1993), p. 634. 173. T. Appelquist and H.D. Politzer, Phys. Rev. Lett. 34, 43 (1975). 174. M.A. Shifman, A.I. Vainshtein, and V.I. Zakharov, Nucl. Phys. B 120, 316 (1977). 175. E. Witten, Nucl. Phys. B 122, 109 (1977). 176. J. Polchinski, Nucl. Phys. B 231, 269 (1984). 177. K. Wilson, Phys. Rev. 179, 1499 (1969); Phys. Rev. D 3, 1818 (1971). 178. W. Zimmermann, Ann. of Phys. 77, 536 (1973); 77, 570 (1973). 179. G. Altarelli and L. Maiani, Phys. Lett. B 52, 351 (1974). 180. M.K. Gaillard and B.W. Lee, Phys. Rev. Lett. 33, 108 (1974). 181. F.J. Gilman and M.B. Wise, Phys. Rev. D 27, 1128 (1983). 182. T. Mannel, W. Roberts and Z. Ryzak, Nucl. Phys. B 368, 204 (1992). 183. J. Soto and R. Tzani, Phys. Lett. B 297, 358 (1992). 184. H. Georgi, B. Grinstein, and M.B. Wise, Phys. Lett. B 252, 456 (1990). 185. N. Isgur and M.B. Wise, Phys. Rev. Lett. 66, 1130 (1991). 186. D. Buskulic et al. (ALEPH collaboration), Phys. Lett. B 311, 425 (1993). 187. A.F. Falk and M. Neubert, Phys. Rev. D 47, 2965 (1993); 47, 2982 (1993). 188. G. Preparata and W.I. Weisberger, Phys. Rev. 175, 1965 (1968). 189. X. Ji and M.J. Musolf, Phys. Lett. B 257, 409 (1991). 190. D.J. Broadhurst and A.G. Grozin, Phys. Lett. B 267, 105 (1991). 191. D.J. Broadhurst, N. Gray, and K. Schilcher, Z. Phys. C 52, 111 (1991). 192. B.S. DeWitt, Phys. Rev. 162, 1195 (1967). 193. G. ’t Hooft, in: Functional and Probabilistic Methods in Quantum Field Theory, Proceedings of the 12th Winter School of Theoretical Physics, Karpacz, Poland, Acta Univ. Wratisl. 38, Vol. 1 (1975). 194. D. Boulware, Phys. Rev. D 23, 389 (1981). 195. L.F. Abbott, Nucl. Phys. B 185, 189 (1981); Acta Phys. Pol. B 13, 33 (1982). 196. H. Kluberg-Stern and J.B. Zuber, Phys. Rev. D 12, 3159 (1975). 197. M. Neubert, Phys. Rev. D 46, 2212 (1992). 198. G.P. Korchemsky and A.V. Radyushkin, Nucl. Phys. B 283, 342 (1987); Phys. Lett. B 279, 359 (1992); G.P. Korchemsky, Mod. Phys. Lett. A 4, 81

199. 200. 201. 202. 203. 204.

205. 206. 207. 208. 209. 210. 211. 212. 213. 214. 215. 216. 217. 218. 219. 220. 221. 222. 223. 224. 225.

1257 (1989). A.F. Falk and B. Grinstein, Phys. Lett. B 247, 406 (1990); 249, 314 (1990). M. Neubert, Nucl. Phys. B 371, 149 (1992). W. Kilian, P. Manakos, and T. Mannel, Phys. Rev. D 48, 1321 (1993). M. Luke and A.V. Manohar, Phys. Lett. B 286, 348 (1992). M. Neubert, Phys. Lett. B 306, 357 (1993); Phys. Rev. D 49, 1542 (1994). J.D. Bjorken, in: Results and Perspectives in Particle Physics, Proceedings of the 4th Rencontres de Physique de la Vall? e d’Aoste, La Thuile, Italy, 1990, edited by M. Greco (Editions Fronti` eres, Gif-sur-Yvette, 1990), p. 583; in: Gauge Bosons and Heavy Quarks, Proceedings of the 18th SLAC Summer Institute on Particle Physics, Stanford, California, 1990, edited by J.F. Hawthorne (SLAC Report No. 378, Stanford, 1991), p. 167. A.F. Falk, Nucl. Phys. B 378, 79 (1992). H.D. Politzer, Phys. Lett. B 250, 128 (1990). T. Mannel, W. Roberts, and Z. Ryzak, Phys. Lett. B 271, 421 (1991). M. Neubert, Nucl. Phys. B 416, 786 (1994). P. Cho and B. Grinstein, Phys. Lett. B 285, 153 (1992). M. Neubert and V. Rieckert, Nucl. Phys. B 382, 97 (1992). T. Mannel, CERN preprint CERN-TH.7162/94 (1994). N. Isgur and M.B. Wise, Phys. Rev. D 43, 819 (1991). J.D. Bjorken, I. Dunietz, and J. Taron, Nucl. Phys. B 371, 111 (1992). N. Isgur, M.B. Wise, and M. Youssefmir, Phys. Lett. B 254, 215 (1991). B. Blok and M. Shifman, Phys.

赞助商链接

- Quark Masses and Mixing Angles
- Light Quark Masses and Mixing Angles
- Quark Masses and Mixing with A4 Family Symmetry
- IFP-725-UNC Horizontal Symmetry for Quark and Squark Masses in
- Quark and Lepton Masses and Mixing Angles from Superstring Constructions
- Semileptonic $B_c$ Decay and Heavy Quark Spin Symmetry
- Heavy Quark Symmetry and the Skyrme Model
- Quark masses and mixing angles in heterotic orbifold models
- Semileptonic bc to cc Baryon Decay and Heavy Quark Spin Symmetry
- Flavour-symmetry breaking of the quark condensate and chiral corrections to the Gell-Mann-O
- Quark Masses and Mixing Angles
- Doubly Heavy Baryons, Heavy Quark-DiQuark Symmetry and NRQCD
- Heavy Quark Symmetry Violation in Semileptonic Decays of D Mesons
- New Relations and Constraints on Quark Spin-flavor Contents in Symmetry-breaking Chiral Qua
- Using Heavy Quark Spin Symmetry in Semileptonic $B_c$ Decays

更多相关文章：
更多相关标签：