The impact of interpolation in high-resolution spectroscopy -- The overlooked role of interpolation in radial velocity extraction

Resumo Técnico: O Impacto da Interpolação na Espectroscopia de Alta Resolução

Definição do Problema
A detecção e caracterização de exoplanetas do tipo Terra requerem precisão de velocidade radial (RV) no nível de cm s$^{-1}$. Embora os espectrógrafos de última geração (ex: ESPRESSO, EXPRES) alcancem alta estabilidade instrumental, sinais espúrios frequentemente surgem dos pipelines de análise de dados, em vez dos próprios instrumentos. Especificamente, os métodos de extração de RV baseados em templates exigem a interpolação de espectros estelares para uma grade de comprimento de onda comum para construir templates de alto SNR e para comparar observações individuais com eles. Este artigo investiga os vieses sistemáticos introduzidos pela escolha do algoritmo de interpolação durante este processo. Os autores postulam que a interpolação assume implicitamente uma forma específica para as linhas espectrais; se essas suposições não forem válidas, o processo pode introduzir resíduos de fluxo, assimetrias de linha e vieses sistemáticos de RV, particularmente quando a variação da Velocidade Radial Terrestre Bariocêntrica (BERV) é pequena em relação ao tamanho do pixel.

Metodologia
O estudo emprega uma abordagem multifacetada combinando geração de dados sintéticos e a análise de observações reais do ESPRESSO:

1. Geração de Dados Sintéticos:

   • Perfis Gaussianos: Os autores geraram conjuntos de dados sintéticos usando funções Gaussianas para isolar os efeitos da interpolação. Eles simularam uma linha Gaussiana única e uma "floresta" de linhas (somando parâmetros derivados de linhas detectadas pelo ARES em $\epsilon$ Indi e HD40307).
   • Variação de Amostragem: Para mimetizar os efeitos do movimento de BERV e derivas instrumentais, os espectros sintéticos foram deslocados Doppler através de uma faixa de 0 a 500 m s$^{-1}$ (cobrindo o tamanho médio do pixel do ESPRESSO). O fluxo foi então interpolado de volta para a grade de comprimento de onda original.
   • Injeção de Ruído: Simulações foram realizadas em condições sem ruído e com ruído de Poisson injetado para avaliar as dependências da Razão Sinal-Ruído (SNR) (variando de SNR=10 a SNR=78).

2. Algoritmos de Interpolação:
   O estudo comparou cinco estratégias de interpolação implementadas no pipeline s-BART (Semi-Bayesian Approach for RVs with Template matching):

   • Spline cúbica (padrão).
   • Spline quadrática.
   • Polinômio de Interpolação de Hermite Cúbica por Partes (PCHIP).
   • Akima.
   • Floater-Hormann.

3. Observações Reais:
   O pipeline s-BART foi aplicado a observações do ESPRESSO de quatro estrelas com distintas estratégias de observação:

   • Alta cadência, noite única: HD40307 (K2) e $\epsilon$ Indi (K6), utilizadas para campanhas asterosísmicas.
   • Multi-noite/longo prazo: Proxima Centauri (M5) e HD39091 (G0), utilizadas para detecção e caracterização de exoplanetas.
   • As RVs foram extraídas usando diferentes algoritmos de interpolação e comparadas contra os resultados padrão de spline cúbica.

Principais Resultados

• Resíduos de Fluxo e Assimetria: A interpolação introduz resíduos no fluxo das linhas espectrais, particularmente quando o centro da linha cai entre os nós dos pixels. A magnitude desses resíduos depende da densidade de amostragem (pixels por FWHM) e do algoritmo utilizado. As estratégias de spline cúbica e Floater-Hormann geralmente produziram resíduos abaixo de 100 ppm para linhas Gaussianas únicas. No entanto, assimetrias no perfil da linha foram observadas, as quais são mais pronunciadas em linhas mal amostradas (baixo SNR ou espaçamento de pixels largo).
• Vieses de RV Simulados (Linha Única & Floresta de Linhas):
  • Em cenários Gaussianos perfeitamente simétricos e sem ruído, a interpolação introduziu um sinal de RV sistemático com uma amplitude pico-a-pico de $\sim$mm s$^{-1}$. Este sinal exibiu uma forma senoidal correlacionada com a variação de BERV.
  • Em cenários de baixo SNR (SNR=10), a amplitude deste viés sistemático aumentou significativamente, atingindo valores pico-a-pico de $\sim$20–40 m s$^{-1}$.
  • A amplitude do viés é inversamente correlacionada com o SNR. À medida que o SNR aumenta, o sinal sistemático diminui, convergindo para o regime de mm s$^{-1}$.
  • O estudo identificou três regimes na resposta de RV ao SNR: comportamento senoidal em baixo e alto SNR, e um perfil quase plano em SNR intermediário, provavelmente devido à interferência destrutiva entre efeitos sistemáticos determinísticos e estocásticos.
• Observações Reais (Noite Única):
  • Para observações de alta cadência com pequena variação de BERV (comparável ao tamanho do pixel), a escolha do algoritmo de interpolação impacta significativamente a série temporal de RV.
  • No caso de menor SNR (HD40307), diferentes algoritmos produziram resíduos pico-a-pico de $\sim$15–25 m s$^{-1}$. O algoritmo Floater-Hormann, que teve bom desempenho em simulações simétricas, na verdade aumentou o sinal sistemático em dados reais.
  • No caso de maior SNR ($\epsilon$ Indi), as diferenças entre algoritmos foram menores, com resíduos pico-a-pico em torno de 1 m s$^{-1}$, concordando dentro de 1-$\sigma$ de incerteza.
• Observações Reais (Multi-noite):
  • Para observações distribuídas sobre uma grande janela de BERV (Proxima Centauri e HD39091), o sinal sistemático causado pela interpolação é significativamente amortecido.
  • O espalhamento de RV pico-a-pico devido à escolha da interpolação foi encontrado como um limite superior de 20 cm s$^{-1}$, o que está bem abaixo da incerteza de 1-$\sigma$ das medições.

Significância e Conclusões
O artigo conclui que a interpolação espectral é uma fonte não negligenciável de erro sistemático na extração de RV, particularmente para casos científicos específicos:

1. Campanhas de Alta Cadência, Baixo SNR: Casos científicos envolvendo observações intensivas de base curta (ex: asterosismologia ou monitoramento de alta frequência) com pequenas variações de BERV são os mais suscetíveis a vieses induzidos pela interpolação, que podem atingir o nível de m s$^{-1}$.
2. Impacto Amplo: Embora o viés seja minimizado em levantamentos de exoplanetas de longo prazo e alto SNR (onde a cobertura de BERV é ampla), os autores argumentam que artefatos de interpolação afetam todos os casos científicos onde os espectros são colocados em uma grade comum. Isso inclui não apenas a detecção de exoplanetas, mas também estudos atmosféricos, análise asterossísmica e determinações de redshift cosmológico.
3. Dependência de Algoritmo: Nenhum algoritmo de interpolação é universalmente ótimo. O desempenho depende fortemente do SNR, da densidade de amostragem e da forma intrínseca das linhas espectrais. O estudo destaca que as suposições sobre a forma da linha inerentes aos métodos de interpolação podem levar à subestimação do fluxo e de vieses sistemáticos nos centros das linhas.

Os autores enfatizam que, embora o viés de nível mm s$^{-1}$ em condições ideais esteja abaixo da precisão instrumental atual, os vieses de nível m s$^{-1}$ observados em cenários de alta cadência e baixo SNR representam um sistema crítico que deve ser contabilizado em futuras análises de RV de alta precisão.