Ciência da Liberação Antecipada do exoplaneta WASP

Notícias

LarLar / Notícias / Ciência da Liberação Antecipada do exoplaneta WASP

May 12, 2023

Ciência da Liberação Antecipada do exoplaneta WASP

Natureza volume 614, páginas

Nature volume 614, páginas 670–675 (2023) Cite este artigo

4559 Acessos

8 Citações

20 Altmétrico

Detalhes das métricas

O exoplaneta de massa de Saturno WASP-39b tem sido objeto de extensos esforços para determinar suas propriedades atmosféricas usando espectroscopia de transmissão1,2,3,4. No entanto, esses esforços foram prejudicados por degenerações de modelagem entre a composição e as propriedades da nuvem causadas pela qualidade limitada dos dados5,6,7,8,9. Aqui apresentamos o espectro de transmissão de WASP-39b obtido usando o modo de espectroscopia sem fendas de objeto único (SOSS) do instrumento Near Infrared Imager and Slitless Spectrograph (NIRISS) no JWST. Este espectro abrange 0,6–2,8 μm de comprimento de onda e mostra várias bandas de absorção de água, o dupleto de ressonância de potássio e assinaturas de nuvens. A precisão e a ampla cobertura de comprimento de onda do NIRISS/SOSS nos permite quebrar as degenerescências do modelo entre as propriedades da nuvem e a composição atmosférica do WASP-39b, favorecendo um aprimoramento de elemento pesado ('metalicidade') de cerca de 10 a 30 vezes o valor solar, uma razão subsolar de carbono para oxigênio (C/O) e uma razão solar para supersolar de potássio para oxigênio (K/O). As observações também são melhor explicadas por nuvens não cinzentas dependentes do comprimento de onda com cobertura não homogênea do terminador do planeta.

Observamos um trânsito de WASP-39 b usando o NIRISS10 no JWST como parte do Transiting Exoplanet Community Early Release Science Program11,12. Nossas observações duraram 8,2 h começando em 26 de julho de 2022 20:45 UTC, cobrindo o trânsito de 2,8 h, bem como 3,0 h antes e 2,4 h após o trânsito para estabelecer uma linha de base do fluxo. Os dados foram obtidos no modo SOSS, que cobre simultaneamente a faixa de comprimento de onda de 0,6 a 2,8 μm em duas ordens espectrais no mesmo detector. A ordem 1 contém a faixa espectral entre 0,6 e 2,8 μm em uma potência de resolução média de R ≣ λ/Δλ = 700, enquanto a ordem 2 fornece a faixa espectral de 0,6–1,4 μm em uma potência de resolução média de R = 1.400. No modo SOSS, os espectros são espalhados por mais de 20 pixels na direção de dispersão cruzada por meio de uma lente de desfocagem cilíndrica (consulte Dados estendidos Fig. 1), permitindo tempos de integração mais longos e reduzindo o impacto das diferenças de nível de pixel na resposta do detector. No entanto, esse desfoque resulta na sobreposição física de ambas as ordens no detector. A observação da série temporal foi composta por 537 integrações de 49,4 s (nove grupos por integração), correspondendo a um ciclo de trabalho de 89%.

Extraímos os espectros estelares das observações de séries temporais usando seis pipelines diferentes para testar o impacto das diferenças no rastreamento de ordem espectral, correção de ruído 1/f, remoção de fundo e metodologia de extração de espectro (consulte Métodos e Dados Estendidos Figs. 2 e 3) . Criamos curvas de luz espectrofotométrica para cada pipeline (Fig. 1) e somamos os dados para criar curvas de luz branca por ordem espectral (Dados Estendidos Fig. 4). As curvas espectrofotométricas e de luz branca são amplamente livres de sistemática instrumental, exceto por uma tendência linear de taxa constante no tempo e um efeito de rampa exponencial nos primeiros 15 minutos da série temporal. As profundidades de trânsito ajustadas foram agrupadas em 80 mudanças de comprimento de onda espectral na ordem 1 e 20 na ordem 2 para criar espectros de transmissão em R ≈ 300. Apresentamos os espectros dos pipelines de redução nirHiss, Supreme-SPOON e espectroscopia de trânsito na Fig. 2. Encontramos resultados consistentes entre as tubulações, com os espectros derivados também concordando com observações anteriores do Telescópio Espacial Hubble (HST) (ver também Dados Estendidos Fig. 5).

Um modelo de trânsito de exoplanetas (linha sólida) foi ajustado a cada curva de luz com ajuste cromático usando uma lei quadrática de escurecimento de membros. Os coeficientes de escurecimento dos membros, razão de raio planeta-estrela (Rp/R*) e fluxo fora do trânsito foram variados em cada canal de comprimento de onda, enquanto todos os outros parâmetros foram fixos. Os resíduos para os modelos de melhor ajuste são mostrados para cada curva de luz. A faixa de comprimento de onda para cada canal é indicada no painel a, enquanto a dispersão de partes por milhão (ppm) nos resíduos é indicada no painel b. Calculamos o ppm como o desvio padrão dos resíduos fora de trânsito. Citamos a razão do ruído do fóton previsto para cada caixa entre parênteses. As reduções são das rotinas nirHiss e chromatic_fitting descritas em Métodos. Definimos nossos erros como as incertezas 1σ extraídas dos espectros estelares. (https://github.com/afeinstein20/wasp39b_niriss_paper/blob/main/scripts/figure1.py).

30σ detection of the molecule (see Methods). Similarly, the potassium doublet at 0.768 μm is detected in the data at 6.8σ. Signatures of CO and/or CO2 are identified because of their contribution to the spectrum past 2.3 μm. We find a 3.6σ significance model preference for CO and no significant preference for CO2 (see Methods)./p> 2 μm data; the high-metallicity models (M/H > 2) expect larger transit depths than that seen in the data. The same reference model is plotted as a thick black line in both panels. We define our errors as the 1σ uncertainties extracted from the 16th and 84th percentiles of the transit depths fit from each pipeline. (https://github.com/afeinstein20/wasp39b_niriss_paper/blob/main/scripts/figure4.py)./p>2 μm; see Figs. 3 and 4b). The preferred range of metallicities provides the best fit to the shape and size of the muted water-vapour features shortward of 2 μm in combination with the larger water and CO/CO2 feature longward of 2 μm, regardless of the assumed cloud treatment in our models./p>1.5σ, which includes the zeroth-order sources. For photutils.Background2D, we used a filter size of (3, 2) pixels and a box size of (2, 2) pixels. Once the background is removed from the median F277W frame, we apply a Gaussian filter with a width of 2 to smooth out any further small-scale background noise. To apply the median F277W frame to the stage 2 science integrations, we scaled it to two isolated zeroth-order sources in the science integrations at x1 ∈ [900, 1,100], y1 ∈ [150, 250] and x2 ∈ [1,800, 2,000], y2 ∈ [150, 250]. We applied the average scaling to all integrations. We found the average F277W background scaling to be 2.81. We apply the scaled background frame to each time-series observation integration./p> 4 using this method. We interpolate over any further bad pixels by taking the median value of the two surrounding pixels along the column. We extract the spectra using a box-extraction routine and ignore any contaminants from overlapping orders or from any potential background orders. We use a box diameter of 24 pixels for both orders 1 and 2./p>5σ outliers in the resulting spectra are then identified and clipped. At present, supreme-SPOON does not explicitly treat contamination from zeroth orders of background stars that intersect the trace./p>0.001 in the jwst_niriss_spectrace_0023.fits, provided by the JWST CRDS) and a column-based 1 × median absolute deviation filter to find the illuminated pixels. Then, we calculated the column-based median of the image—using only the unilluminated pixels—and subtracted it from the image. Finally, we calculated the column-based median of the IMFD (Image-MedianFrame Difference)—using only the unilluminated pixels—and subtracted it from the image. This process is not efficient in subtracting 100% of the background contamination, which was removed during the last analysis step (spectrum decontamination)./p>2 μm, as the spectra at these wavelengths are more sensitive to the treatment of cloud properties (see next subsection for details). The best-fit spectra from PICASO, ATMO and PHOENIX indicate atmospheric metallicities of M/H = 1.7, 1.0 and 2.0, respectively. These models also consistently indicate that the C/O ratio is between 0.229 and 0.389, corresponding to the lowest C/O ratio grid point in each grid (see the main text for why models prefer lower C/O ratios). Thus, the super-solar metallicity and sub-solar C/O ratio of WASP-39b are consistent across the different model interpretations of the NIRISS/SOSS transmission spectrum./p>30σ, K at 6.8σ and CO at 3.6σ, but no notable detections of Na, CH4, CO2, HCN and H2S. The best-fit metallicity across all models is about 10–30 times solar, the best fit K/O ratio 1–2 times solar and C/O ratio 0.2. Taking the average and standard deviation of the best-fit results for all 20 runs (that is, five models on four data resolutions), we find an average M/H = 19 times solar with a standard deviation of 5 times solar and an average K/O ratio 1.5 times solar with a standard deviation of 0.26 times solar./p>

 0.85 μm, which is where the order overlap region begins. Across all pipelines, the shape of the spectra, as well as overall absorption features, cosmic-ray-removal techniques and noise levels are consistent. (https://github.com/afeinstein20/wasp39b_niriss_paper/blob/main/scripts/edfigure3.py)./p>