Blazar flares from plasma blobs crossing the broad-line region

1. Enunciato del Problema

L'articolo affronta la sfida di modellare la variabilità estrema osservata nei Quasar a Spettro Radio Piatto (FSRQ), in particolare il blazar 3C 279. Il 20 dicembre 2013, 3C 279 ha esibito un unico flare "orfano" di raggi $\gamma$ con le seguenti caratteristiche:

• Variabilità Rapida: Una scala temporale di raddoppio del flusso di sole poche ore ($\sim$2 ore).
• Spettro Rigido: Un indice fotonico molto rigido ($\Gamma_\gamma \approx 1.7$).
• Asimmetria: Una curva di luce con una rapida salita e un lento decadimento.
• Assenza di Controparte Ottica: Nessuna variabilità significativa nelle bande ottica o infrarossa, implicando una Compton dominance estrema ($L_\gamma/L_{syn} > 300$).

I modelli leptonici a una zona standard faticano a spiegare questo evento senza invocare parametri non fisici, come campi magnetici estremamente bassi, spettri di iniezione di elettroni insolitamente rigidi o meccanismi di iniezione di particelle dipendenti dal tempo ad hoc. Gli autori cercano un meccanismo fisico che produca naturalmente queste caratteristiche senza regolazioni arbitrarie.

2. Metodologia

Gli autori propongono uno scenario a due zone utilizzando il codice radiativo dipendente dal tempo EMBLEM. Il modello consiste in:

• Zona 1 (Blob Stazionario): Una regione di emissione stazionaria situata al di fuori del toro polveroso. Essa accounts per l'emissione di quiescenza (stato di bassa attività), dominando i componenti di sincrotrone (ottico/IR) e Compton Self-Sincrotrone (SSC, raggi X).
• Zona 2 (Blob Accelerante): Un blob di plasma iniettato alla base del getto che accelera mentre si propaga attraverso la Zona di Accelerazione e Collimazione (ACZ) del getto. Questo blob è responsabile del flare.
  • Legge di Accelerazione: Il blob accelera all'interno di una sezione parabolica del getto ($z \leq R_A$) dove il fattore di Lorentz di bulk scala come $\Gamma \propto z^{1/2}$, raggiungendo un fattore terminale $\Gamma_{max}$ alla transizione verso un getto conico ($z > R_A$).
  • Campi di Fotoni Esterni: Il blob si muove attraverso campi di radiazione esterni stratificati spazialmente: il Disco di Accrescimento, la Regione di Linee Larga (BLR) e il Toro Polveroso.
• Processi Fisici:
  • Processi Radiativi: Sincrotrone, SSC e scattering Compton Inverso Esterno (EIC).
  • Evoluzione delle Particelle: La distribuzione degli elettroni è modellata utilizzando un'equazione di Fokker-Planck, che tiene conto dell'iniezione, del raffreddamento radiativo (sincrotrone + IC), dell'espansione adiabatica e della fuga delle particelle.
  • Assorbimento: Assorbimento interno $\gamma\gamma$ da parte dei fotoni della BLR e assorbimento esterno da parte della Luce di Fondo Extragalattica (EBL).

Ipotesi Chiave: Il flare non è causato da un'improvvisa iniezione di particelle o da un cambiamento nel campo magnetico, ma piuttosto dall'evoluzione geometrica del fattore di Lorentz di bulk ($\Gamma$) del blob e dalla sua posizione relativa alla BLR. Mentre il blob accelera, il fattore Doppler ($\delta$) aumenta e la densità di energia dei fotoni esterni nel riferimento del blob ($U'_{ext}$) cambia, portando a un potenziamento transitorio dell'emissione EIC.

3. Contributi Chiave

• Prima Modellazione Dipendente dal Tempo: Questo è il primo tentativo di modellare con precisione l'evoluzione asimmetrica del flusso di raggi $\gamma$ dell'evento del 2013 utilizzando un trattamento completamente dipendente dal tempo delle popolazioni di elettroni e dei processi radiativi, piuttosto che approssimazioni stazionarie.
• Meccanismo Geometrico per i Flare: L'articolo dimostra che l'accelerazione di bulk combinata con la propagazione attraverso la BLR fornisce un meccanismo naturale per produrre flare ad alta energia e curve di luce asimmetriche senza richiedere iniezioni di particelle dipendenti dal tempo ad hoc.
• Risoluzione del Flare "Orfano": Il modello spiega la mancanza di variabilità ottica mostrando che il blob stazionario (Blob 1) domina l'emissione di sincrotrone, mentre il blob accelerante (Blob 2) domina l'emissione EIC solo durante la fase specifica di attraversamento della BLR.
• Previsione di Firme Ritardate: Il modello prevede specifiche firme multi-lunghezza d'onda ritardate che differiscono dai modelli standard di accelerazione da shock.

4. Risultati

• Riproduzione delle Osservazioni: Il modello riproduce con successo la distribuzione spettrale di energia (SED) osservata e le curve di luce per il flare del 2013:
  • Banda dei Raggi $\gamma$: Un flare breve e intenso con uno spettro rigido e una rapida salita seguita da un lento decadimento. La salita corrisponde al blob che accelera ed entra nella regione ad alta densità di radiazione della BLR; il lento decadimento corrisponde al blob che esce dalla BLR dove la densità di fotoni esterni diminuisce.
  • Banda Ottica/IR: Rimane costante, dominata dal Blob Stazionario 1.
• Parametri Fisici:
  • Il flare si verifica quando il blob raggiunge il raggio interno della BLR ($R_{BLR} \approx 3.1 \times 10^{17}$ cm) mentre accelera.
  • La dimensione della zona di accelerazione è stimata in $R_A \approx 4 \times 10^{17}$ cm.
  • Il fattore di Lorentz massimo è $\Gamma_{max} \approx 30$.
  • La luminosità totale dedotta è $\sim 2.1 \times 10^{46}$ erg/s (circa $1/3$ della luminosità di Eddington), il che è fisicamente plausibile.
• Previsioni:
  • Flare Ritardato nei Raggi X/EUV: Mentre il blob esce dalla BLR, il campo di fotoni esterno si indebolisce, riducendo il tasso di raffreddamento degli elettroni. Ciò porta a un graduale accumulo di densità di energia degli elettroni, prevedendo un potenziamento ritardato nell'emissione di sincrotrone e SSC (banda dei raggi X) dopo il picco dei raggi $\gamma$.
  • Ritardo VHE: Se fosse presente un'emissione ad Altissima Energia (VHE), questa sarebbe ritardata rispetto al flare di Fermi-LAT a causa della riduzione dell'assorbimento interno $\gamma\gamma$ una volta che il blob lascia la regione densa della BLR.

5. Significato

• Cambiamento di Paradigma: Lo studio sfida la visione standard secondo cui i flare dei blazar sono guidati esclusivamente da eventi impulsivi di accelerazione di particelle (ad esempio, riconnessione magnetica o shock). Invece, mette in risalto gli effetti cinematici di bulk (accelerazione e geometria) come driver primari per la variabilità estrema.
• Applicabilità: Gli autori sostengono che questo scenario sia applicabile a una classe più ampia di flare FSRQ, in particolare eventi di raggi $\gamma$ "orfani" con curve di luce asimmetriche (ad esempio, PKS 1510-089, PKS 1222+216).
• Test Osservativi: L'articolo fornisce chiari diagnostici osservativi per future campagne:
  • Cercare potenziamenti ritardati nei raggi X dopo i flare orfani di raggi $\gamma$.
  • Cercare ritardi VHE rispetto ai flare GeV.
  • Il monitoraggio multi-lunghezza d'onda simultaneo è cruciale per distinguere tra meccanismi di flare guidati dall'accelerazione e guidati dalla geometria.

In conclusione, Le Bihan et al. forniscono una spiegazione robusta e fisicamente motivata per il flare del 2013 di 3C 279, attribuendo la variabilità estrema all'interazione tra l'accelerazione di bulk del getto e la struttura spaziale dei campi di radiazione esterni, offrendo un nuovo quadro per comprendere la fenomenologia ad alta energia dei blazar.