27 years of Spaceborne IR Astronomy: An ISO, Spitzer, WISE and NEOWISE Survey for Large-Amplitude Variability in Young Stellar Objects

Resumen Técnico: 27 años de astronomía infrarroja desde el espacio: un estudio de ISO, Spitzer, WISE y NEOWISE para la variabilidad de gran amplitud en objetos estelares jóvenes

Planteamiento del problema
Los Objetos Estelares Jóvenes (YSO, por sus siglas en inglés) exhiben variabilidad en todo el espectro electromagnético, pero la caracterización de la variabilidad de gran amplitud (definida aquí como cambios de flujo $\ge$ factor de 2) en longitudes de onda $\lambda > 3$ $\mu$m es crítica para comprender los procesos de acreción de masa. Si bien estudios previos han identificado ráfagas de acreción episódicas (por ejemplo, tipos FU Ori y EX Lupi), la frecuencia, las amplitudes y las duraciones de estos eventos permanecen insuficientemente restringidas debido a líneas de base temporales limitadas y datos inhomogéneos. Específicamente, la contribución de las ráfagas sostenidas al ensamblaje total de la masa estelar es incierta. Para abordar esto, se requiere un conjunto de datos continuo y de larga línea de base para distinguir entre ráfagas sostenidas, eventos de desvanecimiento y otros modos de variabilidad a través de toda la secuencia evolutiva de los YSO, desde protoestrellas profundamente incrustadas hasta estrellas de la secuencia principal temprana con discos.

Metodología
Los autores construyeron curvas de luz compuestas que abarcan 27 años (1997–2024) para YSO dentro de 500 pc del Sol. El conjunto de datos integra fotometría de tres telescopios infrarrojos espaciales a través de cinco fases de misión:

1. Observatorio Infrarrojo Espacial (ISO): Datos de ISOCAM de 1997 en la banda LW2 (6.7 $\mu$m).
2. Telescopio Espacial Spitzer: Datos de la misión criogénica (2003–2009) y de la misión cálida (2009–2019) utilizando la cámara IRAC en las bandas de 3.6 $\mu$m (I1), 4.5 $\mu$m (I2), 5.8 $\mu$m (I3) y 8.0 $\mu$m (I4).
3. WISE/NEOWISE: Datos de la encuesta inicial de WISE (2010) y de la misión reactivada NEOWISE (2013–2024) en las bandas de 3.4 $\mu$m (W1) y 4.6 $\mu$m (W2).

Ensamblaje y Calibración de Datos:

• Selección de Fuentes: La muestra consiste en 221 YSO identificados por estudios de Spitzer como poseedores de discos o envolventes, con detecciones en ISO, Spitzer y NEOWISE.
• Calibración Cruzada: Los autores desarrollaron un esquema de calibración cruzada riguroso para alinear los perfiles de filtro dispares. Los flujos de ISOCAM LW2 se ajustaron para coincidir con Spitzer I3, I4 y un flujo combinado de I3+I4. Del mismo modo, las magnitudes de WISE/NEOWISE W1 y W2 se ajustaron para alinearse con Spitzer I1 e I2, teniendo en cuenta las diferencias en la forma de los filtros y las distribuciones de energía espectral (SED) de las fuentes.
• Manejo de la Saturación: Se tuvo especial cuidado en identificar y marcar los píxeles saturados en los datos de NEOWISE, comparando los catálogos de IRSA y un-Timely para evaluar los efectos de saturación en fuentes brillantes.
• Construcción de Curvas de Luz: Se generó una curva de luz compuesta por fuente: una curva de 3.6 $\mu$m (I1 + W1) y una de 4.5 $\mu$m (I2 + W2 + ISO LW2). La variabilidad se cuantificó mediante $\Delta m$ (cambio de magnitud relativo a la magnitud mínima observada).

Contribuciones Clave y Resultados
El estudio identificó 39 YSO que exhiben variabilidad de gran amplitud ($\Delta m \ge 0.75$ mag, o un cambio de flujo de un factor de 2) y los clasificó en cuatro categorías provisionales basadas en la morfología de sus curvas de luz de 27 años:

1. Ráfagas Completadas (4 fuentes): Estas fuentes muestran un aumento en el flujo (factor $\ge$ 2), una meseta de flujo alto sostenida que dura $>5$ años y un descenso posterior hacia los niveles previos a la ráfaga. Ejemplos incluyen la protoestrella V1647 Ori (J05461313-0006045) y tres estrellas de la secuencia principal temprana con discos. Las duraciones de estos eventos oscilan entre 6 y 17 años.
2. Ráfagas en Curso (3 fuentes): Estas fuentes muestran un aumento en el flujo seguido de una fase de flujo alto sostenida que continúa a través de la época final de 2024. Las duraciones se estiman en $>13$ a $>20$ años.
3. Desvanecimientos (6 fuentes): Estas fuentes exhiben un descenso sostenido en el flujo de $\ge 0.75$ mag a lo largo de la curva de luz, interpretado como la fase de decaimiento de ráfagas que comenzaron antes de las observaciones de ISO en 1997.
4. Fluctuadores (26 fuentes): La categoría más común, estas fuentes muestran variaciones de gran amplitud en escalas de tiempo de meses a años, pero carecen de la meseta de flujo alto sostenida o del descenso monotónico característico de las otras categorías. Este grupo incluye tanto protoestrellas como estrellas de la secuencia principal temprana con discos.

Características de la Variabilidad:

• Escalas de Tiempo: El análisis de los diagramas $\Delta m$ vs. $\Delta t$ revela que las mayores amplitudes para ráfagas y desvanecimientos ocurren en escalas de tiempo $>1000$ días. En contraste, los fluctuadores muestran una variabilidad significativa incluso en escalas de tiempo $<100$ días, con una amplitud del percentil 90 que permanece relativamente constante en todas las escalas de tiempo muestreadas.
• Evolución del Color: Se utilizaron diagramas de color-magnitud (I1-I2 vs. I2 y W1-W2 vs. W2) para probar la hipótesis de la extinción. Los autores encuentran que, para las ráfagas y desvanecimientos, los cambios de color son inconsistentes con los vectores de extinción de polvo; en su lugar, los colores permanecen constantes o se enrojecen ligeramente a medida que aumenta el brillo. Esto respalda la interpretación de que estos eventos son impulsados por cambios en la luminosidad intrínseca debido a variaciones en la tasa de acreción, en lugar de por la obstrucción. Los fluctuadores muestran una distribución de pendientes más amplia, lo que sugiere una mezcla de variabilidad impulsada por la acreción y alguna contribución de la extinción.
• Dinámica de las Ráfagas: El estudio señala que las fases de ascenso y descenso pueden ser rápidas (1–2 años) o lentas (más de una década). Los datos sugieren que las fases de ascenso y decaimiento de las ráfagas son a menudo más cortas que la fase de meseta sostenida.

Significado y Reivindicaciones
El artículo afirma que la línea de base de 27 años proporciona las mejores restricciones actuales sobre las duraciones de las ráfagas impulsadas por acreción en los YSO. Los autores concluyen que:

• La variabilidad de gran amplitud está dominada por cambios de luminosidad impulsados por la acreción en lugar de la extinción.
• La categoría de "fluctuadores" representa un modo de variabilidad distinto y prevalente, lo que potencialmente indica que los estados de alta acreción son a menudo inestables o que estas fuentes están experimentando ráfagas donde las fases de ascenso/descenso caen fuera de la ventana de 27 años.
• La duración observada de las ráfagas (que oscila entre 6 y $>20$ años) y la prevalencia de los fluctuadores sugieren que la contribución de la acreción episódica al ensamblaje de la masa estelar es compleja y puede involucrar episodios de acreción inestables y no monotónicos.
• Los autores evitan explícitamente clasificar estas fuentes en las categorías estándar de FU Ori o EX Lupi basándose únicamente en las curvas de luz, citando la falta de datos espectroscópicos simultáneos y el potencial de un muestreo incompleto para ocultar la verdadera naturaleza de los eventos.

El estudio enfatiza que, si bien los datos apoyan la imagen de que la mayor parte de la masa estelar se ensambla durante ráfagas sostenidas, los mecanismos específicos que impulsan el comportamiento de los "fluctuadores" y la relación entre las fluctuaciones y las ráfagas completas siguen siendo preguntas abiertas para investigaciones futuras.