Context.
Stellar activity and instrumental signals are the main limitations to the detection of Earth-like planets using the radial-velocity (RV) technique. Recent studies show that the key to ...mitigating those perturbing effects might reside in analysing the spectra themselves, rather than the RV time series and a few activity proxies.
Aims.
The goal of this paper is to demonstrate that we can reach further improvement in RV precision by performing a principal component analysis (PCA) decomposition of the shell time series, with the shell as the projection of a spectrum onto the space-normalised flux versus flux gradient.
Methods.
By performing a PCA decomposition of shell time series, it is possible to obtain a basis of first-order spectral variations that are not related to Keplerian motion. The time coefficients associated with this basis can then be used to correct for non-Dopplerian signatures in RVs.
Results.
We applied this new method on the YARARA post-processed spectra time series of HD 10700 (
τ
Ceti) and HD 128621 (
α
Cen B). On HD 10700, we demonstrate, thanks to planetary signal injections, that this new approach can successfully disentangle real Dopplerian signals from instrumental systematics. The application of this new methodology on HD 128621 shows that the strong stellar activity signal seen at the stellar rotational period and one-year aliases becomes insignificant in a periodogram analysis. The RV root mean square on the 5-yr data is reduced from 2.44 m s
−1
down to 1.73 m s
−1
. This new approach allows us to strongly mitigate stellar activity, however, noise injections tests indicate that rather high signal-to-noise ratio (S/N > 250) is required to correct for the observed activity signal on HD 128621.
Aims.
Even the most precise radial-velocity instruments gather high-resolution spectra that present systematic errors that a data reduction pipeline cannot identify and correct for efficiently by ...simply analysing a set of calibrations and a single science frame. In this paper we aim at improving the radial-velocity precision of HARPS measurements by ‘cleaning’ individual extracted spectra using the wealth of information contained in spectral time series.
Methods.
We developed YARARA, a post-processing pipeline designed to clean high-resolution spectra of instrumental systematics and atmospheric contamination. Spectra are corrected for: tellurics, interference patterns, detector stitching, ghosts, and fibre B contaminations, as well as more advanced spectral line-by-line corrections. YARARA uses principal component analysis on spectral time series with prior information to disentangle contaminations from real Doppler shifts. We applied YARARA to three systems, HD 10700, HD 215152, and HD 10180, and compared our results to the standard HARPS data reduction software and the SERVAL post-processing pipeline.
Results.
We ran YARARA on the radial-velocity dataset of three stars intensively observed with HARPS: HD 10700, HD 215152, and HD 10180. For HD 10700, we show that YARARA enables us to obtain radial-velocity measurements that present an rms smaller than 1 m s
−1
over the 13 years of the HARPS observations, which is 20% and 10% better than the HARPS data reduction software and the SERVAL post-processing pipeline, respectively. We also injected simulated planets into the data of HD 10700 and demonstrated that YARARA does not alter pure Doppler-shifted signals. For HD 215152, we demonstrated that the 1-year signal visible in the periodogram becomes marginal after processing with YARARA and that the signals of the known planets become more significant. Finally, for HD 10180, the six known exoplanets are well recovered, although different orbital parameters and planetary masses are provided by the new reduced spectra.
Conclusions.
The post-processing correction of spectra using spectral time series allows the radial-velocity precision of HARPS data to be significantly improved and demonstrates that for the extremely quiet star HD 10700 a radial-velocity rms better than 1 m s
−1
can be reached over the 13 years of HARPS observations. Since the processing proposed in this paper does not absorb planetary signals, its application to intensively followed systems is promising and will certainly result in advances in the detections of the lightest exoplanets.
ABSTRACT
The detection of Earth-like planets with the radial-velocity (RV) method is currently limited by the presence of stellar activity signatures. On rotational time-scales, spots and plages (or ...faculae) are known to introduce different RV signals, but their corrections require better activity proxies. The best-known chromospheric activity proxies in the visible are the Ca ii H&K lines, but the physical quantities measured by their profiles need to be clarified. We first investigate resolved images of the Sun in order to better understand the spectrum of plages, spots, and the network using the Meudon spectroheliogram. We show that distinct line profiles are produced by plages, spots, and by the network component and we also derived the centre-to-limb variations of the three profiles. Some care is required to disentangle their contributions due to their similarities. By combining disc-integrated spectra from the ISS high-resolution spectrograph with SDO direct images of the Sun, we managed to extract a high-resolution emission spectrum of the different components, which tend to confirm the spectra extracted from the Meudon spectroheliogram datacubes. Similar results were obtained with the HARPS-N Sun-as-a-star spectra. We concluded using a three-component model that the temporal variation of the popular Sindex contains, on average for the 24th solar cycle: $70\pm 12{{\ \rm per\ cent}}$ of plage, $26\pm 12{{\ \rm per\ cent}}$ of network, and $4\pm 4 {{\ \rm per\ cent}}$ of spots. This preliminary investigation suggests that a detailed study of the Ca ii H&K profiles may provide rich information about the filling factor and distribution of different types of active regions.
Context.
The detection of Earth-like planets with the radial velocity (RV) method is extremely challenging today due to the presence of non-Doppler signatures such as stellar activity and ...instrumental signals that mimic and hide the signals of exoplanets. In a previous paper, we presented the YARARA pipeline, which implements corrections for telluric absorption, stellar activity, and instrumental systematics at the spectral level, and then it extracts line-by-line (LBL) RVs with a significantly better precision than standard pipelines.
Aims.
In this paper, we demonstrate that further gains in RV precision can be achieved by performing principal component analysis (PCA) decomposition on the LBL RVs.
Methods.
The mean-insensitive nature of PCA means that it is unaffected by true Doppler shifts, and thus can be used to isolate and correct nuisance signals other than planets.
Results.
We analysed the data of 20 intensively observed HARPS targets by applying our PCA approach on the LBL RVs obtained by YARARA. The first principal components show similarities across most of the stars and correspond to newly identified instrumental systematics for which we can now correct. For several targets, this results in an unprecedented RV root-mean-square of around 90 cm s
−1
over the full lifetime of HARPS. We used the corrected RVs to confirm a previously published 120-day signal around 61 Vir, and to detect a super-Earth candidate (
K
~ 60 ± 6 cm s
−1
,
m
sin
i
= 6.6 ± 0.7
M
⊕
) around the G6V star HD 20794, which spends part of its 600-day orbit within the habitable zone of the host star.
Conclusions.
This study highlights the potential of LBL PCA to identify and correct hitherto unknown, long-term instrumental effects and thereby extend the sensitivity of existing and future instruments towards the Earth analogue regime.
We present the analysis of TESS optical photometry of WASP-121b, which reveals the phase curve of this transiting ultra-hot Jupiter. Its hotspot is located at the sub-stellar point, showing ...inefficient heat transport from the dayside (2870 ± 50 K) to the nightside (<2500 K at 3
σ
) at the altitudes probed by TESS. The TESS eclipse depth, measured at the shortest wavelength to date for WASP-121b, confirms the strong deviation from blackbody planetary emission. Our atmospheric retrieval on the complete emission spectrum supports the presence of a temperature inversion, which can be explained by the presence of VO and possibly TiO and FeH. The strong planetary emission at short wavelengths could arise from an H
−
continuum.
Comparisons of the alignment of exoplanets with a common host star and each other can be used to distinguish among concurrent evolution scenarios for the star and the planets. However, multi-planet ...systems usually host mini-Neptunes and super-Earths, whose sizes make orbital architecture measurements challenging. We introduce the Rossiter-McLaughlin effect Revolutions (RMR) technique, which can access the spin-orbit angle of small exoplanets by exploiting the full extent of information contained in spectral transit time series. We validated the technique through its application to published HARPS-N data of the mini-Neptune HD 3167c (
P
= 29.8 days), refining its high sky-projected spin-orbit angle (−108.9
−5.5
+5.4°
), and we applied it to new ESPRESSO observations of the super-Earth HD 3167 b (
P
= 0.96 days), revealing an aligned orbit (−6.6
−7.9
+6.6°
). Surprisingly different variations in the contrast of the stellar lines occulted by the two planets can be reconciled by assuming a latitudinal dependence of the stellar line shape. In this scenario, a joint fit to both datasets constrains the inclination of the star (111.6
−3.3
+3.1°
) and the 3D spin-orbit angles of HD 3167b (29.5
−9.4
+7.2°
) and HD 3167c (107.7
−4.9
+5.1°
). The projected spin-orbit angles do not depend on the model for the line contrast variations, and so, with a mutual inclination of 102.3
−8.0
+7.4°
, we can conclude that the two planets are on perpendicular orbits. This could be explained by HD 3167b being strongly coupled to the star and retaining its primordial alignment, whereas HD 3167c would have been brought to a nearly polar orbit via secular gravitational interactions with an outer companion. Follow-up observations of the system and simulations of its dynamical evolution are required to search for this companion and explore the likelihood of this scenario. HD 3167 b (
R
= 1.7
R
Earth
) is the smallest exoplanet with a confirmed spectroscopic Rossiter-McLaughlin signal. The RMR technique opens the way to determining the orbital architectures of the super-Earth and Earth-sized planet populations.
Context.
The solar telescope connected to HARPS-N has been observing the Sun since the summer of 2015. Such a high-cadence, long-baseline data set is crucial for understanding spurious ...radial-velocity signals induced by our Sun and by the instrument. On the instrumental side, this data set allowed us to detect sub- m s
−1
systematics that needed to be corrected for.
Aims.
The goals of this manuscript are to (i) present a new data reduction software for HARPS-N, (ii) demonstrate the improvement brought by this new software during the first three years of the HARPS-N solar data set, and (iii) release all the obtained solar products, from extracted spectra to precise radial velocities.
Methods.
To correct for the instrumental systematics observed in the data reduced with the current version of the HARPS-N data reduction software (DRS version 3.7), we adapted the newly available ESPRESSO DRS (version 2.2.3) to HARPS-N and developed new optimised recipes for the spectrograph. We then compared the first three years of HARPS-N solar data reduced with the current and new DRS.
Results.
The most significant improvement brought by the new DRS is a strong decrease in the day-to-day radial-velocity scatter, from 1.27 to 1.07 m s
−1
; this is thanks to a more robust method to derive wavelength solutions, but also to the use of calibrations closer in time. The newly derived solar radial-velocities are also better correlated with the chromospheric activity level of the Sun in the long term, with a Pearson correlation coefficient of 0.93 compared to 0.77 before, which is expected from our understanding of stellar signals. Finally, we also discuss how HARPS-N spectral ghosts contaminate the measurement of the calcium activity index, and we present an efficient technique to derive an index free of instrumental systematics.
Conclusions.
This paper presents a new data reduction software for HARPS-N and demonstrates its improvements, mainly in terms of radial-velocity precision, when applied to the first three years of the HARPS-N solar data set. Those newly reduced solar data, representing an unprecedented time series of 34 550 high-resolution spectra and precise radial velocities, are released alongside this paper. Those data are crucial to understand stellar activity signals in solar-type stars further and develop the mitigating techniques that will allow us to detect other Earths.
ABSTRACT
Stellar magnetic activity produces time-varying distortions in the photospheric line profiles of solar-type stars. These lead to systematic errors in high-precision radial-velocity ...measurements, which limit efforts to discover and measure the masses of low-mass exoplanets with orbital periods of more than a few tens of days. We present a new data-driven method for separating Doppler shifts of dynamical origin from apparent velocity variations arising from variability-induced changes in the stellar spectrum. We show that the autocorrelation function (ACF) of the cross-correlation function used to measure radial velocities is effectively invariant to translation. By projecting the radial velocities on to a subspace labelled by the observation identifiers and spanned by the amplitude coefficients of the ACF’s principal components, we can isolate and subtract velocity perturbations caused by stellar magnetic activity. We test the method on a 5-yr time sequence of 853 daily 15-min observations of the solar spectrum from the HARPS-N instrument and solar-telescope feed on the 3.58-m Telescopio Nazionale Galileo. After removal of the activity signals, the heliocentric solar velocity residuals are found to be Gaussian and nearly uncorrelated. We inject synthetic low-mass planet signals with amplitude K = 40 cm s−1 into the solar observations at a wide range of orbital periods. Projection into the orthogonal complement of the ACF subspace isolates these signals effectively from solar activity signals. Their semi-amplitudes are recovered with a precision of ∼ 6.6 cm s−1, opening the door to Doppler detection and characterization of terrestrial-mass planets around well-observed, bright main-sequence stars across a wide range of orbital periods.
Context.
The study of exoplanet atmospheres is essential for understanding the formation, evolution, and composition of exoplanets. The transmission spectroscopy technique is playing a significant ...role in this domain. In particular, the combination of state-of-the-art spectrographs at low- and high-spectral resolution is key to our understanding of atmospheric structure and composition.
Aims.
We observed two transits of the close-in sub-Saturn-mass planet, WASP-127b, with ESPRESSO in the frame of the Guaranteed Time Observations Consortium. We aim to use these transit observations to study the system architecture and the exoplanet atmosphere simultaneously.
Methods.
We used the Reloaded Rossiter-McLaughlin technique to measure the projected obliquity
λ
and the projected rotational velocity
v
eq
⋅sin(
i
*
). We extracted the high-resolution transmission spectrum of the planet to study atomic lines. We also proposed a new cross-correlation framework to search for molecular species and we applied it to water vapor.
Results.
The planet is orbiting its slowly rotating host star (
v
eq
⋅sin(
i
*
) = 0.53
−0.05
+0.07
km s
−1
) on a retrograde misaligned orbit (
λ
= −128.41
−5.46
+5.60
°). We detected the sodium line core at the 9-
σ
confidence level with an excess absorption of 0.34 ± 0.04%, a blueshift of 2.74 ± 0.79 km s
−1
, and a full width at half maximum of 15.18 ± 1.75 km s
−1
. However, we did not detect the presence of other atomic species but set upper limits of only a few scale heights. Finally, we put a 3-
σ
upper limit on the average depth of the 1600 strongest water lines at equilibrium temperature in the visible band of 38 ppm. This constrains the cloud-deck pressure between 0.3 and 0.5 mbar by combining our data with low-resolution data in the near-infrared and models computed for this planet.
Conclusions.
WASP-127b, with an age of about 10 Gyr, is an unexpected exoplanet by its orbital architecture but also by the small extension of its sodium atmosphere (~7 scale heights). ESPRESSO allows us to take a step forward in the detection of weak signals, thus bringing strong constraints on the presence of clouds in exoplanet atmospheres. The framework proposed in this work can be applied to search for molecular species and study cloud-decks in other exoplanets.
ABSTRACT We present the discovery of three super-Earth candidates orbiting HD 48948, a bright K-dwarf star with an apparent magnitude of $m_V$ = 8.58 mag. As part of the HARPS-N Rocky Planet Search ...programme, we collect 189 high-precision radial velocity measurements using the HARPS-N spectrograph from 2013 October 6, to 2023 April 16. Various methodologies are applied to extract the radial velocities from the spectra, and we conduct a comprehensive comparative analysis of the outcomes obtained through these diverse extraction techniques. To ensure the robustness of our findings, we employ several methods to address stellar variability, with a focus on Gaussian Process regression. To account for the impact of stellar variability and correlated noise in the radial velocity data set, we include activity indicators, such as $\log R^{^{\prime }}_{\mathrm{HK}}$ and bisector span, in the multidimensional Gaussian Process regression. Our analysis reveals three planetary candidates with orbital periods of 7.3, 38, and 151 d, and minimum masses estimated at $4.88 \pm 0.21$ M$_{\oplus }$, $7.27 \pm 0.70$ M$_{\oplus }$, and $10.59 \pm 1.00$ M$_{\oplus }$, respectively. The outermost planet resides within the (temperate) habitable zone, positioned at a projected distance of $0.029\,{\rm arcsec}$ from its star. Given the close proximity of this planetary system, situated at a distance of 16.8 parsecs, HD 48498 emerges as a promising target (closest super-Earth around FGK stars) for future high-contrast direct imaging and high-resolution spectroscopic studies.