JWST/NIRSpec IFU Spectra of HR 8799 c, d, e (2.9–5.2 um; R~2,700)

JWST/NIRSpec IFU Spectra of HR 8799 c, d, eFrom work Ruffio, JB., Xuan, J.W., Chachan, Y. et al. Jupiter-like uniform metal enrichment in a system of multiple giant exoplanets. Nat Astron (2026).https://doi.org/10.1038/s41550-026-02783-z This folder contains the 2.9–5.2 micron, R ~ 2,700 spectra of the three HR 8799 planets (c, d, e) obtained with JWST/NIRSpec IFU on July 29, 2023 (UT) as part of the Cycle 1 GTO program 1188 (PI: Hodapp). See spectra_hpf1dspline_HR8799cde_Flambda.png/pdf for an illustration of the planet spectra. The total on-source integration time was ~2.2 hours. For each planet (c, d, e), four FITS files are provided with filename: HR8799[planet]combined[nrs1|nrs2]hpf1dspline_1dspectrum[Flambda|Fnu].fits where: (nrs1, nrs2) refers to the two NIRSpec detectors corresponding to 2.9-4.1 um and 4.2-5.2 um respectively. Flambda = flux density in W/m2/um Fnu = flux density in muJy (ie, 1e-6 Jy) Each FITS file contains the following extensions: WAVE = wavelength sampling (micron) FLUX = spectrum (W/m2/um or muJy) ERR = 1-sigma uncertainty (same units as flux) COV = covariance matrix CORR = correlation matrix The fits files can be loaded in Python like this: from astropy.io import fits with fits.open(spectrum_filename_Flambda) as hdulist: wv_sampling = hdulist['WAVE'].data # wavelength [micron] spectrum_Flambda = hdulist['FLUX'].data # flux [W/m2/um] speckle_std_Flambda = hdulist['ERR'].data cov_matrix_Flambda = hdulist['COV'].data corr_matrix = hdulist['CORR'].data Finally, an example Jupyter notebook is also provided: 20251105_covariance_chi2_HR8799cde_example.ipynb This notebook shows how to load the spectra and compute a chi-squared statistic using the covariance matrix.