Researcher profile

Nikku Madhusudhan

Nikku Madhusudhan contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author
21works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

21 published item(s)

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Characterising Atmospheres of Cloudy Temperate Mini-Neptunes with JWST

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.

preprint2022arXiv

Constraints on TESS albedos for five hot Jupiters

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.

preprint2022arXiv

Molecular tracers of planet formation in the atmospheres of hot Jupiters

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.

preprint2022arXiv

On Atmospheric Retrievals of Exoplanets with Inhomogeneous Terminators

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.

preprint2020arXiv

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

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.

preprint2020arXiv

Assessment of Supervised Machine Learning for Atmospheric Retrieval of Exoplanets

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.

preprint2020arXiv

Considerations for Atmospheric Retrieval of High-Precision Brown Dwarf Spectra

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.

preprint2020arXiv

Detection of neutral atomic species in the ultra-hot jupiter WASP-121b

The class of ultra-hot Jupiters comprises giant exoplanets undergoing intense irradiation from their host stars. They have proved to be a particularly interesting population for their orbital and atmospheric properties. One such planet, WASP-121~b, is in a highly misaligned orbit close to its Roche limit, and its atmosphere exhibits a thermal inversion. These properties make WASP-121~b an interesting target for additional atmospheric characterization. In this paper, we present analysis of archival high-resolution optical spectra obtained during transits of WASP-121~b. Artifacts from the Rossiter-McLaughlin effect and Center-to-Limb Variation are deemed negligible. However, we discuss scenarios where these effects warrant more careful treatment by modeling the WASP-121 system and varying its properties. We report a new detection of atmospheric absorption from H$α$ in the planet with a transit depth of $1.87\pm0.11\%$. We further confirm a previous detection of the Na I doublet, and report a new detection of Fe I via cross-correlation with a model template. We attribute the H$α$ absorption to an extended Hydrogen atmosphere, potentially undergoing escape, and the Fe I to equilibrium chemistry at the planetary photosphere. These detections help to constrain the composition and chemical processes in the atmosphere of WASP-121~b.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

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

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.

preprint2020arXiv

Statistical Characterization of Hot Jupiter Atmospheres using Spitzer&#39;s Secondary Eclipses

We report 78 secondary eclipse depths for a sample of 36 transiting hot Jupiters observed at 3.6- and 4.5 microns using the Spitzer Space Telescope. Our eclipse results for 27 of these planets are new, and include highly irradiated worlds such as KELT-7b, WASP-87b, WASP-76b, and WASP-64b, and important targets for JWST such as WASP-62b. We find that WASP-62b has a slightly eccentric orbit e cos(omega) = 0.00614+/- 0.00064, and we confirm the eccentricity of HAT-P-13b and WASP-14b. The remainder are individually consistent with circular orbits, but we find statistical evidence for eccentricity increasing with orbital period in our range from 1 to 5 days. Our day-side brightness temperatures for the planets yield information on albedo and heat redistribution, following Cowan and Agol (2011). Planets having maximum day side temperatures exceeding ~ 2200K are consistent with zero albedo and distribution of stellar irradiance uniformly over the day-side hemisphere. Our most intriguing result is that we detect a systematic difference between the emergent spectra of these hot Jupiters as compared to blackbodies. The ratio of observed brightness temperatures, Tb(4.5)/Tb(3.6), increases with equilibrium temperature by 100 +/- 24 parts-per-million per Kelvin, over the entire temperature range in our sample (800K to 2500K). No existing model predicts this trend over such a large range of temperature. We suggest that this may be due to a structural difference in the atmospheric temperature profile between the real planetary atmospheres as compared to models.

preprint2020arXiv

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

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.

preprint2020arXiv

The interior and atmosphere of the habitable-zone exoplanet K2-18b

Exoplanets orbiting M dwarfs present a valuable opportunity for their detection and atmospheric characterisation. This is evident from recent inferences of H$_2$O in such atmospheres, including that of the habitable-zone exoplanet K2-18b. With a bulk density between Earth and Neptune, K2-18b may be expected to possess a H/He envelope. However, the extent of such an envelope and the thermodynamic conditions of the interior remain unexplored. In the present work, we investigate the atmospheric and interior properties of K2-18b based on its bulk properties and its atmospheric transmission spectrum. We constrain the atmosphere to be H$_2$-rich with a H$_2$O volume mixing ratio of $0.02-14.8$%, consistent with previous studies, and find a depletion of CH$_4$ and NH$_3$, indicating chemical disequilibrium. We do not conclusively detect clouds/hazes in the observable atmosphere. We use the bulk parameters and retrieved atmospheric properties to constrain the internal structure and thermodynamic conditions in the planet. The constraints on the interior allow multiple scenarios between rocky worlds with massive H/He envelopes and water worlds with thin envelopes. We constrain the mass fraction of the H/He envelope to be $\lesssim 6$%; spanning $\lesssim 10^{-5}$ for a predominantly water world to $\sim6$% for a pure iron interior. The thermodynamic conditions at the surface of the H$_2$O layer range from the super-critical to liquid phases, with a range of solutions allowing for habitable conditions on K2-18b. Our results demonstrate that the potential for habitable conditions is not necessarily restricted to Earth-like rocky exoplanets.

preprint2019arXiv

HyDRA-H: Simultaneous Hybrid Retrieval of Exoplanetary Emission Spectra

High-resolution Doppler spectroscopy has been used to detect several chemical species in exoplanetary atmospheres. Such detections have traditionally relied on cross correlation of observed spectra against spectral model templates, an approach that is successful for detecting chemical species but not optimised for constraining abundances. Recent work has explored ways to perform atmospheric retrievals on high-resolution spectra (HRS) and combine them with retrievals routinely performed for low-resolution spectra (LRS) by developing a mapping from the cross correlation function to a likelihood metric. We build upon previous studies and report HyDRA-H, a hybrid retrieval code for simultaneous analysis of low- and high- resolution thermal emission spectra of exoplanets in a fully Bayesian approach. We demonstrate HyDRA-H on the hot Jupiter HD 209458b as a case study. We validate our HRS retrieval capability by confirming previous results and report a simultaneous hybrid retrieval using both HRS and LRS data. The LRS data span the HST WFC3 (1.1-1.7 $μ$m) and Spitzer photometry (IRAC 3.6-8 $μ$m) bands, while the HRS data were obtained with CRIRES on VLT at 2.3 $μ$m. The constraints on the composition and temperature profiles for the hybrid retrieval are more stringent than retrievals with either LRS or HRS datasets individually. We retrieve abundances of $\log(\mathrm{H_2O)} = -4.11^{+0.91}_{-0.30}$ and $\log(\mathrm{CO}) = {-2.16}^{+0.99}_{-0.47}$, and $\mathrm{C/O} = 0.99^{+0.01}_{-0.02}$, consistent with previous works. We constrain the photospheric temperature to be $1498^{+216}_{-57}$ K, consistent with the equilibrium temperature. Our results demonstrate the significant advantages of hybrid retrievals by combining strengths of both HRS and LRS observations which probe complementary aspects of exoplanetary atmospheres.

preprint2019arXiv

JexoSim: A time domain simulator of exoplanet transit spectroscopy with JWST

The James Webb Space Telescope (JWST) will perform exoplanet transit spectroscopy in the coming decade promising transformative science results. All four instruments on board can be used for this technique which reconstructs the atmospheric transmission or emission spectrum of an exoplanet from wavelength-dependent light curve measurements. Astrophysical and instrumental noise and systematics can affect the precision and accuracy of the final spectrum, and hence, the atmospheric properties derived from the spectrum. Correlated noise and time-dependent systematics that can bias the measured signal must be accounted for in the final uncertainties. However quantifying these effects can be difficult with real data or simple analytic tools. Existing publicly-available simulators for JWST do not adequately simulate complex time domain processes on exoplanetary transit observations. We report JexoSim, a dedicated time domain simulator for JWST including all four instruments for exoplanet transit spectroscopy. JexoSim models both the astrophysics and the instrument, generating 2-D images in simulated time akin to a real observation. JexoSim can capture correlated noise and systematic biases on the light curve giving it great versatility. Potential applications of JexoSim include performance testing of JWST instruments, assessing science return, and testing data reduction pipelines. We describe JexoSim, validate it against other simulators, and present examples of its utility.

preprint2019arXiv

Mass-Metallicity Trends in Transiting Exoplanets from Atmospheric Abundances of H$_2$O, Na, and K

Atmospheric compositions can provide powerful diagnostics of formation and migration histories of planetary systems. We investigate constraints on atmospheric abundances of H$_2$O, Na, and K, in a sample of transiting exoplanets using latest transmission spectra and new H$_2$ broadened opacities of Na and K. Our sample of 19 exoplanets spans from cool mini-Neptunes to hot Jupiters, with equilibrium temperatures between $\sim$300 and 2700 K. Using homogeneous Bayesian retrievals we report atmospheric abundances of Na, K, and H$_2$O, and their detection significances, confirming 6 planets with strong Na detections, 6 with K, and 14 with H$_2$O. We find a mass-metallicity trend of increasing H$_2$O abundances with decreasing mass, spanning generally substellar values for gas giants and stellar/superstellar for Neptunes and mini-Neptunes. However, the overall trend in H$_2$O abundances, from mini-Neptunes to hot Jupiters, is significantly lower than the mass-metallicity relation for carbon in the solar system giant planets and similar predictions for exoplanets. On the other hand, the Na and K abundances for the gas giants are stellar or superstellar, consistent with each other, and generally consistent with the solar system metallicity trend. The H$_2$O abundances in hot gas giants are likely due to low oxygen abundances relative to other elements rather than low overall metallicities, and provide new constraints on their formation mechanisms. The differing trends in the abundances of species argue against the use of chemical equilibrium models with metallicity as one free parameter in atmospheric retrievals, as different elements can be differently enhanced.

preprint2011arXiv

Possible thermochemical disequilibrium in the atmosphere of the exoplanet GJ 436b

The nearby extrasolar planet GJ 436b--which has been labelled as a &#39;hot Neptune&#39;--reveals itself by the dimming of light as it crosses in front of and behind its parent star as seen from Earth. Respectively known as the primary transit and secondary eclipse, the former constrains the planet&#39;s radius and mass, and the latter constrains the planet&#39;s temperature and, with measurements at multiple wavelengths, its atmospheric composition. Previous work using transmission spectroscopy failed to detect the 1.4-μm water vapour band, leaving the planet&#39;s atmospheric composition poorly constrained. Here we report the detection of planetary thermal emission from the dayside of GJ 436b at multiple infrared wavelengths during the secondary eclipse. The best-fit compositional models contain a high CO abundance and a substantial methane (CH4) deficiency relative to thermochemical equilibrium models for the predicted hydrogen-dominated atmosphere. Moreover, we report the presence of some H2O and traces of CO2. Because CH4 is expected to be the dominant carbon-bearing species, disequilibrium processes such as vertical mixing and polymerization of methane into substances such as ethylene may be required to explain the hot Neptune&#39;s small CH4-to-CO ratio, which is at least 10^5 times smaller than predicted.