Issue |
A&A
Volume 656, December 2021
|
|
---|---|---|
Article Number | A119 | |
Number of page(s) | 8 | |
Section | Planets and planetary systems | |
DOI | https://doi.org/10.1051/0004-6361/202141696 | |
Published online | 08 December 2021 |
Detection of OH in the ultra-hot Jupiter WASP-76b
1
Leiden Observatory, Leiden University,
Postbus 9513,
2300 RA
Leiden,
The Netherlands
e-mail: rlandman@strw.leidenuniv.nl
2
Max-Planck-Institut für Astronomie,
Königstuhl 17,
69117
Heidelberg,
Germany
Received:
2
July
2021
Accepted:
15
October
2021
Context. Ultra-hot Jupiters have dayside temperatures at which most molecules are expected to thermally dissociate. The dissociation of water vapour results in the production of the hydroxyl radical (OH). While OH absorption is easily observed in near-infrared spectra of M dwarfs, which have similar effective temperatures as ultra-hot Jupiters, it is often not considered when studying the atmospheres of ultra-hot Jupiters. Ground-based high-resolution spectroscopy during the primary transit is a powerful tool for detecting molecular absorption in these planets.
Aims. We aim to assess the presence and detectability of OH in the atmosphere of the ultra-hot Jupiter WASP-76b.
Methods. We use high-resolution spectroscopic observations of a transit of WASP-76b obtained using CARMENES. After validating the OH line list, we generate model transit spectra of WASP-76b with petitRADTRANS. The data are corrected for stellar and telluric contamination and cross-correlated with the model spectra. After combining all cross-correlation functions from the transit, a detection map is constructed. Constraints on the planet properties from the OH absorption are obtained from a Markov chain Monte Carlo analysis.
Results. OH is detected in the atmosphere of WASP-76b with a peak signal-to-noise ratio of 6.1. From the retrieval we obtain Kp = 232 ± 12 km s−1 and a blueshift of − 13.2 ± 1.6 km s−1, which are offset from the expected velocities. Considering the fast spin rotation of the planet, the blueshift is best explained with the signal predominantly originating from the evening terminator and the presence of a strong dayside-to-nightside wind. The increased Kp over its expected value (196.5 km s−1) is, however, a bit puzzling. The signal is found to be broad, with a full width at half maximum of 16.8−4.0+4.6 km s−1. The retrieval results in a weak constraint on the mean temperature of 2700–3700 K at the pressure range of the OH signal.
Conclusions. We show that OH is readily observable in the transit spectra of ultra-hot Jupiters. Studying this molecule can provide insights into the molecular dissociation processes in the atmospheres of such planets.
Key words: planets and satellites: atmospheres / planets and satellites: individual: WASP-76b / techniques: spectroscopic / methods: observational
© ESO 2021
1 Introduction
Ultra-hot Jupiters (UHJs) are gas giant planets that orbit extremely close to their host star; as such, they receive large amounts of radiation, resulting in dayside temperatures of >2200 K (Parmentier et al. 2018). They provide a unique opportunity to study atmospheric chemistry under extreme conditions that are not present in any of the planets in the Solar System. Furthermore, because of their short orbital periods and bloated atmospheres, they are well suited for atmospheric characterisation through transmission spectroscopy. At these high dayside temperatures, most molecules are thermally dissociated into their atomic constituents. This has led to a large range of detections of atoms and ions in the atmospheres of UHJs (e.g. Hoeijmakers et al. 2019, 2020; Casasayas-Barris et al. 2019; Yan et al. 2019; Tabernero et al. 2021). However, some diatomic and triatomic molecules are still expected to be present at temperatures typical for UHJs. For example, CO, H2O, OH, FeH, TiO, and VO absorption can be observed in the spectra of M dwarfs (Reiners et al. 2018; Rajpurohit et al. 2018), which have similar effective temperatures as the daysides of UHJs, albeit at higher surface gravities. Furthermore, recombination of these molecules is expected to occur on the limbs or nightside of the planet (Parmentier et al. 2018). However, high spectral resolution searches for these molecules in UHJs have mostly been unsuccessful (e.g. Kesseli et al. 2020; Merritt et al. 2020). A detection of TiO was claimed in the atmosphere of WASP-33b (Nugroho et al. 2017), but not reproduced by Herman et al. (2020) or Serindag et al. (2021). Recently, new evidence for TiO in the dayside spectrum of WASP-33b has been presented (Cont et al. 2021).
Space-based low spectral resolution observations of UHJs with the Wide Field Camera 3 (WFC3) of the Hubble Space Telescope (HST) also show mostly featureless and blackbody-like emission spectra (e.g. Mansfield et al. 2018; Kreidberg et al. 2018; Arcangeli et al. 2018). This has been suggested as being the result of the dissociation of water on the dayside of these planets, in combination with the consequent increase in H− opacity (Lothringer et al. 2018; Parmentier et al. 2018; Arcangeli et al. 2018; Gandhi et al. 2020). Indeed, Parmentier et al. (2018) show that for a subset of UHJs, which have sufficiently high dayside temperatures and sufficiently low surface gravities, it is expected that water is partially or completely thermally dissociated. In the case of WASP-76b, their models show that about half of the water vapour is thermally dissociated at the expected 1.4 μm photosphere (Parmentier et al. 2018, Fig. 13).
One of the products of the dissociation of water is the hydroxyl radical (OH). Still, OH is often not considered when studying atmospheresof UHJs. There have been some attempts at modelling the expected OH abundance in hot Jupiters and rocky exoplanets (Moses et al. 2011; Hu et al. 2012; Miguel & Kaltenegger 2014). These have mostly focused on planets with equilibrium temperatures of < 2000 K. For these planets, the main OH production mechanism is the photolysis of water, which enhances its abundance in the upper atmosphere over the equilibrium case. Indeed, OH has been detected in the atmospheres of the Earth (Meinel 1950), Venus (Piccioni et al. 2008), and Mars (Clancy et al. 2013). However, in the case of UHJs, thermal dissociation of water is expected to be the main OH production mechanism, increasing its abundance and making it potentially easier to detect in these planets. Nugroho et al. (2021) presented the first evidence for OH in the dayside emission spectrum of WASP-33b at a significance of 5.5σ, showing that this is indeed the case. In contrast, water was only marginally detected in WASP-33b, indicating that it is largely thermally dissociated in the upper atmosphere.
WASP-76b is an UHJ orbiting an F7V star with an orbital period of 1.8 days and an equilibrium temperature of ~ 2160 ±40 K (West et al. 2016). Observations of the transmission and emission spectrum with HST/WFC3 revealed the presence of water in the atmosphere of WASP-76b (Tsiaras et al. 2018; Edwards et al. 2020). Furthermore, the WFC3 observations showed hints of TiO and a thermal inversion in its atmosphere (Edwards et al. 2020). High-resolution transmission spectroscopy from the ground using HARPS showed the presence of sodium (Seidel et al. 2019; Žák et al. 2019), which was later confirmed using the HST’s Space Telescope Imaging Spectrograph (von Essen et al. 2020). Furthermore, Na I, Fe I, Li I, Mg I, Ca II, and Mn I were detected using ESPRESSO by Tabernero et al. (2021). However, no molecular absorption from species such as TiO or VO was observed in the optical. Ehrenreich et al. (2020) showed that the iron absorption feature in the ESPRESSO data is asymmetric during the transit, which was later confirmed with HARPS-N (Kesseli & Snellen 2021). The iron absorption was found to increase in strength along the transit and to progressively blueshift from between 0 and −5 km s−1 during ingress up to − 11 km s−1 at the end of the transit. The tidally locked rotation of the planet, with a rotation speed of ~ 5.3 km s−1, results in the absorption on the trailing limb being more blueshifted than the absorption of the leading limb. Furthermore, winds blowing from the hot dayside to the cooler nightside can result in an additional blueshift of several km s−1. The blueshift of − 11 ±0.7 km s−1 of the iron signal during the second half of the transit can thus be explained by the combination of the signal mostly originating from the trailing evening-side limb and the presence of a strong dayside-to-nightside wind. The absence of the signal on the morning side was interpreted as iron condensing across the nightside. This was the first evidence for a chemical gradient in the atmosphere of an exoplanet. Three-dimensional Monte Carlo radiative transfer simulations by Wardenier et al. (2021) show that condensation of iron on the nightside could indeed explain the observations. However, they can also be explained by a large temperature contrast between the leading and trailing limbs of the planet.
In this paper we report on the detection of OH in the atmosphere of WASP-76b using high-resolution transmission spectroscopy with CARMENES (Calar Alto high-Resolution search for M dwarfs with Exoearths with Near-infrared and optical Echelle Spectrographs). This paper is structured as follows: Sect. 2 reports on the observations and data processing. Section 3 describes the modelling of the OH transmission spectrum of WASP-76b. Section 4 presents the results, which are then discussed in Sect. 5.
Parameters of WASP 76b and its host star.
2 Observations and data reduction
2.1 Observations
We used archival data obtained with CARMENES, a high-resolution spectrograph at the 3.5 meter telescope at the Calar Alto Observatory. It covers both the near-infrared (NIR) and the visible with spectral resolutions of R ~ 80 400 and R ~ 94 600, respectively (Quirrenbach et al. 2014, 2018). We downloaded the publicly available NIR observations of the 4 October 2018 transit of WASP-76b from the CAHA Archive1. The retrieved CARMENES spectra had already been reduced using the CARACAL pipeline (Zechmeister et al. 2014; Caballero et al. 2016). The observations consist of 44 exposures of 498 seconds each, with an average signal-to-noise ratio (S/N) per pixel of ~56 in the continuum. We excluded the last eight and the first out-of-transit spectra from our analysis because of a lower S/N. The remaining 35 spectra cover orbital phases between −0.068 and +0.053. The relevant parameters of WASP-76b and its host star are given in Table 1. Throughout this paper, we only use spectral orders of the NIR channel that are not dominated by telluric absorption, as shown in Fig. 2.
2.2 Stellar and telluric corrections
We followed a similar cleaning procedure as previous work with CARMENES data (e.g. Alonso-Floriano et al. 2019; Sánchez-López et al. 2019; Kesseli et al. 2020). All these steps were performed separately for each of the used orders. First, we removed the continuum by fitting a third-order polynomial and dividing by this fit. Subsequently, we flagged wavelength bins that have on average less than 40% flux or 10% excess flux with respect to the continuum. This removes those parts of the spectra that are dominated by telluric absorption or sky emission lines. Finally, we flagged any remaining 5σ outliers. The flagged data points are not considered in the remainder of the analysis.
The spectra are still dominated by telluric and stellar features at this point. We used the SysRem algorithm (Tamuz et al. 2005)to correct for this. SysRem uses an iterative, noise-weighted principal component analysis to remove any shared systematics between light curves. SysRem has been successfully used for removing tellurics in high-resolution spectroscopy data in previous work (e.g. Birkby et al. 2017; Nugroho et al. 2017; Alonso-Floriano et al. 2019; Sánchez-López et al. 2019). Such a blind detrending algorithm can be applied because the radial velocity (RV) of the planet changes by tens of km s−1 during the transit, while stellar and telluric features virtually remain at the same wavelength location. However, when too many SysRem iterations are applied, this leads to the subtraction of the planet signal itself. On the other hand, removing too few modes leads to residual tellurics, which increases the noise and can lead to spurious signals. We used nine SysRem iterations for each of the spectral orders and models in this paper. Other studies have optimised the number of SysRem iterations for each of the spectral orders separately (e.g. Sánchez-López et al. 2019, 2020). This optimisation requires an accurate model of the strength of the planet signal to estimate when self-subtraction of the planet signal occurs, which is not trivial. Therefore, we chose a conservative approach and used the same number of SysRem iterations for each of the spectral orders and models. Appendix B shows the impact of changing the number of SysRem iterations on the S/N of both the observed signal and an injected signal. The pixel uncertainties used for SysRem are obtained from the result of the CARACAL pipeline and are propagated through the preprocessing steps.
![]() |
Fig. 1 Volume mixing ratio of H2O and OH using the chemical equilibrium model from Mollière et al. (2017) and an isothermal P-T profile for different temperatures. |
2.3 Cross-correlation
After removing the stellar and telluric features, individual spectral lines from the planet are still hidden in the noise. The signal from all of the tens to hundreds of individual spectral lines can be combined by cross-correlating with a representative template spectrum of the targeted species. The cross-correlation technique for detecting molecular absorption or emission from exoplanet atmospheres has been extensively used (e.g. Snellen et al. 2010, 2014; Brogi et al. 2012; Giacobbe et al. 2021). This cross-correlation is performed by linearly interpolating the template spectrum to the Doppler-shifted wavelengths of the data and taking the inner product. We did this for all orders simultaneously with an RV range from −300 to +300 km s−1 with a step size of 1.3 km s−1. This was done for each of the spectra, which results in a 2D cross-correlation map. Since the radial component of the orbital velocity of the planet changes as a function of the orbital phase, the planet signal would appear as a slanted line in this cross-correlation map. The RV of the planet is given by:
(1)
Here, vbar is the barycentric velocity of the Earth at the time of the exposure, vsys the systemic velocity of the system, Kp the RV semi-amplitude of the planet, and ϕ the orbital phase of the planet. The cross-correlation functions (CCFs) were then shifted to the planet’s rest frame for each exposure. To further enhance the signal, we then summed the CCFs from all in-transit exposures, obtaining a total CCF for the entire transit. In contrast to some previous studies, we did not employ a weighting scheme for different spectral orders or exposures. The S/N of the cross-correlation peak was then estimated by subtracting the mean value and dividing by the standard deviation. The mean and standard deviation were calculated excluding the central region of ± 50 km s−1.
3 Models
The cross-correlation requires a model of the OH transmission spectrum of the planet. These model spectra are generated using the high-resolution mode of petitRADTRANS (Mollière et al. 2019), a radiative transfer code developed specifically for the characterisation of exoplanet atmospheres. As input it requires pre-computed opacity grids of the desired molecules, the abundance profiles of the molecules, the pressure-temperature (P-T) profile of the atmosphere, and the surface gravity and radius of the planet. The OH line list was retrieved from the HITEMP database (Rothman et al. 2010). This line list was validated on the spectrum of the M dwarf GX And, and was found to be sufficiently accurate for this study (see Appendix A).
We also included continuum opacity due to H2 −H2, H2 −He scattering and H− in our models. Since high-resolution transmission spectroscopy is not very sensitive to absolute abundances and changes in the P-T profile, we simply used an isothermal atmosphere. The OH abundance and the mean molecular weight were determined using the chemical equilibrium model from Mollière et al. (2017), assuming solar abundance ratios and metallicity. The OH abundances for isothermal atmospheres of different temperatures are shown in Fig. 1. It shows that for higher temperatures, OH is expected to be located lower in the atmosphere with a higher volume mixing ratio. This is because at a higher pressure, water vapour thermally dissociates at a higher temperature, as shown in Fig. 1.
The output of petitRADTRANS consists of the effective transit radius of the planet at each sampled wavelength. The transit depth for an isothermal atmosphere at the equilibrium temperature of 2160 K is shown in Fig. 2. This shows that there is some absorption around 1 μm, but most of the OH opacity is located between 1.4 and 1.7 μm. This is in agreement with our findings from the line validation using an M dwarf (Appendix A). We included all spectral orders that are not heavily contaminated by tellurics, as indicated in Fig. 2.
The resulting models were subsequently convolved with a Gaussian kernel that matched the spectral resolution of CARMENES. Furthermore, the RV change of the planet during a single exposure of 498 s is about 3.9 km s−1. This results in a further smearing of the signal. To take this into account, we convolved the model with a box car with this width. Finally, we removed the continuum of the model by fitting a third-order polynomial and divided by this.
![]() |
Fig. 2 Top: telluric transmission as a function of wavelength obtained from ESO SkyCalc (Noll et al. 2012). Middle: transit depth as a function of wavelength, modelled using petitRADTRANS over the entire CARMENES NIR channel and assuming an isothermal atmosphere at the planet equilibrium temperature of 2160 K. The red region indicates the spectral range used for our analysis. Bottom: same but for T = 3200 K. |
4 Results
The cross-correlation map we obtained using the model that assumes an isothermal atmosphere with a temperature of 2160 K is shown in Fig. 3. The slanted white lines indicate where the planet signal is expected to be. This cross-correlation map was converted to a detection significance for each Kp -vrest combinationby co-adding the CCFs in the rest frame of the planet following Eq. (1). The resulting signal-to-noise map is shown in Fig. 4. We detect OH with a peak S/N of 6.1. Also shown are the 1D CCFs for both Kp = 196.5 km s−1, the expected value, and 232.2 km s−1, which is the value we obtain from the retrieval.
We performed a Markov chain Monte Carlo (MCMC) retrieval to constrain the atmospheric conditions in the terminator region of WASP 76b, following the log-likelihood mapping from Gibson et al. (2020):
(2)
Here, N is the number of wavelength bins, α the model scaling parameter, fi the data, mi the model, and σi the error at wavelength bin i. As free parameters we used Kp, the RV of the OH signal in the planet rest frame, vrest, the temperature in the planet’s atmosphere, the full width at half maximum (FWHM) of the Gaussian broadening kernel applied to the transmission spectrum obtained from petitRADTRANS, and the scaling parameter α. We used uniform priors with bounds specified in Table 2 and an isothermal P-T profile, and calculated the OH abundance following the chemical equilibrium model from Mollière et al. (2017). To decrease the computation time, we pre-computed the models on a grid with a temperature resolution of 50 K and interpolated linearly between the two closest temperatures. The reason that we included the FWHM of the signal is because this can point to dynamics in the planet’s atmosphere. Broadening of the signal can, for example, be caused by the rotation of the planet (Snellen et al. 2014; Brogi et al. 2016), equatorial super-rotation, or high-altitude vertical winds (Seidel et al. 2020). Indeed, absorption features from WASP-76b have been observed to be significantly broadened (Seidel et al. 2019; Tabernero et al. 2021).
The retrieval was done using the emcee package for python (Foreman-Mackey et al. 2013). We used 50 walkers, ran the MCMC for 3000 iterations, and discarded the first 1000 iterations as the burn-in. The resulting corner plot of our retrieval is shown in Fig. 5.
![]() |
Fig. 3 Resulting CCF for each of the spectra, with the RV in the Earth’s rest frame on the horizontal axis and the planet orbital phase on the vertical axis. The horizontal dashed white lines indicate the start and end of the transit, while the slanted vertical lines show the expected slope and location of the planetary signal. |
![]() |
Fig. 4 Bottom: detection significance as a function of RV semi-amplitude, Kp, and velocity relative to the planet rest frame, vrest, showing a ~6σ peak. Top: detection significance at the expected Kp of 196.5 km s−1 and at Kp = 232.2 km s−1, which is obtained from the retrieval. |
![]() |
Fig. 5 Resulting posterior distributions for the OH signal in WASP-76b obtained from the MCMC retrieval. The dashed vertical lines indicate the 16th, 50th, and 84th percentiles. |
Parameters included in the retrieval and their priors.
5 Discussion
From the retrieval we obtain Kp = 232.2 ± 12 km s−1; this is offset from the expected planet orbital velocity of 196.5 km s−1, which was calculated from the stellar mass, inclination, and orbital period (Ehrenreich et al. 2020). From the posterior distribution we find that the expected Kp is excluded at a significance of 2.7σ. Furthermore, we find that the signal has a very significant blueshift of − 13.2 ± 1.6 km s−1 with respect to the planet rest frame.
Wardenier et al. (2021) show that rotation and atmospheric dynamics can lead to peaks away from the expected Kp and vrest. The signal from the evening limb is blueshifted due to a combination of planet rotation and a strong dayside-to-nightside wind, leading to a total blueshift of − 11 ± 0.7 km s−1 (Ehrenreich et al. 2020). On the other hand, on the morning limb, the planet rotation and dayside-to-nightside wind are in opposite directions, resulting in a signal between 0 and −5 km s−1. Our retrieved vrest is similar to that found for iron at the end of the transit (Ehrenreich et al. 2020). We therefore argue that the OH signal presented here is dominated by the evening terminator. This may be the result of a chemical asymmetry as the atmosphere at the evening terminator consists of gas emerging from the hotter dayside, where the OH abundance is expected to be higher (see Fig. 1). In addition, photochemical processes may have had the chance to further enrich the atmosphere in OH (Moses et al. 2011; Miguel & Kaltenegger 2014). Alternatively, this may be caused by the larger scale height on the evening terminator or clouds forming on the morning terminator, which both lead to the evening terminator dominating the transmission signal (Wardenier et al. 2021; Savel et al. 2021). However, we cannot rule out contributions from the morning terminator at velocities between 0 and −5 km s−1, especially given the broadness of the signal.
The offset in Kp, if real, is a bit puzzling. The simulations from Wardenier et al. (2021) indicate that 3D effects would mostly lead to lower Kp values than expected. This means that effects not included in their model need to be invoked to explain the OH signal obtained here. The lower Kp values simulated by Wardenier et al. (2021) result from an interplay between planet rotation, variations in viewing angle, and a particular iron-distribution and/or thermal structure in the planet atmosphere. This results in the morning-side and evening-side terminators dominating the transmission spectrum in the first and second half of the transit, respectively, leading to a decrease in the observed Doppler shift of the signal and thus a lower Kp. The marginally higher-than-expected Kp from the planet orbit observed here would imply that the blueshift of the signal decreases along the transit. While the origin of this is unclear, it is interesting that Sánchez-López et al. (2021) find similar increased Kp values for H2O and HCN. Furthermore, Kesseli et al. (2021) also find evidence for increased Kp values for specific species (Li, Na, and H). Higher S/N observations will be needed to validate the offset in Kp observed here and to study the time-resolved velocity shift of the OH signal in more detail. Furthermore, 3D global circulation models of the OH distribution and the subsequent transmission spectrum will be needed to interpret the signals, but this is beyond the scope of this work.
We find a weak constraint on the mean atmospheric temperature of 2700–3700 K (68% confidence interval). While the error margins are large, this is much higher than the equilibrium temperature of 2160 ± 40 K (West et al. 2016) and is similar to the value obtained in Seidel et al. (2021) for Na absorption. This supports the evidence that the signal originates from the hotter parts of the planet. However, a shared retrieval with both low- and high-resolution spectra may provide better constraints on the thermal structure in the atmosphere (Brogi et al. 2017).
We find a value of km s−1 for the FWHM of the signal. This means the absorption features are significantly broader than expected from instrumental and observational effects only. A broadened absorption signal in WASP-76b was also found for iron (Ehrenreich et al. 2020), sodium (Seidel et al. 2019), and many other species (Tabernero et al. 2021). Seidel et al. (2021) show that this broadening may be explained by including high velocity vertical winds in the upper atmosphere. The broad OH signal observed here may again be evidence for strong dynamics in the upper atmosphere of WASP-76b. The scaling parameter
indicates that the strength of the observed signal is close to the expected strength from the models.
An important caveat to the retrieval is that we used models that assume a simple 1D isothermal atmosphere, but including 3D effects isnecessary for interpreting the transmission spectra of UHJs (e.g. Caldas et al. 2019; MacDonald et al. 2020; Wardenier et al. 2021). Furthermore, the procedure of removing the stellar and telluric contamination using SysRem may alter the planetary signal and bias the obtained posteriors (Brogi & Line 2019). It may especially lead to a decrease in the signal strength, decreasing the scaling parameter, and a smearing of the signal, increasing the observed FWHM.
Finally, our detection shows that OH is one of most abundant spectroscopically active molecules in the atmosphere of WASP-76b and that it should be considered when studying similar UHJs. Our detection reinforces the conclusions made in Nugroho et al. (2021), who detect OH in the dayside emission spectrum of WASP-33b. WASP-33b has a dayside temperature even higher than that of WASP-76b, meaning that water is expected to thermally dissociate at higher pressures as compared to WASP-76b. Since emission spectroscopy generally probes higher pressures than transmission spectroscopy, this is consistent with both detections. OH may therefore be used as a probe for the thermal dissociation processes in these UHJs.
Acknowledgements
R.L., A.S.L, A.K, and I.S. acknowledge funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program under grant agreement No 694513. P.M. acknowledges support from the European Research Council under the European Union’s Horizon 2020 research and innovation program under grant agreement No. 832428. This research has made use of the Spanish Virtual Observatory (http://svo.cab.inta-csic.es) supported by the MINECO/FEDER through grant AyA2017-84089.7. We thank the anonymous referee for the valuable comments that have improved the quality of this work.
Appendix A Validation of the OH line list
We validated the HITEMP OH line list (Rothman et al. 2010) on the spectra of an M dwarf. We used CARMENES observations of GX And that were published as part of a survey searching for exoplanets around M dwarfs (Reiners et al. 2018). GX And is an M1 star with a J-band magnitude of 5.25. We used the observation taken on 10 November 2016, the OH lines of which are already denoted in the appendix of Reiners et al. (2018). The cross-correlation template was generated using petitRADTRANS (Mollière et al. 2019), as described in Sect. 3. Figure A.1 shows the resulting cross-correlation S/N at the star’s systemic velocity for each of the spectral orders. We see that the highest S/N is obtained in the reddest orders of the CARMENES data.
![]() |
Fig. A.1 Cross-correlation S/N with the spectrum of GX And with the OH model template for each of the spectral orders from the CARMENES NIR channel. |
Appendix B Impact of SysRem iterations on the signal-to-noise ratio
Figure B.1 shows the effect of changing the number of SysRem iterations on the peak S/N of the observed OH signal. Also shown is the retrieved S/N for an injected signal at the expected strength with an offset of +100 km/s. While the maximum S/N for the injected signal is obtained at five SysRem iterations, we observed some clear residuals in the cross-correlation map due to telluric contamination. We therefore chose to use nine SysRem iterations, which gives the highest S/N of the observed signal and gives a much cleaner CCF map.
![]() |
Fig. B.1 Obtained peak S/N as a function of the number of SysRem iterations applied to the data. |
References
- Alonso-Floriano, F. J., Sánchez-López, A., Snellen, I. A. G., et al. 2019, A&A, 621, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Arcangeli, J., Désert, J.-M., Line, M. R., et al. 2018, ApJ, 855, L30 [Google Scholar]
- Birkby, J. L., de Kok, R. J., Brogi, M., Schwarz, H., & Snellen, I. A. G. 2017, AJ, 153, 138 [NASA ADS] [CrossRef] [Google Scholar]
- Brogi, M., & Line, M. R. 2019, AJ, 157, 114 [Google Scholar]
- Brogi, M., Snellen, I. A. G., de Kok, R. J., et al. 2012, Nature, 486, 502 [Google Scholar]
- Brogi, M., de Kok, R. J., Albrecht, S., et al. 2016, ApJ, 817, 106 [Google Scholar]
- Brogi, M., Line, M., Bean, J., Désert, J. M., & Schwarz, H. 2017, ApJ, 839, L2 [Google Scholar]
- Caballero, J. A., Guàrdia, J., López del Fresno, M., et al. 2016, SPIE Conf. Ser., 9910, 99100E [Google Scholar]
- Caldas, A., Leconte, J., Selsis, F., et al. 2019, A&A, 623, A161 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Casasayas-Barris, N., Pallé, E., Yan, F., et al. 2019, A&A, 628, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Clancy, R., Sandor, B. J., García-Muñoz, A., et al. 2013, Icarus, 226, 272 [CrossRef] [Google Scholar]
- Cont, D., Yan, F., Reiners, A., et al. 2021, A&A, 651, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Edwards, B., Changeat, Q., Baeyens, R., et al. 2020, AJ, 160, 8 [Google Scholar]
- Ehrenreich, D., Lovis, C., Allart, R., et al. 2020, Nature, 580, 597 [Google Scholar]
- Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [Google Scholar]
- Gandhi, S., Madhusudhan, N., & Mandell, A. 2020, AJ, 159, 232 [Google Scholar]
- Giacobbe, P., Brogi, M., Gandhi, S., et al. 2021, Nature, 592, 205 [NASA ADS] [CrossRef] [Google Scholar]
- Gibson, N. P., Merritt, S., Nugroho, S. K., et al. 2020, MNRAS, 493, 2215 [Google Scholar]
- Herman, M. K., de Mooij, E. J. W., Jayawardhana, R., & Brogi, M. 2020, AJ, 160, 93 [Google Scholar]
- Hoeijmakers, H. J., Ehrenreich, D., Kitzmann, D., et al. 2019, A&A, 627, A165 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hoeijmakers, H. J., Seidel, J. V., Pino, L., et al. 2020, A&A, 641, A123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hu, R., Seager, S., & Bains, W. 2012, ApJ, 761, 166 [CrossRef] [Google Scholar]
- Kesseli, A. Y., & Snellen, I. A. G. 2021, ApJ, 908, L17 [Google Scholar]
- Kesseli, A. Y., Snellen, I. A. G., Alonso-Floriano, F. J., Mollière, P., & Serindag, D. B. 2020, AJ, 160, 228 [NASA ADS] [CrossRef] [Google Scholar]
- Kesseli, A. Y., Snellen, I., Casasayas-Barris, N., Mollière, P., & Sanchez-Lopez, A. 2021, AJ, accepted [arXiv: arXiv:2111.09916] [Google Scholar]
- Kreidberg, L., Line, M. R., Parmentier, V., et al. 2018, AJ, 156, 17 [NASA ADS] [CrossRef] [Google Scholar]
- Lothringer, J. D., Barman, T., & Koskinen, T. 2018, ApJ, 866, 27 [NASA ADS] [CrossRef] [Google Scholar]
- MacDonald, R. J., Goyal, J. M., & Lewis, N. K. 2020, ApJ, 893, L43 [NASA ADS] [CrossRef] [Google Scholar]
- Mansfield, M., Bean, J. L., Line, M. R., et al. 2018, AJ, 156, 10 [NASA ADS] [CrossRef] [Google Scholar]
- Meinel, I. A. B. 1950, ApJ, 111, 555 [NASA ADS] [CrossRef] [Google Scholar]
- Merritt, S. R., Gibson, N. P., Nugroho, S. K., et al. 2020, A&A, 636, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Miguel, Y., & Kaltenegger, L. 2014, ApJ, 780, 166 [Google Scholar]
- Mollière, P., van Boekel, R., Bouwman, J., et al. 2017, A&A, 600, A10 [Google Scholar]
- Mollière, P., Wardenier, J. P., van Boekel, R., et al. 2019, A&A, 627, A67 [Google Scholar]
- Moses, J. I., Visscher, C., Fortney, J. J., et al. 2011, ApJ, 737, 15 [Google Scholar]
- Noll, S., Kausch, W., Barden, M., et al. 2012, A&A, 543, A92 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Nugroho, S. K., Kawahara, H., Masuda, K., et al. 2017, AJ, 154, 221 [Google Scholar]
- Nugroho, S. K., Kawahara, H., Gibson, N. P., et al. 2021, ApJ, 910, L9 [CrossRef] [Google Scholar]
- Parmentier, V., Line, M. R., Bean, J. L., et al. 2018, A&A, 617, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Piccioni, G., Drossart, P., Zasova, L., et al. 2008, A&A, 483, L29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Quirrenbach, A., Amado, P. J., Caballero, J. A., et al. 2014, SPIE Conf. Ser., 9147, 91471F [Google Scholar]
- Quirrenbach, A., Amado, P. J., Ribas, I., et al. 2018, SPIE Conf. Ser., 10702, 107020W [Google Scholar]
- Rajpurohit, A. S., Allard, F., Rajpurohit, S., et al. 2018, A&A, 620, A180 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Reiners, A., Zechmeister, M., Caballero, J. A., et al. 2018, A&A, 612, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rothman, L. S., Gordon, I. E., Barber, R. J., et al. 2010, J. Quant. Spectr. Rad. Transf., 111, 2139 [Google Scholar]
- Sánchez-López, A., Alonso-Floriano, F. J., López-Puertas, M., et al. 2019, A&A, 630, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sánchez-López, A., López-Puertas, M., Snellen, I. A. G., et al. 2020, A&A, 643, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sánchez-López, A., Landman, R., Mollière, P., et al. 2021, A&A, submitted [Google Scholar]
- Savel, A. B., Kempton, E. M. R., Malik, M., et al. 2021, AAS J. submitted [arXiv:2109.00163] [Google Scholar]
- Seidel, J. V., Ehrenreich, D., Wyttenbach, A., et al. 2019, A&A, 623, A166 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Seidel, J. V., Ehrenreich, D., Pino, L., et al. 2020, A&A, 633, A86 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Seidel, J. V., Ehrenreich, D., Allart, R., et al. 2021, A&A, 653, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Serindag, D. B., Nugroho, S. K., Mollière, P., et al. 2021, A&A, 645, A90 [EDP Sciences] [Google Scholar]
- Snellen, I. A. G., de Kok, R. J., de Mooij, E. J. W., & Albrecht, S. 2010, Nature, 465, 1049 [Google Scholar]
- Snellen, I. A. G., Brandl, B. R., de Kok, R. J., et al. 2014, Nature, 509, 63 [Google Scholar]
- Tabernero, H. M., Zapatero Osorio, M. R., Allart, R., et al. 2021, A&A, 646, A158 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Tamuz, O., Mazeh, T., & Zucker, S. 2005, MNRAS, 356, 1466 [Google Scholar]
- Tsiaras, A., Waldmann, I. P., Zingales, T., et al. 2018, AJ, 155, 156 [NASA ADS] [CrossRef] [Google Scholar]
- von Essen, C., Mallonn, M., Hermansen, S., et al. 2020, A&A, 637, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wardenier, J. P., Parmentier, V., Lee, E. K. H., Line, M. R., & Gharib-Nezhad, E. 2021, MNRAS, 506, 1258 [NASA ADS] [CrossRef] [Google Scholar]
- West, R. G., Hellier, C., Almenara, J. M., et al. 2016, A&A, 585, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Yan, F., Casasayas-Barris, N., Molaverdikhani, K., et al. 2019, A&A, 632, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Zechmeister, M., Anglada-Escudé, G., & Reiners, A. 2014, A&A, 561, A59 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Žák, J., Kabáth, P., Boffin, H. M. J., Ivanov, V. D., & Skarka, M. 2019, AJ, 158, 120 [Google Scholar]
All Tables
All Figures
![]() |
Fig. 1 Volume mixing ratio of H2O and OH using the chemical equilibrium model from Mollière et al. (2017) and an isothermal P-T profile for different temperatures. |
In the text |
![]() |
Fig. 2 Top: telluric transmission as a function of wavelength obtained from ESO SkyCalc (Noll et al. 2012). Middle: transit depth as a function of wavelength, modelled using petitRADTRANS over the entire CARMENES NIR channel and assuming an isothermal atmosphere at the planet equilibrium temperature of 2160 K. The red region indicates the spectral range used for our analysis. Bottom: same but for T = 3200 K. |
In the text |
![]() |
Fig. 3 Resulting CCF for each of the spectra, with the RV in the Earth’s rest frame on the horizontal axis and the planet orbital phase on the vertical axis. The horizontal dashed white lines indicate the start and end of the transit, while the slanted vertical lines show the expected slope and location of the planetary signal. |
In the text |
![]() |
Fig. 4 Bottom: detection significance as a function of RV semi-amplitude, Kp, and velocity relative to the planet rest frame, vrest, showing a ~6σ peak. Top: detection significance at the expected Kp of 196.5 km s−1 and at Kp = 232.2 km s−1, which is obtained from the retrieval. |
In the text |
![]() |
Fig. 5 Resulting posterior distributions for the OH signal in WASP-76b obtained from the MCMC retrieval. The dashed vertical lines indicate the 16th, 50th, and 84th percentiles. |
In the text |
![]() |
Fig. A.1 Cross-correlation S/N with the spectrum of GX And with the OH model template for each of the spectral orders from the CARMENES NIR channel. |
In the text |
![]() |
Fig. B.1 Obtained peak S/N as a function of the number of SysRem iterations applied to the data. |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.