Apuntes de La Práctica de Astrofisica

Apuntes de la práctica de Astrofísica
6 de diciembre de 2018
Willy!
ma
ble
Los resueltos de
pro
No hay ALF
Sobre estos apuntes Estos apuntes/resueltos que usted está viendo fueron creados por un alumno mientras cursaba la
materia. Es por ello que podrían haber errores de tipeo, errores conceptuales, de interpretación en los resultados, etc. Use
estos apuntes con precaución. Estos apuntes no son oficiales de ninguna cátedra. Lea atentamente el prospecto. En caso de
notar algún efecto adverso suspenda inmediatamente su uso y consulte con su profesor de cabecera.
El alumno autor de estos apuntes cursó la materia durante el segundo cuatrimestre de 2018, este link conduce a la
página oficial del curso.
Encontrá más resueltos de Alf en este link.
Box 1 - ¿Cómo se hacen estos apuntes?
Estos apuntes están hechos usando un programa llamado Lyxa . Para hacer los dibujos se usó Inkscape y después se
insertó las imágenes en formato svgb directamente en Lyx.
En este repositorio de GitHub se encuentra la plantilla (template) que Alf usa actualmente, con todo lo necesario
para compilarla y empezar a divertirse.
a Lyx es una interfaz gráfica para Latex que hace que la escritura se vuelva extremadamente fluida y veloz (al punto de poderse
tomar apuntes en vivo durante una clase).

b svg es el formato nativo de Inkscape.
1
ÍNDICE ÍNDICE
Índice
1. Guía 1 - El universo 5
1.1. Cúmulos estelares . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
2. Guía 2 - Fluidos clásicos 8
3. Guía 3 - Relatividad 13
4. Guía 4 - Atmósferas estelares 16
5. Guía 5 - Medios interestelares 23

5.1. Rankine-Hugoniot . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
6. Guía 6 - Estructura y evolución estelar 28

6.1. Resumen . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 29
7. Guía 7 - Campos magnéticos 32
8. Guía 8 - Galaxias 36
9. Guía 9 - Cosmología 41
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

2 Si ves un error, avisale a Alf , 201812062056 CC0
Índice alfabético
Guía 1 problema 16, 7
3
REFERENCIAS REFERENCIAS
Referencias
[Carroll and Ostlie, 2017] Carroll, B. W. and Ostlie, D. A. (2017). An introduction to modern astrophysics. Cambridge
University Press.
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

1 GUÍA 1 - EL UNIVERSO
1. Guía 1 - El universo
La guía 1 la podemos dividir en cuatro partes:
1. Órdenes de magnitud.
2. Determinación de no sé qué.
3. Radiación de cuerpo negro.
4. Leyes de Kepler.
Guía 1, ejercicio 2
Fotometría La fotometría Consiste en medir la intensidad de la radiación que proviene de algún astro en una determinada
banda. Por ejemplo tenemos la banda del ultravioleta, la del azul, etc. En particular las bandas y sus anchos son
Sistema fotométrico de Johnson-Morgan (∼1950)
Banda Longitud de onda λ Ancho de banda ∆λ
U Ultra violeta 365 nm 66 nm
V Visible 580 nm 88 nm
B Azul 440 nm 94 nm
La elección de estas bandas es histórica. Hoy se las sigue utilizando por cuestiones históricas, pero uno podría mirar
cualquier banda.
Espectrometría Se determina la luminosidad de un astro como función de la longitud de onda. Es como una fotometría
pero con más detalle. En la espectrometría se pueden observar muchas cosas: la composición de una estrella y/o del medio
en el que viajó la luz a partir de las líneas de absorsión, también se puede estudiar el corrimiento doppler lo cual tiene
información sobre el movimiento de la estrella, etc.
Polarimetría Se estudia la polarización de la radiación recibida lo cual brinda información acerca del medio en el que ésta
se propagó. En particular puede indicar la presencia de campos magnéticos en el medio interestelar.
Fácil: E = hν por lo tanto λ = hcE . La longitud de onda es entonces λ ∼ 21 cm. Esta línea espectral fue la primera que
se utilizó para mapear la distribución de hidrógeno en el universo.
Ítem a
En cuanto a la órbita circular de la tierra la podemos justificar, aproximadamente, con el hecho de que siempre vemos
al sol del mismo tamaño con lo cual la distancia tiene que ser más o menos constante. En cuanto a la redondez de la órbita
de Venus, no podemos decir nada. Simplemente es cómodo. Y en cuanto a la coplanaridad1 tampoco. Para poder justificar
estas otras dos hipótesis tendríamos que recurrir a observaciones astronómicas. La resolución del problema es trivial usando
geometría. En particular se obtiene que si d y D son las distancias mínima y máxima respectivamente, se satisface que
d+D
AU =
2
Ítem b Vale 1,497 × 108 km.
Ítem c El dibujito es así:

Observador B
Observador A
Tierra Venus
1 Sorprendente, la palabra “coplanaridad” existe. Pensé que la estaba inventando, ja.
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

1 GUÍA 1 - EL UNIVERSO
donde ` = 10 × 103 km y δ = 49,800 . Ahora sólo hay que hacer trigonometría para encontrar la distancia mínima
d ≈ 4,14 × 107 km
Usando el dato del “tamaño angular” cuando está más lejos se termina encontrando que
AU ∼ 1,491 × 108 km
Vamos a suponer que los objetos emiten radiación de cuerpo negro. Esto no es cierto, pero mal que mal le pega para
tener una idea:
← Cuerpo negro ideal

← Radiación afuera de la atmósfera
← Radiación a nivel del suelo
La densidad espectral de energía en el cuerpo negro (que obtuvo Planck) es
2hν 3 1
Bν (T ) = hν
c2 e kT −1
Como estamos interesados en la longitud de onda vamos a hacer el cambio de variable tal que
dν
Bν dν = −Bλ dλ ⇒ Bλ = −Bν
dλ
c dν
Usando que ν = λ encontramos que dλ = − λc2 y entonces
c 2h c 3 1
Bλ = hc
λ2 c2 λ e λkT −1
Ahora sólo hay que encontrar para qué λ se da el máximo de esta expresión y despejar en función de T . No es más que
derivar y despejar. Se encuentra que
0,00289 m K−1
λmáximo =
T
Si todos colectáramos de alguna forma estos 1,38 kW m−2 fundimos a Edenor. Para encontrar la luminosidad del sol L
lo que tenemos que hacer es integrar este número en toda la superficie esférica que encierra al sol. Es decir multiplicamos
1,38 kW m−2 por la superficie de una esfera de radio 1 AU:
2
L = 1,38 kW m−2 × 4π (1 AU)
1 AU = 1,5 × 1011 m → ≈ 36 × 1022 kW
Según Wikipedia en el año 2005 la “potencia promedio mundial” fue de 15 TW = 1,5 × 109 kW.
Para calcular la temperatura del sol usamos Stephan-Boltzmann:
L 2
2 = σTefectiva
4πR
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

1.1 Cúmulos estelares 1 GUÍA 1 - EL UNIVERSO
Sale usando las leyes de Kepler. En particular la tercera ley dice que
4π 2
T2 = a3
G (m1 + m2 )
con lo cual
4π 2 a3
m1 + m2 =
GT 2
En el caso del sistema sol-tierra se considera que m m⊕ donde m⊕ es la masa de la tierra.
1.1. Cúmulos estelares

Son cúmulos de estrellas que están ligadas en forma gravitacional. Las estrellas de un cúmulo tienen todas la misma edad
y provienen de un mismo lugar. La edad de las estrellas se puede determinar en función de su color.
Hay dos tipos de cúmulos estelares:
Abiertos. Tienen del orden de 100 estrellas y son uniformes. Las estrellas suelen ser azules que son estrellas jóvenes. Aún
continúan formándose nuevos cúmulos abiertos. Se encuentran sobre el disco galáctico.
Globulares. Tienen del orden de 10 × 103 estrellas. Son de forma esférica y poseen gigantes rojas, que son estrellas viejas. Ya
no se siguen formando, los que hay se formaron hace mucho. No se encuentran sobre el disco galáctico sino que
están en el halo.
El movimiento propio es el movimiento que tiene un cuerpo celeste en la dirección perpendicular a r̂ siendo r̂ el versor
radial centrado en la tierra. El dibujito es así:
Punto de
convergencia
Cúmulo
Earth
Las velocidades en realidad no apuntan al punto de convergencia, es sólo una “ilusión óptica” por el hecho de que estamos
mirándolas sobre la esfera celeste. O algo así... La verdad que no entiendo bien qué es cada cosa...
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

2 GUÍA 2 - FLUIDOS CLÁSICOS
vradial
vtangencial
vt µd
tan 30 ◦ = = ⇒ d = 43 pc
vr vr
vt = µd
∆s = αd
α = 0,1100
Para calcular el tamaño simplemente hacemos
T
Earth
T/2
tan 10 ◦ = ⇒ T = 15 pc
d
Lo hicieron pero me tuve que ir a cursar E3.
2. Guía 2 - Fluidos clásicos

Tenemos la ecuación de Navier-Stokes que es

∂v
ρ + (v · ∇) v = −∇ · p + f → Navier-Stokes
∂t
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

donde v es el campo de velocidades, ρ es la densidad, p es el tensor de presiones y f es la fuerza externa por unidad de
volumen (por ejemplo un potencial gravitatorio en el 99 % de los casos). El término (v · ∇) v se conoce como “término
convectivo”.
Para fluidos ideales vale que p = p1. Esto implica que ∇ · p = ∇p con lo cual la ecuación de Navier-Stokes para fluidos
ideales es
∂v ∇p
+ (v · ∇) v = − + f → Navier-Stokes para fluidos ideales
∂t ρ
Obsérvese que se ha redefinido f , ahora es la densidad de fuerza por unidad de masa.
Tenemos un fluido de densidad uniforme ρ 6= ρ (r). La ecuación de Navier-Stokes es entonces
1
0 = − ∇p − ∇φ
ρ

p
Pues ρ 6= ρ (r) → = −∇ +φ
ρ
con lo cual
p
+ φ = constante
ρ
Si evaluamos esto sobre una equipotencial tenemos que φ es constante con lo cual
dpcen una equipotencial = otra constante ≡ p0
es decir que la presión es constante. Si ahora consideramos la ecuación de estado
T = T (p, ρ)
tenemos que evidentemente T es una constante sobre la equipotencial.
LTE significa Local Thermal Equilibrium.
Ítem a Planteamos Navier-Stokes en el caso estacionario y con la fuerza gravitatoria

1 MG
∇p = − 2 r̂
ρ r
Ahora usamos el hecho de que esto es un gas ideal con lo cual la ecuación de estado es
p = nkT
ρ
nm = ρ → = kT
m
y entonces
kT MG
∇ρ = − 2 r̂
mρ r
(se usó además que la atmósfera es isotérmica con lo cual ∇T = 0). Ahora resolvemos, la componente radial del gradiente
es:
1 kT ∂ρ MG
=− 2
ρ m ∂r r
y si ahora integramos desde la superficie del planeta R hasta un r cualquiera obtenemos

mM G 1 1
ρ = ρ0 exp −
kT r R
Ahora podemos calcular fácilmente la presión usando la ecuación de de estado. Es prácticamente lo mismo

mM G 1 1
p = p0 exp −
kT r R
Vamos a considerar ahora el límite de bajas alturas, donde tiene sentido que esto funcione. Vamos a considerar entonces
z R. En este caso queda que
kT
ρ = ρ0 e−z mg
kT
A la constante mg se la conoce como altura de escala.
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Ítem b La densidad de columna se define según

ˆ∞ ˆ∞
z
e− H ρ0 dz = −ρ0 H
def
ρ dz =
0 0
Ítem c Pasamos a un modelo de atmósfera un poco más realista (la isotérmica es re trucha) que es el de una atmósfera
adiabática en la que se satisface
pρ−γ = constante → Adiabática
c
con γ = cvp . El hecho de que la atmósfera sea adiabática implica que no hay intercambio de calor entre los elementos de
fluido. Ahora es repetir lo mismo de antes pero utilizando esta “nueva ecuación de estado”. La solución es
1
GM ργ0
γ−1
1 1
ρ = ρ0 1 + (γ − 1) −
γp0 r R
Tomando ahora el límite z = r − R → 0 (bajas alturas) se obtiene

1
gργ0 (γ − 1) γ−1

−


ρ (z) = ρ 0 1 z
 γp0
γ
γ γ−1


p (z) = ρ0 1 − gρ0 (γ − 1)
 z
γp0
Ahora podríamos buscar también la temperatura usando que T ∝ ρp .
Vamos a meter una perturbación en la atmósfera y ver qué ocurre. Vamos a asumir que hay simetría estacionaria en x e
y (salvo por la columna, claro). Entonces tamos a tener p (z) , T (z) y ρ (z). La ecuación de Navier-Stokes es entonces
dp
= −ρg
dz
o bien
dp ρg
= − dz
p p
ρ
Ahora consideramos un gas ideal con lo cual p = m kT y entonces
dp mg
= − dz
p kT
kT 1
≡H→ = − dz
mg H
A continuación usamos la ecuación de balance de calor que es

3 ∂T
nk + v 0 · ∇T + p∇ · v 0 = 0 → Eq. de balance de calor
2 ∂t
∂T
Reemplazando v 0 = v0 ẑ, usando situación estacionaria ⇒ ∂t = 0, obtenemos
3 ∂T ∂v0
nkv0 +p =0
2 ∂z ∂z
Por último vamos a usar la ecuación de continuidad del fluido
∂ρ
+ ∇ · (ρv 0 ) = 0 → Eq. de continuidad
∂t
En el caso estacionario y con simetría en x e y esto queda
∂ (ρv0 )
= 0
∂z
∂n ∂v0
v0 +n = ← nm = ρ
∂z ∂z
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Reemplazando en la ecuación de balance de calor

3 ∂T v0 ∂n
nkv0 −p =0
2 ∂z n ∂z
p
y ahora usamos n = kT y no sé qué más hacemos y nos queda
 
pv0 ∂n  1 ∂p p ∂T 
− = −kT v0  −
 
n ∂z  kT |{z}
∂z kT 2 ∂z 

p
−H

3 ∂T p p ∂T
nkv0 − kT v0 − − =0
2 ∂z HkT kT 2 ∂z

3 1 ∂T 1 1 ∂T
v0 p + + =0
2 T ∂z H T ∂z
5 k
cp =
2m
mcp ∂T 1
+ =0
kT ∂z H
y esto es lo que queríamos demostrar (?), la ecuación de balance de calor para la columna.
Ítem b La variación de la temperatura es
dT ∂T
+ v · ∇T
= 0
dt ∂t
∂T
= v0
∂z

kT
= v0 −
mcp H
g
= −v0
cp
donde hemos reemplazado la expresión de H.
Ítem c Nosotros tenemos que 

dT −α Afuera de la columna
= g
dz  Adentro de la columna
cp
En el caso en que α < cgp lo que pasa es que los elementos de fluido que suben por la columna se enfrían más rápido que la
atmósfera circundante. En este caso la densidad va a aumentar y el fluido se irá nuevamente hacia abajo. En caso contrario,
i.e. cuando α > cgp , lo que ocurre es que los elementos de la columna están siempre más calientes dentro de la columna que
afuera con lo cual se enfriarán más lento y subirán más rápido. En este caso se desarrolla la inestabilidad convectiva.
Vamos a estudiar la deformación de un planeta producto de la rotación.
Ítem a La fuerza centrífuga es

F centrífuga = −ρω × (ω × r)
con ω = ωẑ la velocidad angular. Expresando en polares

ω = ω cos θ r̂ − sin θ θ̂
nos queda que

F centrífuga = ρω 2 r r̂ sin2 θ + θ̂ sin θ cos θ
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Por otro lado la fuerza gravitatoria es
M (r) G

−ρ
 r̂ para r < R
F gravedad = r2
 − ρ M G r̂

para r ≥ R
r2
La ecuación de Navier-Stokes para un fluido bonito nos queda
∇p = −∇φ
El potencial es entonces 
4 r2 r2
φexterior = πρ2 G − ρω 2 sin2 θ + constante


3 2 2
2
interior ρM G 2 r
sin2 θ

φ
 =− − ρω
r 2
Aplicando la condición de contorno de que φexterior = φinterior para r = R se despeja la constante que termina valiendo
− 32 ρM G
R .
2
Ítem b φ (r = R) = − ρM G 2R 2
R − ρω 2 sin θ
∇ (p + φ) = 0 ⇒ p + φ = constante
dpcsuperficie = p0 = constante
dpcsuperficie + dφcsuperficie = constante

con lo cual
dφcsupefficie = constante
2
Ahora lo que vamos a hacer es evaluar φ (r = R) = − ρM G 2R 2
R − ρω 2 sin θ en un punto de referencia, por ejemplo el polo con
R0 y θ = 0. Obtenemos
ρM G R2 ρM G
− − ρω 2 sin2 θ = −
R 2 R0
de donde podemos despejar
R ω 2 R3 sin2 θ
=1+
R0 2GM
R
Para seguir resolviendo eso habría que integrar en forma numérica. A mano llegamos hasta acá. Lo importante es que R0 ≥1
con lo cual cualquier R es mayor que el R0 de referencia tomado en el polo.
Ítem c Vamos a definir

δR = Recuador − R0 R0
Ahora me perdí (estamos resolviendo esto a las chapas en los últimos 10 minutos de clase)
δR Re − R0
=
R0 R0
ω 2 Re3
=
2M G
3
ω 2 (R0 + δR)
=
2GM
ω 2 R03 + 3R02 δR
≈
2GM
ω 2 R03
=
2GM − 3ω 2 R0
ω 2 R03
= 1+
2GM − 3ω 2 R0
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

3 GUÍA 3 - RELATIVIDAD
M = 5,68 × 1026 kg




ω = 1,64 × 10−4 s−1

Ítem d Saturno es el planeta más achatado del sistema solar. Usando los valores obtenemos
R0 = 5,43 × 107 m



Recuador = 6,02 × 107 m


Recuador

 = 1,068
 R0 teórico

 Recuador

 = 1,108
R0 observado
donde el teórico es con la cuenta del ítem previo mientras que el observado es haciendo el cociente entre las cantidades
numéricas de este ítem.
3. Guía 3 - Relatividad
Vamos a comenzar haciendo un repaso. Consideremos dos sistemas de referencia, S y S 0 , tal que S 0 se mueve con
velocidad v = v x̂ visto desde el sistema S (el caso típico de todo libro). En este caso las coordenadas se relacionan mediante
la transformación de Lorentz 0
 t = γ (x − βct)

x0 = γ (ct − βx)

Lorentz boost →
y 0 = y


 0
x =x
v 1
donde β = c yγ=√ . En forma matricial lo anterior es
1−β 2
 0   
t γ −γβ t
x0  −γβ γ  x
 0 =   
y   1  y 
z0 1 z
y en notación de índices
x0µ = Λµ ν xν → Notación de índices
donde se aplica el criterio de sumación de Einstein.
Se define el cuadrivector momento
E
c
µ
px 
p = 
py 

pz
γmc2

=
γmv
y se obtiene la relación
E 2 = m2 c4 + c2 p2
También tenemos el efecto Doppler relativista
νobservada = νemitida (1 − β cos θ) γ → Doppler relativista
donde θ es el siguiente ángulo
← observador
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Ítem a El factor de Lorentz es γ. La componente temporal del cuadrimomento es p0 = γmc2 ≡ E. Conociendo la energía
(de los rayos cósmicos) 1020 eV, la masa del protón y la velocidad de la luz podemos calcular trivialmente que
γ ∼ 1011
Ítem b La contracción de la longitud es

L0
L =
γ
100 pc
=
1011
= 10−9 pc
18
1 pc = 3 × 10 cm → ∼ 109 cm
= 104 km
Eso es menos que la distancia entre la tierra y la luna.

El tiempo propio que le demanda al protón atravesar la galaxia es
t
t0 =
γ
L0
=
γvp
L
vp ≈ c → ≈
c
104 km
=
108 m s−1
= 0,1 s
El tiempo medido por un observador en la tierra va a ser

100 pc
ten la tierra =
c
≈ 300 años
Ítem c Hay que hacer el camino inverso a lo de recién. No entiendo por qué puso que
t = 2τn
t 0 t
Nn
que lo sacó haciendo Nn (t) = Nn (0) 2− τn = 4 ⇒ 2− τn = 14 . Ahora es despejar la formulita: la distancia recorrida por el
neutrón es
En
d= 2τn
mn c
donde hemos usado que vn ≈ c y En = γn mn c2 y d = cγn 2τn .
Ítem a Queremos estudiar el efecto Compton y relacionar la frecuencia del fotón final con el ángulo de escatereo. Por
conservación de momento sabemos que
pγ − pγ final = pe− final
con lo cual
p2γ + p2γ final − 2pγ pγ final cos θ = p2e− final
En cuanto a la energía tenemos que q

hν + mc2 = hνfinal + m2 c4 + c2 p2e− final
Juntando las dos coas tenemos

2
2h (ν − νfinal ) mc2 + h2 (ν − νfinal ) +
m2
c4 =
m2
c4 + c2 p2γ + c2 p2γ

final − 2c2 pγ pγ final cos θ
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Usando ahora que E = cp para los fotones tenemos que
mc2 2
2 (ν − νfinal ) + (ν − νfinal ) = ν 2 + νfinal
2
+ 2ννfinal − 2ννfinal cos θ
h
mc2
(ν − νfinal ) = ννfinal cos θ
h
y finalmente
1 1 h
− = (1 − cos θ)
νfinal ν mc2
Ítem b Ahora tenemos el efecto Compton inverso, es decir un electrón que se mueve a las chapas y choca a un fotón. Lo
único que hay que hacer es aplicar un boost tal que nos vamos al sistema en reposo del electrón. Luego tenemos el mismo
problema que antes. Entonces tenemos que en el sistema en reposo vale lo mismo que en el ítem a:
0 1
νfinal = → En el sistema en reposo
1
ν0 + h
mc2 1 − cos θ0
Ahora tenemos que aplicar un boost al sistema laboratorio para ver cómo queda esto. En primer lugar transformamos la
frecuencia usando el Doppler relativista
0
γ 1 + β cos θ0

νfinal = νfinal
y entonces
γ 1 + β cos θ0

νfinal = 1 h
1 − cos θ0

ν0 + mc2
Ahora tenemos un ángulo φ que no sé qué representa pero aparece en una fórmula que es
ν 0 = νγ (1 − β cos φ)
Creo que es el ángulo en el que observamos el proceso en el sistema laboratorio, pero no sé con respecto a qué. No hay ningún
dibujo que nos ayude a interpretar qué es. Lo que le quedó es
γ 1 + β cos θ0

νfinal = 1 h 0

νγ(1−β cos θ) + mc2 1 − cos θ
Ítem c Queremos encontrar la máxima energía del fotón dispersado en el ítem b. No entendí nada. Hay dos ángulos φ y θ
que no sé qué son, y dadas las preguntas del curso creo que nadie entiende bien qué son.
El potencial gravitatorio en la superficie de la enana blanca es
GMenanita
dφcSuperficie de la enana blanca = −
Renanita
mientras que
dφcSuperficie de la tierra ≈ 0 → Debido a que MTierra ≪ MEnana y RTierra ∼ REnana
por lo tanto
GM
∆φ =
R
Lo que sigue no le entendí nada...
hc hc
∆E = −
λobservado λ0
∆λg = λobservado − λ0
def ∆λ
z= → Red shift
λ0
∆λg GM
= = 1,49 × 10−4 → Red shift gravitatorio
λ0 Rc2
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

4 GUÍA 4 - ATMÓSFERAS ESTELARES
Tenemos que ver cuál es la energía mínima del γ para que pueda ocurrir la reacción p + γ → n + π + . Vamos a usar por
un lado la conservación de cuadrimomento
pµinicial = pµfinal → Por conservación
y por otro el invariante relativista
p2 = pµ pµ → Invariante relativista
Entonces
(pinicial ) µ (pinicial ) µ = (pfinal ) µ (pfinal ) µ
En particular podemos evaluar el pinicial en el frame en que el protón está en reposo y el pfinal en el sistema centro de masa.
Esto es 2
2
dpinicial cSistema en el que el protón está en reposo = (dpfinal cSistema centro de masa )
Las componentes de cada bicho son
Efotón + mp c2
  


  Efotón 
dpinicial cSistema en el que el protón está en reposo =  c

 



  0 

 0
2 2
 


 mn c + mπ c

  0 
 dpfinal cSistema centro de masa = 

 



 0 
0
donde hemos asumido que la energía del fotón es la mínima tal que las partículas finales están en reposo. Entonces tenemos
que
2
2 Efotón 2
Efotón + mp c2 − ≥ (mn + mπ ) c4
c
4. Guía 4 - Atmósferas estelares

Tenemos una función de distribución de fotones
f (r, p) → Distribución de fotones
que es el número de fotones en la posición r con un impulso p. Además tenemos la intensidad de radiación
Iν → Intensidad específica
Energía radiada
que es la Área×Tiempo×Ángulo sólido×Frecuencia
. Queremos relacionar las dos cantidades. Para ello consideramos el hecho de que
el impulso del fotón es
hν
p=
c
con lo cual (no entiendo)
con lo cual
Iν
= #de fotones
hν
y finalmente
h4 ν 3
Iν = f
Willy!
c2
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Ítem a Creo que no es más que la conservación de la energía. La energía total que se encuentra en la cáscara esférica que
se propaga en forma radial desde la estrella hacia la tierra es constante. Esto es justamente la luminosidad de la estrella.
Entonces tenemos que
(
4πR2 qν (R) Evaluada sobre la superficie de la estrella
Luminosidad → Lν =
4πr2 qν (r) Evaluada en un r arbitrario
entonces 2
R
qν (r) = qν (R)
r
Ítem b La intensidad en la superficie de la estrella es

(
Iν para θ < π/2
Iν (θ) =
0 en otro caso
Gráficamente esto sería así:

Este lado irradia
Ahora sólo hay que integrar

ˆ
qν = Iν ẑ cos θ dΩ
ˆ2π ˆπ/2
= Iν dφ cos θ sin θ dθẑ
0 0
= πIν ẑ
Para calcular las demás cantidades simplemente recurrimos a su definición

ˆ
1
Jν = Iν dΩ
4π
y ˆ
pij
ν = Iν ui uj dΩ
Ítem c Ahora asumimos que

2hν 3 hν −1
Iν = Bν = e kT − 1
c2
Tenemos que la luminosidad por unidad de frecuencia es
Lν = 4πR2 qν (R)
= 4π 2 R2 Iν
3 −1
2 2hν hν
= (2πR) 2
e kT − 1
c
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

La luminosidad total es simplemente la integral de esto anterior sobre todo el espectro, es decir
ˆ∞
L = Lν dν
0
ˆ∞
2h 2 ν3
= (2πR) 2 hν dν
c e kT −1
0
hν 

x=
def
4 ˆ∞
x3

kT →
 2h 2 kT
= (2πR) 2 dx
h  c h ex−1
dx = dν  0
kT
ˆ∞
x3 π4 2π 5 k 4 4
dx = → = 4πR2 T
ex−1 15 15h3 c2
0
= 4πR2 σStephan-Boltzmann T 4
Dada una magnitud patrón m0 (estrella Vega) entonces la magnitud de otra estrella se calculaba usando

q
m = m0 − log2,5 → Definición antigua
q0
Esto viene de una definición que hizo un tal Hiparco en la antigua Grecia. Sin embargo es un bajón trabajar con un logaritmo
en base 2,5. En consecuencia se redefinió usando que
1/5
(100) = 2,512 . . .
con lo cual ahora se definió

q
m = m0 − log(100)1/5
q0

log qq0
= m− 1/5

log (100)
y se termina obteniendo la definición que usamos hoy en día

q
Magnitud aparente → m = m0 − 2,5 log → Definición moderna
q0
Todo esto lo podemos hacer en función de ν, en cuyo caso obtendríamos

qν
mν = mν0 − 2,5 log
qν0
Para pasar de longitud de frecuencia a longitud de onda usamos

dλ
qλ | dλ | = qν | dν | ⇒ qν = qλ

dν
(las barras de módulo son para evitar quilombos con los signos). Usando ahora λν = c encontramos que
λ2
qν = qλ
c
Ahora vamos a trabajar en la banda azul B. Usando el dato de que
mν,Vega ≡ 0 ∀ν → Por consigna
tenemos que
qν,?
mB,? = −2,5 log
qν,Vega
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!


mB,?
qν,? = qν,Vega 10− 2,5
λ2 λ2 mB,?
qν = qλ → = qλ,Vega 10− 2,5
c c
Poniendo todos los datos que figuran en la tablita de la consigna termina quedando que
qB,? = . . .
Ahora sólo habría que repetir este proceso para cada punto de la tablita y hacer el grafiquito.
Ítem b La magnitud absoluta es la magnitud con la que mediría a una estrella si ésta se encontrare a una distancia de 10 pc.
Para encontrar cómo se relaciona la magnitud absoluta M con la magnitud aparente m recordemos en primer lugar que el
flujo medido a distintas distancias se relaciona según
2
R
q (r) = q (R)
r
Entonces la relación que hay entre las magnitudes es

q (10 pc)
M = m − 2,5 log → Magnitud absoluta
q (r)
2
q(10 pc) r R 2

Despejando q(r) = 10 pc a partir de q (r) = q (R) r lo que encontramos es que

r
M = m − 5 log → Magnitud absoluta
10 pc
De aquí podemos despejar

m−M
r = 10 pc × 10 5
Los valores de m los tenemos

en la tablita
de la consigna. Los valores de M los obtenemos usando los q del ítem previo en
la fórmula M = m − 2,5 log q(10 pc)
q(r) .
Por otro lado sabemos que la luminosidad es
Lν = 4πr2 qν → Luminosidad
por lo tanto usando los qν del ítem previo y el r que recién calculamos, podemos encontrar la luminosidad.
Ítem c La luminosidad de Planck es

2 2h ν 3
Lν = (2πR) hν → Planck
c2 e kT −1
Haciendo el ajuste de encuentra que (
T ∼ 60 kK
R ∼ 8 × 108 m
Tenemos una estrella que tiene más o menos el tamaño del sol pero es mucho más caliente.
Ítem d En el ejercicio anterior, o en el anterior anterior, encontramos que para un objeto que sigue la ley de Planck vale
que su luminosidad es
L = 4πR2 σT 4
Teniendo R y T del ítem previo sólo hay que hacer la cuenta.
En cuanto a la magnitud bolométrica usamos

L
MBolométrica = MBolométrica − 2,5 log
| {z } L
Referencia
Sólo hay que hacer la cuentita.
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Ítem a Comencemos calculando la presión de radiación
tr (pij )
pν =
ˆ 3ˆ
1 1
pij = Jν ui uj dΩ → = Jν ûi ûj dΩ
c 3c
4πJν
=
3c
4πIν
Jν = Iν por isotropía → =
3c
1
´
(recordemos que Jν = 4π Iν dΩ = Iν cuando hay isotropía). Ahora recordemos que
(
qν (R) = 4πIν → Flujo en la superficie de la estrella
Lν = 4πR2 qν (R) → Luminosidad
Usando esto podemos ver que

Lν
pν =
12cπR2
por lo tanto la presión de radiación es
ˆ∞
p = pν dν → Presión de radiación
0
L
=
12πcR2
Ahora podemos calcular la fuerza de radiación como
Frad = σp → Fuerza por unidad de volumen
con lo cual
Lσ
r̂ F rad =
12πcR2
Si bien los fotones están haciendo un movimiento aleatorio en este gas, el gradiente radial de presiones termina generando
esta fuerza.
Ítem b Queremos ver quién gana, si la gravedad o la presión electromagnética. Esto es

Lσ
Frad 2
R
12πc
= GM
g
2
R

Lσ
=
12πcGM
L
∝
M
Aquí lo que vemos es que dependiendo de la relación entre L y M vamos a tener la estabilidad de la estrella:
Frad
Si g > 0 entonces la estrella expulsa material.
Frad
Si g < 0 la estrella colapsa sobre si misma.
Ítem c Vamos a considerar el scattering de Thomson que es el límite de baja energía del scattering de Compton. La
característica de este scattering es que es elástico y se conserva la energía de cada uno de los elementos por separado.
Si la sección eficaz por cada electrón es σT (según consigna) lo que tenemos que hacer es encontrar la sección eficaz por
unidad de volumen. Esto es simplemente multiplicarla por la densidad numérica de electrones, es decir
σ = ne σT → Sección eficaz por unidad de volumen
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Entonces
Frad Lσ
= ← Del ítem previo
g 12πcGM
Lne σT
=
12πcGM
Ahora consideramos la condición de cuasi neutralidad que implica que
ne ≈ np → Cuasi neutralidad
Por último simplemente consideramos la densidad de masa:

ρp
np =
mp
Reemplazando obtenemos
g 12πGM mp c
=
Frad σT Lρp
Ítem d El límite de Eddington relaciona la máxima luminosidad de una estrella en función de la masa. Para encontrarlo
vamos a considerar lo siguiente
Fgravedad
La estrella no colapsa ⇐⇒ <1
Fradiación
por lo tanto
gmp
=1
Frad Límite de Eddington
12πcGM mp
L< → Límite de Eddington
σT
Finalmente consideramos 3,5 3,5
L M M
= ⇒ L = L
L M M
por lo tanto
M < despejar de lo de antes
Ítem a y b (?) La ecuación de transferencia radiativa es

1 ∂Iν
+ (u · ∇) Iν = ρ [eν + σν Jν − kν Iν − σν Iν ]
c ∂t
donde eν es la emisividad, σν es el coeficiente de scattering, kν es la absorción. Como estamos en el caso estacionario entonces
∂Iν
∂t ≡ 0.
Theta
En este caso el observador está en (

z=∞
Observador →
τν = 0
y
dτν = −κν ρ dz
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Además se define (igual que en la teórica)

def
µ = cos θ
y la función fuente
eν + σν Jν
Sν = → Función fuente
kν + σν
Entonces
τνdIν τν
e− µ µ = (Iν − Sν ) e− µ
dτν
d − τµν τν
µ Iν e = −Sν e− µ
dτν
Iν (z = ∞) = Iν (τν = 0)
ˆ∞
1
Sν e− /µ dt
t
=
µ
0
Ítem c Asumiendo que

Sν ≈ aτν + b → Eddington-Barbier
podemos hacer
Iν (z = ∞) = Iν (τν = 0)
ˆ∞
1
Sν e− /µ dt
t
=
µ
0
1
bµ + aµ2

=
µ2
= aµ + b = Sν (τν = 1)
Entonces en definitiva lo que tenemos es que
Iν = aµ + b → En el observador
Ahora vamos a usar esto para explicar el oscurecimiento del limbo oscurecimiento del limbo
En este problema tenemos lo mismo de siempre pero al revés:
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

5 GUÍA 5 - MEDIOS INTERESTELARES
m_{\nu,0}
Theta
Tierra
Estamos interesados en conocer cómo afecta la atmósfera así después podemos eliminar su influencia en la observación.
Entonces si mν es la magnitud de la estrella “si la atmósfera no estuviera”, o sea la posta, vamos a tener que
mν,suelo = mν,0 + κν cos θ → Consigna
siendo κν la opacidad de la atmósfera terrestre.

Haciendo lo mismo que en el problema de hace un rato llegamos a
Iν,suelo = Iν,0 e−τν sec θ
donde τν es la profundidad óptica cenital. Además el flujo es
qν,suelo = qν,0 e−τν sec θ

qν,0
mν,0 = mν,ref − 2,5 log
qν,ref

qν,suelo
mν,suelo = mν,ref − 2,5 log
qν,ref
Entonces

qν,suelo
mν,suelo − mν,0 = −2,5 log
qν,0
−2,5 log e−τν sec θ

=
= 2,5τν log10 (e) sec θ → Es lo que pedía el problema!
Para medir mν,0 desde la superficie de la tierra (o sea, midiendo mν,suelo ) se toman muchas mediciones para distintos valores
de θ y luego se ajusta la función
mν,suelo = mν,0 + κν sec θ
a los datos medidos. La ordenada al origen es el parámetro de interés.
5. Guía 5 - Medios interestelares

Esta es una guía re tranqui, como para tener idea de órdenes de magnitud y nada más.
Salen todos con la siguiente formulita: r
8kT
hvi =
πm
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Del problema anterior tomamos que la temperatura de las regiones H I es de TH I ∼ 100 K. Por equipartición de energía
tenemos que
m
2 3
v = kT
2 2
Para la velocidad de rotación es similar
I
2 3
ω = kT
2 2
siendo I el momento de inercia.
Para calcular el momento de inercia I y la masa m de las partículas asumimos que éstas son esferas con las características
que nos da la consigna.
Ítem a Una forma de hacerlo es al estilo de Física 3: tomar un elemento de volumen de la esfera de gas y estudiar cómo
interactúa. Tipo ley de Gauss. Consideremos el siguiente dibujito
Elemento de masa
δm Nube de gas de radio R

r
El delta de energía que aporta este delta de masa es
Gm (r) δm
δE = −
r
Gm (r) ρ (r) δv
=−
r
siendo m (r) la masa encerrada en la esfera de radio r y δv el diferencial de volumen. Integrando lo anterior sobre todo el
volumen de la esfera obtenemos
ˆ
E = δE
ˆ
Gm (r) ρ (r) 3
= − d r
r
ˆ
def R m (r) ρ (r) 3 ΘGM 2
Θ= d r → = −
M2 r R
Si ahora asumimos que la densidad ρ (r) es uniforme ρ0 entonces
ˆR
R 4
Θ = 2 ρ20 4π πr4 dr
M 3
0
4 3
M = πR3 ρ0 → =
3 5
con lo cual
3GM 2
E=− → Para esfera uniforme
Willy!
5R
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Ítem b Tenemos que usar la ecuación de un fluido ideal
∇p + ρ∇φ = 0 → Balance hidrostático
Para encontrar la presión hacemos el siguiente truco: primero multiplicamos escalarmente por r de modo tal que nos quede
una ecuación escalar, y además integramos en todo el volumen así podemos transformar las integrales de volumen en integrales
de superficie donde nos quede la presión sobre la superficie que es lo que buscamos. Esto es
ˆ
(∇p + ρ∇φ) · r d3 r = 0
ˆ
(∇p · r + ρ∇φ · r) d3 r =
Ahora consideramos
ˆ ˆ ˆ ∇·r≡3
z }| {
3 3
∇p · r d r = ∇ · (pr) d r − p ∇ · r d3 r
V V
) ‹ ˆ
Teorema de la divergencia 2
→ = pr · r̂R dΩ − 3kT n d3 r
pV = nkT
∂V V
| {z }
N
dpc∂V ≡ p0 → = p0 4πR3 − 3kT N
= p0 4πR3 − 2cV T M
Por otro lado la integral con φ es

ˆ ˆ
3 Gm (r)
ρ∇φ · r d r = ρ r̂ · r d3 r
r
ΘM 2 G
Θ = ··· → =
R
Reemplazando estas dos expresiones finalmente podemos despejar
cV M T ΘGM 2
p0 = −
2πR3 4πR4
Ítem c La forma funcional es algo así
RJ
Para encontrar el máximo simplemente derivamos:

∂p0 3cV M T 4ΘGM 2
=− − =0
∂R 2πR4 3πR5
De aquí se despeja el radio de Jeans que es
2ΘGM
RJ = −
3cV T
4ΘGM
=−
ρkT
El valor de presión máxima es
37

kT
P0M = 10 3 2 3 → No se llega a ver bien, quizá está mal copiada
2 Θ M G m
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

5.1 Rankine-Hugoniot 5 GUÍA 5 - MEDIOS INTERESTELARES
y de aquí se puede despejar la masa de Jeans que es

3
10125 kT 1
MJ2 =
π1024 mG ρ
Podemos analizar el equilibrio y vamos a encontrar que
Colapso Acá es
gravitatorio estable
Presión
Radio de la nube
RJ
Todo esto parece que está en el Battaner.
5.1. Rankine-Hugoniot
Las relaciones de Rankine-Hugoniot son


 ρ1 u1 = ρ2 u2

p1 + ρ1 u21 = p2 + ρ2 u22

Rankine-Hugoniot →
 u2 5 p1 u2 5 p2
 1+ = 2+


2 2 ρ1 2 2 ρ2
γ
donde el 52 = γ−1 con γ = 1 + n2 siendo n el número de grados de libertad de la molécula. Las relaciones de Rankine-Hugoniot
sólo son válidas en el referencial de la onda de choque!
Tenemos
p2,rho2,T2 p1,rho1,T1
En el referencial del
fluido no chocado
onda de choque
Si nos pasamos al referencial de la onda de choque vamos a tener que el fluido es el que se está moviendo, y es en este
referencial en el cual valen las relaciones de Rankine-Hugoniot.
Vamos a definir, igual que en la teórica, s
 ρ2
ψ =



 ρ1
 u
M= → Número de mach
 c
γp


 c2 = → Velocidad del sonido


ρ
Vamos a expresar todo en términos del número de mach del lado no chocado, es decir

 p2 = p2 (M1 )


T = T (M )
2 2 1

 ρ2 = ρ2 (M1 )

M2 = M2 (M1 )

Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

5.1 Rankine-Hugoniot 5 GUÍA 5 - MEDIOS INTERESTELARES
Vamos a considerar la ecuación de estado de un gas ideal

nkT → Gas ideal p =
ρ
= kT
m
Combinando esta ecuación de estado con la ecuación de la velocidad del sonido podemos obtener
k
T c2 = γ
m
Ahora de la primera ecuación de Rankine-Hugoniot metiendo la ecuación de estado obtenemos
p2 m p1 m
u2 = u1
kT2 kT1
Combinándola con no sé qué obtenemos
2 2
p2 u2 1 p1 u1 1
=
c2 T2 c1 T1
y entonces tenemos que
2
T2 p2 M2
=
T1 p1 M1
Por otro lado, de la tercera relación de Rankine-Hugoniot junto con la ecuación de estado y la de la velocidad del sonido
tenemos que ! !
5 T1 k u21 5 T2 k u22
1+ k
= 1+ k
2 m 5T1 m 2 m 5T2 m
k

3c22 = 5T2

Usando ahora que γ = 53 (pues tenemos gas monoatómico o algo de eso) y que m vamos a llegar a
3c2 = 5T k

1 1
m
M21 M22

T1 1 + = T2 1 +
3 3
Por último consideremos la ecuación segunda de Rankine-Hugoniot. Combinándola con la ecuación de estado y la de la
velocidad del sonido vamos a encontrar que
p2 1 + 53 M21
=
p1 1 + 53 M22
2
Juntando esto con TT21 = pp21 M
M1
2
encontramos
2
1 + 13 M22

p2 M2
=
p1 M1 1 + 13 M21
M21 M22

Si ahora usamos T1 1 + 3 = T2 1 + 3 finalmente llegamos a
2 2
1 + 53 M21 1 + 13 M22

M1
=
1 + 35 M22 M2 1 + 31 M21
Como se puede ver esto es simétrico en M21 y M22 . Esto implica que M1 = ±M2 son soluciones. Si despejamos acá
encontramos
1 5
M21 − M22 1 + M21 + M22 − M21 M22 = 0

3 3
Las soluciones M1 = ±M2 no nos interesan ya que son triviales. Las soluciones interesantes salen de 1 + 13 M21 + M22 −

5 2 2
3 M1 M2 = 0. De aquí sacamos
M21 + 3
M2 = → Solución
5M21 − 1
Sabiendo el valor de M2 ahora podemos calcular todo lo demás

5 2 1

 p 2 = p 1 M −
4 1 4





4M2


ρ2 = ρ1 2 1

 M1 + 3

3 + M21 5M21 − 1



T2 = T1

16M21
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

6 GUÍA 6 - ESTRUCTURA Y EVOLUCIÓN ESTELAR
6. Guía 6 - Estructura y evolución estelar

Bibliografía:
Carol capítulo 11 (el mejor)
Long gate high eergy astro phisicc 2
Bataner cap 5
GM 2 Energía gravitatoria
τKH = = ∼ 107 años
RL Luminosidad
r
R3
τH = ∼ 1600 s
GM
QM c2
τm = ∼ 1011 años
L
El gráfico lo podemos encontrar en el capítulo 4 del Choduli o en el 18 del Kipenjam. Es algo así
56
El motivo por el cual las reacciones se cortan en F e es que necesitamos que el resultado tenga menor energía que el
estado inicial Entonces sólo podemos proceder así:
Sube la energía de ligadura

por lo tanto habremos liberado
energía
La fusión La fisión libera

libera energía energía1
Debido a que a la izquierda del F e los átomos tienen cada vez más energía de ligadura entonces cuando se fusionan dos
átomos el resultado será un átomo más estable que los originales. En cambio del lado derecho del F e lo que pasa es que si
hacemos una fusión vamos a obtener un producto que es más inestable que los originales, entonces se va a fisionar.
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

6.1 Resumen 6 GUÍA 6 - ESTRUCTURA Y EVOLUCIÓN ESTELAR
6.1. Resumen
El siguiente resumen se puede encontrar en la bibliografía [Carroll and Ostlie, 2017, p. 420].
Ec. estructura
• Equilibrio hidrostático relaciona p con r [Carroll and Ostlie, 2017, p. 377]
dp GMr ρ
=− 2
dr r
• Conservación de masa relaciona M con r [Carroll and Ostlie, 2017, p. 378]
dMr
= 4πr2 ρ
dr
• Generación de energía relaciona L con r [Carroll and Ostlie, 2017, p. 398]
dL
= 4πr2 ρε
dr
(
gravitatoria
donde ε → .
nuclear
• Transporte de energía relaciona T con r [Carroll and Ostlie, 2017, p. 406]
dT 3 κR ρ Lr
=−
dr 4ac T 3 4πr2
donde κR es la opacidad de un señor que se pronuncia “rosland”. WARNING: los símbolos pueden estar mal
copiados ya que no se entienden bien del pizarrón.
Relaciones constitutivas 
P (ρ, T, Xi )

ε (ρ, T, Xi )

κR (ρ, T, Xi )

En general se considera un gas ideal con lo cual la relación constitutiva para la presión es
ρkT 1
P = + aT 4 → Gas ideal (?)
µmH 3
Para ε se suele considerar la reacción protón protón, CN0 (o CNO, no sé). Para κR no vamos a considerar nada realista.
Condiciones de contorno
• Para r → 0 tenemos que ME BORRÓ TODO EL PIZARRÓN.
Tenemos 
γ
p = kρ

→ Consigna

1 d dψ

 2 x2 + ψn = 0
x dx dx
donde r = αx, ρ = ρ0 ψ n y γ = 1 + n1 .
Ítem a La solución es del tipo críptica en el sentido de que es una sarta de cuentas a interpretar por el alumno. La solución
del ítem (a) es la siguiente:
1 dp Gm (r)
+ =0
ρ dr r2
ˆr
4πr03 ρ r0 dr0

m (r) =
0
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

ˆr
1 dp G
4πr02 ρ r0 dr0 = 0

⇒ + 2
ρ dr r
0
p = kργ
ρ = ρ0 ψ →n
= kργ0 ψ nγ
= kργ0 ψ n+1
dp dψ
= kργ0 (n + 1) ψ n
dr dr
ˆ r
r2 dψ
kργ (n + 1) ψ n + 4πG r02 ρ r0 dr0 = 0

ρ0 ψ n 0 dr
0
kργ−2

(n + 1) d 2 dψ
0
r + ψn = 0
4πGr2 dr dr
ρ (r = 0) = ρ0 ⇒ ψ (x = 0) = 1
d2 ψ 2 dψ
+ + ψn = 0
dx2 x dx

dψ
=0
dx x=0
Ítem b
kρ0γ−2 (n + 1)
α2 =
4πG
4πGα2 = kργ−2
0 (n + 1)
1 1
γ =1+ → γ − 2 = −1 +
n n
1
−1+ n
4πGα2 = kρ0 (n + 1)
Ítem c
ψ (x) = 1 + b2 x2 + b4 x4 + b6 x6
Ahora simplemente metemos esta solución en la ecuación de Lane-Emden [Carroll and Ostlie, 2017, p. 426] y ver cuánto
valen las bi . Parece que queda que
1


 b2 = −

 6
n


b4 =

 120
2
 b6 = n − n



3024 1890
1/m
ρkT mp mkρ0
(a) p = m →T = ρk = k ψ
16σT 3 dT
q=− r̂
3ρcκR dr
1
ε (r) = ∇·q
ρ
1 1 d 2
= r q
ρ r2 dr
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!


2qc C
ε (x) = − A+B+ ψ 3−2m
αρ0 4
donde
16σ 4−n nK
qc = − ρ n
3cκR 0 k
y A, B y C son otras cosas así de feas que no tiene sentido copiar...
MN RN
,
M? R?
ˆx1
M? = ρ (x) dx con x1 α = R?
0
ˆx1
I (x1 ) = ε (x) dx
0
ˆx2
x2 tal que ε (x) dx = 0,9I (x1 )
0
x2 α = RN
ˆx2
MN = ρ (x) dx
0
Ítem e Masa de un politropa.
ˆR?
M = 4πρ (r) r2 dr
0
ˆx1
3
r = αx → = 4πα ρ (x) x2 dx
0
ˆx1
= 4πα3 ρ0 ψ n x2 dx
0

n 2 d 2 dψ
ψ x =− x
dx dx

dψ
M = −4πα3 ρ0 x2
dx x1
" #3/2
kργ−2

0 (n + 1) 2 dψ
= −4πρ0 x
4πG dx x1
c
Estrellas adiabáticas tenemos gas monoatómico γ = cVp = 35 y entonces n = 1,5.
ρ
 kT = βp
m


 4 1/3
4 3k 1−β
Estrellas con presión de radiación. aT = (1 − β) p → p = am 4 β4 ρ4/3 → n = 3 .



 3
0<β<1

h2 3
2/3 ? 5/3
Gas de Fermi no relativista p = 5me 8π ne → n = 1,5
1 3
1/3 4/3
Gas de Fermi relativista p = 8 π chne → n = 3.
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

7 GUÍA 7 - CAMPOS MAGNÉTICOS
7. Guía 7 - Campos magnéticos

Hasta ahora hemos ignorado completamente los efectos de los campos magnéticos. Pero son importantes.
Actividad solar o magnética Podemos distinguir fenómenos transitorios o periódicos
Transitorios. Tienen que ver con la reconexión magnética de Sweet-Parker.

• Fulguraciones. Son liberaciones súbitas de energía, del orden de 1032 erg en 100 s.
• Eyecciones coronales de masa. Se emiten del orden de 1012 kg en un tiempo de 10 minutos.
Periódicos.
• Ciclo de manchas solares.
Todos estos fenómenos se conocen en detalle para el sol, y se presume que ocurren en todas las estrellas de similares
características.
Un pulsar es una estrella de neutrones que tiene un período de rotación muy corto, típicamente del orden del segundo. El
motivo por el cual rotan tan rápidamente es que tienen un radio muy pequeño, entonces la conservación del impulso angular
impone que su velocidad angular debe ser muy elevada.
Ítem a La energía de rotación de una esfera es

I 2
E = ω
2
2 M R2 2
I= M R2 para esfera → = ω
5 5
Entonces
dE M R2 dω
= 2ω
dt 5 dt
2π 2M R2 2π 2π
ω= → =
P 5 P P2
(estamos poniendo sólo los módulos, no ponemos los signos menos). Usando el dato de la consigna sobre la pérdida de energía
finalmente encontramos que s
6M R2 Ṗ P c3
= |m|
5π µ0
Ahora consideramos que el dipolo está alineado con el eje z con lo cual su campo magnético será

µ0 mẑ 3m cos θ
B (r, θ) = − 3 + r̂
4π r r3
Ítem b Reemplazando se obtiene B (R, 0) ≈ 9,68 × 108 T. Comparando este número con el valor en la superficie de la tierra
B (superficie tierra) ≈ 40 µT vemos que es bastante importante.
mp M G
Ítem c La fuerza de Lorentz es F L = qv × B ∼ 10−4 N mientras que la gravitatoria es F G = R2 ∼ 10−15 N. Como se
puede ver el movimiento estará gobernado por la fuerza de Lorentz.
La consigna está mal, debería decir “ejercicio 3 de la práctica 5”. Ahora la ecuación de equilibrio hidrostático será
1
0 = −∇p − ρ∇φ + J × B.
c
Por Maxwell sabemos que
4π 1 ∂E
∇×B = J+ .
c c ∂t
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

En los casos de interés en astrofísica siempre ocurre que

∂E
≈ 0 → En astro
∂t
por lo tanto lo que nos queda es
1
0 = −∇p − ρ∇φ + (∇ × B) × B.
4π
El último término es
((∇ × B) × B)i = εijk (∇ × B)j Bk

= εikj εj`m ∂` Bm Bk
εijk εj`m = δk` δim − δkm δi` → = Bk ∂k Bi − Bk ∂i Bk
B2
= (B · ∇) B − ∇ = (∇ × B) × B.
2
Volviendo a la ecuación del fluido nos queda
B2

1
0 = −∇p − ρ∇φ − ∇ + (B · ∇) B
8π |4π {z }
| {z }
pmagnética Tensión magnética
donde reconocemos a la presión magnética ya que entra en la ecuación de la misma forma que la presión. En cuanto a la
tensión magnética, el último término, su interpretación no es tan simple.
Ahora integramos esto en el volumen de la nube y lo que habíamos encontrado en el problema 3 de la práctica 5 era
ˆ
3 3
 ∇p · r d r = 4πR p0 − 2M cV T


Problema 3, práctica 5 → ˆ
M 2 ΘG
 ρ∇φ · r d3 r = −


R
2 ´
donde Θ = − 4πR
M2 m (r) ρ dr. Ahora lo que habría que calcular son los términos restantes que incluyen al campo magnético.
El primer término con campo magnético es
ˆ ˆ
 
2 2 2
B 3 B B
∇ ·rd r = ∇ · r − ∇ · r  d3 r
8π 8π 8π | {z }
=3
‹ ˆ
B2 B2 3
Teorema de Gauss → = r · ds − 3 d r.
8π 8π
| {z }
Energía magnética
En cuanto al otro término éste es

ˆ
1
(B · ∇) B · r d3 r = ...
4π
..
.
ˆ Bi Bj rj
ˆ
1 z }| { B2 3
= ∇ · (BB · r) d3 r − d r
4π 4π
| {z } | {z }
2EB
‹
1
= BB · r · ds − 2EB .
4π
Juntando todos los términos obtenemos
‹ B2
M 2 ΘG BB −
4πR3 = 2M cV T + + EB + 2
· r · ds
R 4π
y a veces se suele definir el tensor de Maxwell como
B2
BB −
def 2
M= → Tensor de Maxwell
Willy!
4π
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Ahora hacemos la suposición de congelamiento magnético. Esto consiste en asumir que el plasma ‚ y las líneas de campo
magnético se mueven juntos. Bajo esta hipótesis se puede mostrar que se conserva el flujo magnético B · ds. No lo pude
seguir mucho pero esto implica que BR2 se conserva. Creo que integró el flujo en el borde de la nube.
Acá me distraje, mala mía, y me perdí varias cosas. Lo que termina quedando es que la masa de Jeans es
3/2
vA
MJB = MJ 1+ 2
cs
donde
2 B2
vA = → Velocidad de Alfvén
4πρ
que es una velocidad característica de no sé qué, velocidad de Alfvén, que vamos a ver en la teórica. La conclusión es que si
hay un campo magnético vamos a precisar una masa mayor para producir el colapso gravitatorio.
Este problema es prácticamente igual a lo que hicimos hoy en la teórica. Tenemos que linealizar las ecuaciones MHD bajo
algunas hipótesis. Por un lado tenemos que es incompresible por lo tanto
ρ = ρ0 = constante. → Incompresible
Por otro lado tenemos no sé qué ideales que implica que

p p0
= γ . → Algo que es ideal
ργ ρ0
Asumimos equilibrio estático

u0 = 0 → equilibrio estático
y además un campo magnético que es
B = B0 ẑ
con B0 constante. Las perturbaciones son
(
u = u0 + u1
Perturbaciones → .
B = B0 + B1
Proponemos una solución de la forma
u1 , B 1 ∼ exp (i [k · r − ωt]) → Proponemos
y entonces las ecuaciones son
∂t ρ + ∇ · (ρu) = 0 → Continuidad



1


E + u × B = ηJ → Ohm





 c

4π 1 ∂E


∇×B = J+
 c c ∂t
∂B


∇2

= ∇ × (u × B) + η0 B


∂t



1


· ∇) u] = − ∇p
+ J × B + → Navier-Stokes

ρ [∂ u + (u ρ0
E

t
c
∂E
donde hemos cancelado la ∂t pues ∂E
u 2
∂t
∼ 1
|J | c
y además vamos a asumir Rm 1 que implica que
| ∇ × (u × B) |
Rm = 1.
| η 0 ∇2 B |
Los términos que cancelamos en la ecuación de Navier-Stokes no sé por qué es. El primero (u · ∇) u me parece que es porque
queda cuadrático en las perturbaciones. El ∇p y el último ni idea.
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Haciendo no sé qué nos queda

1 1
J ×B = (∇ × B 1 ) × B 0 .
ρ∂t u =
c 4π
Me resulta imposible seguir las clases de esta profesora. Llega un momento que me pierdo. Anotó lo siguiente
)
∂t u1 = iωu1 1
⇒ −iωρu1 = (k × B 1 ) × B 0
∇ × B 1 = ik × B 1 4π
(k × B 1 ) × B 0 = −B 0 × (k × B 1 )
Vaca meos caballo → = − [k (B 0 · B 1 ) − B 1 (B 0 · B)]
Ya está, es imposible entender lo que está haciendo. Pegó volantazo y se fue al otro lado del pizarrón. Ni idea.
1
ρωu1 = [k (B 0 · B 1 ) − B 1 (B 0 · k)] → Ecuación (3)
4π
4π

∇ × B = − J
 → Ecuación (1)
Esto viene de no sé dónde → c
 ∂B = ∇ × (u × B) → Ecuación (2)

∂t
Ahora hacemos (2)+perturbación
−iωB 1 = ∇ × (u1 × B 0 )
= (B 0 · ∇) u1 − B 0 (∇ · u1 )
Eq. continuidad
∂ρ
+ ∇ · (ρu) = 0
∂t
∂ρ
+ (u · ∇) ρ + ρ∇ · u = 0
|∂t {z }
=0 por incompresibilidad
∇·u=0
Me quedé sin batería /.
Ítem a La energía magnética es
ˆ
B2
EM = dV
8π
B2
= 2 0 πR2 L
8π
B02 R2 L
= .
4
Ítem b Hay que usar el modelo de Sweet-Parker de reconexión magnética que vimos en la teórica. Hay que estimar la
potencia
EM
hP i = .
τ
Tenemos que estimar τ . Para ello procedemos de la forma más cabeza: geometría de dos cilindros que se superponen. Tenemos
que
2R
τ∼ .
Uin
Entonces
B 2 RLUin
hP i ∼ 0 .
8
Lo del número de Lundquist está relacionado con la velocidad de Alfvén
Uin
S − /2 =
1
vA
con vA = √B0 . Finalmente nos queda
4πρ
B 3 RL
hP iSP ∼ √ 0 .
8 4πρS 1/2
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

8 GUÍA 8 - GALAXIAS
Ítem c Usando los números de la consigna se obtiene
hP iSP ∼ 1027 erg s−1
mientras que la potencia observada, según consigna, es
hP iobservada ∼ 1028 erg s−1 .
En consecuencia el modelo de Sweet-Parker no puede explicar esto. Existen modelos más complejos que dan mejor, modelos
tridimensionales (a diferencia del modelo de Sweet-Parker que es bidimensional).
8. Guía 8 - Galaxias
Esto es lo que hicimos hace un rato en la teórica. Recordemos cómo se define el LSR (local standard of rest, estándar
local de reposo). El LSR se define según la siguiente imagen
Las velocidades en el LSR se definen de modo tal que


Π = 0

En el LSR → Z = 0 .

Θ = Θ0

Vamos a asumir que el movimiento del sol alrededor del centro galáctico es circular. Consideremos además las cantidades

N
 1 X
hΠ − Π i = Πi





 N i=1
Esto es lo que medimos →
N

 1 X
hΠ − Π i = Zi



N i=1


donde i es una suma sobre las “estrellas cercanas” al sol. Por otro lado tenemos lo siguiente
(
hΠ − Π i = hΠi − Π
hZ − Z i = hZi − Z
Vamos a asumir como hipótesis que estas cantidades dan cero cuando se hace el promedio en una cantidad representativa,
por una cuestión de simetría/isotropía. Es decir
(
hΠi = 0
Cuando medimos muchas estrellas → .
hZi = 0
Entonces el promedio de las mediciones nos termina dando Π y Z . Me voy a cursar E3.
Gráficamente tenemos lo siguiente
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Theta_0
r
Theta
R0 R
Di
re
cc
ión
Centro
r
galáctico
Haciendo un poco de geometría tenemos que
R R0 R0
= = Eq. (I)
sin ` sin β cos α
R0 cos ` = r + R cos (π/2 − α)

= r + R sin α Eq. (II)
Ahora podemos proyectar la velocidad de la estrella en las componentes de la visual y su perpendicular. Es decir
vradial = Θ cos α − Θ0 cos (π/2 − `)
= Θ cos α − Θ0 sin `.
Ahora podemos usar la Eq. (I) podemos reemplazar cos α = R0/R sin ` por lo tanto

R0
vradial = Θ − Θ0 sin `
R
Ω es velocidad angular → = R0 (Ω − Ω0 ) sin `.
En cuanto a la velocidad tangencial tenemos
vtangencial = Θ sin α − Θ0 sin (π/2 − `)
= Θ sin α − Θ0 cos `

R0 cos ` − r R0 r
Usamos eq. (II) ⇒ sin α = → = Θ − Θ0 cos ` − Θ
R R R
= R0 (Ω − Ω0 ) cos ` − Ωr.
Hasta aquí no hemos hecho ninguna aproximación, esto es completamente general. Ahora vamos a asumir que estamos
mirando una estrella que está cerca del sol, es decir
R0 , R r. → Asumimos esto
Vamos ahora a aproximar la diferencia de velocidades angulares

dΩ
Ω − Ω0 ≈ (R − R0 ) + . . .
dR R0
!
Θ 1 dΘ Θ
Ω= → = − (R − R0 )
R R dR R=R0 R2 R=R0
!
1 dΘ Θ0
= − (R − R0 )
R0 dR R0 R02
!
1 dΘ Θ0
= − − r cos `.
R dR R0 R02
Reemplazando esto en la vradial obtenemos
" #
dΘ Θ0
vradial = − − r |cos {z
` sin `}
dR R0 R0
1
2 sin(2`)
!
1 d (ΩR)
= Ω0 − r sin (2`) .
2 dR Ω0
| {z }
Constante A de Oort
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Listo. Para la velocidad tangencial, haciendo algo similar, obtenemos

" # " #
1 d (ΩR) 1 d (ΩR)
vtangencial = Ω0 − r cos (2`) − Ω0 + r.
2 dR R0 2 dR R0
| {z } | {z }
Constante A de Oort Constante B de Oort
¿Para qué sirve todo esto? Tenemos un sistema

 vradial
 = A sin (2`)
r
 vtangencial = A cos (2`) − B
r
por lo tanto midiendo las velocidades de muchas estrellas cercanas al sol podemos ajustar estas funciones para estimar A y
B, y una vez que conocemos A y B las usamos para conocer la velocidad angular del sol alrededor del centro de la galaxia,
dΩ
Ω0 , y el rate con el que cambia Ω con el radio, dR R0
. Parece que recientemente se han medido estas cantidades y se obtuvo
(
A = 14,8 km s−1 kpc−1
Mediciones → .
B = −12,4 km s−1 kpc−1
Todo esto parece que está disponible en [Carroll and Ostlie, 2017, capítulo 12].
Ítem a Consideramos masa unitaria (m = 1). Para un potencial central sabemos que se conserva el momento angular L.
Además (o en consecuencia) el movimiento está contenido en el plano perpendicular a L. Sabiendo esto, podemos utilizar
coordenadas polares para describir el movimiento. Entonces
r = rr̂



ṙ = ṙr̂ + rθ̇θ̂ .

2

r̈ = r̈ − rθ̇ r̂ + 2ṙθ̇ + rθ̈ θ̂
Las ecuaciones de Newton son, entonces,
r̂) r̈ − rθ̇2 = − dφ = −∇φ

dr .
θ̂) 2ṙθ̇ + rθ̈ = 0

De la segunda ecuación sacamos

d 2
r θ̇ = 0 ⇒ r2 θ̇ = L. → Se conserva!
dt
L
De aquí sacamos que θ̇ = r2 , cantidad que podemos despejar en la ecuación para r̂. Obtenemos
d2 r L2 d2 r L2

dφ d
−r 4 =− ⇒ + φ+ 2 =0
dt2 r dr dt2 dr 2r
| {z }
φefectivo
y entonces
d2 r dφefectivo
+ = 0.
dt2 dr
Ítem b Ahora vamos a hacer pequeñas perturbaciones, para lo cual proponemos que
r = r0 + x (t)
con x (t) r0 , es decir que tenemos una órbita circular de radio r0 con una perturbación x (t). Lo que queremos encontrar
es que a primer orden vale que
d2 x
+ κ2 x = 0.
dt2
d2 r dφefectivo
Para ello vamos a hacer un Taylor de dt2 + = 0 tal que
dr
2
dφe dφe d φe
(r − r0 ) + O x2

= + 2
dr dr r0 dr r0 | {z }
x
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

l k
dφe
y usamos que dr = 0 pues se trata de un punto de equilibrio (con solución r0 ). Entonces tenemos que
r0
d2 x
2
d φe
2
+ x=0
dt dr2 r0
de donde obtenemos
d2 φe

κ2 = .
dr2 r0
Para que hayan oscilaciones es necesario que κ2 > 0.

l k
dΩ2
Ítem c Queremos encontrar que κ2 = 4Ω2 (r0 ) + r0 dr donde Ω es la frecuencia angular de rotación de la galaxia. Para
r0
ello simplemente reemplazamos en la expresión que encontramos antes
2
2 d φe
κ = ← Encontramos antes
dr2 r0
L2
2
L2

d
φe = φ + 2 → = φ +
2r dr2 2r2 r0
2
d φ 3L2
= + .
dr2 r0 r04
Antes por otro lado habíamos encontrado que

dφe
= 0 → Pues es un punto de equilibrio
dr r0
lo cual implica que

L2 Ω2 (r0 ) r04

dφ
= 3 = = Ω2 (r0 ) r0 .
dr r0 r0 r03
Si ahora derivamos con respecto a r0 obtenemos
& %
d2 φ d Ω2

= Ω (r0 ) + .
dr2 r0 dr
r0
& %
2 d Ω2
Ω (r0 ) + r0 + 3Ω2 (r0 ) →?
dr2
r0
L2 = Ω2 (r0 ) r04
κ2 = Ω20 ⇒ elipse
Ítem d
2Ω (r0 )
θ (t) = θ0 + Ω (r0 ) t − A sin (κt)
r0 κ

dθ L L L x
= 2 = 2 ≈ 1 − 2
dt r (r0 + x) r02 r0
ˆθ ˆt !
0 L 2x t0
dθ = 1− dt0
r02 r0
θ0 0
L 2L
θ − θ0 = 2 t− 2 A sin (κt + α)
r0 r0 r0 κ
Entonces tenemos que
x (t) ∼ A cos (κt + ϕ)

.
θ (t) ∼ − 2Ω (r0 ) A sin (κt + α)
r0 κ
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Recordemos que la emisión de HI es la famosa línea de los veinti un centímetros. La órbita de los planetas alrededor del
sol está dada por
mv 2 GM m 1
= ⇒ v∼√
r r2 r
Lo anterior está asumiendo que toda la masa M está adentro de la órbita. En cambio, si pensamos que el planeta está
orbitando “adentro del sol” entonces
mv 2 GM (r) m
= ⇒ v∼r
r r2
donde M (r) = 34 πρr3 es la masa encerrada por la órbita. Entonces lo que esperamos es algo así:
Velocidad de rotación
Distancia del centro

Acá estamos Acá estamos afuera
"adentro del núcleo" del núcleo
Sin embargo se observa lo de siempre de la materia oscura

Curva de rotación observada
Curva de rotación teórica

usando la masa visible
Distance from de center
De esta forma podemos estimar cuál es la distribución de materia que generaría dicha curva. Lo que se obtiene al plantear
simetría esférica es
dMr V2
= 4πr2 ρ =
dr G
y se obtiene una densidad
V2
ρ= .
4πGr2
En 1996 un señor propuso que
ρ0
ρNFW =
r r 2

a 1+ a
donde a es un parámetro de ajuste, y que es una forma funcional de densidad de materia (oscura más bariónica) que ajusta
más o menos bien en todas las galaxias.
Lente gravitacional Estamos haciendo algo medio trucho, planteamos Newton para un fotón (de masa cero, ja). Pensemos
en un objeto que pasa cerca de un objeto masivo de masa M a una distancia R:
GM
m → Newton ma =
R2
donde m es la masa del fotón (jaja). Cancelando tenemos que
GM
a=
R2
es la aceleración (perpendicular) que sufre el fotón. Entonces
∆v = a∆t
2R 

∆t =  2M G
c → = .
GM  Rc
a= 2 
Willy!
R
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

9 GUÍA 9 - COSMOLOGÍA
Parece que si uno hace relatividad general esto es lo mismo sólo que en vez de un 2 aparece un 4:
4M G
∆v = . → Si usáramos relatividad general
Rc2
Candidatos a materia oscura

1. MOND. Se propone que la gravedad es distinta, es decir F = GM r2
m
no es válida. Sin embargo tienen problemas a la
hora de explicar las colisiones de galaxias (por ejemplo la bullet galaxy).
2. Los neutrinos y los neutralinos supersimétricos.
3. MACHOS.
4. WIMPs. En este link buscan WIMPs.
9. Guía 9 - Cosmología
Tenemos que n (r) es la densidad´numérica de galaxias. φ (L) dL es la probabilidad de encontrar una galaxia con lumi-
∞
nosidad L. Entonces necesariamente 0 φ (L) dL = 1.
Ítem a Podemos decir que φ (q) dq es el número de galaxias con flujo entre q y q + dq por unidad de ángulo sólido. Entonces
ˆ
 ∞ 
φ (q) dq dΩ =  n (r) φ (L) r2 dr dq dΩ

ˆr=0
L = 4πr2 q → = n (r) φ (L (q, r)) r2 dr dq dΩ.
| {z }
φ(q)
Acá hubo un leve debate de cosas mal hechas y cosas bien hechas. Ahora parece que vamos a hacer un cambio de variable
r → L usando s
L
r=
4πq
por lo tanto dr = √ 1 dL y entonces
2 L4πq
ˆ∞ s !
L L 1
φ (q) = n φ (L) √ dL
4πq 4πq 2 L4πq
0
ˆ∞ s !
q −3/2 L
φ (L) L /2 dL
1
= 3/2
n
2 (4π) 4πq
0
ˆ √
n0 q −3/2
Universo homogéneo ⇒ n = n0 → = 3/2
Lφ (L) dL.
2 (4π)
Listo, hemos demostrado que φ (q) ∝ q −3/2 .
Ítem b La cuenta está hecha en el Bataner. El flujo total será

ˆ∞
Q = qφ (q) dq
0
√ ∞
φ ∼ q − /2 →
3
= d qc0
= ∞.
Tenemos que desarrollar un modelo de cosmología newtoniano. Se basa en la ecuación de Euler para fluidos, que está
basada en las leyes de Newton. La consigna dice “obtenga la ley de Hubble”, pero entiendo que esto lo metemos de prepo en
base a las observaciones.
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Ítem a El fluido debe satisfacer la ecuación de continuidad

∂ρ
+ ∇ · (ρv) = 0 → Continuidad
∂t
y si asumimos el principio cosmológico entonces
∇ρ = 0 → Principio cosmológico
y entonces
1 ∂ρ
∇·v =− .
ρ ∂t
Sabiendo que v = H (t) r entonces ∇ · v = 3H (t) por lo tanto
1 ∂ρ
H (t) = − .
3ρ ∂t
Ítem b Se define el factor de escala cosmológico R (t) según
Ṙ
H≡ .
R
1 ∂ρ
Sabiendo del ítem previo que H (t) = − 3ρ ∂t entonces vemos que
1 ∂ρ 3 ∂R
+ =0
ρ ∂t R ∂t
de donde encontramos que
ρR3 = constante.
Aplicando condiciones iniciales tenemos que
3
R0
ρ (R) = ρ0 .
R
Obsérvese que en ningún momento dijimos que R tenga algo que ver con el radio del universo, aunque sí se lo puede interpretar
como tal.
Véase que
v ≈ cz
= Hr
1 dR
= r
R dt
R0 − R
En relatividad general esto no es aproximado, es igual → ≈ r
tR
R0
r = ct → = −1 c
R
y entonces
R0
z= − 1 → Redshift
R
Ítem c Usamos la ecuación de Euler: Universo homogéneo
∂v ∇p7
+ (v · ∇) v = − − ∇φg
∂t ρ
(Obsérvese que dijimos que ∇p = 0 porque el universo es homogéneo. En tal caso lo mismo debería ocurrir para ∇φg , sin
embargo si no lo conservamos la cosmología newtoniana falla completamente. Parece que en el Battaner o en el Chodurry lo
dicen que esto es re trucho y que ∇φg debe conservarse “porque sí”.) y si ahora reemplazamos v = H (t) r encontramos que
dH
r + H 2 r = −∇φg .
dt
Si ahora tomamos divergencia a esta ecuación obtenemos

dH
∇· r + H r = − ∇ 2 φg
2
dt | {z }
4πGρ
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

y nos queda
dH
3 + 3H 2 + 4πGρ = 0.
dt
Ṙ
A continuación reemplazamos H = R y lo que nos queda es
3 d2 R
+ 4πGρ = 0
R dR2
o bien
4π Gρ
R̈ + = 0.
3 R
R0 3

Hace un rato vimos que ρ (R) = ρ0 R por lo tanto
4π Gρ0 R03
R̈ = − .
3 R2
Ítem d Para integrar usamos el truquito de siempre:

d dR
R̈ =
dt dt

d dR dR
=
dR dt dt
2 !
d 1 dR
= .
dR 2 dt
Gρ0 R03
Ahora usamos R̈ = − 4π
3 R2 y entonces
1 2 4πGρ0 R03
Ṙ = − k0
2 3R
o bien
8πGρ0 R03
Ṙ2 = −k
3R
donde k 0 y k son constantes de integración.
Para entender por qué sólo nos interesa que k ∈ {−1, 0, 1} recordemos que R está definido a menos de una constante
multiplicativa a través de
Ṙ
H ≡ . → Definición de R invariante ante R → R0 = αR
R
Ítem e Véase que r

8πGρ0 R03
Ṙ = ± − k.
3R
Acá no lo pude seguir... Pero hizo el típico grafiquito
R(t) k=-1 (universo abierto)
k=0 (universo plano)
k=1 (universo cerrado)
Tiempo
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!

Ítem f Se define el parámetro de densidad
def 8πGρ
Ω= . → Parámetro de densidad
3H 2
R0 3 Ṙ

Si usamos ρ = ρ0 R yH= R, entonces
8πGρ0 R03
H 2 R2 = −k
3R

8πGρ
H 2 R2 − 1 =k
3H 2
H 2 R2 (Ω − 1) = k
y acá vemos que 
 Ω0 = 1
 ⇒ k=0
Lo podemos medir → Ω0 > 0 ⇒ k=1

Ω0 < 1 ⇒ k = −1

El valor de Ω0 se puede medir por lo tanto podemos saber si estamos en un universo abierto, plano o cerrado.
No sé de dónde salió el ρcrítico que es
3H02
ρcrítico = .
8πG
Ítem g Consideramos la ecuación de balance de calor y un gas monoatómico formado por hidrógeno

3 ρ ∂T
k + v · ∇T + p∇ · v = 0.
2m ∂t
Como el universo es homogéneo entonces

∇T = 0. → Es homogéneo
Ṙ
Si reemplazamos v = Hr y H = R, y asumimos lo que dice la consigna llegamos a que
(
T ∼ R−2 if pm pr
T ∼ R−1 if pr pm
donde hay que usar

a 4
pr = T → Presión de radiación
3
4σ
donde a = c con σ la constante de Stephan-Boltzmann.
Willy!
ma
ble
Los resueltos de
pro
No hay
ALF ←click me!


Apuntes de La Práctica de Astrofisica

Cargado por

Copyright:

Formatos disponibles

Apuntes de La Práctica de Astrofisica

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Apuntes de La Práctica de Astrofisica

Cargado por

Copyright:

Formatos disponibles

Apuntes de la práctica de Astrofísica

Encontrá más resueltos de Alf en este link.

Box 1 - ¿Cómo se hacen estos apuntes?

tomar apuntes en vivo durante una clase).

2. Guía 2 - Fluidos clásicos 8

4. Guía 4 - Atmósferas estelares 16

5. Guía 5 - Medios interestelares 23

6. Guía 6 - Estructura y evolución estelar 28

7. Guía 7 - Campos magnéticos 32

ALF ←click me!

ALF ←click me!

Ítem b Vale 1,497 × 108 km.

Ítem c El dibujito es así:

ALF ←click me!

← Cuerpo negro ideal

La densidad espectral de energía en el cuerpo negro (que obtuvo Planck) es

ALF ←click me!

1.1. Cúmulos estelares

ALF ←click me!

2. Guía 2 - Fluidos clásicos

ALF ←click me!

Ítem a Planteamos Navier-Stokes en el caso estacionario y con la fuerza gravitatoria

ALF ←click me!

Ítem b La densidad de columna se define según

Tomando ahora el límite z = r − R → 0 (bajas alturas) se obtiene

Ahora podríamos buscar también la temperatura usando que T ∝ ρp .

ALF ←click me!

Reemplazando en la ecuación de balance de calor

Ítem b La variación de la temperatura es

donde hemos reemplazado la expresión de H.

Ítem c Nosotros tenemos que 

Ítem a La fuerza centrífuga es

nos queda que  

ALF ←click me!

Por otro lado la fuerza gravitatoria es

dpcsuperficie + dφcsuperficie = constante

Ítem c Vamos a definir

ALF ←click me!

νobservada = νemitida (1 − β cos θ) γ → Doppler relativista

donde θ es el siguiente ángulo

ALF ←click me!

Ítem b La contracción de la longitud es

Eso es menos que la distancia entre la tierra y la luna.

El tiempo medido por un observador en la tierra va a ser

En cuanto a la energía tenemos que q

Juntando las dos coas tenemos

ALF ←click me!

Usando ahora que E = cp para los fotones tenemos que

ALF ←click me!

4. Guía 4 - Atmósferas estelares

ALF ←click me!

Ítem b La intensidad en la superficie de la estrella es

Gráficamente esto sería así:

Ahora sólo hay que integrar

Para calcular las demás cantidades simplemente recurrimos a su definición

Ítem c Ahora asumimos que

ALF ←click me!

con lo cual ahora se definió

y se termina obteniendo la definición que usamos hoy en día

nos queda que