Study of an interaction between the jet and an interstellar medium of M87 with a spectral analysis by using Chandra

Technische Samenvatting: Studie naar een interactie tussen de jet en een interstellaire materie van M87 met een spectrale analyse met behulp van Chandra

Probleemstelling
Het M87-radiogalactisch centrum, gelegen in het centrum van het Virgocluster, vertoont een prominente jet die zich uitstrekt over 20 boogseconden. Hoewel de oorsprong van de TeV-gamma-straling in M87 over het algemeen wordt toegeschreven aan de kern of specifieke knopen zoals HST-1, blijven de fysieke mechanismen die de interactie tussen de jet en het omringende interstellaire medium (ISM) drijven, onderwerp van onderzoek. Eerdere spectrale analyses van de M87-jet met Chandra-data (bijv. Wilson & Yang 2002; Perlman & Wilson 2005; Sun et al. 2018) pasten voornamelijk X-ray energiespectra aan met power-law modellen, uitgaande van synchrotronemissie. Deze studies stuitten echter op beperkingen:

1. Onvoldoende statistiek: Veel analyses vertrouwden op korte blootstellingstijden (bijv. 13–38 ks), wat leidde tot slechte statistiek waarbij complexe modellen acceptabel leken te zijn simpelweg door de lage datakwaliteit.
2. Fouten in de achtergrondsubtractie: Eerdere werken gebruikten vaak een gemeenschappelijk achtergrondgebied voor de gehele jet. Dit is fysiek incorrect omdat het hete gas van het Virgocluster, dat bijdraagt aan de achtergrond, in intensiteit varieert met de afstand tot de kern.
3. Verwaarlozing van thermische componenten: Er is een mogelijkheid dat de jet interageert met het ISM via schokken, waarbij gas wordt verhit om een thermische component (bremsstrahlung) te produceren naast de niet-thermische synchrotronemissie. Deze component kan in studies met lage statistiek over het hoofd zijn gezien.

Methodologie
Om deze problemen aan te pakken, maakt deze studie gebruik van een aanzienlijk grotere dataset van de Chandra X-ray Observatory, die observaties beslaat van 2000 tot 2018 met een totale blootstellingstijd van ongeveer 800 ks.

• Dataselectie en Pile-up Behandeling: De analyse richt zich op vier specifieke regio's: de kern, knoop HST-1, knoop D en knoop A. Vanwege de hoge helderheid van de kern en HST-1 werden de data zorgvuldig geselecteerd om "pile-up" te beheersen (waarbij meerdere fotonen een enkel pixel raken, wat het spectrum vertekent). De data werden gefilterd om ervoor te zorgen dat de pile-up fracties onder de 5% bleven voor de meeste analyses, hoewel 10% werd geaccepteerd voor de kern en HST-1 in de eerdere dataset om de statistiek te maximaliseren.
• Achtergrondmodellering: In tegen tegenstelling tot eerdere studies werden er voor elke knoop individuele achtergrondregio's gekozen (ten noorden en zuiden van de bron) om rekening te houden met de ruimtelijke variatie van het hete gas van het Virgocluster. Voor knoop A werd een specifieke achtergrondregio gekozen om bekende zachte X-ray dips te vermijden.
• Spectrale Analyse: Energiespectra werden geëxtraheerd met CIAO-software (v4.15) en geanalyseerd met XSPEC. Het energiebereik werd ingesteld van 0,5 keV tot 10 keV om kalibratieproblemen bij lagere energieën te vermijden.
• Model Fitting:
  1. Een eenvoudig power-law model (dat synchrotronemissie vertegenwoordigt) met foto-elektrische absorptie werd getest.
  2. Indien het power-law model statistisch werd afgewezen (geïndiceerd door een hoge $\chi^2$ en een soft excess), werd een gecombineerd model van een power law en een APEC-model (thermische bremsstrahlung) toegepast.
  3. Fysische parameters, waaronder neutrale dichtheid, plasma-dichtheid en versnelde elektrondichtheid, werden afgeleid uit de best-fit modellen.
  4. De flux van niet-thermische bremsstrahlung, voortkomend uit de interactie van versnelde elektronen met het ISM, werd berekend en vergeleken met Fermi gamma-ray observaties.

Belangrijkste Resultaten

1. Spectrale Kenmerken:

   • Kern en Knoop D: De X-ray energiespectra voor zowel de kern als knoop D worden goed beschreven door een eenvoudig power-law model, consistent met synchrotronemissie.
   • Knoop HST-1 en Knoop A: Voor deze regio's wordt een zuiver power-law model statistisch afgewezen.
     • HST-1: In de gecombineerde dataset (inclusief 10% pile-up data) levert een power-law model een kanskans van $2,6 \times 10^{-4}$ op. Het toevoegen van een thermische component (APEC) zorgt voor een redelijke fit.
     • Knoop A: Het power-law model wordt sterk afgewezen ($\chi^2 = 628,76$ voor 355 d.o.f, kans $1,6 \times 10^{-17}$). Een gecombineerd model van een power law en thermische bremsstrahlung past goed bij de data.

2. Eigenschappen van de Thermische Component:

   • De thermische component gedetecteerd in de knopen HST-1 en A wordt gekenmerkt door een temperatuur van 0,2 keV en een metaalovervloed van 0,00 (vrijwel nul).
   • De metaalovervloed is aanzienlijk lager dan die van het omringende hete gas van het Virgocluster (die ongeveer ~1 zonne abundantie is), wat suggereert dat het verhitte gas binnen de jet zelf ontstaat in plaats van in het omringende medium.
   • De thermische fluxratio in knoop A is ongeveer 22% van de totale flux.

3. Fysische Parameters:

   • Plasma Dichtheid: Afgeleid van de APEC normalisatie, is de elektrondichtheid in knoop A ongeveer $11,0 \, \text{cm}^{-3}$.
   • Schokverhitting: De aanwezigheid van een thermische component met een lage metalliciteit en hoge temperatuur relatief aan het omringende medium ondersteunt de hypothese dat gas binnen de jet wordt verhit door een schokgolf die voortvloeit uit de interactie met het ISM.
   • Neutrale Dichtheid: Berekend uit het verschil tussen de X-ray afgeleide kolomdichtheid en 21 cm radio-observaties, wat wijst op de aanwezigheid van koud materiaal rond de jet.

4. Niet-Thermische Bremsstrahlung en Fermi Vergelijking:

   • De studie heeft de verwachte flux van niet-thermische bremsstrahlung berekend van versnelde elektronen die interageren met het ISM.
   • De berekende flux van de kern en knoop HST-1 is orders van grootte lager dan de geobserveerde GeV-flux van Fermi.
   • Conclusie: Er is geen bijdrage van niet-thermische bremsstrahlung in de kern of HST-1 aan de geobserveerde Fermi-flux. De paper merkt op dat knoop A mogelijk een bijdrage levert als het magnetische veld daar aanzienlijk lager is dan aangenomen (<10 $\mu$G), maar dit vereist verder onderzoek.

Betekenis en Claims
De paper claimt dat door gebruik te maken van hoog-statistische data (800 ks) en correcte, ruimtelijk afhankelijke achtergrondsubtractie, het erin is geslaagd een thermische component te identificeren in de X-ray spectra van de M87-jet knopen HST-1 en A. Deze bevinding daagt eerdere interpretaties uit die uitsluitend vertrouwden op power-law modellen.

De auteurs concluderen dat:

1. De thermische emissie waarschijnlijk gas binnen de jet is dat wordt verhit door een schokinteractie met het ISM, wat wordt aangetoond door de lage metaalovervloed en hoge temperatuur.
2. Het mechanisme van niet-thermische bremsstrahlung (elektron-ISM interactie) de geobserveerde GeV gamma-straalflux van M87 niet verklaart, wat de conclusie versterkt dat de oorsprong van de TeV/GeV emissie waarschijnlijk elders ligt (bijv. de kern of andere deeltjesversnellingsmechanismen).
3. Eerdere studies deze thermische kenmerken mogelijk hebben gemist door onvoldoende statistiek en onjuiste achtergrondmodellering.

De studie stelt geen nieuwe experimentele opstellingen voor, maar benadrukt de noodzaak van kwalitatieve analyse van archiefdata om het multi-wavelength energiespectrum en de fysische condities van relativistische jets op te lossen.