Abstract
The time-variable velocity fields of solar-type stars limit the precision of radial-velocity determinations of their planets’ masses, obstructing detection of Earth twins. Since 2015 July, ...we have been monitoring disc-integrated sunlight in daytime using a purpose-built solar telescope and fibre feed to the HARPS-N stellar radial-velocity spectrometer. We present and analyse the solar radial-velocity measurements and cross-correlation function (CCF) parameters obtained in the first 3 yr of observation, interpreting them in the context of spatially resolved solar observations. We describe a Bayesian mixture-model approach to automated data-quality monitoring. We provide dynamical and daily differential-extinction corrections to place the radial velocities in the heliocentric reference frame, and the CCF shape parameters in the sidereal frame. We achieve a photon-noise-limited radial-velocity precision better than 0.43 m s−1 per 5-min observation. The day-to-day precision is limited by zero-point calibration uncertainty with an RMS scatter of about 0.4 m s−1. We find significant signals from granulation and solar activity. Within a day, granulation noise dominates, with an amplitude of about 0.4 m s−1 and an autocorrelation half-life of 15 min. On longer time-scales, activity dominates. Sunspot groups broaden the CCF as they cross the solar disc. Facular regions temporarily reduce the intrinsic asymmetry of the CCF. The radial-velocity increase that accompanies an active-region passage has a typical amplitude of 5 m s−1 and is correlated with the line asymmetry, but leads it by 3 d. Spectral line-shape variability thus shows promise as a proxy for recovering the true radial velocity.
We describe an analytical method for computing the orbital parameters of a planet from the periodogram of a radial velocity signal. The method is very efficient and provides a good approximation of ...the orbital parameters. The accuracy is mainly limited by the accuracy of the computation of the Fourier decomposition of the signal which is sensitive to sampling and noise. Our method is complementary with more accurate (and more expensive in computer time) numerical algorithms (e.g. Levenberg-Marquardt, Markov chain Monte Carlo, genetic algorithms). Indeed, the analytical approximation can be used as an initial condition to accelerate the convergence of these numerical methods. Our method can be applied iteratively to search for multiple planets in the same system.
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.
State-of-the-art radial-velocity (RV) exoplanet searches are currently limited by RV signals arising from stellar magnetic activity. We analyze solar observations acquired over a 3 yr period during ...the decline of Carrington Cycle 24 to test models of RV variation of Sun-like stars. A purpose-built solar telescope at the High Accuracy Radial-velocity Planet Searcher for the Northern hemisphere (HARPS-N) provides disk-integrated solar spectra, from which we extract RVs and log R HK ′ . The Solar Dynamics Observatory (SDO) provides disk-resolved images of magnetic activity. The Solar Radiation and Climate Experiment (SORCE) provides near-continuous solar photometry, analogous to a Kepler light curve. We verify that the SORCE photometry and HARPS-N log R HK ′ correlate strongly with the SDO-derived magnetic filling factor, while the HARPS-N RV variations do not. To explain this discrepancy, we test existing models of RV variations. We estimate the contributions of the suppression of convective blueshift and the rotational imbalance due to brightness inhomogeneities to the observed HARPS-N RVs. We investigate the time variation of these contributions over several rotation periods, and how these contributions depend on the area of active regions. We find that magnetic active regions smaller than 60 Mm2 do not significantly suppress convective blueshift. Our area-dependent model reduces the amplitude of activity-induced RV variations by a factor of two. The present study highlights the need to identify a proxy that correlates specifically with large, bright magnetic regions on the surfaces of exoplanet-hosting stars.
ABSTRACT
With high-resolution spectroscopy, we can study exoplanet atmospheres and learn about their chemical composition, temperature profiles, and presence of clouds and winds, mainly in hot, giant ...planets. State-of-the-art instrumentation is pushing these studies towards smaller exoplanets. Of special interest are the few planets in the ‘Neptune desert’, a lack of Neptune-sized planets in close orbits around their hosts. Here, we assess the presence of water in one such planet, the bloated super-Neptune WASP-166 b, which orbits an F9-type star in a short orbit of 5.4 d. Despite its close-in orbit, WASP-166 b preserved its atmosphere, making it a benchmark target for exoplanet atmosphere studies in the desert. We analyse two transits observed in the visible with ESPRESSO. We clean the spectra from the Earth’s telluric absorption via principal component analysis, which is crucial to the search for water in exoplanets. We use a cross-correlation-to-likelihood mapping to simultaneously estimate limits on the abundance of water and the altitude of a cloud layer, which points towards a low water abundance and/or high clouds. We tentatively detect a water signal blue-shifted ∼5 $\mathrm{km\, s^{-1}}$ from the planetary rest frame. Injection and retrieval of model spectra show that a solar-composition, cloud-free atmosphere would be detected at high significance. This is only possible in the visible due to the capabilities of ESPRESSO and the collecting power of the VLT. This work provides further insight on the Neptune desert planet WASP-166 b, which will be observed with JWST.
ABSTRACT
The hot Neptune desert, a distinct lack of highly irradiated planets in the size range of Neptune, remains one of the most intriguing results of exoplanet population studies. A deeper ...understanding of the atmosphere of exoplanets sitting at the edge or even within the Neptune desert will allow us to better understand if planetary formation or evolution processes are at the origin of the desert. A detection of sodium in WASP-166 b was presented previously with tentative line broadening at the $3.4\, \sigma$ with the HARPS spectrograph. We update this result with two transits observed with the ESPRESSO spectrograph, confirming the detection in each night and the broadened character of the line. This result marks the first confirmed resolved sodium detection within the Neptune desert. In this work, we additionally highlight the importance of treating low-SNR spectral regions particularly where absorption lines of stellar sodium and planetary sodium overlap at mid transit – an important caveat for future observations of the system.
Abstract
Radial velocity (RV) searches for Earth-mass exoplanets in the habitable zone around Sun-like stars are limited by the effects of stellar variability on the host star. In particular, ...suppression of convective blueshift and brightness inhomogeneities due to photospheric faculae/plage and starspots are the dominant contribution to the variability of such stellar RVs. Gaussian process (GP) regression is a powerful tool for statistically modeling these quasi-periodic variations. We investigate the limits of this technique using 800 days of RVs from the solar telescope on the High Accuracy Radial velocity Planet Searcher for the Northern hemisphere (HARPS-N) spectrograph. These data provide a well-sampled time series of stellar RV variations. Into this data set, we inject Keplerian signals with periods between 100 and 500 days and amplitudes between 0.6 and 2.4 m s
−1
. We use GP regression to fit the resulting RVs and determine the statistical significance of recovered periods and amplitudes. We then generate synthetic RVs with the same covariance properties as the solar data to determine a lower bound on the observational baseline necessary to detect low-mass planets in Venus-like orbits around a Sun-like star. Our simulations show that discovering planets with a larger mass (∼0.5 m s
−1
) using current-generation spectrographs and GP regression will require more than 12 yr of densely sampled RV observations. Furthermore, even with a perfect model of stellar variability, discovering a true exo-Venus (∼0.1 m s
−1
) with current instruments would take over 15 yr. Therefore, next-generation spectrographs and better models of stellar variability are required for detection of such planets.
ABSTRACT
In this paper, we present high-resolution spectroscopic transit observations from ESPRESSO of the super-Neptune WASP-166 b. In addition to spectroscopic ESPRESSO data, we analyse photometric ...data from TESS of six WASP-166 b transits along with simultaneous NGTS observations of the ESPRESSO runs. These observations were used to fit for the planetary parameters as well as assessing the level of stellar activity (e.g. spot crossings, flares) present during the ESPRESSO observations. We utilize the reloaded Rossiter McLaughlin (RRM) technique to spatially resolve the stellar surface, characterizing the centre-to-limb convection-induced variations, and to refine the star–planet obliquity. We find WASP-166 b has a projected obliquity of $\lambda = -15.52^{+2.85}_{-2.76}\, ^{\circ }$ and vsin (i) = 4.97 ± 0.09 km s−1 which is consistent with the literature. We were able to characterize centre-to-limb convective variations as a result of granulation on the surface of the star on the order of a few km s−1 for the first time. We modelled the centre-to-limb convective variations using a linear, quadratic, and cubic model with the cubic being preferred. In addition, by modelling the differential rotation and centre-to-limb convective variations simultaneously, we were able to retrieve a potential antisolar differential rotational shear (α ∼ −0.5) and stellar inclination (i* either 42.03$^{+9.13}_{-9.60}\, ^{\circ }$ or 133.64$^{+8.42}_{-7.98}\, ^{\circ }$ if the star is pointing towards or away from us). Finally, we investigate how the shape of the cross-correlation functions change as a function of limb angle and compare our results to magnetohydrodynamic simulations.
ABSTRACT
In this paper, we present observations of two high-resolution transit data sets obtained with ESPRESSO of the bloated sub-Saturn planet WASP-131 b. We have simultaneous photometric ...observations with NGTS and EulerCam. In addition, we utilized photometric light curves from TESS, WASP, EulerCam, and TRAPPIST of multiple transits to fit for the planetary parameters and update the ephemeris. We spatially resolve the stellar surface of WASP-131 utilizing the Reloaded Rossiter McLaughlin technique to search for centre-to-limb convective variations, stellar differential rotation, and to determine the star–planet obliquity for the first time. We find WASP-131 is misaligned on a nearly retrograde orbit with a projected obliquity of $\lambda = 162.4\substack{+1.3 \\
-1.2}^{\circ }$ . In addition, we determined a stellar differential rotation shear of α = 0.61 ± 0.06 and disentangled the stellar inclination ($i_* = 40.9\substack{+13.3 \\
-8.5}^{\circ }$ ) from the projected rotational velocity, resulting in an equatorial velocity of $v_{\rm {eq}} = 7.7\substack{+1.5 \\
-1.3}$ km s−1. In turn, we determined the true 3D obliquity of $\psi = 123.7\substack{+12.8 \\
-8.0}^{\circ }$ , meaning the planet is on a perpendicular/polar orbit. Therefore, we explored possible mechanisms for the planetary system’s formation and evolution. Finally, we searched for centre-to-limb convective variations where there was a null detection, indicating that centre-to-limb convective variations are not prominent in this star or are hidden within red noise.
.
Since 2012, thanks to the installation of the high-resolution echelle spectrograph in the optical range HARPS-N, the Italian telescope TNG (La Palma) became one of the key facilities for the study ...of the extrasolar planets. In 2014 TNG also offered GIANO to the scientific community, providing a near-infrared (NIR) cross-dispersed echelle spectroscopy covering 0.97-2.45μm at a resolution of 50000. GIANO, although designed for direct light-feed from the telescope at the Nasmyth-B focus, was provisionally mounted on the rotating building and connected via fibers to only available interface at the Nasmyth-A focal plane. The synergy between these two instruments is particularly appealing for a wide range of science cases, especially for the search of exoplanets around young and active stars and the characterisation of their atmosphere. Through the funding scheme “WOW” (a Way to Others Worlds), the Italian National Institute for Astrophysics (INAF) proposed to position GIANO at the focal station for which it was originally designed and the simultaneous use of these spectrographs with the aim to achieve high-resolution spectroscopy in a wide wavelength range (0.383-2.45μm) obtained in a single exposure, giving rise to the project called GIARPS (GIANO-B & HARPS-N). Because of its characteristics, GIARPS can be considered the first and unique worldwide instrument providing not only high resolution in a large wavelength band, but also a high-precision radial velocity measurement both in the visible and in the NIR arm, since in the next future GIANO-B will be equipped with gas absorption cells.