Spectroscopic Variability of the Broad H$β$ Emission Line in Sloan Digital Sky Survey Quasars

Sintesi Tecnica: Variabilità Spettroscopica della Riga di Emissione Larga H$\beta$ nei Quasar della Sloan Digital Sky Survey

Definizione del Problema
Sebbene i motori centrali dei quasar siano spazialmente irrisolti, la loro struttura e dinamica possono essere sondate attraverso la variabilità delle loro righe di emissione. Le analisi precedenti sulla variabilità spettroscopica dei quasar sono state ampiamente limitate a piccoli campioni di singoli oggetti a causa della mancanza di dati su larga scala nel dominio del tempo. Sebbene la Sloan Digital Sky Survey abbia ampliato la popolazione nota di quasar, le prime iterazioni si sono concentrate su osservazioni a epoca singola per studi cosmologici. Sebbene SDSS IV abbia introdotto la spettroscopia ripetuta attraverso il Progetto di Mappatura della Reverberazione (RM) e la Time Domain Spectroscopic Survey (TDSS), mancava un catalogo completo delle proprietà di variabilità per l'esteso campione di spettri disponibili tramite la Data Release 16 (DR16). Nello specifico, vi era la necessità di caratterizzare la variabilità della riga larga H$\beta$ (bH$\beta$), inclusi i cambiamenti nella velocità radiale, nella luminosità e nella forma del profilo di riga, attraverso una popolazione statisticamente significativa.

Metodologia
Gli autori hanno costruito un catalogo delle proprietà di variabilità spettroscopica per tutti i quasar SDSS con $z < 0,8$ e con almeno due osservazioni incluse nel Catalogo Quasar DR16. Il campione comprende 10.336 quasar unici e 31.320 spettri, risultando in un dataset dominato da coppie di spettri del progetto RM, in particolare per intervalli temporali brevi.

1. Decomposizione Spettrale: Per isolare la componente larga di H$\beta$, gli autori hanno eseguito una decomposizione spettrale su ogni spettro all'interno della finestra nel frame di riposo 3750–5550 Å. Il modello ha fittato simultaneamente un continuo a legge di potenza, un complesso di Fe II ottico (utilizzando il template I Zw 1), un continuo della galassia ospite (utilizzando i template MILES), un continuo di Balmer e varie righe di emissione. Le righe strette (es. [O III], H$\beta$ stretta) e le righe larghe (H$\beta$ larga, H$\gamma$, H$\delta$) sono state modellate utilizzando profili gaussiani. Per mitigare i minimi locali nell'ipotesi $\chi^2$ senza il costo computazionale di un intero MCMC per l'intero campione, gli autori hanno impiegato una procedura di fitting multi-step utilizzando ottimizzatori $\chi^2$ rapidi (BFGS, Nelder, Powell) e ipotesi iniziali intelligenti dei parametri.
2. Misurazione della Variabilità: Per ogni coppia di spettri, gli autori hanno misurato la variazione della velocità radiale del centroide ($\Delta v_{rad}$) tramite cross-correlazione dei profili dei dati isolati di H$\beta$ larga. Altre proprietà misurate includono la variazione della luminosità del continuo a 5100 Å, la luminosità della H$\beta$ larga, la larghezza a metà altezza (FWHM), l'ampiezza equivalente (EW) e l'asimmetria di Pearson. Le proprietà fisiche a epoca singola, come la massa del buco nero ($M_{BH}$) e la luminosità bolometrica ($L_{bol}$), sono state derivate utilizzando relazioni di scala consolidate.
3. Stima dell'Incertezza: Gli autori hanno sviluppato un metodo migliorato per stimare le incertezze in $\Delta v_{rad}$ che tenga conto sia della decomposizione spettrale che della cross-correlazione. Hanno utilizzato un campione pilota di 170 coppie modellate con MCMC (utilizzando il campionatore emcee) per costruire relazioni di scala. Una regressione lineare multipla (MLR) è stata applicata per predire l'incertezza $\sigma(\Delta v_{rad})$ basandosi sul rapporto segnale-rumore della H$\beta$ larga (SNR(bH$\beta$)) e sulla grandezza stessa di $\Delta v_{rad}$. Questo è stato validato contro una stima non parametrica derivata da coppie con baseline $< 2$ giorni, dove la vera variabilità è assunta come trascurabile.

Contributi Chiave

• Creazione del Catalolo: Il lavoro presenta il primo catalogo su larga scala delle proprietà di variabilità della riga larga H$\beta$ per il campione di quasar SDSS DR16, fornendo proprietà misurate per i singoli spettri e differenze per tutte le coppie spettrali.
• Modellazione Migliorata dell'Incertezza: Gli autori propongono un framework robusto per stimare le incertezze in $\Delta v_{rad}$ quando viene utilizzata la decomposizione spettrale, andando oltre i semplici errori di cross-correlazione per includere le incertezze di decomposizione tramite relazioni di scala derivate da MCMC.
• Caratterizzazione delle Forme di Distribuzione: Lo studio fornisce una caratterizzazione statistica della distribuzione di $\Delta v_{rad}$ e di altre metriche di variabilità in funzione degli intervalli di tempo nel frame di riposo.

Risultati

• Distribuzione di $\Delta v_{rad}$: Gli autori riscontrano che $\Delta v_{rad}$ non è distribuito normalmente. Per baseline nel frame di riposo fino a circa 300 giorni, la distribuzione è meglio rappresentata da una distribuzione esponenziale. Per baseline superiori a circa 1.000 giorni, la distribuzione è consistente con una distribuzione Lorentziana (con una semilarghezza a metà altezza $\gamma = 81$ km s$^{-1}$). Ciò contrasta con le assunzioni precedenti di residui gaussiani negli studi di variabilità a lungo termine.
• Bias del Campione e Cadenza: L'analisi evidenzia una differenza significativa tra il sottocampione RM e il campione non-RM. Sebbene il campione RM costituisca una piccola frazione di quasar unici, esso domina il dataset in termini di coppie spettrali. Il campione RM campiona uniformemente un ampio intervallo di masse di buchi neri e luminosità, mentre il campione non-RM mira preferenzialmente a quasar con maggiore luminosità e buchi neri più massicci.
• Asimmetria del FWHM: È stata osservata una significativa asimmetria nella variazione del FWHM ($\Delta$FWHM) per baseline di 300–2.000 giorni, dove una diminuzione del FWHM era più probabile rispetto a un aumento. Gli autori dimostrano che questa asimmetria è guidata quasi interamente da pochi quasar RM che sono transitati da stati deboli/larghi a stati brillanti/stretti su baseline di diversi anni. Rimuovendo questi oggetti specifici, l'asimmetria viene soppressa, illustrando l'impatto della statistica di piccoli numeri nei campioni dominati da RM.
• Proprietà Fisiche: Le masse dei buchi neri derivate ($10^7$–$10^9 M_\odot$) e le luminosità bolometriche ($10^{44}$–$10^{47}$ erg s$^{-1}$) sono consistenti con le misurazioni precedenti di quasar SDSS a basso redshift.

Significatività
Il lavoro fornisce vincoli osservativi essenziali sulla variabilità a livello di popolazione delle righe di emissione larghe dei quasar. Questi vincoli sono critici per:

1. Testare i Modelli BLR: Le distinte forme di distribuzione di $\Delta v_{rad}$ a diverse scale temporali offrono punti di riferimento per i modelli teorici della regione delle righe larghe (BLR), distinguendo tra effetti di tempo di attraversamento della luce ed evoluzione temporale dinamica.
2. Ricerca di Sistemi Binari di Buchi Neri Supermassicci (SMBHB): Lo studio affina la comprensione del "rumore" nelle misurazioni della velocità del centroide per i quasar a singolo buco nero. Stabilendo che i cambiamenti del centroide su tempi brevi non sono gaussiani, gli autori forniscono vincoli empirici più accurati per aiutare a distinguere le firme monotone o sinusoidali di SMBHB dalla variabilità intrinseca dei normali quasar.
3. Analisi Future di SDSS: La metodologia e i dati servono come fondamento per future analisi nel dominio del tempo, in particolare con SDSS V, evidenziando la necessità di considerare attentamente le differenze di cadenza e i bias di selezione tra i campioni RM e non-RM.