Mena Miguel, Vintimilla Rafael, Peñafiel

James

Universidad San Francisco de Quito

Proyecto QFT
Sección Transversal y Tasas de Decaimiento

Instructor: Ernesto Contreras

December 19, 2023
Contents

1 Sección Transversal y Tasas de Decaimiento 2

1.1 Introducción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.2 Regla de Oro de Fermi . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 2
1.3 Matriz S . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 4
1.4 Tasa de Decaimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 5
1.5 Sección Transversal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 6

2 Ejemplo: Mesón Vectorial Producción y Decaimiento 15

2.1 Producción . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.1.1 e+ e− → µ+ µ− . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.1.2 Matriz S . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
2.1.3 Aproximación no relativista . . . . . . . . . . . . . . . . . . . . . . . . . . . 16
2.2 Mesón vectorial . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 17
2.3 Sección transversal . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 18
2.4 Tasa de Decaimiento . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 19
2.5 Variables Mandelstam . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 20

3 Más Ejemplos 21
3.1 Introduction . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
3.2 Dispersión de Compton . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 21
3.2.1 Suma sobre los estados de polarización . . . . . . . . . . . . . . . . . . . . . 22
3.2.2 Álgebra de trazas . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
3.2.3 Variables de Mandelstam . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 23
3.2.4 Cálculo del elemento de la matriz cuadrada completa . . . . . . . . . . . . . 24
3.3 Evaluación en distintos marcos de referencia . . . . . . . . . . . . . . . . . . . . . . 25
3.3.1 Marco de referencia del centro de masa . . . . . . . . . . . . . . . . . . . . . 25
3.3.2 Evaluación de la integral en el espacio de fases del marco de referencia del
centro de masa . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
3.3.3 Marco de Laboratorio . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
3.4 Gráficos . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 32

1
Chapter 1

Sección Transversal y Tasas de

Decaimiento

1.1 Introducción
En la interpretación de la teorı́a cuántica, las probabilidades de que ocurran ciertos eventos son el
cuadrado (módulo) de las amplitudes cuánticas. En este trabajo, nos enfocaremos en calcular estas
probabilidades, que son conocidas como tasas de desintegración y secciones transversales. Es lo que
medimos y estas mediciones son fundamentales para entender cómo las partı́culas interactúan y se
transforman bajo diferentes condiciones.
Una pequeña sutileza a considerar es que los elementos de la matriz S, denotados como ⟨f |S−1|i⟩,
incluyen un factor de (2π)4 δ (4) (pF − pI ). Este aspecto es crucial ya que terminamos con el cuadrado
de una función delta. Esta particularidad matemática tiene importantes implicaciones fı́sicas, ya
que refleja las condiciones de conservación de energı́a y momento en las interacciones de partı́culas.
Además, debemos tener en cuenta que este enfoque se basa en el trabajo en un espacio infinito,
un concepto fundamental en la fı́sica cuántica de campos. Trabajar en un espacio infinito nos
permite aplicar ciertas simplificaciones y asumir condiciones ideales, lo que facilita el cálculo de
estas probabilidades. Sin embargo, también conlleva ciertas complicaciones teóricas, como el manejo
del cuadrado de la función delta, que requiere un tratamiento cuidadoso y métodos matemáticos
avanzados.
En resumen, esta sección profundiza en la teorı́a cuántica de campos, explorando cómo se pueden
calcular las probabilidades de diversos procesos cuánticos mediante el uso de amplitudes cuánticas
y la matriz S. Este estudio no solo es fundamental para entender la teorı́a cuántica a un nivel más
profundo, sino también para aplicar estos conceptos en investigaciones prácticas y experimentales
en el campo de la fı́sica de partı́culas.

1.2 Regla de Oro de Fermi

La evolución temporal de un estado cuántico en presencia de una interacción está descrita pertur-
bativamente por la fórmula de Dyson. Para dos estados propios de energı́a |m⟩ y |n⟩, con Em ̸= En ,

2
tenemos, al primer orden en la interacción [1],
 Z t 
⟨m|U (t)|n⟩ = −i m dtHI (t) n
0
Z t
′
= −i ⟨m |Hint | n⟩ dt′ eiωt (1.1)
0
iωt
e −1
= − ⟨m |Hint | n⟩
ω
donde ω = Em − En , es la siguiente:
Esto nos da la probabilidad de la transición de |n⟩ a |m⟩ en el tiempo t, como

[sin2 (ωt/2)]
Pn→m (t) = |⟨m|U (t)|n⟩|2 = 4 |⟨m |Hint | n⟩|2 (1.2)
ω2
2
Usando eiθ − 1 = 4 sin2 (θ/2). Esta expresión matemática describe la probabilidad de tran-
sición entre dos estados cuánticos en un tiempo dado, en función de la diferencia de energı́a entre
esos estados y el tiempo transcurrido.

Figure 1.1: Gráfico de sin2 (ωt/2) /ω 2 versus ω para un valor fijo t; ω = (Em − En ) /2.[2]
 

La probabilidad de transición, como se muestra en la Figura 1.1, tiene un patrón de interferencia

y la probabilidad de transición de encontrar el sistema en un estado |m⟩ de energı́a Em es mayor
solo cuando En ≃ Em o cuando ω ≃ 0.
A medida que t → ∞, la función en la figura comienza a acercarse a una función delta.
Aprovechamos la siguiente identidad de la función delta

sin2 (ax)
δ(x) = lim (1.3)
a→∞ πax2

Para ver que a medida que t → ∞

2 sin2 (ωt/2)
→ πtδ(ω) (1.4)
ω2
En la Figura 1.1 hemos asumido que ω es una variable continua; es decir, hemos considerado
un continuo de estados finales. Si estos estados tienen una densidad ρ(E). En el lı́mite t → ∞,

3
tenemos la probabilidad de transición.
2 sin2 (ωt/2)
Z  
2
Pn→m = dEm ρ (Em ) 2 |⟨m |Hint | n⟩|
ω2 (1.5)
2
→ 2π |⟨m |Hint | n⟩| ρ (En ) t
Esto proporciona una probabilidad constante por unidad de tiempo para la transición de estados
alrededor de la misma energı́a En ∼ Em = E.
Ṗn→m = 2π |⟨m |Hint | n⟩|2 ρ(E) (1.6)
Esta es la Regla de Oro de Fermi.

Sin embargo, si primero calculamos la amplitud para el estado |n⟩ en t → −∞ para hacer la
transición al estado |m⟩ en t → +∞. Entonces obtenemos 1.1.
Z ∞
′
lim ⟨m|U (t, −t)|n⟩ = −i ⟨m |Hint | n⟩ dt′ eiωt
t→∞ −∞ (1.7)
= −i ⟨m |Hint | n⟩ δ(ω)
Ahora, al elevar al cuadrado la amplitud para obtener la probabilidad, nos encontramos con el
problema del cuadrado de la función delta:
Pn→m = |⟨m |Hint | n⟩|2 (2π)2 δ(ω)2 . (1.8)
Siguiendo los cálculos anteriores, nos damos cuenta de que la infinito adicional proviene porque
Pm→n es la probabilidad de que la transición ocurra en un tiempo infinito. Entonces hacemos la
correspondencia [1]
(2π)2 δ(ω)2 = (2π)δ(ω)T (1.9)
Donde T representa t → ∞. Al dividir por el intervalo de tiempo infinito T , estamos esencialmente
convirtiendo una probabilidad total a lo largo de un tiempo infinito en una tasa de probabilidad
por unidad de tiempo que nos permite renormalizar la teorı́a.
Ṗn→m = 2π |⟨m |Hint | n⟩|2 δ(ω) (1.10)
Que al integrar sobre la densidad de estados finales nos da nuevamente la Regla de Oro de Fermi
[1].

1.3 Matriz S
En teorı́a cuántica de campos podemos calcular los elementos de la matriz S para ir de un estado
|i⟩ al estado |f ⟩ utilizando la siguiente ecuación [3]
⟨f |S − 1|i⟩ = iMf i (2π)4 δ 4 (Pf − Pi ) (1.11)
Donde hemos descartado el modelo trivial, sin interacción S − 1. Además, los elemento de la
matriz siempre reflejan la conservación del cuadrimomentum (2π)4 δ 4 (Pf − Pi ). Extraemos este
factor para una expresión dependiente únicamente del elemento invariante de la matriz iMf i . Esta
es la amplitud de dispersión y puede calcularse utilizando los diagramas de Feynman. Para lo
expuesto en este proyecto tendremos en mente procesos donde dos partı́culas se dispersan, o una
partı́cula decae. Como es usual en teorı́a cuántica de campos, la probabilidad de que la partı́cula
sea encontrada en el estado final habiendo sido creada en el estado inicial es el cuadrado de la
amplitud cuántica.

4
1.4 Tasa de Decaimiento
Ahora veamos la probabilidad de que una sola partı́cula |i⟩ con momento
PpI (I=inicial) decaiga en
cierto número de partı́culas |f ⟩ con momento pi y momento total pF = i pi . Esto se da por [1]

|⟨f |S − 1|i⟩|2
Pi→f = (1.12)
⟨f | f ⟩⟨i | i⟩

De la fórmula de normalización relativista ⟨p | q⟩ = (2π)3 2Ep⃗ δ (3) (⃗p − ⃗q),

⟨i | i⟩ = (2π)3 2E−
→ δ (3) (0) = 2E−
p I
→V
p I
(1.13)
donde hemos reemplazado δ (3) (0) por el volumen del espacio de 3 dimensiones. De manera
similar, Y
⟨f | f ⟩ = 2Ep⃗i V (1.14)
estados finales

Si colocamos nuestra partı́cula inicial en reposo, de modo que p⃗I = 0 y Ep⃗I = m, obtenemos la
probabilidad de desintegración

|Mf i |2 Y 1
Pi→f = (2π)8 δ (4)2 (pI − pF ) (1.15)
2mV estados finales
2Ep⃗i V

Como en la segunda derivación de la Regla de Oro de Fermi 1.9, intercambiamos una de las funciones
delta por el volumen del espacio-tiempo: (2π)4 δ (4) (pI − pF ) = V T

|Mf i |2 Y 1
Pi→f = (2π)4 δ (4) (pI − pF ) V T (1.16)
2mV estados finales
2Ep⃗i V

Ahora podemos dividir por T para obtener la función de transición por unidad de tiempo.

|Mf i |2 Y 1
Ṗi→f = (2π)4 δ (4) (pI − pF ) (1.17)
2m estados finales
2Ep⃗i V

Hemos llegado a la fórmula análoga a 1.10 en QFT pero todavı́a tenemos que preocuparnos por
sumar sobre todos los estados finales. Integramos sobre todos los momentos posibles de las partı́culas
finales:
V d3 pi
Z Y
Ṗi→f 3
(1.18)
estados finales
(2π)
Los factores de volumen espacial V en esta medida cancelan aquellos en 1.17, mientras que los
factores de 1/2Ep⃗i en 1.17 producen la medida invariante de Lorentz para las integrales en el
espacion de fase invariante de Lorentz o LIP S por sus siglas en inglés [3].
Y d3 pi 1
dLIP S = (2π)4 δ (4) (pF − pI ) (1.19)
estados finales
(2π)3 2E−
→pi

El último paso es sumar sobre todos los estados finales con diferentes números (y posiblemente
tipos) de partı́culas. Esto nos da nuestra expresión final para la probabilidad de desintegración por
unidad de tiempo,. Z
S X
Γ = Ṗi→f = |Mf i |2 dLIP S (1.20)
2m estados finales

5
y es igual al recı́proco de la vida media τ = 1/Γ.
Si hay algunas partı́culas idénticas en el estado final, las contaremos en exceso al integrar sobre todas
las configuraciones de momento. El factor de simetrı́a S se encarga de esto. Si hay ni partı́culas
idénticas del tipo i en el estado final, entonces
Y 1
S= (1.21)
i
ni !

En resumen, la generalización de esta fórmula cumple las siguientes reglas [4]:

- Defina Mf i como en 1.11.

- Construya la medida invariante de Lorentz dLIP S con:

- un factor de (2π)4 δ 4 (pF − pI ),

- un factor de
d3 pi
(2π)3 2E−→
pi

por cada partı́cula en el estado final,

- un factor de 1/2m la partı́cula inicial en reposo,

-un factor de simetrı́a S por las partı́culas idénticas en el estado final.

1.5 Sección Transversal

En este proyecto final, se procederá a obtener la expresión para la sección transversal de dispersión
diferencial utilizando la matriz de dispersión. Esta expresión es fundamental, ya que vincula la
matriz de dispersión con valores que pueden ser medidos en el ámbito fı́sico. El enfoque particular
de esta derivación se centrará en el proceso en el cual dos partı́culas se separan en n partı́culas. Este
proceso no es solo uno de los más frecuentes en la teorı́a cuántica de campos, sino que también es
representativo de lo que comúnmente ocurre en los experimentos de dispersión, donde tı́picamente
se observan dos haces de partı́culas interactuando entre sı́ como observamos en la figura 1.2.

ρB v ρA

lB lA

Figure 1.2: Movimiento de partı́culas con densidades ρB y ρA .

Consideremos un material muy fino, cuyas partı́culas son tan pequeñas que no se superponen
cuando se observan de frente. Alternativamente, imagina una nube esparcida de partı́culas, en lugar
de enfocarse en un objetivo especı́fico. Bajo la suposición de que los objetivos están ampliamente
dispersos, la probabilidad de dispersión de una partı́cula se determina por la fracción del área del
objetivo ocupada por las partı́culas. El área total ocupada por las partı́culas es igual al área de
una partı́cula individual multiplicada por su número total. Dicha área de una partı́cula individual
es conocida como la sección transversal de dispersión total, que representa el tamaño del área que
un proyectil necesita impactar para dispersar la partı́cula objetivo en cualquier dirección. Esencial-
mente, esto representa el tamaño de la partı́cula objetivo. No obstante, en algunos casos se requiere

6
un análisis más detallado, como cuando se utiliza un detector que solo abarca un ángulo sólido
especı́fico alrededor del evento de dispersión, haciendo que la sección transversal de dispersión total
sea menos relevante. En estos casos, resulta más útil determinar el área especı́fica de la partı́cula
objetivo que debe ser impactada para que la dispersión ocurra dentro del ángulo sólido cubierto por
el detector.
Multiplicando este valor por el número de partı́culas en el objetivo y dividiéndolo por el área total
del objetivo, obtenemos la probabilidad de que las partı́culas proyectiles se dispersen especı́ficamente
hacia el detector situado en ese ángulo sólido. Para cualquier ángulo sólido arbitrario, se puede
calcular una medida llamada sección transversal de dispersión diferencial. Integrando esta sección
transversal diferencial sobre el ángulo sólido en cuestión, obtenemos el área necesaria. Si se integra
sobre un ángulo sólido que abarca toda la esfera, resulta en la sección transversal de dispersión
total.
Es importante destacar que en estos experimentos, el objetivo puede ser una única partı́cula,
como otro haz de un acelerador, o puede tener dimensiones reducidas en comparación con la distancia
al detector, haciendo que las diferencias angulares entre las partı́culas del objetivo y el detector sean
mı́nimas. Por lo tanto, se puede asumir que la sección transversal para un ángulo sólido dado es la
misma para todas las partı́culas en el objetivo.
Estos conceptos de secciones transversales de dispersión total y diferencial adquieren capas
adicionales de complejidad que los hacen valiosos en experimentos reales de fı́sica nuclear y de
partı́culas. Primero, las partı́culas cuánticas no tienen áreas definidas claramente, aunque anteri-
ormente se las trató como si las tuvieran, como esferas sólidas. A pesar de no tener áreas bien
definidas, las partı́culas cuánticas tienden a dispersarse en ciertos ángulos sólidos de manera con-
fiable, una tendencia que varı́a según el tipo de partı́cula en el objetivo y en el haz incidente. Por
lo tanto, deben tener un área efectiva, lo que es clave en el concepto de secciones transversales de
dispersión total y diferencial en fı́sica de partı́culas.
Otra capa de complejidad a considerar es que las secciones transversales de dispersión pueden
aplicarse a tipos especı́ficos de dispersión, no solo a la dispersión en general. Cuando las partı́culas
colisionan, pueden transformarse en otros tipos de partı́culas y en diferentes cantidades, no solo
rebotar. Se pueden definir secciones transversales de dispersión no solo para la dispersión en un
ángulo sólido especı́fico, sino también para un tipo especı́fico de dispersión, y luego para que las
partı́culas resultantes caigan en un conjunto especı́fico de ángulos sólidos y rangos de magnitudes de
momento. Incluso sin variar los tipos de partı́culas objetivo y proyectil, el ”área” efectiva asociada
con la dispersión para producir un conjunto especı́fico de partı́culas resultantes cambia dependi-
endo de la reacción de partı́culas considerada. Además, las secciones transversales de dispersión
generalmente se escriben también para rangos especı́ficos de magnitudes de momento de salida. Las
secciones transversales de dispersión diferencial se definen para intervalos de momentos de salida
diferenciales arbitrarios. Por tanto, para obtener el área seccional transversal efectiva para un pro-
ceso de dispersión que produzca partı́culas de salida en un conjunto especı́fico de ángulos sólidos y
momentos, se requiere un análisis detallado.
Hasta este punto, hemos definido la sección transversal de dispersión como el área seccional
transversal efectiva necesaria para producir las probabilidades y tasas de dispersión observadas.
Esto puede ser reformulado de una manera quizás menos intuitiva, pero más conveniente para
cálculos.

7
 Definamos los siguientes términos.

dσ : sección cruzada diferencial de dispersión

σ : sección cruzada del total de dispersión
A : area del objetivo
v : velocidad del haz
ρB : densidad del haz
ρA : densidad del objetivo
ρB v : flujo incidente
NB : número de partı́culas en el pulso del haz = ρB vA
NA : número de partı́culas en el objetivo = ρA V
lA : longitud del haz de partı́cula de A
lA : longitud del haz de partı́cula de B
N : número de eventos
nB : número de partı́culas en el pulso del haz por unidad de área
b : parámetro de impacto
V : Volumen del objetivo

Además, tenemos lo siguiente:

d(dPf ←i )
: Tasa de transición diferencial por unidad de volumen (1.22)
dV dt
dPf ←i = probabilidad de transición (1.23)
La ecuación 1.23 se refiere a los diferenciales del volumen en el espacio de momento saliente so-
bre el que se basa. Estos diferenciales pueden expresarse en términos de coordenadas esféricas,
dividiéndolos en diferenciales de ángulo sólido y de magnitud de momento. Con estas definiciones
en mano, es posible formular una ecuación inicial para la sección transversal de dispersión. Lo más
sencillo es empezar con la sección transversal total y después elaborar la versión diferencial. Para
comenzar, considerando que la sección transversal total representa el área efectiva de la sección
transversal de la partı́cula para un proceso de dispersión especı́fico, se puede formular la siguiente
expresión:
NT σ
Probabilidad del evento = (1.24)
A
Por consiguiente, el número de eventos de dispersión es simplemente como sigue
NA σ
Número de eventos = NB (1.25)
A
Se resuelve para la sección transversal.
number of events
σ= ·A (1.26)
NB NA
Sustituyendo y reacomodando.

números de eventos números de eventos/V t

σ= ·A= (1.27)
ρB · vt · ρA · A · V (ρB · v)ρA

8
Sustituyendo una vez más.
Tasa de transición por unidad de volumen
σ= (1.28)
densidad del objetivo · flujo incidente
La sección transversal diferencial de dispersión serı́a.

1 d(dPf ←i )
dσ = · (1.29)
densidad del objetivo · flujo incidente dV dt
Con la fórmula base para la sección transversal diferencial de dispersión ya definida, ahora se puede
continuar con la deducción de la fórmula explı́cita para el escenario en el que dos partı́culas se
dispersan en el marco de la teorı́a cuántica de campos.
Para este fin primero es necesario establecer paquetes de onda que simbolizan las partı́culas
en su estado inicial. Estos paquetes de ondas deben evolucionar a lo largo de un perı́odo extenso
utilizando el operador de evolución temporal derivado de la teorı́a del campo interactivo. Después,
se debe hacer coincidir el estado final obtenido con paquetes de ondas que correspondan a un grupo
especı́fico de partı́culas en su estado final. Este procedimiento nos permite determinar la amplitud
de probabilidad para alcanzar ese estado final, la cual tiene una relación directa con la sección
transversal. Se descubre que si los paquetes de ondas son muy estrechos en el espectro de momento,
la amplitud de probabilidad dependerá exclusivamente de los momentos de dichos paquetes y no de
la forma especı́fica que tengan.
Un paquete de onda representando un estado deseado |ϕ⟩ puede ser representado de la siguiente
manera [5]
d3 k
Z
1
|ϕ⟩ = 3
√ ϕ(k) |k⟩ (1.30)
(2π) 2Ek
Donde ϕ(k) es la transformada de Fourier de la función de onda especial, y |k⟩√es un estado de
una partı́cula de momento k en la teorı́a de interacción. Y, donde el factor 2E convierte la
normalización relativista del estado de una partı́cula a la normalización convencional en el cual la
suma de todas las probabilidades sea 1.

d3 k
Z
⟨ϕ|ϕ⟩ = 1 si 3
|ϕ(k)|2 = 1. (1.31)
(2π)

La probabilidad que se quiere calcular es

Pi→j = |⟨ϕ1 ϕ2 · · · |ϕA ϕB ⟩|2 . (1.32)

Al establecer el estado cuántico compuesto por |ϕA ϕB ⟩ en un punto lejano en el pasado y después
se considera el caso en que sus respectivos paquetes de ondas se enfocan alrededor de valores de
momentos especı́ficos, se logra delinear un estado entrante para |pA pB ⟩in con momentos iniciales
bien definidos. Este estado puede ser interpretado como la suma de varios estados más simples,
cada uno representando diferentes momentos definidos, combinados de manera lineal.
Por otro lado, es crucial tomar en cuenta cómo el paquete de ondas ϕB se desplaza en relación
a ϕA en términos de su posición (véase figura 1.3). Aunque esta consideración podrı́a omitirse en
la expresión de ϕB (kB ), optamos por una práctica en la que las funciones de onda que describen el
momento están alineadas (es decir, el parámetro de impacto b es cero) y ajustamos la función de
onda ϕB con un término exponencial para reflejar su desplazamiento espacial. Debido a que ϕA y
ϕB se formulan de manera independiente en lugares distintos, la descripción matemática del estado
inicial se realiza a través de la integral 1.33 que combina ambas funciones de onda, considerando la
influencia de su separación espacial.

9
 B
b
A

Figure 1.3: Paquetes de ondas incidentes son uniformemente distribuidas en el parámetro de impacto
b

d3 kA d3 kB ϕA (kA )ϕB (kB )e−ib·kB

Z
|ϕA ϕB ⟩in = p |kA kB ⟩in . (1.33)
(2π)3 (2π)3 (2EA )(2EB )
Similarmente, tenemos lo mismo para el estado final de las partı́culas
!
Y Z d3 pf ϕf (pf )
out ⟨ϕ1 ϕ2 · · · | = 3
p ·out ⟨p1 p2 · · ·| (1.34)
f
(2π) 2E f

Para comenzar con el cálculo de la sección transversal diferencial de dispersión tomando en

cuenta lo mencionado anteriormente, partimos de lo siguiente:
X
4 4
out ⟨p1 p2 . . . |kA kB ⟩in = iMf i (2π) δ (kA + kB − ( pf )) (1.35)

Hay que tomar en cuenta que no está normalizado los estados. Para hallar la probabilidad de ir de
un estado |i⟩ al estado |f ⟩ se realiza lo siguiente:

|out ⟨p1 p2 . . . |kA kB ⟩in |2

Pi→j = (1.36)
out ⟨p1 p2 . . . |p1 p2 . . .⟩out ·in ⟨kA kB |kA kB ⟩in

Y luego,
|iMf i (2π)4 δ 4 (kA + kB − ( pf ))|2
P
Pi→j = (1.37)
out ⟨p1 p2 . . . |p1 p2 . . .⟩out ·in ⟨kA kB |kA kB ⟩in

Si separamos
|(2π)4 δ 4 (kA + kB − ( pf ))|2 |Mf i |2
P
Pi→j = (1.38)
out ⟨p1 p2 . . . |p1 p2 . . .⟩out ·in ⟨kA kB |kA kB ⟩in

Aquı́ el módulo al cuadrado de la amplitud produce una delta al cuadrado sin significado aparente.
Este problema técnico surge porque nuestra amplitud se mide entre estados de onda plana no nor-
malizables, los cuales se extiende a través de todo el espacio-tiempo. Esto implica que el proceso de
dispersión ocurre en todo lugar y en todo momento y no se puede utilizar directamente para calcular
probabilidades en estos procesos. Una solución es construir paquetes de ondas normalizados, los
cuales se encuentran bien separados tanto en el pasado remoto como en el futuro lejano. Optaremos
por una solución en donde nuestro sisema en una caja con un volumen V = L3 . También suponemos
que la interacción se limita a actuar solo durante un perı́odo de tiempo de orden de T. Lo impor-
tante es que las respuestas finales sean independientes de los valores de V y T , reproduciendo ası́
los resultados obtenidos con el método más sofisticado de paquetes de onda.
En un espacio-tiempo finito con paquedes de ondas planas la amplitud de transición de i a f
está dado por la ecuación 1.11. Ahora, para eliminar los delta de dirac al cuadrado tomando en
cuenta una caja de tamaño finito L y tiempo finito T se realiza lo siguiente
Z ∞
f (x) δ 2 (x) dx (1.39)
−∞

10
Podemos expresar la ecuación 1.39 de la siguiente manera
Z ∞
g(x) δ(x) dx (1.40)
−∞

Donde g(x) = f (x)δ(x). La integral anterior resulta en lo siguiente

Z ∞
g(0) = δ(0)f (0) = δ(0) f (x) δ(x) dx (1.41)
−∞

Y como Z ∞ Z ∞
2
f (x) δ (x) dx = f (x) δ(0)δ(x) dx (1.42)
−∞ −∞

De lo anterior podemos ver que

δ 2 (x) = δ(0)δ(x) (1.43)
Por lo tanto,
(2π)4 δ 4 (0) = V T (1.44)
Regresando a la amplitud de probabilidad habiendo obtenido este resultado
(2π)4 δ 4 (kA + kB − ( pf ))|Mf i |2 V T
P
Pi→j = (1.45)
out ⟨p1 p2 . . . |p1 p2 . . .⟩out ·in ⟨kA kB |kA kB ⟩in

Por otro lado, en la teorı́a libre se tiene lo siguiente

|p1 p2 ⟩in = 2Ea 2Eb a†a a†b |0⟩

p p
(1.46)

Y p p
in ⟨p1 p2 · · ·| = 2Ea 2Eb |0⟩ aa ab (1.47)
Entonces,
out ⟨p1 p2 · · · |p1 p2 ⟩out = 2Ea 2Eb ⟨0|aa ab a†a a†b |0⟩ (1.48)
Como los operadores con diferentes ı́ndices conmutan, podemos escribirlo de la siguiente forma

out ⟨p1 p2 |p1 p2 ⟩out = 2Ea 2Eb ⟨0|aa a†a ab a†b |0⟩ (1.49)

Podemos utilizar la relación de conmutación de los operadores ab y a†b el cual es

ab a†b = a†b ab + (2π)3 δ(0) (1.50)

Por lo tanto, y si aplicamos esta relación de conmutación al vacı́o obtenemos

out ⟨p1 p2 |p1 p2 ⟩out = 2Ea 2Eb (2π)3 δ(0)(2π)3 δ(0) ⟨0|0⟩ (1.51)

Y tomando en cuenta la ecuación 1.44

out ⟨p1 p2 |p1 p2 ⟩out = 2Ea 2Eb V 2 (1.52)

Por lo tanto,
(2π)4 δ 4 (kA + kB − ( pf ))|Mf i |2 V T
P
Pi→j = (1.53)
2Ea 2Eb V 2 2Ec 2Ed V 2
Si queremos la probabilidad de transición por unidad de tiempo, entonces
(2π)4 δ 4 (kA + kB − ( pf ))|Mf i |2
P
Pi→j
= (1.54)
T 2Ea 2Eb V 2Ec 2Ed V · V

11
Aquı́ vale la pena aclarar que el componente matricial M guarda una analogı́a con la amplitud
de dispersión f utilizada en la mecánica cuántica de una sola partı́cula. Su utilidad radica en que
facilita la distinción entre los aspectos de la fı́sica que son influenciados por los detalles especı́ficos
del hamiltoniano de interacción, o sea, la dinámica, de aquellos que no están involucrados, conocidos
como cinemática.Es preciso entender cómo se puede reconstruir la sección transversal σ a partir de
M. Para lograr esto, debemos estimar la probabildiad de que el estado inicial |ϕA ϕB ⟩ se disperse y
evolucione hacia un estado final compuesto por n partı́culas, cuyos momentos están contenidos den-
tro de un volumen diferencial especificado por d3 p1 . . . d3 pn , siguiendo la normalización establecida
anteriormente. Entonces,
Y d3 pf 1
2
dPi→j = |out ⟨p1 p2 . . . |kA kB ⟩in | · (1.55)
f
(2π)3 2Ef

Por lo tanto,
(2π)4 δ 4 (kA + kB − ( pf ))|Mf i |2 Y d3 pf 1
P
dPi→j
= · (1.56)
dt 2Ea 2Eb · V f
(2π)3 2Ef

Ahora, dividimos para el flujo incidente.

número de partı́culas
ρB v = flujo incidente = (1.57)
área · tiempo
Aquı́, haremos el siguiente truco.
número de partı́culas l
ρB v = flujo incidente = · (1.58)
área · tiempo l

v
ρB v = flujo incidente = número de partı́culas · (1.59)
V
Donde l es una longitud que corresponde a los haces y, recordando que l/t es una velocidad.
Asimismo, el área por la longitud serı́a un volumen [4]. Entonces,

(2π)4 δ 4 (kA + kB − ( pf ))|Mf i |2 Y d3 pf 1

P
dPi→j
= · (1.60)
dt · flujo incidente 2Ea 2Eb · número de partı́culas · v f (2π)3 2Ef

Finalmente,

(2π)4 δ 4 (kA + kB − ( pf ))|Mf i |2 Y d3 pf 1

P
dPi→j · número de partı́culas
= · (1.61)
dt · flujo incidente 2Ea 2Eb · v f
(2π)3 2Ef

Note que

1 d(dPf ←i ) dPi→j · número de partı́culas

dσ = · = (1.62)
densidad del objetivo · flujo incidente dV dt dt · flujo incidente
Por lo tanto, hemos obtenido la expresión de la sección transversal diferencial de dispersión.
En este caso, v es la velocidad que tiene una de las partı́culas si estamos trabajando en el sistema
de referencia del laboratorio. En el caso de que estemos en el centro de momentos se tendrı́a una
velocidad relativa de las partı́culas.

(2π)4 δ 4 (kA + kB − ( pf ))|Mf i |2 Y d3 pf 1

P
dσ = · (1.63)
2Ea 2Eb · |vA − vB | f
(2π)3 2Ef

12
Este es un resultado obtenido a partir de la definición de la sección transversal y de manera intuitiva.
Podemos llegar al mismo resultado si consideramos que el número de partı́culas dispersadas para
un solo objetivo NA = 1 y para nB partı́culas viniendo por unidad de área transversal es
X Z
N= dPi = d2 b nB dPi→j (b) (1.64)
i
all incident
particles

Si nB es constante, y tomando en cuenta la ecuación 1.25, entonces

Z
N
σ= = d2 b dPi→j (b̄) (1.65)
nB
Aquı́, reemplazamos la forma de dPi→j de la ecuación 1.55 y, también, la amplitud de los estados
de las ecuaciones 1.33 y 1.34.
!Z !
Y d3 pf Y d3 ki ϕ(ki )
dσ = 3
d2 b
f
(2π) 2Ef i=A,B
(2π)3 2Ei
 3 ′ ∗ ′ 
d ki ϕ (ki ) ib·(kb −kb′ ) (1.66)
× e (out ⟨p p
1 2 . . . |k k ⟩
A B in )
(2π)3 2Ei′
× (out ⟨p1 p2 . . . |kA kB ⟩in )∗ .

Recordemos que Z
′
d2 b eib·(kb −kb ) = (2π)2 δ (2) (kb − kb′ ) (1.67)

Si tomamos en cuenta la ecuación 1.35 y su conjugado, podemos observar que al usar estas deltas
junto con la ecuación 1.67, desarrollamos seis de las integrales para k ′ de la ecuación 1.66.
Z Z X X 
3 ′ ′ (4)
d kA d3 kB δ ki − pf δ (2) (kb − kb′ ) (1.68)

Donde (kb⊥ = (kb′ )⊥ ). Ası́,

Z
′ ′ ′ ′
X X
= ki⊥ × dkAz dkBz δ(kAz + kBz − pzf )δ(EA′ + EB′ − Ef ) (1.69)

Utilizando otro de los delta para que nos quede una integral sobre KA
Z q q X 
′z ′2 ′
= dkA δ kA + m2A + kB2 + m2B − Ef (1.70)
′ ′
kBz = pzf −kAz
P

Finalmente, obtenemos
′ ′
kz kz
= A − B = |vA − vB | (1.71)
EA EB
√ ′ ′
Cabe aclarar que hemos utilizado el hecho de que como E = k 2 + m2 , obtenemos kBz = pf −kAz .
P

dE ′ A k′A dE ′ B k′B dE ′ B dE ′ B
= , = , = − (1.72)
dk ′ A E ′A dk ′ B E ′B dk ′ A dk ′ A
También, hemos utilizado el siguiente hecho.
Z
1
dx δ(f (x) − f (x0 )) = (1.73)
|f ′ (x 0 )|

13
Al final, obtenemos la sección transversal diferencial de dispersión [6].
!Z
Y d3 pf d3 kA
Z 3
d kB
dσ = 3 3
f
(2π) 2Ef (2π) (2π)3
(1.74)
2
|M| X
× |ϕA (kA )|2 |ϕB (kB )|2 (2π)4 δ (4) (kA + kB − pf )
2EA 2EB |vA − vB |

Para estados con funciones de onda centrados en un momento dado δ(⃗k − p⃗) obtenemos
1
dσ =
2EA 2EB |vA − vB |
Y d3 pf
!
X (1.75)
2 4 (4)
× |M| (2π) δ (kA + k B − pf ).
f
(2π)3 2Ef

Donde hemos obtenido la misma expresión que la ecuación 1.63, siendo esta la forma final de la
relación entre la matriz S y la sección transversal.

En resumen, la generalización de esta fórmula cumple las siguientes reglas [4]:

- un factor de (2π)4 δ 4 (kA + kB −

P
pf ) /v, donde v es el flujo,

- un factor de
d3 pf
(2π)3 2Ef
por cada partı́cula en el estado final,

- un factor de
1
4EA EB ,
donde EA y EB son las energı́as de las dos partı́culas en el estado inicial.

14
Chapter 2

Ejemplo: Mesón Vectorial Producción y

Decaimiento

2.1 Producción
Hasta ahora hemos considerado que los estados iniciales y finales de los procesos de dispersión sean
estados de partı́culas aisladas individuales. Sin embargo, muy cerca del umbral, la atracción de
Coulomb de los muones deberı́a convertirse en un efecto importante. Justo debajo del umbral,
todavı́a podemos formar pares µ+ µ− en estados ligados electromagnéticos.

2.1.1 e+ e− → µ+ µ−
La reacción e+ e− → µ+ µ− es el proceso más simple de todos los de QED (Electrodinámica
Cuántica), pero también uno de los más importantes en la fı́sica de altas energı́as. Es funda-
mental para entender todas las reacciones en colisionadores de e+ e− y, de hecho, se utiliza para
calibrar dichas máquinas [5].

2.1.2 Matriz S
El diagrama que describe el proceso QED es:

Figure 2.1: Diagrama de Feynman del proceso e+ e− → µ+ µ− .[5]

Empleando las siguientes reglas de feynman para procesos QED:

• Vértice: −ieγ µ .
−igµν
• Propagador fotónico q 2 +iϵ
.

15
 • Lı́neas externas
′
- v̄ s (p′ ) representa un positrón entrante.
- us (p) representa un electrón entrante.
- ūr (k) representa un muón saliente.
′
- v r (k ′ ) representa un anti-muón saliente.

Aplicando las reglas tenemos para la amplitud

 
− ′ − ′ s′ ′ −igµν ′
+ + µ s
ūr (k) (−ieγ ν ) v r (k ′ ) .


iM e (p)e (p ) → µ (k)µ (k ) = v̄ (p ) (−ieγ ) u (p) (2.1)
q2
ie2
= (v̄ (p′ ) γ µ u(p)) (ū(k)γµ v (k ′ )) (2.2)
q2
Dejando implı́citos los ı́ndices de spin. Sin embargo, en este régimen la polarización del estado
ligado, sensible a efectos relativistas, se dificulta mucho.

2.1.3 Aproximación no relativista

Afortunadamente, muchos de los sistemas ligados familiares en la naturaleza pueden ser tratados
(al menos como una buena primera aproximación) como sistemas no relativistas, en los cuales los
movimientos internos son lentos. Sin embargo, el proceso de crear las partı́culas constituyentes del
vacı́o sigue siendo un efecto relativista, requiriendo teorı́a cuántica de campos para su descripción
adecuada. El electrón y el positrón siguen siendo relativistas, por lo que simplificamos la expresión
escogiendo una helicidad definida. Tomamos el electrón diestro, moviéndose en la dirección +z, y
el positrón zurdo, moviéndose en la dirección −z Figura 2.2.

 
1
Figure 2.2: Para el electrón tenemos ξ = , correspondiente a espı́n hacia arriba en la dirección
0  
0
z, mientras que para el positrón tenemos ξ = , también espı́n hacia arriba en la dirección z.
1

Los espinores son:    

0 0
√  0  √  0 
u(p) = 2E 
 1 ;
 v (p′ ) = 2E 
 0 
 (2.3)
0 −1
Entonces, de la Ec. 2.2 tenemos:
 
′ µ µ 1
v̄ (p ) γ u(p) = 2E(0, −1)σ (2.4)
0
= −2E(0, 1, i, 0). (2.5)
En la otra mitad del elemento de la matriz debemos usar las expresiones no relativistas:
√ √
   ′ 
ξ ′ ξ
u(k) = m , v (k ) = m . (2.6)
ξ −ξ ′

16
 Entonces, de la Ec. 2.2 tenemos:
ξ′
  
µ ′ † †

σ̄ µ 0
ū(k)γ v (k ) = m ξ , ξ
0 σµ −ξ ′
( (2.7)
0 para µ = 0,
=
−2mξ † σ i ξ ′ para µ = i.

Para evaluar M, calculamos el producto punto entre 2.5 y 2.7 y multiplicamos por e2 /q 2 =
e2 /4m2 :
e2
M e− + + −
4m2 ξ † σ 1 + iσ 2 ξ ′

 
e
R L → µ µ = − 2
(2.8)
4m    
2 † 0 1 0 −i
= −e ξ +i ξ′ (2.9)
1 0 i 0
 
2 † 0 1
= −2e ξ ξ′. (2.10)
0 0
Dado que no hay dependencia angular en esta expresión, es igualmente probable que los muones
salgan en cualquier dirección. Más precisamente, son emitidos en una onda s; su momento angular
orbital es cero. La conservación del momento angular por lo tanto requiere que el espı́n total del
estado final sea 1, y de hecho el producto de matrices da cero a menos que tanto el muón como el
antimuón tengan el espı́n hacia arriba a lo largo del eje z (ver Fig. 2.3).

Figure 2.3: En el lı́mite relativista podemos asegurar que el espı́n total del sistema se conserva

2.2 Mesón vectorial

Consideremos primero el caso en que los espines del electrón y el positrón ambos apuntan hacia
arriba a lo largo del eje z. El único tipo de estado ligado que podemos producir tendrá un spin
total de 1. La amplitud de esta configuración acorde con la Ec. 2.10 es:

M (↑↑→ k1 ↑, k2 ↑) = −2e2 (2.11)

Las coordenadas centro de masa R y relativas r son
1
R= (r1 + r2 ) , r = r1 − r2 . (2.12)
2
Tienen momento conjugado
1
K = k1 + k 2 , k= (k1 − k2 ) (2.13)
2
17
 Donde K = 0 en el comovil del sistema ligado. El estado ligado es simplemente una super-
posición lineal de estados libres de r o k definidos, ponderados por la función de onda ψ(r) solución
de la ecuación de Schrödinger para un potencial columbiano. Por conveniencia construimos esta
superposición en el espacio de momentos, utilizando la transformada de Fourier de ψ(r):
Z
ψ(k)
e = d3 xeik·r ψ(r) (2.14)

Una expresión explı́cita para un estado ligado con masa M ≈ 2m, momento K = 0 y espı́n 1
orientado hacia arriba es entonces:
√ Z
d3 k 1 1
|B⟩ = 2M ψ̃(k) √ √ |k ↑, −k ↑⟩. (2.15)
(2π)3 2m 2m
√
Donde los factores deberı́an involucrar 2E±k , pero para un estado ligado no relativista, |k| ≪ m.
Dada esta expresión para el estado ligado, podemos escribir inmediatamente la amplitud para su
producción:
√ Z
d3 k e∗ 1 1
M(↑↑→ B) = 2M ψ (k) √ √ M(↑↑→ k ↑, −k ↑) (2.16)
(2π)3 2m 2m
Dado que la amplitud del estado libre de 2.11 es independiente de los momentos de los muones:
r
d3 k e∗
Z
2 2


M(↑↑→ B) = −2e ψ (k) (2.17)
M (2π)3
r
2
= (−2e2 )ψ ∗ (0) (2.18)
M
Esta fórmula muestra que la amplitud para la creación de un mesón vectorial es proporcional a
la carga al cuadrado e2 , y depende del valor de la función de onda del estado ligado ψ en el origen.

2.3 Sección transversal

La Ecuación 2.18 puede convertirse directamente en una sección transversal para la producción de
estados ligados µ+ µ− en la aniquilación e+ e− . Sin embargo, si relajamos la suposición de polar- √
ización del spin del estado ligado mediante vectores de polarización: Estado inicial ϵ+ = (x̂+iŷ)/ 2
de Ec. 2.5, a n un estado ligado spin-1 o spin-0. Entonces la Ec. 2.18 puede reescribirse de la Ec.
2.2 en una forma más general como
r
− + 2
−2e2 (n∗ · ϵ+ ) ψ ∗ (0)



M eR eL → B = (2.19)
M
Y si los electrones están inicialmente no polarizados, la sección transversal para la producción
de B implicará el promedio de polarización
1
1
|n∗ · ϵ+ |2 + |n∗ · ϵ− |2 = (nx )2 + (ny )2


(2.20)
4 4
√
Donde ϵ− = (x̂ − iŷ)/ 2 para una polarización de estados ligados a lo largo del eje de colisión
e+ e− . Todo lo que queda es introducir esta expresión en la fórmula de la sección transversal derivada
en la Sección 1.5. Dado que solo hay dos partı́culas en el estado final y estamos trabajando en el

18
marco de referencia del centro de masas. Para nuestro problema, |vA − vB | = 2 y EA = EB = m,
ası́ que tenemos de Ec. 1.63

(2π)4 δ 4 (p + p′ − K)|Mf i |2 d3 K 1
dσ e+ e− → B = · (2.21)
2m2m · 2 (2π)3 2EK

Asumiendo un haz de electrones no polarizado y sumando 2.20 sobre las tres posibles direcciones
de n.

(2π)4 δ 4 (p + p′ − K) d3 K 1 2 4 1
dσ e+ e− → B = |ψ(0)|2


· 3
· 4e (2.22)
2m2m · 2 (2π) 2EK M 2
Encontramos la siguiente expresión para la sección transversal total para la producción del estado
ligado:
Z 3
+ −

1 1 1 dK 1 2
1
σ e e →B = 3
(2π)4 δ (4) (p + p′ − K) · 4e4 |ψ(0)|2 . (2.23)
2 2m 2m (2π) 2EK M 2

La integral del espacio de fase de 1-cuerpo solo puede eliminar tres de las cuatro funciones delta.
Es convencional reescribir la última función delta usando δ (P 0 − K 0 ) = 2K 0 δ (P 2 − K 2 ). Entonces

|ψ(0)|2
σ e+ e− → B = 64π 3 α2 2 2



δ E cm − M . (2.24)
M3
La última función delta impone la restricción de que la energı́a total del centro de masas debe ser
igual a la masa del estado ligado; por lo tanto, el estado ligado se produce como una resonancia en
la aniquilación e+ e− . Si el estado ligado tiene una vida útil finita, esta función delta se ensanchará
en un pico de resonancia.

2.4 Tasa de Decaimiento

Si el estado ligado B puede ser producido a partir de e+ e− , también puede aniquilarse de vuelta
a e+ e− , o a cualquier otro par de leptones suficientemente ligeros. De acuerdo con 1.20, el ancho
total para este modo de desintegración se da por
Z
+ − 1
|M|2 dLIP S,


Γ B→e e = (2.25)
2M
donde S = 1 y M es simplemente el conjugado complejo del elemento de matriz de la sección
2.3 que usamos para calcular la producción de B. Como ahora el proceso se invierte
Z Z 3 Z 3 ′
dp 1 dp 1
dLIP S = (2π)4 δ (4) (K − p − p′ ) (2.26)
2π 2m 2π 2m

Como p⃗′ = −⃗p en el marco de referencia del centro de masa integramos en p′ las 3 deltas espaciales.
Pasamos a coordenadas polares

(P 0 )2 dP 0 dΩ 1
Z Z
dLIP S = δ(K 0 − P 0 ) (2.27)
(2π)2 4m2
(P 0 )2 dP 0
Z Z
dcos(θ)
= δ(M − 2m) (2.28)
2π (2m)2
 
1 dcos(θ)
= (2.29)
2π 2

19
Ası́
8e4
Z  
1 1 d cos θ
|ψ(0)|2 |n · ϵ|2 + |n · ϵ∗ |2 .


Γ= (2.30)
2M 2π 2 M
Ahora debemos sumar sobre los estados de polarización del electrón y promediar sobre los tres
posibles valores de n. Ası́ obtenemos

16πα2 |ψ(0)|2
Γ B → e+ e− = . (2.31)
3 M2

2.5 Variables Mandelstam

Vemos que en muchas de las amplitudes anteriores (en particular aquellas que incluyen el intercambio
de solo una partı́cula) las mismas combinaciones de momentos aparecen frecuentemente en los
denominadores. Hay nombres estándar para varias sumas y diferencias de momentos: se conocen
como variables de Mandelstam. Son
2
s = (p1 + p2 )2 = (p′1 + p′2 )
2 2
t = (p1 − p′1 ) = (p2 − p′2 )
2 2
u = (p1 − p′2 ) = (p2 − p′1 )

donde, como en los ejemplos anteriores, p1 y p2 son los momentos de las dos partı́culas iniciales, y
p′1 y p′2 son los momentos de las dos partı́culas finales. Podemos definir estas variables ya sea que las
partı́culas involucradas en la dispersión sean las mismas o diferentes. Para entender qué significan
estas variables, supongamos que las cuatro partı́culas son iguales. Nos situamos en el marco del
centro de masas, de modo que las dos partı́culas iniciales tienen cuatrimomentos

p1 = (E, 0, 0, p) y p2 = (E, 0, 0, −p) (2.32)

Luego las partı́culas se dispersan en algún ángulo θ y salen con momentos

p′1 = (E, 0, p sin θ, p cos θ) y p′2 = (E, 0, −p sin θ, −p cos θ) (2.33)

Entonces, de las definiciones anteriores, tenemos que

s = 4E 2 y t = −2p2 (1 − cos θ) y u = −2p2 (1 + cos θ) (2.34)

La variable s mide la energı́a total del centro de masas de la colisión, mientras que las variables
t y u son medidas del intercambio de momento entre partı́culas. (Básicamente son equivalentes,
solo que con las partı́culas salientes intercambiadas). Ahora, las amplitudes que involucran el
intercambio de una sola partı́cula pueden escribirse simplemente en términos de las variables de
Mandelstam. Veremos ejemplos de esto en la siguiente sección con más ejemplos.

20
Chapter 3

Más Ejemplos

3.1 Introduction
Hay muchos procesos fı́sicos que se pueden explicar con la matriz S y las secciones transversales,
vamos a centrarnos en uno de ellos en el cual usamos todos los conceptos estudiados en el capı́tulo
anterior.
En esta sección se va a calcular la sección transversal, o cross section de la dispersión de Compton,
se va a obtener las ecuaciones de Klein - Nishina, además de relacionar distintos valores de α con
distintas funciones de probabilidad en gráficos polares.

3.2 Dispersión de Compton

Uno de los ejemplos más conocidos de dispersión es la dispersión de Compton, en ella fotones de alta
energı́a interactúan con los electrones de las capas externas de un átomo y los arrancan del mismo,
o si son electrones libres los desvı́an de su curso natural y los fotones también se ven desplazados.
Para proceder con los cálculos necesitamos los diagramas de Feynman de esta interacción, mismos
que se muestran en la figura 3.1

Figure 3.1: Diagramas de Feynman asociados a la interacción de Compton [?]

A partir del diagrama, y usando las reglas de Feynman para interpretar los diagramas, se llega

21
a que la matriz S viene dada por la siguiente expresión:
i(p/ + k/ + m)
iS = ū(p′ )(−ieγ µ )ϵ∗µ (k ′ ) 2 2
(−ieγ ν )ϵν (k)u(p)
(p + k) − m
′
′ ′
i(p/ − k/ + m)
ν
+ū(p )(−ieγ )ϵν (k ) (−ieγ µ )ϵ∗µ (k)u(p)
(p − k ′ )2 − m2
′
" #
µ ν ν ν
γ ( p
/ + /
k + m)γ γ ( p
/ − /
k + m)γ
= −ie2 ϵ∗µ (k ′ )ϵν (k)ū(p′ ) + u(p)
(p + k)2 − m2 (p − k ′ )2 − m2

Esta última expresión puede simplificarse aún más si consideramos las condiciones fı́sicas de la
dispersi ’on de Compton: p2 = m2 y k 2 = 0 entonces los denominadores se ven reducidos a
(p + k)2 − m2 = p2 + 2p · k + k 2 − m2 = 2p · k y (p − k ′ )2 − m2 = p2 − 2p · k ′ + k ′2 − m2 = −2p · k ′ .
Para simplificar los numeradores, usamos el álgebra de Dirac

(p/ + m)γ ν u(p) = (2pν − γ ν p/ + γ ν m)u(p)

= 2pν u(p) − γ ν (p/ − m)u(p)
= 2pν u(p)

La expresión simplificada, haciendo la sustitución en el numerador y el denominador queda:

 µ ν ′
γ k/γ + 2γ µ pν γ ν k/ γ µ + 2γ ν pµ

2 ∗ ′ ′
iS = −ie ϵµ (k )ϵν (k)ū(p ) + u(p). (3.1)
2p · k 2p · k ′
Con S simplificada podemos romper la suma en dos matrices S1 y S2 donde S1 y S2 son los sumandos
 µ ν
−γ k/γ + 2γ µ pν

2 ∗ ′ ′
iS1 = −ie ϵµ (k )ϵν (k)ū(p ) u(p). (3.2)
2p · k
Nos centramos en S1 ya que el cálculo para S2 es idéntico, toca encontrar el conjugado de S1
 σ ρ
−γ k/γ + 2γ σ pρ

∗ ′ ∗
2
−iS1 = ie ϵρ (k )ϵσ (k)ū(p) u(p′ ). (3.3)
2p · k
Ahora, multiplicamos las ecuaciónes 3.2 y 3.3 para obtener el cuadrado de S1
 µ ν
γ k/γ + 2γ µ pν
 σ ρ
γ k/γ + 2γ σ pρ
 
4 ∗ ′ ′ ∗ ′
2
|S1 | = e ϵµ (k )ϵρ (k )ϵν (k)ϵσ (k)ū(p ) u(p)ū(p) u(p′ ) (3.4)
2p · k 2p · k

3.2.1 Suma sobre los estados de polarización

Actualmente, la expresión mantiene especificados los estados de espı́n y de polarización para elec-
trones y fotones. Sin embargo, experimentalmente es tı́pico que involucre fotones no polarizados
colisionando con electrones en un medio no polarizado, entonces debemos promediar sobre esos
estados. Igualmente, los detectores de fotones y electrones no detectan la polarización, entonces
debemos sumar sobre los estados finales de espı́n. Estamos interesados en la matriz sin polarizar.
1 X 1 XXX 1 X
|S(s, r → s′ , r′ )|2 = |S|2 = S̄ 2 (3.5)
2 s 2 r s′ r ′ 4 espines

La suma de los espines para los estados electrónicos puede ser realizada usando la relación de
completitud: X
us (p)ūs (p) = p/ + m.
s=1,2

22
Mientras, la suma de las polarizaciones para los fotones puede ser reemplazada por la métrica
X
ϵ∗r r
µ (k)ϵν (k) → −gµν
r=1,2

Hacemos la suma sobre los estados de espı́n y queda:

e4
 µ ν
γ k/γ + 2γ µ pν
 σ ρ
γ k/γ + 2γ σ pρ
 
¯ 2 ′
|S1 | = gµρ gνσ (p/ + m) (/p/ + m) . (3.6)
4 2p · k 2p · k

3.2.2 Álgebra de trazas

La ecuación 3.6 muestra que tiene la forma de traza de un producto de 4 matrices 4 - dimensionales:
3
X 3
X
tr(BCDE) = (bcde)f f = bf i cih dhg egf ,
f =0 g,h,i,f =0

entonces sustituyendo la ecuación 3.6 en lo anterior nos queda

e4 1
|S̄1 |2 = A,
4 2p · k
donde
A = tr[(p/′ + m)(γ µ k/γ ν + 2γ µ pν )(p/ + m)(γ σ k/γ ρ + 2γ σ pρ )]
La expansión de A consta de 16 términos, Para simplificar un poco los cálculos usamos
A = I + II + ... + X V + X VI
Donde I,...,X VI son las trazas individuales. Si en el proceso de multiplicación nos movemos desde
el final hasta el comienzo de la expresión, mientras ciclamos dentro de los paréntesis de izquierda
a derecha, encontraremos que los términos III, IV, VII − X .X III, X IV contienen un número
impar de matrices gamma y por lo tanto, estos términos son nulos. Sumando los términos no nulos
uno encuentra que
A = 32[2m4 + 2m2 (p · k) − m2 (p · p′ ) − m2 (p′ · k) + (p · k)(p′ · k)] (3.7)

3.2.3 Variables de Mandelstam

Para escribir las ecuaciones de una manera Lorentziana, usamos las variables de Mandelstam de las
que se profundizó en el capı́tulo anterior
s = (pµ + k µ )2 = (p′µ + k ′µ )2
t = (k µ − k ′µ )2 = (pµ − p′µ )2
u = (pµ − k ′µ )2 = (k µ − p′µ )2
Como las variables satisfacen s + t + u = 2m2e , ası́ como pµ pµ = p′µ p′µ = m2e y k µ kµ = k ′µ kµ′ = 0
podemos usar las relaciones:
1
pµ kµ = (s − m2 )
2
µ ′ 1
p pµ = (s + u)
2
µ ′ 1 2
k pµ = (m − u).
2

23
Esto da, sustituyendo en 3.7:

A = 8[2m2 (2s + t + u) + (s − m2 )(u − m2 )]

= 8[4m4 + 2m2 (s − m2 ) − (s − m2 )(u − m2 )]

4m4 2m2 u − m2
 
2 4
|S1 | = 2e + − (3.8)
(s − m2 )2 s − m2 s − m2
Ahora calculamos el resto de los términos de |S|2 .

3.2.4 Cálculo del elemento de la matriz cuadrada completa

Simplificamos S2 de la misma manera que S1
′
−γ ν k/ γ µ + 2γ ν pµ
 
iS2 = −ie2 ϵ∗ν (k ′ )ϵµ (k)ū(p′ ) u(p). (3.9)
−2p · k ′

Conjugamos 3.9 de la misma manera que en 3.3

′
−γ ρ k/ γ σ + 2γ ρ pσ
 
−iS2∗ 2
= ie ϵρ (k ′
)ϵ∗σ (k)ū(p) ′
u(p′ ). (3.10)
2p · k

Multiplicamos 3.9 y 3.10 entre sı́

′  σ ′ ρ
γ ν k/ γ µ + 2γ ν pµ γ k/ γ + 2γ σ pρ
  
2
|S2 | = e4 ϵ∗ν (k ′ )ϵρ (k ′ )ϵµ (k)ϵ∗σ (k)ū(p′ ) u(p)ū(p) u(p′ ) (3.11)
2p · k ′ 2p · k ′

Comparando con 3.4 notamos la misma estructura en la ecuación 3.11 bajo el reemplazo k → k ′ (s →
u). Esto es, la evaluación es simétrica e inmediatamente nos permite escribir:

4m4 2m2 s − m2
 
2 4
|S1 | = 2e + − (3.12)
(u − m2 )2 y − m2 u − m2

Ahora todo lo que queda es calcular el término de interferencia

 µ ν  σ ′ ρ
−γ k/γ + 2γ µ pν −γ k/ γ + 20γ σ pρ
 
∗ 4 ∗ ′ ′ ∗ ′
2|S1 ||S2 | = e ϵµ (k )ϵρ (k )ϵν (k)ϵσ (k)ū(p ) u(p)ū(p) u(p′ )
2p · k 2p · k ′
(3.13)
Debido a la simetrı́a entre esta expresión y la de los cuadrados de S1 y S2 podemos ver que esta
evaluación es similar a la anterior:

¯ e4 1
∗
2S1 S2 = ′
A′ (3.14)
2 (2p · k)(−2p · k )

Donde A′ es la traza
′
A′ = tr[(p/′ + m)(γ µ k/γ ν + 2γ µ pν )(p/ + m)(−γµ k/ γν + 2γν pµ )]

Cuando expandamos la matriz va a contar con 16 términos

A′ = I ′ + II ′ + ... + X VI ′

24
Haciendo un tratamiento idéntico al que se comentó en la subsección anterior, los 8 términos nulos
son los mismos y los 8 términos no nulos son:

A′ = 32(k · k ′ )(p · p′ ) − 32(p · k)(p · p′ ) + 16m2 (p′ · k) + 32(p · k ′ )(p · p′ ) − 16m2 (p′ · k ′ ) + 16m2 (p · p′ )
(3.15)
−16m2 (k · k ′ ) + 32m2 (p · k) − 32m2 (p · k ′ ) + 16m4
(3.16)

Sustituyendo las variables de Mandelstam en la ecuación 3.16 nos da:

A′ = 8[4m4 + m2 (u − m2 ) + m2 (s − m2 )]
8m4 2m2 2m2
 
¯ ∗
2S1 S2 = 2e 2
+ −
s − m2 u − m2 s − m2

Finalmente, haciendo juntar los términos queda:

" 2 #
2 2
  
¯ 2 = 2e4 4m4 1 1 1 1 u − m s − m
|S| + + 4m2 + − − (3.17)
s − m2 u − m2 s − m2 u − m2 s − m2 u − m2

Para calcular la sección transversal es necesario especificar el marco de referencia

3.3 Evaluación en distintos marcos de referencia

3.3.1 Marco de referencia del centro de masa
En el marco de referencia del centro de masa el cuadrimomento de las partı́culas se escribirá

k = (ω ∗ , 0, 0, ω ∗ )
k ′ = (ω ∗ , ω ∗ sinθ∗ , 0, ω ∗ cosθ∗ )
p = (E, 0, 0, −ω ∗ )
p′ = (E, −ω ∗ sinθ∗ , 0, −ω ∗ cosθ∗ )

Donde ω ∗ y θ∗ son la energı́a del fotón y el ángulo de dispersión en el marco del centro de masa.
Esto hace que las variables de Mandelstam sean

s = 2ω ∗ (E + ω ∗ ) + m2
t = 2ω ∗2 (cosθ∗ − 1)
u = −2ω ∗ (E + ω ∗ cosθ∗ + m2 )

Entonces la sección transversal diferencial se convierte en

1
dσ = |S̄|2 dΠ2 (3.18)
ω∗
4Eω ∗ |1 − |
E
En términos de invariantes en el marco de referencia:
1
dσ = 2
|S̄|2 dΠ2 (3.19)
2(s − m )

25
3.3.2 Evaluación de la integral en el espacio de fases del marco de
referencia del centro de masa
Podemos evaluar la integral referente al estado final en el espacio de fases en el marco del centro de
masa:
d3 p′ 1
Z Z Z 3 ′
dk 1
3
d Π2 = 3
(2π)4 δ (4) (p + k − p′ − k ′ ). (3.20)
(2π) 2Ep′ (2π)3 2Ek′
Integramos sobre las 3 componentes de p′ , lo cual, dado p = −k, establece p′ = −k ′ , entonces con
cambio de variable a coordenadas esféricas en el espacio de fases tenemos:
dk ′ k ′2 dcos(θ∗ )dϕ∗
Z Z Z Z
dΠ2 = δ(E + ω ∗ − Ep′ (k ′ ) − Ek′ (k ′ )). (3.21)
(2π)2 4Ep′ Ek′
Sea como sea, la colisión es simétrica sobre ϕ entonces recogemos otro factor de 2π. Ahora tenemos
una integral sobre momentos, pero una delta de Dirac en términos de energı́a, que son funciones
de momento via la relación de dispersión. Entonces, si f (k ′ ) = E + ω ∗ − Ep′ (k ′ ) − Ek′ (k ′ ) debemos
evaluar
dk ′ k ′2 dcosθ∗ δ(k ′ − k0′ )
Z Z Z
dΠ2 =
8πEp′ Ek′ ∂f (k0′ )
∂k ′
Z ′2 −1
k dcosθ∗ k ′ k′

= +
8πEp′ Ek′ Ep′ Ek′
Z ′2
k dcosθ∗ k′
 
= . (3.22)
8π Ep′ + Ek′
Donde f (k0′ ) = 0. El valor exacto de k0′ no es necesario, pero integrar sobre esta distribución fuerza
la conservación de la energı́a E + ω ∗ = Ep′ + Ek′ . Entonces,
ω∗
Z Z
dΠ2 = dcosθ∗ (3.23)
8π(E + ω ∗ )
También se debe calcular en el marco del centro de masa dt = 2ω ∗2 dcosθ∗ . Entonces,
Z Z
dt
dΠ2 = ∗
16πω (E + ω∗)
Z
dt
=
8π(s − m2 )
Esto nos da una fórmula invariante de Lorentz para la sección transversal:
dσ 1
= |S̄|2 (3.24)
dt 16π(s − m2 )2
Donde |S̄|2 está dada por la ecuación 3.17. Para expresar esto en términos de variables con más
sentido fı́sico debemos cambiar al marco de laboratorio.

3.3.3 Marco de Laboratorio

En el marco de laboratorio el electrón está en reposo y el cuadrimomento de las partı́culas es:
k = (ω, 0, 0, ω)
′
k = (ω, ωsinθ, 0, ωcosθ)
p = (m, 0)
p = (E, −p′ )
′

26
Esto da las variables de Mandelstam en el marco de laboratorio:

s = 2mw + m2
t = 2ww′ (cosθ − 1)

u = −2mw′ + m2 (3.25)
Notamos que en este marco w′ y θ no son variables independientes. Para encontrar la dependencia
derivamos la fórmula de Compton. Conservando el cuadrimomento p + k = p′ + k ′ de nuevo:

p′2 = (p + k − k ′ )2
= p2 + 2p · (k − k ′ ) − 2k · k ′ + k 2 + k ′2
m2 = m2 + 2m(w − w′ ) + 2ww′ (cosθ − 1)
Entonces,
1 1 1
′
− = (1 − cosθ),
w w m
Reordenando para w′ :
w
w′ = w
1 + (1 − cosθ)
m
Esto nos permite expresar t en términos del cosθ :
2w2 (1 − cosθ)
t= w
1 + (1 − cosθ)
m
dt 2w2
= w = 2w′2
dcosθ (1 + (1 − cosθ)) 2
m
Entonces encontramos la sección transversal diferencial respecto al ángulo en el marco del labora-
torio:
dσ dσ dt 1
= = 2w′2 |S̄|2
dcosθ dt dcosθ 16π(2mw)2
Para expresar |S̄|2 en términos del marco de laboratorio sustituimos las expresiones obtenidas en
3.17 para obtener
" 2 #
′
  
cosθ − 1 cosθ − 1 w w
|S̄|2 = 2e4 4m4 + 4m2 + + ′
2m2 2m2 w w
 ′ 
4 w w 2
= 2e + ′ − sin θ
w w

Lo que finalmente produce la fórmula de Klein - Nishina:

 2  ′
πα2 w′

dσ w w 2
= 2 + ′ − sin θ , (3.26)
dcosθ m w w w

donde α es un parámetro conocido, es la relación que existe entre la energı́a del fotón incidente y la
masa - energı́a del electrón, w′ y w son las frecuencias del fotón saliente e incidente respectivamente.

27
Figure 3.2: Densidad de probabilidad con α muy bajo

28
Figure 3.3: Densidad de probabilidad con α todavı́a bajo

29
Figure 3.4: Densidad de probabilidad con α mediano

30
Figure 3.5: Densidad de probabilidad con α alto

31
3.4 Gráficos
Para ver a efectos prácticos cómo funciona la fórmula obtenida en la ecuación 3.26 vamos a graficar
algunos ejemplos de funciones de probabilidad en gráficas polares, ya que el parámetro del que
depende la función de probabilidad es θ, donde los radios son la densidad de probabilidad, note que
la densidad de probabilidad no tiene por qué estar normalizada. Notamos que entre más alto
es α, menor será el ángulo dispersado y mayor será la densidad de probabilidad máxima, esto es
porque α es la relación entre la energı́a del fotón incidente y la masa del electrón [7].

32
Bibliography

[1] D. Tong, “Quantum field theory, university of cambridge part iii mathematical tripos,” Quantum
Field Theory. University of Cambridge, p. 71, 2012.

[2] N. Zettili, Quantum mechanics concepts and applications. A John Wiley and Sons, Ltd.s, 2009.

[3] A. Belyaev, “Path integral in quantum field theory.”

[4] F. Gross, Relativistic quantum mechanics and field theory. John Wiley & Sons, 1999.

[5] M. E. Peskin, An introduction to quantum field theory. CRC press, 2018.

[6] H. Nastase, Introduction to Quantum Field Theory. Cambridge University Press, 2019.

[7] B. Robson and S. Sutanto, “Relativistic wave equations and compton scattering,” International
Journal of Theoretical Physics, vol. 40, pp. 1491–1499, 2001.

33