A Census of Variable Radio Sources at $3\,$GHz

Technische Samenvatting: Een Cijfer van Variabele Radiobronnen bij 3 GHz

Probleem en Motivatie
Radiovariabiliteit is een alomtegenwoordig verschijnsel over diverse astronomische tijdschalen, variërend van milliseconden bij pulsars tot decennia bij Actieve Galactische Kernen (AGN). Hoewel gerichte monitoring van bekende variabelen en smalle-veldopnames inzichten hebben verschaft, blijft de bredere populatie van variabele radiobronnen slecht gekarakteriseerd vanwege beperkte bemonstering. Eerdere grote oppervlaktesurveys, zoals FIRST en NVSS, identificeerden variabiliteit maar waren vaak beperkt door gevoeligheid of cadans. Bij lagere frequenties is variabiliteit zeldzaam, maar bij hogere frequenties (boven een paar GHz) wordt het algemener, toch ontbreekt er een bijna-hemeldekende telling bij 3 GHz. Dit werk adresseert de noodzaak om het tempo, de schaal en de aard van radiovariabiliteit over de hele hemel te karakteriseren om dynamische processen zoals accretiemechanica, jetvoortplanting en sterrenontwikkeling beter te begrijpen.

Methodologie
De auteurs gebruiken data uit de eerste twee epocha van de Karl G. Jansky Very Large Array Sky Survey (VLASS), een meer-epocha survey die werkt bij $\nu \sim 3$ GHz. De analyse steunt op de VLASS Quick Look-afbeeldingen en catalogi, die, hoewel onderhevig aan bekende kwaliteitsbeperkingen (zoals fouten in fluxdichtheidsschaling en artefacten), het volledige hemelgebied ten noorden van $-40^\circ$ declinatie bestrijken.

1. Dataselectie: De studie richt zich op compacte bronnen die in zowel Epoch 1 (2017–2019) als Epoch 2 (2020–2022) zijn gedetecteerd, gescheiden door een mediaan tijd van 974 dagen. Bronnen werden gefilterd op kwaliteitsvlaggen die echte detecties aangeven en op de status opgelost/onopgelost, gebaseerd op de verhouding van geïntegreerde tot piekfluxdichtheid ($S_{total}/S_{peak}$) versus signaal-ruisverhouding.
2. Systematische Correctie: Om rekening te houden met niet-homogene fluxkalibratie over het survey-oppervlak, analyseerden de auteurs de verhouding van mediaan piekfluxdichtheden tussen epocha per tegel. Ze identificeerden een bimodaliteit in fluxverhoudingen afhankelijk van het sub-epoch van observatie en corrigeerden de Epoch 2 fluxdichtheden om een mediaan verhouding van één binnen elke tegel te waarborgen.
3. Identificatie van Variabiliteit: Kandidaat-variabelen werden geïdentificeerd met behulp van twee parameters: de Variabiliteitsstatistiek ($V_s$), die de absolute fluxverandering ten opzichte van meetonzekerheid vertegenwoordigt, en de modulatie-index ($m$), die de fractieverandering vertegenwoordigt. Bronnen werden geselecteerd als $|V_s| > 4\sigma_{V_s, Tegels}$ en $|m| > 0.26$.
4. Kwaliteitsborging: Een visuele inspectie van 4.124 kandidaten werd uitgevoerd met een drie-panelen afbeeldingsformaat (Epoch 1, Epoch 2 en verschilafbeelding) met geconvolueerde resolutie om variaties in de Point Spread Function (PSF) te mitigeren. Een consensusvlag ($f_{vi}$) werd toegewezen op basis van onafhankelijke classificaties door drie auteurs. De uiteindelijke robuuste steekproef bestaat uit 3.618 bronnen waarbij alle classificeerders het erover eens waren dat de bron een echte variabele was.
5. Multigolfbandkarakterisering: De variabele steekproef werd gekruist met SDSS, LS DR8 (voor roodverschuivingen), AllWISE (voor infraroodkleuren) en de 4FGL Fermi-LAT catalogus (voor $\gamma$-straal tegenhangers) om de aard van de variabele bronnen te bepalen.

Belangrijkste Resultaten

• Variabiliteitspercentage: Ongeveer 3.600 compacte bronnen bleken significant te variëren. Voor bronnen met een gemiddelde fluxdichtheid ($\mu_S$) helderder dan 20 mJy vertonen 4,9% variabiliteit ($|m| > 0.26$). Dit percentage neemt toe met helderheid, tot 9,0% voor bronnen met $\mu_S > 300$ mJy.
• Fluxveranderingen: De studie onderscheidt tussen flaring en fading bronnen. Hoewel het aantal flaring en fading bronnen ongeveer gelijk is, vertonen flaring bronnen iets grotere absolute fluxveranderingen (mediaan $|\Delta S| = 15.6$ mJy) vergeleken met fading bronnen (mediaan $|\Delta S| = 13.4$ mJy), een verschil dat statistisch significant bleek ($p = 4 \times 10^{-4}$).
• Aard van de Bron: Multigolfbandanalyse onthult dat de meerderheid van de variabele radiobronnen extragalactisch is.
  • Blazars: Infraroodkleuranalyse geeft aan dat blazars significant oververtegenwoordigd zijn onder variabele bronnen in vergelijking met de algemene radiopopulatie. Ongeveer 68% van de variabelen met WISE-kleuren is consistent met blazars (56% BL Lac's, 40% FSRQ's), vergeleken met 42% in de oudersteekproef.
  • Galactische Bronnen: De grootste fractieveranderingen in helderheid worden vertoond door galactische bronnen, specifiek radiosterren (bijv. RS CVn-systemen) en microquasars.
  • $\gamma$-straal Correlatie: Variabele radiobronnen hebben een grotere kans op een $\gamma$-straal tegenhanger in de 4FGL catalogus (8%) dan de algemene radiopopulatie (2%). Bovendien vertonen deze tegenhangers hogere $\gamma$-straal variabiliteitsindices, wat een correlatie suggereert tussen radio- en hoog-energetische variabiliteit.
• Roodverschuivingsverdeling: De roodverschuivingsverdeling van variabele bronnen weerspiegelt die van de niet-variabele oudersteekproef, wat suggereert dat variabelen bestaan gedurende dezelfde epocha van de kosmische geschiedenis. Variabelen met $\gamma$-straaldetecties hebben echter de neiging om op lagere roodverschuivingen te zitten ($z \sim 0.64$) dan die zonder ($z \sim 1.22$).
• Extreme Voorbeelden: Het paper identificeert specifieke extreme gevallen, waaronder LB 2449 (een BL Lac met een vijfvoudige helderheidstoename) en PKS 1502+106 (een Flat Spectrum Radio Quasar met een bijna 50% daling). Galactische voorbeelden zijn Cygnus X-3 (een microquasar die een 5,4-voudige toename vertoont) en $\sigma$ Gem (een RS CVn-ster die een grote fractiedaling vertoont).

Betekenis en Beweringen
De auteurs beweren dat dit werk de eerste bijna-hemeldekende telling van variabele radiobronnen bij 3 GHz biedt, waarmee een basislijn wordt gelegd voor de populatie van variabelen op tijdschalen van $\sim$2,5 jaar. De studie toont aan dat:

1. Blazars domineren de variabele populatie: De oververtegenwoordiging van blazars suggereert dat diepe radio-tijdsdomeinsurveys effectieve hulpmiddelen zijn voor het identificeren van zwakke, hoog-roodverschuivende blazars die mogelijk worden gemist door huidige $\gamma$-straalsurveys.
2. Methodologische robuustheid: Ondanks het gebruik van Quick Look-afbeeldingen met bekende systematische fouten, isoleerde de statistische aanpak per tegel en de visuele inspectie succesvol echte variabelen, waarbij vervuiling door valse detecties minimaal leek.
3. Galactische bruikbaarheid: De survey is effectief in het identificeren van galactische radiosterren en microquasars, met name die welke extreme fractievariabiliteit vertonen.
4. Toekomstpotentieel: De auteurs merken op dat de toevoeging van een derde VLASS-epocha statistiek met drie punten voor variabiliteit zal mogelijk maken, wat de detectie van stochastische variabelen verbetert en bronnen met consistente langetermijntrends identificeert (bijv. nieuw ingeschakelde jets). Ze benadrukken ook het potentieel om "quasi-stabiele" zwakke puntbronnen te identificeren voor gebruik als fasekalibratoren voor radio-telescopen van de volgende generatie zoals de ngVLA en SKA.

Het paper concludeert dat radiovariabiliteit-surveys een krachtige, complementaire methode zijn voor het bestuderen van AGN-fysica en sterrenontwikkeling, en in staat zijn populaties van objecten bloot te leggen die zwak zijn of zich op hoge roodverschuivingen bevinden ten opzichte van die gevonden in andere golfbanden.