Free Access
Issue
A&A
Volume 663, July 2022
Article Number A160
Number of page(s) 12
Section Planets and planetary systems
DOI https://doi.org/10.1051/0004-6361/202142559
Published online 28 July 2022

© ESO 2022

Introduction

GJ 436 was the first M dwarf across which the Rossiter-McLaughlin (RM) effect was detected, leading to the measurement of a highly misaligned orbit (Bourrier et al. 2018b) for its evaporating (Kulow et al. 2014; Ehrenreich et al. 2015; Lavie et al. 2017; dos Santos et al. 2019) warm Neptune companion. This peculiar architecture supports a late Kozai migration (Naoz 2016) induced by an outer companion (Bourrier et al. 2018b), a scenario that was originally proposed by Beust et al. (2012) to explain the current nonzero orbital eccentricity of GJ 436b (e = 0.152, Trifonov et al. 2018), which should otherwise have been quickly circularized by tidal interactions with the star (Mardling 2008). In a Kozai scenario, a distant perturbing body traps the inner planet in a secular dynamical resonance far away from the star for several billion years, only to subsequently decouple and migrate inwards on a misaligned and eccentric orbit (e.g., Fabrycky & Tremaine 2007; Anderson et al. 2016). GJ 436b might be a representative example of Neptune-size planets that started to evaporate long after their formation because of this delayed migration (Bourrier et al. 2018b; Owen & Lai 2018; Correia et al. 2020; Attia et al. 2021). Albrecht et al. 2021 recently presented statistical evidence for a preponderance of exoplanets on nearly polar orbits, several of which have noncircular orbits, lie at the edge of the hot Neptune desert, and are observed to be evaporating. Noticeable examples with known planetary companions include HAT-P-11 b (see Yee et al. 2018; Ben-Jaffel et al. 2022) and WASP-107 b (see Spake et al. 2018, Spake et al. 2021; Allart et al. 2019; Kirk et al. 2020; Piaulet et al. 2021). GJ 3470 b is also worth mentioning as it exhibits all of the aforementioned characteristics of a delayed Kozai migration (see Bourrier et al. 2018a; Palle et al. 2020; Ninan et al. 2020; Stefansson et al. 2022), although, like GJ 436b, its putative companion remains undetected. While secular resonance crossing driven by the protoplanetary disk (Petrovich et al. 2020) can also account for the eccentric and misaligned orbits of low-mass planets around the desert, it requires that the planet formed on a close-in orbit and thus does not explain the long-term survival of their evaporating atmospheres. The above systems thus support the critical role that Kozai resonance could play in forming the desert of hot Neptunes, by delaying the arrival at their close-in location and their subsequent evaporation. The case of GJ 436b illustrates the utility of transit spectroscopy in determining the spin–orbit angles of close-in planets (see the review by Triaud 2018), in constraining the mechanisms that altered their orbital architecture, and in unraveling their complex evolution.

Here we analyze the RM effect of GJ 436b in two transits recently observed with ESPRESSO (Sect. 2). We compare the results obtained using the classical analysis of the RV anomaly (Sect. 3) with those obtained using the recently developed Reloaded RM and RM Revolutions techniques, analyzing the stellar surface properties along the transit chord (Sect. 4). Using the latter technique, we perform a joint fit to the ESPRESSO and archival HARPS/HARPS-N data to refine the architecture of the GJ 436 system (Sect. 5).

2 Observations and Data Reduction

We exploit two datasets obtained with ESPRESSO during the transit of GJ 436b on 27 February 2019 (Visit 1) and 29 April 2019 (Visit 2). During each visit, 49 exposures were obtained at about the same orbital phases, with 10 in-transit exposures, 17 before the transit in Visit 1 (respectively 16 in Visit 2), and 22 after the transit in Visit 1 (respectively 23 in Visit 2). GJ 436 was observed on fiber A, while fiber B was used to monitor sky contamination. Because of the faintness of the star in the optical (V = 10.6), it was observed in ESPRESSO high-resolution mode (R ~ 140 000) with a binning of 2 × 1 to reduce readout noise (instrumental mode HR21, for more information see Pepe et al. 2021). Spectra were extracted from the detector images, corrected, and calibrated using version 2.2.8 of the Data Reduction Software (DRS) pipeline (Pepe et al. 2021). One of the DRS corrections concerns the color effect caused by the variability of extinction induced by Earth’s atmosphere (e.g., Bourrier & Hébrard 2014; Bourrier et al. 2018b; Wehbe et al. 2020). The flux balance of the GJ 436 spectra was reset to a M3 stellar spectrum template before they were passed through weighted cross-correlation (Baranne et al. 1996; Pepe et al. 2002) with an M3 numerical mask to compute cross-correlation functions (CCFs) with a step of 0.5 km s−1. This mask is part of a new set that were built with weights more representative of the photonic error on the line positions, as described in Bourrier et al. (2021).

Properties of the GJ 436 system used and derived in our analyses are reported in Table 1. Recent CHEOPS observations of GJ 436b (Maxted 2022) improve the precision on mid-transit times at the epochs of our observations to 8 s, removing potential biases associated with ephemeris (e.g., Casasayas-Barris et al. 2021). The planet-to-star radius ratio was fixed to the error-weighted average of values obtained at different precisions with various instruments by Lanotte et al. (2014), Knutson et al. (2014), Lothringer et al. (2018), and Maxted (2022). None of these studies found any strong evidence for variations in transit depth with wavelength.

Throughout the paper, the posterior probability distributions (PDFs) of free parameters describing models fitted to the data are sampled using emcee MCMC (Foreman-Mackey et al. 2013), adjusting the number of walkers and the burn-in phase based on the degrees of freedom of the considered problem and the convergence of the chains. Best-fit values for these parameters are set to the median of their PDFs, and their 1σ uncertainty ranges are defined using highest density intervals.

Table 1

Properties of the GJ 436 system.

3 Analysis of the Disk-Integrated CCFs

3.1 Stellar Line Properties

The CCFDI (for “disk-integrated”) produced by the DRS correspond to the light coming from the entire star. CCFDI of each individual exposure are aligned by correcting their velocity table for the Keplerian motion of the star. CCFDI outside of the transits are co-added to build master-outs representative of the unocculted star, which display side lobes typical of M dwarf CCFs (Fig. 1). The lobes are thought to arise from the region around each mask line devoid of other stellar lines, while the lower continuum beyond the lobes would arise from the dense forest of lines in these cool star spectra being well-characterized by an average distance between lines. As such, the peaks of the lobes are better estimates of the actual continuum. We use the double-Gaussian model introduced by Bourrier et al. (2018b) to fit the CCFs of GJ 436, as it captures their full profile well and was shown to improve the stability of the measurements. The model consists of the sum of a Gaussian function representing the continuum and side lobes, and an inverted Gaussian function representing the core. Model parameters are the amplitude, FWHM, and RV centroid of the core Gaussian component, and the amplitude ratio, FWHM ratio, and RV centroid difference between the lobe and core Gaussian components.

All CCFDI are shifted to the rest frame of the star using the centroids derived from the fits to the master-outs. The best-fit double-Gaussian model to the master-out CCFDI in each visit is then used as a template (i.e., the amplitude ratio, FWHM ratio, and RV centroid difference between the lobe and core Gaussian components are fixed to its properties) to fit each individual exposure. We performed the same operations on CCFs reduced with the mask for M3-type stars previously used by the ESPRESSO DRS, and used the out-of-transit time-series of RV residuals, contrast, and FWHM derived from the fits to both CCF series to assess the stability of the measurements and the improvements associated with the new CCF mask (Table 2). This mask yields shallower but much narrower CCFDI, with relative precisions on the contrast and FWHM improved by 27% and 44% in Visit 1, 8% and 39% in Visit 2, respectively. The stability of the RV measurements from Visits 1 and 2 is improved by 8% and 19%, respectively. This analysis suggests that the new masks devised in Bourrier et al. (2021) improve the stability of CCF not only from K-type stellar spectra but also from M dwarf spectra. The above comparison and subsequent analyses were performed with CCFDI reduced by the DRS from the sky-subtracted products, as this improves the overall stability of their properties despite an increased dispersion on the contrast in Visit 2.

thumbnail Fig. 1

Master-out CCFDI of GJ 436 in Visit 1 (blue profile) and its bestfit double-Gaussian model (black dashed profile). The vertical dotted line indicates the measured systemic velocity.

3.2 Classical RM Analysis

Figure 2 shows the RV time-series derived from the CCFDI. In contrast to the HARPS/HARPS-N data in Bourrier et al. (2018b) the precision of the ESPRESSO data is sufficient to detect the anomaly induced by the planet, despite its small amplitude (~1 ms–1). The RM signal is consistent between the two visits within the quoted error bars, strengthening its velocimetric detection. We fitted the two visits together using the ARoME code1, which implements the prescription of Boué et al. (2013). Residuals on RV were phase-folded using the revised ephemeris (Table 1). The ARoME code does not model the Power-2 limbdarkening law used by Maxted (2022), and we therefore used the limb-darkening calculator tool2 to derive quadratic coefficients c1 = 0.440 ± 0.008 and c2 = 0.238 ± 0.012 in the wavelength range of ESPRESSO (adopting the Phoenix ACES model atmospheres; Husser et al. 2013). We also adopted a stellar macroturbulence velocity of 1 km s−1 from Fig. 11 in Wende et al. (2009) (although this parameter has little impact in our final solution), the planetary parameters reported in Table 1, and an ESPRESSO instrumental CCF width of 2.1 km s−1 (Pepe et al. 2021).

We ran a MCMC simulation with 4000 samples in order to derive the sky-projected obliquity λ and stellar rotational velocity υ sin i, together with their associated uncertainties. We set wide and uniform priors, U(0,180)${\cal U}{\left( {0,180} \right)^ \circ }$ on λ and U(0,5)km s-1${\cal U}\left( {0,5} \right){\rm{km}}\,{{\rm{s}}^{{\rm{ - 1}}}}$ on υ sin i. The best-fit RM model shown in Fig. 2 is obtained for λ=102.518.5+17.2$\lambda = 102.5{_{ - 18.5}^{ + 17.2}^ \circ }$ and υ sin i=27234+40m s-1${i_ \star } = 272_{ - 34}^{ + 40}{\rm{m}}\,{{\rm{s}}^{{\rm{ - 1}}}}$. These results are consistent within 1σ with the values reported in Bourrier et al. (2018b), although their uncertainties are reduced by up to a factor of two given the higher S/N and larger spectral coverage of the ESPRESSO data (Table 3).

Table 2

Quality assessment of ESPRESSO CCFDI.

Table 3

Comparison of RM analyses.

thumbnail Fig. 2

Top: RM anomaly induced by GJ 436 b in the RV residuals measured in Visits 1 (orange) and 2 (blue). Black dots represent a binning of the RV data with a bin size of 2.5 over the combined transits. The red line shows the classical fit to the anomaly, with associated uncertainty depicted by the gray area. Vertical dashed lines indicate the first and fourth transit contacts. Bottom: RV residuals, which display an overall dispersion of 65 cm s−1.

4 Analysis of the Planet-Occulted CCFs

RVs condense the information contained in disk-integrated CCF profiles into a single measurement, which limits our ability to detect and interpret the occultation of the stellar surface by the planet in the resulting RV anomaly. A classical analysis of the velocimetric RM effect can be further biased if the occulted stellar line profile is not well modeled or varies along the transit chord (Cegla et al. 2016b; Bourrier et al. 2017). Here we analyze the data using the Reloaded (Cegla et al. 2016a) and Revolutions (Bourrier et al. 2021) techniques, which provide direct access to the occultation signal and avoid these biases through the analysis of the planet-occulted starlight rather than the disk-integrated starlight. In both techniques, the disk-integrated CCFs are rescaled to their correct relative flux level, which is critical to extract CCFs from the planet-occulted regions. Whereas the Reloaded RM technique analyzes the stellar line from each planet-occulted region independently to retrieve and interpret the stellar surface RVs along the transit chord, the RM Revolutions technique exploits all available information by analyzing the planet-occulted CCFs together with a joint model. Hereafter, we characterize planet-occulted stellar lines according to their core RV centroid, and also according to the “true” contrast and FWHM calculated using the peaks of their model lobes as an estimate for the continuum.

4.1 Extraction and Analysis of Individual Exposures

CCFDI are scaled to a common flux level outside of the transit, and to the flux expected from the planetary disk absorption during transit using a light curve computed with the batman package (Kreidberg 2015) and the properties in Table 1. CCFs from the planet-occulted regions are retrieved by subtracting the scaled CCFDI from their corresponding master-out, and are then reset to a common flux level to yield intrinsic CCFintr that allow for a more direct comparison of the local stellar lines (Fig. 3).

Following the RM Revolutions approach (Bourrier et al. 2021), we first analyze individual exposures to evaluate the quality of each CCFintr and assess the possibility for line profile variations along the transit chord. Like the disk-integrated line, the average intrinsic stellar line displays side lobes and CCFintr were therefore fitted with a double-Gaussian profile as well. As in Bourrier et al. (2018b), we fixed the core-to-lobe FWHM ratio, amplitude ratio, and centroid difference of the CCFintr model to the values derived from the fit to the master-out CCFdi in each visit. We set uniform priors, U(1,1)km s1$U\left( { - 1,1} \right){\rm{km }}{{\rm{s}}^{ - 1}}$ on the model RV centroid (i.e., about three times the maximum stellar surface RV), U(0,10)km s1$U\left( {0,10} \right){\rm{km }}{{\rm{s}}^{ - 1}}$ on the FWHM (i.e., about three times the width of the CCFDI, which is assumed to be similar to that of the CCFintr given the slow stellar rotation), and U(0,1)$U\left( {0,1} \right)$ on the contrast. We ran a set of 100 walkers for 2000 steps, with a burn-in phase of 500 steps. The model is fitted to the data over [−50,50] km s−1 in the star rest frame.

The stellar line is clearly detected in all individual CCFintr (Fig. 3), with well-defined PDFs for their model parameters (Fig. A.1), except for the first and last exposures in both visits. Their S/N is much lower than the rest of the series, which is due to the darkened flux of the stellar limb and its partial occultation by the planet, and they were excluded from further analysis. Properties of the local stellar lines along the transit chord are shown in Fig. 4. The local line contrast, FWHM, and RVs are consistent between the two visits, as expected from the stability of the host star (Bourrier et al. 2018b). Surface RVs are all positive and remain at about the same value between 200 and 300 m s−1, which is in agreement with the results from Bourrier et al. (2018b) and the classical RM analysis (Sect. 3.2), and confirms that GJ 436b is on a polar orbit across the redshifted stellar hemisphere.

thumbnail Fig. 3

Maps of the CCFintr during the transit of GJ 436b in Visits 1 (upper panel) and 2 (lower panel). Left panels show measured profiles. The core of the stellar line from the planet-occulted regions is clearly visible as a bright streak, with side lobes visible on both sides as dark streaks. Middle panels show the RM Revolutions best-fit profiles. Right panels show residuals between measured and theoretical profiles. Transit contacts are shown as green dashed lines. Values are colored as a function of the flux, and plotted as a function of RV in the star rest frame (in abscissa) and orbital phase (in ordinate). The green solid line shows the stellar surface RV model from the RM Revolutions best fit.

4.2 Reloaded RM fit to Surface RVs

The Reloaded method is restricted to the interpretation of the RV centroids from the planet-occulted CCFs, which requires that they have a strong enough S/N to be fitted with a stellar line model. This is the case here, and we were therefore able to fit the model described in Cegla et al. (2016a) to the surface RV series. We assume solid-body rotation for the star and oversample each exposure by 5 to account for the possible blur induced by the planet motion (Bourrier et al. (2017)). Uniform priors are set on the model parameters, U(0,10) km s1$U\left( {0,10} \right)\,{\rm{km}}\,{{\rm{s}}^{ - 1}}$ for υ sin i*, and over its definition range U(180,180)°$U\left( { - 180,180} \right)^\circ $ for λ. We ran 20 walkers for 2000 steps, with a burn-in phase of 500 steps. The best fit yields vsini*=292.949.9+41.9 m s1$v\,\sin \,{i_*} = 292.9_{ - 49.9}^{ + 41.9}\,{\rm{m}}\,{{\rm{s}}^{ - 1}}$ and λ=107.519.3+23.6°$\lambda = 107.5_{ - 19.3}^{ + 23.6^\circ }$. These results are consistent within 1σ with those derived from the Reloaded analysis of HARPS/HARPS-N transits in Bourrier et al. (2018b) ((vsini*=33066+91m s1;λ=7224+33°)$\left( {v\,\sin \,{i_*} = 330_{ - 66}^{ + 91}{\rm{m}}\,{{\rm{s}}^{ - {\rm{1}}}};\lambda = 72_{ - 24}^{ + 33^\circ }} \right)$), and comparable in terms of precision with the results from the classical RM analysis (Table 3). We note that, despite their precision, the ESPRESSO data on GJ 436b do not allow us to constrain differential rotation, as a fit following the prescription of Cegla et al. (2016a) results in a uniform PDF for the relative differential rotation rate. This is in agreement with the conclusions of Browning (2008) and Morin et al. (2008) that differential rotation in fully convective dwarf stars is quenched by their strong magnetic fields.

4.3 RM Revolutions Fit to Intrinsic Line Maps

The RM Revolutions approach exploits the full information contained in the transit time-series by directly fitting a model of the stellar line to all CCFintr simultaneously (details can be found in Bourrier et al. 2021). The local stellar line is modeled as a double-Gaussian profile, whose centroids are set by the same surface RV model as in Sect. 4.3. Figure 5 shows that the master disk-integrated and intrinsic CCFs have the same profile. Comparisons of the Bayesian information criterion (BIC, Schwarz 1978; Kass & Raftery 1995; Liddle 2007) confirm that it is not justified to model the local stellar line with profiles specific to each visit (BIC = 2276) or common to both visits (BIC = 2259), compared to using the profiles describing the disk-integrated line (BIC = 2239). Properties linking the core and lobe components of the double-Gaussian CCFintr model were therefore fixed to those derived for the master CCFDI in each visit (Sect. 3.1). The line shape is therefore controlled by its contrast and FWHM, which are chosen to be constant along the transit chord but specific to each visit (allowing for center-to-limb variations in the local contrast increased the BIC of the fit by 14). The time-series of theoretical stellar lines was convolved with a Gaussian profile of width equivalent to ESPRESSO resolving power, before being fitted to the CCFintr in both visits. MCMC jump parameters are the line contrast and FWHM, λ, and υ sin i*. Uniform priors are set on all parameters, and over the same range as for the Reloaded fit for λ and υ sin i*, between U(0,10)km s1$U\left( {0,10} \right){\rm{km}}\,{{\rm{s}}^{ - 1}}$ for the local FWHM, and between U(0,1)$U\left( {0,1} \right)$ for the local contrast. We ran 60 walkers for 3000 steps, with a burn-in phase of 800 steps.

The best-fit line model accurately reproduces the local stellar lines along the transit chord, as can be seen from the residual maps in Fig. 6. These maps also highlight the stability of the star over the visits and the overall quality of the data, with no spurious features visible in the out-of-transit residuals. The surface RVs and line properties associated with the best fit are shown in Fig. 4. We derive vsini*=300.557.0+45.9m s1$v\,\sin \,{i_*} = 300.5_{57.0}^{ + 45.9}{\rm{m}}\,{{\rm{s}}^{ - 1}}$ and λ=114.117.8+22.8°$\lambda = 114.1_{ - 17.8}^{ + 22.8^\circ }$. Again these values are consistent within 1σ and comparable in terms of precision with those derived from the classical RM and Reloaded analysis (Table 3). Varying the mid-transit time and transit depth within their 3σ uncertainties and letting the semimajor axis and orbital inclination free to vary with priors set from Maxted (2022) has no impact on the derived values.

thumbnail Fig. 4

Properties of the stellar surface regions occulted by GJ 436b in ESPRESSO Visits 1 (blue) and 2 (red), derived from the fits to individual exposures. Vertical bars indicate their associated 1σ highest density intervals. Horizontal bars indicate the duration of each exposure. Solid lines are derived from the best RM Revolutions fit to the CCFintr maps (contrast and FWHM are specific to each visit and colored accordingly). The dotted line in the top panel is the best Reloaded RM fit to the surface RV series, nearly identical to the model from the RM Revolutions fit. Dashed vertical lines are the transit contacts.

thumbnail Fig. 5

Comparison between the master-out CCFDI averaged over both ESPRESSO visits (blue), and the average of all CCFintr (magenta) after they were shifted to a common rest frame using the best-fit surface RV model for GJ436.

thumbnail Fig. 6

Maps of the out-of-transit residuals, and of the in-transit residuals between CCFintr and their best-fit RM Revolutions model in ESPRESSO Visit 1 (top panel) and Visit 2 (bottom panel). Transit contacts are shown as green dashed lines.

5 Revision of the System Architecture

To refine the orbital architecture of the GJ 436 system as much as possible we performed a combined RM Revolutions fit to the archival HARPS/HARPS-N and new ESPRESSO data.

5.1 Reduction of HARPS/HARPS-N Data

One HARPS transit and two HARPS-N transits of GJ 436b were used by Bourrier et al. (2018b) to measure its obliquity using the reloaded RM method (we refer to this publication for details). Here we reduced and analyzed these datasets following the same approach as with the ESPRESSO data, albeit with some specificities, which are described below.

Older versions of the HARPS/HARPS-N DRS did not correct for the color effect in M dwarf data, which was therefore corrected for manually by Bourrier et al. (2018b). Here we used versions 3.5 and 3.7 of the HARPS and HARPS-N DRS, which apply the color correction automatically. The same template as with ESPRESSO was used. All exposures were well corrected for, except for the last one in each HARPS-N night obtained at low S/N because of high airmass and shorter exposure time. These were excluded from our analysis. CCFs were calculated with a step of 0.82 km s−1 using the same mask as the ESPRESSO data.

We observed significant correlations between the contrast of the CCFdi and the S/N. These correlations are similar between the HARPS-N visits and particulary strong in the HARPS visit (Fig. 7). The origin of these variations is unclear, and we therefore applied an empirical correction obtained by fitting a linear model between the out-of-transit contrast values and the S/N at 550 nm. This polynomial degree was set by minimizing the BIC of the fit. The correction is applied by temporarily subtracting CCFdi for their continuum flux (so that their continuum is set to zero), “stretching” them using the linear contrast model (so that their contrast is corrected to the expected value), and then adding back the continuum flux. We note that no equivalent correlations are observed in the ESPRESSO data.

5.2 Global Revolutions Fit

CCFintr were extracted in the same way as for the ESPRESSO data. GJ 436 is occulted by the planet in 12 exposures in the HARPS data, and in 9 exposures in the HARPS-N data. Applying the same analysis as in Sect. 4.1 (see Figs. A.2 and A.3), we exclude exposures that display uniform RV PDFs and/or contrast and FWHM PDFs that are consistent with null values, which is an indication that the S/N of their CCFintr is too low to detect the stellar line. This concerns exposures at the limbs of the star: the three first and three last ones in the HARPS transit, and the first and last ones in the HARPS-N transits.

We used the RM Revolutions technique to fit the five visits together, using the same procedure and settings as in Sect. 4.3. The double-Gaussian profile of the local line model is set to the template derived for the disk-integrated line in each visit, with the contrast and FWHM of the line let free to vary in each visit. Theoretical line models are convolved with the instrumental response of each instrument before they are compared with the data. Figure 8 shows the surface RVs derived from the individual analysis of each exposure (as in Sect. 4.1), along with the best-fit model derived from the best RM Revolutions fit. We derive vsini*=293.552.2+43.7 m s1$v\,\sin \,{i_*} = 293.5_{ - 52.2}^{ + 43.7}\,{\rm{m}}\,{{\rm{s}}^{ - 1}}$ and λ=113.517.3+23.3°$\lambda = 113.5_{ - 17.3}^{ + 23.3^\circ }$. The correction applied to the HARPS/HARPS-N CCFDI improved the stability of the CCFintr contrast, but has negligible impact on these results.

We note that combining the HARPS/HARPS-N and ESPRESSO data does not substantially improve the precision on the derived properties (Table 3). To understand why, we performed a RM Revolutions fit on the HARPS/HARPS-N data alone, which yielded vsini*=290148+128 m s1$v\,\sin \,{i_*} = 290_{ - 148}^{ + 128}\,{\rm{m}}\,{{\rm{s}}^{ - 1}}$ and λ=10250+57°$\lambda = 102_{ - 50}^{ + 57^\circ }$. Those results remain consistent within 1σ with those obtained by Bourrier et al. (2018b) but υ sin i* is shifted to lower values while λ is shifted to larger values, with both properties having larger uncertainties. First, this explains why the ESPRESSO dataset dominates the global fit of the five visits. Also, these differences with Bourrier et al. (2018b) could be due to biases in the reloaded RM technique these authors applied to the HARPS/HARPS-N datasets. Indeed, Fig. 8 shows that the surface RVs have asymmetrical uncertainty ranges resulting from non-Gaussian PDFs from the MCMC fits to the HARPS/HARPS-N CCFintr. In such a case, the least-square minimization used to fit CCFintr in the traditional reloaded RM approach yields biased estimates for surface RVs and their uncertainties, which then bias the properties derived by fitting those RV series. This highlights the utility of the RM Revolutions technique in fitting the series of CCFintr directly using a M.M. approach rather than fitting their properties in a two-stage approach as with the reloaded RM technique.

Our measurement of υ sin i* is the most precise to date of all those obtained by analysis of the RM effect, with only five other exoplanet host stars having measurements more precise than 50ms−1 (based on the NASA exoplanet archive). Following the approach of Bourrier et al. (2018b), we combined our PDF for υ sin i* with Gaussian PDFs for the stellar radius (0.438 ± 0.013 R, Maxted 2022) and rotation period (44.09 ± 0.08 d, Bourrier et al. 2018b) to build a PDF for the stellar inclination, revising its value to i*=35.77.6+5.9°${i_*} = 35.7_{ - 7.6}^{ + 5.9^\circ }$ or 144.25.9+7.6°$144.2_{ - 5.9}^{ + 7.6^\circ }$. The PDFs for the true 3D obliquity of GJ 436b associated with these two degenerate configurations have similar widths and medians3 and we averaged them to derive ψ=103.211.512.8°${\rm{\psi = 103}}{\rm{.2}}_{ - 11.5}^{12.8^\circ }$ (Table 1). The corresponding architecture for the GJ 436 system is shown in Fig. 9.

thumbnail Fig. 7

Contrast of the HARPS (magenta), HARPS-N (green and orange), and ESPRESSO (blue and red) CCFDI. HARPS/HARPS-N data are shown before (top panels) and after (bottom panels) correcting for their correlation with the S/N at 550 nm. Solid lines show best-fit models to the out-of-transit data (squares; in-transit points are shown as disks). Error bars are not shown for the sake of clarity.

thumbnail Fig. 8

Stellar surface RVs of the regions occulted by GJ 436b in ESPRESSO visits 1 (blue) and 2 (red), HARPS Visit A (magenta), and HARPS-N visits B (orange) and C (green) derived from the fits to individual exposures (same details as Fig. 4). The solid black curve is the best-fit model derived from the joint Revolutions fit to all visits.

thumbnail Fig. 9

Projection of GJ 436 in the plane of the sky for the best-fit orbital architecture, in the configuration where the stellar spin (shown as a black arrow extending from the north pole) is pointing toward the Earth. The south pole is shown as en empty disk. The stellar equator is represented as a solid (dashed) black line on the visible (hidden) stellar hemisphere. The stellar disk is colored as a function of its surface RV field. The normal to the orbital plane of GJ 436b is shown as a green arrow. The thick green solid curve represents the best-fit orbital trajectory. The thin lines surrounding it show orbits obtained for orbital inclination, semi-major axis, and sky-projected obliquity values drawn randomly within 1σ from their probability distributions. The star, planet (black disk), and orbit are to scale.

6 Conclusion

The classical, Reloaded RM, and RM Revolutions analysis of two transits of GJ436b observed with ESPRESSO yield similar results for the sky-projected obliquity λ and stellar rotational velocity υ sin i*, which are also consistent with the values obtained by Bourrier et al. (2018b) (Table 3). This unambiguously confirms the conclusions of these latter authors that GJ 436b is on a polar orbit around its slow-rotating host.

The three analysis techniques yield unbiased results on λ and υ sin i* with comparable precision when the occulted stellar lines are detectable in each exposure, are constant across the transit chord, and solid-body rotation dominates their Doppler shift, as is the case for GJ 436b.

The smaller uncertainties provided by the classical RM fit of GJ 436b likely arise from the underlying assumption of a constant stellar line profile, which removes freedom to the fit compared to the other techniques. However, this assumption is only validated by the independent analysis of the local line profiles in the Reloaded RM and RM Revolutions techniques. In a general case, the classical RM approach should be used with caution because it does not allow the occulted stellar lines to be extracted and studied, making its results prone to bias if these lines vary.

The comparable precision between the Reloaded RM and RM Revolutions fits of GJ 436b likely arises from the detectability of the occulted stellar lines in nearly all exposures. Indeed, the reloaded RM technique relies on the possibility to derive the RVs of planet-occulted regions by fitting the occulted stellar line in individual exposures. This technique thus loses precision as planets get smaller/host stars fainter and fewer lines can be detected, up to the point where it fails entirely. This was the case with the super-Earth HD 3167b, whose signal could only be detected through the RM Revolutions technique (Bourrier et al. 2021).

While the classical, Reloaded RM, and RM Revolutions techniques yield comparable results in ideal cases, the latter is required to fully exploit a RM transit dataset, and assess and avoid the presence of possible biases due to line shape variations, and we encourage its use to measure obliquities from all types of planetary systems.

Acknowledgements

We are grateful for a constructive review from the referee, which led us to go deeper into our analyses. The authors acknowledge the ESPRESSO project team for its effort and dedication in building the ESPRESSO instrument. This work has been carried out in the frame of the National Centre for Competence in Research PlanetS supported by the Swiss National Science Foundation (SNSF). The authors acknowledge the financial support of the SNSF. This project has received funding from the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation programme (project Spice Dune, grant agreement No 947634; project Four Aces; grant agreement No 724427). This work was supported by FCT – Fundação para a Ciências e a Tecnologia through national funds and by FEDER through COMPETE2020 – Programa Operacional Competitividade e Internacionalização by these grants: UID/FIS/04434/2019; UIDB/04434/2020; UIDP/04434/2020; PTDC/FIS-AST/32113/2017 & POCI-01-0145-FEDER-032113; PTDC/FIS-AST/28953/2017 & POCI-01-0145-FEDER-028953; PTDC/FIS-AST/28987/2017 & POCI-01-0145-FEDER-028987. M.R.Z.O. acknowledges funding from the Spanish Ministery of Science and Innovation through project PID2019-109522GB-C51. V.A. acknowledges the support from FCT through Investigador FCT contracts nr. IF/00650/2015/CP1273/CT0001. NJN acknowledges support from the following projects: UIDB/04434/2020 & UIDP/04434/2020, CERN/FIS-PAR/0037/2019, PTDC/FIS-OUT/29048/2017, COMPETE2020: POCI-01-0145-FEDER-028987 & FCT: PTDC/FISAST/28987/2017. ASM, JIGH, CAP, and RR acknowledge financial support from the Spanish Ministry of Science and Innovation (MICINN) project AYA2017-86389-P. JIGH also acknowledges financial support from the Spanish MICINN under 2013 Ramón y Cajal program RYC-2013-14875. R.A. is a Trottier Postdoctoral Fellow and acknowledges support from the Trottier Family Foundation. This work was supported in part through a grant from FRQNT. M.R.Z.O. acknoledges financial support from project PID2019-109522GB-C51 from the Spanish Ministerio de Ciencia e Innovación. DE acknowledges financial support from the Swiss National Science Foundation for project 200021_200726. S.G.S. acknowledges the support from FCT through Investigador FCT contract nr. CEECIND/00826/2018 and POPH/FSE (EC).

Appendix A PDFs of the Fits to Individual Exposures

thumbnail Fig. A.1

PDFs of the contrast (upper panels), FWHM (middle panels), and RVs (lower panels) of the Double-Gaussian line model fitted to individual CCFintr in ESPRESSO Visits 1 (left) and 2 (right). Blue lines indicate the PDF median values, with dashed green lines showing the 1σ highest density intervals (multiple in case of multimodal PDFs). Intransit exposure indexes are shown in each subplot.

thumbnail Fig. A.2

PDFs of the CCFintr properties for the HARPS transit (same details as Fig. A.1).

thumbnail Fig. A.3

PDFs of the CCFintr properties for the first (left) and second (right) HARPS-N transit (same details as Fig. A.1).

References

  1. Albrecht, S.H., Marcussen, M.L., Winn, J.N., Dawson, R.I., & Knudstrup, E. 2021, ApJ, 916, L1 [NASA ADS] [CrossRef] [Google Scholar]
  2. Allart, R., Bourrier, V., Lovis, C., et al. 2019, A&A, 623, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Anderson, K.R., Storch, N.I., & Lai, D. 2016, MNRAS, 456, 3671 [NASA ADS] [CrossRef] [Google Scholar]
  4. Attia, M., Bourrier, V., Eggenberger, P., et al. 2021, A&A, 647, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Baranne, A., Queloz, D., Mayor, M., et al. 1996, A&AS, 119, 373 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Ben-Jaffel, L., Ballester, G.E.G. Munoz, A.G., et al. 2022, Nat. Astron., 6, 141 [NASA ADS] [CrossRef] [Google Scholar]
  7. Beust, H., Bonfils, X., Montagnier, G., Delfosse, X., & Forveille, T. 2012, A&A, 545, A88 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Boué, G., Montalto, M., Boisse, I., Oshagh, M., & Santos, N.C. 2013, A&A, 550, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  9. Bourrier, V., & Hébrard, G. 2014, A&A, 569, A65 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Bourrier, V., Cegla, H.M., Lovis, C., & Wyttenbach, A. 2017, A&A, 599, A33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  11. Bourrier, V., Lecavelier des Etangs, A., Ehrenreich, D., et al. 2018a, A&A, 620, A147 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Bourrier, V., Lovis, C., Beust, H., et al. 2018b, Nature, 553, 477 [NASA ADS] [CrossRef] [Google Scholar]
  13. Bourrier, V., Lovis, C., Cretignier, M., et al. 2021, A&A, 654, A152 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Browning, M.K. 2008, ApJ, 676, 1262 [NASA ADS] [CrossRef] [Google Scholar]
  15. Casasayas-Barris, N., Palle, E., Stangret, M., et al. 2021, A&A, 647, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Cegla, H.M., Lovis, C., Bourrier, V., et al. 2016a, A&A, 588, A127 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  17. Cegla, H.M., Oshagh, M., Watson, C.A., et al. 2016b, ApJ, 819, 67 [NASA ADS] [CrossRef] [Google Scholar]
  18. Correia, A.C.M., Bourrier, V., & Delisle, J.B. 2020, A&A, 635, A37 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. dos Santos, L.A., Ehrenreich, D., Bourrier, V., et al. 2019, A&A, 629, A47 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Ehrenreich, D., Bourrier, V., Wheatley, P.J., et al. 2015, Nature, 522, 459 [NASA ADS] [CrossRef] [Google Scholar]
  21. Fabrycky, D., & Tremaine, S. 2007, ApJ, 669, 1298 [NASA ADS] [CrossRef] [Google Scholar]
  22. Foreman-Mackey, D., Hogg, D.W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [NASA ADS] [CrossRef] [Google Scholar]
  23. Husser, T.O., Wende-von Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Kass, R.E., & Raftery, A.E. 1995, J. Am. Stat. Assoc., 90, 773 [Google Scholar]
  25. Kirk, J., Alam, M.K., Löpez-Morales, M., & Zeng, L. 2020, AJ, 159, 115 [NASA ADS] [CrossRef] [Google Scholar]
  26. Knutson, H.A., Dragomir, D., Kreidberg, L., et al. 2014, ApJ, 794, 155 [NASA ADS] [CrossRef] [Google Scholar]
  27. Kreidberg, L. 2015, PASP, 127, 1161 [Google Scholar]
  28. Kulow, J.R., France, K., Linsky, J., & Loyd, R.O.P. 2014, ApJ, 786, 132 [NASA ADS] [CrossRef] [Google Scholar]
  29. Lanotte, A.A., Gillon, M., Demory, B.-O., et al. 2014, A&A, 572, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  30. Lavie, B., Ehrenreich, D., Bourrier, V., et al. 2017, A&A, 605, L7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Liddle, A.R. 2007, MNRAS, 377, L74 [NASA ADS] [CrossRef] [Google Scholar]
  32. Lothringer, J.D., Barman, T., & Koskinen, T. 2018, ApJ, 866, 27 [NASA ADS] [CrossRef] [Google Scholar]
  33. Mardling, R.A. 2008, ArXiv e-prints [arXiv:0805.1928] [Google Scholar]
  34. Maxted, P.F.L. 2018, A&A, 616, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Maxted, P.F.L., Ehrenreich, D., Wilson, T.G., et al. 2022, MNRAS, 514, 77 [NASA ADS] [CrossRef] [Google Scholar]
  36. Morello, G., Tsiaras, A., Howarth, I.D., & Homeier, D. 2017, AJ, 154, 111 [NASA ADS] [CrossRef] [Google Scholar]
  37. Morin, J., Donati, J.F., Petit, P., et al. 2008, MNRAS, 390, 567 [NASA ADS] [CrossRef] [Google Scholar]
  38. Naoz, S. 2016, ARA&A, 54, 441 [Google Scholar]
  39. Ninan, J.P., Stefansson, G., Mahadevan, S., et al. 2020, ApJ, 894, 97 [NASA ADS] [CrossRef] [Google Scholar]
  40. Owen, J.E., & Lai, D. 2018, MNRAS, 479, 5012 [NASA ADS] [CrossRef] [Google Scholar]
  41. Palle, E., Nortmann, L., Casasayas-Barris, N., et al. 2020, A&A, 638, A61 [EDP Sciences] [Google Scholar]
  42. Pepe, F., Mayor, M., Galland, F., et al. 2002, A&A, 388, 632 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Pepe, F., Cristiani, S., Rebolo, R., et al. 2021, A&A, 645, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  44. Petrovich, C., Munoz, D.J., Kratter, K.M., & Malhotra, R. 2020, ApJ, 902, L5 [NASA ADS] [CrossRef] [Google Scholar]
  45. Piaulet, C., Benneke, B., Rubenzahl, R.A., et al. 2021, AJ, 161, 70 [NASA ADS] [CrossRef] [Google Scholar]
  46. Schwarz, G. 1978, Ann. Statist., 6, 461 [Google Scholar]
  47. Spake, J.J., Sing, D.K., Evans, T.M., et al. 2018, Nature, 557, 68 [NASA ADS] [CrossRef] [Google Scholar]
  48. Spake, J.J., Oklopcic, A., & Hillenbrand, L.A. 2021, AJ, 162, 284 [NASA ADS] [CrossRef] [Google Scholar]
  49. Stefansson, G., Mahadevan, S., Petrovich, C., et al. 2022, ApJ, 931, L15 [NASA ADS] [CrossRef] [Google Scholar]
  50. Triaud, A.H.M.J. 2018, The Rossiter-McLaughlin Effect in Exoplanet Research, 2 [Google Scholar]
  51. Trifonov, T., Kürster, M., Zechmeister, M., et al. 2018, A&A, 609, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  52. Wehbe, B., Cabral, A., Martins, J.H.C., et al. 2020, MNRAS, 491, 3515 [CrossRef] [Google Scholar]
  53. Wende, S., Reiners, A., & Ludwig, H.G. 2009, A&A, 508, 1429 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Yee, S.W., Petigura, E.A., Fulton, B.J., et al. 2018, AJ, 155, 255 [NASA ADS] [CrossRef] [Google Scholar]

2

https://exoctk.stsci.edu/limb_darkening. This tool is part of the Exoplanet Characterization toolkit, https://exoctk. readthedocs.io/en/latest/

3

For transiting planets, the 3D obliquity Ψ = arccos(sin(i) cos(λ) sin(ip) + cos(i) cos(ip)) approximates as arccos(sin(i) cos(λ)).

All Tables

Table 1

Properties of the GJ 436 system.

Table 2

Quality assessment of ESPRESSO CCFDI.

Table 3

Comparison of RM analyses.

All Figures

thumbnail Fig. 1

Master-out CCFDI of GJ 436 in Visit 1 (blue profile) and its bestfit double-Gaussian model (black dashed profile). The vertical dotted line indicates the measured systemic velocity.

In the text
thumbnail Fig. 2

Top: RM anomaly induced by GJ 436 b in the RV residuals measured in Visits 1 (orange) and 2 (blue). Black dots represent a binning of the RV data with a bin size of 2.5 over the combined transits. The red line shows the classical fit to the anomaly, with associated uncertainty depicted by the gray area. Vertical dashed lines indicate the first and fourth transit contacts. Bottom: RV residuals, which display an overall dispersion of 65 cm s−1.

In the text
thumbnail Fig. 3

Maps of the CCFintr during the transit of GJ 436b in Visits 1 (upper panel) and 2 (lower panel). Left panels show measured profiles. The core of the stellar line from the planet-occulted regions is clearly visible as a bright streak, with side lobes visible on both sides as dark streaks. Middle panels show the RM Revolutions best-fit profiles. Right panels show residuals between measured and theoretical profiles. Transit contacts are shown as green dashed lines. Values are colored as a function of the flux, and plotted as a function of RV in the star rest frame (in abscissa) and orbital phase (in ordinate). The green solid line shows the stellar surface RV model from the RM Revolutions best fit.

In the text
thumbnail Fig. 4

Properties of the stellar surface regions occulted by GJ 436b in ESPRESSO Visits 1 (blue) and 2 (red), derived from the fits to individual exposures. Vertical bars indicate their associated 1σ highest density intervals. Horizontal bars indicate the duration of each exposure. Solid lines are derived from the best RM Revolutions fit to the CCFintr maps (contrast and FWHM are specific to each visit and colored accordingly). The dotted line in the top panel is the best Reloaded RM fit to the surface RV series, nearly identical to the model from the RM Revolutions fit. Dashed vertical lines are the transit contacts.

In the text
thumbnail Fig. 5

Comparison between the master-out CCFDI averaged over both ESPRESSO visits (blue), and the average of all CCFintr (magenta) after they were shifted to a common rest frame using the best-fit surface RV model for GJ436.

In the text
thumbnail Fig. 6

Maps of the out-of-transit residuals, and of the in-transit residuals between CCFintr and their best-fit RM Revolutions model in ESPRESSO Visit 1 (top panel) and Visit 2 (bottom panel). Transit contacts are shown as green dashed lines.

In the text
thumbnail Fig. 7

Contrast of the HARPS (magenta), HARPS-N (green and orange), and ESPRESSO (blue and red) CCFDI. HARPS/HARPS-N data are shown before (top panels) and after (bottom panels) correcting for their correlation with the S/N at 550 nm. Solid lines show best-fit models to the out-of-transit data (squares; in-transit points are shown as disks). Error bars are not shown for the sake of clarity.

In the text
thumbnail Fig. 8

Stellar surface RVs of the regions occulted by GJ 436b in ESPRESSO visits 1 (blue) and 2 (red), HARPS Visit A (magenta), and HARPS-N visits B (orange) and C (green) derived from the fits to individual exposures (same details as Fig. 4). The solid black curve is the best-fit model derived from the joint Revolutions fit to all visits.

In the text
thumbnail Fig. 9

Projection of GJ 436 in the plane of the sky for the best-fit orbital architecture, in the configuration where the stellar spin (shown as a black arrow extending from the north pole) is pointing toward the Earth. The south pole is shown as en empty disk. The stellar equator is represented as a solid (dashed) black line on the visible (hidden) stellar hemisphere. The stellar disk is colored as a function of its surface RV field. The normal to the orbital plane of GJ 436b is shown as a green arrow. The thick green solid curve represents the best-fit orbital trajectory. The thin lines surrounding it show orbits obtained for orbital inclination, semi-major axis, and sky-projected obliquity values drawn randomly within 1σ from their probability distributions. The star, planet (black disk), and orbit are to scale.

In the text
thumbnail Fig. A.1

PDFs of the contrast (upper panels), FWHM (middle panels), and RVs (lower panels) of the Double-Gaussian line model fitted to individual CCFintr in ESPRESSO Visits 1 (left) and 2 (right). Blue lines indicate the PDF median values, with dashed green lines showing the 1σ highest density intervals (multiple in case of multimodal PDFs). Intransit exposure indexes are shown in each subplot.

In the text
thumbnail Fig. A.2

PDFs of the CCFintr properties for the HARPS transit (same details as Fig. A.1).

In the text
thumbnail Fig. A.3

PDFs of the CCFintr properties for the first (left) and second (right) HARPS-N transit (same details as Fig. A.1).

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.