A universal scaling between damping time and period of quasi-periodic pulsations from solar EUV brightenings to X-ray stellar flares

Resumen Técnico: Una escala universal entre el tiempo de amortiguamiento y el periodo de las pulsaciones cuasi-periódicas

Planteamiento del problema
Las pulsaciones cuasi-periódicas (QPP, por sus siglas en inglés) son una firma ubicua de la liberación impulsiva de energía en llamaradas solares y estelares. Si bien las QPP han sido estudiadas extensamente en llamaradas solares y estelares de X de gran escala, su presencia en eventos de escala mucho menor, específicamente en los brillos en Ultravioleta Extremo (EUV), ha sido confirmada recientemente. Una pregunta fundamental sigue siendo si las QPP observadas a través de escalas espaciales y energéticas tan dispares —desde brillos EUV del tamaño de una nanoflarea hasta masivas llamaradas estelares— comparten un origen físico común. Un trabajo previo de I. H. Cho et al. (2016) estableció una relación de escala entre el tiempo de amortiguamiento y el periodo de oscilación para llamaradas solares y estelares de rayos X, sugiriendo un mecanismo unificado. Sin embargo, no estaba claro si esta escala se extiende a los brillos EUV más pequeños y fríos.

Metodología
Los autores realizaron un análisis estadístico de las QPP en la fase de decaimiento de eventos impulsivos utilizando datos de imagen de alta resolución.

• Fuentes de datos: El estudio analizó 2,146 brillos EUV observados por el Solar Orbiter/Extreme Ultraviolet Imager (EUI) en 174 Å y 300 llamaradas solares EUV observadas por el Solar Dynamics Observatory/Atmospheric Imaging Assembly (SDO/AIA) en 171 Å.
• Selección de la muestra: Del conjunto de datos inicial, los autores identificaron 185 brillos EUV y 89 llamaradas solares (que comprenden eventos de clases C, M y X) que exhibieron QPP amortiguadas estadísticamente significativas.
• Procesamiento de señales: Para aislar los componentes oscilatorios, las curvas de luz de la fase de decaimiento fueron procesadas mediante la Descomposición de Modos Empíricos (EMD). Se eliminó la función de modo intrínseco (IMF) con la escala de tiempo más larga para eliminar la tendencia de fondo de variación lenta. Las señales residuales se ajustaron con una función armónica amortiguada para extraer el periodo de oscilación ($P$) y el tiempo de amortiguamiento ($\tau$).
• Validación estadística: Las oscilaciones se retuvieron solo si cumplían criterios estrictos: potencia espectral superior al umbral de significancia del 90% (vía aleatorización de Fisher), incertidumbres menores que los valores ajustados, periodos superiores a cinco veces la cadencia y consistencia entre los ajustes en el dominio del tiempo y los espectros de Fourier en el dominio de la frecuencia.
• Análisis de escala: La relación entre $\tau$ y $P$ se analizó en el espacio logarítmico utilizando regresión lineal bayesiana. Este enfoque tuvo en cuenta las incertidumbres de medición y la dispersión intrínseca. Los autores compararon modelos que asumían un exponente de escala universal único frente a modelos con exponentes distintos para diferentes poblaciones de eventos (brillos EUV frente a llamaradas; EUV frente a rayos X) utilizando validación cruzada de dejar uno fuera con muestreo de importancia suavizado por Pareto.

Resultos clave

• Detección de QPP amortiguadas: El estudio confirmó que las QPP amortiguadas de periodo único están presentes en las fases de decaimiento tanto en los brillos EUV de pequeña escala como en las llamaradas EUV solares de gran escala.
• Factores de calidad: Las distribuciones de los factores de calidad ($Q = \tau/P$) para los brillos EUV y las llamaradas EUV resultaron ser estadísticamente distinguibles (prueba Kolmogorov–Smirnov $p=0.002$), mostrando los brillos EUV una cola de alto $Q$ más extendida. Sin embargo, ambas distribuciones exhibieron valores modales similares (0.6 para brillos, 0.5 para llamaradas).
• Escala universal: Se encontró una tendencia monotónica clara entre el tiempo de amortiguamiento y el periodo de oscilación en todos los conjuntos de datos.
  • Para los brillos EUV y las llamaradas EUV combinados, los datos se describen bien mediante una relación de ley de potencia única $\tau \propto P^{\alpha}$ con un índice $\alpha \approx 0.87$.
  • La comparación de modelos bayesianos indicó que permitir pendientes separadas para los brillos y las llamaradas no mejoró significativamente el rendimiento predictivo, lo que respalda una relación de escala compartida dentro de los datos de EUV.
  • Al combinarse con los datos de llamaradas solares y estelares de rayos X reportados previamente (I. H. Cho et al. 2016), el conjunto completo (que abarca los regímenes de EUV y rayos X) continúa siguiendo una escala de ley de potencia única con $\alpha \approx 0.92$. La inclusión de exponentes distintos para las poblaciones de EUV y rayos X no fue estadísticamente requerida.

Significancia y afirmaciones
La afirmación principal de este trabajo es que la escala entre el tiempo de amortiguamiento y el periodo de oscilación para las QPP de fase de decaimiento es universal, persistiendo a través de un vasto rango de energías y escalas espaciales, desde los brillos EUV solares hasta las llamaradas estelares de rayos X.

• Implicaciones físicas: La persistencia de esta escala sugiere que las QPP en estos diversos eventos están gobernadas por un mecanismo físico subyacente común. Los autores discuten dos candidatos potenciales consistentes con el índice de escala observado ($\sim0.87-0.92$): ondas de flexión (kink waves) estacionarias y oscilaciones magnetoacústicas lentas estacionarias. Los resultados son particularmente consistentes con las oscilaciones de modo lento, donde la conducción térmica puede actuar como el mecanismo de amortiguamiento.
• Calentamiento coronal: Los hallazgos apoyan la interpretación de que los brillos EUV de pequeña escala pueden representar eventos de llamarada elementales, reforzando la hipótesis de que podrían desempeñar un papel en el calentamiento coronal.
• Seismología coronal: El estudio destaca el potencial diagnóstico de las QPP. Si las oscilaciones se identifican como modos magnetoacústicos lentos, el periodo medido y la escala de longitud del evento pueden utilizarse para estimar la temperatura del plasma en brillos EUV no resueltos o para inferir longitudes de bucle en llamaradas estelares, extendiendo así la utilidad de la seismología coronal a un rango más amplio de fenómenos astrofísicos.

Los autores mantienen una postura modesta, señalando que, si bien la escala respalda un mecanismo común, la identificación robusta del modo de onda específico (flexión vs. lento) requiere futuras observaciones multi-longitud de onda para distinguir entre las interpretaciones magnetohidrodinámicas competidoras.