Atmosferas Estelares
Atmosferas Estelares
Atmosferas Estelares
Introducción
Las estrellas son gases y conocemos las ecuaciones de estado de los mis-
mos, de modo que podemos obtener los valores de la temperatura, la presión
y la densidad en cualquier punto de la misma y bajo cualquier circunstancia
a través de modelos fı́sicos.
¿ Por qué son tan importantes los modelos estelares? Un modelo nos
describe como se va a comportar una estrella bajo unas determinadas condi-
ciones; entre otros parámeros podemos obtener la temperatura y la presión,
de modo que nos dará la localización de la misma en el diagrama HR. Con-
forme varı́an estos parámetros también lo hacen el radio, la temperatura
efectiva y la luminosidad, de modo que podemos conocer la trayectoria de la
estrella en este diagrama. Y al contrario, conociendo el lugar de una estrella
en el diagrama HR podemos obtener su estructura interna.
1
se comporta como ideal, con una ecuación de estado
k
Pg = ρT
µmH
siendo µ el peso molecular medio definido como µ = (2 X+ 31 Y + Z )
4 2
Las ecuaciones necesarias para estudiar la estructura estelar son las sigu-
ientes:
• Equilibrio hidrostático
dP m(r)G
=− ρ(r)
dr r2
Esta ecuación debe cumplirse para todas las estrellas y nos indica que
en cualquier punto de la misma la presión de radiación debe ser igual
al peso de la columna de material que hay por encima. Es simplemente
una generalización del equilibrio que debe de existir entre la presión y
la gravedad.
• Conservación de la masa
dm(r)
= 4πr 2 ρ(r)
dr
2
Estas expresiones se aplican a los interiores o a las atmósferas de forma
independiente para poder obtener un resultado sobre los mismos.
Interiores estelares
y es independiente de la frecuencia.
Modelos lineales
La suposición más simple que se hace para poder resolver las ecuaciones
3
es que la densidad varı́a linealmente con el radio
r
ρ = ρc (1 − )
R
A partir de aquı́ debemos obtener los valores de la densidad, la presión y la
temperatura en función de la masa, el radio y la luminosidad.
La ecuación de equilibrio hidrostático queda
dP m(r)G r
= 2
ρc (1 − )
dr r R
Utilizando la ecuación de la conservación de la masa se puede obtener el valor
de m(r) y de ρc . Y con ellas, integrando, la expresión de P
2 7 r3 1 r4
P = −Gπρ2c [ r 2 − + ] + Pc
3 9 R 4 R2
2
con Pc = 54 GM
πR4
. La expresión de la temperatura incluso en un modelo tan
sencillo como este, requiere de cálculos tediosos, ası́ como de las expresiones
para ǫ(r) y κ.
Polı́tropos
4
Para n=3 tenemos una estrella con masa independiente de la densidad
central. En este caso, la estrella es soportada por la presión de la radiación
y la presión del gas y se llama ” modelo estándard de Eddington”.
Atmósferas estelares
La atmósfera de una estrella está formada por las capas más externas
y que son lo único observable desde el exterior. Además, es aquı́ donde se
forman las lı́neas espectrales. A pesar de que hay numerosos factores que
complican el cálculo (rotación, campo magnético, etc) se pueden estudiar
utilizando modelos simples de estrellas no magnetizadas y sin rotación. Con-
sideraremos primero modelos de atmósferas completamente radiativas, ya
que esta es la forma más usual de intercambiar energı́a en ellas. En este
caso, el problema es más simple que para los interiores ya que no hay fuentes
de energı́a y se conoce la estructura de las capas inferiores.
Los principales parámetros en el estudio de las atmósferas estelares son
la temperatura efectiva (Tef f ), la aceleración de la gravedad en la superficie
g, y la composición quı́mica. Para un conjunto de parámetros la estructura
es calculable. Sin embargo, es mucho más fácil hacer el proceso inverso. Es
decir, a partir de las condiciones observadas de color, tamaño e intensidad
de las lı́neas espectrales intentar obtener el conjunto (Tef f , g, Z) que mejor
se ajuste. Las lı́neas espectrales contienen mucha información, tanto en la
intensidad como en la forma, sobre la presión, la composición y la presencia
de turbulencia o campos magnéticos en la zona que se originaron.
Una de las complicaciones en el estudio de las atmósferas está en el hecho
de que claramente no cumplen las condiciones de LTE, y además, la radiación
no es plackiana cerca de la superficie.
5
Modelos de atmósferas
Esta ecuación, utilizándose después que las otras dos anteriores, contiene
toda la información necesaria. Si se utiliza una profundidad óptica tal que
τ = cosΘ se obtiene la llamada aproximación de Eddington-Barbier en donde
Iν (Θ) = Bν (T (τ = cosΘ)). Es decir, que la radiación que vemos se origina a
una profundidad óptica de 1 a lo largo de la lı́nea de visión.
Modelo plano-paralelo
6
Supongamos un modelo de capas plano-paralelas de profundidad h. Esta
aproximación es posible porque tanto el radio total como el radio de curvatura
de la estrella son mucho mayores que la anchura de la atmósfera, ya que para
la mayorı́a de los casos h = 0.001R. Además, también consideramos LTE
4
y equilibrio radiativo (Frad = σTef f = cte). El sistema de referencia que
utilizamos es el (h, k), donde h es la profundidad fı́sica, aumentando hacia
dentro y k es la normal, cuyo coseno sobre h se denomina µ.
La ecuación de equilibrio hidrostático se puede integrar suponiendo que
la temperatura es constante y tendremos
ghµmH
P = Po exp[ ]
kT
Al factor de la exponencial que nos darı́a la distancia sobre la cual la
presión cambia un factor e se le llama escala de alturas, H.
La ecuación de transporte radiativo se puede escribir dSν = − dz
µ
siendo
Sν la función fuente que aquı́ se define como
Sν = (1 − Aν )Bν (T ) + Aν Jν
sca
con Aν = κκν y Jν la intensidad especı́fica media.
Utilizando la ecuación de equilibrio térmico en atmósferas se puede obtener
c dPr 4
− = σTef f
κ dh
donde Pr = 0∞ Pν dν es la presión de radiación total. Utilizando la definición
R
4 2
q = CσTef f Pr (S) =
3
Para escribir estas ecuaciones hemos supuesto implicitamente que hay
LTE y que la función fuente es isotrópica. La temperatura se obtiene como
la temperatura efectiva en τ = 2/3.
Atmósfera gris
7
las longitudes de onda, y la ecuación de transporte se puede integrar en
frecuencias obteniéndose
1 dF
− =J −S
4 dτ
De modo que como el flujo es constante tendremos que la función fuente es
igual a la intensidad media, J = S, en primera aproximación.
Adoptamos ahora de nuevo LTE, de modo que la ecuación de transferencia
integrada nos da la intensidad en función de µ (= cosΘ) y de τ ,
dI 1
Z
µ =I −J =I − Idµ
dτ 2
Esta ecuación es conocida como de Mine-Schwarzchild. Se puede resolver
formalmente, I(µ, τ ), aunque no tiene solución analı́tica. Una aproximación
se puede obtener como
3
J(τ ) = F [τ + q(y)]
4π
siendo q(y) la llamada función de Hobbs y vale 0.377 para τ = 0 y 0.71
para τ = ∞, ambos valores muy parecidos a 2/3. Para resolverla se ha
utilizado la aproximación de cierre de Eddington, que es Prad ≈ 4πJ
3c
. Con
esta aproximación se puede explicar el oscurecimiento de los bordes, ya que
cuando nos acercamos a los mismos el ángulo (µ) aumenta y la intensidad
que se observa es menor.
Transiciones atómicas
8
los mismos que está en el estado de excitación apropiado para absorber el
fotón correspondiente. Este número se puede obtener a partir de la llamada
ecuación de Saha, siempre que supongamos que el gas está en equilibrio
térmico. Para obtenerla partimos de la expresión general de Boltzman
Ns gs
= exp(−Ξs /kT )
No go
siendo No el número de átomos en el nivel fundamental, NS el de átomos en
el nivel excitado, gs y go son los pesos estadı́sticos de esos estados atómicos
y ΞS es la energı́a de excitación definida como la diferencia de energı́a entre
el nivel fundamental y el excitado Ξs = Ξion (1 − n1 ). Mediante esta expresión
podemos conocer la excitación del átomo, es decir, las poblaciones de los
niveles energéticos de cada especie. Esta expresión se puede referir al número
total de átomos neutros existentes (N), definido como Σn=0 Ns , en la forma
Ns gs exp(−Ξs /kT )
=
N Σn=0 gs exp(−Ξs /kT )
donde el denominador es la llamada suma de estados o función de partición.
Esta expresión nos da la excitación térmica de los átomos cuantizados.
Ahora debemos de considerar también la excitación de los átomos libres.
Para ello sea N1 el número de átomos una vez ionizados y No el de átomos
neutros. El problema principal que se plantea es calcular el peso estadı́stico
del átomo ionizado en el estado fundamental más el electrón libre. Este se
puede calcular integrando la probabilidad de un estado con momento p y
energı́a cinética m v 2 /2 en el espacio de momentos. Es decir,
Z inf ty p2 /2 m
exp(− ) 4 π p2 dp = (2 π m k T )3/2 h−3 Ne−1
0 kT
Entonces, g = g1,0 2 (2 π m k T )3/2 . De modo que la ecuación de Saha que
nos da la ionización se puede expresar como
9
a desionizarse por lo cual la fracción de átomos activos disminuye y por
tanto dejan de observarse en las atmósferas estelares. Además, la fracción de
ionización aumenta con la profundidad, como cabia esperar.
Estas expresiones se usan cuando el proceso dominante en la población
de los niveles de energı́a son las colisiones.
Perfiles de lı́neas
10
siendo f (v) la fracción de átomos moviendose a la velocidad v. El ensan-
chamiento Doppler es más importante en la longitud de onda central y de-
crece muy rapidamente (de forma exponencial) conforme nos alejamos de
ella. Este ensanchamiento tiene a su vez dos componentes: el movimiento
térmico de los átomos y el movimiento turbulento de la atmósfera (rotación,
movimientos convectivos, salidas de masa, etc). El primero sigue la dis-
tribución de Boltzman mientras que el segundo es mucho más complejo. En
general, la turbulencia es más importante en las atmósferas de las estrellas
gigantes y supergigantes. Cuando la escala de los movimientos turbulentos
es menor que la escala de la atmósfera, entonces se puede suponer que siguen
2 2 2
una distribución
q gaussiana y la velocidad Doppler total serı́a vD = vth + vtur
con vth = 2kT m
. En este caso, generalmente el ensanchamiento térmico dom-
ina sobre todas las demás y se puede tomar como el único responsable de la
foma del perfil. La opacidad total serı́a
ao 1 ∆ν 2
a(ν) = √ exp −( )
π ∆νD ∆νD
con ∆ν = ν − νo .
Otro fenómeno que influye en los perfiles es el de las colisiones entre iones
y átomos neutros. Esto tiene como resultado un amortiguamiento de la lı́nea
pero con un valor mayor de γ.
Por último, para algunos átomos se produce efecto Stark. Este es de-
bido a los campos eléctricos originados por iones lentos distribuidos al azar.
En este caso, las alas de las lı́neas presentarán un ensanchamiento propor-
cional a 1/∆λ5/2 . También la profundidad óptica produce un mayor o menor
ensanchamiento de las alas de las lı́neas.
Cuando se toman en cuentan los ensanchamientos debidos a todos estos
factores se puede escribir que
πe2 1 1 α exp(−(v − u)2 )
Z
a(νo ) = f√ [ du]
me c π ∆νD π u2 − α 2
∆ν ∆ν ′ γ
con v = ∆ν D
, u = ∆ν D
y α = 4π∆ν D
. A la cantidad entre paréntesis se le
llama función de Voigt y al perfil total se le llama perfil de Voigt.
En el centro de la lı́nea la opacidad está siempre dominada por el ensan-
chamiento Doppler, independientemente de la cantidad de amortiguamiento
que exista y este sólo influye en la forma de las lı́neas. Estas tienen normal-
mente unas alas muy extendidas.
Se define la cantidad llamada anchura equivalente de la lı́nea como
Z
EW = (Fc − Fl )/Fc dλ
11
Estrellas reales
Las estrellas reales son mucho más complejas que lo expuesto anterior-
mente. En los modelos, tanto de atmósferas como de interiores, las capas
están claramente delimitadas y cuando se produce un intercambio de ma-
terial se hace de forma instantánea, sobre todo en la zona convectiva. Sin
embargo, en modelos más reales se impone semiconvección en las fronteras
entre capas, es decir, hay un intercambio lento en una región finita. Además,
este intercambio de masa es difusivo, con lo que le lleva un tiempo no cero
en producirse.
Por último, también se impone que las células convectivas penetren dentro
de la zona radiativa estable (overshooting). La distancia de penetración es
el gran tema de discusión, y se han propuesto valores desde 0 hasta 2 Hp ,
siendo Hp la escala de alturas de la presión.
La semiconvección es importante, sobre todo en estrellas poco masivas,
hasta 5 M⊙ , mientras que el overshooting afecta a todas las estrellas que
tengan un núcleo convectivo.
La consecuencia principal de la existencia de semiconvección y overshoot-
ing es el alargamiento de la vida de la estrella en la secuencia principal en casi
un factor 2. Esto es debido a que la región accesible a la combustión nuclear
se hace mayor y por tanto es mayor también el tiempo durante el cual se
pueden producir reacciones nucleares. Sin embargo, la luminosidad también
aumenta y esto hace que la evolución posterior sea mucho más rápida.
12