Search

In this paper, we present observations of two high-resolution transit datasets obtained with ESPRESSO of the bloated sub-Saturn planet WASP-131~b. We have simultaneous photometric observations with NGTS and EulerCam. In addition, we utilised photometric lightcurves 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 utilising 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 $\alpha = 0.61 \pm 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., Comment: 15 Pages, 10 Figures and 4 Tables Accepted for Publication in MNRAS. arXiv admin note: text overlap with arXiv:2207.10127

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 {\sl 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 utilise the Reloaded Rossiter McLaughlin (RRM) technique to spatially resolve the stellar surface, characterising 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 $v\sin(i) = 4.97 \pm 0.09$~kms$^{-1}$ which is consistent with the literature. We were able to characterise centre-to-limb convective variations as a result of granulation on the surface of the star on the order of a few kms$^{-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 anti-solar differential rotational shear ($\alpha \sim$ -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., Comment: 18 Pages, 13 Figures, 4 Tables. Accepted for publication in MNRAS

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-size 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 days. 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 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., Accepted for publication in MNRAS

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-166b 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., Letter, re-submitted to MNRAS after minor referee report; comments welcome

Computing and using opacities is a key part of modeling and interpreting data of exoplanetary atmospheres. Since the underlying spectroscopic line lists are constantly expanding and currently include up to ~ 10^10 - 10^11 transition lines, the opacity calculator codes need to become more powerful. Here we present major upgrades to the HELIOS-K GPU-accelerated opacity calculator and describe the necessary steps to process large line lists within a reasonable amount of time. Besides performance improvements, we include more capabilities and present a toolbox for handling different atomic and molecular data sets: from downloading and pre-processing the data to performing the opacity calculations in a user-friendly way. HELIOS-K supports line lists from ExoMol, HITRAN, HITEMP, NIST, Kurucz and VALD3. By matching the resolution of 0.1 cm^-1 and cutting length of 25 cm^-1 used by the ExoCross code for timing performance (251 seconds excluding data read-in time), HELIOS-K can process the ExoMol BT2 water line list in 12.5 seconds. Using a resolution of 0.01 cm^-1, it takes 45 seconds - equivalent to about 10^7 lines per second. As a wavenumber resolution of 0.01 cm^-1 suffices for most exoplanetary atmosphere spectroscopic calculations, we adopt this resolution in calculating opacity functions for several hundred atomic and molecular species, and make them freely available on the open-access DACE database. For the opacity calculations of the database, we use a cutting length of 100 cm^-1 for molecules and no cutting length for atoms. Our opacities are available for downloading from https://dace.unige.ch/opacityDatabase and may be visualized using https://dace.unige.ch/opacity., Published in The Astrophysical Journal Supplement Series

The solar telescope connected to HARPS-N has been observing the Sun since the summer of 2015. Such 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-\ms\,systematics that needed to be corrected for. The goal of this manuscript is to i) present a new data reduction software for HARPS-N, ii) demonstrate the improvement brought by this new software on 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. 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 optimized recipes for the spectrograph. We then compared the first three years of HARPS-N solar data reduced with the current and new DRS. 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\ms; 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 on 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 present an efficient technique to derive an index free of instrumental systematics. This paper presents a new data reduction software for HARPS-N, and demonstrates its improvements [...], Comment: 20 pages, 15 figures, 2 tables, version accepted for publication in Astronomy and Astrophysics

Funding: A.C.C. acknowledges support from the Science and Technology Facilities Council (STFC) consolidated grant No. ST/R000824/1. 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. Postprint

Transiting ultra-hot Jupiters are ideal candidates to study the exoplanet atmospheres and their dynamics, particularly by means of high-resolution, high signal-to-noise ratio spectra. One such object is KELT-20b, orbiting the fast rotating A2-type star KELT-20. Many atomic species have already been found in its atmosphere, with blueshifted signals that hints at the presence of a day-to-night side wind. We aimed to observe the atmospheric Rossiter-McLaughlin effect in the ultra-hot Jupiter KELT-20b, and to study any variation of the atmospheric signal during the transit. For this purpose, we analysed five nights of HARPS-N spectra covering five transits of KELT-20b. We computed the mean line profiles of the spectra with a least-squares deconvolution, and then we extracted the stellar radial velocities by fitting them with a rotational broadening profile in order to obtain the radial velocity time-series. We used the mean line profile residuals tomography to analyse the planetary atmospheric signal and its variations. We also used the cross-correlation method to study an already known double-peak feature in the FeI planetary signal. We observed both the classical and the atmospheric Rossiter-McLaughlin effect in the radial velocity time-series. The latter gave us an estimate of the radius of the planetary atmosphere that correlates with the stellar mask used in our work: R(p+atmo)/Rp = 1.13 +/- 0.02). We isolated the planetary atmospheric trace in the tomography, and we found radial velocity variations of the planetary atmospheric signal during transit with an overall blueshift of approximatively 10 km/s, along with small variations in the signal's depth and, less significant, in the full width at half maximum (FWHM). We also find a possible variation in the structure and position of FeI signal in different transits., Comment: 13 pages, 20 figures. Accepted for publication in A&A

Context. The existence of hot Jupiters is still not well understood. Two main channels are thought to be responsible for their current location: a smooth planet migration through the protoplanetary disk or the circularization of an initial highly eccentric orbit by tidal dissipation leading to a strong decrease in the semimajor axis. Different formation scenarios result in different observable effects, such as orbital parameters (obliquity and eccentricity) or frequency of planets at different stellar ages. Aims. In the context of the GAPS Young Objects project, we are carrying out a radial velocity survey with the aim of searching and characterizing young hot-Jupiter planets. Our purpose is to put constraints on evolutionary models and establish statistical properties, such as the frequency of these planets from a homogeneous sample. Methods. Since young stars are in general magnetically very active, we performed multi-band (visible and near-infrared) spectroscopy with simultaneous GIANO-B + HARPS-N (GIARPS) observing mode at TNG. This helps in dealing with stellar activity and distinguishing the nature of radial velocity variations: stellar activity will introduce a wavelength-dependent radial velocity amplitude, whereas a Keplerian signal is achromatic. As a pilot study, we present here the cases of two known hot Jupiters orbiting young stars: HD 285507 b and AD Leo b. Results. Our analysis of simultaneous high-precision GIARPS spectroscopic data confirms the Keplerian nature of the variation in the HD 285507 radial velocities and refines the orbital parameters of the hot Jupiter, obtaining an eccentricity consistent with a circular orbit. Instead, our analysis does not confirm the signal previously attributed to a planet orbiting AD Leo. This demonstrates the power of the multi-band spectroscopic technique when observing active stars., With funding from the Spanish government through the "María de Maeztu Unit of Excellence" accreditation (MDM-2017-0737)

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}{\prime }. 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}{\prime } 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.

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 years 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-minute 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 minutes. On longer timescales, 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 days. Spectral line-shape variability thus shows promise as a proxy for recovering the true radial velocity., 19 pages, 15 figures, accepted for publication in MNRAS

Context. Stellar activity is currently challenging the detection of young planets via the radial velocity (RV) technique. Aims. We attempt to definitively discriminate the nature of the RV variations for the young active K5 star BD+20 1790, for which visible (VIS) RV measurements show divergent results on the existence of a substellar companion. Methods. We compare VIS data with high precision RVs in the near infrared (NIR) range by using the GIANO - B and IGRINS spectrographs. In addition, we present for the first time simultaneous VIS-NIR observations obtained with GIARPS (GIANO - B and HARPS - N) at Telescopio Nazionale Galileo (TNG). Orbital RVs are achromatic, so the RV amplitude does not change at different wavelengths, while stellar activity induces wavelength-dependent RV variations, which are significantly reduced in the NIR range with respect to the VIS. Results. The NIR radial velocity measurements from GIANO - B and IGRINS show an average amplitude of about one quarter with respect to previously published VIS data, as expected when the RV jitter is due to stellar activity. Coeval multi-band photometry surprisingly shows larger amplitudes in the NIR range, explainable with a mixture of cool and hot spots in the same active region. Conclusions. In this work, the claimed massive planet around BD+20 1790 is ruled out by our data. We exploited the crucial role of multi- wavelength spectroscopy when observing young active stars: thanks to facilities like GIARPS that provide simultaneous observations, this method can reach its maximum potential., 12 pages, 7 figures

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., accepted to A&A