Resumen Iar

Resumen IAR
Introducción a la Astrofı́sica Relativista

Iara Pintos
1. El espacio-tiempo y la relatividad especial

Hay dos tipos de propiedades de las cosas: intrı́nsecas y relacionales. Las primeras solo dependen
de la cosa en cuestión (por ejemplo, la carga eléctrica de una partı́cula) mientras que las segundas
dependen también de otras cosas (por ejemplo, la velocidad de la partı́cula). Cuando las cosas se
combinan para formar nuevas cosas, las cosas resultantes pueden tener propiedades emergentes,
que las constituyentes no tienen. A su vez, las propiedades emergentes pueden ser intrı́nsecas o
relacionales. Ası́, por ejemplo, un gas puede tener temperatura y presión, propiedades de las que
carecen las moléculas constitutivas.
Dada una cosa x llamaremos P al conjunto de todas sus propiedades:
P = {p/px }.
Los elementos de P pueden ser representados por funciones matemáticas (supuesto metodológico
de la Ciencia). Llamamos espacio de estados S(x) de una cosa x al conjunto de funciones (de dominio
M) que representan a los elementos de P. Una ley es una restricción sobre S(x). Nos dice que las
propiedades de una cosa no pueden tomar cualquier valor. Llamaremos SL (x) al conjunto de estados
legales de x. Estos son los estados accesibles en principio a la cosa x de acuerdo con las restricciones
legales que imperan sobre ella. El estado real de una cosa concreta x es un punto de Sl (x).
Un cambio es un par ordenado de estados de la cosa que cambia:
(s1 , s2 ) ∈ EL (x) = SL (x) × SL (x).
El conjunto de todos los cambios de una cosa es el espacio de eventos EL (x) de esa cosa. Definimos
el espacio-tiempo de la siguiente manera:
El espacio-tiempo es un sistema fı́sico que contiene a la colección de todos los eventos.
Todo lo que ha ocurrido, ocurre o ocurrirá a alguna cosa es un punto (elemento) del espacio-
tiempo. Un proceso (sucesión de cambios) es una lı́nea (o subconjunto) del espacio-tiempo.
Al ser el espacio-tiempo una colección y no un mero conjunto es una entidad fı́sica y no un mero
concepto matemático.
Debemos ahora caracterizar matematicamente al espacio-tiempo si queremos hacer predicciones
precisas sobre ciertos eventos. Postulamos
El espacio-tiempo se representa por una variedad real cuadri-dimensional diferenciable.
Una variedad real se puede pensar como un conjunto (estrictamente un espacio topológico) que
puede ser completamente cubierto por subconjuntos cuyos elementos pueden ser puestos en corres-
pondencia 1 a 1 con subconjuntos de R4 (si la variedad es cuadri-dimensional; Rn si es n-dimensional).
En forma estricta, M es una variedad real n-dimensional diferenciable si y solo si:
M es un conjunto.
1
∃O/O = {Oα ⊂ M }
Todo elemento p ∈ M es tal que ∃Oα ∈ O/p ∈ Oα .
∀α, ∃Φα : Oα → Uα , con Uα subconjunto abierto de Rn .
Si existen dos conjuntos O1 y O2 /O1 ∩ O2 ̸= ∅ =⇒ ∃Φ1 · Φ−1
1 que pone en correspondencia 1
a 1 los puntos de U1 ⊂ R con los de U2 ⊂ R .
n n
De esta definición se desprende el porqué se postula una variedad para representar el espacio-
tiempo: independientemente de la estructura geométrica de este, pueden adoptarse coordenadas
(números reales) para describir procesos que ocurren en él.
Dado un elemento p ∈ M podemos describirlo usando distintos sistemás de coordenadas. Por
ejemplo,
p ⇐⇒ {xµ },
p ⇐⇒ {x′µ },
∃x′µ = x′µ ({xµ }),
donde µ = 0, 1, 2, 3. Adoptamos cuatro coordenadas porque el mundo parece ser 4-dimensional (con
tres dimensiones espaciales y una temporal), pero en principio no hay limitaciones en ese sentido.
Para poder caracterizar objetos fı́sicos sobre nuestra variedad (esto es, para poder describir pro-
cesos reales), necesitamos definir sobre ella objetos matemáticos que puedan ser utilizados luego en
la representación objetos fı́sicos y sus propiedades.
1.1. Objetos y estructuras sobre la variedad

Los objetos sobre la variedad se definen por sus propiedades de transformación frente a cambios
de coordenadas {xµ } ⇐⇒ {x′µ }. El objeto más simple es un escalar:
ϕ(xµ ) = ϕ′ (x′µ ).
El valor del escalar no cambia cuando el sistema coordenado cambia de {xµ } a {x′µ }, pero la
forma sı́ puede cambiar.
Introduzcamos ahora un objeto de cuatro componentes Aµ . Si ante un cambio de coordenadas
{x } → {x′µ }, este se transforma como
µ
∂x′µ
A′µ = Aν ,
∂xν
entonces Aµ es un vector contravariante. Ejemplo: un vector contravariante es el que une dos puntos
arbitrariamente próximos de la variedad
∂x′µ ν
dx′µ = dx .
∂xν
Los vectores covariantes se definen por medio de la siguiente ley de transformación ante cambios
de coordenadas:
∂xν
Bµ′ = Bν .
∂x′µ
Un ejemplo es el gradiente de un campo escalar;
∂ϕ ∂xµ ∂ϕ
=
∂x′ν ∂x′ν ∂xµ
2
A los vectores sobre la variedad se los puede llamar tensores de rango 1. Podemos definir tensores
contravariantes (o covariantes) de rango arbitrario:
n
∂x′µ ∂xσ
z }| {
′ ...µ...
T...ν... = ... ρ ... ′ν ...T...σ....
...ρ...
| {z } ∂x ∂x
m
El tensor T es n veces contravariante y m veces covariante.

Ejemplo de tensor de rango 2: tensor de energı́a-impulso que caracteriza a cualquier sistema fı́sico.
Un campo tensorial definido sobre alguna región de la variedad es una asociación de un tensor de
la misma variedad a cada punto de la región:
...ρ...
p → T...σ.... (p),
donde T...σ....
...ρ...
(p) es el valor del tensor en p. El campo tensorial se llama continuo o diferenciable si las
componentes del tensor lo son.
Debido a sus propiedades de transformación, toda relación entre tensores mantendrá su forma al
cambiar de sistemas de coordenadas. Como las leyes de la fı́sica se espera que tengan la misma forma
en todos los sistemas coordenados, deben ser expresadas en forma de ecuaciones tensoriales. En tal
caso, se dice que son covariantes. La covarianza es validez de las ecuaciones bajo transformaciones
generales de coordenadas.
Aunque hemos definido objetos sobre la variedad que representa al espacio-tiempo, aún no sabe-
mos medir distancias entre puntos (eventos) de la variedad. Para eso debemos asignar una estructura
geométrica a la variedad. Esto se hace introduciendo un tensor de rango 2 llamado tensor métrico.
El mismo nos dice como calcular la distancia ds entre dos eventos arbitrariamente próximos ϵ1 y ϵ2 :
ds2 = gµν dxµ dxν ,
donde gµν es el tensor métrico. En un espacio-tiempo euclı́deo, por ejemplo
+1 si µ=ν
(
gµν = δµν = y ds2 = (dxµ )2 .
X
0 si µ ̸= ν µ
El tensor métrico gµν tiene un contenido empı́rico: depende de la naturaleza del Mundo. Minkowski
propuso que el tensor métrico del espacio-tiempo es un tensor pseudo-euclı́deo de rango 2 y traza -2:
1 0 0 0
 
0 −1 0 0
gµν = ηµν = ,

ηµα η αν = δµν ,
0 0 −1 0 

0 0 0 −1
conocido como tensor de Minkowski. Aquı́, el intervalo resulta
ds2 = ηµν dxµ dxν = (dx0 )2 − (dx1 )2 − (dx2 )2 − (dx3 )2 .
La geometrı́a resultante es pseudoeuclı́dea. Y ηµν tiene el mismo valor sobre toda la variedad.
Las coordenadas i = 1, 2, 3 se denominan espaciales y la coordenada que aparece con el signo
opuesto es llamada temporal x0 = ct. Vemos que c es una constante que permite uniformar las
dimensiones y su valor coincide con el de la velocidad de la luz en el vacı́o.
3
La introducción del intervalo nos permite, en cada evento, dividir la variedad en tres regiones
bien definidas:
ds2 < 0 : región tipo espacio,
ds2 = 0 : región tipo luz,
ds2 > 0 : región tipo tiempo.
Solo las dos ultimas regiones son accesibles a sistemas fı́sicos. La región tipo luz, en partı́cular,
solo a sistemas que se muevan a la velocidad de la luz. Si reescribimos:
dx2 + dy 2 + dz 2
!
ds = dt
2 2 2
c − = dt2 (c2 − v 2 ),
dt2
con v la velocidad del sistema, obtenemos las siguientes equivalencias
ds2 < 0 ←→ v > c,

ds2 = 0 :←→ v = c,
ds2 > 0 :←→ v < c.
Como nunca ocurre que v > c, los eventos en la región tipo espacio no pueden estar causalmente
ligados a eventos de las otras regiones.
Esto se puede representar gráficamente mediante el
cono de luz. Al graficar eliminamos la dimensión z y
c = 1. Las partı́culas materiales (v < c) tienen un estado
ϵ0 coincidente con el vértice del cono solo pueden sufrir
cambios o procesos que los lleven a estados que están
dentro del cono (afuera del cono v > c). Las partı́culas
sin masa se mueven a la velocidad de la luz y evolucionan
sobre la superficie del cono.
La introducción de una métrica sobre la variedad ha
permitido especificar la estructura causal de la misma,
esto es, identificar qué eventos o sucesos pueden estar o
no causalmente relacionados respecto a cualquier punto
de la variedad.
Usando la métrica podemos, además, diferenciar aho-
ra, pasado y futuro de un dado evento. Un vector xµ
dentro del cono señala hacia el futuro si
Figura 1: Diagrama de un cono de luz. ηµν xµ T ν > 0,
con T ν = (1, 0, 0, 0). En forma similar, xµ señala hacia el pasado si
ηµν xµ T ν < 0.
Podemos definir el tiempo propio de un sistema fı́sico que se mueve con velocidad v respecto de
un cierto sistema coordenado como
1
!
v2
dτ = 2 ds2 ,
2
dτ = dt
2 2
1− 2
c c
que es el tiempo que mide un reloj fijo al sistema fı́sico.
4
Introduciendo el factor de Lorentz γ = (1 − β 2 )−1/2 , con β = v/c, reescribimos
dt
dτ = .
γ
Como γ ≥ 1, el tiempo respecto al sistema propio se dilata. Si β = β(t), tenemos

ˆ t2
dt
τ= .
t1 γ(t)
1.2. El grupo de Lorentz

El grupo más general de transformaciones lı́neales y homogénas entre dos sistemas de referencia
inerciales es el de Lorentz, que deja invariante la métrica de Minkowski:
∂x′µ
L = {xµ → x′µ = Lµν xν }, con Lµν = . .
∂xν
Como la invarianza del intervalo exige la invarianza del tensor métrico, se cumple que
ds2 = ηµν dx′µ dx′ν = ηµν Lµα Lνβ = ηαβ dxα dxβ ,
entonces
ηαβ = ηµν Lµα Lνβ ,
y el elemento neutro del grupo de Lorentz es δνµ
∂xµ
δνµ = = η µα ηνα .
∂xν
La matriz Lµν no es singular y por lo tanto invertible y el elemento inverso es la inversa de la
matriz Lµν . Además, las transformaciones son únicas.
El producto escalar entre dos vectores U y V resulta un invariante
Uα′ V ′α = Lβα Uβ Lαµ V µ = Uµ V µ
Hay cuatro casos posibles de conjunto de transformaciones de Lorentz posibles:
1. det(L) = 1 y L00 ≥ 1 es el grupo propio de Lorentz,
2. det(L) = 1 y L00 ≤ 1 son las inversiones espacio temporales,
3. det(L) = −1 y L00 ≥ 1 son las inversiones espaciales y
4. det(L) = −1 y L00 ≤ 1 son las inversiones temporales.
El grupo de Lorentz es un subgrupo propio del más general grupo de Poincaré. Este consis-
te en las transformaciones lı́neales inhomogéneas que dejan la métrica ηµν invariante. Se trata de
transformaciones de Lorentz más una traslación arbitraria en el espacio tiempo
P = {xµ → x′µ = Lµν xν + tµ }
5
1.3. Mecánica relativista
Consideremos una partı́cula moviéndose en el espacio-tiempo. Su trayectoria es
xµ = xµ (τ ).
La acción de una partı́cula es ˆ 2

S= L dt,
1
donde L = T − U es el Lagrangiano de la partı́cula, siendo T y U la energı́a cinética y potencial,
respectivamente. El Lagrangiano tiene unidades de energı́a y la acción, de energı́a por tiempo. Su
interpretación es que es una medida de la cantidad de cambios que ocurren al sistema correspondiente
al Lagrangiano L en un determinado lapso. Sobre trayectorias reales se cumple δS = 0. Esto refleja
un principio general de la naturaleza: el Principio de Mı́nima Acción. La evolución dinámica de un
sistema fı́sico, ir de un estado a otro, ocurre de tal forma que el número de cambios requeridos es el
mı́nimo posible.
La partı́cula de nuestro ejemplo se mueve libremente sobre una geodésica y su acción es
ˆ τ2
c2 dt2
S = −α ds, ds2 = ηµν dxµ dxν = 2 ,
τ1 γ
y vemos que ˆ τ2 √ √
S = −α c2 − v 2 dt =⇒ L = −α c2 − v 2 .
τ1
Para velocidades bajas, hacemos un desarrollo en Taylor alrededor de v = 0 y vemos que
L ≈ −αc + 21 αc v 2 + ... expresión relativista,

)
=⇒ α = mc.
L = 12 mv 2 expresión newtoniana
Definiciones:
Cuadrivector velocidad
dxµ
vµ ≡ = (γc, γ⃗v ) = ẋµ .
dτ
Cuadri-aceleración:
d 2 xµ
v = aµ ≡ = ẍµ .
dτ 2
Cuadri-impulso
P⃗ = γm⃗v ,
(
∂L
pµ = − µ =⇒
∂ ẋ p0 = γmc
Energı́a (usando el Hamiltoniano H)
E = H = P⃗ · ⃗v − L = γmc2
Con la definición de energı́a encontramos

E ⃗

P ≡ (γmc, γm⃗v ) =
µ
,P , P µ = η µν Pν = mv µ .
c
6
Notemos que
P µ Pµ = m2 c2 =⇒ E 2 = m2 c4 + c2 p2 ,
que es la expresión para la energı́a de una partı́cula libre y tiene un termino de energı́a en reposo y
otro dinámico.
Para hallar la ecuación de movimiento variamos la acción y encontramos que
dv α dP µ
= 0 o también = 0.
dτ dτ
Si hay una fuerza externa
dP µ
fµ = .
dτ
Vemos que cuando no hay fuerzas, el momento se conserva.
Siempre que el espacio-tiempo tenga métrica Minkoskiana podremos obtener reglas de transfor-
mación para todas nuestras cantidades. Pero la geometrı́a real del espacio tiempo parece ser pseudo-
Riemanniana (esto es, solo el espacio tangente en cada punto es Minkoskiano y la curvatura global
es no nula, dependiendo del contenido de energı́a e impulso de la materia).
1.4. Elementos de la relatividad general

La distribución de materia existente en el Universo determina la geometrı́a del mismo, a través
del campo tensorial métrico gµν (x). El potencial gravitacional se representa a través de gµν (x), que
se obtiene resolviendo las ecuaciones de Einstein-Maxwell:
1 8πG mat
Rµν − Rgµν + Λgµν = 4 Tµν + Tµν
EM
,
2 c
con Λ la constante cosmológica, Rµν = g λσ Rλσµν y R = g µν Rµν son el tensor y el escalar de Ricci y
Rλσµν es el tensor de Riemann.
Las fuentes de la curvatura del espacio-tiempo son la materia y los campos electromagnéticos
que esta produce si está cargada. Es decir, las fuentes de las ecuaciones de E-M son el tensor de
energı́a-impulso de la materia, que en el caso de un fluido ideal tiene la forma
mat
Tµν = (ρ + p)uµ uν + pgµν ,
1 1

Tµν =
EM α αβ
Fµα Fν − Fαβ F gµν ,
4π 4
con Fµν = ∂m uAν − ∂ν Aµ .
Las ecuaciones de conservación de energı́a y momento angular se pueden derivar de
∇ν T µν = 0.
1.4.1. Agujeros negros

Los agujeros negros son regiones del espacio-tiempo limitadas por un horizonte de eventos, tal
que las cosas que crucen el horizonte hacia el agujero negro no podrán volver a cruzarlo en el sentido
opuesto.
Las principales soluciones analı́ticas conocidas de las ecuaciones de Einstein-Maxwell que descri-
ben agujeros negros son
7
Métrica de Schwarzschild: corresponde a un espacio-tiempo determinado por una masa M que
no rota y con carga neta nula.
Métrica de Kerr: corresponde a un espacio-tiempo determinado por una masa M en rotación
y con carga neta nula.
Métrica de Reissner-Nordstrom: corresponde a un espacio-tiempo determinado por una masa
M cargada pero con momento angular nulo.
Métrica de Kerr-Newman: corresponde a un espacio-tiempo determinado por una masa M en
rotación y cargada. Es la solución más general.
Solo la solución de Kerr se piensa que corresponde a agujeros negros que existen en la naturaleza.
2. Partı́culas elementales <3

Las ideas actuales acerca de los elementos básicos que forman todas las cosas que hay en el
Universo están expresadas en el llamado Modelo Estándar de las partı́culas elementales. Según el,
hay 12 tipos diferentes de partı́culas que forman todas las cosas. Estas partı́culas se dividen en dos
grupos: quarks y leptones, que se dividen en seis y seis. Existen además cuatro bosónes de gauge,
que son las partı́culas mediadoras de las interacciónes electromagnética, débil y fuerte. La ultima
partı́cula que compone el modelo es el bosón de Higgs.
Figura 2: Partı́culas que componen el modelo estándar: tres generaciones de fermiones (seis quarks
y seis leptones), los cuatro bosones de gauge y el bosón de Higgs.
8
Los seis tipos (usualmente llamados sabores) de quarks se denominan, en orden de masa creciente,
up (u), down (d), strange (s), charm (c), bottom (b) y top (t). Los seis tipos de leptones son el electrón
(e− ), el neutrino electrónico (νe ), el muón (µ− ), el neutrino muónico (νµ ), el tau (τ − ) y el neutrino
tau (ντ ).
Las partı́culas elementales tienen, a la vez, propiedades elementales. Estas propiedades consis-
ten en la capacidad de cambiar el estado de otras partı́culas a través de fuerzas o interacciónes
elementales. Existen cuatro interacciónes elementales:
gravedad,
electromagnetismo,
fuerza fuerte y
fuerza débil.
Todas estas interacciónes actúan a distancia pero con diferente rango o alcance:
La gravedad tiene rango infinito pero es muy débil. Actúa sobre todas las partı́culas y es una
fuerza atractiva. Es la interacción que determina la estructura a gran escala del Universo.
El electromagnetismo también tiene rango infinito pero es una interacción más fuerte. Actúa
sobre aquellas partı́culas que tienen carga eléctrica (que puede ser negativa o positiva). El
número total de cargas positivas es igual al número total de cargas negativas, por lo que sobre
distancias grandes la fuerza electromagnética se anula. El electromagnetismo es la interacción
dominante a escala humana, pero no tiene influencia sobre la estructura a gran escala del
Universo.
La interacción fuerte es de corto alcance (∼ 10−15 m ≡ 1 fermi). Solo actúa sobre aquellas
partı́culas que tienen carga fuerte, usualmente llamada color. Los quarks tienen color, pero
no ası́ los leptones. La interacción fuerte cohesiona a los quarks para formar partı́culas más
complejas llamadas hadrones.
La interacción débil es de rango aún más corto que la fuerte (∼ 10−18 m) y actúa tanto sobre los
quarks como sobre los leptones (además de la gravedad es la única fuerza que actúa sobre los
neutrinos), aunque con ciertas restricciones. Los quarks pueden cambiar su sabor únicamente
por acción de la interacción débil. La propiedad cuántica de las partı́culas que determina su
acoplamiento con la interacción débil es el isoespı́n débil; es el análogo a la carga eléctrica y al
color para las interacciónes electromagnética y fuerte, respectivamente.
Respecto a la intensidad relativa de las distintas interacciónes, consideramos los efectos sobre dos
protones en contacto
Fuerza fuerte 1
Fuerza electromagnética 10−2
Fuerza débil 10−7
Fuerza gravitacional 10−39
Como la interacción fuerte no actúa sobre los leptones, estos no forman partı́culas estables más
complejas.
9
Los leptones y quarks pueden dividirse, además,
en tres generaciones, como muestra el modelo estan-
dar (Figura 2). Cada generación agrupa dos lepto-
nes (uno de carga -1 y un neutrino) y dos quarks
(uno de carga 2/3 y otro -1/3). Respecto a la ma-
sa, cada partı́cula en una dada generación tiene una
masa mayor que la de la partı́cula correspondiente
de la generación anterior. Debido a esta diferencia
en las masas, las partı́culas de las generaciones II y
III decaen en las de la generación I, que son estables
y componen la materia ordinaria.
El modo de interacción débil depende de la ge-
neración a la que pertenezcan las partı́culas. En el
caso de los leptones, solo actúa entre ellos que pertenezcan a la misma generación. Respecto de los
quarks, es la única interacción que puede cambiar su sabor. Las transiciones más probables son aque-
llas entre quarks de la misma generación. Las transiciones entre generaciones están permitidas pero
son improbables (están suprimidas).
2.1. Leptones
Los leptones son las partı́culas elementales que no se acoplan a la interacción fuerte. Los seis
leptones se pueden clasificar de acuerdo con sus propiedades: la carga eléctrica Q (el electrón, el
muón y el tau tienen carga eléctrica Q = −1 y los neutrinos, Q = 0) y el número leptónico L, que es
una propiedad intrı́nseca de este tipo de partı́culas. Hay tres tipos diferentes: Le , Lµ y Lτ , cada uno
caracterı́stico de una generación. La masa de los neutrinos es muy pequeña comparada con la de los
otros leptones, pero no nula.
Toda interacción que involucre leptones debe conservar tanto la carga eléctrica como el número
leptónico.
Generación leptón Q Le Lµ Lτ mc2

I e−1 -1 1 0 0 0.511 MeV
νe 0 1 0 0 < 2 eV
II µ−1 -1 0 1 0 105.66 MeV
νµ 0 0 1 0 <0.19 MeV
III τ −1 -1 0 0 1 1776.86 MeV
ντ 0 0 0 1 <18.2 MeV
2.1.1. Antipartı́culas
La teorı́a de la relatividad implica que. por cada leptón ℓ con número leptónico Lℓ = 1 y carga
eléctrica Qℓ = −1 debe existir un antileptón ℓ con número leptónico Lℓ = −1 y carga eléctrica
Qℓ = −1
El concepto de antipartı́cula surge naturalmente en relatividad especial a partir de la relación
entre energı́a y momento para una partı́cula libre
1/2
E = ± m2 c4 + c2 p2
10
2.2. Hadrones
Los quarks (q) se combinan para formar partı́culas no elementales llamadas hadrones. Tres quarks
forman un barión (qqq), tres antiquarks forman un antibarión (qqq) y un quark y un antiquark forman
un mesón (qq).
El protón y el neutrón, llamados en conjunto núcleones por componer el núcleo atómico, son
ejemplos de bariónes. El protón está formado por dos quarks u y un quark d (p ≡ uud) y el neutrón,
por dos quarks d y un quark u (n ≡ udd).
El protón es el barión más ligero y es estable. Todos los demás bariónes son inestables y decaen.
Como el protón tiene carga eléctrica Qp = +1 y el neutrón es neutro, de lo anterior se deduce que
los quarks tienen carga eléctrica fraccionaria (respecto de la del electrón). En partı́cular,
2 1
Qu = + , Qd = −
3 3
La fuerza débil no puede cambiar un quark en un leptón. De aquı́ que los quarks tengan su
propio número llamado número bariónico. Vale 1/3 para los quarks y 0 para los leptones. El número
bariónico se conserva en todas las interacciónes. Los antiquarks tienen carga y número bariónico
opuesto a los quarks.
Combinando los tres quarks más ligeros se pueden formar distintos bariónes. El llamado octeto
bariónico esta formado por las siguientes partı́culas, donde señalamos la extraneza S (la extrañeza
de un hadrón decrece en una unidad por cada quark s que contenga) y las masas en GeV.
Nucleones p ≡ uud (0.938 GeV), S=0

n ≡ udd (0.940 GeV)
Σ+ ≡ uus (1.189 GeV)
Sigma Σ0 , Λ ≡ uds (1.192 GeV) S = −1
Σ− ≡ dds (1.197 GeV)
Xi Ξ0 ≡ uss (1.314 GeV) S = −2
Ξ− ≡ dss (1.32 GeV)
A cualquier hadrón que tenga extrañeza no nula se la llama partı́cula extraña. A los hadrones
extraños que decaen que no contengan quarks c, t ni b se los suele llamar hiperiones.
Los quarks dentro de los bariónes tienen una variedad de niveles energéticos, por lo que los
estados excitados dan lugar a partı́culas más masivas, llamadas resonancias. Las más conocidas son
las partı́culas ∆, de extrañeza S = 0. También tenemos partı́culas Σ∗ con resonancias más masivas y
S = −1. Además, podemos formar un barión solo con quarks extraños, la partı́cula Ω− con S = −3,
Ω− ≡ sss.
Los mesones son hadrones formados por un quark y un antiquark, que no tienen por qué ser del
mismo sabor. Los más relevantes (o ligeros) son:
Pion: π − ≡ du, π 0 ≡ uu, π + ≡ ud, con S = 0.
Eta: η ≡ dd, η ′ ≡ ss, con S = 0.
Kaon:
• Con extrañeza S = 1: K 0 ≡ ds, K + ≡ us.

0
• Con extrañeza S = −1: K − ≡ su, K ≡ sd.
11
2.3. Partı́culas de interacción
Las diferentes cargas (carga eléctrica, sabor, color) corresponden a excitaciones de campos que
existen en el espacio-tiempo. Las perturbaciones de estos campos son discretas y forman las partı́culas
llamadas bosones de gauge. Estas partı́culas son las mediadoras de las interacciones fundamentales.
Sus propiedades son
Partı́cula Q mc2 [GeV]

Fotón (γ) 0 0
Gluón (g) 0 0
W± ±1 80.379 ± 0.012
Z0 0 91,1876 ± 0,0021
El cuanto de interacción del campo electromagnético es el fotón. El fotón, a diferencia de las

partı́culas cargadas que producen el campo, no tiene carga eléctrica. La masa del fotón es cero.
El cuanto de interacción del campo fuerte es el gluón. El gluón tiene masa nula pero es portador
de color, por lo que experimenta la fuerza.
Las perturbaciones del campo débil son los bosónes cargados W + y W − , que tienen carga eléctrica
Q ± 1. Además, hay un bosón neutro, el Z 0 , que se manifiesta en condiciones de unificación con el
electromagnetismo. A diferencia del fotón y el gluón, los mediadores de la fuerza débil son muy
masivos (casi 100 veces más masivos que el protón).
Finalmente, el hipotético cuanto de interacción del campo gravitatorio es el gravitón. El gravitón,
de existir, serı́a una partı́cula sin carga eléctrica ni masa y que sentirı́a los efectos del campo gra-
vitatorio. La relatividad general, sin embargo, no describe a la gravitación como el resultado de un
campo, sino como efecto de la curvatura del espacio-tiempo. De allı́ que la relatividad general no
sea compatible con un bosón de gauge. El gravitón deberı́a aparecer en una teorı́a que tenga como
lı́mite de bajas energı́as a la relatividad general, pero que sea en sı́ una teorı́a de campo y no del
espacio-tiempo.
2.4. El bosón de Higgs

El bosón de Higgs es el cuanto del campo de Higgs. Este campo permea el espacio-tiempo, tiene
un valor de vacı́o no nulo e interaccióna con las partı́culas elementales otorgándoles masa a través de
un mecanismo conocido como mecanismo de Higgs. Explica la masa no nula de los fotones, W ± y Z
y, de manera algo distinta, da cuenta de la masa de los quarks y leptones.
La mayor parte de la masa de las partı́culas no elementales (y de los cuerpos macroscópicos) no
tiene su origen en la interacción con el campo de Higgs. Se debe, en cambio, a la energı́a cinética de
las partı́culas elementales que la componen y a la energı́a de ligadura que las mantiene unidas.
El bosón de Higgs es una partı́cula inestable, de carga nula, spin cero y de masa no nula, por
lo que es su propia antipartı́cula. La masa del bosón de Higgs no es predicha por el modelo, es un
parámetro libre del mismo y debe medirse experimentalmente.
2.5. Partı́culas más allá del Modelo Estándar

Hay varias razones para pensar que el Modelo Estándar es incompleto y existen partı́culas no
previstas en el mismo. La principal de estas razones es que las observaciones astronómicas revelan
la existencia de anomalı́as gravitacionales: la medición de las curvas de rotación de las galaxias a
través de las lı́neas del hidrógeno neutro revelan que las mismas no obedecen a la ley de gravitación
12
Newtoniana (ni a la de Einstein) a menos que exista materia oscura. Esta materia es materia que
no emite gravitacionalmente ni sufre interacciónes hadrónicas. No se trata de neutrinos, ya que las
observaciones indican que si existe su ecuación de estado es no relativista.
Entonces, se postuló la existencia de partı́culas no incluidas en el Modelo Estándar y que solo
interactúan débilmente con el resto de la materia. Estas partı́culas suelen denominarse WIMPs
(Weakly Interacting Massive Particles). Si esas partı́culas existen, deberı́an ser mucho más númerosas
que los bariones, tal ası́ que se estima que esta materia constituirı́a el 26 % del universo, mientras
que los bariones, apenas el 4 %. El resto se deberı́a a un campo de densidad de energı́a negativa.
Alternativamente, se ha sugerido que la gravedad podria comportarse en forma diferente de lo
postulado por Einstein y Newton, pero no está clara la situacion.
2.6. Interacciones entre partı́culas

Las partı́culas elementales interactúan cuando se acercan lo suficiente como para que las fuerzas
fundamentales puedan actuar modificando su estado. Las interacciónes entre partı́culas involucran
el intercambio entre ellas de bosones de gauge, que son los mensajeros o portadores de las fuerzas
fundamentales.
Una interacción puede tener dos resultados:
La trayectoria de las partı́culas en el espacio-tiempo se modifica.
Las partı́culas cambian y/o nuevas partı́culas aparecen.
El resultado de la interacción queda completamente determinado por las condiciones iniciales,

las fuerzas fundamentales y las leyes de conservación. Cuanto más energéticas sean las
partı́culas, mayor será el rango de posibilidades para la actuación de las diferentes fuerzas.
Todas las interacciónes deben conservar el cuadri-impulso y la carga eléctrica. Las interacciónes
que involucran hadrones conservan el número bariónico y las que involucran leptones, el número
leptónico. El sabor no se conserva, ya que la fuerza débil puede cambiarlo.
Consideremos una reacción permitida (que conserva todos los números cuánticos necesarios)
A + B → C + D.
Entonces, cualquier reacción donde una de las partı́culas pase al otro lado como su antipartı́cula
también esta permitida, por ejemplo
A → B + C + D, (2.1)
A + C → B + D, (2.2)
y la reacción inversa también esta permitida
C + D → A + B.
Sin embargo, hay que tener en cuenta que algunas de estas reacciones puede no ocurrir ya que no
se conserva el cuadri-momento. Se dice entonces que están cinemáticamente prohibidas. Por ejemplo,
si la masa de A es menor que la de la suma de B, C y D la reacción (2.1) no es posible. De la misma
manera, si A y C son menos masivas que B y D, la reacción de la ecuación (2.2) no podrá ocurrir a
menos que la energı́a cinética inicial supere cierto valor umbral.
Reacciones importantes que involucran leptones:
13
El muón puede decaer en un electrón y dos neutrinos
µ− → e− + νµ + ν e .
Decaimiento β del neutron:

n → p + e− + ν e .
Conversion de un proton en un neutron:
p → n + e+ + νe ,
pero para protones libres esta reacción esta cismáticamente prohibida porque la masa del
neutrón es mayor que la del protón, pero puede ocurrir dentro del núcleo atómico.
También podemos encontrar las reacciones de aniquilación
e− + e+ → γ + γ,
e− + e+ → µ+ + µ− ,
e− + e+ → τ + + τ − .
Consideremos la interacción de dos protones. La reacción más simple es la dispersión o scattering
p + p → p + p.
Cuando se incrementa la energı́a son posibles nuevos canales:
p + p → p + p + π0,
p + p → p + p + π+ + π−
y analicemos la primera en términos de quarks
p + p → p + p + π0,
       
u u u u !
u
u + u → u + u +
       
u
d d d d
Parte de la energı́a fue usada para crear un par quark-antiquark (formando un mesón). Notemos
que la carga y el número bariónico se conservan.
Al aumentar la energı́a de los protones incidentes, puede aumentar la multiplicidad de los piones
producidos. En general,
p + p → p + p + aπ 0 + b(π + + π − ), a, b ∈ N0 .
Si aumenta aun más la energı́a pueden aparecer partı́culas extrañas
p + p → p + Σ+ + K 0 ,
0
p + p → p + n + K+ + K .
La colisión suele ser tal que la mayor parte de la energı́a se usa para crear un estado excitado de
uno de los protones, que al decaer produce un leading pión que transfiere la energı́a a rayos gamma.
p + p → p + ∆+ + aπ 0 + b(π + + π − ), a, b ∈ N0 .,
→ p + p + π 0 + aπ 0 + b(π + + π − ),
A los productos mesónicos diferentes del leading pión se les suele llamar fireball: F = aπ 0 +b(π + +
π−)
14
2.6.1. Decaimiento de partı́culas
Un decaimiento es el proceso por el cual la energı́a de una partı́cula aislada se transforma en un
cierto número de partı́culas con menor energı́a. Un decaimiento también puede ocurrir cuando una
partı́cula pasa a un nivel de energı́a menor. Podemos enumerar diferentes tipos de decaimientos:
Decaimiento electromagnético
Dos hadrones pueden estar formados por los mismos quarks pero en distintos niveles de energı́a.
El decaimiento se produce por la emisión de ese exceso de energı́a en forma de un fotón. Ejemplo:
Σ0 → Λ + γ.
uds uds
Otro tipo:
∆+ → p + γ,
∆0 → n + γ.
En algunos casos el decaimiento electromagnético opera por aniquilación directa de un par quark-
antiquark. Un decaimiento de este tipo muy relevante en astrofı́sica es el del pión neutro en dos
fotones
π 0 → γ + γ.
uu
Los decaimientos mediados por la interacción electromagnética ocurren sobre los tiempos tı́picos
∼ 10−10 − 10−16 s. Por ejemplo, la vida media del π 0 es τπ0 = (8,52 + 0,18) × 10−17 s.
Decaimientos fuertes
Se trate de decaimientos en los cuales hay creación de quarks. Ejemplo:
∆+ → p + π 0
uud uud uu
Aquı́, el par uu es creado por gluones.
Los decaimientos mediados por la fuerza fuerte ocurren sobre tiempos extremadamente cortos
(∼ 10−25 s) en comparación con los decaimientos electromagnéticos.
Decaimientos débiles
Son decaimientos en los cuales hay creación de quarks por acción de los bosones W ± . Por ejemplo,
cuando un quark s (Q = 1/3) se convierte en un u (Q = 2/3) la diferencia de masa y carga escapa
como un W − . Cuando un quark u cambia en un d (Q = −1/3) se emite un W + .
La energı́a de los bosones W ± puede ser usada para crear pares quark-antiquark.
Decaimiento del neutrón
Un neutrón aislado tiene una masa mayor que la del protón y es por lo tanto inestable. Decae
según la siguiente reacción
n → p + e + νe
0 = +1 + (−1) + 0 (carga)
1 = 1 + 0 + 0 (número bariónico)
0 = 0 + 1 − 1 (número leptónico)
15
Este es un decaimiento débil en el que se producen leptones. La vida media del neutrón aislado
es τn = (880,2 ± 1,0) s.
Decaimiento de piones cargados
En el caso de los piones cargados, el decaimiento está mediado por la interacción débil. El modo
más probable es
π − → µ − + ν µ , π + → µ + + νµ ,
que procede con intercambio de bosones W ± . La vida media de los piones cargados es de τπ± =
2,60 × 10−8 s, varios ordenes de magnitud mayor que los valores tı́picos de decaimientos fuertes y
electromagnéticos.
Decaimiento de leptones
El muón y el tau (y sus antipartı́culas) decaen débilmente en leptones de generación más liviana
con la emisión de un W ± .
µ− → e− + νµ + ν e , µ + → e + + ν µ + νe .
El tiempo de vida media del muón es τµ ∼ 2,2 × 10−6 s, mientras que la del tau, ττ ∼ 2,90 × 10−13
s. El tau, sin embargo, es suficientemente masivo como para decaer en hadrones. Estos decaimientos
son, en conjunto, más probables que los puramente leptónicos.
2.7. Propiedades intrı́nsecas de las partı́culas

2.7.1. El spin
El spin es un momento angular intrı́nseco que tienen todas las partı́culas. Es una propiedad
cuántica que no tiene análogo clásico. El spin se mide en múltiplos de ℏ = h/2π. Las partı́culas
materiales (tanto los quarks como leptones) tienen spin ℏ/2 y son llamadas fermiones (pues tienen
spin semientero). Los fermiones obedecen el principio de exclusión de Pauli: dos fermiones idénticos
no pueden ocupar el mismo estado cuántico.
El spin de todas los bosones de gauge es un múltiplo entero de ℏ: los fotones, gluones, W ± y Z 0
tienen spin ℏ. De existir, los gravitones tendrı́an spin 2ℏ. Las partı́culas con spin entero se llaman
bosones y no obedecen al principio de Exclusion de Pauli.
2.8. El color
Los quarks sienten la fuerza fuerte mientras que los leptones no. Esto quiere decir que hay una
propiedad que tienen los quarks que los leptones no: el color, que es la carga asociada a la fuerza
fuerte. Hay tres tipos de carga de color: rojo (r), verde (g) y azul (b).
Una caracterı́stica fundamental de la fuerza fuerte es que determina que para poder existir como
partı́cula libre, un hadrón no debe tener color neto. Esto se logra con cualquier combinación de
quarks que incluya los tres colores o un color y su anticolor, lo que explica por que los hadrones solo
se observan en combinaciones de quarks de la forma qqq, qqq y qq.
Por ejemplo, en términos de la carga de color de los quarks que lo forman, un protón puede ser
p ≡ ur ub dg o p ≡ ur ug db .
La propiedad anterior, conocida como confinamiento de color, implica que no pueden existir los
quarks libres: están confinados dentro de los hadrones. Sin embargo, dentro de los hadrones, los
quarks pueden considerarse esencialmente partı́culas libres. Esto se relaciona con otra propiedad de
la fuerza fuerte llamada libertad asintótica: su intensidad decrece a medida que decrecen las escalas
16
de distancia (y energı́a) involucradas en el proceso. Cuando en una interacción se intenta extraer
un quark de un hadrón, la densidad de energı́a en el campo fuerte aumenta hasta que supera el
umbral de creación de un par quark-antiquark. Estos se combinan en partı́culas sin color que son las
que se observan como las partı́culas libres producto de la interacción. Este proceso se conoce como
hadronización.
Aunque a fuerza fuerte entre quarks no puede describirse en términos clásicos, se la suele repre-
sentar de forma aproximada a través del siguiente potencial:
4 αs
V (color) = − + kr,
3 r
donde r es la separación entre los quarks, αs ∼ 1 J·m y k = 1,36 × 105 J m−1 .
Una tercera propiedad importante de la fuerza fuerte es que su partı́cula mensajera, el gluón, tiene
carga de color y por lo tanto también siente la fuerza fuerte. Esto significa que, además de interactuar
con los quarks, los gluones interactúan también entre ellos mediante la fuerza fuerte (mediante el
intercambio de gluones). Al tener color los gluones están sujetos al confinamiento de los hadrones;
al intentar separar gluones se producen los mismos efectos de hadronización que en el caso de los
quarks.
En una interacción fuerte un quark puede cambiar su color (pero no su sabor). Como el color debe
conservarse, la diferencia de color debe llevársela un gluón. Los gluones tienen entonces dos colores:
un color y un anticolor. A pesar de que r, g y b y sus anticolores permiten formar 9 combinaciones,
solo 8 permitidas por la teorı́a. No es posible tener gluones con carga de color de tipo rr, gg o bb, es
decir, sin color neto.
3. Aceleración de partı́culas
Una partı́cula cargada que se mueve por una región del espacio donde hay campo electromagnético
experimenta la fuerza de Lorentz dada por
 
⃗
d⃗p ⃗ + ⃗v × B  ,
= e E
dt c
donde ⃗v es la velocidad de la partı́cula y e su carga eléctrica.

El segundo término del lado derecho de la ecuación corresponde a la fuerza ejercida por el campo
magnético, que es perpendicular tanto a la velocidad de la partı́cula como al campo magnético, puesto
que es el resultado de un producto vectorial. Al ser perpendicular a la velocidad, actúa como fuerza
centrı́peta: modifica la dirección del movimiento de la partı́cula manteniendo constante el modulo de
la velocidad. Esto significa que el campo magnético no realizará trabajo sobre la carga y no modifica
entonces su energı́a cinética. La variación en la energı́a se debe solo a la acción del campo eléctrico:
dE ⃗ = e d⃗r · ∇V,
= e⃗v · E ⃗
dt dt
con V el potencial escalar. En el caso unidimensional,
∆E ≃ e∆V.
Ası́, la forma más sencilla de acelerar una carga es someterla a una diferencia de potencial.
17
3.1. Aceleradores artificiales
El primer acelerador de partı́culas moderno fue el ciclotrón de Lawrence, desarrollado entre 1928
y 1931. Un diagrama se muestra en la Figura 3.
El instrumento consta de dos regiones semi-circulares donde un campo magnético deflecta las
partı́culas cargadas. El radio de giro de las partı́culas es
(E/eV)
rg = cm,
300Z (B/G)
donde Z es el número atómico.

Separando las dos regiones semi-circulares, hay una
brecha (gap) donde se aplica una diferencia de potencial.
Las partı́culas se inyectaban en el gap (por ejemplo, co-
locando un material radioactivo) y eran aceleradas por la
diferencia de potencial. Debido a los efectos de campo, las
partı́culas describı́an una semi-circunferencia y volvı́an a
cruzar el gap, donde la diferencia de potencial habı́a sido
invertida durante el tiempo de vuelo. El proceso se repetı́a
una y otra vez, generándose más y más energı́a, hasta que
el giroradio de la partı́cula se hacia igual al máximo espacio
disponible.
Para lograr mayores energı́as se desarrollaron los acele-
radores tipo sincrotrones. En estos, las partı́culas se mue-
ven a lo largo de una trayectoria circular fija y se utilizan
campos magnéticos variables. Las diferencias de potencial
se aplican en diferentes tramos del circuito. Se utilizan ima-
Figura 3: Esquema de un acelerador ci-
nes cuadrupolares para mantener el haz de partı́culas bien
clotrón.
colimado e imanes superconductores para desviar el haz y
mantenerlo dentro del anillo de circulación.
La energı́a máxima que puede obtener una partı́cula en un acelerador sincrotrón esta determinada
por
la intensidad máxima de los campos magnéticos,
el tamaño del anillo y
las pérdidas de energı́a de las partı́culas (por radiación sincrotrón).
Un elemento esencial de todo experimento con aceleradores de partı́culas, independiente de la

forma del acelerador y del hecho de que se colisionen haces de partı́culas entre sı́ o con blancos fijos,
lo constituyen los detectores. Estos pueden ser de muchos tipos, pero esencialmente todos se basan
en el proceso de ionización. La forma en la cual los iones que se forman son usados para rastrear el
movimiento de las partı́culas viene dado por el tipo especifico de detector.
3.2. Rayos cósmicos

Los rayos cósmicos son partı́culas energéticas que llegan a la Tierra desde el espacio exterior
(fuera de la atmósfera).
18
Las partı́culas primarias que forman los rayos cósmicos son protones (86 %), partı́culas α (núcleos
de He, 11 %), electrones (2 %) y núcleos más pesados (1 %). Hay, además, pequeñas proporciones de
positrones y antiprotones que se creen de origen secundario (esto es, que se generan debido a las
interacciónes de los rayos cósmicos primarios con el medio interestelar).
El espectro de energı́a de los
rayos cósmicos se extiende desde
energı́as del orden de 106 eV (rango
dominado por contribuciones locales
del Sistema Solar) hasta energı́as por
encima de 1020 . En el rango más al-
to de energı́a, la composición de los
rayos cósmicos es incierta.
A energı́as ∼ 1011−12 eV el flu-
jo de rayos cósmicos en la Tierra
es de ∼ 1 partı́cula por m2 s−1 . A
energı́as ∼ 1015−16 eV, donde el es-
pectro cambia de N (E) ∝ E −2,7 a
N (E) ∝ E −3 , el flujo es de 1 partı́cu-
la por m2 año−1 . A energı́as muy al-
tas, arriba de 1018 eV, donde el es-
pectro parece volver a endurecerse,
el flujo es de ∼ 1 por km2 año−1 .
Vemos el espectro en la Figura 4.
Entre unos pocos GeV y la llama-
da rodilla (knee) del espectro, a unos
1014−15 eV, el espectro está bien des-
cripto por
N (E)dE ∝ E −2,7 dE.
A energı́as mayores a ∼ 30 GeV,

donde los efectos de modulación de
los campos magnéticos solar y te-
Figura 4: Espectro observado de rayos cósmicos. rrestre dejan de ser importantes, la
radiación cósmica parece ser comple-
tamente isotrópica. Esto es de esperar ya que, independientemente de la naturaleza de las fuentes,
el campo magnético galáctico destruye la anisotropı́as de partı́culas cargadas. Solo en el caso de las
energı́as más altas podrı́a esperarse observar anisotropı́a en las direcciónes de arribo de los rayos
cósmicos desde fuentes cercanas. Las partı́culas más energéticas deben ser de origen extragaláctico,
ya que su giroradio no puede ser contenido en la Galaxia.
Los efectos de los campos magnéticos hacen que una astronomı́a galáctica de rayos cósmicos sea
imposible, ya que en las partı́culas detectadas no guardan memoria de las fuentes en su dirección de
arribo. Solo en el caso de fuentes extragalácticas cercanas y para partı́culas con energı́as por arriba
de 1019 eV tal astronomı́a es factible.
La densidad de energı́a de rayos cósmicos en regiones lejos de la influencia solar es
eV
ωRC ∼ 1 .
cm3
19
Esta densidad es comparable a la de la luz estelar (∼ 0,6 eV/cm3 ), a la del fondo cósmico de
radiación (∼ 0,26 eV/cm3 ) y a la de campo magnético galáctico (∼ 0,6 eV/cm3 , para B ∼ 3 µG).
En base a la densidad medida de rayos cósmicos se puede calcular la potencia inyectada por los
mismos en la Galaxia. Consideremos que tiene un disco de radio RG = 15 pc y espesor hG = 200 pc.
Su volumen será
VG = πR2 hG ∼ 4 × 1066 cm3 .
El tiempo de residencia de los rayos cósmicos en el disco vendrá determinado por la escala temporal
de difusión de las partı́culas
3r2
td ∼ ,
D
donde r es la distancia a recorrer por las partı́culas y D es el coeficiente de difusión. Para el caso de
la Galaxia, r = hG y D ∼ 1028 cm2 s−1 , por lo que resulta td ∼ 1014 s. Luego,
VG ωRC
WRC = ∼ 4,1 × 1040 erg s−1
td
Esta es la potencia total de los rayos cósmicos en la Galaxia. Cualquier acelerador o conjunto de
aceleradores que los produzca debe satisfacer este supuesto energético.
A fin de postular posibles aceleradores galácticos de rayos cósmicos, el primer paso es explorar
si existen sistemas astrofı́sicos capaces de generar grandes diferencias de potencial que permitan
acelerar en forma electrostática partı́culas cargadas.
3.2.1. El criterio de Hillas

Si una partı́cula escapa de la región donde se estaba acelerando, no podrá obtener más energı́a.
Esta situación impone un lı́mite en la energı́a máxima
Emax = qBR,
donde q es la carga eléctrica de la partı́cula acelerada, B es el campo magnético y R, el tamaño del
acelerador. Esta ecuación se obtiene exigiendo que el radio de Larmor de la partı́cula, Rg = E/(qB),
no exceda el tamaño de la región de aceleración. Este es un criterio geométrico general conocido
como criterio de Hillas y es útil para seleccionar posibles sitios de aceleración.
Si la carga de la partı́cula es Ze, la energı́a máxi-
ma es
Emax = qZBR.
El diagrama de Hillas (Figura 5) es un gráfi-
co que muestra la relación entre el tamaño de una
fuente y el campo magnético para los posibles esce-
narios de aceleración de rayos cósmicos. Las lı́neas
diagonales indican diferentes energı́as máximas de
las partı́culas. Estas energı́as son posibles para una
cierta fuente si la lı́nea pasa por la región del dia-
grama ocupada por la fuente. Se puede ver que
las fuentes más compactas requieren grandes cam-
pos magnéticos para lograr la aceleración de rayos
cósmicos más energéticos. En cambio, los objetos
extendidos pueden acomodar a estas partı́culas en
principio con campos magnéticos mucho más débi-
les.
Figura 5
20
3.2.2. Púlsares
Los púlsares son estrellas de neutrones: estrellas colapsadas donde la presión de la gravedad
es sostenida por la presión de degeneración de los nucleones. El tamaño tı́pico de las estrellas de
neutrones es R∗ = 106 cm y su masa M + 1,4M⊙ . Esto hace que su densidad sea
m∗ g
ρ∗ ∼ ∼ 6,6 × 1014 3 .
(4/3)πR∗
3 cm
Al colapsar, las estrellas arrastran su campo magnético junto con la materia, por lo que las
estrellas de neutrones resultantes están fuertemente magnetizadas y en rotación rápida (debido a la
conservación del momento angular). El campo magnético resultante es, en primera aproximación,
dipolar.
El campo magnético superficial tı́pico B ≈ 1012 G. Sean Ω la velocidad angular de la estrella y

P = 2π/Ω su periodo de rotación. A una distancia
P

RL = c/Ω ∼ 4,8 × 9
cm,
1s
las lı́neas de campo magnético dipolar deberı́an moverse más
rápido que la luz, y con ellas las partı́culas que por ellas se mue-
van. De aquı́ que a partir de RL las lı́neas se abren cerrándose
en el infinito. Las partı́culas que se mueven sobre estas lı́neas
de campo pueden entonces escapar del sistema.
Si una esfera magnetizada rota, las cargas sobre ella experi-
mentarı́an una fuerza
e ⃗
F⃗ = Ω × ⃗r × B.⃗
c
Esta fuerza hará que las cargas de signo opuesto se sepa-
ren originando un campo eléctrico tal que se anule la fuerza de
Lorentz. En primera aproximación
⃗ = −1 Ω

E ⃗ × ⃗r × B.
⃗
c
Como la fuerza de Lorentz no es la única fuerza que actúa,
la zona inmediatamente cercana al púlsar presenta una gran complejidad.
21
La rotación de la esfera separa la carga y esta no se modifica a menos que cambie Ω. Por lo tanto,
el potencial eléctrico es estático. El potencial electrostático generado en la superficie (r = R∗ ) de la
estrella de neutrones será
ΩBR∗2
−1
B R∗ P

V ∼ ER∗2 = ∼ 6 × 1016 V.
c 1012 G 106 cm 1s
Vemos que las partı́culas arrancadas de la superficie y que se mueven por las lı́neas de campo
pueden ser aceleradas hasta grandes energı́as.
Sin embargo, como las partı́culas que escapan del sistema lo hacen a través de las lı́neas abiertas,
el flujo generado quedara determinado por r = a, donde a es el radio donde se originan las lı́neas
abiertas. Sabiendo que las lı́neas de un dipolo quedan definidas por sen2 θ/r, debe ser que a ∼ =
R∗ sen θ0 . Partiendo de
!1/2 s
sen2 θ0 ∼ Ω R∗3 Ω
3 −2
R∗3 B R∗ P

= =⇒ a ∼
= ∼
= =⇒ V ≈ 1013 V.
R∗ c RL c 1012 G 106 cm 1s
Vemos que los púlsares pueden acelerar partı́culas hasta energı́as de ∼ 1013 eV = 10 TeV. Si el
púlsar rota muy rápidamente (P ∼ 0,1 − 0,01), Emax ∼ 1015−16 eV. El lı́mite efectivo es mejor que
este valor ua que las partı́culas sufrirán pérdidas radiativas durante el proceso de aceleración.
Para saber si los púlsares contribuyen al grueso de los rayos cósmicos galácticos, necesitamos saber
cuantos púlsares con periodos P < 1s hay en la galaxia y cuanta energı́a en partı́culas relativistas
deposita un púlsar en el medio interestelar.
La tasa de nacimiento de púlsares es de aproximadamente 1 cada 80 años. Como el campo
magnético decae con el tiempo, la vida media de ellos es ∼ 10 Myr. Entonces el número de púlsares
es
107
N∼ ∼ 105 .
80
Por otro lado, tenemos el cociente σ entre las densidad de energı́a asociadas es
WPoynting
σ= , (3.1)
Wviento
que depende de la distancia al púlsar y, en cierta medida, de la dirección. Se estima que vale σ ∼ 1000
cerca de la superficie del púlsar y cae hasta un valor de σ ∼ 0,1 en unos cuantos radios estelares
3 −2
B R∗ P

σ ∼ 10 .
1012 G 106 cm 1s
Para P = 0,1 s tenemos σ ∼ 1000 y por lo tanto, para una potencia tı́pica en un púlsar, se obtiene
que la potencia inyectada en rayos cósmicos por un púlsar es
Ė 104
WRC ∼ ≈ 6 × 1031 erg s−1 × 3 ∼ 6 × 1032 erg s−1 .
σ 10
Si hay del orden de 105 púlsares en la Galaxia
púlsares
WRC ∼ 6 × 1037 erg s−1 .
Este valor es mucho menor que el necesario para explicar el origen de los rayos cósmicos galácticos.
Debe existir entonces otro tipo de fuente astrofı́sica en Galaxia capaz de acelerar rayos cósmicos.
22
3.2.3. Remanentes de supernova
Cuando ocurre una explosión de supernova, el material que formaba la estrella es eyectado y el
medio interestelar es perturbado en la región donde ocurrió la explosión.
Una explosión de supernova tı́pica libera una energı́a E ∼ 1051 erg. Esto se produce en una escala
temporal muy corta, mucho más corta que cualquier otra escala temporal involucrada. La velocidad
de la materia eyectada por la explosión veyec , se relaciona con la energı́a a través de la expresión
!1/2 !−1/2
E M
veyec ≈ 104 km s−1
1051 erg M⊙
!1/2 !−1/2
E M
≈ 10 −2
pc yr−1
1051 erg M⊙
La velocidad del material eyectado tı́picamente es mayor que la velocidad del sonido en el medio
interestelar, por lo que origina ondas de choque al propagarse. Una onda de choque es una
perturbación que produce una discontinuidad en los parámetros que caracterizan un
medio. El espacio queda dividido por la onda de choque en dos regiones: la región chocada y la
no chocada. Para que se genere un choque es necesario que la perturbación tenga una velocidad
supersónica, puesto que de lo contrario la perturbación se disiparı́a en forma de ondas sonoras. Sean
Vch la velocidad de la onda de choque y Cs la velocidad del sonido en el medio. El número de Mach
se define como M ≡ Vch /Cs . En el caso de las supernovas, la onda de choque se mueve hacia afuera
del sitio de la explosión, viajando delante de la superficie de separación entre el medio interestelar y
el material eyectado.
Figura 6: Choque plano visto en el sistema de referencia del medio 1 (izquierda) y en el sistema de
referencia de la onda de choque (derecha). Notar que V2 = Uch − U2 < V1 = Uch .
El medio interestelar puede modelarse a través de una ecuación de estado politrópica tal que
P V γ = cte, T V γ−1 = cte T P (1−γ)/γ = cte,
donde γ = CP /CV es el ı́ndice adiabático del gas. Los calores especı́ficos a presión y volumen constante
son ! ! ! ! !
dQ dU dV V dQ dU
CP = = +P , y CV = = ,
dT P dT P dT P dT V dT V
respectivamente. Para un gas monoatómico, γ = 5/3.
23
Las propiedades del medio perturbado se relacionan con las del medio sin perturbar por medio
de las relaciones de Rankine-Hugoniot (también llamadas relaciones de salto)
ρ1 V1 = ρ2 V2 ,
+ P1 = ρ2 V22 + P2 ,
ρ1 V12
γP1 1 γP2 1
+ V12 = + V22 ,
(γ − 1)ρ1 2 (γ − 1)ρ2 2
las cuales se deducen a partir de las ecuaciones de conservación de la energı́a, momento y masa.
En el caso de que una onda de choque fuerte (M1 ≫ 1) y no relativista (Vch < 0,1c) que se
propaga en un gas ideal no conductor, en estado estacionario y con campo magnético despreciable,
la relación entre las densidades a ambos lados del frente de choque resulta
ρ2 γ+1
ξ= = ,
ρ1 γ−1
donde los subı́ndices 1 y 2 designan a las propiedades del medio no chocado y chocado, respectiva-
mente. Para un gas monoatómico γ = 5/3 y entonces
ρ2
= 4.
ρ1
Las velocidades del material se relacionan como
V2 1
= ξ −1 = ,
V1 4
mientras que el cociente de temperaturas es
T2 2γ(γ − 1)M12 5
= = M12 .
T1 (γ + 1) 2 16
Como M1 puede ser muy grande, la región chocada puede estar muy caliente respecto a la no
chocada.
En el caso de incorporar un campo magnético perpendicular al frente de choque, se encuentra
una relación adicional dada por
B2 ρ2
= = 4.
B1 ρ1
Estos resultados muestran que la temperatura, presión, densidad y campo magnético del medio
chocado aumentan, mientras que la velocidad del fluido decrece. En consecuencia, el material detrás
del choque ya no es supersónico.
Etapas de un remanente de supernova
Debido a que la inyección de una explosión de SN es puntual en el tiempo, la onda de choque
producida se diluye a medida que se propaga en el medio interestelar. Durante la propagación, las
propiedades del choque van cambiando con el tiempo. Podemos separar tres etapas principales
Fase inicial o de expansión libre: en una etapa inicial, luego de la explosión de SN, el material
eyectado se mueve a una velocidad uniforme
r ∝ t,
24
donde r es el radio de la onda de choque. Esta fase termina cuando la masa del medio barrido
por la onda de choque es iguala la masa eyectada en la explosión:
4 −1/3
3
ρMIE rI,f = Meyec =⇒ rI,f ∝ ρMIE Meyec
1/3
,
3
siendo ρMIE la densidad del medio interestelar. Considerando una densidad de partı́culas del
medio interestelar nH = ρMIE /mp , puede escribirse rI,f en unidades convenientes como
!1/3 −1/3
Meyec nH
rI,f ≃ 2 pc.
M⊙ 1 cm−3
Esta distancia corresponde a un tiempo

!5/6 !−1/2 −1/3
rI,f Meeyc E nH
tI, f = ≃ 200 yr.
veyec M⊙ 1051 erg 1 cm−3
Fase Sedov o adiabática: la siguiente fase de la evolución del remanente de supernova es la

fase de Sedov. Aquı́, la energı́a es aproximadamente constante ya que las pérdidas radiativas se
pueden despreciar. La evolución es adiabática, es decir,
1 4π

E≈ ρMIE r3 veyec
2
∼ ρMIE r3 ṙ2 .
2 3
Luego,
!1/5
3/2 dr −1/2 E
r ≃ ρMIE E 1/2
=⇒ r ∝ t2/5 .
dt ρMIE
Asumiendo equipartición entre la energı́a interna y la cinética

3 1 2
nH kT ∼ Meyec veyec ,
2 2
se obtiene que la temperatura del material eyectado por la SN y calentado por la onda de
choque es
T ∝ v 2 ∝ r−3 E51 n−1
H .
En unidades convenientes
!−3 !−6/5
r 2/5 −2/5 t
T ≃ 6 × 10 8 −1
E51 nH K = 10 6
E51 nH K.
2 pc 3 × 104 yr
Cuando las pérdidas radiativas empiezan a afectar la dinámica del remanente, este abandona
la fase de Sedov. Invirtiendo la ecuación para la temperatura obtenemos la edad de la fase de
Sedov
−5/6 1/3 −1/3
tSedov ≃ 3 × 104 T6 E51 nH yr,
donde T6 = T /(106 K).
25
Fase radiativa: cuando la escala temporal de enfriamiento radiativo del gas, tenf , se hace menor
que tSedov , el remanente entra en la llamada fase radiativa. Dado que
3/2
T
tenf =≃ 4 × 10 6 yr,
4
nH
esta condición se cumple cuando
ṙ ≤ 200km s−1 (E51 n2H )1/14 .
Como la dependencia en E y nH es débil, esto ocurre para velocidades v ∼ 200 km s−1 .

En la nueva fase, el interior del remanente está formado por gas caliente rodeado por una cascara
de gas frı́o. A medida que la cascara avanza, va barriendo el medio e incrementando su masa.
De aquı́ que, a pesar de ir desacelerándose, su momento radial permanezca aproximadamente
constante
d 4π

ρr ṙ ≈ 0.
3
dt 3
Si la cascara se formo en t0 con radio r = r0 y velocidad v = v0 , entonces
4π 3 4π 3
ρr ṙ = ρr v0 =⇒
3 3 0

v0 1/4
r = r0 1 + 4 t − t0 ,
r0

v0 3/4
ṙ = v0 1 + 4 t − t0 ,
r0
Para t grandes, r ∝ t1/4 , luego

!−3/4
t
ṙ ∝ t −3/4
≃ 200 km s −1
3 × 104 yr
Cuando t ∼ 105 yr, la velocidad cae a un valor ∼ 10 km s−1 y el remanente comienza a disiparse.
Las explosiones de supernova inyectan en el medio interestelar cantidades enormes de energı́a,

del orden de ∼ 1051 . Es razonable preguntarse, entonces, su juegan algún papel en la producción de
rayos cósmicos. La tasa de explosiones de supernova es de 1 cada 50 años en la Galaxia. La potencia
que inyectan en el medio interestelar es
1051 erg
WSN = ∼ 6,3 × 1041 erg s−1 .
50 × 3,15 × 107 s
Entonces, si las supernovas pudieran transformar menos del 10 % de su energı́a en rayos cósmicos,
serı́an capaces de explicar el origen de estos.
El hecho de que se observe emisión de origen no térmico en remanentes de supernovas implica
que en ellos hay partı́culas relativistas. La detección de rayos X de origen sincrotrónico indica la
presencia de electrones con energı́as por arriba del TeV y la detección de radiación gamma podrı́a
señalar la presencia de hadrones con energı́as similares. El mecanismo que acelera estas partı́culas
parecerı́a ser difusivo (lo que lleva fácilmente a un espectro del tipo ley de potencias), pero existen
26
también otros mecanismos posibles de aceleración. Lo que todos tienen en común es que la presencia
de un campo magnético juega un papel crucial en la aceleración de las partı́culas.
La relevancia de un proceso de aceleración de partı́culas esta directamente relacionada con la efi-
ciencia de aceleración del mismo, ηacel . Para determinar que tan rápidamente se aceleran las partı́culas
resulta conveniente definir un tiempo caracterı́stico de aceleración parametrizado como
qcB
acel = ηacel
t−1 ,
E
con q y E la carga y energı́a de la partı́cula, respectivamente. La máxima eficiencia posible de
aceleración corresponde a ηacel = 1, mientras que valores ηacel ≪ 1 implican aceleración poco eficiente.
3.3. Mecanismo de aceleración difusivo

Se arroja una pelota perfectamente elástica contra un muro. Si la velocidad de la pelota es −V⃗
en el sistema de referencia del laboratorio, rebotará con velocidad +V⃗ .
Si ahora la pared se mueve con una velocidad ⃗u, luego del choque la velocidad de la pelota será
V⃗ + 2⃗u (la velocidad no cambia en el sistema de referencia de la pared).
Si en vez de una se tienen dos paredes, una moviéndose hacia la otra, la partı́cula ganara velocidad
con cada choque hasta que
1. las paredes se detengan
2. o la pelota se haga tan energética que atraviese una de las paredes - en cuyo caso se dice que
la pared se ha hecho transparente a la pelota.
Reemplacemos ahora las paredes por nubes en movimiento que contienen campos magnéticos
turbulentos. La turbulencia puede actuar como un espejo magnético creando regiones con mayor
densidad de energı́a magnética donde las partı́culas cargadas son deflectadas. Esta situación es fı́si-
camente improbable y no es posible generar de esta manera los rayos cósmicos galácticos.
Figura 7: Colisión elástica de una partı́cula con una pared en movimiento. Se muestra la velocidad
de la partı́cula antes (izquierda) y después (derecha) de la colisión.
El mismo mecanismo básico puede operar, sin embargo, en remanentes de supernova o en otros
sistemas astrofı́sicos con choques. Consideremos una onda de choque que genera una supernova como
la que ya conocemos. En ambas regiones puede establecerse turbulencia magnética: en la región
chocada por la compresión del fluido y el desarrollo de inestabilidades, y en la región no chocada por
los efectos de los propios rayos cósmicos que atraviesan el frente de choque y perturban el plasma
generando inestabilidades.
27
Supongamos que el choque se mueve con ve-
locidad Vs en la dirección del eje x. El plasma en
la región chocada se mueve con velocidad Vp . La
velocidad del choque y la del gas detrás de él se
relacionan por el factor de compresión ξ como
Vs
= ξ.
Vp
Para un choque fuerte no relativista, ξ = 4.

Usando transformaciones relativistas puede es-
tablecerse que la variación neta de energı́a de una partı́cula con energı́a inicial E que completa un
ciclo yendo de la región chocada a la no chocada y de vuelta a la chocada es
4 ξ−1
!
⟨∆E⟩ Vs
≈ .
E 3 ξ c
Para obtener el espectro de energı́a tenemos que hallar la probabilidad de que una dada partı́cula
cruce el choque un número arbitrario de veces. Para ello consideremos el proceso en el sistema de
referencia del choque.
El flujo neto de partı́culas que se pierden en la región
chocada es
Vs
rper = nRC cm−2 s−1
ξ
donde nRC es la densidad de partı́culas relativistas.
En la región no chocada, un rayo cósmico con velocidad
v (en el sistema de referencia del laboratorio) que forma un
ángulo θ, se acercara al choque con velocidad Vs + v cos θ
en el sistema de referencia del choque. A fin de que pueda
cruzar el choque,
Vs
cos θ > − .
v
Entonces, suponiendo que los rayos cósmicos en la región no chocada son isotrópicos, la razón a
la cual las partı́culas cruzan desde la región no chocada a la chocada es
ˆ 1
1 v
rcruce = nRC Vs + v cos θ) · 2πd(cos θ) =⇒ rcruce ≈ nRC cm−2 s−1 .
4π − Vvs 4
La probabilidad de cruzar y luego escapar en la región chocada es

rper Vs
Pescape = ≈4 .
rcruce ξv
La probabilidad de que la partı́cula sea deflectada en la región chocada y vuelva al choque es
Pregreso = 1 − Pescape .
Por lo tanto, la probabilidad de que cruce el choque n o más veces es

n
P(cruce>n) = 1 − Pescape .
28
Como el incremento fraccional de energı́a por cruce es ∆E/E, la energı́a luego de cruzar n veces
es !n
∆E
E = E0 1+ ,
E
donde E0 es la energı́a inicial. Si tomamos logaritmo a ambos lados,

∆E ln E/E0
!
E

ln = n ln 1 + , donde n = .
E0 E ln 1 + (∆E/E)
Luego, el espectro de partı́culas aceleradas y con energı́as mayores a E será

n
J(> E) ∝ 1 − Pescape ,
ln(E/E0 )/ ln(1+(∆E/E))
J(> E) = k 1 − Pescape , k = cte,

ln E/E0
ln J(> E) = k ′ + ln 1 − Pescape , k ′ = ln(k) = cte
ln 1 + (∆E/E)
Finalmente,
ln J(> E) = k ′′ − (Γ − 1) ln(E),
donde
ln 1 − Pescape ln 1 − 4Vs
ξv
Γ=1− =1− .
ln 1 + (∆E/E) ln 1 + 4 (ξ−1)Vs
3 ξc
Haciendo un desarrollo en serie de Taylor para x ∼ 0, encontramos que

ξ+2
Γ= .
ξ−1
Por lo tanto, el espectro diferencial de rayos cósmicos que emerja de la fuente luego de n cruces
será
J(E) ∝ E −Γ .
Para una onda de choque fuerte

ξ = 4 =⇒ Γ = 2.
Si bien el espectro observado de rayos cósmicos es de la forma J(E) ∝ E −2,7 , la diferencia en el
ı́ndice espectral se explica por efectos de difusión de las partı́culas en el medio interestelar.
Para calcular la energı́a máxima que pueden alcanzar las partı́culas, calculamos la tasa de ganancia
de energı́a de las partı́culas e igualamos la tasa de pérdidas radiativas. Ahı́ podemos despejar la
energı́a máxima, siempre y cuando esta satisfaga que el giroradio de las partı́culas con esa energı́a
sea menor que el tamaño de la región de aceleración.
La tasa de aceleración es
1 dE 4 (ξ − 1)Vs −1
t−1
acel ≈ ≈ tciclo .
E dt 3 ξc
Aquı́, tciclo ∼ D/Vs2 depende del coeficiente de difusión D del medio. Este coeficiente se supone
dependiente de la energı́a D = D(E) ∝ E δ , con δ > 0. El coeficiente de difusión a lo largo del campo
29
magnético es un cierto número de veces el valor del coeficiente de difusión mı́nimo, conocido como
coeficiente de difusión de Bohm:
1 E
DBohm = rg c, rg ≈ ,
3 qB
donde aproximamos el giroradio rg considerando que las partı́culas son relativistas (v ∼ c).
Definimos choque paralelo si la normal al frente de choque es paralela a B.
⃗ Denominamos choque
perpendicular si Vs ⊥ B. Entonces, el coeficiente de difusión paralelo es
⃗ ⃗
D∥ = ζDBohm .
Para valores tı́picos de Vs = 0,1c y ζ = 10, obtenemos
dE
≈ 1,5 × 10−4 ec B.
dt
acel
Por otro lado, el coeficiente de difusión perpendicular es
D∥ dE
D⊥ ≈ =⇒ ≈ 4 × 10−2 ec B.
c+ζ 2 dt
acel
Las ondas de choque de los remanentes de supernova permanecen fuertes como para acelerar
partı́culas durante unos ∼ 103 años. Las energı́as máximas alcanzables (considerando B ∼ 10−5 G y
pérdidas sincrotrón) son
SN
Emax ∼ 1014 Z eV para choques paralelos,
SN
Emax ∼ 1016 Z eV para choques perpendiculares.
Para el caso de que el choque sea oblicuo, con un ángulo θ entre la normal al choque y el campo
magnético, el coeficiente de difusión es
D = D∥ cos2 θ + D⊥ sen2 θ.
En esta situación el valor de la energı́a máxima esta comprendido entre los valores máximos dados
anteriormente.
Efectos adicionales importantes a tener en cuenta en diversas situaciones astrofı́sicas son
1. La modificación del choque por efecto de la presión de los rayos cósmicos que están siendo
acelerados, lo que introduce no linealidades.
2. Ondas de choque relativistas. En este caso el indica adiabático del gas es γ = 4/3, lo que
lleva a ξ = 7 y por lo tanto a espectros con Γ ∼ 1,5 si se tienen en cuenta efectos de anisotropı́a.
En choques relativistas se tiene tı́picamente 1,5 ≤ Γ ≤ 2.
3. Efectos termodinámicos producidos por la radiación que escapa de la región chocada, afec-
tando los centros de dispersión magnéticos.
4. Inestabilidades introducidas por los propios rayos cósmicos que pueden llevar a una magni-
ficaron de los campos magnéticos (inestabilidades de Bell).
30
Hay una gran variedad de fuentes astrofı́sicas en las que se forman ondas de choque y pueden
acelerarse partı́culas cargadas hasta velocidades relativistas. Estas fuentes incluyen:
jets en AGNs y microcuasares,
ondas de choque producidas por eruptores de rayos gamma,
lóbulos y manchas calientes (hot spots) en radio galaxias,
asociaciones de estrellas masivas,
estrellas tempranas con vientos poderosos,
cúmulos de galaxias donde colisiones pueden formar ondas de choque,
sistemas binarios con colisión de vientos estelares y
choques de nubes de alta velocidad con el medio interestelar.
3.4. Magnetohidrodinámica
La teorı́a de la magnetohidrodinámica (MHD) es una posible aproximación al problema de un
fluido en presencia de un campo electromagnético. Para que esta aproximación sea valida, se deben
cumplir dos condiciones básicas
el fluido debe ser eléctricamente casi neutro (es decir, con una densidad de carga muy pequeño)
y
debe formar un plasma (estar completamente ionizado).
En MHD no relativista, el campo relevante es el campo magnético B.

⃗ De hecho, si la velocidad
del flujo es v ≪ c, en todas partes, se puede demostrar que |E| ⃗ Esto también significa que la
⃗ ≪ |B|.
corriente puede ignorarse, de modo que
⃗ = 4π J,
⃗ ×B
∇ ⃗ (3.2)
c
donde J⃗ es el vector de corriente eléctrica. Teniendo en cuenta, además, que la densidad de carga es
muy pequeña, la expresión de la fuerza de Lorentz se simplifica a
J⃗ × B
⃗ 1 ⃗
F⃗L = == ⃗ × B.
∇×B ⃗
c 4π
Otra de las hipótesis de la MHD es que el fluido satisface la ley de Ohm en la forma simple
!
⃗ + ⃗v × B
J⃗ = σ E ⃗ , (3.3)
c
donde σ es la conductividad eléctrica del plasma.

A estas ecuaciones debemos agregarles el resto de las ecuaciones de Maxwell:
⃗ = − 1 ∂B ,
⃗
⃗ = 0,
∇·B ⃗ ×E
∇ ⃗ = 4πρe .
⃗ ·E
∇
c ∂t
31
La densidad de carga ρe no es un parámetro relevante y si es necesario se puede calcular a partir
de la ley de Gauss una vez que encontremos E.⃗
Combinando la ley de Ampére-Maxwell (3.2), la ley de Ohm (3.3) y la ley de Faraday, obtenemos
la ecuación central de la MHD, llamada ecuación de inducción para el campo magnético
⃗
∂B 2
⃗ + c ∇2 B.

=∇
⃗ × ⃗v × B ⃗
∂t 4πσ
La ecuación de inducción establece que el campo magnético en un punto dado en el espacio varı́a
en el tiempo porque es advectado con el flujo (primer término en el lado derecho) y porque se difunde
(segundo término). Una suposición que se suele hacer en las aplicaciones astrofı́sicas de la MHD es
que la conductividad del plasma es muy grande. Entonces, el termino difusivo puede despreciarse.
Esa aproximación se conoce como MHD ideal. En este caso, las lineas de campo magnético están
congeladas en el fluido: si el fluido se mueve la lı́nea lo acompaña conservando su topologı́a.
El conjunto de ecuaciones de MHD debe completarse con las ecuaciones para la conservación de
la masa y el momento (incluida la fuerza de Lorentz, pero no la viscosidad):
∂ρ ⃗
+ ∇ · ρ⃗v = 0,
∂t
∂⃗v ⃗ ⃗ + 1 ∇

ρ + ⃗v · ∇⃗v = −∇P
⃗ − ρ∇Φ ⃗ ×B
⃗ × B,
⃗
∂t 4π
más una ecuación para la conservación de energı́a y una ecuación de estado. Aquı́, como de costumbre,
P es la presión, ρ es la densidad de masa y Φ es el potencial gravitacional.
La transformación de energı́a magnética requiere romper la topologı́a de las lı́neas: MHD resistiva.
En un fluido de conductividad imperfecta, el campo magnético puede moverse a través del fluido
siguiendo una ley de difusión con la resistividad del plasma sirviendo como constante de difusión.
Además, las lı́neas de campo magnético pueden moverse por el fluido transfiriendo energı́a al mismo.
3.5. Reconexión magnética

La energı́a almacenada en campos magnéticos puede liberarse súbitamente en eventos de reco-
nexión magnética, en los que se produce un campo en la topologı́a de las lineas de campo. Esta
energı́a es transferida a las partı́culas del plasma en escalas de tiempo cortas y puede dar lugar a
una población de partı́culas no térmicas. Para que este proceso puede ocurrir son necesarios:
Un plasma con una región de conductividad finita. Caso contrario, no es posible que la topologı́a
del campo magnético se rompa, pues las lı́neas están congeladas.
Flujos convergentes de material. Si las lı́neas de fluido de los plasmas se desplazan en una
dirección convergente, las lı́neas de campo magnético también se acercan entre ellas.
Regiones con polaridades de campo magnético opuestas. Cuando las lı́neas de campo magnético
paralelas y de polaridad opuesta se acercan a una corta distancia, las corrientes de plasma no
pueden sostenerlas si la conductividad es finita y las lı́neas reconectan.
Supongamos un campo magnético en el plano xy y un campo eléctrico en la dirección del eje z.
En la denominada zona neutral tipo-X, el campo magnético es débil y las partı́culas cargadas son
aceleradas por un campo eléctrico inducido en la dirección z, prácticamente uniforme alrededor de
la región tipo-X. Si los flujos convergentes tienen su movimiento en la dirección del eje y, entonces
el plasma luego de la reconexión es eyectado en la dirección del eje x debido a la tensión magnética.
32
La velocidad con el que el plasma es acelerado
puede ser del orden de la velocidad del Alfvén del
plasma:
B
Va = √ ,
ρµ0
donde B es la intensidad del campo magnético, ρ
es la densidad y µ0 , la permeabilidad magnética
del plasma previo a la reconexión.
La eficiencia de aceleración de este proceso depende de la tasa de reconexión del campo magnético.
Si la misma es alta, la reconexión es rápida y el proceso es eficiente. No obstante, en el modelo
propuesto, la longitud de la región de reconexión es grande y puede demostrarse que en ese caso la
tasa de aceleración se vuelve lenta. Por este motivo hace falta agregar un ingrediente más para que
la reconexión magnética sea un proceso viable de aceleración de rayos cósmicos: la turbulencia.
3.5.1. Reconexión magnética turbulenta

Si el fluido es turbulento, entonces la reconexión magnética puede ocurrir como un gran número
de reconexiónes de pequeña escala. La reconexión sobre estas escalas pequeñas determinan la tasa de
reconexión local, mientras que la tasa de reconexión global es mucho mayor debido a las múltiples
reconexiónes.
En este escenario, las partı́culas relativistas son aceleradas en múltiples islas magnéticas. La
eficiencia de aceleración resulta del orden
vrec 2

ηacel ∼ 0,3 ,
c
donde vrec ∼ VA es la velocidad de reconexión. Dependiendo de las condiciones en el plasma, el ı́ndice
espectral del espectro de partı́culas puede ser p < 2,5.
Vale destacar que existen muchas fuentes astrofı́sicas en las que la reconexión magnética puede
ser un fenómeno importante. Algunas de ellas son
El Sol (los procesos de reconexión se originan en la corona solar debido al movimiento del
plasma subyacente).
Sistemas con discos de acreción tales como MQs y los flujos magnéticos en AGNs.
Estrellas altamente magnetizadas, como las T-Tauri.
Magnetares
Fast radio bursts y GRB.
Nebulosas creadas por púlsares, como la del Cangrejo.
4. Ecuación de transporte
4.1. Difusión y convección
La difusión consiste en el movimiento propio, caótico, de las partı́culas (en este caso, rayos cósmi-
cos) respecto a un medio. Este movimiento es estocástico debido a las múltiples interacciónes que un
rayo cósmico sufre con los átomos del medio y/o las irregularidades magnéticas.
33
El movimiento de las partı́culas puede modelar-
se como una caminata al azar, por lo que, en prome-
dio, el desplazamiento ∆r que una partı́cula realiza
al cabo de un intervalo de tiempo ∆t por difusión es
nulo, es decir, ⟨∆r⟩ = 0. No obstante, se tiene que
⟨∆r2 ⟩ = 6D∆t ̸= 0. El coeficiente de difusión D es
una medida de la facilidad con la que la partı́cula
se propaga en ese medio: cuanto mayor es D, más
lejos difunde en el mismo intervalo de tiempo.
Otro proceso de transporte de partı́culas en un
medio es la convección (o advección). En este ca-
Figura 8: Dos partı́culas con distinto camino
so las partı́culas son arrastradas por el movimiento
libre medio l atravesando un gas con tamaño li-
macroscópico del medio, como granos de arena en el
neal caracterı́stico R. Para la partı́cula de arriba
viento. Lo que ancla a los rayos cósmicos al fluido es
el medio es transparente - no interactúa -, pero
el campo magnético, el cual, a su vez, esta congelado
la de abajo difunde pues el medio es opaco. Si
en las lı́neas de fluido en un plasma de conductivi-
la velocidad de la particula es V , entonces en el
dad infinita. Notemos que el tiempo de convección
primer caso el tiempo de cruce es tcruce ∼ R/V ,
no depende de las energı́as de las partı́culas, su no
mientras que el segundo tcruce ∼ R2 /D ≫ R/V .
de la velocidad del fluido y de la escala espacial del
problema.
4.2. Ecuación de transporte

Nos interesa estudiar cómo evoluciona dicha población de rayos cósmicos debido a la interacción
con los blancos en la región (campo magnético, materia y radiación) y a los fenómenos de transporte.
Caracterizamos la población de partı́culas relativistas a través de la función de distribución
dN
n(⃗r, E, t) = ,
dE dV
es decir, el número de partı́culas a tiempo t en la posición ⃗r por unidad de volumen dV con energı́a
en el intervalo dE alrededor de E.
La evolución de esta distribución viene determinada por la ecuación de transporte o ecuación
cinética, cuya forma es

∂n ⃗ ⃗ ∂ bn n
− ∇ · D∇n + + = Q(E, ⃗r, t),
∂t ∂E T
donde D = D(E, ⃗r) es el coeficiente de difusión, b = b(E, ⃗r, t) = dE/dt ≤ 0 es la tasa total de pérdida
de energı́a de las partı́culas por distintos procesos, T = T (E, ⃗r) es la escala temporal de escape de
las partı́culas de la región de interés y Q(E, ⃗r, t) es el término fuente de inyección.
La ecuación de transporte describe el balance entre los procesos que inyectan y remueven partı́cu-
las por unidad de tiempo y por unidad de volumen en el espacio de fases de las posiciones (y even-
tualmente direcciones de movimiento) y la energı́a. El segundo término del lado izquierdo da cuenta
del transporte espacial de partı́culas por difusión. En el caso más general el coeficiente de difusión
D es un tensor; si el medio es uniforme, D = D(E). En general, se asume que la dependencia en la
energı́a es una ley de potencias
D(E) = D0 E δ δ > 0,
34
pero esto puede variar significativamente en diversos medios. En ausencia de cualquier otra infor-
mación se suele adoptar el coeficiente de difusión mı́nimo, conocido como coeficiente de difusión de
Bohm
1
DB = rg c,
3
donde rg es el giroradio de las partı́culas. En algunos casos, la difusión compite con la convección
como mecanismo de transporte de las partı́culas (por ejemplo, en regiones donde hay vientos fuertes).
El tercer término representa las pérdidas de energı́a continuas, es decir, aquellas interacciónes que
producen cambios pequeños en la energı́a de las partı́culas respecto de su energı́a inicial. Aquı́, b es
la tasa (total) de pérdida de energı́a, si sobre las partı́culas actúan varios procesos de enfriamiento
i = 1, ..., N , resulta
N N
dE
b(E, ⃗r, t) = bi =
X X
.
i=1 i=4 dt i
El cuarto término del lado izquierdo da cuenta del escape de partı́culas de la región de interés; el
tiempo de escape es T . En este termino pueden introducirse de forma efectiva el escape por convección
y también por difusión - de manera de simplificar aun más la ecuación de transporte eliminando el
segundo término. En ese caso T −1 = t−1 dif + tconv , donde tdif y tconv son tiempos tı́picos de escape por
−1
difusión y convección, respectivamente. En el termino de escape también pueden incluirse efectos de

pérdidas de energı́a catastróficos, en los que la partı́cula pierde tanta energı́a que es removida del
intervalo de energı́a de interés. Si las partı́culas son inestables, también puede tenerse en cuenta en
este termino la desaparición de partı́culas por decaimiento; en este caso, T = tdec es el tiempo de
decaimiento.
Finalmente, el término de la derecha es el término fuente. Aquı́ se incluyen todos los fenómenos
que inyectan partı́culas en la región de interés. Por ejemplo, la inyección de partı́culas mediante algún
mecanismo de aceleración puede introducirse de manera efectiva en este termino utilizando alguna
función inyección adecuada.
4.3. Soluciones de la ecuación de transporte

4.3.1. Solucion general
La ecuación de transporte es una ecuación diferencial en derivadas parciales de tipo inhomogeneo.
Puede resolverse utilizando el metodo de la función de Green. Para D = D(E), b = b(E) y T = T (E),
la función de Green G(E, ⃗r, t) satidface la ecuación de transporte si el término de fuente lo podemos
representar como Q(E, ⃗r, t) = δ(E − E0 )δ(⃗r − ⃗r0 )δ(t − t0 ), donde E0 , ⃗r0 y t0 son los valores iniciales.
La solucion general es
1 (⃗r − ⃗r0 )2
" #
τ
G(E, ⃗r, t; E0 , ⃗r0 , t0 ) = exp − − δ(t − t0 − τ ),
|b(E)|(4πλ)3/2 T 4λ
con ˆ ˆ
E E
dE D(E)
τ ≡ τ (E, E0 ) = , λ ≡ λ(E, E0 ) = dE.
E0 b(E) E0 b(E)
Entonces, la solucion general de la ecuación de transporte es
ˆ +∞ ˆ ˆ t
n(E, ⃗r, t) = 3
d r0 dE0 dt0 Q(E0 , ⃗r0 , t0 )G(E, ⃗r, t; E0 , ⃗r0 , t0 ).
−∞ ∞
35
4.3.2. Caso estacionario, homogéneo y sin escape
Consideremos el caso sencillo de un sistema de partı́culas en estado estacionario en una región ho-
mogénea, donde puede despreciarse el escape (T → ∞) y en la que se inyectan partı́culas relativistas
con un espectro Q(E) = KE −p , con p > 0. La ecuación para n(E) es en este caso
d h i
b(E)n(E) = Q(E). (4.1)
dE
Para p ̸= 1 y suponiendo que b(E)n(E) → 0 para E → ∞, la solución resulta
K E −(p−1)
n(E) =
(p − 1) |b(E)|
4.3.3. Casos no estacionarios

Supongamos que en una región actúa, durante un intervalo acotado de tiempo, un mecanismo
de inyección continua de partı́culas que puede describirse mediante una función de la forma Q(E) =
KE −p (p > 0) para t ≤ t0 y Q(E) = 0 para T > t0 . Supongamos además que las pérdidas de energı́a
son del tipo AE 2 , con A una constante. Resolviendo la ecuación de transporte
∂n(E, t) ∂ h i
=− b(E)n(E, t) + Q(E, t),
∂t ∂E
y considerando el caso de una inyección instantánea en t = 0 del tipo Q(E, t) = KE p δ(t), la
distribución de partı́culas resulta
KE p (1 − AEt)p−2 si AEt < 1

(
n(E, t) =
0 si AEt ≥ 1
Notar que para p = 2 el espectro no se modifica en su forma, solo va cambiando la energı́a máxima
de las partı́culas de acuerdo a Emax (t) = 1/(At).
5. Procesos radiativos
5.1. Conceptos basicos
El concepto básico en el estudio de los procesos radiativos es el de intensidad:
Intensidad: número de partı́culas incidentes por unidades de área por unidad de ángulo sólido
y por unidad de tiempo sobre un cierto detector no identificado.
Designaremos a la intensidad de partı́culas tipo i por Ii . A partir de la intensidad, se define el

flujo de partı́culas como ˆ
F = I cos θ dΩ,
Ω
donde θ es el ángulo entre la dirección de movimiento de las partı́culas y la normal al área sobre
la que inciden. La integración se hace sobre todo el ángulo solido subtendido por las partı́culas. En
36
particular, para un flujo isotrópico de radiación F = πI y el número de partı́culas por unidad de
volumen es n = (4π/v)I, donde v es la velocidad de las partı́culas. Como nos interesan partı́culas
relativistas, v ≈ c.
Por lo general, las partı́culas que se detectan no tienen una única energı́a, ası́ que es conveniente
introducir una función de distribución en energı́as n(E):
número de partı́culas h i
n(E) = , n(E) = cm−3 erg−1
unidad de volumen × unidad de energı́a
El número de partı́culas por unidad de volumen con energı́as entre E y E + dE es entonces
n(E)dE, por lo que el número total de partı́culas por unidad de volumen se puede calcular como
ˆ ∞
n= n(E)dE, [n] = cm−3 .
0
El número de partı́culas por unidad de volumen con energı́as mayores a E y el correspondiente

flujo integrado son ˆ ∞ ˆ ∞
n(> E) = n(E )dE, F (> E) =
′
F (E ′ )dE
E E
La densidad de energı́a de las partı́culas se puede calcular como
ˆ ∞
U= E n(E)dE, [U ] = erg cm−3 ,
0
por lo que el flujo de energı́a será

c h i
S= U, S = erg cm−2 s−1
4π
Finalmente, consideremos una fuente localizada a una distancia d y que radia isotrópicamente.
En ese caso, su luminosidad puede definirse como
ˆ ∞ h i
L(> E) = 4πd 2
F (E ′ )dE, L = erg s−1
E
Para una fuente de radio R que radia isotrópicamente, L = 4πR2 S.
5.1.1. Sección eficaz

Consideremos un flujo de partı́culas de clase a con velocidad va que interaccionan con un blanco
de volumen dV = dA dx formado por partı́culas de tipo b.
El número dNi de interacciónes de tipo i que ocurren en el intervalo de tiempo dt en el volumen
dV será proporcional a
el número de partı́culas de tipo b que haya en dV ,
el número de partı́culas incidentes de tipo a que atraviesan dA en dt.
Si nb es la densidad de partı́culas de tipo b en el blanco y na es la densidad de partı́culas incidentes
de tipo a en el sistema de referencia del blanco, entonces
dNi = dσi (nb dV )(na va dt).
La constante de proporcionalidad dσi es la sección eficaz diferencial que caracteriza la probabilidad
de interacción de las partı́culas de tipo a y b a través del mecanismo i.
37
La sección eficaz total σi para la interacción i es la suma
sobre todos los posibles estados (valores de energı́a, momento y
dirección de movimiento) de las partı́culas después de la inter-
acción. Tiene unidades de área, [σ] = cm2 .
Tanto dσi como σi son invariantes relativistas.
La unidad de medida tı́pica para la sección eficaz de inter-
acción entre partı́culas es el barn (b): 1 b = 10−24 cm2 .
La sección eficaz total para la interacción entre dos tipos de
partı́culas se obtiene sumando las secciones eficaces de todos
los posibles procesos que pueden ocurrir durante la interacción
entre partı́culas
σtot =
X
σi .
i
La probabilidad relativa de que ocurra un cierto canal de reacción i es entonces

σi
Pi = .
σtot
A partir de la sección eficaz se pueden definir otras cantidades útiles. La primera es el camino
libre medio de una partı́cula en un medio de densidad numérica de blancos n
λ = (nσ)−1 .
El camino libre medio puede interpretarse como la distancia promedio que recorre la partı́cula
en el medio entre dos interacciones. Una cantidad relacionada es el tiempo de enfriamiento de una
partı́cula
λ
tenf = ,
κv
donde v es la velocidad de la partı́cula y κ es la ineslasticidad, que es la fracción de su energı́a inicial
E0 que pierde la partı́cula en una interacción
|∆E|
κ= .
E0
Para un proceso de pérdida de energı́a que pueda considerarse continuo, el tiempo de enfriamiento
se define a partir de la tasa de pérdida de energı́a como
1 dE
enf = −
t−1 .
E dt
A t−1
enf se lo suele llamar tasa de enfriamiento. Si las partı́culas pierden energı́a por más de un
proceso, entonces
−1
= t−1
X
tenf enf,i
i
El proceso de pérdida de energı́a más relevante será aquel que tenga el tiempo de enfriamiento
más pequeño. Notemos que esto dependerá tanto de la probabilidad de que la partı́cula interaccione
por ese mecanismo (representada por la sección eficaz) como de cuánta energı́a pierda la partı́cula
por interacción (inelasticidad).
38
En el caso general, tanto la sección eficaz como el camino libre medio y la inelasticidad dependerán
de la energı́a (y eventualmente de la dirección de movimiento) de las partı́culas. La sección eficaz
dependerá también en general de las propiedades de los blancos y en muchos casos se necesitarı́ntegrar
sobre estas variables.
La radiación γ se produce por la interacción de partı́culas de algún tipo a con distribución
na (Ea , Ωa , ⃗r) con un blanco de tipo b, formado por partı́culas materiales o por fotones, distribuidos
según nb (Eb , Ωb , ⃗r), La emisividad de la radiación γ resultante será
˛ ˆ ˛ ˆ
dσ
qγ (Eγ , Ωγ , ⃗r) = dΩa dEa dΩb dEb na (Ea , Ωa , ⃗r)nb (Eb , Ωb , ⃗r) × vr (1 − β⃗a · β⃗b )
dEγ dΩγ
donde vr es la velocidad relativa entre partı́culas y blancos, β⃗a · β⃗b = βa βb cos ψ , ψ es el ángulo de
colisión (ángulo entre las direcciones de movimiento de las partı́culas que interaccionan) y la sección
eficaz (doble) diferencial depende, en general, de Ea,b , ψ, Eγ y Ωγ . La emisividad representa el número
de fotones por unidad de energı́a, volumen y ángulo solido emitidos en ⃗r a tiempo t; en el sistema
cgs sus unidades son [qγ ] =erg−1 s−1 cm−3 sr−1 . Notar que la ecuación anterior es simplemente una
generalización de la tasa de colisiones (5.1.1), integrada sobre todas las variables que describen el
estado de partı́culas y blancos.
A partir de la emisividad definimos la luminosidad especifica Lγ integrando sobre el volumen de
la región de emisión ˆ
Eγ Lγ (Eγ ) = 4πEγ
2
qγ (Eγ , ⃗r)dV,
V
donde se supuso que la emisividad es isotrópica. La luminosidad integrada (en algún rango de
energı́as) tieneh unidades
i de potencia, [L] = erg s , mientras que las unidades de la luminosidad
−1
especifica son Lγ (Eγ ) = erg s−1 erg−1 .
5.2. Radiación térmica

La forma más sencilla de producir radiación electromagnética es calentando un gas. Al hacerlo, la
agitación térmica del gas aumenta y las partı́culas cargadas (generalmente electrones) son aceleradas
en interacciones coulombianas (generalmente con núcleos atómicos) en el gas. Ası́, estas partı́culas
producen emisión libre-libre (o Bremsstrahlung); como el gas esta en equilibrio termodinámico, la
emisión tiene un espectro térmico1 . Un caso particular a muchos objetos astronómicos (por ejemplo, el
CMB) es el de un cuerpo negro, que es un absorbente radiativo perfecto en equilibrio termodinámico
a una temperatura T (o sea, un plasma ópticamente grueso). El espectro de emisión de un cuerpo
negro en función de la energı́a esta dado por la distribución de Planck:
8π 2
" #
Eph h i
nBB (Eph ) = 3 3 , nBB = erg−1 cm−3.
h c exp(Eph /kT ) − 1
Aquı́, Eph es la energı́a de los fotones de frecuencia ν y k es la constante de Boltzman. El máximo

de la distribución de Planck ocurre para
T

Eph, max ≈ 1,59 kT ≈ 1,37 × 10 −10
MeV,
K
1
Aunque un espectro térmico no es sinónimo de un espectro de cuerpo negro. De hecho, el espectro termino de un
gas óptimamente delgado usualmente presenta lı́neas de emisión que dependen de las abundancias quı́micas del gas,
mientras que en un espectro de cuerpo negro no guarda esta información puesto que las lı́neas son reprocesadas.
39
que es la Ley de desplazamiento de Wien y Eph, max ∝ T .
La energı́a media ⟨Eph ⟩ de los fotones emitidos por un cuerpo negro puede obtenerse dividiendo
la densidad de energı́a ˆ ∞
UBB = Eph nBB (Eph )dEph ,
0
por el número total del fotones por unidad de volumen
ˆ ∞
NBB = nBB (Eph )dEph .
0
De aquı́ resulta que

T

⟨Eph ⟩ = 2,7kT ≈ 2,3 × 10−10 MeV
K
Para obtener fotones con energı́as medias de 1
GeV por radiación de cuerpo negro, las temperatu-
ras deben ser del orden de T ∼ 1013 K. Estas tem-
peraturas no pueden hallarse en objetos astrofı́sicos
usuales. Solo son posibles durante breves lapsos en
eventos explosivo, como las erupciones de rayos γ
(GRB) o el Big Bang.
La temperatura tı́pica en un GRB es del orden
de 1013 K. La densidad de fotones de una fuente
tipo cuerpo negro con esa temperatura es extrema-
damente alta: NBB ∼ 4,6 × 1046 cm−3 . El camino
libre medio de un fotón en un medio de densidad
nγ es λγ ∼ (nγ σγγ )−1 , donde σγγ es la sección efi-
caz para la aniquilación de dos fotones en un par
electrón-positrón. Para un valor tı́pico de la sección
eficaz σγγ ∼ 10−24 cm2 ,
λγ ≪ 1 cm,
luego, la fuente será auto-absorbida y los fotones no podrán escapar de ella.
De todo esto se concluye que las fuentes de rayos γ que se observan en el continuo son de origen
no térmico.
5.3. Radiación sincrotrón

5.3.1. Radiación sincrotrón de una partı́cula
El movimiento de una partı́cula de carga e y masa m en un campo magnético queda determinado
por la fuerza de Lorentz
d(γm⃗v ) e
= ⃗v × B⃗ .
dt c
Como la fuerza es perpendicular a ⃗v , entonces v = cte. Como la fuerza también es perpen-
dicular a B,
⃗ la componente de la velocidad paralela al campo permanecerá constante, v∥ = cte.
Entonces, v⊥ = v 2 − v∥2 = cte. De aquı́ que la partı́cula se mueva describiendo una hélice con su
q
eje paralelo al campo magnético. Al ángulo α que forman la velocidad de la partı́cula y el cam-
po magnético se lo conoce como pitch angle. La frecuencia de giro en el plano normal a B ⃗ es
40
eB mc2 eB
ωB = = ,
mc E γmc
donde se introdujo el factor de Lorentz γ de la partı́cula, tal
que E = γmc2 .
La expresión relativista para potencia total radiada por una
carga acelerada esta dada por
2e4
P = γ 2 v⊥
2 2
B .
3m2 c5
Usando que β = v/c = v⊥ /c sen α e introduciendo el radio clásico del electrón re = e2 /me c2 , la
potencia total puede escribirse como
2 me
2
P = cre2 β 2 γ 2 B 2 sen2 α.
3 m
Si la partı́cula radiada es ultra relativista, β ∼ 1 y la potencia total radiada resulta
2 me
2
P = cre2 γ 2 B 2 sen2 α.
3 m
Gran parte de la radiación es emitida en un cono de ángulo de apertura θ ∼ 1/γ en la dirección
del momento de la partı́cula, por lo que la emisión será tanto más colimada cuando más energética
sea la partı́cula.
La pérdida de energı́a de la partı́cula por unidad de tiempo se obtiene directamente a partir de
la expresión para la potencia
2 me
2
dE

Psincr =− =− cre2 γ 2 B 2 sen2 α. (5.1)
dt 3 m
sincr
Para un dado valor del factor de Lorentz, las pérdidas dependen de la masa de la partı́cula como
m−2 , por lo que son ∼ 106 veces más importantes para electrones que para protones.
Definiendo la sección eficaz de Thompson σT = (8π/3)re2 ∼ 6,65×10−25 cm2 y promediando sobre
el ángulo α suponiendo una distribución isotrópica (P (α)dα = 1/2 sen α dα, α ∈ [0, π]) se obtiene
que
dE 4 me 2 me 2 B 2 2

=− cσT Umag γ ∼ −6,66 × 10−4 γ eV s−1 ,
dt 3 m m G
sincr
con Umag densidad de energı́a magnética.
La ecuación (5.1) da la pérdida total de energı́a por radiación sincrotrón de una partı́cula con
factor de Lorentz γ. La emisión no es monoenergética; su distribución en energı́a esta dada por
√ 3 ˆ
3 e B sen α Eph ∞
P (γ, Eph , α) = K5/3 (ζ)dζ, (5.2)
hmc2 Ec Eph /Ec
con K5/3 la función de Bessel modificada de segunda especia y de orden 5/3. La energı́a caracterı́stica
Ec se define como
3 ehB sen α 2
Ec = γ . (5.3)
4π mc
41
La función de la ecuación (5.2) tiene un máximo agudo para Eph,max ∼ 0,29Ec . En unidades
correspondientes, es aproximadamente,
me B 2

Eph,max ≈ 5,1 × 10 −15
γ sen α MeV.
m G
Entonces, para que un electrón con γ ∼ 103 emita rayos γ de energı́a 5 GeV hace falta un campo
muy alto, B ∼ 1012 G.
5.3.2. Radiación sincrotrón de una distribución de partı́culas

Consideremos ahora una distribución en energı́a de partı́culas n(E, α) en algún intervalo de energı́a
E min ≤ E ≤ E max . En ese caso, el espectro sincrotrón puede calcularse integrando la ecuación (5.2)
en la energı́a de las partı́culas y en el pitch angle
ˆ ˆ E max
P (Eph ) = P (E, Eph , α)n(E, α)dE dΩα .
Ωα E min
Las unidades de esta ecuación son de potencial por unidad de energı́a por unidad de volumen,
[P (Eph )] =s−1 cm−3 . Supongamos que la distribución de partı́culas es isotrópica y su forma funcional
es la de una ley de potencias en la energı́a
n(E)dE = K0 E −p dE, [n] = erg−1 cm−3 , (5.4)
donde K0 es una constante y p es el ı́ndice espectral.

El punto importante es que el espectro sincrotrón de una distribución de partı́culas tipo ley de
potencias es otra ley de potencias en la energı́a de los fotones
p−1
P (Eph ) ∝ Eph
−δ
, δ= .
2
Otra caracterı́stica de la radiación sincrotrón es que esta intrı́nsecamente polarizada. En el caso

de un campo magnético homogéneo el grado de polarización lı́neal
p+1
Π0 (p) = .
p + 7/3
Esto da valores de 69−75 % para p = 2−3. Si el campo magnético tiene una componente aleatoria
el nivel de polarización será menor,
!
B02
Π(p) = Π0 (p) ,
B02 + Br2
donde Br es la componente aleatoria del campo.
5.3.3. Absorción de la radiación sincrotrón

Los fotones sincrotrón pueden ser absorbidos por los propios electrones en presencia de un campo
magnético. Esto lleva a una modificación del espectro sincrotrón a bajas frecuencias.
42
Figura 9: Espectro sincrotrón generado por una distribución de partı́culas N (Ee ) ∝ Ee−2 . La curva
solida es la distribución espectral en energı́a (SED), mientras que las curvas a rayas son la contribución
de electrones de distintas energı́as y las curvas punto-raya y punto-punto-raya son ajustes lineales en
distintas porciones del espectro. La variable graficada es ϵL(ϵ) = ϵP (ϵ).
Para una distribución en energı́a de las partı́culas tipo ley de potencias como la de la ecuación
(5.4), el coeficiente de absorción de la radiación es
αν = Aν −(p+4)/2 ,
donde ν es la frecuencia y A una función complicada de p; A = A(p).

Luego, la intensidad de la radiación resultante será
jν
I(ν) = 1 − exp(−αν l ,
αν
donde jν es la emisividad de la fuente (cantidad de energı́a emitida por unidad de tiempo, de fre-
cuencia, de ángulo solido y de volumen) y l es su dimensión lineal.
Cuando τ ≡ αν l ≪ 1, la fuente es transparente a su propia radiación y decimos que es ópticamente
delgada. En ese caso,
jν
I(ν) = lı́m 1 − exp(−αν l ≈ jν l ∝ ν −δ .
αν →0 αν
Si en cambio τν ≫ 1, la fuente es opaca u ópticamente gruesa y

jν ν −(p−1)/2
I(ν) ∼ ∝ −(p+4)/2 = ν 5/2 .
αν ν
La condición τ ≈ 1 divide entonces el espectro de la fuente en dos regiones. Si llamamos νa a la
frecuencia critica tal que τ (νa ) = 1, entonces
1. ν < νa =⇒ τν > 1: la fuente es ópticamente gruesa y I(ν) ∝ ν 5/2 ,
2. ν > νa =⇒ τν < 1: la fuente es ópticamente gruesa y I(ν) ∝ ν δ ,
43
Como τν = αν l la frecuencia crı́tica depende del tamaño
de la fuente. Una fuente en expansión, como un jet, presen-
ta regiones de diferentes tamaños que se hacen ópticamente
delgadas a diferentes frecuencias. El resultado es un espectro
chato, combinación de los distintos espectros emitidos en dis-
tintas zonas.
La presencia de un plasma térmico también puede modifi-
car el espectro sincrotrón, ya que el plasma puede absorber la
radiación. El coeficiente de absorción en este caso es αν ∝ ν −2 ,
por lo que
jν
I(ν) ∼ ∝ ν −(α−2)
αν
5.3.4. Lı́mite cuántico

Si el campo magnético donde son inyectadas las partı́culas relativistas es muy intenso, la apro-
ximación clásica para la radiación sincrotrón de las secciones anteriores deja de ser válida. En este
régimen, además, comienzan a ser probables la aniquilación y creación de pares en un campo magnéti-
co B (en ausencia de campos, los procesos están prohibidos)
e+ + e− + B → γ + B,
γ + B → e+ + e− + B.
Para electrones, la probabilidad de que sucedan estos fenómenos cuánticos es apreciable cuando
la energı́a caracterı́stica de los fotones sincrotrón en la aproximación clásica Ec. (5.3) es del orden de
la energı́a Ee del electrón. Esto ocurre para Ee ≳ Ecrit con
3
Bc m2e

Ecrit = me c 2
, Bc eℏ ≈ 4,4 × 1013 G.
B c
con Bc el campo magnético critico. Ası́, por ejemplo, en un campo magnético B = 107 G, electrones
con energı́a Ee ≳ 2,2×1012 eV serán capaces de radiar fotones sincrotrón con energı́as suficientemente
altas para crear pares, iniciando cascadas electromagnéticas.
5.3.5. Efecto sobre el espectro de electrones

Las pérdidas radiativas sufridas por las partı́culas relativistas modifican su distribución en energı́a.
Para determinar como, consideremos la solución de la ecuación (4.1).
Para pérdidas por radiación sincrotrón b(E) ∝ E 2 , luego n(E) ∝ E −(p+1) . O sea, que el espectro
de los electrones se hace más blando respecto de la inyección, incrementándose en uno la potencia de
E. (Ejemplo: Figura 10).
5.4. Radiación de curvatura

El movimiento de una partı́cula relativista de carga q = Ze y energı́a E en un campo magnético
B es una hélice de pitch angle α y radio de curvatura
⃗
E
rg = ,
qB sen α
que sigue la curvatura de la lı́neas de campo.
44
Luego, la curvatura de la trayectoria de la partı́cula (y por
ende su espectro radiativo) se debe tanto a la curvatura de
la hélice (determinada por los valores de E, B y α) como a
la curvatura de las lı́neas de campo Rc . Cuando la contribu-
ción más importante a la curvatura de la trayectoria es la del
movimiento helicoidal (rg ≪ Rc ), el espectro radiativo queda
determinado por la intensidad de B; ⃗ este es el caso de la ra-
diación sincrotrón. En cambio, si la hélice esta muy estirada
(rg ≫ Rc ), la mayor contribución a la curvatura de la trayec-
toria se debe a la de las lı́neas de campo; en este caso se dice
que la partı́cula emite radiación de curvatura. Se suele afir-
mar que la radiación de curvatura es aquella que emiten las
partı́culas cuando se mueven sobre las lı́neas de campo. Esto Figura 10: Espectro de partı́culas re-
no es estrictamente correcto, pero es una buena aproxima- lativistas en estado estacionario con
ción para B muy grandes ya que en ese régimen las intensas perdidas b(E) = AE 2 . A energı́as al-
pérdidas de energı́a radiativas producen que el pitch angle se tas (donde dominan las perdidas por
reduzca considerablemente. radiación sincrotrón) la distribución
Si el valor de B es suficientemente alto como para poder se quiebra haciéndose mas blanda.
aproximar que la curvatura de la trayectoria de las partı́culas
es igual a la de las lı́neas de campo, entonces la radiación se
producirá puramente en el régimen de curvatura. En este caso, el espectro de emisión de un electrón
con energı́a Ee = γe me c2 es el mismo que en el caso de la radiación sincrotrón - excepto que el
giroradio se reemplaza por el radio de curvatura Rc de la lı́nea de campo
√ 2 ˆ
3 e γe Eph ∞ 3 hc 3
P (Ee , Eph ) = K5/3 (ζ)dζ, Ec = γ . (5.5)
h Rc Ec Eph /Ec 4π Rc e
La energı́a perdida por el electrón es radiada en la dirección del movimiento dentro de un cono de
ángulo 1/γe . La tasa total de pérdida de energı́a se calcula integrando el espectro de emisión sobre
la energı́a de los fotones
ˆ
dEe 2 e2 c 4
= − P (Ee , Eph )dEph = − γ .
dt 3 Rc2 e
curv
Notar que, en esta aproximación, la tasa de pérdidas no depende de la intensidad del campo
magnético sino solo de su radio de curvatura.
Si se tiene una distribución de electrones N (Ee ) en un campo con curvatura Rc , el espectro
total emitido por radiación de curvatura se obtiene integrando la ecuación (5.5) en la energı́a de los
electrones ˆ Eemax
P (Eph ) = P (Ee , Eph )N (Ee )dEe .
Eemin
Nuevamente, en esta aproximación el resultado no dependerá de B sino de Rc . Cuanto menor sea
Rc , mayor será la potencia radiada.
Si una partı́cula se mueve alrededor de una lı́nea de campo con curvatura finita, en general emitirá
tanto radiación de curvatura como sincrotrón. Se pueden igualar las pérdidas sincrotrón con las de
curvatura para definir un ángulo critico αcrit , de tal forma que la radiación de curvatura domina para
valores del pitch angle α < αcrit :
γme c2 rg
sen(αcrit ) = = .
eBRc Rc
45
Figura 11: Diagrama de la interacción Compton inversa en el sistema del laboratorio (izquierda) y
el sistema en reposo del electrón (derecha).
Cuando α ∼ αcrit ambos procesos deben tenerse en cuenta. Si el campo magnético es regular, la
radiación de curvatura estará polarizada como la sincrotrón.
La radiación de curvatura es un proceso importante en la región polar de los púlsares. Allı́ el
radio de curvatura de las lı́neas de campo magnético vale Rc ∼ (cR∗ /Ω)1/2 , donde R∗ es el radio de
la estrella de neutrones y Ω su velocidad de rotacion.
5.5. Radiación Compton inversa (IC)

Cuando un fotón de energı́a Eph es dispersado por un electrón de energı́a Ee
e− + γ → e− + γ,
el electrón puede ceder energı́a al fotón y producir rayos gamma. Este proceso se denomina efecto
Compton inverso.
En la Figura 11 mostramos un esquema de la interacción en el sistema propio del electrón K ′
(donde estaba inicialmente en reposo) y en el sistema de laboratorio K, donde el electrón está en
movimiento antes de la colisión. A partir de la conservación de la energı́a y del momento puede
hallarse fácilmente la energı́a del fotón luego de la colisión en el sistema propio del electrón
′
Eph
Eγ′ = ,
1 + (Eph
′
/me c2 )(− cos θ)
donde θ′ es el ángulo de dispersión (aquel que forman las direcciones de movimiento del fotón antes
y después de la interacción). Notemos que Eγ′ < Eph ′
, por lo que el electrón le cede energı́a al fotón.
La energı́a final del fotón en el sistema de referencia del laboratorio es entonces
Eγ = γe Eγ′ (1 + βe cos θ2′ ),
q
donde γe es el factor de Lorentz del electrón y βe = 1 − γe−2 .
Si Eph
′
≪ me c2 la interacción ocurre en el llamado lı́mite de Thomson; la condición equivalente en
el sistema de laboratorio es Ee Eph ≪ m2e c4 . En este régimen, la colisión es casi elástica en el sistema
en reposo del electrón y Eγ′ ≈ Eph
′
. La máxima energı́a que puede alcanzar el fotón dispersado es
Eγ,max ≈ 4γe2 Eph ,

que corresponde al caso de una colisión frontal. Aunque la energı́a caracterı́stica de los fotones
dispersados es grande, aun es mucho menor que la del electrón - que en este lı́mite solo pierde una
pequeña fracción de su energı́a en cada interacción.
46
Consideremos un fotón de energı́a Eph que se mueve formando un ángulo θ respecto de la dirección
de la velocidad del electrón. En el sistema en reposo del electrón este ángulo vale
sen θ 1
!
θ
tan θ =
′
≈ − cot ,
γe (cos θ − βe ) γe 2
si el electrón es muy energético (β ≈ 1). Entonces, el electrón ve a los fotones incidir sobre el en
dirección contraria a la de su movimiento y formando un cono de semiapertura ∼ 1/γe .
Algo parecido ocurre con el ángulo de dispersión de los fotones en el sistema de laboratorio. Su
valor es
cos θ2′ + βe
cos θ2 = ,
1 + βe cos θ2′
donde θ2′ es el ángulo de dispersión en el sistema en reposo del electrón. Cuando βe ∼ 1, cos θ2 ∼ 1,
por lo que los fotones sin dispersados en la dirección del movimiento del electrón antes de la colisión,
dentro de un cono de semi apertura pequeña ∼ 1/γe .
5.5.1. La sección eficaz IC

La sección eficaz diferencial exacta para la interacción Compton está dada por la fórmula de
Klein-Nishina
!2
Eγ′ ′
Eγ′ ′
! !
dσKN re2 Eph Eph
= + − sen 2 ′
θ δ E ′
− ,
dEγ′ dΩ′γ 2 ′
Eph Eγ′ ′
Eph γ
1 + (Eph
′
/me c2 )(1 − cos θ′ )
donde re es el radio clásico del electrón. En el lı́mite de Thompson, Eγ′ ≈ Eph

′
y la sección eficaz se
reduce a
dσKN re2
≈ 1 + cos 2 ′
θ δ E ′
− E ′
ph .
dEγ′ dΩ′γ 2 γ
La sección eficaz total σIC es un invariante relativista, ası́ que podemos integrar la sección eficaz
de Klein-Nishina. Si definimos x = Ee Eph /m2e c4 , la sección eficaz total (promediada en ángulo) en el
sistema de laboratorio resulta
3σT 2 2 1 4 1
" #

σIC = 1 − − 2 ln 1 + 2x + + − ,
8x x x 2 x 2(1 + 2x)2
donde σT es la sección eficaz de Thomson. En la Figura 12 se grafica la sección eficaz total. Se

observa claramente que existen dos regı́menes distintos de interacción. En el régimen de Thomson,
para x ≪ 1,
σIC ≈ σT (1 − 2x), x ≪ 1.
Para x ≫ 1 la sección eficaz total decrece abruptamente,
3 11
σIC ≈ σT x−1 ln 4x − , x ≫ 1.
8 6
Este lı́mite se conoce como régimen de Klein-Nishina. En este régimen el electrón transfiere casi
toda su energı́a al fotón en una única colisión y Ee ∼ Eγ . Sin embargo, como la sección eficaz
se hace mucho más pequeña que la de Thomson, el flujo γ resultante (proporcional al número de
interacciones) se reduce considerablemente.
47
Figura 12: Sección eficaz de la interacción Compton inversa en función de la variable x = Ee Eph /m2e c4
(panel superior) y de la energı́a del electrón Ee para una energı́a inicial del fotón fija ϵ = 1 eV (panel
inferior).
48
5.5.2. Espectro de emisión
Para una distribución de energı́a de electrones ne (Ee , Ωe , ⃗r) que interactúa con una distribución de
fotones nph (Eph , Ωph , ⃗r), ambas en unidades de erg−1 cm−3 sr−1 , la emisividad de fotones, con unidades
de erg−1 s−1 cm−3 sr−1 , puede escribirse el caso más general como
ˆ ∞ ˛ ˆ Eemax ˛
dσKN
qγ (Eγ , Ωγ , ⃗r) = dEph dΩph dEe dΩe c(−βe cos ψ) ×ne (Ee , Ωe , ⃗r)nph (Eph , Ωph , ⃗r).
0 Eemin dEγ dΩγ
Aquı́, ψ es el ángulo de colisión (aquel entre las direcciones de movimiento del fotón y el electrón
antes de la dispersión). Esta expresión puede simplificarse considerablemente bajo ciertas suposicio-
nes. Si γe ≫ 1, entonces la dirección de movimiento del fotón dispersado será aproximadamente la
misma que la del electrón antes de la colisión: Ωγ ∼ Ωe . Entonces,
dσKN dσKN
∼ δ(Ωγ − Ωe ), (5.6)
dEγ dΩγ dEγ
y podemos integrar sobre Ωγ . A esta aproximación (ecuación (5.6)) para la sección eficaz doble
diferencial se la conoce como aproximación head-on. La expresión entonces nos queda
4πre2 2z 2z 2
" #
dσKN z2
= 1+ − + 2 ,
dEγ bθEe 2(1 − z) bθ (1 − z) bθ (1 − z)2
con z = Eγ /Ee < bθ /(1 + bθ ) y bθ = 2(1 − βe cos ψ)Eph Ee /(me c2 )2 . En la aproximación head-on,
entonces
ˆ E max ˛ ˆ Eemax
ph dσKN
qγ (Eγ , Ωγ , ⃗r) = c dEph dΩph dEe (1 − βe cos ψ) × ne (Ee , Ωe , ⃗r)nph (Eph , Ωph , ⃗r).
min
Eph Eemin dEγ
(5.7)
Los lı́mites para Eph permitidos por la cinemática son
z 1 2Eγ

min
Eph = , max
Eph = .
1−z 2γe (1 − βe cos ψ) 1 − βe cos ψ
Un caso importante es aquel en que el campo de fotones blanco y la distribución de electrones
son isotrópicos
1 1
nph (Eph , Ωph , ⃗r) = nph (Eph , ⃗r), ne (Ee , Ωe , ⃗r) = ne (Ee , ⃗r)
4π 4π
Bajo esta aproximación el ángulo de colisión ψ también estará distribuido isotrópicamente, por
lo que es posible promediar la ecuación (5.7) sobre Ωph e integrar en Ωγ . La expresión final para la
emisividad (ahora en unidades de erg−1 s−1 cm−3 ) es
ˆ Eemax ˆ E max
ph
qγ (Eγ , ⃗r) = 8πre c
2
dEe ne (Ee , ⃗r) dEph nph (Eph , ⃗r)f (Eγ , Ee , Eph ), (5.8)
Eemin min
Eph
con f (Eγ , Ee , Eph ) una función muy complicada donde aparecen z = Eγ /Ee y b = 4Eph Ee /(me c2 )2 .
Los lı́mites de integración sobre Eph quedan determinados por la condición
Eph b
≤z≤ .
Ee 1+b
49
Un caso sencillo para el que, bajo ciertas aproximaciones, es posible resolver la ecuación (5.8) es
el de una distribución isotrópica tipo ley de potencias de electrones,
ne (Ee ) = K0 Ee−p , Eemin < Ee < Eemax ,
que interacciona en el régimen de Thomson con una distribución de fotones blanco nph (Eph ). En un
intervalo de energı́a de los fotones dispersados para el que se satisfaga que
1 1
s s
Eγ Eγ
γemin < , γemax ≫ ,
2 Eph 2 Eph
podemos integrar la ecuación (5.8) para dar

ˆ
(p−1)/2
gγ (Eγ ) ∝ Eγ−(p+1)/2 dEph Eph nph (Eph ).
O sea que el espectro de fotones dispersados es también una ley de potencias en la energı́a Eγ de
ı́ndice espectral δ = (p+)/2.
5.5.3. Tasa de pérdida de energı́a

El calculo de la tasa de pérdida de energı́a para un electrón de energı́a Ee por interacciones IC
con un campo de fotones puede calcularse de manera similar a lo hecho en la sección anterior. Se
debe tener en cuenta la energı́a perdida por el electrón en cada colisión e integrar sobre todos los
posibles estados finales del fotón dispersado. De manera general tenemos
ˆ ˛ ˆ ˛
dEe dσKN
− (Ee ) = c Eγ dΩγ dEph dΩph (Eγ − Eph )nph (Eph , Ωph ) × (1 − βe cos ψ) .
dt dEγ dΩγ
El resultado toma una forma sencilla en el régimen de Thomson cuando el campo de fotones
blanco es isotrópico
ˆ
dEe 4
= − cσT Uph γe2 , Uph = dEph Eph nph (Eph ),
dt 3
Th
con Uph la densidad de energı́a del campo de fotones blanco. Esta expresión es formalmente idéntica
a la tasa de pérdida de energı́a de un electrón por radiación sincrotrón si la densidad de energı́a
del campo magnético Umag se reemplaza por la del campo de fotones. En el régimen de Thomson,
la importancia relativa de las pérdidas de energı́a sincrotrón e IC queda fijada por la densidad de
energı́a de los campos blanco
tIC, Th Umag
= .
tsincr Uph
En el régimen de Klein-Nishina, para un campo de fotones monocromático de energı́a ⟨Eph ⟩, la
tasa de pérdidas es
!3 "
3 4γe ⟨Eph ⟩ 11
! #
dEe me c2
≈ − cσT Uph ln − ,
dt 8 ⟨Eph ⟩ me c 2 6
KN
50
donde Uph = n⟨Eph ⟩ y n es la densidad numérica de fotones blanco. Mientras que en el régimen de
Thomson la tasa de pérdidas depende fuertemente de la energı́a de los electrones como ∝ γe2 en el
régimen de Klein-Nishina la dependencia es muy débil, ∝ ln γe .
Entonces, las pérdidas radiativas sufridas por las partı́culas relativistas modifican su distribución
en energı́a en estado estacionario respecto del espectro de inyección. Consideremos la ecuación (4.1)
para analizar el efecto de las pérdidas de energı́a por IC. En el régimen de Thomson, b(E) ∝ Ee2 ,
luego n(Ee ) ∝ Ee−(p+1) . O sea, que el espectro de los electrones se hace más ”blandorespecto de la
inyección, incrementándose en uno el ı́ndice espectral. En el régimen de Klein-Nishina, por otro lado
b(E) ≈ cte y entonces n(Ee ) ∝ Ee−(p−1) . El ı́ndice espectral disminuye en uno y el espectro de los
electrones se hace más duro respecto del de inyección. Si el mecanismo de enfriamiento dominante
es el IC y la energı́a de los electrones abarca un intervalo suficientemente amplio como para que la
interacción ocurre en ambos regı́menes, la distribución en energı́a de los electrones presentara un
quiebre debido al cambio en el ı́ndice espectral aproximadamente para la energı́a de transición.
5.6. Radiación por producción de foto-mesones

La interacción de protones energéticos con fotones puede resultar en la creación de piones a través
de las reacciones
p + γ → p + π 0 , y también p + γ → n + π + .
Para que esto ocurra la energı́a del fotón en el sistema de reposo del protón debe superar el valor
umbral !
m π
ϵ′th = mπ 1 + = 144,7 MeV.
2mp
Para energı́as mayores también es posible la reacción de múltiples piones.
p + γ → p + aπ 0 + b(π + + π − ),
p + γ → n + π + + aπ 0 + b(π + + π − ) a, b ∈ N0 .
La vida media de los piones cargados es de 2,6×10−8 s. Decaen produciendo, con una probabilidad
del 99,9877 %, un neutrino y un muón, que a su vez decae en un electrón/positrón y otro neutrino
π + → µ + + νµ µ+ → e+ + νe ,
π − → µ− + νµ µ+ → e+ + νe .
Los piones neutros tienen una vida media mucho más corta, ∼ 8,4 × 10−17 s. Decaen el 98,798 %
de las veces en dos fotones
π 0 → 2γ.
En la Figura 13 se muestra un gráfico de la sección eficaz de creación de piones en colisiones pγ.
La sección eficaz para el proceso p + γ → p + π 0 tiene un pico de ∼ 0,25mb para ϵ′ ∼ 300 MeV y
luego decrece a menos de 30 µb cerca de ϵ′ ∼ 1 GeV. Algo similar ocurre para la sección eficaz del
proceso p + γ → n + π + , pero con un segundo pico de ∼ 0,1 mb para ϵ ∼ 700 MeV. El valor medio
de la sección eficaz total, para todos los canales de interacción, es del orden de ⟨σpγ ⟩ ∼ 0,1 mb.
Se propuso la siguiente parametrización
340 µb si 200MeV ≤ ϵ′ ≤ 500MeV,

(
σpγ ≈
120 µb si ϵ′ ≥ 500MeV
51
Figura 13: Sección eficaz para la creación de piones por interacciones protón-fotón. Los puntos son
datos experimentales y la curva solida la seccion eficaz programada.
La tasa de pérdida de energı́a para un protón de factor de Lorentz γp que interacciona con un
campo de fotones con una distribución de energı́a nph (ϵ) esta dada por
pγ ˆ ˆ
dEp mp c3 ∞
nph (ϵ) 2ϵγp
− (Ee ) = dϵ 2 ϵ′th dϵ′ σpγ (ϵ′ )κpγ (ϵ′ )ϵ′ .
dt 2γp2 ϵ′th /2γp ϵ
π
Aquı́, κpγ (ϵ′ ) es la inelasticidad, que se define como la fracción de su energı́a inicial que el protón
pierde en la interacción y también se puede aproximar en forma sencilla por una función de tipo
escalón.
En el caso de un campo de radiación termino de densidad de energı́a ωph con una energı́a media
por fotón ⟨kT ⟩ se obtiene
pγ ! !−2
dEp ωph ⟨kT ⟩
− (Ee ) = 1,8 × 10 10
eV s−1 .
dt erg cm−3 eV
π
Para un campo de radiación con una distribución en energı́as del tipo ley de potencias nph ∝ ϵ−2 ,
en cambio,
pγ
Ep 2
!
dEp ωph

− (Ee ) = 4,3 × 10 −7
eV s−1 .
dt erg cm−3 GeV
π
Para protones con una distribución en energı́a np (Ep , ⃗r) en el intervalo Ep,min ≤ Epp,max que
interaccionan con una distribución de fotones nph (ϵ, ⃗r) la emisividad de piones resultante puede
obtenerse en forma similar a lo visto para la producción de fotones a través de la interacción Compton
inversa, pero con las secciones eficaces adecuadas.
ˆ Epmax ˆ ∞
dσpγ
qπ (Eπ , ⃗r) = c
pγ
dEp np (Ep , ⃗r) dϵ (Eπ , Ep , ϵ)nph (ϵ, ⃗r).
Epmin 0 dEπ
Aproximadamente 1/3 de los piones creados serán π 0 , que decaerán en rayos γ. Esta emisividad,
producto del decaimiento de piones neutros creados en interacciones fotohadrónicas, puede calcularse
como
ˆ
dEp
qγ (Eγ ) = np (Ep )nph (ϵ)Φ(η, x) dϵ.
Ep
52
En esta expresión, np (Ep ) son las distribuciones en energı́a de los protones y fotones, respecti-
vamente, η = 4p /m2p c4 y x = Eγ /Ep . La función Φ(η, x) es una función partida que depende de los
valores de x y el valor de energı́a umbral η0 .
Un tratamiento alternativo más sencillo para calcular la emisividad de rayos γ utiliza la aproxi-
mación de la funcional δ. En este formalismo se supone que la energı́a de cada tipo e partı́cula que
se crea producto de la interacción puede tener un único valor. En el canal de producción de un único
pión por colisión la inelasticidad vale κ1 ∼ 0,2. Entonces, si la energı́a del protón es Ep , la energı́a
de los piones neutros es Eπ0 ≈ 0,2Ep y la de cada fotón Eγ ≈ 0,5Eπ0 = 0,1Ep . En el régimen de
producción de múltiples piones la inelasticidad vale κ2 ∼ 0,6. La energı́a pérdida por el protón se
divide casi completamente entre tres piones (π 0 , π − y π + ) llamados leading pions; cada uno tendrá
una energı́a Eπ ≈ 0,2Ep . La energı́a de cada rayo γ será nuevamente Eγ ≈ 0,1Ep . Teniendo, además,
en cuenta que en cada colisión pγ existe una probabilidad ξ ≈ 0,5 de que el protón se convierta en
un neutrón con la emisión de un π + , la emisividad de rayos γ resulta
qγ (Eγ ) ≈ 20(1 − ξP1 )νpγ (10Eγ )np (10Eγ ),
donde νpγ es la tasa de colisiones y P1 es la probabilidad de que la interacción proceda a través del
canal de creación de un único pión.
5.7. Creación de pares por interacciones foto-hadrónicas

′
Para energı́as del fotón en el sistema de referencia del protón mayores que Ephth
= 2me c2 = 1,022
MeV, de la interacción entre un protón relativista y un fotón pueden crearse pares
p + γ → p + e+ e− .
La sección eficaz diferencial es la de Bethe-Heitler

4α re2 2 2Ee′ − Ee′ + 1
" ! #
dσe±

= Ee′2− +Ee′2+ + Ee′ − Ee′ + ln − ,
′
dEe− ′3
Eph 3 me c2 Eph′
2
con α = 1/137 la constante de estructura fina y las energı́as

primadas son medidas en el sistema de referencia en repo-
so del protón. Se puede aproximar la sección eficaz total y
parametrizar la inelasticidad κe± .
Si bien la sección eficaz para la producción de pares en
colisiones protón-fotón es ∼ 100 veces más grande que la de
la creación de piones, la inelasticidad del proceso de creación
de pares es muy pequeña. κe± ≤ 2me /mp ≪ 1. Esto significa
que, apenas se supera el umbral para la creación de piones,
este ultimo proceso para a dominar las pérdidas radiativas de
los protones por colisiones con fotones.
La tasa de pérdida de energı́a para un protón con factor
de Lorentz γp que interacciona con una distribución de foto-
nes nph se calcula de forma análoga al caso de producción de Figura 14: Tasa de pérdida de energı́a
fotomesones de un protón debido a la interacción
pγ ˆ ˆ con el CMB.
dEp mp c3 ∞ nph (Eph ) 2γp Eph
− = dEph ′ ′
dEph σe± κe± Eph .
dt ± 2γp2 Eph′ th
/2γp Eph2 ′ th
Eph
e
53
La figura 14 muestra un gráfico de la tasa de pérdida de energı́a de un protón en función de su
energı́a debido a la interacción con fotones del fondo cósmico de radiación. Para energı́as mayores
a Ep ∼ ×1019 eV, las pérdidas por creación de mesones dominan completamente el enfriamiento.
Debido a la interacción con los fotones del CMB, no se espera que lleguen a la Tierra protones con
energı́as mayores a ∼ 1020 eV provenientes de fuentes a distancias mayores a unos 50 Mpc. Deberı́a
aparecer, por lo tanto, un quiebre en el espectro de rayos cósmicos a muy altas energı́as, efecto
conocido como de Greisen-Zarsepin-Kusmin (efecto GZK).
5.8. Producción de pares en tripletas

O triplet pair production. Un electrón o un positrón pueden producir un par electrón-positrón al
interactuar con un fotón a través de las reacciones
e− + γ → e− + e− + e+ ,
e+ + γ → e+ + e− + e+ .
Para que el proceso sea posible, la energı́a de las partı́culas debe superar cierto valor umbral.
Medida en el sistema en reposo del electrón (positrón), la energı́a del fotón debe ser Eph
′
> 4me c2 .
La sección eficaz para esta interacción se puede aproximar como
σTPP (ϵ′ph ) = ασT f ϵ′ph ),
donde α = 1/137, ϵ′ph = Eph ′
/me c2 es la energı́a (adimensional) del fotón en el sistema de referencia
en reposo del electrón y f es una función que vale ∼ 1 para energı́as cercanas al umbral. Por ejemplo,
para un electrón de 10 GeV que interacciona con un fotón de 100 keV, σTPP ∼ = 10 mb.
La tasa de pérdida de energı́a para un electrón de energı́a Ee que interacciona con una distribución
de fotones n(ϵph ) es
√ ˆ ˆ
dEe 2 αcσT ϵmax n(ϵph 2ϵph Ee ′ ′ 1/2 ′
− ≈ dϵ 2 dϵph ϵph f (ϵph ).
dt Ee ϵmin ϵph 4
TPP
En el caso particular de una distribución monoenergética de fotones de energı́a ϵ y densidad nph

y aproximando f (ϵ′ ) ∼ 1,
8
1/2
dEe γe
− ≈ αcσT nph ,
dt 3 ϵ
TPP
donde γe = Ee /(me c ).
2
6. Procesos radiativos II
Estos procesos son producidos por la interacción de partı́culas relativistas con materia.
6.1. Radiación Cherenkov

La radiación Cherenkov ocurre cuando una partı́cula cargada viaja a través de un medio dieléctrico
con ı́ndice de refracción n(ω) > 1 con una velocidad que localmente excede la velocidad de la luz
en ese medio. Cuando la partı́cula pasa por el medio dieléctrico interactúa con las moléculas locales
induciéndoles una polarización que desaparece cuando la partı́cula se aleja. En el proceso se emite
un pulso electromagnético, en la mayorı́a de los medios en el rango del UV/visible.
54
La radiación se emite en un cono cuyo eje es paralelo a la velo-
cidad ⃗v de la partı́cula y que tiene un ángulo de semiapertura θ que
depende del ı́ndice de refracción n(ω) del medio,
c
cos θ = ,
v n(ω)
donde ω es la frecuencia de la radiación emitida y c/n, la velocidad

de la luz en el medio.
La energı́a radiada por unidad de frecuencia dω por unidad de
longitud dx recorrida por la partı́cula en el material esta dada por
2πe2 2πe2 v
" # " #
dE c2 dE c2
= 2 ω 1− 2 =⇒ = ω 1 − ,
dωdx c n (ω)v 2 dωdt c2 n2 (ω)v 2
Cher Cher
donde usamos que dx = vdt.

Una aplicación muy importante de la radiación Cherenkov es en la detección de rayos γ de muy
alta energı́a. Cuando estos entran en la atmósfera terrestre, disparan una cascada electromagnética.
Los electrones y positrones que se crean pueden tener energı́as tan altas que su velocidad excede la
de la luz en la atmósfera. La radiación Cherenkov que se emite puede detectarse al nivel del suelo,
proveyendo información sobre los rayos γ primarios.
Es importante enfatizar que la radiación Cherenkov es un proceso macroscópico, ya que involucra
caracterı́sticas globales del medio como el ı́ndice de refracción y que la radiación no es emitida por
la partı́cula sino por el medio.
6.2. Interacciones de electrones relativistas con materia

6.2.1. Bremsstrahlung relativista
Es la radiación producida cuando un electrón relativista es acelerado en el campo electrostático
de un núcleo u otra partı́cula cargada. Consideremos la producción de un fotón de energı́a Eγ por un
electrón relativista con energı́a Ee ≫ me c2 que es dispersado por un núcleo de carga Ze. La sección
eficaz diferencial para la interacción es
dσBr 4αre2 Z 2
(Ee , Eγ ) = ϕ Ee , Eγ ,
dEγ Eγ
donde la función ϕ depende del estadio de ionización de los blancos. Para un núcleo desnudo (es
decir, un gas totalmente ionizado) viene dada por
2 2Ee (Ee − Eγ ) 1
" 2 # ( " # )
Eγ Eγ

ϕ Ee , Eγ = 1+ 1− − 1− × ln − ,
Ee 3 Ee 2
me c E γ 2
mientras que para un núcleo que esta completamente apantallado por todos los electrones (por
ejemplo, hidrógeno neutro del medio interestelar) viene dada por
2 191 1
" 2 #
Eγ Eγ Eγ

ϕ Ee , Eγ = 1+ 1− − 1− ln + 1−
Ee 3 Ee Z 1/3 9 Ee
55
La emisividad de la radiación producida puede obtenerse como
ˆ
cna (⃗r ∞ dσBr h i
qγ (Eγ , ⃗r) =
Br
(Ee , Eγ )ne (Ee , ⃗r)dEe erg−1 cm−3 s−1 sr−1 ,
4π Eγ dEγ
con na la densidad de blancos y ne el espectro de electrones relativistas. Si consideramos que para

Ee ≫ me c2 , la distribución de electrones es
h i
ne (Ee , ⃗r) = Ke Ee−p erg−1 cm−3 ,
y su emisividad Bremsstrahlung al interactuar con atomos con apantallamietno completo resulta
c mm na (⃗r)
!
Ke A
qγBr (Eγ , ⃗r) = Eγ−p , χ0 = g cm−2
4π χ0 p−1 4α N0 Z 2 re2 ln(191Z −1/3 )
donde N0 es el número de Avogadro y A es el peso atómico. Es importante destacar que el espectro

de rayos γ resultante tiene la misma forma que el espectro original de electrones relativistas.
Las pérdidas de energı́a por Bremsstrahlung para los electrones pueden ser muy grandes, pudiendo
radiar fotones con energı́a Eγ ∼ Ee . Las pérdidas para un electrón individual vienen dadas por
ˆ Ee −me c2
dEe dσBr
− = cna Eγ (Ee , Eγ )dEγ .
dt 0 dEγ
Br
En el caso de un plasma ionizado, estas resultan
dEe h i
− = 4na Z 2 re2 αc ln(γe ) + 0,36 Ee .
dt
Br
En el caso de apantallamiento total, tenemos

 
dEe 1
− = 4na Z 2 re2 αc ln(183Z −1/3 ) − Ee .
dt 18
Br
En la Figura 15 se muestra a modo de ejemplo una SED por Bremsstrahlung relativista para una
distribución de electrones de tipo ley de potencias con ı́ndice p = 2,2, Ee,min = 1 MeV, Ee,max = 20
GeV y normalización arbitraria. El espectro para energı́as Eγ < Ee,min tiene pendiente +1, mientras
que para energı́as Ee,min < Eγ < Ee,max tiene pendiente α = −(p − 2). Notemos que la energı́a
máxima de los fotones es del orden de la energı́a máxima de los electrones.
56
Figura 15: SED de partı́culas en un campo de materia homogéneo. Las curvas punteadas y rayada
son ajustes lineales en distintas porciones del espectro.
6.3. Interacciones de protones relativistas con materia

6.3.1. Radiación por decaimiento de piones neutros
En la colisión inelástica de protones relativistas con protones de baja energı́a (colisiones pp) se
producen mesones. Los canales con energı́a umbral más baja corresponden a la creación de piones
p + p → p + p + aπ 0 + b(π + + π − ),
p + p → p + n + π + + aπ 0 + b(π + + π − ),
p + p → n + n + 2π + + aπ 0 + b(π + + π − ),
donde a, b ∈ N. El umbral de energı́a cinética de los protones relativistas para que sea posible la
creación de un único pión neutro es
!
mπ0
Tpth = 2mπ0 c 2
1+ ≈ 280 MeV,
4mp
donde mπ0 ≈ 135 MeV es la energı́a en reposo del pión neutro.
La vida media del pión neutro es τπ0 ≈ 8,5 × 10−17 s. Dicho pión
decae, con una probabilidad del 98,8 %, en dos fotones
π 0 → γ + γ.
Los piones cargados tienen una vida media de τπ± ≈ 2,6 × 10−8
s y decaen, con una probabilidad de ∼ 99,99 %, como
π + → µ + νµ , π − → µ− ν µ
Consideremos el decaimiento de los piones neutros. Si el π 0 decae en reposo, cada uno de los
fotones tendrá una energı́a
1
Eγ = mπ0 c2 ≈ 67,5 MeV.
2
57
Sin embargo, los piones neutros rara vez son creados en reposo y por tanto decaen en vuelo.
Consideremos el decaimiento de un pión que en el sistema de laboratorio tiene una velocidad de
vπ = βπ c. En el sistema de referencia del pión los dos fotones salen formando ángulos α y π + α con
la dirección de la velocidad del pión, afı́n de que se conserve el momento.
La distribución de los fotones es isotrópica en el sistema propio de los piones. El número de
fotones emitidos entre α y α + dα
dNγ
dα = sen α dα.
dα
Si la energı́a de uno (cualquiera) de los fotones en el sistema
CM de Eγcm = mπ0 c2 /2, entonces en el sistema de laboratorio
su energı́a es
E cm (1 − βπ cos α)
Eγ = γ q .
1 − βπ2
Luego, derivando y reemplazando obtendremos que
dNγ 2
=q ,
dEγ Eπ2 − m2π0 c4
que es independiente de la energı́a Eγ . El rango permitido de

energı́as de los fotones es
!1/2 !1/2
1 − βπ 1 + βπ
Eγcm ≤ Eγ ≤ Eγcm (6.1)
1 + βπ 1 − βπ
La distribución espectral de fotones inyectados en el decaimiento de piones neutros con diferentes
energı́as Eπ , donde en cada valor Eπ el espectro de fotones es una constante (rectángulo) que se
extiende entre los lı́mites dados por la ecuación (6.1). El centro de cada rectángulo se ubica en
Eγcm = 67,5 MeV; la altura varia con Eπ ya que el área bajo la curva debe ser igual a 2 (pues se crean
dos fotones por decaimiento). El espectro total es la suma de los espectros individuales.
Supongamos ahora que, en una dada región, protones relativistas con una distribución en energı́a
Np (Ep ) interaccionan con un blanco de protones de baja energı́a y densidad numérica nH . La emisi-
vidad inyectada de piones neutros por colisiones inelásticas pp es
ˆ Epmax
dσπ
qπ (Eπ ) = cnH Np (Ep ) (Ep , Eπ )dEp ,
Epmin dEπ
donde tenemos que dσπ /dEπ es la sección eficaz diferencial para la producción de un pi0 de energı́a
Eπ por un protón de energı́a Ep .
Existen diferentes aproximaciones adecuadas de la sección eficaz. En la aproximación de la funcio-
nal delta, se asume que, dado un valor de energı́a cinética Tp del protón relativista, todos los piones
neutros se crean con la misma energı́a Eπ = κTp . Entonces,
dσπ
(Ep , Eπ ) = δ(Eπ − κTp )σpp (Ep ),
dEπ
donde σpp es la sección eficaz inelástica total para interacciones pp. En esta aproximación, la emisi-
vidad de piones neutros resulta
cnH Ep i Ep i

qπ (Eπ ) = Np mp c2 + σpp mp c2 + .
κ κ κ
58
Figura 16: SED generada por una distribución de protones en un campo de materia homogéneo.
La curva a rayas es un ajuste lineal en la porción del espectro que se comporta como una ley de
potencias.
Experimentalmente, sobre un amplio rango de energı́as (GeV a TeV) la inelasticidad vale aproxi-
madamente κ ∼ 0,17. Una aproximación precisa para la sección eficaz es
 !4 2
Epph
σpp (Ep ) = 34,3 + 1,88L + 0,25L2 1 −  mb,
Ep
donde L = ln(Ep /1 TeV) y Epph = Tpph + mp c2 ≈ 1,22 GeV.

Finalmente, una vez conocida la emisividad de piones neutros qπ (Eπ ), la emisividad de fotones
producto de su decaimiento puede calcularse como
ˆ ∞
qπ (Eπ ) m2 0 c4
qγ (Eγ ) = 2 dEπ , Eπmin = Eγ + π .
4Eγ
q
Eπmin Eπ2 − m2π0 c4
En la Figura 16 se muestra a modo de ejemplo una SED por pp producida por la interacción
de una población homogénea y estacionaria de protones relativistas con ı́ndice espectral p = 2 y
Epmax = 20 TeV con un campo de materia. Notemos que el espectro pp queda limitado al rango de
rayos γ y se extiende desde Eγ ≳ 0,1Eth hasta Eγ ≈ 0,1Epmax . La SED en general es plana, con un
ı́ndice α ≈ 0.
La tasa de pérdida de energı́a por colisiones inelásticas pp, para un protón de energı́a Ep que
interacciona con protones blancos de densidad numérica nH esta dada por
dEp
− = nnH Kpp Ep σpp (Ep )Θ(Ep − Epth ),
dt
pp
donde Θ es la función de Heaviside, Kpp ≈ 0,5 es la inelasticidad total y Epth es el umbral de energı́a
de los protones relativistas para la producción de piones.
59
6.3.2. Radiación por aniquilación protón-antiprotón
La reacción básica es
p + p → ξπ,
pero el canal p + p → π 0 esta prohibido por conservación del momento, mientras que p + p → π 0 + π 0
y p + p → π + + π − están prohibidas por conservación de la paridad. La primera reacción permitida
es
p + p → π+ + π− + π0.
La sección eficaz para esta reacción es del orden de σpp ∼ 2,4 × 10−26 cm2 . En general, la
multiplicidad es ξ > 3. Fotones γ se inyectan por decaimiento de π 0 .
El proceso de aniquilación con creación directa de fotones esta permitido pero la sección eficaz es
muy pequeña
p + p → γ + γ, σpp→γγ ∼ 3 × 10−30 cm2 .
6.3.3. Perdidas por ionización

Cuando un protón o un núcleo (de carga eZ y masa M ) atraviesa un medio, lo ioniza y por lo
tanto pierde energı́a. La tasa de pérdida de energı́a por ionización esta dada por
dE E

− = 7,62 × 10−9 Z 2 nβ −1 22,2 + 4 ln + 2 ln β 2 − 2β 2 eV s−1 ,
dt me c2
ion
donde E ≪ (M/me )M c2 es la energı́a de la partı́cula, β = v/c es su velocidad y n es la densidad

numérica del medio.
Aunque las pérdidas por ionización no conducen a radiación significativa por parte del núcleo,
su calculo puede ser valioso cuando se desea considerar interacciones de un núcleo que primero debe
atravesar un medio. A altas energı́as las pérdidas por creación de piones son dominantes.
Las pérdidas por cada interacción individual no son grandes en general, pero su efecto acumulativo
puede ser significativo. La energı́a necesaria para ionizar un átomo de H en el estado fundamental es
de solo 13,6 eV.
6.3.4. Interacciones pión-núcleo y pión-pión

Si la densidad de hadrones es alta pueden ocurrir interacciones de piones con núcleos o entre
piones.
En el primer caso la sección eficaz presenta un pico para energı́as del pión de ∼ 190 MeV con
valores σπp ∼ 200 mb. Luego, la sección eficaz cae rápidamente hasta unos 40 mb a 1 GeV. A altas
energı́as, la sección eficaz se vuelve suave y se puede parametrizar como
b
σπp = a + ,
p
donde p es el momento en GeV/c del pión. Los valores de a y b para π = , π − andan en el orden de 23
mb y 25 mb GeV/c, respectivamente.
Los π 0 tienen una vida media demasiado corta como para interactuar antes de decaer.
La sección eficaz para la interacción π − π es, a bajas energı́as, σπ−π ∼ 35 mb. A energı́as altas
la sección eficaz es incierta.
60
6.3.5. Interacción neutrón-protón
Se pueden producir neutrones (relativistas) mediante las interacciones
p + γ → n + π,
p + p → p + n + π + + aπ 0 + b(π + + π − ),
p + p → n + n + 2π + + aπ 0 + b(π + + π − ).
Dichos neutrones pueden interaccionar con protones antes de decaer
p + n → p + p + π − + aπ 0 + b(π + + π − ).
Para energı́as cinéticas del neutrón incidente mayores que 290 MeV la sección eficaz para colisiones
np permanece constante en
σnp ∼
= 39,5 ± 1,0 mb.
A energı́as menores la sección eficaz se incrementa hasta llegar a 1 b a energı́as del MeV. Por
debajo de 100 MeV la sección eficaz se comporta como E −1 .
6.4. Aniquilación de electrones y positrones

La aniquilación de electrones con positrones puede ser una fuente importante de rayos γ a través
de la reacción
e+ + e− → γ + γ.
Cuando dos partı́culas se aniquilan en reposo, la energı́a de los fotones resultantes es simplemente
Eγ = me c2 = 0,511 MeV.
61
La radiación de lı́nea de esta energı́a se suele llamar radiación de aniquilación. Si uno de los
leptones se mueve a gran velocidad, uno de los fotones emergentes tendrá gran energı́a, mientras que
el otro tendrá una energı́a ∼ me c2 .
La sección eficaz para la aniquilación de un positrón de energı́a Ee = γme c2 con un electrón en
reposo es
πre2 γ 2 + 4γ + 1 γ+3
" q #

σe± = ln γ + γ − 1 − √ 2
2 .
γ+1 γ2 − 1 γ −1
Para los casos
πre2 h i
σe± ≈ ln(2γ) − 1 , γ ≫ 1,
γ
πr2 v
σe± ≈ e, β = ≪ 1.
β c
Para un plasma de electrones y positrones, la tasa de aniquilacion es
3 σT c cm3

Re± = ln(x) + x−1/2 , Re± = ,
8 x s
donde x = γe+ γe− y la exactitud de esta aproximación es mejor que el 14 %.
Si se tiene una distribución no térmica ne± (Ee± ) de pares en el plasma, el número de aniquilaciones
por unidad de tiempo es
Ṅ = ne+ ne− Re± dEe+ de− dV
La aniquilación de pares e± puede ocurrir también con la emisión de un único fotón, pero en este
caso el electrón debe estar ligado a un átomo. El átomo asegura la conservación del momento. La
sección eficaz del proceso es
 5 4 2
 ≈ 4πZ γα re para γ ≫ 1,
σe1±ph = 5 4 2
≈ 4πZ 3α re β para β ≪ 1,
donde la energı́a del positrón es γme c2 y el núcleo atómico tiene carga eZ.
Las reacciones e+ + e− → iγ también son posibles, pero la sección eficaz cae con un factor
1
i−2
q i−2 ∼ .
137
Para tres fotones, entonces, la reacción es 137 veces menos probable que la más usual con dos. Si
queremos 4 fotones, es 18769 veces menos probable!!!
Otras formás de aniquilación de pares e± son posibles:
e+ + e− → µ+ + µ− : que es una reacción electromagnética, no débil. La sección eficaz es

σe→µ ∼
= 87/s mb, donde s es el cuadrado de la energı́a en el centro de masas.
e+ + e− → ν + ν: con una sección eficaz σe→ν ∼
= 10−10 /(6π)s b. Aquı́, interactúa el bosón Z 0 .
e+ +e− → q+q → hadrones: esta es una interacción electromagnética. Su sección eficaz depende
del tipo de quark. Si queremos uu, dd y ss, la sección eficaz es el doble de la σe→µ .
62
7. Absorción
Los rayos γ creados por partı́culas relativistas en una fuente astrofı́sica pueden ser absorbidos por
campos de radiación o materia en la fuente misma o en su trayecto al observador. Los principales
mecanismos de absorción son la creación de pares en el campo Coulombiano de un núcleo y la creación
de pares por aniquilación de fotones:
γ + γ → e+ + e− .
A bajas energı́as, el efecto Compton directo puede ser relevante y a energı́as muy altas y en
presencia de campos magnéticos intensos los fotones γ pueden crear pares
γ + B → e+ + e− .
Supongamos que la intensidad original de los rayos γ es Iγ0 (Eγ ) y que se los inyecta en un medio
de densidad de blancos n. Luego de atravesar una distancia x, la intensidad será
Iγ (Eγ ) = Iγ0 (Eγ )e−τ ,
donde la cantidad τ = σnx es la profundidad óptica del medio y σ es la sección eficaz para interacción
de los rayos γ con los blancos. En general, la profundidad óptica es una integral de lı́nea que tiene
en cuenta el efecto acumulado de la absorción a lo largo de la visual.
El camino libre medio del fotón en un medio uniforme de densidad n es
1
λγ = .
σn
La probabilidad de que el fotón sea absorbido luego de haber atravesado una distancia L es
!
L
PL = 1 − exp − .
λγ
7.1. Procesos de absorción de energı́a

7.1.1. Creación de pares en un campo Coulombiano
El umbral de energı́a para que un rayo γ cree un par e± en el campo electrostático de un nucleón
N de carga eZ es 2me c2 ∼ 1,022 MeV.
La sección eficaz para el proceso es
7 2Eγ 109 2Ee+ Ee− me c2

σγN →e± = 4α Z 2 re2 ln − para ≪ ,
9 me c2 54 Eγ α Z2
7 183 1 2Ee+ Ee− me c2

σγN →e± = 4α Z 2 re2 ln − para ≫ .
9 Z 1/3 54 Eγ α Z2
Aquı́, Ee± es la energı́a de los leptones producidos. La primera ecuación se usa para núcleos sin
apantallamiento, mientras que la segunda, para casos con apantallamiento.
Notemos que como re = 2,8 × 10−13 cm, la sección eficaz es pequeña y el camino libre medio
en el medio interestelar (donde tı́picamente nMIE ∼ 0,01 − 1 cm−3 ) es extremadamente largo. Esto
significa que la Galaxia es esencialmente transparente a los rayos γ.
63
Figura 17: Sección eficaz total para la creación de pares por aniquilación de dos fotones en un campo
de radiación blanco isotrópico.
7.1.2. Absorción por creación de pares en un campo de radiación

Un fotón de energı́a Eγ puede producir un par e± en una interacción con un fotón de energı́a Eph
si
Eγ Eph (1 − cos θ) > 2(me c2 )2 ,
donde θ es el ángulo entre las direcciones de movimiento de los fotones.
La sección eficaz total para la aniquilación γ + γ → e+ + e− esta dada por
3 1+β 2m2e c4
" !#

σγγ (Eγ , Eph ) = σT 1 − β 2 2β β 2 − 2 + 3 − β 4 ln , 1 − β2 = ,
16 1−β Eγ Eph (1 − cos θ)
para 0 ≤ β < 1.
Si asumimos isotropı́a en la distribución de los fotones blancos, podemos aproximar la sección efi-
caz σγγ (s0 ), donde el parámetro s0 = (Eγ Eph )/(me c2 )2 determina la energı́a de los leptones inyectados
en el sistema de referencia del centro de masas:
3σT 1 1 1 4 1 1
s
√
" #
√
σγγ (s0 ) = 2 s0 + ln s0 − + ln s0 + s0 − 1 − s0 + − 1−
2s0 2 6 2s0 9 9s0 s0
En la Figura 17 se muestra la sección eficaz total para una distribución de fotones blanco isotrópica
en función del parámetro s0 . En el limite s0 → 1, la sección eficaz se anula, mientras que es similar
a la del efecto Compton inverso para s0 ≫ 1. Para una dada energı́a Eγ de los fotones γ, la energı́a
umbral de los fotones blanco es Eth = (me c2 )2 /Eγ . La sección eficaz de la interacción tiene un máximo
pronunciado σγγ ≈ 0,2σT para s0 ≳ 3; en consecuencia, fotones con una dada energı́a Eγ son más
eficientemente absorbidos por fotones blanco con energı́as cercanas a la energı́a umbral en un rango
angosto alrededor de Eph ≳ eEth .
64
Si un fotón de energı́a Eγ debe atravesar una región del tamaño R con un campo de fotones
isotrópico de densidad nph , la profundidad óptica será
ˆ ∞ ˆ R
τγ (Eγ ) = nph (Eph , r)σγγ (Eγ , Eph )dEph dr.
Eth (Eγ ) 0
En medios astrofı́sicos nph puede ser muy alta en cercanı́a de fuentes, por lo que la absorción en
campos de radiación es un fenómeno común. Por otro lado, la densidad de fotones media interestelar
de nuestra galaxia es demasiado baja como para que la opacidad sea significativa. En cambio, sobre
escalas cosmológicas, la opacidad del fondo cósmico de radiación a la propagación de rayos γ puede
ser importante.
En el caso de tener una fuente de rayos γ con una luminosidad intrı́nseca Lγ , la atenuacion por
γ + γ → e+ + e− intrı́nseca será determinada por la cantidad de fotones absorbente y el tamaño de
la fuente. La mayor eficiencia de absorción γγ se da para fotones blanco cuya energı́a es cercana a la
umbral de la interacción, y por lo tanto inversamente proporcional a la de los rayos γ
(me c2 )2
Eph ∼ .
Eγ
El parámetro de compacticidad l para una fuente se define como

Lph hi
l≡ , l = erg s−1 cm−1 ,
R
donde R es el tamaño tı́pico y Lph , la luminosidad del campo de radiación absorbente de la fuente.
Aproximando Lph ≈ 4πcnph R2 Eph es posible estimar la opacidad τγγ según
σγγ
τγγ (Eγ ) ≈ σγγ nph R ≈ l.
4πcEph
La región central de un AGN puede tener luminosidades isotrópicas de Lγ ∼ 1048 erg s−1 y
dimensiones R ∼ 1015 cm, lo que da una compacticidad l ∼ 1033 erg s−1 cm−1 . En el caso particular
en que Eγ ∼ 1 MeV ∼ Eph y usando que cerca del umbral σγγ ∼ = σT /2, resulta τγγ ∼ 102 . Por
consiguiente, fuentes como esta pueden estas fuertemente auto-absorbidas a menos que la emisión
no sea isotrópica. De hecho, el que se observe emisión γ indica que debe existir una importante
anisotropı́a o beaming en cual es consistente con la idea de que la radiación γ se origina en jets -
chorros de partı́culas eyectadas por la fuente central.
7.1.3. Absorción en campos magnéticos

Un fotón γ puede ser convertido en un par e± en presencia de un campo magnético
γ + B → e+ + e− .
La razón a la que se produce esta transformación es muy pequeña a menos que
Eγ B m2e c3

ξ= ∼ 0,1, Bcr = ∼ 4,4 × 1013 G.
2me c2 Bcr eℏ
Se necesitan, pues, campos magnéticos muy fuertes para que este fenómeno domine la atenuación
de los rayos γ. Estos campos se dan en los púlsares.
65
El camino libre medio de un fotón con energı́a Eγ > me c2 en un campo B es
44 4
!
ℏ Bcr
λγ = exp ,
10(e ℏ) me c B sen(θ)
2 3ξ
donde θ es el ángulo entre la dirección de propagación del fotón y B.

⃗ Para θ = 0 no hay producción
de pares. Como la creación magnética de pares es muy sensible a la intensidad del campo magnético,
el criterio para que ocurra es
B sen(θ) 1
!
Eγ
≥ .
2me c 2 Bcr 15
7.1.4. Interacción Compton directa

La interacción Compton directa puede ser una causa importante de absorción para fotones de
energı́as ∼ 1 MeV. Al colisionar el fotón con el electrón, este gana energı́a cinética que le entrega el
fotón.
La sección eficaz para la interacción Compton directa depende de la polarización de los fotones.
Si el vector eléctrico de los fotones incidentes forma un ángulo Θ con el de los fotones emergentes, la
sección eficaz puede expresarse como (en unidades de cm2 electrón−1 )
!2
Eγ1 Eγ0 Eγ1
!
r2
dσC Eγ1 , Eγ0 , Θ = e + − 2 + 4 cos2 Θ dΩ,
4 Eγ0 Eγ1 Eγ0
donde Eγ1 es la energı́a del fotón dispersado y Eγ0 , la energı́a original. Si el fotón incidente no esta
polarizado,
!2
re2 Eγ1 Eγ0 Eγ1
!

dσC Eγ , Eγ , Θ =
1 0
+ − sen Θ dΩ.
2
4 Eγ0 Eγ1 Eγ0
A bajas energı́as, e integrando sobre dΩ, esta ultima expresión se transforma en la sección eficaz
Thomson
8π 2
σT = r .
3 e
7.1.5. Debilitamiento de rayos γ por efectos Doppler y gravitacional

Los rayos γ pueden ser afectados por el estado de movimiento de la fuente o por la presencia de
campos gravitacionales intensos en la fuente.
Si la fuente se mueve hacia el detector a una velocidad cβ formando un ángulo θ con la visual,
entonces la energı́a de la radiación emitida se ve modificada por efecto Doppler
h i−1
Eγobs = δEγ0 , δ = γ 1 − β cos θ ,
con δ el factor Doppler, Eγ0 la energı́a del fotón en el sistema de referencia de la fuente y Eγobs
la energı́a en el sistema de referencia del observador. Cuando la fuente se aleja del observador,
δ[γ(1 + β| cos θ|)]−1 < 1 y la energı́a observada resulta menor que la emitida.
La expansión cosmológica del Universo también modifica la energı́a de los fotones. El corrimiento
al rojo cosmológico se define como
∆λ λ − λ0
z≡ = ,
λ0 λ0
66
donde λ es la longitud de onda de los fotones recibidos y λ0 es la que tenı́an al ser emitidos. De aquı́
que la variación de energı́a del fotón debida a la expansión es
Eγ0
Eγobs = .
1+z
El corrimiento al rojo puede expresarse en función de los factores de escala del Universo al tiempo
actual y cuando el fotón fue emitido
R(tactual )
1+z = ,
R(temitido )
siendo la diferencia R(tactual ) − R(temitido ) el tiempo en que la radiación viajo hasta el observador. El
factor de escala se determina a través de un modelo cosmológico relativista.
Si la fuente de los fotones esta sometida a un campo gravitacional fuerte, la energı́a de un fotón
emitido en la fuente sufrirá un corrimiento al rojo de origen gravitacional. Si la fuente tiene una masa
M y un radio R, la variación de la energı́a de un fotón emitido con energı́a Eγ0 será
GM 0
∆Eγ = E ,
c2 R γ
de tal forma que Eγobs = Eγ0 − ∆Eγ .
7.2. Cascadas electromagnéticas

La inyección de rayos γ de alta energı́a formados por materia (plasma, gas, solidos), radiación
o campos magnéticos resulta en la formación de pares e± si la profundidad óptica a la propagación
de los rayos γ es τ ≳ 1. Si estos pares tienen suficiente energı́a, pueden producir, a su vez. más
radiación γ por Bremsstrahlung relativista, efecto Compton inverso y/o radiación sincrotrón. En estas
condiciones, una cascada electromagnética se desarrolla en el medio. El resultado de esta cascada
es degradar la energı́a de los fotones originales y multiplicar el número de leptones. El espectro
emergente depende del espectro original de inyección y de las caracterı́sticas del medio.
Una vez iniciada, la cascada se desarrollará hasta que las escalas temporales de los diferentes
procesos radiativos en competición que producen fotones fuera del rango γ sea menor que la de
los procesos que resultan en rayos γ. La cascada también se detendrá si el tiempo de enfriamiento
radiativo de las partı́culas excede el tiempo de las pérdidas no radiativas o si la profundidad óptica
a la propagación de los rayos γ cae por debajo de 1.
Cascadas electromagnéticas se desarrollan en numerosas situaciones astrofı́sicas. Por ejemplo, la
inyección de rayos γ en binarias de rayos X o AGNs puede iniciar cascadas. En la magnetósfera
de los púlsares, se espera el desarrollo de cascadas en los gaps electrostáticos donde los electrones
son acelerados. Los rayos γ originados en fuentes extragalácticas distantes pueden inicial cascadas
electromagnéticas en el fondo cósmico de radiación.
A energı́as para las cuales τ > 1, el efecto sobre el espectro inyectado es una disminución en
un factor exp(−τeff ), donde la opacidad efectiva es menor que la original. Dependiendo de la dureza
(pendiente) del espectro inyectado, puede ser más o menos notable la acumulación de fotones con
energı́as por debajo de la condición τ = 1.
67
Figura 18: Efecto de las cascadas electromagnéticas en el espectro de un microcuásar. La curva rosa
punteada muestra el espectro intrı́nseco y las otras curvas, el espectro reprocesado. El efecto de las
cascadas varı́a con la fase orbital, siendo mas relevante cuando el objeto compacto esta detrás de la
estrella (curva solida).
7.2.1. Cascadas electromagnéticas en la materia

Consideremos una cascada electromagnética en un medio donde en Bremsstrahlung relativista es
el principal mecanismo de enfriamiento. Supongamos que un fotón de energı́a Eγ ≫ me c2 incida en
tal medio. Su camino libre medio será
1
λγ ∼ ,
σγp n
donde σγp es la sección eficaz para la creación de un par e± en el medio y n es su densidad.
El par creado a una profundidad λγ = R dentro del medio tendrá, a su vez, un camino libre medio
de
1
λe± ∼ ,
σB n
donde σB es la sección eficaz para Bremsstrahlung relativista. La energı́a de cada elemento del par
(1)
e± es Ee± ∼ E0 /2, donde E0 = Eγ . Después de una distancia λB (que es ∼ λγ = R) cada partı́cula
radia un nuevo fotón con energı́a
E0
Eγ(1) = .
4
Estos fotones, al cabo de una distancia R, crearán nuevos pares. A medida que la cascada se
desarrolla, el número total de fotones y leptones aumenta pero su energı́a media disminuye. A una
profundidad L = z/R, el número de partı́culas será
N (L) ∼ eL , Eγ ∼ E0 e−L ,
con Eγ la energı́a por partı́cula a esa profundidad. Es decir, que la energı́a decae aproximadamente
en forma exponencial con la profundidad.
La cascada cesa cuando se alcanza una energı́a
E0

Ec ∼ E0 e−Lmax
, Lmax ∼ ln ,
Ec
68
donde Lmax es la distancia a la cual el Bremsstrahlung deja de dominar las pérdidas. El número de
fotones a esa profundidad es 1/3 del número total de partı́culas
2 E0 2 E0

Ne± Lmax ∼ exp ln = , Ntotal ∼ e−Lmax .
3 Ec 3 Ec
7.3. Cascadas hadrónicas

Cuando un protón de muy alta energı́a impacta en un núcleo, interacciona con algún nucleón
produciendo piones. También pueden producirse partı́culas extrañas y antinucleones. Si la energı́a
es suficiente, cada nuevo hadrón puede también interaccionar dentro del núcleo dando lugar a una
mini-cascada hadrónica. Los nucleones que interaccionan con el protón primario por lo general son
removidos del núcleo dejándolo en un estado que puede ser inestable; como consecuencia el núcleo se
fragmenta y se emiten núcleos más livianos en un proceso conocido como spallion. Los núcleos ligeros
inyectados se suelen llamar astillas o fragmentos de astillado. Estos fragmentos son emitidos más o
menos isotrópicamente en el sistema de laboratorio. En cambio, los productos de la cascada hadrónica
tienen una componente importante de su momento en la dirección de movimiento del protón original.
Normalmente, neutrones son eyectados tanto por el núcleo original como por los fragmentos.
Los piones producto de la cascada hadrónica decaen según
π 0 → 2γ o π + → µ+ + νµ o π + → µ− + ν µ .
Los muones, a su vez, decaen como
µ+ → e+ νe + ν µ , µ− → e− + ν e + νµ .
Tanto los rayos γ como los e± producidos en estos decaimientos dan lugar a cascadas electro-
magnéticas.
Los muones tienen vidas medias relativamente largas, por lo que antes de decaer pueden perder
energı́a por los mismos mecanismos radiativos que afectan a los electrones y positrones. Los muones
resultado de una cascada iniciada por un rayo cósmico en la atmósfera pueden alcanzar la superficie
terrestre. Los detectores de rayos cósmicos, precisamente, operan detectando radiación Cherenkov
que producen en la atmósfera o la que generan en tanques cerrados llenos de agua.
8. Detectores
Los problemas y desafı́os que plantea la detección y medición de rayos γ de origen cósmico son
formidables. De hecho, la situación de la astronomı́a a estas energı́as es única, en el sentido de que
solo en esta banda el ruido de fondo es mucho mayor que las señales a detectar.
Para cuantificar este problema, observemos que el flujo medio de rayos X de 1 keV incidente sobre
un detector en orbita es de ∼ 15 cm−2 s−1 , lo que es unas 100 veces más que el flujo de rayos cósmicos
en las mismas condiciones. En cambio, el flujo medio de rayos γ con Eγ > 100 MeV es de ∼ 2 × 10−4
cm−2 s−1 sr−1 - unas mil veces mejor que el correspondiente flujo de rayos cósmicos.
Al interaccionar con un detector, los rayos cósmicos producen partı́culas secundarias que a su vez
producen rayos γ localmente. La eliminación de estas contribuciones locales al ”fondo”de radiación
ha sido y es el principal problema técnico en la astronomı́a de rayos γ.
Dado que los rayos γ interaccionan con la materia por diferentes procesos con distinta probabilidad
según su energı́a, las técnicas de detección varı́an con la energı́a de los fotones a detectar. La relativa
importante de las distintas formas de interacción para los rayos γ se ven en la Figura 19.
69
Figura 19: Importancia relativa de los diferentes mecanismos de perdida de energıa de un rayo γ en
funcion de la energıa del foton y del numero atomico del material blanco.
Los distintos instrumentos harán uso de estos efectos a las distintas energı́as. A energı́as muy
altas (E ≥ 100 MeV), los rayos γ producen cascadas en la atmósfera y sus caracterı́sticas pueden ser
inferidas a partir de la reconstrucción de las cascadas.
8.1. Astronomı́a gamma desde tierra: telescopios Cherenkov

Cuando un rayo γ ingresa a la atmósfera inicia
una cascada electromagnética (extended air shower,
EAS). Si Eγ > 100 GeV, la cascada se desarrolla
hasta una altitud de unos pocos km sobre el ni-
vel del mar (∼ 10 km para Eγ ≈ 100 GeV-1 TeV)
donde se alcanza el número máximo de partı́culas
secundarias, luego se desvanece más profundo en la
atmósfera. Los leptones secundarios son relativistas
y producen luz Cherenkov en el UV (con un pico en
∼ 300 − 350 nm a nivel del suelo) al propagarse en
la atmósfera. La emisión Cherenkov se produce en
forma de un pulso de duración de algunos ns y esta
contenida en un cono cuyo eje coincide con la direc-
ción de arribo del rayo γ. El área iluminada en tierra
(light pool) tiene unos 250 m de diámetro, aproxi-
madamente, y la densidad de fotones que arriban
durante el pulso es de ∼ 100 cm−2 a 1 TeV.
La luz Cherenkov es colectada en telescopios for-
mados por arreglos de espejos que la enfocan en una
cámara que opera como detector. Las caracterı́sti-
cas de la imagen obtenida permiten determinar la dirección de arribo y energı́a del rayo γ. La forma
de la misma permite rechazar eventos que no se deben a cascadas iniciadas por fotones. En particular,
las imágenes debidas a cascadas electromagnéticas son elipses angostas y elongadas en el plano de
la cámara. El eje mayor de la elipse corresponde a la extensión vertical de la cascada y apunta en
la dirección de la fuente en el campo visual del telescopio. Una sola imagen no permite reconstruir
70
apropiadamente la geometrı́a espacial de la cascada, por lo que se utilizar arreglos de telescopios que
proveen una visión estereoscópica: si varios telescopios detectan la misma cascada, la posición de la
fuente queda determinada por la intersección de los ejes mayores de todas las imágenes.
La energı́a total emitida durante la cascada en forma de luz Cherenkov es proporcional a la energı́a
del rayo γ primario. La intensidad de la imagen de la cascada combinada con el dato de la distancia
de cada telescopio al punto de inicio de la cascada permiten estimar la energı́a del rayo γ. El rango
de energı́as potencialmente detectables es de ∼ 50 GeV a ∼ 100 Tev, dependiendo del instrumento,
el tiempo de exposición y la fuente observada. La resolución en energı́a tı́pica es de 15 − 20 % para
Eγ ≳ 100 GeV.
El detector consiste en un espejo que colecta y refleja la luz Cherenkov sobre un foco donde hay
un tubo fotomultiplicador (PMT). Si el tiempo de integración del PMT es mayor que la duración del
destello Cherenkov, la señal detectada es
ˆ E2
C= C(E)ζ(E)AdE,
E1
donde C(E) es le flujo Cherenkov dentro de los lı́mites de sensibilidad en energı́a E1 y E2 del PMT,
ζ(E) es la curva de respuesta de PMT y A es el área colectora. A su vez, C(E) depende de la trans-
misión atmosférica T (E) y del espectro de emisividad de la cascada ϵ(E): C(E) = κϵ(E)T (E), κ =
cte.
La señal Cherenkov se detecta sobre un ruido producido por el cielo nocturno. El brillo del cielo
es ˆ E2
B= B(E)ζ(E)τ AΩdE,
E1
donde τ es el intervalo de integración, B(E) la dependencia del brillo con la energı́a de los fotones y
Ω el ángulo solido cubierto por el telescopio.
La relación señal-ruido es entonces
ˆ E2 ( )1/2
S S ζ(E)A
=√ = C(E) dE.
N B E1 B(E)Ωτ
La mı́nima señal detectable es ∝ (S/N )−1 . De aquı́ que la energı́a mı́nima detectable sea
)1/2
1
(
B(E)Ωτ
Eγmin ∝
ζ(E)A C(E)
El brillo del cielo esta constituido por la luz de las estrellas, el Sol, la luna, el brillo nocturno,
meteoros, rayos, nubes y fuentes de origen artificial. Para minimizar estos efectos, se observa en sitios
aislados durante noches sin luna.
El principal problema para las observaciones Cherenkov, sin embargo, es el ruido de fondo produ-
cido por los rayos cósmicos. Las cascadas hadrónicas iniciadas en la atmósfera por los rayos cósmicos
tienen componentes electromagnéticas que emiten luz Cherenkov.
Hay diferencias ente las cascadas iniciadas por rayos γ y aquellas iniciadas por rayos cósmicos
que permiten diferenciarlas. La tarea no es fácil: las cascadas iniciadas por rayos cósmicos, a una
energı́a dada, son 103 veces más numerosas que las iniciadas por rayos γ.
71
La diferenciación efectiva de ambos tipos de
eventos se hace sobre la base de una técnica cono-
cida como Cherenkov imaging, por lo que este tipo
de detectores se conoce como Imaging Atmospheric
Cherenkov Telescopes. La técnica consiste en pixe-
lar la imagen de la lluvia Cherenkov utilizando un
arreglo de muchos PMT que forman una cámara
colocada en el foco del telescopio. Al obtenerse una
imagen de la lluvia Cherenkov es posible descartar
eventos de fondo en base a la forma y orientación de
Figura 20: Imágenes en la cámara de un telesco-
la imagen. Como dijimos, las imagenes de cascadas
pio Cherenkov de una cascada iniciada por un
iniciadas por rayos γ son elipses angostas y elon-
rayo γ (izquierda) y un protón (derecha). La es-
gadas. Las cascadas hadrónicas iniciadas por rayos
cala de colores corresponde a la intensidad de
cósmicos producen partı́culas secundarias cargadas
la luz Cherenkov detectada en cada pı́xel.
con una importante componente transversal de su
momento, que a su vez disparan sub-cascadas. Co-
mo consecuencia, las imágenes de las cascadas iniciadas por rayos cósmicos son más anchas e irre-
gulares que las de las cascadas iniciadas por rayos γ. Para fuentes puntuales, como AGNs distantes,
los modernos telescopios tipo IACT pueden rechazar hasta el 99,999 % de los eventos disparados por
rayos cósmicos. El rechazo es menos eficiente para fuentes extendidas, pero aun se puede aplicar la
técnica par obtener mapas de emisión en rayos γ. La resolución angular depende de la energı́a, pero
es tı́picamente de 0,1◦ para Eγ ≳ 100GeV .
Los detectores de tipo IATC de generacion actual y futura incluyen los siguientes instrumentos:
H.E.S.S (High Energy Stereoscopic System): 4 telescopios de 12 m de diámetro dispuestos
en un cuadrado y uno de 28 m de diametro en el centro, que operan en modo estereoscopico.
Puede detectar rayos γ con energı́as por encima de ∼ 50 GeV con una resolución angular de
≲ 0,1◦ .
VERITAS (Very Energetic Radiation Imaging Telescope Array System): Es un arre-

glo de 4 telescopios de 12 m de diametro separados unos 100 m entre ellos, cada uno formado
por 350 espejos individuales. Detecta rayos γ con sensibilidad máxima en el intervalo 0,1 − 10
TeV, con una resolución en energı́a del 17 % y una resolución angular de 0,08◦ a 1 TeV.
MAGIC (Major Atmospheric Gamma Imaging Cherenkov): consta de dos telescopios

de 17 m de diámetro que observan en modo estereoscópico. Detecta rayos γ con energı́as
entre∼ 50 GeV y ∼ 20 TeV. Para fuentes puntuales Eγ ≳ 200 GeV, la resolución angular es de
0,07◦ y la resolución en energı́a del 16 %.
CTA (Cherenkov Telescope Array): El arreglo norte será mas pequeno y esta pensado
especialmente para detectar fotones entre 20 GeV y 20 TeV (rango de energıas bajo a medio),
mientras que el arreglo sur operara en el rango completo de energıas entre 20 GeV y 300 TeV. Se
espera que CTA tenga una sensitividad 10 veces mejor que los instrumentos actuales. Debido
a su gran area colectora y ubicación en los dos hemisferios, se podra observar practicamente
todo el cielo, con una tasa de deteccion 10 veces mayor a la actual.
Además de los telescopios Cherenkov basados en cámaras y en la técnica de imágenes, existe otro
tipo de instrumentos donde lo que se detecta es la luz Cherenkov emitida por las partı́culas secun-
darias cargadas creadas en la cascada, dentro de enormes tanques llenos de agua. Eventualmente,
72
también los rayos γ pueden interactuar dentro del tanque para producir pares electrón/positrón que
luego emiten luz Cherenkov.
La emisión Cherenkov en agua es eficiente debido al alto ı́ndice de refracción del medio. Como
vimos, la luz Cherenkov se emite dentro de un cono en la dirección de movimiento de la partı́cula
relativista; el ángulo de apertura del cono depende del ı́ndice de refracción del medio. En aire (naire =
1,0003), el ángulo de apertura del cono de emisión es de ∼ 1◦ , mientras que en agua (nagua = 1,33) es
de ∼ 41◦ . Que el cono de emisión sea tan grande en agua permite detectar la emisión de prácticamente
todas las partı́culas que entran al tanque, mediante alguno de los varios PMTs instalados en su
interior. Uno de los observatorios de rayos γ que opera basado en esta tecnica es HAWC (High-
Altitude Water Cherenkov Gamma-Ray Observatory).
8.2. Astronomı́a gamma espacial

8.2.1. 30 MeV ≤ Eγ ≤ 300 Gev
En este rango de energı́a se utilizar detectores en orbita basados en la creación de pares electrón-
positrón en interacciones de rayos γ con materia. El detector es el telescopio mismo y la superficie
del detector es el área colectora.
El principal tipo de telescopio por produccion de pares es la camara de chispas (spark chamber).
Sus elementos basicos son
El rastreador (tracker): determina la dirección de llegada del rayo γ a partir de las tra-
yectorias de un par creado por el. En general, consiste en una serie de placas metálicas (lo
suficientemente gruesas para que absorba el rayo γ y cree el par pero no tanto como para que
afecte su trayectoria) colocadas en paralelo en una cámara sellada y llena de gas a presión, de
modo que el paso de un protón o electrón deje una estela de ionización. Las placas están conec-
tadas de forma tal que cuando una partı́cula cargada se mueve en la cámara hay una diferencia
de potencial entre placas sucesivas. Al pasar la partı́cula cargada e ionizar el gas (generalmente
mezcla de neón y etano), se produce una chispa eléctrica entre placas. Las posiciones en las
que se producen las chispas permiten determinar el movimiento de la partı́cula.
73
El disparados (trigger): activa la diferencia de potencial entre las placas cuando atravie-
sa una partı́cula cargada. Su eficacia es la que básicamente determina el campo visual del
instrumento.
El calorı́metro: absorbe los e± para medir su energı́a (por lo tanto, la del rayo γ). Su espesor
debe ser de varios caminos libres medios del electrón en el material usado.
El escudo anti-coincidencia (anticoincidence shield): rodea al receptor e indica el paso

de partı́culas cargadas. Tiene una sección eficaz muy baja para interacciones con rayos γ. Gene-
ralmente es un centelleador plástico monitoreado por fotomultiplicadores; si estos se disparan
(indicando el paso de un rayo cósmico) no se aplica la diferencia de potencial a las placas y
el evento no se registra. Esto permite eliminar el ruido de fondo producido por rayos cósmicos
cargados (aunque no por neutrones).
El instrumento más importante que uso una cámara de destellos fue el EGRET (Energetic
Gamma-Ray Experiment Telescope).
Los telescopios actuales no usan cámaras con gas sino detectores de silicio, que son semiconduc-
tores que permiten medir el punto de interacción de una partı́cula cargada con gran precisión. La
interacción con la partı́cula cargada genera corriente en el semiconductor y permite una inmediata
digitalización de la lectura del instrumento. Observatorios satélites que incorporan esta tecnologia
son Fermi (con el instrumendo LAT) y AGILE (Astro-Rivelatore Gamma a Immagini Leggero).
8.2.2. 1 MeV ≤ Eγ ≤ 30 MeV

El coeficiente de absorción de un rayo γ en un medio
material alcanza un máximo en el rango 1 a 10 MeV. El
mecanismo de interacción dominante es el de la interacción
Compton directa: un rayo γ es dispersado por un electrón
y le transfiere parte de su energı́a a este. El rayo γ original
pierde energı́a y altera su dirección de propagación. La
energı́a E1e transmitida al electrón y la energı́a Eγ del fotón
incidente se relacionan con el ángulo de dispersión del fotón
como
Eγ2 (1 − cos θ)
E1e = .
Eγ (1 − cos θ) + m3 c2
Un esquema de un detector consiste de dos planos de
detección P1 y P2 . En el primer plano hay un detector dis-
persor (donde el rayo γ interacciona por efecto Compton)
y en el segundo un detector absorbente (donde el fotón dis-
persado es absorbido). Las posiciones de las interacciones y la pérdida de energı́a del fotón en cada
plano detector se miden para estimar la energı́a y dirección de arribo del rayo γ. Ambos detectores
están separados por una distancia d ∼ 2 m. El detector P1 esta hecho de un material de bajo Z
mientras que P2 esta hecho con materiales de alto Z a fin de facilitar la absorción. La determinación
de E1e y E2e permite obtener Eγ y θ, pero no es suficiente para determinar la dirección de llegada del
rayo γ original, sino que hay todo un cono de posibilidades (no muy buena resolución angular).
74
8.3. Eγ ≤ 1 MeV
A bajas energı́as el efecto fotoeléctrico es la interacción dominante de fotones con materia. Para
detectarlos se utilizan distintos tipos de centelladores en los telescopios. Un centellador consiste de
un material que convierte en luz parte de la energı́a pérdida por ionización por una partı́cula cargada.
Se suele acoplan con un instrumento fotoeléctrico que convierte la luz en una señal eléctrica.
El principal problema de estos detectores es la contaminación de la señal por eventos disparados
por rayos cósmicos. Para aliviarlo, usamos escudos activos formados por centelladores.
La dirección de arribo de los rayos γ suele determinarse por el uso de colimadores, opacos a los
rayos γ y su geometrı́a permite determinar un rago de ángulos de llegada para los rayos γ.
8.4. Emisión difusa y detección de fuentes puntuales

La interacción de los rayos cósmicos son el medio interestelar de la Galaxia produce emisión γ
difusa. Toda fuente discreta debe ser detectada sobre este fondo. Se aplican técnicas básicas para
remover la contaminación del fondo difuso:
Tecnica on/off: consiste en observar alternativamente en dirección a la supuesta fuente y

luego a una región cercana libre de fuentes. El resultado de esta ultima observación se sustrae
de la primera a fin de determinar el flujo real de la fuente. El problema de esta técnica es que la
radiación de fondo dista de ser uniforme y una mala elección de la región off puede introducir
errores considerables.
Modelo teórico del fondo que se sustrae a las observaciones: el modelo tiene como
entradas las distribuciones de gas y polvo en la Galaxia, ası́ como campos de radiación. Se
calcula, usando todos los procesos radiativos relevantes, la distribución de emisión γ resultante,
función de la energı́a de los fotones, longitud y latitud galáctica. El modelo de emisión difusa
se sustrae sistemáticamente de todas las observaciones realizadas con el instrumento a fin de
identificar fuentes discretas. Si el modelo es incompleto entonces se pueden generar fuentes
espúreas.
Técnica de ocultación: se utiliza un ocultador hecho de algún material pesado para tapar
la fuente del campo visual para medir los fotones que vienen del fondo. Luego, la fuente se
descubre y se vuelve a medir. El espectro de la fuente se obtiene al sustraer ambos espectros.
De los cambios en la modulación (si el ocultador se mueve hacia delante y hacia atrás) pueden
inferirse cambios en el fondo - técnica conocida como ”chopper technique”.
8.5. Astronomı́a de neutrinos

Los neutrinos son partı́culas elusivas debido a su pequeña sección eficaz de interacción con la
materia. Esa misma caracterı́stica les da importancia astrofı́sica, ya que escapan de las diferentes
fuentes sin ser afectados por problemas de absorción o pérdidas. Los neutrinos son partı́culas muy
ligeras, eléctricamente neutras y solo interactúan débilmente.
Los antineutrinos electrónicos producidos en reactores nucleares con energı́as de ∼ 1 MeV tienen
una sección eficaz de solo ∼ 10−44 cm2 . Esto significa que solo un neutrino entre 1011 interactuara
cuando viaje a lo largo del diámetro de la Tierra. Aunque la energı́a sea mayor, la sección eficaz sigue
siendo pequeñı́sima. Esto implica que los detectores de neutrinos sean gigantescos para lograr tasas
de eventos apreciables.
75
El principio básico de un telescopio de neutrinos
es el de una matriz formada por detectores ópticos
en agua o hielo, lo que ofrece un gran volumen ab-
sorbente para los neutrinos, mientras que al mis-
mo tiempo proporciona protección contra los rayos
cósmicos secundarios. Las partı́culas cargadas rela-
tivistas producidas por la interacción de los neutri-
nos emiten radiación Cherenkov en el medio trans-
parente. Un detector debe medir con gran precisión
el número y el tiempo de llegada de estos fotones en
una matriz tridimensional de PMT, cuya informa-
ción permite reconstruir algunas de las propiedades
del neutrino original.
Figura 21: Sección eficaz para neutrinos y anti-
Los neutrinos de altas energı́as solo pueden ser
neutrinos muónicos según una función particu-
producidos por procesos hadrónicos (decaimiento de
lar de distribución de patrones.
π ± u otros mesones). Esos mesones cargados pueden
ser el resultado de la interacción de protones u otros
núcleos con campos de materia en el entorno de las fuentes o con campos de radiación. En esas
interacciones se produce prácticamente el mismo número de cada especie de pión. El π 0 decae en dos
γ mientras que los π ± → µ± + ν µ (νµ ) y luego µ± → e± + ν e (νe ), por lo tanto hay tres neutrinos por
cada pión y seis neutrinos por cada dos rayos γ.
La relación entre la energı́a del neutrino producido en estas reacciones y el protón original es
Eν ≈ 0,05E, mientras que la relación entre las luminosidades es Lν /Lγ = 1/3.
8.5.1. Telescopios de neutrinos

Los telescopios deben ser capaces de detectar neutrinos individuales, determinar su energı́a y
direcciones de llegada. Esta es una tarea formidable no solo por la pequeña sección eficaz del neu-
trino, sino también por el enorme ruido producido en los instrumentos por los muones y neutrinos
atmosféricos que son creados por los rayos cósmicos que bombardean la Tierra.
El principio básico de detección es que un neutrino interactúa con un átomo en una gran masa
de agua o hielo produciendo un muón. El muón es una partı́cula cargada muy penetrante que se
mueve por ese medio más rápido que la luz, produciendo luz Cherenkov en el rango óptico. Esa
luz es captada por arreglos 3D de tubos fotomultiplicadores permitiendo reconstruir la dirección de
propagación y energı́a del neutrino original.
Los neutrinos de alta energı́a interactúan con un nucleón N de un núcleo atómico a través de
interacciones débiles de corriente cargada (CC), con i = e, µ, τ
νi + N → i + X,
o interacciones débiles de corriente neutra (CN)
νi + N → νi + X,
donde X es el sistema hadrónico que se forma en la interacción y que se lleva parte de la energı́a del
neutrino incidente.
Las partı́culas cargadas producidas por la interacción de neutrinos viajan a través del medio hasta
que decaigan o interactúen. El camino libre medio depende de la energı́a de la partı́cula y las pérdidas
en el medio.
76
Los neutrinos de fuentes astrofı́sicas discretas deben detectarse sobre diferentes fondos difusos
debido a interacciones de rayos cósmicos con la atmósfera y fuentes no resueltas.
Además de los fondos debidos a neutrinos, existe una fuerte contaminación de los detectores por
muones creados en la atmósfera por los rayos cósmicos.
Los telescopios de neutrinos, al contrario de todos los demás instrumentos dedicados a la as-
tronomı́a, están ”mirando hacia abajo”, para evitar la contaminación del enorme flujo de muones
atmosféricos.
El número de eventos detectados sobre el intervalo de tiempo T es
ˆ
Nν dϕν eff
= dEν A (Eν ),
T dEν ν
donde dϕν /dEν es el flujo diferencial de neutrinos y Aeff

ν el area efectiva del detector, que se suele
medir en cm . Suponiendo un área de sección eficaz 1 km2 y una eficiencia de detección del 10 %,
2
el área efectiva a energı́as de 1 TeV corresponde aproximadamente a Aeffν (1 TeV)a ≈ 10 cm y el

3 2
número de eventos esperados es ∼ 1,5 por año. Por lo tanto, se requiere un detector del tamaño de
1 km3 para detectar algunas interacciones de neutrinos por año de una fuente galáctica, de acuerdo
con modelos de flujo de neutrinos estimados por las observaciones de rayos γ de TeV.
9. Fuentes de rayos γ
Llamamos fuente de rayos γ a un sistema astrofı́sico que emite una fracción significativa de
su luminosidad electromagnética a energı́as mayores que 0.5 MeV. Las fuentes de rayos γ pueden
clasificarse en 2 grandes grupos: pasivas y activas.
Las fuentes pasivas son simplemente ”blancos”para flujos de partı́culas relativistas originadas en
una región diferente. Las fuentes activas, por el contrario, aceleran partı́culas hasta velocidades rela-
tivistas y la interacción de estas partı́culas con los diferentes campos (materiales o electromagnéticos)
locales da lugar a la radiación γ.
9.1. Fuentes pasivas

Podemos separar en dos tipos: difusas o discretas.
Las fuentes difusas están formadas por un medio extendido y de baja densidad que es atravesado
por rayos cósmicos. Observacionalmente, la fuente difusa más importante es el medio interestelar
acumulado sobe el plano de la Galaxia. La irradiación de este material (básicamente H, más trazas
de elementos más pesados como He y material molecular como CO y polvo) produce emisión γ de
origen hadrónico a través de interacciones pp. Su intensidad dependerá de la densidad del medio y de
la densidad local de los rayos cósmicos. La emisión muestra un fuerte gradiente de intensidad hacia
bajas latitudes galácticas, donde se concentra la mayorı́a del gas. Esta emisión es la contribución más
importante al fondo difuso detectado por varios instrumentos y debe sustraerse a fin de detectar las
fuentes galácticas discretas.
Si escribimos ϵγ (⃗r) = qγ (⃗r)/n(⃗r) donde n es la densidad del medio y qγ (⃗r) la emisividad por
decaimientos de π 0 en la dirección ⃗r, el flujo total recibido a una distancia d será
ˆ
1
Fγ = n(⃗r)ϵγ (⃗r)d3 r,
4πd2
77
donde la integral se extiende sobre toda la región donde se distribuye el gas. Si la emisividad local
(en la vecindad del Sol) es ϵγ,0
ϵγ ωcr
∼ = κS ,
ϵγ,0 ωcr, 0
donde ωcr es la densidad de energı́a de los rayos cósmicos, que localmente vale ωcr, 0 ∼ 1 eV cm−3 .
Luego ˆ
1
Fγ ∼ κS ϵγ,0 n(⃗r)d3 r.
4πd2
De aquı́ que, si se conoce la distribución de densidad del medio y su distancia pueden hacerse
inferencias sobre la distribución de los rayos cósmicos de la Galaxia.
Sin embargo, la contribución del decaimiento de piones neutros no es la única a la emisión di-
fusa total. La interacción de la componente leptónica de los rayos cósmicos con el gas interestelar
también puede producir radiación Bremsstrahlung relativista. Además, la interacción de estas mis-
mas partı́culas con la radiación cósmica de fondo a 2.73 K puede resultar en la producción de rayos
γ. Entonces, deben desarrollarse modelos muy detallados para hacer una reproducción confiable de
la emisión difusa de la Galaxia, que deben contar con una descripción cuantitativa detallada de la
distribución de materia en la Galaxia y con modelos adecuados de los diferentes campos radiativos
(CMB, campos infrarrojos y luz estelar).
Las fuentes pasivas discretas están formadas básicamente por medio interestelar altamente es-
tructurado. Pueden ser nubes moleculares masivas, nubes de polvo o incluso nubes más pequeñas
ubicadas localmente o próximas a aceleradores de rayos cósmicos.
Si la densidad de una de tales nubes es constante, formada principalmente por H, podemos escribir
! !−2
Fγ Mcl d ωcr,0

∼ 10−9 κS ηA
cm s
−2 −1 1000M⊙ kpc eVcm−3
donde ηA es un factor que tiene en cuenta la presencia en el medio de elementos más pesados que el
H.
Vemos que si κS ∼ 1 la mayorı́a de las nubes moleculares no serı́an detectables para instrumentos
como EGRET a menos que o estén muy cerca o haya un acelerador cerca de manera tal que κS > 1.
Este ultimo caso se puede dar si la nube esta cerca de una fuente activa como un remanente
de supernova o un microcuasar, objetos capaces de producir rayos cósmicos que se difunden en el
medio. Es interesante notar que hay una acumulación de fuentes de rayos γ en regiones de formación
estelar conocidas como asociaciones OB. Aquı́ hay objetos jóvenes que pueden acelerar partı́culas
(púlsares, MQ, SNR) y nubes que pueden ofrecer blancos pasivos a las partı́culas relativistas que en
estas regiones tienen una densidad de energı́a mayor que la usual (κS ≫ 1).
Como los campos magnéticos en las nubes moleculares son mayores que en el medio interestelar,
una vez que los rayos cósmicos penetran en ellos pueden quedar atrapados allı́, con tiempos de
difusión muy largos, lo que aumenta la densidad de partı́culas relativistas en su interior.
9.2. Fuentes activas

Las fuentes activas de rayos γ son aquellas que son acelerados eficientes de partı́culas relativistas
y, al mismo tiempo, ofrecen campos adecuados para que estas partı́culas interaccionen y emitan
radiación γ. Podemos dividir a las fuentes activas en dos grandes clases: acretantes y no acretantes.
En las no acretantes galácticas podemos incluir a
púlsares aislados,
78
remanentes de supernova y
sistemas binarios de estrellas tempranas con vientos en colisión.
En las fuentes extragalacticas tenemos
galaxias eruptivas (starbursts),
cúmulos de galaxias y
galaxias normales.
Las galaxias normales son fuentes débiles. Emiten radiación difusa que debido a su distancia puede
aparecer como emisión discreta en el campo visual (como el caso de la Gran Nube de Magallanes,
detectada por EGRET). Como las galaxias normales son fuentes débiles, solo las más cercanas pueden
ser detectadas con la tecnologı́a actual.
Las starbursts son galaxias con episodios de formación estelar masiva. Están llenas de gas y se
espera que la densidad de rayos cósmicos sea alta.
9.2.1. Púlsares
⃗ ·B
Si Ω ⃗ < 0, el campo eléctrico sobre el polar cap del púlsar se dirige siempre hacia afuera de
la estrella. A diferencia de los leptones, los iones pueden estar firmemente aferrados a la superficie
formando una red cristalina. En cambio, los positrones pueden ser libremente acelerados. Estos se
mueven a lo largo de las lı́neas de campo sufriendo pérdidas por radiación de curvatura. Si el factor
de Lorentz de los positrones es γe , entonces
2 e2
2 1/2 1/2
dγe c R∗ P

= γe4 , Rc ∼ 10
4,9
dt 3 m e c2 Rc 1cm 1s
con Rc el radio de curvatura. Esta energı́a se emite en fotones γ que inician una cascada electro-
magnética a través de la reacción
γ + B → e+ + e− .
Los positrones son acelerados a su vez hacia fuera del púlsar, produciendo nuevos fotones γ y
más pares e± . Los electrones son acelerados hacia la superficie donde impactan produciendo un
calentamiento que resulta en la emisión de rayos X térmicos. La interacción de estos rayos X con
los fotones γ producidos por los electrones al moverse hacia el púlsar sobre las lı́neas de campo da
lugar a pares sobre la superficie, de donde salen los positrones originalmente acelerados. Los rayos γ
finalmente se escapan formando un cono con centro en el momento dipolar magnético. Además de
esta emisión no isotrópica, puede haber emisión de rayos γ en el viento del púlsar y en la llamada
nebulosa sincrotrón.
⃗ no es paralelo a µ
Como el púlsar esta en rotación, si Ω ⃗ en el sistema del observador, la emisión
aparece pulsada. El perı́odo de pulsación es igual al periodo de rotación P .
La cascada electromagnética produce un aumento del número de pares que son eyectados a lo
largo de las lı́neas de campo abierto formando el viento del púlsar. Una parte significativa de la
energı́a, sin embargo, escapa en forma de rayos γ.
79
El número de pares creados depende de las ca-
racterı́sticas del sistema y los parámetros particu-
lares asumidos, pero tı́picamente la multiplicidad
(número de secundarios producidos por cada leptón
primario inyectado) es κ ∼ 102 − 103 .
La luminosidad total es
! −2
m e c2 κ BS P

Ltotal ∼ 2,7×10 γ
33
erg s−1 ,
erg 103 1012 G 1s
donde para γ ∼ 1000, BS ∼ 1012 G y P ∼ 0,1 s se

tiene Ltotal ∼ 1032 erg s−1 .
La fracción que se emite en rayos γ la podemos
encontrar si consideramos que hay una dependencia
de γ con la distancia al centro de la estrella y que
las partı́culas escapan más allá del cilindro de luz
c γ0 − γfinal
γfinal = γ(Rcl ), Rcl = ∼ Ltotal .
Ωγ γ0
Esto resulta del orden de 1031 erg s−1 para púlsa-

res con P ∼ 0,1. Debido a la dependencia con P −2 ,
púlsares más rápidos pueden ser significativamente
más luminosos.
En la Figura 22 mostramos los pulsos de la emi-
sión a diferentes longitudes de onda. En rayos γ Figura 22: Existen tres regiones de radiacion
es más común una estructura doble del pulso que no termica asociadas con laemision por perdida
puede deberse a un efecto de apantallamiento de de energıa rotacional en una estrella de neutro-
la radiación dentro del cono de emisión en el polar nes. La primera region corresponde al pulsar y
cap o a efectos asociados a la presencia de un gap su magnetosfera, hasta el cilindro de luz. Es-
exterior. ta emision es periodica (pulsada). La segunda
La existencia de estas pulsaciones indica que la region corresponde al viento del pulsar, cuya
fuente es un objeto compacto en rotación. Sobre emision es debida principalmente a dispersiones
escalas de tiempo largas, sin embargo, los púlsares Compton inverso de los e relativistas. Por ulti-
no parecen ser variables. mo la nebulosa sincrotron, donde las partıculas
Los espectros de los púlsares de rayos γ mues- son re-aceleradas en el choque entre el viento y
tran diferentes contribuciones: emisión no termina el medio interestelar.
de origen coherente en radio, emisión de rayos X
térmica (proveniente de la superficie) y emisión no
térmica incoherente en rayos γ, aunque la mayor parte de la luminosidad es producida en rayos γ.
9.2.2. Remanentes de supernova

Pueden acelerar partı́culas cargadas hasta velocidades relativistas por medio del mecanismo de
Fermi. Es un proceso de aceleración difusivo que lleva a una ley de potencias en la distribución de
las partı́culas h i
n(E)dE = KE −Γ dE, n = cm−3 ,
con un flujo isotrópico I(E) = c/(4π)n(E).
80
Estas partı́culas (electrones o protones) pueden interaccionar con el material barrido por la ex-
plosión de supernova produciendo rayos γ.
Consideremos un remanente de unos 10 pc de radio que se expande en un medio que originalmente
tenı́a una densidad n = 0,1 cm−3 . El material desplazado por la onda de choque es N ∼ 1059
partı́culas. Si es un medio formado por H, M ∼ 100M⊙ . Si el remanente esta a una distancia d
! !
M d
Fπ0 →2γ (E > 100 MeV) ∼ 10 −10
κS cm−2 s−1 .
100M⊙ kpc
Para un remanente a 2 kpc se precisa κS > 200 a fin de que sea detectable por un instrumento
como EGRET. Por el contrario, si el remanente se encuentra en un medio denso (n > 1 cm−3 ), esta
en interacción con alguna nube molecular o es muy cercano, puede ser detectable.
En remanentes no muy brillantes en radio (emisión producida por la radiación sincrotrón de lep-
tones relativistas), se espera en general que la emisión de origen hadrónico domine sobre el Bremss-
trahlung relativista. En remanentes jóvenes, por otro lado, con alta densidad de rayos cósmicos que
se expanden en un medio no muy denso, la interacción Compton inversa con fotones del CMB puede
ser dominante.
9.2.3. Binarias de estrellas tempranas

Estos sistemas están formados por estrellas tempranas como Wolf-Rayet o estrellas O, que tienen
fuertes vientos. Estos colisionan formando una región delimitada por dos ondas de choque que pueden
acelerar partı́culas hasta velocidades relativistas.
Los electrones se enfrı́an a través de radiación
sincrotrón en radio (puede llegar hasta el IR o el
óptico). Como las estrellas emiten mucha radiación
electromagnética (con pico en UV), los electrones
se enfrı́an también por interacciones Compton in-
verso. Estos últimos fotones caen en el rango X-γ.
Los protones acelerados en la región de colisión de
vientos también pueden producir rayos γ por de-
caimiento de π 0 ; este proceso puede ser eficiente si
los vientos son lo suficientemente densos, o si los
protones logran difundirse hasta la base del viento.
La única binaria con colisión de vientos que se ha confirmado que es una fuente de rayos γ y rayos
X no térmicos es η-Carinae, aunque este sistema no tiene emisión sincrotrón observable puesto que
dicha radiación es autoabsorbida en la fuente.
9.2.4. Núcleos galácticos activos (AGN)

Fueron las primeras fuentes extragalácticas de rayos γ identificadas.
Los AGNs presentan emisión electromagnética a lo largo de todo el espectro. De hecho, los blazares
-un tipo de AGN con jets apuntando hacia el observador - son las fuentes de rayos γ más númerosas
del tercer catalogo Fermi. Las distancias determinadas a partir del corrimiento cosmológico al rojo
de sus lı́neas espectrales y los flujos observados permiten inferir luminosidades gigantescas (≥ 1044
erg s−1 ). Al mismo tiempo, suelen presentar variabilidad muy rápida, lo que implica que la región
de emisión es extremadamente compacta. Resulta natural, entonces, suponer que la acreción de
materia sobre un objeto compacto (por ejemplo, un agujero negro supermasivo con M ∼ 107−10 M⊙
es responsable de la emisión.
81
Para que la acreción esférica de un gas sobre un
objeto se mantenga, la atracción gravitatoria debe
superar a la presión de radiación
σT L GM mp 4πGcmp
|F⃗rad | ≤ |F⃗grav | =⇒ ≤ =⇒ L M
4πcr 2 r 2 σT
donde M es la masa del objeto compacto y L, su
luminosidad. Además, asumimos un gas de H.
La luminosidad a la cual la acreción se detiene
se denomina luminosidad de Eddington
4πGcmp
!
M
LEdd = M∼
= 1,3×1046 erg s−1 .
σT 108 M⊙
Asociada a esta luminosidad puede definirse la
tasa de acreción de Eddington
!
LEdd M
ṀEdd = 2 ∼ 0,23 M⊙ yr−1 .
c 10 M⊙
8
Algunos AGNs presentan luminosidades mayores que la de Eddington. Esto se interpreta como
un indicio de la anisotropı́a de la emisión. De hecho, observaciones con interferómetros de radio
muestran que la emisión no térmica forma jets que emanan de la fuente central.
La tasa a la cual la energı́a potencial de la materia en acreción es convertida en radiación es
GM m dU GM dm GM Ṁ
U= , L≈ = = ,
r dt r dt r
donde Ṁ es la tasa de acreción sobre el objeto compacto. Si consideramos que la ultima órbita estable
alrededor del agujero negro esta a 3RS , donde RS = 2GM/c2 es el radio de Schwarszchild, la energı́a
liberada por la materia acretante es
GM Ṁ 1
L= ∼ Ṁ c2 .
3RS 6
La eficiencia del proceso, por tanto, es mucho mayor que la de procesos termonucleares. De hecho,
la eficiencia radiativa real es un poco menor ya que por el teorema del Virial, la mitad de la energı́a
va a calentar el gas y la otra mitad es radiada. Luego
L ∼ 0,1Ṁ c2 .
La energı́a radiada se emite según la ley de Wien

!1/4
GM Ṁ GM Ṁ
= 2πr2 σSB T 4 =⇒ T (r) = .
2r 4πσSB r3
Si la materia tiene momento angular formará un disco alrededor del objeto compacto. Entonces,
la temperatura variará como
1/2 #)1/4
3GM Ṁ
( "
Rin

T (r) = 1− ,
8πσSB r 3 r
82
donde Rin es el borde interno del disco de acreción. A cada r el disco radia como un cuerpo negro
ν3
Bν (T ) ∝ .
exp(hν/kT ) − 1
El espectro total se obtiene integrando sobre la superficie del disco
ˆ Rout
I(ν) ∝ Bν (T (r))r dr.
Rin
Vemos un espectro tı́pico en la Figura 23.

A energı́as más altas la emisión puede tener una componente debida a Comptonizacion térmica
pasando luego a ser no térmica y dominada por jets.
El mecanismo de producción de los jets no esta cla-
ramente establecido. Es muy probable que los campos
magnéticos asociados al disco de acreción jueguen un
papel importante en el lanzamiento y la colimación del
plasma relativista. El fluido puede comenzar como un
plasma térmico que se mueve a velocidades relativistas
macroscópicas. Ondas de choque pueden entonces trans-
formar la energı́a cinética macroscópica en energı́a de las
partı́culas, que al volverse relativistas se enfrı́an por ra-
diación sincrotrón e interacciones Compton inverso. Los
fotones blanco para estas ultimas interacciones pueden
ser los propios fotones sincrotrón o fotones de fuentes
externas como el disco o radiación reprocesada en nubes
que orbitan cerca del objeto compacto. Las nubes más Figura 23: Espectro tı́pico del disco de
cercanas se encuentran en la llamada región de lı́nea an- acreción delgado, en función de la frecuen-
chas (broad line region, BLR). Se trata de nubes que se cia.
mueven a velocidades ∼ 5000 km s−1 en el campo gravi-
tacional del objeto compacto. Más lejos hay una región de lı́neas angostas (narrow line region, NLR),
más lentas.
Si el jet tiene contenido hadrónico, rayos γ se pueden producir por reacciones del tipo
p + γ → p + π 0 , π 0 → 2γ,
p + γ → n + π + , π + → µ+ + νµ ,
p + γ → n + π+ + π−,
donde los fotones blanco para la reacción original pueden venir de una corona formada por un plasma
muy caliente (T ∼ 109 ) alrededor del agujero negro. Este plasma, que probablemente es calentado
por procesos de reconexión magnética, se enfrı́a por Comptonizacion de los fotones más frı́os del
disco produciendo rayos X que pueden interaccionar con protones relativistas para producir piones.
El decaimiento de estos piones lleva a rayos γ y a cascadas electromagnéticas. También es posible
que un haz de protones interaccione con una nube molecular y se produzcan rayos γ a través de
reacciones
p + p → p + p + π 0 , π 0 → γ + γ.
El contenido de materia de los jets de los AGNs, sin embargo, no esta claramente estableci-
do. Podrı́a tratarse de un fluido formado por e± relativistas, electrones y positrones relativistas y
electrones y protones frı́os, o por una mezcla de estas posibilidades.
83
9.2.5. Microcuasares (Mqs)
Se trata de sistemas binarios forma-
dos por una estrella y un objeto compac-
to (agujero negro o estrella de neutrones).
Si la estrella es de la secuencia principal
de baja masa, hablamos de Mqs de baja
masa. Si la compañera es, en cambio, una
estrella temprana de gran masa, se habla
de Mqs de alta masa.
Las caracterı́sticas definitorias de es-
tos sistemás son que la materia de la es-
trella se acreta sobre el objeto compacto
formando un disco de acreción que emite
en rayos X y que sobre el eje de rotación
del sistema se forman y coliman jets re-
lativistas que pueden ser detectados en
radio debido a su emisión sincrotrón. En
su morfologı́a general, los Mqs son ver-
siones a escala reducida de los quásares
extragalácticos.
Los Mqs se observan en distintos estados espectrales en rayos X. Los dos estados más caracterı́sti-
cos que pueden observarse en la mayorı́a de las fuentes son
1. Estado high-soft: estado dominado ”térmicamente”. Se caracteriza porque la distribución

espectral de energı́a tiene un claro pico a energı́as de ∼ 1 keV. Esta emisión se interpreta como
originada en el disco de acreción y puede aproximarse por un cuerpo negro con temperaturas
tales que kT ∼ 1 keV. En este estado, además, se observa una contribución no térmica débil
en forma de ley de potencia con un espectro muy blando (ı́ndice espectral nph ∝ E −Γ con
Γ < 2). Esta componente se cree originada en una corona de plasma caliente alrededor del
objeto compacto. En este estado no se observan jets. La emisión esta dominada por el disco de
acreción, que se extiende hasta la última órbita estable.
2. Estado low-hard: estado dominado por la emisión no térmica. Se caracteriza por una dis-
tribución espectral de energı́as en la cual la componente térmica está disminuida y corrida
hacia energı́as más bajas. Hay ahora una fuerte contribución en forma de ley de potencia con
un espectro mucho más duro (Γ < 1,6). En este estado se observan jets que emiten radiación
sincrotrón en radio, que podrı́a extenderse hasta los rayos X duros. También es posible que a
estas energı́as haya una contribución Compton inverso (E ∼ 100 keV). El disco de acreción, en
esta configuración, parece estar truncado a una cierta distancia del objeto compacto y el papel
de al corona ahora es más importante.
Los Mqs permanecen la mayor parte del tiempo en el estado low-hard. El campo al estado high-
soft suele ser rápido y va acompañado de la eyección de componentes o plasmones relativistas cuyas
velocidades inferidas son mayores de lo que se espera sea la velocidad macroscópica del jet en el
estado low-hard. El movimiento de estas componentes, cuando es proyectado en el plano del cielo,
puede aparecer como supralumı́nico, debido al efecto de aberración relativista.
84
En el estado low-hard, la potencia total
del jet Lj esta relacionada con la tasa de
acreción sobre el objeto compacto
Lj = qj Ṁ c2 , qj ∼ 0,1 − 0,01.
Una parte de esta potencia se halla en

forma de partı́culas relativistas, por lo que
el jet puede generar rayos γ por medio de
interacciones Compton inverso. Los campos
de fotones blanco para los leptones pueden
ser el propio campo sincrotrón, el campo de
disco, la corona o la estrella compañera (no
presente en el caso de los AGN). Si el jet
tiene hadrónes relativistas, entonces inter-
acciones con protones del viento de la es-
trella compañera pueden dar lugar a rayos
γ en el caso de Mqs de alta masa.
Debido a que presentan orbitas excéntri-
cas y cambios de estado, los Mqs deberı́an
ser fuentes de rayos γ variables. La prece-
sión del jet también puede jugar un papel
importante en la generación de la variabili-
dad, cambiando el ángulo de la emisión con
la lı́nea de la visual.
9.3. Fuentes transitorias de

rayos γ: Erupciones de rayos
gamma
Las erupciones de rayos γ (Gamma-Ray Bursts, GRB) son un fenómeno astronómico caracterizado
por un rápido incremento de la radiación γ que llega a la Tierra desde el espacio. Este incremento
puede ser tal que llegue a superar a toda otra fuente γ del Universo. La duración de este fenómeno es
muy corta, usualmente yendo de unos pocos segundos a algunas decenas de ellos. Se han observado
GRB con duraciones extremas del ms y de decenas de minutos. En promedio se detecta entre 1 y 2
GRB por dı́a.
9.3.1. Caracteristicas fenomenologicas

La morfologı́a de los perfiles temporales de los GRB es extremadamente variada. Algunos tienen
una estructura sencilla caracterizada por un solo pico, mientras que otros presentan picos múltiples
y subestructura el nivel de milisegundo.
La distribución temporal de los GRBs parece ser bimodal, existiendo dos clases bien definidas
GRB cortos con duración T < 2 s y GRB largos con T ≥ 2 s. Se supone que esta dicotomı́a refleja
una diferencia intrı́nseca en el mecanismo que genera ambas clases de eventos.
85
9.3.2. Propiedades espectrales
La energı́a de los fotones que caracterizan a los GRB esta tı́picamente en el rango que va de
algunas decenas de keV a unos pocos MeV. En algunos casos excepcionales se han observado fotones
de hasta 10 GeV. El espectro es claramente no térmico y puede representarse por
(
AE −α e−(E/E0 ) , E < E0
n(E)dE = n0
BE β , E > E0
con α ∈ (∼, ∼ 1), β ∈ (∼ −2, −3) y E0 (∼ 0,1, ∼ 1) MeV.

La intensidad integrada en el tiempo T es del orden de
F ∼ (0,1 − 10) × 10−6 erg cm−2 .
Distribución espacial: el instrumento BATSE mostró claramente que la distribución de los

GRBs es altamente isotrópica, lo que sugiere un origen extragaláctico.
La distribución de intensidades, sin embargo, muestra una falta de eventos débiles respecto a lo
que se espera de una población homogéneamente distribuida en el espacio. Esto sugiere que hay un
lı́mite espacial máximo hasta el cual se producen los GRB y que estamos observando algunos de esos
eventos extremadamente lejanos.
Contrapartidas a bajas energı́as: satélites de sayos X como Beppo SaX han detectado con-
trapartidas a energı́as más bajas de ciertos GRBs. Estas también han sido halladas con telescopios
ópticos y de radio. Las duraciones de los eventos son mayores a energı́as más bajas, llegando a unas
cuantas semanas y meses en radio. Las contrapartidas son variables y su intensidad decae como leyes
de potencia: Fν ∝ t−α , (α = 1,1 − 2,1). Los espectros de las contrapartidas también son leyes de
potencia o en algunos casos leyes de potencia quebradas.
Galaxias anfitrionas: debido a las localizaciones precisas obtenidas primero por Beppo Sax y
Swift, ha sido posible determinar el corrimiento al rojo e incluso el tipo de galaxia en las que ocurren
los GRBs. Parece existir una tendencia a que los GRB de larga duración ocurran en galaxias ricas
en formación estelar.
9.3.3. Modelos para Gamma-Ray Bursts

Modelo fenomenológico básico: la ”bola de fuego”. La escala temporal de variabilidad de los GRBs
es en general muy corta δT ∼ 1 ms, por lo que a primer orden
R ≲ c δT ∼ 3 × 102 km,
o sea que la fuente inicial debe ser muy compacta. El objeto más compacto que existe es un agujero
negro, para el cual R = 2GM/c2 , lo que para un GRB implica
c3 δT
M≥ ∼ 102 M⊙ .
2G
Si el objeto no es un agujero negro, su masa debe ser aun menor. De aquı́ vemos que, sea lo que
sea, un GRB debe ser un fenómeno de naturaleza estelar.
La energı́a liberada en un GRB es enorme
! !2
F d
E0 = F (4 ) ≡ 10
2 51
erg.
10 erg cm−2
−6 3Gpc
86
Si esta cantidad enorme de energı́a es liberada en una región muy pequeña de volumen V ∼ 1023
cm3 , la fuente será opaca a su propia emisión γ, ya que debido a interacciones γ + γ → e− + e+
resulta τγγ ≫ 1. Especı́ficamente,
! !2 !−2
fp σT F D2 F d δT
τγγ = ≡ 1014 fp ,
2
R me c 2 10 erg cm−2
−7 3Gpc 10ms
donde fp es la fracción de fotones con energı́as por arriba del umbral de creación de fotones. Para
valores tı́picos τγγ ∼ 1014 , la fuente es opaca. Esto presenta dos grandes problemás:
1. la radiación deberı́a ser térmica, mientras que la que se observa es no térmica,
2. para tal opacidad no deberı́an observarse fotones de más de 0,511 MeV, mientras que claramente
se observan muchos fotones con Eph > 1 MeV.
La solución a este problema consiste en suponer expansión relativista del plasma anterior. La
enorme presión de la radiación dentro de la fuente hace que esta se expanda relativisticamente con
un factor de Lorentz γ. Los fotones observados, pues están corridos hacia el azul, por lo que el factor
fp debe sufrir una expansión relativista.
La opacidad corregida es
!2 !−2
1014
!
F d δT
τγγ ≡ 2+2α .
γ 10 erg cm−2
−7 3Gpc 10ms
Como α ∼ 2, para tener τγγ < 1 se debe cumplir que γ ≥ 102 . Cuando la fuente llega a expandirse
a esta velocidad se hace transparente permitiendo que se observe la emisión no térmica de su interior.
Esta expansión ultrarelativista, a su vez, limita el contenido de materia de la bola de fuego. En
efecto, su γ > 100 entonces M < 10−5 M⊙ . Este es el famoso problema de la contaminación bariónica
que deben evitar los modelos de GRB.
En el modelo de ”bola de fuego”, la energı́a interna del sistema se transforma en energı́a cinética
de expansión macroscópica. Cuando la cascara en expansión de la bola de fuego colisiona con el medio
interestelar produce una onda de choque. Esta transforma parte de la energı́a cinética de expansión
en energı́a cinética de las partı́culas microscópicas, que entonces es radiada por mecanismo sincrotrón
o Compton inverso.
Ondas de choque pueden ser también producidas por colisiones entre cascaras eyectadas con dife-
rentes velocidades. Estos shocks son llamados internos y dan lugar a la variada morfologı́a temporal
de los GRB. Los shocks externos producen las contrapartidas en rayos X y frecuencias más bajas.
9.3.4. Colimación
En algunos GRB para los que se ha medido el corrimiento al rojo, las energı́as internas inferidas
bajo la suposición de emisión isotrópica son tan grandes como 1054 erg. Esto es más de M⊙ c2 , por
lo que se supone que la emisión no debe ser isotrópica. Los requerimientos energéticos disminuyen
considerablemente si se asume que la emisión esta colimada en forma de jets. La colimación, además,
ayuda a explicar ciertas peculiaridades de la emisión temporal de las contrapartidas a bajas energı́as.
87
9.3.5. Mecanismo central
El mecanismo que produce la liberación de la energı́a en forma de radiación esta oculto por la
”bola de fuego”, que es opaca a su estadio inicial.
Se suele suponer que los GRB cortos son el resultado de fusiones de estrellas de neutrones con
agujeros negros en sistemás binarios de periodo ultra-corto, mientras que los GRB largos podrı́an ser
el resultado de la implosión de estrellas muy masivas.
Este ultimo modelo, llamado collapsar, consiste en la formación de un agujero negro en el centro
de una estrella masiva. El agujero es rodeado por un disco de acreción ultradenso, produciéndose
jets relativistas en el proceso. Parece haber abundante evidencia que asocia GRB con explosiones de
supernova.
9.3.6. Implicaciones cosmológicas

El Universo es transparente a los GRBs debido a que su emisión γ es de baja energı́a.
Los GRBs pueden dar información importante sobre la formación de las primeras estrellas y (a
través de las contrapartes de baja energı́a) del medio que las rodeaba.
En principio, debido a las restricciones existentes sobre la energı́a liberada, los GRB pueden dar
información sobre la energı́a oscura y los parámetros cosmológicos, en un dominio de z donde las
supernovas ya no son identificadas.
9.4. Fuentes no identificadas de rayos γ

La gran mayorı́a de las fuentes de rayos γ con Eγ > 100 MeV detectadas por EGRET no han
sido identificadas. Entre las identificadas hay 70 AGN, 6 púlsares, 1 fulguración solar, la Gran Nube
de Magallanes y una radio-galaxia.
Las fuentes no identificadas pueden dividirse en 3 grandes grupos:
1. Fuentes distribuidas isotrópicamente, de un claro origen extragaláctico.
2. Fuentes a latitudes galácticas medias que parecen formar un halo alrededor del centro galáctico.
Estas fuentes deben ser de origen galáctico, pero viejas (edades ∼ 108−9 anos).
3. Fuentes sobre el plano galáctico: están muy bien correlacionadas con trazadores de objetos de
Población I, por lo que deben ser jóvenes (edades 106−7 años) y se acumulan en las regiones
de formación estelar y en los brazos espirales de la Galaxia. Estas fuentes a su vez se pueden
dividir en variables y no variables. Las variables pueden ser objetos compactos como los Mqs,
mientras que las otras pueden incluir púlsares, nubes moleculares, remanentes de supernova,
entre otros.
10. Aspectos cosmológicos

El Universo esta lleno de radiación difusa conocida como radiación extragaláctiva difusa (RED).
La distribución de esta RED es isotrópica, pero en sus inhomogeneidades y en su espectro guarda
información sobre la composición y evolución del Universo.
Se han sugerido varias posibilidades sobre el origen de RED:
1. que se deba a la superposición de fuentes discretas no resueltas por los instrumentos,
88
2. que se origine por mecanismos difusos que operan a gran escala en el Universo o
3. una mezcla de los dos casos anteriores.
La multitud de objetos y procesos que puedan manifestarse en la RED solamente puede des-
entrañarse si se tienen determinaciones precisas de como se distribuye la radiación a lo largo del
espectro electromagnético.
La emisión difusa en microondas, infrarrojo, óptico y ultravioleta está casi con seguridad dominada
por procesos térmicos. La radiación de fondo de microondas es única en el sentido de que tiene un
origen cosmológico claro: es la radiación que escapo durante el periodo de recombinación del Universo
cuando la temperatura cayo, debido a la expansión, a niveles tales que los electrones pudieron ser
capturados por los núcleos - cambiando entonces súbitamente la opacidad y haciéndose el Universo
transparente a los fotones. La radiación del fondo cósmico de microondas (CMB) es el más perfecto
ejemplo de espectro de cuerpo negro conocido. Corresponde a un cuerpo negro con una temperatura
T = 2,725 K. No obstante, existen pequeñas anisotropı́as en el CMB del orden de 10−3 K. El estudio
de estas anisotropı́as en el CMB sirve para establecer las fluctuaciones originales de densidad a partir
de las cuales se formaron las galaxias.
En los rangos infrarrojo, óptico y ultravioleta, la emisión difusa esta dominada por radiación
térmica debida a galaxias, estrellas individuales y polvo. Las relaciones entre las intensidades de la
RED en las distintas bandas guarda relación con los contenidos y la evolución del Universo. Sin
embargo, estas componentes son influenciadas más por la evolución de las galaxias que por el modelo
cosmológico subyacente.
Las longitudes de onda a las cuales la RED tiene contribuciones no térmicas dominante son radio,
rayos X y rayos γ. Lo que parece dominar la emisión aquı́ es la contribución de diversos tipos de
núcleos galácticos activos.
La densidad de energı́a y de fotones de la RED están dominadas por la contribución del CMB. En
la actualidad, nuestro conocimiento de la radiación de fondo infrarroja y ultravioleta es demasiado
pobre como para arrojar valores confiables de densidad.
Longitud de onda Densidad de energı́a (eV m−3 ) Densidad de fotones (m−3 )

Radio ∼ 5 × 10−2 ∼ 106
Microondas 3 × 105 5 × 108
Infrarrojo - -
Optico ∼ 2 × 103 ∼ 103
Ultravioleta - -
Rayos X 75 3 × 10−3
Rayos γ 8.2 1,3 × 10−6
Recién con la llegada del satélite Compton y sus instrumentos COMPTEL y EGRET se pudo
medio con razonable exactitud el fondo difuso gamma. En el rango 100 MeV −820 GeV, de acuerdo
con las observaciones de Fermi-LAT, la forma del espectro puede ajustarse por una ley de potencias
con un corte exponencial
I(E) ∝ E −Γ exp − E/Ecorte ,
con Γ2,30 ± 0,02 y Ecorte 366 ± 71 GeV. A energı́as mayores de 100 GeV los rayos γ pueden ser
absorbidos por el fondo infrarrojo. La detección de fuentes discretas de rayos γ de energı́as de TeV
es siempre de objetos cercanos. De hecho, los cambios de su distribución espectral con el corrimiento
al rojo de las detecciones puede utilizarse para imponer restricciones al fondo infrarrojo.
Los objetos discretos que contribuyen a la RED en rayos γ son:
89
Galaxias normales: es difı́cil de establecer la contribución ya que no se conoce la función de
luminosidad de las galaxias normales. Su emisividad a estas energı́as depende de su contenido
de polvo y gas, ası́ como de la población de rayos cósmicos que exista en ellas.
AGNs: probablemente aportan la contribución dominante. Como los principales AGNs que
emiten en γ son fuentes intensas en radio, deberı́a haber una correlación entre los fondos
difusos en ambas bandas. Sin embargo, no es posible cuantificar la contribución de AGNs
radio-silenciosos y débiles. El espectro promedio en rayos γ de los AGNs detectados es de la
forma I(E) ∝ E −2,13±0,31 , lo que los hace candidatos firmás para explicar el fondo.
Galaxias ultra-luminosas: son galaxias con mucho gas y formación estelar, por lo que se
espera también que sean ricas en explosiones de supernova y, por ende, en rayos cósmicos.
Supernovas: las extragalácticas pueden contribuir a la RED en energı́as bajas debido a las
desexitaciones nucleares a las que dan lugar, como son 56 Ni →56 Co →56 Fe. Esta cadena produce
lı́neas de emisión entre 0,85 MeV y 3.25 MeV. Perdidas por dispersión Compton de estos fotones
pueden resultar en una distribución continua.
Agujeros negros primordiales: se podrı́an haber formado agujeros negros de baja masa
(< 1015 g) en el Big Bang, que se enfrı́an por radiación Hawking y se desintegran liberando
rayos γ. Sin embargo, el espectro esperado de la desintegración de una población primordial de
agujeros negros es una ley de potencia quebrada:
E −1 si E < 100 MeV
(
I(E) ∝
E −3 si E ≥ 100 MeV
que entra en conflicto con las determinaciones observacionales de la RED. Por lo tanto, queda
descartado que hay habido una formación masiva de agujeros negros primordiales de baja masa
en el Big Bang.
Además, se sugirieron fuentes intrı́nsecamente difusas:
Aniquilacion materia-antimateria: si el Universo estuviese lleno de antimateria uniforme-
mente distribuida, entonces las aniquilaciones
p + p → π+ + π− + π0,
producirı́an emisión γ difusa. Esta deberı́a exhibir un espectro que, promediando sobre todos
los redshifts, serı́a del tipo I(E) ∝ E −2,8 , lo que esta en conflicto con las observaciones. En
consecuencia, el Universo parece estar formado por materia - al menos en su mayor parte.
Partı́culas exoticas: observaciones astrofı́sicas de distintos tipos indican que debe existir
materia oscura en el Universo. Una porción considerable de esta materia parece ser no bariónica
y podrı́a estar formada por neutrinos u otras partı́culas más exóticas que no aparecen en el
modelo estándar. Partı́culas exóticas como WIMPS aparecen en extensiones supersimétricas
del modelo estándar. Estas partı́culas pueden aniquilarse produciendo rayos γ y contribuyendo
a la RED.
Halo galáctico extendido: es posible que exista una emisión difusa asociada a un halo
galáctico. partı́culas relativistas en este halo podrı́an contribuir al fondo difuso de rayos γ
a altas latitudes galácticas. Sin embargo, no es claro como las partı́culas se difundirı́an y
acelerarı́an en el halo, en el caso de existir.
90
Al presente, la opinión dominante es que la contribución dominante a la RED en rayos γ proviene
de AGNs no resueltos.
91

Resumen Iar

Cargado por

Copyright:

Formatos disponibles

Resumen Iar

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Resumen Iar

Cargado por

Copyright:

Formatos disponibles

Resumen IAR

Introducción a la Astrofı́sica Relativista

1. El espacio-tiempo y la relatividad especial

(s1 , s2 ) ∈ EL (x) = SL (x) × SL (x).

El espacio-tiempo es un sistema fı́sico que contiene a la colección de todos los eventos.

El espacio-tiempo se representa por una variedad real cuadri-dimensional diferenciable.

1.1. Objetos y estructuras sobre la variedad

El tensor T es n veces contravariante y m veces covariante.

ds2 = gµν dxµ dxν ,

donde gµν es el tensor métrico. En un espacio-tiempo euclı́deo, por ejemplo

conocido como tensor de Minkowski. Aquı́, el intervalo resulta

ds2 = ηµν dxµ dxν = (dx0 )2 − (dx1 )2 − (dx2 )2 − (dx3 )2 .

con v la velocidad del sistema, obtenemos las siguientes equivalencias

ds2 < 0 ←→ v > c,

Figura 1: Diagrama de un cono de luz. ηµν xµ T ν > 0,

con T ν = (1, 0, 0, 0). En forma similar, xµ señala hacia el pasado si

que es el tiempo que mide un reloj fijo al sistema fı́sico.

Como γ ≥ 1, el tiempo respecto al sistema propio se dilata. Si β = β(t), tenemos

1.2. El grupo de Lorentz

Uα′ V ′α = Lβα Uβ Lαµ V µ = Uµ V µ

Hay cuatro casos posibles de conjunto de transformaciones de Lorentz posibles:

1. det(L) = 1 y L00 ≥ 1 es el grupo propio de Lorentz,

2. det(L) = 1 y L00 ≤ 1 son las inversiones espacio temporales,

3. det(L) = −1 y L00 ≥ 1 son las inversiones espaciales y

4. det(L) = −1 y L00 ≤ 1 son las inversiones temporales.

P = {xµ → x′µ = Lµν xν + tµ }

La acción de una partı́cula es ˆ 2

Para velocidades bajas, hacemos un desarrollo en Taylor alrededor de v = 0 y vemos que

L ≈ −αc + 21 αc v 2 + ... expresión relativista,

Energı́a (usando el Hamiltoniano H)

Con la definición de energı́a encontramos

1.4. Elementos de la relatividad general

1.4.1. Agujeros negros

2. Partı́culas elementales <3

Generación leptón Q Le Lµ Lτ mc2

Nucleones p ≡ uud (0.938 GeV), S=0

Pion: π − ≡ du, π 0 ≡ uu, π + ≡ ud, con S = 0.

Eta: η ≡ dd, η ′ ≡ ss, con S = 0.

• Con extrañeza S = 1: K 0 ≡ ds, K + ≡ us.

Partı́cula Q mc2 [GeV]

El cuanto de interacción del campo electromagnético es el fotón. El fotón, a diferencia de las

2.4. El bosón de Higgs

2.5. Partı́culas más allá del Modelo Estándar

2.6. Interacciones entre partı́culas

El resultado de la interacción queda completamente determinado por las condiciones iniciales,

Decaimiento β del neutron:

2.7. Propiedades intrı́nsecas de las partı́culas

donde ⃗v es la velocidad de la partı́cula y e su carga eléctrica.

donde Z es el número atómico.

la intensidad máxima de los campos magnéticos,

el tamaño del anillo y

las pérdidas de energı́a de las partı́culas (por radiación sincrotrón).

Un elemento esencial de todo experimento con aceleradores de partı́culas, independiente de la

3.2. Rayos cósmicos

N (E)dE ∝ E −2,7 dE.

A energı́as mayores a ∼ 30 GeV,

3.2.1. El criterio de Hillas

El campo magnético superficial tı́pico B ≈ 1012 G. Sean Ω la velocidad angular de la estrella y

P V γ = cte, T V γ−1 = cte T P (1−γ)/γ = cte,

Esta distancia corresponde a un tiempo