U NIVERSIDAD C OMPLUTENSE DE M ADRID

A PUNTES DEL CURSO

Campos Cuánticos

Antonio L. Maroto Ignazio Scimemi

D EPARTAMENTO DE F ÍSICA T EÓRICA
I NSTITUTO DE F ÍSICA DE PARTÍCULAS Y DEL C OSMOS IPARCOS-UCM

C URSO 2022-23

September 6, 2023
Contenidos

1 Introducción 4
1.1 Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4

2 Transformaciones de Lorentz 11
2.1 Definición . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 11
2.2 Grupos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
2.3 Representaciones de grupos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2.4 Grupos de Lie . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 14
2.5 El álgebra de Lorentz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.6 Representaciones espinoriales . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.6.1 Representaciones tensoriales . . . . . . . . . . . . . . . . . . . . . . . . . 21

3 Teoría clásica de campos 23

3.1 Campos clásicos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
3.2 Principio de mínima acción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 28
3.3 El teorema de Noether . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 30
3.4 Invariancia bajo traslaciones: tensor energía-momento . . . . . . . . . . . . . 34
3.5 Invariancia Lorentz . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 35

4 Cuantización de campos escalares libres 37

4.1 El campo escalar real clásico . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 37
4.2 Cuantización en Mecánica Cuántica . . . . . . . . . . . . . . . . . . . . . . . . . 39
4.3 Cuantización del campo escalar real . . . . . . . . . . . . . . . . . . . . . . . . . 40
4.3.1 Espacio de Fock . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 44
4.4 Causalidad . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 46
4.5 El propagador de Klein-Gordon . . . . . . . . . . . . . . . . . . . . . . . . . . . 47
4.6 Cuantización del campo escalar complejo: antipartículas . . . . . . . . . . . . 50
4.6.1 El campos escalar complejo clásico . . . . . . . . . . . . . . . . . . . . . 50

1
Campos Cuánticos 2022-23 UCM

4.6.2 Cuantización del campo escalar complejo . . . . . . . . . . . . . . . . . 52

4.7 Propagador . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 54

5 Cuantización de campos espinoriales libres 55

5.1 Campos espinoriales clásicos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 55
5.1.1 El Lagrangiano de Weyl . . . . . . . . . . . . . . . . . . . . . . . . . . . . 56
5.1.2 El Lagrangiano de Dirac . . . . . . . . . . . . . . . . . . . . . . . . . . . . 58
5.1.3 Soluciones de la ecuación de Dirac . . . . . . . . . . . . . . . . . . . . . 61
5.1.4 Corrientes, simetrías y transformaciones de Lorentz . . . . . . . . . . . 63
5.2 Cuantización del campo de Dirac . . . . . . . . . . . . . . . . . . . . . . . . . . 65
5.3 El propagador de Dirac . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 68
5.4 Fermiones de Weyl sin masa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 70

6 Cuantización del campo electromagnético libre 72

6.1 El campo electromagnético clásico . . . . . . . . . . . . . . . . . . . . . . . . . 72
6.1.1 Simetría gauge . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 73
6.1.2 Tensor energía-momento . . . . . . . . . . . . . . . . . . . . . . . . . . . 74
6.1.3 Interacciones electromagnéticas . . . . . . . . . . . . . . . . . . . . . . . 75
6.2 Cuantización del campo electromagnético libre . . . . . . . . . . . . . . . . . . 76
6.2.1 Cuantización en el gauge de radiación . . . . . . . . . . . . . . . . . . . 77
6.2.2 Propagador del fotón en el gauge de radiación . . . . . . . . . . . . . . . 80
6.2.3 Cuantización covariante . . . . . . . . . . . . . . . . . . . . . . . . . . . . 80
6.2.4 Propagador del fotón en la formulación covariante . . . . . . . . . . . . 85

7 Procesos de interacción 86
7.1 Imagen de interacción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 87
7.2 La matriz S . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 90
7.2.1 Ejemplo: desintegración de mesones . . . . . . . . . . . . . . . . . . . . 91
7.3 El teorema de Wick y los diagramas de Feynman . . . . . . . . . . . . . . . . . 94
7.3.1 Ejemplo: dispersión de nucleones . . . . . . . . . . . . . . . . . . . . . . 96
7.3.2 Diagramas de Feynman . . . . . . . . . . . . . . . . . . . . . . . . . . . . 97
7.3.3 Reglas de Feynman . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 98
7.4 Secciones eficaces y ritmos de desintegración . . . . . . . . . . . . . . . . . . . 100
7.4.1 Ritmo de desintegración . . . . . . . . . . . . . . . . . . . . . . . . . . . . 102
7.4.2 Secciones eficaces . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 103
7.4.3 Estados finales de dos partículas . . . . . . . . . . . . . . . . . . . . . . . 103
7.4.4 Ejemplo: Anchura de desintegración de mesones . . . . . . . . . . . . . 105

Antonio L. Maroto, Ignazio Scimemi 2

Campos Cuánticos 2022-23 UCM

8 Electrodinámica cuántica 107

8.1 Reglas de Feynman para QED . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 107
8.1.1 Ejemplo: e + e − → µ+ µ− . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 108
8.2 Procesos elementales . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 112
8.2.1 Sección eficaz de aniquilación e + e − → µ+ µ− . . . . . . . . . . . . . . . . 112
8.2.2 Dispersión Compton: e − γ → e − γ . . . . . . . . . . . . . . . . . . . . . . . 115

Antonio L. Maroto, Ignazio Scimemi 3

1. Introducción

Notación
Unidades naturales: ℏ=c =1
Métrica: g µν = η µν = diag(1, −1, −1, −1)

Cuadrivectores: p µ = (p 0 , p
⃗), (⃗
p corresponde al trivector con los índices arriba).

Productos escalares: p · x ≡ p µ x µ = η µν p µ x ν = p 0 x 0 − p
⃗ ·⃗
x

Productos vectoriales: ⃗=⃗

C A ×B ⃗ , C i = ϵi j k A j B k con ϵ123 = 1.
i jk ilm
ϵ ϵ = δ j l δkm − δ j m δkl

Derivadas: ∂i = ∂/∂x i , (las componentes de ⃗

∇ son ∂i )
↔
f ∂µ g = f ∂µ g − g ∂µ f
R n i px
Delta de Dirac: d xe = (2π)n δ(n) (p)

Matrices de Pauli:
µ ¶ µ ¶ µ ¶
1 0 1 2 0 −i 3 1 0
σ = , σ = , σ =
1 0 i 0 0 −1

[σi , σ j ] = 2i ϵi j k σk

1.1 Introducción
La Teoría Cuántica de Campos surge históricamente del intento de contruir una teoría
cuántica para sistemas de partículas relativistas. La Mecánica Cuántica que habéis es-
tudiado en cursos anteriores es fundamentalmente una teoría no relativista que permite

4
Campos Cuánticos 2022-23 UCM

describir sistemas con un número fijo de partículas que pueden interaccionar con fuerzas
externas. Con el desarrollo de la Física Nuclear y de Partículas a partir de los años 20, fue
necesario encontrar un formalismo cuántico válido para estos sistemas que involucraban
partículas relativistas y muy especialmente para los fotones asociados a la interacción elec-
tromagnética.
Los primeros intentos realizados para encontrar una extensión relativista de la ecuación
de Schrödinger se deben a Vladimir Fock, Oskar Klein y Walter Gordon en 1926. A contin-
uación realizaremos un breve recorrido histórico de las dificultades encontradas en este
proceso que finalmente conduciría a la construcción de la Teoría Cuántica de Campos.
La ecuación de Schrödinger ordinaria puede obtenerse de la expresión no relativista
para la energía de un partícula libre (sin interacciones) de masa m

⃗2
p
E= (1.1)
2m
mediante la identificación con los operardores

∂
E → iℏ (1.2)
∂t
⃗ → −i ℏ⃗
p ∇ (1.3)

de forma que sustituyendo obtenemos la ecuación de onda

∂Ψ ℏ2 2
µ ¶
iℏ = − ∇ Ψ (1.4)
∂t 2m

donde Ψ es la función de onda de la partícula que tiene una interpretación probabilística.

Esta ecuación es invariante bajo transformaciones de Galileo1 , pero no bajo transforma-
ciones de Lorentz.
Para partículas relativistas podemos utilizar en su lugar la relación de la energía total
relativista

E 2 = (⃗
p 2c 2 + m2c 4) (1.5)

de forma que utilizando los operadores anteriores llegamos a

∂2 Ψ ¡ 2 2 2
−ℏ2 2 4
Ψ
¢
= −ℏ c ∇ + m c (1.6)
∂t 2
1
La función de onda adquiere una fase en la transformación

Antonio L. Maroto, Ignazio Scimemi 5

Campos Cuánticos 2022-23 UCM

que en unidades naturales puede escribirse como

(□ + m 2 )Ψ = 0 (1.7)

que es la llamada ecuación de Klein-Gordon. Esta ecuación, es invariante bajo transfor-

maciones de Lorentz, puesto que □Ψ = ηµν ∂µ ∂ν Ψ es un escalar Lorentz si Ψ lo es. Sin
embargo, si queremos interpretar Ψ como función de onda de una partícula individual,
la extensión relativista de la Mecánica Cuántica encuentra importantes dificultades. En
efecto, puesto que ahora la relación entre la energía y el momento depende de E 2 tenemos
que las soluciones en forma de ondas planas pueden tomar dos formas distintas

Ψ+ (t , x) = C p e −i E p t +i p⃗·⃗x (Modos de frecuencia positiva) (1.8)

Ψ− (t , x) = C p e i E p t +i p⃗·⃗x (Modos de frecuencia negativa) (1.9)

donde C p es la constante de normalización y

q
Ep = p ⃗2 + m 2 (1.10)

con E p ≥ 0. Los modos de frecuencia positiva corresponden a funciones de onda de partícu-

las con momento p ⃗ y energía E p , ya que

−i ⃗ ⃗ Ψ+
∇ Ψ+ = p (1.11)
∂
i Ψ+ = E p Ψ+ (1.12)
∂t
Sin embargo, los modos con frecuencia negativa corresponderían a partículas con mo-
mento p ⃗ pero energía negativa −E p .
Por tanto, a diferencia del caso no relativista en el que la energía de las partículas es
siempre positiva, el espectro de la ecuación de Klein-Gordon correspondería a |E | ≥ m (ver
Fig. 1.1). Esto es, existen soluciones con energía arbitrariamente negativa (el espectro no
está acotado inferiormente). En el caso de una teoría sin interacciones, esto no plantea
ningún problema, pero si incluímos algún tipo de interacción, el estado con E = m (estado
con momento cero que sería el estado fundamental en el caso no relativista) podría decaer
en un estado con energía negativa generando una cantidad arbitraria de radiación con
energía positiva, es decir no habría estados estables para la partícula.
Una posible solución a este problema fue propuesta por Paul Dirac en 1928. Si supone-
mos que las partículas son fermiones, éstas deberían cumplir el principio de exclusión de
Pauli que dice que cada estado del espectro sólo puede estar ocupado por una partícula.
Ahora bien, si además suponemos que todos los estados con energía E < −m están ya ocu-
pados, entonces no serían posibles las transiciones a dichos estado y por tanto el estado

Antonio L. Maroto, Ignazio Scimemi 6

Campos Cuánticos 2022-23 UCM

Figure 1.1: Espectro de la ecuación de Klein-Gordon

con E = m sería estable. Esta configuración en la que los estados con energía negativa están
ocupados es lo que se conoce como mar de Dirac. Dirac obtuvo además la versión rela-
tivista de la ecuación de Schrödinger para partículas de spin 1/2 (ecuación de Dirac) que
como veremos no coincide con la ecuación de Klein-Gordon vista anteriormente. Pero lo
que es más importante, la solución de Dirac al problema de los estados de energía negativa
llevó a una predicción sorprendente, la existencia de antipartículas. En efecto, considere-
mos el caso en el que un fotón con energía suficientemente alta E > 2m colisiona con un
fermión del mar. En este caso, como vemos en la Fig. 1.1, podría excitarlo a un estado con
E > m y dejaría por tanto un hueco en el mar de Dirac. Estos huecos se podrían entender
como antipartículas, es decir partículas con la misma masa, pero carga opuesta, lo que
llevó a Dirac a postular la existencia del positrón (antipartícula del electrón) que de hecho
fue posteriormente descubierto por Carl Anderson en 1932.
A pesar del éxito de la ecuación de Dirac en la interpretación de los estados de energía
negativa, otros nuevos problemas surgieron inmediatamente. Por un lado, esta solución
sólo es válida para partículas fermiónicas y, por otra, una posible inconsistencia aparece
en el caso de interacción de una partícula relativista con un potencial escalón. Ésta es la

Antonio L. Maroto, Ignazio Scimemi 7

Campos Cuánticos 2022-23 UCM

Figure 1.2: Paradoja de Klein

llamada paradoja de Klein descubierta por Oskar Klein en 1929. 2 .

Consideremos una partícula libre de masa m descrita por la ecuación de Klein-Gordon
con energía positiva E que incide desde la izquierda sobre una barrera de potencial de
altura V0 (ver Fig. 1.2). La solución en el lado izquierdo de la barrera corresponderá a
una onda plana incidente propagándose hacia la derecha más una onda reflejada que se
propaga hacia la izquierda de forma que

ΨI (t , x) = e −i E t +i p 1 x + Re −i E t −i p 1 x (1.13)

Por otra parte, a la derecha de la barrera, tendremos sólo la onda transmitida

ΨI I (t , x) = Te −i E t +i p 2 x (1.14)

donde
p
p1 = E 2 − m2
q (1.15)
p2 = (E − V0 )2 − m 2

Nótese que todas las soluciones tienen energía positiva. Imponiendo ahora continuidad de
la función de onda y su derivada sobre la barrera, es decir, ΨI (t , 0) = ΨI I (t , 0) y ∂x ΨI (t , 0) =
2
La paradoja de Klein fue originalmente deducida para la ecuación de Dirac, pero está presente también
en la ecuación de Klein-Gordon

Antonio L. Maroto, Ignazio Scimemi 8

Campos Cuánticos 2022-23 UCM

∂x ΨI I (t , 0) obtenemos para los coeficientes de reflexión y transmisión

2p 1
T=
p1 + p2
p1 − p2 (1.16)
R=
p1 + p2

Formalmente la expresión para los coeficientes de reflexión y transmisión es la misma que

en el caso de Schrödinger, pero las expresiones de p 1 y p 2 son distintas.
En el caso no relativista3 la solución del problema implicaba que si:

• E > V0 tanto p 1 como p 2 son reales de forma que parte de la onda se refleja y parte
se transmite.

• E < V0 tenemos que p 2 es imaginario puro, la onda no se transmite a través de la

barrera y la solución es exponencialmente decreciente dentro de la barrera con |R|2 =
1, tenemos por tanto reflexión total.

En el caso relativista, tenemos de nuevo dos casos

• E > V0 + m (la energía de la partícula es mayor que la energía en reposo más la de la

barrera), p 1 y p 2 son reales y tenemos de nuevo parte transmitida y parte reflejada.

• E < V0 + m ahora sin embargo hay dos posibilidades:

– E > V0 − m entonces p 2 es imaginario puro y tenemos de nuevo reflexión total.

– E < V0 − m (nótese que E > m siempre para tener p 1 real y por tanto esta condi-
ción implica V0 > 2m) en este caso encontramos sorprendentemente que tanto
p 1 como p 2 son reales4 y tenemos tanto onda reflejada como transmitida, a pe-
sar de no tener energía suficiente para superar la barrera (Paradoja de Klein).

Esta inconsistencia de nuevo se presenta al intentar interpretar la función Ψ(t , x) como

una función de onda para una partícula individual. Una interpretación de Ψ(t , x) en tér-
minos de múltiples partículas permite describir el proceso como una reflexión de una
partícula en una barrera que es seguido por la producción de un par partícula-antipartícula
a partir de la energía de la barrera. Nótese que este caso patológico corresponde a V0 > 2m,
es decir, la barrera de potencial tiene energía suficiente para producir el par partícula-
antipartícula.
3
p p
En el caso no relativista p 1 = 2mE , p 2 = 2m(E − V0 )
4
Para poder tener p 2 real hay dos posibilidades E − V0 > m o V0 − E > m

Antonio L. Maroto, Ignazio Scimemi 9

Campos Cuánticos 2022-23 UCM

A la vista de esta paradoja, resulta fácil de entender por qué no es posible cuantizar
partículas relativistas como en el caso no relativista, es decir construyendo una ecuación
para las funciones de onda de partículas individuales. Esto se debe a que no tiene sentido
suponer que en un proceso relativista se conserve el número de partículas. En efecto, la
equivalencia entre masa y energía de la Relatividad Especial permite la creación de pares
reales partícula-antipartícula. Incluso si la energía no es suficiente para la creación de
los pares, el principio de indeterminación de Heisenberg nos dice que la energía de un
sistema tiene una indeterminación dada por ∆E ∆t ≥ h. Por tanto dicha energía puede
sufrir fluctuaciones importantes durante un intervalo de tiempo suficientemente corto y
por tanto permite también la creación de pares virtuales partícula-antipartícula.
La solución a estos problemas surgió unos años más tarde. A partir del año 1930, el pro-
pio P. Dirac, junto con V. Fock, P. Jordan o E. Wigner, encontraron que sería posible dar sen-
tido a las ecuaciones relativistas, no como generalizaciones de la ecuación de Schrödinger
para partículas, sino como ecuaciones del movimiento de campos. Las excitaciones de
estos campos son las que podrían interpretarse como partículas. Este proceso se conoce
también como segunda cuantización. Así por ejemplo, los fotones serían las excitaciones
del campo electromagnético. En este formalismo, el número de partículas no se conser-
varía (las partículas podrían crearse y destruirse) y por tanto la noción de partícula como
objeto fundamental con identidad física propia se pierde. En su lugar el campo se convierte
en el objeto fundamental. Nótese que si las partículas se consideran las excitaciones de un
cierto campo, podemos tener tantas partículas como queramos y por tanto el número de
grados de libertad de nuestro sistema será en general infinito. El espacio de Hilbert aso-
ciado a la Teoría Cuántica de Campos se conoce como espacio de Fock. En resumen, la
extensión relativista consistente de la Mecánica Cuántica no es una Mecánica Cuántica
Relativista sino una Teoría Cuántica de Campos.
La Teoría Cuántica de Campos es a día de hoy la herramienta fundamental para enten-
der la Física de los consituyentes fundamentales de la materia. Proporciona un formalismo
potente y preciso que permite describir procesos en un rango de energías desde los MeV
a los TeV alcanzados en el LHC y permite realizar predicciones con un enorme grado de
precisión.
Sin embargo y a pesar de que los sistemas relativistas son su campo de aplicación por
excelencia, los métodos de la Teoría Cuántica de Campos tienen una gran generalidad y
pueden aplicarse en otras muchas ramas de la Física en las que nos enfrentemos a fenó-
menos colectivos que involucren un gran número de grados de libertad. Así por ejemplo,
en Física de la Materia Condensada, las excitaciones en sólidos pueden entenderse tam-
bién como partículas asociadas a campos y pueden estudiarse utilizando las técnicas de la
teoría de campos.

Antonio L. Maroto, Ignazio Scimemi 10

2. Transformaciones de Lorentz

En este capítulo repasaremos la definición y las propiedades de las transformaciones de

Lorentz entendidas como grupo de Lie. No entramos en demasiados detalles, porque son
objeto de otros cursos, sino simplemente repasaremos los conceptos básicos necesarios
para introducir la teoría clásica de campos.

2.1 Definición
Las transformaciones de Lorentz se definen como las transformaciones de las coordenadas
espacio-temporales de un suceso
µ ν
x ′µ = Λ νx (2.1)

que dejan invariante la cantidad

x · x = η µν x µ x ν = t 2 − x 2 − y 2 − z 2 (2.2)

es decir,

x ′ · x ′ =η µν x ′µ x ′ν = η µν x µ x ν = x · x (2.3)
µ
Por tanto la condición que debe satisfacer la matriz Λ ν es
µ ρ
η µν x ′µ x ′ν = η µν (Λ ρx )(Λν λ x λ ) = η ρλ x ρ x λ (2.4)

Puesto que esta condición se debe cumplir para un x ρ arbitrario tenemos que
µ ν
η µν Λ ρΛ λ = η ρλ (2.5)

que en notación matricial se escribe

ΛT ηΛ = η, (2.6)

11
Campos Cuánticos 2022-23 UCM

Tomando el determinante de la ecuación anterior encontramos

(det Λ)2 = 1 =⇒ (det Λ) = ±1. (2.7)

Sólo las transformaciones con det Λ = 1 pueden conectarse de forma continua con la iden-
tidad y se conocen como transformaciones de Lorentz propias. 1 .
Por otra parte, si tomamos la componente (00) de la ecuación (2.5) tenemos

(Λ0 0 )2 − (Λi 0 )2 = 1
X
(2.8)
i

y por tanto (Λ0 0 )2 ≥ 1. Es decir, el grupo de Lorentz propio tiene dos componentes de-
sconectadas, una con Λ0 0 ≥ 1 y otra con Λ0 0 ≤ −1, conocidas como ortocrona y no or-
tocrona. Sólo las transformaciones ortocronas pueden conectarse de forma continua con
la identidad. Por tanto nos limitaremos a las transformaciones de Lorentz propias or-
tocronas con Λ0 0 ≥ 1. Estas consisten en las rotaciones y los boosts dados por

t ′ = γ(t + ⃗
v ·⃗
x)
x ∥′ = γ(⃗
⃗ x∥ + ⃗
vt)
′
⃗
x⊥ =⃗
x⊥ (2.9)

donde γ = (1 − v 2 )−1/2 con v 2 < 1. 2

2.2 Grupos
Dado un conjunto de transformación, podemos comprobar si consituyen un grupo en un
sentido matemático estricto. Si este es el caso, podemos aprovechar todo el conjunto de
propiedades de este tipo de objetos. Por ejemplo, hemos comenzado a definir las transfor-
maciones de Lorentz en el espacio de coordendas espacio-temporales (4 dimensiones). Sin
embargo, podemos encontrar un conjunto de transformaciones con las mismas propiedades
1
El grupo de matrices ortogonales de dimensión n = 4 que cumplen A T A = 1 se denomina O(4). El grupo
que cumple ΛT ηΛ = η es el grupo O(3, 1). Si además tienen determinante positivo se denota SO(3, 1).
2
El conjunto de transformaciones de un espacio vectorial 4D que deja invariante la forma s 2 (x) = (x 0 )2 −⃗
x2
también se denomina grupo de Lorentz homogéneo (como hemos visto incluye rotaciones, boosts, inversión
espacial I s = diag(1, −1, −1, −1) e inversión temporal I t = diag(−1, 1, 1, 1)). El conjunto de transformaciones
de un espacio vectorial 4D que deja invariante la forma s 2 (x, y) = (x 0 − y 0 )2 − (⃗ x −⃗
y )2 se denomina grupo
de Lorentz no homogéneo o grupo de Poincaré (incluye rotaciones, boosts, inversión espacial, inversión
temporal y traslaciones).

Antonio L. Maroto, Ignazio Scimemi 12

Campos Cuánticos 2022-23 UCM

en espacios vectoriales de diferentes dimensiones. Esta es la base de la clasificación de

magnitudes físicas en escalares, fermiónicas, vectoriales, etc.
Un grupo G está definido por las siguientes reglas genéricas. Necesitamos un conjunto
de objetos (a, b, c, . . . ) con una ley de composición a · b tal que

• ∀ a, b ∈ G tenemos a · b ∈ G ;

• Existe el elemento identidad e dentro del grupo, tal que a · e = e · a = a, ∀a ∈ G ;

• Existe un inverso para cada elemento tal que si a ∈ G ⇒ a −1 ∈ G y a · a −1 = a −1 · a = e

• Se cumple la propiedad asociativa, a · (b · c) = (a · b) · c, ∀a, b, c ∈ G .

Si a · b = b · a, ∀a, b ∈ G el grupo se dice abeliano.

El grupo de Lorentz ortocrono propio forma por tanto un grupo no abeliano.

2.3 Representaciones de grupos

Una representación (lineal) R de un grupo es una operación que asigna a cada elemento
génerico abstracto g un operador lineal D R (g ) (una matriz en el caso de dimensión finita)
definido sobre un espacio lineal

g → D R (g ) (2.10)

con las siguientes propiedades

1) D R (e) = 1, donde 1 es el operador identidad.

2) D R (g 1 )D R (g 2 ) = D R (g 1 g 2 ), es decir, la estructura de grupo se preserva.

El espacio sobre el que actúan los operadores D R (g ) se denomina espacio base. El ejemplo
más simple de representación son las representaciones matriciales en las que asignamos
una matriz a cada elemento del grupo. En este caso la base es un espacio vectorial de
dimensión finita n y por tanto D R (g ) son matrices n × n con componentes D R (g )ij , i , j =
1...n. La dimensión de la representación es simplemente la dimensión n del espacio base.
Por tanto podemos asignar distintas representaciones a un mismo grupo.
Supongamos ahora un elemento del espacio base (φ1 , φ2 , ...φn ). Un elemento g del
grupo induce por tanto una transformación de los elementos del espacio base dada por

φi → D R (g )ij φ j . (2.11)

Antonio L. Maroto, Ignazio Scimemi 13

Campos Cuánticos 2022-23 UCM

Por ejemplo, en nuestra definición de las transformaciones de Lorentz (2.1), el espacio

base sería el espacio vectorial R4 de las coordenadas espacio-temporales, y por tanto n =
4. Podemos por tanto tomar (φ1 , φ2 , ...φn ) = (x 0 , x 1 , x 2 , x 3 ) y las D R (g )ij serían una repre-
µ
sentación de dimensión n = 4 dada por las matrices Λ ν . Sin embargo, podríamos con-
struir representaciones de otras dimensiones que actuaran sobre otros objetos (espinores,
tensores, etc). Finalmente:

1) Una representación se llama reducible si hay un subespacio invariante, es decir, todas

las transformaciones D R (g ) envían los elementos del subespacio en el mismo sube-
spacio.

2) Dos representaciones son equivalentes si hay un cambio de base que lleva una a la
otra. Así R, R ′ son equivalentes si hay una matriz para el cambio de base S tal que
∀g , D R ′ (g ) = S −1 D R (g )S.

3) La representación (no trivial) de dimensión más baja se suele denominar representación

fundamental.

2.4 Grupos de Lie

Un conjunto particular de grupos muy útiles en física son los llamados grupos de Lie. Un
grupo de Lie es un grupo cuyos elementos g (θ) dependen de forma continua y diferen-
ciable de un conjunto de parámetros reales θ a , a = 1, ..N . Por ejemplo, las traslaciones y
rotaciones son grupos de Lie. Para las traslaciones tenemos

x ′µ = x µ + a µ (2.12)

de modo que a µ es un conjunto de 4 parámetros. Para las rotaciones tenemos 3 ángulos,

x ′0 = x 0
x ′i = R ij (θ1 , θ2 , θ3 )x j (2.13)

y θk definen la rotación con respecto a cada eje.

En general, se suelen elegir los parámetros de forma que θ a = 0 corresponda a la iden-
tidad del grupo, esto es g (0) = e.
El grupo de Lorentz es también un grupo de Lie que depende de 6 parámetros (3 ángu-
los correspondientes a las rotaciones y las 3 componenes del vector velocidad correspon-
dientes a los boosts).

Antonio L. Maroto, Ignazio Scimemi 14

Campos Cuánticos 2022-23 UCM

Transformaciones infinitesimales, generadores del grupo y álgebra de Lie

Consideremos una cierta representación matricial de un elemento de un grupo de Lie g (θ)
que denotaremos simplemente como D R (θ). Por las hipótesis iniciales de continuidad y
diferenciabilidad de los grupos de Lie respecto a los parámetros, para θ a infinitesimales,
es decir, cerca del elemento identidad tendremos
∂D R ¯¯
¯
a a
D R (θ) ≃ 1 + i θa TR con TR ≡ −i (2.14)
∂θ ¯ a θ=0
los TRa , con a
= 1, . . . , N , se llaman los generadores del grupo que son por tanto matrices de
la misma dimensión que la representación. Nótese que tenemos tantos generadores como
parámetros independientes. Es decir, para las traslaciones tendríamos 4 generadores, para
las rotaciones 3, etc. El factor i se elige de modo que, si los generadores son hermíticos
entonces D R es unitario, es decir, tenemos una representación unitaria.
Nótese que siempre podemos escribir una transformación finita como límite de m
transformaciones infinitesimales,
θa TRa m
µ ¶
a
D R (θ) = lim 1 + i ≡ e i θa TR (2.15)
m→∞ m
y esta es la forma que generalmente se considera en muchas aplicaciones. Antes de es-
tablecer la forma de los TRa para nuestros casos relevantes, hay una propiedad importante
de los generadores que tenemos que tratar: el álgebra del grupo. En general, el producto
de dos transformaciones del grupo debe dar una tercera transformación del grupo,
a b c
e i α a T R e i βb T R = e i δc T R (2.16)
Debemos por tanto determinar δc = δc (α, β). Nótese que para matrices A y B , en general
e A e B ̸= e A+B y por tanto δa ̸= αa + βa . El truco para derivar δc es el siguiente. Tomamos el
logaritmo en ambos lados
a b
i δc TRc = ln(e i αa TR e i βb TR ) (2.17)
Desarrollamo el miembro derecho hasta segundo orden en los parámetros tenemos
½ i¾
c a 1h a 2 a 2 a b
i δc TR = ln 1 + i (αa + βa )TR − (αa TR ) + (βa TR ) + 2αa βb TR TR (2.18)
2
Usando ahora que ln(1 + x) ≃ x − x 2 /2 y teniendo en cuenta que los generadores no con-
mutan, obtenemos hasta segundo orden,
1h i 1
i δc TRc = i (αa + βa )TRa − (αa TRa )2 + (βa TRa )2 + 2αa βb TRa TRb + (αa + βa )(αb + βb )TRa TRb
2 2
(2.19)

Antonio L. Maroto, Ignazio Scimemi 15

Campos Cuánticos 2022-23 UCM

de donde

αa βb [TRa , TRb ] = −2i (δc (α, β) − αc − βc )TRc ≡ i γc (α, β)TRc (2.20)

Puesto que (2.20) debe ser cierta ∀α y ∀β, γc debe ser lineal en α y β. Es decir podemos
escribir en general γc (α, β) = αa βb f abc con f abc constantes y por tanto

[T a , T b ] = i f abc T c (2.21)

Estas relaciones de conmutación definen el álgebra de Lie del grupo y las f abc se denomi-
nan constantes de estructura.
Nótese que las constantes de estructura son las mismas para todas las representaciones
del grupo. Esto se debe a que el producto de dos elementos del grupo siempre da el mismo
tercer elemento y por tanto las relaciones de conmutación no pueden depender de la rep-
resentación. Estos resultados se han obtenido hasta segundo orden en el desarrollo de
Taylor pero se puede comprobar que son ciertos a cualquier orden.

2.5 El álgebra de Lorentz

En esta sección aplicaremos los resultados anteriores a las transformaciones de Lorentz.
Consideremos por tanto una transformación de Lorentz infinitesimal
µ µ µ
Λ ν =δ ν +ω ν (2.22)

de forma que sustituyendo en (2.5) tenemos que

ωµν = −ωνµ (2.23)

es decir, se trata de una matriz antisimétrica 4 × 4 que tiene por tanto 6 elementos inde-
pendientes. El grupo de Lorentz por tanto tiene 6 parámetros independientes. Estos 6
parámetros independientes, como comentamos anteriormente, se pueden identificar con
los ángulos de 3 rotaciones y las 3 componentes de la velocidad ⃗ v del boost.
Denotaremos a los 6 generadores por tanto como J con J = −J νµ (en vez se usar un
µν µν

único índice a usamos un par de índices µν) y por tanto un elemento genérico del grupo
se puede escribir como
i µν
Λ = e − 2 ωµν J (2.24)

El factor 1/2 se introduce para no contar dos veces el mismo generador.

Antonio L. Maroto, Ignazio Scimemi 16

Campos Cuánticos 2022-23 UCM

Por tanto un conjunto de objetos φi con i = 1, . . . , n se transformarán en una repre-

sentación R de dimensión n si bajo transformaciones de Lorentz cambian como
h i µν ii
φi → φ′ = e − 2 ωµν J R
i
φj (2.25)
j

µν
donde J R son los generadores en la representación R. Por tanto, para una transformación
infinitesimal, tendremos que

i i µν
δφi = φ′ − φi = − ωµν [J R ]ij φ j (2.26)
2
µν
Nótese que en [J R ]ij , los índices µ, ν = 0, 1, 2, 3 etiquetan al generador e i , j = 1, 2, ...n son
los índices matriciales de la representación.
Una observación importante es que todas las cantidades físicas pueden clasificarse de
acuerdo con su transformación con respecto al grupo de Lorentz. En efecto:
Un escalar es una cantidad que no cambia bajo transformaciones de Lorentz

φ′ = φ (2.27)

y por tanto J µν = 0. La masa de una partícula por ejemplo es un escalar Lorentz.

Un cuadrivector contravariante se transforma como
µ
V µ → Λ νV ν

donde Λ aquí es una transformación de Lorentz estándar en 4 dimensiones. Un cuadrivec-

tor covariante se transforma como

Vµ → ΛµνVν

con
ρ
Λµν = η µρ ηνσ Λ σ , Vµ = η µνV ν .
Las coordenadas del espacio-tiempo x µ , los cuadrimomentos p µ = (E , p
⃗) son ejemplos de
cuadrivectores.
Para vectores los generadores deben ser matrices 4 × 4 matrices, por tanto podemos
usar para i , j también índices griegos de forma que

i ρ
δV ρ = − ωµν [J µν ]σV σ (2.28)
2

Antonio L. Maroto, Ignazio Scimemi 17

Campos Cuánticos 2022-23 UCM

Ahora bien usando (2.22), tenemos que

ρ
δV ρ = ω σV σ (2.29)

y por tanto la forma antisimetrizada de los generadores en (µν) es directamente

ρ µ
[J µν ] σ = i (ηµρ δνσ − ηνρ δ σ ) (2.30)

Esta es la forma de los generadores en la representación de 4 dimensiones. A partir de esta

expresión podemos obtener el álgebra de Lie simplemente calculando los conmutadores.
En efecto,

[J µν , J ρσ ] = i (ηνρ J µσ − ηµρ J νσ − ηνσ J µρ + ηµσ J νρ ). (2.31)

Este es el álgebra de Lorentz o también álgebra de SO(3, 1). Nótese que aunque hemos
obtenido la expresión del álgebra para la representación de 4 dimensiones, esta expresión
es válida para cualquier representación. Como veremos esto resultará muy útil más ade-
lante para construir las representaciones espinoriales.
Para tener algo de intuición sobre la acción de los generadores, vamos a separar los 6
generadores J µν en dos grupos de 3 generadores,

1
J i = ϵi j k J j k , K i = Ji0 . (2.32)
2
Por tanto las relaciones de commutación se escriben

[J i , J j ] = i ϵi j k J k (2.33)
i j i jk k
[J , K ] = i ϵ K (2.34)
[K i , K j ] = −i ϵi j k J k (2.35)

Vemos que la primera ecuación corresponde al álgebra del momento angular. De hecho,
los J i son los generadores del grupo de rotaciones en 3 dimensiones SO(3), que tiene el
mismo álgebra que el grupo SU(2). Por otra parte, K i son los generadores de los boosts.
Comprobemos esta interpretación. Para ello vamos a descomponer una transformación
de Lorentz genérica de la siguiente forma

1
θ i = ϵi j k ω j k
2
η = ωi 0
i
(2.36)

Antonio L. Maroto, Ignazio Scimemi 18

Campos Cuánticos 2022-23 UCM

de forma que seguimos teniendo 6 parámetros. Ahora podemos escribir

1 3
ωµν J µν = ω12 J 12 + ω13 J 13 + ω23 J 23 + ωi 0 J i 0 = ⃗
θ ·⃗ ⃗
J −⃗η · K
X
(2.37)
2 i =1

donde se ha usado ωi 0 = −ωi 0 = −ηi , ω12 = ω12 = θ 3 , etc. Por tanto una transformación de
Lorentz genérica se puede escribir como
⃗⃗ ⃗
Λ = e −i (θ· J −⃗ηK ) (2.38)

Consideremos como ejemplo la transformación generada por el parámetro θ 3 = ω12

sobre el cuadrivector x µ . Es decir:
µ
δx µ = −i ω12 [J 12 ] ν x ν = θ 3 (η1µ δ2ν − η2µ δ1ν )x ν (2.39)

de donde δx = −θ 3 y y δy = θ 3 x, que correponde a una rotación de ángulo θ 3 alrededor del

eje Z . Luego efectivamente los parámetros θ i corresponden a rotaciones alrededor del eje
x i en sentido antihorario. De forma análoga puede verse que ηi corresponde a un boost a
lo largo del eje x i .

2.6 Representaciones espinoriales

Hasta ahora hemos visto las representaciones escalar (dimensión 1) y vector (dimensión
4) del grupo de Lorentz. Como veremos es posible construir representaciones de dimen-
sión 2, conocidas como representaciones espinoriales que constituyen la representación
fundamental del grupo de Lorentz.
Los generadores J i y K i definidos anteriormente se mezclan entre sí por las propiedades
del álgebra como puede verse en (2.34) y (2.35). Es posible, sin embargo, definir dos con-
juntos de generadores que no se mezclen. En concreto,
1
M i = (J i + i K i )
2
1
N = (J i − i K i )
i
(2.40)
2
cumplen

[M i , M j ] = i ϵi j k M k
[N i , N j ] = i ϵi j k N k
[M i , N j ] = 0 (2.41)

Antonio L. Maroto, Ignazio Scimemi 19

Campos Cuánticos 2022-23 UCM

Vemos que tanto los generadores M i como los de N i cumplen el álgebra de SU(2). Se dice
que el grupo de Lorentz es isomorfo al grupo SU(2)×SU(2). Por tanto, si consideramos
transformaciones que sólo involucren a los generadores M i , éstas cumplirán el álgebra de
Lorentz (2.31) sin involucrar a los generadores N i y viceversa. Por tanto, tanto M i como
N i pueden generar representaciones independientes del grupo de Lorentz.
En efecto, consideremos el caso
i
Mi = 0 J i = σ2
⇒ (2.42)
σi σi
Ni = 2
Ki =i 2

con σi las matrices de Pauli 2 × 2. Puede verse que esta elección cumple el álgebra de
Lorentz. Se trata por tanto de una representación de dimensión 2 que actúa en un espa-
cio base de dimensión 2. Denotemos por ψα L con α = 1, 2 los elementos de este espacio
vectorial. Los elementos del espacio base se transforman por tanto como
⃗ ⃗
σ
ψL → ΛL ψL = e (−i θ−⃗η)· 2 ψL (2.43)

Los objetos que se transforman de esta forma se conocen como espinores de Weyl (levó-
giros) y como veremos representan a partículas de spin 1/2 sin masa (neutrinos sin masa).
De forma totalmente análoga podríamos construir la representación
i i
M i = σ2 J i = σ2
⇒ (2.44)
i i σi
N =0 K = −i 2

En este caso, los elementos del espacio base se denominan espinores de Weyl dextrógiros
ψα
R con α = 1, 2 y

⃗ ⃗
σ
ψR → ΛR ψL = e (−i θ+⃗η)· 2 ψR (2.45)

Dado que el grupo de Lorentz puede entenderse como dos copias3 del grupo de rota-
ciones asociado al momento angular, podemos asociar un valor de momento angular a
cada una de ellas. De esta forma podemos etiquetar cada representación como ( j L , j R ). La
3
x ) → (t , −⃗
Bajo paridad (t ,⃗ x ) los generadores de los boosts cambian de signo puesto que la velocidad
cambia de signo K⃗ → −K ⃗ , mientras que los generadores de las rotaciones no cambian ⃗ J →⃗ J . Es decir, los
generadores en las dos representaciones se intercambian entre sí M ⃗ ↔N⃗ , por lo que los espinores L pasan a
transformarse como R y viceversa.

Antonio L. Maroto, Ignazio Scimemi 20

Campos Cuánticos 2022-23 UCM

dimensión total de la representación será por tanto n = (2 j L + 1)(2 j R + 1). En efecto, la rep-
resentación escalar, de dimensión n = 1 correponde a (0, 0). Las dos representaciones de
Weyl de dimensión n = 2, corresponden a los momentos angulares (1/2, 0) y (0, 1/2). Puede
verse que la representación (1/2, 1/2) de dimensión n = 4 es la representación vectorial,
etc.
Existe además una representación reducible muy interesante que como veremos per-
mite describir fermiones con masa que se conoce como representación de Dirac. Esta
L
corresponde a (1/2, 0) (0, 1/2) y tiene dimensión n = 4. La forma de las transformaciones
es simplemente:

ΛL 0
µ ¶
ΛD = (2.46)
0 ΛR

y los elementos del espacio base de dimensión 4, se conocen como espinores de Dirac

ψL
µ ¶
ψD = (2.47)
ψR

que como vemos podemos escribir en términos de los espinores de Weyl4 . Como veremos
los espinores de Dirac describen a los leptones (electrones, muones y taus) y a los quarks.

2.6.1 Representaciones tensoriales

Es posible construir representaciones de dimensión superior. En concreto, un tensor T µν
dos veces contravariante se transforma bajo transformaciones de Lorentz como
µ
T µν → Λ ρ Λνσ T ρσ (2.48)

Nótese que T µν tiene 16 componentes independientes que se mezclan entre sí bajo una
transformación de Lorentz. Por tanto, la representación anterior tiene dimensión n = 16.
En general un tensor con un cierto número de índices arriba y cierto número de índices
µ ρ
abajo se transformará con un factor Λ ρ por cada índice arriba y un factor Λµ para cada
índice abajo. En general las representaciones tensoriales como la anterior son representa-
ciones reducibles. El campo gravitatorio en Relatividad General se representa a través de
un campo tensorial simétrico g µν .
Hasta aquí nos hemos limitado a estudiar representaciones de dimensión finita que
son adecuadas para describir magnitudes de partículas individuales. Sin embargo, en el
4
Más concretamente esta es la expresión de un espinor de Dirac en la llamada base quiral

Antonio L. Maroto, Ignazio Scimemi 21

Campos Cuánticos 2022-23 UCM

próximo tema veremos que también es posible estudiar las propiedades de transforma-
ción de campos. Los campos son funciones y por tanto están definidas en un espacio base
de dimensión infinita. Veremos por tanto cómo construir estas representaciones de di-
mensión infinita que vendrán dadas por operadores diferenciales en vez de matrices.

Antonio L. Maroto, Ignazio Scimemi 22

3. Teoría clásica de campos

El principal objetivo de este tema es estudiar los campos clásicos y sus propiedades de
transformación bajo transformaciones de Lorentz. Esto nos permitirá construir acciones
invariantes Lorentz que a su vez nos proporcionarán la dinámica de dichos campos.

3.1 Campos clásicos

Un campo clásico es simplemente una función de las coordenadas espacio-temporales
φ(x). Esta función puede tener varias componentes φ(x) = (φ1 (x), . . . , φn (x)). Bajo trans-
formaciones de Lorentz
µ ν
x ′µ = Λ νx (3.1)

el campo φ(x) tendrá ciertas propiedades de transformación, de forma que

φ(x) → φ′ (x ′ ) (3.2)

Definir cómo se transforma un campo es por tanto determinar φ′ (x ′ ) en términos de φ(x).

Nótese que la transformación afecta tanto a las componentes del campo i = 1, . . . , n, que se
podrán mezclar entre sí, como a las coordenadas. En efecto, consideremos una transfor-
mación infinitesimal de forma que podemos escribir

δφ(x) = φ′ (x ′ ) − φ(x) = φ′ (x + δx) − φ(x) = φ′ (x) − φ(x) + δx ρ ∂ρ φ(x) + O (δx 2 ). (3.3)

En esta expresión hemos separado la parte de la transformación que depende de δx de la

parte que no depende del cambio en las coordenadas. De hecho podemos definir

δ0 φ(x) = φ′ (x) − φ(x) (3.4)

de forma que

δφ(x) = δ0 φ(x) + δx ρ ∂ρ φ(x) + O (δx 2 ). (3.5)

23
Campos Cuánticos 2022-23 UCM

Utilizando la forma infinitesimal de las transformaciones Lorentz (2.22) podemos es-

cribir
ρ i ρ
δx ρ = ω σx
σ
= − ωµν [J µν ]σ x σ (3.6)
2
con los generadores en la respresentación vectorial (2.30). Por tanto, la parte de la trans-
formación dependiente de δx se lee
i ρ i
δx ρ ∂ρ φ(x) = − ωµν [J µν ]σ x σ ∂ρ φ(x) = ωµν L µν φ(x) (3.7)
2 2
donde hemos definido el operador

L µν = i (x µ ∂ν − x ν ∂µ ) (3.8)

Estos operadores cumplen las relaciones de conmutación(2.31) del álgebra de Lorentz y

por tanto constituyen una representación del grupo. Nótese que en este caso el espacio
base es un espacio de funciones y por tanto la representación es de dimensión infinita.
Puesto que podemos escribir L i j = x i p j − x j p i (i ∂i = p i ) y si definimos L i = ϵi j k L j k /2,
tenemos L i = ϵi j k x j p k . Es decir, podemos identificar L i con el momento angular orbital.
De nuevo las componentes K i = L i 0 corresponden a los generadores de los boosts.
Podemos por tanto escribir
i
δφ(x) = δ0 φ(x) + ωµν L µν φ(x) (3.9)
2
La parte δ0 φ(x) dependerá del tipo de campo que tengamos (escalar, vector, espinor, etc),
mientras que el segundo término es igual para todos ellos. Los distintos tipos de campos
se definen de la siguiente forma

Campos escalares
Se define un campo escalar como aquél que es invariante bajo transformaciones de Lorentz,
esto es

φ′ (x ′ ) = φ(x) (3.10)

y por tanto

δφ(x) = 0 (3.11)

En este caso, sustituyendo en (3.9) obtenemos

i
δ0 φ(x) = − ωµν L µν φ(x) (3.12)
2

Antonio L. Maroto, Ignazio Scimemi 24

Campos Cuánticos 2022-23 UCM

Figure 3.1: Transformación (δφ) de campos escalares y vectoriales bajo una rotación. La
región indicada corresponde a la zona donde el campo es no nulo. La rotación afecta de
dos formas distintas al campo: por una parte, cambia la zona donde el campo es no nulo y
por otra (en el caso vectorial) cambia la orientación del campo.

Campos vectoriales
Un campo vectorial es el que se transforma como
µ µ
V ′ (x ′ ) = Λ νV ν (x) (3.13)

de forma que

i ρ
δV ρ (x) = − ωµν [S µν ]σV σ (3.14)
2
donde hemos denotado ahora los generadores en la representación vectorial (2.30) como
S µν , es decir 1 ,
ρ µ
[S µν ] σ = i (ηµρ δνσ − ηνρ δ σ ) (3.15)

Por tanto, tenemos que

i ρ ρ
δ0V ρ (x) = − ωµν (L µν δ σ + [S µν ]σ )V σ (3.16)
2
1
Este cambio en la notación obedece a que ahora los generadores de las transformaciones δ0 V que pode-
mos denotar J tienen dos contribuciones, la parte de momento angular orbital L y la parte de momento
angular de spin S, es decir, podemos escribir δ0 V = − 2i ωµν J µν V = − 2i ωµν (L µν + S µν )V

Antonio L. Maroto, Ignazio Scimemi 25

Campos Cuánticos 2022-23 UCM

Campos de Weyl
Un campo de Weyl (levógiro) es aquél que bajo transformaciones de Lorentz cambia de la
siguiente forma

ψ′L (x ′ ) = ΛL ψL (x) (3.17)

donde ΛL viene dado por (2.45) que podemos escribir como

i µν
ΛL = e − 2 ωµν S (3.18)

donde los generadores son

1
S i j = ϵi j k σk (3.19)
2
i
S i 0 = σi (3.20)
2
Por tanto
i
δψL (x) = − ωµν S µν ψL (x) (3.21)
2
de donde
i
δ0 ψL (x) = − ωµν (L µν + S µν )ψL (x) (3.22)
2
En el caso dextrógiro, tendremos

1
S i j = ϵi j k σk (3.23)
2
i
S i 0 = − σi (3.24)
2
y el mismo L µν .

Campos de Dirac
Finalmente en el caso de Dirac, los campos se transforman de acuerdo con

ψ′D (x ′ ) = ΛD ψL (x) (3.25)

Antonio L. Maroto, Ignazio Scimemi 26

Campos Cuánticos 2022-23 UCM

donde ΛD está dado en (2.47) y por tanto

i
δψD (x) = − ωµν S µν ψD (x) (3.26)
2
donde

1 i j k σk 0
µ ¶
ij
S = ϵ (3.27)
2 0 σk
µ i
i σ
¶
i0 0
S = (3.28)
2 0 −σi

y de forma totalmente análoga al caso de Weyl

i
δ0 ψD (x) = − ωµν (L µν + S µν )ψD (x) (3.29)
2

Campos tensoriales
De forma totalmente análoga es posible construir campos tensoriales de orden arbitrario
con la correspondiente representación. Así por ejemplo un campo tensorial dos veces con-
travariante es aquel que cumple
µν µ
T′ (x ′ ) = Λ ρ Λνσ T ρσ (x) (3.30)

En resumen, vemos que en general, bajo transformaciones de Lorentz, un campo se

transforma como
i
δφ(x) = − ωµν S µν φ(x) (3.31)
2
i
δ0 φ(x) = − ωµν (L µν + S µν )φ(x) (3.32)
2
donde el operador de momento angular L µν es universal, independiente del tipo de campo,
y la parte de spin S µν depende de la representación concreta (escalar, vector, espinor, etc).
Nótese que en el caso escalar S µν = 0.
Una vez que hemos definido nuestros objetos básicos, los campos, y determinado sus
propiedades de transformación, vamos a construir acciones invariantes Lorentz que nos
permitan obtener sus ecuaciones de movimiento relativistas.

Antonio L. Maroto, Ignazio Scimemi 27

Campos Cuánticos 2022-23 UCM

3.2 Principio de mínima acción

El principio de mínima acción (o principio de Hamilton) se ha estudiado en Mecánica
Clásica y permite obtener las ecuaciones de movimiento de un sistema mecánico con un
número finito n de grados de libertad. Estas ecuaciones determinan las trayectorias del
sistema que hacen estacionaria la integral de acción. En este apartado veremos que es
posible extender este principio al caso de campos. Comenzaremos repasando el caso de la
Mecánica Clásica.
Consideremos por tanto un sistema mecánico con n coordendas generalizadas inde-
pendientes q i (t ), i = 1, . . . , n. El Lagrangiano del sistema L(q, q̇) en su forma más sencilla
tomará la forma
1X 2 2
L(q, q̇) = m i q̇ i − V (q)
2 i
y la acción se define como Z
S= d t L(q, q̇)

El principio de mínima acción establece que si fijamos los extremos de la trayectoria: q(t i n ) =
q i n y q(t f ) = q f , la trayectoria clásica del sistema es un extremo de la integral de acción, es
decir, cerca de la trayectoria clásica
Z tf
δ d t L(q, q̇) = 0
ti

La variación se realiza manteniendo fijos los extremos, es decir,

tf tf X ∂L ∂L
Z Z µ ¶
δS = d t δL = dt δq i + δq̇ i = 0
ti ti i ∂q i ∂q̇ i

d
Suponiendo que las variaciones y las derivadas conmutan δq̇ i = d t δq i e integrando por
partes tenemos Z tf X ∂L d ∂L
µ ¶
dt − δq i = 0
ti i ∂q i d t ∂q̇ i
donde hemos usado que δq i (t i n ) = δq i (t f ) = 0. Puesto que el resultado debe ser válido
∀δq, obtenemos finalmente las ecuaciones de Euler-Lagrange

∂L d ∂L
µ ¶
− = 0 i = 1, . . . , n
∂q i d t ∂q̇ i

Antonio L. Maroto, Ignazio Scimemi 28

Campos Cuánticos 2022-23 UCM

Un sistema de ecuaciones equivalente puede obtenerse con el formalismo Hamiltoniano.

El Hamiltoniano del sistema se define como
X
H (p, q) = p i q̇ i − L(q, q̇(p, q))
i

donde los momentos canónicos conjugados son

∂L
pi ≡
∂q̇ i

Hemos visto que un campo es una función de las coordenadas espacio-temporales,

por tanto podemos extender el formalismo Lagrangiano de funciones del tiempo q i (t ) a
funciones del espacio tiempo φi (x)

q i (t ) → φi (t ,⃗
x)

Podemos ahora definir una densidad Lagrangiana2

L (φ, ∂µ φ) (3.33)

y definimos el Lagrangiano del campo como

Z
L = d 3 x L (φ, ∂µ φ) (3.34)

La acción será por tanto

Z Z
S[φ] = dt L = d 4 x L (φ, ∂µ φ) (3.35)

donde la integración se realiza sobre todo el espacio-tiempo y las condiciones de contorno

implican que las variaciones de los campos δφi se anulan en el infinito. Podemos por
tanto calcular la variación de la integral de acción de forma totalmente análoga al caso de
la Mecánica Clásica
X ∂L ∂L
Z µ ¶
4
δS = d x δφi + δ(∂µ φi )
i ∂φi ∂(∂µ φi )
X ∂L ∂L
Z µ ¶
4
= d x − ∂µ δφi = 0 (3.36)
i ∂φi ∂(∂µ φi )
2
Consideraremos únicamente teorías locales en las que la densidad Lagrangiana depende del campo y de
sus derivadas en un único punto.

Antonio L. Maroto, Ignazio Scimemi 29

Campos Cuánticos 2022-23 UCM

y por tanto, las ecuaciones de Euler-Lagrange son

∂L ∂L
µ ¶
− ∂µ = 0, i = 1, . . . n (3.37)
∂φi ∂(∂µ φi )

Nótese que dos densidades Lagrangianas que difieren en una derivada total conducen
a las mismas ecuaciones de Euler-Lagrange. En efecto, consideremos

L ′ = L + ∂µ K µ con K µ = K µ (φ). (3.38)

En un volumen finito V limitado por una superficie Σ el teorema de Stokes implica

Z Z
4 µ
d x ∂µ K = d A n µ K µ (3.39)
V Σ

con d A el elemento de superficie y n µ la dirección normal saliente a la superficie. Si K µ

decrece suficientemente rápido en infinito este término va a cero en el límite de volumen
infinito y por tanto S ′ = S.
Siguiendo con la analogía con la Mecánica Clásica, podemos definir el momento canónico
conjugado del campo φi (x) como
∂L
πi (x) = (3.40)
∂(∂0 φi (x))
y la densidad Hamiltoniana

H = πi (x)∂0 φi (x) − L
X
(3.41)
i

El Hamiltoniano del campo será

Z
H= d 3x H (3.42)

3.3 El teorema de Noether

Este teorema relaciona las simetrías de la integral de acción con la existencia de leyes de
conservación. Consideremos una transformación infinitesimal que puede afectar tanto a
las coordenadas como a los campos, esto es

x µ → x ′µ = x µ + δx µ
φi (x) → φ′i (x ′ ) = φi (x) + δφi (x) (3.43)

Antonio L. Maroto, Ignazio Scimemi 30

Campos Cuánticos 2022-23 UCM

Se dice que estas transformaciones definen una simetría de la acción si dejan S invariante
para todo φi (x) (sea o no solución de las ecuaciones de movimiento).
Una observación importante es que una simetría requiere que la acción S sea invari-
ante, pero no necesariamente la densidad Lagrangiana L . Es decir, la condición de simetría
es S ′ [φ′ ] = S[φ] y por tanto
Z Z
δS[φ] = d x L (φ (x ), ∂µ φ (x )) − d 4 x L (φ(x), ∂µ φ(x)) = 0
4 ′ ′ ′ ′ ′
(3.44)
V′ V

donde V ′ denota el mismo volumen de integración V pero expresado en coordendas x ′ .

Haciendo el cambio de variable a las coordenadas x en la primera integral, tendremos que
usar el determinante jacobiano que a primer orden es

d 4 x ′ = d 4 x(1 + ∂µ δx µ ) (3.45)

y por tanto
Z
δS[φ] = d 4 x(1 + ∂µ δx µ ) L (φ′ (x), ∂µ φ′ (x)) + ∂µ L (φ(x), ∂µ φ(x)) δx µ
¡ ¢

ZV
− d 4 x L (φ(x), ∂µ φ(x)) (3.46)
V

Usando la definición (3.5)

φ′i (x) = φi (x) + δ0 φi (x) (3.47)

con

δ0 φi (x) = δφi (x) − δx λ ∂λ φi (x) (3.48)

podemos escribir

∂L ∂L
µ ¶
′ ′
L (φ (x), ∂µ φ (x)) = L (φ(x), ∂µ φ(x)) + δ0 φi + δ0 (∂µ φi ) (3.49)
∂φi ∂(∂µ φi )

Ahora usando δ0 (∂µ φi ) = ∂µ (δ0 φ) tenemos

∂L ∂L ∂L
µ
µ ¶¶ µ ¶
′ ′
L (φ (x), ∂µ φ (x)) = L (φ(x), ∂µ φ(x)) + − ∂µ δ0 φi + ∂µ δ0 φi
∂φi ∂(∂µ φi ) ∂(∂µ φi )
(3.50)

Antonio L. Maroto, Ignazio Scimemi 31

Campos Cuánticos 2022-23 UCM

Sustituyendo en (3.46) y quedándonos a primer orden obtenemos

Z
δS[φ] = d 4 x (∂µ δx µ ) L + (∂µ L ) δx µ
£

µV
∂L ∂L ∂L
µ ¶¶ µ ¶¸
+ − ∂µ δ0 φi + ∂µ δ0 φi = 0 (3.51)
∂φi ∂(∂µ φi ) ∂(∂µ φi )

que se puede reescribir como

∂L ∂L ∂L
Z ·µ µ ¶¶ µ ¶¸
4 µ
δS[φ] = d x − ∂µ δ0 φi + ∂µ δ0 φi + L δx =0 (3.52)
V ∂φi ∂(∂µ φi ) ∂(∂µ φi )

Puesto que esta variación debe ser cero para cualquier valor del volumen de integración V ,
el integrando debe ser cero por lo que
∂L ∂L ∂L
µ µ ¶¶ µ ¶
µ
− ∂µ δ0 φi + ∂µ δ0 φi + L δx = 0 (3.53)
∂φi ∂(∂µ φi ) ∂(∂µ φi )

Ahora bien, el primer paréntesis no son más que las ecuaciones de Euler-Lagrange (3.37),
y por tanto sobre soluciones de las ecuaciones clásicas de movimiento φi (x) = φcl
i
(x) lleg-
amos a

∂µ j µ (φcl ) = 0 (3.54)

donde
∂L
µ ¶
µ
j = δ0 φi + L δx µ (3.55)
∂(∂µ φi )

Supongamos ahora que las transformaciones dependen de una serie de parámetros

a
θ (x) con a = 1, . . . , N de forma que
µ
δx µ = θ a (x)X a (x)
δφi (x) = θ a (x)Yi ,a (φ, ∂µ φ) (3.56)

Estas transformaciones pueden corresponder a cada uno de los parámetros de un grupo

de Lie. En este caso, (3.48) se escribirá

δ0 φi (x) = θ a (x)(Yi ,a (x) − X aλ (x)∂λ φi (x)) (3.57)

y por tanto la corriente (3.55) se escribirá

∂L
µ ¶
µ λ µ
j = (Yi ,a (x) − X a (x)∂λ φi (x)) + L X a (x) θ a (x) (3.58)
∂(∂µ φi )

Antonio L. Maroto, Ignazio Scimemi 32

Campos Cuánticos 2022-23 UCM

Puesto que los parámetros θ a son independientes, podemos definir una corriente conser-
vada para cada uno de ellos:

µ ∂L µ
ja = (Yi ,a (x) − X aλ (x)∂λ φi (x)) + L X a (x), a = 1, . . . , N (3.59)
∂(∂µ φi )

de forma que sobre soluciones de las ecuaciones del movimiento

µ
∂µ j a (φcl ) = 0, a = 1, . . . , N (3.60)

Es decir, si la acción es invariante bajo las transformaciones de simetría (3.56) exis-

µ
ten N corrientes de Noether j a conservadas cuando se evalúan sobre las soluciones clási-
cas. Es decir, hay una corriente conservada por cada generador de las transformaciones de
simetría. Este es el Teorema de Noether.
Para cada corriente conservada tenemos una carga conservada dada por:
Z
Q a ≡ d 3 x j a0 (x) (3.61)

Efectivamente, usando (3.60) tenemos

Z Z
∂0Q a = d x ∂0 j a (x) = − d 3 x ∂i j ai (x) = 0
3 0
(3.62)

donde en la última integral hemos usado el teorema de Stokes y asumido que los campos
se anulan en el infinito.

Clasificación de las simetrías

Dependiendo de la forma de las transformaciones (3.56), podemos distinguir distintos
tipos de simetrías:

Simetrías globales: los parámetros θ a son independientes de las coordenadas; las simetrías
globales son un caso especial de las simetrías locales.

Simetrías locales: los parámetros θ a (x) dependen de las coordenadas

µ
Simetrías internas: X a (x) = 0, es decir, las coordenadas no cambian.
µ
Simetrías espacio-temporales: X a (x) ̸= 0, es decir, involucran cambios en las coorde-
nadas.

Antonio L. Maroto, Ignazio Scimemi 33

Campos Cuánticos 2022-23 UCM

3.4 Invariancia bajo traslaciones: tensor energía-momento

El ejemplo más sencillo es la simetría bajo traslaciones, dadas en el caso infinitesimal por

x ′µ = x µ + θ µ
φ′i (x ′ ) = φi (x) (3.63)

En este caso, el índice a = 1, . . . , 4 podemos sustituirlo por un índice griego ν, de forma que
µ µ
Xν = δ ν
Yi ,ν = 0 (3.64)

La corriente conservada se lee en este caso

∂L
θ µν = ∂ν φi − ηµν L . (3.65)
∂(∂µ φi )
µ
donde hemos usado θ µν = ηνρ θ ρ . Esta corriente conservada se conoce como tensor energía-
momento.
La conservación de la corriente implica

∂µ θ µν (φcl ) = 0 (3.66)

y las cargas son

Z
ν
P ≡ d 3 x θ 0ν (φcl ) (3.67)

que es el cuadrimomento del campo.

Nótese que θ µν no es simétrico en sus índices. Sin embargo, podemos construir un
tensor simétrico

T µν = θ µν + ∂ρ A ρµν (3.68)

con A un tensor antisimétrico en los índices ρ, µ elegido adecuadamente. Precisamente

debido a la antisimetría de A, el nuevo tensor también es conservado

∂µ T µν = ∂µ ∂ρ A ρµν = 0.

Para µ = 0 tenemos ∂ρ A ρ0ν = ∂i A i 0ν que es una divergencia espacial y por tanto no cambia
las cargas conservadas en (3.67). En otras palabras, θ µν y T µν son físicamente equivalentes
y podemos trabajar con cualquiera de ellos. T µν se conoce como tensor energía-momento
simetrizado de Belinfante-Rosenfeld.

Antonio L. Maroto, Ignazio Scimemi 34

Campos Cuánticos 2022-23 UCM

3.5 Invariancia Lorentz

Consideremos ahora el caso en el que la acción es invariante bajo transformaciones de
Lorentz dadas anteriormente en (3.6) y (3.31)
µ
x ′µ = x µ + ω ν x ν
i
φ′i (x ′ ) = φi (x) − ωµν [S µν ] i φ j (x)
j
(3.69)
2
de forma que ahora el índice a se convierte en un par de índices griegos ν, ρ

1
X µνρ = (ηµν x ρ − ηµρ x ν )
2
νρ i
Yi = − [S νρ ] i φ j (x)
j
(3.70)
2
j
donde los generadores [S νρ ] i se escribirán en la representación correspondiente a los
campos φi .
La correspondiente corriente de Noether será

∂L i νρ j
· µ ¶ ¸
µ λνρ µνρ
j = − [S ] i φ j (x) − X (x)∂λ φi (x) + L X (x) ωνρ (3.71)
∂(∂µ φi ) 2

Esta corriente se puede escribir como

1
j µ = M µνρ ωνρ (3.72)
2
donde M µνρ se puede escribir en función del tensor energía-momento como

M µνρ = x ν θ µρ − x ρ θ µν + S µνρ (3.73)

con
∂L
S µνρ = −i [S νρ ] i φ j
j
(3.74)
∂(∂µ φi )

La conservación de la corriente implica

∂µ M µνρ = 0 (3.75)

Antonio L. Maroto, Ignazio Scimemi 35

Campos Cuánticos 2022-23 UCM

Por tanto, tenemos 6 corrientes conservadas, M µνρ una por cada uno de los 6 parámetros
independientes de la transformación de Lorentz ωνρ . Nótese que por construcción M µνρ
es antisimétrico en (νρ). Las correspondientes cargas conservadas son:
Z Z
νρ
M = d xM 3 0νρ
= d 3 x(x ν p ρ − x ρ p ν + S 0νρ ) (3.76)

con p ν = θ 0ν la densidad de cuadrimomento.

Nótese que las 3 cargas conservadas M i j , asociadas a las rotaciones, proporcionan el
momento angular del campo como

1
J i = ϵi j k M j k (3.77)
2

que como vemos tienen la contribución del momento angular orbital (x i p j − x j p i ) y de

spin S 0i j . Las M 0i corresponderían a las 3 cargas conservadas asociadas a los boosts. 3

3
La interpretación física de estas cargas conservados no es tan directa y puede verse que corresponden a
la posición del centro de momentos del campo en el instante t = 0.

Antonio L. Maroto, Ignazio Scimemi 36

4. Cuantización de campos escalares libres

En este tema obtendremos la acción clásica invariante Lorentz para un campo escalar
(tanto real como complejo) y extenderemos los métodos de cuantización de la Mecánica
Cuántica al caso de campos.

4.1 El campo escalar real clásico

Para construir una teoría relativista para un campo escalar real φ(x), debemos partir de
una acción invariante bajo transformaciones de Lorentz. Una acción que describa una
dinámica no trivial debe contener términos con derivadas de los campos ∂µ φ. Para poder
tener una densidad Lagrangian invariante Lorentz el índice µ debe aparecer contraído para
formar un escalar. La forma más simple es contraerlo con otro factor ∂µ φ. Por otra parte,
la medida de integración d 4 x es invariante Lorentz gracias a la condición det Λ = 1. Por
tanto, la acción más simple es de la forma1
µ ¶
1 µν 1 2 2
Z
4
S= d x η ∂µ φ∂ν φ − m φ (4.1)
2 2
donde también hemos incluido un término de energía potencial que corresponde al lla-
mado término de masa. Las ecuaciones de movimiento son
∂L ∂L
µ ¶
− ∂µ = 0,
∂φ ∂(∂µ φ)
−m 2 φ − ∂µ (ηµν ∂ν φ) = 0 (4.2)

que proporciona la ecuación de Klein-Gordon

(□ + m 2 )φ = 0 (4.3)
1
En unidades naturales la acción S[φ] debe ser adimensional, por lo que la dimensión del campo escalar
debe ser [φ] = E

37
Campos Cuánticos 2022-23 UCM

Una solución a esta ecuación es la onda plana e ±i px con p 2 = p 02 −⃗ p 2 = m 2 que es la llamada

relación de dispersión relativista con m la masa del campo. Teniendo en cuenta que el
campo φ debe ser real la solución más general se puede escribir como
d 3p
Z ¯
−i px ∗ i px ¯
φ(x) = (a p⃗ e + a p⃗ e )¯ (4.4)
(2π)3 2E p
p
p 0 =E p

donde los coeficientes a p⃗ (complejos en general) dependen del momento p ⃗ y el factor

3
p
1/((2π) 2E p ) se introduce por conveniencia. Nótese, como ya vimos en la Introducción,
que aunque p 0 = E p > 0, siempre tenemos modos de energía positiva y negativa en la
solución (esto es modos con e −i p 0 t y e i p 0 t respectivamente).
El momento conjugado de φ es
∂L
πφ = = ∂0 φ (4.5)
∂(∂0 φ)
y la densidad Hamiltoniana
1
H = (π2φ + (⃗
∇φ)2 + m 2 φ2 ) (4.6)
2
El tensor energía-momento es

θ µν = ∂µ φ∂ν φ − ηµν L (4.7)

El Hamiltoniano viene dado por

Z Z
3 00
H= d xθ = d 3 xH (4.8)

y el trimomento es por tanto

Z Z
i
P = 3
d xθ = 0i
d 3 x πφ ∂i φ (4.9)

Por otra parte las corrientes de Noether asociadas a la invariancia Lorentz en esta caso
corresponden a:

M µνρ = x ν θ µρ − x ρ θ µν (4.10)

puesto que en el caso escalar S µνρ = 0, los campos escalares representan partículas de
spin cero y al momento angular total del campo sólo contribuye el momento angular or-
bital. Las cargas conservadas son por tanto
Z Z
νρ
M = d xM 3 0νρ
= d 3 x(x ν p ρ − x ρ p ν ) (4.11)

Antonio L. Maroto, Ignazio Scimemi 38

Campos Cuánticos 2022-23 UCM

y el momento angular del campo escalar es

1
Z
J i = ϵi j k M j k → ⃗
J =− x ×⃗
d 3 x πφ (⃗ ∇φ) (4.12)
2
El Lagrangiano que hemos introducido sólo incluye un término cinético y el término
de masa. Es posible incluir autointeracciones a través de un término de potencial V (φ)
adicional
1 m2 2
L = ∂µ φ∂µ φ − φ − V (φ). (4.13)
2 2

4.2 Cuantización en Mecánica Cuántica

Para cuantizar el campo vamos a extender la cuantización canónica de sistemas mecáni-
cos. Consideremos un sistema mecánico con n grados de libertad con coordenadas q i y
momentos canónicos conjugados p i = ∂L/∂q̇ i , i = 1..n. El Hamiltoniano como vimos está
dado por X
H (p, q) = p i q̇ i − L(q, q̇).
i
Para cuantizar el sistema las variables dinámicas q, p, H , etc se promocionan a operadores
hermíticos q̂, p̂, Ĥ que actúan en el espacio de Hilbert de los estados |Ψ〉 ∈ H . El proceso
de cuantización canónica impone que los operadores asociados a las variables canónicas
conjugadas satisfacen las relaciones de conmutación

[q̂ i , p̂ j ] = i δi j (4.14)

[q̂ i , q̂ j ] = 0 (4.15)

[p̂ i , p̂ j ] = 0 (4.16)

La evolución del estado del sistema viene dada en la imagen de Schrödinger por

∂ ¡ ¢
iℏ |Ψ (t )〉 = Ĥ p̂, q̂ |Ψ (t )〉 (4.17)
∂t
El estado del sistema en el instante t , |Ψ (t )〉 viene univocamente determinado por el es-
tado inicial |Ψ (t 0 )〉 de forma que

|Ψ (t )〉 = U (t , t 0 )|Ψ (t 0 )〉 (4.18)

Antonio L. Maroto, Ignazio Scimemi 39

Campos Cuánticos 2022-23 UCM

donde el operarador de evolución U (t , t 0 ) es un operador unitario que en el caso en el que

el operador Hamiltoniano no dependa explícitamente del tiempo viene dado por

Û (t , t 0 ) = e −i Ĥ (t −t0 ) (4.19)

En esta imagen los estados evolucionan en el tiempo pero los observables no lo hacen.
También es posible describir la dinámica del sistema en la imagen de Heisenberg en la
que el estado del sistema se fija al estado inicial |Ψ (t 0 )〉 mientras que toda la evolución
está en los operadores que vienen dados por

 H (t ) = Û † (t , t 0 )  S Û (t , t 0 ) (4.20)

donde S y H denotan la imagen de Schrödinger y Heisenberg respectivamente. La ecuación

de evolución de los operadores en la imagen de Heisenberg es (suponiendo ∂ Â S /∂t = 0)

∂ Â H £ ¤
i = Â H , Ĥ H (4.21)
∂t
De esta forma los valores esperados obtenidos en ambas imágenes son iguales

A (t ) = 〈Ψ (t ) ¯ Â S ¯ Ψ (t )〉 = 〈Ψ (t 0 ) ¯ Â H (t )¯ Ψ (t 0 )〉
¯ ¯ ¯ ¯
(4.22)

Nótese que en la imagen de Heisenberg las relaciones de conmutación se traducen en rela-

ciones de conmutación a tiempos iguales, esto es

[q̂ i ,H (t ), p̂ j ,H (t )] = i δi j (4.23)

4.3 Cuantización del campo escalar real

Consideremos por tanto el Lagrangiano del campo escalar real con masa libre (sin interac-
ciones)
1 1
L = ηµν ∂µ φ∂ν φ − m 2 φ2 (4.24)
2 2
De forma que como hemos visto en (4.5), el momento canónico conjugado del campo
φ(t ,⃗
x ) es

∂L
πφ (t ,⃗
x) = = ∂0 φ (4.25)
∂(∂0 φ)

Antonio L. Maroto, Ignazio Scimemi 40

Campos Cuánticos 2022-23 UCM

Podemos por tanto cuantizar el campo escalar por analogía con la cuantización en ima-
gen de Heisenberg. Así podríamos identificar el campo φ(t ,⃗ x ) con los operadores q̂ i ,H (t ),
de forma que el índice discreto i pasaría a convertirse en el vector posición ⃗ x continuo y
de forma análoga p̂ i ,H (t ) → πφ (t ,⃗
x ). De esta forma las reglas de cuantización canónica
supone promover φ y π a operadores e imponer la regla de conmutación a tiempos iguales

φ (t ,⃗
x ) , πφ t ,⃗ = i δ(3) ⃗
£ ¡ ¢¤ ¡ ¢
y x −⃗
y (4.26)
φ (t ,⃗
x ) , φ t ,⃗
£ ¡ ¢¤
y = 0 (4.27)
πφ (t ,⃗x ) , πφ t ,⃗
£ ¡ ¢¤
y = 0 (4.28)

Nótese que la delta de Kronecker ahora se ha convertido en delta de Dirac por el carácter
continuo de la etiqueta ⃗
x . El campo φ(t ,⃗
x ) se convierte ahora en un operador hermítico,
de forma que la solución general de la ecuación de Klein-Gordon se escribe 2

d 3p
Z
φ(x) = (a p⃗ e −i px + a p⃗† e i px )|p 0 =E p (4.29)
(2π)3 2E p
p

⃗2 + m 2 y donde a p⃗ se promueve a operador y a p⃗∗ se convierte en el operador

p
con E p = + p
adjunto a p⃗† de forma que φ es hermítico. De forma totalmente análoga tenemos
s
3 Ep
d p
Z
πφ (x) = ∂0 φ(x) = −i (a p e −i px − a p† e i px )|p 0 =E p (4.30)
(2π)3 2

Nótese que la única evolución temporal aparece a través de las exponenciales e ±i p 0 t de

forma que los operadores a p y a p† son independientes del tiempo y por tanto son los oper-
adores en la imagen de Schrödinger.
Las relaciones de conmutación anteriores son equivalentes a las siguientes relaciones
para los operadores a y a †
h i
a p , a q† = (2π)3 δ(3) p⃗ −⃗
¡ ¢
q (4.31)
£ ¤
ap , aq = 0 (4.32)
h i
a p† , a q† = 0 (4.33)

2
p
La normalización 1/ 2E p puede diferir respecto a algunos libros que usan 1/(2E p )

Antonio L. Maroto, Ignazio Scimemi 41

Campos Cuánticos 2022-23 UCM

En efecto, podemos ver que tomando t = 0 sin pérdida de generalidad podemos escribir3

d 3p
Z
φ(0,⃗
x) = †
(a p + a −p )e i p⃗·⃗x (4.35)
3
p
(2π) 2E p

s
3 Eq
d q
Z
πφ (0,⃗
y ) = −i †
(a q − a −q )e i ⃗q ·⃗y (4.36)
(2π)3 2

Por tanto
s
−i d 3 pd 3 q E q ³h
Z i h i´
φ (0,⃗
x ) , π 0,⃗ † +
e i (p⃗·⃗x +⃗q ·⃗y )
£ ¡ ¢¤
y = a −p , a q − a p , a −q
2 (2π)6 Ep
s
−i d 3 pd 3 q E q ¡
Z
3 (3) 3 (3)
¢¢ i (p⃗·⃗x +⃗q ·⃗y )
δ ⃗ ⃗ δ ⃗ ⃗
¡ ¢ ¡
= − (2π) q + p − (2π) p + q e
2 (2π)6 Ep
d 3 p i p⃗·(⃗x −⃗y )
Z
= i e = i δ(3) (⃗
x −⃗y) (4.37)
(2π)3

Podemos ahora obtener la expresión para el operador Hamiltoniano (4.6) en térmi-

nos de los operadores creación y destrucción (la expresión coincide tanto en la imagen
de Schrödinger como Heisenberg) por lo que calcularemos de nuevo con t = 0.
· ¸
1 2 1 1 2 2
Z Z
H = d x H = d x πφ + (⃗
3 3 2
∇φ) + m φ
2 2 2
Z 3 3 Ã p
d pd q i (p⃗+⃗q )·⃗x Ep Eq ³
Z ´³ ´
= d 3x e − a p − a †
−p a q − a †
−q
(2π)6 4
!
p ·⃗
−⃗ q + m2 ³ †
´³
†
´
+ p a p + a −p a q + a −q (4.38)
4 Ep Eq

p +⃗
Realizando ahora la integral en x obtenemos una fator(2π)3 δ(3) (⃗ q ) que nos permite
3 ⃗y
Donde hemos usado que E p es una función par en p
Z ∞ Z ∞
d pe i px = d pe −i px (4.34)
−∞ −∞

Antonio L. Maroto, Ignazio Scimemi 42

Campos Cuánticos 2022-23 UCM

integrar en q obteniendo
d 3p Ep ³ ´ p ⃗2 + m 2 ³
Z µ ´³ ´³ ´¶
† † † †
H = − a p − a −p a −p − a p + a p + a −p a −p + a p
(2π)3 4 4E p
d 3 p E p h³
Z ´³ ´ ³ ´³ ´i
† † † †
= a p + a −p a −p + a p − a p − a −p a −p − a p
(2π)3 4
d 3p Ep ³
Z ´
† †
= a a
p p + a a
p p (4.39)
(2π)3 2
Podemos ahora reescribir este expresión como
d 3p
· ¸
1
Z
† †
H= E p a p a p + [a p , a p ] (4.40)
(2π)3 2
El último término es proporcional a δ(3) (0) y es por tanto divergente. Esta divergencia es
esperable pues corresponde a la suma de la energía de vacío de los infinitos osciladores
armónicos. Este es el primer ejemplo de la aparición de divergencias en teoría cuántica
de campos. Estas divergencias requieren de un proceso de renormalización para poder
eliminarlas. En espacio-tiempo plano, la renormalización de la energía de vacío puede
realizarse directamente eliminando este último término, dado que los procesos físicos sólo
son sensibles a diferencias de energías pero no a sus valores absolutos. En presencia de
campos gravitatorios, este proceso no puede llevarse a cabo de forma directa, puesto que
la densidad de energía es fuente de campo gravitatorio. Esto plantea uno de los problemas
fundamentales de la Teoría Cuántica de Campos que es la renormalización de la energía
de vacío en espacio-tiempo curvo.
Es posible formalizar este procedimiento de eliminación de la energía de vacío uti-
lizando la llamada ordenación normal de operadores. Dado un operador A denotamos
su forma normalmente ordenada como : A : que se obtiene colocando todos lo operadores
de creación a la izquierda y los de destrucción a la derecha. Así, por ejemplo : a p a p† := a p† a p .
Podemos también obtener la expresión para el operador momento total,
d 3p i †
Z Z Z
i 3 0i 3
P = d x θ = − d x πφ ∂i φ = p ap ap (4.41)
(2π)3
En este caso no es necesaria la ordenación normal, porque el término del vacío se cancela
al integrar en d 3 p. Como puede verse inmediatamente H y P ⃗ conmutan, por lo que confir-
mamos que tanto P ⃗ como H tienen las mismas expresiones en la imagen de Schrödinger o
Heisenberg.
De forma totalmente análoga al oscilador armónico podemos definir el operador número
de partículas con momento p ⃗ como
N p = a p† a p (4.42)

Antonio L. Maroto, Ignazio Scimemi 43

Campos Cuánticos 2022-23 UCM

4.3.1 Espacio de Fock

Las relaciones de conmutación (4.31), (4.32) y (4.33) para cada modo de Fourier coinciden
con las de los operadores de creación y destrucción del oscilador armónico en Mecánica
Cuántica. Podemos por tanto definir el estado de vacío |0〉 como el estado que es aniquilado
por todos los a p , esto es

a p⃗ |0〉 = 0 ∀⃗
p (4.43)

Normalizamos el estado de vacío de forma que 〈0|0〉 = 1. El resto de estados del espacio de
Fock pueden obtenerse actuando sobre el vacío con los distintos operadores de creación.
Los estados de una partícula se escriben por tanto
¡ ¢1/2 †
p 〉 = 2E p
|⃗ a p |0〉 (4.44)

Los estados de varias partículas serán

¢1/2 ¡ ¢1/2 † †
⃗2 , . . . p
⃗n 〉 = 2E p 1 a p 1 a p 2 . . . a p† n |0〉
¡
p1, p
|⃗ . . . 2E p n (4.45)
¡ ¢1/2
Los factores 2E p se introducen por conveniencia.
Usando las definiciones anteriores y las reglas de conmutación obtenemos

p 2 〉 = 〈p 1 |p 2 〉〈0|a p 1 a p† 2 |0〉
p 1 |⃗
〈⃗
¢1/2 ¡ ¢1/2
〈0|[a p 1 , a p† 2 ]|0〉
¡
= 2E p 1 2E P 2
= 2E p 1 (2π)3 δ(3) (⃗ ⃗2 )
p1 − p (4.46)

La combinación E p 1 δ(3) (⃗ p1 − p ⃗2 ) puede demostrarse que es invarante Lorentz. Nótese que

estados con distinto número de partículas son ortogonales entre sí. Por ejemplo
¢1/2 ¡ ¢1/2 ¡ ¢1/2
⃗3 〉 = 2E p 1 〈0|a p 1 a p† 2 a p† 3 |0〉
¡
p 1 |⃗
〈⃗ p2, p 2E P 2 2E p 3
¢1/2 ¡ ¢1/2 ¡ ¢1/2 ³ ´
〈0|a p† 2 a p 1 a p† 3 |0〉 − (2π)3 δ(3) (⃗ ⃗2 )〈0|a p† 3 |0〉 = 0
¡
= 2E p 1 2E P 2 2E p 3 p1 − p
(4.47)

El espacio de Fock estará formado por tanto por la suma directa de los espacios de Hilbert
de una partícula, dos partículas, etc.
Utilizando la expresión de H en términos de los operadores de creación y destrucción
es posible calcular los conmutadores
h i
H , a p† = E p a p† (4.48)
£ ¤
H , a p = −E p a p (4.49)

Antonio L. Maroto, Ignazio Scimemi 44

Campos Cuánticos 2022-23 UCM

Con lo que podemos obtener fácilmente la energía de un estado genérico del espacio de
Fock

H |⃗ ⃗2 , . . . p
p1, p ⃗n 〉 = (E p 1 + E p 2 + . . . E p n )|⃗ ⃗2 , . . . p
p1, p ⃗n 〉 (4.50)

de forma totalmente análoga obtenemos

⃗ |⃗
P ⃗2 , . . . p
p1, p ⃗n 〉 = (⃗ ⃗2 + . . . p
p1 + p ⃗n )|⃗ ⃗2 , . . . p
p1, p ⃗n 〉 (4.51)

Resulta por tanto natural denotar estas excitaciones del campo como partículas puesto
que son excitaciones discretas con la energía y el momento correcto. Así el estado |⃗ ⃗2 , . . . p
p1, p ⃗n 〉
correspondería a un conjunto de partículas con energías E p 1 , E p 2 , . . . E p n y momentos
⃗1 , p
p ⃗2 , . . . p
⃗n . Nótese que la expresión del operador campo (4.29) contiene tanto modos
de frecuencia positiva (multiplicados por operadores de destrucción) como modos de fre-
cuencia negativa (multiplicados por operadores de creación) aunque como hemos visto los
estados de Fock tienen energía positiva. Es decir, en la interpretación como excitaciones
de campos los modos de frecuencia negativa no implican la aparición de partículas con
energías negativas. También es importante notar que estas partículas no están localizadas
en el espacio sino que a p† crea autoestados del momento.
Este formalismo nos permite determinar la estadística de nuestras partículas. Consid-
¢1/2 ¡ ¢1/2 † †
⃗2 〉 = 2E p 1 a p 1 a p 2 |0〉. Puesto que a p† 1 y
¡
eremos el estado de dos partículas |⃗ p1, p 2E p 2
a p† 2 conmutan, este estado es idéntico al estado
¢1/2 ¡ ¢1/2 † †
⃗1 〉 = 2E p 2
¡
p2, p
|⃗ 2E p 1 a p 2 a p 1 |0〉

De hecho un modo p ⃗ puede contener un número arbitrariamente alto de partículas. Por

tanto concluimos que las partículas de Klein-Gordon satisfacen la estadística de Bose.
Consideremos ahora la interpretación del estado φ(0,⃗ x )|0〉. Actúando con (4.29) sobre
el vacío tenemos

d 3p d 3p
Z Z
⃗·⃗
φ(0,⃗
x )|0〉 = e −i p x †
a p |0〉 = 3
e −i p⃗·⃗x |⃗
p〉 (4.52)
3
p
(2π) 2E p (2π) (2E p )

Nótese que salvo por el factor 1/(2E p ) que luego estudiaremos, esta es la misma expre-
x 〉 en términos de autoestados de
sión que la del desarrollo del autoestado de posición |⃗
momento en Mecánica Cuántica. Por tanto podemos interpretar que el operador φ(0,⃗ x)
actuando sobre el vacío crea una partícula en la posición ⃗ x . Podemos confirmar esta in-

Antonio L. Maroto, Ignazio Scimemi 45

Campos Cuánticos 2022-23 UCM

terpretación calculando directamente

d 3p
Z
x )|⃗
〈0|φ(t ,⃗ q 〉 = 〈0| (a p e −i (E p t −⃗p ·⃗x ) + a p† e i (E p t −⃗p ·⃗x ) )|⃗
q〉
3
p
(2π) 2E p
s
d 3p Eq
Z
= 〈0| 3
(a p e −i (E p t −⃗p ·⃗x ) + a p† e i (E p t −⃗p ·⃗x ) )a q† |0〉
(2π) Ep
s
d 3p E q −i (E t −⃗p ·⃗x )
Z
= 3
e p
〈0|a p a q† |0〉
(2π) Ep
s
d 3p E q −i (E t −⃗p ·⃗x )
Z
= 3
e p
〈0|[a p , a q† ]|0〉 = e −i (E q t −⃗q ·⃗x ) (4.53)
(2π) Ep

q 〉 ∝ e i ⃗q ·⃗x , vemos que

x |⃗
de forma totalmente análoga al resultado en Mecánica Cuántica 〈⃗
incluyendo la dependencia temporal obtenemos la función de onda de una partícula de
momento ⃗ q y energía positiva.

4.4 Causalidad
A diferencia de la Mecánica Cuántica (no relativista), en el caso relativista puesto que ninguna
señal puede propagarse a una velocidad mayor que la velocidad de la luz, se requiere que
dos medidas realizadas en regiones cuasalmente desconectadas del espacio-tiempo no
puedan tener ningún tipo de relación, es decir, deben ser independientes. En términos
matemáticos, esto significa que cualquier par de observables A(x) y B (y) asociados a me-
didas realizadas en puntos del espacio-tiempo separados espacialmente deben conmutar,
esto es

[A(x), B (y)] = 0, (x − y)2 < 0 (4.54)

Consideremos por tanto el caso más sencillo en el que los operadores A y B se identifican
con el propio campo

d 3p d 3q
Z Z h i
−i px † i px −i q y † iqy
[φ(x), φ(y)] = ap e + ap e , aq e + aq e
(2π)3 2E p (2π)3 2E q
p p

d 3p
Z ³ ´
−i p(x−y) i p(x−y)
= e − e (4.55)
(2π)3 (2E p )

En primer lugar, vemos que el conmutador de los campos libres es un número no un oper-
ador (c-número, que conmuta con cualquier operador). Nótese ahora que tanto la medida

Antonio L. Maroto, Ignazio Scimemi 46

Campos Cuánticos 2022-23 UCM

de integración como las exponenciales son invariantes Lorentz y por tanto la expresión
anterior también es invariante Lorentz. Ahora bien, para sucesos con separación espacial
tenemos (x − y)2 < 0 y siempre es posible encontrar un sistema de referencia en el que los
dos sucesos sean simultáneos x 0 = y 0 y por tanto, dado que el conmutador es invariante
Lorentz, podemos evaluar la integral en este sistema

d 3p
Z ³ ´
⃗(⃗
ip x −⃗
y) ⃗(⃗
−i p x −⃗
y)
[φ(x), φ(y)] = e − e (4.56)
(2π)3 (2E p )

Usando (4.34) las integrales de los dos términos se cancelan y por tanto

[φ(x), φ(y)] = 0, (x − y)2 < 0 (4.57)

Nótese que esta conclusión sólo es válida para puntos con separación espacial y que para
llegar a este resultado es necesaria la presencia de modos de frecuencia positiva y negativa
en el desarrollo de los campos.

4.5 El propagador de Klein-Gordon

Consideremos el conmutador [φ(x), φ(y))] en más detalle en el caso general cuando la sep-
aración de los sucesos no es necesariamente de género espacio. La integral en (4.55) en
general no está bien definida y es necesario dar una prescripción para poder dar sentido
a esta expresión. Puesto que el conmutador es un c-número podemos escribirlo también
como [φ(x), φ(y))] = 〈0|[φ(x), φ(y))]|0〉, es decir

d 3p
Z ³ ´
−i p(x−y) i p(x−y)
〈0|[φ(x), φ(y))]|0〉 = e − e
(2π)3 (2E p )
à !
d 3p
¯ ¯
1 −i p (x−y ) ¯¯ 1
Z
e −i p (x−y ) ¯¯
¯
= e +
(2π)3 2E p ¯
p 0 =E
−2E p p 0 =−E p
3 0
d p dp −1
Z Z
= 3 2 2
e −i p (x−y ) , (x 0 > y 0 ) (4.58)
(2π) 2πi p − m

donde en el último paso hemos supuesto que x 0 > y 0 y hemos usado integración por resid-
uos en la variable p 0 con el contorno en Fig. 4.1 cerrado por abajo por lo que aparece el
signo menos global. Hemos desplazado los polos en p 0 = ±E p del eje real añadiendo −i ϵ
de forma que la integral por residuos puede hacerse. El resultado final se obtiene tomando
ϵ → 0.

Antonio L. Maroto, Ignazio Scimemi 47

Campos Cuánticos 2022-23 UCM

Figure 4.1: Propagador retardado: contorno de integración. En el caso x 0 > y 0 el contorno

se cierra por abajo con lo que los polos se recorren en sentido horario por lo que aparece
un signo menos en la integral. En el caso x 0 < y 0 se cierran por la parte superior.

En el caso x 0 < y 0 podemos cerrar el contorno por la parte superior por lo que la inte-
gral es cero. Por tanto la integral (4.58) junto con la prescripción anterior para recorrer los
polos nos permiten definir la función que denotaremos como

D R x − y = θ x 0 − y 0 〈0 ¯ φ (x) , φ y ¯ 0〉
¡ ¢ ¡ ¢ ¯£ ¡ ¢¤¯
(4.59)

Podemos por tanto escribir esta función de forma más compacta como

d 4p i
Z
D R (x − y) = 4 p 2 − m2
e −i p (x−y ) (4.60)
(2π)
con la prescripción sobre los contornos que hemos descrito antes. Si actuamos ahora con
el operador de Klein-Gordon sobre esta función

d 4p
Z
2
e −i p (x−y ) = −i δ(4) x − y
¡ ¢ ¡ ¢ ¡ ¢
□+m x D R x − y = −i 4
(4.61)
(2π)

Antonio L. Maroto, Ignazio Scimemi 48

Campos Cuánticos 2022-23 UCM

Figure 4.2: Propagador de Feynman: contorno de integración. En el caso x 0 > y 0 el con-

torno se cierra por abajo con lo que el polo +E p se recorren en sentido horario por lo que
aparece un signo menos en la integral. En el caso x 0 < y 0 se cierran por la parte superior y
el polo −E p se recorre en sentido antihorario y tiene signo positivo.

Vemos por tanto que D R (x − y) es una función de Green del operador de Klein-Gordon.
Puesto que se anula para x 0 < y 0 se denomina la función de Green retardada o propagador
retardado.
Sin embargo es posible utilizar cuatro prescripciones distintas para realizar la integral
por residuos. En concreto una particularmente importante consiste en utilizar el contorno
de la Fig.4.2. De nuevo si x 0 > y 0 el contorno se cierra por abajo con lo que el polo +E p se
recorren en sentido horario por lo que aparece un signo menos en la integral. En el caso
x 0 < y 0 se cierran por la parte superior y el polo −E p se recorre en sentido antihorario y
tiene signo positivo. Una forma de implementar esta prescripición en la propia integral es
escribirla de la forma
d 4p i
Z
∆(x − y) = 4 p 2 − m2 + i ϵ
e −i p (x−y ) (4.62)
(2π)
puesto que ahora los polos se encuentran p 0 = ±(E p − i ϵ), de forma que el polo −E p se

Antonio L. Maroto, Ignazio Scimemi 49

Campos Cuánticos 2022-23 UCM

desplaza hacia arriba y el E p hacia abajo del eje real. De esta forma tenemos

d3p
 R

(2π)3 (2E p )
e −i p(x−y) , x0 > y 0
∆(x − y) = R d3p
(4.63)

(2π)3 (2E p )
e −i p(y−x) , x0 < y 0

pero

d 3p
Z
〈0 ¯φ (x) φ y ¯ 0〉 = e −i p(x−y)
¯ ¡ ¢¯
(4.64)
(2π)3 (2E p )

luego podemos escribir

∆(x − y) = θ x 0 − y 0 〈0 ¯φ (x) φ y ¯ 0〉 + θ y 0 − x 0 〈0|φ y φ (x) |0〉

¡ ¢ ¯ ¡ ¢¯ ¡ ¢ ¡ ¢

= 〈0|T φ (x) φ y |0〉

¡ ¢
(4.65)

que es el llamado propagador de Feynman. En la última línea se ha introducido el llamado

símbolo de ordenación temporal T que ordena los operadores de forma que los últimos
temporalmente aparecen a la izquierda. El propgador de Feynman es también una función
de Green del operador de Klein-Gordon y en espacio de momentos se escribe

i
∆(p) = (4.66)
p 2 − m2 + i ϵ

y como veremos juega un papel fundamental en la construcción de las llamadas reglas de

Feynman.

4.6 Cuantización del campo escalar complejo: antipartícu-

las

4.6.1 El campos escalar complejo clásico

Un campo escalar complejo puede escribirse como combinación de dos campos escalares
reales
1
φ = p (φ1 + i φ2 ). (4.67)
2

Antonio L. Maroto, Ignazio Scimemi 50

Campos Cuánticos 2022-23 UCM

La acción puede escribirse como

µ ¶
1 µν 1 2 2 1 µν 1 2 2
Z
4
S[φ] = d x η ∂µ φ1 ∂ν φ1 − m φ1 + η ∂µ φ2 ∂ν φ2 − m φ2
2 2 2 2
Z
= d 4 x(ηµν ∂µ φ∗ ∂ν φ − m 2 φ∗ φ) (4.68)

de forma que las ecuaciones de movimiento serán las mismas que en el caso real

(□ + m 2 )φ = 0, (□ + m 2 )φ∗ = 0 (4.69)

o en términos de los campos reales

(□ + m 2 )φ1,2 = 0 (4.70)

Podemos por tanto describir el campo escalar complejo en términos de dos campos inde-
pendientes φ1 y φ2 o bien φ y φ∗
Las soluciones de la ecuación de Klein-Gordon serán
d 3p
Z
φ(x) = (a p e −i px + b p∗ e i px )|p 0 =E p (4.71)
3
p
(2π) 2E p
A diferencia del caso real, ahora a p , b p son independientes.
Observando la acción S, vemos que ésta tiene una simetría adicional con respecto al
campo real dada por

φ′ (x ′ ) = e i θ φ(x)
φ′ (x ′ ) = e −i θ φ∗ (x)
∗
(4.72)

donde θ es un parámetro real. Se trata por tanto de una simetría interna global del tipo
U (1) con X µ = 0, Yφ = i φ, Yφ∗ = −i φ∗ y por tanto la corriente Noether será
∂L
jµ =
X
Yi , (4.73)
i =φ,φ∗ ∂(∂µ φi )

de donde

j µ = i (φ∂µ φ∗ − φ∗ ∂µ φ) (4.74)

y por tanto la carga conservada será

↔
Z Z
Q = d 3 x j 0 = i d 3 x φ∗ ∂0 φ. (4.75)

que como veremos puede identificarse con la carga eléctrica cuando el campo escalar se
acopla al campo electromagnético.

Antonio L. Maroto, Ignazio Scimemi 51

Campos Cuánticos 2022-23 UCM

4.6.2 Cuantización del campo escalar complejo

Podemos por tanto promover los coeficientes a p y b p a operadores y desarrollar por tanto
el campo y su hermítico como

d 3p
Z
φ= (a p e −i px + b p† e i px )|p 0 =E p (4.76)
(2π)3 2E p
p

d 3p
Z
†
φ = (b p e −i px + a p† e i px )|p 0 =E p (4.77)
(2π)3 2E p
p

donde los operadores a p , b p son independientes.

Los momentos canónicos conjugados serán
∂L
πφ = = ∂0 φ† (4.78)
∂(∂0 φ)
∂L
π†φ = = ∂0 φ (4.79)
∂(∂0 φ† )
Por otra parte la densidad Hamiltoniana será

H = π†φ πφ + ⃗
∇φ†⃗
∇φ + m 2 φ† φ (4.80)

De esta forma la relaciones de conmutación a tiempos iguales serán

h ¡ ¢i
φ (t ,⃗
x ) , πφ t ,⃗ = φ† (t ,⃗x ) , π†φ t ,⃗
y = i δ(3) ⃗
£ ¡ ¢¤ ¡ ¢
y x −⃗
y (4.81)
h ¡ ¢i
φ (t ,⃗
x ) , φ t ,⃗ = φ† (t ,⃗x ) , φ† t ,⃗
£ ¡ ¢¤
y y =0 (4.82)
h ¡ ¢i
πφ (t ,⃗x ) , πφ t ,⃗ = π†φ (t ,⃗
x ) , π†φ t ,⃗
£ ¡ ¢¤
y y =0 (4.83)

y el resto de conmutadores nulos. Las relaciones de conmutación para los operadores a p y

b p serán
h i h i
a p , a q† = b p , b q† = (2π)3 δ(3) p⃗ −⃗
¡ ¢
q (4.84)

y el resto de conmutadores nulos.

El espacio de Fock se construirá a partir del estado de vacío que se define como

a p |0〉 = b p |0〉 = 0, ∀⃗
p (4.85)

Vemos que el campo φ crea partículas de tipo b y destruye partículas de tipo a, mientras
que φ† hace lo contrario.

Antonio L. Maroto, Ignazio Scimemi 52

Campos Cuánticos 2022-23 UCM

El Hamiltoniano del sistema después de la ordenación normal, será

d 3p
Z h i
† †
H= E p a a
p p + b b
p p (4.86)
(2π)3
y el operador momento

d 3p i †
Z Z Z
i 3 0i 3 † †
P = d x θ = − d x (πφ ∂i φ + πφ ∂i φ ) = p (a p a p + b p† b p ) (4.87)
(2π)3
¡ ¢1/2 † ¡ ¢1/2 †
Los estados de una partícula tienen degeneración doble, puesto que 2E p a p |0〉 y 2E p b p |0〉
tienen la misma energía y el mismo trimomento. El Lagrangiano como vimos tiene una
simetría U (1) cuya carga asociada viene dada por
↔
Z Z
Q = d 3 x j 0 = i d 3 x φ† ∂0 φ

d 3p d 3q
Z ³³ ´ ³ ´
= i d 3x a p† e i px + b p e −i px ∂0 a q e −i q x + b q† e i q x
(2π)3 2E p (2π)3 2E q
p p
³ ´³ ´´
− ∂0 a p† e i px + b p e −i px a q e −i q x + b q† e i q x
d 3p d 3q
Z ³³ ´ ³ ´
= d 3x a p† e i px + b p e −i px E q a q e −i q x − b q† e i q x
(2π)3 2E p (2π)3 2E q
p p
³ ´³ ´´
+ E p a p† e i px − b p e −i px a q e −i q x + b q† e i q x (4.88)

donde en todo momento suponemos p 0 = E p y q 0 = E q . La integración en d 3 x produce en

p −⃗
los términos que involucran un creador y un aniquilador un término (2π)3 δ(3) (⃗ q ), mien-
tras que en los términos que involucran dos creadores o dos aniquiladores (2π)3 δ(3) (⃗ p +⃗
q ).
En ambos casos, imponen |⃗ p | = |⃗
q | y por tanto E p = E q . Consideremos por ejemplo los
† †
términos que involucran a p b q

d 3p d 3q
Z
a p† b q† (E p − E q )e i (E p +E q )t (2π)3 δ(3) (⃗
p +⃗
q) = 0 (4.89)
(2π)3 2E p (2π)3 2E q
p p

y lo mismo ocurre con los términos que involucran a p b q , por lo que finalemente podemos
escribir
d 3p d 3q
Z ³ ´
Q= a p† a q (E p + E q )e i (E p −E q )t − b p b q† (E p + E q )e −i (E p −E q )t
(2π)3 2E p (2π)3 2E q
p p

d 3p ³ †
Z ´
3 (3) †
× (2π) δ (⃗ p −⃗
q) = a a p − b b
p p (4.90)
(2π)3 p

Antonio L. Maroto, Ignazio Scimemi 53

Campos Cuánticos 2022-23 UCM

Utilizando ahora la ordenación normal4 de los operadores tenemos,

d 3p ³ †
Z ´
†
: Q := a a p − b p p = N a − Nb
b (4.91)
(2π)3 p

donde N a,b es el número total de partículas de tipo a y b respectivamente. Vemos por tanto
que a p† crea partículas con masa m, momento p ⃗, energía E p y carga +1, mientras que b p†
crea particulas con la misma masa, momento y energía pero con carga −1. Podemos por
tanto denominar a la partícula de tipo b como la anipartícula de la de tipo a. En el caso del
campo escalar real, vimos que a p = b p y por tanto la partícula es su propia antipartícula.
Nótese que [Q, H ] = 0 y por tanto de acuerdo con (4.21), la carga total N a − Nb se con-
serva, aunque los números de partículas N a y Nb podrían no conservarse separadamente.
Por otra parte, como hemos comentado antes

Q|0〉 = 0 (4.92)

el vacío tiene carga total nula.

4.7 Propagador
Siguiendo los mismos pasos que en caso del campo escalar real, puede verse que la canti-
dad

∆(x − y) = 〈0|T φ (x) φ† y |0〉

¡ ¢
(4.93)

es una función de Green del operador de Klein-Gordon donde de nuevo

d 4p i
Z
∆(x − y) = 4 p 2 − m2 + i ϵ
e −i p (x−y ) (4.94)
(2π)

Nótese que el operador φ† crea una partícula de carga +1 o bien destruye una partícula
de carga -1, mientras que φ al contrario, crea una partículas de carga -1 o destruye una
partícula de carga +1, por lo que el producto φ(x)φ† (y) no cambia la carga total.

4
Nótese que en el caso de la carga, la ordenación normal a priori no está justificada puesto que a diferencia
de la energía la carga total sí es observable. En este caso la ordenación normal se justifica porque en la
definición de la carga hay una ambigüedad en el orden de los operadores φ† y ∂0 φ. Podemos elegir φ† ∂0 φ o
∂0 φφ† . En este caso la ordenación normal garantiza que la carga del vacío es cero

Antonio L. Maroto, Ignazio Scimemi 54

5. Cuantización de campos espinoriales li-
bres

En este tema construimos la acción clásica invariante Lorentz para campo espinoriales
tanto de Weyl como de Dirac y aplicamos el método de cuantización canónica para campos
fermiónicos.

5.1 Campos espinoriales clásicos

Para construir una teoría relativista para campos espinoriales debemos partir de una ac-
ción invariante Lorentz que, como en el caso escalar, involucre derivadas de los campos.
Estas derivadas deben contraerse con objetos vectoriales construidos a partir de los pro-
pios campos espinoriales. Veamos cómo hacerlo a partir de los espinores de Weyl levógiros
y dextrógiros introducidos anterioremente.
Una propiedad importante de las matrices de Pauli es

σ2 σi σ2 = −σi ∗ i = 1, 2, 3

de forma que las matrices de transformación de los espinores de Weyl levógiros y dextró-
giros ΛL y ΛR (2.45) se pueden relacionar por

σ2 Λ∗L σ2 = ΛR . (5.1)

Esto signfica que que el espinor (σ2 ψ∗L ) se comporta como un espinor dextrógiro, es decir,
bajo transformaciones de Lorentz

σ2 ψ∗L → σ2 (ΛL ψL )∗ = σ2 Λ∗L σ2 σ2 ψ∗L = ΛR (σ2 ψ∗L ). (5.2)

Definamos ahora el espinor conjugado de ψL como

ψcL = i σ2 ψ∗L (5.3)

55
Campos Cuánticos 2022-23 UCM

Tomando ahora el complejo conjugado de la ecuación anterior y denotando por ψR al es-

pinor dextrógiro ψcL tenemos que ψL = −i σ2 ψ∗R que podemos identificar con el espinor
conjugado de ψR , es decir

ψcR = −i σ2 ψ∗R (5.4)

Por tanto, la operación de conjugación satisface

(ψcL )c = (i σ2 ψ∗L )c = −i σ2 (i σ2 ψ∗L )∗ = ψL ,

Es decir, si la aplicamos dos veces recuperamos el espinor inicial.

Consideremos ahora dos espinores de Weyl independientes ψL y ξR y definamos sus
conjugados como
ξL ≡ −i σ2 ξ∗R
ψR ≡ i σ2 ψ∗L
Denotando ahora

σµ = (1, σi ); σ̄µ = (1, −σi ) (5.5)

Se pueden construir los siguientes vectores contravariantes complejos

ξ†R σµ ψR , (5.6)
ξ†L σ̄µ ψL . (5.7)

5.1.1 El Lagrangiano de Weyl

Consideremos un único campo espinorial levógiro ψL , podemos por tanto construir un
Lagrangiano invariante Lorentz de la siguiente forma

L L = i ψ†L σ̄µ ∂µ ψL (5.8)

que es el Lagrangiano de Weyl. El factor i garantiza que la acción es real usando que las
matrices de Pauli son hermíticas.
Puesto que el campo de Weyl es complejo, las ecuaciones de movimiento se obtienen
variando tanto con respecto a ψL como a su complejo conjugado ψ∗L . Puesto que ∂µ ψ∗ no
aparece en el Lagrangiano, las ecuaciones de movimiento son simplemente

∂L L
=0 (5.9)
∂ψ∗L

Antonio L. Maroto, Ignazio Scimemi 56

Campos Cuánticos 2022-23 UCM

de donde

i σ̄µ ∂µ ψL = 0 (5.10)

que es la ecuación de Weyl y que puede escribirse como

(∂0 − σi ∂i )ψL = 0 (5.11)

Se trata por tanto de una ecuación de primer orden. Derivando respecto al tiempo y us-
ando de nuevo la ecuación llegamos a

∂20 ψL = σ j σi ∂ j ∂i ψL (5.12)

pero usando σi σ j = δi j + i ϵi j k σk , llegamos finalmente a

(∂20 − ∂2i )ψL = 0 (5.13)

es decir, la ecuación de Klein-Gordon para un campo sin masa

□ψ L = 0 (5.14)

Vemos que la ecuación de Weyl describe espinores sin masa. Por tanto las soluciones con
frecuencia positiva podrán escribirse como

ψL (x) = u L e −i px con u L = biespinor constante (5.15)

donde p 0 = E p = |⃗
p |.
Sustituyendo en la ecuación de Weyl tenemos

(∂0 − σi ∂i )u L e −i px = −i (p 0 + p
⃗ ·⃗
σ)u L = 0 (5.16)

y por tanto

⃗ ·⃗
p σ
u L = −u L (5.17)
p0

Por tanto, teniendo en cuenta que para espinores de Weyl, ⃗

J =⃗ ⃗/p 0 = p̂
σ/2 y usando que p
es un vector unitario podemos escribir

1
(p̂ · ⃗
J )u L = − u L (5.18)
2

Antonio L. Maroto, Ignazio Scimemi 57

Campos Cuánticos 2022-23 UCM

es decir, un espinor de Weyl levógiro tiene helicidad h = −1/2. El espinor dextrógiro tendrá
helicidad h = 1/2. .
El tensor energía-momento puede verse que tiene la siguiente forma

θ µν = i ψ†L σ̄µ ∂ν ψL . (5.19)

donde se ha utilizado que sobre soluciones de las ecuaciones de movimiento el Lagrangiano

se anula.
La acción de Weyl es invariante bajo transformaciones de fase globales ψL → e i θ ψL .
Existe por tanto una simetría U (1) con su corriente y correspondiente carga conservada
µ
J U(1) = ψ†L σ̄µ ψL (5.20)
Z
QU (1) = d 3 x ψ†L ψL (5.21)

De forma totalmente análoga puede construirse un Lagrangiano para fermiones de Weyl

dextrógiros dado por

L R = i ψ†R σµ ∂µ ψR (5.22)

Las ecuaciones de movimiento son

i σµ ∂µ ψR = i (∂0 + σi ∂i )ψR = 0, (5.23)

Un ejemplo de fermiones de Weyl serían los neutrinos si su masa fuese nula.

5.1.2 El Lagrangiano de Dirac

Como vimos anteriormente, los espinores de Dirac se construyen usando espinores de
Weyl levógiros y dextrógiros. Usando ambos tipos de espinores es posible construir los
siguientes escalares Lorentz
ψ†L ψR y ψ†R ψL
Bajo una transformación de Lorentz tenemos

ψ†L ψR → ψ†L Λ†L ΛR ψR = ψ†L ψR

puesto que Λ†L ΛR =1. Sin embargo bajo una transformación de paridad un espinor levógiro
se transforma en un dextrógiro y viceversa de forma que

ψ†L ψR + ψ†R ψL par bajo paridad, combinación escalar

ψ†L ψR − ψ†R ψL impar bajo paridad, combinación pseudo-escalar. (5.24)

Antonio L. Maroto, Ignazio Scimemi 58

Campos Cuánticos 2022-23 UCM

El Lagrangiano de Dirac se construye para preservar la paridad, de forma que

L D = i ψ†L σ̄µ ∂µ ψL + i ψ†R σµ ∂µ ψR − m(ψ†L ψR + ψ†R ψL ) (5.25)

donde m es una constante escalar que como veremos se identifica con la masa del campo.
Las ecuaciones de movimiento se obtienen como en el caso de Weyl considerando ψL,R
y ψ†L,R como independientes de forma que

i σ̄µ ∂µ ψL = mψR
i σµ ∂µ ψR = mψL . (5.26)

Si multiplicamos la primera ecuación por i σν ∂ν y usando la segunda obtenemos

−σµ σ̄ν ∂µ ∂ν ψL = m 2 ψL

y usando el hecho de que el producto de derivadas es simétrico y que σµ σ̄ν + σν σ̄µ = 2ηµν
tenemos

(□ + m 2 )ψL = 0. (5.27)

y de forma similar

(□ + m 2 )ψR = 0. (5.28)

por tanto ambos campos satisfacen la ecuación de Klein-Gordon con masa.

Usando la notación de cuadriespinores de Dirac

ψL
µ ¶
ψ= (5.29)
ψR

Esta es la llamada representación quiral de los espinores de Dirac.

Podemos ahora definir las matrices gamma de Dirac

σi
µ ¶ µ ¶
0 0 1 i 0
γ = ; γ = ; (5.30)
1 0 −σi 0

de forma que

0 σµ
µ ¶
µ
γ = ; (5.31)
σ̄µ 0

Antonio L. Maroto, Ignazio Scimemi 59

Campos Cuánticos 2022-23 UCM

de forma que el Lagrangiano de Dirac se escribe1

L D = ψ† γ0 (i γµ ∂µ − m)ψ ≡ ψ̄(i ∂/ − m)ψ (5.32)

donde
ψ̄ ≡ ψ† γ0
y hemos definido γµ ∂µ ≡ ∂.
/ La ecuación de movimiento correspondiente es la ecuación de
Dirac

(i ∂/ − m)ψ = 0. (5.33)

que para el espinor ψ̄ implica

−i ∂µ ψ̄γµ − m ψ̄ = 0 (5.34)

Podemos definir también la matriz γ5 como

µ ¶
5 0 1 2 3 −1 0
γ = iγ γ γ γ = . (5.35)
0 1

que permite definir los proyectores left y right como P L = 12 (1 − γ5 ) y P R = 12 (1 + γ5 ) que

seleccionan la parte levógira o dextrógira del espinor. Una propiedad fundamental de la
matriz γ5 es que anticonmuta con el resto de matrices gamma

{γµ , γ5 } = 0. (5.36)

La forma en la que hemos escrito los campos y las matrices de Dirac no es única 2 . Sin
embargo, en cualquier representación se cumple

γµ γν + γν γµ = {γµ , γν } = 2ηµν . (5.39)

1
Puesto que en unidades naturales la densidad Lagrangiana tiene dimensiones [L ] = E 4 , la dimensión del
campo de Dirac debe ser [ψ] = E 3/2
2
Otra representación de los campos y las matrices es la llamada representación estándar o de Dirac
à 1
p (ψL + ψR )
!
ψ= 2 , (5.37)
p1 (ψL − ψR )
2

y las matrices gamma se escriben como

σi
µ ¶ µ ¶ µ ¶
1 0 0 0 1
γ0 = ; γi = ; γ5 = i γ0 γ1 γ2 γ3 = . (5.38)
0 −1 −σi 0 1 0

Antonio L. Maroto, Ignazio Scimemi 60

Campos Cuánticos 2022-23 UCM

Propiedades de las matrices gamma

γλ γλ = 4
γλ γµ γλ = −2γµ
γλ γµ γν γλ = 4ηµν
tr γµ = tr γ5 = tr (γ5 γµ ) = 0
tr(número impar de γµ ) = 0
tr (γµ γν ) = 4ηµν
tr (γµ γν γρ γσ ) = 4(ηµν ηρσ − ηµρ ηνσ + ηµσ ηνρ )
tr (γ5 γµ γν γρ γσ ) = 4i ϵµνρσ con ϵ0123 = +1

5.1.3 Soluciones de la ecuación de Dirac

La solución general de la ecuación de Dirac puede escribirse como combinación de modos
de frecuencia positiva

ψ(x) = u(p)e −i px (5.40)

y modos de frecuencia negativa

ψ(x) = v(p)e i px (5.41)

con u(p) y v(p) cuadri-espinores constantes y donde las ecuaciones de Klein-Gordon (5.27)
y (5.28) imponen p 2 = m 2 con p 0 = E p . Sustituyendo en la ecuación de Dirac obtenemos

(p/ − m)u(p) = 0 (5.42)

(p/ + m)v(p) = 0 (5.43)

En la representación quiral podemos escribir

µ ¶
u L (p)
u(p) = . (5.44)
u R (p)

En el caso m ̸= 0, si nos colocamos en el sistema en reposo en el que p µ = (m,⃗0) la ecuación

se simplifica

(p/ − m)u(p) = m(γ0 − 1)u(p) = 0 ↔ u L (p) = u R (p), (5.45)

Antonio L. Maroto, Ignazio Scimemi 61

Campos Cuánticos 2022-23 UCM

Una normalización conveniente para los bi-espinores u L y u R es

p
uL = uR = mξ con ξ† ξ = 1. (5.46)

Una base ortonormal para los bi-espinores ξ será por tanto

µ ¶ µ ¶
1 1 2 0
ξ = , ξ = . (5.47)
0 1

Para obtener la solución con un momento genérico, simplemente hacemos un boost so-
bre la solución anterior usando (2.46). Suponiendo que p ⃗ está en la dirección z, obten-
emos que la solución general de la ecuación (5.42) será una combinación lineal de los
cuadriespinores:
 hp 3 3
i 
3 1−σ + 3 1+σ ξs
p
E + p 2 E − p 2 i
u s (p) =  hp 1+σ 3
1−σ3
 , s = 1, 2 (5.48)
E + p 3 2 + E − p 3 2 ξs
p

En el límite ultrarrelativista p µ → (E , 0, 0, E ) obtenemos

p p ξ2
µ ¶ µ ¶
1 0 2
u (p) → 2E , u (p) → 2E , (5.49)
ξ1 0

es decir, recuperamos los espinores de Weyl.

Para los cuadri-espinores v(p) procediendo de forma similar llegamos a
 hp i 
1−σ3 1+σ3 s
η
p
E +p 3 + E −p 3
2 2
v s (p) =  hp 1+σ 3
1−σ 3
i  (5.50)
s
η
p
− E +p 3 + E −p 3
2 2

donde al igual que antes, los bi-espinores son

µ ¶ µ ¶
1 1 2 0
η = , η = . (5.51)
0 1

Resumimos a continuación algunas de las propiedades de los espinores. Usamos la no-

tación ū = u † γ0 , v̄ = v † γ0 y la normalización

ξ†r ξs = δr s
η†r ηs = δr s (5.52)

Antonio L. Maroto, Ignazio Scimemi 62

Campos Cuánticos 2022-23 UCM

que implica

ū r (p)u s (p) = 2mδr s (5.53)

r s rs
v̄ (p)v (p) = −2mδ (5.54)
†r s rs
u (p)u (p) = 2E p δ (5.55)
†r s rs
v (p)v (p) = 2E p δ (5.56)
r s r s
ū (p)v (p) = v̄ (p)u (p) = 0. (5.57)

De aquí se sigue

u s (p)ū s (p) = p/ + m
X
(5.58)
s=1,2
v s (p)v̄ s (p) = p/ − m
X
(5.59)
s=1,2

De la ecuación de Dirac tenemos

ū(p/ − m) = 0
v̄(p/ + m) = 0 (5.60)

p µ ū r (p)γµ u s (p) = m ū r (p)u s (p) = 2m 2 δr s → ū r (p)γµ u s (p) = 2p µ δr s (5.61)

De forma similar

v̄ r (p)γµ v s (p) = 2p µ δr s . (5.62)

5.1.4 Corrientes, simetrías y transformaciones de Lorentz

Una base completa de matrices en el espacio de cuadri-espinores es

i
1, γ5 , γµ , γµ γ5 , σµν = [γµ , γν ],
2
y por tanto, podemos construir las siguientes cantidades (escalar, pseudo-escalar, corri-
ente vectorial, corriente axial, tensor) bilineales

ψ̄ψ, ψ̄γ5 ψ, ψ̄γµ ψ, ψ̄γµ γ5 ψ, ψ̄σµν ψ.

Antonio L. Maroto, Ignazio Scimemi 63

Campos Cuánticos 2022-23 UCM

En términos de las matrices γ se puede ver que la transformación de Lorentz de los

espinores de Dirac puede escribirse como
i µν
ΛD = e − 4 ωµν σ (5.63)

y por tanto,
1
J µν = σµν
2
es el generador de las transformaciones de Lorentz.
La acción de Dirac es invariante bajo transformaciones de fase globales ψ → e i θ ψ. Ex-
iste por tanto una simetría U (1) cuya corriente vectorial conservada es
µ
J V = ψ̄γµ ψ (5.64)

En el caso m = 0, se puede ver que la acción de Dirac también es invariante bajo transfor-
5
maciones axiales ψ → e i θγ ψ cuya corriente axial conservada es
µ
J A = ψ̄γµ γ5 ψ (5.65)

La acción de Dirac es además invariante bajo una serie de tranformaciones discretas.

En concreto:

Conjugación de carga (C)

x ) → ψc (t ,⃗
ψ(t ,⃗ x ) = −i γ2 ψ∗ (t ,⃗
x)
Como vemos transforma los espinores L y R en sus correspondientes conjugados de carga.

Paridad (P)

x ) → (t , −⃗
(t ,⃗ x)
ψ(t ,⃗ x ′ ) = γ0 ψ(t , −⃗
x ) → ψ′ (t ′ ,⃗ x)
Intercambia los espinores L y R.

Inversión temporal (T)

x ) → (−t ,⃗
(t ,⃗ x)
ψ(t ,⃗ x ′ ) = −γ1 γ3 ψ(−t ,⃗
x ) → ψ′ (t ′ ,⃗ x)

Antonio L. Maroto, Ignazio Scimemi 64

Campos Cuánticos 2022-23 UCM

5.2 Cuantización del campo de Dirac

Consideremos el Lagrangiano para un espinor de Dirac con masa m

L = ψ̄ i γµ ∂µ − m ψ
¡ ¢
(5.66)

donde ψ̄ = ψ† γ0 y la correspondiente ecuación de movimiento es

(i ∂/ − m)ψ = 0. (5.67)

con ∂/ = γµ ∂µ . La ecuación para el espinor ψ̄ se puede escribir como

−i ∂µ ψ̄γµ − m ψ̄ = 0 (5.68)

que se expresa habitualmente como ψ̄(i ∂/ + m) = 0. La ecuación de Dirac puede por tanto
entenderse como la ecuación para un campo complejo ψ. El correspondiente momento
canónico conjugado será3
∂L
πψ = = i ψ̄γ0 = i ψ† (5.69)
∂(∂0 ψ)
y la densidad Hamiltoniana se puede escribir como

γ ·⃗
H = πψ ∂0 ψ − L = i ψ† ∂0 ψ − ψ̄ i γ0 ∂0 + i⃗ γ ·⃗
∇ − m ψ = ψ̄ −i⃗ ∇+m ψ
¡ ¢ ¡ ¢
(5.70)

A diferencia de los campos escalares anteriores que cumplen la estadística bosónica,

en el caso de fermiones las reglas de conmutación se establecen con anticonmutadores en
vez de conmutadores. En concreto, si usáramos conmutadores puede verse que el Hamil-
toniano no estaría acotado inferiormente (ver Peskin para una discusión detallada). Por
tanto las reglas de (anti)conmutación a tiempos iguales serían:
n ¡ ¢o
x ) , ψ†b t ,⃗
ψa (t ,⃗ = δ(3) ⃗ y δab
¡ ¢
y x −⃗
n ¡ ¢o
ψa (t ,⃗
x ) , ψb t ,⃗ x ) , ψ†b t ,⃗
= ψ†a (t ,⃗
© ¡ ¢ª
y y =0 (5.71)

donde los índices a, b = 1 . . . 4 corresponden a las 4 componentes del espinor de Dirac.

Como vimos anteriormente, el desarollo de una solución general de la ecuación de
Dirac se escribe
d 3p
Z X ³ s
¡ ¢ −i px † s
¡ ¢ i px ´
ψ (x) = a p,s u p e + b p,s v p e (5.72)
(2π)3 2E p s=1,2
p

3
El momento canónico conjugado de ψ† es nulo.

Antonio L. Maroto, Ignazio Scimemi 65

Campos Cuánticos 2022-23 UCM

d 3p
Z X ³ s
¡ ¢ −i px † s
¡ ¢ i px ´
ψ̄ (x) = b p,s v̄ p e + a p,s ū p e (5.73)
(2π)3 2E p s=1,2
p

donde estamos suponiendo p 0 = E p y los operadores de creación y destrucción tienen

tanto índices de momento como de spin. De nuevo, podemos ver que las siguientes rela-
ciones de anticonmutación satisfacen las relaciones (5.71)
n o n o
† †
= (2π)3 δ(3) p ⃗ −⃗
q δr s
¡ ¢
a p,r , a q,s = b p,r , b q,s
© ª © ª n † †
o n
† †
o
a p,r , a q,s = b p,r , b q,s = a p,r , a q,s = b p,r , b q,s =0 (5.74)

y el resto de anticonmutadores cero.

Definimos de nuevo el vacío como

a p,s |0〉 = b p,s |0〉 = 0, ∀⃗

p, s (5.75)

Puesto que los operdores de creación anticonmutan, los estados correspondientes a partícu-
las múltiples son antisimétricos en el intercambio de dos partículas como corresponde a
la estadística de Fermi-Dirac. Los estados de una partícula serán por tanto
q
†
p , s〉 =
|⃗ 2E p a p,s |0〉 (5.76)

que cumplen por tanto

〈⃗ p , s〉 = 2E p (2π)3 δ(3) p⃗ −⃗
q δr s
¡ ¢
q , r |⃗ (5.77)

y de forma totalmente análoga para las antipartículas b.

Sustituyendo en la expresión del Hamiltoniano, obtenemos depués de la ordenación
normal

d 3p X
Z Z ³ ´
3
: H := d x : ψ̄ −i⃗ ⃗
γ · ∇ + m ψ := † †
¡ ¢
E p a p,s a p,s + b p,s b p,s (5.78)
(2π)3 s=1,2

y de forma totalmente análoga en la expresión para el operador del momento lineal total

d 3p X
Z Z ³ ´
⃗= d x ψ −i ⃗
3 †
∇ ψ= †
⃗ a p,s †
¢ ¡
P p a p,s + b p,s b p,s (5.79)
(2π)3 s=1,2

Antonio L. Maroto, Ignazio Scimemi 66

Campos Cuánticos 2022-23 UCM

Como vimos, el Lagrangiano de Dirac también tenía una simetría U (1) cuya carga viene
dada por

d 3p X ³
Z Z ´
3 † † †
: Q := d x : ψ ψ := a p,s a p,s − b p,s b p,s (5.80)
(2π)3 s

Nótese que en el caso fermiónico, debido a la anticonmutación de los operadores, la orde-

† †
nación normal lleva un signo menos : a p,s a q,r := −a q,r a p,s por cada intercambio de un par
de operadores.
Por tanto, el operador a † crea fermiones con carga fermiónica +1 (partículas) mientras
que el operador b † crea partículas con carga fermiónica −1 (antipartículas). Veremos que
cuando acoplemos el campo del electrón al campo electromagnético encontraremos que
esta carga carga fermiónica se relacionará con la carga eléctrica por un factor de normal-
ización.
Finalmente, podemos también calcular el operador de spin a partir de la corriente de
Noether obtenida en (3.76). El momento angular total para el campo de Dirac puede es-
cribirse como
µ ¶
1
Z
⃗J = d 3 x ψ† ⃗ ∇) + ⃗
x × (−i ⃗ Σ ψ (5.81)
2

donde

σi
µ ¶
i 0
Σ = (5.82)
0 σi

†
Usando la expresión anterior podemos calcular el spin del estado a 0,s |0〉. Por simplicidad
µ
nos limitaremos al caso de una partícula en reposo p = (m, 0) en un estado de "spin" s.
Para la componente J 3 obtenemos

† 1 †
J 3 a 0,1 |0〉 = a 0,1 |0〉
2
† 1 †
J 3 a 0,2 |0〉 = − a 0,2 |0〉 (5.83)
2
†
Es decir el estado a 0,1 |0〉 es autoestado de la tercera componente del momento angular con
†
autovalor +1/2 y a 0,2 |0〉 con autovalor −1/2, que es lo que esperaríamos para una partícula
de spin 1/2 como el electrón.

Antonio L. Maroto, Ignazio Scimemi 67

Campos Cuánticos 2022-23 UCM

Estado H ⃗
P J3 Q

†
a 0,1 |0〉 Ep ⃗
p +1/2 +1
†
a 0,2 |0〉 Ep ⃗
p -1/2 +1
†
b 0,1 |0〉 Ep ⃗
p -1/2 -1
†
b 0,2 |0〉 Ep ⃗
p +1/2 -1

Table 5.1: Estados de una partícula de un fermión de Dirac

†
Podríamos haber repetido el cálculo con los estados b 0,s |0〉, el resultado tiene el signo
cambiado, es decir

† 1 †
J 3 b 0,1 |0〉 = − b 0,1 |0〉
2
† 1 †
J 3 b 0,2 |0〉 = b 0,2 |0〉 (5.84)
2
Por tanto, la antipartícula con s = 2 es la que tiene momento angular 1/2. Vemos en la Tabla
† †
las propiedades de los estados de una partícula a 0,s |0〉 y b 0,s |0〉 que efectivamente permite
identificarlas como fermión y antifermión.
Vemos por tanto que al igual que ocurría con los escalares, podemos interpretar el es-
tado ψs (0,⃗x )|0〉 como el correspondiente a un antifermión (positrón) en la posición ⃗ x con
la polarización de spin corresondiente a s (seleccionamos una de las dos componentes)
mientras que ψ̄s (0,⃗ x )|0〉 correspondería a un fermión (electrón) en la posición ⃗ x con spin
correspondiente a s.

5.3 El propagador de Dirac

Calculemos en primer lugar el valor esperado

d 3p d 3q
Z
〈0 ¯ψa (x) ψ̄b y ¯ 0〉 = 〈0|
¯ ¡ ¢¯
(2π)3 2E p (2π)3 2E q
p p
X³ ´³ ´
a p,s u as p e −i px + b p,s
†
v as p e i px b q,r v̄ br q e −i q y + a q,r
†
ū br q e i q y |0〉
¡ ¢ ¡ ¢ ¡ ¢ ¡ ¢
(5.85)
s,r

Antonio L. Maroto, Ignazio Scimemi 68

Campos Cuánticos 2022-23 UCM

Puesto que sólo los términos de la forma aa † sobreviven y usando la regla de anticon-
mutación llegamos a
d 3p d 3q
Z X s ¡ ¢ r ¡ ¢ −i px i q y †
〈0 ¯ψa (x) ψ̄b y ¯ 0〉 =
¯ ¡ ¢¯
u a p ū b q e e 〈0|a p,s a q,r |0〉
3 3
p p
(2π) 2E p (2π) 2E q s,r
d 3p d 3q
Z X s ¡ ¢ r ¡ ¢ −i px i q y
= u a p ū b q e e (2π)3 δ(3) (⃗p −⃗
q )δr s 〈0|0〉
(2π)3 2E p (2π)3 2E q s,r
p p

d 3p d 3p
Z X s ¡ ¢ s ¡ ¢ −i p(x−y) Z
= u p ū p e = (p/ + m)e −i p(x−y)
(2π)3 (2E p ) s a b
(2π)3 (2E p ) ab
d 3p
Z
= (i ∂/x + m)ab 3
e −i p(x−y) (5.86)
(2π) (2E p )
De forma totalmente análoga obtenemos
d 3p
Z X s ¡ ¢ s ¡ ¢ −i p(y−x)
〈0 ψ̄b y ψa (x) 0〉 =
¯ ¡ ¢ ¯
¯ ¯ v p v̄ b p e
(2π) (2E p ) s a
3

d 3p d 3p
Z Z
−i p(y−x)
= 3
(p/ − m)ab e = −(i ∂/x + m)ab 3
e −i p(y−x)
(2π) (2E p ) (2π) (2E p )
(5.87)
Como hicimos en el caso de Klein-Gordon podemos definir la función de Green retardada4
S Rab x − y = θ x 0 − y 0 〈0 ¯ ψa (x) , ψ̄b (y) ¯ 0〉
¡ ¢ ¡ ¢ ¯© ª¯
(5.88)
Podemos comprobar inmediatamente a partir de (4.55) que
S R x − y = (i ∂/x + m)D R (x − y)
¡ ¢
(5.89)
con D R (x − y) la función de Green retardada de Klein-Gordon (4.60). Podemos entonces
escribir la función de Green en espacio de Fourier como
d 4p i
Z
S R x − y = (i ∂/x + m)D R (x − y) = (i ∂/x + m) e −i p (x−y )
¡ ¢
4 p 2 − m2
(2π)
4
d p i (p/ + m) −i p (x−y )
Z
= e (5.90)
(2π)4 p 2 − m 2
4
Puede verse que si en vez de anticonmutadores hubiéramos usado conmutadores el propagador retar-
dado no se cancelaría para sucesos con separación espacial, violando por tanto causalidad. Estos resultados
se pueden resumir en el llamado teorema de conexión spin-estadística: campos con spin semientero deben
cuantizarse como fermiones (es decir, con anti-conmutadores), mientras que los campos con spin entero
deben cuantizarse como bosones (es decir, con conumtadores). Si una teoría de campos se cuantiza con
la conexión incorrecta, la teoría se vuelve no local, con violaciones de causalidad y/o no tiene un estado
fundamental, o contiene estados en su espectro con norma negativa.

Antonio L. Maroto, Ignazio Scimemi 69

Campos Cuánticos 2022-23 UCM

donde la integral en p 0 se ha realizado de acuerdo con el contorno de la Fig. 4.1. Nótese

que S R es una matriz 4 × 4 en los índices de Dirac. Podemos demostrar fácilemente que S R
es función de Green del operador de Dirac.

(i ∂/x − m)S R (x − y) = i δ(4) (x − y) (5.91)

De forma totalmente análoga podríamos haber hecho la integral usando el contorno

de la Fig. 4.2, de forma que obtenemos el propagador de Feynman fermiónico.

d 4 p i (p/ + m) −i p (x−y )
Z
¡ ¢
SF x − y = e
(2π)4 p 2 − m 2 + i ϵ
( ¯
〈0 ¯ψ (x) ψ̄ y ¯ 0〉 para x 0 > y 0
¡ ¢¯
=
−〈0 ¯ψ̄ y ψ (x)¯ 0〉 para x 0 < y 0
¯ ¡ ¢ ¯

= 〈0 ¯T ψ (x) ψ̄ y ¯ 0〉
¯ ¡ ¢¯
(5.92)

Donde nótese que la ordenación temporal de operadores fermiónicos tiene un signo menos
cuando se intercambia el orden de los operadores.
Podemos escribir el propagador de Feynman para fermiones de Dirac en espacio de
Fourier como
i (p/ + m)
S F (p) = (5.93)
p 2 − m2 + i ϵ

5.4 Fermiones de Weyl sin masa

Consideraremos ahora la cuantización de fermiones de Weyl sin masa. Podemos utilizar
la notación con cuatro componentes de Dirac. Así podemos escribir la expresión para un
fermión de Weyl levórigo en la representación quiral como

ψL
µ ¶
ΨL = (5.94)
0

y usando el desarrollo en operadores de creación y destrucción

d 3p
Z X ³ s
¡ ¢ −i px † s
¡ ¢ i px ´
ΨL (x) = a p,s u L p e + b p,s L p e
v (5.95)
(2π)3 2E p s=1,2
p

donde u L y v L son espinores con cuatro componentes de la forma

µ s¶ µ s¶
s UL s V
uL = , vL = L (5.96)
0 0

Antonio L. Maroto, Ignazio Scimemi 70

Campos Cuánticos 2022-23 UCM

Estado (ΨL ) H ⃗
P J3 Q Estado (ΨR ) H ⃗
P J3 Q

† †
a p,2 |0〉 Ep ⃗
p -1/2 +1 a p,1 |0〉 Ep ⃗
p 1/2 +1
† †
b p,2 |0〉 Ep ⃗
p +1/2 -1 b p,1 |0〉 Ep ⃗
p -1/2 -1

Table 5.2: Estados de una partícula para fermiones de Weyl levógiros y dextrógiros

Puesto que en este caso m = 0 podemos comparar con el límite ultrarrelativista de los es-
pinores de Dirac, de forma que vemos que para s = 1, UL1 = VL1 = 0 y por tanto sólo con-
tribuye s = 2. Es decir,

d 3p
Z ³
2
¡ ¢ −i px † 2
¡ ¢ i px ´
ΨL (x) = a u
p,2 L p e + b v
p,2 L p e (5.97)
(2π)3 2E p
p

y por tanto

d 3p
Z ³
2
¡ ¢ −i px † 2
¡ ¢ i px ´
Ψ̄L (x) = b p,2 v̄ L p e + a p,2 L p e
ū (5.98)
(2π)3 2E p
p

† †
La Tabla 5.2 indica que b p,2 crea un antifermión con helicidad +1/2 y a p,2 un fermión de
helicidad −1/2. Es decir, un fermión de Weyl levórigo sólo describe partículas con helici-
dad −1/2 y antipartículas de helicidad 1/2.
Por su parte, para un fermión de Weyl dextrógiro,
µ ¶
0
ΨR = (5.99)
ψR

Ahora sólo el término s = 1 sobrevive y por tanto describe partículas de helicidad 1/2 y
antipartículas de helicidad −1/2. Si ignoramos la pequeña masa de los neutrinos del Mod-
elo Estándar, estos corresponderían a fermiones de Weyl levórigos. La carga (5.21) en este
caso no puede interpretarse como carga eléctrica puesto que los neutrinos no se acoplan
al campo electromagnético, sino como número léptónico.

Antonio L. Maroto, Ignazio Scimemi 71

6. Cuantización del campo electromagnético
libre

6.1 El campo electromagnético clásico

El campo electromagnético se describe a través de un campo vectorial A µ . El tensor de
⃗ = (E 1 , E 2 , E 3 ) y magnéticos B
Faraday que contiene los campos eléctricos E ⃗ = (B 1 , B 2 , B 3 )
está dado por
F µν = ∂µ A ν − ∂ν A µ (6.1)
F 0i = ∂0 A i − ∂i A 0 = ∂0 A i + ∂i A 0 = −E i (6.2)
F i j = ∂i A j − ∂ j A i = −ϵi j k B k (6.3)
El Lagrangiano invariante Lorentz se puede construir contrayendo el tensor de Faraday de
la siguiente forma 1
1 1 2
⃗ −B
L = − F µν F µν = (E ⃗2) (6.4)
4 2
Las ecuaciones de movimiento correspondientes son las ecuaciones de Maxwell libres
∂µ F µν = 0 ↔ ⃗ ⃗ = 0; ⃗
∇·E ⃗ = ∂0 E
∇×B ⃗ (6.5)
Definiendo el tensor dual del tensor de Faraday
1 ∗ µν
= ϵµνρσ F ρσ
F
2
y derivando obtenemos las llamadas identidades de Bianchi que conducen a las otras dos
ecuaciones de Maxwell
0 = ∂µ ∗ F µν = ϵµνρσ ∂µ ∂ρ A σ ↔ ⃗ ⃗ = 0; ⃗
∇·B ⃗ = −∂0 B
∇×E ⃗. (6.6)
1
Puesto que en unidades naturales la densidad Lagrangiana tiene dimensiones [L ] = E 4 , la dimensión del
campo debe ser [A µ ] = E

72
Campos Cuánticos 2022-23 UCM

6.1.1 Simetría gauge

Una simetría fundamental del Lagrangiano de Maxwell, aparte de la simetría Lorentz, es la
simetría gauge dada por

A µ (x) → A ′µ (x) = A µ (x) − ∂µ θ(x) (6.7)

donde θ(x) es una función arbitraria de x. Se trata por tanto de una simetría interna local.
Como puede verse directamente, el tensor de Faraday es invariante bajo transformaciones
gauge y por tanto el Lagrangiano también lo es.
La invariancia gauge implica que de los cuatro campos A µ algunos son redundantes.
Así por ejemplo, podemos encontrar una transformación gauge que haga A ′0 = 0. En efecto,
Z t
Aµ → A ′µ = A µ − ∂µ d t ′ A 0 (⃗
x , t ′ ).

Nótese que aún podemos hacer transformaciones gauge que mantengan la condición A ′′0 =
0, en concreto todas aquéllas en las que el parámetro θ = θ(⃗
x ) no dependa de t

A ′µ → A ′′µ = A ′µ − ∂µ θ(⃗
x ),

Escogiendo
d3y ∂ ′i
Z
θ(⃗
x) = − A (⃗
y , t ).
y | ∂y i
x −⃗
4π|⃗
x ) definido de esta forma es independiente del tiempo ∂0 ∂i A ′i = ∂i (−E i +
vemos que θ(⃗
∂ A ) = −⃗
i ′0 ⃗ = 0, donde hemos usado A ′0 = 0. Usando ahora
∇·E
µ ¶
1
⃗
∇x2
= −δ3 (⃗
x −⃗
y)
x −⃗
4π|⃗ y|

obtenemos
∇· ⃗
⃗ ∇· ⃗
A ′′ = ⃗ A′ − ⃗
∇2 θ = 0.
Por tanto, usando la simetría gauge hemos encontrado que podemos fijar el llamado gauge
de radiación

∇· ⃗
A 0 = 0, ⃗ A = 0. (6.8)

El gauge de radiación implica el llamado gauge de Lorentz

∂µ A µ = 0, (6.9)

Antonio L. Maroto, Ignazio Scimemi 73

Campos Cuánticos 2022-23 UCM

Nótese sin embargo que imponer el gauge de Lorentz no garantiza estar en el gauge de
radiación.
En el gauge de radiación, las ecuaciones de movimiento se reducen a

∂µ F µν = ∂µ (∂µ A ν − ∂ν A µ ) = □ A ν = 0 (6.10)

Es decir, las componentes del campo vectorial se comportan como cuatro campos es-
calares sin masa.
Para los campo electromagnéticos, la simetría gauge y el hecho de que no tienen masa
implica que sólo dos de las cuatro componentes son independientes (los fotones sólo tienen
dos grados de libertad). Podemos por tanto escribir las soluciones reales de las ecuaciones
(6.10) como

A µ (x) = εµ (k)e −i kx + ε∗µ (k)e i kx (6.11)

con k 2 = 0 y εµ k µ = 0. En el gauge de radiación ε0 = 0 y ⃗ε · ⃗

k = 0 que corresponde a dos
grados de libertad.
Nótese que si añadimos un término de masa al Lagrangiano

1 M2
L = − F µν F µν + Aµ Aµ (6.12)
4 2
la teoría no respeta la invariancia gauge. Este es el Lagrangiano de Proca que describe
campos vectoriales con masa como pueden ser los bosones vectoriales W, Z . En este caso
el campo tiene tres grados de libertad.

6.1.2 Tensor energía-momento

Podemos calcular el tensor energía-momento a partir del Lagrangiano de Maxwell
1
θ µν = −F µρ ∂ν A ρ + ηµν (F ρσ F ρσ ) (6.13)
4
Vemos que esta expresión no es explícitamente invariante gauge puesto que depende de
los campos A ρ . Sin embargo, bajo transformaciones gauge y usando las ecuaciones del
movimiento ∂µ F µν = 0 tenemos

θ µν → θ µν + F µρ ∂ν ∂ρ θ = θ µν + ∂ρ (F µρ ∂ν θ) (6.14)

y las cargas conservadas cambian como

Z Z
P → P + d x ∂ρ (F ∂ θ) = P + d 3 x ∂i (F 0i ∂ν θ)
ν ν 3 0ρ ν ν
(6.15)

Antonio L. Maroto, Ignazio Scimemi 74

Campos Cuánticos 2022-23 UCM

y por tanto el último término se integra a cero. Es decir, la energía y el momento son in-
variantes gauge.
Podemos obtener una expresión explícitamente invariante gauge añadiendo un tér-
mino ∂ρ (F µρ A ν ) a θ µν que como vimos conduce a las mismas cargas conservadas. De esta
forma se obtiene
1
T µν = F µρ F ρν + ηµν F 2 (6.16)
4

con F 2 = F αβ F αβ . Esta expresión conduce directamente a la energía

1 2
Z Z
0
P = 3
d xT 00
= d 3x ⃗ +B
(E ⃗2) (6.17)
2

y al vector de Poynting para el trimomento del campo

Z Z
⃗
P = d x T = d 3x E
3 0i ⃗ ×B
⃗
¡ ¢
(6.18)

6.1.3 Interacciones electromagnéticas

En electrodinámica clásica podemos incluir una corriente externa a las ecuaciones de Maxwell
de forma que
∂µ F µν = j ν
y puesto que ∂ν ∂µ F µν = 0, se debe cumplir que ∂ν j ν = 0, es decir, la corriente se conserva.
La acción correspondiente es
µ ¶
1 µν
Z
4 µ
S = d x − F F µν − j A µ
4

La conservación de la corriente garantiza que la acción es invariante gauge. En efecto, bajo

transformaciones gauge tendremos
Z Z
S → S ′ = S + d 4 x j µ ∂µ θ = S − d 4 x θ(∂µ j µ ) = S

Podemos ahora obtener la corriente electromagnética j µ en términos de los campos

asociados a las partículas cargadas, por ejemplo, en términos de los campos espinoriales
de los electrones. Para ello, vamos a ver que la simetría gauge establece la forma en la que
deben acoplarse los campos electromagnético y los campos fermiónicos. Recordemos que
la acción de Dirac era invariante bajo transformaciones internas globales (pero no locales).

Antonio L. Maroto, Ignazio Scimemi 75

Campos Cuánticos 2022-23 UCM

Podemos intentar hacer el Lagrangiano de Dirac invariante bajo tranformaciones locales

(gauge) de la forma
ψ → e i qθ(x) ψ,
donde q es un número que podremos identificar con la carga del campo espinorial y θ(x)
es la función de la transformación gauge. Para ello tenemos que modificar el término con
derivada e introducir la llamada derivada covariante gauge

D µ = ∂µ + i q A µ (6.19)

de forma que el Lagrangiano de Dirac con acoplamiento electromagnético se escribe

L D = ψ̄(i D
/ − m)ψ (6.20)

Puede verse fácilmente que

D µ ψ → D µ′ ψ′ = e i qθ(x) D µ ψ

y por tanto el Lagrangiano es invariante gauge. Vemos que la derivada covariante gauge
introduce automáticamente un acoplo entre el campo electromagnético y los campos es-
pinoriales (acoplamiento mínimo). La corriente electromagnética corresponde de hecho
a la corriente vectorial j µ = q ψ̄γµ ψ 2 .
El mismo acoplamiento mínimo puede emplearse con los campos escalares complejos

L = D µ φ∗ D µ φ − m 2 φ∗ φ (6.21)

6.2 Cuantización del campo electromagnético libre

La cuantización del campo electromagnético, a pesar de ser la primera en llevarse a cabo
históricamente, presenta dificultades adicionales con respecto a las de partículas de spin
0 y 1/2 vistas anteriormente. Estas dificultades aparecen debido a la simetría gauge que
hacen que algunas de las componentes del campo A µ sean redundantes. En efecto, con-
sideremos el Lagrangiano del campo electromagnético libre
1
L = − F µν F µν (6.22)
4
donde F µν = ∂µ A ν − ∂ν A µ . Las correspondientes ecuaciones de movimiento son como vi-
mos

□ A µ − ∂µ (∂ν A ν ) = 0 (6.23)
2
En esta notación la carga del electrón sería q = e con e < 0

Antonio L. Maroto, Ignazio Scimemi 76

Campos Cuánticos 2022-23 UCM

El momento canónico conjugado de A µ será

∂L
πµ = = F µ0 (6.24)
∂(∂0 A µ )

Encontramos aquí el primer obstáculo para la cuantización canónica y es que el campo A 0

no tiene momento conjugado puesto que π0 = 0.
Como veremos a continuación hay fundamentalmente dos formas de evitar este prob-
lema que nos llevarán a la cuantización en el gauge de radiación y a la llamada cuanti-
zación covariante.

6.2.1 Cuantización en el gauge de radiación

Fijemos por tanto el gauge de radiación de forma que

∇· ⃗
A 0 = 0, ⃗ A=0 (6.25)

En este gauge tenemos que ∂µ A µ = 0 y por tanto las ecuaciones de movimiento se reducen
a

□ Ai = 0 (6.26)

cuya solución como vimos puede escribirse como

d 3p
Z X ³ ¡ ¢ † i px ´
⃗ ⃗
ε ⃗ λ −i px
ε ∗
⃗ λ
¢ ¡
A (x) = p , a p,λ e +⃗ p , a p,λ e (6.27)
(2π)3 2ωp λ=1,2 p 0 =ωp
p

donde usamos la notación ωp = |⃗ ∇· ⃗

p |. La condición de gauge ⃗ ⃗ ·⃗ε(⃗
A = 0 impone p p , λ) = 0.
Por tanto λ corresponde a las dos polarizaciones independientes transversas a p ⃗. Las dos
componentes físicas del campo electromagnético están dadas por los vectores de polar-
ización unitarios ⃗ε(⃗
p , 1) y ⃗ε(⃗ ⃗. Puede de-
p , 2) ortogonales entre sí y perpendiculares a p
mostrarse que

X³ i ¡ ¢´ pi p j
ε p ⃗, λ ε j ∗ p⃗, λ = δi j − 2 ≡ P i j (⃗
¢ ¡
p) (6.28)
λ ⃗
p

El momento canónico conjugado de las componentes físicas A i será

πi = F i 0 = ∂0 A i = E i (6.29)

y puesto que A 0 = 0 no tenemos ahora el problema con su momento conjugado.

Antonio L. Maroto, Ignazio Scimemi 77

Campos Cuánticos 2022-23 UCM

En principio, las reglas de conmutación deberían ser [A i (t ,⃗

x ), π j (t ,⃗
y )] = −i δi j δ(3) (⃗
x−
⃗
y ). Sin embargo esta regla de conmutación no es compatible con la condición de gauge.
Efectivamente, si tomamos la divergencia con respecto a ⃗ x , tendremos que el miembro
de la izquierda se anula mientrar que el miembro de la derecha no lo hace. La forma de
resolver este problema es sustituir el miembro de la derecha por la llamada delta de Dirac
transversa que viene dada por

d 3k ki k j
Z µ ¶
ij i⃗
k (⃗ y)
x −⃗ ij
δt r (⃗
x −⃗
y) = e δ − (6.30)
(2π)3 ⃗
k2
ij y ij
x −⃗
que efectivamente cumple ∂ix δt r (⃗ x −⃗
y ) = ∂i δt r (⃗ y ) = 0. Por tanto, la regla de con-
mutación correcta será
ij
[A i (t ,⃗
x ), π j (t ,⃗
y )] = −i δt r (⃗
x −⃗
y) (6.31)

y el resto de conmutadores nulos. Podemos por tanto, al igual que en los casos anteriores,
imponer las reglas de conmutación
† 3 (3)
[a p,λ , a q,λ′ ] = δλλ′ (2π) δ p −⃗
(⃗ q) (6.32)
† †
[a p,λ , a q,λ′ ] = [a p,λ , a q,λ′] = 0 (6.33)

Podemos ahora construir el espacio de Fock definiendo el vacío como

p , λ = 1, 2
a p,λ |0〉 = 0, ∀⃗ (6.34)

El operador Hamiltoniano se puede obtener a partir de la definición canónica

1 2
⃗ +B
H = πi ∂0 A i − L = (E ⃗2) (6.35)
2
y por tanto usando la ordenación normal tenemos

1 d 3p X
Z Z
: H := ⃗2 + B
d 3x : E ⃗ 2 := †
ωp a p,λ a p,λ (6.36)
2 (2π)3 λ=1,2

De forma análoga, el operador momento corresponde con

d 3p X
Z Z
⃗ :=
:P 3 ⃗ ×B
d x :E ⃗ := †
⃗ a p,λ
p a p,λ (6.37)
3
(2π) λ=1,2

Antonio L. Maroto, Ignazio Scimemi 78

Campos Cuánticos 2022-23 UCM

p †
Por tanto vemos que los estados de una partícula 2ωp a p,λ |0〉 corresponden a partículas
con energía ωp y momento p ⃗. Para determinar el spin de los estados, debemos de nuevo,
como hicimos en el caso fermiónico construir la corriente de Noether correspondiente al
momento angular.
En este caso las expresión para el operador momento angular total es
Z ³ ¡ ´
⃗
J = d 3 x : πi ⃗x ×⃗ π× ⃗
∇ A i +⃗
¢
A : (6.38)

Podemos de nuevo identificar el segundo término con el operador de spin de forma que
Z
⃗ π× ⃗
S = d 3x : ⃗
¡ ¢
A : (6.39)

y podemos por tanto determinar el spin de los estados de una partícula.

Consideremos el caso S 3 con p ⃗ = (0, 0, p). Por tanto los vectores de polarización pode-
mos tomarlos como ⃗ε(⃗p , 1) = (1, 0, 0), ⃗ε(⃗ ⃗·
p , 2) = (0, 1, 0). En este caso la helicidad h = p
⃗ p| = S3
S/|⃗

† †
S 3 a p,1 |0〉 = i a p,2 |0〉 (6.40)

y de forma totalmente análoga tenemos

† †
S 3 a p,2 |0〉 = −i a p,1 |0〉 (6.41)

Es decir, las polarizaciones lineales no son autoestados de la helicidad. Sin embargo las
polarizaciones circulares
† †
a p,1 ± i a p,2
†
a p,± = p (6.42)
2

sí lo son, de forma que

† †
S 3 a p,+ |0〉 = a p,+ |0〉
S 3 a p,−
† †
|0〉 = −a p,− |0〉 (6.43)

La conclusión por tanto es que los estados 2ωp a p,± describen partículas con energía ωp ,
p

momento p ⃗, spin 1 y helicidad ±1. Estas partículas son los fotones. Nótese que en este
caso, como en el de campo escalar real, los fotones son sus propias antipartículas.

Antonio L. Maroto, Ignazio Scimemi 79

Campos Cuánticos 2022-23 UCM

6.2.2 Propagador del fotón en el gauge de radiación

Como en los casos anteriores comencemos calculando
d 3p d 3q
Z
i j †
X i¡ ¢ j∗
〈0|A (x)A (y)|0〉 = 3
p 3
p ε ⃗
p , λ ε (⃗ q , λ′ )e −i px e i q y 〈0|a p,λ a q,λ′ |0〉
(2π) 2ωp (2π) 2ωq λ,λ′
d 3q d 3p
Z X i¡ Z
¢ j∗ −i p(x−y)
= ε ⃗
p , λ ε p
(⃗ , λ)e = P i j (⃗
p )e −i p(x−y) (6.44)
(2π)3 2ωp λ (2π)3 2ωp

Podemos por tanto calcular ahora el propagador retardado

d 3q
Z ³ ´
ij
D R (x − y) = θ(x 0 − y 0 )〈0|[A i (x), A j (y)]|0〉 = θ(x 0 − y 0 ) P ij
p
(⃗ ) e −i p(x−y)
− e −i p(y−x)
(2π)3 2ωq
d 4p i j i
Z
= 4
P (⃗ p ) 2 e −i p(x−y) (6.45)
(2π) p

donde el contorno en la integración en p 0 corresponde a la Fig. 4.1. De forma totalmente

análoga, la expresión para el propagador de Feynman es

d 4p i j i
Z
ij i j
D (x − y) = 〈0|T A (x), A (y)|0〉 = P (⃗
p) 2 e −i p(x−y) (6.46)
(2π) 4 p +iϵ
donde el contorno de integración es el de la Fig. 4.2. Por tanto en espacio de Fourier el
propagador de Feynman del fotón en el gauge de radiación es
i
D i j (p) = P i j (⃗
p) (6.47)
p2 + i ϵ

6.2.3 Cuantización covariante

Como hemos visto antes, el problema fundamental a la hora de cuantizar la teoría de
Maxwell es que no hay un término ∂0 A 0 en el Lagrangiano que permita obtener un mo-
mento π0 no nulo. La cuantización en el gauge de radiación resuelve el problema elim-
inando la componente A 0 como campo físico a través de la condición de gauge A 0 = 0.
Sin embargo esta condición no es invariante Lorentz con lo que se pierda la covariancia
manifiesta de la teoría.
Una forma de evitar ambos problemas es el llamado método de cuantización covari-
ante. El primer paso es la modifición del Lagrangiano de la siguiente forma
1 1
L ′ = − F µν F µν − (∂µ A µ )2 (6.48)
4 2

Antonio L. Maroto, Ignazio Scimemi 80

Campos Cuánticos 2022-23 UCM

En principio esta teoría es muy distinta de la teoría de Maxwell (no es ni siquiera invariante
gauge) sin embargo las ecuaciones de movimiento para todas las componentes son

□ Aµ = 0 (6.49)

y por tanto coinciden para los modos transversos con las del gauge de radiación.
Podemos ahora calcular los momentos canónicos conjugados de las cuatro compo-
nentes A µ

∂L ′
πµ = (6.50)
∂(∂0 A µ )

es decir

π0 = −∂µ A µ
πi = F i 0 = ∂0 A i − ∂i A 0 = E i (6.51)

Como vemos coincide con los anteriores para las componentes espaciales πi pero ahora
π0 es no nulo. Podemos intentar imponer las reglas de conmutación canónicas de forma
manifiestamente covariante

x ), πν (t ,⃗
[A µ (t ,⃗ y )] = i δνµ δ(3) (⃗
x −⃗
y) (6.52)

De forma que con todos los índice arriba se leería

[A µ (t ,⃗
x ), πν (t ,⃗
y )] = i ηµν δ(3) (⃗
x −⃗
y) (6.53)

y el resto de conmutadores nulos.

Puesto que todas las componentes cumplen ahora una ecuación de ondas, la solución
más general podrá escribirse como

d 3p 3 ³
Z
µ µ −i px µ∗
¢ † i px ´
ε ⃗ λ ε ⃗ λ
X ¡ ¢ ¡
A (x) = p , a p,λ e +⃗ p , a p,λ e (6.54)
(2π)3 2ωp λ=0 p 0 =ωp
p

donde ahora tenemos cuatro posible polarizaciones ϵµ (⃗ p , λ), λ = 0, . . . 3 y puesto que todas
las componentes satisfacen la misma ecuación de ondas tenemos para todos los modos
p 2 = 0. En el frame en el que p µ = (p, 0, 0, p), los vectores de polarización se puede elegir
como

ϵµ (⃗
p , 0) = (1, 0, 0, 0), ϵµ (⃗
p , 1) = (0, 1, 0, 0),
ϵµ (⃗
p , 2) = (0, 0, 1, 0), ϵµ (⃗
p , 3) = (0, 0, 0, 1) (6.55)

Antonio L. Maroto, Ignazio Scimemi 81

Campos Cuánticos 2022-23 UCM

µ
o de forma más compacta ϵµ (⃗p , λ) = δ λ . En cualquier otro inercial pueden obtenerse real-
izando la correspondiente transformación de Lorentz.
En general se cumple que

η µν ϵµ (⃗
p , λ)ϵν∗ (⃗
p , λ′ ) = ηλλ
′
(6.56)

ϵµ (⃗
p , λ)ϵν∗ (⃗
p , λ)
= ηµν
X
µ ν∗
(6.57)
λ η µν ϵ (⃗ p , λ)ϵ (⃗ p , λ)

independientemente del inercial3 .

Las relaciones de conmutación anteriores puede obtenerse imponiendo las siguientes
relaciones de conmutación sobre los operadores de creación y destrucción

† 3 (3)
[a p,λ , a q,λ′ ] = −η λλ′ (2π) δ p −⃗
(⃗ q) (6.59)
† †
[a p,λ , a q,λ′ ] = [a p,λ , a q,λ′ ] = 0, λ, λ′ = 0 . . . 3 (6.60)

El punto fundamental aquí es que el conmutador correspodiente a λ = λ′ = 0 tiene el

signo opuesto al de los espaciales y esto hace que los estados de una partícula correspon-
dientes tengan normas negativas, lo cual podría ser catastrófico para la consistencia de la
teoría. En efecto consideremos los estados
q
†
p , λ〉 = 2ωp a p,λ
|⃗ |0〉 (6.61)

por tanto
† †
p , λ|⃗
〈⃗ p , λ〉 = 2ωp 〈0|a p,λ a p,λ |0〉 = 2ωp 〈0|[a p,λ , a p,λ ]|0〉 = −η λλ′ 2ωp (2π)3 δ(3) (0) (6.62)

Puesto que los productos escalares en Mecánica Cuántica tienen una interpretación en tér-
minos de probabilidades, un espacio de Fock con un producto escalar no definido positivo
no tiene una interpretación probabilística. Por otra parte, hay que notar que ni los estados
† †
creados por a p,0 ni los creados por a p,3 correspondían a estados físicos.
3
Nótese que esta expresión también puede escribirse
X λλ µ
η ϵ (⃗ p , λ)ϵν∗ (⃗
p , λ) = ηµν (6.58)
λ

donde ηλλ es simplemente un signo.

Antonio L. Maroto, Ignazio Scimemi 82

Campos Cuánticos 2022-23 UCM

La solución a estos problemas puede encontrarse en el hecho de que hemos añadido un

término (∂µ A µ )2 al Lagrangiano. La idea de la cuantización covariante es limitar el espacio
de Fock de forma que sobre los estados físicos |phys〉 tengamos

〈phys’|(∂µ A µ )|phys〉 = 0 (6.63)

que es lo que se conoce como condición de Gupta-Bleuler. Por tanto, sobre estados físicos
se recuperarán las ecuaciones de Maxwell usuales. Veamos cómo implementar esta condi-
ción de forma más concreta. Para ello vamos a separar el operador ∂µ A µ en su parte de
frecuencia positiva y negativa.

∂µ A µ = (∂µ A µ )+ + (∂µ A µ )− (6.64)

donde

d 3p 3 ³
Z ´
µ + µ −i px
ε ⃗ λ
X ¡ ¢
(∂µ A ) = −i p µ p , a p,λ e
(2π)3 2ωp λ=0
p

d 3p 3 ³
Z
µ − µ∗
¢ † i px ´
ε ⃗ λ
X ¡
(∂µ A ) = i p µ p , a p,λ e (6.65)
(2π)3 2ωp λ=0
p

Puesto que (∂µ A µ )− = (∂µ A µ )+ tenemos que la condición de Gupta-Bleuler se cumple si

†

(∂µ A µ )+ |phys〉 = 0 (6.66)

Puesto que es una condición lineal, tendremos que el espacio de Fock de los estados físi-
cos también lo será, es decir si |phys1 〉 y |phys2 〉 son estados físicos entonces α|phys1 〉 +
β|phys2 〉 también lo será.
Veamos el efecto de esta condición sobre los estados de una partícula. Para ello con-
p †
sideremos un estado de una partícula general |ψ〉 = 3λ′ =0 c λ′ 2ωk a k,λ
P
′ |0〉. Tomando por
simplicidad k µ = (k, 0, 0, k), tenemos que

d 3p 3 ³ 3
Z ´X
(∂µ A µ )+ |ψ〉 = −i µ
⃗ −i px †
p
ε λ
X ¡ ¢
3
p p µ p , a p,λ e c λ′ 2ωk a k,λ′ |0〉
(2π) 2ωp λ=0 λ′ =0
3 ³
µ
´
k µ δλ e −i kx c λ′ η λ,λ′ |0〉 = i k(c 0 + c 3 )e −i kx |0〉 = 0
X
=i (6.67)
λλ′ =0

es decir, la condición de Gupta-Bleuler requiere c 0 + c 3 = 0. Es decir:

† †
|φ〉 = (a k,0 − a k,3 )|0〉 (6.68)

Antonio L. Maroto, Ignazio Scimemi 83

Campos Cuánticos 2022-23 UCM

sería un estado físico. Por otra parte, nótese que los estados de fotones transversales de la
† †
forma a k,1,2 |0〉 también son físicos. Sin embargo, ni los fotones temporales a k,0 |0〉 ni los
†
fotones longitudinales a k,3 |0〉 son físicos, sólo la combinación anterior lo es. Por tanto, la
forma más general de los estados de una partícula sería

|ψ〉 = |ψT 〉 + c|φ〉 (6.69)

con |ψT 〉 un estado de fotones transerversos.

Calculemos ahora la norma del estado físico |φ〉,
† † † †
〈φ|φ〉 = 〈0|(a k,0 − a k,3 )(a k,0 − a k,3 )|0〉 = 〈0|(a k,0 a k,0 + a k,3 a k,3 )|0〉
† †
= 〈0|([a k,0 , a k,0 ] + [a k,3 , a k,3 ])|0〉 = 0 (6.70)

debido a la diferencia de signo en los conmutadores. Es decir la norma de los estados

|φ〉 es cero y son por tanto ortogonales a todos los estados físicos. Por tanto, todos los
productos escalares de estados |ψT 〉 + c|φ〉 son idénticos a los productos escalares de los
estados transversos |ψT 〉.
Calculemos ahora la energía y el momento de los nuevos estados. La expresión del
Hamiltoniano es ahora diferente
à !
d 3p
Z
† †
ωp −a p,0
X
: H := a p,0 + a p,λ a p,λ (6.71)
(2π)3 λ=1,2,3

y de forma análoga
à !
d 3p
Z
⃗ := †
⃗ −a p,0 a p,0 + †
X
:P p a p,λ a p,λ (6.72)
(2π)3 λ=1,2,3

vemos por tanto que los fotones temporales contribuyen a la energía y al momento con un
signo opuesto al de las componentes transversas. De hecho puede verse que puesto que
los estados físicos contienen el mismo número de fotones temporales y longitudinales
sus contribuciones se cancelan y por tanto sobre los estados físicos únicamente los fotones
transversos contribuyen a la energía y el momento.
En conclusión, a pesar de que en el formalismo covariante aparecen nuevos estados,
estos no son físicos, puesto que no contribuyen a la energía ni al momento, ni a los pro-
ductos escalares. Dado su carácter manifiestamente covariante, este es el formalismo ha-
bitualmente empleado en los cálculos en electrodinámica cuántica.

Antonio L. Maroto, Ignazio Scimemi 84

Campos Cuánticos 2022-23 UCM

6.2.4 Propagador del fotón en la formulación covariante

Como en los casos anteriores comencemos calculando

d 3p d 3q
Z
〈0|A µ (x)A ν (y)|0〉 =
X µ¡ ¢ ν∗ †
3
p 3
p ε ⃗
p , λ ε (⃗ q , λ′ )e −i px e i q y 〈0|a p,λ a q,λ′ |0〉
(2π) 2ωp (2π) 2ωq λ,λ′
d 3p d 3p
Z X λλ µ ¡ Z
¢ ν∗
=− 3
η ε ⃗
p , λ ε p
(⃗ , λ)e −i p(x−y)
= − 3
ηµν e −i p(x−y) (6.73)
(2π) 2ωp λ (2π) 2ωp

Podemos por tanto calcular ahora el propagador retardado

d 3q
Z
µν
³ ´
D R (x − y) = θ(x 0 − y 0 )〈0|[A µ (x), A ν (y)]|0〉 = −θ(x 0 − y 0 ) ηµν e −i p(x−y) − e −i p(y−x)
(2π)3 2ωq
d 4 p ηµν −i p(x−y)
Z
=− i e (6.74)
(2π)4 p 2

donde el contorno en la integración en p 0 corresponde a la Fig. 4.1. De forma totalmente

análoga, la expresión para el propagador de Feynman es

d 4 p −i ηµν −i p(x−y)
Z
D µν (x − y) = 〈0|T A µ (x)A ν (y)|0〉 = e (6.75)
(2π)4 p 2 + i ϵ

donde el contorno de integración es el de la Fig. 4.2. Por tanto en espacio de Fourier el

propagador de Feynman covariante del fotón es

−i ηµν
D µν (p) = (6.76)
p2 + i ϵ

Esta expresión es la que usaremos como parte de las reglas de Feynman en QED.

Antonio L. Maroto, Ignazio Scimemi 85

7. Procesos de interacción

En los temas anteriores hemos estudiado cómo cuantizar teorías libres. En una teoría li-
bre, la evolución de los estados es trivial. En efecto, consideremos por simplicidad un es-
p 〉. En imagen de Schrödinger, el estado simplemente evolu-
tado inicial de una partícula |⃗
cionará con una fase, esto es

e −i H (t −t0 ) |⃗
p 〉 = e −i E p (t −t0 ) |⃗
p〉

donde hemos usado que los estados de Fock son autoestados del Hamiltoniano, y por tanto
el estado final será también un estado de una partícula con el mismo momento p⃗. Es decir,
en la teoría libre no es posible que ocurran procesos de dispersión en los que cambie el
momento de la partícula, ni de creación o aniquilación de partículas.
Los procesos físicamente interesantes aparecen cuando hay interacciones entre los
campos. Desde el punto de vista matemático, las interacciones se describen a través de
términos que involucren más de dos campos en el Lagrangiano. Dichos términos pueden
contener potencias mayores de un mismo campo, como

λ 4
φ
4!
en el caso de una teoría escalar, o productos que involucren distintos tipos de campos
como
−q ψ̄γµ ψA µ
en electrodinámica.
En este tema presentaremos el formalismo para calcular las probabilidades de transi-
ción entre distintos estados en presencia de interacciones y cómo, a partir de esas probabil-
idades, calcular las cantidades que pueden medirse experimentalmente (ritmos de desin-
tegración y secciones eficaces).

86
Campos Cuánticos 2022-23 UCM

7.1 Imagen de interacción

Consideremos una teoría con interacción, en la que podemos separar el Hamiltoniano del
sistema en dos partes

H = H0 + Hi nt (7.1)

donde H0 corresponde al Hamiltoniano de la teoría libre que hemos visto en los temas
anteriores y Hi nt es el Hamiltoniano de las interacciones. Así, por ejemplo en la teoría
escalar con interacción
1 1 λ
L = ηµν ∂µ φ∂ν φ − m 2 φ2 − φ4 (7.2)
2 2 4!
podemos definir

λ 4
Z Z
Hi nt = − d 3 x L i nt = d 3x φ (7.3)
4!

donde hemos usado que el término de interacción no involucra derivadas de los campos.
Supondremos que Hi nt es una pequeña perturbación con respecto a H0 .
Como vimos en temas anteriores, en Mecánica Cuántica es posible trabajar en la

Imagen de Schrödinger

• A S . Los operadores no dependen del tiempo

• |Ψ (t )〉S = e −i H (t −t0 ) |Ψ (t 0 )〉S . Los estados evolucionan.

Imagen de Heisenberg

• A H (t ) = e i H (t −t0 ) A H (t 0 )e −i H (t −t0 ) . Los operadores evolucionan.

• |Ψ〉H . Los estados no evolucionan

En el caso en el que tengamos interacciones es posible definir una tercera imagen en la

que tanto los estados como los operadores dependen del tiempo

Antonio L. Maroto, Ignazio Scimemi 87

Campos Cuánticos 2022-23 UCM

Imagen de interacción

• A I (t ) = e i H0 (t −t0 ) A I (t 0 )e −i H0 (t −t0 ) . Los operadores evolucionan con la parte libre del

Hamiltoniano, es decir, como los operadores libres.

• |Ψ (t )〉I = e i H0 (t −t0 ) e −i H (t −t0 ) |Ψ (t 0 )〉I . Los estados evolucionan con el hamiltoniano

libre y el total.

Nótese que A I (t 0 ) = A H (t 0 ) = A S y que |Ψ (t 0 )〉I = |Ψ (t 0 )〉S = |Ψ〉H . Nótese también que si

Hi nt = 0, las imágenes de Heisenberg y de interacción coinciden.
Podemos, por tanto, definir el operador de evolución en la imagen de interacción

|Ψ (t )〉I = U I (t , t 0 )|Ψ (t 0 )〉I (7.4)

como

U I (t , t 0 ) = e i H0 (t −t0 ) e −i H (t −t0 ) (7.5)

que cumple U I (t 0 , t 0 ) = 1 y satisface la ecuación

∂
i U I (t , t 0 ) = H I (t )U I (t , t 0 ) (7.6)
∂t

donde H I (t ) = e i H0 (t −t0 ) Hi nt e −i H0 (t −t0 ) . Así, en el ejemplo anterior

λ
Z
H I (t ) = d 3 x φ4I (7.7)
4!
La solución de la ecuación (7.6) puede escribirse formalmante como
Z t
U I (t , t 0 ) = 1 + (−i ) d t ′ H I (t ′ )U I (t ′ , t 0 ) (7.8)
t0

utilizando el hecho de que H I es una pequeña perturbación, podemos escribir una solu-
ción en serie de la forma
Z t Z t Z t1
2
U I (t , t 0 ) = 1 + (−i ) d t 1 H I (t 1 ) + (−i ) d t1 d t 2 H I (t 1 )H I (t 2 )
t0 t0 t0
Z t Z t1 Z t2
+ (−i )3 d t1 d t2 d t 3 H I (t 1 )H I (t 2 )H I (t 3 ) + . . . (7.9)
t0 t0 t0

Se puede comprobar que al derivar cada término se obtiene el anterior multiplicado por
−i H I (t ) y por tanto es solución perturbativa de la ecuación (7.6). Nótese que en la integral

Antonio L. Maroto, Ignazio Scimemi 88

Campos Cuánticos 2022-23 UCM

Figure 7.1: Integración temporal en (7.10).

anterior los factores H I (t i ) aparecen en orden temporal (los últimos a la izquierda) puesto
que por ejemplo t 2 ≤ t 1 . Esto implica que por ejemplo el término H I2
Z t Z t1 1
Z t Z t
d t1 d t 2 H I (t 1 )H I (t 2 ) = d t1 d t 2 T {H I (t 1 )H I (t 2 )} (7.10)
t0 t0 2 t0 t0

donde T denota el producto cronológico de los operadores, como puede verse en la Fig.7.1.
Expresiones similares pueden encontrarse para el resto de términos de forma que pode-
mos escribir

(−i )2 t
Z t Z Z t
U I (t , t 0 ) = 1 + (−i ) d t 1 H I (t 1 ) + d t1 d t 2 T {H I (t 1 )H I (t 2 )} + . . .
t0 2! t0 t0
µ Z t ¶
′ ′
= T exp −i d t H I (t ) (7.11)
t0

que es la llamada fórmula de Dyson.

Antonio L. Maroto, Ignazio Scimemi 89

Campos Cuánticos 2022-23 UCM

7.2 La matriz S
Hemos visto en la sección anterior que debido a la presencia de la interacción, la evolu-
ción de los estados es no trivial. Consideremos por tanto un sistema que inicialmente se
encuentra en un cierto estado |i 〉I y queremos determinar la probabilidad de que se trans-
forme en un estado final | f 〉I debido a la interacción. Es decir, si el estado de nuestro
sistema viene descrito en la imagen de interacción por |ψ(t )〉I , tendremos

lim |ψ(t )〉I = |i 〉I (7.12)

t →−∞

Por tanto, la amplitud de probabilidad de que este estado se transforme en un cierto estado
| f 〉I en t → ∞ vendrá dado por el elemento de matriz S

S f i = lim I 〈 f |ψ(t )〉I = I 〈 f |S|i 〉I (7.13)

t →∞

Expresado en términos del operador de evolución en la imagen de interacción tendremos

S f i = lim lim I〈f |U I (t 2 , t 1 )|i 〉I (7.14)

t 2 →∞ t 1 →−∞

es decir el operador S puede escribirse como

µ Z ∞ ¶ µ Z ¶
′ ′ 4
S = U I (∞, −∞) = T exp −i d t H I (t ) = T exp −i d xH I (x) (7.15)
−∞

Un aspecto fundamental de este resultado es que puesto que los campos en imagen de
interacción evolucionan como campos libres, podremos usar la descomposición en oper-
adores de creacion y destrucción vista en los temas anteriores cuando calculemos H I .
El operador S es un operador unitario S † S = SS † = 1 lo que representa la conservación
de la probabilidad. Es conveniente introducir la llamada matriz T de la forma

S = 1+iT (7.16)

En los procesos de dispersión suponemos que en los estados iniciales y finales las partícu-
las están suficientemente separadas como para ignorar las interacciones. Es decir, los es-
tados inicial y final |i 〉I y | f 〉I se supondrán por simplicidad como estados libres, es decir
serán autoestados del Hamiltoniano libre H0 . Esto plantea dos problemas:

• Este formalismo no permite trabajar con estados ligados. Por ejemplo, no podremos
describir el proceso en el que un electrón y un protón colisionan y forma un átomo
de hidrógeno.

Antonio L. Maroto, Ignazio Scimemi 90

Campos Cuánticos 2022-23 UCM

• Aún más importante es que incluso una única partícula alejada de otras partículas
nunca puede considerarse de forma estricta como libre en teoría de campos. Así
por ejemplo, un electrón siempre aparecerá rodeado por un nube de fotones vir-
tuales. Esto conduce a la noción de renormalización que no veremos en este curso.
La forma rigurosa de tratar con este problema se basa en la llamada fórmula de re-
ducción de Lehamnn-Symanzik-Zimmermann (LSZ).

7.2.1 Ejemplo: desintegración de mesones

Consideremos como ejemplo la teoría de Yukawa escalar1 para un campo escalar real φ
(mesón neutro) de masa m acoplado a un escalar complejo ϕ (nucleón cargado) de masa
M . El Lagrangiano total será
1 1
L = ηµν ∂µ φ∂ν φ − m 2 φ2 + ηµν ∂µ ϕ∗ ∂ν ϕ − M 2 ϕ∗ ϕ − g ϕ† ϕφ (7.17)
2 2
donde el Lagrangiano de interacción será

L i nt = −g ϕ† ϕφ

Con g una constante con dimensiones de masa. Por tanto el Hamiltoniano de interacción
será Z
Hi nt = g d 3 x ϕ† ϕφ

y el operador de la matriz S
µ Z ¶
4 †
S = T exp −i g d xϕ ϕφ (7.18)

donde hemos quitado el subíndice I de los campos por simplicidad. Suponiendo que la
constante de acoplo g es pequeña comparadas con m y M , podemos quedarnos con el
primer término en el desarrollo perturbativo, es decir
Z
S = 1 − i g d 4 xϕ† ϕφ + . . . (7.19)

Consideremos ahora el proceso de desintegración de un mesón φ en un par nucleón-

antinucleón.
φ → ϕϕ̄
1
Se trata de una versión simplificada de la teoría original de Yukawa que acoplaba escalares reales con
fermiones de Dirac y por tanto no representa procesos reales.

Antonio L. Maroto, Ignazio Scimemi 91

Campos Cuánticos 2022-23 UCM

Por tanto, los estados inicial y final serán:

q
|i 〉 = 2E p a p† |0〉 (7.20)
q
|f 〉 = 4E q1 E q2 b q† 1 c q† 2 |0〉 (7.21)

donde para los campos en la imagen de interacción podemos utilizar la descomposición

usual en operadores de creación y destrucción libres φ ∼ a + a † , ϕ ∼ b +c † y ϕ† ∼ c +b † . Por
tanto, la amplitud de transición será
Z
S f i = −i g 〈 f | d 4 xϕ† (x)ϕ(x)φ(x)|i 〉 (7.22)

Nótese que puesto que |i 〉 ̸= | f 〉 el primer término del desarrollo es cero. Al sustituir el
campo φ en términos de operadores de creación y destrucción, el operador a † actuando
sobre |i 〉 dará lugar a un estado con dos mesones. Puesto que el estado final no contiene
mesones, el producto dará cero. Luego nos quedamos sólo con el término de aniquilación
p
d 3k 2E p
Z Z
S f i = −i g 〈 f | 4 †
d xϕ (x)ϕ(x) 3
p a k a p† e −i kx |0〉
(2π) 2E k
Z
= −i g 〈 f | d 4 xϕ† (x)ϕ(x)e −i px |0〉 (7.23)

Ahora sólo los términos con b † y c † contribuirán al actuar sobre el vacío, puesto que | f 〉
contiene un nucleón y un antinucleón, por tanto:
p
d 3 k 1 d 3 k 2 4E q1 E q2
Z Z
S f i = −i g 〈0| d 4x c q2 b q1 c k† b k† |0〉e i (k1 +k2 −p)x
(2π)6
p 1 2
4E k1 E k2
= −i g (2π)4 δ(4) (p − q 1 − q 2 ) (7.24)

que proporciona la amplitud de desintegración de un mesón con momento p en un par

nucleón-antinucleón con momentos q 1 y q 2 . Vemos que la delta de Dirac impone la con-
servación del cuadrimomento. Más adelante veremos cómo convertir esta amplitud en
una cantidad medible como el tiempo de vida media del mesón.
Una representación gráfica de la amplitud anterior está dada en la Fig. 7.2 donde rep-
resentamos:

• El mesón inicial φ con una línea discontinua.

• El nucleón final ϕ como una línea continua con una flecha saliente.

Antonio L. Maroto, Ignazio Scimemi 92

Campos Cuánticos 2022-23 UCM

Figure 7.2: Diagrama de Feynman correspondientes a la desintegración de un mesón.

• El antinucleón final ϕ̄ como una línea continua con una flecha entrante.

Las flechas en las líneas de las partículas finales indican el sentido del flujo de la carga (la
línea del mesón neutro no lleva flecha). Este gráfico consitituye el diagrama de Feynman
del proceso de desintegración.
Este diagrama es el caso más sencillo en el que simplemente aparece el vértice de tres
partículas que contiene L i nt . Este diagrama contribuye a la amplitud con un factor −i g
por el vértice y un factor (2π)4 δ(4) (p − q 1 − q 2 ) donde aparece la suma de todos los cuadri-
momentos entrantes menos los cuadrimomentos salientes del vértice.
Aunque este cálculo lo hemos realizado en un ejemplo concreto, el procedimiento es
completamente general (para interacciones que no involucren derivadas). Es decir, dada
una densidad Lagrangiana de interacción L i nt podemos calcular la contribución del vér-
tice de interacción de la siguiente forma:

• Tomar un estado inicial |i 〉 y un estado final | f 〉 que involucren en total tantas partícu-
las y/o antipartículas como campos aparecen en L i nt . En el ejemplo anterior ten-
emos una partícula φ en el estado inicial y un partícula y una antipartícula ϕ en el
estado final.

• Calcular 〈 f |i d 4 xL i nt |i 〉
R

• El resultado se podrá escribir como i Mi f (2π)4 δ(4) (P i − P f ) donde P i es la suma de

los cuadrimomentos de las partículas iniciales y P f el de las partículas finales.

Antonio L. Maroto, Ignazio Scimemi 93

Campos Cuánticos 2022-23 UCM

7.3 El teorema de Wick y los diagramas de Feynman

La fórmula de Dyson requiere el cálculo de términos que en general tienen la forma

〈 f |T {H I (x 1 )H I (x 2 ) . . . H I (x n )}|i 〉

Su cálculo podría simplificarse si los operadores de aniquilación que aparecen en los H I

pudieran moverse a la derecha como en el caso de la ordenación normal y actuar sobre el
estado |i 〉. Precisamente el teorema de Wick nos dice cómo pasar de los productos orde-
nados temporalmente a productos ordenados normalmente.
Consideremos un ejemplo sencillo con campos escalares reales:

T {φ(x)φ(y)} (7.25)

donde hemos suprimido el subíndice I por simplicidad.

Para calcular el producto temporal separamos el operador φ(x) en sus partes de fre-
cuencia positiva y negativa

φ(x) = φ+ (x) + φ− (x) (7.26)

donde
d 3p
Z
+
φ (x) = a p e −i px (7.27)
(2π)3 2E p
p

d 3p
Z
−
φ (x) = a p† e i px (7.28)
3
p
(2π) 2E p

Así, si x 0 > y 0 tendremos

T {φ(x)φ(y)} = φ(x)φ(y) (7.29)

+ − + −
= (φ (x) + φ (x))(φ (y) + φ (y))
= φ+ (x)φ+ (y) + φ− (x)φ− (y) + φ− (x)φ+ (y) + φ+ (x)φ− (y) (7.30)

donde todos los términos salvo el último están ordenados normalmente. Podemos en-
tonces reescribir el último usando el conmutador, de forma

T {φ(x)φ(y)} =: φ(x)φ(y) : +[φ+ (x), φ− (y)] (7.31)

De forma totalmente análoga si y 0 < x 0 tendremos

T {φ(x)φ(y)} =: φ(x)φ(y) : +[φ+ (y), φ− (x)] (7.32)

Antonio L. Maroto, Ignazio Scimemi 94

Campos Cuánticos 2022-23 UCM

y por tanto, en general

T {φ(x)φ(y)} =: φ(x)φ(y) : +θ(x 0 − y 0 )[φ+ (x), φ− (y)] + θ(y 0 − x 0 )[φ+ (y), φ− (x)]
=: φ(x)φ(y) : +∆(x, y) (7.33)

donde ∆(x, y) no es más que el propagador de Feynman (4.63). Es decir, el producto cronó-
logico de dos campos puede escribirse como el producto normal más el propagador de
Feynman.
Podemos extender este resultado al caso general con n campos. Consideremos la ca-
dena φ1 φ2 φ3 φ4 . Definimos la contracción de dos operadores como la operación que con-
vierte el par de operadores en un propagador de Feynman y deja los demás operadores
igual. Esto es

φ1 φ2 φ3 φ4 = φ2 φ4 ∆(x 1 − x 3 ) (7.34)

En el caso de campos complejos φ1 φ†2 = ∆(x 1 − x 2 ), con φ1 φ2 = 0. En el caso fermiónico

tendríamos ψ1 ψ̄2 = S F (x 1 − x 2 ), con un signo menos cuando se necesite intercambiar dos
operadores fermiónicos para poner los dos operadores a contraer juntos. Así por ejemplo
ψ1 ψ2 ψ̄3 = −S F (x 1 − x 3 )ψ2 y ψ1 ψ2 = 0.
El teorema de Wick2 afirma que

T (φ1 φ2 . . . φn ) =: φ1 φ2 . . . φn : + : todas las contracciones posibles : (7.35)

Así por ejemplo:

T (φ1 φ2 φ3 ) =: φ1 φ2 φ3 : + : φ1 φ2 φ3 : + : φ1 φ2 φ3 : + : φ1 φ2 φ3 :
=: φ1 φ2 φ3 : +φ3 ∆(x 1 − x 2 ) + φ1 ∆(x 2 − x 3 ) + φ2 ∆(x 1 − x 3 ) (7.36)

En el caso de cuatro campos, serían posibles contracciones simples del tipo φ1 φ2 φ3 φ4 y

dobles como φ1 φ2 φ3 φ4 .
2
El teorema de Wick también puede formularse directamente para operadores de creación de destrucción,
de forma que
ABC D · · · =: ABC D · · · : + : todas las contracciones posibles :

donde A, B, ... son operadores de creación o destrucción. La contracción en este caso se define como AB =
AB − : AB :. Así, para el caso escalar a p a †q = [a p , a q† ] = (2π)3 δ(3) (⃗ q ), pero a †q a p = 0
p −⃗

Antonio L. Maroto, Ignazio Scimemi 95

Campos Cuánticos 2022-23 UCM

7.3.1 Ejemplo: dispersión de nucleones

Consideremos de nuevo la interacción de Yukawa. Estamos interesados en calcular la am-
plitud de probabilidad para la dispersión de nucleones, esto es

ϕϕ → ϕϕ

Por tanto, los estados iniciales y finales ahora serán

q
|i 〉 = 4E p 1 E p 2 b p† 1 b p† 2 |0〉 = |⃗
p1, p ⃗2 〉 (7.37)
q
| f 〉 = 4E p ′ E p ′ b p† ′ b p† ′ |0〉 = |⃗
p 1′ , p
⃗2′ 〉 (7.38)
1 2 1 2

De nuevo queremos calcular 〈 f |S|i 〉 aunque como no nos interesa la parte sin dispersión,
en realidad vamos a calcular 〈 f |S − 1|i 〉 = 〈 f |i T |i 〉. Nótese que el primer término del de-
sarrollo del operador de Dyson O (g ) no permite conectar un estado de dos nucleones ini-
ciales con un estado de dos nucleones finales, puesto que eso requiere al menos cuatro
operadores ϕ. Por tanto, este proceso, sólo puede tener lugar al siguiente orden en teoría
de perturbaciones, esto es O (g 2 ). El término a considerar será entonces:

(−i g )2
Z
iTf i = 〈f | d 4 x 1 d 4 x 2 T {ϕ† (x 1 )ϕ(x 1 )φ(x 1 )ϕ† (x 2 )ϕ(x 2 )φ(x 2 )}|i 〉 (7.39)
2
Usando el teorema de Wick, podemos desarrollar el producto ordenado temporalmente y
encontramos que el único término que permite conectar el estado inicial y final es

: ϕ† (x 1 )ϕ(x 1 )ϕ† (x 2 )ϕ(x 2 ) : φ(x 1 )φ(x 2 ) (7.40)

Ahora bien recordemos que en el producto normal tenemos los aniquiladores a la derecha
y que ϕ ∼ b + c † y ϕ† ∼ c + b † por lo que podemos escribir (ignorando por ahora la contrac-
ción de los dos últimos operadores que sale del producto)

〈⃗ ⃗2′ | : ϕ† (x 1 )ϕ(x 1 )ϕ† (x 2 )ϕ(x 2 ) : |⃗

p 1′ , p ⃗2 〉 = 〈⃗
p1, p ⃗2′ |ϕ† (x 1 )ϕ† (x 2 )|0〉〈0|ϕ(x 1 )ϕ(x 2 )|⃗
p 1′ , p ⃗2 〉
p1, p
(7.41)

p 〉 = e −i px obtenemos
Usando ahora (4.53) 〈0|ϕ(x)|⃗

〈p 1′ , p 2′ | : ϕ† (x 1 )ϕ(x 1 )ϕ† (x 2 )ϕ(x 2 ) : |p 1 , p 2 〉

³ ′ ′ ′ ′
´³ ´
= e i (p 1 x1 +p 2 x2 ) + e i (p 1 x2 +p 2 x1 ) e −i (p 1 x1 +p 2 x2 ) + e −i (p 1 x2 +p 2 x1 )
′ ′ ′ ′
= e i x1 (p 1 −p 1 )+i x2 (p 2 −p 2 ) + e i x1 (p 2 −p 1 )+i x2 (p 1 −p 2 ) + (x 1 ↔ x 2 ) (7.42)

Antonio L. Maroto, Ignazio Scimemi 96

Campos Cuánticos 2022-23 UCM

Por tanto,

(−i g )2 d 4 k i e i k(x1 −x2 )

Z Z
4 4
iTf i = d x1 d x2 [ ] (7.43)
2 (2π)4 k 2 − m 2 + i ϵ

donde [ ] es la expresión (7.42) y la última integral es el término de la contracción φ(x 1 )φ(x 2 ) =

∆(x 1 −x 2 ) donde m es la masa del mesón φ. Realizando las integrales en x 1 y x 2 obtenemos
deltas de Dirac y finalmente llegamos a

d 4k i (2π)8
Z
2
£ (4)
i T f i = (−i g ) δ (p 1 − p 1′ − k)δ(4) (p 2 − p 2′ + k)
(2π) k − m + i ϵ
4 2 2

+δ(4) (p 1 − p 2′ − k)δ(4) (p 2 − p 1′ + k)
¤
(7.44)
i i
· ¸
= (−i g )2 ′ 2 2
+ ′ 2 (2π)4 δ(4) (p 1 + p 2 − p 1′ − p 2′ ) (7.45)
(p 1 − p 1 ) − m (p 1 − p 2 ) − m 2

El elemento de matriz T que acabamos de obtener puede entenderse diagramáticamente

como la suma de los dos diagramas en la Fig. 7.3, si les asociamos las siguientes reglas:

• Por cada vértice un factor (−i g )(2π)4 δ(4) ( i k i ) donde aparece la suma de todos los
P

cuadrimomentos entrantes en el vértice menos los salientes.

• Por cada línea interna de mesones un factor

d 4k i
Z

(2π) k − m 2 + i ϵ
4 2

Vemos efectivamente que aplicando estas reglas a los diagramas anteriores se obtiene la
expresión (7.44).

7.3.2 Diagramas de Feynman

Aunque hemos considerado únicamente un par de ejemplos concretos, se puede demostrar
de forma totalmente general que existe una relación uno a uno entre los distintos térmi-
nos del desarrollo perturbativo de la amplitud de transición y los diagramas de Feyn-
man.
Por tanto, para calcular i T f i = 〈 f |i T |i 〉 tendremos que dibujar todos los posibles dia-
gramas de Feynman (conectados) que conecten las partículas de los estados iniciales con
las finales. Para ello tendremos que seguir los siguientes pasos:

• Identificar los distintos vértices de interacción que aparecen en L i nt . En el ejemplo

anterior sólo había un vértice de tres partículas (mesón-(anti)barión-(anti)barión).

Antonio L. Maroto, Ignazio Scimemi 97

Campos Cuánticos 2022-23 UCM

Figure 7.3: Diagramas de Feynman correspondientes a la dispersión nucleón-nucleón.

• Dibujar una línea externa por cada partícula del estado inicial |i 〉 y por cada partícula
del estado final | f 〉. En el ejemplo anterior se usaban líneas discontinuas para los
mesones y líneas continuas para los bariones y antibariones.

• Asignar un momento con sentido entrante o saliente a cada línea.

• En las líneas correspondiente a partículas con carga conservada (no necesariamente

eléctrica) asignar una flecha que indique el sentido de la carga. Entrante (saliente)
para partículas (antipartículas) en el estado inicial y al revés en el estado final.

• Conectar de todas las formas posibles las líneas de las partículas iniciales con las fi-
nales utilizando los vértices y líneas internas. Asignar también un momento entrante
o saliente del vértice a cada línea interna.

Una vez dibujados todos los diagramas de Feynman que contribuyen a un cierto pro-
ceso, calcularemos su contribución a la amplitud de transición usando las siguientes re-
glas:

7.3.3 Reglas de Feynman

Hemos visto que en los elementos de matriz i T f i que hemos calculado aparecen factores
δ(4) (P i − P f ), donde P i = k j y P f = p k son los cuadrimomentos totales iniciales y fi-
P P

nales, que aseguran la conservación del cuadrimomento. Es conveniente por tanto extraer

Antonio L. Maroto, Ignazio Scimemi 98

Campos Cuánticos 2022-23 UCM

estos factores y escribir los elementos de matriz T como

i T f i = (2π)4 δ(4) (P i − P f )i M f i (7.46)

donde se ha definido el llamado elemento de matriz invariante M f i . Podemos finalmente

escribir unas reglas de Feynman aún más sencillas para calcular las amplitudes invariantes
i M f i . Consideraremos el caso de la teoría escalar de Yukawa.

• Escribir un factor (−i g ) por cada vértice.

• Imponer conservación del cuadrimomento en cada vértice.

• Escribir un propagador de Feynman por cada línea interna, esto es,

i
k 2 − m2 + i ϵ
por cada línea interna de mesones y

i
k2 − M 2 + i ϵ
por cada línea interna de nucleones.

• Integrar sobre el momento interno que fluya en líneas cerradas (loops)

d 4k
Z

(2π)4

Respecto al último punto, los diagramas que hemos visto son diagramas a nivel árbol (tree-
level) que no incluyen loops. Típicamente las contribuciones de los loops son divergentes
y requieren un proceso de regularización y renormalización para proporcionar resultados
finitos que no veremos en este curso.
Podemos considerar ahora el calculo de la amplitud invariantes i M a orden g 2 para la
dispersión nucleón-antinucleón ϕϕ̄ → ϕϕ̄ usando las reglas de Feynman anteriores. En la
Fig. 7.4 pueden verse los diagramas correspondientes que como vemos difieren de los de
la dispersión de nucleones. La amplitud invariante se obtiene inmediatamente

i i
· ¸
2
i M f i = (−i g ) + (7.47)
(p 1 − p 1′ )2 − m 2 (p 1 + p 2 )2 − m 2

Antonio L. Maroto, Ignazio Scimemi 99

Campos Cuánticos 2022-23 UCM

Figure 7.4: Diagramas de Feynman correspondientes a la dispersión nucleón-antinucleón.

7.4 Secciones eficaces y ritmos de desintegración

En las secciones anteriores hemos visto cómo calcular elementos de matriz entre estados
iniciales y finales. En esta sección calcularemos las cantidades que pueden ser directa-
mente medidas en experimentos, es decir, las secciones eficaces y los ritmos de desinte-
gración.
En la sección anterior hemos visto que la matriz S puede escribirse como

S = 1+iT (7.48)

Considerando un estado inicial |i 〉 = |⃗

k 1⃗
k 2 . . .⃗
k m 〉 y un estado final | f 〉 = |⃗ ⃗2 . . . p
p1p ⃗n 〉, el
término T se podía escribir en términos del elemento de matriz invariante M f i de forma
que

S f i = 〈 f |S|i 〉 = 〈 f |i 〉 + (2π)4 δ(4) (P i − P f )i M f i (7.49)

P P
donde P i = k j y P f = p k son los cuadrimomentos totales iniciales y finales.
La probabilidad de la transición del estado |i 〉 al | f 〉 vendrá dada por tanto por el mó-
dulo al cuadrado del elemento de matriz S, es decir

|〈 f |S|i 〉|2 (7.50)

Sin embargo, el cálculo de esta cantidad presenta algunos problemas. En primer lugar, para
interpretar (7.50) como una probabilidad, los estados deben estar correctamente normal-
izados, es decir 〈i |i 〉 = 〈 f | f 〉 = 1. Sin embargo, hemos visto en (4.46) que la normalización

Antonio L. Maroto, Ignazio Scimemi 100

Campos Cuánticos 2022-23 UCM

relativista de los estados de Fock no es ésta. Por lo que deberemos dividir entre la norma de
dichos estados. Por otra parte, puesto que los valores de los momentos de los estados son
continuos, la expresión anterior debe entenderse como una densidad de probabilidad en
el espacio de estados finales. Por tanto podemos definir la probabilidad de transición de
un estado inicial |i 〉 a un conjunto de estados finales con tri-momentos entre p⃗j y p
⃗ j +d p
⃗j ,
j = 1 . . . n como

|〈 f |S|i 〉|2
dP = dΠ (7.51)
〈 f | f 〉〈i |i 〉

donde d Π denota el volumen de estados de las partículas finales. El cálculo de esta ex-
presión plantea de nuevo otras dificultades debido a que la normalización relativista (4.46)
es divergente. Debemos por tanto dar sentido a esta expresión. Para ello es conveniente
pasar a realizar los cálculos en una caja de volumen finito V y luegoR tomar el límite V → ∞.
Nótese en primer lugar que de la definición de la delta de Dirac d 3 x e i p⃗⃗x = (2π)3 δ(3) (⃗
p ),
podemos escribir tomando p ⃗ = 0 en la expresión anterior y haciendo la integral en un vol-
umen finito V

V = (2π)3 δ(3) (0) (7.52)

Por tanto, los estados de Fock en un volumen finito estarán normalizados como

p |⃗
〈⃗ p 〉 = 2E p V (7.53)

Calculemos ahora el número de estados d Π con tri-momentos entre p ⃗j y p

⃗j + d p
⃗j , j =
1 . . . n. Consideremos por simplicidad una caja cúbica de lado L de forma que V = L 3 . Para
partículas en una caja, los momentos están cuantizados de forma que p ⃗ = (2π/L)(n x⃗ ex +
n y⃗ e y + n z⃗
e z ) con n x , n y , n z = 0, ±1, ±2 . . . (en unidades naturales h = 2π, condiciones de
contorno periódicas). Por tanto, la densidad de estados en espacio de momentos será n e =
(L/(2π))3 . Por lo tanto, el número de estados en una región del espacio de momentos con
volumen d 3 p será n e d 3 p = V /(2π)3 d 3 p para cada partícula, y por tanto para la n partículas
en el estado final, tendremos
n V
dΠ = d 3p j
Y
3
(7.54)
j =1 (2π)

Por otra parte, en el caso en el que los estados inicial y final sean distintos tendremos que

|〈 f |S|i 〉|2 = |(2π)4 δ(4) (P i − P f )i M f i |2 (7.55)

Antonio L. Maroto, Ignazio Scimemi 101

Campos Cuánticos 2022-23 UCM

Nótese que de forma análoga al caso de volumen espacial finito, si también trabajamos en
un intervalo temporal finito T , podemos escribir V T = (2π)4 δ(4) (0) y por tanto la expresión
anterior que involucra el producto de dos deltas de Dirac puede reescribirse como

|〈 f |S|i 〉|2 = (2π)4 δ(4) (P i − P f )V T |M f i |2 (7.56)

donde hemos usado (2π)8 δ(4) (P i − P f )2 = V T (2π)4 δ(4) (P i − P f ) en un volumen espacio-

temporal finito.

7.4.1 Ritmo de desintegración

Consideremos el caso de una partícula en el estado inicial con momento p ⃗ que se desin-
tegra en una serie de partículas con momentos p ⃗ j , j = 1 . . . n, de forma que |i 〉 = |⃗
p〉 y
podemos escribir la probabilidad de desintegración como

|〈 f |S|i 〉|2 (2π)4 δ(4) (P i − P f )V T |M f i |2 Y n V

dP = dΠ = Qn 3
d 3p j (7.57)
〈i |i 〉〈 f | f 〉 2E p V j =1 (2E j V ) j =1 (2π)

Podemos por tanto obtener para el ritmo de desintegración es decir, la probabilidad de

desintegración por unidad de tiempo d Γ = d P /T la siguiente expresión
|M f i |2 Y
à !
4 (4)
n d 3p j
d Γ = (2π) δ
X
p − pj (7.58)
j 2E p j =1 (2π)3 2E j

Se suele definir el llamado espacio de fases invariante Lorentz (LIPS) de n partículas como
la cantidad
à !
n d 3p j
d Φ(n) = (2π)4 δ(4) p − p j
X Y
3
(7.59)
j j =1 (2π) 2E j

De forma que el ritmo de desintegración puede escribirse como

|M f i |2
dΓ = d Φ(n) (7.60)
2E p
Γ se conoce también como anchura de desintegración y su inverso τ = 1/Γ es el tiempo
de vida media de la partícula.
Nótese que tanto M como el LIPS son invariantes Lorentz, sin embargo E p depende
del sistema de referencia y tiene en cuenta el efecto de dilatación temporal en las desinte-
graciones al cambiar de inercial. Finalmente, si en el estado final m de las partículas son
idénticas sus permutaciones serían indistinguibles y por tanto habría que dividir el LIPS
por m!.

Antonio L. Maroto, Ignazio Scimemi 102

Campos Cuánticos 2022-23 UCM

7.4.2 Secciones eficaces

Consideremos ahora el caso de dos partículas en el estado inicial, esto es |i 〉 = |⃗ k 1 ,⃗
k 2 〉 que
colisionan dando lugar a una serie de partículas con momentos p ⃗ j , j = 1 . . . n, en el estado
3
final. y nos interesa medir la sección eficaz del proceso, esto es, el número de colisiones
que tienen lugar por unidad de tiempo normalizado por el flujo de partículas que colisio-
nan. Esto es, la sección eficaz del proceso vendrá dada por

dP
dσ = (7.61)
TΦ
En el sistema LAB en el que la partícula 2 está en reposo, el flujo vendría dado simplemente
por la velocidad de la partícula 1 dividido entre el volumen Φ = |⃗ v 1 |/V . En el sistema centro
de masas (CM), las dos partículas contribuyen al flujo y por tanto Φ = |⃗ v1 − ⃗v 2 |/V donde
⃗
v =p⃗/E p . Por tanto, sustituyendo en (7.61)

dP |〈 f |S|i 〉|2 (2π)4 δ(4) (P i − P f )V T |M f i |2 n V

dσ = dΠ = d 3p j
Y
= Qn (7.62)
T Φ 〈i |i 〉〈 f | f 〉T Φ v1 − ⃗
|⃗ v 2 |(2E 1 )(2E 2 )V T j =1 (2E j V ) j =1 (2π)3

que podemos escribir en términos del LIPS como

|M f i |2
dσ = d Φ(n) (7.63)
v1 − ⃗
|⃗ v 2 |(2E 1 )(2E 2 )

Nótese que la sección eficaz no es invariante Lorentz puesto que la cantidad I = E 1 E 2 |⃗

v1 −
⃗
v 2 | cambia bajo transformaciones de Lorentz. Es decir, distintos observadores medirán
secciones eficaces diferentes para el mismo proceso. Sin embargo, se puede ver que si las
dos partículas son colineales, I es invariante bajo boosts a lo largo de la dirección de col-
isión4 . Por tanto, las secciones eficaces medidas en los sistemas CM y LAB son las mismas.

7.4.3 Estados finales de dos partículas

Consideremos ahora el caso particular de colisiones 2 → 2 entre partículas con momentos

p1 + p2 → p3 + p4 (7.64)
3
En realidad, en los experimentos de colisión tenemos dos haces de partículas que colisionan. En este
caso habría que multiplicar numerador y denominador por las densidades de partículas de los haces n 1 y n 2
pero el resultado sería el mismo
q
4
De hecho, en el caso colineal I coincide con la cantidad I col l = (p 1 µ p 2 µ )2 − (m 1 m 2 )2 que es invariante
Lorentz.

Antonio L. Maroto, Ignazio Scimemi 103

Campos Cuánticos 2022-23 UCM

donde las masas de las partículas iniciales son m 1 y m 2 y las de las finales m 3 y m 4 . Es útil
introducir las llamadas variables de Mandelstam

s = (p 1 + p 2 )2 (7.65)
t = (p 1 − p 3 )2 (7.66)
2
u = (p 1 − p 4 ) (7.67)

Las variables son invariants Lorentz y cumplen

s + t + u = m 12 + m 22 + m 32 + m 42 (7.68)

como puede comprobarse inmediatamente usando conservando del cuadrimomento total

p1 + p2 = p3 + p4.
En el sistema CM tendremos p ⃗1 = −⃗ ⃗3 = −⃗
p2 y p p 4 con E 1 + E 2 = E 3 + E 4 = EC M y por
p
tanto se puede ver inmediatamente que s = EC M es la energía total en el sistema centro
de masas.
Calculemos ahora el LIPS (7.59) d Φ(2) en este caso. Puesto que en el sistema CM p ⃗1 +
⃗2 = 0 tendremos
p 5

d 3p3 d 3p4
d Φ(2) = (2π)4 δ(3) −⃗ ⃗4 δ(EC M − E 3 − E 4 )
¡ ¢
p3 − p (7.69)
(2π)3 2E 3 (2π)3 2E 4
q
e integremos sobre los momentos finales. Tenemos por tanto, E 3 = p ⃗32 + m 32 y la integral
q
⃗4 impone E 4 = p
sobre p ⃗32 + m 42 de forma que

d 3p3
Z
d Φ(2) = (2π)δ(EC M − E 3 − E 4 ) (7.70)
(2π)3 (2E 3 )(2E 4 )

p 3 | tendremos
separando ahora la parte angular del módulo |⃗

d |⃗ p 3 |2
p 3 ||⃗
Z
(2)
dΦ = dΩ (2π)δ(EC M − E 3 − E 4 ) (7.71)
(2π)3 (2E 3 )(2E 4 )

ahora haciendo el cambio de variable x = EC M − E 3 − E 4 tenemos que

¯ ¯
¯ d x ¯ |⃗
¯ = p 3 | + |⃗
p3| EC M
µ ¶
¯
¯ d |⃗ = p
|⃗ 3 | (7.72)
p3| ¯ E3 E4 E3E4
5 ⃗3 + p
Nótese que la delta es la que impone p ⃗4 = 0

Antonio L. Maroto, Ignazio Scimemi 104

Campos Cuánticos 2022-23 UCM

y haciendo la integral en x llegamos a

p3|
1 |⃗
d Φ(2) = 2
dΩ (7.73)
16π EC M
Por otra parte, en el sistema CM tenemos que

⃗2 ¯¯
⃗1 p
¯ ¯
¯p EC M
v1 − ⃗
|⃗ v 2 | = ¯ − ¯ = |⃗
¯ p1| (7.74)
E1 E2 E1E2
de forma que podemos escribir finalmente para la sección eficaz diferencial 7.63, en el
sistema CM, integrando sobre estados finales

1 p3|
2 |⃗
d σC M = |M f i | dΩ (7.75)
64π2 s p1|
|⃗

donde hemos usado s = EC2 M y |⃗ p 1 | y |⃗

p 3 | se escriben como función de s usando E 1 + E 2 =
p
E 3 + E 4 = s. Si la colisión es elástica con m 1 = m 3 y m 2 = m 4 , podemos además usar que
p 3 | = |⃗
|⃗ p 1 |.

7.4.4 Ejemplo: Anchura de desintegración de mesones

Consideremos de nuevo la teoría escalar de Yukawa y calculemos el ritmo de desinte-
gración de mesones en el sistema CM.

φ → ϕ ϕ̄

Tenemos dos partículas en el estado final por lo que

|M f i |2 |M f i |2 p3|
1 |⃗
dΓ = d Φ(2) = 2
dΩ (7.76)
2E p 2E p 16π EC M

La amplitud de desintegración (7.24) corresponde a una amplitud invariante

i M f i = −i g (7.77)

y en el sistema CM, el mesón estaría en reposo y por tanto p µ = (m, 0) con E p = EC M = m.

µ µ
⃗3 ) y el antinucleón p 4 = (m/2, −⃗
Por otra parte, el nucleón saliente tendría p 3 = (m/2, p p 3 ).
µ 2
Por lo que, usando p 3 p 3µ = M llegamos a
s
m2
p3| =
|⃗ − M2
4

Antonio L. Maroto, Ignazio Scimemi 105

Campos Cuánticos 2022-23 UCM

de donde
s
g2 M2
dΓ = 1−4 dΩ (7.78)
64π2 m m2

Por tanto, integrando en el ángulo sólido llegamos al ritmo de desintegración total o an-
chura total
s
2
g M2
Γ= 1−4 2 (7.79)
16πm m

El tiempo de vida media correspondiente es

¶−1/2
M2
µ
16πm
τ= 1 − 4 (7.80)
g2 m2

Antonio L. Maroto, Ignazio Scimemi 106

8. Electrodinámica cuántica

La electrodinámica cuántica (QED) describe la interacción entre electrones (o cualquier

otra partícula cargada de spin 1/2 como los muones, quarks, etc) y los fotones, así como las
interacciones de las partículas cargadas entre sí. La teoría fue desarrollada por F. Dyson,
R. Feynman, J. Schwinger y S. Tomonaga a finales de los años 40 y es la primera teoría
cuántica de campos que permitió hacer cálculos de forma consistente. La constante de
acoplo de la teoría viene dada por la carga del electrón e, o más concretamente, por la
constante de estructura fina α = e 2 /4πℏc ≃ 1/137 que es una cantidad pequeña y por tanto
la teoría permite un desarrollo perturbativo como el que hemos visto en el tema anterior.
En este tema estudiaremos las reglas de Feynman para QED y los principales procesos
de interacción.

8.1 Reglas de Feynman para QED

La densidad Lagrangiana que describe la interacción entre fermiones de Dirac con masa
m y carga e con el campo electromagnético es

1 1
LQE D = ψ̄(i ∂/ − m)ψ − F µν F µν − (∂µ A µ )2 − e ψ̄γµ ψA µ (8.1)
4 2
donde usaremos la cuantización covariante del campo electromagnético. La densidad La-
grangiana de interacción está dado por L i nt = −e ψ̄γµ ψA µ y por tanto
Z
Hi nt = e d 3 x ψ̄γµ ψA µ (8.2)

Vemos por tanto que el único tipo de vértice que aparece en la teoría es un vértice de tres
partículas que involucra dos campos fermiónicos y un fotón.
Cuando tenemos varios tipos de fermiones interaccionando electromagnéticamante,

107
Campos Cuánticos 2022-23 UCM

la densidad Lagrangiana será

1 1
LQE D = ψ̄a (i ∂/ − m a )ψa − F µν F µν − (∂µ A µ )2 − |e| Q a ψ̄a γµ ψa A µ
X X
(8.3)
a 4 2 a

donde m a es la masa y Q a es la carga del fermión correspondiente en unidades de la carga

del electrón. Es decir, Q a = −1 para electrones, muones y taus, mientras que Q a = +2/3
para quarks tipo "up" y Q a = −1/3 para quarks tipo "down".
Para derivar las reglas de Feynman de QED, consideraremos el proceso más sencillo, la
aniquilación de un electrón y un positrón en un par muón-antimuón. La carga del muón
es la misma que la del electrón y su masa es m µ ≃ 200m e .

8.1.1 Ejemplo: e + e − → µ+ µ−
Tenemos en este caso dos familias de partículas cargadas con la misma carga por tanto
Z
Hi nt = e d 3 x ψ̄a γµ ψa A µ
X
(8.4)
a=e,m

donde ψe y ψm denotan los campos del electrón y del muón respectivamente.

Ahora tenemos que el Hamiltoniano de interacción contiene dos tipos de vértices: uno
de electrones con fotones y otros de muones con fotones, con la misma constante de
acoplo e. El desarrollo del campo de electrones en operdores de creación y destrucción
irá como ψe ∼ au + b † v, ψ̄e ∼ b v̄ + a † ū, mientras que el de los muones ψm ∼ cu + d † v,
ψ̄m ∼ d v̄ + c † ū. Por otra parte para los fotones tendremos A µ ∼ f ϵµ + f † ϵ∗µ . Nótese que
hemos denotado de forma distinta los operadores de creación y destrucción de cada tipo
de partícula para evitar confundirlas.
El estado inicial del sistema contiene un electrón con momento p ⃗1 y estado de spin r y
un positrón con p ⃗2 y estado de spin s. El final, un muón con momento p⃗1′ y estado de spin
′
r y un antimuón con p ⃗2 y s , es decir
′ ′

q
|i 〉 = 4E p 1 E p 2 a p† 1 ,r b p† 2 ,s |0〉 (8.5)
q
|f 〉 = 4E p ′ E p ′ c p† ′ ,r ′ d p† ′ ,s ′ |0〉 (8.6)
1 2 1 2

Para poder conectar estos dos estados necesitamos de nuevo ir a segundo orden en el de-
sarrollo de Dyson (necesitamos dos aniquiladores de electrones y dos creadores de muones),

Antonio L. Maroto, Ignazio Scimemi 108

Campos Cuánticos 2022-23 UCM

por tanto:

(−i e)2
Z
d 4 x1 d 4 x2 T { ψ̄a (x 1 )γµ ψa (x 1 )A µ (x 1 ) ψ̄b (x 2 )γν ψb (x 2 )A ν (x 2 )}|i 〉
X X
iTf i = 〈f |
2 a b
(8.7)

Utilizando el Teorema de Wick, el único término que proporciona un resultado no nulo

es el que tiene los fotones contraidos (puesto que los estados inicial y final no contienen
fotones, cualquier término que contenga campos fotones dará cero al evaluarlo entre los
estados)

ψ̄a (x 1 )γµ ψa (x 1 ) ψ̄b (x 2 )γν ψb (x 2 ) : A µ (x 1 )A ν (x 2 )

X X
: (8.8)
a=e,m b=e,m

Ahora bien, de estos términos, sólo el que contenga dos campos electrónicos y dos muóni-
cos dará una contribución no nula y puesto que los campos electrónicos y muónicos anti-
conmutan, tendremos que los dos posibles términos son idénticos:
Z
i T f i = (−i e)2 d 4 x 1 d 4 x 2 〈 f | : ψ̄m (x 1 )γµ ψm (x 1 )ψ̄e (x 2 )γν ψe (x 2 ) : |i 〉A µ (x 1 )A ν (x 2 ) (8.9)

Ahora bien

〈 f | : ψ̄m (x 1 )γµ ψm (x 1 )ψ̄e (x 2 )γν ψe (x 2 ) : |i 〉 = 〈 f | : ψ̄m (x 1 )γµ ψm (x 1 ) : |0〉〈0| : ψ̄e (x 2 )γν ψe (x 2 ) : |i 〉

(8.10)

donde

〈0| : ψ̄e (x 2 )γν ψe (x 2 ) : |i 〉

d 3 k1 d 3 k2 q
Z
〈0|b k2 ,s ′′ a k1 ,r ′′ a p† 1 ,r b p† 2 ,s |0〉v̄(k 2 , s ′′ )γν u(k 1 , r ′′ )e −i (k1 +k2 )x2
X
= 4E p 1 E p 2
6
p
(2π) 4E k1 E k2 r ′′ ,s ′′

= v̄(p 2 , s)γν u(p 1 , r )e −i (p 1 +p 2 )x2 (8.11)

y de forma análoga para el término de los muones

〈 f | : ψ̄m (x 1 )γµ ψm (x 1 ) : |0〉

d 3 k1 d 3 k2 q
Z
〈0|d p ′ ,s ′ c p ′ ,r ′ c k† ,r ′′ d k† ,s ′′ |0〉ū(k 1 , r ′′ )γµ v(k 2 , s ′′ )e i (k1 +k2 )x1
X
= 4E p ′ E p ′
6
p 1 2 2 1
(2π) 4E k1 E k2 ′′
r ,s ′′
1 2

= ū(p 1′ , r ′ )γµ v(p 2′ , s ′ )e i (p 1 +p 2 )x1

′ ′
(8.12)

Antonio L. Maroto, Ignazio Scimemi 109

Campos Cuánticos 2022-23 UCM

Figure 8.1: Diagrama de Feynman correspondientes a la aniquilación electrón-positrón en

muón-antimuón.

y podemos escribir
Z
i T f i = (−i e)2 d 4 x 1 d 4 x 2 ū(p 1′ , r ′ )γµ v(p 2′ , s ′ )e i (p 1 +p 2 )x1 v̄(p 2 , s)γν u(p 1 , r )e −i (p 1 +p 2 )x2
′ ′

d 4 p −i η µν −i p(x1 −x2 )
Z
× e
(2π)4 p 2 + i ϵ
−i η µν
= ū(p 1′ , r ′ )(−i eγµ )v(p 2′ , s ′ ) v̄(p 2 , s)(−i eγν )u(p 1 , r )(2π)4 δ(4) (p 1 + p 2 − p 1′ − p 2′ )
(p 1 + p 2 )2 + i ϵ
(8.13)

Es decir, el elemento de matriz invariante es

−i η µν
i M f i = ū(p 1′ , r ′ )(−i eγµ )v(p 2′ , s ′ ) v̄(p 2 , s)(−i eγν )u(p 1 , r ) (8.14)
(p 1 + p 2 )2 + i ϵ

De forma totalmente análoga al tema anterior, podemos representar gráficamente este

proceso con un único diagrama de Feynman (ver Fig. 8.1) y asociando las siguientes reglas
de Feynman de QED. Esta misma amplitud podría emplearse para la aniquilación en pares
quark-antiquark simplemente cambiando la carga que aparece en el acoplo.

Antonio L. Maroto, Ignazio Scimemi 110

Campos Cuánticos 2022-23 UCM

Reglas de Feynman de QED

En general, para otros fermiones tendremos que reemplazar e → Q|e| en el vértice

Antonio L. Maroto, Ignazio Scimemi 111

Campos Cuánticos 2022-23 UCM

como hemos comentado antes. Cuando tenemos que sumar varios diagramas es impor-
tante determinar el signo relativo entre ellos. Para ello debemos ver si la ordenación de
los espinores en las amplitudes es la misma o no. Así por ejemplo, dos amplitudes con es-
pinores ūu v̄ v y v̄u ūv tendrán signos opuestos puesto que requieren intercambiar el orden
de operadores de creación o destrucción fermiónicos.

8.2 Procesos elementales

Una vez obtenidas la reglas de Feynman podemos utilizarlas para calcular todos los proce-
sos que involucran la interacción electromagnética de fermiones cargados.
Consideremos en más detalle la aniquilación electrón-positrón estudiada anteriormente
y obtengamos la correspondiente sección eficaz

8.2.1 Sección eficaz de aniquilación e + e − → µ+ µ−

A partir de la amplitud invariante calculada anteriormente en (8.14) vamos a calcular la
sección eficaz de aniquilación en el sistema centro de masas usando (7.75). Necesitamos
por tanto evaluar |M |2 . Para ello necesitamos obtener los complejos conjugados de térmi-
nos de la forma v̄γµ u, que usando (γµ )† = γ0 γµ γ0 , se escribe

(v̄γµ u)∗ = u † (γµ )† (γ0 )† v = u † (γµ )† γ0 v = u † γ0 γµ v = ūγµ v (8.15)

Por tanto, podemos escribir (por simplicidad hemos quitado la etiqueta de spin)

e4 µ ν
|M |2 = ′ ′ ′ ′
¡ ¢¡ ¢
ū(p 1 )γ v(p 2 ) v̄(p 2 )γ u(p 1 ) v̄(p 2 )γµ u(p 1 )ū(p 1 )γν v(p 2 ) (8.16)
(p 1 + p 2 )4

En la mayor parte de los casos, los haces de electrones y positrones iniciales no están polar-
izados, por tanto la sección eficaz que medimos es un promedio sobre los estados de spin r
y s de los electrones y positrones iniciales. Por otra parte, los detectores de muones no son
sensibles al spin de los mismos, por tanto la sección eficaz es la suma sobre los estados de
spin finales r ′ y s ′ . Es decir, tenemos que calcular

1X1XX
|M (r, s → r ′ , s ′ )|2 (8.17)
2 r 2 s r ′ ,s ′

Antonio L. Maroto, Ignazio Scimemi 112

Campos Cuánticos 2022-23 UCM

Estas sumas como veremos simplifican enormemente el cálculo de la sección eficaz. En

efecto, podemos emplear para la suma en estados de spin las relaciones (5.59)
X
u(p, s)ū(p, s) = p/ + m (8.18)
s=1,2
X
v(p, s)v̄(p, s) = p/ − m (8.19)
s=1,2

Así, cogiendo el segundo factor de (8.16) tenemos que escribiendo explícitamente los índices
espinoriales
X X
v̄(p 2 , s)γµ u(p 1 , r )ū(p 1 , r )γν v(p 2 , s) = v̄ a (p 2 , s)(γµ )ab u b (p 1 , r )ū c (p 1 , r )(γν )cd v d (p 2 , s)
r,s r,s
£ ¤
=(p/ 2 − m e )d a (γµ )ab (p/ 1 + m e )bc (γν )cd = tr (p/ 2 − m e )γµ (p/ 1 + m e )γν (8.20)

Procedemos de forma análoga con el primer factor y obtenemos

1 X e4
|M |2 = tr (p/ ′1 + m µ )γµ (p/ ′2 − m µ )γν tr (p/ 2 − m e )γµ (p/ 1 + m e )γν (8.21)
£ ¤ £ ¤
4 spi ns 4(p 1 + p 2 )4

Utilizando ahora las propiedades de las traza de las matrices gamma tenemos
h i
tr (p/ 2 − m e )γµ (p/ 1 + m e )γν = 4 p 1 µ p 2 ν + p 1 ν p 2 µ − η µν (p 1 · p 2 + m e2 )
£ ¤
(8.22)

y
′µ ′ν ′ν ′µ
h i
µ ′ ν µν ′ 2
£ ′ ′
tr (p/ 1 + m µ )γ (p/ 2 − m µ )γ = 4 p 1 p 2 + p 1 p 2 − η (p 1 · p 2 + m µ )
¤
(8.23)

Multiplicando los dos factores obtenemos directamente

1 X 8e 4
|M |2 = (p 1 · p 1′ )(p 2 · p 2′ ) + (p 1 · p 2′ )(p 2 · p 1′ ) − (p 1′ · p 2′ )(p 1 · p 2 + m e2 )
£
4 spi ns (p 1 + p 2 ) 4
i
−(p 1 · p 2 )(p 1′ · p 2′ + m µ2 ) + 2(p 1 · p 2 + m e2 )(p 1′ · p 2′ + m µ2 ) (8.24)

Podemos ahora considerar la siguiente aproximación. Puesto que la masa del muón es
mucho mayor que la masa del electrón, la producción de muones requiere electrones ul-
trarrelativistas, por lo que podemos despreciar la masa del electrón y suponer m e = 0. En
ese caso la expresión anterior se simplifica de forma que

1 X 8e 4 h i
|M |2 = (p 1 · p ′
1 )(p 2 · p ′
2 ) + (p 1 · p ′
2 )(p 2 · p ′
1 ) + (p 1 · p 2 )m 2
µ (8.25)
4 spi ns (p 1 + p 2 )4

Antonio L. Maroto, Ignazio Scimemi 113

Campos Cuánticos 2022-23 UCM

Figure 8.2: Colisión en el sistema centro de masas con p̂ 1 · p̂ 1′ = p̂ 2 · p̂ 2′ = cos θ

En el sistema centro de masa (ver Fig.8.2), tenemos que en el estado inicial p ⃗1 = −⃗ p2 y

por tanto como ambas partículas tienen la misma masa E 1 = E 2 . De la misma forma en el
p
estado final p⃗1′ = −⃗
p 2′ y E 1′ = E 2′ . La energía total en el centro de masas será s = E 1 + E 2 =
2E 1 y por conservación de la energía E 1 = E 1′ . Por tanto

p 1 · p 1′ = p 2 · p 2′ = E 12 − E 1 |⃗
p 1′ | cos θ

donde hemos usado que m e = 0. Por otra parte,

p 1 · p 2′ = p 2 · p 1′ = E 12 + E 1 |⃗
p 1′ | cos θ

y finalmente
p 1 · p 2 = 2E 12
El momento transferido que aparece en el propagador será

(p 1 + p 2 )2 = 4E 12

con lo que podemos escribir finalmente

1 X 8e 4 h 2 i
|M |2 = (E 1 − E p
|⃗
1 1
′
| cos θ) 2
+ (E 2
1 + E p
|⃗
1 1
′
| cos θ)2
+ 2E 2 2
m
1 µ
4 spi ns 16E 14
2 2
"Ã ! Ã ! #
m µ m µ
= e 4 1 + 2 + 1 − 2 cos2 θ (8.26)
E1 E1

Antonio L. Maroto, Ignazio Scimemi 114

Campos Cuánticos 2022-23 UCM

Por tanto, la sección eficaz diferencial de aniquiliación en el sistema centro de masas será
2 2
"Ã ! Ã ! #
d σC M 1 |⃗ p 1′ | 1 X e 4 p
|⃗ ′
1 | m µ m µ
= |M f i |2 = 1 + 2 + 1 − 2 cos2 θ
dΩ 64π2 s |⃗
p 1 | 4 spi n 256π2 E 12 E 1 E1 E1
vÃ
2 m µ2 m µ2
! "Ã ! Ã ! #
α2 u m
u
µ
= t 1− 1 + 2 + 1 − 2 cos2 θ (8.27)
16E 12 E 12 E1 E1

que como vemos es de orden α2 . Integrando la sección eficaz diferencial en d Ω = d φd cos θ

obtenemos la sección eficaz total en función de la energía total en el sistema centro de
masas
vÃ
4m µ2 2m µ2
!Ã !
2u
u
4πα t
σC M = 1− 1+ (8.28)
3s s s

Vemos que la sección eficaz es nula para s < 4m µ2 puesto que en ese caso, el par electrón-
positrón no tiene energía suficiente para producir el par muón-anitmuón.

8.2.2 Dispersión Compton: e − γ → e − γ

Vamos a considerar ahora un tipo de proceso diferente, la dispersión electrón-fotón que
involucra la suma de los dos diagramas de Feynman en la Fig.8.3. Usando las reglas de
Feynman de QED, la amplitud invariante del proceso está dada por
i (p/ 1 + p/ 2 + m)
i M = ū(p 1′ , r ′ )(−i eγµ )ϵ∗µ (p 2′ , λ′ ) (−i eγν )ϵν (p 2 , λ)u(p 1 , r )
(p 1 + p 2 )2 − m 2
ν i (p/ 1 − p/′ 2 + m)
′ ′
+ ū(p 1 , r )(−i eγ )ϵν (p 2 , λ) ′ 2 2
(−i eγµ )ϵ∗µ (p 2′ , λ′ )u(p 1 , r )
(p 1 − p 2 ) − m
· µ
γ (p/ 1 + p/ 2 + m)γν γν (p/ 1 − p/′ 2 + m)γµ
¸
2 ′ ′ ∗ ′ ′
= −i e ū(p 1 , r )ϵµ (p 2 , λ )ϵν (p 2 , λ) + u(p 1 , r )
(p 1 + p 2 )2 − m 2 (p 1 − p 2′ )2 − m 2
(8.29)

Nótese que en las dos amplitudes el orden de los espinores es el mismo luego deben sumarse.
Esta expresión puede simplificarse, teniendo en cuenta que p 12 = m 2 y p 22 = p 2′2 = 0 con lo
que (p 1 + p 2 )2 − m 2 = 2p 1 · p 2 y (p 1 − p 2′ ) − m 2 = −2p 1 · p 2′ . Por otra parte, en el numerador
de los propagadores podemos usar que

(p/ 1 + m)γν u(p 1 ) = (2p 1ν − γν p/ 1 + γν m)u(p 1 ) = 2p 1ν u(p 1 ) − γν (p/ 1 − m)u(p 1 ) = 2p 1ν u(p 1 )

(8.30)

Antonio L. Maroto, Ignazio Scimemi 115

Campos Cuánticos 2022-23 UCM

Figure 8.3: Diagramas de Feynman correspondientes a la dispersión Compton electrón-

fotón.

Por tanto,
µ#
γµ p/ 2 γν + 2γµ p 1ν −γµ p/ ′2 γν + 2γν p 1
"
2
i M = −i e ū(p 1′ , r ′ )ϵ∗µ (p 2′ , λ′ )ϵν (p 2 , λ) + u(p 1 , r )
2p 1 · p 2 −2p 1 · p 2′
(8.31)

De forma totalmente análoga al ejemplo anterior, ahora podemos promediar sobre el spin
y la polarización de las partículas iniciales y sumar sobre las finales. Para los fermiones
usamos el mismo truco, mientras que para los fotones puede demostrarse que
X ∗
ϵµ (p, λ)ϵν (p, λ) → −η µν (8.32)
λ

Esta expresión no es una identidad estricta, pero se puede demostrar que es válida siempre
que aparezca en el interior de una amplitud de QED. Utilizando las expresiones anteriores
llegamos a
p 1 · p 2′ p 1 · p 2
· µ ¶ µ ¶2 ¸
1 X 2 4 2 1 1 4 1 1
|M | = 2e + + 2m − +m −
4 spi ns p 1 · p 2 p 1 · p 2′ p 1 · p 2 p 1 · p 2′ p 1 · p 2 p 1 · p 2′
(8.33)

Antonio L. Maroto, Ignazio Scimemi 116

Campos Cuánticos 2022-23 UCM

que permite calcular la sección eficaz del proceso. En el sistema laboratorio, la sección
eficaz diferencial viene dada por la llamada fórmula de Klein-Nishina
µ ¶2 µ ′
d σL AB α2 ω′ ω ω
¶
2
= + − sin θ (8.34)
dΩ 2m 2 ω ω ω′

donde ω y ω′ son las energías inicial y final del fotón relacionadas por la fórmula de Comp-
ton
1 1 1
− = (1 − cos θ)
ω ω m
′

con θ el ángulo entre la dirección del fotón incidente y el saliente.

Antonio L. Maroto, Ignazio Scimemi 117