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 ...
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