Primary Beam Chromaticity in HIRAX: I. Characterization from Simulations and Power Spectrum Implications

Sintesi Tecnica: Cromaticità del Fascio Primario in HIRAX: I. Caratterizzazione da Simulazioni e Implicazioni sullo Spettro di Potenza

Enunciato del Problema
L'esperimento Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) mira a vincolare l'energia oscura mappando le Oscillazioni Acustiche dei Barioni (BAO) tramite la mappatura dell'intensità a 21 cm. Una sfida critica in questo sforzo è la presenza di foreground (primari), sia galattici che extragalattici, che sono di ordini di grandezza più luminosi del segnale cosmologico. Sebbene i foreground siano spettralmente regolari, la risposta dipendente dalla frequenza (cromatica) del fascio del radiotelescopio causa la "fuga" di questi foreground regolari nella finestra del segnale pulito nello spazio di Fourier. Questa fuga è esacerbata dalle strutture spaziali e spettrali complesse dei lobi secondari e del lobo principale. Studi precedenti hanno modellato i fasci fino a due o tre lobi secondari, ma una contaminazione significativa può derivare da sorgenti in lobi secondari più lontani, e i formalismi esistenti spesso mancano della capacità di caratterizzare la cromaticità in interferometri a campo largo che si estendono profondamente nei lobi secondari.

Metodologia
Gli autori presentano una tecnica di modellazione del fascio motivata fisicamente, progettata per rappresentare analiticamente l'intero fascio Stokes I di HIRAX, estendendosi oltre il lobo principale.

1. Simulazioni Elettromagnetiche: Lo studio utilizza simulazioni elettromagnetiche (EM) eseguite nel CST Studio Suite. Il modello semplifica la configurazione disco-feed di HIRAX (rimuovendo gambe e cavi per ridurre il costo computazionale) pur mantenendo il diametro effettivo di circa 4 anni e un rapporto focale f/0.21. Le simulazioni coprono la larghezza di banda 400–800 MHz con 81 canali.
2. Base Analitica (Trasformate di Zernike): Per modellare il fascio, gli autori impiegano una base flessibile derivata dalla trasformata di Fourier dei polinomi di Zernike, denominata "Trasformate di Zernike" (ZT). A differenza della modellazione tradizionale nello spazio delle immagini che si concentra sul lobo principale, questa base permette di separare le strutture radiali e azimutali. Il modello è definito come una somma di funzioni di base reali $X^m_n(\theta, \phi)$, dove $n$ rappresenta l'ordine radiale e $m$ l'ordine azimutale.
3. Procedura di Fitting: I fasci simulati vengono adattati canale per canale utilizzando una routine di ottimizzazione numerica. Gli autori utilizzano un approccio a minimi quadrati pesati, dando maggior peso ai lobi secondari per evitare l'over-fitting sul lobo principale. Il fitting determina i coefficienti $\alpha^m_n$ e un parametro di scala dipendente dalla frequenza $\theta_0$.
4. Analisi dello Spettro di Potenza: Per quantificare l'impatto degli errori di modellazione del fascio, gli autori simulano visibilità mock per un array HIRAX a 256 elementi utilizzando il Global Sky Model (GSM) e un campo casuale gaussiano per il segnale a 21 cm. Effettuano la sottrazione del modello utilizzando varie approssimazioni del fascio (Gaussiana ingenua, Gaussiana adattata e modelli ZT con diverse troncature radiali e azimutali) per recuperare lo spettro di potenza a 21 cm. La fedeltà del recupero è valutata utilizzando una "funzione di trasferimento dello spettro di potenza ingenua" ($T$), definita come il rapporto tra lo spettro di potenza del segnale recuperato e lo spettro di potenza del segnale vero.

Contributi Chiave e Risultati

• Caratterizzazione del Fascio: Lo studio caratterizza la cromaticità del lobo principale e dei lobi secondari di HIRAX. Il lobo principale presenta una struttura a "ripple del fascio" (variazione sinusoidale in larghezza del fascio e ampiezza) causata dalle riflessioni tra il feed e la superficie del disco, che introduce una struttura spettrale che si mescola con i foreground.
• Efficacia della Modellazione: Il modello delle Trasformate di Zernike cattura con successo le complesse strutture dei lobi secondari. Gli autori riscontrano che la modellazione fino all'ordine azimutale ottupolare (variazione angolare del quarto ordine, $m_{max}=4$) è sufficiente per minimizzare la fuga al di fuori del wedge con un bias minimo.
• Impatto sullo Spettro di Potenza:
  • Modelli solo del lobo principale: I modelli che ignorano i lobi secondari o che si affidano a semplici adattamenti gaussiani (anche se catturano la cromaticità del lobo principale) producono una significativa fuga di foreground e strutture a ripple residue nello spettro di potenza, particolarmente alle scale $k \sim 0.25$ e $0.5 \, \text{Mpc}^{-1}$.
  • Modi Radiali vs. Azimutali: Aumentare l'ordine radiale ($n$) è essenziale per catturare l'estensione dei lobi secondari. Tuttavia, lo studio indica che catturare fino a $m_{max}=2$ (modi quadrupolari) è ampiamente sufficiente per evitare problemi maggiori riguardanti l'asimmetria azimutale, sebbene ordini superiori contribuiscano al fitting.
  • Mitigazione della Fuga: Man mano che la complessità del modello del fascio aumenta (avvicinandosi al fascio simulato "ground truth"), la fuga di potenza dei foreground nella finestra del segnale pulito diminuisce e la funzione di trasferimento si avvicina all'unità.
• Limitazioni: Gli autori notano che le attuali simulazioni assumono visibilità prive di rumore, ridondanza perfetta e risposte del fascio uniformi. Dichiarano esplicitamente che gli effetti di cross-coupling, la fuga di polarizzazione (Stokes Q, U, V) e l'impatto di piatti strettamente impacchettati sulla forma del fascio non sono affrontati in questo lavoro.

Significatività e Rivendicazioni
Il documento rivendica di essere il primo di una serie in due parti riguardante la cromaticità del fascio primario in HIRAX. La sua significatività primaria risiede nell'aver stabilito un formalismo per modellare l'intero fascio, inclusi i lobi secondari lontani e la loro cromaticità, il che è vitale per la mappatura dell'intensità a 21 cm ad alta precisione.

Gli autori rivendicano con modestia che il loro lavoro:

1. Fornisce una base analitica necessaria per comprendere il comportamento del fascio che si estende profondamente nei lobi secondari, una lacuna nella letteratura precedente sui 21 cm.
2. Dimostra che trascurare la cromaticità dei lobi secondari e le asimmetrie azimutali porta a una significativa fuga di foreground e bias nella stima dello spettro di potenza.
3. Stabilisce obiettivi specifici per le future campagne di mappatura del fascio tramite droni, suggerendo che la modellazione fino all'ordine azimutale ottupolare e la cattura di una sostanziale estensione radiale siano necessari per minimizzare gli errori sistematici.
4. Pone le basi per un secondo articolo che applicherà questo framework di modellazione a dati reali del fascio mappati da droni e proporrà specifiche strategie di mitigazione dei foreground.

Il documento non pretende di aver risolto interamente il problema della sottrazione dei foreground, ma identifica la modellazione del fascio come un sistema critico che deve essere controllato per raggiungere gli obiettivi cosmologici del sondaggio. Enfatizza che, sebbene le funzioni di trasferimento possano teoricamente debiasare i risultati, la caratterizzazione accurata della risposta del fascio rimane il prerequisito per qualsiasi strategia di pulizia o evitamento dei foreground di successo.