A Spitzer Space Telescope Exploration Science Program to Search for Y Dwarf Variability

Resumen Técnico: Un Programa de Ciencia de Exploración del Telescopio Espacial Spitzer para Buscar la Variabilidad en Enanas Y

Planteamiento del Probleza
Las atmósferas de las enanas ultrafrías (tipos espectrales M7 y posteriores) están dominadas por la opacidad molecular y de condensados en lugar de la opacidad atómica. Si bien los modelos sugieren que las estructuras de nubes heterogéneas moduladas por la rotación impulsan la variabilidad fotométrica en las enanas L y T, las enanas Y más frías han sido difíciles de estudiar debido a su debilidad intrínseca en las longitudes de onda óptica e infrarrojo cercano. Estudios previos en el infrarrojo medio para enanas Y se vieron limitados por tamaños de muestra pequeños o grandes incertidumbres fotométricas. Además, aunque la variabilidad en las enanas L y T se ha vinculado a la evolución de las nubes, el mecanismo que impulsa la variabilidad en las enanas Y más frías —donde los flujos de calor interno son bajos y los tiempos de disipación son largos— permanece menos restringido. El objetivo principal de este trabajo es buscar y caracterizar la variabilidad fotométrica en el infrarrojo medio en una muestra exhaustiva de enanas Y para determinar las fracciones de variabilidad, los periodos de rotación y las tendencias de amplitud, y para probar la hipótesis de que las variaciones en la estructura de las nubes son el principal motor de esta variabilidad a lo largo de la secuencia de enanas marrones.

Metodología
Los autores llevaron a cabo un Programa de Ciencia de Exploración del Telescopio Espacial Spitzer (Programa #90015) dirigido a 14 enanas Y confirmadas espectroscópicamente en el momento de la propuesta (2012). Para crear una muestra completa de todas las enanas Y observadas para variabilidad con Spitzer, se incluyeron dos objetos adicionales, WD 0806−661B y WISE J085510.83−071442.5, elevando el total de la muestra a 16 objetos.

• Observaciones: Los objetivos fueron observados en "modo de observación fija" (staring mode) utilizando la Cámara de Matriz Infrarroja (IRAC) tanto en [3.6] como en [4.5] $\mu$m. Cada objeto fue monitoreado continuamente durante 24 horas (12 horas por banda) para cubrir potenciales periodos de rotación. Las observaciones se repitieron en una segunda época, separada por 122 a 464 días, para investigar la estabilidad a largo plazo y los cambios en la morfología de la curva de luz.
• Reducción de Datos: Los fotogramas de Datos Calibrados Básicos (BCD) se convirtieron a electrones totales. El centrado se realizó utilizando funciones gaussianas 2D, y la fotometría de apertura se llevó a cabo con un radio de 2 píxeles. Los niveles de fondo se determinaron a partir de regiones libres de fuentes en lugar de utilizar anillos para evitar la contaminación. Se eliminaron los valores atípicos que excedían 50 veces la desviación absoluta de la mediana (MAD).
• Análisis de Variabilidad: Se emplearon dos métodos estadísticos distintos para identificar fuentes variables:
  1. Comparación de Modelos Bayesianos: Utilizando el muestreo de Cadenas de Markov Monte Carlo (MCMC), los autores compararon tres modelos: un flujo constante, una curva de seno simple y una curva de seno doble. Se utilizó el Criterio de Información Bayesiano (BIC) para seleccionar el mejor modelo de ajuste. Una fuente se consideró variable si el $\Delta$BIC entre el modelo constante y el de seno (o entre el de seno y el de seno doble) excedía 6.
  2. Análisis de Periodograma: Se generaron periodogramas de Lomb-Scargle para cada curva de luz. Las fuentes con potencia que excedía el Nivel de Falsa Alarma del 5% (confianza del 95%) fueron clasificadas como variables.
• Completitud y Confianza: Para derivar fracciones de variabilidad robustas con intervalos de confianza, los autores emplearon un formalismo bayesiano (Lafrenière et al. 2007). Calcularon la completitud de la encuesta ($p_i$) para cada objetivo generando 10,000 curvas de luz sintéticas con fases, amplitudes y periodos aleatorios, determinando la fracción detectable por su canal de análisis.

Resultos Clave

• Fracciones de Variabilidad: El estudio encontró que una fracción significativa de las enanas Y exhibe variabilidad en el infrarrojo medio.
  • Época 1: Las fracciones de variabilidad fueron $59\% \pm 15\%$ en [3.6] y $64^{+10}_{-13}\%$ en [4.5] (el análisis bayesiano arrojó fracciones más bajas de 19% y 56%, respectivamente, resaltando la dependencia del método).
  • Época 2: Las fracciones de variabilidad fueron $35^{+17}_{-11}\%$ en [3.6] y $75^{+8}_{-15}\%$ en [4.5].
  • La fracción de variabilidad es consistentemente más alta en [4.5] que en [3.6], probablemente debido a la mayor relación señal-ruido en [4.5], donde las enanas Y son más brillantes.
• Estabilidad de la Curva de Luz: Aunque muchas curvas de luz permanecieron estables durante escalas de tiempo de meses, tres objetos (WISE 0713−29, WISE 1405−55 y WISE 1639−68) mostraron cambios claros en la amplitud y/o la forma entre las épocas. Notablemente, WISE 1405−55 exhibió una transición de una variabilidad de bajo nivel en la primera época a una variabilidad significativa en la segunda, con un cambio dramático en la forma de la curva de luz en [4.5].
• Periodos de Rotación: Se identificaron cinco enanas Y con periodos de rotación seguros que oscilan entre 2.44 y 8.42 horas, derivados de periodos consistentes a través de múltiples bandas y épocas. Tres objetos adicionales tienen periodos tentativos.
• Tendencias de Amplitud: El estudio confirma que la tendencia de aumento de la amplitud máxima de variabilidad con la disminución del tipo espectral (observada en enanas L y T) se extiende a las enanas Y de tipo temprano. Las variaciones de gran amplitud (>3%) son más comunes entre las enanas Y que entre las L y T.
• Análisis de Fase: Para WISE 0359−54, que mostró una fuerte variabilidad en ambas bandas y ambas épocas, no se detectó diferencia de fase entre [3.6] y [4.5], lo cual es consistente con que estas bandas sondean capas atmosféricas de baja presión similares.

Significado y Reivindicaciones
El artículo afirma que las fracciones de variabilidad derivadas para las enanas L, T y Y (combinando este trabajo con el estudio de S. A. Metchev et al. 2015 de enanas L/T) apoyan débilmente el paradigma estándar de que la variabilidad de las enanas marrones es causada por variaciones en la estructura horizontal y/o vertical de los condensados de nubes. Específicamente, el patrón observado de alta variabilidad en enanas L, una caída en enanas T (asociada con la limpieza de nubes) y un retorno a la alta variabilidad en enanas Y (asociada con la formación de nuevos condensados como KCl, Na$_2$S y hielo de agua) se alinea con la evolución "nublado $\to$ libre de nubes $\to$ nublado" de la secuencia espectral.

Los autores enfatizan que, si bien la detección de la variabilidad depende fuertemente de la técnica de análisis (Bayesiano vs. periodograma), la conclusión general de altas fracciones de variabilidad en las enanas Y es robusta. Los resultados proporcionan una línea base para futuros estudios fotométricos y espectroscópicos de la variabilidad de las enanas Y utilizando el Telescopio Espacial James Webb (JWST). El artículo no pretende demostrar definitivamente el mecanismo físico de la variabilidad, pero señala que las tendencias observadas son consistentes con los modelos impulsados por nubes.