21.01.2015 Views

Summary of Lecture Notes with Questions - Institute for Particle ...

Summary of Lecture Notes with Questions - Institute for Particle ...

Summary of Lecture Notes with Questions - Institute for Particle ...

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Durham University,<br />

Ephiphany Term 2011,<br />

Version <strong>of</strong> July 4, 2011<br />

Flavour Physics and Effective Theories (FPE)<br />

A brief summary <strong>of</strong> basic concepts and central results will be given.<br />

self-contained script!<br />

This is not a<br />

Dr Thorsten Feldmann<br />

IPPP Durham


Contents<br />

1 Quark Flavour in the Standard Model 1<br />

2 The Weak Effective Hamiltonian 4<br />

2.1 Quantum Corrections to b → csū . . . . . . . . . . . . . . . . . . . . . . 5<br />

2.2 Wilson Coefficients in Perturbation Theory . . . . . . . . . . . . . . . . . 9<br />

2.3 Penguins and Boxes . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12<br />

3 Flavour Transitions (on the hadron level) 16<br />

3.1 Example: ¯B → DK decays . . . . . . . . . . . . . . . . . . . . . . . . . . 16<br />

3.2 Example: B → πK decays . . . . . . . . . . . . . . . . . . . . . . . . . . 20<br />

4 Mixing and CP Violation 22<br />

4.1 K 0 - ¯K 0 Mixing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 22<br />

4.2 B 0 - ¯B 0 Mixing . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30<br />

5 Semi-Leptonic Decays 36<br />

5.1 |V cb | from exclusive decays and HQET . . . . . . . . . . . . . . . . . . . 36<br />

5.2 |V cb | from inclusive decays . . . . . . . . . . . . . . . . . . . . . . . . . . 43<br />

5.3 |V ub | from inclusive decays and SCET . . . . . . . . . . . . . . . . . . . . 47


1 Quark Flavour in the Standard Model<br />

• Quark multiplets in the SM, and the Yukawa Sector <strong>of</strong> the Lagrangian<br />

L yuk = −Y ij<br />

D<br />

¯Q i LHD j R − Y ij<br />

U<br />

¯Q i L ˜HU j R + h.c.<br />

coupling left-handed doublets and right-handed singlets <strong>with</strong> the Higgs field.<br />

• Diagonalization <strong>of</strong> the Yukawa matrix by bi-unitary trans<strong>for</strong>mations<br />

and m u = y u 〈H 0 〉 etc.<br />

Y U = V uL diag(y u , y c , y t ) V †<br />

u R<br />

, Y D = V dL diag(y d , y s , y b ) V †<br />

d R<br />

• Flavour mixing in charged weak currents (weak basis → mass basis)<br />

Ū L gW µ γ µ D L → Ū ′ LV †<br />

u L<br />

gW µ γ µ V DL D ′ L<br />

≡ (V CKM ) ij (Ū ′ L) i gW µ γ µ (D ′ L) j<br />

<strong>with</strong> unitary Cabbibo-Kobayashi-Maskawa matrix, V CKM = V †<br />

u L<br />

V dL .<br />

• CKM parameter counting:<br />

– 36=18+18 parameters in Y U and Y D<br />

– broken flavour symmetries <strong>for</strong> quark multiplets:<br />

makes 26 symmetry generators.<br />

U(3) QL × U(3) UR × U(3) DR /U(1) B<br />

– leaves 10 physical parameters: 6 masses, 3 angles, 1 CP-violating phase<br />

• interactions <strong>with</strong> γ and Z 0 unaffected by unitary rotations: no flavour-changing<br />

neutral currents (FCNCs) in the SM (at tree-level)<br />

• FCNCs can be induced through quantum loops. Example: s → dγ (diagram, see<br />

lecture), which involves the CKM factors V is Vid ∗ <strong>with</strong> i = u, c, t.<br />

• W ± only couple to left-handed fields: Flavour transitions <strong>with</strong> right-handed quarks<br />

q R suppressed by m q /m W .<br />

• GIM mechanism (Glashow-Iliopoulos-Maiani): For equal quark masses in FCNC<br />

loops, the contributions would add up to zero, since ∑ i V isVid ∗ = 0 as a consequence<br />

<strong>of</strong> the CKM-unitarity. For hierarchical quark masses, the heaviest quark usually<br />

dominates (e.g. top-quark in b → sγ).<br />

1


• Standard CKM parametrization (PDG) in terms <strong>of</strong> angles and phase.<br />

• Wolfenstein parametrization, based on θ 13 ≪ θ 23 ≪ θ 12 ≪ 1<br />

λ ≡ sin θ 12 =<br />

¯ρ + i¯η ≡ − V udV ∗<br />

ub<br />

V cd V ∗<br />

cb<br />

|V us |<br />

√<br />

|Vud | 2 + |V us | 2 ,<br />

≃ sin θ 13e iδ<br />

Aλ 3 .<br />

Aλ2 ≡ sin θ 23 = λ |V cb|<br />

|V us | ,<br />

• Unitary triangle(s) in complex plane: Standard triangle from<br />

0 = V ∗<br />

ub V ud + V ∗<br />

cb V cd + V ∗<br />

tb V td<br />

V ∗<br />

cb V cd<br />

<strong>with</strong> tip <strong>of</strong> the triangle at ¯ρ + i¯η. Area <strong>of</strong> any CKM triangle given by Jarlskog<br />

determinant J/2.<br />

• Sides and angles <strong>of</strong> the triangle from absolute values and arguments <strong>of</strong> ratios <strong>of</strong><br />

CKM elements. Aim <strong>of</strong> flavour physics: Overconstrain the CKM triangle to precisely<br />

measure (¯ρ, ¯η) and test the CKM mechanism in the SM.<br />

– |V ij | from semi-leptonic decays (i ≠ t, no top bound states)<br />

– |V tb | ≃ 1 from t → bW +<br />

– |V td | and |V ts | from loop induced rare decays or meson-mixing involving b → d,<br />

b → s or s → d transitions.<br />

– CKM angles from interference effects, leading to CP asymmetries between<br />

particle and anti-particle decay.<br />

Requires precision experiments <strong>with</strong> kaons and B-mesons and theoretical control<br />

on strong interaction effects.<br />

<strong>Questions</strong>/Exercises:<br />

• Calculate<br />

V CKM = R 23 (θ 23 ) · diag(1, 1, e iδ ) · R 13 (θ 13 ) · diag(1, 1, e −iδ ) · R 12 (θ 12 )<br />

<strong>with</strong> R ij =rotations in the i − j plane.<br />

• Disentangle angles and phases in the Yukawa matrices and in the broken<br />

flavour symmetry generators, and find that there are 3 physical angles and<br />

1 physical CP phase left. Repeat the analysis <strong>for</strong> (2,4,n) generations.<br />

• What is the GIM mechanism<br />

2


Further Reading:<br />

• [1, page 13-16]; [2, page 3-11]; [3];<br />

3


2 The Weak Effective Hamiltonian<br />

Example: b → csū decay<br />

• Can be approximated by 4-quark interactions, similar to the Fermi model <strong>for</strong> muon<br />

decay (µ − → e −¯ν e ν µ )<br />

• In this approximation, the decay amplitude can be obtained from an effective<br />

operator in the interaction Hamiltonian:<br />

<strong>with</strong> two left-handed currents,<br />

O eff = G F<br />

√ J α<br />

(b→c) g αβ J (u→s)<br />

β<br />

(1)<br />

2<br />

J (b→c)<br />

α = V cb [¯cγ α (1 − γ 5 )b] ,<br />

J (u→s)<br />

β<br />

and the usual definition <strong>of</strong> the Fermi constant,<br />

= V ∗<br />

us [¯sγ β (1 − γ 5 )u] , (2)<br />

G F<br />

√<br />

2<br />

≡<br />

g2<br />

8m 2 W<br />

≃ 1.166 · 10 −5 GeV −2 .<br />

• We encounter one <strong>of</strong> the simplest examples <strong>for</strong> an Effective Theory, <strong>with</strong> the<br />

following general features:<br />

– operators in the effective interaction Hamiltonian have dimension > 4 (in the<br />

above case dim=6), and the corresponding coupling constants (here G F ) have<br />

negative mass dimension, and the theory – in general – is non-renormalizable,<br />

i.e. only valid up to some cut-<strong>of</strong>f scale (here m W );<br />

– effective operators are contain light (low-energetic) field modes, only<br />

(here light quarks);<br />

– effect <strong>of</strong> large energy/mass scales (here m W ) absorbed into coupling constants.<br />

4


<strong>Questions</strong>/Exercises:<br />

• Write down the Feynman diagram and the effective operator <strong>for</strong> the decay<br />

c → su ¯d in the Fermi model.<br />

2.1 Quantum Corrections to b → csū<br />

A typical QCD loop correction involves the exchange <strong>of</strong> an additional gluon between the<br />

two hadronic currents:<br />

• The loop momentum q µ can take any values between −∞ and +∞, and there<strong>for</strong>e<br />

a simple expansion in |q|/m W as in the Fermi model is not possible.<br />

• The idea <strong>of</strong> “Factorization” <strong>with</strong>in the context <strong>of</strong> Effective Field Theories is<br />

to systematically separate the cases <strong>of</strong><br />

– |q| ≥ m W (“short-distance” dynamics)<br />

– |q| ≪ m W (“long-distance” dynamics)<br />

An intuitive approach (which can be mathematically realized <strong>with</strong>in e.g. dimensional<br />

regularization) is to consider a Feynman integral in the “full theory” (which, itself can<br />

be an effective field theory) and divide it into its different momentum regions by<br />

expanding the integrand according to the scaling <strong>with</strong> the large mass/energy parameter:<br />

• In our case, the full theory integral <strong>for</strong> the above Feynman diagram would yield a<br />

generic function<br />

(<br />

I α s , m b<br />

, m )<br />

c<br />

, . . .<br />

m W m b<br />

where we assumed the function I to be dimensionless by factoring out an overall<br />

factor <strong>of</strong> G F , and there<strong>for</strong>e it can only depend on dimensionless quantities, i.e.<br />

the strong coupling constant and ratios <strong>of</strong> masses (and external momenta) <strong>of</strong> the<br />

particles being involved. 1<br />

1 We assume here, that potential UV divergences have been renormalized by appropriate counter<br />

terms and absorbed into running coupling/mass parameters.<br />

5


• In the infra-red (IR) region <strong>of</strong> “s<strong>of</strong>t” loop momenta |q µ | ≪ m W , we can still approximate<br />

the W -exchange by a point-like 4-fermion interaction, which correspond<br />

to taking the limit m W → ∞ in the integrand, and the above diagram would reduce<br />

to:<br />

– It represents a diagram <strong>with</strong>in the effective theory, i.e. a contribution to the<br />

vertex correction to the effective Hamiltonian.<br />

– Since one propagator has shrunk to a point, the diagram behaves differently<br />

<strong>for</strong> large momenta. In particular – in the general case – additional<br />

UV-divergences appear which are not compensated by the counter terms inherited<br />

from the full theory. One there<strong>for</strong>e has to introduce an additional<br />

UV-regularization (e.g. by an explicit cut-<strong>of</strong>f <strong>of</strong> order m W , or in terms <strong>of</strong><br />

dimensional regularisation).<br />

– The corresponding (dimensionally renormalized) integral will have the generic<br />

<strong>for</strong>m<br />

(<br />

I IR = I IR α s , ln µ , m )<br />

c<br />

, . . . ,<br />

m b m b<br />

i.e. it does not depend on the heavy mass m W , but instead it depends on an<br />

arbitrary UV-renormalization scale µ.<br />

• On the other hand, in the UV-region, we have<br />

|q| m W ≫ m b , m c , . . .<br />

and it can be obtained from the full theory integral by neglecting the masses (as<br />

well as external momenta) <strong>of</strong> the light degrees <strong>of</strong> freedom, m b , m c , . . . → 0.<br />

– The UV behaviour is as <strong>for</strong> the full theory (i.e. regularized by the same counter<br />

terms).<br />

– But now, one encounters new IR divergences, because the light masses have<br />

been set to zero.<br />

6


There<strong>for</strong>e, the UV-region <strong>of</strong> the integral (after renormalization) will have the<br />

generic <strong>for</strong>m<br />

(<br />

I UV = I UV α s , ln µ )<br />

m W<br />

depending now on an IR-regulator scale µ, but being manifestly independent <strong>of</strong> the<br />

low-energy (long-distance) parameters <strong>of</strong> the theory. It thus represents a genuine<br />

“short-distance” effect which can be interpreted as a quantum correction to socalled<br />

Wilson Coefficients multiplying the effective operators in the interaction<br />

Hamiltonian.<br />

• Notice that, due to the colour exchange, the operator basis generally has to be<br />

extended.<br />

To summarize, a generic 1-loop Feynman integral in the full theory is reproduced by the<br />

sum <strong>of</strong> a 1-loop matrix element <strong>of</strong> an effective operator <strong>with</strong>in the effective theory (the<br />

IR region), and a 1-loop short-distance coefficient (from the UV region) multiplying the<br />

tree-level matrix element <strong>of</strong> (another) operator in the effective theory:<br />

I<br />

(<br />

α s , m b<br />

m W<br />

, m c<br />

m b<br />

, . . .<br />

)<br />

= 〈O〉 1−loop<br />

(<br />

Effective Operators <strong>for</strong> b → csū<br />

α s , ln µ m b<br />

, m c<br />

m b<br />

, . . .<br />

)<br />

+ C ′ (<br />

α s , ln µ<br />

m W<br />

)<br />

〈O ′ 〉 tree (3)<br />

Use symmetry considerations (colour, flavour, Lorentz, . . . ) to constrain the operator<br />

basis<br />

• short-distance QCD corrections preserve chirality (Wilson coefficients calculated<br />

<strong>with</strong> m q → 0)<br />

⇒ only (V − A) ⊗ (V − A) structure involving left-handed quark fields<br />

• there are two independent possibilities to construct a colour singlet from two quark<br />

and two anti-quark fields. Two independent colour structures can be chosen as<br />

O 1 = (¯s ) i<br />

Lγ α u j j<br />

L<br />

(¯c<br />

L γα bL) i ,<br />

O 2 = (¯s ) )<br />

i<br />

Lγ α u i j<br />

L<br />

(¯c<br />

L γα b j L<br />

(4)<br />

• The decays b → csū are thus described by two current-current operators in the<br />

effective Hamiltonian and the corresponding Wilson coefficients<br />

7


H eff = − 4G F<br />

√<br />

2<br />

V cb V ∗<br />

∑<br />

C i (µ) O i + h.c. (5)<br />

us<br />

i=1,2<br />

Here, the Wilson coefficients C i (µ) contain all in<strong>for</strong>mation on short-distance physics<br />

above the scale µ, whereas the matrix elements <strong>of</strong> the operator reproduce the IR<br />

behaviour <strong>of</strong> the theory.<br />

<strong>Questions</strong>/Exercises:<br />

– How are the different momentum regions <strong>of</strong> Feynman integrals in the full<br />

SM theory reproduced <strong>with</strong>in the effective theory based on the weak effective<br />

Hamiltonian<br />

– Write down the two current-current operators <strong>for</strong> the decay c → su ¯d. Explain<br />

the appearance <strong>of</strong> the two colour structures.<br />

– Why does the electroweak Hamiltonian in the SM does not contain currentcurrent<br />

operators <strong>with</strong> right-handed quarks<br />

– Consider the toy integral<br />

∫ ∞<br />

0<br />

dk<br />

M − m<br />

(k + M)(k + m)<br />

∗ Calculate the integral exactly, and then determine the leading term <strong>for</strong><br />

m/M ≪ 1.<br />

∗ Calculate the UV-region <strong>of</strong> the integral by setting m → 0 in the integrand.<br />

Use a regulator (k/µ) δ in the integrand <strong>with</strong> δ → 0 + .<br />

∗ Calculate the IR-region <strong>of</strong> the integral by considering M → ∞ in the<br />

integrand, and use the same regulator (k/µ) δ <strong>with</strong> δ → 0 − .<br />

∗ Verify that the sum <strong>of</strong> the IR and UV region reproduce the full integral<br />

in the limit m/M ≪ 1.<br />

Further Reading:<br />

– [2, chapter 2.1]<br />

8


2.2 Wilson Coefficients in Perturbation Theory<br />

Following the above lines, the 1-loop result <strong>for</strong> the Wilson coefficients C 1,2 (µ) (renormalized<br />

in the MS scheme) reads<br />

{ } 0<br />

C i (µ) =<br />

1<br />

+ α s(µ)<br />

4π<br />

(ln µ2<br />

m 2 W<br />

+ 11 ) { } 3<br />

6 −1<br />

(6)<br />

As indicated, the result depends on the arbitrary renormalization scale µ. The advantage<br />

<strong>of</strong> the effective-theory construction precisely lies in the fact that the Wilson coefficients<br />

only depend on the physical short-distance scale m W . There<strong>for</strong>e, the natural procedure<br />

is to consider µ ∼ O(m W ), which in this context is called the “Matching Scale”:<br />

• the perturbative expansion <strong>of</strong> C i (µ ∼ m W ) does not contain coefficients enhanced<br />

by large logarithms<br />

• the Wilson coefficients can be reliably calculated in fixed-order perturbation<br />

theory, since α s (m W )/π ≪ 1.<br />

Anomalous Dimension<br />

In order to compare <strong>with</strong> experiment, the matrix elements <strong>of</strong> the effective-theory operators<br />

are needed at low-energy (hadronic) scales. For instance, in B-meson decays<br />

we typically need µ ∼ m b in order that no large logarithms ln µ m b<br />

appear. From the<br />

effective-theory construction, it is clear that only the combination<br />

∑<br />

C i (µ)〈O i 〉(µ)<br />

i<br />

is µ-independent in perturbation theory. The advantage <strong>of</strong> the effective theory is that it<br />

also allows to calculate the scale dependence <strong>of</strong> the operator matrix elements, respectively<br />

<strong>of</strong> the Wilson coefficients, <strong>with</strong>in the effective theory. The scale dependence is quantified<br />

by the Anomalous Dimension Matrix γ ij which itself has an expansion in α s (µ),<br />

∂<br />

∂ ln µ C i(µ) ≡ C j (µ) γ ji (µ) , (7)<br />

γ(µ) = γ(α s (µ)) = α s(µ)<br />

4π γ(1) + . . . (8)<br />

where γ (1,2,...) are numerical coefficient matrices that can be calculated from the renormalization<br />

constants <strong>for</strong> the 〈O i 〉. In our case<br />

9


γ (1) =<br />

( ) −2 6<br />

6 −2<br />

(9)<br />

<strong>with</strong> eigenvectors C ± = 1 √<br />

2<br />

(C 2 ± C 1 ) and eigenvalues γ (1)<br />

± = +4, −8. The differential<br />

equation (7) can then <strong>for</strong>mally be solved by separation <strong>of</strong> variables, leading to<br />

ln C ∫ ln µ<br />

±(µ)<br />

C ± (M) = d ln µ ′ γ ± (µ ′ ) (10)<br />

ln M<br />

Substituting the QCD β-function,<br />

d ln µ ′ = dα s<br />

2β , 2β = −2β 0<br />

4π α2 s + . . . (11)<br />

and keeping only the leading coefficients β 0 and γ (1)<br />

± , one ends up <strong>with</strong> the solution<br />

( ) (1) −γ αs<br />

±<br />

(µ)<br />

/2β 0<br />

C ± (µ) ≃ C ± (m W )<br />

α s (m W )<br />

(12)<br />

α<br />

Re-expanding the ratio s(µ)<br />

α s(m W<br />

, one realizes that the above <strong>for</strong>mula resums the socalled<br />

“leading logarithms”. Starting at the matching scale m W we can thus use<br />

)<br />

the ”renormalization-group running” <strong>of</strong> the Wilson coefficients in the effective theory<br />

to obtain “renormalization-group improved” predictions at the low-energy scale<br />

(as compared to the fixed-order result in (6)).<br />

Numerical Values (<strong>for</strong> our example):<br />

C 1 (m b ) C 2 (m b ) approximation<br />

−0.514 1.026 leading-log<br />

−0.303 1.008 next-to-leading-log<br />

Remark: Possible effect <strong>of</strong> New Physics contributions:<br />

• change value <strong>of</strong> Wilson coefficients at matching scale<br />

• change operator basis (right-handed, scalar, or tensor currents)<br />

10


<strong>Questions</strong>/Exercises:<br />

• Pro<strong>of</strong> the Fierz identity between colour matrices,<br />

t a ikt a jl = − 1<br />

3N C<br />

δ ik δ jl + 1 2 δ ilδ jk ,<br />

to reproduce the relative weights between the 1-loop corrections to C 1 and<br />

C 2 in (6).<br />

• Verify the eigenvectors and eigenvalues <strong>of</strong> the anomalous dimension matrix<br />

in (9).<br />

• Work out the details in the derivation <strong>of</strong> (12).<br />

• Re-expand (12) in fixed-order perturbation theory, using the one-loop approximation<br />

<strong>for</strong> the ratio α s (µ)/α s (m W ). Compare <strong>with</strong> (6).<br />

Further Reading:<br />

• [2, chapter 2.2]; [4, chapter 18.2]; [1, page 14-41]<br />

11


2.3 Penguins and Boxes<br />

Example: b → sq¯q decays<br />

We are now considering hadronic b → s transitions <strong>with</strong> an arbitrary (flavour-singlet) q¯q<br />

pair in the final state (analogously b → dq¯q or s → dq¯q).<br />

• As be<strong>for</strong>e, we obtain effective current-current operators, but now there are two<br />

different flavour structures contributing<br />

V ub V ∗<br />

us [ū L γ µ b L ] [¯s L γ µ u L ] ≡ λ u O (u)<br />

2 ,<br />

V cb V ∗<br />

cs [¯c L γ µ b L ] [¯s L γ µ c L ] ≡ λ c O (c)<br />

2 . (13)<br />

By gluon exchange, we also generate the corresponding current-current operators<br />

<strong>with</strong> differing colour structure<br />

O (u)<br />

1 , O (c)<br />

1 .<br />

• The same final state can also be obtained by so-called penguin diagrams, where<br />

the q¯q-pair is produced from gluon radiation <strong>of</strong>f a loop diagram:<br />

– The produced q¯q pair now can refer to left- or right-handed fields. This gives<br />

rise to four new types <strong>of</strong> “strong penguin operators” O 3−6 <strong>of</strong> the <strong>for</strong>m<br />

∑<br />

[¯s L γ µ b L ] [¯q L γ µ q L ] ,<br />

q=u,d,s,c,b<br />

∑<br />

q=u,d,s,c,b<br />

[¯s L γ µ b L ] [¯q R γ µ q R ] , (14)<br />

each coming in two colour variants. The Wilson coefficients <strong>for</strong> these operators<br />

are suppressed by the strong coupling constant and a loop factor.<br />

– Actually, we also have to allow <strong>for</strong> operators describing the partonic b → sg<br />

process, because it will generate the same flavour quantum numbers. Using<br />

QCD gauge invariance and exploiting the equations <strong>of</strong> motion <strong>for</strong> the quark<br />

12


fields, one realizes that the only relevant operators have to have the <strong>for</strong>m <strong>of</strong> a<br />

“chromomagnetic penguin” coupling left- and right-handed quark fields,<br />

g s<br />

O g 8 ≡<br />

16π m 2 b ¯s L σ µν G µν b R ,<br />

O g 8 ′ ≡<br />

g s<br />

16π m 2 s ¯s R σ µν G µν b L . (15)<br />

Notice that the operators include an additional quark-mass factor, and there<strong>for</strong>e<br />

<strong>for</strong>mally are counted as dim-6. The second operator is mostly negligible<br />

since m s ≪ m b .<br />

– The CKM factor <strong>for</strong> the penguin operators can be determined by exploiting<br />

the GIM mechanism: The result <strong>of</strong> any penguin diagram gets contributions<br />

from up, charm and top-quarks in the loop. As discussed be<strong>for</strong>e, in the<br />

matching calculation we can set m u = m c = 0, and there<strong>for</strong>e, generically, one<br />

ends up <strong>with</strong> functions<br />

∑<br />

V ib Vis ∗ f(m 2 i /m 2 W ) = λ u f(0) + λ c f(0) + λ t f(m 2 t /m 2 W )<br />

i=u,c,t<br />

= λ t<br />

(<br />

f(m<br />

2<br />

t /m 2 W ) − f(0) ) , (16)<br />

where in the last step, we have used the CKM unitarity.<br />

The weak effective Hamiltonian <strong>for</strong> the b → sq¯q transitions (including strong-interaction<br />

effects) there<strong>for</strong>e reads<br />

H eff = − 4G F<br />

√<br />

2<br />

{ ∑<br />

+ h.c.<br />

i=1,2<br />

(<br />

C i (µ) λ u O (u)<br />

i<br />

(<br />

)<br />

6∑<br />

)}<br />

+ λ c O (c)<br />

i + λ t C i (µ)O i + C8(µ)O g g 8<br />

i=3<br />

(17)<br />

Electroweak penguin operators<br />

Electroweak corrections to the weak effective Hamiltonian can become important <strong>for</strong><br />

• precision observables, in particular isospin asymmetries,<br />

• weak decays <strong>with</strong> leptons or photons in the final state.<br />

To this end, penguin and box diagram <strong>with</strong> additional photon or Z-exchange have to be<br />

included. This yields:<br />

13


• electroweak corrections to the matching and running <strong>of</strong> the Wilson coefficients in<br />

(17),<br />

• new isospin-violating effects, to be accounted <strong>for</strong> by electroweak penguin operators,<br />

O 7 =<br />

∑<br />

e q [¯s L γ µ b L ] [¯q L γ µ q L ] , O 8 = ∑<br />

]<br />

i<br />

e q<br />

[¯s L γ µ b j j<br />

L<br />

[¯q<br />

L γ µqL] i ,<br />

O 9 =<br />

q=u,d,s,c,b<br />

∑<br />

q=u,d,s,c,b<br />

q=u,d,s,c,b<br />

e q [¯s L γ µ b L ] [¯q R γ µ q R ] , O 10 = ∑<br />

q=u,d,s,c,b<br />

]<br />

i<br />

e q<br />

[¯s L γ µ b j j<br />

L<br />

[¯q<br />

R γ ]<br />

µqR<br />

i ,<br />

(18)<br />

which differ from the strong penguin operators by an additional factor e q refering<br />

to the fractional quark charge, and by the fact that their Wilson coefficients are<br />

suppressed by an additional factor α em /α s .<br />

• a new electromagnetic operator, describing b → sγ,<br />

O γ 7 ≡<br />

where we neglected O γ 7 ′ <strong>for</strong> m s → 0.<br />

Semi-leptonic operators<br />

e<br />

16π 2 m b ¯s L σ µν b R F µν , (19)<br />

For decays into charged leptons, b → sl + l − , another set <strong>of</strong> semi-leptonic operators has<br />

to be added,<br />

O 9,V = [¯s L γ µ b L ] ¯lγ µ l ,<br />

O 10,A = [¯s L γ µ b L ] ¯lγ µ γ 5 l , (20)<br />

For decays into neutrinos, b → sl + l − , only the corresponding 4-fermion operator has to<br />

be taken into account, as the neutrinos can only arise from short-distance processes, i.e.<br />

penguin diagrams <strong>with</strong> Z-boson radiation or box-diagrams <strong>with</strong> two W -boson exchange.<br />

Operators <strong>for</strong> neutral meson mixing<br />

In neutral meson mixing, we consider flavour transitions <strong>of</strong> the <strong>for</strong>m<br />

b ¯d ↔ d¯b (|∆B| = 2, ∆S = 0) ,<br />

b¯s ↔ s¯b (|∆B| = |∆S| = 2) ,<br />

s ¯d ↔ d¯s (|∆S| = 2, ∆B = 0) , (21)<br />

This requires box diagrams <strong>with</strong> two W -exchanges again involving only left-handed<br />

quarks,<br />

14


H |∆B|=2,∆S=0<br />

eff<br />

= G2 F<br />

4π 2 m2 W (VtbV ∗<br />

td ) 2 C(m t /m W , µ) [¯bL ] ]<br />

γ µ d L<br />

[¯bL γ µ d L + h.c. (22)<br />

Because the field content <strong>of</strong> the two currents in brackets is identical:<br />

• there is only one independent colour structure,<br />

• the CKM factor is the simple square (not the absolute square) <strong>of</strong> the CKM elements.<br />

<strong>Questions</strong>/Exercises:<br />

• What would be the CKM factor <strong>of</strong> the penguin operators in s → dq¯q and in<br />

c → uq¯q decays<br />

• What are the crucial differences between electroweak and strong penguin<br />

operators<br />

• Draw an example <strong>for</strong> an effective-theory diagram contributing to b → sγ and<br />

involving one <strong>of</strong> the strong penguin operators O 3−6 (or the operator O g 8).<br />

• Draw examples or box and penguin diagrams contributing to b → sl + l − and<br />

b → sν¯ν.<br />

Further Reading:<br />

• [1, part two (p. 43ff)]<br />

15


3 Flavour Transitions (on the hadron level)<br />

The analysis <strong>of</strong> quark-flavour transitions is complicated by the fact that quark are confined<br />

in hadronic bound states. As an example, we will study the different flavour<br />

topologies that can contribute to 2-body non-leptonic decays.<br />

3.1 Example: ¯B → DK decays<br />

As an example, we will study the quark-level transition b → csū (see previous section)<br />

<strong>with</strong>in the non-leptonic decay modes B → DK. Depending on the spectator quark in<br />

the B-meson, we may consider the following hadronic decays: 2<br />

(a) ¯B0 d → D + K − <strong>with</strong> (b ¯d) → (c ¯d)(sū) , (23)<br />

(b)<br />

(c)<br />

¯B0 d → D 0 ¯K0<br />

B − → D 0 ¯K−<br />

<strong>with</strong> (b ¯d) → (cū)(s ¯d) , (24)<br />

<strong>with</strong> (bū) → (cū)(sū) , (25)<br />

(c) ¯B0 s → D + s ¯K − <strong>with</strong> (b¯s) → (c¯s)(sū) . (26)<br />

Let us consider the “flavour topologies” <strong>for</strong> each decay mode in turn<br />

(a) ¯B 0 d → D+ K − : Here the ¯d-quark simply acts as a spectator quark in the B → D<br />

transition, while the kaon flavour stems from the quarks emitted from the weak<br />

decay.<br />

This case is sometimes refered to as a “Tree Topology” (class-I).<br />

(b) ¯B 0 d → D0 ¯K0 : Here the ¯d-quark acts as a spectator quark in the B → K transition,<br />

while the D-meson flavour stems from the quarks emitted from the weak decay.<br />

2 We may have also considered B → DK decays induced by b → us¯c decays.<br />

16


This case is sometimes refered to as a “Tree Topology” (class-II).<br />

(c) B − → D 0 K − : Here the ū-quark can act as a spectator quark in either the B → D<br />

or the B → K transition.<br />

There<strong>for</strong>e, this decay is induced by “Tree Topologies” <strong>of</strong> class-I and class-II.<br />

(d) ¯B s<br />

0 → D s + K − : This can be induced by a class-I tree topology, but also by an<br />

“Annihilation Topology” (class-III), where both quarks from the initial-state<br />

meson participate in the weak decay.<br />

Each decay topology gives rise to qualitatively different hadronic bound state effects:<br />

17


Class-I Topology<br />

In the “naive“ picture, we would neglect interactions between the emitted quarks <strong>for</strong>ming<br />

the kaon and the quarks involved in the B → D transition. Starting from the hadronic<br />

matrix element <strong>of</strong> the weak effective Hamiltonian,<br />

〈D + K − |H eff | ¯B 0 d〉 ∼ 4G F<br />

√<br />

2<br />

V cb V ∗<br />

cs C i (µ) 〈D + K − |O i | ¯B 0 d〉<br />

we obtain the so-called ”naive factorization approximation“,<br />

〈D + K − |O 1,2 | ¯B 0 d〉 ∼ 〈D + |J (b→c) |B 0 d〉<br />

} {{ } × 〈K− |J (u→s) |0〉<br />

} {{ } .<br />

B → D <strong>for</strong>m factor<br />

kaon decay constant<br />

which only involves to (relatively) simple hadronic objects:<br />

• The B → D transition <strong>for</strong>m factor F B→D (q 2 = m 2 K ).<br />

• The kaon decay constant f K .<br />

In order to determine the individual contributions from the two operators O 1,2 , we have<br />

to project onto colour singlet currents:<br />

• The operator O 2 is already given as the product <strong>of</strong> two colour-singlet currents<br />

J (b→c) and J (u→s) .<br />

• For the operator O 1 we make use <strong>of</strong> the Fierz identity (see above)<br />

(¯c i Lγ µ b j L )(¯sj L γµ u i L) = 1 3 J (b→c) J (u→s) + octet × octet<br />

Consequently, <strong>for</strong> the class-I topologies we obtain<br />

∑<br />

C i (µ) 〈D + K − |O i | ¯B d〉 0 ∼<br />

i=1,2<br />

(<br />

C 2 (µ) + 1 )<br />

3 C 1(µ) · F B→D (m 2 K) · f K<br />

However, we immediately see that the above <strong>for</strong>mula cannot be correct, as the transition<br />

<strong>for</strong>m factors and the decay constants do not compensate the scale-dependence <strong>of</strong> the<br />

Wilson coefficients.<br />

18


The correct scale-dependence <strong>of</strong> the hadronic matrix elements actually comes from the<br />

”non-factorizing” gluon cross-talk between the quarks in the kaon and the quarks<br />

in the B → D transitions. For the naive factorization hypothesis to yield a good approximation<br />

to the real case, we would have to show that the non-factorizing effects<br />

are<br />

• suppressed by α s (m b ) in perturbation theory<br />

(i.e. dominated by short-distances <strong>of</strong> order δx ∼ 1/m b ),<br />

• or suppressed by Λ QCD /m b (sub-leading non-perturbative effects).<br />

An intuitive explanation involves the so-called colour-transparency argument:<br />

• The light kaon (in the B-meson rest frame) is rather energetic (E ∼ 2 GeV<br />

∼ O(m b /2)). There<strong>for</strong>e the non-factorizing gluons effectively “see” a rather small<br />

colour dipole <strong>of</strong> a size δx ∼ 1/E. The coupling constant <strong>of</strong> the gluons to this dipole<br />

is thus given by α s (E) ≪ 1.<br />

• Couplings to higher multipoles in the hadronic transitions will be suppressed by<br />

Λ QCD /E.<br />

A diagrammatic pro<strong>of</strong> that the colour-transparency argument <strong>for</strong> a certain class <strong>of</strong> nonleptonic<br />

B decays is correct has been given in [5]. Moreover, the perturbative corrections<br />

can be systematically calculated, giving rise to the so-called “QCD(-improved) Factorization<br />

Approach”. The <strong>for</strong>mula <strong>for</strong> the hadronic matrix elements is modified as<br />

∑<br />

C i (µ) 〈D + K − |O i | ¯B d〉<br />

0<br />

i=1,2<br />

∼ F B→D · f K ·<br />

∫ 1<br />

0<br />

(<br />

du C 2 (µ) + 1 3 C 1(µ) + α )<br />

s(µ)<br />

4π t(u, µ) φ K(u, µ) + . . .<br />

• The function t(u, µ) can be calculated <strong>with</strong>in QCD perturbation theory and depends<br />

on the energy/momentum fractions u and ū = 1 − u <strong>for</strong> a quark–antiquark<br />

state in the kaon.<br />

• There<strong>for</strong>e, the final result <strong>for</strong> the transition amplitude is obtained by a convolution<br />

integral, involving a new hadronic input function, φ K (u, µ), characterizing the<br />

probability amplitude <strong>for</strong> the energy/momentum fraction u in the kaon, <strong>with</strong><br />

∫ 1<br />

0<br />

du φ K (u, µ) = 1 .<br />

• The calculation <strong>of</strong> Λ QCD /m b corrections <strong>with</strong>in the QCDF approach, on the other<br />

hand, is difficult and yields an irreducible source <strong>of</strong> theoretical uncertainties.<br />

19


Class-II Topology<br />

• Naively, we would expect a similar factorization into a B → K <strong>for</strong>m factor and the<br />

D-meson decay constant.<br />

• Concerning the colour structure, the role <strong>of</strong> O 1 and O 2 is now reversed, and the<br />

amplitude would be proportional to<br />

C 1 + 1 3 C 2<br />

which happens to be a rather small number (see Table in the previous chapter).<br />

The class-II decay are there<strong>for</strong>e also refered to as “colour-suppressed”.<br />

• Finally, the colour transparency argument does not apply <strong>for</strong> the class-II topologies,<br />

because the ratio <strong>of</strong> E/m D in the decay is moderate and there<strong>for</strong>e the emitted D-<br />

meson corresponds to a rather large colour dipole.<br />

Fortunately, in this case, the transition <strong>of</strong> a B-meson at rest into an energetic kaon<br />

itself is a sub-leading effect, since the spectator quark has to be accelerated from<br />

essentially at rest to a constituent <strong>of</strong> a fast kaon. There<strong>for</strong>e the class-II topologies<br />

<strong>for</strong> B → DK are suppressed by 1/E K <strong>with</strong> respect to the class-I topology.<br />

Annihilation Topology<br />

Again, to produce the additional fast (back-to-back) s¯s pair from the vacuum, one receives<br />

a suppression <strong>of</strong> the decay amplitude <strong>with</strong> 1/E ∼ 1/m b , and the annihilation<br />

topology is sub-leading compared to class-I decays.<br />

3.2 Example: B → πK decays<br />

We briefly comment on the modifications – compared to the previous subsection – that<br />

arise when both final-state mesons are light:<br />

• Both final-state mesons are energetic, E ∼ m b /2, and there<strong>for</strong>e class-I and class-II<br />

decays can be treated on the same footing, and the colour transparency argument<br />

applies in both cases.<br />

• The perturbative calculation <strong>of</strong> the non-factorizable gluon exchange diagrams, now<br />

depends on the momentum fraction <strong>of</strong> all the various quarks involved in the decay.<br />

As a consequence, one has to specify the momentum distributions φ π (u) and φ K (u)<br />

as well as the distribution φ B (ω), where ω is the energy <strong>of</strong> the spectator quark in<br />

the B-meson.<br />

• As the fundamental quark transitions involve b → sq¯q, one has to extend the<br />

analysis to include hadronic matrix elements <strong>of</strong> the penguin operators (O 3−6,8 g ,...).<br />

The give rise to flavour topologies like the following,<br />

20


We thus encounter many contributions, involving many hadronic input parameters, and<br />

different kind <strong>of</strong> suppression factors (CKM elements, Wilson coefficients, colour factors,<br />

α s , 1/m b , . . . ). The analysis <strong>of</strong> systematic theoretical uncertainties is thus very complicated<br />

and difficult. As a general rule, we may conclude:<br />

• The colour-allowed amplitudes (class-I) are generally well-predicted and theoretical<br />

estimates are in line <strong>with</strong> experimental measurements.<br />

• The colour-suppressed amplitudes (class-II) have large theoretical errors because<br />

<strong>of</strong> numerical cancellations.<br />

• The annihilation topologies cannot be predicted reliably (but are expected to be<br />

sub-leading).<br />

<strong>Questions</strong>/Exercises:<br />

• Draw the possible flavour topologies <strong>for</strong> B → DK decays, but now <strong>with</strong> the<br />

decay currents <strong>for</strong> b → us¯c.<br />

• Which decay operators and topologies appear in B → Dπ decays<br />

• What is naive factorization, and how can it be improved<br />

• Why are the class-II and class-III topologies in B → DK decays suppressed<br />

<strong>with</strong> respect to class-I <br />

Further Reading:<br />

• [6]<br />

21


4 Mixing and CP Violation<br />

The ∆S, B = 2 operators specified in chapter 2, induce mixing between neutral mesons<br />

and their anti-mesons.<br />

4.1 K 0 - ¯K 0 Mixing<br />

• The neutral kaons can be specified by their flavour content<br />

K 0 ∼ (¯sd) ,<br />

¯K0 ∼ (s ¯d)<br />

These flavour eigenstates are not the physical (mass) eigenstates.<br />

• Chosing a phase convention <strong>for</strong> CP-trans<strong>for</strong>mations,<br />

we can construct CP eigenstates,<br />

CP |K 0 〉 = −| ¯K 0 〉 , CP | ¯K 0 〉 = −|K 0 〉 ,<br />

|K 1 〉 ≡ 1 √<br />

2<br />

(<br />

|K 0 〉 − | ¯K 0 〉 ) , CP |K 1 〉 = +|K 1 〉 ,<br />

|K 2 〉 ≡ 1 √<br />

2<br />

(<br />

|K 0 〉 + | ¯K 0 〉 ) , CP |K 2 〉 = −|K 2 〉 . (27)<br />

They would be the physical eigenstates, if CP were conserved.<br />

• The neutral kaons can decay into pions, <strong>with</strong> |ππ〉 being CP-even, and |πππ〉 being<br />

CP-odd. If CP were a good symmetry, we expect<br />

K 1 → 2π<br />

large phase space → large width → short life-time<br />

K 2 → 3π small phase space → smaller width → larger life-time (28)<br />

• In reality, we observe small corrections from CP-violation in weak interactions, i.e.<br />

small admixtures <strong>of</strong> the opposite CP states in physical mass eigenstates and decay<br />

amplitudes. We classify the physical kaon eigenstates by their lifetime, and define<br />

K S = K 1 + ¯ɛK<br />

√ 2<br />

= 1 (1 + ¯ɛ)K 0 − (1 − ¯ɛ)<br />

√ ¯K 0<br />

√ ,<br />

1 + |¯ɛ|<br />

2 2 1 + |¯ɛ|<br />

2<br />

K L = K 2 + ¯ɛK<br />

√ 1<br />

= 1 (1 + ¯ɛ)K 0 + (1 − ¯ɛ)<br />

√ ¯K 0<br />

√ , (29)<br />

1 + |¯ɛ|<br />

2 2 1 + |¯ɛ|<br />

2<br />

where ¯ɛ is a (convention-dependent) small complex parameter.<br />

22


• The experimentally measured mass difference amounts to<br />

∆M K = M(K L ) − M(K S ) ∼ 3.5 · 10 −15 GeV<br />

and the width difference in the kaon system satisfies to a good approximation<br />

Theoretical Description<br />

∆Γ K ≈ −2∆M K .<br />

The masses and decay widths <strong>of</strong> the neutral kaons describe the time-evolution <strong>of</strong> the<br />

K 0 − ¯K 0 system in quantum mechanics,<br />

i dψ<br />

( ) |K<br />

dt = Ĥ ψ(t) , ψ(t) ∼ 0 (t)〉<br />

| ¯K 0 (30)<br />

(t)〉<br />

where the Hamilton operator contains a hermitian and an anti-hermitian part, related<br />

to virtual and real intermediate states, respectively,<br />

Ĥ = ˆM − i ˆΓ<br />

2<br />

(31)<br />

Here ˆM and ˆΓ are hermitian matrices,<br />

M 12 = M12, ∗ Γ 21 = Γ ∗ 12 .<br />

The CPT theorem further relates particle and anti-particle properties, such that<br />

M 11 = M 22 ≡ M , Γ 11 = Γ 22 ≡ Γ ,<br />

and the Hamiltonian can be written as<br />

( )<br />

M − iΓ/2 M12 − iΓ<br />

Ĥ =<br />

12 /2<br />

M12 ∗ − iΓ ∗ . (32)<br />

12/2 M − iΓ/2<br />

The eigenvalues <strong>of</strong> the matrix yield the mass and width differences:<br />

∆M = 2 Re Q , ∆Γ = −4 Im Q ,<br />

<strong>with</strong> Q = √ (M 12 − iΓ 12 /2)(M ∗ 12 − iΓ ∗ 12/2) . (33)<br />

23


The eigenvectors determine ¯ɛ,<br />

√<br />

1 − ¯ɛ<br />

1 + ¯ɛ = M12 ∗ − iΓ ∗ 12/2<br />

M 12 − iΓ 12 /2<br />

≡ r exp(iκ) . (34)<br />

The deviation <strong>of</strong> the parameter r from unity (or ¯ɛ from zero) is a measure <strong>for</strong> CP violation<br />

in kaon mixing,<br />

r = 1 +<br />

2|Γ 12 | 2<br />

4|M 12 | 2 + |Γ 12 | 2 Im M 12<br />

Γ 12<br />

≈ 1 + Im M 12<br />

Γ 12<br />

. (35)<br />

CP asymmetry in semi-leptonic kaon decays<br />

We consider the CP asymmetry <strong>for</strong> semi-leptonic K L decays into π ± ,<br />

a SL (K L ) = Γ(K L → π − e + ν e ) − Γ(K L → π + e −¯ν e )<br />

Γ(K L → π − e + ν e ) − Γ(K L → π + e −¯ν e )<br />

(36)<br />

Realizing that K 0 ( ¯K 0 ) can only decay into π − (π + ), we obtain<br />

Γ(K L → π − e + ν e ) ∝ ∣ ∣ (1 + ¯ɛ) A(K 0 → π − e + ν e ) ∣ ∣ 2 ,<br />

Γ(K L → π + e −¯ν e ) ∝ ∣ ∣ (1 − ¯ɛ) A( ¯K0 → π + e −¯ν e ) ∣ ∣ 2 (37)<br />

Since the strong interactions are CP-conserving (and there are no interference effects),<br />

the absolute values <strong>for</strong> the amplitudes cancel in the ratio, and one obtains<br />

a SL (K L ) = |1 + ¯ɛ|2 − |1 − ¯ɛ| 2 2Re ¯ɛ<br />

=<br />

|1 + ¯ɛ| 2 + |1 − ¯ɛ|<br />

2<br />

1 + |¯ɛ| = 1 − r2<br />

2 1 + r 2<br />

≈ 2Re ¯ɛ ≈ Im Γ 12<br />

M 12<br />

(38)<br />

24


This kind <strong>of</strong> CP asymmetry only depends on the K 0 - ¯K 0 mixing parameters and is thus<br />

dubbed<br />

“CP-Violation in Mixing”<br />

Experimentally one finds Re ¯ɛ ∼ 10 −3 .<br />

Direct and indirect CP violation in kaon decays<br />

CP violation can also manifest itself in the “wrong” decays K L → ππ (or similarly <strong>for</strong><br />

K S → 3π).<br />

K L ∼ K 2 + ¯ɛK 1<br />

↓<br />

direct CP-viol.<br />

↓<br />

ππ<br />

↓<br />

indirect CP-viol.<br />

↓<br />

ππ<br />

<strong>with</strong> “indirect CP violation“ from the small K 1 admixture in K L , and ”direct CP<br />

violation” through the CP-violating effects in the weak decay K 2 → ππ itself.<br />

• As a measure <strong>for</strong> the indirect CP violation, one defines the parameter<br />

ɛ ≡ A(K L → (ππ) I=0 )<br />

A(K S → (ππ) I=0 )<br />

(39)<br />

which compares the amplitudes <strong>for</strong> decays into the isospin-zero (I = 0) projection<br />

<strong>of</strong> the two pions in the final state (the intereference effects <strong>with</strong> the I = 2 projection<br />

will be relevant <strong>for</strong> the direct CP violation discussed below). The parameter ɛ can<br />

be expressed as<br />

ɛ = ¯ɛ + iξ ≈ exp(iπ/4) √<br />

2 ∆MK<br />

(Im M 12 + 2ξ Re M 12 ) (40)<br />

25


where the parameter ξ = Im A 0 /Re A 0 measures the weak phase <strong>of</strong> the I = 0 decay<br />

amplitude A 0 as defined below. Notice that ɛ is phase-convention independent, and<br />

Re ɛ = Re ¯ɛ.<br />

The decay K → (ππ) I=0 proceeds dominantly via s → udū and thus receives a CKM<br />

factor V us Vud ∗ which is real in the standard parametrization. There<strong>for</strong>e, ξ ≈ 0, and the<br />

parameter ɛ is basically sensitive to Im M 12 .<br />

• For direct CP-violation, we need the interference <strong>of</strong> two independent decay amplitudes.<br />

In K → ππ decays these are given by the two alternative isospin projections,<br />

3 √<br />

3<br />

A(K + → π + π 0 ) =<br />

2 A 2 e iδ 2<br />

,<br />

A(K 0 → π + π − ) =<br />

√<br />

2<br />

3 A 0 e iδ 0<br />

+<br />

√<br />

1<br />

3 A 2 e iδ 2<br />

,<br />

A(K 0 → π 0 π 0 ) =<br />

√<br />

2<br />

3 A 0 e iδ 0<br />

− 2<br />

√<br />

1<br />

3 A 2 e iδ 2<br />

. (41)<br />

The pre-factors are obtained from the SU(2) I Clebsch-Gordan coefficients, and<br />

δ 0,2 denote the so-called “strong phases” (which arise as a consequence <strong>of</strong> nonfactorizable<br />

QCD rescattering effects between the final-state quarks). The weak<br />

phases are still contained in the two independent amplitudes A 0 and A 2 , such that<br />

the CP-conjugate decays will proceed via A ∗ 0,2 (while δ 0,2 will not change sign).<br />

As a measure <strong>for</strong> indirect CP violation one then considers a particular combination<br />

<strong>of</strong> decay amplitudes where the contribution from the indirect CP violation drops<br />

out,<br />

ɛ ′ ≡ √ 1 ( A(KL → (ππ) I=2 )<br />

2 A(K S → (ππ) I=0 ) − A(K S → (ππ) I=2 )<br />

A(K S → (ππ) I=0 )<br />

)<br />

A(K L → (ππ) I=0 )<br />

A(K S → (ππ) I=0 )<br />

= 1 √<br />

2<br />

Im A 2<br />

A 0<br />

exp(iΦ ɛ ′) , <strong>with</strong> Φ ɛ ′ = π 2 + δ 2 − δ 0 (42)<br />

– Im A 2<br />

A 0<br />

comes from the different CKM factors <strong>of</strong> the operators in H eff contributing<br />

to the I = 0 or I = 2 final state.<br />

– The phase difference δ 2 − δ 0 comes from strong-interaction effects, and can be<br />

independently measured in ππ scattering, leading to Φ ɛ ′ ≈ π/4.<br />

3 The combination <strong>of</strong> two isospin-1 pion triplets yields 3 × 3 = 1 + 3 + 5, refering to I = 0, 1, 2.<br />

However, the I = 1 projection does not contribute because <strong>of</strong> Bose symmetry.<br />

26


• From the experimental decay widths, one can then extract the value <strong>of</strong> Re (ɛ ′ /ɛ)<br />

from the ratios<br />

leading to<br />

η 00 = A(K L → π 0 π 0 )<br />

A(K S → π 0 π 0 ) , η +− = A(K L → π + π − )<br />

A(K S → π + π − )<br />

( )<br />

∣ ɛ<br />

′<br />

Re ≈ 1 −<br />

η 00 ∣∣∣<br />

2<br />

ɛ ∣η +−<br />

(43)<br />

ɛ from H eff<br />

The value <strong>of</strong> ɛ can be related to the matrix element <strong>of</strong> the weak effective Hamiltonian<br />

<strong>for</strong> kaon mixing,<br />

2m K M ∗ 12 = 〈 ¯K 0 |H eff (∆S = 2)|K 0 〉 , (44)<br />

where the factor 2m K comes from the normalization <strong>of</strong> the hadronic states.<br />

• In the case <strong>of</strong> K 0 - ¯K 0 mixing, we have to take into account the specific CKM<br />

hierarchy <strong>for</strong> s → d transitions:<br />

– The top contribution yields a CKM factor (V td V ∗<br />

ts) 2 ∼ λ 10 .<br />

– Although the charm contribution is suppressed by factors <strong>of</strong> m 2 c/m 2 t , this is<br />

compensated by a much larger CKM factor, (V cd V ∗<br />

cs) 2 ∼ λ 2 .<br />

As a consequence, also the box diagrams <strong>with</strong> 2 charm-quark propagators or <strong>with</strong><br />

one charm-quark and one top-quar propagator will be important.<br />

This effect also contributes to the relative phase between M 12 and Γ 12 .<br />

• The hadronic matrix element M 12 can then be expressed in terms <strong>of</strong> function<br />

S 0 arising from the calculation <strong>of</strong> the box diagrams <strong>with</strong> different masses x c,t =<br />

m 2 c,t/MW 2 Ṁ 12 = G2 F<br />

12π 2 F 2 K ˆB K m K M 2 W<br />

× [ (λ ∗ c) 2 η 1 S 0 (x c ) + (λ ∗ t ) 2 η 2 S 0 (x t ) + 2λ ∗ cλ ∗ t η 3 S 0 (x c , x t ) ] (45)<br />

27


Figure 1: Determination <strong>of</strong> ¯ρ and ¯η from ɛ K (green band), ∆m d,s (yellow/orange) and<br />

sin 2β (blue).<br />

Here η i are short-distance QCD corrections, and F 2 K ˆB K parametrizes the relevant<br />

hadronic matrix element<br />

〈 ¯K 0 |(¯sd) (V −A) (¯sd) (V −A) |K 0 〉 ≡ 8 3 m K F 2 K ˆB K (1 + QCD corr.) (46)<br />

where F K is the kaon decay constant, and the normalization factors are chosen in<br />

such a way that in the naive factorization approximation, 〈 ¯K 0 |(¯sd) (V −A) (¯sd) (V −A) |K 0 〉 ≈<br />

〈 ¯K 0 |J|0〉〈0|J|K 0 〉, the quantity ˆB K is unity,<br />

ˆB K = 1 + non-factorizable corrections (47)<br />

• The actual numerical estimate <strong>for</strong> ˆB K ≈ 80% is obtained from non-perturbative<br />

methods. The theoretical <strong>for</strong>mula <strong>for</strong> M 12 can then be turned into a theoretical<br />

prediction <strong>for</strong> ɛ. Comparison <strong>with</strong> the experimental results provides a constraint on<br />

the CP-phase δ entering the CKM element V td (in the standard parametrization).<br />

The constraint takes the <strong>for</strong>m <strong>of</strong> a hyperbola in the ¯ρ-¯η plane <strong>for</strong> the standard<br />

CKM triangle <strong>with</strong> an error band dominated by the uncertainties on ˆB K , see the<br />

green band in Fig. 1.<br />

28


• Finally, the experimental measurement <strong>for</strong> direct CP violation in kaon decays yields<br />

( ) ɛ<br />

′<br />

Re ∼ 2 · 10 −3 ,<br />

ɛ<br />

i.e. the effect is further suppressed by 3 orders <strong>of</strong> magnitude compared to the<br />

indirect CP violation. A theoretical prediction requires the quantitive knowledge<br />

<strong>of</strong> the hadronic matrix elements responsible <strong>for</strong> the I = 0 and I = 2 amplitudes,<br />

and takes the generic <strong>for</strong>m<br />

ɛ ′<br />

ɛ ∼ Im λ t [〈· · · 〉 I=0 − 〈· · · 〉 I=2 ]<br />

Due to the cancellation effects occuring in the brackets, the theoretical predictions<br />

are thus rather uncertain.<br />

<strong>Summary</strong> <strong>of</strong> the different types <strong>of</strong> CP-violation observed in the neutral kaon system:<br />

• Re ɛ = Re ¯ɛ: CP violation in mixing<br />

• Im ɛ: CP-violation in the interference between mixing ang decay<br />

• Re ɛ ′ : direct CP violation in decay, where<br />

Re ɛ ′ = − 1 √<br />

2<br />

∣ ∣∣∣ A 2<br />

A 0<br />

∣ ∣∣∣<br />

sin(φ 2 − φ 0 ) sin(δ 2 − δ 0 ) (48)<br />

i.e. direct CP violation requires a non-vanishing difference φ 2 − φ 0 ≠ 0 <strong>of</strong> weak<br />

phases and a non-vanishing difference δ 2 − δ 0 ≠ 0 <strong>of</strong> strong phases in the I = 0<br />

and I = 2 amplitudes.<br />

<strong>Questions</strong>/Exercises:<br />

• Verify the eigenvalues and eigenvectors <strong>for</strong> the matrix Ĥ in (32).<br />

• What is the difference between direct and indirect CP violation in kaon<br />

decays<br />

• Why do we have to take into account the box diagrams <strong>with</strong> charm quarks<br />

<strong>for</strong> kaon mixing<br />

Further Reading:<br />

• [2, chapter 3.1-3.4]; [7]<br />

29


4.2 B 0 - ¯B 0 Mixing<br />

The description <strong>of</strong> mixing between neutral B d or B s mesons is similar to the <strong>for</strong>malism<br />

<strong>for</strong> kaon mixing. However, there are certain differences related to the different size <strong>of</strong><br />

the relevant CKM element and the different hadronic decay channels. Especially, in the<br />

B 0 - ¯B 0 system, one has<br />

∆Γ ≪ ∆M<br />

and the eigenstates will be classified by masses (heavy and light) rather than lifetimes,<br />

B H = pB 0 + q ¯B 0 ,<br />

B L = pB 0 − q ¯B 0 , (49)<br />

where we have introduced a new convenient parametrisation in terms <strong>of</strong> two complex<br />

parameters p, q which can be related to the ¯ɛ notation via<br />

p, q =<br />

1 ± ¯ɛ B<br />

√<br />

2(1 + |¯ɛB | 2 ) .<br />

As a consequence <strong>of</strong> Γ 12 ≪ M 12 in the B q -system (q = d, s), we can approximate<br />

and<br />

∆M q = M Bq<br />

H<br />

∆Γ q = Γ Bq<br />

H<br />

− M Bq<br />

L<br />

− ΓBq L<br />

(q)<br />

≃ 2|M 12 | ,<br />

(q)<br />

Re (M 12 Γ (q)∗<br />

12 )<br />

≃ 2<br />

|M (q)<br />

12 |<br />

(50)<br />

q<br />

p ≃ M (<br />

12<br />

∗ 1 − 1 |M 12 | 2 Im Γ )<br />

12<br />

≈ −arg M 12 . (51)<br />

M 12<br />

As a consequence,<br />

• the semi-leptonic asymmetry is tiny, a SL (B) ≈ 10 −4 , and there<strong>for</strong>e it is difficult<br />

measure this observable precisely and to quantify CP violation in B-meson mixing.<br />

• The ratio q/p is to a very good approximation a pure phase, i.e. |q/p| ≃ 1.<br />

Since, in the B-system, the top-quark always dominates the box diagrams, the<br />

weak phases are obtained from<br />

30


(M ∗ 12) (d) ∝ (V td V ∗<br />

tb) 2 ,<br />

(M ∗ 12) (s) ∝ (V ts V ∗<br />

tb) 2 (52)<br />

such that<br />

( q<br />

p)<br />

d,s<br />

≃ e 2iφd,s M<br />

<strong>with</strong> φ d M = −β , φ s M = −β s ≈ 0 (53)<br />

where β is one <strong>of</strong> the angles in the standard unitarity triangle, and β s the according<br />

tiny angle in the (squashed) triangle <strong>for</strong> b → s transitions.<br />

• ∆M d and ∆M s are directly proportional to the relevant CKM elements |V td | 2 and<br />

|V ts | 2 , respectively. The main uncertainties again stem from the hadronic matrix<br />

element <strong>of</strong> the operator in the ∆B = 2 weak effective Hamiltonian, which are<br />

parametrized analagously as <strong>for</strong> the kaon case by f 2 B q<br />

ˆBq .<br />

– ∆M d alone restricts the length <strong>of</strong> the right side <strong>of</strong> the UT, |V tb V ts /V cb V cs .<br />

Comparison <strong>with</strong> experiment gives the yellow band in Fig. 1, <strong>with</strong> a large<br />

uncertainty from f 2 B ˆB d .<br />

– A sizeable part <strong>of</strong> the uncertainties drops out in the ratio ∆M d /∆M s which<br />

again basically constrains the same side <strong>of</strong> the UT, since |V td /V ts | ≈ |V td /V cb |.<br />

Comparison <strong>with</strong> experiment yields the (narrower) orange band in Fig. 1.<br />

From Fig. 1, we see that ∆M d /∆M s and ɛ K alone are not sufficient to precisely determine<br />

¯ρ and ¯η.<br />

Time-dependent CP asymmetries at B-factories<br />

The new idea is to exploit the fact that the B d mesons have a relatively large life-time. In<br />

so-called asymmetric “B-factories”, one produces B- ¯B mesons <strong>with</strong> a large Lorentz<br />

boost which allows to identify the spatial displacement between the production and decay<br />

vertex <strong>of</strong> the B-mesons. This can be turned into in<strong>for</strong>mation <strong>of</strong> the time-evolution <strong>of</strong> (a<br />

statistical sample <strong>of</strong>) B-mesons. The idea <strong>of</strong> asymmetric B-factories is illustrated below:<br />

31


• e + e − collisions at the Υ(4S) resonance, M 2 e + e − ≃ M 2 Υ(4S) .<br />

• Υ(4S) almost exclusively decays coherently into B- ¯B mesons.<br />

• The decay <strong>of</strong> one <strong>of</strong> the mesons into a flavour eigenstate can be used to “tag” the<br />

charge <strong>of</strong> the b-quark at a time τ 1 defined by the corresponding displacement <strong>of</strong><br />

the vertex.<br />

• Exploiting the EPR paradoxon, this fixes the charge <strong>of</strong> the other b-quark at that<br />

time to be opposite.<br />

• Measuring the decay vertex <strong>of</strong> the other B-meson, we can thus test the quantum<br />

mechanical time evolution,<br />

|B meas. (τ 2 )〉 = U(τ 2 , τ 1 ) |B meas (τ 1 )〉<br />

In particular, we consider decays into CP-eigenstates |f〉, and define the time-dependent<br />

CP asymmetries <strong>for</strong> t = (τ 2 − τ 1 ) as<br />

a CP (t, f) = Γ(B0 (t) → f) − Γ( ¯B 0 (t) → f)<br />

Γ(B 0 (t) → f) + Γ( ¯B 0 (t) → f)<br />

(54)<br />

Now, the time-evolution <strong>of</strong> the flavour eigenstates follows from the evolution <strong>of</strong> the mass<br />

eigentstates,<br />

32


|B 0 (t)〉 = f + (t) |B 0 〉 + q p f −(t) | ¯B 0 〉 ,<br />

| ¯B 0 (t)〉 = p q f −(t) |B 0 〉 + f + (t) | ¯B 0 〉 , (55)<br />

<strong>with</strong><br />

f ± (t) = 1 2<br />

(<br />

e<br />

−iM H t ± e −iM Lt ) . (56)<br />

With this, the time-dependent CP asymmetry takes a relatively simple <strong>for</strong>m<br />

a CP (t, f) = A decay<br />

CP<br />

(B → f) cos(∆Mt) + Ainterf. CP (B → f) sin(∆Mt) (57)<br />

If we define<br />

ξ f ≡ q p<br />

A( ¯B 0 → f)<br />

A(B 0 → f) ≃ A( ¯B 0 → f)<br />

e2iφ M<br />

A(B 0 → f)<br />

(58)<br />

one obtains<br />

A decay<br />

CP<br />

= 1 − |ξ f| 2<br />

1 + |ξ f | , 2 Ainterf CP = 2Im ξ f<br />

(59)<br />

1 + |ξ f | 2<br />

The Golden Decay Mode B → J/ψK s<br />

The expression are particularly simple if we consider decay modes which are dominated<br />

by decay amplitudes <strong>with</strong> a single CKM factor. The well-known example is the decay<br />

B 0 / ¯B 0 → J/ψK S which is dominated by the operator inducing b → c¯cs transitions,<br />

involving V cb Vcs ∗ (being real in the standard parametrisation). As a consequence, <strong>for</strong><br />

such decays we have<br />

A( ¯B 0 → f)<br />

A(B 0 → f) = −η f e −2iφ decay<br />

where η f = ±1 refers to the CP-parity <strong>of</strong> the final state. From this, we have |ξ f | = 1,<br />

and there<strong>for</strong>e A decay<br />

CP<br />

= 0, while the interference term simply yields<br />

33


a interf<br />

CP (t, f) = η f sin 2β sin(∆Mt) (60)<br />

We thus have a simple theoretical prediction <strong>for</strong> a time-dependent oscillation, <strong>with</strong> a<br />

frequency set by the mass difference, and an amplitude yielding a direct measure <strong>for</strong> one<br />

<strong>of</strong> the CKM angles. Fitting the theoretical prediction to the experimental measurments<br />

(like the following one from the BaBar experiment),<br />

one obtains a precise value<br />

sin 2β = 0.681 ± 0.025<br />

which translates into the blue stripe in Fig. 1. Together <strong>with</strong> ɛ K and ∆M d /∆M s this<br />

already gives a satisfactory determination <strong>of</strong> ¯ρ and ¯η.<br />

Comment on sub-leading decay topologies in B → JψK<br />

• The main decay mechanism is through the current-current operator <strong>for</strong> b → c¯cs, yielding<br />

a “tree” amplitude<br />

T c V cb V ∗<br />

cs ∼ T c λ 2 .<br />

• The final state can also be reached via penguin operators <strong>for</strong> b → q¯qs (including q = c)<br />

P t V tb V ∗<br />

ts ∼ P t λ 2 .<br />

As the penguin operators have smaller Wilson coefficients, we have P t ≪ T c .<br />

• Finally, we also encounter current-current operators <strong>for</strong> b → uūs transitions,<br />

T u V ub V ∗ us ∼ T u λ 4 .<br />

As the uū quarks have to be turned into c¯c quarks by strong interactions, T u ≪ T c<br />

34


In total, we thus have<br />

A tot = T c V cb V ∗<br />

cs + P t V tb V ∗<br />

ts + T u V ub V ∗<br />

us = (T c − P t ) V cb V ∗<br />

cs + (T u − P t )V ub V ∗<br />

us<br />

≃ (T c − P t ) V cb V ∗<br />

cs # (61)<br />

where the neglected term is doubly suppressed by a relative factor λ 2 and by (T u −P t )/T c .<br />

<strong>Questions</strong>/Exercises:<br />

• Derive the expressions <strong>for</strong> A decay,interf.<br />

CP<br />

in (59).<br />

Further Reading:<br />

• [2, chapter 3.4-3.7,4.1,4.4]<br />

35


5 Semi-Leptonic Decays<br />

Semi-leptonic decay amplitudes q → q ′ lν are proportional to V qq ′. From semi-leptonic<br />

B-meson decays, we can extract |V cb | and |V ub |.<br />

(a) exclusive B → π(ρ)lν yields |V ub |; exclusive B → D(D ∗ )lν yields |V cb |<br />

(b) inclusive B → X u lν yields |V ub |, inclusive B → X c lν yields |V cb |<br />

(c) exclusive B + → τ + ν τ yields |V ub |; (exclusive B c → τ + ν τ yields |V cb |)<br />

5.1 |V cb | from exclusive decays and HQET<br />

In order to determine |V cb | we need, <strong>for</strong> instance,<br />

〈D ∗ lν|Heff<br />

b→clν |B〉 = G F<br />

√ 〈l¯ν|¯lγ µ (1 − γ 5 )ν|0〉〈D ∗ |¯cγ µ (1 − γ 5 )b|B〉 (62)<br />

2<br />

where the leptonic matrix element is straight-<strong>for</strong>ward, while the hadronic matrix element<br />

involves a B → D ∗ transition <strong>for</strong>m factor which parametrizes the non-perturbative QCD<br />

dynamics. The idea is to exploit the hierarchy Λ QCD ≪ m b,c . In the limit m b,c → ∞,<br />

the heavy quarks in the transitions act as “static sources <strong>of</strong> a colour field”, and one<br />

expects analogous simplifications as <strong>for</strong> the transitions between hydrogen-like atoms in<br />

non-relativistic quantum mechanics.<br />

Heavy-Quark Effective Theory (HQET)<br />

Within relativistic quantum field theory, we can <strong>for</strong>malize the heavy quark limit by<br />

constructing an appropriate effective theory. In contrast to the case <strong>of</strong> the effective weak<br />

Hamiltonian, the heavy quark field is not completely removed from the theory. To this<br />

end, one first per<strong>for</strong>ms a field redefinition <strong>for</strong> (Q = b, c)<br />

Q(x) = e −im Q v·x (h v (x) + H v (x)) (63)<br />

<strong>with</strong><br />

h v = e im Q v·x 1 + /v<br />

2<br />

H v = e im Q v·x 1 − /v<br />

2<br />

Q(x) ,<br />

Q(x) (64)<br />

Here, one has factored out the plane-wave behaviour corresponding to the on-shell motion<br />

<strong>of</strong> a heavy quark (i.e. the static limit in the heavy quark’s rest frame). Furthermore, v µ<br />

36


efers to the 4-velocity <strong>of</strong> the heavy quark (v 2 = 1), and (1±/v)/2 are projectors in Dirac<br />

space.<br />

Inserting the new parametrization into the QCD Dirac-Lagrangian, and using the projection<br />

properties, one obtains<br />

L = (¯hv + ¯H v<br />

)<br />

e<br />

im Q v·x (i /D − m Q ) e −im Qv·x (h v + H v )<br />

= ¯h v (iv · D)h v − ¯H v (iv · D + 2m Q )H v + ¯h v i /D ⊥ H v + ¯H v i /D ⊥ h v . (65)<br />

where /D ⊥ = /D − /v (v · D). Now, only the kinetic term <strong>for</strong> the projection H v contains<br />

the heavy-quark mass m Q , and we can thus eliminate this field by making use <strong>of</strong> the<br />

classical equations <strong>of</strong> motion (which corresponds to tree-level matching <strong>of</strong> the effective<br />

theory),<br />

∂L<br />

∂ ¯H v<br />

= −(iv · D + 2m Q )H v + i /D ⊥ h v = 0<br />

⇒H v =<br />

1<br />

iv · D + 2m Q<br />

i /D ⊥ h v ∼ O<br />

(<br />

ΛQCD<br />

)<br />

h v (66)<br />

m Q<br />

where the power counting is obtained under the assumption that the residual modulation<br />

<strong>of</strong> the heavy-quark field <strong>with</strong> respect to the plane-wave factor is generated by interactions<br />

<strong>with</strong> s<strong>of</strong>t gluons, having momenta/field strength <strong>of</strong> order Λ QCD . There<strong>for</strong>e<br />

H v : “small component” <strong>of</strong> the heavy-quark Dirac spinor,<br />

h v : “large component” <strong>of</strong> the heavy-quark Dirac spinor.<br />

If we now expand the non-local term 1/(iv·D+2m Q ) in powers <strong>of</strong> Λ QCD /m Q , we obtain an<br />

effective theory where the field H v does not propagate anymore. As in the example <strong>of</strong> the<br />

weak effective Hamiltonian, this will change the UV properties <strong>of</strong> Feynman integrals, and<br />

the resulting Lagrangian corresponds to an effective theory <strong>for</strong> heavy quarks interacting<br />

<strong>with</strong> s<strong>of</strong>t gluons and light quarks.<br />

“Heavy Quark Effective Theory” (HQET)<br />

The quantum corrections from hard momenta (<strong>of</strong> order m b ), again, will be absorbed into<br />

Wilson coefficients <strong>of</strong> HQET operators,<br />

Let us look at the leading-order terms (at tree-level) by inserting H v ≃ i /D ⊥<br />

2m Q<br />

h v ,<br />

37


L HQET = ¯h v (iv · D) h v + ¯h v i /D ⊥<br />

1<br />

2m Q<br />

i /D ⊥ h v + . . . (67)<br />

• The leading term in this expansion<br />

– is independent <strong>of</strong> the heavy quark mass<br />

⇒ (approximate) Heavy Quark Flavour Symmetry<br />

– is independent <strong>of</strong> heavy quark spin orientation (no magnetic interaction <strong>with</strong><br />

gluon field)<br />

⇒ (approximate) Heavy Quark Spin Symmetry<br />

• The O(1/m Q ) term can be rewritten, using<br />

leading to<br />

/D ⊥ /D ⊥ = (D ⊥ ) 2 + 1 4 [γ µ, γ ν ] [D µ ⊥ , Dν ⊥] = (D ⊥ ) 2 + g 2 σ µνG µν<br />

L HQET<br />

∣<br />

∣1/mQ = ¯h v<br />

{ (iD⊥ ) 2<br />

2m Q<br />

− g<br />

4m Q<br />

σ µν G µν }<br />

h v (68)<br />

– The first term represents the (generalization <strong>of</strong> the) non-relatistic kinetic energy.<br />

It breaks the HQET flavour symmetry. Its Wilson coefficient does not<br />

receive radiative corrections because it is protected by Lorentz symmetry.<br />

– The second term represents the (generalization <strong>of</strong> the) non-relatistic Pauli<br />

term. It breaks the HQET flavour and spin symmetries. Its Wilson coefficient<br />

receives radiative corrections, and in the leading-log-approximation, it is given<br />

by<br />

C mag (µ) ∣ ( ) 9/(33−2nf )<br />

αs (m b )<br />

LLA<br />

=<br />

α s (µ)<br />

38


Consequences <strong>for</strong> the the Hadron Spectrum<br />

Consider pseudoscalar or vector mesons containing one heavy quark<br />

H ≡ (Q¯q) spin−0 ,<br />

H ∗ ≡ (Q¯q) spin−1<br />

In the heavy-mass limit, they would correspond to a fundamental doublet <strong>of</strong> the heavyquark<br />

spin symmetry. The masses <strong>of</strong> the meson can thus be written as an expansion in<br />

1/m Q :<br />

m H = m Q + ¯Λ − λ 1<br />

2m Q<br />

− 3λ 2<br />

2m Q<br />

+ . . .<br />

m H ∗ = m Q + ¯Λ − λ 1<br />

2m Q<br />

+ λ 2<br />

2m Q<br />

+ . . . (69)<br />

• m Q comes from the heavy-quark mass itself,<br />

• ¯Λ parametrizes the expectation value <strong>of</strong> the leading HQET operator ¯h v (ivD)h v ,<br />

• λ 1 parametrizes the 1/m Q correction from the kinetic energy operator,<br />

• λ 2 parametrizes the 1/m Q correction from the magnetic term, where the different<br />

coefficients <strong>for</strong> pseudoscalar and vector mesons are given by Clebsch-Gordan<br />

coefficients <strong>of</strong> SU(2).<br />

As a consequence, one obtains simple approximate relations, like<br />

or<br />

m 2 H ∗ − m2 H = const.<br />

m 2 B ∗ − m2 B ≃ 0.49 GeV 2 ,<br />

m 2 D ∗ − m2 D ≃ 0.55 GeV 2 , (70)<br />

m Hs − m Hq<br />

= ¯Λ s − ¯Λ q = const.<br />

Feynman rules in HQET<br />

m Bs − m B ≃ m Ds − m D ≃ 100 MeV . (71)<br />

From the leading HQET Lagrangian, ¯h v (ivD)h v , one reads <strong>of</strong>f the Feynman rules <strong>for</strong><br />

• the heavy-quark propagator:<br />

i<br />

v·k+iɛ<br />

• heavy-quark gluon vertex: ig s v µ t A 39


Alternatively, we can obtain the HQET Feynman rules from the corresponding QCD<br />

Feynman rules. We decompose the heavy quark momentum as<br />

p µ = m Q v µ + k µ<br />

where k µ is the “residual momentum” <strong>of</strong> the heavy quark. Then the heavy-quark<br />

propagator can be approximated as<br />

/p + m Q<br />

p 2 − m 2 Q + iɛ = m Q(1 + /v) + /k<br />

2m Q v · k + k 2 + iɛ ≃ 1 1 + /v<br />

v · k + iɛ 2<br />

(72)<br />

which reproduces the HQET propagator times the heavy-quark projector (which is implicit<br />

in HQET).<br />

The heavy-quark gluon vertex is reproduced by realizing that the vertex always appears<br />

between heavy quark propagator or on-shell heavy quark spinors, which both provide a<br />

projector 1+/v<br />

. One then has<br />

2<br />

1 + /v<br />

2<br />

ig s γ µ t a 1 + /v<br />

2<br />

= ig s v µ t a 1 + /v<br />

2<br />

. (73)<br />

• The HQET Feynman rules thus reproduce the s<strong>of</strong>t region <strong>of</strong> Feynman integrals,<br />

where the residual momentum <strong>of</strong> the heavy quark k µ is small compared to its mass.<br />

• The hard momentum region are accounted <strong>for</strong> by Wilson coefficients <strong>for</strong> the subleading<br />

operators in the HQET Lagrangian and external decay currents (see below).<br />

Heavy-Quark Currents<br />

In terms <strong>of</strong> full QCD fields, the semi-leptonic b → clν transitions are induced by the<br />

currents<br />

J QCD = ¯c(x) γ µ (1 − γ 5 ) b(x) ≡ ¯cΓb<br />

At leading-order in the HQET expansion, this can be replaced by<br />

J QCD ≃ e −imbv·x imcv′·x<br />

e ¯h(c)<br />

v<br />

(x)Γh (b)<br />

′ v (x) ≡ e −imbv·x e imcv′·x J HQET (74)<br />

Similarly, we may replace the hadronic states by the corresponding states in the heavymass<br />

limit,<br />

|B(p)〉 → √ m B |B(v)〉 etc.<br />

where we have factored out the mass-dependent normalization <strong>for</strong> convenience.<br />

40


To derive the consequences <strong>of</strong> the heavy quark symmeties, we first study the representation<br />

<strong>of</strong> heavy-meson states under Lorentz trans<strong>for</strong>mations<br />

<strong>with</strong><br />

M (Q)<br />

v<br />

(Q α¯q β ) ∼ ( )<br />

M (Q) αβ<br />

v<br />

→ D(Λ)M (Q)<br />

v D(Λ) −1<br />

under Lorentz trans<strong>for</strong>mations, where D(Λ) determines the trans<strong>for</strong>mations <strong>of</strong> the fundamental<br />

Dirac spinors.<br />

• As a consequence <strong>of</strong> the HQET spin symmetry, we can combine pseudoscalar and<br />

vector meson states into a common (momentum-space) representation matrix <strong>for</strong><br />

the spin doublet,<br />

M (Q)<br />

v = 1 + /v [ɛ/(v, λ) a V (v, λ) + iγ 5 a P (v, λ)] (75)<br />

2<br />

which includes the HQ-projector, the polarization vector ɛ <strong>for</strong> a vector meson, as<br />

well as the annihilation operators <strong>for</strong> vector or pseudoscalar meson states. Notice<br />

that M (Q)<br />

v /v = −M (Q)<br />

v .<br />

• In terms <strong>of</strong> this representation matrix, the hadronic matrix element <strong>of</strong> the decay<br />

current can be written as<br />

〈D(v ′ )|¯h (c)<br />

v<br />

Γh (b)<br />

′ v<br />

[<br />

|B(v)〉 ∝ tr Ξ(v, v ′ )<br />

]<br />

(c) ¯M<br />

v<br />

ΓM (b)<br />

′ v<br />

(76)<br />

Here the Dirac matrix Ξ(v, v ′ ) parametrizes the (unspecified) bound-state dynamics<br />

<strong>of</strong> the light degrees <strong>of</strong> freedom, whereas the remaining terms are dictated by the<br />

HQET and Lorentz symmetries.<br />

• Now, Ξ can only be a function <strong>of</strong> v and v ′ , and due to the heavy-quark projection<br />

properties and the parity conservation <strong>of</strong> strong interactions, all allowed Dirac<br />

structures can be reduced to the unit matrix,<br />

Ξ(v, v ′ ) ∼ ξ(v · v ′ ) 1<br />

With this, all transition matrix element between pseudoscalar and vector B or D-mesons,<br />

can be expressed by a unique “Isgur-Wise“ Form Factor ξ(ω = v · v ′ ),<br />

41


〈D(v ′ )|¯c v ′γ µ b v |B(v)〉 = ξ(ω) ( v µ + v ′ µ)<br />

,<br />

〈D ∗ (v ′ , ɛ)|¯c v ′γ µ γ 5 b v |B(v)〉 = −iξ(ω) ( (1 + ω)ɛ ∗ µ − (v · ɛ ∗ )v ′ µ)<br />

,<br />

〈D ∗ (v ′ , ɛ)|¯c v ′γ µ b v |B(v)〉 = ξ(ω) ε µναβ ɛ ∗ν v ′α v β .<br />

We may also consider the case where Q = Q ′ = b and v = v ′ (”Isgur-Wise limit“, ω = 1).<br />

This corresponds to the matrix element <strong>of</strong> a conserved current,<br />

〈B(v)|¯h (b)<br />

v γ µ h (b)<br />

v |B(v)〉 = 2ξ(ω = 1) v µ<br />

≃ 〈B(p)|¯bγ µ b|B(p)〉/m B = 2F (q 2 = 0)p µ /m B (77)<br />

such that<br />

ξ(1) = F (0) = 1<br />

In semi-leptonic decays, the Isgur-Wise limit ω = 1 corresponds to the maximal momentum<br />

transfer to the lepton pair,<br />

q 2 = (p B − p D ) 2 = (m B − m D ) 2 = q 2 max .<br />

Then, the CKM element |V cb | can be extracted from<br />

|V cb | ≃ |F (q2 → q 2 max) V cb | exp.<br />

|F (v · v ′ = 1)| theor.<br />

(78)<br />

where<br />

|F (v · v ′ = 1)| theor. = 1 + O(1/m Q ) + O(α s ) (79)<br />

In particular, it turns out that <strong>for</strong> B → D ∗ decays, there are no 1/m Q corrections at<br />

ω = 1 (”Luke’s theorem“), i.e. the power corrections start at O(1/m b · 1/m c ). The<br />

relevant <strong>for</strong>m factor including the corrections is then estimated as<br />

F D ∗(ω = 1) = 0.91 ± 0.05<br />

42


which means a 10% deviation from unity, <strong>with</strong> a relative hadronic uncertainty <strong>of</strong> 50%,<br />

implying a 5% theoretical uncertainty on the extraction <strong>of</strong> |V cb |.<br />

<strong>Questions</strong>/Exercises:<br />

• Prove that (1 ± /v)/2 are projectors in Dirac space. Verify that in the heavyquark’s<br />

rest frame, the projections h v and H v correspond to the particle and<br />

anti-particle components <strong>of</strong> the Dirac spinor.<br />

• Calculate<br />

1 ± /v<br />

2<br />

γ µ<br />

1 ± /v<br />

2<br />

and 1 ∓ /v<br />

2<br />

γ µ<br />

1 ± /v<br />

2<br />

• Verify the derivation <strong>of</strong> the HQET Lagrangian and derive the Feynman rules<br />

<strong>for</strong> the HQET propagator, the leading heavy-quark gluon vertex, and the<br />

sub-leading quark-gluon vertices from the kinetic and chromomagnetic operator.<br />

• Verify (77).<br />

Further Reading:<br />

• [8, chapter 2.5, 2.6, 2.9, 3.4, 4.3]<br />

.<br />

5.2 |V cb | from inclusive decays<br />

In the inclusive decay B → X c l¯ν, we do not specify the final state completely, except<br />

<strong>for</strong> demanding open charm in the final state. The situation is somewhat similar to<br />

deep inelastic scattering: Instead <strong>of</strong> elastic parton-electron scattering, the partonic subprocess<br />

is given by the weak b → cl¯ν decay.<br />

The differential decay rate can be written in terms <strong>of</strong> a ”hadronic tensor W µν “ and a<br />

”leptonic tensor L µν “,<br />

d 3 Γ<br />

dq 2 dE l dE ν<br />

= G2 F<br />

4π 3 |V cb| 2 W µν L µν (80)<br />

<strong>with</strong> q µ L µν = q ν L µν = 0 <strong>for</strong> m l = 0 (i.e. <strong>for</strong> electrons and muons).<br />

• The computation <strong>of</strong> L µν is straight<strong>for</strong>ward<br />

L µν = 1 2 tr [/p νγ µ /p l γ ν (1 − γ 5 )]<br />

• The hadronic tensor contains the hadronic matrix element and phase-space factors<br />

43


W µν = (2π)3<br />

2m B<br />

∑ ∫<br />

X c<br />

δ (4) (p B − q − p X )<br />

× 〈 ¯B|¯bγ µ P L c|X c 〉〈X c |¯cγ ν P L b| ¯B〉 (81)<br />

and it can be decomposed into the most general Lorentz structures<br />

W µν ≡ −g µν W 1 + v µ v ν W 2 + iɛ µναβ v α q β W 3<br />

+ q µ q ν W 4 + (v µ q ν + q µ v ν )W 5 (82)<br />

The W i are real function <strong>of</strong> q 2 and v · q. W 4,5 do not contribute when contracted<br />

<strong>with</strong> L µν (<strong>for</strong> m l = 0).<br />

• The aim is now to calculate the W i in the heavy-mass limit in terms <strong>of</strong> HQET parameters<br />

and α s (m b ). The general procedure is to make use <strong>of</strong> the optical theorem<br />

to relate the decay rate to the imaginary part <strong>of</strong> the <strong>for</strong>ward-scattering amplitude,<br />

where<br />

W i = − 1 π Im T i<br />

T µν =<br />

−i ∫<br />

2m B<br />

d 4 x e −iq·x 〈 ¯B|T [ J † µ(x)J ν (0) ] | ¯B〉 (83)<br />

Now, the time-ordered product <strong>of</strong> the two currents can be computed in perturbation<br />

theory, by contracting the quark fields to propagators. Formally, this amounts to<br />

an ”Operator Product Expandion“ (OPE), where, in momentum space, the<br />

product <strong>of</strong> two operators at different space-time points is expanded into a series <strong>of</strong><br />

local operators <strong>with</strong> the corresponding Wilson coefficients.<br />

∫<br />

−i<br />

d 4 x e −iqx T [ J † µ(x)J ν (0) ] ≡ ∑ i<br />

C i (q) O i (84)<br />

The operator <strong>with</strong> the lowest dimension, contributing to the imaginary part <strong>of</strong> T µν<br />

has the <strong>for</strong>m ¯bΓb and gets a tree-level contribution from contracting the charmquark<br />

fields to a charm propagator. The Wilson coefficient is obtained from match-<br />

44


ing the expression<br />

ū(p b )γ µ P L<br />

/p b − q/<br />

(p b − q) 2 − m 2 c + iɛ γ νP L u(p b )<br />

From this one immediately sees that Im T i corresponds to an on-shell charm quark,<br />

(p b − q) 2 − m 2 c = 0, in other words at leading order one has<br />

LO:<br />

dΓ(B → X c lν) ∼ dΓ(b → clν)<br />

i.e. the hadronic decay rate can be approximated by the partonic one. The matrix<br />

elements <strong>of</strong> the local operators, 〈 ¯B|¯bΓb| ¯B〉, are simply determined by the trivial<br />

probability to find a b-quark in the B-meson, and do not introduce a new hadronic<br />

parameter (at this order), apart from the b-quark mass (see below).<br />

Working out the Dirac and Lorentz algebra, the final result <strong>for</strong> the LO differential<br />

rate then reads<br />

d 3 Γ<br />

dE e dE ν dq = G2 F |V [<br />

] (<br />

cb| 2<br />

− q2 E ν<br />

+ E 2 2π 3 e E ν δ E e + E ν − q2 + m 2 b − )<br />

m2 c<br />

2m b 2m b<br />

(85)<br />

• For instance, we can study the electron-energy spectrum by integrating over E ν and<br />

q 2 . Introducing the dimensionless variables y = 2E e /m b and ρ = m 2 c/m 2 b ∼ 0.08,<br />

one obtains<br />

dΓ<br />

dy = G2 F |V [<br />

]<br />

cb| 2 m 5 b<br />

3y 2 − 2y 3 − 3y 2 ρ − 3y2 ρ 2<br />

96π 3 (1 − y) + (3y2 − y 3 )ρ 3<br />

θ(1 − ρ − y)<br />

2 (1 − y) 3<br />

(86)<br />

45


– Notice the strong dependence <strong>of</strong> the decay width on the b-quark mass. Remember<br />

that m b is a scheme- and scale-dependent quantity, which means that<br />

the related uncertainties can only be resolved by calculating the higher orders<br />

in perturbation theory.<br />

• The experimental procedure to extract a value <strong>for</strong> |V cb | from the comparison <strong>of</strong><br />

experimental data and theoretical predictions goes as follows:<br />

– Compute the α s corrections to the partonic rate.<br />

– Include higher orders in the OPE, <strong>with</strong> operators <strong>of</strong> the <strong>for</strong>m<br />

¯bΓ(iD⊥ ) · · · (ivD) · · · (iD ⊥ )b<br />

which can be related to HQET parameters ¯Λ, λ 1 , λ 2 etc.<br />

– measure moments 〈y n 〉 <strong>with</strong> n = 0..3, <strong>of</strong> the electron-energy spectrum<br />

– compare <strong>with</strong> the theoretical <strong>for</strong>mula and simultaneously fit the HQET parameters<br />

and |V cb |.<br />

– Notice that the final result will depend on the renormalization scheme.<br />

• The result <strong>of</strong> such analysis on average leads to [9]<br />

|V cb | incl. = (41.5 ± 0.7) · 10 −3 (87)<br />

which should be compared <strong>with</strong> the number obtained from exclusive B → D/D ∗ lν<br />

decays<br />

|V cb | incl. = (38.7 ± 1.1) · 10 −3 (88)<br />

The slight discrepancy between the two numbers can be taken as a statistical fluctuation<br />

in the data, an underestimate <strong>of</strong> perturbative and non-perturbative theoretical<br />

uncertainties, or even an effect <strong>of</strong> new-physics contributions (<strong>for</strong> instance,<br />

from right-handed weak currents).<br />

<strong>Questions</strong>/Exercises:<br />

• Calculate the (normalized) moments 〈y n 〉/〈y 0 〉 (n = 1..3) <strong>of</strong> the electronenergy<br />

spectrum, approximating ρ = 0.<br />

46


Further Reading:<br />

• [8, chapter 6.1-6.2]; [9]<br />

5.3 |V ub | from inclusive decays and SCET<br />

• The determination <strong>of</strong> |V ub | from exclusive decay modes like B → π(ρ)lν heavily<br />

relies on non-perturbative methods to estimate the <strong>for</strong>m factors. The PDG quotes<br />

<strong>with</strong> an uncertainty 10 − 15%.<br />

|V ub | excl. = 3.38 · 10 −3<br />

• From the theoretical point <strong>of</strong> view, the determination from hte inclusive decays<br />

¯B → X u l¯ν seems more favourable, as the same theoretical framework (OPE) as<br />

<strong>for</strong> B → X c lν can be used.<br />

• However, because |V cb | 2 ≫ |V ub | 2 , the b → clν decays actually represent a large<br />

background <strong>for</strong> the small b → ulν signal.<br />

• The usual experimental strategy is to introduce kinematic cuts, e.g.<br />

M 2 X < M 2 X c,min ≃ M 2 D .<br />

• Un<strong>for</strong>tunately, <strong>for</strong> such kinematic cuts, most <strong>of</strong> the remaining events <strong>for</strong> B → X u lν<br />

correspond to the final state X u <strong>for</strong>ming hadronic jets, i.e. the energy E X is <strong>of</strong><br />

order m b /2 while the invariant mass M 2 X ≪ m2 b .<br />

(The particles <strong>for</strong>ming the jet move almost collinear to each other – in the context <strong>of</strong> the effective<br />

theory to be constructed, we sometimes call these particle modes ”hard-collinear“.)<br />

m b v + s<strong>of</strong>t<br />

¯ν l<br />

l<br />

B<br />

s<strong>of</strong>t<br />

hard-collinear (jet)<br />

X u<br />

To get further insight, let us consider the LO diagram <strong>for</strong> the functions T i again:<br />

47


– The momentum <strong>of</strong> the (now up-quark) propagator in the <strong>for</strong>ward-scattering<br />

amplitude is given by<br />

∆ µ = (m b v µ + k µ ) − q µ = (m b v µ + k µ ) − (m B v m u − p µ X ) = pµ X + (kµ − ¯Λv µ )<br />

(89)<br />

The virtuality <strong>of</strong> the propagator follows as (<strong>for</strong> m u ≃ 0)<br />

∆ 2 = p 2 X + 2k · p X − 2¯Λv · p X + O(Λ 2 QCD) (90)<br />

– The OPE case discussed be<strong>for</strong>e in the context <strong>of</strong> B → X c lν corresponds to<br />

hadronic final states which are dominated by momenta satisfying p 2 X ∼ m2 b<br />

(notice, that the larger |p X |, the larger the available phase space, and there<strong>for</strong>e events <strong>with</strong><br />

smaller p 2 X are less important)<br />

In this case, we can approximate ∆ 2 ≃ p 2 X , and the partonic decay rate does<br />

not depend on the residual momentum k µ which justifies the OPE in terms<br />

<strong>of</strong> local operators.<br />

– In contrast, the jet configuration corresponds to events, where<br />

p 2 X ∼ m b Λ ≪ m 2 b<br />

”hard-collinear“<br />

where <strong>for</strong> the considered case, we can consider Λ ∼ Λ QCD . In this case, the<br />

terms depending on the residual heavy-quark momentum in (90) cannot be<br />

neglected anymore.<br />

As a consequence, the decay rate will depend on the distribution <strong>of</strong> the b-<br />

quark momentum k µ in the B-meson, which can be described by a parton<br />

distribution function (as in DIS), dubbed ”Shape Function“ in the context<br />

<strong>of</strong> inclusive B decays.<br />

Phenomenology <strong>for</strong> |V ub | incl.<br />

• One can construct models <strong>for</strong> the shape functions, implementing theory constraints<br />

from the OPE. The variance in ”reasonable“ model functions gives an estimate <strong>for</strong><br />

the related systematic theoretical uncertainties.<br />

48


• Alternatively, one can extrac the shape-function from other inclusive decays, noteably<br />

the radiative decay B → X s γ, which is based on the partonic rate <strong>for</strong> b → sγ,<br />

and the shape-function region corresponds to large photon energies.<br />

The PDG quotes a value<br />

|V ub | incl. = 4.27 ± 0.38 · 10 −3<br />

which is somewhat higher than the the exclusive value.<br />

The effect on the CKM triangle fit is via the ratio |V ub |/|V cb | which enters the length <strong>of</strong><br />

the left side <strong>of</strong> the triangle, see the dark green band in the figure below:<br />

As one can see, the determination <strong>of</strong> |V ub /V cb | from semi-leptonic decays (in particular<br />

from the inclusive rates) favour slightly larger values than the ones corresponding to the<br />

overall CKM triangle fit.<br />

S<strong>of</strong>t-collinear effective theory (SCET)<br />

The theoretical description <strong>of</strong> the shape-function region can be <strong>for</strong>mulated in terms <strong>of</strong><br />

an effective theory which simultaneously reproduces the dynamics <strong>of</strong> s<strong>of</strong>t and (hard-<br />

)collinear particle fields, ”s<strong>of</strong>t-collinear effective theory“ (SCET).<br />

Starting point is the decomposition <strong>of</strong> Lorentz vectors according to ”light-cone kinematics“,<br />

introducing light-like vectors (given in the B-meson rest frame)<br />

n µ ± = (1, 0, 0, ±1) , n + · n − = 2 , n 2 + = n 2 − = 0 (91)<br />

which allow to decompose any vector as<br />

p µ = (n + p) nµ −<br />

2 + pµ ⊥ + (n −p) nµ +<br />

2 , p2 = (n + p)(n − p) + p 2 ⊥ (92)<br />

49


We distinguish:<br />

• (hard-)collinear particles (in jets), <strong>with</strong><br />

|n + p X | ∼ m b ≫ |p X,⊥ | ∼ √ m b Λ ≫ |n − p X | ∼ Λ ,<br />

p 2 X ∼ m b Λ , (93)<br />

• s<strong>of</strong>t particles (from B-meson), <strong>with</strong>in<br />

|n + k| ∼ |k ⊥ | ∼ |n − k| ∼ Λ , k 2 ∼ Λ 2 . (94)<br />

The virtuality <strong>of</strong> the up-quark propagator in the leading diagram <strong>for</strong> the T i in b → ulν<br />

is then given by<br />

∆ 2 ≃ p 2 X + (n − k − ¯Λ)(n + p X )<br />

which shows that we actually need the distribution <strong>of</strong> the light-cone component (n − k)<br />

<strong>of</strong> the residual b-quark momentum in the B-meson, f (B)<br />

b<br />

(n − k), defining the relevant<br />

shape-function (or pdf).<br />

The effective theory again should reproduce the low-energy (or, as we will see, more<br />

precisely low-virtuality) region(s) <strong>of</strong> loop integrals in perturbation theory, whereas the<br />

hard momentum region will be absorbed into coefficients (more precisely, coefficient<br />

functions). As an example, consider the 1-loop correction to the b → u decay vertex:<br />

• In full QCD, the resulting Feynman integral has the structure:<br />

∫<br />

d D l<br />

(2π) [numerator] 1<br />

D l 2 + iɛ ·<br />

1<br />

(m b v + k − l) 2 − m 2 b + iɛ ·<br />

1<br />

(∆ − l) 2 + iɛ<br />

(95)<br />

where l is the loop-momentum carried by the gluon, and the result will be a nonanalytic<br />

function <strong>of</strong> ∆ 2 /m 2 b and (v · k)/m b.<br />

50


• The hard momentum region is obtained by expanding the integrand under the assumption<br />

|l µ | ∼ m b ∼ (n + ∆), and putting all sub-leading momentum components<br />

(∆ ⊥ , (n − ∆), k) to zero.<br />

hard:<br />

∫<br />

× 1<br />

l 2 + iɛ ·<br />

d D l<br />

(2π) D [numerator]<br />

1<br />

(m b v − l) 2 − m 2 b + iɛ ·<br />

1<br />

l 2 − (n + ∆)(n − l) + iɛ<br />

(96)<br />

As usual, we will absorb the contribution <strong>of</strong> the hard region into a Wilson coefficient<br />

<strong>for</strong> the current operator in the effective theory. An important difference to the cases<br />

discussed so far is, that the Wilson coefficient now depends on the large momentum<br />

component (n + ∆) <strong>of</strong> the u-quark jet,<br />

C i = C i (µ, m b , n + ∆) (97)<br />

• The hard-collinear region corresponds to the case<br />

which leads to<br />

|n + l| ∼ m b ∼ |n + ∆| ≫ |l ⊥ | ∼ |∆ ⊥ | ≫ |n − l| ∼ |n − ∆| ∼ |k|<br />

hard-collinear:<br />

× 1<br />

l 2 + iɛ ·<br />

∫<br />

d D l<br />

(2π) D [numerator]<br />

1<br />

−m b (n + l) + iɛ ·<br />

1<br />

(∆ − l) 2 + iɛ<br />

(98)<br />

It should be reproduced by the Feynman rules in SCET (including the unconventional<br />

factor 1/(n + l)). The result <strong>for</strong> the hard-collinear region is independent <strong>of</strong><br />

the heavy-quark residual momentum, and only depends on the virtuality ∆ 2 <strong>of</strong> the<br />

u-quark jet. It thus represents a universal perturbative correction to the propagator<br />

<strong>of</strong> the hard-collinear u-quark. In SCET, these effects are collected in the<br />

so-called ”Jet Function“<br />

51


(hard-collinear) Jet-Function: J = J(µ, ∆ 2 )<br />

• Finally, the s<strong>of</strong>t region <strong>of</strong> the integral corresponds to |l µ | ∼ |k µ ∼ Λ QCD , which<br />

gives rise to<br />

s<strong>of</strong>t:<br />

∫<br />

d D l<br />

(2π) D [numerator]<br />

× 1<br />

l 2 + iɛ ·<br />

1<br />

2v · (k − l) + iɛ ·<br />

1<br />

−(n + ∆)(n − l) + iɛ<br />

(99)<br />

which contains the familiar HQET propagator <strong>for</strong> a quasi-static b-quark, and again<br />

an unconventional 1/(n − l) factor to be reproduced by the SCET Feynman rules.<br />

The result <strong>for</strong> the s<strong>of</strong>t integral depends only on the residual momentum k (apart<br />

from a global kinematic 1/(n + ∆) factor) and thus represents an inherent property<br />

<strong>of</strong> the B-meson which can be attributed to a non-trivial ”Shape Function“<br />

(s<strong>of</strong>t) Shape-Function: S = S(µ, n − k)<br />

describing the distribution <strong>of</strong> the residual b-quark momentum in the B-meson.<br />

Putting everything together, one ends up <strong>with</strong> a ”Facotrization Theorem“, which<br />

takes the schematic <strong>for</strong>m<br />

dΓ ∝ C · J ⊗ S (100)<br />

which takes the <strong>for</strong>m <strong>of</strong> a convolution in the residual-momentum variable (n − k) (every<br />

allowed value <strong>of</strong> (n − k) <strong>for</strong> fixed external kinematics (q 2 , E e , E ν ) corresponds to another<br />

value <strong>of</strong> ∆ µ ). Here the Wilson coefficient C can be calculated in full QCD, the jet<br />

function J can be calculated in SCET, and the shape-function S can be caluclated<br />

in HQET. SCET/HQET also describes the dependence <strong>of</strong> the various functions under<br />

changes <strong>of</strong> the factorization scale µ, and there<strong>for</strong>e the above <strong>for</strong>mula can be exploited to<br />

calculate the decay rate in renormalization-group improved perturbation theory, where<br />

52


• The hard function is matched at µ ∼ m b and evolved down to µ ∼ √ ∆ 2 .<br />

• The jet function is matched at µ ∼ √ ∆ 2 .<br />

• The s<strong>of</strong>t function is modelled at µ ∼ 1−2 GeV and evolved (upwards) to µ ∼ √ ∆ 2 .<br />

<strong>Questions</strong>/Exercises:<br />

• Compare the analysis <strong>for</strong> the b → u vertex in the shape-function region in<br />

SCET, <strong>with</strong> the corresponding vertex correction <strong>for</strong> the b → c current in<br />

HQET. Verify that the required Wilson coefficient does not depend on the<br />

residual momentum <strong>of</strong> the heavy quark.<br />

Further Reading:<br />

• [10] and refs. therein.<br />

53


References<br />

[1] G. Buchalla, A. J. Buras, M. E. Lautenbacher, Rev. Mod. Phys. 68 (1996) 1125-<br />

1144. [hep-ph/9512380].<br />

[2] A. J. Buras, “Flavor dynamics: CP violation and rare decays,” [hep-ph/0101336].<br />

[3] The CKM Quark-Mixing Matrix, in C. Amsler et al., Phys. Lett. B667, 1 (2008),<br />

Available online at:<br />

http://pdg.lbl.gov/2010/reviews/rpp2010-rev-ckm-matrix.pdf.<br />

[4] M. E. Peskin, D. V. Schroeder, Reading, USA: Addison-Wesley (1995) 842 p.<br />

[5] M. Beneke, G. Buchalla, M. Neubert, C. T. Sachrajda, Nucl. Phys. B591 (2000)<br />

313-418. [hep-ph/0006124].<br />

[6] M. Neubert, “<strong>Lecture</strong>s on the theory <strong>of</strong> nonleptonic B decays,” [hep-ph/0012204].<br />

[7] CP violation in meson decays, in C. Amsler et al., Phys. Lett. B667, 1 (2008),<br />

Available online at:<br />

http://pdg.lbl.gov/2010/reviews/rpp2010-rev-cp-violation.pdf.<br />

[8] A. V. Manohar, M. B. Wise, Camb. Monogr. Part. Phys. Nucl. Phys. Cosmol. 10<br />

(2000) 1-191.<br />

[9] Determination <strong>of</strong> V cb and V ub , in C. Amsler et al., Phys. Lett. B667, 1 (2008),<br />

Available online at:<br />

http://pdg.lbl.gov/2010/reviews/rpp2010-rev-vcb-vub.pdf<br />

[10] M. Neubert, [hep-ph/0512222].<br />

54

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!