Bayesian Evidence for Inspiraling Hotspot Motion in the Galactic Center

Sintesi Tecnica: Evidenza Bayesiana per il Moto di Inspiral di un Hotspot nel Centro Galattico

Definizione del Problema
Le osservazioni polarimetriche del buco nero supermassiccio Sagittarius A* (Sgr A*) durante gli episodi di flare rivelano traiettorie evolutive a forma di loop nel piano degli Stokes $Q-U$, note come "loop $Q-U$". Precedenti analisi dell'osservazione dell'Atacama Large Millimeter/submillimeter Array (ALMA) dell'11 aprile 2017 hanno modellato con successo il loop esterno primario utilizzando hotspot transitori su orbite circolari (kepleriane) fisse. Tuttavia, tali modelli non riescono a riprodurre il successivo loop interno osservato nello stesso dataset, poiché la parte successiva dei dati è incoerente con qualsiasi traiettoria strettamente periodica. Mentre lavori precedenti hanno testato le deviazioni dal moto kepleriano all'interno di un framework di orbita circolare, era necessaria una prescrizione cinematica più generale e non circolare per spiegare l'intera evoluzione temporale del segnale polarizzato, in particolare la transizione tra il loop primario e quello interno.

Metodologia
Gli autori sviluppano un framework bayesiano per inferire la cinematica dipendente dal tempo di un hotspot emissivo da curve di luce polarizzate. Il nucleo della metodologia consiste in:

• Modello Cinematico Generalizzato: La traiettoria della regione emissiva è modellata come una sequenza continua e a tratti di $K$ segmenti cinematici. Ogni segmento è descritto da parametri di velocità radiale ($\beta_r$) e azimutale ($\beta_\phi$) rispetto a un flusso geodetico circolare (kepleriano). Questi parametri interpolano tra un flusso kepleriano generalizzato e una caduta libera dall'infinito con momento angolare nullo.
  • $\beta_r < 1$ introduce una componente di inspiral radiale.
  • $\beta_\phi < 1$ riduce la velocità angolare al di sotto del valore kepleriano.
• Trasferimento Radiativo: Il modello calcola l'emissione polarizzata nello spaziotempo di Kerr utilizzando il codice di ray-tracing analitico adattivo AART. Esso incorpora un trattamento di "luce lenta" per tenere conto degli effetti del tempo di percorrenza della luce, mappando i tempi di emissione sugli tempi di arrivo all'osservatore. La polarizzazione è calcolata tramite il trasporto parallelo del quadrivettore di polarizzazione usando la costante di Penrose-Walker, tenendo conto del lensing gravitazionale e del redshift.
• Inferenza Statistica: L'analisi impiega un algoritmo Markov Chain Monte Carlo a temperatura parallela (PT-MCMC) (utilizzando reddemcee) per campionare la distribuzione a posteriori dei parametri del modello. Il vettore dei parametri include la posizione iniziale, i parametri di keplerianità per ogni segmento, i pesi della durata del segmento, gli offset di polarizzazione di fondo, un angolo di rotazione del cielo (PA), le componenti del campo magnetico e gli indici spettrali.
• Confronto dei Modelli: Diversi profili cinematici vengono confrontati utilizzando l'evidenza bayesiana (calcolata tramite integrazione termodinamica e il metodo dei stepping-stones) e la statistica chi-quadro ridotta ($\chi^2_{\text{eff}}$). Lo studio analizza due sottoinsiemi di dati: una finestra più breve ($t_{\text{obs}} \sim 230 M$) che copre solo il loop primario, e una finestra estesa ($t_{\text{obs}} \sim 350 M$) che include il loop interno.

Contributi Chiave e Risultati
Il paper presenta un confronto sistematico tra modelli con $K=1, 2$ e $3$ segmenti cinematici rispetto ai dati ALMA del 2017:

1. Limitazioni dei Modelli a Segmento Singolo:

   • Per il solo loop primario ($t_{\text{obs}} \sim 230 M$), i dati mostrano una lieve preferenza per un inspiral non kepleriano ($\beta_r < 1$) rispetto a un'orbita strettamente kepleriana. Tuttavia, estendere questi modelli a segmento singolo oltre il loop primario risulta in un marcato degrado della qualità del fit.
   • Nessun modello a segmento singolo fornisce un fit adeguato per l'intero dataset ($t_{\text{obs}} \sim 350 M$).

2. Evidenza per un Moto di Inspiral a Multi-Segmento:

   • Quando viene incluso il loop interno, i modelli che permettono alla cinematica di cambiare lungo la traiettoria ($K=2$) sono fortemente preferiti.
   • Il modello con la massima evidenza ($E_2$) presenta due segmenti con parametri radiali e azimutali liberi ($\beta_r, \beta_\phi$) ma un singolo indice spettrale condiviso. Questo modello è favorito rispetto al miglior controparte kepleriana a due segmenti di $\Delta \ln Z \approx 46$.
   • I dati non supportano un modello a tre segmenti ($K=3$); la complessità aggiuntiva non è giustificata dall'evidenza, poiché la distribuzione a posteriori collassa verso la soluzione a due segmenti.

3. Cinematica Inferita:

   • Segmento 1 (Loop Primario): Il flusso esibisce un lieve inspiral ($\beta_{r,1} \approx 0.95$) ed è coerente con una rotazione kepleriana ($\beta_{\phi,1} \approx 0.96$).
   • Segmento 2 (Loop Interno): Oltre un tempo di rottura a $t_{\text{break}} \approx 275 M$ (tempo di emissione), il flusso subisce un significativo cambiamento cinematico. L'afflusso radiale accelera ($\beta_{r,2} \approx 0.86$) e il supporto rotazionale scende bruscamente ($\beta_{\phi,2} \approx 0.24$), indicando una transizione verso un flusso dominato dall'advezione e sub-kepleriano.
   • La traiettoria termina ben al di fuori dell'Orbita Circolare Ultima Stabile (ISCO).

4. Parametri Fisici:

   • Gli offset di polarizzazione di fondo e l'angolo di rotazione del cielo inferiti sono coerenti con misurazioni indipendenti precedenti.
   • L'indice spettrale ($\alpha_\nu \approx 0.35$) è coerente con i valori misurati immediatamente dopo il flare di raggi X.
   • La geometria del campo magnetico è inferita avere una componente verticale dominante, coerente con argomentazioni precedenti ma derivata qui senza forti a priori.

Significatività e Rivendicazioni
Il paper rivendica un'evidenza statisticamente significativa che la regione emissiva in Sgr A* segua una traiettoria di inspiral non geodesica piuttosto che un'orbita strettamente kepleriana. Questa preferenza è debole quando viene considerato solo il loop primario, ma si rafforza drasticamente quando viene modellato l'intero dataset, incluso il loop interno. I risultati suggeriscono che l'hotspot non sta semplicemente orbitando, ma viene trasportato verso l'interno da un flusso di accrezione dominato dall'advezione, perdendo il supporto rotazionale man mano che si avvicina all'orizzonte degli eventi.

Gli autori posizionano questo lavoro come un passo necessario per caratterizzare il sistema Sgr A* e comprendere la dinamica di accrezione magnetizzata. Notano che, sebbene il loro modello descriva con successo i dati, l'assunzione di cinematica costante a tratti sia un'approssimazione; i lavori futuri dovrebbero esplorare traiettorie di accelerazione fluida e includere la variabilità di spin e inclinazione. Lo studio sottolinea che questi modelli teorici migliorati sono essenziali per interpretare le future osservazioni ad alta risoluzione angolare del next-generation Event Horizon Telescope (ngEHT) e del Black Hole Explorer (BHEX).