We report high-resolution spectroscopic detection of TiO molecular signature in the day-side spectra of WASP-33b, the second hottest known hot Jupiter. We used the High Dispersion Spectrograph (HDS; ...R ∼ 165,000) on the Subaru telescope in the wavelength range of 0.62-0.88 m to obtain the day-side spectra of WASP-33b. We suppress and correct the systematic effects of the instrument and the telluric and stellar lines using the SYSREM algorithm after the selection of good orders based on Barnard's star and other M-type stars. We detect a 4.8 signal at an orbital velocity of km s−1 and systemic velocity of km s−1, which agree with the derived values from a previous analysis of the primary transit. Our detection with the temperature inversion model implies the existence of a stratosphere in its atmosphere; however, we were unable to constrain the volume mixing ratio of the detected TiO. We also measure the stellar radial velocity and use it to obtain a more stringent constraint on the orbital velocity, km s−1. Our results demonstrate that high-dispersion spectroscopy is a powerful tool to characterize the atmosphere of an exoplanet, even in the optical wavelength range, and shows a promising potential in using and developing similar techniques with high-dispersion spectrograph on current 10 m class and future extremely large telescopes.
ABSTRACT
High-resolution Doppler-resolved spectroscopy has presented new opportunities for studying the atmospheres of exoplanets. While the ‘classical’ cross-correlation approach has proven to be ...efficient at finding atmospheric species, it is unable to perform direct atmospheric retrievals. Recent work has shown that retrievals are possible using a direct likelihood evaluation or likelihood ‘mappings’. The unique aspect of high-resolution retrievals is that the data-processing methods required to remove the stellar and telluric lines also distort the underlying exoplanet’s signal and therefore the forward model must be pre-processed to match this filtering. This was the key remaining limitation in our previously published framework. This paper directly addresses this by introducing a simple and fast model-filtering technique that can replicate the processing performed by algorithms such as SysRem and PCA. This enables retrievals to be performed without having to perform expensive injection and pre-processing steps for every model. We show that we can reliably constrain quantitative measures of the atmosphere from transmission spectra including the temperature–pressure profile, relative abundances, planetary velocities, and rotational broadening parameters. Finally, we demonstrate our framework using UVES transmission spectroscopy of WASP-121b. We constrain the temperature–pressure profile and relative abundances of Fe, Cr, and V to be log10(χFe/χCr) = 1.66 ± 0.28, log10(χFe/χV) = 3.78 ± 0.29, and log10(χFe/χMg) = −1.26 ± 0.60. The relative abundances are consistent with solar values, with the exception of Fe/Mg, where the large Mg abundance is probably explained by the escaping atmosphere of WASP-121b that is not accounted for in our atmospheric model.
ABSTRACT
High-resolution Doppler-resolved spectroscopy has opened up a new window into the atmospheres of both transiting and non-transiting exoplanets. Here, we present VLT/UVES observations of a ...transit of WASP-121b, an ‘ultra-hot’ Jupiter previously found to exhibit a temperature inversion and detections of multiple species at optical wavelengths. We present initial results using the blue arm of UVES (≈3700–5000 Å), recovering a clear signal of neutral Fe in the planet’s atmosphere at >8$\, \sigma$, which could contribute to (or even fully explain) the temperature inversion in the stratosphere. However, using standard cross-correlation methods, it is difficult to extract physical parameters such as temperature and abundances. Recent pioneering efforts have sought to develop likelihood ‘mappings’ that can be used to directly fit models to high-resolution data sets. We introduce a new framework that directly computes the likelihood of the model fit to the data, and can be used to explore the posterior distribution of parametrised model atmospheres via MCMC techniques. Our method also recovers the physical extent of the atmosphere, as well as account for time- and wavelength-dependent uncertainties. We measure a temperature of $3710^{+490}_{-510}$ K, indicating a higher temperature in the upper atmosphere when compared to low-resolution observations. We also show that the Fe i signal is physically separated from the exospheric Fe ii. However, the temperature measurements are highly degenerate with aerosol properties; detection of additional species, using more sophisticated atmospheric models, or combining these methods with low-resolution spectra should help break these degeneracies.
Abstract
We report on Fe
i
in the dayside atmosphere of the ultra-hot-Jupiter WASP-33b, providing evidence for a thermal inversion in the presence of an atomic species. We also introduce a new way to ...constrain the planet’s brightness variation throughout its orbit, including its day–night contrast and peak phase offset, using high-resolution Doppler spectroscopy alone. We do so by analyzing high-resolution optical spectra of six arcs of the planet’s phase curve, using Echelle SpectroPolarimetric Device for the Observation of Stars (ESPaDOnS) on the Canada–France–Hawaii telescope and High Dispersion Spectrograph on the Subaru telescope. By employing a likelihood mapping technique, we explore the marginalized distributions of parameterized atmospheric models, and detect Fe
i
emission at high significance (>10.4
σ
) in our combined data sets, located at
K
p
= 222.1 ± 0.4 km s
−1
and
v
sys
= −6.5 ± 0.3 km s
−1
. Our values agree with previous reports. By accounting for WASP-33b’s brightness variation, we find evidence that its nightside flux is <10% of the dayside flux and the emission peak is shifted westward of the substellar point, assuming the spectrum is dominated by Fe
i
. Our ESPaDOnS data, which cover phases before and after the secondary eclipse more evenly, weakly constrain the phase offset to +22 ± 12 degrees. We caution that the derived volume-mixing ratio depends on our choice of temperature-pressure profile, but note it does not significantly influence our constraints on day–night contrast or phase offset. Finally, we use simulations to illustrate how observations with increased phase coverage and higher signal-to-noise ratios can improve these constraints, showcasing the expanding capabilities of high-resolution Doppler spectroscopy.
Abstract
We analyze the high-resolution emission spectrum of WASP-33b taken using the High Dispersion Spectrograph (
R
≈ 165,000) on the 8.2 m Subaru telescope. The data cover
λ
≈ 6170–8817 Å, ...divided over 30 spectral orders. The telluric and stellar lines are removed using a de-trending algorithm, S
ys
R
em
, before cross-correlating with planetary spectral templates. We calculate the templates assuming a 1D plane-parallel hydrostatic atmosphere including continuum opacity of bound–free H
−
and Rayleigh scattering by H
2
with a range of constant abundances of Fe
i
. Using a likelihood-mapping analysis, we detect an Fe
i
emission signature at 6.4
σ
located at
of 226.0
km s
−1
and
v
sys
of −3.2
km s
−1
—consistent with the planet’s expected velocity in the literature. We also confirm the existence of a thermal inversion in the dayside of the planet, which is very likely to be caused by the presence of Fe
i
and previously detected TiO in the atmosphere. This makes WASP-33b one of the prime targets to study the relative contributions of both species to the energy budget of an ultra-hot Jupiter.
ABSTRACT
Ultrahot Jupiters (UHJs) present excellent targets for atmospheric characterization. Their hot dayside temperatures (T ≳ 2200 K) strongly suppress the formation of condensates, leading to ...clear and highly inflated atmospheres extremely conducive to transmission spectroscopy. Recent studies using optical high-resolution spectra have discovered a plethora of neutral and ionized atomic species in UHJs, placing constraints on their atmospheric structure and composition. Our recent work has presented a search for molecular features and detection of Fe i in the UHJ WASP-121b using Very Large Telescope (VLT)/UV–Visual Echelle Spectrograph (UVES) transmission spectroscopy. Here, we present a systematic search for atomic species in its atmosphere using cross-correlation methods. In a single transit, we uncover potential signals of 17 atomic species that we investigate further, categorizing five as strong detections, three as tentative detections, and nine as weak signals worthy of further exploration. We confirm previous detections of Cr i, V i, Ca i, K i, and exospheric H i and Ca ii made with the High Accuracy Radial velocity Planet Searcher (HARPS) and the Echelle SPectrograph for Rocky Exoplanets and Stable Spectroscopic Observations (ESPRESSO), and independently re-recover our previous detection of Fe i at 8.8σ using both the blue and red arms of the UVES data. We also add a novel detection of Sc ii at 4.2σ. Our results further demonstrate the richness of UHJs for optical high-resolution spectroscopy.
ABSTRACT
We analyse the transmission spectra of KELT-20b/MASCARA-2b to search for possible thermal inversion agents. The data consist of three transits obtained using HARPSN and one using CARMENES. ...We removed stellar and telluric lines before cross-correlating the residuals with spectroscopic templates produced using a 1D plane-parallel model, assuming an isothermal atmosphere and chemical equilibrium at solar metallicity. Using a likelihood-mapping method, we detect Fe i at > 13σ, Ca ii H$\&$K at > 6σ and confirm the previous detections of Fe ii, Ca ii IR Triplet, and Na i D. The detected signal of Fe i is shifted by −3.4 ± 0.4 km s−1 from the planetary rest frame, which indicates a strong day–night wind. Our likelihood-mapping technique also reveals that the absorption features of the detected species extend to different altitudes in the planet’s atmosphere. Assuming that the line lists are accurate, we do not detect other potential thermal inversion agents (NaH, MgH, AlO, SH, CaO, VO, FeH, and TiO) suggesting that non-chemical equilibrium mechanisms (e.g. a cold-trap) might have removed Ti- and V-bearing species from the upper atmosphere. Our results, therefore, show that KELT-20b/MASCARA-2b cannot possess an inversion layer caused by a TiO/VO-related mechanism. The presence of an inversion layer would therefore likely be caused by metal atoms such as Fe i and Fe ii. Finally, we report a double-peak structure in the Fe i signal in all of our data sets that could be a signature of atmospheric dynamics. However, further investigation is needed to robustly determine the origin of the signal.
Abstract Transmission spectroscopy presents one of the most successful approaches for investigating the atmospheres of exoplanets. We analyzed the near-infrared high-resolution transmission spectrum ...of a hot Saturn, HD 149026 b, taken using CARMENES spectrograph ( R ∼ 80,400 ). We found evidence of H 2 O at a signal-to-noise ratio (S/N) of ∼4.8. We also performed grid search using a Bayesian framework and constrained the orbital velocity K p and rest velocity V rest to 158.17 − 7.90 + 8.31 km s −1 and 2.57 − 0.57 + 0.54 km s −1 , respectively. While the retrieved K p value is consistent with theoretical prediction, the retrieved V rest value is highly redshifted (>3 σ ). This might be an indication of either anomalous atmospheric dynamics at play or an orbit with nonzero eccentricity. Additionally, we searched for HCN but no successful detection has been made, possibly due to the relatively low-S/N data set. The detection of H 2 O and subsequent retrieval of its abundance, coupled with analysis of other species such as CO in the K band, for example, might help us to get some information about the atmospheric C/O ratio and metallicity, which in turn could give us some insight into the planet’s formation scenario.
Abstract
We report the first detection of a hydroxyl radical (OH) emission signature in the planetary atmosphere outside the solar system, in this case, in the dayside of WASP-33b. We analyze ...high-resolution near-infrared emission spectra of WASP-33b taken using the InfraRed Doppler spectrograph on the 8.2 m Subaru telescope. The telluric and stellar lines are removed using a detrending algorithm,
SysRem
. The residuals are then cross-correlated with OH and H
2
O planetary spectrum templates produced using several different line lists. We check and confirm the accuracy of OH line lists by cross-correlating with the spectrum of GJ 436. As a result, we detect the emission signature of OH at
K
p
of
km s
−1
and
v
sys
of −0.3
km s
−1
with a signal-to-noise ratio (S/N) of 5.4 and a significance of 5.5
σ
. Additionally, we marginally detect H
2
O emission in the
H
-band with an S/N of 4.0 and a significance of 5.2
σ
using the POKAZATEL line list. However, no significant signal is detected using the HITEMP 2010, which might be due to differences in line positions and strengths, as well as the incompleteness of the line lists. Nonetheless, this marginal detection is consistent with the prediction that H
2
O is mostly thermally dissociated in the upper atmosphere of the ultra-hot Jupiters. Therefore, along with CO, OH is expected to be one of the most abundant O-bearing molecules in the dayside atmosphere of ultra-hot Jupiters and should be considered when studying their atmospheres.
Context.
Efficient absorption of stellar ultraviolet and visible radiation by TiO and VO is predicted to drive temperature inversions in the upper atmospheres of hot Jupiters. However, very few ...inversions or detections of TiO or VO have been reported, and results are often contradictory.
Aims.
Using the improved ExoMol T
OTO
line list, we searched for TiO emission in the dayside spectrum of WASP-33b using the same data in which the molecule was previously detected with an older line list at 4.8
σ
. We intended to confirm the molecular detection and quantify the signal improvement offered by the ExoMol T
OTO
line list.
Methods.
Data from the High Dispersion Spectrograph on the Subaru Telescope was extracted and reduced in an identical manner to the previous study. Stellar and telluric contamination were then removed. High-resolution TiO emission models of WASP-33b were created that spanned a range of molecular abundances using the radiative transfer code
petitRADTRANS
, and were subsequently cross-correlated with the data.
Results.
We measure a 4.3
σ
TiO emission signature using the ExoMol T
OTO
models, corresponding to a WASP-33b orbital velocity semi-amplitude of
K
p
=252.9
−5.3
+5.0
km s
-1
and a system velocity of
v
sys
=−23.0
−4.6
+4.7
km s
-1
. Injection-recovery tests using models based on the new and earlier line lists indicate that if the new models provide a perfect match to the planet spectrum, the significance of the TiO detection should have increased by a factor of ~2.
Conclusions.
Although the TiO signal we find is statistically significant, comparison with previous works makes our result too ambiguous to claim a clear-cut detection. Unexpectedly, the new ExoMol T
OTO
models provide a weaker signal than that found previously, which is offset in
K
p
-
v
sys
space. This sheds some doubt on both detections, especially in light of a recently published TiO non-detection using a different dataset.