Ejercicios entregables de Astrofísica Miguel Ángel Aloy
Tercer curso, Grado en Física, UV Departamento de Astronomía y Astrofísica, UV
Tema 8.– Estructura estelar. Ejercicio entregable
8.1 Modelos estelares politrópicos
Figura 1: Soluciones de la ecuación de Lane-Emden para n = 0, 1, 2, 3, 4, 5.
Introdución: los modelos politrópicos son modelos en simetría esférica y equilibrio hidrostático
en los que la presión y la densidad están relacionados mediante una función polítropo:
p = Kρ1+1/n ,
donde K y n son constantes reales positivas, y n es el índice del polítropo. La estructura de
estos modelos se describe mediante las ecuaciones de equilibrio hidrostático y conservación de
la masa
dp G Mr ρ
=− , (1)
dr r2
dMr
= 4πr2 ρ, (2)
dr
que pueden combinarse en una única ecuación diferencial, de segundo orden, dada por
1 d r2 dp
= −4πGρ. (3)
r2 dr ρ dr
Para su resolución, se introducen las variables adimensionales θ, ξ
ρ = ρc θ n , p = pc θn+1 , r = αξ
donde ρc y pc son la densidad y la presión centrales, respectivamente, y α es una constante con
dimensiones de longitud definida por
(1−n)/n
K(n + 1)ρc
α2 = .
4πG
1
Ejercicios entregables de Astrofísica Miguel Ángel Aloy
Tercer curso, Grado en Física, UV Departamento de Astronomía y Astrofísica, UV
Sustituyendo las variables adimensionales en la ecuación (3), obtenemos la ecuación en variables
adimensionales
1 d 2 dθ
ξ = −θn , (4)
ξ 2 dξ dξ
ecuación de Lane-Emden de índice n para estrellas politrópicas. Las condiciones de contorno
para esta ecuación son:
dθ
θ = 1, = 0, en ξ = 0.
dξ
La solución tiene sentido físico mientras θ ≥ 0, de manera que la superficie de la estrella se
obtiene para ξ = ξ1 , donde θ = 0 y, de acuerdo con la definición, la presión y la densidad
también se anulan.
Una vez resuelta la ecuación de Lane-Emden y obtenida θn (ξ), θnn da directamente el perfil
de densidad en unidades de la densidad central, ρn (ξ)/ρc , y θnn+1 , el perfil de presiones en
unidades de la presión central, pn (ξ)/pc . Considerando una ecuación de estado de gas ideal,
para la que p ∝ ρ T , θn representa el perfil de temperatura en unidades de la temperatura
central, Tn (ξ)/Tc .
La ecuación de Lane-Emden tiene solución analítica para n = 0, 1, 5. El caso n = 0 se co-
rresponde con un fluido incompresible homogéneo. El caso n = 5 se corresponde con un modelo
de radio infinito. Se puede probar que todos los modelos con n > 5 tienen radio infinito. Hay
dos casos especialmente interesantes para los que la ecuación de Lane-Emden debe resolverse
numéricamente. Se trata de los casos n = 1.5 y n = 3.
El caso n = 1.5 describe objetos tan distintos como enanas blancas de masa pequeña so-
portadas por la presión de degeneración de los electrones en el límite no relativista, o estrellas
convectivas de secuencia principal. El caso n = 3, por su parte, describe tanto enanas blancas
de masa grande o estrellas de neutrones soportadas por la presión de degeneración (de los elec-
trones, en el caso de las enanas blancas; de los neutrones en el caso de las estrellas de neutrones)
en el límite relativista. El caso n = 3 describe también estrellas radiativas de secuencia principal
y puede utilizarse para construir un modelo bastante aproximado del Sol.
Objetivo: el objetivo de la práctica es familiarizarse con las soluciones de la ecuación de Lane-
Emden y resolverla numéricamente al menos para los casos n = 1.5 y n = 3.
Metodología: puedes resolver la ecuación usando cualquier estrategia. Puedes utilizar progra-
mas de cálculo simbólico o programar tú el algoritmo de integración. En las líneas que siguen,
describimos la estrategia para resolver la ecuación de Lane-Emden con el método de Euler.
La ecuación de Lane-Emden es una ecuación diferencial de segundo orden que podemos
reescribir como:
d2 θn 2 dθn
+ + θnn = 0
dξ 2 ξ dξ
(θn hace referencia a la función θ(ξ) solución de la ecuación de Lane-Emden de índice n).
Definiendo
dθn
φn (ξ) =
dξ
obtenemos
dφn 2
= −θnn − φn .
dξ ξ
2
Ejercicios entregables de Astrofísica Miguel Ángel Aloy
Tercer curso, Grado en Física, UV Departamento de Astronomía y Astrofísica, UV
Las dos ecuaciones acopladas anteriores pueden ahora resolverse para obtener θn (ξ).
En el método de Euler definimos
dθn
θn (ξi+1 ) = θn (ξi ) + h (ξi ) = θn (ξi ) + h φn (ξi ),
dξ
dφn n 2
φn (ξi+1 ) = φn (ξi ) + h (ξi ) = φn (ξi ) − h θn (ξi ) + φn (ξi ) .
dξ ξi
h = ξi+1 − ξi es el paso de integración.
Inicia la integración en ξ = 0 con las condiciones de contorno θn (0) = 1, φn (0) = 0. Usa un
paso de integración h = 0.0001. Continúa la integración hasta que θn (ξ) < 0 o ξ > 10.
Resultados:
1. Representa las funciones θn , θnn , θnn+1 en función del radio adimensional ξ para n =
0, 1, 1.5, 3, 5.
2. Deduce la expresión para la masa incluida hasta un cierto radio en unidades de la masa
total: 2
Mn (ξ) ξ φn (ξ)
= .
Mn (ξ1 ) ξ1 φn (ξ1 )
3. Representa la masa incluida en unidades de la masa total en función de ξ para los índices
del polítropo del punto 1.
8.2 Un modelo para el Sol (Opcional)
Objetivo: nos planteamos la construcción de un modelo sencillo para representar la estructura
interna del Sol basado en un polítropo de índice n = 3.
Metodología: considera el Sol gobernado por una ecuación de estado de gas ideal totalmente
ionizado
k
p= ρT,
µmu
1
≈ 2X + 3Y /4 + Z/2
µ
(donde X, Y y Z son las fracciones de masa promedio de hidrógeno, helio y metales, respecti-
vamente, del Sol), y densidad y temperatura centrales
ρc = 1.58 × 102 g cm−3 , Tc = 1.57 × 107 K.
Bajo estas hipótesis, se puede obtener la constante K del polítropo y, a partir de ella, la
constante de proporcionalidad α que permitirá obtener el radio del Sol según nuestro modelo
politrópico,
R = αξ1 ,
3
Ejercicios entregables de Astrofísica Miguel Ángel Aloy
Tercer curso, Grado en Física, UV Departamento de Astronomía y Astrofísica, UV
y su masa
M = −4πα3 ρc ξ12 φ3 (ξ1 ).
Para el cálculo de la luminosidad, consideraremos la producción de energía nuclear por
unidad de masa correspondiente a las cadenas PP y que viene dada por
−2/3 −1/3
PP = 2.41 × 106 ρX 2 T6 exp (−33.80T6 ) erg s−1 g−1 ,
donde T6 representa la temperatura en unidades de 106 K y la densidad, ρ, se expresa en g/cm3 .
Ahora, la luminosidad del modelo se obtendrá integrando numéricamente la ecuación para
el transporte radiativo
Z R Z ξ1
2 6 3 2 2 −2/3 16/3 −1/3 −1/3
L = 4πr ρPP dr = (9.64 × 10 )πα ρc X Tc, 6 ξ 2 θ3 exp (−33.80Tc, 6 θ3 )dξ.
0 0
Resultados:
1. Justifica que el peso molecular medio de un gas compuesto por una mezcla de hidrógeno,
helio y elementos más pesados (metales) todos totalmente ionizados, verifica aproxima-
damente la expresión
1
≈ 2X + 3Y /4 + Z/2
µ
2. Obtén los valores del radio y la masa de nuestro modelo, R , M , y compáralos con los
valores correctos del Sol.
3. Representa los valores físicos de la temperatura en función de la distancia en unidades
de r/R , T (r/R ). Haz lo mismo con la densidad, ρ(r/R ), y la distribución de masa,
M(r/R ).
4. Calcula el valor de la constante que acompaña a la integral de la luminosidad,
−2/3
L0 = (9.64 × 106 )πα3 ρ2c X 2 Tc, 6 ,
y comprueba que vale
L0 = 3.20 × 1040 erg s−1 .
5. Integra numéricamente la integral de la luminosidad y representa su valor en función de
la distancia en unidades de r/R , L(r/R ). Finalmente, determina la luminosidad de
nuestro modelo, L , y compárala con el valor correcto para el Sol.
