Letter 1
Caracterización de Supernovas con IRAF
A NDRÉS A LEJANDRO R OJAS G ARCÍA1,* AND A NGEL U RIEL M OLINA S ERRATOS2,*
* Departamento de física, Centro Universitario de Ciencias Exactas e Ingenierías - CUCEI, Universidad de Guadalajara
1 andres.rojas4426@[Link]
2 angel.molina4448@[Link]
20 de mayo de 2024
1. INTRODUCCIÓN Tipo Ib
Las supernovas, fenómenos astronómicos de extrema Origen: Se forman a partir del colapso del núcleo
importancia, representan el explosivo final de la vida de de estrellas masivas (mayores de 8 masas solares)
una estrella masiva. Estos eventos cósmicos son crucia- que han perdido su envoltura de hidrógeno, proba-
les para comprender la evolución estelar, así como para blemente debido a fuertes vientos estelares o a la
investigar la composición del universo y sus procesos interacción con una estrella compañera [1].
fundamentales.
Espectro:
Se clasifican principalmente en dos grandes categorías,
basadas en sus espectros y en sus mecanismos de explo- 1. Características notables: Ausencia de líneas de
sión: las supernovas de tipo I y las supernovas de tipo II hidrógeno.
[1]. Cada una de estas categorías se divide en subtipos
2. Presencia de líneas de helio (He I) en el espectro.
adicionales según características específicas observadas
en sus espectros. 3. También pueden mostrar líneas de otros elemen-
tos como oxígeno, magnesio, y calcio.
Supernovas de Tipo I
Tipo Ic
Las supernovas de tipo I se caracterizan por la ausencia
de líneas de hidrógeno en sus espectros. Se dividen en Origen: Similar a las de tipo Ib, pero las estrellas
tres subtipos principales: progenitoras también han perdido su capa de helio,
además de la de hidrógeno [1].
Tipo Ia Espectro:
Origen: Se cree que resultan de la explosión termo-
nuclear de una enana blanca en un sistema binario, 1. Características notables: Ausencia de líneas de
cuando la enana blanca acumula suficiente masa de hidrógeno y helio.
su compañera para alcanzar el límite de Chandra- 2. Dominadas por líneas de elementos pesados co-
sekhar (1̃.4 masas solares) [2]. mo el oxígeno, calcio y magnesio.
3. La falta de helio en el espectro distingue a las
Espectro:
supernovas Ic de las Ib.
1. Características notables: Ausencia de líneas de
hidrógeno. Supernovas de Tipo II
Las supernovas de tipo II se caracterizan por la presen-
2. Presentan líneas de absorción del silicio ioniza- cia de líneas de hidrógeno en sus espectros. Resultan del
do (Si II) a 615 nm (más prominente cerca del colapso del núcleo de estrellas masivas que han retenido
máximo brillo). su envoltura de hidrógeno [1].
3. A lo largo del tiempo, los espectros de las su-
pernovas Ia muestran líneas de elementos inter- Tipo II-P (Platea)
medios como el oxígeno, magnesio y calcio, que Origen: Colapso del núcleo de estrellas masivas que
evolucionan hacia líneas de hierro y cobalto a han retenido su envoltura de hidrógeno.
medida que la supernova envejece y el material
expulsado se vuelve más transparente. Espectro:
Letter 2
1. Características notables: Presencia de líneas de Cada tipo y subtipo de supernova muestra un espectro
hidrógeno. distintivo que refleja su origen y el proceso de explosión,
2. Fase de meseta en su curva de luz, donde la proporcionando a los astrónomos información sobre la
luminosidad se mantiene casi constante durante evolución estelar y los mecanismos de las explosiones de
varias semanas debido a la recombinación del supernovas[6].
hidrógeno en la envoltura expandida. [3]
Tipo II-L (Lineal) Efecto Doppler
Origen: Similar a las II-P, pero pueden tener envoltu- El efecto Doppler se basa en el principio de que cuando
ras menos masivas. una fuente de luz se mueve hacia ti, las ondas de luz se
comprimen, lo que hace que la longitud de onda aparente
Espectro:
se vuelva más corta, lo que se llama corrimiento al azul.
1. Características notables: Presencia de líneas de Por otro lado, cuando una fuente se aleja de ti, las ondas
hidrógeno. de luz se estiran, lo que resulta en una longitud de onda
aparente más larga, conocida como corrimiento al rojo [7].
2. Curva de luz que muestra un declive lineal en
lugar de una meseta. [4] En el contexto de las supernovas, este fenómeno se usa
para determinar la velocidad a la que las estrellas están
Tipo IIb viajando, ya sea hacia nosotros o alejándose de nosotros,
así como para estudiar la expansión del universo. Cuan-
Origen: Estrellas que han perdido la mayor parte,
do se observan las líneas espectrales de una supernova
pero no toda, de su envoltura de hidrógeno[1]
en el espectro electromagnético, si estas líneas aparecen
Espectro: desplazadas hacia el rojo, se infiere que la supernova y su
entorno están alejándose de nosotros. Si están desplazadas
1. Características notables: Inicialmente presentan hacia el azul, se deduce que se están acercando[8].
líneas de hidrógeno. El estudio del corrimiento al rojo de las supernovas ha
2. Con el tiempo, el espectro evoluciona para pare- sido fundamental en la cosmología moderna. Por ejemplo,
cerse más al de una supernova de tipo Ib, mos- la observación de supernovas tipo Ia con corrimientos
trando líneas de helio a medida que el material al rojo mayores de lo esperado llevó a la conclusión de
externo se dispersa. que el universo no solo está en expansión, sino que esta
expansión se está acelerando. Esto llevó al descubrimiento
Tipo IIn (Interacción con el medio circunestelar) de la energía oscura, una misteriosa forma de energía que
Origen: Colapso del núcleo de estrellas masivas ro- parece estar impulsando esta expansión acelerada [8].
deadas por densas envolturas de material expulsado
en episodios previos.
Espectro: En este contexto, el presente proyecto se centra en la
observación y análisis de una supernova cercana, con
1. Características notables: Líneas estrechas de emi- el objetivo principal de obtener su espectro utilizando
sión de hidrógeno superpuestas sobre líneas an- datos proporcionados por el Observatorio Astronómico
chas de emisión y absorción. Nacional en San Pedro Mártir, México (Figura 1).
2. Las líneas estrechas indican interacción con ma- La obtención del espectro de una supernova cercana
terial circunestelar denso y lentamente en expan- ofrece una valiosa oportunidad para estudiar sus caracte-
sión. [4] rísticas distintivas y líneas espectrales. A través del aná-
lisis con IRAF, buscamos determinar la composición quí-
Caracterización Espectral de Supernovas mica de los materiales expulsados durante la explosión,
La caracterización espectral de las supernovas es cru- así como clasificar el tipo de supernova en cuestión. Este
cial para entender la física de las explosiones, la naturaleza proceso involucra la identificación y catalogación de las
de las estrellas progenitoras y la interacción con su en- líneas espectrales características, cuya presencia y intensi-
torno. Los espectros permiten identificar la composición dad proporcionan información crucial sobre la naturaleza
química del material expulsado, la velocidad de expan- física y química del evento.
sión de la supernova, y la energía liberada en la explosión. La colaboración entre observatorios y equipos de in-
Estas observaciones se realizan a través de telescopios vestigación es fundamental para abordar estos desafíos
equipados con espectrógrafos, que descomponen la luz astronómicos. En este proyecto, aprovecharemos los re-
de la supernova en sus componentes espectrales, revelan- cursos del Observatorio de San Pedro Mártir para llevar a
do detalles que no son visibles en observaciones de luz cabo observaciones precisas y detalladas de la supernova
integrada [5]. seleccionada.
Letter 3
todo el mundo.
Este estudio pretende arrojar luz sobre los misterios
que rodean a las supernovas, proporcionando una visión
más profunda de su naturaleza y contribuyendo al avance
del campo de la astrofísica. La investigación aquí propues-
ta representa un paso significativo hacia la comprensión
de los procesos fundamentales que dan forma al cosmos
y su diversidad cósmica. Al profundizar en el estudio de
las supernovas, no solo estamos explorando los límites
Fig. 1. Observatorio San Pedro Mártir [9]. de nuestra comprensión del cosmos, sino que también
estamos contribuyendo al avance del campo de la astro-
física. Cada supernova es un laboratorio cósmico único
IRAF, ”Image Reduction and Analysis Faci-
que nos ofrece una ventana privilegiada hacia los proce-
lity”(Facilidad de Reducción y Análisis de Imágenes),
sos que dan forma al universo. Este proyecto representa
es un poderoso software utilizado en astronomía para
un esfuerzo dedicado para desentrañar los secretos de
el procesamiento, análisis y visualización de datos
estos eventos explosivos y nos acerca un paso más hacia
astronómicos. Desarrollado por el National Optical
la comprensión completa de los fenómenos astrofísicos
Astronomy Observatory (NOAO) en la década de 1980,
más poderosos y fascinantes del universo.
IRAF se ha convertido en un estándar en la comunidad
astronómica debido a su flexibilidad, robustez y amplia
gama de herramientas especializadas. 2. METODOLOGÍA
No obstante, en la actualidad, IRAF ya no es amplia- A. Adquisición de Datos
mente utilizado debido a la falta de actualizaciones en los
El profesor nos proporcionó los datos obtenidos por el
últimos años. Esta falta de mantenimiento ha llevado a
observatorio San Pedro Martir [1] en las noches del 24 y
una disminución en su relevancia, ya que no se adapta
25 de febrero del año 2014; en las cuales se encuentran los
a los avances tecnológicos y a las necesidades actuales
archivos con los cuales vamos a trabajar.
de los astrónomos. Las nuevas generaciones de software
astronómico ofrecen interfaces más modernas, mejor inte-
B. Procesamiento de Datos
gración con tecnologías recientes y un soporte continuo,
factores cruciales para el trabajo eficiente y preciso en la Para comenzar con este proyecto el primer paso fue ins-
investigación astronómica. A pesar de esto, IRAF sigue talar el sistema operativo Linux, utilizando la distribución
siendo recurrido por muchos astrónomos gracias a su le- Ubuntu 22.04.3 LTS en nuestra computadora de trabajo,
gado de confiabilidad y su extensa documentación, lo que ya que se hizo uso del software IRAF para procesar los
lo hace una opción familiar y aún útil para ciertas tareas datos proporcionados.
específicas.
C. Extracción de Espectros
Con una arquitectura modular y extensible, IRAF per-
mite a los astrónomos realizar una variedad de tareas, El primer paso fue identificar el archivo que contenía la
desde la reducción de datos brutos de observación hasta imagen de la supernova en la noche del 24 de febrero del
la elaboración de imágenes finales y el análisis de datos 2014 para después, mediante el comando apextract, extraer
espectroscópicos. Entre sus características principales se el espectro de la misma. Para esto, una vez dentro de la
incluyen algoritmos avanzados para el procesamiento de interfaz del apall:
imágenes, la manipulación de datos espectrales, la calibra-
ción fotométrica y la creación de mosaicos de imágenes.
El software IRAF se ejecuta en sistemas Linux y se
distribuye junto con un extenso conjunto de paquetes y
utilidades, lo que lo hace adecuado para su uso en una
amplia variedad de plataformas informáticas. Además,
IRAF ofrece una interfaz de línea de comandos altamente
flexible que permite a los usuarios personalizar sus flujos
de trabajo según sus necesidades específicas.
A lo largo de los años, IRAF ha seguido evolucionando
con actualizaciones regulares y la contribución de la
comunidad astronómica, manteniendo su relevancia
y utilidad en la era moderna de la astronomía digital.
Su robusta funcionalidad y su capacidad para manejar
grandes volúmenes de datos astronómicos lo convierten
en una herramienta indispensable para investigadores en Fig. 2. Interfaz apall.
Letter 4
Lo primero fue seleccionar el archivo con el que se
va a trabajar colocanto el nombre de este en el apartado
llamado “input”, después en la sección APERTURE PA-
RAMETERS se colocan los parametros de estracción; en
el apartado “Lower aperture limit relative to bright”se co-
locó el valor de -5 y en el apartado “Upper aperture limit
relative to bright”se colocó el valor de 5, esto quiere decir
que el espectro se extraera en un rango de [-5,5] pixeles
relativos al pico del brillo. En la sección BACKGROUND
PARAMETERS se selecciona una región del fondo de la
imagen para que el software compara dicha región con la
extraída de la supernova y pueda eliminar algo de “rui-
do”causado por la galaxia, para esto se selecciono el rango
[-120:-90,90,120] en el apartado “b sampl”.
Una vez colocados los parametros salimos de la interfaz
del apextract y ejecutamos el comando apall para que se
extraiga el espectro, una vez que el software nos arroja Fig. 4. Espectro de la lámpara en pixeles.
la gráfica presionamos ”q” para avanzar y ENTER a to-
do lo que sigue hasta que llegamos a una parte donde
debemos hacer un ajuste de curva manual, eliminando en el apartado “input ”se coloca el nombre del archivo de
tdodos los puntos que estén relativamente lejos de la línea la lampara y en el apartado “referen ”se coloca el nombre
de tendencia central que se nos muestra ahí mismo, para del archivo creado anteriormente con el espectro de la
eliminar los puntos se coloca el cursor sobre uno de estos supernova, después se modifican los siguientes apartados
y se presiona la tecla ”d”, repetimos hasta terminar con
todos los puntos que se encuentren alejados de la linea de interac = no
tendencia y nuevamente presionamos ”q” para avanzar y recente = no
ENTER a lo que nos indique.
Una vez terminado esto, se crea un nuevo archivo el cuál trace = no
contiene la gráfica del espectro pero no está en las uni-
dades de medición correctas, para esto, por cada imagen EXTRACTION PARAMETERS
estelar tomada por el observatorio se hace una toma a una
bacgro = none
lampara cuya composición es conocida para posterior-
mente calibrar el espectro obtenido de la toma de datos Después de modificar esto, salimos de la interfaz y ejecuta-
anterior, lo cual fue lo que hicimos como siguiente paso. mos con el comando apall y se crea un nuevo archivo con
el espectro de la lámpara el cual tampoco está calibrado
pero ese es el siguiente paso a realizar
D. Calibración del espectro de la lámpara
Para cambiar los valores en los que se encuentra la
longitud de onda de la lámpara se utiliza el comando
identify, una vez dentro de la interfaz de este comando:
Fig. 3. Espectro de la supernova con longitud de onda
en píxeles.
Ahora es mediante la misma interfaz de apextract, se
extrajo el espectro de la lámpara correspondiente a la toma
de datos de la supernova, dentro de la interfaz: Fig. 5. Interfaz identify.
Letter 5
Se coloca el nombre del archivo del espectro de la lám-
para en el apartado “images ”mientras que en el apartados
“coordi ”se coloca el nombre del archivo de una lista, tam-
bién proporcionada por el profesor, la cual contiene todas
las longitudes de onda correspondientes a los valores ya
establecidos en el espectro actual de la lámpara: Cu, He,
Ar, y Ne. Al salir de la interfaz se ejecuta el comando
“identify ”.
E. Calibración del espectro de la supernova
Lo siguiente fue utilizar el archivo del espectro de la
lámpara ya calibrado para calibrar el espectro de la super-
nova, para esto se utiliza el comando hedit, dentro de la
interfaz se modifican los siguientes parametros:
Fig. 7. Espectro de la supernova con longitud de onda
en angstroms (Å)
F. Determinación de la composición de la supernova
Una vez obtenido el espectro con los valores correctos
de longitud de onda, por último se toman los valores don-
de se encuentran los valles del espectro y se comparan
con los valores de los posibles elementos que se podrían
encontrar en ella para así obtener una buena aproxima-
ción y determinar cuáles son los aquellos elementos que
coinciden con estos y por lo tanto seguramente compo-
nen la supernova que se observó. Al tener estos datos es
momento en el cual podemos clasificar a la supernova y
Fig. 6. Interfaz hedit. saber de cuál se trata.
image = “Nombre del archivo del espectro de la su- 3. RESULTADOS
pernova ”
A. Determinación de Composición
fields = REFSPEC1 Resultados del análisis espectroscópico de SN2014J:
value = “Nombre del archivo del espectro de la lám- Línea espectral Longitud de onda [Å] Línea espectral Longitud de onda [Å]
para ”
1 4476 8 6039
2 4744 9 6138
add = yes
3 4818 10 6283
4 5074 11 6477
show = yes
5 5326 12 6489
6 5718 13 6873
Una vez modificados estos valores, salimos de la interfaz y
7 5895 14 7083
se ejecuta el comando con “hedit ”y después presionamos
ENTER en todo lo siguiente, lo siguiente fue utilizar el
Cuadro 1. Lineas de absorción detectadas en nuestro
comando dispcor, dentro de su interfaz se modifican los
espectro.
siguientes parametros:
input = “Nombre del espectro de la supernova ” Tras comparar las líneas espectrales observadas en
SN2014J con las líneas espectrales ya conocidas (Véase
output = “Nuevo nombre para el archivo que se va a imagen 8) de elementos y compuestos, se han obtenido
crear ” los siguientes resultados:
confirm = yes Silicio (Si): Se observaron líneas espectrales caracte-
rísticas de la presencia de silicio en el espectro de
Se ejecuta con “:go ”y ENTER SN2014J.
Letter 6
Sodio (Na): Aunque se observó una línea espectral
pronunciada de sodio, se asume que esta es el resul-
tado de la materia gaseosa interestelar presente en la
galaxia M82. Es decir, no es propia de la supernova
SN2014J.
Hierro (Fe): Se identificaron líneas espectrales aso-
ciadas con el hierro, indicando su presencia en el
espectro de la supernova.
Estos resultados sugieren la presencia predominante de
hierro y silicio con breves apariciones de sodio en la com-
posición química de SN2014J. Esto proporciona informa-
ción valiosa sobre la evolución estelar y las condiciones Fig. 8. Líneas espectrales conocidas de SN2014J. Recupe-
físicas en el entorno de la explosión supernova. rado de [Link]
[11]
B. Clasificación de Supernovas
Basándonos en los elementos detectados en el análisis
espectroscópico de SN2014J, podemos determinar su cla- y 66 días (22 de marzo) después de su explosión. Esos
sificación como una supernova de tipo Ia. Las supernovas perfiles han sido indicados con líneas blancas y amarillas
de tipo Ia son conocidas por la presencia de líneas espec- respectivamente. Podemos estimar la velocidad de cierre
trales de compuestos ferrosos y silicatos entre otros. Estas radial (llámese expansión) de la explosión de superno-
supernovas se producen en sistemas binarios donde una va a partir de las dos longitudes de onda observadas de
enana blanca acumula material de una estrella compañe- silicio ionizado (Si II) a 6355 Å y silicio neutro (Si I) a
ra hasta alcanzar una masa crítica, desencadenando una 5948 Å. A continuación la tabla [2] indica las velocidades
explosión termonuclear. de expansión calculadas a partir de los desplazamientos
La detección de silicio en particular es característica hacia el azul de estas dos longitudes de onda de líneas de
de las supernovas de tipo Ia, ya que se forma durante la absorción observados.
explosión termonuclear cuando el material de la enana Los perfiles tienen una diferencia en las distribuciones
blanca se fusiona. La presencia de sodio y hierro también Doppler de los átomos de silicio ionizados en 6355 Å. En
es consistente con este tipo de supernova. el perfil del 22 de marzo, la absorción en 6140 Å parece
Por lo tanto, la clasificación de SN2014J como una su- más débil que la del 22 de febrero, mientras que se puede
pernova de tipo Ia está respaldada por los resultados del observar la absorción alrededor de 6355 Å, longitud de
análisis espectroscópico, lo que sugiere que se trata de onda original sin el desplazamiento Doppler [12]. Parece
una explosión termonuclear en una enana blanca. que la absorción de Si II se ha dividido en dos partes. Se
considera que la materia interna con menor velocidad de
C. Validación de Resultados expansión se ha vuelto observable debido a la dilución de
Un estudio en particular que respalde la clasificación la supernova.
de SN2014J como una supernova de tipo Ia basada en la La velocidad de recesión (VR ) de la galaxia debido a
detección de elementos como silicio, sodio y hierro es el la expansión cósmica (es decir, la velocidad de la galaxia
trabajo de Li, W. et al.,”Observations of Type Ia supernova M82), se puede estimar a partir de los desplazamientos
2014J for nearly 900 days and constraints on its progenitor hacia el rojo de las líneas (resultados estimados en la tabla
system”. publicado en The Astrophysical Journal letters [2]). El valor de VR según su corrimiento al rojo, es de
en 2019. Este estudio proporciona una extensa cobertura unos 400 km/s, la velocidad de expansión del SN2014J,
espectroscópica de SN2014J durante un período de casi incluida la corrección, se puede calcular entre 10.000 y
900 días después de su descubrimiento [10]. 12.000 km/s.
En dicho artículo analizan detalladamente el espectro
de SN2014J y examinan la evolución temporal de las ca- Línea de absorción Si (I) Si (II)
racterísticas espectrales, incluidas las líneas de elementos
como silicio, sodio y hierro. Su estudio confirma la presen- L. de onda observada (λi ) 6140 Å 5720 Å
cia de estos elementos en el espectro de la supernova, lo L. de onda teórica (λ0 ) 6355 Å 5948 Å
que respalda su clasificación como una supernova de tipo Corrimiento al azul (z) 3,38 × 10−2 3,83 × 10−2
Ia.
Vel. de acercamiento (V) 9,98 × 103 km/s 1,13 × 104 km/s
D. Efecto Doppler, expansión y velocidad de acercamien- Vel. de expansión (V ′ ) 1,04 × 104 km/s 1,17 × 104 km/s
to
La figura [8] muestra los espectrogramas y perfiles es- Cuadro 2. Valores calculados para corrimiento al azul,
pectrales observados de SN2014J 38 días (22 de febrero) velocidad de acercamiento y velocidad de expansión
Letter 7
Los valores calculados fueron obtenidos mediante las 9. N. Triedo, «San Pedro Mártir y los confines del
siguientes ecuaciones: universo,» [Link]
[Link] (2019). Re-
∆λ λ − λi cuperado el dd/mm/aaaa.
z= = 0 (1) 10. W. Li, X. Wang, M. Hu, et al., The Astrophys. J. 882, 30 (2019).
λ0 λ0 11. «M82 - Supernova SN 2014J,» Recuperado el 17/04/2024. https://
[Link]/object_e/m82_2014j.htm.
(1 + z )2 − 1
V = c· (2) 12. National Institute of Standards and Technology (NIST), «Atomic spectra
(1 + z )2 + 1 database,» (s.f.). Recuperado el 10/05/2024.
V ′ = V + VR = V + 400 km/s (3)
en donde
z es el corrimiento al azul
λi es la longitud de onda observada
λ0 es la longitud de onda teórica
V es la velocidad de acercamiento
V ′ es la velocidad de expansión
VR es la velocidad de la galaxia M82
4. CONCLUSIONES
El análisis realizado utilizando IRAF ha proporcionado
una comprensión más profunda de la composición y cla-
sificación de la supernova SN2014J. Nuestros resultados
confirman que SN2014J es una supernova de tipo Ia, lo
que sugiere un origen en un sistema binario que involucra
una enana blanca acumulando material de una estrella
compañera.
Además, hemos identificado la presencia de elementos
clave en el espectro de la supernova, incluyendo hierro
(Fe), silicio (Si) y sodio (Na), aunque este último no sea
propio de la supernova. Estos hallazgos respaldan la teo-
ría de que las supernovas de tipo Ia son el resultado de
la detonación termonuclear de una enana blanca, lo que
libera grandes cantidades de energía y produce elementos
más pesados a través de la nucleosí[Link] otro lado,
Logramos utilizar nuestros datos recabados para calcular
más detalles acerca de este gran fenómeno astronómico
como la velocidad de expansión y el corrimiento al azul,
mismos que concuerdan con datos de otras investigacio-
nes.
REFERENCIAS
1. M. Turatto, S. Benetti, and E. Cappellaro, Handbook of Supernovae
(Springer International Publishing, 2017).
2. D. Branch and A. M. Khokhlov, Phys. Reports 256, 53 (2001).
3. S. E. Woosley, A. Heger, and T. A. Weaver, Rev. Mod. Phys. 74, 1015
(2002).
4. C. Gall, J. Hjorth, and A. C. Andersen, Astron. Astrophys. 528, A13
(2012).
5. A. V. Filippenko, Annu. Rev. Astron. Astrophys. 35, 309 (1997).
6. S. Blondin, T. Matheson, R. P. Kirshner, et al., The Astron. J. 143, 126
(2012).
7. A. G. Riess, A. V. Filippenko, P. Challis, et al., The Astron. J. 116, 1009
(1998).
8. E. Hubble, Proc. National Acad. Sci. 15, 168 (1929).
