ASTROFÍSICA DE ALTAS ENERGÍAS
Actividad Guiada:
Variabilidad y curvas de luz en rayos X
Máster en Astronomía y Astrofísica
Profesor VIU: Pere Blay
e-mail: [Link]@[Link]
Máster en Astronomía y Astrofísica/VIU
Índice
1.- Objetivo 3
2.- Introducción 3
3.- Análisis temporal 4
4.- Period04 5
4.1.- Uso de Period04 5
5.- Actividad guiada 7
6.- Tarea para el Portafolios 8
PARTE OBLIGATORIA 8
PARTE VOLUNTARIA 8
Memoria 9
Página 2
Máster en Astronomía y Astrofísica/VIU
1.- Objetivo
Aprender analizar una curva de luz de una fuente de rayos X.
Obtener información relevante sobre el comportamiento temporal del objeto a estudiar.
Familiarizarse con herramientas de análisis de series temporales (Period04)
2.- Introducción
Las fuentes de rayos X pueden presentar variabilidad temporal en un rango grande
de escalas temporales, desde milisegundos, en los periodos de pulsación más cortos,
hasta centenares de días, en los periodos orbitales más largos. A partir de una lista de
eventos, generada con los datos almacenados por nuestro detector, se puede estudiar
dicha variación temporal para poder investigar las propiedades físicas del objeto
observado.
Hay muchas herramientas de análisis temporal. Cada investigador acaba
familiarizándose con alguna en particular o incluso desarrollando las suyas propias.
Desde programas en IDL (Interactive Data Language), Matlab o R, hasta pequeños
códigos en C, o FORTRAN. En esta actividad utilizaremos Period04, porque ya os es
conocido.
3.- Análisis temporal
Los instrumentos de altas energías nos proporcionan una lista de eventos. Es
decir, una lista de detecciones en nuestro plano detector que incluye información sobre la
posición del detector o el lugar donde se ha producido la detección, en qué momento se
ha producido, y también sobre cuál era la energía del fotón incidente. Si nos quedamos
con el número de fotones para cada energía, tendremos un espectro, como los que vimos
en la actividad anterior. Sin embargo, si la información que escogemos es el número de
fotones detectados en un intervalo temporal determinado, estamos construyendo una
curva de luz.
Página 3
Máster en Astronomía y Astrofísica/VIU
Hay muchos fenómenos que pueden producir cambios en el brillo en rayos X que
medimos de una fuente deterinada. Entre las muchas fuentes posibles de variabilidad
podemos destacar:
Movimientos orbitales: Producirán variabilidad si las condiciones de producción
de rayos X cambian a lo largo del movimiento orbital, en caso de sistemas binarios
(órbitas elípticas, distribución irregular de material a o largo de la órbita, etc)
Rotación de estrellas de neutrones: Dan lugar a las conocidas pulsaciones en
rayos X. Cuando cae material hacia una estrella de neutrones no lo hace
directamente, más bien es atrapado por las líneas de campo magnético de la
estrella de neutrones y, de esta forma, conducida hacia uno de los polos
magnéticos de ésta. De esta manera los choques de material contra la superficie
de la estrella de neutrones se producen únicamente en sus polos magnéticos. Si el
eje magnético no está alineado con el de rotación, como resultado observaremos
de forma intermitente la fuente (pulsos) como si fuese la luz de un faro.
Cambio en las condiciones físicas del material emisor: Como ocurre en el
material en expansión expulsado en una explosión de supernova, se va enfriando
con el tiempo, disminuyendo la cantidad de energía emitida en e rango de los rayos
X.
La variabilidad puede ser periódica, quasi-periódica o aperiódica. Es periódica
cuando la curva de luz presenta un patrón que se repite continuamente, puede ser debido
a movimientos de rotación u orbitales, por ejemplo. En estos casos un periodograma
(herramienta que nos da información sobre qué frecuencias o periodos dominan en una
curva de luz), mostrará picos bien definidos y altos en las frecuencias correspondientes al
periodo de dichos movimientos. En el caso de variabilidad quasi-periódica, el
periodograma nos mostrará picos muy anchos y no tan significativos (no tan altos). Por
útimo si la variabilidad es periódica (como en una explosión de supernova, por ejemplo) el
periodograma sólo nos mostrará ruido. Vemos que el periodograma es la herramienta
esencial para estudiar qué tipo de variabilidad presenta nuestra fuente. El aspecto de la
curva de luz de un objeto, además, nos estará indicando que tipo de fenómeno o
propiedad del objeto es el causante de dicha variabilidad. Por eso es importante un
análisis concienzudo de las curvas de luz de fuentes de rayos X.
Página 4
Máster en Astronomía y Astrofísica/VIU
4.- Period04
Aunque hay muchas librerías y código para analisis temporal, utilizaremos
Period04, ya que os es conocido. Dentro de las aplicaciones específicas de la astrofísica
de altas energías, podemos destacar las herramientas efsearch y powerspec, del paquete
de las FTOOLS llamado Xronos. Son herramientas muy básicas y con muchas
limitaciones. Por esta razón generalmente son utilizadas como una herramienta de
diagnóstico o primer vistazo, para luego pasar a utilizar herramientas más especializadas.
El programa está disponible en el link:
[Link]
4.1.- Uso de Period04
Podéis instalar el programa en cualquier PC, recomendamos el uso de la máquina
virtual Linux que ya conocéis. Se ejecuta con el comando 'period04', o haciendo click en el
icono correspondiente si la habéis instalado en windows o mac.
La primera acción a llevar a cabo es cargar la curva de luz, y desmarcar la casilla
que indica que el eje Y representa magnitudes. Para ello buscamos el botón 'Import Time
String', que nos abrirá un diálogo para poder cargar la cruva de luz. Nos mostrará un
cuadro donde podemos seleccionar a qué corresponde cada columna. La primera
indicaremos que es “Time” (estará en días para todas las curvas de luz menos las que
tienen 'ppuls' en el nombre, que estará en segundos), la segunda columna la marcamos
como 'Observed' y el resto de columnas como 'Ignore’. Seguidamente buscamos, debajo
de la información de la curva de luz, la casilla a desmarcar titulada como ”Time string is in
magnitudes” y que viene marcada por defecto.
Página 5
Máster en Astronomía y Astrofísica/VIU
Desde la ventana principal es importante ver la curva de luz para ver si a simple
vista nos da alguna pista del comportamiento de la fuente. Podemos ampliar zonas de la
cruva de luz para investigar los detalles más finos.
Pasamos a la pestaña de Fourier. Allí nos aseguramos que la frecuencia máxima
para analizar sea de 0.1 para las curvas de luz en las que buscaremos periodos orbitales
(ficheros terminados en d1_n.dat), y la frecuencia de Nyquist, calculada por el propio
Period04, para las curvas de luz con las que buscaremos periodos de pulsación (con un
'ppuls_n.dat' en el nombre del archivo).
Después de realizar la descomposición de Fourier de la curva de luz, podemos ver
a gráfica con el resultado (el periodograma) y determinar si alguna frecuencia domina o
no. Es posible que en algunos casos sea necesario dividir la curva de luz en trozos más
pequeños para trabajar con más comodidad, para ello, desde la pestaña ’Time String’
vamos al menú 'Specials' de la parte superior, y seleccionamos 'Subdivide Time String'.
Tendremos que indicar con que criterios dividimos en trozos la curva de luz para,
posteriormente, hacer un análisis de Fourier en cada trozo. Para ello, en la ventana
emergente seleccionamos a opción 'Subdivide by fixed time intervals'. Un intervalo de
unos 2000 días suele ser interesante, pero la elección dependerá de la curva de luz con la
que estemos trabajando.
Los detalles para el uso de Period04 comando se encuentran en la ayuda disponible
online:
[Link]
Página 6
Máster en Astronomía y Astrofísica/VIU
5.- Actividad guiada
La práctica consiste en los siguientes pasos:
1. Leer y tener a mano el manual de Period04.
2. Os daremos un conjunto de curvas de luz para trabajar con ellas.
3. Hay que ver qué aspecto general tiene la curva de luz, si aparentemente tiene
variaciones periódicas o no, si tiene épocas de emisión más intensa o si más o menos
siempre se mantiene alrededor de un valor medio, etc. Tomar nota de la posibilidad de
que sea interesante más adelante dividir la curva de luz en trozos y calcular el tamaño
temporal de dichos trozos.
4. Posteriormente se pasa a realizar el análisis de Fourier. Hay que representar el
periodograma e inspeccionar su aspecto. Si alguna frecuencia destaca la apuntamos y
calculamos el periodo correspondiente (la inversa de la frecuencia).
Para determinar que un pico sea válido para una buena medida de periodicidad, ha de
cumplir que su relación señal ruido sea mayor que 4. Para calcular la señal ruido de un
pico, hemos de acudir al menú “Special” y seleccionar la opción “Calculate Noise”.
Tenemos que seleccionar o escribir la frecuencia para la cual queremos calcular la
relación señal ruido y el “Box Size”, que será el número de puntos alrededor del pico que
consideramos para la medida. Éste puede ser de ~6 para picos estrechos y ~10 para
picos anchos. La relación señal ruido es el valor “snr” que nos mostrará Period04
después de pulsar “Calculate” en ésta ventana. Recordemos que ha de ser mayor que 4
para que el pico sea válido y nos permita calcular periodicidades.
Los períodos los estaremos midiendo en días y las frecuencias en 1/día. Éste último valor,
la frecuencia, es el que nos da el periodograma de Period04. Para las curvas de luz con
un 'ppuls' en el nombre, las frecuencias están en Hercios (Hz o 1/segundo) y los periodos
calculados estarán en segundos. En resumen, para una frecuencia dada “f”, el periodo
correspondiente se calcula de esta forma:
P = 1f
Se estima el error de medida del periodo obteniendo la significancia de la detección del
pico a la frecuencia dada con la opción de “Calculate LS uncertainties” dentro de la
lengüeta “Goodness of fit”, en el apartado de ajustes (“Fit”), una vez tenemos el error de la
frecuancia (segunda columna), el error del periodo se obtiene de la siguiente manera:
Ef
EP =
f2
Página 7
Máster en Astronomía y Astrofísica/VIU
Donde Ef es el error en la frecuencia. El valor final que tenemos que dar es el periodo y su
error de la forma P ± E P .
5. Opcionalmente: usando la base de datos de SIMBAD1, obtener información adicional
sobre la fuente para determinar a que corresponde la variabilidad que hemos encontrado.
6. Si fuese preciso, dividir la curva de luz en trozos y realizar el análisis de la curva de luz
en cada trozo, comentando los resultados obtenidos.
6.- Tarea para el Portafolios
La tarea a realizar en esta práctica será estudiar la curva de luz de un conjunto de fuentes
con Period04. Sugerimos comentar el aspecto general de la curva de luz. Antes de
comenzar a calcular el periodograma. Se toma nota del aspecto general del periodograma
y si alguna frecuencia destaca se calcula el periodo correspondiente (os saldrá en días
para los ficheros acabados en “d1_n.dat” y en segundos para los ficheros acabados en
“[Link]”). Se valorará el hacer comentarios sobre la posible naturaleza de dicho
periodo. Valorar brevemente si es interesante o no dividir la curva de luz en trozos y en
caso positivo, realizar el análisis de Fourier en cada trozo.
PARTE OBLIGATORIA
Como tarea obligatoria hay que realizar el análisis descrito con anterioridad para unos de
los ficheros terminados en “[Link]” (tiempo en segundos) y para al menos 5 de los
ficheros acabados en “d1_n.dat” (tiempo en días).
PARTE VOLUNTARIA
Realizar el análisis para todos los ficheros y/o intentar contestar alguna de estas
preguntas:
1.- ¿Cómo crees que afecta el movimiento orbital de la estrella de neutrones al periodo de
pulsación que detectamos?
2.- ¿Crees que puede haber factores geométricos (orientación relativa a nosotros de la
órbita) que puedan resultar en la no detección de periodicidades?
3.- En algunas fuentes se habla de periodos ‘super-orbitales’. Busca referencias y cuenta
brevemente a qué se refieren.
1
[Link]
Página 8
Máster en Astronomía y Astrofísica/VIU
Memoria
La tarea se ha de entregar en formato PDF. Los criterios de evaluación se pueden consultar en la
sección Guía Didáctica del CampusVIU. Es importante mantener una determinada estructura en el
documento. Debe contener los siguientes apartados
1.-Introducción
Donde se detalle qué se quiere conseguir y cómo (qué pasos se van a seguir)
2.- Desarrollo de la actividad
Donde se detalle el trabajo realizado con cada espectro.
3.- Discusión de los resultados/conclusiones
Donde se detalle el resultado obtenido para cada espectro y se comente cualquier cosa de
interés relacionada con el proceso, dificultades, etc. Aquí se incluirán los puntos 2, 3, y 4 si
procede, de la parte voluntaria.
4.- Bibliografía
Es importante, sobre todo para facilitar la evaluación de la actividad, una redacción clara y
coherente.
Página 9
