Showing 1–50 of 179 results for author: Madhusudhan, N

The Hycean Paradigm in the Search for Life Elsewhere

Authors: Nikku Madhusudhan

Abstract: The search for habitable conditions and signs of life on exoplanets is a major frontier in modern astronomy. Detecting atmospheric signatures of Earth-like exoplanets is challenging due to their small sizes and relatively thin atmospheres. Recently, a new class of habitable sub-Neptune exoplanets, called Hycean worlds, has been theorized. Hycean worlds are planets with H2-rich atmospheres and plan… ▽ More The search for habitable conditions and signs of life on exoplanets is a major frontier in modern astronomy. Detecting atmospheric signatures of Earth-like exoplanets is challenging due to their small sizes and relatively thin atmospheres. Recently, a new class of habitable sub-Neptune exoplanets, called Hycean worlds, has been theorized. Hycean worlds are planets with H2-rich atmospheres and planet-wide oceans with thermodynamic conditions similar to those in the Earth's oceans. Their large sizes and extended atmospheres, compared to rocky planets of similar mass, make Hycean worlds significantly more accessible to atmospheric observations. These planets open a new avenue in the search for planetary habitability and life elsewhere using spectroscopic observations with the James Webb Space Telescope (JWST). We observed the transmission spectrum of a candidate Hycean world, K2-18 b, recently with JWST in its first year of operations. The spectrum reveals multiple spectral features of carbon bearing molecules in the planetary atmosphere, leading to the first detections of methane (CH4) and carbon dioxide (CO2) in a habitable-zone exoplanet. We discuss inferences of the atmospheric chemical composition and its implications for the atmospheric, interior and surface conditions on the planet, along with the possibility of a habitable ocean underneath the atmosphere. We discuss new observational and theoretical developments in this emerging frontier and their implications for exoplanetary habitability and search for life elsewhere. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: To be published in "Astrophysics: The James Webb Space Telescope: from first light to new world views", proceedings of the Pontifical Academy of Sciences workshop February 27-29 2024, Vatican City, E.F. van Dishoeck and G. Consolmagno (eds)

The PLATO Mission

Authors: Heike Rauer, Conny Aerts, Juan Cabrera, Magali Deleuil, Anders Erikson, Laurent Gizon, Mariejo Goupil, Ana Heras, Jose Lorenzo-Alvarez, Filippo Marliani, Cesar Martin-Garcia, J. Miguel Mas-Hesse, Laurence O'Rourke, Hugh Osborn, Isabella Pagano, Giampaolo Piotto, Don Pollacco, Roberto Ragazzoni, Gavin Ramsay, Stéphane Udry, Thierry Appourchaux, Willy Benz, Alexis Brandeker, Manuel Güdel, Eduardo Janot-Pacheco , et al. (801 additional authors not shown)

Abstract: PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observati… ▽ More PLATO (PLAnetary Transits and Oscillations of stars) is ESA's M3 mission designed to detect and characterise extrasolar planets and perform asteroseismic monitoring of a large number of stars. PLATO will detect small planets (down to <2 R_(Earth)) around bright stars (<11 mag), including terrestrial planets in the habitable zone of solar-like stars. With the complement of radial velocity observations from the ground, planets will be characterised for their radius, mass, and age with high accuracy (5 %, 10 %, 10 % for an Earth-Sun combination respectively). PLATO will provide us with a large-scale catalogue of well-characterised small planets up to intermediate orbital periods, relevant for a meaningful comparison to planet formation theories and to better understand planet evolution. It will make possible comparative exoplanetology to place our Solar System planets in a broader context. In parallel, PLATO will study (host) stars using asteroseismology, allowing us to determine the stellar properties with high accuracy, substantially enhancing our knowledge of stellar structure and evolution. The payload instrument consists of 26 cameras with 12cm aperture each. For at least four years, the mission will perform high-precision photometric measurements. Here we review the science objectives, present PLATO's target samples and fields, provide an overview of expected core science performance as well as a description of the instrument and the mission profile at the beginning of the serial production of the flight cameras. PLATO is scheduled for a launch date end 2026. This overview therefore provides a summary of the mission to the community in preparation of the upcoming operational phases. △ Less

Submitted 8 June, 2024; originally announced June 2024.

Exoplanet Transit Spectroscopy with JWST NIRSpec: Diagnostics and Homogeneous Case Study of WASP-39 b

Authors: Subhajit Sarkar, Nikku Madhusudhan, Savvas Constantinou, Måns Holmberg

Abstract: The JWST has ushered in a new era of exoplanet transit spectroscopy. Among the JWST instruments, the Near-Infrared Spectrograph (NIRSpec) has the most extensive set of configurations for exoplanet time series observations. The NIRSpec Prism and G395H grating represent two extremes in NIRSpec instrument modes, with the Prism spanning a wider spectral range (0.6-5.3 $μ$m at lower resolution (R… ▽ More The JWST has ushered in a new era of exoplanet transit spectroscopy. Among the JWST instruments, the Near-Infrared Spectrograph (NIRSpec) has the most extensive set of configurations for exoplanet time series observations. The NIRSpec Prism and G395H grating represent two extremes in NIRSpec instrument modes, with the Prism spanning a wider spectral range (0.6-5.3 $μ$m at lower resolution (R$\sim$100) compared to G395H (2.87-5.14 $μ$m; R$\sim$2700). In this work, we develop a new data reduction framework, JexoPipe, to conduct a homogeneous assessment of the two NIRSpec modes for exoplanet spectroscopy. We use observations of the hot Saturn WASP-39 b obtained as part of the JWST Transiting Exoplanets ERS program to assess the spectral quality and stability between the two instrument modes at different epochs. We explore the noise sources, effect of saturation, and offsets in transmission spectra between the different instrument modes and also between the two G395H NRS detectors. We find an inter-detector offset in G395H of $\sim$ 40-50 ppm, consistent with recent studies. We find evidence for correlated noise in the Prism white light curve. We find the G395H spectrum to be of higher precision compared to the Prism at the same resolution. We also compare the JexoPipe spectra with those reported from other pipelines. Our work underscores the need for robust assessment of instrument performance and identification of optimal practices for JWST data reduction and analyses. △ Less

Submitted 10 May, 2024; originally announced May 2024.

Comments: 27 pages, 23 figures

Revealing H$_2$O dissociation in WASP-76~b through combined high- and low-resolution transmission spectroscopy

Authors: Siddharth Gandhi, Rico Landman, Ignas Snellen, Luis Welbanks, Nikku Madhusudhan, Matteo Brogi

Abstract: Numerous chemical constraints have been possible for exoplanetary atmospheres thanks to high-resolution spectroscopy (HRS) from ground-based facilities as well as low-resolution spectroscopy (LRS) from space. These two techniques have complementary strengths, and hence combined HRS and LRS analyses have the potential for more accurate abundance constraints and increased sensitivity to trace specie… ▽ More Numerous chemical constraints have been possible for exoplanetary atmospheres thanks to high-resolution spectroscopy (HRS) from ground-based facilities as well as low-resolution spectroscopy (LRS) from space. These two techniques have complementary strengths, and hence combined HRS and LRS analyses have the potential for more accurate abundance constraints and increased sensitivity to trace species. In this work we retrieve the atmosphere of the ultra-hot Jupiter WASP-76~b, using high-resolution CARMENES/CAHA and low-resolution HST WFC3 and Spitzer observations of the primary eclipse. As such hot planets are expected to have a substantial fraction of H$_2$O dissociated, we conduct retrievals including both H$_2$O and OH. We explore two retrieval models, one with self-consistent treatment of H$_2$O dissociation and another where H$_2$O and OH are vertically-homogeneous. Both models constrain H$_2$O and OH, with H$_2$O primarily detected by LRS and OH through HRS, highlighting the strengths of each technique and demonstrating the need for combined retrievals to fully constrain chemical compositions. We see only a slight preference for the H$_2$O-dissociation model given that the photospheric constraints for both are very similar, indicating $\log(\mathrm{OH/H_2O}) = 0.7^{+0.3}_{-0.3}$ at 1.5~mbar, showing that the majority of the H$_2$O in the photosphere is dissociated. However, the bulk O/H and C/O ratios inferred from the models differs significantly, and highlights the challenge of constraining bulk compositions from photospheric abundances with strong vertical chemical gradients. Further observations with JWST and ground-based facilities may help shed more light on these processes. △ Less

Submitted 3 May, 2024; originally announced May 2024.

Comments: 11 pages, 6 figures, 2 tables, accepted for publication in MNRAS

Feasibility of High-Resolution Transmission Spectroscopy for Low-Velocity Exoplanets

Authors: Connor Cheverall, Nikku Madhusudhan

Abstract: In recent years, high-resolution transmission spectroscopy in the near-infrared has led to detections of prominent molecules in several giant exoplanets on close-in orbits. This approach has traditionally relied on the large Doppler shifts of the planetary spectral lines induced by the high velocities of the close-in planets, which were considered necessary for separating them from the quasi-stati… ▽ More In recent years, high-resolution transmission spectroscopy in the near-infrared has led to detections of prominent molecules in several giant exoplanets on close-in orbits. This approach has traditionally relied on the large Doppler shifts of the planetary spectral lines induced by the high velocities of the close-in planets, which were considered necessary for separating them from the quasi-static stellar and telluric lines. In this work we demonstrate the feasibility of high-resolution transmission spectroscopy for chemical detections in atmospheres of temperate low-mass exoplanets around M dwarfs with low radial velocity variation during transit. We pursue this goal using model injection and recovery tests with H- and K- band high-resolution spectroscopy of the temperate sub-Neptune TOI-732 c, observed using the IGRINS spectrograph on Gemini-S. We show that planetary signals in transit may be recovered when the change in the planet's radial velocity is very small, down to sub-pixel velocities. This is possible due to the presence of the planetary signal in only a subset of the observed spectra. A sufficient number of out-of-transit spectra can create enough contrast between the planet signal and telluric/stellar contaminants that the planet signal does not constitute a principal component of the time-series spectra and can therefore be isolated using PCA-based detrending without relying on a significant Doppler shift. We additionally explore novel metrics for finding such signals, and investigate trends in their detectability. Our work extends the scope of high-resolution transmission spectroscopy and creates a pathway towards the characterisation of habitable sub-Neptune worlds with ground-based facilities. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: Accepted for publication in AJ

High-resolution Spectroscopic Reconnaissance of a Temperate Sub-Neptune

Authors: Samuel H. C. Cabot, Nikku Madhusudhan, Savvas Constantinou, Diana Valencia, Johanna M. Vos, Thomas Masseron, Connor J. Cheverall

Abstract: The study of temperate sub-Neptunes is the new frontier in exoplanetary science. A major development in the past year has been the first detection of carbon-bearing molecules in the atmosphere of a temperate sub-Neptune, K2-18 b, a possible Hycean world, with the James Webb Space Telescope (JWST). The JWST is poised to characterise the atmospheres of several other such planets with important impli… ▽ More The study of temperate sub-Neptunes is the new frontier in exoplanetary science. A major development in the past year has been the first detection of carbon-bearing molecules in the atmosphere of a temperate sub-Neptune, K2-18 b, a possible Hycean world, with the James Webb Space Telescope (JWST). The JWST is poised to characterise the atmospheres of several other such planets with important implications for planetary processes in the temperate regime. Meanwhile, ground-based high-resolution spectroscopy has been highly successful in detecting chemical signatures of giant exoplanets, though low-mass planets have remained elusive. In the present work, we report the atmospheric reconnaissance of a temperate sub-Neptune using ground-based high-resolution transmission spectroscopy. The long orbital period and the low systemic velocity results in a low planetary radial velocity during transit, making this system a valuable testbed for high-resolution spectroscopy of temperate sub-Neptunes. We observe high-resolution time-series spectroscopy in the H- and K-bands during the planetary transit with the IGRINS instrument (R$\sim$45,000) on Gemini-South. Using observations from a single transit we find marginal evidence (2.2$σ$) for the presence of methane (CH$_4$) in the atmosphere and no evidence for ammonia (NH$_3$) despite its strong detectability for a cloud-free H$_2$-rich atmosphere. We assess our findings using injection tests with different atmospheric scenarios, and find them to be consistent with a high CH$_4$/NH$_3$ ratio and/or the presence of high-altitude clouds. Our results demonstrate the capability of Gemini-S/IGRINS for atmospheric characterization of temperate sub-Neptunes, and the complementarity between space- and ground-based facilities in this planetary regime. △ Less

Submitted 27 March, 2024; originally announced March 2024.

Comments: Accepted for publication in ApJ Letters

VIRA: An Exoplanet Atmospheric Retrieval Framework for JWST Transmission Spectroscopy

Authors: Savvas Constantinou, Nikku Madhusudhan

Abstract: JWST observations are leading to important new insights into exoplanetary atmospheres through transmission spectroscopy. In order to harness the full potential of the broad spectral range and high sensitivity of JWST, atmospheric retrievals of exoplanets require a high level of robustness and accuracy in the underlying models. We present the VIRA retrieval framework which implements a range of mod… ▽ More JWST observations are leading to important new insights into exoplanetary atmospheres through transmission spectroscopy. In order to harness the full potential of the broad spectral range and high sensitivity of JWST, atmospheric retrievals of exoplanets require a high level of robustness and accuracy in the underlying models. We present the VIRA retrieval framework which implements a range of modelling and inference capabilities motivated by early JWST observations of exoplanet transmission spectra. This includes three complementary approaches to modelling atmospheric composition, three atmospheric aerosol models, including a physically-motivated Mie scattering approach, and consideration of correlated noise. VIRA enables a cascading retrieval architecture involving a sequence of retrievals with increasing sophistication. We demonstrate VIRA using a JWST transmission spectrum of the hot Saturn WASP-39 b in the $\sim$1-5 $μ$m range. In addition to confirming prior chemical inferences, we retrieve molecular abundances for H$_2$O, CO, CO$_2$, SO$_2$ and H$_2$S, resulting in super-solar elemental abundances of log(O/H)=$-2.0\pm0.2$, log(C/H)=$-2.1\pm0.2$ and log(S/H)=$-3.6\pm0.2$, along with C/O and S/O ratios of $0.83^{+0.05}_{-0.07}$ and $0.029^{+0.012}_{-0.009}$, respectively, in the free chemistry case. The abundances correspond to $20.1^{+10.5}_{-8.1}\times$, $28.2^{+16.3}_{-12.1}\times$ and $20.8^{+10.3}_{-7.5}\times$ solar values for O/H, C/H and S/H, respectively, compared to C/H $= 8.67\pm0.35 \times$ solar for Saturn. Our results demonstrate how JWST transmission spectroscopy combined with retrieval frameworks like VIRA can measure multi-elemental abundances for giant exoplanets and enable comparative characterisation with solar system planets. △ Less

Submitted 7 March, 2024; originally announced March 2024.

Comments: Accepted for publication in MNRAS

Possible Hycean conditions in the sub-Neptune TOI-270 d

Authors: Måns Holmberg, Nikku Madhusudhan

Abstract: The JWST has ushered in a new era in atmospheric characterisations of temperate low-mass exoplanets with recent detections of carbon-bearing molecules in the candidate Hycean world K2-18 b. We investigated JWST observations of the TOI-270 system, with two sub-Neptunes simultaneously transiting the nearby M dwarf during the visit. We report our atmospheric characterisation of the outer planet TOI-2… ▽ More The JWST has ushered in a new era in atmospheric characterisations of temperate low-mass exoplanets with recent detections of carbon-bearing molecules in the candidate Hycean world K2-18 b. We investigated JWST observations of the TOI-270 system, with two sub-Neptunes simultaneously transiting the nearby M dwarf during the visit. We report our atmospheric characterisation of the outer planet TOI-270 d, a candidate Hycean world, with JWST transmission spectroscopy using the NIRSpec G395H instrument in the 2.7-5.2 $μ$m range, combined with previous observations obtained with the HST WFC3 spectrograph (1.1-1.6 $μ$m). The spectrum reveals strong signatures of CH4 and CO2 at 3.8-4.9$σ$ and 2.9-3.9$σ$ confidence, respectively, and no evidence of NH3. The abundant CH4 and CO2, at ~0.1-1% mixing ratios, and the non-detection of NH3 are similar to the findings reported for K2-18 b and consistent with predictions for a Hycean world with a planet-wide ocean under a H2-rich atmosphere. We also report evidence of CS2 at a 2.3-3.0$σ$ confidence and a potential inference of H2O at 1.6-4.4$σ$, depending on the data analysis approach, and discuss possible interpretations of these results. The spectrum does not provide strong constraints on the presence of clouds or hazes in the observable atmosphere, nor any evidence for the effects of stellar heterogeneities, which is consistent with previous studies. For the smaller inner planet TOI-270 b, we find that the spectrum is inconsistent with a featureless spectrum at ~3$σ$, showing some preference for an H2-rich atmosphere in a super-Earth. We discuss the implications of our findings and future prospects. △ Less

Submitted 5 March, 2024; originally announced March 2024.

Comments: Accepted for publication in Astronomy & Astrophysics Letters

On the Ocean Conditions of Hycean Worlds

Authors: Frances E. Rigby, Nikku Madhusudhan

Abstract: Recent studies have suggested the possibility of Hycean worlds, characterised by deep liquid water oceans beneath H$_2$-rich atmospheres. These planets significantly widen the range of planetary properties over which habitable conditions could exist. We conduct internal structure modelling of Hycean worlds to investigate the range of interior compositions, ocean depths and atmospheric mass fractio… ▽ More Recent studies have suggested the possibility of Hycean worlds, characterised by deep liquid water oceans beneath H$_2$-rich atmospheres. These planets significantly widen the range of planetary properties over which habitable conditions could exist. We conduct internal structure modelling of Hycean worlds to investigate the range of interior compositions, ocean depths and atmospheric mass fractions possible. Our investigation explicitly considers habitable oceans, where the surface conditions are limited to those that can support potential life. The ocean depths depend on the surface gravity and temperature, confirming previous studies, and span 10s to $\sim$1000 km for Hycean conditions, reaching ocean base pressures up to $\sim$6$\times$10$^4$ bar before transitioning to high-pressure ice. We explore in detail test cases of five Hycean candidates, placing constraints on their possible ocean depths and interior compositions based on their bulk properties. We report limits on their atmospheric mass fractions admissible for Hycean conditions, as well as those allowed for other possible interior compositions. For the Hycean conditions considered, across these candidates we find the admissible mass fractions of the H/He envelopes to be $\lesssim$10$^{-3}$. At the other extreme, the maximum H/He mass fractions allowed for these planets can be up to $\sim$4-8$\%$, representing purely rocky interiors with no H$_2$O layer. These results highlight the diverse conditions possible among these planets and demonstrate their potential to host habitable conditions under vastly different circumstances to the Earth. Upcoming JWST observations of candidate Hycean worlds will allow for improved constraints on the nature of their atmospheres and interiors. △ Less

Submitted 19 February, 2024; originally announced February 2024.

Comments: Accepted for publication in MNRAS

Carbon-bearing Molecules in a Possible Hycean Atmosphere

Authors: Nikku Madhusudhan, Subhajit Sarkar, Savvas Constantinou, Måns Holmberg, Anjali A. A. Piette, Julianne I. Moses

Abstract: The search for habitable environments and biomarkers in exoplanetary atmospheres is the holy grail of exoplanet science. The detection of atmospheric signatures of habitable Earth-like exoplanets is challenging due to their small planet-star size contrast and thin atmospheres with high mean molecular weight. Recently, a new class of habitable exoplanets, called Hycean worlds, has been proposed, de… ▽ More The search for habitable environments and biomarkers in exoplanetary atmospheres is the holy grail of exoplanet science. The detection of atmospheric signatures of habitable Earth-like exoplanets is challenging due to their small planet-star size contrast and thin atmospheres with high mean molecular weight. Recently, a new class of habitable exoplanets, called Hycean worlds, has been proposed, defined as temperate ocean-covered worlds with H2-rich atmospheres. Their large sizes and extended atmospheres, compared to rocky planets of the same mass, make Hycean worlds significantly more accessible to atmospheric spectroscopy with the JWST. Here we report a transmission spectrum of the candidate Hycean world, K2-18 b, observed with the JWST NIRISS and NIRSpec instruments in the 0.9-5.2 $μ$m range. The spectrum reveals strong detections of methane (CH4) and carbon dioxide (CO2) at 5$σ$ and 3$σ$ confidence, respectively, with high volume mixing ratios of ~1% each in a H2-rich atmosphere. The abundant CH4 and CO2 along with the non-detection of ammonia (NH3) are consistent with chemical predictions for an ocean under a temperate H2-rich atmosphere on K2-18 b. The spectrum also suggests potential signs of dimethyl sulfide (DMS), which has been predicted to be an observable biomarker in Hycean worlds, motivating considerations of possible biological activity on the planet. The detection of CH4 resolves the long-standing missing methane problem for temperate exoplanets and the degeneracy in the atmospheric composition of K2-18 b from previous observations. We discuss possible implications of the findings, open questions, and future observations to explore this new regime in the search for life elsewhere. △ Less

Submitted 4 October, 2023; v1 submitted 11 September, 2023; originally announced September 2023.

Comments: Accepted for publication in ApJ Letters, updated with proof edits

Chemical Conditions on Hycean Worlds

Authors: Nikku Madhusudhan, Julianne I. Moses, Frances Rigby, Edouard Barrier

Abstract: Traditionally, the search for life on exoplanets has been predominantly focused on rocky exoplanets. Hycean worlds are a class of habitable sub-Neptunes with planet-wide oceans and H2-rich atmospheres. Their broad range of possible sizes and temperatures lead to a wide habitable zone and high potential for discovery and atmospheric characterization using transit spectroscopy. Over a dozen candidat… ▽ More Traditionally, the search for life on exoplanets has been predominantly focused on rocky exoplanets. Hycean worlds are a class of habitable sub-Neptunes with planet-wide oceans and H2-rich atmospheres. Their broad range of possible sizes and temperatures lead to a wide habitable zone and high potential for discovery and atmospheric characterization using transit spectroscopy. Over a dozen candidate Hycean planets are already known to be transiting nearby M dwarfs, making them promising targets for atmospheric characterization with the James Webb Space Telescope (JWST). In this work, we investigate possible chemical conditions on a canonical Hycean world, focusing on (a) the present and primordial molecular composition of the atmosphere, and (b) the inventory of bioessential elements for the origin and sustenance of life in the ocean. Based on photochemical and kinetic modeling for a range of conditions, we discuss the possible chemical evolution and observable present-day composition of its atmosphere. In particular, for reduced primordial conditions the early atmospheric evolution passes through a phase that is rich in organic molecules that could provide important feedstock for prebiotic chemistry. We investigate avenues for delivering bioessential metals to the ocean, considering the challenging lack of weathering from a rocky surface and the ocean separated from the rocky core by a thick icy mantle. Based on ocean depths from internal structure modelling and elemental estimates for the early Earth's oceans, we estimate the requirements for bioessential metals in such a planet. We find that the requirements can be met for plausible assumptions about impact history and atmospheric sedimentation, and supplemented by other steady state sources. We discuss the observational prospects for atmospheric characterisation of Hycean worlds. △ Less

Submitted 23 June, 2023; originally announced June 2023.

Comments: To be published in Faraday Discussions 2023, Astrochemistry at High Resolution

Exoplanet Spectroscopy with JWST NIRISS: Diagnostics and Case Studies

Authors: Måns Holmberg, Nikku Madhusudhan

Abstract: The James Webb Space Telescope (JWST) is ushering in a new era in remote sensing of exoplanetary atmospheres. Atmospheric retrievals of exoplanets can be highly sensitive to high-precision JWST data. It is, therefore, imperative to characterise the instruments and noise sources using early observations to enable robust characterisation of exoplanetary atmospheres using JWST-quality spectra. The pr… ▽ More The James Webb Space Telescope (JWST) is ushering in a new era in remote sensing of exoplanetary atmospheres. Atmospheric retrievals of exoplanets can be highly sensitive to high-precision JWST data. It is, therefore, imperative to characterise the instruments and noise sources using early observations to enable robust characterisation of exoplanetary atmospheres using JWST-quality spectra. The present work is a step in that direction, focusing on the NIRISS SOSS instrument mode, with a wavelength coverage of 0.6 - 2.8 μm and R ~ 700. Using a custom-built pipeline, JExoRes, we investigate key diagnostics of NIRISS SOSS with observations of two giant exoplanets, WASP-39 b and WASP-96 b, as case studies. We conduct a detailed evaluation of the different aspects of the data reduction and analysis, including sources of contamination, 1/f noise, and system properties such as limb darkening. The slitless nature of NIRISS SOSS makes it susceptible to contamination due to background sources. We present a method to model and correct for dispersed field stars which can significantly improve the accuracy of the observed spectra. In doing so, we also report an empirically determined throughput function for the instrument. We find significant correlated noise in the derived spectra, which may be attributed to 1/f noise, and discuss its implications for spectral binning. We quantify the covariance matrix which would enable the consideration of correlated noise in atmospheric retrievals. Finally, we conduct a comparative assessment of NIRISS SOSS spectra of WASP-39 b reported using different pipelines and highlight important lessons for exoplanet spectroscopy with JWST NIRISS. △ Less

Submitted 7 June, 2023; originally announced June 2023.

Comments: Accepted for publication in MNRAS

Robustness Measures for Molecular Detections using High-Resolution Transmission Spectroscopy of Exoplanets

Authors: Connor J. Cheverall, Nikku Madhusudhan, Måns Holmberg

Abstract: Ground-based high-resolution transmission spectroscopy has emerged as a promising technique for detecting chemicals in transiting exoplanetary atmospheres. Despite chemical inferences in several exoplanets and previous robustness studies, a robust and consistent detrending method to remove telluric and stellar features from transmission spectra has yet to be agreed upon. In this work we investigat… ▽ More Ground-based high-resolution transmission spectroscopy has emerged as a promising technique for detecting chemicals in transiting exoplanetary atmospheres. Despite chemical inferences in several exoplanets and previous robustness studies, a robust and consistent detrending method to remove telluric and stellar features from transmission spectra has yet to be agreed upon. In this work we investigate the robustness of metrics used to optimise PCA-based detrending for high-resolution transmission spectra of exoplanets in the near-infrared. As a case study, we consider observations of the hot Jupiter HD 189733 b obtained using the CARMENES spectrograph on the 3.5 m CAHA telescope. We confirm that optimising the detrending parameters to maximise the S/N of a cross-correlation signal in the presence of noise has the potential to bias the detection significance at the planetary velocity of optimisation. However, we find that optimisation using the difference between a signal-injected cross-correlation function and the direct cross-correlation function (CCF) is more robust against over-optimisation of noise and spurious signals. We additionally examine the robustness of weighting the contribution of each order to the final CCF, and of S/N calculations. Using a prescribed robust methodology, we confirm H2O in the atmosphere of HD 189733 b (S/N = 6.1). We then investigate two further case studies, of exoplanets HD 209458 b and WASP-76 b, confirming OH in the atmosphere of WASP-76 b (S/N = 4.7), and demonstrating how non-robust methods may induce false positive or inflated detections. Our findings pave the way towards a robust framework for homogeneous characterisation of exoplanetary atmospheres using high-resolution transmission spectroscopy in the near-infrared. △ Less

Submitted 2 March, 2023; originally announced March 2023.

Comments: Accepted for publication in MNRAS

A precise blue-optical transmission spectrum from the ground: Evidence for haze in the atmosphere of WASP-74b

Authors: Petros Spyratos, Nikolay K. Nikolov, Savvas Constantinou, John Southworth, Nikku Madhusudhan, Elyar Sedaghati, David Ehrenreich, Luigi Mancini

Abstract: We report transmission spectroscopy of the bloated hot Jupiter WASP-74b in the wavelength range from 4000 to 6200 Å. We observe two transit events with the Very Large Telescope FOcal Reducer and Spectrograph (VLT FORS2) and present a new method to measure the exoplanet transit depth as a function of wavelength. The new method removes the need for a reference star in correcting the spectroscopic li… ▽ More We report transmission spectroscopy of the bloated hot Jupiter WASP-74b in the wavelength range from 4000 to 6200 Å. We observe two transit events with the Very Large Telescope FOcal Reducer and Spectrograph (VLT FORS2) and present a new method to measure the exoplanet transit depth as a function of wavelength. The new method removes the need for a reference star in correcting the spectroscopic light curves for the impact of atmospheric extinction. It also provides improved precision, compared to other techniques, reaching an average transit depth uncertainty of 211 ppm for a solar-type star of V=9.8 mag and over wavelength bins of 80 Å. The VLT transmission spectrum is analysed both individually and in combination with published data from Hubble Space Telescope (HST) and Spitzer. The spectrum is found to exhibit a mostly featureless slope and equilibrium chemistry retrievals with PLATON favour hazes in the upper atmosphere of the exoplanet. Free chemistry retrievals with AURA further support the presence of hazes. While additional constraints are possible depending on the choice of atmospheric model, they are not robust and may be influenced by residual systematics in the data sets. Our results demonstrate the utility of new techniques in the analysis of optical, ground-based spectroscopic data and can be highly complementary to follow-up observations in the infrared with JWST. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Comments: 18 pages, 12 figures, 4 tables, Accepted for publication in MNRAS

Early Insights for Atmospheric Retrievals of Exoplanets using JWST Transit Spectroscopy

Authors: Savvas Constantinou, Nikku Madhusudhan, Siddharth Gandhi

Abstract: We have entered the era of the James Webb Space Telescope (JWST). We use the first JWST transmission spectrum of the hot Saturn-mass exoplanet, WASP-39 b, obtained with the NIRSpec instrument in the 3-5 $μ$m range to investigate (a) what atmospheric constraints are possible with JWST-quality data in this spectral range, (b) requirements for atmospheric models used in retrievals, (c) effect of diff… ▽ More We have entered the era of the James Webb Space Telescope (JWST). We use the first JWST transmission spectrum of the hot Saturn-mass exoplanet, WASP-39 b, obtained with the NIRSpec instrument in the 3-5 $μ$m range to investigate (a) what atmospheric constraints are possible with JWST-quality data in this spectral range, (b) requirements for atmospheric models used in retrievals, (c) effect of differences between data reduction pipelines on retrieved atmospheric properties, and (d) complementarity between JWST data in the 3-5 $μ$m range and HST observations at shorter wavelengths. JWST spectra in the 3-5 $μ$m range provide a promising avenue for chemical detections while encompassing a window in cloud opacity for several prominent aerosols. We confirm recent inferences of CO$_2$, SO$_2$, H$_2$O, and CO in WASP-39 b, report tentative evidence for H$_2$S, and retrieve elemental abundances consistent with Saturn's metallicity. We retrieve molecular abundances with $\sim$0.3-0.6 dex precision with this relatively limited spectral range. When considering the 3-5 $μ$m data alone, reported differences in spectra with different reduction pipelines can affect abundance estimates by up to $\sim$1 dex and the detectability of less prominent species. Complementing with data at shorter wavelengths, e.g. with other JWST instruments or HST WFC3 ($\sim$0.8-1.7 $μ$m), can significantly improve the accuracy and precision of the abundance estimates. The high data quality enables constraints on aerosol properties, including their composition, modal size and extent, motivating their consideration in retrievals. Our results highlight the promise of JWST exoplanet spectroscopy, while underscoring the importance of robust data reduction and atmospheric retrieval approaches in the JWST era. △ Less

Submitted 6 January, 2023; originally announced January 2023.

Comments: Accepted for publication in ApJ Letters

On the Application of Bayesian Leave-One-Out Cross-Validation to Exoplanet Atmospheric Analysis

Authors: Luis Welbanks, Peter McGill, Michael Line, Nikku Madhusudhan

Abstract: Over the last decade, exoplanetary transmission spectra have yielded an unprecedented understanding about the physical and chemical nature of planets outside our solar system. Physical and chemical knowledge is mainly extracted via fitting competing models to spectroscopic data, based on some goodness-of-fit metric. However, current employed metrics shed little light on how exactly a given model i… ▽ More Over the last decade, exoplanetary transmission spectra have yielded an unprecedented understanding about the physical and chemical nature of planets outside our solar system. Physical and chemical knowledge is mainly extracted via fitting competing models to spectroscopic data, based on some goodness-of-fit metric. However, current employed metrics shed little light on how exactly a given model is failing at the individual data point level and where it could be improved. As the quality of our data and complexity of our models increases, there is an urgent need to better understand which observations are driving our model interpretations. Here we present the application of Bayesian leave-one-out cross-validation to assess the performance of exoplanet atmospheric models and compute the expected log pointwise predictive density (elpd$_\text{LOO}$). elpd$_\text{LOO}$ estimates the out-of-sample predictive accuracy of an atmospheric model at data point resolution providing interpretable model criticism. We introduce and demonstrate this method on synthetic HST transmission spectra of a hot Jupiter. We apply elpd$_\text{LOO}$ to interpret current observations of HAT-P-41b and assess the reliability of recent inferences of H$^-$ in its atmosphere. We find that previous detections of H$^{-}$ are dependent solely on a single data point. This new metric for exoplanetary retrievals complements and expands our repertoire of tools to better understand the limits of our models and data. elpd$_\text{LOO}$ provides the means to interrogate models at the single data point level, a prerequisite for robustly interpreting the imminent wealth of spectroscopic information coming from JWST. △ Less

Submitted 7 December, 2022; originally announced December 2022.

Comments: Accepted for publication in The Astronomical Journal

Hubble Space Telescope transmission spectroscopy for the temperate sub-Neptune TOI-270d: a possible hydrogen-rich atmosphere containing water vapour

Authors: Thomas Mikal-Evans, Nikku Madhusudhan, Jason Dittmann, Maximilian N. Guenther, Luis Welbanks, Vincent Van Eylen, Ian J. M. Crossfield, Tansu Daylan, Laura Kreidberg

Abstract: TOI-270d is a temperate sub-Neptune discovered by the Transiting Exoplanet Survey Satellite (TESS) around a bright (J=9.1mag) M3V host star. With an approximate radius of 2RE and equilibrium temperature of 350K, TOI-270d is one of the most promising small exoplanets for atmospheric characterisation using transit spectroscopy. Here we present a primary transit observation of TOI-270d made with the… ▽ More TOI-270d is a temperate sub-Neptune discovered by the Transiting Exoplanet Survey Satellite (TESS) around a bright (J=9.1mag) M3V host star. With an approximate radius of 2RE and equilibrium temperature of 350K, TOI-270d is one of the most promising small exoplanets for atmospheric characterisation using transit spectroscopy. Here we present a primary transit observation of TOI-270d made with the Hubble Space Telescope (HST) Wide Field Camera 3 (WFC3) spectrograph across the 1.126-1.644 micron wavelength range, and a 95% credible upper limit of $8.2 \times 10^{-14}$ erg s$^{-1}$ cm$^{-2}$ A$^{-1}$ arcsec$^{-2}$ for the stellar Ly-alpha emission obtained using the Space Telescope Imaging Spectrograph (STIS). The transmission spectrum derived from the TESS and WFC3 data provides evidence for molecular absorption by a hydrogen-rich atmosphere at 4-sigma significance relative to a featureless spectrum. The strongest evidence for any individual absorber is obtained for H2O, which is favoured at 3-sigma significance. When retrieving on the WFC3 data alone and allowing for the possibility of a heterogeneous stellar brightness profile, the detection significance of H2O is reduced to 2.8-sigma. Further observations are therefore required to robustly determine the atmospheric composition of TOI-270d and assess the impact of stellar heterogeneity. If confirmed, our findings would make TOI-270d one of the smallest and coolest exoplanets to date with detected atmospheric spectral features. △ Less

Submitted 28 November, 2022; originally announced November 2022.

Comments: Accepted for publication in AAS journals on November 22, 2022 (received July 5, 2022; revised October 30, 2022)

A First Look at CRIRES+: Performance Assessment and Exoplanet Spectroscopy

Authors: Måns Holmberg, Nikku Madhusudhan

Abstract: High-resolution spectroscopy has proven to be a powerful avenue for atmospheric remote sensing of exoplanets. Recently, ESO commissioned the CRIRES+ high-resolution infrared spectrograph at VLT. CRIRES+ is a cross-dispersed spectrograph with high throughput and wide wavelength coverage across the near-infrared (0.95-5.3 $μ$m), designed to be particularly suited for atmospheric characterisation of… ▽ More High-resolution spectroscopy has proven to be a powerful avenue for atmospheric remote sensing of exoplanets. Recently, ESO commissioned the CRIRES+ high-resolution infrared spectrograph at VLT. CRIRES+ is a cross-dispersed spectrograph with high throughput and wide wavelength coverage across the near-infrared (0.95-5.3 $μ$m), designed to be particularly suited for atmospheric characterisation of exoplanets. In this work, we report early insights into the performance of CRIRES+ for exoplanet spectroscopy and conduct a detailed assessment of the data reduction procedure. Because of the novelty of the instrument, we perform two independent data reduction strategies, using the official CR2RES pipeline and our new custom-built ExoRES pipeline. Using science verification observations we find that the spectral resolving power of CRIRES+ can reach $R \gtrsim 100,000$ for optimal observing conditions. Similarly, we find the signal-to-noise ratio (S/N) to be consistent with expected and empirical estimates for the observations considered. As a case study, we perform the first application of CRIRES+ to the atmospheric characterisation of an exoplanet - the ultra-hot Jupiter MASCARA-1 b. We detect CO and H$_2$O in the atmosphere of MASCARA-1 b at a S/N of 12.9 and 5.3, respectively, and a temperature inversion revealed through the CO and H$_2$O emission lines, the first for an exoplanet. We find a combined S/N of 13.8 for CO and H$_2$O together, with a preference for lower H$_2$O abundance compared to CO. Our findings demonstrate the scientific potential of CRIRES+ and highlight the excellent opportunity for high-resolution atmospheric spectroscopy of diverse exoplanets. △ Less

Submitted 25 August, 2022; v1 submitted 21 June, 2022; originally announced June 2022.

Comments: Accepted version updated with proof edits

The Origin and Evolution of Saturn: A Post-Cassini Perspective

Authors: Sushil K. Atreya, Aurélien Crida, Tristan Guillot, Cheng Li, Jonathan I. Lunine, Nikku Madhusudhan, Olivier Mousis, Michael H. Wong

Abstract: The Saturn System has been studied in detail by the Cassini-Huygens Mission. A major thrust of those investigations has been to understand how Saturn formed and evolved and to place Saturn in the context of other gas giants and planetary systems in general. Two models have been proposed for the formation of the giant planets,the core accretion model and the disk instability model. The heavy elemen… ▽ More The Saturn System has been studied in detail by the Cassini-Huygens Mission. A major thrust of those investigations has been to understand how Saturn formed and evolved and to place Saturn in the context of other gas giants and planetary systems in general. Two models have been proposed for the formation of the giant planets,the core accretion model and the disk instability model. The heavy element enrichment, core size, and internal structure of Saturn, compared to Jupiter strongly favor the core accretion model as for Jupiter. Two features of the core accretion model that are distinct from the disk instability model are the growth of a core with a mass several times that of the Earth, followed by runaway collapse of gas onto the core once a mass threshold is reached. The heavy element core grows slowly over millions of years through accretion of cm-m sized pebbles, even larger bodies, and moon sized embryos in the gaseous disk. The abundance pattern of heavy elements is thus a key constraint on formation models. C, N, S, and P at Saturn are presently known to varying degree of uncertainty. The He to H ratio in the atmosphere is crucial for understanding heat balance, interior processes, and planetary evolution, but present values at Saturn range from low to high, allowing for a wide range of possibilities. While the very low values are favored to explain excess luminosity, high values might indicate presence of layered convection in the interior, resulting in slow cooling. Additional insight into Saturn's formation comes from the unique data on the rings from Cassini's Grand Finale orbits. While the solar system is the only analog for the extra solar systems, detection of the alkali metals and water in giant exoplanets is useful for understanding the formation and evolution of Saturn, where such data are presently lacking. △ Less

Submitted 13 May, 2022; originally announced May 2022.

Comments: 32 pages, 3 tables, 9 figures

Characterising Atmospheres of Cloudy Temperate Mini-Neptunes with JWST

Authors: Savvas Constantinou, Nikku Madhusudhan

Abstract: The upcoming James Webb Space Telescope (JWST) promises a generational shift in the study of temperate mini-Neptune atmospheres using transit spectroscopy. High-altitude clouds however threaten to impede their atmospheric characterisation by muting spectral features. In this study, we systematically investigate JWST instrument configurations for characterising cloudy mini-Neptune atmospheres, asse… ▽ More The upcoming James Webb Space Telescope (JWST) promises a generational shift in the study of temperate mini-Neptune atmospheres using transit spectroscopy. High-altitude clouds however threaten to impede their atmospheric characterisation by muting spectral features. In this study, we systematically investigate JWST instrument configurations for characterising cloudy mini-Neptune atmospheres, assessing the importance of instrument choice and wavelength coverage, focusing on NIRISS and NIRSpec. We consider two temperate mini-Neptunes orbiting nearby M dwarfs, K2-18 b and TOI-732 c, with equilibrium temperatures below 400 K, as case studies and assess observations using different instrument configurations with one transit per instrument. We find that their JWST transmission spectra with modest observing time and adequate wavelength coverage can provide precise abundance constraints of key molecules H2O, CH4, and NH3 even in the presence of clouds at significantly high altitudes. The best constraints are obtained by combining all three high-resolution NIRSpec gratings (G140H+G235H+G395H) that together span the ~1-5 $μ$m range. Single-transit observations with this three-instrument configuration allow precise abundance constraints for cloud-top pressures as low as 3 mbar and 0.1 mbar for K2-18 b and TOI-732 c, respectively, assuming a nominal 10x solar metallicity. The constraints vary with instrument combinations. We find that NIRSpec G235H+G395H is the optimal two-instrument configuration, while NIRISS or NIRSpec G235H is optimal for single-instrument observations. Absent high-altitude clouds, even single-instrument observations can provide good abundance constraints for these planets. Our findings underscore the promise of JWST transmission spectroscopy for characterising temperate mini-Neptunes orbiting nearby M dwarfs. △ Less

Submitted 3 May, 2022; originally announced May 2022.

Comments: Accepted for publication in MNRAS, 19 pages, 12 figures

A survey of sodium absorption in ten giant exoplanets with high-resolution transmission spectroscopy

Authors: Adam B. Langeveld, Nikku Madhusudhan, Samuel H. C. Cabot

Abstract: The alkali metal sodium (Na) is one of the most commonly detected chemical species in the upper atmospheres of giant exoplanets. In this work we conducted a homogeneous survey of Na in a diverse sample of ten highly irradiated giant exoplanets using high-resolution transmission spectroscopy. Our sample includes nine planets with previous Na detections and one new detection. We confirm previous det… ▽ More The alkali metal sodium (Na) is one of the most commonly detected chemical species in the upper atmospheres of giant exoplanets. In this work we conducted a homogeneous survey of Na in a diverse sample of ten highly irradiated giant exoplanets using high-resolution transmission spectroscopy. Our sample includes nine planets with previous Na detections and one new detection. We confirm previous detections and assess multiple approaches for deriving Na line properties from high-resolution transmission spectra. The homogeneously measured sodium line depths were used to constrain the atmospheric heights ($H_{\text{Na}}$) with respect to the planetary radii ($R_{\text{p}}$). We assess an empirical trend describing the relative atmospheric height ($H_{\text{Na}}/R_{\text{p}}$) as a function of planetary equilibrium temperature ($T_{\text{eq}}$) and surface gravity ($g$), in which $H_{\text{Na}}/R_{\text{p}}$ decreases exponentially with $ξ\propto gT_{\text{eq}}$, approaching a constant at large $ξ$. We also report the sodium D2/D1 line ratios across our sample and find that seven targets have line ratios that are consistent with unity. Finally, we measured net blueshifted offsets of the sodium absorption lines from their rest frame wavelengths for all ten planets, corresponding to day-night wind velocities of a few km s$^{-1}$. This suggests that the broad sample of exoplanets share common underlying processes that govern atmospheric dynamics. Our study highlights a promising avenue for using high-resolution transmission spectroscopy to further our understanding of how atmospheric characteristics vary over a diverse sample of exoplanets. △ Less

Submitted 4 July, 2022; v1 submitted 3 May, 2022; originally announced May 2022.

Comments: Accepted in MNRAS. 22 pages, 10 figures

Constraints on TESS albedos for five hot Jupiters

Authors: Martin Blažek, Petr Kabáth, Anjali A. A. Piette, Nikku Madhusudhan, Marek Skarka, Ján Šubjak, David R. Anderson, Henri M. J. Boffin, Claudio C. Cáceres, Neale P. Gibson, Sergio Hoyer, Valentin D. Ivanov, Patricio M. Rojo

Abstract: Photometric observations of occultations of transiting exoplanets can place important constraints on the thermal emission and albedos of their atmospheres. We analyse photometric measurements and derive geometric albedo ($A_\mathrm{g}$) constraints for five hot Jupiters observed with TESS in the optical: WASP-18 b, WASP-36 b, WASP-43 b, WASP-50 b and WASP-51 b. For WASP-43 b, our results are compl… ▽ More Photometric observations of occultations of transiting exoplanets can place important constraints on the thermal emission and albedos of their atmospheres. We analyse photometric measurements and derive geometric albedo ($A_\mathrm{g}$) constraints for five hot Jupiters observed with TESS in the optical: WASP-18 b, WASP-36 b, WASP-43 b, WASP-50 b and WASP-51 b. For WASP-43 b, our results are complemented by a VLT/HAWK-I observation in the near-infrared at $2.09~μ$m. We derive the first geometric albedo constraints for WASP-50 b and WASP-51 b: $A_\mathrm{g}<0.445$ and $A_\mathrm{g}<0.368$, respectively. We find that WASP-43 b and WASP-18 b are both consistent with low geometric albedos ($A_\mathrm{g}<0.16$) even though they lie at opposite ends of the hot Jupiter temperature range with equilibrium temperatures of $\sim1400$ K and $\sim2500$ K, respectively. We report self-consistent atmospheric models which explain broadband observations for both planets from TESS, \HST, \Spitzer and VLT/HAWK-I. We find that the data of both hot Jupiters can be explained by thermal emission alone and inefficient day-night energy redistribution. The data do not require optical scattering from clouds/hazes, consistent with the low geometric albedos observed. △ Less

Submitted 7 April, 2022; originally announced April 2022.

Comments: 21 pages, 7 figures

Aura-3D: A Three-dimensional Atmospheric Retrieval Framework for Exoplanet Transmission Spectra

Authors: Matthew C. Nixon, Nikku Madhusudhan

Abstract: Atmospheric retrievals of exoplanet transmission spectra allow constraints on the composition and structure of the day-night terminator region. Such retrievals in the past have typically assumed one-dimensional temperature structures which were adequate to explain extant observations. However, the increasing data quality expected from exoplanet spectroscopy with JWST motivates considerations of mu… ▽ More Atmospheric retrievals of exoplanet transmission spectra allow constraints on the composition and structure of the day-night terminator region. Such retrievals in the past have typically assumed one-dimensional temperature structures which were adequate to explain extant observations. However, the increasing data quality expected from exoplanet spectroscopy with JWST motivates considerations of multidimensional atmospheric retrievals. We present AURA-3D, a three-dimensional atmospheric retrieval framework for exoplanet transmission spectra. AURA-3D includes a forward model that enables rapid computation of transmission spectra in 3D geometry for a given atmospheric structure and can, therefore, be used for atmospheric retrievals as well as for computing spectra from General Circulation Models (GCMs). In order to efficiently explore the space of possible 3D temperature structures in retrievals, we develop a parametric 3D pressure-temperature profile which can accurately represent azimuthally-averaged temperature structures of a range of hot Jupiter GCMs. We apply our retrieval framework to simulated JWST observations of hot Jupiter transmission spectra, obtaining accurate estimates of the day-night temperature variation across the terminator as well as the abundances of chemical species. We demonstrate an example of a model hot Jupiter transmission spectrum for which a traditional 1D retrieval of JWST-quality data returns biased abundance estimates, whereas a retrieval including a day-night temperature gradient can accurately retrieve the true abundances. Our forward model also has the capability to include inhomogeneous chemistry as well as variable clouds/hazes. This new retrieval framework opens the field to detailed multidimensional atmospheric characterisation using transmission spectra of exoplanets in the JWST era. △ Less

Submitted 23 August, 2022; v1 submitted 10 January, 2022; originally announced January 2022.

Comments: 27 pages, 16 figures. Accepted for publication in ApJ

Journal ref: Nixon, M. C. and Madhusudhan, N., 2022, The Astrophysical Journal, Volume 935, Issue 2, id.73, 20 pp

On Atmospheric Retrievals of Exoplanets with Inhomogeneous Terminators

Authors: Luis Welbanks, Nikku Madhusudhan

Abstract: The complexity of atmospheric retrieval models is largely data-driven and one-dimensional models have generally been considered adequate with current data quality. However, recent studies have suggested that using 1D models in retrievals can result in anomalously cool terminator temperatures and biased abundance estimates even with existing transmission spectra of hot Jupiters. Motivated by these… ▽ More The complexity of atmospheric retrieval models is largely data-driven and one-dimensional models have generally been considered adequate with current data quality. However, recent studies have suggested that using 1D models in retrievals can result in anomalously cool terminator temperatures and biased abundance estimates even with existing transmission spectra of hot Jupiters. Motivated by these claims and upcoming high-quality transmission spectra we systematically explore the limitations of 1D models using synthetic and current observations. We use 1D models of varying complexity, both analytic and numerical, to revisit claims of biases when interpreting transmission spectra of hot Jupiters with inhomogeneous terminator compositions. Overall, we find the reported biases to be resulting from specific model assumptions rather than intrinsic limitations of 1D atmospheric models in retrieving current observations of asymmetric terminators. Additionally, we revise atmospheric retrievals of the hot Jupiter WASP-43b ($T_{\rm eq}=1440$ K) and the ultra-hot Jupiter WASP-103b ($T_{\rm eq}=2484$ K ) for which previous studies inferred abnormally cool atmospheric temperatures. We retrieve temperatures consistent with expectations. We note, however, that in the limit of extreme terminator inhomogeneities and high data quality some atmospheric inferences may conceivably be biased, although to a lesser extent than previously claimed. To address such cases, we implement a 2D retrieval framework for transmission spectra which allows accurate constraints on average atmospheric properties and provides insights into the spectral ranges where the imprints of atmospheric inhomogeneities are strongest. Our study highlights the need for careful considerations of model assumptions and data quality before attributing biases in retrieved estimates to unaccounted atmospheric inhomogeneities. △ Less

Submitted 3 May, 2022; v1 submitted 16 December, 2021; originally announced December 2021.

Comments: Replaced with version accepted for publication in The Astrophysical Journal

HyDRo: Atmospheric Retrieval of Rocky Exoplanets in Thermal Emission

Authors: Anjali A. A. Piette, Nikku Madhusudhan, Avi M. Mandell

Abstract: Emission spectroscopy is a promising technique to observe atmospheres of rocky exoplanets, probing both their chemistry and thermal profiles. We present HyDRo, an atmospheric retrieval framework for thermal emission spectra of rocky exoplanets. HyDRo does not make prior assumptions about the background atmospheric composition, and can therefore be used to interpret spectra of secondary atmospheres… ▽ More Emission spectroscopy is a promising technique to observe atmospheres of rocky exoplanets, probing both their chemistry and thermal profiles. We present HyDRo, an atmospheric retrieval framework for thermal emission spectra of rocky exoplanets. HyDRo does not make prior assumptions about the background atmospheric composition, and can therefore be used to interpret spectra of secondary atmospheres with unknown compositions. We use HyDRo to assess the chemical constraints which can be placed on rocky exoplanet atmospheres using JWST. Firstly, we identify the best currently-known rocky exoplanet candidates for spectroscopic observations in thermal emission with JWST, finding >30 known rocky exoplanets whose thermal emission will be detectable by JWST/MIRI in fewer than 10 eclipses at R~10. We then consider the observations required to characterise the atmospheres of three promising rocky exoplanets across the ~400-800 K equilibrium temperature range: Trappist-1 b, GJ 1132 b, and LHS 3844 b. Considering a range of CO_2- to H_2O-rich atmospheric compositions, we find that as few as 8 eclipses of LHS 3844 b or GJ 1132 b with MIRI will be able to place important constraints on the chemical compositions of their atmospheres. This includes confident detections of CO_2 and H_2O in the case of a cloud-free CO_2-rich composition, besides ruling out a bare rock scenario. Similarly, 30 eclipses of Trappist-1 b with MIRI/LRS can allow detections of a cloud-free CO_2-rich or CO_2-H_2O atmosphere. HyDRo will allow important atmospheric constraints for rocky exoplanets using JWST observations, providing clues about their geochemical environments. △ Less

Submitted 9 December, 2021; originally announced December 2021.

Comments: Accepted for publication in MNRAS. 21 pages, 10 figures, 4 tables

Molecular tracers of planet formation in the atmospheres of hot Jupiters

Authors: Richard Hobbs, Oliver Shorttle, Nikku Madhusudhan

Abstract: The atmospheric chemical composition of a hot Jupiter can lead to insights into where in its natal protoplanetary disk it formed and its subsequent migration pathway. We use a 1-D chemical kinetics code to compute a suite of models across a range of elemental abundances to investigate the resultant abundances of key molecules in hot jupiter atmospheres. Our parameter sweep spans metallicities betw… ▽ More The atmospheric chemical composition of a hot Jupiter can lead to insights into where in its natal protoplanetary disk it formed and its subsequent migration pathway. We use a 1-D chemical kinetics code to compute a suite of models across a range of elemental abundances to investigate the resultant abundances of key molecules in hot jupiter atmospheres. Our parameter sweep spans metallicities between 0.1x and 10x solar values for the C/H, O/H and N/H ratios, and equilibrium temperatures of 1000K and 2000K. We link this parameter sweep to the formation and migration models from previous works to predict connections between the atmospheric molecular abundances and formation pathways, for the molecules \ce{H2O}, \ce{CO}, \ce{CH4}, \ce{CO2}, \ce{HCN} and \ce{NH3}. We investigate atmospheric \ce{H2O} abundances in eight hot Jupiters reported in the literature. All eight planets fall within our predicted ranges for various formation models, however six of them are degenerate between multiple models and, hence, require additional molecular detections for constraining their formation histories. The other two planets, HD 189733~b and HD 209458~b, have water abundances that fall within ranges expected from planets that formed beyond the \ce{CO2} snowline. Finally, we investigate the detections of \ce{H2O}, \ce{CO}, \ce{CH4}, \ce{CO2}, \ce{HCN} and \ce{NH3} in the atmosphere of HD 209458~b and find that, within the framework of our model, the abundances of these molecules best match with a planet that formed between the \ce{CO2} and \ce{CO} snowlines and then underwent disk-free migration to reach its current location. △ Less

Submitted 25 July, 2022; v1 submitted 9 December, 2021; originally announced December 2021.

Comments: 14 pages, 11 figures

Atmospheric Modelling and Retrieval

Authors: Jonathan J. Fortney, Joanna K. Barstow, Nikku Madhusudhan

Abstract: This brief review focuses on methods and applications of modeling exoplanetary atmospheres. We discuss various kinds of state of the art self-consistent and retrieval models in 1D and multi-D with a focus on open questions and short- and long-term goals in the field. Expertise previously developed in modeling cool stellar atmospheres and in modeling solar system planetary atmospheres has proven va… ▽ More This brief review focuses on methods and applications of modeling exoplanetary atmospheres. We discuss various kinds of state of the art self-consistent and retrieval models in 1D and multi-D with a focus on open questions and short- and long-term goals in the field. Expertise previously developed in modeling cool stellar atmospheres and in modeling solar system planetary atmospheres has proven valuable to the field, and will continue to do so. We described upcoming opportunities for making progress in our understanding of atmospheres, and close with what we see as the field's challenges. △ Less

Submitted 4 September, 2021; originally announced September 2021.

Comments: Accepted in Autumn of 2020. To appear as a book chapter in "ExoFrontiers: Big questions in exoplanetary science", Ed. N Madhusudhan (Bristol: IOP Publishing Ltd) AAS-IOP ebooks https://iopscience.iop.org/bookListInfo/aas-iop-astronomy

JexoSim 2.0: End-to-End JWST Simulator for Exoplanet Spectroscopy -- Implementation and Case Studies

Authors: Subhajit Sarkar, Nikku Madhusudhan

Abstract: The recently developed JWST Exoplanet Observation Simulator (JexoSim) simulates transit spectroscopic observations of exoplanets by JWST with each of its four instruments using a time-domain approach. Previously we reported the validation of JexoSim against Pandexo and instrument team simulators. In the present study, we report a substantially enhanced version, JexoSim 2.0, which improves on the o… ▽ More The recently developed JWST Exoplanet Observation Simulator (JexoSim) simulates transit spectroscopic observations of exoplanets by JWST with each of its four instruments using a time-domain approach. Previously we reported the validation of JexoSim against Pandexo and instrument team simulators. In the present study, we report a substantially enhanced version, JexoSim 2.0, which improves on the original version through incorporation of new noise sources, enhanced treatment of stellar and planetary signals and instrumental effects, as well as improved user-operability and optimisations for increased speed and efficiency. A near complete set of instrument modes for exoplanet time-series observations is now included. In this paper we report the implementation of JexoSim 2.0 and assess performance metrics for JWST in end-member scenarios using the hot Jupiter HD 209458 b and the mini-Neptune K2-18 b. We show how JexoSim can be used to compare performance across the different JWST instruments, selecting an optimal combination of instrument and subarray modes, producing synthetic transmission spectra for each planet. These studies indicate that the 1.4 μm water feature detected in the atmosphere of K2-18 b using the Hubble WFC3 might be observable in just one transit observation with JWST with either NIRISS or NIRSpec. JexoSim 2.0 can be used to investigate the impact of complex noise and systematic effects on the final spectrum, plan observations and test the feasibility of novel science cases for JWST. It can also be customised for other astrophysical applications beyond exoplanet spectroscopy. JexoSim 2.0 is now available for use by the scientific community. △ Less

Submitted 1 September, 2021; originally announced September 2021.

Comments: Accepted for publication in MNRAS, 21 pages, 9 figures, 4 tables

Habitability and Biosignatures of Hycean Worlds

Authors: Nikku Madhusudhan, Anjali A. A. Piette, Savvas Constantinou

Abstract: We investigate a new class of habitable planets composed of water-rich interiors with massive oceans underlying H2-rich atmospheres, referred to here as Hycean worlds. With densities between those of rocky super-Earths and more extended mini-Neptunes, Hycean planets can be optimal candidates in the search for exoplanetary habitability and may be abundant in the exoplanet population. We investigate… ▽ More We investigate a new class of habitable planets composed of water-rich interiors with massive oceans underlying H2-rich atmospheres, referred to here as Hycean worlds. With densities between those of rocky super-Earths and more extended mini-Neptunes, Hycean planets can be optimal candidates in the search for exoplanetary habitability and may be abundant in the exoplanet population. We investigate the bulk properties (masses, radii, and temperatures), potential for habitability, and observable biosignatures of Hycean planets. We show that Hycean planets can be significantly larger compared to previous considerations for habitable planets, with radii as large as 2.6 Earth radii (2.3 Earth radii) for a mass of 10 Earth masses (5 Earth masses). We construct the Hycean habitable zone (HZ), considering stellar hosts from late M to sun-like stars, and find it to be significantly wider than the terrestrial-like HZ. While the inner boundary of the Hycean HZ corresponds to equilibrium temperatures as high as ~500 K for late M dwarfs, the outer boundary is unrestricted to arbitrarily large orbital separations. Our investigations include tidally locked `Dark Hycean' worlds that permit habitable conditions only on their permanent nightsides and `Cold Hycean' worlds that see negligible irradiation. Finally, we investigate the observability of possible biosignatures in Hycean atmospheres. We find that a number of trace terrestrial biomarkers which may be expected to be present in Hycean atmospheres would be readily detectable using modest observing time with the James Webb Space Telescope (JWST). We identify a sizable sample of nearby potential Hycean planets that can be ideal targets for such observations in search of exoplanetary biosignatures. △ Less

Submitted 1 September, 2021; v1 submitted 24 August, 2021; originally announced August 2021.

Comments: Published in The Astrophysical Journal

Transmission spectroscopy with VLT FORS2: a featureless spectrum for the low-density transiting exoplanet WASP-88b

Authors: Petros Spyratos, Nikolay Nikolov, John Southworth, Savvas Constantinou, Nikku Madhusudhan, Aarynn L. Carter, Ernst J. W. de Mooij, Jonathan J. Fortney, Neale P. Gibson, Jayesh M. Goyal, Christiane Helling, Nathan J. Mayne, Thomas Mikal-Evans

Abstract: We present ground-based optical transmission spectroscopy of the low-density hot Jupiter WASP-88b covering the wavelength range 4413-8333 Å with the FORS2 spectrograph on the Very Large Telescope. The FORS2 white light curves exhibit a significant time-correlated noise which we model using a Gaussian Process and remove as a wavelength-independent component from the spectroscopic light curves. We a… ▽ More We present ground-based optical transmission spectroscopy of the low-density hot Jupiter WASP-88b covering the wavelength range 4413-8333 Å with the FORS2 spectrograph on the Very Large Telescope. The FORS2 white light curves exhibit a significant time-correlated noise which we model using a Gaussian Process and remove as a wavelength-independent component from the spectroscopic light curves. We analyse complementary photometric observations from the Transiting Exoplanet Survey Satellite and refine the system properties and ephemeris. We find a featureless transmission spectrum with increased absorption towards shorter wavelengths. We perform an atmospheric retrieval analysis with the AURA code, finding tentative evidence for haze in the upper atmospheric layers and a lower likelihood for a dense cloud deck. Whilst our retrieval analysis results point toward clouds and hazes, further evidence is needed to definitively reject a clear-sky scenario. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Comments: 19 pages, 12 figures, 3 tables, Accepted for publication in MNRAS

How Deep Is the Ocean? Exploring the phase structure of water-rich sub-Neptunes

Authors: Matthew C. Nixon, Nikku Madhusudhan

Abstract: Understanding the internal structures of planets with a large H$_2$O component is important for the characterisation of sub-Neptune planets. The finding that the mini-Neptune K2-18b could host a liquid water ocean beneath a mostly hydrogen envelope motivates a detailed examination of the phase structures of water-rich planets. To this end, we present new internal structure models for super-Earths… ▽ More Understanding the internal structures of planets with a large H$_2$O component is important for the characterisation of sub-Neptune planets. The finding that the mini-Neptune K2-18b could host a liquid water ocean beneath a mostly hydrogen envelope motivates a detailed examination of the phase structures of water-rich planets. To this end, we present new internal structure models for super-Earths and mini-Neptunes that enable detailed characterisation of a planet's water component. We use our models to explore the possible phase structures of water worlds and find that a diverse range of interiors are possible, from oceans sandwiched between two layers of ice to supercritical interiors beneath steam atmospheres. We determine how the bulk properties and surface conditions of a water world affect its ocean depth, finding that oceans can be up to hundreds of times deeper than on Earth. For example, a planet with a 300 K surface can possess H$_2$O oceans with depths from 30-500 km, depending on its mass and composition. We also constrain the region of mass-radius space in which planets with H/He envelopes could host liquid H$_2$O, noting that the liquid phase can persist at temperatures up to 647 K at high pressures of $218$-$7\times10^4$ bar. Such H/He envelopes could contribute significantly to the planet radius while retaining liquid water at the surface, depending on the planet mass and temperature profile. Our findings highlight the exciting possibility that habitable conditions may be present on planets much larger than Earth. △ Less

Submitted 3 June, 2021; originally announced June 2021.

Comments: 19 pages, 15 figures. Accepted for publication in MNRAS

Ground-Based Transmission Spectroscopy with VLT FORS2: Evidence for faculae and clouds in the optical spectrum of the warm Saturn WASP-110b

Authors: Nikolay Nikolov, Gracjan Maciejewski, Savvas Constantinou, Nikku Madhusudhan, Jonathan J. Fortney, Barry Smalley, Aarynn L. Carter, Ernst J. W. de Mooij, Benjamin Drummond, Neale P. Gibson, Christiane Helling, Nathan J. Mayne, Thomas Mikal-Evans, David K. Sing, Jamie Wilson

Abstract: We present a ground-based optical transmission spectrum for the warm Saturn-mass exoplanet WASP-110b from two transit observations made with the FOcal Reducer and Spectrograph (FORS2) on the Very Large Telescope (VLT). The spectrum covers the wavelength range from 4000 to 8333Å, which is binned in 46 transit depths measured to an averaged precision of 220 parts per million (ppm) over an averaged 8… ▽ More We present a ground-based optical transmission spectrum for the warm Saturn-mass exoplanet WASP-110b from two transit observations made with the FOcal Reducer and Spectrograph (FORS2) on the Very Large Telescope (VLT). The spectrum covers the wavelength range from 4000 to 8333Å, which is binned in 46 transit depths measured to an averaged precision of 220 parts per million (ppm) over an averaged 80Å~bin for a Vmag=12.8 star. The measured transit depths are unaffected by a dilution from a close A-type field dwarf, which was fully resolved. The overall main characteristic of the transmission spectrum is an increasing radius with wavelength and a lack of the theoretically predicted pressure-broadened sodium and potassium absorption features for a cloud-free atmosphere. We analyze archival high-resolution optical spectroscopy and find evidence for low to moderate activity of the host star, which we take into account in the atmospheric retrieval analysis. Using the AURA retrieval code, we find that the observed transmission spectrum can be best explained by a combination of unocculted stellar faculae and a cloud deck. Transmission spectra of cloud-free and hazy atmospheres are rejected at a high confidence. With a possible cloud deck at its terminator, WASP-110b joins the increasing population of irradiated hot-Jupiter exoplanets with cloudy atmospheres observed in transmission. △ Less

Submitted 13 May, 2021; originally announced May 2021.

Comments: Accepted for publication in The Astronomical Journal

Aurora: A Generalised Retrieval Framework for Exoplanetary Transmission Spectra

Authors: Luis Welbanks, Nikku Madhusudhan

Abstract: Atmospheric retrievals of exoplanetary transmission spectra provide important constraints on various properties such as chemical abundances, cloud/haze properties, and characteristic temperatures, at the day-night atmospheric terminator. To date, most spectra have been observed for giant exoplanets due to which retrievals typically assume H-rich atmospheres. However, recent observations of mini-Ne… ▽ More Atmospheric retrievals of exoplanetary transmission spectra provide important constraints on various properties such as chemical abundances, cloud/haze properties, and characteristic temperatures, at the day-night atmospheric terminator. To date, most spectra have been observed for giant exoplanets due to which retrievals typically assume H-rich atmospheres. However, recent observations of mini-Neptunes/super-Earths, and the promise of upcoming facilities including JWST, call for a new generation of retrievals that can address a wide range of atmospheric compositions and related complexities. Here we report Aurora, a next-generation atmospheric retrieval framework that builds upon state-of-the-art architectures and incorporates the following key advancements: a) a generalised compositional retrieval allowing for H-rich and H-poor atmospheres, b) a generalised prescription for inhomogeneous clouds/hazes, c) multiple Bayesian inference algorithms for high-dimensional retrievals, d) modular considerations for refraction, forward scattering, and Mie-scattering, and e) noise modeling functionalities. We demonstrate Aurora on current and/or synthetic observations of hot Jupiter HD209458b, mini-Neptune K218b, and rocky exoplanet TRAPPIST1d. Using current HD209458b spectra, we demonstrate the robustness of our framework and cloud/haze prescription against assumptions of H-rich/H-poor atmospheres, improving on previous treatments. Using real and synthetic spectra of K218b, we demonstrate the agnostic approach to confidently constrain its bulk atmospheric composition and obtain precise abundance estimates. For TRAPPIST1d, 10 JWST NIRSpec transits can enable identification of the main atmospheric component for cloud-free CO$_2$-rich and N$_2$-rich atmospheres, and abundance constraints on trace gases including initial indications of O$_3$ if present at enhanced levels ($\sim$10-100x Earth levels). △ Less

Submitted 15 March, 2021; originally announced March 2021.

Comments: Accepted for publication in ApJ

Sulfur Chemistry in the Atmospheres of Warm and Hot Jupiters

Authors: Richard Hobbs, Paul Rimmer, Oliver Shorttle, Nikku Madhusudhan

Abstract: We present and validate a new network of atmospheric thermo-chemical and photo-chemical sulfur reactions. We use a 1-D chemical kinetics model to investigate these reactions as part of a broader HCNO chemical network in a series of hot and warm Jupiters. We find that temperatures approaching $1400\,\mathrm{K}$ are favourable for the production of H2S and HS around $\mathrm{10^{-3}\,bar}$, the atmo… ▽ More We present and validate a new network of atmospheric thermo-chemical and photo-chemical sulfur reactions. We use a 1-D chemical kinetics model to investigate these reactions as part of a broader HCNO chemical network in a series of hot and warm Jupiters. We find that temperatures approaching $1400\,\mathrm{K}$ are favourable for the production of H2S and HS around $\mathrm{10^{-3}\,bar}$, the atmospheric level where detection by transit spectroscopy may be possible, leading to mixing ratios of around $10^{-6}$. At lower temperatures, down to $1000\,\mathrm{K}$, the abundance of S2 can be up to a mixing ratio of $10^{-5}$ at the same pressure, at the expense of H2S and HS, which are depleted down to a mixing ratio of $10^{-6}$. We also investigate how the inclusion of sulfur can manifest in an atmosphere indirectly, by its effect on the abundance of non-sulfur-bearing species. We find that in a model of the atmosphere of HD 209458 b, the inclusion of sulfur can lower the abundance of NH3, CH4 and HCN by up to two orders of magnitude around $\mathrm{10^{-3}\,bar}$. In the atmosphere of the warm Jupiter 51 Eri b, we additionally find the inclusion of sulphur depletes the peak abundance of CO2 by a factor of five, qualitatively consistent with prior models. We note that many of the reactions used in the network have poorly determined rates, especially at higher temperatures. To obtain a truly accurate idea of the impact of sulfur chemistry in hot and warm Jupiter atmospheres, new measurements of these reaction rates must take place. △ Less

Submitted 25 June, 2021; v1 submitted 20 January, 2021; originally announced January 2021.

Comments: 17 pages, 9 Figures

Assessing telluric correction methods for Na detections with high-resolution exoplanet transmission spectroscopy

Authors: Adam B. Langeveld, Nikku Madhusudhan, Samuel H. C. Cabot, Simon T. Hodgkin

Abstract: Using high-resolution ground-based transmission spectroscopy to probe exoplanetary atmospheres is difficult due to the inherent telluric contamination from absorption in Earth's atmosphere. A variety of methods have previously been used to remove telluric features in the optical regime and calculate the planetary transmission spectrum. In this paper we present and compare two such methods, specifi… ▽ More Using high-resolution ground-based transmission spectroscopy to probe exoplanetary atmospheres is difficult due to the inherent telluric contamination from absorption in Earth's atmosphere. A variety of methods have previously been used to remove telluric features in the optical regime and calculate the planetary transmission spectrum. In this paper we present and compare two such methods, specifically focusing on Na detections using high-resolution optical transmission spectra: (1) calculating the telluric absorption empirically based on the airmass, and (2) using a model of the Earth's transmission spectrum. We test these methods on the transmission spectrum of the hot Jupiter HD 189733 b using archival data obtained with the HARPS spectrograph during three transits. Using models for Centre-to-Limb Variation and the Rossiter-McLaughlin effect, spurious signals which are imprinted within the transmission spectrum are reduced. We find that correcting tellurics with an atmospheric model of the Earth is more robust and produces consistent results when applied to data from different nights with changing atmospheric conditions. We confirm the detection of sodium in the atmosphere of HD 189733 b, with doublet line contrasts of -0.64 $\pm$ 0.07 % (D2) and -0.53 $\pm$ 0.07 % (D1). The average line contrast corresponds to an effective photosphere in the Na line located around 1.13 R$_p$. We also confirm an overall blueshift of the line centroids corresponding to net atmospheric eastward winds with a speed of 1.8 $\pm$ 1.2 km/s. Our study highlights the importance of accurate telluric removal for consistent and reliable characterisation of exoplanetary atmospheres using high-resolution transmission spectroscopy. △ Less

Submitted 21 April, 2021; v1 submitted 13 January, 2021; originally announced January 2021.

Comments: Published in MNRAS, 13 pages, 9 figures. Updated with proof corrections

The Hubble PanCET Program: A Metal-rich Atmosphere for the Inflated Hot Jupiter HAT-P-41b

Authors: Kyle B. Sheppard, Luis Welbanks, Avi Mandell, Nikku Madhusudhan, Nikolay Nikolov, Drake Deming, Gregory W. Henry, Michael H. Williamson, David K. Sing, Mercedes López-Morales, Jegug Ih, Jorge Sanz-Forcada, Panayotis Lavvas, Gilda E. Ballester, Thomas M. Evans, Antonio García Muñoz, Leonardo A. Dos Santos

Abstract: We present a comprehensive analysis of the 0.3--5\,$μ$m transit spectrum for the inflated hot Jupiter HAT-P-41b. The planet was observed in transit with Hubble STIS and WFC3 as part of the Hubble Panchromatic Comparative Exoplanet Treasury (PanCET) program, and we combine those data with warm \textit{Spitzer} transit observations. We extract transit depths from each of the data sets, presenting th… ▽ More We present a comprehensive analysis of the 0.3--5\,$μ$m transit spectrum for the inflated hot Jupiter HAT-P-41b. The planet was observed in transit with Hubble STIS and WFC3 as part of the Hubble Panchromatic Comparative Exoplanet Treasury (PanCET) program, and we combine those data with warm \textit{Spitzer} transit observations. We extract transit depths from each of the data sets, presenting the STIS transit spectrum (0.29--0.93\,$μ$m) for the first time. We retrieve the transit spectrum both with a free-chemistry retrieval suite (AURA) and a complementary chemical equilibrium retrieval suite (PLATON) to constrain the atmospheric properties at the day-night terminator. Both methods provide an excellent fit to the observed spectrum. Both AURA and PLATON retrieve a metal-rich atmosphere for almost all model assumptions (most likely O/H ratio of $\log_{10}{Z/Z_{\odot}} = 1.46^{+0.53}_{-0.68}$ and $\log_{10}{Z/Z_{\odot}} = 2.33^{+0.23}_{-0.25}$, respectively); this is driven by a 4.9-$σ$ detection of H$_2$O as well as evidence of gas absorption in the optical ($>$2.7-$σ$ detection) due to Na, AlO and/or VO/TiO, though no individual species is strongly detected. Both retrievals determine the transit spectrum to be consistent with a clear atmosphere, with no evidence of haze or high-altitude clouds. Interior modeling constraints on the maximum atmospheric metallicity ($\log_{10}{Z/Z_{\odot}} < 1.7$) favor the AURA results. The inferred elemental oxygen abundance suggests that HAT-P-41b has one of the most metal-rich atmospheres of any hot Jupiters known to date. Overall, the inferred high metallicity and high inflation make HAT-P-41b an interesting test case for planet formation theories. △ Less

Submitted 19 October, 2020; originally announced October 2020.

Comments: Resubmitted to AAS Journals after revisions

Search for radio emission from the exoplanets Qatar-1b and WASP-80b near 150 MHz using the Giant Metrewave Radio Telescope

Authors: D. A. Green, N. Madhusudhan

Abstract: We present radio observations made towards the exoplanets Qatar-1b and WASP-80b near 150~MHz with the Giant Meterwave Radio Telescope. These targets are relatively nearby irradiated giant exoplanets, a hot Jupiter and a hot Saturn, with sizes comparable to Jupiter but different masses and lower densities. Both the targets are expected to host extended H/He envelopes like Jupiter, with comparable o… ▽ More We present radio observations made towards the exoplanets Qatar-1b and WASP-80b near 150~MHz with the Giant Meterwave Radio Telescope. These targets are relatively nearby irradiated giant exoplanets, a hot Jupiter and a hot Saturn, with sizes comparable to Jupiter but different masses and lower densities. Both the targets are expected to host extended H/He envelopes like Jupiter, with comparable or larger magnetic moments. No radio emission was detected from these exoplanets, with 3sigma limits of 5.9 and 5.2 mJy for Qatar-1b and WASP-80b, respectively, from these targeted observations. These are considerably deeper limits than those available for exoplanets from wide field surveys at similar frequencies. We also present archival VLA observations of a previously reported radio source close to Vir 61 (which has three exoplanets). The VLA observations resolve the source, which we identify as an extragalactic radio source, i.e. a chance association with Vir 61. Additionally, we cross-match a recent exoplanet catalogue with the TIFR GMRT Sky Survey ADR1 radio catalogue, but do not find any convincing associations. △ Less

Submitted 9 November, 2020; v1 submitted 13 October, 2020; originally announced October 2020.

Comments: MNRAS in press

On the Temperature Profiles and Emission Spectra of Mini-Neptune Atmospheres

Authors: Anjali A. A. Piette, Nikku Madhusudhan

Abstract: Atmospheric observations of mini-Neptunes orbiting M-dwarfs are beginning to provide constraints on their chemical and thermal properties, while also providing clues about their interiors and potential surfaces. With their relatively large scale heights and large planet-star contrasts, mini-Neptunes are currently ideal targets towards the goal of characterising temperate low-mass exoplanets. Under… ▽ More Atmospheric observations of mini-Neptunes orbiting M-dwarfs are beginning to provide constraints on their chemical and thermal properties, while also providing clues about their interiors and potential surfaces. With their relatively large scale heights and large planet-star contrasts, mini-Neptunes are currently ideal targets towards the goal of characterising temperate low-mass exoplanets. Understanding the thermal structures and spectral appearances of mini-Neptunes is important to understand various aspects of their atmospheres, including radiative/convective energy transport, boundary conditions for the interior, and their potential habitability. In the present study, we explore these aspects of mini-Neptunes using self-consistent models of their atmospheres. We begin by exploring the effects of irradiation, internal flux, metallicity, clouds and hazes on the atmospheric temperature profiles and thermal emission spectra of temperate mini-Neptunes. In particular, we investigate the impact of these properties on the radiative-convective boundary and the thermodynamic conditions in the lower atmosphere, which serves as the interface with the interior and/or a potential surface. Building on recent suggestions of habitability of the mini-Neptune K2-18 b, we find a range of physically-motivated atmospheric conditions that allow for liquid water under the H_2-rich atmospheres of such planets. We find that observations of thermal emission with JWST/MIRI spectrophotometry can place useful constraints on the habitability of temperate mini-Neptunes such as K2-18 b, and provide more detailed constraints on the chemical and thermal properties of warmer planets such as GJ 3470 b. Our results underpin the potential of temperate mini-Neptunes such as K2-18 b as promising candidates in the search for habitable exoplanets. △ Less

Submitted 23 September, 2020; originally announced September 2020.

Comments: Accepted for publication in ApJ, 25 pages, 11 figures

A planetary system with two transiting mini-Neptunes near the radius valley transition around the bright M dwarf TOI-776

Authors: R. Luque, L. M. Serrano, K. Molaverdikhani, M. C. Nixon, J. H. Livingston, E. W. Guenther, E. Pallé, N. Madhusudhan, G. Nowak, J. Korth, W. D. Cochran, T. Hirano, P. Chaturvedi, E. Goffo, S. Albrecht, O. Barragán, C Briceño, J. Cabrera, D. Charbonneau, R. Cloutier, K. A. Collins, K. I. Collins, K. D. Colón, I. J. M. Crossfield, Sz. Csizmadia , et al. (41 additional authors not shown)

Abstract: We report the discovery and characterization of two transiting planets around the bright M1 V star LP 961-53 (TOI-776, J = 8.5 mag, M = 0.54+-0.03 Msun) detected during Sector 10 observations of the Transiting Exoplanet Survey Satellite (TESS). Combining the TESS photometry with HARPS radial velocities, as well as ground-based follow-up transit observations from MEarth and LCOGT telescopes, we mea… ▽ More We report the discovery and characterization of two transiting planets around the bright M1 V star LP 961-53 (TOI-776, J = 8.5 mag, M = 0.54+-0.03 Msun) detected during Sector 10 observations of the Transiting Exoplanet Survey Satellite (TESS). Combining the TESS photometry with HARPS radial velocities, as well as ground-based follow-up transit observations from MEarth and LCOGT telescopes, we measured for the inner planet, TOI-776 b, a period of 8.25 d, a radius of 1.85+-0.13 Re, and a mass of 4.0+-0.9 Me; and for the outer planet, TOI-776 c, a period of 15.66 d, a radius of 2.02+-0.14 Re, and a mass of 5.3+-1.8 Me. The Doppler data shows one additional signal, with a period of 34 d, associated with the rotational period of the star. The analysis of fifteen years of ground-based photometric monitoring data and the inspection of different spectral line indicators confirm this assumption. The bulk densities of TOI-776 b and c allow for a wide range of possible interior and atmospheric compositions. However, both planets have retained a significant atmosphere, with slightly different envelope mass fractions. Thanks to their location near the radius gap for M dwarfs, we can start to explore the mechanism(s) responsible for the radius valley emergence around low-mass stars as compared to solar-like stars. While a larger sample of well-characterized planets in this parameter space is still needed to draw firm conclusions, we tentatively estimate that the stellar mass below which thermally-driven mass loss is no longer the main formation pathway for sculpting the radius valley is between 0.63 and 0.54 Msun. Due to the brightness of the star, the TOI-776 system is also an excellent target for the James Webb Space Telescope, providing a remarkable laboratory to break the degeneracy in planetary interior models and to test formation and evolution theories of small planets around low-mass stars. △ Less

Submitted 30 November, 2020; v1 submitted 17 September, 2020; originally announced September 2020.

Comments: 26 pages, 15 figures, 8 tables. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 645, A41 (2021)

A Featureless Infrared Transmission Spectrum for the Super-Puff Planet Kepler-79d

Authors: Yayaati Chachan, Daniel Jontof-Hutter, Heather A. Knutson, Danica Adams, Peter Gao, Björn Benneke, Zachory Berta-Thompson, Fei Dai, Drake Deming, Eric Ford, Eve J. Lee, Jessica E. Libby-Roberts, Nikku Madhusudhan, Hannah R. Wakeford, Ian Wong

Abstract: Extremely low density planets ('super-puffs') are a small but intriguing subset of the transiting planet population. With masses in the super-Earth range ($1-10$ M$_{\oplus}$) and radii akin to those of giant planets ($>4$ R$_{\oplus}$), their large envelopes may have been accreted beyond the water snow line and many appear to be susceptible to catastrophic mass loss. Both the presence of water an… ▽ More Extremely low density planets ('super-puffs') are a small but intriguing subset of the transiting planet population. With masses in the super-Earth range ($1-10$ M$_{\oplus}$) and radii akin to those of giant planets ($>4$ R$_{\oplus}$), their large envelopes may have been accreted beyond the water snow line and many appear to be susceptible to catastrophic mass loss. Both the presence of water and the importance of mass loss can be explored using transmission spectroscopy. Here, we present new HST WFC3 spectroscopy and updated Kepler transit depth measurements for the super-puff Kepler-79d. We do not detect any molecular absorption features in the $1.1-1.7$ $μ$m WFC3 bandpass and the combination of Kepler and WFC3 data are consistent with a flat line model, indicating the presence of aerosols in the atmosphere. We compare the shape of Kepler-79d's transmission spectrum to predictions from a microphysical haze model that incorporates an outward particle flux due to ongoing mass loss. We find that photochemical hazes offer an attractive explanation for the observed properties of super-puffs like Kepler-79d, as they simultaneously render the near-infrared spectrum featureless and reduce the inferred envelope mass loss rate by moving the measured radius (optical depth unity surface during transit) to lower pressures. We revisit the broader question of mass loss rates for super-puffs and find that the age estimates and mass loss rates for the majority of super-puffs can be reconciled if hazes move the photosphere from the typically assumed pressure of $\sim 10$ mbar to $\sim 10 \; μ$bar. △ Less

Submitted 19 August, 2020; v1 submitted 12 August, 2020; originally announced August 2020.

Comments: Awaiting publication in AJ. Small updates in Table 6 and Fig 11. Table 3 and 6 will be provided in MRT format upon publication

Journal ref: AJ, 160, 201 (2020)

Considerations for Atmospheric Retrieval of High-Precision Brown Dwarf Spectra

Authors: Anjali A. A. Piette, Nikku Madhusudhan

Abstract: Isolated brown dwarfs provide remarkable laboratories for understanding atmospheric physics in the low-irradiation regime, and can be observed more precisely than exoplanets. As such, they provide a glimpse into the future of high-SNR observations of exoplanets. In this work, we investigate several new considerations that are important for atmospheric retrievals of high-quality thermal emission sp… ▽ More Isolated brown dwarfs provide remarkable laboratories for understanding atmospheric physics in the low-irradiation regime, and can be observed more precisely than exoplanets. As such, they provide a glimpse into the future of high-SNR observations of exoplanets. In this work, we investigate several new considerations that are important for atmospheric retrievals of high-quality thermal emission spectra of sub-stellar objects. We pursue this using an adaptation of the HyDRA atmospheric retrieval code. We propose a parametric pressure-temperature (P-T) profile for brown dwarfs consisting of multiple atmospheric layers, parameterised by the temperature change across each layer. This model allows the steep temperature gradient of brown dwarf atmospheres to be accurately retrieved while avoiding commonly-encountered numerical artefacts. The P-T model is especially flexible in the photosphere, which can reach a few tens of bar for T-dwarfs. We demonstrate an approach to include model uncertainties in the retrieval, focusing on uncertainties introduced by finite spectral and vertical resolution in the atmospheric model used for retrieval (~8\% in the present case). We validate our retrieval framework by applying it to a simulated data set and then apply it to the HST/WFC3 spectrum of the T-dwarf 2MASS J2339+1352. We retrieve sub-solar abundances of H2O and CH4 in the object at ~0.1 dex precision. Additionally, we constrain the temperature structure to within ~100 K in the photosphere. Our results demonstrate the promise of high-SNR spectra to provide high-precision abundance estimates of sub-stellar objects. △ Less

Submitted 29 July, 2020; originally announced July 2020.

Comments: Accepted for publication in MNRAS, 21 pages, 15 figures

Ground-Based Transmission Spectroscopy with FORS2: A featureless optical transmission spectrum and detection of H$_2$O for the ultra-hot Jupiter WASP-103b

Authors: J. Wilson, N. P. Gibson, N. Nikolov, S. Constantinou, N. Madhusudhan, J. Goyal, J. K. Barstow, A. L. Carter, E. J. W. de Mooij, B. Drummond, T. Mikal-Evans, C. Helling, N. J. Mayne, D. K. Sing

Abstract: We report ground-based transmission spectroscopy of the highly irradiated and ultra-short period hot-Jupiter WASP-103b covering the wavelength range $\approx$ 400-600 nm using the FORS2 instrument on the Very Large Telescope. The light curves show significant time-correlated noise which is mainly invariant in wavelength and which we model using a Gaussian process. The precision of our transmission… ▽ More We report ground-based transmission spectroscopy of the highly irradiated and ultra-short period hot-Jupiter WASP-103b covering the wavelength range $\approx$ 400-600 nm using the FORS2 instrument on the Very Large Telescope. The light curves show significant time-correlated noise which is mainly invariant in wavelength and which we model using a Gaussian process. The precision of our transmission spectrum is improved by applying a common-mode correction derived from the white light curve, reaching typical uncertainties in transit depth of $\approx$ 2x10$^{-4}$ in wavelength bins of 15 nm. After correction for flux contamination from a blended companion star, our observations reveal a featureless spectrum across the full range of the FORS2 observations and we are unable to confirm the Na absorption previously inferred using Gemini/GMOS or the strong Rayleigh scattering observed using broad-band light curves. We performed a Bayesian atmospheric retrieval on the full optical-infrared transmission spectrum using the additional data from Gemini/GMOS, HST/WFC3 and Spitzer observations and recover evidence for H$_2$O absorption at the 4.0$σ$ level. However, our observations are not able to completely rule out the presence of Na, which is found at 2.0$σ$ in our retrievals. This may in part be explained by patchy/inhomogeneous clouds or hazes damping any absorption features in our FORS2 spectrum, but an inherently small scale height also makes this feature challenging to probe from the ground. Our results nonetheless demonstrate the continuing potential of ground-based observations for investigating exoplanet atmospheres and emphasise the need for the application of consistent and robust statistical techniques to low-resolution spectra in the presence of instrumental systematics. △ Less

Submitted 27 July, 2020; originally announced July 2020.

Comments: 17 pages, 8 Figures, Accepted for publication in MNRAS

Detection of Na in WASP-21b's lower and upper atmosphere

Authors: G. Chen, N. Casasayas-Barris, E. Palle, L. Welbanks, N. Madhusudhan, R. Luque, F. Murgas

Abstract: Optical transmission spectroscopy provides crucial constraints on the reference pressure levels and scattering properties for hot Jupiter atmospheres. For certain planets, where alkali atoms are detected in the atmosphere, their line profiles could serve as a good probe to link upper and lower atmospheric layers. WASP-21b is a Saturn-mass hot Jupiter orbiting a thick disc star, with a low density… ▽ More Optical transmission spectroscopy provides crucial constraints on the reference pressure levels and scattering properties for hot Jupiter atmospheres. For certain planets, where alkali atoms are detected in the atmosphere, their line profiles could serve as a good probe to link upper and lower atmospheric layers. WASP-21b is a Saturn-mass hot Jupiter orbiting a thick disc star, with a low density and an equilibrium temperature of 1333 K, which makes it a good target for transmission spectroscopy. Here, we present a low-resolution transmission spectrum for WASP-21b based in one transit observed by the OSIRIS spectrograph at the 10.4 m Gran Telescopio Canarias (GTC), and a high-resolution transmission spectrum based in three transits observed by HARPS-N at Telescopio Nazinale Galileo (TNG) and HARPS at the ESO 3.6 m telescope. We performed spectral retrieval analysis on GTC's low-resolution transmission spectrum and report the detection of Na at a confidence level of $>$3.5-$σ$. The Na line exhibits a broad line profile that can be attributed to pressure broadening, indicating a mostly clear planetary atmosphere. The spectrum shows a tentative excess absorption at the K D$_1$ line. Using HARPS-N and HARPS, we spectrally resolved the Na doublet transmission spectrum. An excess absorption at the Na doublet is detected during the transit, and shows a radial velocity shift consistent with the planet orbital motion. We proposed a metric to quantitatively distinguish hot Jupiters with relatively clear atmospheres from others, and WASP-21b has the largest metric value among all the characterized hot Jupiters. The detection of Na at both lower and upper atmosphere of WASP-21b reveals that it is an ideal target for future follow-up observations, providing the opportunity to understand the nature of its atmosphere across a wide range of pressure levels. △ Less

Submitted 27 July, 2020; originally announced July 2020.

Comments: 16 pages, 14 figures, accepted for publication in A&A on July 27, 2020

Journal ref: A&A 642, A54 (2020)

TESS unveils the optical phase curve of KELT-1b. Thermal emission and ellipsoidal variation from the brown dwarf companion, and activity from the star

Authors: C. von Essen, M. Mallonn, A. Piette, N. B. Cowan, N. Madhusudhan, E. Agol, V. Antoci, K. Poppenhaeger, K. G. Stassun, S. Khalafinejad, G. Tautvaišienė

Abstract: We present the detection and analysis of the phase curve of KELT-1b at optical wavelengths, analyzing data taken by the Transiting Exoplanet Survey Satellite (TESS). With a mass of ~27 M_J, KELT-1b is a low-mass brown dwarf. Due to the high mass and close proximity of its companion, the host star has a TESS light curve which shows clear ellipsoidal variations. We model the data with a six-componen… ▽ More We present the detection and analysis of the phase curve of KELT-1b at optical wavelengths, analyzing data taken by the Transiting Exoplanet Survey Satellite (TESS). With a mass of ~27 M_J, KELT-1b is a low-mass brown dwarf. Due to the high mass and close proximity of its companion, the host star has a TESS light curve which shows clear ellipsoidal variations. We model the data with a six-component model: secondary eclipse, phase curve accounting for reflected light and thermal emission, Doppler beaming, ellipsoidal variations, stellar activity and the primary transit. We determine the secondary eclipse depth in the TESS bandpass to be 304 +/- 75 parts-per-million (ppm), the most accurate eclipse depth determined so far for KELT-1b. We measure the amplitude of the phase curve to be 128 +/- 27 ppm, with a corresponding eastward offset between the region of maximum brightness and the substellar point of 19.2 +/- 9.6 degrees, in good agreement with Spitzer measurements. We determine day and night brightness temperatures of 3201 +/- 147 K and 1484 +/- 110 K, respectively, slightly higher than those from Spitzer 3.6 and 4.5 micrometer data. A one-dimensional self-consistent atmospheric model can explain the TESS and Spitzer day side brightness temperatures with thermal emission alone and no reflected light. The night side data can be explained by a model with an internal temperature of ~1100 K, which may be related to the inflated radius. The difference between the TESS and Spitzer brightness temperatures can be explained by stronger molecular opacity in the Spitzer bands. On the night side, this opacity is due primarily to CH4 and CO while on the day side it is due to H2-H2 and H2-He collision-induced absorption. △ Less

Submitted 15 March, 2021; v1 submitted 17 June, 2020; originally announced June 2020.

Comments: 18 pages, 12 figures

Neutral Cr and V in the atmosphere of ultra hot jupiter WASP-121 b

Authors: Maya Ben-Yami, Nikku Madhusudhan, Samuel H. C. Cabot, Savvas Constantinou, Anjali Piette, Siddharth Gandhi, Luis Welbanks

Abstract: Ultra hot jupiters (UHJs), giant exoplanets with equilibrium temperatures above 2000 K, are ideal laboratories for studying metal compositions of planetary atmospheres. At these temperatures the thermal dissociation of metal-rich molecules into their constituent elements makes these atmospheres conducive for elemental characterisation. Several elements, mostly ionized metals, have been detected in… ▽ More Ultra hot jupiters (UHJs), giant exoplanets with equilibrium temperatures above 2000 K, are ideal laboratories for studying metal compositions of planetary atmospheres. At these temperatures the thermal dissociation of metal-rich molecules into their constituent elements makes these atmospheres conducive for elemental characterisation. Several elements, mostly ionized metals, have been detected in UHJs recently using high resolution transit spectroscopy. Even though a number of neutral transition metals (e.g., Fe, Ti, V, Cr) are expected to be strong sources of optical/NUV opacity and, hence, influence radiative processes in the lower atmospheres of UHJs, only Fe I has been detected to date. We conduct a systematic search for atomic species in the UHJ WASP-121 b. Using theoretical models we present a metric to predict the atomic species likely to be detectable in such planets with high resolution transmission spectroscopy. We search for the predicted species in observations of WASP-121 b and report the first detections of neutral transition metals Cr I and V I in an exoplanet at 3.6 $σ$ and 4.5 $σ$, respectively. We confirm previous detections of Fe I and Fe II. Whereas Fe II was detected previously in the NUV, we detect it in the optical. We infer that the neutral elements Fe I, V I, and Cr I are present in the lower atmosphere, as predicted by thermochemical equilibrium, while Fe II is a result of photoionisation in the upper atmosphere. Our study highlights the rich chemical diversity of UHJs. △ Less

Submitted 10 June, 2020; originally announced June 2020.

Comments: Accepted for publication in ApJL

Assessing Spectra and Thermal Inversions due to TiO in Hot Jupiter Atmospheres

Authors: Anjali A. A. Piette, Nikku Madhusudhan, Laura K. McKemmish, Siddharth Gandhi, Thomas Masseron, Luis Welbanks

Abstract: Recent detections of thermal inversions in the dayside atmospheres of some hot Jupiters are motivating new avenues to understand the interplay between their temperature structures and other atmospheric conditions. In particular, TiO has long been proposed to cause thermal inversions in hot Jupiters, depending on other factors such as stellar irradiation, C/O, and vertical mixing. TiO also has spec… ▽ More Recent detections of thermal inversions in the dayside atmospheres of some hot Jupiters are motivating new avenues to understand the interplay between their temperature structures and other atmospheric conditions. In particular, TiO has long been proposed to cause thermal inversions in hot Jupiters, depending on other factors such as stellar irradiation, C/O, and vertical mixing. TiO also has spectral features in the optical and near-infrared that have been detected. However, interpretations of TiO signatures rely on the accuracy of TiO opacity used in the models. The recently reported Toto TiO line list provides a new opportunity to investigate these dependencies, which is the goal of the present work. First, we investigate how the Toto line list affects observable transmission and emission spectra of hot Jupiters at low and high resolution. The improvement in the Toto line list compared to a previous line list results in observable differences in the model spectra, particularly in the optical at high resolution. Secondly, we explore the interplay between temperature structure, irradiation and composition with TiO as the primary source of optical opacity, using 1D self-consistent atmospheric models. Among other trends, we find that the propensity for thermal inversions due to TiO peaks at C/O$\sim$0.9, consistent with recent studies. Using these models, we further assess metrics to quantify thermal inversions due to TiO, compared to frequently-used Spitzer photometry, over a range in C/O, irradiation, metallicity, gravity and stellar type. △ Less

Submitted 8 June, 2020; originally announced June 2020.

Comments: Accepted for publication in MNRAS, 19 pages, 14 figures

An Unusual Transmission Spectrum for the Sub-Saturn KELT-11b Suggestive of a Sub-Solar Water Abundance

Authors: Knicole D. Colón, Laura Kreidberg, Luis Welbanks, Michael R. Line, Nikku Madhusudhan, Thomas Beatty, Patrick Tamburo, Kevin B. Stevenson, Avi Mandell, Joseph E. Rodriguez, Thomas Barclay, Eric D. Lopez, Keivan G. Stassun, Daniel Angerhausen, Jonathan J. Fortney, David J. James, Joshua Pepper, John P. Ahlers, Peter Plavchan, Supachai Awiphan, Cliff Kotnik, Kim K. McLeod, Gabriel Murawski, Heena Chotani, Danny LeBrun , et al. (8 additional authors not shown)

Abstract: We present an optical-to-infrared transmission spectrum of the inflated sub-Saturn KELT-11b measured with the Transiting Exoplanet Survey Satellite (TESS), the Hubble Space Telescope (HST) Wide Field Camera 3 G141 spectroscopic grism, and the Spitzer Space Telescope (Spitzer) at 3.6 $μ$m, in addition to a Spitzer 4.5 $μ$m secondary eclipse. The precise HST transmission spectrum notably reveals a l… ▽ More We present an optical-to-infrared transmission spectrum of the inflated sub-Saturn KELT-11b measured with the Transiting Exoplanet Survey Satellite (TESS), the Hubble Space Telescope (HST) Wide Field Camera 3 G141 spectroscopic grism, and the Spitzer Space Telescope (Spitzer) at 3.6 $μ$m, in addition to a Spitzer 4.5 $μ$m secondary eclipse. The precise HST transmission spectrum notably reveals a low-amplitude water feature with an unusual shape. Based on free retrieval analyses with varying molecular abundances, we find strong evidence for water absorption. Depending on model assumptions, we also find tentative evidence for other absorbers (HCN, TiO, and AlO). The retrieved water abundance is generally $\lesssim 0.1\times$ solar (0.001--0.7$\times$ solar over a range of model assumptions), several orders of magnitude lower than expected from planet formation models based on the solar system metallicity trend. We also consider chemical equilibrium and self-consistent 1D radiative-convective equilibrium model fits and find they too prefer low metallicities ($[M/H] \lesssim -2$, consistent with the free retrieval results). However, all the retrievals should be interpreted with some caution since they either require additional absorbers that are far out of chemical equilibrium to explain the shape of the spectrum or are simply poor fits to the data. Finally, we find the Spitzer secondary eclipse is indicative of full heat redistribution from KELT-11b's dayside to nightside, assuming a clear dayside. These potentially unusual results for KELT-11b's composition are suggestive of new challenges on the horizon for atmosphere and formation models in the face of increasingly precise measurements of exoplanet spectra. △ Less

Submitted 20 October, 2020; v1 submitted 11 May, 2020; originally announced May 2020.

Comments: Accepted to The Astronomical Journal. 31 pages, 20 figures, 7 tables

Assessment of Supervised Machine Learning for Atmospheric Retrieval of Exoplanets

Authors: Matthew C. Nixon, Nikku Madhusudhan

Abstract: Atmospheric retrieval of exoplanets from spectroscopic observations requires an extensive exploration of a highly degenerate and high-dimensional parameter space to accurately constrain atmospheric parameters. Retrieval methods commonly conduct Bayesian parameter estimation and statistical inference using sampling algorithms such as Markov Chain Monte Carlo (MCMC) or Nested Sampling. Recently seve… ▽ More Atmospheric retrieval of exoplanets from spectroscopic observations requires an extensive exploration of a highly degenerate and high-dimensional parameter space to accurately constrain atmospheric parameters. Retrieval methods commonly conduct Bayesian parameter estimation and statistical inference using sampling algorithms such as Markov Chain Monte Carlo (MCMC) or Nested Sampling. Recently several attempts have been made to use machine learning algorithms either to complement or replace fully Bayesian methods. While much progress has been made, these approaches are still at times unable to accurately reproduce results from contemporary Bayesian retrievals. The goal of our present work is to investigate the efficacy of machine learning for atmospheric retrieval. As a case study, we use the Random Forest supervised machine learning algorithm which has been applied previously with some success for atmospheric retrieval of the hot Jupiter WASP-12b using its near-infrared transmission spectrum. We reproduce previous results using the same approach and the same semi-analytic models, and subsequently extend this method to develop a new algorithm that results in a closer match to a fully Bayesian retrieval. We combine this new method with a fully numerical atmospheric model and demonstrate excellent agreement with a Bayesian retrieval of the transmission spectrum of another hot Jupiter, HD 209458b. Despite this success, and achieving high computational efficiency, we still find that the machine learning approach is computationally prohibitive for high-dimensional parameter spaces that are routinely explored with Bayesian retrievals with modest computational resources. We discuss the trade offs and potential avenues for the future. △ Less

Submitted 22 April, 2020; originally announced April 2020.

Comments: 13 pages, 14 figures. Accepted for publication in MNRAS

H- and Dissociation in Ultra-hot Jupiters: A Retrieval Case Study of WASP-18b

Authors: Siddharth Gandhi, Nikku Madhusudhan, Avi Mandell

Abstract: Atmospheres of a number of ultra-hot Jupiters (UHJs) with temperatures $\gtrsim$2000~K have been observed recently. Many of these planets show largely featureless thermal spectra in the near-infrared observed with the HST WFC3 spectrograph (1.1-1.7~$μ$m) even though this spectral range contains strong H$_2$O opacity. Recent works have proposed the possibility of H- opacity masking the H$_2$O featu… ▽ More Atmospheres of a number of ultra-hot Jupiters (UHJs) with temperatures $\gtrsim$2000~K have been observed recently. Many of these planets show largely featureless thermal spectra in the near-infrared observed with the HST WFC3 spectrograph (1.1-1.7~$μ$m) even though this spectral range contains strong H$_2$O opacity. Recent works have proposed the possibility of H- opacity masking the H$_2$O feature and/or thermal dissociation of H$_2$O causing its apparent depletion at the high temperatures of UHJs. In this work we test these hypotheses using observations of the exoplanet WASP-18b as a case study. We report detailed atmospheric retrievals of the planet using the HyDRA retrieval code, extended to include the effects of H- opacity and thermal dissociation. We report constraints on the H$_2$O, CO and H- abundances as well as the pressure-temperature profile of the dayside atmosphere for retrievals with and without H-/dissociation for each dataset. We find that the H$_2$O and H- abundances are relatively unconstrained given the featureless WFC3 spectra. We do not conclusively detect H- in the planet contrary to previous studies which used equilibrium models to infer its presence. The constraint on the CO abundance depends on the combination of WFC3 and Spitzer data, ranging from solar to super-solar CO values. We additionally see signs of a thermal inversion from two of the datasets. Our study demonstrates the potential of atmospheric retrievals of UHJs including the effects of H- and thermal dissociation of molecules. △ Less

Submitted 15 April, 2020; originally announced April 2020.

Comments: 19 pages, 10 figures, accepted for publication in The Astronomical Journal

Molecular Cross Sections for High Resolution Spectroscopy of Super Earths, Warm Neptunes and Hot Jupiters

Authors: Siddharth Gandhi, Matteo Brogi, Sergei N. Yurchenko, Jonathan Tennyson, Phillip A. Coles, Rebecca K. Webb, Jayne L. Birkby, Gloria Guilluy, George A. Hawker, Nikku Madhusudhan, Aldo S. Bonomo, Alessandro Sozzetti

Abstract: High resolution spectroscopy (HRS) has been used to detect a number of species in the atmospheres of hot Jupiters. Key to such detections is accurately and precisely modelled spectra for cross-correlation against the R$\gtrsim$20,000 observations. There is a need for the latest generation of opacities which form the basis for high signal-to-noise detections using such spectra. In this study we pre… ▽ More High resolution spectroscopy (HRS) has been used to detect a number of species in the atmospheres of hot Jupiters. Key to such detections is accurately and precisely modelled spectra for cross-correlation against the R$\gtrsim$20,000 observations. There is a need for the latest generation of opacities which form the basis for high signal-to-noise detections using such spectra. In this study we present and make publicly available cross sections for six molecular species, H$_2$O, CO, HCN, CH$_4$, NH$_3$ and CO$_2$ using the latest line lists most suitable for low- and high-resolution spectroscopy. We focus on the infrared (0.95-5~$μ$m) and between 500-1500~K where these species have strong spectral signatures. We generate these cross sections on a grid of pressures and temperatures typical for the photospheres of super Earth, warm Neptunes and hot Jupiters using the latest H$_2$ and He pressure broadening. We highlight the most prominent infrared spectral features by modelling three representative exoplanets, GJ~1214~b, GJ~3470~b and HD~189733~b, which encompass a wide range in temperature, mass and radii. In addition, we verify the line lists for H$_2$O, CO and HCN with previous high resolution observations of hot Jupiters. However, we are unable to detect CH$_4$ with our new cross sections from HRS observations of HD~102195~b. These high accuracy opacities are critical for atmospheric detections with HRS and will be continually updated as new data becomes available. △ Less

Submitted 8 April, 2020; originally announced April 2020.

Comments: 16 pages, 11 figures, accepted for publication in MNRAS