Probing jet base emission of M87* with the 2021 Event Horizon Telescope observations

Technische Samenvatting: Het onderzoeken van de emissie van de jet-basis van M87 met de 2021 Event Horizon Telescope-observaties*

Probleemstelling
Eerdere Event Horizon Telescope (EHT) observaties van M87* in 2017 en 2018 hebben de compacte ringstructuur (de schaduw van het zwarte gat) succesvol opgelost bij 230 GHz. Echter, deze datasets leden aan een "missing flux"-probleem (ontbrekende flux), waarbij de totale fluxdichtheid gemeten op triviale (co-gelokaliseerde) baselines (bijv. ALMA/APEX, JCMT/SMA) ongeveer 1,0–1,4 Jy was, aanzienlijk hoger dan de ~0,5–1,0 Jy die door compacte ring-alleen modellen werd teruggevonden. Deze discrepantie suggereerde de aanwezigheid van uitgebreide emissie op grote ruimtelijke schalen, waarschijnlijk geassocieerd met de basis van de jet. De 2017 en 2018 arrays beschikten niet over de noodzakelijke korte en intermediaire baselines om deze uitgebreide structuur te beperken of te imageren, wat leidde tot onopgeloste degeneraties in de reconstructie van de jet-basis.

Methodologie
Deze studie maakt gebruik van de aanzienlijk verbeterde $(u, v)$-dekking van de 2021 EHT-observatiecampagne, specifiek de integratie van de 12-m Kitt Peak (KP) telescoop (VS) en de NOEMA-array (Frankrijk). Deze nieuwe stations introduceerden twee cruciale intermediaire baselines van gemiddelde lengte:

1. PV-NOEMA: 1100 km, gevoelig voor ruimtelijke schalen van ~250 $\mu$as (0,02 pc).
2. SMT-KP: 100 km, gevoelig voor ruimtelijke schalen van ~2500 $\mu$as (0,2 pc).

Vanwege de schaarste van deze intermediaire baselines is directe imaging van de zwakke, uitgebreide jet-basis niet robuust. In plaats daarvan hanteert de studie een geometrische model-fitting benadering met behulp van closure phases (sluitingsfasen). Closure phases zijn robuust tegen station-gebaseerde gain-kalibratiefouten en zijn gevoelig voor asymmetrieën in de helderheidsverdeling. De analyse richt zich op closure-triangles gevormd met de zeer gevoelige ALMA-station (ALMA-PV-NOEMA en ALMA-SMT-KP).

De auteurs vergelijken geobserveerde closure phases met modellen bestaande uit:

• Een compacte ringcomponent (gereconstrueerd via DoG-HIT en andere algoritmen).
• Een extra, zwakke, offset Gaussische component die de emissie van de jet-basis representeert.

De modelparameters (positie, fluxdichtheid en full width at half maximum) worden geoptimaliseerd met behulp van Markov chain Monte Carlo (MCMC) sampling om de gereduceerde $\chi^2_\nu$ residuen tussen het model en de geobserveerde closure phases te minimaliseren. Alternatieve modellen, inclusief asymmetrische Gaussians, dubbele Gaussians en GMVA-afgeleide jet-structuren, worden ook getest.

Belangrijkste Resultaten

• Detectie van Residuen: De compacte ring-alleen modellen slagen er niet in om de geobserveerde closure phases op de ALMA-PV-NOEMA en ALMA-SMT-KP triangles volledig te reproduceren, waarbij significante residuen (een "bump" en een "offset" van enkele graden) optreden.
• Gaussische Component Fit: Het introduceren van een enkele, offset Gaussische component verbetert de fit aanzienlijk. Het best-fit model plaatst een zwakke Gaussische component op een geprojecteerde scheiding van $\approx 5500$ AU (gecentreerd bij $\Delta$R.A. $\approx 320$ $\mu$as, $\Delta$Dec. $\approx 60$ $\mu$as, relatief aan het ringcentrum) met een geschatte fluxdichtheid van ~60 mJy en een grootte (FWHM) van ~180 $\mu$as.
• Statistische Verbetering: De inclusie van deze component vermindert de $\chi^2_\nu$ voor de ALMA-PV-NOEMA triangle van 2,73 naar 1,29 en voor de ALMA-SMT-KP triangle van 2,52 naar 1,12. De globale fit over alle triangles verbetert van 1,50 naar 1,39.
• Modelbeperkingen: Complexere modellen (asymmetrische Gaussians of een tweede Gaussische component) bieden geen statistisch significante verbetering ten opzijde van het enkelvoudige Gaussische model. De data bevoordelen geen helderdere emissie op deze schalen; een aanzienlijk helderdere component zou residuen produceren die groter zijn dan geobserveerd.
• Fysieke Consistentie: De locatie van de teruggevonden Gaussian komt overeen met de zuidelijke/zuidwestelijke arm van de edge-brightened jet geobserveerd in 86 GHz GMVA-beelden, wat consistent is met de verwachte richting van de grootschalige jet.
• Impact op Vorige Data: De studie bevestigt dat deze zwakke component ondetecteerbaar zou zijn geweest in de 2017 en 2018 data vanwege het gebrek aan intermediaire baselines, wat de robuustheid van de eerdere ring-alleen reconstructies voor die tijdperken valideert.

Betekenis en Claims
Het artikel claimt de eerste beperkingen op 230 GHz te bieden voor zwakke, asymmetrische emissie op intermediaire schalen nabij M87*. De primaire bijdrage is de identificatie van een plausibele jet-basis component die een deel van de voorheen onverklaarde "missing flux" verklaart.

De auteurs behouden echter een bescheiden en conservatieve houding met betrekking tot de definitieve aard van deze detectie:

1. Morfologische Beperkingen: Vanwege het beperkte aantal intermediaire baselines (slechts twee), kunnen de auteurs de gedetailleerde morfologie van de jet-basis niet robuust reconstrueren. De teruggevonden Gaussian wordt behandeld als een bovengrens (upper limit) op de gereconstrueerde jet-basis fluxdichtheid bij deze specifieke ruimtelijke schalen.
2. Systematische Onzekerheden: Hoewel polarisatie-lekken en andere bekende systematieken zijn uitgesloten als de bron van het signaal, erkennen de auteurs dat onbekende systematische fouten niet volledig uitgesloten kunnen worden, gezien de afhankelijkheid van single-baseline/kleine-triangle metingen.
3. Toekomstige Vereisten: Het artikel concludeert dat definitieve bevestiging en hogere-precisie beperkingen op de structuur en aard van de jet-basis toekomstige EHT-observaties zullen vereisen met een uitgebreidere intermediaire-baseline dekking (meer stations) en potentieel uitgebreide frequentiebereiken.

Samenvattend toont het werk aan dat de 2021 EHT-data bewijs levert voor een zwakke, naar het zuidwesten verschoven emissiecomponent die consistent is met de jet-basis, maar de huidige datadichtheid beperkt het vermogen om de structuur ervan volledig op te lossen, wat verdere observaties noodzakelijk maakt voor een volledige fysieke karakterisering.