The Shape of the Vertical Action Distribution Locates the Scatterers that Heat the Galactic Disc

Technische Samenvatting: De vorm van de verticale actieverdeling lokaliseert de verstrooiers die de Galactische schijf verhitten

Probleemstelling
De verticale snelheidsdispersie van sterren in de Melkweg neemt toe met de leeftijd, een fenomeen dat bekend staat als schijfverhitting (disc heating). Hoewel het bestaan van deze verhitting goed gevestigd is, blijft de identiteit en ruimtelijke verdeling van de "verstrooiers" (scatterers) die hiervoor verantwoordelijk zijn, onopgelost. Er bestaan twee primaire klassen modellen:

1. Volume-vullende verstrooiers: Perturbaties zoals spiraalstructuren, de Galactische staaf of donkere materie-subhalo's die op een ster inwerken gedurende haar gehele baan.
2. In het midvlak geconcentreerde verstrooiers: Perturbaties, zoals reusachtige moleculaire wolken (GMC's), die geconcentreerd zijn in een dunne laag rond het Galactische midvlak en alleen interageren met sterren tijdens hun verticale kruisingen.

Eerdere pogingen om tussen deze mechanismen te onderscheiden met behulp van de leeftijd–snelheidsdispersie-relatie of de vorm van de snelheidsellipsoïde waren inconclusief vanwege degeneraties tussen de abundantie, massa en de geschiedenis van de verhitting. Het artikel betoogt dat de verticale actie ($J_z$), een adiabatische invariant, een superieure diagnostiek biedt omdat het de baan labelt in plaats van de instantane snelheid, en ongevoelig is voor de groei van het potentiaal van de schijf in de loop van de tijd.

Methodologie
De auteurs leiden de theoretische distributie van de verticale actie af voor een coëvale populatie sterren die verhit zijn door verstrooiers met een willekeurige verticale hoogteverdeling.

• Theoretische Afleiding:

  • Het verhittingsproces wordt gemodelleerd als een diffusie in de verticale actieruimte met behulp van de Fokker–Planck-vergelijking.
  • De auteurs leiden de diffusiviteit $D(J_z)$ af als een functie van het hoogteprofiel $q(z)$ van de verstrooiers. Ze tonen aan dat de diffusiviteit schaalt als een machtswet $D \propto J_z^b$, waarbij de exponent $b$ uitsluitend wordt bepaald door de geometrie van de verstrooiers.
  • Volume-vullend geval: Als verstrooiers het volume vullen, groeit de diffusiviteit lineair met de actie ($b=1$), wat resulteert in een exponentiële distributie $p(J_z) \propto \exp(-J_z/J_0)$.
  • Midvlak-laag geval: Als verstrooiers geconcentreerd zijn in een dunne laag, neemt de interactietijd af naarmate de kruisingssnelheid van de ster toeneemt. Dit vermindert de groei van de diffusiviteit tot $D \propto J_z^{1/2}$ ($b=1/2$), wat resulteert in een "gestrekte exponentiële" distributie $p(J_z) \propto \exp[-(J_z/J_0)^{3/2}]$.
  • De algemene oplossing voor een coëvale populatie is $p(J_z) \propto \exp[-(J_z/J_0)^{2-b}]$, waarbij de vorm (staartexponent) enkel afhangt van de geometrie ($b$), terwijl de schaal ($J_0$) de amplitude en geschiedenis van de verhitting codeert.

• Data en Fitting:

  • Het model wordt gefit aan een steekproef van 7.589 low-$\alpha$ rode clump sterren uit de APOGEE DR14 survey, gekruist met Gaia DR2 astrometrie.
  • De steekproef is beperkt tot sterren met effectieve radii $5 \le R_{\text{eff}} \le 10$ kpc en leeftijden $2 \le \tau \le 8$ Gyr.
  • Een cruciaal onderdeel van de methodologie is de rigoureuze behandeling van de selectiefunctie, die sterren met een hoge actie bij de magnitudegrens van de survey bevooroordeeld mist. De auteurs passen een per-ster selectiefunctie, afgeleid van Ting & Rix (2019), direct toe op de model-likelihood.
  • De fittingprocedure marginaliseert over meetfouten in leeftijd, actie en straal met behulp van Gaussische kwadratuur. Cruciaal is dat de schaalparameter $J_0$ bij elke leeftijd vrij wordt gelaten (geparameteriseerd door vier knopen en een radiale helling). Dit maakt het mogelijk om de vormexponent $b$ te meten zonder een specifieke geschiedenis van verhitting of geboortedispersie aan te nemen, waardoor de geometrie effectief wordt ontkoppeld van de verhittingssnelheid.

Belangrijkste Resultaten

1. Gemeten Exponent: De fit levert een diffusiviteitshelling op van $b = 0.51^{+0.06}_{-0.7}$.

   • Deze waarde is consistent met de voorspelling voor een in het midvlak geconcentreerde laag ($b=0.5$), maar is sterk inconsistent met de volume-vullende voorspelling ($b=1$).
   • Synthetische tests bevestigen dat de methode de correcte exponent onbevooroordeeld kan terugvinden en de twee geometrieën scheidt met een significantie van $5.6\sigma$.
   • Het resultaat is robuust over verschillende radiale bins en leeftijdscategorieën, en is bestand tegen variaties in de aangenomen actiefouten.

2. Effectieve Massa van de Verstrooiers:

   • Gebruikmakend van de gemeten vrije schalen (en de afgeleide exponent $b \approx 0.5$), berekenen de auteurs de verhittingsamplitude.
   • Voor de leeftijdscategorie 2–4 Gyr, bij het vergelijken van de verhittingsamplitude met de huidige moleculaire oppervlaktedichtheid, wordt een effectieve massa van de verstrooiers berekend van $2.7 \times 10^6 M_\odot$.
   • Deze waarde valt binnen het bereik ($1.4\text{--}3.4 \times 10^6 M_\odot$) dat is afgeleid van onafhankelijke wolkencatalogi en massafuncties (bijv. Miville-Deschênes et al. 2017).
   • Voor oudere sterren neemt de schijnbare effectieve massa toe wanneer deze wordt vergeleken met de huidige gasdichtheid. Echter, wanneer gecorrigeerd wordt voor de afnemende sterformatiegeschiedenis (en dus de afnemende gasdichtheid) over de afgelopen vele gigajaren, convergeert de effectieve massa voor alle leeftijdscategorieën naar het bereik $1.9\text{--}2.8 \times 10^6 M_\odot$.

3. Verhittingsgeschiedenis:

   • De vrije schalen groeien met de leeftijd als $J_0 \propto \tau^{1.12}$. Dit is sneller dan de $\tau^{0.67}$ die verwacht wordt bij een constante verhittingssnelheid met een koude geboorte, wat aangeeft dat de schijf in het verleden sterker is verhit.

Betekenis en Claims
Het artikel beweert de langlopende discussie over de locatie van de verstrooiers in de Galactische schijf te hebben opgelost door gebruik te maken van de vorm van de verticale actieverdeling in plaats van enkel het gemiddelde of de snelheidsdispersie.

• Lokalisatie van de Verstrooiers: De gemeten exponent $b \approx 0.5$ levert direct bewijs dat het primaire mechanisme voor verticale verhitting gelokaliseerd is in het Galactische midvlak. Dit sluit volume-vullende mechanismen (zoals donkere materie-subhalo's, buigingsgolven of directe spiraalperturbaties) uit als de dominante bron van diffuse verhitting. De auteurs schatten dat een eventuele volume-vullende component minder dan 7% aan de verhitting bijdraagt.
• Identificatie van de Verstrooiers: De combinatie van de midvlak-geometrie en de afgeleide effectieve massa ($\sim 10^6 M_\odot$) identificeert reusachtige moleculaire wolken (GMC's) als de verstrooiers. De massa is consistent met de bovenkant van de GMC-massafunctie, en de geometrie komt overeen met de dunne laag moleculair gas.
• Resolutie van Eerdere Discrepanties: Het artikel verklaart waarom eerdere studies (bijv. Binney & Lacey 1988) $b=1$ (volume-vullend) afleidden. Die studies analyseerden gemengde populaties of snelheidsdispersies zonder de vorm bij een vaste leeftijd te resolveren. De auteurs demonstreren dat het mengen van populaties met verschillende verhittingsschalen de staart van de distributie kunstmatig afvlakt, waardoor een volume-vullende geometrie ($b \approx 1$) wordt nagebootst, zelfs als de werkelijke onderliggende geometrie een dunne laag is ($b=0.5$).
• Mechanisme versus Energiebron: De auteurs verduidelijken dat hoewel de verstrooiing plaatsvindt in het midvlak (via wolken), de energie mogelijk afkomstig is van bewegingen binnen het vlak (bijv. aangedreven door spiraalstructuren) die tijdens wolk-kruisingen verticaal worden omgezet. De exponent meet de locatie van de omzetting, niet de uiteindelijke bron van de energie.

Concluderend stelt het artikel vast dat de Melkwegschijf wordt verhit nabij het vlak door een evoluerende populatie van objecten met de massa van reusachtige moleculaire wolken, en dat de vorm van de verticale actieverdeling een robuuste, modelonafhankelijke diagnostiek is voor de ruimtelijke verdeling van de verstrooiingsagenten.