Particle acceleration signatures in the time-dependent one-zone synchrotron self-Compton model of blazar flares

Planteamiento del Problema
Los estallidos (flares) de los blázares exhiben variaciones rápidas de flujo a través del espectro electromagnético, que a menudo ocurren en escalas de tiempo de minutos a días. Aunque el origen físico de esta variabilidad no se comprende completamente, generalmente se atribuye ya sea a cambios macroscópicos en la región de emisión (por ejemplo, geometría, fuerza del campo magnético) o a cambios en la distribución de energía de la población de partículas no térmicas. Si bien a menudo se requieren modelos de múltiples zonas sofisticados para describir la emisión estocástica completa de los blázares, el estudio de los estallidos rápidos específicamente puede aproximarse mediante una región de emisión homogénea de una sola zona. Sin embargo, la comprensión sistemática de cómo diferentes mecanismos de inyección y aceleración de partículas se manifiestan en las curvas de luz (LCs) multi-longitud de onda observables sigue siendo incompleta. Los estudios actuales suelen centrarse en el comportamiento estocástico a largo plazo o en eventos individuales específicos sin un marco unificado para distinguir entre procesos físicos genéricos (por ejemplo, inyección frente a aceleración por choque o aceleración estocástica) basándose únicamente en la morfología del estallido.

Metodología
Los autores emplean el código EMBLEM (Evolutionary Modeling of BLob EMission) para simular modelos de sincrotrón auto-Compton (SSC) dependientes del tiempo de una sola zona, adaptados para objetos BL Lac. El estudio se centra en la evolución de la distribución de electrones $N_e(\gamma, t)$ dentro de un único lóbulo de plasma, gobernada por una ecuación de Fokker-Planck que incluye inyección de partículas, enfriamiento radiativo (sincrotrón e inverso-Compton), pérdidas adiabáticas, escape de partículas y mecanismos de aceleración.

El estudio explora sistemáticamente cinco escenarios distintos para modelar la fase de estallido:

1. Escenario 1a (Inyección): Inyección adicional de una distribución de partículas de ley de potencia durante un breve periodo, equilibrada por el enfriamiento y el escape.
2. Escenario 1b (Inyección + Expansión Adiabática): Similar al 1a, pero incluye los efectos de la expansión adiabática del lóbulo.
3. Escenario 2a (Aceleración Fermi-I): Aceleración de choque difusiva (DSA) actuando sobre una población de electrones fríos o ligeramente relativistas, caracterizada por un factor de Lorentz máximo dependiente del tiempo.
4. Escenario 2b (Reaceleración Fermi-I): DSA actuando sobre una población de electrones relativistas existente responsable del estado de quietud.
5. Escenario 3a (Aceleración Fermi-II): Aceleración estocástica (Fermi II) actuando sobre una población fría/ligeramente relativista mediante la dispersión resonante en ondas MHD.
6. Escenario 3b (Reaceleración Fermi-II): Aceleración estocástica actuando sobre la población relativista existente.

Las simulaciones asumen una dispersión de "esfera dura" (tiempos de aceleración y escape independientes de la energía) para mantener un marco simplificado. Las curvas de luz multi-longitud de onda resultantes se analizan en cuatro bandas de energía: Óptica, Rayos X, Gamma de Alta Energía (HE) y Gamma de Muy Alta Energía (VHE). Las métricas clave extraídas incluyen la amplitud de variabilidad (VA), los tiempos de ascenso y decaimiento, la asimetría, la presencia de mesetas (plateaus) y los retrasos temporales entre bandas.

Contribuciones Clave y Resultados
El estudio identifica firmas características en los perfiles de las curvas de luz que permiten la discriminación entre los principales procesos físicos que impulsan los estallidos de los blázares:

• Morfología de la Curva de Luz y Simetría:

  • Escenarios de Inyección (1a/1b): Producen curvas de luz relativamente simétricas con un ascenso cóncavo. A menudo se alcanza una fase de meseta si la duración del estallido es suficiente. En el Escenario 1b, la expansión adiabática induce un descenso sistemático del flujo durante la meseta.
  • Aceleración Fermi-I (2a/2b): La forma del ascenso transiciona de convexa a cóncava dependiendo de la escala de tiempo de aceleración. Una aceleración ineficiente (escalas de tiempo largas) conduce a perfiles asimétricos donde el tiempo de ascenso supera al tiempo de decaimiento. El Escenario 2b (reaceleración) muestra una diferencia marcada en las formas de ascenso entre las bandas de baja energía (óptica/HE) y las de alta energía (Rayos X/VHE).
  • Aceleración Fermi-II (3a/3b): La parte ascendente de la curva de luz es consistentemente convexa. En el Escenario 3b (reaceleración), se observa un retraso temporal distintivo donde los picos de Rayos X y VHE ocurren significativamente antes que los picos de óptica/HE. Este escenario también exhibe un régimen de "alta dominancia de Compton" (CD) donde el pico SSC se desplaza hacia energías más altas.

• Retrasos Temporales Dependientes de la Energía:

  • En los escenarios de Inyección y Reaceleración, el inicio del estallido ocurre aproximadamente de forma simultánea en todas las bandas.
  • En los escenarios de aceleración inicial Fermi-I (2a) y Fermi-II (3a), el inicio del ascenso del flujo se retrasa en las bandas de mayor energía con respecto a las de menor energía. Este retraso es más pronunciado en el escenario Fermi-I.
  • El Escenario 3b exhibe una firma única donde los picos de Rayos X/VHE preceden a los picos de óptica/HE, como resultado directo de la rápida aceleración de la población existente seguido de un fuerte enfriamiento por Compton inverso.

• Histéresis:

  • Los escenarios de inyección muestran una histéresis mínima en la Distribución de Energía Espectral (SED).
  • Los escenarios de aceleración y reaceleración (2a, 2b, 3a, 3b) muestran una histéresis significativa, donde los picos espectrales se desplazan hacia energías más altas durante el ascenso y hacia energías más bajas durante el decaimiento. Este efecto es más pronunciado en el escenario de reaceleración Fermi-II (3b).

• Caso de Estudio (Mrk 421):
  Los autores aplicaron estos escenarios a un estallido de 2013 de Mrk 421. El escenario simple de inyección (1a) proporcionó la mejor reproducción de las curvas de luz observadas en Rayos X y VHE. El escenario de reaceleración Fermi-I (2b) también podría ajustarse a los datos en el límite de una aceleración muy eficiente, mientras que el escenario de reaceleración Fermi-II (3b) no pudo reproducir la forma de ascenso de los Rayos X bajo las restricciones dadas.

Significancia y Reivindicaciones
El artículo argumenta que, si bien los modelos SSC de una sola zona son una simplificación, son suficientes para identificar firmas observacionales genéricas que distingan entre diferentes mecanismos de aceleración de partículas. Los autores afirman que, con datos multi-longitud de onda de suficiente calidad y bien muestreados, se puede discriminar entre inyección, aceleración por choque y aceleración estocástica basándose en:

1. El tiempo relativo de los inicios y picos de los estallidos a través de las bandas de energía.
2. La forma de las partes ascendentes y descendentes de las curvas de luz (convexa vs. cóncava, simétrica vs. asimétrica).
3. La presencia y naturaleza de las mesetas.
4. El grado de histéresis en la SED.

El estudio enfatiza que comparar curvas de luz de al menos dos bandas de energía (específicamente muestreando tanto el ascenso de sincrotrón como el de SSC) es esencial para desentrañar estos escenarios. Los autores señalan que el próximo Observatorio del Telescopio de Cherenkov Array (CTAO) será particularmente sensible a estas firmas debido a su resolución temporal mejorada en la banda VHE. El artículo concluye que expandir tanto la cobertura observacional como el modelado teórico (por ejemplo, incorporando campos de fotones externos para FSRQ o efectos de múltiples zonas) es necesario para resolver completamente los procesos físicos en los estallidos de los blázares.