SKYSURF-11: A New Zodiacal Light Model Optimized for Optical Wavelengths

1. Declaração do Problema

A luz zodiacal (LZ), causada pela dispersão da luz solar por poeira interplanetária (IPD), é o contaminante de primeiro plano dominante em observações espaciais do ultravioleta ao infravermelho, representando mais de 90% do sinal total do céu. A subtração precisa da LZ é crítica para o estudo da Luz Galáctica Difusa (DGL) e da Luz de Fundo Extragaláctica (EBL).

• Limitações dos Modelos Existentes: Os modelos padrão utilizados atualmente (Kelsall et al. 1998 e Wright 1998) foram otimizados usando dados do instrumento COBE/DIRBE, que operou principalmente no infravermelho (1,25–240 µm).
• A Lacuna Óptica: Esses modelos falham em extrapolar com precisão para comprimentos de onda ópticos ($\lambda \lesssim 1,0$ µm). Eles dependem de parâmetros de função de fase não lineares ajustados apenas a bandas de IR, levando a superestimações significativas do brilho do céu no regime óptico.
• Consequências: A subtração imprecisa da LZ leva a erros sistemáticos nas medições da EBL, potencialmente mimetizando um sinal excessivo de "luz difusa" que poderia ser mal interpretado como nova física (por exemplo, matéria escura em decaimento) ou populações de galáxias ausentes.

2. Metodologia

Os autores apresentam o ZodiSURF, um modelo revisado de luz zodiacal 3D otimizado para a faixa óptica e de infravermelho próximo (0,3–1,6 µm).

Fonte de Dados

• Projeto SKYSURF: O modelo é calibrado usando mais de 5.000 medições de brilho de superfície do céu (sky-SB) de alta precisão a partir de dados arquivais do Telescópio Espacial Hubble (HST).
• Controle de Qualidade: O conjunto de dados inclui sinalização rigorosa para excluir imagens com alta contaminação por estrelas/galáxias, luz parasita ou proximidade da Terra/Lua. Apenas "apontamentos independentes" (separados por >10 arcmin e 2 dias) são utilizados para evitar viés.
• Geometria: Os dados abrangem ângulos solares ($\epsilon$) de $80^\circ$ a $180^\circ$, fornecendo uma ampla gama de ângulos de dispersão cruciais para restringir a função de fase, particularmente no regime de retrodispersão (gegenschein).

Estrutura de Modelagem

Os autores mantiveram a estrutura geométrica do modelo de Kelsall (nuvem suave, faixas de poeira, anel circunsolar), mas substituíram a física de dispersão por novas formas analíticas:

1. Irradiância Solar: Atualizada usando o Espectro de Referência Solar Híbrido (HSRS) do TSIS-1.
2. Função de Fase de Dispersão: Substituiu a formulação não linear de Kelsall por uma soma ponderada de três funções Henyey-Greenstein (HG):
   • $i=1$: Dispersão para frente.
   • $i=2$: Dispersão para trás.
   • $i=3$: Componente gegenschein.
   • Os parâmetros ($g_i, w_i$) são permitidos a variar linearmente com o comprimento de onda.
3. Albedo: Modelado como um parâmetro dependente do comprimento de onda.
4. Subtração de DGL: Um novo estimador de DGL foi implementado, correlacionando a DGL óptica com mapas de poeira do Planck de 350 µm e 550 µm (seguindo Postman et al. 2024), em vez da correlação menos precisa de 100 µm.

Procedimento de Ajuste

Um pipeline de ajuste de Cadeia de Markov Monte Carlo (MCMC) em três etapas foi utilizado:

1. Ajuste Conjunto: Ajuste simultâneo dos parâmetros de albedo e função de fase para cada filtro do HST independentemente.
2. Ajuste Refinado: Fixar o albedo a uma tendência linear de comprimento de onda derivada no Passo 1 e reajustar os parâmetros da função de fase.
3. Dependência do Comprimento de Onda: Ajustar tendências lineares aos parâmetros da função de fase em todos os filtros para criar um modelo contínuo válido para qualquer comprimento de onda entre 0,3–1,6 µm.

3. Contribuições Principais

• Modelo ZodiSURF: O primeiro modelo de luz zodiacal 3D especificamente otimizado para comprimentos de onda ópticos (0,3–1,6 µm), disponível publicamente como pacotes Python e IDL.
• Função de Fase Analítica: Introdução de uma função de fase Henyey-Greenstein de 3 componentes dependente do comprimento de onda que captura com precisão a dispersão para frente, para trás e gegenschein no regime óptico.
• Melhoria na Estimativa de DGL: Implementação de um método de subtração de DGL baseado em mapas de poeira submilimétricos (350/550 µm), reduzindo incertezas no regime óptico.
• Modelo de Brinquedo para Poeira Isotrópica: Uma proposta de um componente de poeira esférico, isotrópico e tênue para explicar a luz difusa residual.

4. Resultados Principais

• Precisão do Modelo: O ZodiSURF reduz a incerteza nas previsões de luz zodiacal para $\sim 4,5\%$ em comprimentos de onda ópticos.
• Comparação com Dados:
  • Diferentemente de modelos anteriores (Kelsall, Wright, ZodiPy) que superestimam o brilho do céu por margens significativas em $\lambda < 1,0$ µm, o ZodiSURF corresponde estreitamente às observações do HST.
  • Os resíduos (Céu HST – DGL – ZodiSURF) são planos em função do ângulo solar e da latitude eclíptica, indicando que a física de dispersão e a geometria 3D estão corretamente modeladas.
• Albedo e Função de Fase:
  • Albedo: Encontrado para aumentar linearmente com o comprimento de onda: $A(\lambda) \approx 0,113\lambda + 0,082$.
  • Dispersão: A dispersão para frente torna-se mais proeminente em comprimentos de onda mais longos, enquanto a dispersão para trás diminui. O componente gegenschein é relativamente constante com o comprimento de onda.
  • Implicações de Composição: O albedo geométrico derivado ($p_V \sim 0,07$) e as propriedades de dispersão sugerem que a IPD é dominada por grãos de poeira cometária de baixo albedo, provavelmente maiores que 2 µm.
• Resíduo de Luz Difusa:
  • Após subtrair o ZodiSURF e a DGL, um sinal residual de luz difusa de $0,013 \pm 0,006$ MJy sr$^{-1}$ permanece.
  • Este sinal é isotrópico e não pode ser explicado pela Luz de Galáxias Integrada (IGL) ou modelos padrão de EBL.
  • Modelo de Brinquedo: Os autores propõem que este resíduo poderia ser uma casca de poeira esférica muito tênue (componente isotrópico) localizada a $\sim 4\text{--}5$ UA com um perfil de densidade $n \propto R^{-1,78}$. Adicionar este componente ao modelo melhora significativamente o ajuste aos resíduos.

5. Significado

• Cosmologia: Ao fornecer um modelo de primeiro plano mais preciso, o ZodiSURF permite medições mais precisas da EBL, resolvendo discrepâncias entre medições diretas e contagens de luz de galáxias integradas.
• Ciência de Exoplanetas: A modelagem precisa da LZ é essencial para futuras missões (por exemplo, Observatório de Mundos Habitáveis) que visam detectar exoplanetas, pois informa a modelagem da luz exozodiacal ao redor de outras estrelas.
• Planejamento de Missões: O modelo permite cálculos mais eficientes de tempo de exposição para observatórios atuais e futuros (JWST, Euclid, Roman, SPHEREx) operando no óptico/IR próximo, potencialmente reduzindo pela metade os tempos de exposição necessários ao reduzir erros de subtração de primeiro plano.
• Trabalho Futuro: O artigo prepara o terreno para o SKYSURF-IR (usando dados do JWST) para estender essas restrições a comprimentos de onda mais longos e utilizar polarimetria para restringir ainda mais o componente de poeira isotrópica.