Curso de Radiación - (Alberto Carramiñana)
Curso de Radiación - (Alberto Carramiñana)
Curso de Radiación - (Alberto Carramiñana)
Ondas electromagnéticas
1
Alternativamente, estas ecuaciones se pueden expresar en el sistema in-
ternacional (SI),
~ = ρ/ε0 ,
∇·E ~ = 0,
∇·B
~ ~
~ = − ∂B ,
∇×E ~ = µ0~ + µ0 ε0 ∂ E .
∇×B (1.5)
∂t ∂t
Estas unidades emplean diferentes unidades para los campos eléctrico y
magnético, con ESI = cEcgs , e introducen las constantes ε = 1/4πc, la
permitividad del vacı́o, y µ0 = 4π/c, la permeabilidad magnética del vacı́o,
√
de las cuales se calcula la velocidad de la luz, c = 1/ µ0 ε0 .
Las dos formulaciones anteriores presuponen campos en el vacı́o, o en
términos microscópicos. La formulación macroscópica para un medio de per-
mitividad ε y permeabilidad magnética µ es,
∇·D~ = ρ, ~ = 0,
∇·B
~ ~
~ = − ∂B ,
∇×E ~ = ~ + ∂ D ,
∇×H (1.6)
∂t ∂t
~ = εE,
con D ~ el campo de desplazamiento, y H ~ = B/µ,
~ el campo magnético,
~ siendo ahora la inducción magnética.
con B
En lo que sigue se emplearán las expresiones gaussianas, ecs. ?? a ??.
2
Figura 1.1: Las lı́neas de campo eléctrico divergen de una carga puntual
positiva. Para una carga negativa se invierten las flechas.
3
trostático queda como,
Z ρ(r~0 ) ~r − r~0
~ r) =
E(~ d3 r0 . (1.10)
|~r − r~0 |3
4
Ley de Ampere - caso magnetostático
En el caso estacionario, sin dependencia temporal, la ecuación ?? se escribe
como
∇×B ~ = 4π ~j,
c
~ ~
con j = j(~r). Esta ecuación proviene de la ley de Ampere, que se recupera
al integrar sobre una superficie normal a ~j,
~ = 4π I ~ · ~` = 4π I .
Z I
∇×B ⇒ B (1.12)
A c C c
Si consideramos una superficie acotada por un cı́rculo de radio R, normal n̂
y vector tangente ϕ̂ la ley de Ampere resulta en
~ = 2I ϕ̂,
B
Rc
al dividir por 2πR, la circunferencia de C. El campo magnético correspon-
diente a una corriente I que lleva un cable lineal infinito decae como 1/R,
siendo R la distancia a la recta.
Maxwell generalizó la ley de Ampere al introducir la corriente de despla-
zamiento en la ecuación ??. El término añadido por Maxwell corresponde a
la inducción de un campo magnético por un campo eléctrico variable en el
tiempo - fuera del caso magnetostático.
5
Figura 1.3: Un campo magnético con dependencia armónica en el tiempo
induce un voltaje con misma dependencia temporal al atravesar un circuito.
~ + ~v × B
F~ = q E ~ . (1.15)
c
La fuerza de Lorentz proporciona una aceleración paralela al campo eléctrico
y perpendicular al campo magnético. En el caso de campos con una direc-
~ y girar
ción fija la tendencia de una carga es moverse en la dirección de E
~
alrededor de las lı́neas de campo magnético B.
El teorema de Poynting proporciona una relación análoga a la ecuación
de continuidad, que se deriva de estimar el trabajo que realiza la fuerza de
Lorentz sobre un conjunto de cargas discretas. Definiendo la densidad de
6
energı́a del campo electromagnético,
1 2
u≡ E + B2 , (1.16)
8π
y el flujo de energı́a electromagnética, o vector de Poynting:
~≡ c E
S ~ ×B
~ , (1.17)
4π
obtenemos de las ecuaciones de Maxwell la relación:
∂u ~ = −~ · E
~.
+∇·S (1.18)
∂t
La diferencia principal entre la ecuación ?? y la ecuación ?? radica en el
término de la derecha: mientras que la carga eléctrica es una cantidad que
se conserva, la energı́a electromagnética puede no conservarse y ser disipada
o producida. El principio de conservación de la energı́a se mantiene ya que
la energı́a electromagnética disipada se convierte en otro tipo de energı́a,
~ corresponde con la ley de Ohm.
como mecánica o térmica. El término ~ · E
7
Estas ecuaciones pueden simplificarse aprovechando que los potenciales no
están definidos de manera unı́voca. Dada una función escalar arbitraria ψ,
podemos reemplazar el potencial vectorial A ~ por A ~0 = A~ + ∇ψ y el po-
0
tencial escalar ϕ por ϕ = ϕ − (1/c)(∂ψ/∂t) , sin dejar de cumplir con las
ecuaciones homogéneas de Maxwell. El seleccionar la función ψ de manera
a cumplir ciertas condiciones corresponde a la elección de una norma. En la
práctica hay dos normas que permiten simplificar las ecuaciones ?? y ??: la
norma de Coulomb y la norma de Lorentz.
La norma de Coulomb:
~ = 0,
∇·A (1.23)
es conveniente para tratar problemas de electrostática. La primera de las
ecuaciones de Maxwell queda como una ecuación de Poisson, ∇2 ϕ = −4πρ,
de solución general
ρ(~r0 , t) 3 0
Z
ϕ(~r, t) = d r . (1.24)
|~r − ~r0 |
En el caso estacionario (∂ϕ/∂t = 0, ∂ A/∂t~ = 0), el potencial vectorial
~ = −4π~/c.
también satisface una ecuación de Poisson, ∇2 A
La norma de Lorentz:
∇·A ~ + 1 ∂ϕ = 0 , (1.25)
c ∂t
permite escribir las dos ecuaciones inhomogéneas de Maxwell en forma de
ecuaciones de onda con fuentes:
1 ∂2ϕ 1 ∂2A~ 4π
∇2 ϕ − = −4πρ , ~−
∇2 A = − ~ . (1.26)
c2 ∂t2 2
c ∂t2 c
~ en términos de las
Las ecuaciones ?? pueden integrarse, quedando ϕ y A
densidades de carga y corriente:
ρ(~r0 , t − |~r − ~r0 |/c) 3 0
Z
ϕ(~r, t) = d r , (1.27)
|~r − ~r0 |
~(~r0 , t − |~r − ~r0 |/c) 3 0
Z
~
A(~r, t) = d r , (1.28)
|~r − ~r0 |
Las fuentes de los campos deben ser evaluadas en tret = t − |~r − ~r0 |/c,
manifestando el retrazo temporal causado por la propagación a velocidad
finita de la interacción electromagnética (cap 4).
(Shu 11; Ryb 2.1; Jackson 6)
8
1.2. Ondas electromagnéticas en el vacı́o
1.2.1. Ecuación de onda
En ausencia de cargas y corrientes las ecuaciones de Maxwell dan lugar a
una ecuación de onda para los campos eléctrico y magnético:
1 ∂2E~ 2~
~−
∇2 E = 0 , ∇2 B~ − 1 ∂ B = 0. (1.29)
2
c ∂t2 c ∂t2
2
9
Figura 1.4: El espectro electromagnético: arriba se muestran los rangos es-
pectrales, en términos de energı́a y longitud de onda del fotón. En medio se
ilustra la absorción de atmósfera mediante una lı́nea horizontal interrumpida
por las ventanas de transparencia, marcando encima de la lı́nea telescopios
espaciales y por debajo instrumentos terrestres. Abajo se ilustra el rango de
emisión térmica y no térmica.
10
Figura 1.5: Transparencia atmosférica en el sitio del Caltech Submillimeter
Observatory en Hawái, para una densidad de columna de vapor de agua de
0.5mm en la atmósfera.
11
Infrarrojo
La banda infrarroja se describe comúnmente en términos de longitudes de
onda, siendo raro el uso de frecuencias por encima de 1 THz. El infrarrojo
lejano cubre aproximadamente 1 mm > ∼λ> ∼ 20 µm; el infrarrojo medio entre
20 y 2 micras, descrito por las bandas atmosféricas LMNQ; y el infrarrojo
cercano entre 2 y 0.8 µm=8000Å, descrito por las bandas JHK. El infrarro-
jo lejano es accesible sólamente desde el espacio, empleando observatorios
espaciales como Spitzer o JWST, y con telescopios en aviones, como SOFIA;
el infrarrojo medio y cercano es parcialmente accesible desde sitios secos a
gran altura.
Visible
El visible está definido por la respuesta del ojo humano, cuya respuesta
abarca los colores del arcoı́ris, contenidos en el intervalo 8000 Å > ∼ λ >∼
4000 Å, que corresponde también a fotones entre 1.5 y 3 eV. Observaciones
fotométricas emplean tradicionalmente las bandas UBVRI, del sistema de
Johnson, que corresponden aproximadamente a los colores ultravioleta (U),
azul (B), visual (V) -entre el verde y el amarillo-, rojo (R) e infrarrojo
cercano (I). Este sistema se ha expandido al infrarrojo para definir la bandas
JHKLMNQ. La tabla 1 muestra los rangos de estas bandas y el flujo de un
objeto de magnitud cero, expresado en Janskys1 .
Visible
Filtro U B V R I
λ0 (Å) 3600 4400 5500 7000 9000
∆λ0 (Å) 68 98 89 220 240
F0 (Jy) 1880 4650 3950 2870 2240
IR cercano IR medio
Filtro J H K L2 M N Q
λ0 (µm) 1.25 1.65 2.20 3.40 5.0 10.2 21.0
∆λ0 (µm) 0.30 0.35 0.40 0.55 0.3 5 8
F0 (Jy) 1770 636 629 312 183 43 10
1
1 Jy = 10−23 erg cm−2 s−1 Hz−1
12
Ultravioleta
El ultravioleta cercano es observable desde la Tierra, en la banda U, hasta
alrededor de λ >∼ 3200Å, donde el ozono atmosférico y otras moléculas impo-
nen un corte al espectro accesible. Observaciones desde el espacio permiten
estudiar la banda de 3000Å a λ ' 911 Å, lı́mite de ionización del hidrógeno,
donde el medio interestelar mismo se vuelve opaco a la radiación por medio
del efecto fotoeléctrico. Esto convierte al ultravioleta lejano en una banda
mayormente inaccesible, en el que observaciones desde el espacio sólo son
posibles en direcciones privilegiadas con baja densidad de columna de gas
interestelar. Esto se extiende hasta alrededor de λ ∼ 120 Å, o E ∼ 100 eV,
frontera con la banda de rayos X.
Rayos X
Los rayos X se describen en términos de energı́a con hν > ∼ 100 eV, has-
ta alrededor de 511 keV, la masa en reposo del electrón (E = mc2 ). Son
observables exclusivamente desde el espacio, con el factor adicional de la ab-
sorción por efecto fotoeléctrico del medio interestelar en la banda de rayos
X blandos, 0.1–1 keV. Los telescopios espaciales de rayos X, como Chandra
X-ray Observatory y XMM-Newton, son instrumentos sofisticados de altı́si-
ma precisión basados en el principio de enfocar haces mediante reflección
por incidencia rasante, funcional entre 0.1 y 10 keV, . Los rayos X duros,
de decenas o cientos de keVs, se observan normalmente sin colimación, em-
pleando el efecto fotoeléctrico. La introducción de telescopios con máscara
codificada ha permitido mejorar considerablemente la resolución angular de
telescopios como BeppoSAX e INTEGRAL. Más recientemente, NuSTAR ha
extendido el principio de la óptica de incidencia rasante hasta energı́aspor
encima de 10 keV.
Rayos gamma
En principio la radiación γ comprende toda radiación electromagnética con
energı́a mayor a 0.511 MeV, hasta valores arbitrariamente altos. Histórica-
mente fue descubierta como uno de los tres productos de la radioactividad,
siendo los otros rayos α, constituidos por núcleos de helio, y rayos β, consti-
tuidos por positrones e+ . Los rayos γ producto de la radioactividad corres-
ponden a emisiones mono-energéticas debido a transiciones entre niveles de
energı́a en núcleos atómicos que tienen valores entre decenas de keV hasta al-
gunos MeV. Fotones de mayores energı́as son producidos por aniquilación de
13
partı́culas y antipartı́culas, e− e+ → γγ, o procesos que involucran partı́cu-
las relativistas aceleradas previamente. Rayos gamma por debajo de 0.1 a
1 MeV pueden ser estudiados mediante el efecto fotoeléctrico; para energı́as
entre 1-30 MeV se emplea el efecto Compton para su detección con teles-
copios como COMPTEL, instrumento a bordo del Compton Gamma-Ray
Observatory (CGRO) entre 1991 y 2000, que realizó el único mapeo del cielo
a la fecha en la banda 0.7-30 MeV. Rayos γ de alta energı́a, entre 30 MeV y
unos 100 GeV, son estudiados con telescopios de producción de pares, como
CGRO-EGRET en 1991-2000, y más recientemente con Fermi-LAT.
Es posible la detección de fotones con energı́a por encima de 100 GeV
desde la Tierra a través de la detección de cascadas electromagnéticas em-
pleando telescopios de tipo Cherenkov atmosférico, o con detectores de su-
perficie, originalmente diseñados para la detección de rayos cósmicos me-
diante el mismo principio. Los telescopios Cherenkov atmosféricos, como
VERITAS, HESS y MAGIC, registran la emisión debida a partı́culas carga-
das que se mueven más rápido que la velocidad de la luz en el aire, v > c/n
con n el ı́ndice de refracción. La detección directa de partı́culas secunda-
rias de estas cascadas se puede hacer con detectores Cherenkov de agua,
como el High Altitude Water Cherenkov (HAWC) γ-ray observatory insta-
lado en Sierra Negra. HAWC es capaz de detectar fotones con energı́as de
E> 14
∼ 100 TeV = 10 eV, El observatorio de rayos cósmicos Pierre Auger ha
puesto cotas para el flujo de fotones con energı́as E > 2 × 1018 eV.
~ med = R(E)
E ~ = x̂ x cos(ωt − φx ) + ŷ y cos(ωt − φy ) . (1.36)
14
Figura 1.6: Modos de polarización elı́ptica para ∆φ entre 0 y 2π, y conside-
rando x = y ⇒ Q = 0. A la izquierda, ∆φ = nπ/2, obtenemos polarización
lineal, V = 0, o circular, U = 0.
15
inferior izquierda y superior derecha de la caja. Este es el caso de
polarización lineal.
16
Polarización ∆φ x , y Q U V χ β
Circular izquierda π/2 x = y 0 0 +I 0/0 +π/4
Circular derecha 3π/2 x = y 0 0 −I 0/0 −π/4
Elı́ptica izquierda 0 < ∆φ < π x 6= y 6= 0 >0 6= 0 >0 6= 0
Elı́ptica derecha π < ∆φ < 2π x 6= y 6= 0 <0 6= 0 <0 6= 0
Lineal (x) 0 − 2π y = 0 +I 0 0 0 0
Lineal (y) 0 − 2π x = 0 −I 0 0 0 0
Lineal (xy) 0 x = y 0 +I 0 +π 0
Lineal (−xy) π x = y 0 −I 0 −π 0
17
1.3. Ondas electromagnéticas en un plasma
El medio interestelar no es estrictamente un vacı́o y es posible medir efectos
de la propagación de las ondas electromágneticas por la presencia de car-
gas y corrientes. Se observan varios efectos importantes que se describen al
considerar estas cargas.
Consideremos la densidad de carga y de corriente de un plasma formado
por electrones de carga −e y densidad numérica ne , e iones de carga +Ze y
densidad numérica ni :
siendo ~ve y ~vi la velocidad tı́pica de cada especie2 , donde ~ve 6= ~vi . Plantea-
mos ρ = ρe + ρi = 0, y proponemos un comportamiento oscilatorio para la
~
densidad de corriente3 , ~ = ~0 ei(k·~r−ωt) . Generalizando el mismo comporta-
miento para los campos electromagnéticos, podemos escribir las ecuaciones
de Maxwell correspondientes como,
~k · E
~ = ~k · B ω~ 4π ω~
~ = 0, i ~k × E
~ =i B, i ~k × B
~ = ~ − i E . (1.43)
c c c
El comportamiento oscilatorio de ~ es una respuesta del medio a los campos.
De las ecuaciones ?? se desprende (B ~ ⊥ ~k) y (B
~ ⊥ E).
~
18
movilidad de los electrones. Bajo las suposiciones anteriores, el movimiento
de los electrones puede restringirse al plano xy,
x̂ ± iŷ −iωt
~v = v0 √ e ,
2
de donde
eB −ieE~
~ ∓i
−iωm~v = −eE ~v ⇒ ~v = , (1.45)
c m(ω ± ω` )
siendo ω` ≡ eB/mc la frecuencia de Larmor. Podemos obtener la expresión
para la densidad de corriente,
ine e2 ~
~ = −ne e~ve + ni Zi e~vi ' E, (1.46)
m(ω ∓ ω` )
ω
i~k × E
~ = −i ~.
εE
c
donde introducimos la constante dieléctrica,
ωp2
ε=1− , (1.47)
ω(ω ∓ ω` )
1/2
y definimos la frecuencia de plasma, ωp ≡ 4πne e2 /me ' 5.6 × 104 Hz · ne .
p
19
El número de onda k es imaginario para ω < ωp , lo que implica un comporta-
miento exponencial del campo, E ~ ∝ e±|k|z , y la atenuación correspondiente
de la onda, con una longitud caracterı́stica
εω 2 − k 2 c2 = 0 ,
√
calculamos la velocidad de grupo de la onda vg ≡ dω/dk = c ε ≤ c, ya que
para el plasma ε < 1 (ec. ??). La onda se propaga más lentamente que en
el vacı́o, con una velocidad dependiente de la frecuencia. De hecho, si ne es
variable, las ondas pueden no viajar en lı́nea recta.
Una aplicación relevante al medio interestelar proviene del estudio de
pulsares, donde pulsos de distintas frecuencias ω viajan a distintas veloci-
dades, ya que ε = ε(ω), y llegan con distintos retrasos con respecto a una
propagación en el vacı́o. Supóngase un pulsar a una distancia d y ω ωp .
El tiempo requerido para que una onda de frecuencia ω llegue a la Tierra
está dado por:
( )Z
4πe2 d 2πe2 D
Z d d
dz d 1
t= ' + ne dz = + . (1.50)
0 vg c 2 mω 2 c 0 c mc ω 2
20
relación ?? nos da el valor de D, la densidad de electrones integrada a lo largo
de la visual. En la práctica, se grafica el desfasamiento temporal del pulso en
21
distintas frecuencias (figura ??), y se ajusta conforme a ??. El valor prome-
dio de la densidad de electrones, deducido a partir de mediciones en pulsares
con distancias conocidas, da hne i ≈ 0.03 cm−3 . Existe un modelo bastante
completo de la distribución de electrones libres en la Galaxia, obtenido a
partir de mediciones a pulsares (Taylor & Cordes 1993). El valor promedio
de la densidad de electrones es frecuentemente usado para tener una prime-
ra estimación de distancias a otros pulsares, dDM = D/(0.03 cm−3 ). Esta
estimación debe tomarse con cierta precaución ya que puede estar afectada
por sobre-densidades locales de electrones.
k± = k0 ∓ ∆k , (1.52)
ωp2 ω`
∆k = , (1.53)
2ωc
es la diferencia en números de onda con respecto a luz blanca asociada a dos
distintas polarizaciones (derecha o izquierda). Si uno considera una onda ini-
~ = x̂E0 e−iωt ,
cialmente polarizada linealmente a lo largo del eje x̂, es decir E
que recorre una distancia d a lo largo del eje ẑ, el eje de polarización res-
pectivo rota un ángulo ψ dado por,
2πe3
ψ= RM , (1.54)
m2e c2 ω 2
22
que no es posible medir DM , es posible estimar Bk conociendo la distancia
D E
d al objeto y suponiendo una densidad de electrones, Bk ≈ RM/(hne i d).
Es a través de éstas mediciones que los radioastrónomos han concluido que
el campo magnético promedio en la Galaxia es de 3 microgauss.
23
Capı́tulo 2
Transferencia radiativa
Versión 1 de noviembre de 2016
2.1. Definiciones
2.1.1. Intensidad y cantidades derivadas
La variable fundamental para la descripción del transporte de energı́a ra-
diativa es la intensidad especı́fica, Iν , que se refiere conceptualmente a un
haz de radiación. Se define a partir del flujo de energı́a (E) por unidad de
tiempo (t), área (A), frecuencia (ν) y ángulo sólido (Ω),
dE
= Iν (k̂, ~r, t) k̂ · n̂ dΩ dν , (2.1)
dA dt
con n̂ la normal al elemento de área dA que atraviesa el haz, el cual se propa-
ga en la dirección k̂. Iν es función de la posición, del tiempo, de la frecuencia
y del vector de propagación1 Iν contiene información tanto de la distribución
espacial y espectral del campo de radiación, como de su nivel de isotropı́a.
Las unidades (cgs) de la intensidad especı́fica son [erg cm−2 s−1 Hz−1 Sr−1 ].
1
Formalmente la combinación (k̂, ν) corresponde con el vector ~k = (2πν/c)k̂. El vector
de propagación, k̂ = ẑ cos θ +(x̂ cos φ+ ŷ sin φ) sin θ se relaciona directamente con el ángulo
sólido, dΩ = sin θdθdφ.
1
Figura 2.1: Todos los haces de radiación que salen de (1) con una apertura
dada por el ángulo sólido dΩ1 atraviesan la diferencial de área dA2 = r2 dΩ1 ,
siendo r la distancia entre ambos puntos. De misma manera, aquellos que
pasan por el punto (2) habiendo salido de (1) abarcan un ángulo sólido
dΩ2 = dA1 /r2 . Se desprende que dA1 dΩ1 = dA2 dΩ2 .
donde dΩ1 es el ángulo sólido que abarca el receptor (2) visto desde (1). Por
otro lado, la cantidad de energı́a recibida en (2) está dada por
2
la recibida dentro del cono de apertura dΩ2 coincide con la emitida dentro
del área dA1 . En virtud de que dA1 dΩ1 = dA2 dΩ2 (figura 2.1), tenemos
teniendo unidades de [erg cm−2 s−1 Hz−1 ]. La intensidad media tiene la mis-
ma forma funcional que la densidad de energı́a, uν , que se define también
integrando la intensidad sobre el ángulo sólido,
1 4π
Z
uν ≡ Iν dΩ = Jν ; (2.5)
c c
el factor c nos permite pasar de flujo a densidad. La presión de radiación
tiene unidades de [erg cm−3 Hz−1 ]; representa la cantidad de energı́a por
unidad de volumen y de rango espectral en forma de radiación.
Flujo de energı́a
El flujo de energı́a es una cantidad que se mide frecuentemente de manera
observacional. Corresponde al primer momento de la distribución angular
de la intensidad, Z
Fν ≡ Iν cos θ dΩ , (2.6)
3
El flujo de energı́a tiene la particularidad de ser cero en el caso isotrópico.
Considerando una superficie con normal el eje ẑ, podemos definir dos compo-
nentes del flujo: la de incidencia paralela a ẑ y la de incidencia anti-paralela.
En coordenadas esféricas,
Z 2π Z π/2 Z 2π Z π/2
Fν+ = Iν cos θ sin θdθdϕ , Fν− = Iν cos θ sin θdθdϕ ,
0 0 0 π
(2.7)
de manera que Fν = Fν+ − Fν− .
Es importante distinguir entre el flujo emitido, Fν1 , y el recibido, Fν2 ;
bajo la conservación de la intensidad se relacionan aproximadamente como
Iν ' Fν1 /∆Ω1 = Fν2 /∆Ω2 , siendo ∆Ω1 el ángulo sólido dentro del cual
se da la emisión y ∆Ω2 el ángulo sólido que abarca la fuente vista por el
observador. Generalmente se tiene ∆Ω2 ∆Ω1 , y el flujo emitido es mucho
mayor que el observado, Fν1 Fν2 .
La forma más común de medir el flujo es mediante el sistema de mag-
nitudes, definidas de manera logarı́tmica en el óptico y bandas aledañas a
partir de:
mν = −2.5 log(Fν /Fν0 ),
siendo Fν0 un flujo de referencia para la frecuencia ν. El sistema de magni-
tudes se define de manera que la estrella Vega tenga magnitud cero en las
distintas bandas fotométricas. En radio-astronomı́a es más común el uso del
Jansky, que tiene la bondad de ser una unidad lineal, definida como
Presión de radiación
La presión de radiación corresponde al segundo momento de la intensidad,
1
Z
Pν ≡ Iν cos2 θ dΩ , (2.8)
c
4
La intensidad especı́fica y sus momentos caracterizan las componentes
espectrales de un haz de radiación; pueden integrarse sobre frecuencias,
Z ∞ Z ∞ Z ∞
I= Iν dν , u= uν dν , F = Fν dν , ... (2.9)
0 0 0
para obtener las cantidades globales.
5
Figura 2.2: (a) Campo isotrópico de radiación; (b) campo de una fuente
isotrópica; (c) esfera radiando semi-isotropicamente.
6
que modificar esta relación para dar lugar a la ecuación de transferencia
radiativa. Los procesos a considerar son: emisión espontánea e inducida,
absorción y dispersión.
dE = ν dV dΩ dt dν ,
εν = 4π ν /ρ ,
αν = n σν ,
7
siendo n la densidad numérica de partı́culas. La sección eficaz tiene unidades
de área, [cm2 ]. Alternativamente, se usa el coeficiente de opacidad, κν , que
se relaciona con la densidad de masa,
αν = ρ κν ,
8
donde Sν ≡ ν /αν es la función fuente (emisión / absorción). Expresada de
esta forma, la ecuación de transferencia (2.20) tiene como solución formal:
Z τν
0
Iν (τν ) = Iν (0)e−τν + Sν (τν0 )e−(τν −τν ) dτν0 . (2.21)
0
9
Figura 2.3: Absorción (arriba) y emi-
sión (abajo) de radiación dentro de un
cono de apertura dΩ en la dirección k̂.
Una fracción aλ de la radiación inciden-
te de longitud de onda λ es absorbida.
eλ (T )
eλ (T ) = aλ (T )Kλ (T ) ⇒ = Kλ (T ) , (2.23)
aλ (T )
10
Figura 2.4: Cuerpo negro emitiendo de manera isotrópica (rojo) después de
termalizar radiación anisotrópica incidente (en azul).
cλ−5
Kλ = ,
ea/λT − 1
la cual ajustó correctamente los datos experimentales, y cumple con el prin-
cipio de máxima entropı́a. Esta expresión requiere que los intercambios de
energı́a de los resonadores de Planck con el campo de radiación se divida
11
en un número entero de elementos, o cuantos, de dimensiones ε = hν. En
1905, Rayleigh y Jeans mostraron por separado que el postulado cuántico
es incompatible con los principios clásicos, y que bajo el principio clásico de
equipartición de la energı́a la radiación en equilibrio termodinámico deberı́a
cumplir uλ (T ) = 8πkT /λ4 , la distribución de Rayleigh-Jeans.
2hν 3 /c2
Bν (T ) ≡ , (2.24)
ehν/kT − 1
Flujo de energı́a
Al ser la función de Planck isotrópica, tenemos un flujo de energı́a nulo,
Fν = 0. Sin embargo, al integrar sobre medio hemisferio, de acuerdo a la
12
relación 2.7, tenemos Fν+ = πBν (T ). Esto coincide con el cálculo para una
fuente puntual isotrópica. Si integramos sobre frecuencias obtenemos,
F (T ) = πB(T ) = σT 4 . (2.26)
Presión de radiación
Al tratarse de una emisión isotrópica, tenemos que la presión de radiación
está dada por, pν = uν /3 , de donde p = aT 4 /3 . Esta es formalmente la
ecuación de estado de un gas de radiación.
Bν dν = Bλ dλ .
13
Figura 2.5: La función de Planck, Bν (arriba) o Bλ (abajo). A la izquierda
en unidades lineales y a la derecha en logarı́timicas. Las lı́neas punteadas
denotan las aproximaciones de Rayleigh-Jeans y de Wien.
Lı́mite de Rayleigh-Jeans
Para hν/kT 1 la Planckiana tiende a la función de Rayleigh-Jeans,
2ν 2
Iν(RJ) (T ) = kT , (2.30)
c2
14
derivada por Rayleigh y Jeans en 1905 bajo el principio clásico de equi-
partición de la energı́a. Si la función de Rayleigh-Jeans fuera correcta, un
objeto a temperatura ambiente serı́a extremadamente brillante en el ultra-
violeta
R ∞ RJ
extremo, llevando a la divergencia de su integral sobre frecuencias,
0 Iν dν → ∞, conocida en su momento como la catástrofe ultravioleta.
Lı́mite de Wien
Para frecuencias grandes, definidas por hν/kT 1, la función de Planck
tiende a la función de Wien,
2hν 3 −hν/kT
Iν(w) (T ) = e , (2.31)
c2
historicamente la primera descripción de la emisión de cuerpo negro. La
función de Wien es integrable y cumple cualitativamente con la ley de des-
plazamiento, λmax T = constante; pero falla a frecuencias pequeñas donde
predice incorrectamente Iν ∝ ν 3 . Para frecuencias altas la intensidad de la
radiación de cuerpo negro decrece exponencialmente con ν.
Temperaturas astrofı́sicas
La función de Planck permite asignar un valor de temperatura a una me-
dición de flujo o intensidad, particularmente para fuentes resueltas angu-
larmente, y valorar las condiciones fı́sicas del emisor. Los tres principales
estimadores son:
La temperatura de brillo, que se define a partir de la relación
Iν = Bν (Tb ) .
15
La temperatura efectiva se define a partir del flujo integrado en el
emisor,
F = σ Te4 .
El flujo medido se relaciona con el emitido mediante el tamaño angu-
lar de la fuente, Fobs = Fem ∆Ωobs /∆Ωem . Para una fuente isotrópica
∆Ωem = π, mientras que ∆Ωobs = πδθ2 , siendo δθ el radio angular de
la fuente observada. En el caso del Sol, Fobs = σ Te4 δθ2 corresponde
con la constante solar, s ' 1367W/m2 .
La temperatura efectiva es una medida del flujo integrado y se usa
frecuentemente para describir tipos de estrellas. La luminosidad de
una estrella está dada por
L∗ = 4πR∗2 σTe4 .
Para el Sol Te ' 5770 K; una estrella roja tiene Te ' 3000 K, mientras
que una estrella azul puede alcanzar Te ' 30 000 K.
ν /αν = Bν (T ) , (2.32)
16
En el caso de radiación térmica en el régimen de Rayleigh-Jeans, la
ecuación de transferencia puede escribirse como
dIν dTb
= −Iν + Bν (T (τν )) ⇒ = −Tb + T . (2.33)
dτν dτν
Por un lado, se tiene que la emisión térmica es igual a la emisión de cuerpo
negro si dIν = 0; lo cual se da para medios ópticamente gruesos, τν → ∞.
Para un medio en estricto equilibrio termodinámico, T es independiente de
τν , y la solución de la ecuación de transferencia puede expresarse como una
relación entre la temperatura de brillo y la del medio,
17
la emisión espontánea, caracterizada por el coeficiente A21 que mide
la probabilidad de transición por unidad de tiempo (s−1 ).
A21 /B21
n1 B12 J¯ = n2 B21 J¯ + n2 A21 ⇒ J¯ = , (2.34)
n1 B12 /n2 B21 − 1
siendo ésta una relación fundamental entre las poblaciones de los niveles y
la intensidad de la radiación. Si además se tiene el equilibrio termodinámico,
las poblaciones de los niveles se relacionan a través de la ley de Boltzmann,
18
de donde
A21 /B21
J¯ = . (2.35)
(g1 B12 /g2 B21 ) e−hν/kT − 1
Identificando esta expresión con la ley de Kirchhoff, Jν = Bν , se obtienen
las relaciones entre los coeficientes de Einstein:
2hν 3
g1 B12 = g2 B21 , A21 = B21 . (2.36)
c2
Estas son propiedades inherentes a los sistemas cuánticos, independientes
de las poblaciones de los niveles de energı́a, que relacionan de manera fun-
damental los procesos de emisión y de absorción. Este tipo de relaciones
se denominan relaciones de balance detallado y pueden emplearse fuera de
equilibrio termodinámico. Los coeficientes de emisión y absorción de estos
sistemas vienen dados por
hν hν
ν = φ(ν)n2 A21 , αν = φ(ν) (n1 B12 − n2 B21 ) . (2.37)
4π 4π
El cociente entre los coeficientes de emisión y absorción es la función fuente,
que en el caso mas general viene dada por
2.5. Dispersión
El proceso de dispersión pura (scattering en Inglés) consiste en la absorción
y re-emisión isotrópica y coherente de radiación. El término coherente indica
que las distribuciones espectrales de la radiación incidente y dispersada son
19
iguales. La dispersión se describe mediante un coeficiente de dispersión3 ,
(d)
αν , que cumple:
ν(d) = αν(d) Jν , (2.39)
donde Jν es la intensidad media de la radiación. En el caso particular de
dispersión sin absorción y sin emisión espontánea tenemos Sν = Jν , y la
ecuación de transferencia queda como:
dIν
= −αν(d) (Iν − Jν ) . (2.40)
ds
Dado que Jν es una integral de Iν , la ecuación 2.40 es de tipo integro-
(d)
diferencial y no es fácilmente soluble. El término αν Iν denota la fracción
de radiación “interceptada” por el dispersor, la cual es re-emitida isotrópi-
(d) (d)
camente a través del término αν Jν . Nótese que si αν depende de ν, la
radiación de distribución espectral localmente, Iν (s) 6= Iν (0), aunque no lo
haga en promedio, Jν (s) = Jν (0). Un ejemplo es la dispersión de Rayleigh
en la atmósfera terrestre, función creciente con la frecuencia, que resulta en
un Sol enrojecido y un cielo azul.
En el caso de radiación térmica podemos considerar la acción conjunta de
dispersión, absorción y emisión espontánea, distinguiendo el coeficiente de
(a)
absorción como αν y suponiendo válida la ley de Kirchhoff para el término
(a) (a)
de absorción, ν = αν Bν . La ecuación de transferencia radiativa queda
como:
dIν
= − αν(a) + αν(d) (Iν − Sν ) , (2.41)
ds
quedando la función fuente redefinida como:
(a) (d)
αν Bν + αν Jν
Sν ≡ (a) (d)
. (2.42)
αν + αν
(a)
La emisión espontánea está contenida en el término αν Bν , mientras que
(a)
la absorción pura es αν Iν . Los términos de dispersión son aquellos propor-
cionales a αν (d). A la suma de los coeficientes de dispersión y absorción,
(a) (d)
αν + αν se le denomina coeficiente de extinción. Un ejemplo concreto de
su empleo es la corrección de datos fotométricos o espectroscópicos en el
visible por extinción debida a la atmósfera terrestre.
3 (d)
Rybicki emplea σν , ν , aunque esta ya se introdujo para denotar la sección eficaz.
20
2.6. Atmósferas plano paralelas
Una atmósfera plano-paralela está estructurada a lo largo de una sola direc-
ción, que denominaremos ẑ (figura 2.7). Esta aproximación es útil para el
estudio de atmósferas delgadas en estrellas o planetas. Los tres casos que se
presentan son: la deducción de la ecuación de difusión radiativa; la aproxi-
mación de Eddington, incluyendo el caso de una atmósfera gris; y una solu-
ción para la transferencia de energı́a radiativa en atmósferas plano-paralelas
presentada por Simonneau.
21
El inverso de µ se denomina masa de aire. Podemos despejar Iν planteando
una relación propicia para un proceso de iteración,
µ ∂Iν
Iν = Sν − a . (2.44)
αν + ανs ∂z
A orden cero despreciamos el cambio de la intensidad con la profundidad,
(0) (0) (0)
por lo que tenemos Iν (z, µ) = Sν (z). Dado que Sν no depende de µ,
(0) (0)
es isotrópica ⇒ Jν = Sν . Esto es particularmente válido para un medio
(0)
con emisión térmica, Sν = Bν . La aproximación de Rosseland consiste en
suponer que la función fuente está dada por una superposición de funciones
de Planck con T = T (z). Substituyendo en la ecuación 2.44 obtenemos, a
primer orden,
µ dBν
(1)
Iν = Bν − a . (2.45)
αν + ανs dz
Podemos entonces calcular el flujo de radiación
Z +1 Z +1
2π dBν
Fν (z) = 2π Iν (z, µ) µdµ = − a s
µ2 dµ
−1 α ν + αν dz −1
4π dBν dT
= − a s
. (2.46)
3(αν + αν ) dT dz
Si integramos sobre frecuencias obtenemos la ecuación de difusión radiativa,
16σT 3 dT
F (z) = − , (2.47)
3αR dz
donde, empleando las propiedades de la función de Planck, definimos la
opacidad de Rosseland como
Z ∞ Z ∞
1 1 dBν dBν
≡ dν dν . (2.48)
αR 0 ανa + ανs dT 0 dT
La ecuación de difusión radiativa describe el flujo de energı́a originado por un
gradiente de temperatura. Esta relación de particular utilidad en el estudio
de interiores estelares.
22
Figura 2.8: Trayectorias de rayos de luz provenientes del centro del Sol o del
limbo.
1 +1
Z +1
c 4π
Z
Jν = Iν dµ = uν = aν , Fν = 2π Iν µdµ = bν ,
2 −1 4π −1 3
Z +1
2π 4π aν
pν = Iν µ2 dµ = .
c −1 3 c
En particular se tiene que campos semi-isotrópicos cumplen la misma rela-
ción entre densidad y presión de radiación que en el caso isotrópico,
1
pν = uν . (2.50)
3
Al substituir la aproximación (2.49) en la ecuación de transferencia radiativa,
µ(∂Iν /∂τν ) = Iν − Sν , se relacionan las funciones aν y bν con Sν .
Entre las aplicaciones más conocidas de la aproximación de Eddington
está la atmósfera gris, en la que la opacidad es independiente de la frecuencia,
τν = τ . Trabajando con las funciones a y b integradas sobre frecuencia,
y suponiendo un medio térmico, S = B = (σ/π)T 4 , podemos relacionar
T = T (τ ) con la temperatura efectiva imponiendo la condición de que en la
23
frontera externa de la atmósfera, τ = 0, el flujo es solamente emergente e
igual a F = σTe4 . Se obtiene
3σ 4 2 3 2
I(τ, µ) = Te µ + τ + , 4
T (τ ) = Te4 τ + . (2.51)
4π 3 4 3
Vemos que T (0) ' 0.841Te y que T (τ ) = Te en τ = 2/3. De la forma
completa para la intensidad (2.51) obtenemos la ley de oscurecimiento al
limbo,
I(0, µ) 3 2
= µ+ . (2.52)
I(0, 1) 5 3
Bajo las aproximaciones de Eddington y de atmósfera gris, la intensidad en
el limbo solar debe ser 0.4 veces la observada en el centro del disco del Sol.
24
Procedemos a integrar sobre λ, tomando en cuenta las definiciones de
opacidad media de Rosseland (ec. 2.48) y de Planck,
Z ∞
1
κp = κλ Bλ dλ , (2.57)
B(T ) 0
4π κp Te4
F = πB(0) + B(0) ⇒ T (0)4 = . (2.58)
3 κR 1 + 4κp /3κR
25
Capı́tulo 3
Termodinámica y mecánica
estadı́stica
3.1. Termodinámica
3.1.1. Las leyes de la termodinámica
La termodinámica estudia el comportamiento de sistemas macroscópicos
descritos por variables como energı́a, presión, volumen, temperatura o in-
cluso tensión, magnetización, etc.... La termodinámica no considera el estado
microscópico del sistema, dado por los estados de las partı́culas individuales
que lo componen.
El formalismo de la termodinámica se basa en el concepto de tempera-
tura, definido a partir del concepto de equilibrio termodinámico:
un sistema está en equilibrio termodinámico si las variables que
lo describen no varı́an con el tiempo.
De manera equivalente, dos sistemas están en equilibrio termodinámico
mútuo si al entrar en contacto las variables termodinámicas respectivas no
varı́an. Por definición la termodinámica se aplica a sistemas en equilibrio y
se basa en las leyes brevemente descritas a continuación.
1
Ley cero - definición de temperatura
Si de tres sistemas A, B y C, los sistemas A y B se encuentran cada uno en
equilibrio termodinámico con C, entonces A y B se encuentran en equilibrio
uno con el otro.
La ley cero define el concepto de temperatura a través de los conjuntos
de sistemas en equilibrio mutuo (clases de equivalencia): si A y B están en
equilibrio, entonces TA = TB .
dU = d−Q − P dV + µ dN +
X
Xı dxı , (3.1)
ı
d−Q = T dS
dS ≥ 0 . (3.2)
2
Como consecuencia un estado de equilibrio es un estado de máxima entropı́a.
Procesos que implican dS = 0 en un sistema cerrado son reversibles, mientras
que aquellos que implican dS > 0 son irreversibles.
Los principios de máxima entropı́a y mı́nima energı́a establecen que el
equilibrio termodinámico entre dos sistemas requiere no sólo TA = TB (equi-
librio térmico), sino también PA = PB (equilibrio mecánico). El requerimien-
to extra de equilibrio quı́mico se traduce en µA = µB . El equilibrio requiere
una relación entre las variables del sistema denominada ecuación de estado.
3
respectı́vamente. Para el gas ideal se tiene, empleando d−Q = dU + P dV ,
!
d−Q dU dU dV
CV = = , CP = +P = CV + N k , (3.5)
dT V
dT dT dT P
4
termodinámicas del sistema macroscópico basándose en una descripción es-
tadı́stica. Hay cuatro requisitos básicos para dar una descripción estadı́stica
de un sistema fı́sico: (1) especificar el estado del sistema; (2) especificar el
espacio de estados accesibles, o ensemble estadı́stico; (3) especificar la pro-
babilidad de cada estado; (4) hacer el cálculo estadı́stico.
# 1 - Estado de un sistema
La descripción del estado de un sistema puede darse conforme las leyes de
la mecánica clásica o la mecánica cuántica, dependiendo de qué formalis-
mo sea más exacto o incluso más práctico. En mecánica clásica el estado
de una partı́cula en un instante t puede describirse con su posición y mo-
mento, {~x(t), p~(t)}, suponiendo que tenemos el conocimiento necesario para
calcular las fuerzas que actúan sobre la partı́cula. En el caso de un siste-
ma de N partı́culas, requerirı́amos el conjunto de posiciones y momentos,
{~xj , p~j ; j = 1, ..., N }, que basta para especificar el (micro)-estado del siste-
ma; requerimos también conocer la interacción entre las partı́culas, que da
lugar a la fuerza que siente la i-ésima partı́cula, F~i (~x1 , ...~xN ) y de ahı́ cal-
cular cómo evolucionan posiciones y momentos. Esto es posible en principio
pero no siempre en la práctica.
En mecánica cuántica el estado de un sistema queda definido por la
función de onda, ψ(~x1 , . . . , ~xN ), que depende de las 3N coordenadas de las
partı́culas. Frecuentemente la función de onda puede escribirse en términos
de otras variables, concretamente los números cuánticos. Por ejemplo, un sis-
tema formado por dos osciladores armónicos tridimensionales está plenamen-
te definido si especificamos los seis números cuánticos {n1x , n1y , n1z , n2x , n2y , n2z }.
Las funciones de onda ψ{n1x ,n1y ,n1z ,n2x ,n2y ,n2z } son enumerables y podemos
hacer la estadı́stica del sistema. Nótese que hay varios estados con misma
energı́a, E = (n1x + n1y + n1z + n2x + n2y + n2z + 3)h̄ω .
En general, el localizar los estados dentro del espacio-fase de posiciones
y momentos es particularmente ventajoso. Esta descripción no especifica el
estado de cada partı́cula, sino cuantas partı́culas tienen mismo estado defi-
nido por una posición ~x y un momento p~. Cada par {~x, p~} define un punto
en el espacio-fase y la descripción del sistema está completa si se especifica
la ocupación del espacio-fase. El principio de exclusión de Pauli prohibe que
dos partı́culas tengan el mismo estado, dentro de la indeterminación inhe-
rente a la mecánica cuántica, que fija el tamaño de las celdas del espacio
fase a ∆x∆p = h/2π, (Figura 3.1). El lı́mite clásico corresponde a h → 0,
el paso de celdas discretas a continuas. Una descripción cuántica completa
considera el espı́n de las partı́culas, pudiendo ser entero -bosones- o semien-
5
Figura 3.1: La división del espacio
fase en celdas de tamaño ∆x∆p = h̄.
# 2 - Ensemble estadı́stico
Podemos especificar un estado arbitrario de nuestro sistema, pero no todos
los estados son consistentes con una medición macroscópica. En particular,
al fijar la energı́a total del sistema a un valor E restringe fuertemente el
conjunto de microestados posibles a aquellos que cumplen esta restricción,
o estados accesibles. Supongamos un sistema de cuatro partı́culas, pudiendo
cada una de ellas tomar como energı́a un valor entero entre 1 y 10. El número
total de estados que uno puede construir es 104 = 10000, de los cuales 35
son accesibles si la energı́a total es E = 8.
La restricción que se impone define el tipo de ensemble estadı́stico en
cuestión. Los dos ensembles que veremos son el canónico, donde la restricción
al sistema está en asignar a la energı́a a un valor E, y el gran canónico, que
restringe energı́a y número de partı́culas, {E, N }.
# 3 - Probabilidad de un estado
La probabilidad de un estado dado está dada por el siguiente postulado
fundamental:
en un sistema aislado en equilibrio todos los estados accesibles
tienen misma probabilidad.
6
En el ejemplo de la sección anterior, si el sistema tiene una energı́a total
E = 8 la probabilidad de que las cuatro partı́culas tengan la misma energı́a
es P (2222) = 1/35.
Este postulado trae implı́cito el principio de máxima entropı́a. El sistema
puede adquirir de manera estadı́stica todos los estados accesibles, maximi-
zando el desorden microscópico del mismo. La entropı́a se define a partir de
Ω, el número de estados accesibles,
# 4 - Cálculos estadı́sticos
Sea un sistema con energı́a total E y Ω(E) estados accesibles. De estos, sea
Ω(E, y ) el número de aquellos que cumplen una restricción adicional, y = y .
La probabilidad de encontrar al sistema con y = y está dada entonces por
P (y ) = Ω(E; y )/Ω(E) y el valor promedio de la variable y es:
1 X
ȳ = y Ω(E, y ) . (3.9)
Ω(E)
7
3.2.2. Relación entre variables macro y microscópicas
Existen varias relaciones entre variables termodinámicas que se derivan de
las leyes de la termodinámica. En particular la primera ley implica2 :
∂E ∂E ∂E
T = , P =− , µ= . (3.10)
∂S ∂V ∂N
El número de estados accesibles es función de E y las variables extensivas
x, como V y N ; es decir Ω = Ω(E, x). A partir de 3.8 y la primera ley de
la termodinámica, escrita como d ln Ω = βdE − β Xdx, con β = 1/kT ,
P
∂ ln Ω0
ln Ω0 (E 0 − Er ) ' ln Ω0 (E (0) ) − Er . (3.13)
∂E 0 0
∂ ln Ω0 ∂ ln Ω0
0
= ≡ β, (3.14)
∂E ∂E 0
donde por definición β el inverso de la temperatura de A0 , igual a la de
A en equilibrio. Se tiene entonces que Pr ∝ e−βEr y podemos escribir la
probabilidad normalizada del estado r como
e−βEr
Pr = P −βEr . (3.15)
e
r
8
Figura 3.2: Intercambio de
energı́a entre un sistema
canónico y su entorno, hasta
alcanzar el equilibrio termo-
dinámico.
de donde
∂ ln Z 1 ∂ ln Z 1 ∂ ln Z 1 ∂ ln Z
E=− ,X=− → P =+ , µ=− .
∂β β ∂x β ∂V β ∂N
(3.17)
La expresión 3.17 viene dando la ecuación de estado del sistema. La relación
entre la función de partición y la entropı́a se obtiene identificando 3
∂ ln Z ∂ ln Z
d ln Z = dβ + dx = −Edβ − βXdx = −Edβ − β(dE − T dS) ,
∂β ∂x
9
3.3.2. Gas ideal clásico
El gas ideal consiste en un conjunto de N partı́culas libres y la energı́a del
sistema está dada por la suma de las energı́as cinéticas individuales:
N
X p |2
|~
E= . (3.18)
=1
2m
3 2πmk
S = N k ln V + ln T + σ con σ = 3/2 ln e . (3.21)
2 h20
10
donde σs = σ+1. Esta expresión satisface SB = 2SA , resolviendo la paradoja
de Gibbs.
El problema irresoluble es el no cumplir la tercera ley, ya que S → −∞
para T → 0. Como veremos más abajo, el lı́mite de bajas temperaturas es
un caso cuántico.
El lı́mite clásico
Un gas ideal puede describirse empleando la mecánica clásica siempre y
cuando la ocupación de las celdas de espacio-fase sea baja. Esta condición
puede plantearse como
[∆x]3 [∆p]3 N h̄3 ,
11
donde [∆x]3 y [∆p]3 indican el rango de posiciones y momentos ocupados por
las N partı́culas. Claramente [∆x]3 = V , mientras que [∆p]3 ' (3mkT )3/2 .
La condición para la aplicabilidad de la estadı́stica clásica viene siendo,
h̄2
n−2/3 kT , (3.26)
3m
con n = N/V , la densidad numérica. El régimen cuántico se aplica para
temperaturas bajas o densidades altas.
12
donde w = exp {−β(ε − µ)}. La primera igualdad indica que la función
de partición es la sumatoria sobre todos los estados r del producto sobre
los numeros de ocupación correspondientes a cada estado, definido por un
(r) (r) (r)
conjunto de números de ocupación {n1 , n2 , n3 , . . .}. El término de la
derecha indica un producto sobre todas las celdas de espacio fase de la
sumatoria sobre todos los números de ocupación posibles. Las energı́as ε
dependen de la ubicación de la -ésima celda de espacio-fase pero no de
(r)
su estado de ocupación (es decir de n , por tanto de r), lo cual permite
permutar el producto y la suma sobre estados. Consideremos un ejemplo
simple, el de un sistema con tres microestados posibles representados por
tres celdas de espacio fase, = {a, b, c} que pueden ser ocupadas por cero o
una partı́cula (n = 0, 1), resultando en ocho estados r con distintos valores
de Er y Nr . En este sistema los dos términos de la igualdad (3.31) se escriben
como:
111
X 1
Y
wn = 1 + wa + wb + wc + wa wb + wa wc + wb wc + wa wb wc
r=000 =0
Y 1
X
wn = (1 + wa )(1 + wb )(1 + wc ) , (3.32)
=a,b,c n=0
13
3.4.2. Gas de fotones
En un gas de bosones no hay restricción en cuanto al número de partı́cu-
las que pueden ocupar un estado dado. La función de partición se obtiene
directamente de la sumatoria del lado derecho de la expresión 3.31
∞
1 Y 1
X
wn = ⇒ Z= . (3.33)
n =0
1 − w 1 − e−(ε −µ)/kT
14
Figura 3.3: Arriba: ocupación
media de estados n̄ en un
gas de fermiones como fun-
ción de la energı́a del estado,
para µ/kT = 1.0, 5, 10, 20, 50.
La función tiende a un es-
calón centrado en ε = µ, de
ancho ∆E ∼ kT , a medida
que µ/kT aumenta, como se
ilustra en la figura de abajo.
Gas degenerado
Consideremos primero un gas degenerado de electrones (m = me ), corres-
pondiente a materia contenida en una estrella enanas blancas. Tomando en
cuenta que los núcleos son quienes le dan la masa a la estrella, escribimos
15
la densidad de masa como ρ = Γne mH , siendo ne la densidad de electrones,
mH la unidad de masa atómica y Γ el cociente de nucleones a electrones,
el cual tiende a 2 para material procesado (para Fe Γ = 56/26). El número
total de partı́culas es igual al número de estados contenidos en una esfera
de Fermi, multiplicado por el factor de dos debido al espı́n,
Z pf 3
d p 2V 4π 3
Z
3
N = 2N = 2 d x = 3 p , (3.39)
0 h3 h 3 f
2 pf q 2 2
Z
u = p c + m 2 c4 − m c2 4πp2 dp,
e e
h3 0
" #
2 pf 1 p2 c2
Z
P = 4πp2 dp . (3.43)
h3 0 3
p
p2 c2 + m2e c4
16
relativista y el segundo (B) al caso relativista. Si las partı́culas degeneradas
son electrones, la transición se da para pf = me c2 que implica
8π (ΓmH ) m3e c3
ρ = ρc ≡ ' 0.97 × 106 Γ g cm−3 . (3.44)
3h3
Esta es la densidad tı́pica de una estrella enana blanca. Los dos casos son:
(A) no relativista: ⇒ ρ ρc :
Equivale a despreciar el término pc frente a me c2 en las integrales de
las ecs. 3.43, dando u = 3P/2, igual que para un gas ideal clásico, y,
2/3 5/3
3 h2 ρ
P = . (3.45)
π 20me ΓmH
P = Kρ5/3 es la ecuación de estado de un gas de fermiones no relati-
vista.
(B) relativista: ⇒ ρ ρc :
En cuyo caso el término de momento, pc, domina en las integrales
de las ecuaciones 3.43, resultando en u = 3P , como para un gas de
fotones, y la ecuación de estado,
1/3 4/3
3 hc ρ
P = , (3.46)
π 8 ΓmH
es decir P ∝ ρ4/3 para un gas degenerado de partı́culas relativistas.
Estrellas degeneradas
La estructura de una estrella degenerada puede investigarse de manera cua-
litativa y semi-cuantitativa con un modelo simplificado. Se puede calcular
la energı́a total de la estrella como la suma de la energı́a interna (térmica)
y gravitacional,
Z R Z R
2 GM (r)
Etot = u(r) 4πr dr − ρ(r) 4πr2 dr , (3.47)
0 0 r
siendo u = Kργ /(γ − 1) la densidad de energı́a, R el radio de la estrella
y M (r) la masa contenida dentro de un radio r. Un ejemplo relativamente
simple es el de una estrella con densidad constante, para la cual encontramos
la expresión
γ−1
KM γ 3 3 GM 2
E(R) = R−3(γ−1) − , (3.48)
(γ − 1) 4π 5 R
17
siendo la configuración de equilibrio determinada por la condición de mı́nima
energı́a, dE/dR = 0. Distinguimos los dos casos:
3 (5hc/8G)3/2
M = Mch = ' 3.42 × 1033 g ' 1.72 M . (3.51)
2π (Γmh )2
m4H c3
ρc = ' 6.1 × 1015 g cm−3 , (3.52)
3π 2 h̄3
18
obteniéndose para el caso no relativista
4/3 2 −1/3
3 h M
R= 8/3
' 4.9 km (M/M )−1/3 , (3.53)
2π 4GmN
3 (5hc/8G)3/2
Mns = ' 6.9 M . (3.54)
2πm2h
19
partición de m especies de partı́culas puede escribirse de la forma,
m m N
Y Y ζ
Z= Z = (3.55)
=1 =1
N !
condición como: 0
m m
ν 0
[ζ /n V ]ν = ζ`0 /n0` V
Y Y
` , (3.57)
=1 `=1
20
número de estados de energı́a E (estados degenerados de misma energı́a),
por lo que obtenemos la ley de Boltzmann:
n2 g2
= e(E1 −E2 )/kT . (3.59)
n1 g1
Para un sistema con un número arbitrario de niveles tenemos:
n g e−E /kT
= P −E /kT . (3.60)
n g` e `
`
H(0) *
) H(1) + e , (3.61)
ζ0 /V ζ1 /V ζe /V
= . (3.62)
n0 n1 ne
Tanto el electrón como el protón (H(1)) son partı́culas libres, por lo que la
ley de acción de masas, al considerar las distintas funciones de partición, da:
21
de ionización del átomo considerado6 . Los factores de 2 vienen de los dos
estados de espı́n posibles para el electrón y el protón. En estos términos
obtenemos:
n1 2 2
= e−χ/kT . (3.64)
n0 ne λ3e g0
Electrones no degenerados tienen una ocupación ∆x∆p h, por lo que
ne λ3e 1. De hecho la densidad es generalmente varios órdenes de magnitud
menor a λ−3 e , por lo que es común expresar este término mediante una
exponencial, ne λ3e /2 = e−ξ , con ξ tı́picamente entre 10 y 30. La ec. 3.64
puede generalizarse a otros estados de ionización ı 6= 0 y reescribirse como:
nı+1 gı+1 ξ−χı /kT
= e , (3.65)
nı gı
donde gı+1 es el número de estados para la especie ionizada debidos al espı́n
y gı es el número de estados de la especie no-ionizada.
Para que haya un grado de ionización apreciable se requiere que el
término en la exponencial no sea significativamente mayor que uno, lo cual
sucede para kT < ∼ χı /ξ, correspondiente a temperaturas mas de un orden de
magnitud por debajo de χı /k. La dificultad en ionizar un átomo es compen-
sada con creces por la vastedad del espacio fase disponible para el electrón
libre. Como resultado, mientras que el potencial de ionización del hidrógeno
corresponde a 13.6 eV/k ≈ 158, 000 K, es en las estrellas de tipo A, donde
Te ≈ 10, 000 K y ξ ≈ 15, donde la secuencia de Balmer es mas pronunciada.
La ecuación de Saha es básica para obtener la relación entre la intensidad de
las lı́neas espectrales y la temperatura, correspondiendo a la relación entre
tipo espectral y temperatura (fig 3.4).
6
la energı́a del estado ligado es E = −χ y el término de partición es e−E/kT = eχ/kT
22
Figura 3.4: Intensidad de las principales lı́neas espectrales en función de
la temperatura. La secuencia de tipos espectrales se define a partir de las
intensidades de las distintas lı́neas y por tanto la relación entre intensidad
de lı́neas y temperatura define la secuencia espectral como una secuencia de
temperaturas.
23
Capı́tulo 4
~r0 = ~r − ~v t, (4.1)
~r = ~r0 + ~v t0 , (4.2)
1
partı́cula, ~u, medida por los dos observadores,
d~x0 d~x dt d~x
~u0 = = 0 − ~v 0 = − ~v = ~u − ~v , (4.3)
dt0 dt dt dt
empleando dt0 = dt. Al derivar una segunda vez encontramos que la acela-
ración de la partı́cula es la misma, ~a0 = ~a, de acuerdo al principio de inercia
inscrito en la primera ley de Newton.
ct0 = A ct + B x , x0 = D ct + E x .
2
Figura 4.1: Dos observadores O y O’ con movimiento relativo ~v = βc x̂
viendo la propagación de una onda luminosa emitida en (x = 0, t = 0) y
(x0 = 0, t0 = 0).
ct0
γ −γβ 0 0 ct
x0 −γβ γ 0 0 x
= .
y0 0 0 1 0 y
z0 0 0 0 1 z
3
Figura 4.2: Marcos de referencia O y O’ con sistemas de ejes paralelos y en
movimiento relativo ~v = βc x̂
ct0
γ 0 −γβ 0 ct
x0 0 1 0 0 x
= .
y0 −γβ 0 γ 0 y
z0 0 0 0 1 z
4
que se representa también mediante la matriz
γ −γβx −γβy −γβz
−γβx 1 + (γ − 1)βx2 /β 2 (γ − 1)βx βy /β 2 (γ − 1)βx βz /β 2
.
−γβy (γ − 1)βy βx /β 2 1 + (γ − 1)βy2 /β 2 (γ − 1)βy βz /β 2
−γβz (γ − 1)βz βx /β 2 (γ − 1)βz βy /β 2 1 + (γ − 1)βz2 /β 2
(4.7)
Esta expresión representa cualquier movimiento rectilı́neo uniforme entre
sistemas con ejes paralelos. Las rotaciones entre ejes de coordenadas son
también transformaciones de Lorentz y la forma mas general de estas inclu-
ye combinaciones de rotaciones y movimientos rectilı́neos uniformes.
Jackson §11.7.
Transformación inversa
La transformación inversa se obtiene al invertir la matriz que lleva de (ct0 , x0 )
a (ct, x), o -mas simplemente- reemplazando β por −β. Empleando notación
matricial y omitiendo las coordenadas (y, z), la transformación inversa de
un movimiento a lo largo del eje x̂ es
! ! ! ! ! !
ct0 γ −γβ ct ct γ γβ ct0
= ⇐⇒ = .
x0 −γβ γ x x γβ γ x0
Suma de velocidades
Al combinar dos transformaciones sucesivas paralelas, con v1 = β1 c y v2 =
β2 c, obtenemos una transformación con
γ = γ1 γ2 (1 + β1 β2 ) , β = (β1 + β2 )/(1 + β1 β2 ) ,
Simultaneidad
De acuerdo al observador en el marco de referencia O, los eventos (t1 =
0, x1 = 0) y (t2 = 0, x2 = X 6= 0) son simultáneos, al ocurrir al mismo
tiempo. De la ecuación (4.5) vemos que los dos mismos eventos no son si-
multáneos para el observador en O’, quien los ubica en (t01 = 0, x01 = 0) y
(t02 = −γβ X/c, x02 = γ X). El concepto de simultaneidad depende del marco
de referencia del observador, y por tanto no es absoluto.
5
Dilatación temporal
Dos eventos que ocurren en un mismo lugar a tiempos distintos definen un
intervalo de tiempo. Podemos pensar que estos corresponden con dos tics
de un reloj. Los describimos como (t1 = 0, x1 = 0) y (t2 = T, x2 = 0) de
acuerdo al observador en O. En el marco de referencia O’ los mismos eventos
quedan descritos por (t01 = 0, x01 = 0) y (t02 = γT, x02 = −γβ cT ). La dila-
tación temporal se manifiesta en el marco de referencia O’, en movimiento
con respecto al reloj, en el cual estos estan separados por un lapso mayor
de tiempo: t02 − t01 ≥ t2 − t1 . Fraseado alternativamente, el lapso de tiempo
entre dos eventos es mı́nimo en el marco de referencia del reloj, en el cual
estos suceden en una sola posición. De aquı́ se deriva la definición del tiempo
propio, que es aquél medido en el marco de referencia donde el intervalo entre
eventos es puramente temporal (§4.1.2).
Contracción espacial
La contracción espacial, también contracción de Lorentz, no es evidente a
primera vista en las transformaciones de Lorentz. La medición de una ba-
rra en reposo de longitud L que realiza O puede describirse con dos eventos,
(t1 = 0, x1 = 0) y (t2 = 0, x2 = L) correspondientes a la medición simultánea
de los extremos de la barra. O’ ve la barra en movimiento y, según él, estos
dos eventos ocurren en (t01 = 0, x01 = 0) y (t02 = −γβ L/c, x02 = γ L). La pri-
mera impresión es que dado que x02 − x01 ≥ x2 − x1 , deberiamos hablar de
una expansión espacial, y no de una contracción.
La sutileza reside en que para obtener en O’ la longitud de la barra
en movimiento, la medición de la posición de ambos extremos debe hacer-
se simultáneamente. Como vimos antes, t01 6= t02 , por lo que la medición
simultánea de las dos puntas de la barra que hace O no son eventos si-
multáneos para O’. La solución consiste en dejar t2 libre, siendo que no se
requiere que la medición en O sea simultánea, y exigir t01 = t02 . Insistimos
que t2 6= t1 no representa ningún problema ya que, justamente, la barra está
en reposo para O, quien por tanto no requiere una medición simultánea.
Consideramos la medición de la barra en O con dos eventos (t1 = 0, x1 = 0)
y (t2 , x2 = L), los cuales corresponden para O’ en la medición simultánea
(t01 = 0, x01 = 0) y (t02 = 0, x02 = L0 ). La transformación de Lorentz relevante
implica, (
ct02 = ct2 − β L = 0 ⇒ t2 = βL/c
(4.8)
L = γ (L0 − β ct2 ) ⇒ L0 = L/γ ≤ L .
Al exigir la simultaneidad en O’ se obtiene la contracción de la barra, L0 ≤ L.
6
4.1.3. Cuadrivectores
El formalismo de la relatividad especial generaliza la mecánica de Newton
mediante el uso de cuadrivectores. Los cuadrivectores son generalizaciones de
vectores espaciales, añadiendo una componente temporal fı́sicamente asocia-
da, y permiten modificar las leyes de la mecánica y electrodinámica clásica
de manera consistente con los postulados de la relatividad. En particular,
el requerimiento de que la velocidad de la luz sea independiente del mar-
co de referencia se cumple por la condición de invariancia implı́cita en la
construcción de los cuadrivectores.
El cuadrivector elemental es el de tiempo y posición el cual ubica un even-
to en el espacio y el tiempo. Este cuadrivector tiene componentes x0 , x1 , x2 , x3 ,
abreviadas como xα . El ı́ndice α toma siempre el valor cero para la compo-
nente temporal y 1,2,3 para las tres componentes espaciales respectivas, las
cuales a veces se abrevian con algún ı́ndice latino, . El cuadrivector define
eventos mediante el tiempo y las coordenadas del mismo,
x0
ct !
α
x1 x ct
x = = = . (4.9)
x2 y ~x
x3 z
s2 ≡ −c2 t2 + x2 + y 2 + z 2 = 0 ,
7
Podemos describir la constancia de la velocidad de la luz como ||x||2 = 0.
La construcción de las transformaciones de Lorentz exige que la norma del
cuadrivector de posición sea la misma en diferentes marcos de referencia,
||x|| = ||x0 ||, independientemente de que su valor sea nulo o no. Por defini-
ción las transformaciones de Lorentz conservan la norma del cuadrivector
de posición y si podemos construir otros cuadrivectores cuya norma sea una
cantidad independiente del marco de referencia, estos conservaran su signi-
ficado fı́sico en todo marco de referencia inercial. Las transformaciones de
Lorentz se representan con matrices de la forma Λβα , quedando la transfor-
mación de O a O’ como
x0β = Λβα xα , (4.11)
donde introducimos la convención de que todo ı́ndice repetido como su-
perı́ndice y subı́ndice en una fórmula indica una sumatoria sobre él, es decir
x0β = Λβ0 x0 + Λβ1 x1 + Λβ2 x2 + Λβ3 x3 .
Los cuadrivectores con super-ı́ndice, aα , se denominan contravariantes y
a partir de ellos se construyen los cuadrivectores covariantes,
aα = ηαβ aβ ,
8
dx = dy = dz = 0, por lo que ds2 = −c2 dt2 . Definimos el tiempo propio, τ ,
como aquél medido por el observador “en reposo con un reloj”, es decir3
r
1p 1 dt
dτ ≡ −ds2 = dt2 − 2
(dx2 + dy 2 + dz 2 ) = . (4.14)
c c γ
El tiempo propio es una variable fı́sica unidimensional, es decir escalar, cuya
expresión es la misma en cualquier marco de referencia inercial; es por tanto
un invariante. Con él podemos construir cuadrivectores de norma invariante
al derivar con respecto a dτ . El ejemplo inmediato es el cuadrivector de
velocidades !
α dxα γc
u ≡ = , (4.15)
dτ γ ~v
9
4.1.4. Fotones, efecto Doppler y distorsión de haces
El fotón tiene asociado el cuadrivector de energı́a-momento,
! ! !
α E/c h̄ω/c h̄ω 1
k = = = , (4.18)
p~ h̄~k c k̂
10
la distribución uniforme del haz en ángulo sólido dΩ0 se transforma cómo
dΩ0 = dΩ/γ 2 (1 − β cos θ)2 . Al aumentar β el haz se cierra y la radiación se
corre al azul en la dirección de movimiento, dando lugar a una amplificación
de la intensidad observada, en el efecto conocido como “beaming”. (§4.2.3)
F αβ ≡ ∂ α Aβ − ∂ β Aα , (4.25)
11
representable como
0 Ex Ey Ez
−Ex 0 Bz −By
F αβ = . (4.26)
−Ey −Bz 0 Bx
−Ez By −Bx 0
~ ⊥ + ~v × B
~0 = E
E ~0 = γ E
~k , E ~⊥ ,
k ⊥
c
~v
~ 0 ~ ~ 0 ~ ~
Bk = Bk , B⊥ = γ B⊥ − × E⊥ . (4.29)
c
12
donde es necesario especificar el cuadrivector de fuerza, F, para definir el
problema. Desarrollamos el término de aceleración, notando que dγ/dt =
~˙ separando en componentes paralela y perpendicular a la velocidad,
γ 3 (β~ · β),
~a = ak β̂ + ~a⊥ , para obtener el cuadrivector de energı́a-momento
γ 4 m ~a · β~
γ d(E/c)/dt !
γ 4 m ak β
= = ,
γ 2 m ~a⊥ + γ 4 mak β̂
γ d~
p/dt γ 2 m ~a + γ 2 β~ (~a · β)
~
(4.31)
donde E es la energı́a de la partı́cula y
13
acuerdo a ~re = x̂ vt, visto por un observador en reposo en ~r = ŷ b. En un
marco de referencia en reposo con el electrón tenemos
0 0 0
~ 0 (~r0 , t0 ) = −e ~r = −e x x̂ + y ŷ = −e (−γvtx̂ + bŷ) ,
E ~0 = 0,
B
|~r0 |3 (x02 + y 02 )3/2 [b2 + γ 2 v 2 t2 ]3/2
que con las relaciones (4.29) da,
ρ(~r, t) = −e δ [~r − ~re (t)] , ~(~r, t) = −e ~ve (t) δ [~r − ~re (t)] . (4.35)
14
Figura 4.3: Electrón en movi-
miento arbitrario visto desde
el marco de referencia O. El
vector ~r denota la posición del
observador en el instante t y
~re (te ) la del electrón en el ins-
tante te . El vector R ~ une am-
bos eventos.
−c2 (t − te )2 + |~r(t) − ~re (te )|2 = 0 ⇒ te = t − |~r(t) − ~re (te )|/c . (4.37)
de donde φ0 (~r0 , t0 ) puede escribirse también como φ0 = −e/c (t0 − t0e ). Hacien-
do el producto escalar de Rα con el cuadrivector de velocidad de la partı́cula
medido en t = te , obtenemos el escalar invariante:
15
Estos potenciales son invariantes y toman la forma (4.36) en el marco de
referencia en reposo con el electrón en t0e (~ve0 (t0e ) = 0). La expresión (4.38)
corresponde a los potenciales de Liénard-Wiechert
−e ~ r, t) = −e β~e
φ(~r, t) = , A(~ . (4.39)
R 1 − R̂ · β~e R 1 − R̂ · β~e
te te
~˙
~ ~ × β)
" #
~ r, t) = e (R̂ − β) e R̂ × ((R̂ − β)
E(~ +
R2 γ 2 (1 − R̂ · β)
~ 3
te
Rc ~ 3
(1 − R̂ · β)
te
~ r, t) = R̂ × E(~
B(~ ~ r, t) , (4.40)
16
Figura 4.4: Origen de referencias O elegido dentro de la zona donde se mueve
la carga y tal que re r. El observador está indicado del lado derecho, y el
~ va de la posición de la carga en te a la del observador en t.
vector R
˙ ˙ ˙
dP (t) e2 a2 1 2(R̂ · β̂)(β~ · β̂) (1 − β 2 )(R̂ · β̂)2
= + − ,
dΩ 4πc3 (1 − R̂ · β)
~ 4 (1 − R̂ · β)~ 5 ~ 6
(1 − R̂ · β)
te
(4.44)
˙~
siendo ~a = βc la aceleración de la carga. La potencia radiada por la carga
en su propio marco de referencia se describe con el tiempo te = t − R/c, y
está dada por
dP (te ) dP (t)
= (1 − R̂ · β) . (4.45)
dΩ dΩ
17
Figura 4.5: Patron de emisión para un mo-
vimiento no relativista, dado por la fórmu-
la dipolar, ∝ sin2 θ. Las flechas indican la
dirección del vector de aceleración, con res-
pecto a la cual se mide el ángulo θ.
Caso no relativista
Para una carga en movimiento no-relativista (β 1), los campos de radia-
ción están dados por:
dP (t) e2 a2 e 2 a2
˙
≈ 1 − (R̂ · β̂)2 = sin2 θ , (4.47)
dΩ 4πc3 4πc3
siendo θ el ángulo entre la lı́nea de visión y la aceleración; el patrón de
emisión es independiente de la dirección del vector de velocidad. Al integrar
sobre dΩ hallamos la potencia total radiada en función del tiempo:
2 e2 a2
P (t) = . (4.48)
3 c3
Tres puntos a notar: (i) una carga acelerada radı́a proporcionalmente a a2 e2 ;
(ii) la emisión tiene un patrón de distribución geométrica ∝ sin2 θ (fig. 4.4);
(iii) siendo E ~ rad ∝ R̂ × (R̂ × ~a), si ~a tiene una dirección fija (movimiento
rectilı́neo), entonces Ê también es fija, y la polarización es 100 % lineal.
18
Caso relativista
La expresión de Larmor no relativista (4.48) se generaliza directamente a
una expresión covariante:
2e2
P = 3 aµ aµ , (4.49)
3c
que en un marco de referencia general tiene la forma
!2
~
2e2 2 d(γ β)
dγ
2
2e2 4 2
6 2
P (t) = γ − = γ a ⊥ + γ ak , (4.50)
3c dt dt 3c3
19
Figura 4.6: Patrones de emisión para una carga en movimiento relativista:
izquierda: ~a k ~v , para β ≈ 0 y β ≈ 0.866. La figura de la derecha del recuadro
está reducida un factor de 100. derecha: ~a ⊥ ~v para β ≈ 0 y β = 0.5. La
figura de la derecha del recuadro está reducida un factor de 3.
~
con RE(t) ~
y RE(ω) relacionados mediante la transformada de Fourier,
Z +∞ Z +∞
~ 1 ~ iωt ~ 1 ~ −iωt
RE(ω) =√ RE(t)e dt, RE(t) =√ RE(ω)e dω .
2π −∞ 2π −∞
(4.53)
20
~ real se cumple E(−ω)
Para E ~ =E ~ ∗ (ω).
De la expresión (4.41) para el campo de radiación se tiene
2
~˙
dE e2
Z +∞
R̂ × {( R̂ − ~ × β}
β)
eiωt dt .
= 2
dΩdω 4π c −∞ ~ 3
(1 − R̂ · β)
te
Para avanzar usamos la aproximación t ' te + r/c − r̂ ·~re (te )/c para cambiar
la variable de integración, notando dt = (1 − β~ · r̂)dte ,
2
~˙
dP (ω) dE e2
Z +∞
r̂ × {(r̂ − ~ × β}
β)
iω(t −r̂·~
r /c)
e e
= = 2 e dte ,
dΩ dΩdω 4π c −∞ ~ 2
(1 − r̂ · β)
~˙
~ × β} ~
" #
r̂ × {(r̂ − β) =
d r̂ × (r̂ × β)
,
~ 2
(1 − r̂ · β) dt 1 − r̂ · β~
que permite una integración por partes y reducir la expresión para la distri-
bución espectral y angular a
2
e2 ω 2
Z +∞
dE ~ iω(t−r̂·~re (t)/c) dt ,
= r̂ × (r̂ × β)e (4.54)
dΩdω 4π 2 c −∞
21
4.3. Sistemas de cargas en movimiento
4.3.1. Aproximación dipolar
La aproximación no relativista de la ecuación de Larmor (4.48) puede gene-
ralizarse a un sistema de cargas como
~¨ 2
2|d(t)| ~ =
X
P (t) = , con d(t) qı~rı (t), (4.57)
3c3 ı
dP (ω) dE ω4 ~ ~
4ω 4 |d(ω)| 2
= = |d(ω)|2 sin2 θ , P (ω) = . (4.60)
dΩ dω dΩ 2πc3 3c 3
1 α (~x, t0 )
Z
α
A (~r, t) = δ(t0 − t + |~r − ~x|/c) d3 x dt0 . (4.61)
c |~r − ~x|
22
Aquı́ ~x denota la posición de cada elemento de la distribución de carga
y corriente, siendo t0 el tiempo retardado correspondiente a ~x. Nuevamen-
te, empleamos la aproximación del campo lejano, en la cual la región de
densidad de carga y corriente está localizada y tiene dimensiones mucho
menores que la distancia a ella. El origen de coordenadas se sitúa en el cen-
tro de la distribución de cargas y corrientes, lo que nos permite aproximar
R = |~r − ~x| ≈ r − ~r · ~x en el tiempo retardado, o simplemente R ≈ r en el
término inversamente lineal,
1
Z
α
A (~r, t) = α (~x, t − r/c + r̂ · ~x/c) d3 x .
rc
Afuera de la región de carga y corriente, el campo de radiación se comporta
como una superposición de ondas planas, con E ~ =B~ × r̂. Por tanto, basta
determinar A ~ para calcular ambos campos, siendo B ~ = ∇ × A. ~ Se puede
−2
verificar que despreciando términos en r , se tiene B ~ = (1/c)A ~˙ × r̂, y por
tanto E~ = (1/c)(A~˙ × r̂) × r̂.
~
Al considerar la transformada de Fourier de A(t), obtenemos:
−iωr/c
~ r, ω) = e
Z
A(~ ~(~x, ω)eiωr̂·~x/c d3 x , (4.62)
rc
~
con B(ω) ~
= iω r̂ × A(ω)/c, ~
y E(ω) ~
= −iω r̂ × (r̂ × A(ω))/c. El desarrollo en
serie de la exponencial lleva directamente a la expansión multipolar:
∞
e−iωr/c X 1
Z
~
A(~r, ω) = (iωr̂ · ~x/c)n ~(~x, ω)d3 x , (4.63)
rc n=0 n!
e−iωr/c
Z
~ dip (~r, ω) =
A ~(~x, ω) d3 x , (4.64)
rc
23
e−iωr/c iω
Z
~ cuad (~r, ω) =
A r̂ · ~x ~ d3 x , (4.65)
rc c
e−iωr/c iω 2
Z
~ oct (~r, ω) =
A (r̂ · ~x)2 ~ d3 x . (4.66)
rc c
Utilizando los términos dipolar y cuadrupolar puede escribirse el potencial
en términos de derivadas temporales de los momentos de la distribución de
carga y corriente como
˙
~ d~ ~˙ × r̂ r̂ · Q̈
m
A(~r, t) = + + + ... , (4.67)
rc rc 6rc2
donde m
~ es el momento dipolar magnético y Q el tensor (3D) de momento
cuadrupolar eléctrico, dados por:
1 X X
m
~ = qi ~xi × ~vi , Qab = qi 3xa xb − |~x|2 δab , (4.68)
2c i i
donde Q ~ = r̂ · Q = P
a ra Qab . La distribución angular de la potencia radiada
está dada por
dP c ~ 2
= |rB| .
dΩ 4π
La expansión multipolar resulta en los términos correspondientes en la ecua-
ción de Larmor,
¨~ 2 ...
2|d| ~¨ 2 || Q ||2
2|m|
P (t) = + + + ... , (4.71)
3c3 3c3 180 c5
Las dependencias en frecuencias se obtienen al modificar las derivadas tem-
porales por términos iω. (Shu 13, 14; Rybicki §3.1 a 3.3).
24
4.4. Bremsstrahlung
La radiación debida a la aceleración entre cargas se conoce como bremss-
trahlung, radiación de frenado o radiación libre-libre. El problema consiste
en encontrar las caracterı́sticas de la emisión para dos cargas que interaccio-
nan a través de la fuerza Coulombiana. El caso se trata de manera clásica
si la energı́a radiada en la interacción es mucho menor que la energı́a inicial
del sistema; el caso opuesto debe tratarse con un formalismo cuántico. Se
puede emplear la aproximación dipolar en el caso de movimientos no rela-
tivistas (v c). Nótese que el dipolo de dos partı́culas de misma carga y
misma masa es nulo, y por tanto no hay radiación dentro de la aproximación
dipolar. En consecuencia, dentro de la aproximación dipolar, la interacción
entre dos electrones no produce radiación.
25
En el caso de una deflexión tipo bremsstrahlung, es útil plantear el pro-
blema en términos de la velocidad inicial, v0 , y el parámetro de impacto, b,
parámetros equivalentes a E = mv02 /2, y L = mbv0 :
q √ p
v0 = 2E/m = a/τ , b = L/ 2mE = a ε2 − 1 . (4.76)
26
Figura 4.7: Trayectoria y aceleración en una interaccin con deflexión pe-
queña.
16Z 2 e6
Z bmax
dE(v) dE bmax
= ni ne v (v, b) 2πb db = ne ni ln , (4.83)
dtdV dω bmin dω 3c3 m2 v bmin
27
donde el rango finito de valores de b asegura la convergencia de la integral.
El lı́mite superior puede tomarse como bmax ≈ v/ω, generalmente menor
−1/3
a ne , excepto en casos extremos (T ∼ 10 K, ω ∼ 10 MHz). Por otro lado,
tenemos ya identificados dos lı́mites inferiores: el de la validez de la aproxi-
(1)
mación de deflexiones pequeñas, b bmin ≡ Ze2 /mv 2 ; y el lı́mite cuántico,
(2)
b bmin ≡ h̄/mv. El lı́mite cuántico es superior al de deflexiones pequeñas
si v/c >∼ Zα, donde α ' 1/137 es la constante de estructura fina. En es-
te caso, podemos tomar bmin = Ze2 /mv 2 dentro de la expresión; en caso
contrario, debemos considerar de alguna forma las correcciones cuánticas.
Finalmente, notamos que la condición de que energı́a radiada sea mucho
menor que E viene siendo (dE/dω)(v/b) mv 2 /2, lo cual se traduce en
1/3
Ze2 v 16
b . (4.84)
mv 2 c 3πZ
Este caso se cumple automáticamente al estar considerando deflexiones pe-
queñas en un movimiento no relativista.
El resultado exacto puede referirse a la aproximación aquı́ tratada defi-
niendo el factor de Gaunt, gf f (v, ω), tal que
dE(v) 16πZ 2 e6
= √ 3 2 ne ni gf f (v, ω) . (4.85)
dtdV dω 3 3c m v
√
En el caso considerado aquı́, gf f (v, ω) = ( 3/π) ln(bmax /bmin ).
28
Figura 4.8: Movimiento de una
carga atraı́da por un campo cen-
tral Coulombiano de acuerdo con
la parametrización (4.86). Se tomó
ε = 2 y los puntos están separados
por ∆ψ = 0.1, cubriendo un ran-
go en ψ de ±3, y un movimiento
de abajo hacia arriba. El centro
de atracción está en el origen,√ el
parámetro de impacto es b = 3.
Se muestra el punto mas cercano
de la ası́ntota ψ = +∞, separado
del origen una distancia b.
~
d(ω) = m~r(ω)(q1 /m1 − q2 /m2 ), siendo m la masa reducida y anulándose
d~ para una interacción entre partı́culas idénticas o con mismo cociente carga
masa. Aquı́ usamos (q1 /m1 −q2 /m2 ) ≈ −e/m, y calculamos la transformada
de Fourier de {x, y} empleando la parametrización, de donde se obtiene
√
a (1)0 a ε2 − 1 (1)
x(ω) = − Hiωτ (iωτ ε), y(ω) = Hiωτ (iωτ ε),
2ω 2ω ε
siendo Z +∞
1
Hp(1) = epψ−ix sinh ψ dψ ,
iπ ∞
dE 4π 2 Z 2 e6 0
= ωτ |Hiωτ (iωτ )|Hiωτ (iωτ ) . (4.88)
dV dωdt 3c3 m2 v
29
Figura 4.9: Movimiento
de una carga atraı́da por
un campo central Coulom-
biano. En el panel inferior
se muestran las posiciones
como función del tiempo
{x(t), y(t)} y en el panel
superior las aceleraciones
{ẍ(t), ÿ(t)}, para el mo-
vimiento ilustrado en la
figura 4.8 (ε = 2).
Este resultado concuerda con el resultado para ángulos pequeños (4.85), con
un factor de Gaunt
√ !
π 3 0
gf f (v, ω) = ωτ |Hiωτ (iωτ )|Hiωτ (iωτ ) . (4.89)
4
30
donde ḡf f es un promedio sobre velocidades del factor de Gaunt y ν la
emisividad del proceso, en este caso isotrópico. Integrando sobre frecuencias
encontramos la potencia radiada por unidad de volumen,
s
dE(T ) 32πe6 2πkT
= ne ni Z 2 ḡB
dtdV 3hmc3 3m
= 1.426 × 10−27 erg s−1 cm−3 T 1/2 ne ni Z 2 ḡB , (4.92)
4e6
r
2π
αν = ne ni Z 2 T −1/2 ν −3 1 − e−hν/kT ḡf f ,
3mhc 3km
= 3.7 × 108 cm−1 ne ni Z 2 T −1/2 ν −3 1 − e−hν/kT ḡf f . (4.93)
31
de igual manera. En el caso relativista, el tratamiento no cuántico implica
deflexiones pequeñas, ya que
dE 16 ni Z 2 e6 −1/3
= ln 192 Z , (4.96)
dtdV dω 3m2 c4
donde se pierde el factor γ al retornar al marco del observador. El factor
numérico incluye la constante de estructura fina y v/c → 1. La pérdida total
de energı́a resulta de integrar sobre frecuencias, desde cero hasta E/h̄,
dE 16 Z 2 e6 E −1/3
= ln 192 Z ni dV ,
dt 3h̄m2 c4
32
que corresponde a una pérdida exponencial de energı́a. Esta se puede expre-
sar en términos de la longitud o de la densidad de columna, X,
dE E 716 g cm−2 A
=− , donde X0 = , (4.97)
dX x0 Z(Z + 1.3) ln(183 Z −1/3 ) + 1/8
con A la masa atómica del material involucrado. Esta expresión incorpora
un cálculo más preciso realizado originalmente por Bethe y Heitler. Para
el aire X0 = 36.5 g cm−2 . Un electrón de alta energı́a en la atmósfera, de
espesor total un poco superior a 1000 g cm−2 emite fotones de alta energı́a
al atravesar una densidad de columna del orden de X0 .
33
Larmor, r` ≡ v sin α/ωB . El movimiento del electrón es una superposición de
un movimiento rectilı́neo uniforme en el eje ẑ, con velocidad vk = v cos α, y
un movimiento circular uniforme en el plano xy, con velocidad v⊥ = v sin α,
frecuencia ωB y una dirección de giro determinada por la dirección de ω ~ B.
2e2 µ 2e2 4 2 2
P (t) = 3
aµ a = 3
γ a⊥ = re2 c γ 2 β⊥
2 2
B , (4.100)
3c 3c 3
donde re = e2 /mc2 ' 2.8178 × 10−13 cm, es el radio clásico del electrón.
Promediando la componente de velocidad perpendicular al campo magnéti-
co, ∝ sin α, sobre ángulo sólido, tenemos < β 2 sin2 α2 >= 2β 2 /3, con lo que
podemos expresar la potencia promedio radiada por un electrón como
4
P (t) = σT (γ 2 − 1) uB c . (4.101)
3
σT = (8π/3)re2 es la sección eficaz de Thompson, y uB = B 2 /8π es la densi-
dad de energı́a del campo magnético. Esta expresión es análoga a la potencia
radiada por efecto Compton por un electrón en un campo de radiación de
densidad de energı́a urad .
En el tratamiento cinemático se despreciaron las pérdidas de energı́a al
suponer E = γmc2 constante. Un electrón en un campo magnético pierde
energı́a en una escala de tiempo
−1 !
1 dγ 3 mc γ
τ≡ = 2 2 2
.
γ dt 2 re B⊥ γ −1
34
6
~ = B ẑ
B
5
3
θ
2 x̂
0
0 2 4 6 8 10
˙
donde R̂ ≈ r̂ y el término β~ · β~ en la expresión (4.44) se anula. El tiempo
retardado está dado por
dP (t) 2 β 2 (β − cos ω t )2
e2 ωB B e
= , (4.103)
dΩ 4πc (1 − β cos ωB te )6
35
Figura 4.11: Pulso de la emi-
sión sincrotrón visto en térmi-
nos del tiempo retardado, te y
del tiempo medido por el ob-
servador, t. En este ejemplo se
tomó β = 0.8 ⇒ γ ' 1.67.
El pulso en te tiene duración
∆te ∝ 1/γ, mientras que en t
tiene duración ∆te ∝ 1/γ 3 . Pa-
ra valores grandes de γ, los pul-
sos aparecen como lı́neas.
y aproximando,
2 2
dP (t) e2 ωB
2 1 − γ 2 ωB
2 t2
e 4e2 ωB
2 6 2 2
8 1 − γ ωB t
≈ 16γ 8 2 t 2 )6 ≈ γ 2 t2 )6 , (4.104)
dΩ 4πc (1 + γ 2 ωB e πc (1 + γ 6 ωB
para ωB te 1. El comportamiento del pulso de sincrotrón se ilustra en la
figura 4.11. La emisión está formada por pulsos de ancho ∆t ≈ 2π/γ 3 ωB ,
repitiéndose con periodo T = 2π/ωB . Esto resulta en un espectro de emisión
que cubre frecuencias desde ωB hasta γ 3 ωB . Formalmente el espectro es
discreto y la emisión se da en frecuencias ωn = nωB con n = 1 → γ 3 .
e2 ω 2 +∞ 2
Z
dP (ω) dE ~ iω(t−r̂·~re (t)/c) dt
= = r̂ × (r̂ × β)e
dΩ dωdΩ 4π 2 c −∞
2 2
e ω 2 2
= 2
β ê⊥ A⊥ (θ, ω) − êk Ak (θ, ω) , (4.105)
4π c
36
siendo Ak , A⊥ las amplitudes de las componentes de polarización orientadas
a lo largo de los vectores êk = ŷ, ê⊥ = x̂ sin θ − ẑ cos θ, respectivamente,
ambos perpendiculares a r̂. Las componentes de A(θ, ~ ω) están dadas por,
Z +∞
Ak (θ, ω) = sin ωB t · exp {iωt(1 − β cos θ sin(ωB t)/ωB t)} dt
−∞
Z +∞ n o
≈ ωB t · exp iωt 1/2γ 2 + θ2 /2 + ωB
2 2
t /6 dt,
−∞
1 1
≈ ωB 2
+ θ2 √ K2/3 (x),
γ 3
Z +∞
A⊥ (θ, ω) = sin θ exp {iωt(1 − β cos θ sin(ωB t)/ωB t)} dt
−∞
Z +∞ n o
≈ θ exp iωt 1/2γ 2 + θ2 /2 + ωB
2 2
t /6 dt
−∞
1/2
1 1
≈ ωB θ + θ2 √ K1/3 (x), (4.106)
γ2 3
ωc ≡ 3γ 3 ωB /2 .
37
Figura 4.12: Izquierda: distribución angular de la emisión sincrotrón. Dere-
cha: distribución espectral de la emisión sincrotrón dada por las funciones
F (x) -lı́nea sólida- y G(x) -lı́nea punteada.
38
Esta fórmula describe la emisión por un sólo electrón. En general tendre-
mos una población de electrones y su distribución de energı́a determinará el
espectro total y la polarización correspondiente.
(Rybicki §6; Jackson §14)
4.5.4. Auto-absorción
De acuerdo con el formalismo de transferencia radiativa, a un proceso de
emisión le corresponde un proceso de absorción recı́proco. Para un proceso
en equilibrio, el coeficiente de absorción, αν se deduce de la emisividad, jν ,
39
a través de la ley de Kirchhoff, jν = αν Bν (T ), con Bν la función de Planck.
El caso de radiación sincrotrón no corresponde al de un emisor en equilibrio,
pero es posible obtener el coeficiente de absorción a través de jν = αν Sν .
Consideremos la función fuente como el producto del número de estados por
unidad de volumen de espacio fase (2ν 2 /c3 ), por el flujo de energı́a, dado
por la velocidad de la radiación (c) y la energı́a promedio de cada fotón
Sν = 2 ν 2 /c3 c<E > .
~ + ~v × B
F~ = −e E ~ .
c
40
La dispersión de Thomson se restringe al caso no relativista (v c), en el
que se puede despreciar el efecto del campo magnético, menor al eléctrico
por un factor ∼ v/c. El movimiento del electrón está dado por
¨ e = −eE ε̂ ei(kz−ωt) ,
m~x
para una onda monocromática propagándose con k̂ = ẑ, y ε̂ el vector de
polarización del campo eléctrico. Si ż(0) = 0, la partı́cula restringe su mo-
vimiento al plano xy; y el término eikz representa una fase, al ser z fija.
La potencia radiada, incluyendo el patrón de emisión, se estima de la
aproximación dipolar para la ecuación de Larmor
¨~ 2
dP |r̂ × (r̂ × d)|
= , (4.117)
dΩ 4πc3
siendo d~ = −e~xe , el momento dipolar para el electrón situado en ~xe , y r̂ el
vector unitario que une al electrón con el observador. La potencia promedio
radiada por el electrón por unidad de ángulo sólido queda como:
!2
dP c e2
= |E|2 1 − |r̂ · ε̂0 |2 . (4.118)
dΩ 8π mc2
Podemos cuantificar el efecto de la dispersión Thomson en términos de la
comparación entre la “entrada” (onda plana) y la “salida” (emisión dipolar),
cuantificada por la sección eficaz total, σ, o diferencial, dσ/dΩ:
!2
dP dσ dσ 1 e2
=< S > ⇒ = 1 − |r̂ · ε̂0 |2 . (4.119)
dΩ dΩ dΩ 2 mc2
Usando r̂ = ẑ cos θ + (x̂ cos φ + ŷ sin φ) sin θ , para polarización lineal (ε̂0 =
x̂), obtenemos un patrón angular ∝ (1 − sin2 θ cos2 φ). Para una onda con
polarización circular o para una onda no polarizada, se obtiene
!2
dσ e2 1
= (1 + cos2 θ) . (4.120)
dΩ mc2 2
La integración sobre ángulo sólido nos da la sección eficaz de Thomson,
!2
8π e2 8π 2
σT = = r , (4.121)
3 mc2 3 e
siendo re ≡ e2 /mc2 ≈ 2.82 × 10−13 cm el radio clásico del electrón.
No hay alteración del espectro incidente: la luz dispersada tiene la misma
frecuencia que la onda incidente. Es al considerar el caso relativista -el efecto
Compton- que obtenemos un cambio en la frecuencia de la radiación.
41
4.6.2. El efecto Compton
El efecto Compton es la dispersión de un electrón por un fotón. En
este caso el fotón es tratado como una partı́cula y la interacción como una
colisión,
γe → γe . (4.122)
Este tratamiento es mas apropiado que el de dispersión Thomson para foto-
nes con energı́as comparables con la de la masa del electrón en reposo, mc2 .
En este caso se considera la colisión en el marco de referencia del electrón,
empleando la conservación del cuadrivector de energı́a y momento,
El electrón está inicialmente en reposo, por lo que pe(i) = (mc; 0), mientras
que pe(f ) = (γmc; γmβc). ~ El fotón tiene pγ(i) = (h̄ω0 /c)(1; k̂0 ) antes de la
colisón, y pγ(f ) = (h̄ω1 /c)(1; k̂1 ) después. Se puede manipular las ecuacio-
nes para obtener, independientemente del movimiento final del electrón, la
ecuación de Compton que relaciona las energı́as inicial y final del fotón:
h̄ω0
h̄ω1 = . (4.124)
1 + (h̄ω0 /mc2 ) 1 − k̂0 · k̂1
42
Para ω1 = ω0 recobramos la sección eficaz diferencial de Thomson. Pode-
mos expresar dσ/dΩ como función de cos θ y x = h̄ω0 /mc2 , substituyendo
ec.(4.124) en ec.(4.125), y de ahı́ ver el comportamiento de dσ/dΩ(θ), que se
ilustra en la figura 4.13. Al aumentar la energı́a del fotón, el efecto del mo-
mento que recibe el electrón se manifiesta en que dσ/dΩ es menos isotrópica,
favoreciendo la dispersión con θ pequeñas.
3 1+x 2x(1 + x) 1 1 + 3x
σ = σT − ln(1 + 2x) + ln(1 + 2x) − ,
4 x3 1 + 2x 2x (1 + 2x)2
(4.126)
mas sencillamente expresada como σ = σT · f (h̄ω0 /mc2 ). En el régimen no
relativista, x 1, tenemos
26
σ ≈ σT 1 − 2x + x2 + . . . , (4.127)
5
43
mientras que en el extremo relativista, x 1, la sección eficaz tiende lenta-
mente a cero:
3 1
σ ≈ σT + ln 2x . (4.128)
8x 2
44
Este resultado, junto con la conservación de energı́a, implican
de donde
γ02 ω0 (1 − β~0 · k̂0 )(1 + β~0 · k̂10 )
ω1 ' , (4.136)
1 + 2γ0 ω0 (1 − β~0 · k̂0 )
despreciando términos en ω0 /γ0 1. El resultado final depende de si
γ0 ω0 (1 − β~0 · k̂0 ) -la frecuencia del fotón en el centro de momento- es mucho
mayor o mucho menor que uno:
4
Esta expresión corresponde con la fórmula de Klein-Nishima, pero en un marco de
referencia distinto.
45
γ0 ω0 (1 − β~0 · k̂0 ) 1:
Despreciando el término en (4.136), obtenemos
γ0 ω0 (1 − β~0 · k̂0 ) 1:
El término domina en el denominador y se elimina, quedando
1
ω1 ' γ0 (1 + β~0 · k̂10 ).
2
Esta colisión corresponde al caso ω 1 en el centro de momento, por
lo que el “rebote” es dominante en la sección eficaz y k̂10 → −k̂00 . Un
análisis de la transformación de k̂00 indica que
1 + 2 −→ 3 + 4 , (4.137)
q1 + q2 + q3 + q4 = 0 , (4.138)
46
Figura 4.14: Proceso básico de interac-
ción electrodinámica: la absorción de un
fotón (“viborita”) por un electrón (fle-
cha). Rotando este diagrama podemos
representar la creación de un par e+ e−
por un fotón, con la convención de que
una flecha yendo en reversa representa
una anti-partı́cula o la aniquilación e+ e−
en un fotón. Estos procesos no pueden
conservar energı́a y momento simultánea-
mente.
Otros ejemplos son la interacción entre dos electrones cuyos canales son
e− e− → e− e− , e+ e− → e+ e− , e+ e+ → e+ e+ ,
e − µ− → e − µ − , e+ e− → µ+ µ− , µ+ µ− → e+ e− .
47
por un electrón (representado por la flecha yendo de izquierda a derecha).
Usamos la convención de que el proceso (o el tiempo) procede de izquierda a
derecha. Una flecha yendo en sentido opuesto representa una antipartı́cula.
Un detalle importante es que las interacciones representadas por este diagra-
ma no pueden darse de manera aislada ya que no pueden conservar energı́a
y momento simultáneamente. Por ejemplo en el caso del diagrama (4.14), si
~ se deberı́a cumplir
el electrón incidente tiene velocidad βc,
Figura 4.15: Los tres canales de un mismo proceso: a la izquierda dos dia-
gramas representando la dispersión Compton; en el centro la producción de
un par por dos fotones y a la derecha la aniquilación una partı́cula con su
anti-partı́cula para producir dos fotones.
48
la reacción corresponden con puntas “sueltas” mientras que las partı́culas
que unen dos vértices son denominadas partı́culas “virtuales”.
Los diagramas de Feynman son de utilidad en el cálculo de las secciones
eficaces en cuestión, habiendo una relación cercana entre las secciones efica-
ces de los distintos canales de una reacción genérica. Al ser estos procesos
manifestaciones diversas de un mismo tipo de proceso, las secciones eficaces
son tı́picamente del mismo orden (∼ re2 ) y su cálculo se hace bajo el mismo
formalismo de la electrodinámica cuántica.
(Griffiths; Berestetskii)
Producción de pares
La producción de un par electrón-positrón por dos fotones, denotada
por la reacción γγ → e+ e− , corresponde al diagrama central mostrado en la
figura 4.15. Empleando unidades tales que h̄ = 1, mc2 = 1, la conservación
de energı́a-momento se escribe como
! !
ω0 + ω1 γ0 + γ1
= ,
ω0 k̂0 + ω1 k̂1 γ0 β~0 + γ1 β~1
con la invariancia
49
Figura 4.16: Sección eficaz diferencial (izquierda) y total (derecha) para la
producción de pares γγ → e+ e− . La sección eficaz diferencial adquiere un
comportamiento bi-polar a medida que ω aumenta. La sección eficaz total
es nula si ω ≤ 1.0 y máxima para ω ≈ 1.25.
50
Figura 4.17: Sección eficaz para la aniquilación de pares e+ e− → γγ, en
términos de la velocidad o del factor de Lorentz de los electrones en el
marco de referencia de momento nulo. La sección eficaz diverge para para
electrones lentos, γ → 1.
La sección eficaz diverge para γ → 1, como puede verse en la fig. 4.17, con
σ → πre2 /β. Este comportamiento favorece la aniquilación de pares con baja
velocidad relativa y, por tanto, la formación de una lı́nea de aniquilación
delgada a mc2 ' 0.511MeV, observada por OSSE en el centro de la Galaxia.
Cuando se encuentran dos electrones con espines paralelos, el momento
angular de espı́n es h̄ en total. Dado que un fotón transporta un momen-
to angular igual a h̄, no es posible la aniquilación del par produciendo dos
51
fotones. Si los electrones son suficientemente lentos es posible entonces la for-
mación de un sistema ligado análogo al átomo de hidrógeno conocido como
positronio. Dado que, a diferencia de lo que sucede con el hidrógeno, las dos
partı́culas tienen misma masa, la masa reducida del sistema es µ = me /2
y los niveles de energı́a estan a la mitad que los del hidrógeno (es decir
En ' −6.8 eV/n2 ). Una vez formado, el positronio puede adquirir una con-
figuración de momento angular nulo y decaer en dos fotones de 0.511 MeV.
También es posible el decaimiento en tres fotones, e+ e− → γγγ, cuyas ener-
gı́as son (en promedio) menores a 0.511 MeV, pero pueden tomar cualquier
valor. Este tipo de decaimiento se manifiesta en la emisión, observada por
OSSE en el centro de la Galaxia, de un continuo a energı́as menores a la
de la lı́nea de aniquilación. La detección de este continuo de positronio es
evidencia de la presencia de e+ e− lentos en el centrode la Galaxia.
e− → γe− , γ → e− e+ , e+ e− → γ,
E1 B
E1 αB 2 λ3c > 2 2
∼ (mc ) ⇒ x≡ >
∼ 1.
mc2 Bc
donde Bc ≡ m2 c3 /eh̄ ' 4 × 1013 Gauss. Una reacción de tres vértices rele-
vante en la vecindad de estrellas de neutrones con B ∼ Bc es el “photon
splitting”, γ → γγ.
La reacción e− → γe− en presencia de un campo electrostático corres-
ponde al bremsstrahlung. La consideracón para el canal de producción de
un par e− e+ en la vecindad de un núcleo cargada viene dada en términos
52
del parámetro de impacto y la carga del núcleo,
Ze2
∼ mc2 ⇒ σ ∼ αb2 ∼ αZ 2 re2 .
b
La cálculo de la sección eficaz para la producción de pares γ → e− e+ da
28 2h̄ω 109
σ = α Z 2 re2 log − para h̄ω/mc2 1 .
9 mc2 42
53
Figura 4.18: Ilustración de
una cascada electromagnéti-
ca, con los fotones en azul y
los e± en rojo. El desarrollo
de la cascada es el que corres-
ponde a la sección eficaz de
cada proceso, ignorando las
fluctuaciones inherentes a la
cascada.
54
Capı́tulo 5
d2~r
m = F~ (~r, ~r˙ , t) . (5.1)
dt2
1
Figura 5.1: La posición de la masa del
péndulo rı́gido tiene tres coordenadas
fı́sicas (x, y, z) pero, por las constric-
ciones del punto fijo de soporte y de
la rigidez del péndulo, ` =constante,
es más conveniente su descripción en
términos de una sola coordenada ge-
neralizada, θ(t).
equivalentes a (5.1).
2
5.1.2. Cuatro ejemplos clásicos
Ej.1. Partı́cula libre
El Hamiltoniano de una partı́cula libre en tres dimensiones es
x = ` sin θ, y = ` − ` cos θ,
p2θ
H(θ, pθ ) = + mg`(1 − cos θ) ,
2m `2
con las ecuaciones de movimiento
3
del cual se obtienen las ecuaciones de movimiento
ẋ = p/m, ṗ = −kx .
2
p se obtiene la ecuación del oscilador armónico, ẍ − ω x = 0,
Derivando ẋ
donde ω = k/m. La solución general es de la forma,
x(t) = Aeiωt + Be−iωt , p(t) = mẋ = mω Aeiωt − Be−iωt ,
p2y
! ! !
p2x 1 1 p2z 1
H= + kx x2 + + ky y 2 + + kz z 2 .
2m 2 2m 2 2m 2
con las seis constantes determinadas de las condiciones iniciales, ~r(0) y p~(0).
En el caso particular kx = ky = kz ≡ k el problema adquiere simetrı́a esférica
y puede tratarse, tanto en el formalismo clásico como cuántico, como un caso
particular del potencial central U (r) = kr2 /2.
Expresando la velocidad como ~r˙ = ṙr̂ + rθ̇θ̂ + r sin θφ̇φ̂ tenemos la energı́a
cinética
1 1 1
Ecin = mṙ2 + mr2 θ̇2 + mr2 sin θ2 φ̇2 ,
2 2 2
y podemos aplicar la definición (5.2) para obtener el momento generalizado,
4
y la expresión del Hamiltoniano,
p2r p2θ p2φ
H= + + + U (r) . (5.7)
2m 2mr2 2mr2 sin2 θ
Las ecuaciones de Hamilton conducen a seis ecuaciones de movimiento,
, ṗr = p2θ /mr3 + p2φ /mr3 sin2 θ − dU/dr ,
ṙ = pr /m
θ̇ = pθ /mr2 , ṗθ = p2φ cos θ/mr2 sin3 θ , (5.8)
φ̇ = pφ /mr2 sin2 θ
, ṗφ = 0 .
5
Figura 5.2: Potencial efectivo para el problema de Kepler, tomando GM =
1, m = 1 y considerando L = 0.5, 1.0, 1.5, en colores azul, verde y rojo,
respectivamente.
6
Una propiedad importante de esta función en el caso de una partı́cula
puntual, es que su norma representa la probabilidad de hallar a la
partı́cula en una región dada,
Z b
Prob(a ≤ x ≤ b)(t) = |ψ(x, t)|2 dx . (5.12)
a
Oψ(~
e r , t) = ϕ(~
r, t) , Oψ(~
e p, t) = ϕ(~
p, t) .
x
e φ(~
r, t) = x φ(~r, t), pex φ(~
p, t) = px φ(~
p, t).
q
y de su “incertidumbre”, ∆O ≡ hO2 i − hOi2 .
[A, e ≡A
e B] eBe −B
eAe.
[x
ei , pej ] = ih̄ δij ⇒ [x
e, p
fx ] = [ye, p
fy ] = [ze, p
fz ] = ih̄ , (5.13)
7
Las descripciones en términos de coordenadas o momentos son alternativas y
complementarias: se dice que se sigue la representación de coordenadas, ψ =
ψ(~r, t), o la representación de momentos, ψ = ψ(~ p, t). El cuarto postulado
se satisface definiendo los operadores de posición y momento como
h̄2 2 ∂ψ(~r, t)
− ∇ ψ(~r, t) + U (~r, t)ψ(~r, t) = ih̄ . (5.17)
2m ∂t
Hϕ = E ϕ . (5.18)
8
estacionaria de la ecuación de Schrödinger. En general, hay una familia de
funciones ϕn que satisfacen la ecuación de Schrödinger,
Hϕn = En ϕn , (5.19)
cada una asociada a una energı́a particular En . Es posible tener varios esta-
dos estacionarios con una misma energı́a, en cuyo caso hablamos de estados
degenerados. En general, la ecuación de Schrödinger dependiente del tiempo
(5.16), tiene como solución
ak ϕk (~r)e−iEk t/h̄ ,
X
ψ(~r, t) = (5.20)
k
la energı́a del sistema está dada3 por < ψ|H|ψ >= |ak |2 Ek .
P
h̄2 d2 ϕ(x)
− = E ϕ(x) , (5.21)
2m dx2
de la forma ϕ00 +k 2 ϕ = 0, con k 2 = 2mE/h̄. Tiene dos soluciones linealmente
independientes, dadas por los dos posibles signos de k,
donde el factor (2π)−1/2 asegura que ϕ(x) tiene la misma norma que la
función A(k): la normalización de ϕ(x) requiere la de A(k).
3
empleando la notanción de bra (|ψ >) y kets (< ψ|).
9
En la representación de momentos la ecuación (5.21) tiene la forma,
" #
p2
− E ϕ(p) = 0 , (5.24)
2m
h̄2 d2 ϕ 1
− + mω 2 x2 ϕ = E ϕ . (5.27)
2m dx2 2
Dado que el Hamiltoniano del oscilador tiene la misma dependencia cuadráti-
ca en p que en x, las soluciones en ambas representaciones provienen de la
misma familia de funciones. Sabiendo que la Gaussiana tiene la propiedad de
que su transformada de Fourier es otra Gaussiana, se propone una solución
de la forma ϕ(x) = A exp(−α x2 /2): esta satisface la ecuación de Schrödin-
ger si α = mω/h̄, y E = h̄ω/2; su normalización implica A = (mω/πh̄)1/4 .
10
Yendo un paso más adelante, se puede proponer como solución general
una familia de funciones ϕn (x) = An (α1/2 x) exp(−α x2 /2), siendo An un po-
linomio de grado n, para mostrar que existe un conjunto infinito numerable
de soluciones al oscilador armónico,
1
Hϕn = En ϕn con En = n + h̄ω. (5.28)
2
Las funciones de onda son el producto de un polinomio por la Gaussiana,
1/4 r
1 mω mω
2
ϕn (x) = √ Hn x e−mωx /2h̄ , (5.29)
2n n! πh̄ h̄
11
Figura 5.3: Soluciones del oscilador armónico para n = 0, en azul; n = 1, en
verde; n = 2 en rojo; n = 3 en azul claro. Las gráficas tienen misma escala,
pero con un desplazamiento vertical para mostrarlas juntas.
12
y la ecuación de Schrödinger queda como,
( )
h̄2 1 ∂ 2 h̄2 1 ∂ ∂ϕ 1 ∂2ϕ
− (rϕ)− sin θ+ +U (r)ϕ = E ϕ .
2m r ∂r2 2mr2 sin θ ∂θ ∂θsin2 θ ∂φ2
(5.38)
La parte angular puede resolverse para el caso general mediante el método
de separación de variables, sin requerir conocer el potencial,
donde Y`m (θ, φ) son los armónicos esféricos, que satisfacen las ecuaciones de
valores propios
L2 Y`m (θ, φ) = `(` + 1)h̄2 Y`m (θ, φ) , Lz Y`m (θ, φ) = mh̄Y`m (θ, φ) , (5.40)
q q q
5
3 cos2 θ − 1 15
sin θ cos θ e±iφ 15
sin2 θ e±2iφ
`=2 16π ∓ 8π 32π
13
Dada la solución para la parte angular, al especificar el potencial U (r)
se tiene la ecuación radial para Rn` (r):
( )
h̄2 1 d dRn` `(` + 1)h̄2
− 2
r2 + + U (r) Rn` = E Rn` ; (5.42)
2m r dr dr 2mr2
Método variacional
El método variacional permite obtener aproximaciones a la energı́a del es-
tado base de sistemas para los cuales se desconoce la función de onda. Las
soluciones a la ecuación de Schrödinger para un potencial dado
Hϕj = Ej ϕj , j = 0, . . . (5.43)
son linealmente independientes, < ϕi |ϕj >= δij , formando una base de
manera que cualquier función arbitraria Φ puede escribirse como una com-
binación lineal de las funciones ϕj ,
∞
X
Φ= cj ϕj , (5.44)
j=0
|cj |2 = 1
P
con los coeficientes unı́vocamente determinados y cumpliendo
para Φ normalizada. El método variacional consiste en calcular
14
Método de perturbaciones independientes del tiempo
El método de perturbaciones es una herramienta aplicada históricamente
al cálculo de órbitas planetarias por Laplace en el siglo XIX. Este método
se usa para obtener soluciones aproximadas de la ecuación estacionaria de
Schrödinger, Hϕ = Eϕ, cuando el Hamiltoniano puede separarse en un
término soluble, H(0) , y un término perturbativo, H(1) , cuando los valores
propios de las soluciones correspondientes satisfacen E (1) E (0) . El método
provee términos correctivos para las energı́as y las funciones de onda.
Se considera un Hamiltoniano H = H(0) + H(1) , con funciones propias
(0) (1) (0) (1)
ϕk = ϕk + ϕk , y valores propios Ek = Ek + Ek , donde conocemos
(0) (0) (0) (0)
la familia de soluciones {ϕk ; Ek } de H(0) ϕk = E (0) ϕk . La ecuación
de Schrödinger se desarrolla en cuatro términos, donde podemos eliminar
(0) (0) (0)
H(0) ϕk = Ek ϕk , para obtener
(1) (0) (1) (0) (1) (1) (0) (1) (1)
H(0) ϕk + H(1) ϕk + H(1) ϕk = Ek ϕk + Ek ϕk + Ek ϕk . (5.46)
(0)
Las funciones {ϕk } forman una base que permite desarrollar las funciones
(1)
perturbativas {ϕk }, E X E
(1) (0)
ϕk = ck` ϕ` , (5.47)
`
15
5.3. Estructura de átomos y moléculas
5.3.1. Átomos hidrogenoides
El estudio de los átomos hidrogenoides es una de las aplicaciones más im-
portantes del potencial central, dado por la atracción Coulombiana entre un
núcleo de carga Ze y un electrón de carga −e:
Ze2
U (r) = − . (5.51)
r
La solución a la ecuación (5.42) es una familia infinita y numerable de fun-
ciones Rn` (r) que corresponden a una energı́a
! ! !
me4 Z Ze2 1 mc2 1
En = − = − = − α2 Z , (5.52)
2h̄2 n2 2a0 n2 2 n2
con a0 = h̄2 /me2 ' 5.292 × 10−9 cm el radio de Bohr, α = e2 /h̄c ' 1/137 la
constante de estructura fina, y la masa m es la masa reducida del sistema,
m = me mn /(me + mn ), para un núcleo de masa mn . Los números (n, `)
tienen valores enteros positivos. En este caso, la energı́a es independiente de
`, por lo que existen varias funciones de onda para una energı́a dada, es decir
una degeneración de estados. La función de onda completa es de la forma,
con degeneración gn = 2n2 , determinada por los valores que pueden adquirir
los números cuánticos de momento angular, {`, m}, donde |m| ≤ ` < n, y los
dos valores posibles del número cuántico de espı́n, ms = ±1/2. Las funciones
Rn` (r) se relacionan con los polinomios de Laguerre
s 3 `
2 (n − ` − 1)! 2r
Rn` (r) = − e−r/naz Ln+`
2`+1
(2r/naz ) , (5.54)
naz [(n + `)!]3 naz
donde az = a0 /Z, y Lqp (x) son los polinomios generalizados de Laguerre, que
se calculan a través de los polinomios ordinarios de Laguerre, Lp (r),
dq Lp dp
Lqp (x) = Lp (x) = ex xp e−x ,
(x) , (5.55)
dxq dx p
o directamente de la expresión,
p! dp−q
Lqp (x) = (−1)q ex x−q p−q e−x xp .
(5.56)
(p − q)! dx
16
E=0 n→∞ La parte radial de las funciones de onda
de los primeros niveles de sistemas hidro-
E = −0.85 eV n=4 genoides son,
E = −1.51 eV n=3
R1,0 = 2 a−3/2
z e−r/az ,
r
R2,0 = 2 (2az )−3/2 e−r/2az ,
1−
E = −3.40 eV n=2 2az
1 r
R2,1 = 2 (2az )−3/2 √ e−r/2az .
3 2az
De interés particular para la espectros-
copı́a es la estructura de los niveles de
energı́a de los átomos hidrogenoides. Pa-
ra un electrón alrededor de un núcleo de
masa infinita, las dimensiones están es-
tablecidas por la constante de Rydberg,
1 me e4
R∞ = ' 10 973 731.568 549 m−1 ,
hc 2h̄2
la escala de energı́as,
y de longitudes de onda,
E = −13.60 eV n=1
1
' 911.267 050 550Å .
Niveles de energı́a R∞
átomo hidrogenoide.
17
una energı́a E1 = −R∞ hc/(1 + mp /me ) ' −13.598 296 eV, y puede
ser ionizado por fotones de longitud de onda λ1→∞ ' 911.763 342 Å.
Las principales series espectroscópicas del hidrógeno están dadas por
transiciones a distintos niveles:
18
El FeXXVI es un núcleo de relevancia en rayos X, siendo que la serie
de Lyman cae en ∆E21 ' 6.38 keV. Nótese que para núcleos altamente
ionizados, con un deficit A de electrones, el estado base está a una
energı́a E1 ≈ 13.6 A2 eV.
19
Cuadro 5.2: Series espectroscópicas de algunos sistemas hidrogenoides
20
La corección relativista al momento del electrón es de orden p/mc ∼
αZ 1/2 ' Z 1/2 /137, pequeña para valores comunes de Z; puede descri-
birse mediante una expansión en términos de orden p/mc,
p2 3p4
Ecin = (γ − 1)mc2 ≈ − + ... , (5.58)
2m 8m3 c2
siendo la primera corrección igual a −3p4 /8m3 c2 . El cálculo de per-
turbaciones da
* +
3p4 E1z 2n 3
∆Erel = n`m − 3 2 n`m = α2 2 − . (5.59)
8m c 2n ` + 1/2 2
E1z 2n 3
2
∆Erel + ∆Eeo =α − . (5.62)
2n2 j + 1/2 2
Existe una degeneración en j con el número cuántico ms : por ejemplo,
el nivel 2S1/2 (n = 2, ` = 0, ms = +1/2 ⇒ j = 1/2) tiene misma
energı́a que el nivel 2P1/2 (n = 2, ` = 1, ms = −1/2 ⇒ j = 1/2).
21
El desplazamiento de Lamb es la corección radiativa debida a la
cuantización del campo electromagnético, la cual considera a este como
un conjunto de fotones virtuales que viven un tiempo ∆t < E/h̄,
creando pares e± que originan una distribución espacial de “carga
virtual”, denominada polarización del vacı́o. La corrección es,
E1z 1 − δ`=0
∆Elamb = α3 K(n, `) ± , (5.63)
2n3 π(j + 1/2)(` + 1/2)
Estructura hiperfina
~p ,
El protón tiene un momento magnético asociado a su espı́n S
~p
eS
µ
~ p = γp , (5.64)
mp c
El efecto Zeeman
El momento magnético del átomo, tanto la parte debida al momento an-
~ como la originada en el espı́n, interacciona con un campo
gular orbital, L,
magnético externo B~ resultando en un término adicional al Hamiltoniano
Hz = (eB/2mc) · J, donde J~ = L
~ ~ ~ + S.
~
22
Empleando el método de perturbaciones se encuentra que los niveles de
energı́a se desdoblan de forma proporcional al campo magnético,
j(j + 1) + s(s + 1) − `(` + 1) eh̄B
∆Ez = mj 1 + . (5.66)
2j(j + 1) 2mc
incluyendo el espı́n de cada electrón mediante χms,j . Ası́, el estado del átomo
queda especificado con 4Z números cuánticos,
23
El problema fundamental en la descripción de estos sistemas radica en el
término de repulsiones entre electrones que contiene una doble sumatoria y
no puede ser despreciado5 ,
X X X
e2 /rij / Ze2 /ri ' Z(Z − 1)/2Z 2 → 1/2 para Z 1 .
(5.69)
El método de perturbaciones resulta inadecuado para el estudio de átomos
con varios electrones, incluso en el caso del helio, y sólo puede aportar una
descripción cualitativa del sistema.
Un método más apropiado para tratar este problema consiste en apro-
ximar el potencial mediante “algún” potencial central, Vc (r), y un término
residual que pueda ser despreciado, o tratado como una perturbación,
Z
" 2 #
X p j
H= + Vc (rj ) + W . (5.70)
j
2m
24
Figura 5.4: Un potencial que cumple
con los lı́mites asintóticos, indicados
por las lı́neas punteadas: Vc → −1/r
para r 1,Vc → −Z/r para r 1.
En` > En0 ` para n > n0 , En` > En`0 para ` > `0 . (5.73)
25
En esta expresión se puede ver que dos electrones no pueden ocupar el
mismo estado, ya que si a = b entonces ϕ = 0, de acuerdo con el principio
de exclusión de Pauli, inherente a los fermiones.
La configuración de un átomo se define especificando los valores {n, `}
de los orbitales que ocupan sus electrones: n`x indica que x electrones están
ocupando el estado definido por {n, `}, donde x debe ser menor o igual
a 2(2` + 1), la degeneración del estado de momento angular `, incluyendo
el espı́n. Las configuraciones de estados base de los primeros elementos se
muestran en la tabla 5.3 y en la figura 5.5, transcrita de Cohen-Tannoudji.
El llenado de las capas superiores no se da en secuencia con n a partir de
niveles con n ≥ 3, de forma que generalmente el nivel 4s se ocupa antes que
el 3d. Otro punto relevante es el potencial de primera ionización, mostrados
también en la tabla 5.3: un potencial de ionización bajo se relaciona con la
presencia de lı́neas espectrales en el visible; por ejemplo en Na y Mg.
26
Elemento Z χ (eV) n ` Configuración Capa externa
H 1 13.6 1 0 1s K
He 2 24.6 1 0 1s2 K
Núcleo de He: 2 electrones
Li 3 5.39 2 0 2s L
Be 4 9.32 2 0 2s2 L
B 5 8.30 2 1 2s2 2p L
C 6 11.3 2 1 2s2 2p2 L
N 7 14.5 2 1 2s2 2p3 L
O 8 13.6 2 1 2s2 2p4 L
F 9 17.4 2 1 2s2 2p5 L
Ne 10 21.6 2 1 2s2 2p6 L
Núcleo de Ne: 10 electrones
Na 11 5.14 3 0 3s M
Mg 12 7.64 3 0 3s2 M
Al 13 5.98 3 1 3s2 3p M
Si 14 8.15 3 1 3s2 3p2 M
P 15 10.6 3 1 3s2 3p3 M
S 16 10.4 3 1 3s2 3p4 M
Cl 17 13.0 3 1 3s2 3p5 M
Ar 18 15.8 3 1 3s2 3p6 M
Núcleo de Ar: 18 electrones
K 19 4.34 4 0 4s N
h̄2 2 Z e2
A ZB e2 ZA e2 ZB e2 ZA ZB e2 e2
H12 = − ∇1 + ∇22 − − − − + + .
2me r1A r1B r2A r2B R r12
Aproximación de Born-Oppenheimer
La aproximación de Born-Oppenheimer separa la función de onda en una
parte electrónica y una nuclear,
27
Figura 5.5: Esquema de niveles energéticos de distintos átomos de acuerdo
a un modelo de potencial central. Entre paréntesis se indica la degeneración
de cada nivel, g` = 2(2` + 1). La posición relativa de los niveles superiores
puede variar entre diferentes átomos. Basada en §XIV, Cohen-Tannoudji.
Hel ϕel = Eel (R)ϕel , HAB ϕAB + Eel (R)ϕAB = E ϕAB . (5.78)
28
Movimiento translacional e interno
La componente nuclear del Hamiltoniano puede separarse en dos términos
2 /2M + P 2 /2µ, permitiendo describir el movimiento de
cinéticos, HAB = Pcm R
A y B en términos de la translación del centro de masa,
~ cm = MA ~xA + MB ~xB ,
X M = M A + MB ; (5.79)
MA + MB
y el movimiento relativo entre las dos moléculas, descrito por la coordenada
relativa y la masa reducida,
~ = ~xA − ~xB , 1 1 1
R = + . (5.80)
µ MA M B
~ cm )ϕint (R),
La función de onda de los núcleos se separa, ϕAB = ϕtrans (X ~ con
la parte translacional obedeciendo una ecuación de partı́cula libre, mientras
que la interna, fı́sicamente más relevante, sigue un potencial central:
h̄2 2
− ∇ ϕint + Eel (R)ϕint = Eint ϕint . (5.81)
2µ R
Estados rotacionales
La simetrı́a esférica en la ecuación (5.81) nos permite separar la función de
onda en una parte radial y una de armónicos esféricos,
~ = 1 Zint (R)Ym (θ, φ) ,
ϕint (R) (5.82)
R
donde L2 Ym = ( + 1)h̄2 Ym , Lz Ym = mh̄ Ym , y una configuración niveles
rotacionales asociada,
h̄2
Erot = ( + 1) . (5.83)
2µR2
Estados vibracionales
La parte radial de la ecuación de estructura interna,
!
h̄2 d2 Z ( + 1)h̄2
− + Eel (R)Z = Eint − Z, (5.84)
2µ dR2 2µR2
puede resolverse cuando consideramos la estructura de la molécula en un
estado de equilibrio electrónico, alrededor de un mı́nimo de Eel (R). La ex-
pansión en serie omite el término lineal,
1
Eel (R) = Eel (R0 ) + µω02 (R − R0 )2 + . . . , (5.85)
2
29
dando lugar a una ecuación de oscilador armónico
h̄2 d2 Z(x) 1 2 2
− + µω0 x Z(x) = Evib Z(x) , (5.86)
2µ dx2 2
donde x = R − R0 , ω02 = (1/µ) (d2 Eel /dR2 )0 , Evib = Eint − Eel (R0 ) − Erot ;
por lo que tenemos:
1
vib
Ev = v + h̄ω0 . (5.87)
2
Recapitulando, la energı́a de la molécula se separa en una parte trans-
lacional y otra interna, E = Etrans + Eint , con la parte interna incluyendo
componentes electrónica, vibracional y rotacional: Eint = Eel + Evib + Erot .
Al considerar un punto de equilibrio de la estructura electrónica, que define
una distancia internuclear R0 , la estructura vibracional se aproxima a la de
un oscilador armónico. La estructura rotacional está dada por los armónicos
esféricos, siendo los niveles de energı́a de la forma Erot = ( + 1)(h̄2 /2µR0 ).
Las estructuras vibracional y rotacional de la molécula dependen entonces de
los parámetros R0 y d2 Eel /dR02 , determinados por la estructura electrónica.
donde u(r) = π −1/2 e−r es la función de onda del estado base del hidrógeno,
{n = 1, ` = 0, m = 0}, y c1 , c2 son coeficientes a determinar. Existen √ dos
soluciones: √
una simétrica y una antisimétrica, correspondientes a c1 = 1/ 2,
y c2 = ±1/ 2, que dan lugar a dos formas para Eel (R),
! !
e2 e2 α±β
± ±
Eel (R) =− + Emol = −1 + 2 , (5.89)
a0 a0 1±∆
siendo α, β y ∆ funciones de R:
!
1 ∆ R2
−2R
α=e 1+ , β= − e−R (1 + R) , ∆=e −R
1+R+ ,
R R 3
30
Figura 5.6: Energı́a de enlace en un modelo de molécula H2 . La curva su-
perior, en verde, corresponde a E − , positiva para toda R y, por tanto, a
un sistema no ligado. La solución E + corresponde al enlace molecular con
equilibrio en R ' 2.49 a0
±
siendo −e2 /a0 la energı́a asociada a dos átomos de hidrógeno y Emol repre-
senta la energı́a de enlace, mostrada en la figura 5.6.
31
5.4.1. Método de perturbaciones dependientes del tiempo
El formalismo de perturbaciones para resolver la ecuación estacionaria de
Schrödinger es particularmente útil para el estudio de la estructura de sis-
temas cuánticos, pero no aporta información acerca del comportamiento
temporal de estos, por ejemplo la interacción con un campo de radiación.
Un desarrollo perturbativo dependiente del tiempo aborda el problema
de Hamiltoniano de la forma H = H(0) + H(1) , donde el primer término es
independiente del tiempo y sus estados estacionarios, {ϕk , Ek }, soluciones
de H(0) ϕk = Ek ϕk , son conocidos. Estos representan la estructura básica
del sistema cuántico y su conjunto constituye una base, de forma que es
posible desarrollar cualquier función de onda en términos de la familia de
estados estacionarios,
∞
ck (t) ϕk (~x)e−iEk t/h̄ ,
X
ψ(~x, t) = (5.90)
k=0
32
proporcional a |ċf |2 . La probabilidad promedio por unidad de tiempo de que
ocurra la transición i → f es:
Z 2
1 1 T
(1) iωf i t0 0
χf i = hf | H |ii e dt , (5.94)
T h̄ 0
2
1 q~
H= p
~ − A(~
x , t) + qφ(~
x, t) . (5.95)
2m c
e ~ ~ e2 ~ ~
Hint = p~ · A + A · p~ + A·A , (5.97)
2mc 2mc2
representa la interacción
del sistema con el campo electromagnético. El
~ ~
término p~ · A + A · p~ representa la interacción del sistema con un fotón,
mientras que el término (A ~ · A)
~ representa la emisión o absorción simulta-
nea de dos fotones. Al igual que en el caso clásico, tenemos la libertad de
fijar la norma en que trabajamos, siendo más conveniente aquı́ la norma de
Coulomb, ∇ · A ~ = 0, equivalente a p~ · A
~ = 0.
33
En entornos astrofı́sicos comunes el campo electromagnético incidente es
muy débil comparado al potencial Coulombiano, |E| ~ e/a2 , y el fotón emi-
0
tido en el proceso tiene una energı́a comparable al potencial de ionización,
e2 /a0 ∼ h̄ω, por lo que podemos tratar al campo externo como una pertur-
bación dependiente del tiempo, y despreciar el término de dos fotones. La
importancia relativa de los dos términos radiativos en Hint y el estructural
Hest puede evaluarse usando A ∼ B/k ∼ (c/ω)E, y p ∼ αmc,
|(e/2mc) A ~ · p~| |(e2 /2mc2 )(A~ · A)|
~ e A
e2 /a0
∼ ∼ ∼ 1,
|p2 /2m| |(e/2mc) A ~ · p~ | c p h̄ω
Para calcular el efecto del término dominante en Hint expresamos el
~ x, t) = A(t)ei~k·~x êα , donde êα denota el vector
potencial vectorial como A(~
de polarización para calcular la integral en la ecuación (5.94),
i T
e
Z
~ · p~ i eiωf i t0 dt0 ,
ċf (ω) ≡ − f A
h̄ 0 2mc
2π e
D
~
E
= −i A(ωf i ) f eik·~x êα · p~ i , (5.98)
h̄ 2mc
donde A(ω) es la transformada de Fourier de A(t). El cálculoD de las probabili- E
~
dades de transición se remite al de los elementos de matriz, f | eik·~x êα · p~ |i .
34
5.4.4. Coeficientes de Einstein
Los coeficientes de Einstein se relacionan directamente con las tasas de tran-
sición. Para absorción o emisión estimulada tenemos,
hχf i i = Bf i Jν , (5.100)
~
ω 2 |A(ω)|2
Jν = . (5.101)
2 cT
Podemos expresar los coeficientes de Einstein en términos de la matriz del
dipolo,
8π 2 2 4ωf3i 2
Bf i = 2 d~f i = Bif , Af i = d~f i . (5.102)
3h̄ c 3h̄c3
(Rybicki §10.3).
35
Transiciones cuadrupolares eléctricas
Estas son relevantes cuando las transiciones dipolares no son posibles. Se
deben evaluar los elementos de matriz del tensor de cuadrupolo eléctrico,
hf |Qαβ | ii, donde Qαβ = −e(3xα xβ − r2 δαβ ). Las reglas de selección co-
rrespondientes nos dicen que las transiciones cuadrupolares son permitidas
entre estados con ∆` = 0 y ∆` = ±2, con ∆m = 0, ±1, ±2.
2
donde Ef = D h̄ωf i − Ze /2a
E z es la energı́a cinética del electrón liberado. Del
i~k·~
x
cálculo de f |e êα · p~|i se obtiene la sección eficaz de fotoionización,
1 8π Z 4 me e10
σbf (ω) ' √ gbf (ω) , (5.105)
n5 3 3 ch̄3 (h̄ω)3
36
911Å hasta h̄ω <
∼ 1 keV, abarcando la mayor parte del ultravioleta, hasta
los rayos X blandos. Por encima de estas energı́as se tiene un decaimien-
to ∝ ω −3 que permite estimar la densidad de columna en la dirección de
observación.
(Ryb 10; Shu 21, 22, 23, 24)
Figura 5.7: Sección eficaz para transiciones entre estados ligados y libres -
expresión (5.105). Los picos en σ(ω) corresponden al acceso a las distintas
capas de la configuración electrónica, K, L, M, en un átomo hipotético.
37
Bibliografı́a
[1] Cohen-Tannoudji...
[2] “Classical mechanics”, Goldstein, H., 2nd edition, ed. Addison Wesley.
[4] Shu...
38
Radiación: Tarea 1
Martes 27-sep-2016 → lunes 10-oct-2016
1
c) que la expresión anterior, junto con la condición de que el flujo en
τ = 0 solo tiene una componente emergente (i.e. parte “incidente”
es cero), implican
3 4 2
4
T = Te τ + .
4 3
4. Dos fermiones con tres posibles estados individuales (2 puntos)
Para un sistema de dos fermiones con tres posibles estados individuales,
verifique que:
xn j
X Y Y X nj
j ≡ x j ,
{nj } j j nj
2
Radiación: solución a la tarea 1
Jueves 6-oct-2016
1
de dónde I1 + I2 + I3 = 4πq.
~ = 2qa 3x~r − x̂ .
E
r5 r3
Solución:
Se hace la expansión a primer orden:
!−3/2
2ax a2 3ax
−3/2
−3
2
r ± 2ax + a 2
=r 1± 2 + 2 ' r−3 1 ∓ ,
r r r2
de dónde
~r 3ax 3ax ax̂ 3ax 3ax
~ '
E +q 3 1+ 2 − 1− 2 −q 3 1+ 2 + 1− 2
r r r r r r
~r 6ax ax̂ 3x~r x̂
' +q 3 − 2q 3 = 2qa − 3 . (2)
r r2 r r5 r
3. Atmósfera gris...
2
La presión de radiación:
Z +1 Z +1 Z +1
1 2π 2π 4π a0 (τ )
P = I 2πµ2 dµ = a0 µ2 dµ+ a1 µ3 dµ = ,
c −1 c −1 c −1 c 3
de donde P = u/3.
b) Para estimar la función fuente calculamos la intensidad media,
Z +1
1
J= I dµ = a0 ,
2 −1
3 2
T (τ ) = Te4 τ +
4
.
4 3
3
0, 1; y los macro-estados, {nj } = {000, 100, 010, 001, 110, 101, 011, 111} .
De donde, el término de la izquierda:
Y nj
xna a xnb b xnc c
X X
x =
j
{nj } j {nj }
= 1 + xa + xb + xc + xa xb + xa xc + xb xc + xa xb xc .
El término de la derecha:
Y X nj Y
xj = [1 + xj ] = (1 + xa )(1 + xb )(1 + xc )
j nj j
= 1 + xa + xb + xc + xa xb + xa xc + xb xc + xa xb xc .
4
Radiación: Tarea 2
Jueves 17-nov-2016 → lunes 28-nov-2016
1
objetos se movı́a más rápido que la velocidad de la luz1 . Es posible de-
mostrar que este es un efecto aparente.
Supóngase que un objeto emite una masa luminosa con una velocidad
~v con ángulo θ con respecto a la lı́nea de visión, como se ilustra en la
figura 1.
a) Mostrar que la velocidad trasnversa aparente de la masa luminosa
inferida por el observador (inferida a partir de la distancia y la
velocidad angular) es
v sin θ
vapp = .
1 − (v/c) cos θ
b) Mostrar que vapp puede ser mayor que c. Encontrar el ángulo para
el cual vapp es máximo y mostrar que el valor máximo es vapp = γv.
3. Colisión relativista
Considerese la colisión de dos partı́culas relativistas, a y b, cuyos cua-
drivectores de energı́a momento son2
! !
Ea0 Eb0
pαa0 = , pαb0 = ,
p~a0 p~b0
con pα0 = pαa0 + pαb0 .
1
¡estrı́ctamente prohibido!
2
Por comodidad usamos unidades con h̄ = c = me = 1.
2
a) Mostrar que la transformación de Lorentz que nos lleva al marco
de referencia de momento nulo está dada por β~ = p~0 /E0 .
b) Aplicando esta expresión al caso general del efecto Compton, en
el cual un fotón de cuadri-momento k0 interactúa con un electrón
con p0 , dados por,
! !
ω0 γ0
k0 = , p0 = ,
ω0 k̂0 γ0 β~0
3
Radiación: Solución a la tarea 2
Solución:
La combinación de 75 % de N2 , de peso molecular 2×14, y 25 % de
O2 , de peso molecular 2 × 16, tiene un peso molecular µ = 29. Al
reemplazar la ecuación de estado en la de equilibrio hidrostático,
1 dρ 1 dT µmh g
+ =− . (1)
ρ dz T dz kT
con H0 = kT0 /µmh g ' 8.57049×105 cm, o 8570m, para los datos
especificados y los valores de las constantes mh , g, k.
Solución:
La ecuación de conservación de energı́a, en el caso adiabático, es
dU = T dS − P dV = −P dV . Un gas de moléculas diatómicas
1
tiene una energı́a interna dada por U = (5/2)N kT , siendo N el
número total de partı́culas. De donde
5 5 dV
d(N kT ) + P dV = 0 ⇒ N k dT + N kT = 0,
2 2 V
es decir
5 dT dV dn dρ
=− = = ,
2 T V n ρ
siendo n la densidad numérica, n = N/V = ρ/µmh , donde toma-
mos en cuenta que n y V son variables con z. Reemplazando en
la ec. (1), tenemos
7 1 dT µmh g dT 2 µmh g
=− ⇒ =− = constante ;
2 T dz kT dz 7 k
con los datos obtenemos
dT γ−1 µmh g
=− .
dz γ k
Solución:
Si T (z) = T0 − θz, con dT /dz = −θ constante, tenemos T (z) = 0
en z = T0 /θ ' 4.51 × 106 cm = 45.1 km, para T0 = 20◦ C y el gra-
diente del modelo estándar. El modelo es inaplicable para valores
mayores de z.
Solución:
2
El modelo adiabático falla a 30 km y el estándar a 45 km. Mues-
tro gráficas hasta 20 km. Para expresar las funciones soluciones,
podemos escribir
1 mg
θ= ,
η k
con η = γ/(γ − 1) = 1 + m/2 = 3.5, para atmósfera adiabática de
partı́culas con m grados de libertad; η ' 4.5 para caso estándar
(⇒ γ ' 1.29); η → ∞ caso isotérmico (⇒ γ → 1). Usando este
coeficiente, la solución general para una atmósfera de gradiente
dT /dz constante queda como
T (z)
= T0 (1 − z/ηH0 ) ,
ρ(z) = ρ0 (1 − z/ηH0 )η−1 ,
P (z) = P (1 − z/ηH )η ,
0 0
2. Movimiento superlumı́nico
Al poco tiempo del descubrimiento de los cuasares se hicieron obser-
vaciones que indicaban movimientos más rápidos que la luz, prohibido
por las leyes de la fı́sica. Este es un efecto aparente: supóngase una
masa luminosa con una velocidad ~v con ángulo θ con respecto a la
lı́nea de visión, como se ilustra en la figura 1.
3
a) Mostrar que la velocidad transversa aparente de la masa luminosa
inferida por el observador, es
v sin θ
vapp = .
1 − (v/c) cos θ
Solución:
Consideramos un primer fotón emitido en x1 = d, y1 = 0, que
llega al observador en t1 = x1 /c; un segundo fotón es emitido en
un tiempo posterior t, desde x2 = d − v cos θt, y2 = v sin θt, lle-
gando al observador en t2 = t + x2 /c, de manera que la velocidad
aparente según el observador será
∆y y2 − y1 v sin θt
vapp = = = ,
∆t t + x2 /c − x1 /c t + (d − v cos θt − d)/c
v sin θ
= .
1 − (v/c) cos θ
b) Mostrar que vapp puede ser mayor que c. Encontrar el ángulo para
vapp máxima y mostrar que el máximo es vapp = γv.
Solución:
vapp > c ⇒ (v/c) sin θ > 1−(v/c) cos θ ⇒ (v/c) > (sin θ + cos θ)−1 ,
que puede cumplirse para sin θ + cos θ > 1, que a su vez se cum-
ple para 0 < θ < π/2. El ángulo para el cual vapp es máxima se
determina derivando con θ,
dvapp v cos θ v sin θ(v/c) sin θ
= −
dθ 1 − (v/c) cos θ (1 − (v/c) cos θ)2
v cos θ(1 − (v/c) cos θ) − v sin θ(v/c) sin θ
=
(1 − (v/c) cos θ)2
v cos θ − v 2 /c
= = 0,
(1 − (v/c) cos θ)2
4
3. Colisión relativista
Considerese la colisión de dos partı́culas relativistas, a y b, cuyos cua-
drivectores de energı́a momento son1
! !
Ea0 Eb0
pαa0 = , pαb0 = ,
p~a0 p~b0
con pα0 = pαa0 + pαb0 .
a) Mostrar que la transformación de Lorentz que nos lleva al marco
de referencia de momento nulo está dada por β~ = p~0 /E0 .
Solución
La transformación de Lorentz del momento es,
p~00 = −γ βE
~ 0 /c + p~0 + (γ − 1)(β̂ · p~0 )β̂ ,
Solución
E0 = ω0 + γ0 , p~0 = ω0 k̂0 + γ0 β~0 , de donde
ω0 k̂0 + γ0 β~0
β~ = p~0 /E0 = .
ω0 + γ0
Siendo k̂0 unitario y |β~0 | < 1, se asegura |β|
~ < 1.
1
Unidades con h̄ = c = me = 1.
5
q
c) Mostrar que en el centro de momento se cumple γ00 = 1 + ω002 .
Solución
En el centro de momento p~00 = ω00 k̂00 + γ00 β~00 = 0, de donde
q
0 0 0 0
ω00 k̂00 = −γ00 β~00 ⇒ ω02 = γ02 β02 = γ02 − 1 ⇒ γ00 = 1 + ω002 .
Solución: q
Siendo el marco de momento tenemos también γ10 = 1 + ω102 .
Por lo que la conservación de energı́a queda como
q q
ω00 + 1 + ω002 = ω10 + 1 + ω102 ,
que se cumple solamente si ω00 = ω10 , y por tanto γ00 = γ10 . No-
tamos que no se requiere k̂00 = k̂10 , por lo que el electrón y fotón
pueden modificar la dirección de sus vectores de momento. En el
centro de momento la interacción consiste en la rotación conjunta
de los vectores k̂ 0 y β~ 0 .
z=np.linspace(0.0,25.0,1000)
t0=293.15; h0=8.570; p0 = 1013.25; rho_0 = 1.206e-3
#modelo isotermico
6
t_iso = np.linspace(t0,t0,1000)
rho_iso = rho_0 * np.exp(-z/h0)
p_iso = p0 * np.exp(-z/h0)
#modelo adiabatico
eta=3.5
f = (1.0-z/(eta*h0))
t_ad = t0 * f
rho_ad = rho_0 * f**(eta-1.0)
p_ad = p0 * f**eta
#modelo estandar
eta=4.5
f = (1.0-z/(eta*h0))
t_st = t0 * f
rho_st = rho_0 * f**(eta-1.0)
p_st = p0 * f**eta
plt.ion()
plt.axis([0.0,20.0,-200.0,30.0],fontsize=16)
plt.grid()
plt.xlabel("z (km)",fontsize=16)
plt.ylabel("T(C)",fontsize=16)
plt.plot(z,t_iso-273.15); plt.plot(z,t_ad-273.15); plt.plot(z,t_st-273.15)
plt.figure()
plt.axis([0.0,20.0,0.0,1050.0],fontsize=16)
plt.grid()
plt.xlabel("z (km)",fontsize=16)
plt.ylabel("P(mbar)",fontsize=16)
plt.plot(z,p_iso); plt.plot(z,p_ad); plt.plot(z,p_st)
plt.figure()
plt.axis([0.0,20.0,0.0,1.5e-3],fontsize=16)
plt.grid()
plt.xlabel("z (km)",fontsize=16)
plt.ylabel(r"$\rho$(g cm$^{-3}$)",fontsize=16)
plt.plot(z,rho_iso); plt.plot(z,rho_ad); plt.plot(z,rho_st)
7
Figura 2: Perfiles de presión y densidad. Azul corresponde a modelo isotérmi-
co, verde adiabático y rojo estándar. 8
Figura 3: Perfiles de temperatura. Azul corresponde a modelo isotérmico,
verde adiabático y rojo estándar.