Search

We report JWST NIRSpec/G395H observations of TOI-1685 b, a hot rocky super-Earth orbiting an M2.5V star, during a full orbit. We obtain transmission and emission spectra of the planet and characterize the properties of the phase curve, including its amplitude and offset. The transmission spectrum rules out clear H$_2$-dominated atmospheres, while secondary atmospheres (made of water, methane, or carbon dioxide) cannot be statistically distinguished from a flat line. The emission spectrum is featureless and consistent with a blackbody-like brightness temperature, helping rule out thick atmospheres with high mean molecular weight. Collecting all evidence, the properties of TOI-1685 b are consistent with a blackbody with no heat redistribution and a low albedo, with a dayside brightness temperature 0.98$\pm$0.07 times that of a perfect blackbody in the NIRSpec NRS2 wavelength range (3.823-5.172 um). Our results add to the growing number of seemingly airless M-star rocky planets, thus constraining the location of the "Cosmic Shoreline". Three independent data reductions have been carried out, all showing a high-amplitude correlated noise component in the white and spectroscopic light curves. The correlated noise properties are different between the NRS1 and NRS2 detectors - importantly the timescales of the strongest components (4.5 hours and 2.5 hours, respectively) - suggesting the noise is from instrumental rather than astrophysical origins. We encourage the community to look into the systematics of NIRSpec for long time-series observations., Comment: 26 pages, 19 figures, 8 tables. Submitted to AAS Journals

We investigate the impact of using multipoint p-T profiles of varying complexity on the retrieval of synthetically generated hot Jupiter transmission spectra modelled after state-of-the-art observations of the hot Jupiter WASP-39~b with JWST. We perform homogenised atmospheric retrievals with the TauREx retrieval framework on a sample of synthetically generated transmission spectra, accounting for varying cases of underlying p-T profiles, cloud-top pressures, and expected noise levels. These retrievals are performed using a fixed-pressure multipoint p-T prescription with increasing complexity, ranging from isothermal to an eleven-point profile. We evaluate the performance of the retrievals based on the Bayesian model evidence, and the accuracy of the retrievals compared to the known input parameters. We find that performing atmospheric retrievals using an isothermal prescription for the pressure-temperature profile consistently results in wrongly retrieved atmospheric parameters when compared to the known input parameters. For an underlying p-T profile with a fully positive lapse rate, we find that a two-point profile is sufficient to retrieve the known atmospheric parameters, while under the presence of an atmospheric temperature inversion, we find that a more complex profile is necessary. Our investigation shows that, for a data quality scenario mirroring state-of-the-art observations of a hot Jupiter with JWST, an isothermal p-T prescription is insufficient to correctly retrieve the known atmospheric parameters. We find a model complexity preference dependent on the underlying pressure-temperature structure, but argue that a p-T prescription on the complexity level of a four-point profile should be preferred. This represents the overlap between the lowest number of free parameters and highest model preference in the cases investigated in this work., Comment: 16 pages, 14 figures, accepted for publication in A&A

The distribution of different types of atmospheres and surfaces on rocky planets is one of the major questions in exoplanet astronomy, but there are currently no published unambiguous detections of atmospheres on any rocky exoplanets. The MIRI instrument on JWST can measure thermal emission from tidally locked rocky exoplanets orbiting small, cool stars. This emission is a function of their surface and atmospheric properties, potentially allowing the detection of atmospheres. One technique is to measure day-side emission to search for lower thermal emission than expected for a black-body planet due to atmospheric absorption features. Another technique is to measure phase curves of thermal emission to search for night-side emission due to atmospheric heat redistribution. In this work we compare strategies for detecting atmospheres on rocky exoplanets using these techniques. We simulate secondary eclipse and phase curve observations in the MIRI F1500W and F1280W filters, for a range of surfaces and atmospheres on thirty exoplanets selected for their F1500W signal-to-noise ratio. Our results show that secondary eclipse observations are highly degenerate between surfaces and atmospheres, given the wide range of potential surface albedos. We also show that thick atmospheres can support emission consistent with a black-body planet in these filters. These two results make it difficult to unambiguously detect or rule out atmospheres using their photometric day-side emission, except in a subset of CO$_{2}$-dominated atmospheres. We suggest that an F1500W phase curve could instead be observed for a similar sample of planets, allowing the unambiguous detection of atmospheres by night-side emission., Comment: Updated version for journal submission

We present ground-based photometric observations of secondary eclipses of the hottest known planet KELT-9b using MuSCAT2 and Sinistro. We detect secondary eclipse signals in $i$ and $z_{\rm s}$ with eclipse depths of $373^{+74}_{-75}$ ppm and $638^{+199}_{-178}$, respectively. We perform an atmospheric retrieval on the emission spectrum combined with the data from HST/WFC3, Spitzer, TESS, and CHEOPS to obtain the temperature profile and chemical abundances, including TiO and VO, which have been thought to produce temperature inversion structures in the dayside of ultra-hot Jupiters. While we confirm a strong temperature inversion structure, we find low abundances of TiO and VO with mixing ratios of $\rm{log(TiO)}=-7.80^{+0.15}_{-0.30}$ and $\rm{log(VO)}=-9.60^{+0.64}_{-0.57}$, respectively. The low abundances of TiO and VO are consistent with theoretical predictions for such an ultra-hot atmosphere. In such low abundances, TiO and VO have little effect on the temperature structure of the atmosphere. The abundance of ${\rm e}^{-}$, which serves as a proxy for ${\rm H}^{-}$ ions in this study, is found to be high, with $\rm{log(e^-)}=-4.89\pm{0.06}$. These results indicate that the temperature inversion in KELT-9 b's dayside atmosphere is likely not caused by TiO/VO, but rather by the significant abundance of ${\rm H}^{-}$ ions. The best-fit model cannot fully explain the observed spectrum, and chemical species not included in the retrieval may introduce modeling biases. Future observations with broader wavelength coverage and higher spectral resolution are expected to provide more accurate diagnostics on the presence and abundances of TiO/VO. These advanced observations will overcome the limitations of current data from HST and photometric facilities, which are constrained by narrow wavelength coverage and instrumental systematics., Comment: 16 pages, 9 figures, 4 tables, Accepted for publication in PASJ

Conventional atmospheric retrieval codes are designed to extract information, such as chemical abundances, thermal structures and cloud properties, from fully "reduced" spectra obtained during transit or eclipse. Reduced spectra, however, are assembled by fitting a series of simplified light-curves to time series observations, wavelength-by-wavelength. Thus, spectra are a post-processed summary statistics of the original data, which by definition does not encode all the available information (i.e., astrophysical signal, model covariance, instrumental noise). Here, we explore an alternative inversion strategy where the atmospheric retrieval is performed on the light-curve directly -- i.e., closer to the data. This method is implemented in ExPLOR (EXoplanet Panchromatic Light-curve Observation and Retrieval), a novel atmospheric retrieval code inheriting from the TauREx project. By explicitly considering time in the model, ExPLOR naturally handles transit, eclipse, phase-curve and other complex geometries for transiting exoplanets. In this paper, we have validated this new technique by inverting simulated panchromatic light-curves. The model was tested on realistic simulations of a WASP-43b like exoplanet as observed with the James Webb Space Telescope (JWST) and Ariel telescope. By comparing our panchromatic light-curve approach against conventional spectral retrievals on mock scenarios, we have identified key breaking points in information and noise propagation when employing past literature techniques. Throughout the paper, we discuss the importance of developing "closer-to-data" approaches such as the method presented in this work, and highlight the inevitable increase in model complexity and computing requirements associated with the recent JWST revolution., Comment: Accepted to AJ, 22 pages, 11 figures

We present a comprehensive atmospheric retrieval study of the hot Jupiter WASP-77A\,b using eclipse observations from the Hubble Space Telescope (HST) and JWST. Using atmospheric retrievals, the spectral features of H$_2$O, CO, and TiO are identified, with volume mixing ratios estimated at log$_{\rm 10}$(VMR) = -4.40$^{+0.14}_{-0.11}$, -4.44$^{+0.34}_{-0.28}$, and -6.40$^{+0.22}_{-0.23}$, respectively. We derive the atmospheric carbon-to-oxygen ratio -- a key planetary formation tracer -- to be C/O = 0.54$\pm$0.12, which is consistent with both the stellar host value and previous studies of the planet's atmosphere, suggesting a relatively close-in formation. Computing other elemental ratios (i.e., C/H, O/H, and Ti/H), we conclude that the general enrichment of the atmosphere (i.e., metallicity) is sub-stellar, is depleted in C and O, but that Ti appears slightly super-stellar. A low C and O content could be obtained, in combination with a stellar C/O ratio, if the planet formed outside of the CO$_2$ snow line before migrating inwards. Meanwhile, a super-stellar Ti/H could be obtained by late contamination from refractory rich planetesimals. While broadly in agreement with previous works, we do find some differences and discuss these while also highlighting the need for homogeneous analyses when comparative exoplanetology is conducted., Comment: Accepted for publication in ApJL (December 2023)

We present a comprehensive analysis of the Hubble Space Telescope observations of the atmosphere of WASP-121 b, a ultra-hot Jupiter. After reducing the transit, eclipse, and phase-curve observations with a uniform methodology and addressing the biases from instrument systematics, sophisticated atmospheric retrievals are used to extract robust constraints on the thermal structure, chemistry, and cloud properties of the atmosphere. Our analysis shows that the observations are consistent with a strong thermal inversion beginning at ~0.1 bar on the dayside, solar to subsolar metallicity Z (i.e., -0.77 < log(Z) < 0.05), and super-solar C/O ratio (i.e., 0.59 < C/O < 0.87). More importantly, utilizing the high signal-to-noise ratio and repeated observations of the planet, we identify the following unambiguous time-varying signals in the data: i) a shift of the putative hotspot offset between the two phase-curves and ii) varying spectral signatures in the transits and eclipses. By simulating the global dynamics of WASP-121 b atmosphere at high-resolution, we show that the identified signals are consistent with quasi-periodic weather patterns, hence atmospheric variability, with signatures at the level probed by the observations (~5% to ~10%) that change on a timescale of ~5 planet days; in the simulations, the weather patterns arise from the formation and movement of storms and fronts, causing hot (as well as cold) patches of atmosphere to deform, separate, and mix in time., Comment: Accepted for publication in ApJS. 43 pages, 31 figures, 2 animations (available online at the journal)

The Ariel Space Mission aims to observe a diverse sample of exoplanet atmospheres across a wide wavelength range of 0.5 to 7.8 microns. The observations are organized into four Tiers, with Tier 1 being a reconnaissance survey. This Tier is designed to achieve a sufficient signal-to-noise ratio (S/N) at low spectral resolution in order to identify featureless spectra or detect key molecular species without necessarily constraining their abundances with high confidence. We introduce a P-statistic that uses the abundance posteriors from a spectral retrieval to infer the probability of a molecule's presence in a given planet's atmosphere in Tier 1. We find that this method predicts probabilities that correlate well with the input abundances, indicating considerable predictive power when retrieval models have comparable or higher complexity compared to the data. However, we also demonstrate that the P-statistic loses representativity when the retrieval model has lower complexity, expressed as the inclusion of fewer than the expected molecules. The reliability and predictive power of the P-statistic are assessed on a simulated population of exoplanets with H2-He dominated atmospheres, and forecasting biases are studied and found not to adversely affect the classification of the survey., Comment: 51 pages, 15 figures, accepted in Experimental Astronomy

Aims. We analyse unpublished Spitzer observations of the thermal phase-curve of WASP-121 b, a benchmark ultra-hot Jupiter. Methods. We adopted the wavelet pixel-independent component analysis technique to remove challenging instrumental systematic effects in these datasets and we fit them simultaneously with parametric light-curve models. We also performed phase-curve retrievals to better understand the horizontal and vertical thermal structure of the planetary atmosphere. Results. We measured planetary brightness temperatures of $\sim$2700\,K (dayside) and $\sim$700--1100\,K (nightside), along with modest peak offsets of 5.9$^{\circ} \pm$1.6 (3.6\,$\mu$m) and 5.0$^{\circ}$$_{-3.1}^{+3.4}$ (4.5\,$\mu$m) after mid-eclipse. These results suggest inefficient heat redistribution in the atmosphere of WASP-121 b. The inferred atmospheric Bond albedo and circulation efficiency align well with observed trends for hot giant exoplanets. Interestingly, the measured peak offsets correspond to a westward hot spot, which has rarely been observed. We also report consistent transit depths at 3.6 and 4.5\,$\mu$m, along with updated geometric and orbital parameters. Finally, we compared our Spitzer results with previous measurements, including recent JWST observations. Conclusions. We extracted new information on the thermal properties and dynamics of an exoplanet atmosphere from an especially problematic dataset. This study probes the reliability of exoplanet phase-curve parameters obtained from Spitzer observations when state-of-the-art pipelines are adopted to remove the instrumental systematic effects. It demonstrates that Spitzer phase-curve observations provide a useful baseline for comparison with JWST observations, and shows the increase in parameters precision achieved with the newer telescope., Comment: 14 pages, 10 figures

We present an atmospheric analysis of LTT 9779 b, a rare planet situated in the hot Neptune desert, that has been observed with HST WFC3 G102 and G141. The combined transmission spectrum, which covers 0.8 - 1.6 $\mu$m, shows a gradual increase in transit depth with wavelength. Our preferred atmospheric model shows evidence for H$_{\rm 2}$O, CO$_{\rm 2}$ and FeH with a significance of 3.1 $\sigma$, 2.4 $\sigma$ and 2.1 $\sigma$, respectively. In an attempt to constrain the rate of atmospheric escape for this planet, we search for the 1.083 $\mu$m Helium line in the G102 data but find no evidence of excess absorption that would indicate an escaping atmosphere using this tracer. We refine the orbital ephemerides of LTT 9779 b using our HST data and observations from TESS, searching for evidence of orbital decay or apsidal precession, which is not found. The phase-curve observation of LTT 9779 b with JWST NIRISS should provide deeper insights into the atmosphere of this planet and the expected atmospheric escape might be detected with further observations concentrated on other tracers such as Lyman $\alpha$., Comment: Accepted for publication in AJ

Studying chemistry and chemical composition is fundamental to go back to formation history of planetary systems. We propose here to have another look at five targets to better determine their composition and the chemical mechanisms that take place in their atmospheres. We present a re-analysis of five Hot Jupiters, combining multiple instruments and using Bayesian retrieval methods. We compare different combinations of molecules present in the simulated atmosphere, different chemistry types as well as different clouds parametrization. As a consequence of recent studies questioning the detection of Na and K in the atmosphere of HD 209458b as being potentially contaminated by stellar lines when present, we study the impact on other retrieval parameters of misinterpreting the presence of these alkali species. We use spatially scanned observations from the grisms G102 and G141 of the WFC3 on HST, with a wavelength coverage of $\sim$0.8 to $\sim$1.7 microns. We analyse these data with the publicly available Iraclis pipeline. We added to our datasets STIS observations to increase our wavelength coverage from $\sim$0.4 to $\sim$1.7 microns. We then performed a Bayesian retrieval analysis with the open-source TauREx using a nested sampling algorithm. We explore the influence of including Na and K on the retrieval of the molecules from the atmosphere. Our data re-analysis and Bayesian retrieval are consistent with previous studies but we find small differences in the retrieved parameters. After all, Na and K has no significant impact on the properties of the planet atmospheres. Therefore, we present here our new best-fit models, taking into account molecular abundances varying freely and equilibrium chemistry. This work is a preparation for a future addition of more sophisticated representation of chemistry taking into account disequilibrium effects such as vertical mixing and photochemistry., Comment: 19 pages, 14 figures

Retrieval methods are a powerful analysis technique for modelling exoplanetary atmospheres by estimating the bulk physical and chemical properties that combine in a forward model to best-fit an observed spectrum, and they are increasingly being applied to observations of directly-imaged exoplanets. We have adapted TauREx3, the Bayesian retrieval suite, for the analysis of near-infrared spectrophotometry from directly-imaged gas giant exoplanets and brown dwarfs. We demonstrate TauREx3's applicability to sub-stellar atmospheres by presenting results for brown dwarf benchmark GJ 570D which are consistent with previous retrieval studies, whilst also exhibiting systematic biases associated with the presence of alkali lines. We also present results for the cool exoplanet 51 Eri b, the first application of a free chemistry retrieval analysis to this object, using spectroscopic observations from GPI and SPHERE. While our retrieval analysis is able to explain spectroscopic and photometric observations without employing cloud extinction, we conclude this may be a result of employing a flexible temperature-pressure profile which is able to mimic the presence of clouds. We present Bayesian evidence for an ammonia detection with a 2.7$\sigma$ confidence, the first indication of ammonia in an exoplanetary atmosphere. This is consistent with this molecule being present in brown dwarfs of a similar spectral type. We demonstrate the chemical similarities between 51 Eri b and GJ 570D in relation to their retrieved molecular abundances. Finally, we show that overall retrieval conclusions for 51 Eri b can vary when employing different spectral data and modelling components, such as temperature-pressure and cloud structures.

The chromatic contamination that arises from photospheric heterogeneities e.g. spots and faculae on the host star presents a significant noise source for exoplanet transmission spectra. If this contamination is not corrected for, it can introduce substantial bias in our analysis of the planetary atmosphere. We utilise two stellar models of differing complexity, StARPA and ASteRA, to explore the biases introduced by stellar contamination in retrieval under differing degrees of stellar activity. We use the retrieval framework TauREx3 and a grid of 27 synthetic, spot-contaminated transmission spectra to investigate potential biases and to determine how complex our stellar models must be in order to accurately extract the planetary parameters from transmission spectra. The input observation is generated using the more complex model (StARPA), in which the spot latitude is an additional, fixable parameter. This observation is then fed into a combined stellar-planetary retrieval which contains a simplified stellar model (ASteRA). Our results confirm that the inclusion of stellar activity parameters in retrieval minimises bias under all activity regimes considered. ASteRA performs very well under low to moderate activity conditions, retrieving the planetary parameters with a high degree of accuracy. For the most active cases, characterised by larger, higher temperature contrast spots, some minor residual bias remains due to ASteRA neglecting the interplay between the spot and the limb darkening effect. As a result of this, we find larger errors in retrieved planetary parameters for central spots (0 degrees) and those found close to the limb (60 degrees) than those at intermediate latitudes (30 degrees)., Comment: 34 pages, 20 figures, accepted for publication in ApJ

In this paper, we present YunMa, an exoplanet cloud simulation and retrieval package, which enables the study of cloud microphysics and radiative properties in exoplanetary atmospheres. YunMa simulates the vertical distribution and sizes of cloud particles and their corresponding scattering signature in transit spectra. We validated YunMa against results from the literature. When coupled to the TauREx 3 platform, an open Bayesian framework for spectral retrievals, YunMa enables the retrieval of the cloud properties and parameters from transit spectra of exoplanets. The sedimentation efficiency ($f_{\mathrm{sed}}$), which controls the cloud microphysics, is set as a free parameter in retrievals. We assess the retrieval performances of YunMa through 28 instances of a K2-18 b-like atmosphere with different fractions of H$_2$/He and N$_2$, and assuming water clouds. Our results show a substantial improvement in retrieval performances when using YunMa instead of a simple opaque cloud model and highlight the need to include cloud radiative transfer and microphysics to interpret the next-generation data for exoplanet atmospheres. This work also inspires instrumental development for future flagships by demonstrating retrieval performances with different data quality., Comment: 24 pages, 12 figures, accepted in ApJ

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

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

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

We present the analysis of the atmospheres of 70 gaseous extrasolar planets via transit spectroscopy with Hubble's Wide Field Camera 3 (WFC3). For over half of these, we statistically detect spectral modulation which our retrievals attribute to molecular species. Among these, we use Bayesian Hierarchical Modelling to search for chemical trends with bulk parameters. We use the extracted water abundance to infer the atmospheric metallicity and compare it to the planet's mass. We also run chemical equilibrium retrievals, fitting for the atmospheric metallicity directly. However, although previous studies have found evidence of a mass-metallicity trend, we find no such relation within our data. For the hotter planets within our sample, we find evidence for thermal dissociation of dihydrogen and water via the H$^-$ opacity. We suggest that the general lack of trends seen across this population study could be due to i) the insufficient spectral coverage offered by HST WFC3 G141, ii) the lack of a simple trend across the whole population, iii) the essentially random nature of the target selection for this study or iv) a combination of all the above. We set out how we can learn from this vast dataset going forward in an attempt to ensure comparative planetology can be undertaken in the future with facilities such as JWST, Twinkle and Ariel. We conclude that a wider simultaneous spectral coverage is required as well as a more structured approach to target selection., Comment: Accepted for publication in ApJS

We explore how finite integration time or temporal binning can affect the analysis of exoplanet phase-curves. We provide analytical formulae to account for this effect or, if neglected, to estimate the potential biases in the retrieved parameters. As expected, due to their smoother variations over longer time-scales, phase curves can be binned more heavily than transits without causing severe biases. In the simplest case of a sinusoidal phase curve with period $P$, the integration time $\Delta t$ reduces its amplitude by the scaling factor $\text{sinc}{ \left ( \pi \Delta t / P \right ) }$, without altering its phase or shape. We also provide formulae to predict reasonable parameter error bars from phase-curve observations. Our findings are tested with both synthetic and real datasets, including unmodelled astrophysical signals and/or instrumental systematic effects. Tests with the Spitzer data show that binning can affect the best-fitting parameters beyond predictions, due to the correction of high-frequency correlated noise. Finally, we summarize key guidelines for speeding up the analysis of exoplanet phase curves without introducing significant biases in the retrieved parameters., Comment: Published in MNRAS (here with extended acknowledgments)

We introduce a new Python 1D chemical kinetic code FRECKLL (Full and Reduced Exoplanet Chemical Kinetics distiLLed) to evolve large chemical networks efficiently. FRECKLL employs `distillation' in computing the reaction rates, which minimizes the error bounds to the minimum allowed by double precision values ($\epsilon \leq 10^{-15}$). Compared to summation of rates with traditional algorithms like pairwise summation, distillation provides a tenfold reduction in solver time for both full and reduced networks. Both the full and reduced Venot2020 networks are packaged in FRECKLL as well as a TauREx 3.1 plugin for usage in forward modelling and retrievals of exoplanet atmospheres. We present TauREx retrievals performed on a simulated HD189733 JWST spectra using the full and reduced Venot2020 chemical networks and demonstrate the viability of total disequilibrium chemistry retrievals and the ability for JWST to detect disequilibrium processes., Comment: 21 pages, 11 figures

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

The study of extra-solar planets, or simply, exoplanets, planets outside our own Solar System, is fundamentally a grand quest to understand our place in the Universe. Discoveries in the last two decades have re-defined our understanding of planets, and helped us comprehend the uniqueness of our very own Earth. In recent years the focus has shifted from planet detection to planet characterisation, where key planetary properties are inferred from telescope observations using Monte Carlo-based methods. However, the efficiency of sampling-based methodologies is put under strain by the high-resolution observational data from next generation telescopes, such as the James Webb Space Telescope and the Ariel Space Mission. We are delighted to announce the acceptance of the Ariel ML Data Challenge 2022 as part of the NeurIPS competition track. The goal of this challenge is to identify a reliable and scalable method to perform planetary characterisation. Depending on the chosen track, participants are tasked to provide either quartile estimates or the approximate distribution of key planetary properties. To this end, a synthetic spectroscopic dataset has been generated from the official simulators for the ESA Ariel Space Mission. The aims of the competition are three-fold. 1) To offer a challenging application for comparing and advancing conditional density estimation methods. 2) To provide a valuable contribution towards reliable and efficient analysis of spectroscopic data, enabling astronomers to build a better picture of planetary demographics, and 3) To promote the interaction between ML and exoplanetary science. The competition is open from 15th June and will run until early October, participants of all skill levels are more than welcomed!, Comment: 13 pages, Accepted in the Conference on Neural Information Processing Systems 2022 (NeurIPS 2022), Competition hosted on https://www.ariel-datachallenge.space/

This is an exciting era for exo-planetary exploration. The recently launched JWST, and other upcoming space missions such as Ariel, Twinkle and ELTs are set to bring fresh insights to the convoluted processes of planetary formation and evolution and its connections to atmospheric compositions. However, with new opportunities come new challenges. The field of exoplanet atmospheres is already struggling with the incoming volume and quality of data, and machine learning (ML) techniques lands itself as a promising alternative. Developing techniques of this kind is an inter-disciplinary task, one that requires domain knowledge of the field, access to relevant tools and expert insights on the capability and limitations of current ML models. These stringent requirements have so far limited the developments of ML in the field to a few isolated initiatives. In this paper, We present the Atmospheric Big Challenge Database (ABC Database), a carefully designed, organised and publicly available database dedicated to the study of the inverse problem in the context of exoplanetary studies. We have generated 105,887 forward models and 26,109 complementary posterior distributions generated with Nested Sampling algorithm. Alongside with the database, this paper provides a jargon-free introduction to non-field experts interested to dive into the intricacy of atmospheric studies. This database forms the basis for a multitude of research directions, including, but not limited to, developing rapid inference techniques, benchmarking model performance and mitigating data drifts. A successful application of this database is demonstrated in the NeurIPS Ariel ML Data Challenge 2022., Comment: 21 pages, Accepted in RASTI, challenge accepted in NeurIPS 2022. Database link: https://doi.org/10.5281/zenodo.6770103

Current endeavours in exoplanet characterisation rely on atmospheric retrieval to quantify crucial physical properties of remote exoplanets from observations. However, the scalability and efficiency of the technique are under strain with increasing spectroscopic resolution and forward model complexity. The situation becomes more acute with the recent launch of the James Webb Space Telescope and other upcoming missions. Recent advances in Machine Learning provide optimisation-based Variational Inference as an alternative approach to perform approximate Bayesian Posterior Inference. In this investigation we combined Normalising Flow-based neural network with our newly developed differentiable forward model, Diff-Tau, to perform Bayesian Inference in the context of atmospheric retrieval. Using examples from real and simulated spectroscopic data, we demonstrated the superiority of our proposed framework: 1) Training Our neural network only requires a single observation; 2) It produces high-fidelity posterior distributions similar to sampling-based retrieval and; 3) It requires 75% less forward model computation to converge. 4.) We performed, for the first time, Bayesian model selection on our trained neural network. Our proposed framework contribute towards the latest development of a neural-powered atmospheric retrieval. Its flexibility and speed hold the potential to complement sampling-based approaches in large and complex data sets in the future., Comment: 16 pages, 4 figures, Published in ApJ

Population studies of exoplanets are key to unlocking their statistical properties. So far the inferred properties have been mostly limited to planetary, orbital and stellar parameters extracted from, e.g., Kepler, radial velocity, and GAIA data. More recently an increasing number of exoplanet atmospheres have been observed in detail from space and the ground. Generally, however, these atmospheric studies have focused on individual planets, with the exception of a couple of works which have detected the presence of water vapor and clouds in populations of gaseous planets via transmission spectroscopy. Here, using a suite of retrieval tools, we analyse spectroscopic and photometric data of 25 hot Jupiters, obtained with the Hubble and Spitzer Space Telescopes via the eclipse technique. By applying the tools uniformly across the entire set of 25 planets, we extract robust trends in the thermal structure and chemical properties of hot Jupiters not obtained in past studies. With the recent launch of JWST and the upcoming missions Twinkle, and Ariel, population based studies of exoplanet atmospheres, such as the one presented here, will be a key approach to understanding planet characteristics, formation, and evolution in our galaxy., Comment: 66 pages, 23 figures, 7 tables. Published in The Astrophysical Journal Supplement Series

The discovery of super-Earth and mini-Neptune exoplanets means that atmospheric signals from low-mass, temperate exoplanets are being increasingly studied. The signal acquired as the planet transits its host star, known as the transit depth, is smaller for these planets and, as such, more difficult to analyze. The launch of the space telescopes James Webb (JWST) & Ariel will give rise to an explosion in the quality and quantity of spectroscopic data available for an unprecedented number of exoplanets in our galaxy. Accurately extracting the information content, thereby permitting atmospheric science, of such data-sets will require robust models and techniques. We present here the analysis of simulated transmission spectra for water-rich atmospheres, giving evidence for non-negligible differences in simulated transit depths when self-broadening of H$_2$O is correctly accounted for, compared with the currently typically accepted standard of using H$_2$ and He-broadened cross-sections. Our case-study analysis is carried out on two super-Earths, focusing on water-based atmospheres, ranging from H$_2$-rich to H$_2$O-rich. The transit depth is considerably affected, increasing values by up to 60 ppm, which is shown to be detectable with JWST and Ariel. The differences are most pronounced for the lighter (i.e. $mu$ $\sim$ 4) atmospheres. Our work illustrates that it is imperative that the field of exoplanet spectroscopy moves toward adapted cross-sections, increasingly optimized for high-mu atmospheres for studies of super-Earths and mini-Neptunes.

The Ariel mission, due to launch in 2029, will obtain spectroscopic information for 1000 exoplanets, providing an unprecedented opportunity for comparative exoplanetology. Retrieval codes - parameteric atmospheric models coupled with an inversion algorithm - represent the tool of choice for interpreting Ariel data. Ensuring that reliable and consistent results can be produced by these tools is a critical preparatory step for the mission. Here, we present the results of a retrieval challenge. We use five different exoplanet retrieval codes to analyse the same synthetic datasets, and test a) the ability of each to recover the correct input solution and b) the consistency of the results. We find that generally there is very good agreement between the five codes, and in the majority of cases the correct solutions are recovered. This demonstrates the reproducibility of retrievals for transit spectra of exoplanets, even when codes are not previously benchmarked against each other., Comment: 28 pages, 14 figures. Accepted in Experimental Astronomy (2022)

This work presents a re-analysis of the spectroscopic phase-curve observations of the ultra hot Jupiter WASP-103 b obtained by the Hubble Space Telescope (HST) and the Spitzer Telescope. Traditional 1D and unified 1.5D spectral retrieval techniques are employed, allowing to map the thermal structure and the abundances of trace gases in this planet as a function of longitude. On the day-side, the atmosphere is found to have a strong thermal inversion, with indications of thermal dissociation traced by continuum H- opacity. Water vapor is found across the entire atmosphere but with depleted abundances of around 1e-5, consistent with the thermal dissociation of this molecule. Regarding metal oxide and hydrides, FeH is detected on the hot-spot and the day-side of WASP-103 b, but TiO and VO are not present in detectable quantities. Carbon-bearing species such as CO and CH4 are also found, but since their detection is reliant on the combination of HST and Spizer, the retrieved abundances should be interpreted with caution. Free and Equilibrium chemistry retrievals are overall consistent, allowing to recover robust constraints on the metallicity and C/O ratio for this planet. The analyzed phase-curve data indicates that the atmosphere of WASP-103 b is consistent with solar elemental ratios., Comment: Accepted in AJ, 21 pages, 15 figures, 2 animations

Thermochemical equilibrium is one of the most commonly used assumptions in current exoplanet retrievals. As the James Webb Space Telescope (JWST) and Ariel launch draw near, assessing the underlying biases and assumptions made when applying self-consistent chemistry into spectral retrievals is crucial. Here we use the flexibility of TauREx 3.1 to cross-compare three state of the art chemical equilibrium codes: ACE, FastChem and GGchem. We simulate JWST spectra for ACE, FastChem, GGchem and GGchem+condensation containing only C, H, O, N elements and spectra for FastChem, GGchem and GGchem+condensation with a more extensive range of elements giving seven simulated JWST spectra in total and then cross-retrieve giving a total of 49 retrievals. Our analysis demonstrates that like-for-like, all chemical codes retrieve the correct parameters to <1% of the truth. However, retrievals, where the contained elements do not match the truth, parameters such as metallicity deviate by 20% while maintaining extremely low uncertainties <1% giving false confidence. This point is of major importance for future analyses on JWST and Ariel, highlighting that self-consistent chemical scheme that do not employ the proper assumptions (missing species, fixed elemental ratios, condensation) are at risk of confidently biasing interpretations. Free chemistry retrievals employing parametric descriptions of the chemical profiles can provide alternative unbiased explorations., Comment: 31 pages, 15 figures

Most of the exoplanets detected up to now transit in front of their host stars, allowing for the generation of transmission spectra; the study of exoplanet atmospheres relies heavily upon accurate analysis of these spectra. Recent discoveries mean that the study of atmospheric signals from low-mass, temperate worlds are becoming increasingly common. The observed transit depth in these planets is small and more difficult to analyze. Analysis of simulated transmission spectra for two small, temperate planets (GJ 1214 b and K2-18 b) is presented, giving evidence for significant differences in simulated transit depth when the water vapor continuum is accounted for when compared to models omitting it. These models use cross-sections from the CAVIAR lab experiment for the water self-continuum up to 10,000 cm$^{-1}$; these cross-sections exhibit an inverse relationship with temperature, hence lower-temperature atmospheres are the most significantly impacted. Including the water continuum strongly affects transit depths, increasing values by up to 80 ppm, with the differences for both planets being detectable with the future space missions Ariel and JWST. It is imperative that models of exoplanet spectra move toward adaptive cross-sections, increasingly optimized for H$_2$O-rich atmospheres. This necessitates including absorption contribution from the water vapor continuum into atmospheric simulations.

The ARIEL Space Telescope will provide a large and diverse sample of exoplanet spectra, performing spectroscopic observations of about 1000 exoplanets in the wavelength range $0.5 \to 7.8 \; \mu m$. In this paper, we investigate the information content of ARIEL's Reconnaissance Survey low resolution transmission spectra. Among the goals of the ARIEL Reconnaissance Survey is also to identify planets without molecular features in their atmosphere. In this work, (1) we present a strategy that will allow to select candidate planets to be reobserved in a ARIEL's higher resolution Tier; (2) we propose a metric to preliminary classify exoplanets by their atmospheric composition without performing an atmospheric retrieval; (3) we introduce the possibility to find other methods to better exploit the data scientific content., Comment: 31 pages, 14 figures, submitted to ApJ

We present the transmission spectrum of the inflated hot-Jupiter WASP-17 b, observed with the STIS and WFC3 instruments aboard the Hubble Space Telescope, allowing for a continuous wavelength coverage from ~0.4 to ~1.7 um. Observations taken with IRAC channel 1 and 2 on the Spitzer Space Telescope are also included, adding photometric measurements at 3.6 and 4.5 um. HST spectral data was analysed with Iraclis, a pipeline specialised in the reduction of STIS and WFC3 transit and eclipse observations. Spitzer photometric observations were reduced with the TLCD-LSTM method, utilising recurrent neural networks. The outcome of our reduction produces incompatible results between STIS visit 1 and visit 2, which leads us to consider two scenarios for G430L. Additionally, by modelling the WFC3 data alone, we can extract atmospheric information without having to deal with the contrasting STIS datasets. We run separate retrievals on the three spectral scenarios with the aid of TauREx 3, a fully Bayesian retrieval framework. We find that, independently of the data considered, the exoplanet atmosphere displays strong water signatures and potentially, the presence of aluminium oxide (AlO) and titanium hydride (TiH). A retrieval that includes an extreme photospheric activity of the host star is the preferred model, but we recognise that such a scenario is unlikely for an F6-type star. Due to the incompleteness of all STIS spectral lightcurves, only further observations with this instrument would allow us to properly constrain the atmospheric limb of WASP-17 b, before JWST or Ariel will come online., Comment: 24 pages, 11 figures, 7 tables, accepted to AJ

Ariel, the Atmospheric Remote-sensing Infrared Exoplanet Large-survey, was adopted as the fourth medium-class mission in ESA's Cosmic Vision programme to be launched in 2029. During its 4-year mission, Ariel will study what exoplanets are made of, how they formed and how they evolve, by surveying a diverse sample of about 1000 extrasolar planets, simultaneously in visible and infrared wavelengths. It is the first mission dedicated to measuring the chemical composition and thermal structures of hundreds of transiting exoplanets, enabling planetary science far beyond the boundaries of the Solar System. The payload consists of an off-axis Cassegrain telescope (primary mirror 1100 mm x 730 mm ellipse) and two separate instruments (FGS and AIRS) covering simultaneously 0.5-7.8 micron spectral range. The satellite is best placed into an L2 orbit to maximise the thermal stability and the field of regard. The payload module is passively cooled via a series of V-Groove radiators; the detectors for the AIRS are the only items that require active cooling via an active Ne JT cooler. The Ariel payload is developed by a consortium of more than 50 institutes from 16 ESA countries, which include the UK, France, Italy, Belgium, Poland, Spain, Austria, Denmark, Ireland, Portugal, Czech Republic, Hungary, the Netherlands, Sweden, Norway, Estonia, and a NASA contribution., Comment: Ariel Definition Study Report, 147 pages. Reviewed by ESA Science Advisory Structure in November 2020. Original document available at: https://www.cosmos.esa.int/documents/1783156/3267291/Ariel_RedBook_Nov2020.pdf/

We present a study on the spatially scanned spectroscopic observations of the transit of GJ 1132 b, a warm ($\sim$500 K) Super-Earth (1.13 R$_\oplus$) that was obtained with the G141 grism (1.125 - 1.650 $\mu$m) of the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope. We used the publicly available Iraclis pipeline to extract the planetary transmission spectra from the five visits and produce a precise transmission spectrum. We analysed the spectrum using the TauREx3 atmospheric retrieval code with which we show that the measurements do not contain molecular signatures in the investigated wavelength range and are best-fit with a flat-line model. Our results suggest that the planet does not have a clear primordial, hydrogen-dominated atmosphere. Instead, GJ 1132 b could have a cloudy hydrogen-dominated envelope, a very enriched secondary atmosphere, be airless, or have a tenuous atmosphere that has not been detected. Due to the narrow wavelength coverage of WFC3, these scenarios cannot be distinguished yet but the James Webb Space Telescope may be capable of detecting atmospheric features, although several observations may be required to provide useful constraints., Comment: 17 pages, 11 figures. Accepted for publication in AJ

The analysis of exoplanetary atmospheres often relies upon the observation of transit or eclipse events. While very powerful, these snapshots provide mainly 1-dimensional information on the planet structure and do not easily allow precise latitude-longitude characterisations. The phase curve technique, which consists of measuring the planet emission throughout its entire orbit, can break this limitation and provide useful 2-dimensional thermal and chemical constraints on the atmosphere. As of today however, computing performances have limited our ability to perform unified retrieval studies on the full set of observed spectra from phase curve observations at the same time. Here, we present a new phase curve model that enables fast, unified retrieval capabilities. We apply our technique to the combined phase curve data from the Hubble and Spitzer space telescopes of the hot-Jupiter WASP-43 b. We tested different scenarios and discussed the dependence of our solution to different assumptions in the model. Our more comprehensive approach suggests that multiple interpretation of this dataset are possible but our more complex model is consistent with the presence of thermal inversions and a metal rich atmosphere, contrasting with previous data analyses, although this likely depends on the Spitzer data reduction. The detailed constraints extracted here demonstrate the importance of developing and understanding advanced phase curve techniques, which we believe will unlock access to a richer picture of exoplanet atmospheres., Comment: Accepted in ApJ, 57 pages, 38 figures

Ariel will mark the dawn of a new era as the first large-scale survey characterising exoplanetary atmospheres with science objectives to address fundamental questions about planetary composition, evolution and formation. In this study, we explore the detectability of atmospheres vaporised from magma oceans on dry, rocky Super-Earths orbiting very close to their host stars. The detection of such atmospheres would provide a definitive piece of evidence for rocky planets but are challenging measurements with currently available instruments due to their small spectral signatures. However, some of the hottest planets are believed to have atmospheres composed of vaporised rock, such as Na and SiO, with spectral signatures bright enough to be detected through eclipse observations with planned space-based telescopes. In this study, we find that rocky super-Earths with a irradiation temperature of 3000 K and a distance from Earth of up to 20 pc, as well as planets hotter than 3500 K and closer than 50 pc, have SiO features which are potentially detectable in eclipse spectra observed with Ariel., Comment: 12 pages, 8 figures, accepted for publication in Experimental Astronomy, Ariel Special Issue

Recent studies of ultra hot-Jupiters suggested that their atmospheres could have thermal inversions due to the presence of optical absorbers such as TiO, VO, FeH and other metal hydride/oxides. However, it is expected that these molecules would thermally dissociate at extremely high temperatures, thus leading to featureless spectra in the infrared. KELT-9b, the hottest exoplanet discovered so far, is thought to belong to this regime and host an atmosphere dominated by neutral Hydrogen from dissociation and atomic/ionic species. Here, we analysed the eclipse spectrum obtained using the Hubble Space Telescope's Wide Field Camera 3 (WFC3) and, by utilising the atmospheric retrieval code TauREx3, found that the spectrum is consistent with the presence of molecular species and is poorly fitted by a simple blackbody. In particular, we find that a combination of TiO, VO, FeH and H- provides the best-fit when considering HST, Spitzer and TESS datasets together. Aware of potential biases when combining instruments, we also analysed the HST spectrum alone and found that TiO and VO only were needed in this case. These findings paint a more complex picture of the atmospheres of ultra-hot planets than previously thought, Comment: 11 pages, 4 figures, Accepted in ApJL on the 2nd of January 2021

The composition of giant planets is imprinted by their migration history and the compositional structure of their hosting disks. Studies in recent literature investigate how the abundances of C and O can constrain the formation pathways of giant planets forming within few tens of au from the star. New ALMA observations, however, suggest planet-forming regions possibly extending to hundreds of au. We explore the implications of these wider formation environments through n-body simulations of growing and migrating giant planets embedded in planetesimal disks, coupled with a compositional model of the protoplanetary disk where volatiles are inherited from the molecular cloud and refractories are calibrated against extrasolar and Solar System data. We find that the C/O ratio provides limited insight on the formation pathways of giant planets that undergo large-scale migration. This limitation can be overcome thanks to nitrogen and sulphur. Jointly using the C/N, N/O and C/O ratios breaks any degeneracy in the formation and migration tracks of giant planets. The use of elemental ratios normalized to the respective stellar ratios supplies additional information on the nature of giant planets, thanks to the relative volatility of O, C and N in disks. When the planetary metallicity is dominated by the accretion of solids C/N* $>$ C/O* $>$ N/O* (* denoting this normalized scale), otherwise N/O* $>$ C/O* $>$ C/N*. The S/N ratio provides an additional independent probe into the metallicity of giant planets and their accretion of solids., Comment: 27 pages, 10 figures and 7 tables. Accepted version of the article published on The Astrophysical Journal (includes 1 extra figure and additional text with respect to previous preprint)

Deep learning algorithms are growing in popularity in the field of exoplanetary science due to their ability to model highly non-linear relations and solve interesting problems in a data-driven manner. Several works have attempted to perform fast retrievals of atmospheric parameters with the use of machine learning algorithms like deep neural networks (DNNs). Yet, despite their high predictive power, DNNs are also infamous for being 'black boxes'. It is their apparent lack of explainability that makes the astrophysics community reluctant to adopt them. What are their predictions based on? How confident should we be in them? When are they wrong and how wrong can they be? In this work, we present a number of general evaluation methodologies that can be applied to any trained model and answer questions like these. In particular, we train three different popular DNN architectures to retrieve atmospheric parameters from exoplanet spectra and show that all three achieve good predictive performance. We then present an extensive analysis of the predictions of DNNs, which can inform us - among other things - of the credibility limits for atmospheric parameters for a given instrument and model. Finally, we perform a perturbation-based sensitivity analysis to identify to which features of the spectrum the outcome of the retrieval is most sensitive. We conclude that for different molecules, the wavelength ranges to which the DNN's predictions are most sensitive, indeed coincide with their characteristic absorption regions. The methodologies presented in this work help to improve the evaluation of DNNs and to grant interpretability to their predictions., Comment: 32 pages, 22 figures. Submitted to AJ. Accepted July 20, 2021

Atmospheric characterisation of temperate, rocky planets is the holy grail of exoplanet studies. These worlds are at the limits of our capabilities with current instrumentation in transmission spectroscopy and challenge our state-of-the-art statistical techniques. Here we present the transmission spectrum of the temperate Super-Earth LHS 1140b using the Hubble Space Telescope (HST). The Wide Field Camera 3 (WFC3) G141 grism data of this habitable zone (T$_{\rm{eq}}$ = 235 K) Super-Earth (R = 1.7 $R_\oplus$), shows tentative evidence of water. However, the signal-to-noise ratio, and thus the significance of the detection, is low and stellar contamination models can cause modulation over the spectral band probed. We attempt to correct for contamination using these models and find that, while many still lead to evidence for water, some could provide reasonable fits to the data without the need for molecular absorption although most of these cause also features in the visible ground-based data which are nonphysical. Future observations with the James Webb Space Telescope (JWST) would be capable of confirming, or refuting, this atmospheric detection., Comment: Accepted for publication in AJ on 30th October 2020

We present an atmospheric characterization study of two medium sized planets bracketing the radius of Neptune: HD 106315 c (R$_{\rm{P}}$=4.98 $\pm$ 0.23 R$_{\oplus}$) and HD 3167 c (R$_{\rm{P}}$=2.740$_{-0.100}^{+0.106}$ R$_{\oplus}$). We analyse spatially scanned spectroscopic observations obtained with the G141 grism (1.125 - 1.650 $\mu$m) of the Wide Field Camera 3 (WFC3) onboard the Hubble Space Telescope. We use the publicly available Iraclis pipeline and TauREx3 atmospheric retrieval code and we detect water vapor in the atmosphere of both planets with an abundance of $\log_{10}[\mathrm{H_2O}]=-2.1^{+0.7}_{-1.3}$ ($\sim$5.68$\sigma$) and $\log_{10}[\mathrm{H_2O}]=-4.1^{+0.9}_{-0.9}$ ($\sim$3.17$\sigma$) for HD 106315 c and HD 3167 c, respectively. The transmission spectrum of HD 106315 c shows also a possible evidence of ammonia absorption ($\log_{10}[\mathrm {NH_3}]=-4.3^{+0.7}_{-2.0}$, $\sim$1.97$\sigma$ -even if it is not significant-), whilst carbon dioxide absorption features may be present in the atmosphere of HD 3167 c in the $\sim$1.1-1.6~$\mu$m wavelength range ($\log_{10}[\mathrm{CO_{2}}]= -2.4^{+0.7}_{-1.0}$, $\sim$3.28$\sigma$). However the CO$_2$ detection appears significant, it must be considered carefully and put into perspective. Indeed, CO$_2$ presence is not explained by 1D equilibrium chemistry models, and it could be due to possible systematics. The additional contribution of clouds, CO and CH$_4$ are discussed. HD 106315 c and HD 3167 c will be interesting targets for upcoming telescopes such as the James Webb Space Telescope (JWST) and the Atmospheric Remote-Sensing Infrared Exoplanet Large-Survey (Ariel)., Comment: Accepted for publication in AJ

In the past decade, the analysis of exoplanet atmospheric spectra has revealed the presence of water vapour in almost all the planets observed, with the exception of a fraction of overcast planets. Indeed, water vapour presents a large absorption signature in the wavelength coverage of the Hubble Space Telescope's (HST) Wide Field Camera 3 (WFC3), which is the main space-based observatory for atmospheric studies of exoplanets, making its detection very robust. However, while carbon-bearing species such as methane, carbon monoxide and carbon dioxide are also predicted from current chemical models, their direct detection and abundance characterisation has remained a challenge. Here we analyse the transmission spectrum of the puffy, clear hot-Jupiter KELT-11 b from the HST WFC3 camera. We find that the spectrum is consistent with the presence of water vapor and an additional absorption at longer wavelengths than 1.5um, which could well be explained by a mix of carbon bearing molecules. CO2, when included is systematically detected. One of the main difficulties to constrain the abundance of those molecules is their weak signatures across the HST WFC3 wavelength coverage, particularly when compared to those of water. Through a comprehensive retrieval analysis, we attempt to explain the main degeneracies present in this dataset and explore some of the recurrent challenges that are occurring in retrieval studies (e.g: the impact of model selection, the use of free vs self-consistent chemistry and the combination of instrument observations). Our results make this planet an exceptional example of chemical laboratory where to test current physical and chemical models of hot-Jupiters' atmospheres., Comment: 24 pages, 14 figures, Accepted in AJ

The study of exoplanetary atmospheres relies on detecting minute changes in the transit depth at different wavelengths. To date, a number of ground and space based instruments have been used to obtain transmission spectra of exoplanets in different spectral band. One common practice is to combine observations from different instruments in order to achieve a broader wavelength coverage. We present here two inconsistent observations of WASP-96 b, one by Hubble Space Telescope (HST) and the other by the Very Large Telescope (VLT). We present two key findings in our investigation: 1.) a strong water signature is detected via the HST WFC3 observations. 2.) A notable offset in transit depth (> 1100 ppm) can be seen when the ground-based and space-based observations are combined together. The discrepancy raises the question of whether observations from different instruments could indeed be combined together. We attempt to align the observations by including an additional parameter in our retrieval studies but are unable to definitively ascertain that the aligned observations are indeed compatible. The case of WASP-96 b signals that compatibility of instruments should not be assumed. While wavelength overlaps between instruments can help, it should be noted that combining datasets remains a risky business. The difficulty in combining observations also strengthens the need for next generation instruments which will possess broader spectral coverage., Comment: 18 pages, 13 figures. Submitted to ApJ, currently under review

We present spectral analysis of the transiting Saturn-mass planet WASP-117b, observed with the G141 grism of Wide Field Camera 3 (WFC3) on the Hubble Space Telescope. We reduce and fit the extracted spectrum from the raw transmission data using the open-source software Iraclis before performing a fully Bayesian retrieval using the publicly available analysis suite TauREx 3.0. We detect water vapour alongside a layer of fully opaque cloud, retrieving a terminator temperature of 833 K. In order to quantify the statistical significance of this detection, we employ the Atmospheric Detectability Index (ADI), deriving a value of 2.30, which provides positive but not strong evidence against the flatline model. Due to the eccentric orbit of WASP-117b, it is likely that chemical and mixing timescales oscillate throughout orbit due to the changing temperature, possibly allowing warmer chemistry to remain visible as the planet begins transit, despite the proximity of its point of ingress to apastron. We present simulated spectra of the planet as would be observed by the future space missions Ariel and JWST and show that, despite not being able to probe such chemistry with current HST data, these observatories should make it possible in the not too distant future., Comment: Accepted for publication in AJ

We report follow-up observations of four transiting exoplanets, TRES-2b, HAT-P-22b, HAT-P-36b and XO-2b, as part of the Original Research By Young Twinkle Students (ORBYTS) programme. These observations were taken using the Las Cumbres Observatory Global Telescope Network's (LCOGT) robotic 0.4 m telescopes and were analysed using the HOlomon Photometric Software (HOPS). Such observations are key for ensuring accurate transit times for upcoming telescopes, such as the James Webb Space Telescope (JWST), Twinkle and Ariel, which may seek to characterise the atmospheres of these planets. The data have been uploaded to ExoClock and a significant portion of this work has been completed by secondary school students in London., Comment: Published in the Research Notes of the American Astronomical Society (RNAAS)

In the recent years, retrieval analysis of exoplanet atmospheres have been very successful, providing deep insights on the composition and the temperature structure of these worlds via the transit and eclipse methods. Analysis of spectral phase curve observations, which in theory provides even more information, are still limited to a few planets. In the next decade, new facilities such as NASA-JWST and ESA-Ariel will revolutionise the field of exoplanet atmospheres and we expect that a significant time will be spent on spectral phase curve observations. Most current models are still limited in their analysis of phase curve data as they do not consider the planet atmosphere as a whole or they require large computational resources. In this paper we present a semi-analytical model that will allow to compute exoplanet emission spectra at different phase angles. Our model provides a way to simulate a large number of observations while being only about 4 times slower than the traditional forward model for plane parallel primary eclipse. This model, which is based on the newly developed TauREx3 framework (Al-Refaie 2020), will be further developed to allow for phase curve atmospheric retrievals., Comment: Accepted in ApJ, 18 pages, 10 figures