Introduction to Quantum Chromodynamics

Similar documents
1 Running and matching of the QED coupling constant

Theory of Elementary Particles homework VIII (June 04)

Ultraviolet Divergences

5 Infrared Divergences

REVIEW REVIEW. Quantum Field Theory II

Quantum Field Theory II

Loop corrections in Yukawa theory based on S-51

Light-Cone Quantization of Electrodynamics

Unitarity, Dispersion Relations, Cutkosky s Cutting Rules

Review of scalar field theory. Srednicki 5, 9, 10

A Renormalization Group Primer

Introduction to Quantum Chromodynamics (QCD)

Finite Temperature Field Theory

Lecture notes Particle Physics II. Quantum Chromo Dynamics. 6. Asymptotic Freedom. Michiel Botje Nikhef, Science Park, Amsterdam

PHY 396 L. Solutions for homework set #20.

QFT PS7: Interacting Quantum Field Theory: λφ 4 (30/11/17) The full propagator in λφ 4 theory. Consider a theory of a real scalar field φ

8.324 Quantum Field Theory II. Problem Set 7 Solutions. d λ dλ

Review and Preview. We are testing QED beyond the leading order of perturbation theory. We encounter... IR divergences from soft photons;

Summary of free theory: one particle state: vacuum state is annihilated by all a s: then, one particle state has normalization:

Regularization Physics 230A, Spring 2007, Hitoshi Murayama

Vacuum Energy and Effective Potentials

DR.RUPNATHJI( DR.RUPAK NATH )

6. QED. Particle and Nuclear Physics. Dr. Tina Potter. Dr. Tina Potter 6. QED 1

Introduction to particle physics Lecture 6

Physics 444: Quantum Field Theory 2. Homework 2.

The 1/N expansion method in quantum field theory

PHYS 611, SPR 12: HOMEWORK NO. 4 Solution Guide

QCD and Rescattering in Nuclear Targets Lecture 2

Manifestly diffeomorphism invariant classical Exact Renormalization Group

Renormalization of the Yukawa Theory Physics 230A (Spring 2007), Hitoshi Murayama

Towards a manifestly diffeomorphism invariant Exact Renormalization Group

using D 2 D 2 D 2 = 16p 2 D 2

Quantum Field Theory 2 nd Edition

Loop Corrections: Radiative Corrections, Renormalization and All

21 Renormalization group

A Modern Anatomy of Electron Mass. Abstract

6.1 Quadratic Casimir Invariants

221B Lecture Notes Scattering Theory II

The Non-commutative S matrix

Introduction to the physics of highly charged ions. Lecture 12: Self-energy and vertex correction

2P + E = 3V 3 + 4V 4 (S.2) D = 4 E

NTNU Trondheim, Institutt for fysikk

Introduction to Elementary Particle Physics I

Effective Field Theory for Nuclear Physics! Akshay Vaghani! Mississippi State University!

Quantum Field Theory. and the Standard Model. !H Cambridge UNIVERSITY PRESS MATTHEW D. SCHWARTZ. Harvard University

10 Thermal field theory

Lecture 3 (Part 1) Physics 4213/5213

Effective Field Theory for Cold Atoms I

Scattering Amplitudes

Units, dimensions and regularization

Feynman Rules of Non-Abelian Gauge Theory

Overview of low energy NN interaction and few nucleon systems

Quantum Field Theory II

SM, EWSB & Higgs. MITP Summer School 2017 Joint Challenges for Cosmology and Colliders. Homework & Exercises

Kern- und Teilchenphysik II Lecture 1: QCD

Evaluation of Triangle Diagrams

The Strong Interaction and LHC phenomenology

Lecture 10. September 28, 2017

Maxwell s equations. electric field charge density. current density

Quantum Field Theory Spring 2019 Problem sheet 3 (Part I)

Finite-temperature Field Theory

Donoghue, Golowich, Holstein Chapter 4, 6

Lecture 6:Feynman diagrams and QED

Particle Physics I Lecture Exam Question Sheet

Landau s Fermi Liquid Theory

Week 3: Renormalizable lagrangians and the Standard model lagrangian 1 Reading material from the books

2. HEAVY QUARK PRODUCTION

PAPER 51 ADVANCED QUANTUM FIELD THEORY

4 4 and perturbation theory

Effective Field Theory

Removing Infrared Divergences

Renormalization Scheme Dependence

QCD β Function. ǫ C. multiplet

3 Quantization of the Dirac equation

Introduction to perturbative QCD and factorization

1 Introduction. 1.1 The Standard Model of particle physics The fundamental particles

arxiv:hep-ph/ v1 12 Sep 1996

On the renormalization-scheme dependence in quantum field theory.

Generalized Parton Distributions in PT

arxiv:hep-th/ v1 20 Jul 2004

Electron-positron production in kinematic conditions of PrimEx

A Brief Introduction to AdS/CFT Correspondence

Konishi Anomaly. 32π 2 ǫκλµν F κλ F µν. (5)

Richard Williams C. S. Fischer, W. Heupel, H. Sanchis-Alepuz

Derivation of Electro Weak Unification and Final Form of Standard Model with QCD and Gluons 1 W W W 3

lattice QCD and the hadron spectrum Jozef Dudek ODU/JLab

Lecture 5 The Renormalization Group

Quantum Fields in Curved Spacetime

QFT Dimensional Analysis

8.821 F2008 Lecture 12: Boundary of AdS; Poincaré patch; wave equation in AdS

QFT. Chapter 14: Loop Corrections to the Propagator

lattice QCD and the hadron spectrum Jozef Dudek ODU/JLab

Lecture notes for QFT I (662)

LOOP-TREE DUALITY AND QUANTUM FIELD THEORY IN 4D

752 Final. April 16, Fadeev Popov Ghosts and Non-Abelian Gauge Fields. Tim Wendler BYU Physics and Astronomy. The standard model Lagrangian

The boundary state from open string fields. Yuji Okawa University of Tokyo, Komaba. March 9, 2009 at Nagoya

Zhong-Zhi Xianyu (CMSA Harvard) Tsinghua June 30, 2016

d 1 µ 2 Θ = 0. (4.1) consider first the case of m = 0 where there is no azimuthal dependence on the angle φ.

TWISTOR DIAGRAMS for Yang-Mills scattering amplitudes

Zhong-Bo Kang Los Alamos National Laboratory

Transcription:

Introduction to Quantum Chromodynamics Michal Šumbera Nuclear Physics Institute ASCR, Prague November 26, 2009 Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 1 / 78

Ultraviolet renormalization basic ideas and techniques 1 Introduction 2 Renormalization of φ 4 theory 2 Elements of dimensional regularization Dimensional regularization of φ 4 theory 3 Renormalization in QED Electric charge renormalization Pauli Villars Renormalization group Dimensional transmutation Technique of counterterms Feynman parametrization, momentum integrals and related formulae Dimensional regularization in QED and QCD 3 Renormalization in QCD 4 Exercises Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 2 / 78

Literature Our discussion is based on Quarks, partons and Quantum Chromodynamics by Jiří Chýla Available at http://www-hep.fzu.cz/ chyla/lectures/text.pdf Additional material comes from Quantum Field Theory A Modern Introduction by Michio Kaku Oxford University Press 1993 Additional material comes from Quantum Field Theory in a Nutshell by A. Zee Princenton University Press 2003 Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 3 / 78

Vacuum as a dielectric In QED: any charged particle is surrounded by a dense cloud of e e + virtual pairs that tend to screen the charge of a particle. At large distances effective coupling constant is reduced by presence of this screening charge. At smaller distances a probe can penetrate through this virtual cloud, and hence the QED coupling constant gets larger as we increase the energy of the probe. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 4 / 78

Vacuum as a dielectric Classically, we can think of this in terms of the dielectric constant of the vacuum. Place a charge into a dielectric. Electric field of the dielectric will cause the dipoles within the dielectric medium to line up around the charge and decrease its value. the medium will acquire dielectric constant greater than one. QCD: we have color charges and color coupling constants. Each gluon carries both a color charge and an anti-color magnetic moment. The net effect of polarization of virtual gluons in the vacuum is not to screen the field, but to increase it and affect its color. This is sometimes called antiscreening. In QCD: At large distances (low energy) presence of the cloud of virtual particles creates an antiscreening effect and the net coupling constant gets larger. At smaller distances a probe that comes near a colored particle feels the coupling constant decrease at high energies. Thus, the dielectric constant of the vacuum is less than one for an asymptotically free theory. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 5 / 78

Infinities in perturbative calculations Transition from a finite number of degrees of freedom in QM to an infinite number in QFT we must continually sum over an infinite number of internal modes in loop integrations, leading to divergences. Divergent nature of QFT reflects the fact that the UV region is sensitive to the infinite number of degrees of freedom of the theory. We ll discus UV renormalization of perturbative φ 4, QED and QCD. Renormalization procedure is not merely a technique of removing unpleasant UV infinities, but first of all an effective description of quantum phenomena. Essence of UV renormalization is common to QED and QCD we start with QED. Here the basic quantity governing the strength of interactions of charged particles are their electric charges. It is common to use fine structure constant α e 2 /4π = 1/137 (classically). L QED leads to the pqft, in which physical quantities are expressed as power expansions in α. It turns out that at higher orders the coefficients of these expansions, calculated according to the standard Feynman rules, come out formally infinite. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 6 / 78

Divergences in integration over the loop momenta A: Ultraviolet divergences (UV), coming from integration over the large values of the loop momenta. B: Mass divergences, coming from integration over the region of small virtualities. These small virtualities appear in two different situations: For vanishingly small energy and momentum of the virtual particles. These so called infrared (IR) singularities occur for the massless photon, independently of the electron mass. When two of the three particles in the QED vertex e-γ-e are parallel to each other. This can happen only if m e = m γ = 0. parallel singularities are thus absent in QED. Trace of them shows up potentially large logarithms like α k ln k (Q/m e ) (where Q is some external momentum) at k th order of perturbation theory. Both types of mass singularities can be regularized by introduction of a fictitious mass of the photon. Here we shall analyze UV divergences which lead to the important concept of the renormalized electric charge. The mass singularities are dealt with in a completely different way and are discussed in the text of J. Chýla. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 7 / 78

UV example: φ 4 theory Consider second order diagrams in φ 4. In addition to standard tree diagram (Fig.1, left) there is also so called loop diagram (Fig.1, right). Compare their propagators associated with internal line(s): ( iλ) 2 i q 2 m 2 + iε = i ( iλ)2 (k 4 + k 5 + k 6 ) 2 m 2 + iε 1 2 ( iλ)2 d 4 k (2π) 4 i i k 2 m 2 + iε (k 1 + k 2 k) 2 m 2 + iε (1) d 4 k k 4 (2) Figure 1: Feynman diagrams occuring in order g 2 in perturbation theory. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 8 / 78

Regularization: φ 4 theory Tree diagram (Fig.1, left): internal line is associated with a virtual particle with q 2 = (k 4 + k 5 + k 6 ) 2 m 2. The farther the momentum of the virtual particle is from the mass shell the smaller is the amplitude. Virtual particle is penalized for not being real! Loop diagram integrand is large only if one or the other or both of the virtual particles associated with the two internal lines are close to being real. Once again, there is a penalty for not being real. However for very large virtualities (k ) the integral in (2) diverges logarithmically! Thus, in evaluating d 4 k/(2π) 4 in (2) we should integrate only up to Λ, known as a cutoff. The integral is said to have been regularized. Result is 2iC ln Λ2 K where K k 2 1 + k 2 and C is some numerical constant. Modern view of divergences problem: given model should be regarded as an effective low energy theory, valid up to some energy scale Λ. Going to higher and higher momentum scales may eventually bring us to completely new physics region where our simple φ 4 model is not anymore true! Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 9 / 78

(*) From now on [A] denotes dimension of A in mass units. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 10 / 78 Renormalization and dimensional analysis Since = c = 1 S = d 4 xl is dimensionless: [S] = 0 [L] = 4. In this notation [x] = 1 and [ ] = 1. [ In scalar φ 4 theory L = 1 2 ( φ) 2 m 2 φ 2] λφ 4. For term ( φ) 2 to have the same dimension as L i.e. [ ( φ) 2] = 4 we must have [φ] = 1 implying that [λ] = 0, i.e. coupling is dimensionless. For fermion field ψ we find its dimension from the free field lagrangian L = ψiγ µ µ ψ +... to be [ψ] = 3 2. Looking at coupling f φψψ Yukawa coupling f is also dimensionless. From the Maxwell Lagrangian L = F µν F µν we see that [F µν ] = [F µν ] = 2 [A µ ] = 1 i.e. dimensions of [A µ ] and [ µ ] are the same. Consider QED interaction term ea µ ψγ µ ψ: [A µ ψγ µ ψ] = 4 [e] = 0. N.B. The same conclusion can be also deduce from Coulomb s law written in natural units V (r) = α/r, with the fine structure constant α = e 2 /4π.

Fermi theory of the weak interaction: M if blows up In contrast 4-fermion interaction Lagrangian L = Gψψψψ has [G] = 2. Recall dσ = (2π)4 1 1 n M if 2 d p i n v 1 v 2 2E 1 2E 2 (2π) 3 δ 4 (p 1 + p 2 p i ) SF 2E i i=3 i=3 [σ] = 2, [v] = 0, [E] = 1 [M if ] = 0. Calculate ν + ν ν + ν scattering at s Λ. Since all masses and energies are by definition small compared to the cutoff Λ, we can simply set them equal to zero. In the lowest order: M if G and is finite. In the next order: M if G + G 2 Λ 2 [M if ] = 0 or by looking directly at the Feynman diagram (left) which goes as G 2 Λ d 4 p(1/p)(1/p) G 2 Λ 2. For Λ the theory is sick: M if 2. Infinity is predicted value for a physical quantity! Fermi theory and all other non-renormalizable theories have the ability to announce its own eventual failure and hence their domains of validity. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 11 / 78

Renormalization UV renormalization procedure consists of two distinct steps: 1) Regularization: Formally divergent integrals are regularized, i.e. basically cut off at high values of the loop momenta. There are numerous regularization techniques, for instance: Simple cut-off. Analytical regularization. Pauli Villars and its modifications: gives a clear physical interpretation of the renormalized electric charge while maintaining gauge invariance. In multiloop calculations this technique is, however, rather cumbersome. Dimensional regularization: is much simpler and in nonabelian gauge theories, like QCD, it preserves gauge invariance. The renormalization scale, on which the renormalized electric charge does depend, enters in a rather artificial way and its physical interpretation is therefore rather obscure. 2) Renormalization: Divergent terms are absorbed in the newly defined renormalized quantities (couplant, mass and wave functions). Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 12 / 78

Example: Potential of infinite charged line Consider infinitely long wire carrying constant charge density ρ. Potential at at point P at distance R from the wire: ρ(r) + dx V (R) = dx = ρ (3) r R2 + x 2 is logarithmically divergent potential is ill defined. Contrary to this electric field: Regularizing (3) we obtain: E = V = V (R) ρ V Λ (R) = ρ +Λ Λ + leading to slight redefinition of electric field as: [ E = lim ] V Λ (R) = lim Λ Λ dx <. (4) (R 2 + x 2 ) 3/2 [ ] dx R2 + x = ρ ln Λ2 + R 2 + Λ 2 Λ2 + R 2 Λ ˆR 2ρ R (5) Λ Λ2 + R 2ρ 2 R ˆR (6) Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 13 / 78

Example: Potential of infinite charged line We had to introduce new variable Λ with dimension of length i.e. [Λ] = 1. Notice that the difference: = lim [V Λ(r 2 ) V Λ (r 1 )] = ρ ln r 1 2 Λ r2 2 is well defined. So we can use it to renormalize potential by subtracting V (R) at some fixed value of R = R 0 and taking the limit Λ : (7) V (R) V (R) V (R 0 ) = ρ ln R2 0 R 2 (8) Introduction of dimensionful parameter R 0 has caused non-physical infinities present in V (R) and V (R 0 ) to cancel each other, leaving a finite result with a non-trivial R-dependence. The cutoff Λ has disappeared! This example suggest a strategy for dealing with divergencies: 1 Identify an appropriate way to regularize infinite integrals. 2 Absorb the divergent terms into a redefinition of fields or parameters e.g. via substractions. This step is usually called renormalization. 3 Make sure the procedure is consistent, by checking that the physical results do not depend on the regularization prescription. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 14 / 78

Renormalization of φ 4 theory Divergent structure of any graph can be analyzed in terms of: E = number of external legs I = number of internal lines V = number of vertices L = number of loops Degree of its divergence D comes from power counting: each internal propagator contributes 1/p 2, each loop contributes d 4 p D = 4L 2I. Vertex has 4 lines connecting to it. Each of these lines, in turn, either ends on an external leg, or on one end of an internal leg, which has 2 ends. 4E = 2I + E. Expression for the loop number L = I V + 1 can be obtained as follows: # of independent momenta = # of internal lines I - the constraints coming from momentum conservation. There are V such momentum constraints, minus the overall momentum conservation from the entire graph. # of independent momenta in a Feynman graph is also equal to the # of loop momenta. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 15 / 78

Renormalization of φ 4 theory Inserting these graphical rules into our expression for D we have: D = 4 E. Degree of divergence of any graph in 4 dimensions depends only on the number of external lines, which is a necessary condition for renormalizability. D is thus independent on number of internal loops in the graph 1. N.B. Situation in QED where D = 4 3 2 E ψ E A and E ψ,a is the number of external electron and photon legs, respectively, is very similar. In φ 4 theory only two-point and four-point graphs are divergent. We need to renormalize only two physical quantities: the mass and the coupling constant. Thus, by using only power counting arguments, in principle we can renormalize the entire theory with only two redefinitions corresponding to two physical parameters. (1) In d dimensions D = d + (1 d/2)e + (d 4)V. Because D increases with # of internal vertices, there are problems with renormalizing the theory in higher dimensions. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 16 / 78

Renormalization in a nuttshell: What is actually measured Introduce cutt-off Λ into integral (2) and assume m 2 (k 1 + k 2 ) 2 (so that we can neglect m 2 in the integrand). The scattering amplitude up to the second order in λ reads: M if = iλ + icλ 2 L(s, t, u) + O(λ 3 ) (9) where [ ( ) ( ) ( )] Λ 2 Λ 2 Λ 2 L(s, t, u) ln + ln + ln s t u (10) and s K 2 = (k 1 + k 2 ) 2, t (k 1 k 3 ) 2 and u (k 1 k 4 ) 2. M if is supposed to be an actual (physical) scattering amplitude. M if should not depend on Λ. Any change of Λ must be always compensated by the shift of λ in such a way so that M if does not change. Define physical coupling as: iλ P = iλ + icλ 2 L(s 0, t 0, u 0 ) + O(λ 3 ) (11) and solve it for λ: iλ = iλ P icλ 2 L(s 0, t 0, u 0 )+O(λ 3 ) = iλ P icλ 2 PL(s 0, t 0, u 0 )+O(λ 3 P) (12) where the second equality is allowed to the order of approximation indicated. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 17 / 78

Renormalization in a nuttshell: What is actually measured Now plug this into (9): M if = iλ + icλ 2 L(s, t, u) + O(λ 3 ) = = iλ P icλ 2 PL(s 0, t 0, u 0 ) + icλ 2 PL(s, t, u) + O(λ 3 P) (13) where, once again, all manipulations are legitimate up to the order of approximation indicated. Now in the scattering amplitude M if we have the combination L(s, t, u) L(s 0, t 0, u 0 ) = [ln(s 0 /s) + ln(t 0 /t) + ln(u 0 /u)] and so the scattering amplitude comes out as: [ ( M if = iλ P + icλ 2 s0 ) ( t0 ) ( u0 )] P ln + ln + ln + O(λ 3 s t u P) (14) i.e. likewise (8) M if is now expressed in terms of measurable quantity in case of (14) it is the physical coupling constant λ P. The unmeasurable cutoff Λ has completely disappeared! The lesson: express physical quantities not in terms of fictitious theoretical quantities such as λ, but in terms of measurable quantities such as λ P. Later is in the literature often denoted by λ R and for historical reasons called the renormalized coupling constant. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 18 / 78

Elements of dimensional regularization Consider again φ 4 theory in n=4 dimensional Minkowski space time. L φ 1 2 ( µφ(x)) 2 m 2 φ 2 λ 4! φ4 (x), λ > 0. (15) Simplest one loop diagram Fig. 2b describes order λ 2 correction to the leading order diagram, depicted in Fig. 2a. Since basic four particle vertex is associated with iλ the amplitude corresponding to digram in Fig. 2b is: I (p, n = 4) λ2 2 d 4 k 1 (2π) 4 [k 2 m 2 + iɛ] 1 [(k p) 2 m 2 + iɛ], (16) where p p 1 + p 2. (16) is the same as (2) but in the other variables. Figure 2: Lowest order Feynman diagrams for the elastic scattering of two scalar particles with mass m. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 19 / 78

Dimensional regularization of φ 4 theory I (p, n) = Basic idea of DR: in n dimensions the integral in (16) d n k k n 1. While in 4d it is logarithmically divergent, it converges in n = 1, 2, 3 dimensions. In n = 1, 2, 3 dimensions and in the rest frame of p = (M, 0, 0, 0) (16) equals λ2 2(2π) 4 dk 0 d n 1 1 k [ k0 2 ] k 2 m 2 + iɛ 1 [ (k 0 M) 2 k 2 m 2 + iɛ The integrand depends now merely on k of (n 1) dimensional vector k. Integral over d n 1 k can therefore be performed in polar coordinates: ]. (17) d n k = r n 1 dr sin n 2 θ n 1 sin n 3 θ n 2 sin θ 2 dθ 2 dθ 1, (18) where r k 2 and 0 θ 1 2π is the azimuthal angle, all other angles θ i (0, π) are generalizations of well known polar angles in 3 dimensions. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 20 / 78

Dimensional regularization of φ 4 theory: n < 4 Applying (18) to (17) and denoting ω k2 we get: I (p, n) = λ 2 2(2π) 4 dk 0 dωω n 2 1 1 0 [k0 2 ω2 m 2 + iɛ] [(k 0 M) 2 ω 2 m 2 + iɛ] W (n), (19) where W (n) contains the integral over the angular coordinates W (n) 2π 0 π π dθ 1 dθ 2 sin θ 2 dθ 3 sin 2 θ 3 0 0 π 0 dθ n 2 sin n 3 θ n 2. (20) Angular part of f (x) depending only on r = x 2 of the n dimensional vector x, can be carried out using the standard formulae π ( πγ m+1 ) dθ sin m 2 θ = 0 Γ ( ) m+2 d n xf (r) = 2πn/2 Γ( n 2 2 ) drr n 1 f (r) (21) with the result: I (p, n) = λ 2 2π (n 1)/2 2(2π) 4 Γ( n 1 2 ) dk 0 dωω n 2 1 1 0 [k0 2 ω2 m 2 + iɛ] [(k 0 M) 2 ω 2 m 2 + iɛ]. (22) Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 21 / 78

Dimensional regularization of φ 4 theory: n < 4 For n = 4 (22) reduces to (16), but it has also a well defined meaning for all real n (1, 4)! With I (p, n) defined on an open interval we can continue it analytically to the whole complex n plane. The singularity at n = 4 is a simple pole! To continue I (p, n) below n = 1 we rewrite (22) using l times a simple per partes integration and taking into account that / ω = 2ω / ω 2 I (p, n) = λ2 2π (n 1)/2 2(2π) 4 Γ( n 1 2 + l) ( dk 0 dωω n 2+2l ) l 1 1 0 ω 2 [k0 2 ω2 m 2 + iɛ] [(k 0 M) 2 ω 2 m 2 + iɛ], which converges for all integer 1 2l < n < 4. (23) Rewrite I (p, n) so that it coincides with original form (22) in the region n (1, 4) but has also a well defined meaning in a larger interval. Per partes integration does not change the behaviour of the integrand at large ω, but improves the convergence property at the origin ω = 0, thereby extending the convergence of (23) to lower values of n. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 22 / 78

Dimensional regularization of φ 4 theory: n > 4 To continue analytically above n = 4 we insert into (22) the identity 1 = 1 ( dk0 + dω ). (24) 2 dk 0 dω Performing the integrals over k 0 and ω again using per partes we arrive at: I (p, n) = 1 2 ( n + 6)I (p, n) + I (p, n) I (p, n) = 2 n 4 I (p, n), (25) where I (p, n) = λ2 2π (n 1)/2 2(2π) 4 dk 0 dωω n 2 u(k 0, ω, M), (26) u(k 0, ω, M) Γ( n 1 2 ) 2m 2 (k 2 0 m2 ω 2 + iɛ) 2 [(k 0 M) 2 ω 2 m 2 + iɛ] + 2(k 0 M)M + 2m 2 (k 2 0 m2 ω 2 + iɛ) [(k 0 M) 2 ω 2 m 2 + iɛ] 2. (27) N.B. I (p, n) exists for all 1 < n < 5 I (p, n) in (25), as a function of complex n, has a simple pole at n = 4! Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 23 / 78

Dimensional regularization of φ 4 theory Consider the diagram in Fig. 2c, where the bubble appears in the t channel. First we combine the two propagators in (16) into one, using the identity 1 ab = 1 0 1 dx [ax + b(1 x)] 2, (28) which is a special case of Feynman parametrization. For (16) it gives I (p, n) = λ2 d n k 1 dx 2 (2π) n 0 [(k 2 m 2 )x + ((k p) 2 m 2 )(1 x)] 2 1 = λ2 d n l 1 dx 2 0 (2π) n [l 2 a 2 (x) iɛ] 2, (29) where we have changed the order of integrations and subsituted l k p(1 x), a 2 (x) m 2 p 2 x(1 x), p 2 (p 1 p 3 ) 2. (30) The inner integral in (29) yields: I (p, n) = λ2 4 i (4π) 2 (4π)ε Γ(ε) 1 0 dx [ a 2 (x) ] ε, ε 2 n 2. (31) Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 24 / 78

Dimensional regularization of φ 4 theory: mass parameter µ In n dimensions the coupling λ is no longer dimensionless, as in 4d, but has dimension [λ] = 4 n = 2ε. The same dimension has also I (p, n). For n close to the physical value n = 4 in (31) we can use: Γ(ε) = 1 ε Γ(1 + ε) = 1 ε γ ( ε) n 1 E + ζ(n), (32) n! (4π) ε = 1 + ε ln 4π +, (33) n=2 where ζ(n) is the Riemann zeta function and γ E = 0.5772 We cannot, however, expand in a similar way directly a 2ε because [a] = 2. Physically well defined expansion requires introducing some mass parameter µ that provides a scale for the expansion of the logarithm ( a 2 ) ε (x) µ 2 = 1 ε ln a2 (x) µ 2 +. (34) Inserting identity 1 = µ 2ε µ 2ε into (31) we get [ I (p, n) = i λ2 1 1 2 (4π) 2 µ 2ε ε γ E + ln 4π 1 0 ] dx ln a2 (x) µ 2 + O(ε). (35) Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 25 / 78

Coupling renormalization in φ 4 theory: λ B and λ R New parameter µ is an integral part of DR. I (p, n) does not depend on it, provided we keep all the terms in expansion (34). A truncated version of (31), like the first two terms in the Laurent series in (35), which is the only part of (34) that survives in the limit ε 0, still contains µ despite the absence of the factor µ 2ε to cancel this dependence. This apparent inconsistency disappears in the process of defining the µ-dependent renormalized coupling. Summing the contributions of diagrams in Fig. 2a,b we get : [ iλ B 1 λ B µ 2ε ( 1 1 2 (4π) 2 ε γ E + ln 4π 0 dx ln a2 (x) µ 2 which to the order considered can be rewritten as [ i λ B 1 λ B µ 2ε ( )] [ 1 2 (4π) 2 ε γ E + ln 4π 1 + λ B µ 2ε 1 ] 2 (4π) }{{} 2 dx ln a2 (x) 0 µ }{{ 2 } µ 2ε λ R (µ) finite terms (*) From now on λ B λ to indicate its later interpretation as bare coupling. )], (36) (37) Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 26 / 78

Coupling renormalization in φ 4 theory: λ B and λ R Renormalized coupling λ R (µ) is defined in terms of the bare coupling λ B : [ λ R (µ) µ 2ε µ 2ε ( )] 1 λ B 1 λ B 32π 2 ε γ E + ln 4π. (38) While [λ B ] = 2ε and is held fixed when µ is varied, [λ R (µ)] = 0 and absorbs the singular term 1/ε, eventually plus some finite terms, like γ E + ln 4π in (38). The convention defined in (38) is called MS, that corresponding to vanishing finite terms is called MS (minimal subtraction). Evaluating derivative of λ R (µ) we find a finite result as ε 0: dλ R (µ) d ln µ = µdλ R(µ) dµ = 2ελ B µ 2ε + 1 8π 2 λ2 Bµ 4ε + O(λ 3 B) (39) = 2ελ R + 1 16π 2 λ2 R + O(λ 3 R) β φ (λ R, ε) N.B. In the limit ε 0 term proportional to λ B vanishes but must be kept if the renormalization procedure is carried out in terms of renormalized quantities and so called counterterms. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 27 / 78

Coupling renormalization in φ 4 theory The second, finite, term comes from the product λ 2 B 1 1 32π 2 ε d ln µ = 1 1 λ2 B 32π 2 ε ( 4ε)µ 4ε = 1 8π 2 λ2 Bµ 4ε = λ2 R 8π 2. (40) dµ 4ε where replacement λ 2 B λ2 R in (40) is legal as these two expressions start to differ first at order λ 3 B. (40) determines implicit dependence of λ R (µ) on µ, which cancels explicit dependence on µ of the finite terms in (37). This cancellation holds, however, only for physical quantities, like the cross sections and only if perturbation theory is summed to all orders. Leading order β function coefficient 1/16π 2, obtained within DR in (40), is actually independent of the regularization method used in its derivation! For instance, in the conventional Pauli Villars type of regularization we have instead of (38) λ R (µ) λ B [ 1 λ B 1 32π 2 (ln M2 + finite terms µ 2 where M is large cut off to be sent to infinity. In this case 1/16π 2 comes from derivative with respect to ln µ of logarithmic term ln(m/µ). )] (41) Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 28 / 78

Electric charge renormalization in QED e + µ scattering: in the lowest order process is described by diagram in Fig. 3a, while the order O(α 2 ) corrections correspond to diagrams in Figs. 3b-e. While diagrams on Fig. 3a,e give finite results the other are UV divergent. We ll work in the Feynman covariant gauge, i.e. set α G = 1 in d µν = g µν + (1 α G ) kµkν k 2. Quantity standing in each vertex on Fig. 3: e B bare electric charge. Figure 3: Feynman diagrams for electron muon scattering in a leading (a) and next to leading order (b-d). The loop in c) can equally well be on any other of the four fermion legs while that in d) can circumvent the lower vertex as well. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 29 / 78

Renormalization in QED Electric charge renormalization Leading order (LO) diagram Fig. 3a gives a contribution to the invariant amplitude M if (see Fig. 3 for notation) ( [u(p )( iγ µ )u(p)] [u(k )( iγ ν )u(k)] ig ) µν }{{} q 2 eb, 2 (42) }{{} W µν (k, p, q) D µν (q) where we have separated out the square of the bare electric charge and D µν. Contribution of Fig. 3b has a similar structure: ( ) ( ) i i W µν (k, p, q) q 2 I µν (q, m B ) q 2 eb, 2 (43) where the tensor: d 4 [ ] k I µν (q, m B ) ( 1) (2π) 4 Tr i ( ie B γ µ ) k m B + iɛ ( ie i Bγ ν ) k q m B + iɛ (44) and factor ( 1) comes from Pauli principle in the closed electron loop. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 30 / 78

Renormalization in QED Electric charge renormalization m B bare electron mass, will be throughout this subsection treated as a fixed number, although it also must be renormalized. Structure of I µν is dictated by the gauge invariance: I µν (q, m B ) = ( g µν q 2 + q µ q ν )I (q 2, m B ) (45) with I (q 2, m B ) containing the divergent integral. Naive counting of the powers of k in the integrand of (44) suggests quadratic divergence of the integral, but because of (45) the integrand does not behave like 1/k 2, where k is the loop momentum, but rather like q 2 /k 4, thereby causing only the logarithmic divergence of I (q 2, m B ). N.B. In the following we shall drop the term in (45) proportional to q µ q ν as it vanishes after contraction with the tensor W µν in (43). Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 31 / 78

Renormalization in QED Electric charge renormalization Regularization of diagram in Fig. 3b uses 4 technical ingredients: 1 Schwinger parameterization i i( k + m) = k m + iɛ k 2 m 2 = ( k + m) + iɛ dze iz(k2 m 2 +iɛ) 0 for both fermion propagators forming the closed loop adds two more integration variables z 1, z 2. 2 Changing the original loop variable k to (46) l k qz 2 z 1 + z 2 = k q + qz 1 z 1 + z 2 (47) and using the following identities to perform integrals over d 4 l [ 1, 0, d 4 l (2π) 4 [1, l µ, l µ l ν ] e il 2 (z 1+z 2+iɛ) = 1 16π 2 i(z 1 + z 2 ) 2 ig µν 2(z 1 + z 2 ) ]. (48) Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 32 / 78

Renormalization in QED Electric charge renormalization ( dt 3 Using the identity 1 = 0 t δ 1 z ) 1 + z 2 t to trade one of the integrations over z 1 or z 2 for that over t. (49) 4 The identity 0 dx x ( e iax e ibx) = ln b a. (50) Performing the first two steps we get: [ dz dz 2 1 (z 1+z 2) exp i 2 I µν (q) = ie2 b 4π 2 0 0 ( q 2 z1z2 z 1+z 2 (m 2 b iɛ)(z 1 + z 2 ) where the tensor J µν has two terms: [ ] J µν (q) 2(g µν q 2 z 1 z 2 i q µ q ν ) (z 1 + z 2 ) 2 + g µν q2 z 1 z 2 z 1 + z 2 (z 1 + z 2 ) 2 + m2 B )] J µν (q), First term satisfies the transversality condition q µ I µν (q, m B ) = 0 that follows from gauge invariance, while the second term does not. Closer examination of the integral standing by g µν shows that it actually vanishes! (51) (52) Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 33 / 78

QED Electric charge renormalization: Pauli Villars Using in the transverse term of (52), the identity (49) and carrying out the integral over z 1 or z 2, we get the following (still UV divergent!) expression: I (q 2, m B ) = ie2 B 2π 2 1 0 dz z(1 z) 0 dt t exp [ it ( q 2 z(1 z) m 2 B + iɛ )], (53) where the original logarithmic singularity of the integral over d 4 k has been transformed into the same type of logarithmic singularity over the variable t. Pauli Villars technique for regularization of integrals (53) consists in replacement I (q 2, m B ) I (q 2, m B, M) I (q 2, m B ) I (q 2, M), (54) where the subtraction is understood to be carried using the identity (50). i.e. first subtracting the integrands and then performing the integration. Mass M acts as the UV regulator and at the end of the renormalization procedure should be sent to infinity. Usage of fermion mass M in the subtracted term I (q 2, M) guarantees that this regularization technique preserves gauge invariance, which is of crucial importance for the proof of the full renormalizability of the theory. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 34 / 78

Electric charge renormalization: Pauli Villars Using (50) the regularized expression (54) equals I (q 2, m B, M) = ie2 B 2π 2 1 0 ( dzz(1 z) ln M 2 m 2 B q2 z(1 z) Isolating the logarithmically divergent term this can be rewritten as ii (q 2, m B, M) = e2 B M2 ln 12π2 µ 2 e2 B 2π 2 1 0 µ 2 ). (55) ( m 2 dzz(1 z) ln B q 2 ) z(1 z). (56) N.B. Integral over dimensionless parameter z is the trace of the original integration over the loop momentum k. Note that for any nonzero q 2 expression (55) is regular for m B 0 and can thus be used even for massless fermions. In (56) new dimensional scale µ was introduced, which scales the cut off parameter M in the singular term ln(m 2 /µ 2 ), but obviously I (q 2, m B, M) is in fact independent of µ! Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 35 / 78

Renormalization in QED Electric charge renormalization Once the regularization method has been chosen we can add contributions of diagrams in Fig. 3a,b obtaining (dropping the tensor W µν ) ig [ µν q 2 e2 B 1 e2 B M2 ln 12π2 µ 2 + e2 B 2π 2 1 Rewrite now (57) as a product of two terms Z 1 0 ( m 2 dz z(1 z) ln B q 2 )] z(1 z) µ 2 (57) 3 (A, eb, 2 M, µ) { ig [ }} ]{ [ 1 ( µν q 2 e2 B 1 e2 B M2 ln 12π2 µ }{{ 2 1+ e2 B m 2 2π } 2 dzz(1 z) ln B q 2 ) z(1 z) ] 0 µ 2, }{{} er(a, 2 µ) C(A, q 2, µ, m B ) (58) where UV divergence was absorbed into into newly defined renormalized electric charge er 2 (A, µ) with the rest of the original O(eB 2 ) term left in the finite contribution C(A, q 2, µ, m B ). (*) Label A etc. is used to distinguish between different definitions of e R. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 36 / 78

Renormalization in QED Electric charge renormalization In terms of this new parameter (58) can be written as ig [ ] µν q 2 e2 R(A, µ) 1 + e2 B 2π 2 C(A, q2, µ, m B ) ig µν q 2 e2 R(A, µ), (59) where the second expression, differing from the the first one by replacement [ 1 + e R(A, µ) 2π 2 C(A, q 2, µ, m B ) eb 2 e2 R in the brackets, results if the renormalization procedure is carried out to order O(eB 6 ). In two regions the integral over z in (58) can be performed analytically. For large Q 2 q 2, i.e. for q 2 /mb 2 1 we find: ig [ µν q 2 e2 R(A, µ) 1 + e2 R (A, µ) 12π 2 while for q 2 /mb 2 1 we have: ig [ µν q 2 e2 R(A, µ) 1 + e2 R (A, µ) 12π 2 ln m2 B (ln q2 µ 2 + 5 3 ] )], (60) µ 2 + e2 R (A, µ) q 2 ] 60π 2. (61) m 2 B Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 37 / 78

Renormalization in QED Electric charge renormalization e 2 B 1 In (58) singularity of the loop integral in the term proportional to ln(m 2 /µ 2 ) was isolated and included in definition of renormalized electric charge e 2 R (µ). This procedure is not unique we can include in er 2 (µ) arbitrary finite terms as well e.g. we can separate integral over z in (57) into two parts dzz(1 z) ln M 2 m 2 B + µ2 z(1 z) 2π 2 } 0 {{ } 1 Z 1 3 (B, eb, 2 M, m B, µ) + e2 B 2π 2 1 and include in e 2 R (µ) e2 B Z 1 3 (denoted e 2 R This term can be written as: 1 Z 1 3 (B, eb, 2 M, m B, µ) = e2 B 12π 2 ln M2 µ 2 e2 B 2π 2 0 dzz(1 z) ln m 2 B + µ2 z(1 z) m 2 B + ( q2 )z(1 z) } {{ } C(B, q 2, µ, m B ) (62) 1 0 (B, µ)) the whole first term. ( ) m 2 dzz(1 z) ln B + z(1 z). µ 2 (63) Definition e 2 R (µ) e2 B Z 1 3 is motivated by the requirement that at µ 2 = q 2 finite correction C(B, q 2 = µ 2, µ) = 0, i.e. for µ 2 = q 2 all the effects of the electron loop correction are included in er 2 (B, µ). Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 38 / 78

Renormalization in QED Electric charge renormalization For large µ, i.e. when µ 2 mb 2, we find 1 Z 1 3 (B) e2 B 12π 2 (ln M2 µ 2 + 5 ), (64) 3 which differs from previous case only by presence of additional constant 5/3 in definition of er 2 (B, µ), while in region µ m B we find a very different behavior: 1 Z 1 3 (B) e2 B M2 ln 12π2 mb 2 e2 B 60π 2. µ2 mb 2. (65) For µ 0, er 2 (B, µ) approaches a finite value e2 R (B, 0) ( er(b, 2 µ, m B ) = eb 2 1 e2 B M2 ln 12π2 mb 2 + e2 B µ 2 ) 60π 2 mb 2 ( ) eb 2 1 e2 B M2 ln 12π2 mb 2 er(b, 2 0). (66) Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 39 / 78

Electric charge renormalization: e R (A, µ) Basic difference between definitions A and B of e R (µ) is how the effects of electron mass are taken into account. For m e = 0 they obviously coincide. For this case all the dependence on m B resides in (finite) correction C(A, q 2, µ, m B ). Note in particular the presence in (61) of the term proportional to ln(m 2 B /µ2 ), which would be arbitrarily large and therefore dangerous for µ 0. er 2 (A, µ) contains no information on m B and therefore behaves exactly in the same way as for the massless electron. We shall therefore refer to er 2 (A, µ) as mass independent definition of renormalized charge. Setting µ 2 = q 2 and approaching the limit q 2 0 the expansion parameter e 2 R (A, q2 ) 0, but due to the presence of the term ln( q 2 /m B ) the coefficient function C(A, q 2, µ 2 = q 2, m B ) and so do in fact all higher order coefficients as well. (*) For small µ e 2 R(A, µ) as defined in (58) turns negative and eventually blows up to as µ 0. However, taking into account higher order terms makes the couplant e 2 R(A, µ) indeed vanish at µ = 0. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 40 / 78

Electric charge renormalization: e R (B, µ) Situation is different as e 2 R (B, µ, m B) absorbs part of the dependence on m B and is therefore called mass dependent renormalized charge. For µ m B (apart from 5/3 in (64)) er 2 (B, µ, m B) behaves as in the case A, but below m B µ dependence changes dramatically, leading to the constant limit er 2 (B, 0). Usual practice: all the mass effects are included in er 2 (B, µ2 = q 2 ), which has a finite limit for q 2 0 and so has the coefficient function C(B, q 2, µ 2 = q 2, m B )! Nevertheless in principle both of these definitions of e R (µ) are equally legal. µ can thus be interpreted as the lower bound of the logarithmically divergent integral of the type M dk k = ln M µ. (67) µ Heisenberg uncertainty relations e R (µ) charge (original plus the induced one) inside a sphere of the radius r = 1/µ around the center of an electron. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 41 / 78

Scale independence and Renormalization group Arbitrary scale parameter µ appearing in both e 2 R (µ) and C(q2, µ, m B ) is an unphysical parameter results for physical quantities cannot depend on µ. If we change µ, other parameters, like masses and coupling constants, must also change in order to compensate for this effect. To keep physics invariant, changing the subtraction point µ must be offset by changes in the renormalized physical parameters as a function of the energy. This natural requirement is the essence of the renormalization group (RG), introduced in the early fifties by Gell-Mann, Low, Bogoljubov and others. RG expresses nothing else but internal consistency of renormalized pqft. Let R 1 represent some (unspecified) renormalization scheme. If Γ B is unrenormalized quantity and Γ R1 is same quantity renormalized by the scheme R 1, then: Γ R1 = Z(R 1 )Γ B (68) where Z(R 1 ) is renormalization constant under renormalization scheme R. Let us now choose a different renormalization scheme R 2. Since the unrenormalized quantity Γ B is independent of the renormalization scheme, then: Γ R2 = Z(R 2 )Γ B (69) Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 42 / 78

The renormalization group Relationship between these two renormalized quantities is given by: Γ R2 = Z(R 2, R 1 )Γ B where Z(R 2, R 1 ) = Z(R 1 )/Z(R 2 ) (70) Trivially, this satisfies a group multiplication law: where identity element is Z(R, R) = 1 Z(R 3, R 2 )Z(R 2, R 1 ) = Z(R 3, R 1 ) (71) For example, in φ 4 theory, we have the following relationship between unrenormalized and renormalized quantities: Γ (n) B (p i, g B, m B ) = Z n/2 φ Γ (n) R (p i, g, m, µ) (72) Since unrenormalized bare quantity is independent of µ, the derivative acting on the unrenormalized quantity must, by construction, be zero: 0 = µ ( µ Γ(n) B = µ ) ( µ Z n/2 φ Γ (n) R + Z n/2 φ µ ) µ Γ(n) R (73) Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 43 / 78

The renormalization group Using the chain rule for independent variables µ, g and m: d dµ = µ + g µ g + m µ m we can rewrite (74) as: ( where µ µ + β(g) g nγ(g) + mγ m(g) m β(g) µ g µ (74) ) Γ (n) (p i, g, m, µ) = 0 (75) γ(g) µ µ ln Z φ mγ m (g) µ m µ (76) These are the RG equations which express how renormalized vertex functions change when we make a change in scale µ. We can formally solve the expression (76) for the β function. d dµ = dg β(g) ln µ = µ 0 where µ 0 is some arbitrary reference point. g(µ) g(µ 0) dg β(g) (77) Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 44 / 78

Scale independence and Renormalization group Cancelation mechanism of µ dependence in expansion parameter er 2 (µ) and finite coefficient function C(q 2, µ, m B ) works in such a way that although the full sum of perturbation expansion is µ-independent, its approximation to order er 4 (µ) (in fact, to any finite one) is not! This observation is the first signal of the inevitable ambiguities which appear when renormalized QED and QCD are considered to finite orders. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 45 / 78

Electric charge renormalization Additional comments In any practical calculation we therefore have to choose some value of µ! If the finite correction term, containing in both definitions the logarithm ln( q 2 /µ 2 ) is to be reasonably small we should choose µ 2 q 2. Proof that renormalization procedure can be carried out systematically to all orders of perturbation theory, so that infinities at all orders can be absorbed in the redefined renormalized couplant, masses and wave functions, is rather nontrivial... Although the intermediate steps in the renormalization procedure depend on regularization technique used, the structure of the final results does not! There will always be some free scale parameter denoted µ, on which the redefined couplant as well as the new expansion coefficients will depend Once the renormalized perturbation expansions are expressed in terms of these free parameters and the mentioned invariants, they carry no trace of the regularization technique used. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 46 / 78

Electric charge renormalization Additional comments In this way the divergent term, containing the logarithm ln(m 2 /µ 2 ), is absorbed into definition of renormalized couplants er 2 (A, µ) or e2 R (B, µ) [ ] er(a, 2 µ) eb 2 1 e2 B M2 ln 12π2 µ 2 (78) and analogously for er 2 (B, µ). In the final step M should be taken in (78). For fixed bare charge eb 2 this is, however, impossible. However, since eb 2 appears always in combination with the divergent logarithm ln(m/µ) we can forget about their separate existence and consider the renormalized electric charge e R as the basic QED parameter. This well defined algorithm for evaluation of finite coefficients of perturbative expansions is clearly unsatisfactory on physical grounds. Many theorists, prominent among them Landau and his school, had addressed this problem asking following question: Can one define such a dependence eb 2 (M) of the bare charge on the cut-off M, which would compensate in (78) the term proportional to ln M, coming from the divergent integral, and allow the construction of the finite limit for the renormalized charge er 2 (µ) when M? Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 47 / 78

Electric charge renormalization Higher orders Let us now improve a little bit the approximation we have so far worked in. Instead of a single loop in Fig. 3b consider the sum of terms corresponding to the subset of all Feynman diagrams shown in Fig. 4. Figure 4: Subset of Feynman diagrams giving rise to a geometric series δ(a, x) = 0, Bubbles are connected via single photon lines and each separately yields exactly the same result as the bubble in Fig. 3b. Summing this subset of Feynman diagrams leads to a geometric series, the leading term of which is just (78). This series is easily summed with the result: α R (µ) = α B (M) [ 1 + α B (M) ( β 0 ln µ M + δ(x)) + ] = α B (M) 1 α B (M) (β 0 ln(µ/m) + δ(x)), (79) where β 0 = 2/3π and δ(x) is: δ(b, x) = 2 1 ( ) 1 dzz(1 z) ln + z(1 z) π x 2 0, x µ m B. (80) Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 48 / 78

Electric charge renormalization Higher orders For finite µ and m B 0, i.e. for x, δ(b, x) 5/3 and thus er(b, 2 µ) ebz 2 1 3 (B) = eb [1 2 + β 0 α B (ln µ2 M 2 5 )] 3 (81) reflecting fact that massless fermions contribute to the renormalized e 2 R (µ) at all scales µ with the same strength. For fixed µ and m B (80) diverges like ln m B /µ. This looks strange and unphysical as one would expect that infinitely heavy electrons should not have any influence on the theory at finite energy scales µ. However, resuming all the bubbles in Fig. 4a the resulting α R (B, µ, m B ) has the property we wish it to have: In the limit m B a heavy electron decouples from it. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 49 / 78

Electric charge renormalization The infrared limit A closely related feature of (79) i.e. of α R (µ) = concerns its behavior as µ 0. α B (M) 1 α B (M) [β 0 ln(µ/m) + δ(x)] A For the mass independent definition of the renormalized couplant (79) implies α R (A, µ) 0. B For the mass dependent definition the situation is differrent. Keeping the UV cut off parameter M fixed we find β 0 ln µ M + δ(b, x) β 0 ln m B M, (82) for m B 0 the infrared limit of α R (B, µ) is determined by the mass m B and coincides with that of massless electron at the scale µ = m B! This is a manifestation of the fact that in QED (and in general in theories without confinement of elementary fermions) the IR limit of the gauge coupling is determined essentially by the fermion masses. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 50 / 78

Dimensional transmutation Inverting (79) we have: α B (M) = α R (µ) 1 α R (µ) [β 0 ln(m/µ) + δ(x)]. (83) Thus for any finite α R the bare couplant α B (M) blows up to infinity at some finite value of M, corresponding to the pole of this expression. Note that even if we forget about α B (M) and work directly with α R (µ) (79) implies that for any dependence of α B (M) on M which has a finite limit as M the renormalized couplant α R (µ) vanishes! Inconsistency: a dimensionless quantity α(µ) is expressed as a function of a dimensional one (µ or M). Solution: µ enters in the ratio with some other dimensional parameter, which will be denoted Λ. We should thus write more correctly α R = α R (µ/λ), α B = α B (M/Λ). (84) Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 51 / 78

Dimensional transmutation Λ has nothing to do with the cut-off and represents a fundamental scale parameter, which appears in the theory entirely due to the renormalization procedure. This phenomenon is called dimensional transmutation and is typical quantum phenomenon valid for theories with dimensionless couplants. It is caused by the logarithmic divergences in one-loop diagrams which imply that this constant actually depends on the typical energy scale of the processes under considerations. The running is determined by the β-function and renormalization group. Appearance of Λ is an inevitable consequence of the renormalization procedure but its numerical value is not fixed by these considerations and must be determined from experimental data. Consequently, the strength of the interaction may be described by a dimensionful parameter the energy scale where the interaction strength reaches the value 1. In the case of QCD, this energy scale is called the QCD scale and its value 150 MeV replaces the original dimensionless coupling constant. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 52 / 78

Dimensional transmutation The standard way to define Λ unambiguously is to write down and solve the differential equations which follow from (79) and (83). A we have simply dα R (A, µ/λ) d ln µ = β 0 α 2 R(A, µ/λ) α R (A, µ/λ) = 1 β 0 ln(λ/µ), (85) dα B (M/Λ) d ln M = β 0αB(M/Λ) 2 1 α B (M/Λ) = β 0 ln(λ/m), (86) B the equation for α R (B, µ) is more complicated ( ) dα R (B, µ/λ) dδ(b, x) 1 = β 0 + αr 2 = β 0 α 2 d ln µ d ln µ R(B, µ/λ) dz 6x 2 z 2 (1 z) 2 0 1 + x 2. z(1 z) }{{} h(x) (87) Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 53 / 78

Dimensional transmutation Function h(x) in 87 is: h(x) = 1 6 x 2 + 12 x 3 4 + x 2 ln. = x 2 5 + x 2 4 + x 2 + x 4 + x 2 x This is illustrated in Fig. 5. We see that h(x) approaches unity for x 1 but vanishes like x 2 as x 0. Figure 5: Shape of the function h(x) together with its step approximation (dashed line) and the analytical approximation x 2 /(5 + x 2 ) (dotted curve) Consequently the solutions of (87) are essentially same in A and B in first region, but differ substantially in the second, where α R (A, µ/λ) 0 while α R (B, µ/λ) flattens to a constant value, (see Fig. 6). Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 54 / 78

Dimensional transmutation In both cases Λ labels the elements of an infinite set of solutions of these equations. Note that after taking the derivative of (79) with respect to ln µ the cut off M has completely disappeared from the resulting expression (85) and (87)! Figure 6: Sketch of dependence of α R (µ) on µ for both A and B definitions. Taking (85) or (87) as definitions of the renormalized couplant α R (µ/λ) is the crucial step of the renormalization procedure. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 55 / 78

Effective charge Since α R 0 it makes sense only for µ < Λ and should be trusted only for µ Λ. So even if we start from (85) i.e. α R (A, µ/λ) = [β 0 ln(λ/µ)] 1 we encounter severe problems when we approach short distances. The renormalized couplant diverges at µ = Λ and we can, knowing α R at some µ, evaluate Λ from the any of two expressions [ ] [ ] 1 1 Λ = µ exp = M exp. (88) β 0 α R (A, µ/λ) β 0 α B (M/Λ) For case B same relations hold for x 1, while for x 1 (88) are replaced by more complicated formulae. Thus we can alternatively determine Λ from the knowledge of the bare couplant α at some cut off M. Dependence of α R (µ/λ) on µ has also a very simple and intuitive physical interpretation: er 2 (µ/λ) is the charge contained inside the sphere of the radius r = 1/µ around the core of an electron. It can be viewed as the effective charge at a distance defined by the scale µ. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 56 / 78

Large distance behavior of QED: q 2 0 Use of α R (B, µ) as an expansion parameter is appropriate even in this region as it gives a finite value for the contribution of the sum of diagrams in Fig. 3a,b, equal simply to er 2 (B, 0). Including the (finite) contribution of the sum of other three diagrams in this figure leads to the final result proportional to er 2 (B, 0)(1 + r 1eR 2 (B, 0)), where r 1 is some number. Recalling the way e R (B, µ) was introduced, it is obvious that we could add to the quantity B and simultaneously subtract from C(B, q 2, µ, m B ) in (62) such finite term that the resulting r 1 actually vanishes! Denoting the corresponding renormalized charge e R (C, µ), the q 2 0 limit of the electron muon scattering amplitude, evaluated using this definition of renormalized electric charge, would be proportional to er 2 (C, 0), with no higher order corrections. Repeating this procedure for the case of Compton scattering on an electron leads to another renormalized charge: er 2 (D, µ = 0). According to the conventional definition of electric charge, this quantity coincides with the fine structure constant α cl [ er 2 (D, 0)] 2 /4π = 1 137. Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 57 / 78

Technique of counterterms In QED at one loop all UV divergencies can be absorbed in the definition of the renormalized electron charge e 2 R (µ), electron mass m R (µ), electron and photon fields A µ R (x, µ) and ψ R(x, µ), which replace their bare analogues appearing in the original QED lagrangian. In terms of these quantities all physical observables become finite. At one loop this claim follows if we add the contributions from diagrams in Fig. 3a-d, with the lower vertex amputated. In terms of the renormalization factors Z 1, Z 2F, Z 3, Z m, and writing out only the singular parts of the respective contributions we had from Fig. 3a: ie B γ µ, Fig. 3b: ie B γ µ (Z ( 3 1), ) Fig. 3c: ie B γ µ 1 Z 1 2F, ( Fig. 3d: ie B γ µ Z 1 1 1 ). Michal Šumbera (NPI ASCR, Prague) Introduction to QCD November 26, 2009 58 / 78