Search

With sizable volatile envelopes but smaller radii than the solar system ice giants, sub-Neptunes have been revealed as one of the most common types of planet in the galaxy. While the spectroscopic characterization of larger sub-Neptunes (2.5-4R$_\oplus$) has revealed hydrogen-dominated atmospheres, smaller sub-Neptunes (1.6--2.5R$_\oplus$) could either host thin, rapidly evaporating hydrogen-rich atmospheres or be stable metal-rich "water worlds" with high mean molecular weight atmospheres and a fundamentally different formation and evolutionary history. Here, we present the 0.6--2.8$\mu$m JWST NIRISS/SOSS transmission spectrum of GJ 9827 d, the smallest (1.98 R$_\oplus$) warm (T$_\mathrm{eq, A_B=0.3} \sim 620$K) sub-Neptune where atmospheric absorbers have been detected to date. Our two transit observations with NIRISS/SOSS, combined with the existing HST/WFC3 spectrum, enable us to break the clouds-metallicity degeneracy. We detect water in a highly metal-enriched "steam world" atmosphere (O/H of $\sim 4$ by mass and H$_2$O found to be the background gas with a volume mixing ratio of >31%). We further show that these results are robust to stellar contamination through the transit light source effect. We do not detect escaping metastable He, which, combined with previous nondetections of escaping He and H, supports the steam atmosphere scenario. In water-rich atmospheres, hydrogen loss driven by water photolysis happens predominantly in the ionized form which eludes observational constraints. We also detect several flares in the NIRISS/SOSS light-curves with far-UV energies of the order of 10$^{30}$ erg, highlighting the active nature of the star. Further atmospheric characterization of GJ 9827 d probing carbon or sulfur species could reveal the origin of its high metal enrichment., Comment: 37 pages, 18 figures, Accepted for publication in ApJL

One of the main endeavors of the field of exoplanetary sciences is the characterization of exoplanet atmospheres on a population level. The current method of choice to accomplish this task is transmission spectroscopy, where the apparent radius of a transiting exoplanet is measured at multiple wavelengths in search of atomic and molecular absorption features produced by the upper atmosphere constituents. To extract the planetary radius from a transit light curve, it is necessary to account for the decrease in luminosity away from the center of the projected stellar disk, known as the limb darkening. Physically-motivated parametrizations of the limb darkening, in particular of the quadratic form, are commonly used in exoplanet transit light-curve fitting. Here, we show that such parametrizations can introduce significant wavelength-dependent biases in the transmission spectra currently obtained with all instrument modes of the JWST, and thus have the potential to affect atmospheric inferences. To avoid such biases, we recommend the use of standard limb-darkening parametrizations with wide uninformative priors that allow for non-physical stellar intensity profiles in the transit fits, and thus for a complete and symmetrical exploration of the parameter space. We further find that fitting the light curves at the native resolution results in errors on the measured transit depths that are significantly smaller compared to light curves that are binned in wavelength before fitting, thus potentially maximizing the amount of information that can be extracted from the data., Comment: Accepted for publication in The Astronomical Journal

Even though sub-Neptunes likely represent the most common outcome of planet formation, their natures remain poorly understood. In particular, planets near 1.5-2.5$\,R_\oplus$ often have bulk densities that can be explained equally well with widely different compositions and interior structures, resulting in grossly divergent implications for their formation. Here, we present the full 0.6-5.2$\,\mu \mathrm{m}$ JWST NIRISS/SOSS+NIRSpec/G395H transmission spectrum of the 2.2$\,R_\oplus$ TOI-270d ($4.78\,M_\oplus$, $T_\mathrm{eq}$=350-380 K), delivering unprecedented sensitivity for atmospheric characterization in the sub-Neptune regime. We detect five vibrational bands of CH$_4$ at 1.15, 1.4, 1.7, 2.3, and 3.3$\,\mu$m (9.4$\sigma$), the signature of CO$_2$ at 4.3$\,\mu$m (4.8$\sigma$), water vapor (2.5$\sigma$), and potential signatures of SO$_2$ at 4.0$\,\mu \mathrm{m}$ and CS$_2$ at 4.6$\,\mu\mathrm{m}$. Intriguingly, we find an overall highly metal-rich atmosphere, with a mean molecular weight of $5.47_{-1.14}^{+1.25}$. We infer an atmospheric metal mass fraction of $58_{-12}^{+8}\%$ and a C/O of $0.47_{-0.19}^{+0.16}$, indicating that approximately half the mass of the outer envelope is in high-molecular-weight volatiles (H$_2$O, CH$_4$, CO, CO$_2$) rather than H$_2$/He. We introduce a sub-Neptune classification scheme and identify TOI-270d as a "miscible-envelope sub-Neptune" in which H$_2$/He is well-mixed with the high-molecular-weight volatiles in a miscible supercritical metal-rich envelope. For a fully miscible envelope, we conclude that TOI-270d's interior is $90_{-4}^{+3}\,$wt$\,\%$ rock/iron, indicating that it formed as a rocky planet that accreted a few wt % of H$_2$/He, with the overall envelope metal content explained by magma-ocean/envelope reactions without the need for significant ice accretion. TOI-270d may well be an archetype of the overall population of sub-Neptunes., Comment: 25 pages, 12 figures

The hot-Neptune desert is one of the most sparsely populated regions of the exoplanet parameter space, and atmosphere observations of its few residents can provide insights into how such planets have managed to survive in such an inhospitable environment. Here, we present transmission observations of LTT 9779 b, the only known hot-Neptune to have retained a significant H/He-dominated atmosphere, taken with JWST NIRISS/SOSS. The 0.6-2.85$\mu$m transmission spectrum shows evidence for muted spectral features, rejecting a perfectly flat line at >5$\sigma$. We explore water and methane-dominated atmosphere scenarios for LTT 9779 b's terminator, and retrieval analyses reveal a continuum of potential combinations of metallicity and cloudiness. Through comparisons to previous population synthesis works and our own interior structure modelling, we are able to constrain LTT 9779 b's atmosphere metallicity to 20-850x solar. Within this range of metallicity, our retrieval analyses prefer solutions with clouds at mbar pressures, regardless of whether the atmosphere is water- or methane-dominated -- though cloud-free atmospheres with metallicities >500x solar cannot be entirely ruled out. By comparing self-consistent atmosphere temperature profiles with cloud condensation curves, we find that silicate clouds can readily condense in the terminator region of LTT 9779 b. Advection of these clouds onto the day-side could explain the high day-side albedo previously inferred for this planet and be part of a feedback loop aiding the survival of LTT 9779 b's atmosphere in the hot-Neptune desert., Comment: Accepted for publication in ApJL

Ultra-cool dwarf stars are abundant, long-lived, and uniquely suited to enable the atmospheric study of transiting terrestrial companions with JWST. Amongst them, the most prominent is the M8.5V star TRAPPIST-1 and its seven planets. While JWST Cycle 1 observations have started to yield preliminary insights into the planets, they have also revealed that their atmospheric exploration requires a better understanding of their host star. Here, we propose a roadmap to characterize the TRAPPIST-1 system -- and others like it -- in an efficient and robust manner. We notably recommend that -- although more challenging to schedule -- multi-transit windows be prioritized to mitigate the effects of stellar activity and gather up to twice more transits per JWST hour spent. We conclude that, for such systems, planets cannot be studied in isolation by small programs, but rather need large-scale, jointly space- and ground-based initiatives to fully exploit the capabilities of JWST for the exploration of terrestrial planets.

We present the first analysis of JWST near-infrared spectroscopy of stellar flares from TRAPPIST-1 during transits of rocky exoplanets. Four flares were observed from 0.6--2.8 $\mu$m with NIRISS and 0.6--3.5 $\mu$m with NIRSpec during transits of TRAPPIST-1b, f, and g. We discover P$\alpha$ and Br$\beta$ line emission and characterize flare continuum at wavelengths from 1--3.5 $\mu$m for the first time. Observed lines include H$\alpha$, P$\alpha$-P$\epsilon$, Br$\beta$, He I $\lambda$0.7062$\mu$m, two Ca II infrared triplet (IRT) lines, and the He I IRT. We observe a reversed Paschen decrement from P$\alpha$-P$\gamma$ alongside changes in the light curve shapes of these lines. The continuum of all four flares is well-described by blackbody emission with an effective temperature below 5300 K, lower than temperatures typically observed at optical wavelengths. The 0.6--1 $\mu$m spectra were convolved with the TESS response, enabling us to measure the flare rate of TRAPPIST-1 in the TESS bandpass. We find flares of 10$^{30}$ erg large enough to impact transit spectra occur at a rate of 3.6$\substack{+2.1 \\ -1.3}$ flare d$^{-1}$, $\sim$10$\times$ higher than previous predictions from K2. We measure the amount of flare contamination at 2 $\mu$m for the TRAPPIST-1b and f transits to be 500$\pm$450 and 2100$\pm$400 ppm, respectively. We find up to 80% of flare contamination can be removed, with mitigation most effective from 1.0--2.4 $\mu$m. These results suggest transits affected by flares may still be useful for atmospheric characterization efforts., Comment: 29 pages, 17 figures, 3 tables, accepted to The Astrophysical Journal

Recent work on the characterization of small exoplanets has allowed us to accumulate growing evidence that the sub-Neptunes with radii greater than $\sim2.5\,R_\oplus$ often host H$_2$/He-dominated atmospheres both from measurements of their low bulk densities and direct detections of their low mean-molecular-mass atmospheres. However, the smaller sub-Neptunes in the 1.5-2.2 R$_\oplus$ size regime are much less understood, and often have bulk densities that can be explained either by the H$_2$/He-rich scenario, or by a volatile-dominated composition known as the "water world" scenario. Here, we report the detection of water vapor in the transmission spectrum of the $1.96\pm0.08$ R$_\oplus$ sub-Neptune GJ9827d obtained with the Hubble Space Telescope. We observed 11 HST/WFC3 transits of GJ9827d and find an absorption feature at 1.4$\mu$m in its transit spectrum, which is best explained (at 3.39$\sigma$) by the presence of water in GJ9827d's atmosphere. We further show that this feature cannot be caused by unnoculted star spots during the transits by combining an analysis of the K2 photometry and transit light-source effect retrievals. We reveal that the water absorption feature can be similarly well explained by a small amount of water vapor in a cloudy H$_2$/He atmosphere, or by a water vapor envelope on GJ9827d. Given that recent studies have inferred an important mass-loss rate ($>0.5\,$M$_\oplus$/Gyr) for GJ9827d making it unlikely to retain a H-dominated envelope, our findings highlight GJ9827d as a promising water world candidate that could host a volatile-dominated atmosphere. This water detection also makes GJ9827d the smallest exoplanet with an atmospheric molecular detection to date., Comment: Published in ApJL, 11 pages, 6 figures

TRAPPIST-1 is a nearby system of seven Earth-sized, temperate, rocky exoplanets transiting a Jupiter-sized M8.5V star, ideally suited for in-depth atmospheric studies. Each TRAPPIST-1 planet has been observed in transmission both from space and from the ground, confidently rejecting cloud-free, hydrogen-rich atmospheres. Secondary eclipse observations of TRAPPIST-1 b with JWST/MIRI are consistent with little to no atmosphere given the lack of heat redistribution. Here we present the first transmission spectra of TRAPPIST-1 b obtained with JWST/NIRISS over two visits. The two transmission spectra show moderate to strong evidence of contamination from unocculted stellar heterogeneities, which dominates the signal in both visits. The transmission spectrum of the first visit is consistent with unocculted starspots and the second visit exhibits signatures of unocculted faculae. Fitting the stellar contamination and planetary atmosphere either sequentially or simultaneously, we confirm the absence of cloud-free hydrogen-rich atmospheres, but cannot assess the presence of secondary atmospheres. We find that the uncertainties associated with the lack of stellar model fidelity are one order of magnitude above the observation precision of 89 ppm (combining the two visits). Without affecting the conclusion regarding the atmosphere of TRAPPIST-1 b, this highlights an important caveat for future explorations, which calls for additional observations to characterize stellar heterogeneities empirically and/or theoretical works to improve model fidelity for such cool stars. This need is all the more justified as stellar contamination can affect the search for atmospheres around the outer, cooler TRAPPIST-1 planets for which transmission spectroscopy is currently the most efficient technique., Comment: 26 pages, 11 figures, accepted for publication in The Astrophysical Journal Letters

The Near Infrared Imager and Slitless Spectrograph instrument (NIRISS) is the Canadian Space Agency (CSA) contribution to the suite of four science instruments of JWST. As one of the three NIRISS observing modes, the Single Object Slitless Spectroscopy (SOSS) mode is tailor-made to undertake time-series observations of exoplanets to perform transit spectroscopy. The SOSS permits observing point sources between 0.6 and 2.8 um at a resolving power of 650 at 1.25 um using a slit-less cross-dispersing grism while its defocussing cylindrical lens enables observing targets as bright as J=6.7 by spreading light across 23 pixels along the cross-dispersion axis. This paper officially presents the design of the SOSS mode, its operation, characterization, and its performance, from ground-based testing and flight-based Commissioning. On-sky measurements demonstrate a peak photon conversion efficiency of 55% at 1.2 um. The first time-series on the A-type star BD+60o1753 achieves a flux stability close to the photon-noise limit, so far tested to a level of 20 parts per million on 40-minute time-scales after simply subtracting a long-term trend. Uncorrected 1/f noise residuals underneath the spectral traces add an extra source of noise equivalent to doubling the readout noise. Preliminary analysis of a HAT-P-14b transit time-series indicates that it is difficult to remove all the noise in pixels with partially saturated ramps. Overall, the SOSS delivers performance at the level required to tackle key exoplanet science programs such as detecting secondary atmospheres on terrestrial planets and measuring abundances of several chemical species in gas giants., Comment: 41 pages, 28 figures, accepted for publication in PASP

The future is now - after its long-awaited launch in December 2021, JWST began science operations in July 2022 and is already revolutionizing exoplanet astronomy. The Early Release Observations (ERO) program was designed to provide the first images and spectra from JWST, covering a multitude of science cases and using multiple modes of each on-board instrument. Here, we present transmission spectroscopy observations of the hot-Saturn WASP-96b with the Single Object Slitless Spectroscopy (SOSS) mode of the Near Infrared Imager and Slitless Spectrograph, observed as part of the ERO program. As the SOSS mode presents some unique data reduction challenges, we provide an in-depth walk-through of the major steps necessary for the reduction of SOSS data: including background subtraction, correction of 1/f noise, and treatment of the trace order overlap. We furthermore offer potential routes to correct for field star contamination, which can occur due to the SOSS mode's slitless nature. By comparing our extracted transmission spectrum with grids of atmosphere models, we find an atmosphere metallicity between 1x and 5x solar, and a solar carbon-to-oxygen ratio. Moreover, our models indicate that no grey cloud deck is required to fit WASP-96b's transmission spectrum, but find evidence for a slope shortward of 0.9$\mu$m, which could either be caused by enhanced Rayleigh scattering or the red wing of a pressure-broadened Na feature. Our work demonstrates the unique capabilities of the SOSS mode for exoplanet transmission spectroscopy and presents a step-by-step reduction guide for this new and exciting instrument., Comment: MNRAS, in press. Updated to reflect published version

The Kepler and TESS missions revealed a remarkable abundance of sub-Neptune exoplanets. Despite this abundance, our understanding of the nature and compositional diversity of sub-Neptunes remains limited, to a large part because atmospheric studies via transmission spectroscopy almost exclusively aimed for low-density sub-Neptunes and even those were often affected by high-altitude clouds. The recent TESS discovery of the hot, dense TOI-824b ($2.93\,R_\oplus$ and $18.47\,M_\oplus$) opens a new window into sub-Neptune science by enabling the study of a dense sub-Neptune via secondary eclipses. Here, we present the detection of TOI-824b's hot day side via Spitzer secondary eclipse observations in the $3.6$ and $4.5\,\mathrm{\mu m}$ channels, combined with a reanalysis of its interior composition. The measured eclipse depths (142$^{+57}_{-52}$ and 245$^{+75}_{-77}$ ppm) and brightness temperatures (1463$^{+183}_{-196}$ and 1484$^{+180}_{-202}$ K) indicate a poor heat redistribution ($f>$ 0.49) and a low Bond albedo (A$_{B}<$ 0.26). We conclude that TOI-824b could be an "exposed Neptune mantle": a planet with a Neptune-like water-rich interior that never accreted a hydrogen envelope or that subsequently lost it. The hot day-side temperature is then naturally explained by a high-metallicity envelope re-emitting the bulk of the incoming radiation from the day side. TOI-824b's density is also consistent with a massive rocky core that accreted up to 1% of hydrogen, but the observed eclipse depths favor our high-metallicity GCM simulation to a solar-metallicity GCM simulation with a likelihood ratio of 7:1. The new insights into TOI-824b's nature suggest that the sub-Neptune population may be more diverse than previously thought, with some of the dense hot sub-Neptunes potentially not hosting a hydrogen-rich envelope as generally assumed for sub-Neptunes., Comment: Accepted for publication in The Astrophysical Journal, 14 pages, 5 figures

Transmission spectroscopy provides insight into the atmospheric properties and consequently the formation history, physics, and chemistry of transiting exoplanets. However, obtaining precise inferences of atmospheric properties from transmission spectra requires simultaneously measuring the strength and shape of multiple spectral absorption features from a wide range of chemical species. This has been challenging given the precision and wavelength coverage of previous observatories. Here, we present the transmission spectrum of the Saturn-mass exoplanet WASP-39b obtained using the SOSS mode of the NIRISS instrument on the JWST. This spectrum spans $0.6 - 2.8 \mu$m in wavelength and reveals multiple water absorption bands, the potassium resonance doublet, as well as signatures of clouds. The precision and broad wavelength coverage of NIRISS-SOSS allows us to break model degeneracies between cloud properties and the atmospheric composition of WASP-39b, favoring a heavy element enhancement ("metallicity") of $\sim 10 - 30 \times$ the solar value, a sub-solar carbon-to-oxygen (C/O) ratio, and a solar-to-super-solar potassium-to-oxygen (K/O) ratio. The observations are best explained by wavelength-dependent, non-gray clouds with inhomogeneous coverage of the planet's terminator., Comment: 48 pages, 12 figures, 2 tables. Under review at Nature

Measuring the metallicity and carbon-to-oxygen (C/O) ratio in exoplanet atmospheres is a fundamental step towards constraining the dominant chemical processes at work and, if in equilibrium, revealing planet formation histories. Transmission spectroscopy provides the necessary means by constraining the abundances of oxygen- and carbon-bearing species; however, this requires broad wavelength coverage, moderate spectral resolution, and high precision that, together, are not achievable with previous observatories. Now that JWST has commenced science operations, we are able to observe exoplanets at previously uncharted wavelengths and spectral resolutions. Here we report time-series observations of the transiting exoplanet WASP-39b using JWST's Near InfraRed Camera (NIRCam). The long-wavelength spectroscopic and short-wavelength photometric light curves span 2.0 - 4.0 $\mu$m, exhibit minimal systematics, and reveal well-defined molecular absorption features in the planet's spectrum. Specifically, we detect gaseous H$_2$O in the atmosphere and place an upper limit on the abundance of CH$_4$. The otherwise prominent CO$_2$ feature at 2.8 $\mu$m is largely masked by H$_2$O. The best-fit chemical equilibrium models favour an atmospheric metallicity of 1-100$\times$ solar (i.e., an enrichment of elements heavier than helium relative to the Sun) and a sub-stellar carbon-to-oxygen (C/O) ratio. The inferred high metallicity and low C/O ratio may indicate significant accretion of solid materials during planet formation or disequilibrium processes in the upper atmosphere., Comment: 35 pages, 13 figures, 3 tables, Nature, accepted

Measuring the abundances of carbon and oxygen in exoplanet atmospheres is considered a crucial avenue for unlocking the formation and evolution of exoplanetary systems. Access to an exoplanet's chemical inventory requires high-precision observations, often inferred from individual molecular detections with low-resolution space-based and high-resolution ground-based facilities. Here we report the medium-resolution (R$\sim$600) transmission spectrum of an exoplanet atmosphere between 3-5 $\mu$m covering multiple absorption features for the Saturn-mass exoplanet WASP-39b, obtained with JWST NIRSpec G395H. Our observations achieve 1.46x photon precision, providing an average transit depth uncertainty of 221 ppm per spectroscopic bin, and present minimal impacts from systematic effects. We detect significant absorption from CO$_2$ (28.5$\sigma$) and H$_2$O (21.5$\sigma$), and identify SO$_2$ as the source of absorption at 4.1 $\mu$m (4.8$\sigma$). Best-fit atmospheric models range between 3 and 10x solar metallicity, with sub-solar to solar C/O ratios. These results, including the detection of SO$_2$, underscore the importance of characterising the chemistry in exoplanet atmospheres, and showcase NIRSpec G395H as an excellent mode for time series observations over this critical wavelength range., Comment: 44 pages, 11 figures, 3 tables. Resubmitted after revision to Nature

Carbon dioxide (CO2) is a key chemical species that is found in a wide range of planetary atmospheres. In the context of exoplanets, CO2 is an indicator of the metal enrichment (i.e., elements heavier than helium, also called "metallicity"), and thus formation processes of the primary atmospheres of hot gas giants. It is also one of the most promising species to detect in the secondary atmospheres of terrestrial exoplanets. Previous photometric measurements of transiting planets with the Spitzer Space Telescope have given hints of the presence of CO2 but have not yielded definitive detections due to the lack of unambiguous spectroscopic identification. Here we present the detection of CO2 in the atmosphere of the gas giant exoplanet WASP-39b from transmission spectroscopy observations obtained with JWST as part of the Early Release Science Program (ERS). The data used in this study span 3.0 to 5.5 {\mu}m in wavelength and show a prominent CO2 absorption feature at 4.3 {\mu}m (26{\sigma} significance). The overall spectrum is well matched by one-dimensional, 10x solar metallicity models that assume radiative-convective-thermochemical equilibrium and have moderate cloud opacity. These models predict that the atmosphere should have water, carbon monoxide, and hydrogen sulfide in addition to CO2, but little methane. Furthermore, we also tentatively detect a small absorption feature near 4.0 {\mu}m that is not reproduced by these models., Comment: 27 pages, 6 figures, Accepted for publication in Nature, data and models available at https://doi.10.5281/zenodo.6959427

After a successful launch, the James Webb Space Telescope is preparing to undertake one of its principal missions, the characterization of the atmospheres of exoplanets. The Single Object Slitless Spectroscopy (SOSS) mode of the Near Infrared Imager and Slitless Spectrograph (NIRISS) is the only observing mode that has been specifically designed for this objective. It features a wide simultaneous spectral range (0.6--2.8\,\micron) through two spectral diffraction orders. However, due to mechanical constraints, these two orders overlap slightly over a short range, potentially introducing a ``contamination'' signal in the extracted spectrum. We show that for a typical box extraction, this contaminating signal amounts to 1\% or less over the 1.6--2.8\,\micron\ range (order 1), and up to 1\% over the 0.85--0.95\,\micron\ range (order 2). For observations of exoplanet atmospheres (transits, eclipses or phase curves) where only temporal variations in flux matter, the contamination signal typically biases the results by order of 1\% of the planetary atmosphere spectral features strength. To address this problem, we developed the Algorithm to Treat Order ContAmination (ATOCA). By constructing a linear model of each pixel on the detector, treating the underlying incident spectrum as a free variable, ATOCA is able to perform a simultaneous extraction of both orders. We show that, given appropriate estimates of the spatial trace profiles, the throughputs, the wavelength solutions, as well as the spectral resolution kernels for each order, it is possible to obtain an extracted spectrum accurate to within 10\,ppm over the full spectral range., Comment: Submitted to PASP. 22 pages, 12 figures

$\texttt{Eureka!}$ is a data reduction and analysis pipeline for exoplanet time-series observations, with a particular focus on JWST data. Over the next 1-2 decades, JWST will pursue four main science themes: Early Universe, Galaxies Over Time, Star Lifecycle, and Other Worlds. Our focus is on providing the astronomy community with an open source tool for the reduction and analysis of time-series observations of exoplanets in pursuit of the fourth of these themes, Other Worlds. The goal of $\texttt{Eureka!}$ is to provide an end-to-end pipeline that starts with uncalibrated FITS files and ultimately yields precise exoplanet spectra. The pipeline has a modular structure with six stages, and each stage uses a "Eureka! Control File" (ECF) to allow for easy control of the pipeline's behavior. Stage 5 also uses a "Eureka! Parameter File" (EPF) to control the fitted parameters. We provide template ECFs for the MIRI, NIRCam, NIRISS, and NIRSpec instruments on JWST and the WFC3 instrument on the Hubble Space Telescope (HST). These templates give users a good starting point for their analyses, but $\texttt{Eureka!}$ is not intended to be used as a black box tool, and users should expect to fine-tune some settings for each observation in order to achieve optimal results. At each stage, the pipeline creates intermediate figures and outputs that allow users to compare $\texttt{Eureka!}$'s performance using different parameter settings or to compare $\texttt{Eureka!}$ with an independent pipeline. The ECF used to run each stage is also copied into the output folder from each stage to enhance reproducibility. Finally, while $\texttt{Eureka!}$ has been optimized for exoplanet observations (especially the latter stages of the code), much of the core functionality could also be repurposed for JWST time-series observations in other research domains thanks to $\texttt{Eureka!}$'s modularity., Comment: 6 pages, 1 figure, published in JOSS