Study of an interaction between the jet and an interstellar medium of M87 with a spectral analysis by using Chandra

Resumen Técnico: Estudio de una interacción entre el chorro y un medio interestelar de M87 mediante un análisis espectral utilizando Chandra

Planteamiento del problema
La galaxia de radio M87, ubicada en el centro del cúmulo Virgo, exhibe un prominente chorro (jet) que se extiende 20 segundos de arco. Aunque el origen de los rayos gamma TeV en M87 se atribuye generalmente al núcleo o a nudos específicos como HST-1, los mecanismos físicos que impulsan la interacción entre el chorro y el medio interestelar (ISM) circundante siguen bajo investigación. Los análisis espectrales previos del chorro de M87 utilizando datos de Chandra (por ejemplo, Wilson & Yang 2002; Perlman & Wilson 2005; Sun et al. 2018) ajustaron principalmente los espectros de energía de rayos X con modelos de ley de potencia, asumiendo emisión sincrotrón. Sin embargo, estos estudios enfrentaron limitaciones:

1. Estadísticas insuficientes: Muchos análisis dependieron de tiempos de exposición cortos (por ejemplo, 13–38 ks), lo que resultó en estadísticas deficientes donde modelos complejos podrían parecer aceptables simplemente debido a la baja calidad de los datos.
2. Errores de sustracción del fondo: Los trabajos anteriores a menudo utilizaron una región de fondo común para todo el chorro. Esto es físicamente incorrecto porque el gas caliente del cúmulo Virgo, que contribuye al fondo, varía en intensidad con la distancia al núcleo.
3. Negligencia de componentes térmicos: Existe la posibilidad de que el chorro interactúe con el ISM mediante choques, calentando el gas para producir un componente térmico (bremsstrahlung) además de la emisión sincrotrón no térmica. Este componente puede haber sido omitido en estudios de baja estadística.

Metodología
Para abordar estos problemas, este estudio utiliza un conjunto de datos significativamente mayor del Observatorio de Rayos X Chandra, que cubre observaciones desde el año 2000 hasta 2018 con un tiempo de exposición total de aproximadamente 800 ks.

• Selección de datos y tratamiento de "pile-up": El análisis se centra en cuatro regiones específicas: el núcleo, el nudo HST-1, el nudo D y el nudo A. Debido al alto brillo del núcleo y de HST-1, los datos se seleccionaron cuidadosamente para gestionar el "pile-up" (donde múltiples fotones impactan en un solo píxel, distorsionando el espectro). Los datos fueron filtrados para asegurar que las fracciones de pile-up estuvieran por debajo del 5% en la mayoría de los análisis, aunque se aceptó un 10% para el núcleo y HST-1 en el conjunto de datos anterior para maximizar las estadísticas.
• Modelado del fondo: A diferencia de estudios previos, se seleccionaron regiones de fondo individuales para cada nudo (al norte y al sur de la fuente) para dar cuenta de la variación espacial del gas caliente del cúmulo Virgo. Para el nudo A, se eligió una región de fondo específica para evitar caídas conocidas en los rayos X blandos.
• Análisis espectral: Los espectros de energía se extrajeron utilizando el software CIAO (v4.15) y se analizaron con XSPEC. El rango de energía se estableció de 0.5 keV a 10 keV para evitar problemas de calibración en energías bajas.
• Ajuste de modelos:
  1. Se probó un modelo simple de ley de potencia (que representa la emisión sincrotrón) con absorción fotoeléctrica.
  2. Si el ajuste de la ley de potencia era estadísticamente rechazado (indicado por un $\chi^2$ alto y un exceso blando), se aplicó un modelo combinado de ley de potencia y un modelo APEC (bremsstrahlung térmico).
  3. Se derivaron parámetros físicos, incluyendo densidad neutra, densidad de plasma y densidad de electrones acelerados, a partir de los mejores modelos de ajuste.
  4. El flujo de bremsstrahlung no térmico resultante de la interacción de electrones acelerados con el ISM se calculó y se comparó con las observaciones de rayos gamma de Fermi.

Resultados clave

1. Características espectrales:

   • Núcleo y Nudo D: Los espectros de energía de rayos X tanto para el núcleo como para el nudo D se describen bien mediante un modelo simple de ley de potencia, lo cual es consistente con la emisión sincrotrón.
   • Nudo HST-1 y Nudo A: Para estas regiones, un modelo de pura ley de potencia es estadísticamente rechazado.
     • HST-1: En el conjunto de datos combinado (incluyendo datos con 10% de pile-up), un modelo de ley de potencia arroja una probabilidad de azar de $2.6 \times 10^{-4}$. Añadir un componente térmico (APEC) proporciona un ajuste razonable.
     • Nudo A: El modelo de ley de potencia es fuertemente rechazado ($\chi^2 = 628.76$ para 355 grados de libertad, probabilidad $1.6 \times 10^{-17}$). Un modelo combinado de ley de potencia y bremsstrahlung térmico ajusta bien los datos.

2. Propiedades del componente térmico:

   • El componente térmico detectado en los nudos HST-1 y A se caracteriza por una temperatura de 0.2 keV y una abundancia metálica de 0.00 (esencialmente cero).
   • La abundancia metálica es significativamente menor que la del gas caliente circundante del cúmulo Virgo (que es de aproximadamente ~1 solar), lo que sugiere que el gas calentado se origina dentro del propio chorro en lugar del medio ambiente.
   • La relación de flujo térmico en el nudo A es aproximadamente el 22% del flujo total.

3. Parámetros físicos:

   • Densidad de plasma: Derivada de la normalización APEC, la densidad electrónica en el nudo A es de aproximadamente $11.0 \, \text{cm}^{-3}$.
   • Calentamiento por choque: La presencia de un componente térmico con baja metalicidad y alta temperatura en relación con el medio ambiente respalda la hipótesis de que el gas dentro del chorro es calentado por una onda de choque resultante de la interacción con el ISM.
   • Densidad neutra: Calculada a partir de la diferencia entre la columna de densidad derivada de rayos X y las observaciones de radio de 21 cm, indicando la presencia de material frío alrededor del chorro.

4. Bremsstrahlung no térmico y comparación con Fermi:

   • El estudio calculó el flujo esperado de bremsstrahlung no térmico de los electrones acelerados interactuando con el ISM.
   • El flujo calculado del núcleo y del nudo HST-1 es órdenes de magnitud inferior al flujo observado en GeV por Fermi.
   • Conclusión: No hay contribución del bremsstrahlung no térmico en el núcleo o en HST-1 al flujo observado de Fermi, lo que refuerza que el origen de la emisión TeV/GeV probablemente se encuentra en otro lugar (por ejemplo, el núcleo o diferentes mecanismos de aceleración de partículas).

Significancia y Reivindicaciones
El artículo afirma que, al utilizar datos de alta estadística (800 ks) y una sustracción de fondo espacialmente dependiente y correcta, ha identificado con éxito un componente térmico en los nudos del chorro de M87, HST-1 y A, en los espectros de rayos X. Este hallazgo desafía las interpretaciones previas que dependían únicamente de modelos de ley de potencia.

Los autores concluyen que:

1. La emisión térmica es probablemente gas dentro del chorro calentado por una interacción de choque con el ISM, evidenciado por la baja abundancia metálica y la alta temperatura.
2. El mecanismo de bremsstrahlung no térmico (interacción electrón-ISM) no explica el flujo de rayos gamma en el rango de GeV observado de M87, reforzando que el origen de la emisión TeV/GeV probablemente reside en otra parte (por ejemplo, el núcleo o diferentes mecanismos de aceleración de partículas).
3. Los estudios previos pueden haber pasado por alto estas características térmicas debido a estadísticas insuficientes y un modelado de fondo inadecuado.

El estudio no propone nuevos montajes experimentales, sino que enfatiza la necesidad de un análisis de datos de archivo de alta calidad para resolver el espectro de energía multi-longitud de onda y las condiciones físicas de los chorros relativistas.