The impact of numerical derivatives on radial velocity extraction

Technische Samenvatting: De impact van numerieke afgeleiden op de extractie van radiale snelheid

Probleemstelling
Metingen van radiale snelheid (RV) zijn cruciaal voor het detecteren van exoplaneten, met name voor de zogenaamde 'Earth-twins', waarbij de amplitude van het doelsignaal ongeveer 9 cm s⁻¹ bedraagt. Hoewel template matching (TM) een standaardtechniek is voor RV-extractie, vertrouwen veel moderne implementaties — specifiek line-by-line (LBL) methoden — op de numerieke berekening van de afgeleide van een stellaire template om de optimale snelheidsverschuiving te bepalen. Dit artikel onderzoekt een potentieel systematische fout die inherent is aan deze aanpak: de contaminatie van de afgeleide van de stellaire template door ruis. Hoewel de stellaire template in theoretische afleidingen vaak als ruisvrij wordt aangenomen, wordt deze in de praktijk geconstrueerd uit gestapelde observaties. De auteurs stellen dat het berekenen van een numerieke afgeleide (bijv. via eindige verschillen) op een ruizige template biases en spreiding introduceert in de RV-extractie, wat potentieel de precisie die nodig is voor het detecteren van planeten met een lage massa kan beperken.

Methodologie
De studie maakt gebruik van een vergelijkende analyse tussen twee RV-extractieparadigma's toegepast op dezelfde data:

1. Klassieke Template Matching (TM): Een $\chi^2$-minimalisatieprocedure die het geobserveerde spectrum uitlijnt met een stellaire template, waarbij rekening wordt gehouden met continuümvariaties via een eerste-graads polynoom.
2. Afgeleide-gebaseerde methode: Een formulering gebaseerd op een eerste-graads Taylor-expansie (Connes 1985; Bouchy et al. 2001), waarbij de RV-verschuiving wordt berekend met behulp van de ratio van het fluxverschil tot de afgeleide van de template ($\partial T/\partial \lambda$).

De auteurs evalueren deze methoden met behulp van twee verschillende benaderingen:

• Gesimuleerde data: Zij genereren Gaussische profielen met geïnjecteerde Poisson-ruis om spectrale lijnen te simuleren. Zij vergelijken de RV's afgeleid uit de analytische afgeleide van de Gaussische functie (de grondwaarheid) met de RV's afgeleid uit een numerieke afgeleide (eindige verschillen). De simulaties variëren in lijn-diepte, breedte, bemonstering (pixels per lijnbreedte) en signaal-ruisverhouding (SNR) van $10^2$ tot $10^6$. Ook testen zij "bossen" van Gaussische profielen (tot 1000 lijnen) om het effect van lijnmultipliciteit te beoordelen.
• Echte observaties: De auteurs passen beide methoden toe op ESPRESSO-observaties van de M3-ster K2-129. Zij construeren een stellaire template door sequentieel observaties toe te voegen (van 1 tot 78) om de SNR van de template te variëren. Zij extraheren vervolgens de RV's voor de resterende observaties met zowel de klassieke als de afgeleide methode om de residuen en onzekerheidsschattingen als functie van de template-SNR te analyseren.

Belangrijkste resultaten

• Simulatiebevindingen (Enkele lijn): Bij gesimuleerde Gaussische profielen introduceert het gebruik van een numerieke afgeleide RV-residuen op het niveau van meters per seconde wanneer de per-pixel SNR 100 is. De omvang van deze impact correleert sterk met de SNR en de lijn-diepte; ondiepere lijnen vertonen grotere spreiding. Naarmate de SNR toeneemt, neemt het verschil tussen de numerieke en analytische methoden af, en valt dit onder de ruisvloer van state-of-the-art spectrografen voor een per-pixel SNR > 1000.
• Simulatiebevindingen (Meerdere lijnen): Het opnemen van meerdere spectrale lijnen in de simulaties leidt tot een algehele afname van de contaminatie over alle SNR-regimes, wat suggereert dat de ruis gemiddeld wordt wanneer er veel lijnen worden gebruikt.
• Toepassing op echte data: Analyse van K2-129 bevestigt de trends uit de simulaties. Naarmate de template-SNR toeneemt (door meer observaties toe te voegen), convergeren de RV's die door de afgeleide methode worden afgeleid met die van de klassieke methode. Echter, zelfs met een template geconstrueerd uit 78 observaties (wat een template-SNR oplevert die 7,5 keer hoger is dan die van individuele observaties), blijft er een residu van 12,5 cm s⁻¹ bestaan tussen de twee methoden.
• Onzekerheidsschattingen: De studie merkt op dat de afgeleide methode ook de berekende RV-onzekerheden beïnvloedt. In simulaties heeft de onjuiste bepaling van de template-gradiënt invloed op de onzekerheidsschattingen bij lagere SNR's. In echte data neemt de ratio van onzekerheden tussen de twee methoden toe met de template-SNR, hoewel de auteurs opmerken dat het trekken van definitieve conclusies complex is omdat de twee methoden onzekerheden via verschillende wiskundige wegen schatten.

Betekenis en claims
Het artikel identificeert de numerieke berekening van stellaire afgeleiden in ruizige regimes als een bron van systematische fout in de RV-extractie. De auteurs beweren dat voor lage-SNR-regimes (typisch voor individuele observaties of sterren met minder beschikbare spectra), deze numerieke ruis RV-spreiding kan introduceren op de orde van grootte van meters per seconde of, in geoptimaliseerde scenario's met meerdere lijnen, residuen rond de 12,5 cm s⁻¹.

De betekenis van dit werk ligt in de demonstratie dat zelfs met hoogwaardige data, de methode van afgeleideberekening biases kan introduceren die vergelijkbaar zijn met het signaal van een 'Earth-twin' planeet. De auteurs suggereren dat hoewel het verhogen van het aantal lijnen en de template-SNR het effect verzacht, het de effectiviteit voor alle observatiecampagnes mogelijk niet volledig kan elimineren. Zij stellen voor dat toekomstige mitigatiestrategieën het gladstrijken (smoothing) van de spectra of de numerieke afgeleide kunnen omvatten om de impact van hoogfrequente fluxvariaties op de snelheidschattingen te verminderen. Het artikel concludeert dat deze systematische fout in overweging moet worden genomen bij het streven naar een RV-precisie richting de 10 cm s⁻¹ grens.