Radial velocity and atmospheric parameter calculations for the GaiaNIR spectrograph

Technische samenvatting: Berekeningen van radiale snelheid en atmosferische parameters voor de GaiaNIR-spectrograaf

Probleemstelling
De Gaia-missie heeft een revolutie teweeggebracht in de astrometrie en spectroscopie in het optische regime, maar wordt ernstig beperkt door interstellaire stofextinctie, wat een volledige mapping van het binnenste deel van de Melkweg, de bulge en de verre zijde van de Galactische schijf verhindert. De voorgestelde GaiaNIR-missie beoogt dit te overwinnen door te opereren in het nabij-infrarood (NIR, 800–2300 nm), waar de extinctie aanzienlijk lager is. Hoewel GaiaNIR van plan is astrometrie en fotometrie uit te voeren, zou de wetenschappelijke opbrengst exponentieel kunnen worden versterkt door een spectrograaf met een matig hoge resolutie ($R > 10.000$) toe te voegen om radiale snelheden ($v_{rad}$), atmosferische parameters ($T_{eff}$, $\log g$, $[M/H]$) en chemische abundanties van honderden miljoenen sterren te meten. Echter, technische beperkingen (detectorgrootte, telemetrie en koeleisen) maken het onmogelijk voor een enkele spectrograaf om het volledige bereik van 800–2300 nm op hoge resolutie te bestrijken. Bijgevolg is de cruciale uitdaging het identificeren van het optimale golflengtevenster dat een balans vindt tussen hoge precisie in radiale snelheid en het vermogen om nauwkeurige atmosferische parameters en abundanties af te leiden, terwijl het voordeel van verminderde interstellaire extinctie wordt gemaximaliseerd.

Methodologie
De auteurs hebben een systematische simulatiestudie uitgevoerd om optimale spectrale vensters te identificeren zonder gebruik te maken van een uitgebreide empirische bibliotheek van hoog-resolutie NIR-spectra, die momenteel niet bestaat.

• Generatie van synthetische data: De studie maakte gebruik van de BOSZ synthetische spectrale bibliotheek, waarbij 10.000 mock-spectra werden gegenereerd over een breed scala aan stellaire parameters ($T_{eff}$: 2800–10.000 K, $\log g$: 0–5,5 dex, $[M/H]$: -1,5 tot 0,5 dex). De modellen incorporeerden de LTE-benaderingen met MARCS- en ATLAS9-atmosferen, bijgewerkte continuïteits-opaciteiten en 23 moleculaire lijsten.
• Simulatieparameters: Spectra werden gegenereerd bij vier resoluties ($R = 5.000, 10.000, 15.000, 20.000$) en onderworpen aan willekeurige Gaussische ruis om signaal-ruisverhoudingen (S/N) tussen 5 en 100 te simuleren.
• Analyse van radiale snelheid: Mock-spectra werden Doppler-verschoven en kruisgecorreleerd met ideale templates met behulp van de crosscorrRV-functie. De analyse scande 100 nm brede vensters over het bereik van 850–2250 nm in stappen van 5 nm. De prestaties werden geëvalueerd met behulp van de spreiding ($\sigma$), de median absolute deviation (MAD) en het 95e percentiel (P95) van de snelheidresiduen.
• Afleiding van atmosferische parameters: De FERRE-code werd toegepast om $T_{eff}$, $\log g$ en $[M/H]$ te fitten aan de mock-spectra via $\chi^2$-minimalisatie tegen vooraf berekende grids.
• Stratificatie van de steekproef: De analyse werd gestratificeerd in vijf submonsters op basis van temperatuur en oppervlaktezwaartekracht (Rode Dwergen, Rode Reuzen, FGK dwergen/reuzen en Hele A-type sterren) om rekening te houden met de variërende dominantie van atomaire versus moleculaire lijnen.
• Validatie: De methodologie werd gevalideerd door de afgeleide precisies te vergelijken met de empirische prestaties van Gaia Data Release 3 (DR3) in het optische Ca II triplet-gebied, wat bevestigt dat de simulatiebenadering realistische onzekerheidsschattingen oplevert voor doelwitten met een hoge S/N.

Belangrijkste resultaten
De studie identificeerde zes kandidaat-vensters voor verdere verfijning, waarbij de focus uiteindelijk werd vernauwd tot twee primaire regio's:

1. De K-band (1926–2000 nm):

   • Prestaties: Deze regio biedt de hoogste precisie in radiale snelheid voor koele sterren (FGKM-typen), met geschatte precisies van 160–260 m s$^{-1}$ bij $R \approx 16.000–20.000$ voor de helderste doelwitten.
   • Atmosferische parameters: Het biedt een $T_{eff}$-precisie van beter dan 100 K en een metallischeitprecisie rond de 0,06–0,08 dex, vergelijkbaar met de huidige mogelijkheden van Gaia. De $\log g$-precisie is echter iets lager dan bij optische resultaten.
   • Extinctie: De K-band biedt een cruciaal voordeel met een interstellaire roodmaking die ongeveer 7,4 keer lager is dan in het optische spectrum, wat het in kaart brengen van stof-obscure regio's mogelijk maakt.
   • Abundanties: Tien soorten kunnen worden afgeleid (O, Na, Mg, Al, Si, Ca, Ti, V, Cr, Mn, Ni), hoewel M-dwergen voornamelijk beperkt zijn tot Ca vanwege H$_2$O-blanketing.
   • Beperkingen: De prestaties verslechteren aanzienlijk voor sterren warmer dan 7000–8000 K door een gebrek aan metallische lijnen, hoewel de Br-$\delta$ waterstoflijn beperkt nut biedt.

2. De J/H-band transitie (1158–1202 nm):

   • Prestaties: Deze regio biedt een iets lagere precisie in radiale snelheid ($\approx 340$ m s$^{-1}$) vergeleken met de K-band, maar biedt een breder effectief temperatuurbereik, reikend tot 12.000–16.000 K.
   • Atmosferische parameters: De $T_{eff}$- en metallischeitprecisies zijn competitief, hoewel de bepaling van $\log g$ een uitdaging blijft.
   • Extinctie: De reductie in extinctie is marginaal vergeleken met het optische spectrum, wat de strategische waarde voor het verkennen van de meest obscure galactische regio's beperkt.
   • Abundanties: Tien soorten zijn beschikbaar (Mg, Si, K, Ca, Ti, V, Cr, Mn, Co, Ni). Een duidelijk voordeel is de mogelijkheid om abundanties te meten voor M-dwergen onder de 4000 K (Mg, Si, K, Ca, Ti, Cr), waarbij H$_2$O-lijnen minder dominant zijn dan in de K-band.

Conclusies en Betekenis
Het artikel concludeert dat het 1926–1968 nm venster (gecentreerd rond $\lambda_C \approx 1947$ nm) de optimale strategische compromis vormt voor de GaiaNIR-spectrograaf. Hoewel het smallere 1984–2000 nm venster de hoogste precisie in radiale snelheid biedt, offert het de nauwkeurigheid van atmosferische parameters op. De 1926–1968 nm regio, werkend bij resoluties tot $R \approx 20.000$, levert radiale snelheidsprecisies van 160–260 m s$^{-1}$ en atmosferische parameterprecisies die vergelijkbaar zijn met Gaia, terwijl het profiteert van de superieure transmissie van de K-band door stof.

De auteurs benadrukken dat deze selectie GaiaNIR in staat stelt haar primaire mandaat te vervullen: het in kaart brengen van de stof-obscure binnenkant van de Melkweg en de verre zijde van de schijf. De paper merkt echter ook belangrijke kanttekeningen op:

• Crowding: De K-band kan last hebben van ernstige spectrale crowding in het Galactische vlak, wat mogelijk een lagere resolutie of kleinere venstergroottes vereist, hetgeen verdere studie behoeft.
• M-dwergen: Het onvermogen om robuuste abundanties af te leiden voor de koudste M-dwergen (door H$_2$O) is een opgemerkt nadeel, aangezien deze sterren een groot deel van de doelpopulatie vormen.
• Toekomstig werk: De auteurs stellen dat het creëren van een nieuwe, hoogwaardige lijst voor de K-band essentieel is, aangezien huidige observaties vanaf de aarde zwaar gecontamineerd zijn door tellurische lijnen. Toekomstige studies moeten de afwegingen tussen resolutie, venstergrootte en stellaire crowding adresseren om het ontwerp van de spectrograaf definitief te bepalen.

De studie claimt niet alle ontwerpoplossingen te hebben gevonden, maar biedt een op data gebaseerd kader voor het selecteren van de golflengte-regio die het wetenschappelijk rendement van de GaiaNIR-missie maximaliseert binnen haar technische beperkingen.