Open Access
Issue
A&A
Volume 665, September 2022
Article Number A104
Number of page(s) 9
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/202243854
Published online 16 September 2022

© G. Guilluy et al. 2022

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe-to-Open model. Subscribe to A&A to support open access publication.

1 Introduction

Ground-based, high-resolution (resolving power R ≥ 20000) spectroscopy (HRS) in the near-infrared (NIR) represents an effective approach to investigate exoplanet atmospheres. Indeed, at high spectral resolution, thousands of resolved molecular lines can be disentangled, as they are Doppler shifted by tens of km s-1 due to the planet motion along its orbit from the telluric (Earth’s atmosphere) and stellar lines that are (quasi-)stationary in wavelength.

Hot (Teq > 1000 K) and warm (Teq ≤ 1000 K) inflated giant planets are the most suitable targets for atmospheric characterisation (Seager & Deming 2010). Recently, Giacobbe et al. 2021) reported the simultaneous detection of six molecules in the atmosphere of one of the best-studied hot-Jupiters (HJs), HD 209458b. As only two molecules (CO and H2O) had been detected at the same time in exoplanetary atmospheres in the past (e.g. Tsiaras et al. 2018), the work presented in Giacobbe et al. (2021) revealed a previously unknown chemical complexity. Thus this discovery raises the question of whether the complexity of HD 209458b’s atmosphere is unique or if other exoatmospheres show such a rich molecular composition. Such molecular diversity enables one to derive atmospheric abundances and thus elemental ratios, such as the carbon-to-oxygen (C/O) ratio, can then be obtained, allowing us to derive important clues on the formation and migration histories of hot giant planets (e.g. Madhusudhan 2012; Turrini et al. 2021). For example, Giacobbe et al. (2021) retrieved a C/O ratio ≥1, meaning that HD 209458b formed far from its present location and subsequently migrated inwards.

With a warm temperature and low density (see Table 1), the close-in transiting planet WASP-69b represents a suitable laboratory for atmospheric characterisation studies. Furthermore, the low equilibrium temperature Teq ≤ 1000 K places the planet in a different regime compared to that of the HJ HD 209458b (Teq ~ 1500 K): the equilibrium chemistry shifts in favour of CH4 rather than CO and disequilibrium phenomena (i.e. vertical mixing and photochemical processes) are more likely to alter the molecular abundances compared to the atmospheres of hotter planets (e.g. Moses 2014).

The atmosphere of WASP-69b has been studied with both space-based low-resolution (LR) and ground-based highresolution (HR) in-transit (transmission spectroscopy) observations. LR observations performed with WFC3 on board the Hubble Space Telescope (HST) unveiled the presence of H2O and possibly aerosols (Tsiaras et al. 2018; Fisher & Heng 2018). Sodium has also been detected through HR optical observations (e.g. Casasayas-Barris et al. 2017). Moreover, HR and LR observations covering the NIR He I (1083.3 nm) revealed that the planet has an extended and escaping atmosphere, where the escape is driven by the intense X-ray and extreme ultraviolet radiation received from the active K5V host star (e.g. Nortmann et al. 2018; Vissapragada et al. 2020).

In this work, we present transmission spectroscopy of WASP-69b using HR observations in the NIR, searching for the absorption of multiple molecular species. We observed the WASP-69 system within the long-term observing programme at the Telescopio Nazionale Galileo (TNG) telescope ‘GAPS2: the origin of planetary systems’ – awarded to the Italian Global Architecture of Planetary System (GAPS) Collaboration (Poretti et al. 2016). Part of this programme is devoted to the characterisation of exoplanetary atmospheres at HR. Briefly, we are carrying out transmission and emission spectroscopy of twenty-six hot Jupiters and four hot Neptunes and/or sub-Neptunes in a wide range of Teq (625–2500 K), and hosted by stars of different spectral types (A to M with K < 10.5 mag), to accomplish the following: (i) identify atomic species to estimate the abundance of refractory elements (Pino et al. 2020) and/or to probe atmospheric expansion or escape, for instance from the Ha (Borsa et al. 2021) and/or metastable helium triplet (Guilluy et al. 2019); (ii) detect molecular compounds and estimate the atmospheric C/O ratio and metallicity, which in turn provide constraints on the planet formation and evolution scenarios (Giacobbe et al. 2021); (iii) derive temperature-pressure profiles (Pino et al. 2020; Borsa et al. 2022); and (iv) detect the atmospheric Rossiter-McLaughlin effect (Borsa et al. 2019; Rainer et al. 2021).

Figure 1 places WASP-69b into the broader context of the GAPS atmospheric survey. The symbol size is proportional to the expected atmospheric signal in the K band: S/N=2HRPR2*F,$S/N = {{2{\rm{H}}{{\rm{R}}_{\rm{P}}}} \over {R_ \star ^2}}*\sqrt F ,$(1)

where we account for the stellar magnitude in K band (F=10magK2.5)$\left( {F = {{10}^ - }{{magK} \over {2.5}}} \right)$, the atmospheric scale height (H), the planetary (RP), and the stellar radius (R*). The most suitable targets to perform atmospheric characterisation studies are labelled. As this plot shows, WASP-69b represents one of the most suitable candidates to perform atmospheric studies (2RPH/R*2$R_*^2$ ~ 283 ppm1, Brown et al. 2001) in the temperature range (≤ 1000 K) we aimed to probe within this study.

This paper is organised as follows: in Sect. 2 we present the data we collected and in Sect. 3 we describe the analysis that we performed to extract the planetary transmission spectrum and, therefore, the molecules responsible for absorption from the raw GIANO-B data. Finally, we discuss our results and draw the conclusions in Sect. 4.

Table 1

Physical and orbital parameters of the WASP-69 system adopted in this work.

thumbnail Fig. 1

Mass-equilibrium temperature relation for the planets within the GAPS atmospheres survey. The symbol size is proportional to the expected atmospheric signal (for one scale height) in K band. The six more suitable targets to perform atmospheric characterisation studies are labelled.

thumbnail Fig. 2

Average S/N across the entire GIANO-B spectral range as a function of the orbital phase, with the marker colour indicating the airmass at which the observations have been gathered. The ingress and egress of WASP-69b (ϕ1 in Table 1) are marked with vertical dashed lines. Black stars on the last night (28 October 2021) indicate the discarded AB couples for low S/N.

2 Data sample

We observed the WASP-69 system with the TNG telescope in GIARPS configuration by simultaneously measuring HR spectra in the optical (0.39−0.69 µm) and NIR (0.95−2.45 µm) with the HARPS-N (R ≈ 115 000) and GIANO-B (R ~ 50000) spectrographs. Our data encompass a total of three primary transits of WASP-69 b, which were obtained at UT 24 July 2019, UT 09 August 2020, and UT 28 October 2021. The GIANO-B observations were taken according to an ABAB nodding pattern (Claudi et al. 2017) for an optimal subtraction of the thermal background noise and telluric emission lines. The observations consist of spectra taken before, during, and after the planetary transit; Fig. 2 shows the signal-to-noise ratio (S/N) averaged over the entire GIANO-B spectral coverage as a function of the planet’s orbital phase and of the airmass for each of the three nights. The observations during the last night (28 October 2021) were affected by several thin clouds (cirri), so we decided to discard the AB couples of observations that exhibit a very low S/N compared to the others if the lowest S/N in the couple is less than 15, see Fig. 2). A log of the observations is reported in Table A.1.

3 Data analysis

Here we report the steps we performed in order to recover the planetary signal from the GIANO-B raw images.

(i) Basic data reduction

We processed the raw GIANO-B spectra with the GOFIO pipeline (Rainer et al. 2018) to perform flat-fielding, bad pixel correction, background subtraction, optimal extraction of the 1D spectra, and a preliminary wavelength calibration using a uranium-neon lamp. We corrected small temporal variations in the wavelength solution due to the imperfect stability of the spectrograph (Giacobbe et al. 2021), and we refined the initial wavelength solution by employing a HR transmission spectrum of the Earth’s atmosphere generated via the ESO Sky Model Calculator2. As some spectral regions have too many or too few telluric lines, this calibration refinement was not possible for all the orders (e.g. Brogi et al. 2018; Guilluy et al. 2019; Fossati et al. 2022). In Table A.2, we report the discarded orders for each individual night of observation.

(ii) Telluric and stellar lines’ removal

We then removed the telluric contamination and the stellar spectrum by employing the principal component analysis (PCA) and linear regression techniques (Giacobbe et al. 2021). We selected the optimal number of components Nopt following the procedure described in Giacobbe et al. (2021). More precisely, Nopt varies with the considered order depending on the root mean square (rms) variation. We selected the number for which the first derivative between the last and penultimate component decreases by less than σwhiteNopt1/2${\sigma _{{\rm{white}}}}N_{{\rm{opt}}}^{ - {1 \mathord{\left/ {\vphantom {1 2}} \right. \kern-\nulldelimiterspace} 2}}$, where σwhite is the standard deviation of the full matrix assuming pure white noise. Since in two nights out of three (i.e. UT 24 July 2019 and UT 09 August 2020) the relative velocity shift between telluric and planetary lines was close to zero, we masked the telluric components to avoid contamination in our final spectra. Following Brogi et al. (2018), we thus masked the data assigning zero flux to those spectral channels corresponding to telluric lines deeper than 20%. Finally, we visually inspected the outcome for spurious telluric residuals.

(iii) Search for molecules through cross-correlation

For each considered molecule, we combined the contribution of thousands of lines in the spectral coverage of GIANO-B by performing the cross-correlation (CC) of the residual spectra with isothermal transmission models generated using the GENESIS code (Gandhi & Madhusudhan 2017) adapted for transmission spectroscopy (Pinhas et al. 2018; see Table 2 for the adopted line list). Prior to CC, the models were convolved to the GIANO-B instrument profile (a Gaussian profile with FWHM 5.4 km s−1). Models were calculated between the 100 and 10−8 bar in pressure, at a Teq = 963(18)K, between 0.9 and 2.6 µm, with a constant spacing of 0.01 cm−1. We included collision-induced absorption from H2–H2 and H2–He interactions (Richard et al. 2012). For each night we performed the CC for an optimal selection of GIANO-B orders (Giacobbe et al. 2021) and every phase over a lag vector corresponding to planet radial velocities (RVs) in the range −252 ≤ RV ≤ 252 km s−1, in steps of 3 km s−1. This range was chosen in order to cover all possible RVs of the planet. We then co-added the CC functions over the selected orders, nights, and the orbital phase after shifting them in the planet rest frame by assuming a circular orbit (e.g. Brogi et al. 2012, 2014, 2018). In doing so, we considered a range of planet RV semi-amplitudes 0 ≤ Kp ≤ 200 kms-1, in steps of 1.5 km s-1. We then quantified the confidence level of our detection by computing the S/N detection map (e.g. Brogi et al. 2012, 2014, 2018). For each investigated molecule, we divided the total CC matrix by its standard deviation (calculated by excluding the CC peak). We also used the Welch t-test to compare the ‘in-trail’ values (CC values that carry signal) of the 2D CC matrix aligned in the planet rest frame to those ‘out-of-trail’ (CC values away from the WASP-69b’s RV). We then converted the corresponding t value into the statistical significance at which these two distributions are not drawn from the same parent distribution. We applied this method for the same range of rest-frame velocity Vrest and the planetary semi-amplitude KP adopted before. The final S/N and significance maps for the three nights combined together are shown in the upper and in the middle panels of Fig. 3, respectively. These results show the advantage of a multi-night approach; indeed even if the data quality does not allow one to always detect the signal in each individual night (see Fig. A.1), the detection significance increases when transits are co-added.

The detections in CC have been obtained with isothermal models of one molecule at a time with a volume mixing ratio (VMR) large enough to reveal a given molecule, if present in the atmosphere; it has indeed been demonstrated (e.g. Gandhi et al. 2020) that the CC is only weakly dependent on the absolute VMR. As a good compromise, we found a VMR ~ 10−4.

Table 2

CC and LH results.

(iv) Confidence intervals via likelihood

The CC framework as applied in the previous analysis implicitly assumes that the signal amplitude and the S/N are uniform across the considered spectral range. However, this assumption is only valid at a first-order level. To investigate this aspect more thoroughly, we converted the CC values into LikeliHood (LH) mapping values using the approach proposed by Brogi & Line (2019) and used by Giacobbe et al. (2021). In this formalism, the LH is indeed set up by taking into account the model’s line depth and the S/N order-by-order and spectrum-by-spectrum. Retrieving a peak in the LH at the correct position of the (Vrest, KP) space is hence an additional confirmation of our detection. We thus computed a log-LH function for each selected order, spectrum, and RV shift of the model across the KP versus Vrest space. We used a lag vector corresponding to planet RVs in the range −20 ≤ RV ≤ 20 km s−1, in steps of 1 km s−1. To take into account any possible modification of the planetary signal due to the PCA process, on each theoretical model, we performed the exact telluric removal process applied to the observations. The final log-LH for each (Vrest, KP) couple was obtained by adding all the log-LHs over each order, night, and observed spectrum.

Subsequently, as in Giacobbe et al. (2021), we introduced the line-intensity scaling factor S (−1.1 < log S < 1.1, in steps of 0.1) to understand how much the used models have to be scaled to match the observed data. For a model perfectly matching the data, S is exactly 1 (log S = 0), that is to say the model and the data have the same average line amplitude compared to the local continuum. Scaling factors smaller than 1 indicate that spectral lines are too strong compared to the observations and, for example, the presence of a grey-cloud model can be invoked to account for that, as in Giacobbe et al. (2021). Furthermore, as in the single species models, the continuum can be approximate, and the introduction of the log S allows us to account for the muting effects of other species. On the contrary, in the CC framework, since the line depth does not enter into the calculation, the log S is not necessary. We then converted the LH values into confidence intervals by using the LH-ratio test between the maximum log-LH and each LH value in the (Vrest, KP, and log S) space. The LH confidence intervals’ maps corresponding to the log S, maximising the detection, are shown in the bottom panel of Fig. 3.

We then evaluated the significance of our LH detection by performing a LH-ratio test between the maximum log-LH and the median of the 2D log-LH matrix values for the best-fit scaling factor which exhibited a statistical confidence greater than 1. For each considered molecule, Table 2 reports the results of our analysis in both the CC and the LH frameworks. We consider a molecule as detected if altogether the S/N, the detection significance, and the LH significance are higher than 3σ.

4 Results and discussion

We characterise the atmosphere of WASP-69b showing a multi-species molecular detection in the planet transmission spectrum. Our results indicate the presence of CH4, NH3, CO, C2H2, and H2O, at more than a 3.3σ level (Table 2). We did not detect the presence of HCN or CO2 with a confidence level higher than 3σ (see Table 2). The retrieved Kp values for all molecules are compatible within 1σ with the nominal Kp value in the S/N and t-test maps (see Table 2 and Fig. 3, top and middle panels). We find that the same thing occurs for NH3, CO, and C2H2 when switching to the LH framework, while the peaks of the CH4 and H2O signals are slightly shifted, though consistent with the expected Kp at the 2σ level (see Fig. 3, bottom panel). Since the LH function is more sensitive to the variance of the spectra, these small shifts may be due to the contribution of the other molecules, which are not modelled simultaneously when considering one molecule at a time.

We also checked whether our CO detection might be spurious, as stellar CO lines can still be an important contaminant of the planetary transmission spectrum (Brogi et al. 2016). To resolve the possible ambiguity, we evaluated the expected KP at which a spurious signal of stellar origin should have been maximised as KPstar~vsinisin(2πϕ1)~${K_{P{\rm{star}}}}\~{{v\sin i} \over {\sin \left( {2\pi {\phi _1}} \right)}}\~$29.2 ± 5.3 km s−1 (by using the values in Table 1). While the above equation provides an estimate of the expected position of stellar residuals in the (Vrest, KP) map, the actual position will depend on the interplay between the stellar and the planetary signal, and it might also be affected by any partial removal of stellar residuals via PCA. Thus, one should not expect stellar residuals to arise exactly at the KPstar computed above. In fact, there are cases in the literature where stellar CO causes the planet signal to be detected at higher KP than the orbit of the planet allows (e.g. Chiavassa & Brogi 2019, Fig. 7), which seems counter-intuitive at first.

With the LH framework applied, we could also infer to which extent the model is a plausible representation of the data. However, we measured low log S values for several molecules, as single-species models have a higher variance (deeper lines) than an actual spectrum, because they do not contain the additional opacity from other species. We remark that the assumption of a fixed VMR composition has an impact on the line strengths S, and thus it can lead to ambiguities in the scaling factors’ interpretation. Consequently, the comparison between the S values for the different molecules is not straightforward, meaning that the relative abundances between the different species cannot be estimated easily. The scaling factors obtained here (Table 2) thus must be considered representative and taken with caution.

The large number of individual species’ detections (and non-detections) on WASP-69b provides important constraints on the physical processes that govern the planet’s atmosphere. Using the PYRAT BAY (Cubillos & Blecic 2021) modelling framework (see Appendix B), we explored physically-plausible atmospheric scenarios in radiative and thermochemical equilibrium that could be consistent with the observed molecular detections. Given the striking difference in composition that is expected for scenarios with elemental ratios of C/O < 1 or C/O > 1 (see, e.g. Tsuji 1973; Madhusudhan 2012; Moses et al. 2013), we adopted the C/O ratio and the atmospheric metallicity as the main proxies to explore different plausible compositions for WASP-69b. We found that very high-metallicity regimes (>10× solar) appear disfavoured by the observations because, for C/O < 1 scenarios, the non-detected CO2 molecule would present detectable features since its abundance increases faster with metallicity than for the other species. For C/O > 1 scenarios at high metallicities, the observed H2O features would not be detected since they would be obscured by the CH4 absorption, due to the much lower abundance of H2O relative to CH4. When comparing solar metallicity models with C/O ratios lower and greater than one (Fig. 4), the observations seem to be more consistent with the latter. First, while transmission models for both regimes show clear spectral features of H2O, CH4, and CO, the stronger detection of CH4 than H2O may favour the C/O > 1 scenario (in which CH4 is predominant). Then, the NH3 molecule presents a wide-range yet weak features across the spectrum in both C/O regimes, thus, it does not particularly favour any of the scenarios. However, the C2H2 is hard to explain in equilibrium conditions since its estimated abundance does not yield detectable features in any of the above scenarios. A way to resolve this conundrum is by including disequilibrium chemistry (see e.g. Moses 2014, and references therein), which would significantly enhance the abundances of NH3 and C2H2 by dragging material from the deep interior to the upper atmospheric layers (vertical quenching), by photochemical production, or both. Certainly, the much greater abundance of C2H2 for the C/O > 1 case makes the detection of this molecule more plausible for this scenario than for C/O < 1 (a difference of ~5 orders of magnitude for the cases shown in Fig. 4), as long as disequilibrium processes help to bridge the gap between the WASP-69b detections and the equilibrium theoretical predictions.

To give a more quantitative interpretation of the physical and chemical conditions of WASP-69b’s atmosphere, we computed and tested three different scenarios that could be a possible representation of WASP-69b’s atmosphere:

  • (i) solar metallicity and solar C/O (C/O = 0.55);

  • (ii) solar metallicity and carbon rich (C/O = 2.0); and

  • (iii) solar metallicity, carbon rich (C/O = 2.0), and quenched C2H2 and NH3 to VMRs of 10−6 −10−7.

As opposed to the single-molecule analyses from Sect. 3, here we generated more realistic ‘full spectrum’ models combining opacities from all species. To quantify which scenario describes the observations better, we applied the Wilks theorem (Wilks 1938) by computing the LH-ratio test relative to the LH-maxima (following the procedure of Giacobbe et al. 2021). Our analysis (see Table A.3) indicates that scenario (iii) is marginally (~1σ) preferred compared to scenario (ii), and it is significantly favoured (~3.5σ) over scenario (i). The atmosphere of WASP-69 thus seems to be better described by a scenario where quenching of NH3 and CH4 is taking place. However, we need to stress two aspects: (1) this result does not imply that (iii) is the best-fitting scenario, but it is the preferred one within the tested models; and (2) we assumed a fixed solar metallicity, which might not be the most representative of the planet’s atmosphere. Dealing with (1) and (2) implies a retrieval analysis that goes beyond the scope of this paper.

The emerging picture provides some preliminary indications on WASP-69b’s formation history once put in the context of the host star. WASP-69’s metallicity is super-solar (Anderson et al. 2014), with recent works suggesting metallicity values 2–3 times higher than that of the Sun (Magrini et al. 2022). The atmospheric scenario favoured by our modelling efforts points to a solar metallicity for WASP-69b’s atmosphere, meaning that the atmospheric metallicity is sub-stellar. The combination of substellar planetary metallicity and planetary C/O > 1 is the signpost of giant planets whose content of heavy elements is dominated by the accretion of gas with a limited contribution of solids (Turrini et al. 2022, 2021). Furthermore, this combination suggests that the chemical structure of WASP-69b’s native circumstellar disc was inherited from the parent molecular cloud (Pacetti et al. 2022).

We also note that cloud formation processes offer an alternative interpretation of these results. In an atmosphere with prominent cloud formation, condensation locks and transports heavy elements such as carbon or oxygen across an atmosphere, which can significantly modify the local gas-phase carbon-to-oxygen ratios. For example, for a solar-like composition, oxygen depletion due to condensation can be strong enough to raise C/O ratios up to ≈0.7 (Bilger et al. 2013). For carbon-rich atmospheres, on the other hand, carbon depletion can make an atmosphere locally more oxygen rich, approaching values of C/O ≈ 1.0 at the regions with the strongest carbon depletion (Helling et al. 2017).

While the picture described above will need to be validated by future studies, the simultaneous presence of several chemical species constitutes a very important step forward, as so far only a few molecules have been firmly detected simultaneously in the atmosphere of warm giant planets (e.g. Tsiaras et al. 2018; Fisher & Heng 2018). Our findings reinforce the notion that a rich atmospheric chemistry is not the sole dominion of hotter planets such as HD 209458b (Giacobbe et al. 2021). The frontier in the characterisation of exoplanetary atmospheres is expanding further, with even better prospects in the future when stronger constraints on molecular abundances, and hence on metallicity and elemental ratios (e.g. Brogi et al. 2017; Brogi & Line 2019; Guilluy et al. 2022), will be placed based on systematic joint analyses of HR ground-based observations gathered with GIANO-B, CARMENES, Spirou, NIRPS, and CRIRES+, for instance, and LR space-borne spectra such as those obtained with HST, JWST, and Ariel.

thumbnail Fig. 3

CC both in the S/N (top panel) and in the significance (middle panel) framework, and LH confidence intervals (bottom panel) maps for each investigated molecule as a function of the planet’s maximum RV (KP) and the planet’s rest-frame velocity (Vrest). The dotted lines denote the expected RV semi-amplitude of WASP-69b (Kp=127.111.52+1.49km s1)$\left( {{K_{\rm{p}}} = 127.11_{ - 1.52}^{ + 1.49}{\rm{km }}{{\rm{s}}^{ - 1}}} \right)$ in the planet rest frame (Vrest = 0 km s−1). The detection significance increased when transits were co-added.

thumbnail Fig. 4

Sample atmospheric models for WASP-69b in radiative and thermochemical equilibrium. Left panel: equilibrium temperature profile for two models with solar elemental metallicity and a solar C/O ratio (solid green curve) or a super-solar ratio (dashed blue curve). Middle panel: equilibrium composition for the respective solar C/O model (solid curves) and super-solar C/O model (dashed curves), colour coded for each relevant species (see legend). Right panel: post-produced HR transmission spectra over the GIANO-B wavelength range, for the solar C/O model (green) and super-solar C/O model (blue). The horizontal black lines denote the wavelength ranges where each of the labelled species significantly impact the spectra.

Acknowledgements

We thank the referee, Dr. Michael Line for his insightful comments which helped to improve the quality of our work. We acknowledge financial contributions from PRIN INAF 2019, and from the agreement ASI-INAF number 2018-16-HH.0 (THE StellaR PAth project). P.C. was funded by the Austrian Science Fund (FWF) Erwin Schroedinger Fellowship J4595-N. D.T., E.S. and E.P. acknowledge the support of the Italian National Institute of Astrophysics (INAF) through the INAF Main Stream project “Ariel and the astrochemical link between circumstellar discs and planets” (CUP: C54I19000700005), and of the Italian Space Agency through the ASI-INAF contract no. 2021-5-HH.0. E.P. also acknowledges support from the European Research Council via the Horizon 2020 Framework Programme ERC Synergy “ECOGAL” Project GA-855130. M.B. ackowledges support from from the UK Science and Technology Facilities Council (STFC) research grant ST/T000406/1.

Appendix A Additional figures and tables

Table A.1

Log of the WASP-69b GIANO-B observations.

Table A.2

Discarded orders in the basic data reduction step of the analysis, i.e. alignment in the telluric rest frame and refinement of the wavelength solution.

Table A.3

Comparison of atmospheric ‘full models’.

thumbnail Fig. A.1

S/N maps for the detected molecules for each night considered individually. The CC peaks are indicated by a black cross. As discussed in Giacobbe et al. (2021), species have weaker signatures and are not always detected in each individual transit. Nevertheless, they are firmly detected when more transits are co-added. These results show the advantage of our multi-night approach.

Appendix B The Pyrat Bay code

To compute atmospheric profiles in radiative and thermochemical equilibrium, we employed the Pyrat Bay atmospheric modelling code (Cubillos & Blecic 2021). To evaluate the layer-by-layer fluxes, we implemented the two-stream radiative-transfer scheme of Heng et al. (2014, Appendix B). This calculation does not consider scattering or feedback from the formation of condensates.

Following Malik et al. (2017), the code attains radiative equilibrium by iteratively updating the temperature profile, such that the divergence of net flux traversing each layer of the atmosphere converges to a negligible value. After each temperature update, the PYRAT BAY code re-evaluates the atmospheric composition using the thermochemical-equilibrium abundances code TEA (Blecic et al. 2016). This is done for a reduced chemical network consisting of H-, He-, C-, N-, O-, Na-, and K-bearing species. At the top of the atmosphere, we imposed an incident stellar irradiation from a Kurucz model (Castelli & Kurucz 2003) according to the physical properties of WASP-69, assuming zero Bond albedo and full day-night energy redistribution. At the bottom of the atmosphere, we imposed an internal radiative heat corresponding to a 100 K blackbody.

The atmospheric model spans over a fixed pressure range from 100 to 10−9 bar, and a wavelength grid ranging from 0.3 to 30 µm sampled with at a resolving power of R = 15 000, which is sufficient to contain the majority of the stellar optical radiation and planetary infrared radiation. The opacity sources include line-list data for the seven molecules searched in Section 3: CO, CO2, and CH4 from HITEMP (Rothman et al. 2010; Li et al. 2015; Hargreaves et al. 2020); and H2O, HCN, NH3, and C2H2 from ExoMol (Polyansky et al. 2018; Chubb et al. 2020; Yurchenko et al. 2011; Harris et al. 2006, 2008; Coles et al. 2019). Because of the large size of the ExoMol opacities (several billions of line transitions), we employed the REPACK algorithm (Cubillos 2017) to extract only the dominant transitions, thus reducing the number of transitions by a factor of ~100 without a significant impact on the resulting opacities. Additionally, the model included Na and K resonance-line opacities (Burrows et al. 2000); Rayleigh opacity for H, H2, and He (Kurucz 1970); and collision-induced absorption for H2–H2 and H2–He (Borysow et al. 2001; Borysow 2002; Richard et al. 2012). We considered a few broad regimes by varying two defining properties: elemental metallicities sub-solar to super-solar versus C/O ratios lower and greater than one. Once we estimated these equilibrium models, we post-produced transmission spectra, generated at a high resolving power (R = 80 000), and then convolved to the GIANO-B instrumental resolution (R ≈ 50 000), which we inspected to determine which species show detectable features over the observed wavelength range.

References

  1. Anderson, D. R., Collier Cameron, A., Delrez, L., et al. 2014, MNRAS, 445, 1114 [NASA ADS] [CrossRef] [Google Scholar]
  2. Barber, R. J., Strange, J. K., Hill, C., et al. 2013, MNRAS, 437, 1828 [Google Scholar]
  3. Bilger, C., Rimmer, P., & Helling, C. 2013, MNRAS, 435, 1888 [NASA ADS] [CrossRef] [Google Scholar]
  4. Blecic, J., Harrington, J., & Bowman, M. O. 2016, ApJS, 225, 4 [CrossRef] [Google Scholar]
  5. Bonomo, A. S., Desidera, S., Benatti, S., et al. 2017, A&A, 602, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Borsa, F., Rainer, M., Bonomo, A. S., et al. 2019, A&A, 631, A34 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  7. Borsa, F., Lanza, A. F., Raspantini, I., et al. 2021, A&A, 653, A104 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Borsa, F., Giacobbe, P., Bonomo, A. S., et al. 2022, A&A, 663, A141 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Borysow, A. 2002, A&A, 390, 779 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Borysow, A., Jorgensen, U. G., & Fu, Y. 2001, J. Quant. Spectr. Rad. Transf., 68, 235 [NASA ADS] [CrossRef] [Google Scholar]
  11. Brogi, M., & Line, M. R. 2019, AJ, 157, 114 [Google Scholar]
  12. Brogi, M., Snellen, I. A. G., de Kok, R. J., et al. 2012, Nature, 486, 502 [Google Scholar]
  13. Brogi, M., de Kok, R. J., Birkby, J. L., Schwarz, H., & Snellen, I. A. G. 2014, A&A, 565, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Brogi, M., de Kok, R. J., Albrecht, S., et al. 2016, ApJ, 817, 106 [Google Scholar]
  15. Brogi, M., Line, M., Bean, J., Désert, J.-M., & Schwarz, H. 2017, ApJ, 839, L2 [NASA ADS] [CrossRef] [Google Scholar]
  16. Brogi, M., Giacobbe, P., Guilluy, G., et al. 2018, A&A, 615, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Brown, T. M., Charbonneau, D., Gilliland, R. L., Noyes, R. W., & Burrows, A. 2001, ApJ, 552, 699 [NASA ADS] [CrossRef] [Google Scholar]
  18. Burrows, A., Marley, M. S., & Sharp, C. M. 2000, ApJ, 531, 438 [NASA ADS] [CrossRef] [Google Scholar]
  19. Casasayas-Barris, N., Palle, E., Nowak, G., et al. 2017, A&A, 608, A135 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Castelli, F., & Kurucz, R. L. 2003, IAU Symp., 210, A20 [Google Scholar]
  21. Chiavassa, A., & Brogi, M. 2019, A&A, 631, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Chubb, K. L., Tennyson, J., & Yurchenko, S. N. 2020, MNRAS, 493, 1531 [NASA ADS] [CrossRef] [Google Scholar]
  23. Claudi, R., Benatti, S., Carleo, I., et al. 2017, Euro. Phys. J. Plus, 132, 364 [NASA ADS] [CrossRef] [Google Scholar]
  24. Coles, P. A., Yurchenko, S. N., & Tennyson, J. 2019, MNRAS, 490, 4638 [CrossRef] [Google Scholar]
  25. Cubillos, P. E. 2017, ApJ, 850, 32 [CrossRef] [Google Scholar]
  26. Cubillos, P. E., & Blecic, J. 2021, MNRAS, 505, 2675 [NASA ADS] [CrossRef] [Google Scholar]
  27. Fisher, C., & Heng, K. 2018, MNRAS, 481, 4698 [Google Scholar]
  28. Fossati, L., Guilluy, G., Shaikhislamov, I. F., et al. 2022, A&A, 658, A136 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Gandhi, S., & Madhusudhan, N. 2017, MNRAS, 472, 2334 [Google Scholar]
  31. Gandhi, S., Brogi, M., & Webb, R. K. 2020, MNRAS, 498, 194 [Google Scholar]
  32. Giacobbe, P., Brogi, M., Gandhi, S., et al. 2021, Nature, 592, 205 [NASA ADS] [CrossRef] [Google Scholar]
  33. Guilluy, G., Sozzetti, A., Brogi, M., et al. 2019, A&A, 625, A107 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Guilluy, G., Sozzetti, A., Giacobbe, P., Bonomo, A. S., & Micela, G. 2022, Exp. Astron., 53, 655 [NASA ADS] [CrossRef] [Google Scholar]
  35. Hargreaves, R. J., Gordon, I. E., Rey, M., et al. 2020, ApJS, 247, 55 [NASA ADS] [CrossRef] [Google Scholar]
  36. Harris, G. J., Tennyson, J., Kaminsky, B. M., Pavlenko, Y. V., & Jones, H. R. A. 2006, MNRAS, 367, 400 [Google Scholar]
  37. Harris, G. J., Larner, F. C., Tennyson, J., et al. 2008, MNRAS, 390, 143 [NASA ADS] [CrossRef] [Google Scholar]
  38. Helling, C., Tootill, D., Woitke, P., & Lee, E. 2017, A&A, 603, A123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  39. Heng, K., Mendonça, J. M., & Lee, J.-M. 2014, ApJS, 215, 4 [NASA ADS] [CrossRef] [Google Scholar]
  40. Høg, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27 [Google Scholar]
  41. Huang, X., Schwenke, D., Freedman, R., & Lee, T. 2017, J. Quant. Spectr. Rad. Transf., 203, 224 [CrossRef] [Google Scholar]
  42. Kurucz, R. L. 1970, SAO Special Report, 309 [Google Scholar]
  43. Li, G., Gordon, I. E., Rothman, L. S., et al. 2015, ApJS, 216, 15 [NASA ADS] [CrossRef] [Google Scholar]
  44. Madhusudhan, N. 2012, ApJ, 758, 36 [NASA ADS] [CrossRef] [Google Scholar]
  45. Magrini, L., Danielski, C., Bossini, D., et al. 2022, A&A, 663, A161 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Malik, M., Grosheintz, L., Mendonça, J. M., et al. 2017, AJ, 153, 56 [Google Scholar]
  47. Moses, J. I. 2014, Phil. Trans. R. Soc. A Math. Phys. Eng. Sci., 372, 20130073 [CrossRef] [Google Scholar]
  48. Moses, J. I., Madhusudhan, N., Visscher, C., & Freedman, R. S. 2013, ApJ, 763, 25 [NASA ADS] [CrossRef] [Google Scholar]
  49. Nortmann, L., Pallé, E., Salz, M., et al. 2018, Science, 362, 1388 [Google Scholar]
  50. Pacetti, E., Turrini, D., Schisano, E., et al. 2022, ApJ, accepted ArXiv e-prints [arXiv:2206.14685] [Google Scholar]
  51. Pinhas, A., Rackham, B. V., Madhusudhan, N., & Apai, D. 2018, MNRAS, 480, 5314 [Google Scholar]
  52. Pino, L., Désert, J.-M., Brogi, M., et al. 2020, ApJ, 894, L27 [Google Scholar]
  53. Polyansky, O. L., Kyuberis, A. A., Zobov, N. F., et al. 2018, MNRAS, 480, 2597 [NASA ADS] [CrossRef] [Google Scholar]
  54. Poretti, E., Boccato, C., Claudi, R., et al. 2016, Mem. Soc. Astron. It., 87, 141 [NASA ADS] [Google Scholar]
  55. Rainer, M., Harutyunyan, A., Carleo, I., et al. 2018, SPIE Conf. Ser., 10702, 1070266 [Google Scholar]
  56. Rainer, M., Borsa, F., Pino, L., et al. 2021, A&A, 649, A29 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  57. Richard, C., Gordon, I., Rothman, L., et al. 2012, J. Quant. Spectr. Rad. Transf., 113, 1276 [NASA ADS] [CrossRef] [Google Scholar]
  58. Rothman, L. S., Gordon, I. E., Barber, R. J., et al. 2010, J. Quant. Spectr. Rad. Transf., 111, 2139 [Google Scholar]
  59. Seager, S., & Deming, D. 2010, ARA&A, 48, 631 [NASA ADS] [CrossRef] [Google Scholar]
  60. Tsiaras, A., Waldmann, I. P., Zingales, T., et al. 2018, ApJ, 155, 156 [Google Scholar]
  61. Tsuji, T. 1973, A&A, 23, 411 [NASA ADS] [Google Scholar]
  62. Turrini, D., Schisano, E., Fonte, S., et al. 2021, ApJ, 909, 40 [Google Scholar]
  63. Turrini, D., Codella, C., Danielski, C., et al. 2022, Exp. Astron., 53, 225 [NASA ADS] [CrossRef] [Google Scholar]
  64. Vissapragada, S., Knutson, H. A., Jovanovic, N., et al. 2020, AJ, 159, 278 [Google Scholar]
  65. Wilks, S. S. 1938, Ann. Math. Stat., 9, 60 [Google Scholar]
  66. Yurchenko, S. N., Barber, R. J., & Tennyson, J. 2011, MNRAS, 413, 1828 [Google Scholar]

1

For comparison, HD 209458 b investigated in Giacobbe et al. (2021) has a 2RPH/R*2${{2{R_P}H} \mathord{\left/ {\vphantom {{2{R_P}H} {R_*^2}}} \right. \kern-\nulldelimiterspace} {R_*^2}}$ ~ 175 ppm.

All Tables

Table 1

Physical and orbital parameters of the WASP-69 system adopted in this work.

Table 2

CC and LH results.

Table A.1

Log of the WASP-69b GIANO-B observations.

Table A.2

Discarded orders in the basic data reduction step of the analysis, i.e. alignment in the telluric rest frame and refinement of the wavelength solution.

Table A.3

Comparison of atmospheric ‘full models’.

All Figures

thumbnail Fig. 1

Mass-equilibrium temperature relation for the planets within the GAPS atmospheres survey. The symbol size is proportional to the expected atmospheric signal (for one scale height) in K band. The six more suitable targets to perform atmospheric characterisation studies are labelled.

In the text
thumbnail Fig. 2

Average S/N across the entire GIANO-B spectral range as a function of the orbital phase, with the marker colour indicating the airmass at which the observations have been gathered. The ingress and egress of WASP-69b (ϕ1 in Table 1) are marked with vertical dashed lines. Black stars on the last night (28 October 2021) indicate the discarded AB couples for low S/N.

In the text
thumbnail Fig. 3

CC both in the S/N (top panel) and in the significance (middle panel) framework, and LH confidence intervals (bottom panel) maps for each investigated molecule as a function of the planet’s maximum RV (KP) and the planet’s rest-frame velocity (Vrest). The dotted lines denote the expected RV semi-amplitude of WASP-69b (Kp=127.111.52+1.49km s1)$\left( {{K_{\rm{p}}} = 127.11_{ - 1.52}^{ + 1.49}{\rm{km }}{{\rm{s}}^{ - 1}}} \right)$ in the planet rest frame (Vrest = 0 km s−1). The detection significance increased when transits were co-added.

In the text
thumbnail Fig. 4

Sample atmospheric models for WASP-69b in radiative and thermochemical equilibrium. Left panel: equilibrium temperature profile for two models with solar elemental metallicity and a solar C/O ratio (solid green curve) or a super-solar ratio (dashed blue curve). Middle panel: equilibrium composition for the respective solar C/O model (solid curves) and super-solar C/O model (dashed curves), colour coded for each relevant species (see legend). Right panel: post-produced HR transmission spectra over the GIANO-B wavelength range, for the solar C/O model (green) and super-solar C/O model (blue). The horizontal black lines denote the wavelength ranges where each of the labelled species significantly impact the spectra.

In the text
thumbnail Fig. A.1

S/N maps for the detected molecules for each night considered individually. The CC peaks are indicated by a black cross. As discussed in Giacobbe et al. (2021), species have weaker signatures and are not always detected in each individual transit. Nevertheless, they are firmly detected when more transits are co-added. These results show the advantage of our multi-night approach.

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.