Astrofísica estelar · UNIR
Análisis de estrellas pulsantes: una δ Scuti observada por Kepler
Transformada de Fourier y prewhitening con Period04 sobre una curva de luz de Kepler: identificación de las frecuencias significativas, ajuste multifrecuencia, descarte de frecuencias espurias del satélite y clasificación del tipo de pulsador según Aerts et al. (2010).
Introducción
Las estrellas no son objetos estáticos: algunas muestran variaciones periódicas de brillo causadas por oscilaciones internas que modifican su densidad, su presión y su temperatura. A este grupo pertenecen las estrellas pulsantes, cuyo estudio permite inferir su estructura interna, algo imposible de observar directamente.
En este trabajo se analizan pulsaciones radiales, aquellas en las que la estrella se expande y se contrae manteniendo su simetría esférica. Estas oscilaciones se generan por el mecanismo κ, un proceso que tiene lugar en capas parcialmente ionizadas y que actúa como un motor térmico: durante la contracción la opacidad aumenta y la capa retiene energía; durante la expansión la opacidad disminuye y la energía se libera, lo que produce las variaciones de luminosidad que se observan.
La asterosismología estudia estas variaciones de brillo del mismo modo en que la sismología analiza las vibraciones terrestres. Gracias a misiones como Kepler se dispone de curvas de luz de alta precisión que permiten identificar modos de pulsación. En este informe se aplica la transformada de Fourier con Period04 para detectar las frecuencias significativas, realizar un ajuste multifrecuencia y, finalmente, proponer el tipo de pulsación siguiendo las clasificaciones de Aerts.
Ejercicio 1. Identificación de la primera frecuencia significativa
A partir de una serie temporal de datos reales observados por el satélite Kepler, se busca la primera frecuencia significativa, se determinan su amplitud, su fase y su relación señal-ruido (SNR), y se representa el espectro de amplitudes.
Cualquier señal que cambia con el tiempo, incluida la curva de luz de una estrella, puede descomponerse como la suma de muchas ondas sinusoidales (Aerts, Christensen-Dalsgaard y Kurtz, 2010), es decir, ondas del tipo:
A · sen(2π f t + φ)
donde:
- A es la amplitud (con qué intensidad vibra).
- f es la frecuencia (cuántas veces vibra por día).
- φ es la fase (en qué punto del ciclo está la onda cuando empieza la medida).
De esta manera, aunque la señal parezca compleja, irregular o ruidosa, puede analizarse como una mezcla de ondas simples. Para ello se aplica primero la transformada de Fourier a la curva de luz original, que da el espectro de amplitudes: la descomposición de la variabilidad del brillo de la estrella en sus componentes sinusoidales.
La Figura 1 muestra la variación del flujo observada por Kepler a lo largo de casi 30 días. Se aprecia una modulación periódica claramente detectable, superpuesta a ruido instrumental y a un hueco en los datos asociado a pausas de adquisición típicas de la misión.
Para identificar la primera frecuencia significativa se sigue un proceso secuencial que comienza con el cálculo del espectro de amplitudes de los datos originales.
Cálculo del espectro. Se calcula el espectro de amplitudes de la serie temporal original dentro del rango permitido por el muestreo. En estos datos, la frecuencia de Nyquist, la máxima que puede detectarse sin ambigüedad, es de 717,169 d⁻¹. A partir de este espectro se identifica el pico dominante, que corresponde a la primera frecuencia significativa.
Identificación. La frecuencia con la amplitud más alta es f₁ = 22,790 d⁻¹, que representa el modo de oscilación dominante de la estrella. En el espectro de amplitudes (Figura 3), cada pico representa una oscilación sinusoidal presente en la curva de luz y su altura corresponde a la amplitud de dicha oscilación (A = 2,534 mmag).
Desde el punto de vista conceptual, identificar la primera frecuencia implica detectar el modo oscilatorio que más energía aporta al brillo observado. Además del máximo central, se identifican otros picos de menor amplitud distribuidos a lo largo del espectro que superan el umbral de significancia utilizado (SNR ≥ 4), y que se discuten más adelante. Esto indica que, aunque la mayor parte de la variabilidad se concentra en la frecuencia dominante, la estrella presenta multiplicidad de modos de oscilación que también contribuyen a la señal observada.
Ajuste. Se ajustan simultáneamente la frecuencia, la amplitud y la fase mediante mínimos cuadrados no lineales, para obtener los valores más precisos y poder reconstruir la curva de luz.
Significancia. Antes de buscar una nueva frecuencia es necesario calcular la SNR de la frecuencia recién identificada. Este cálculo no se realiza sobre los datos originales, sino sobre los residuales, es decir, sobre la curva de luz a la que ya se le ha sustraído la frecuencia ajustada.
Siguiendo el criterio estándar para los datos fotométricos de Kepler, una frecuencia se considera significativa únicamente si SNR ≥ 4. La altísima relación señal-ruido obtenida, SNR = 376,599, indica que este modo es físicamente real y domina la variabilidad de la estrella; en una estrella pulsante, suele corresponder a un modo radial o a un modo dominante. La Tabla 1 recoge los parámetros fundamentales de esta frecuencia.
| Parámetro | Valor | Unidad |
|---|---|---|
| Frecuencia (f₁) | 22,790 | ciclos/día (d⁻¹) |
| Amplitud (A₁) | 2,533 | mmag |
| Fase (φ₁) | 0,909 | fracción de ciclo |
| Relación señal-ruido (SNR₁) | 376,599 | — |
Tabla 1. Parámetros ajustados de la primera frecuencia f₁ (modo dominante).
Estos valores resumen la información clave del proceso: la posición del pico más prominente en el dominio de la frecuencia, su amplitud fotométrica, la fase que define el desplazamiento temporal de la señal y la elevada relación señal-ruido que confirma su carácter físico. Con ellos se cierra la caracterización del modo principal, responsable de la mayor parte de la variabilidad de la curva de luz.
Ejercicio 2. Análisis de la segunda frecuencia significativa
Como la primera frecuencia supera ampliamente el criterio de significancia, se repite el proceso para la segunda frecuencia significativa, determinando su amplitud, su fase y su SNR, y se compara su espectro con el anterior.
La dinámica general del método es la misma, pero a partir de este punto la transformada de Fourier ya no se calcula sobre los datos originales, sino sobre los residuales. Esto permite revelar los picos que estaban parcialmente ocultos bajo la potencia del modo principal.
| Parámetro | Valor | Unidad |
|---|---|---|
| Frecuencia (f₂) | 21,714 | ciclos/día (d⁻¹) |
| Amplitud (A₂) | 0,250 | mmag |
| Fase (φ₂) | 0,284 | fracción de ciclo |
| Relación señal-ruido (SNR₂) | 73,108 | — |
Tabla 2. Parámetros ajustados de la segunda frecuencia.
La segunda frecuencia significativa presenta una amplitud considerablemente menor que la dominante, pero su elevada SNR confirma que se trata de un modo real de oscilación.
Al eliminar previamente el modo principal, este segundo pico se vuelve claramente visible en el espectro de los residuales, donde se convierte en la nueva señal más prominente. Se observa un pico estrecho y bien definido en 21,714 d⁻¹, con un nivel general de potencia notablemente más bajo que en el caso anterior, y aparecen picos secundarios que antes quedaban parcialmente enmascarados por la amplitud de f₁. Este comportamiento es el esperado en un proceso de prewhitening: cada sustracción sucesiva limpia el espectro y facilita la detección de modos más débiles.
Los dos espectros estudiados presentan, por tanto, características propias:
- Cambio en la escala de amplitudes. El espectro de la segunda frecuencia presenta valores de amplitud mucho menores que el de f₁.
- Mayor visibilidad relativa de los picos secundarios. Señales que estaban parcialmente ocultas bajo la potencia del modo principal emergen ahora con claridad.
- Nuevo pico dominante. La frecuencia f₂ = 21,714 d⁻¹ se convierte en el máximo del espectro, algo que no podía apreciarse antes por el predominio absoluto del primer modo.
Estas diferencias confirman que el procedimiento es correcto: al eliminar cada modo significativo, el espectro residual revela la estructura multiperiódica real de la estrella.
Ejercicio 3. Detección del conjunto completo de frecuencias
A continuación se determinan todas las demás frecuencias de la estrella siguiendo el criterio de significancia SNR ≥ 4. El procedimiento se repite de forma sistemática: aplicar una nueva transformada de Fourier sobre los residuales, ajustar la frecuencia detectada, calcular su SNR y verificar si supera el umbral. Este proceso iterativo revela progresivamente todos los modos de oscilación presentes en la estrella, y se detiene cuando SNR < 4.
El criterio SNR ≥ 4 (Breger et al., 1993; Aerts et al., 2010) es el estándar en fotometría de alta precisión para decidir si un pico del espectro de Fourier corresponde a una oscilación real de la estrella o es simplemente una fluctuación estadística del ruido.
| Iteración | Frecuencia (d⁻¹) | Amplitud (mmag) | Fase | Ruido (mmag) | SNR |
|---|---|---|---|---|---|
| F1 | 22,7904 | 2,5327 | 0,9094 | 0,0067 | 376,5990 |
| F2 | 21,7145 | 0,2503 | 0,2844 | 0,0034 | 73,1078 |
| F3 | 23,0222 | 0,0384 | 0,7307 | 0,0031 | 12,2776 |
| F4 | 20,5752 | 0,0317 | 0,9335 | 0,0018 | 17,7458 |
| F5 | 0,0534 | 0,0286 | 0,7511 | 0,0030 | 8,9366 |
| F6 | 22,0148 | 0,0302 | 0,8581 | 0,0032 | 8,5519 |
| F7 | 21,9647 | 0,0247 | 0,4722 | 0,0026 | 9,4226 |
| F8 | 42,5917 | 0,0219 | 0,5318 | 0,0022 | 9,8152 |
| F9 | 22,8137 | 0,0180 | 0,5764 | 0,0023 | 9,4848 |
| F10 | 10,3285 | 0,0139 | 0,2907 | 0,0021 | 6,5611 |
| F11 | 22,7620 | 0,0129 | 0,8319 | 0,0020 | 6,6163 |
| F12 | 11,9982 | 0,0111 | 0,0240 | 0,0020 | 5,4162 |
| F13 | 0,3319 | 0,0110 | 0,1559 | 0,0026 | 4,1530 |
| F14 | 43,5975 | 0,0107 | 0,5519 | 0,0020 | 5,2529 |
| F15 | 0,1101 | 0,0096 | 0,8080 | 0,0024 | 4,0391 |
| F16 | 42,2797 | 0,0089 | 0,2416 | 0,0018 | 4,8291 |
| F17 | 9,9015 | 0,0086 | 0,9808 | 0,0020 | 4,4897 |
| F18 | 43,3690 | 0,0087 | 0,6921 | 0,0019 | 4,6758 |
| F19 | 9,4878 | 0,0084 | 0,9800 | 0,0018 | 4,5506 |
| F20 | 679,5884 | 0,0083 | 0,6932 | 0,0020 | 4,2288 |
| F21 | 326,5817 | 0,0078 | 0,2204 | 0,0019 | 3,9897 |
Tabla 3. Frecuencias obtenidas mediante el proceso iterativo de prewhitening, con sus amplitudes, fases y relaciones señal-ruido.
Esta tabla es el resultado final del proceso de prewhitening: cada frecuencia con SNR ≥ 4 se identificó, se ajustó y se sustrajo del conjunto de datos, lo que permitió revelar modos cada vez más débiles. La última frecuencia significativa aparece en la iteración 21, que marca el final del proceso.
Ejercicio 4. Ajuste global de la curva de luz
Con todas las frecuencias validadas, se realiza un ajuste de la curva de luz con todas las frecuencias significativas a la vez, reconstruyendo la señal observada mediante la combinación simultánea de todos los modos detectados. De este modo, el modelo final incorpora únicamente las oscilaciones físicas de la estrella y excluye el ruido residual.
La superposición del modelo y los datos muestra una correspondencia casi perfecta. Esto confirma que los modos extraídos durante el prewhitening reproducen fielmente la variabilidad periódica de la estrella y que no es necesario añadir más frecuencias para describir su comportamiento pulsante.
Además, el ajuste permite extrapolar el comportamiento de la estrella en intervalos sin medidas por huecos en la adquisición de datos. En la Figura 8 se observa un ejemplo claro: aunque la serie temporal presenta un corte, el modelo reproduce de forma continua la variabilidad esperada a partir de los modos de pulsación identificados. Esto es posible porque el ajuste combina todas las frecuencias significativas, de modo que la señal resultante queda completamente definida incluso donde faltan observaciones.
Por último, las pequeñas dispersiones alrededor del ajuste se deben al nivel de ruido fotométrico propio de las observaciones, no a la presencia de oscilaciones adicionales.
Ejercicio 5. Descarte de frecuencias espurias e identificación de modos
Para seleccionar las frecuencias más representativas de la estrella y descartar las frecuencias espurias del satélite Kepler, se comparó la columna de frecuencias de la Tabla 3 con la lista de frecuencias espurias en cadencia corta publicada para la misión (SC spurious frequency summary; Christiansen et al., 2016). La comparación se realizó buscando coincidencias numéricas dentro de la anchura indicada como narrow o broad en dicha tabla. Se observa lo siguiente:
- Modos reales. Ninguna de las frecuencias en el rango de ~20 a 23 d⁻¹ (F1, F2, F3, F4, F6, F7, F9 y F11) coincide con las líneas espurias típicas del satélite (7,50; 16,98; 20,95; 25,06; 29,38; 31–32 d⁻¹, etc.), por lo que pueden considerarse modos reales de pulsación. F4 está en la vecindad de una región contaminada por el satélite, pero no coincide con la frecuencia instrumental listada, por lo que se trata como un modo de pulsación.
- Frecuencias dudosas. Las bajas frecuencias F5, F13 y F15 (≈ 0,05–0,33 d⁻¹) no aparecen en la tabla de artefactos, pero, por su cercanía a 0 d⁻¹, son candidatas a ser tendencias de largo plazo o efectos instrumentales, y se marcan como dudosas.
- Frecuencia instrumental. La frecuencia F20 = 679,59 d⁻¹, muy alta, cae en el mismo rango que la línea espuria U7 de Kepler (~679,54 d⁻¹).
| Iteración | Frecuencia (d⁻¹) | Amplitud (mmag) | Fase | Ruido (mmag) | SNR |
|---|---|---|---|---|---|
| F20 | 679,5884 | 0,0083 | 0,6932 | 0,0020 | 4,2288 |
Tabla 4. Parámetros ajustados de la frecuencia F20 y su relación señal-ruido, comparados con la frecuencia espuria U7 del satélite Kepler (Christiansen et al., 2016).
Puede concluirse, por tanto, que F20 = 679,588 d⁻¹ es una frecuencia instrumental, asociada al muestreo del satélite, y se descarta como modo pulsante de la estrella. F21 ya había sido descartada por no cumplir el criterio SNR ≥ 4.
Para determinar qué tipo de estrella se observa, se comparan las frecuencias representativas con las tablas de Aerts et al. (2010) y se discute de forma razonada qué tipos de pulsación son posibles. En asterosismología, cuando una frecuencia se expresa en ciclos por día (d⁻¹), el periodo se obtiene simplemente como:
P = 1 / f
donde P es el periodo (en días) y f la frecuencia (en ciclos por día). Tomando f₁ = 22,790 d⁻¹, se obtiene P = 0,0438 días; multiplicando por 24, P = 1,0526 h.
Extrapolando este cálculo a todas las frecuencias no descartadas, pueden formarse grupos según su periodo y calcular para cada uno el periodo medio de pulsación en horas, para proponer finalmente con qué tipo de estrella presenta la mayor compatibilidad según la tabla de Aerts, Christensen-Dalsgaard y Kurtz (2010):
| Grupo | Frecuencia (d⁻¹) | Iteraciones | Periodo medio (h) | Compatibilidad (Aerts) |
|---|---|---|---|---|
| I | 20 – 23 | F1, F2, F3, F4, F6, F7, F9, F11 | 1,0821 | δ Sct |
| II | 42 – 43 | F8, F14, F16, F18 | 0,5588 | δ Sct, ¿sdBV de modo g, GW Vir? |
| III | 9 – 12 | F10, F12, F17, F19 | 2,3194 | δ Sct, ¿β Cep, T Tauri pulsantes, Herbig Ae/Be pulsantes? |
Tabla 5. Agrupación de frecuencias representativas y comparación de los periodos medios con los rangos de pulsación típicos según Aerts et al. (2010).
Los modos detectados presentan periodos de entre aproximadamente 0,5 y 3 horas, lo que cae de lleno en el rango característico de las estrellas δ Scuti, cuyos periodos típicos van de 15 minutos a 8 horas (Aerts et al., 2010). Este comportamiento corresponde a modos de presión (p-modes), coherentes con las oscilaciones rápidas y de alta frecuencia que muestra el espectro.
Además, las amplitudes obtenidas (entre ~2,53 mmag y ~0,008 mmag) también son compatibles con las δ Scuti, que presentan variaciones luminosas muy inferiores a 0,3 mag (300 mmag). Por lo tanto, tanto los periodos como las amplitudes sitúan a la estrella dentro del régimen típico de esta clase.
Al comparar los periodos observados con las demás familias de la tabla de Aerts se advierten algunas coincidencias parciales. Sin embargo, estas otras familias pueden descartarse: aunque algunas tienen rangos de periodo que se superponen parcialmente con los valores observados, poseen propiedades físicas incompatibles con la señal analizada. Por ejemplo, las β Cep tienen periodos típicos de entre unas 2 y 8 h (Stankov y Handler, 2005; Aerts et al., 2010), lo que solo coincide con el grupo III; los grupos I y II, más rápidos, quedan por debajo de su límite inferior. β Cep no puede explicar los periodos más cortos, y solo una δ Scuti reproduce de forma simultánea los tres grupos.
En resumen, los periodos, las amplitudes, el número de modos detectados y la estructura del espectro coinciden de forma consistente con una estrella δ Scuti. El resto de candidatos aparecen como coincidencias parciales en uno u otro grupo, pero no pueden reproducir todo el patrón de pulsaciones observado. Esta clasificación es la que mejor se ajusta al conjunto de propiedades derivadas del análisis asterosísmico.
Conclusión
En este trabajo se identificaron 20 frecuencias significativas (SNR ≥ 4) en la curva de luz observada por Kepler, de las que 16 corresponden a modos de pulsación fiables (F20 es instrumental y F5, F13 y F15 quedan como dudosas), tras aplicar un análisis iterativo de prewhitening y un ajuste multifrecuencia preciso. La comparación con el catálogo de artefactos instrumentales permitió descartar las líneas espurias y aislar únicamente los modos reales de la estrella. A partir de la conversión de frecuencias en periodos y de su contraste con las tablas de clasificación de Aerts et al. (2010), se determinó que todos los valores obtenidos se encuentran dentro del rango típico de pulsación de las estrellas δ Scuti. Las amplitudes ultrabajas medidas, junto con la riqueza modal y la estructura final del espectro, refuerzan esta clasificación. Por tanto, el análisis asterosísmico confirma que la estrella estudiada es un claro ejemplo de pulsadora δ Scuti de baja amplitud observada con la precisión fotométrica del satélite Kepler.
Referencias
- Aerts, C., Christensen-Dalsgaard, J., & Kurtz, D. W. (2010). Asteroseismology. Springer.
- Bracewell, R. N. (1986). The Fourier Transform and Its Applications. McGraw-Hill.
- Breger, M., Stich, J., Garrido, R., et al. (1993). Nonradial pulsation of the delta Scuti star BU Cancri in the Praesepe cluster. Astronomy & Astrophysics, 271, 482–486.
- Christiansen, J. L., Van Cleve, J. E., Jenkins, J. M., Caldwell, D. A., Barclay, T., Bryson, S. T., et al. (2016). Kepler Data Characteristics Handbook (KSCI-19040-005). NASA Ames Research Center; Space Telescope Science Institute.
- Lenz, P., & Breger, M. (2005). Period04 User Guide: An analysis tool for astronomical time series. Communications in Asteroseismology, 146, 53–136.
- NASA Kepler Mission (2013). Archivo de datos de fotometría. Material docente de la asignatura Estructura y Evolución Estelar e Interestelar, Universidad Internacional de La Rioja.
- Stankov, A., & Handler, G. (2005). Catalog of Galactic β Cephei stars. The Astrophysical Journal Supplement Series, 158(2), 193–216.