RASSINE: Interactive tool for normalising stellar spectra I. Description and performance of the code

Aims: We provide an open-source code allowing an easy, intuitive, and robust normalisation of spectra. Methods: We developed RASSINE, a Python code for normalising merged 1D spectra through the concepts of convex hulls. The code uses six parameters that can be easily fine-tuned. The code also provides a complete user-friendly interactive interface, including graphical feedback, that helps the user to choose the parameters as easily as possible. To facilitate the normalisation even further, RASSINE can provide a first guess for the parameters that are derived directly from the merged 1D spectrum based on previously performed calibrations. Results: For HARPS spectra of the Sun that were obtained with the HELIOS solar telescope, a continuum accuracy of 0.20% on line depth can be reached after normalisation with RASSINE. This is three times better than with the commonly used method of polynomial fitting. For HARPS spectra of $\alpha$ Cen B, a continuum accuracy of 2.0% is reached. This rather poor accuracy is mainly due to molecular band absorption and the high density of spectral lines in the bluest part of the merged 1D spectrum. When wavelengths shorter than 4500 \AA are excluded, the continuum accuracy improves by up to 1.2%. The line-depth precision on individual spectrum normalisation is estimated to be 0.15%, which can be reduced to the photon-noise limit (0.10%) when a time series of spectra is given as input for RASSINE. Conclusions: With a continuum accuracy higher than the polynomial fitting method and a line-depth precision compatible with photon noise, RASSINE is a tool that can find applications in numerous cases, for example stellar parameter determination, transmission spectroscopy of exoplanet atmospheres, or activity-sensitive line detection.Comment: 13 pages, 9 pages appendix, 9 figure

Modelling stellar variability in archival HARPS data: I -- Rotation and activity properties with multi-dimensional Gaussian Processes

Although instruments for measuring the radial velocities (RVs) of stars now routinely reach sub-meter per second accuracy, the detection of low-mass planets is still very challenging. The rotational modulation and evolution of spots and/or faculae can induce variations in the RVs at the level of a few m/s in Sun-like stars. To overcome this, a multi-dimensional Gaussian Process framework has been developed to model the stellar activity signal using spectroscopic activity indicators together with the RVs. A recently published computationally efficient implementation of this framework, S+LEAF 2, enables the rapid analysis of large samples of targets with sizeable data sets. In this work, we apply this framework to HARPS observations of 268 well-observed targets with precisely determined stellar parameters. Our long-term goal is to quantify the effectiveness of this framework to model and mitigate activity signals for stars of different spectral types and activity levels. In this first paper in the series, we initially focus on the activity indicators (S-index and Bisector Inverse Slope), and use them to a) measure rotation periods for 49 slow rotators in our sample, b) explore the impact of these results on the spin-down of middle-aged late F, G & K stars, and c) explore indirectly how the spot to facular ratio varies across our sample. Our results should provide valuable clues for planning future RV planet surveys such as the Terra Hunting Experiment or the PLATO ground-based follow-up observations program, and help fine-tune current stellar structure and evolution models.Comment: Accepted for publication in MNRA

Modelling stellar variability in archival HARPS data:I -- Rotation and activity properties with multi-dimensional Gaussian Processes

Although instruments for measuring the radial velocities (RVs) of stars now routinely reach sub-meter per second accuracy, the detection of low-mass planets is still very challenging. The rotational modulation and evolution of spots and/or faculae can induce variations in the RVs at the level of a few m/s in Sun-like stars. To overcome this, a multi-dimensional Gaussian Process framework has been developed to model the stellar activity signal using spectroscopic activity indicators together with the RVs. A recently published computationally efficient implementation of this framework, S+LEAF 2, enables the rapid analysis of large samples of targets with sizeable data sets. In this work, we apply this framework to HARPS observations of 268 well-observed targets with precisely determined stellar parameters. Our long-term goal is to quantify the effectiveness of this framework to model and mitigate activity signals for stars of different spectral types and activity levels. In this first paper in the series, we initially focus on the activity indicators (S-index and Bisector Inverse Slope), and use them to a) measure rotation periods for 49 slow rotators in our sample, b) explore the impact of these results on the spin-down of middle-aged late F, G & K stars, and c) explore indirectly how the spot to facular ratio varies across our sample. Our results should provide valuable clues for planning future RV planet surveys such as the Terra Hunting Experiment or the PLATO ground-based follow-up observations program, and help fine-tune current stellar structure and evolution models

Modelling stellar variability in archival HARPS data: I - rotation and activity properties with multi-dimensional Gaussian processes

Although instruments for measuring the radial velocities (RVs) of stars now routinely reach sub-meter per second accuracy, the detection of low-mass planets is still very challenging. The rotational modulation and evolution of spots and/or faculae can induce variations in the RVs at the level of a few m/s in Sun-like stars. To overcome this, a multi-dimensional Gaussian Process framework has been developed to model the stellar activity signal using spectroscopic activity indicators together with the RVs. A recently published computationally efficient implementation of this framework, S+LEAF 2, enables the rapid analysis of large samples of targets with sizeable data sets. In this work, we apply this framework to HARPS observations of 268 well-observed targets with precisely determined stellar parameters. Our long-term goal is to quantify the effectiveness of this framework to model and mitigate activity signals for stars of different spectral types and activity levels. In this first paper in the series, we initially focus on the activity indicators (S-index and Bisector Inverse Slope), and use them to a) measure rotation periods for 49 slow rotators in our sample, b) explore the impact of these results on the spin-down of middle-aged late F, G & K stars, and c) explore indirectly how the spot to facular ratio varies across our sample. Our results should provide valuable clues for planning future RV planet surveys such as the Terra Hunting Experiment or the PLATO ground-based follow-up observations program, and help fine-tune current stellar structure and evolution models

TOI-837 b is a young Saturn-sized exoplanet with a massive 70 M_⊕ core

We present an exhaustive photometric and spectroscopic analysis of TOI-837, a F9/G0 35 Myr young star, hosting a transiting exoplanet, TOI-837 b, with an orbital period of ∼8.32 d. Utilizing data from the Transiting Exoplanet Survey Satellite and groundbased observations, we determine a planetary radius of 0.818+0.034−0.024 RJ for TOI-837 b. Through detailed High Accuracy Radial Velocity Planet Searcherspectroscopic time series analysis, we derive a Dopplersemi-amplitude of 34.7+5.3−5.6 m s−1, corresponding to a planetary mass of 0.379+0.058−0.061 MJ. The derived planetary properties suggest a substantial core of approximately 70 M⊕, constituting about 60 per cent of the planet’s total mass. This finding poses a significant challenge to existing theoretical models of core formation. We propose that future atmospheric observations with JWST could provide insights into resolving ambiguities of TOI-837 b, offering new perspectives on its composition, formation, and evolution

Investigating stellar activity through eight years of Sun-as-a-star observations

Stellar magnetic activity induces both distortions and Doppler-shiftsin the absorption line profiles of Sun-like stars. Those effects produce apparent radial velocity (RV) signals which greatly hamper the search for potentially habitable, Earth-like planets. In this work, we investigate these distortions in the Sun using cross-correlation functions (CCFs), derived from intensive monitoring with the high-precision spectrograph HARPS-N. We show that the RV signal arising from line-shape variations on time-scales associated with the Sun’s rotation and activity cycle can be robustly extracted from the data, reducing the RV dispersion by half. Once these have been corrected, activity-induced Doppler-shifts remain, that are modulated at the solar rotation period, and that are most effectively modelled in the time domain, using Gaussian processes (GPs). Planet signatures are still best retrieved with multidimensonal GPs, when activity is jointly modelled from the raw RVs and indicators of the line width or of the Ca II H & K emission. After GP modelling, the residual RVs exhibit a dispersion of 0.6–0.8 m s−1, likely to be dominated by signals induced by supergranulation. Finally, we find that the statistical properties of the RVs evolve significantly over time, and that this evolution is primarily driven by sunspots, which control the smoothness of the signal. Such evolution, which reduces the sensitivity to long-period planet signatures, is no longer seen in the activity-induced Doppler-shifts, which is promising for long term RV monitoring surveys such as the Terra Hunting Experiment or the PLATO follow-up campaign

TOI-837 b is a young Saturn-sized exoplanet with a massive 70 M_⊕ core

We present an exhaustive photometric and spectroscopic analysis of TOI-837, a F9/G0 35 Myr young star, hosting a transiting exoplanet, TOI-837 b, with an orbital period of ∼8.32 d. Utilizing data from the Transiting Exoplanet Survey Satellite and groundbased observations, we determine a planetary radius of 0.818+0.034−0.024 RJ for TOI-837 b. Through detailed High Accuracy Radial Velocity Planet Searcherspectroscopic time series analysis, we derive a Dopplersemi-amplitude of 34.7+5.3−5.6 m s−1, corresponding to a planetary mass of 0.379+0.058−0.061 MJ. The derived planetary properties suggest a substantial core of approximately 70 M⊕, constituting about 60 per cent of the planet’s total mass. This finding poses a significant challenge to existing theoretical models of core formation. We propose that future atmospheric observations with JWST could provide insights into resolving ambiguities of TOI-837 b, offering new perspectives on its composition, formation, and evolution

Investigating stellar activity through eight years of Sun-as-a-star observations

Funding: BK, SA, OB, HY, and NKOS acknowledge funding from the European Research Council under the European Union’s Horizon 2020 research and innovation programme (grant agreement no 865624, GPRV). MC acknowledges the SNSF support under grant P500PT_211024. FR is funded by the University of Exeter’s College of Engineering, Maths and Physical Sciences, UK. ACC acknowledges support from STFC consolidated grant number ST/V000861/1, and EPSRC grant number EP/Z000181/1 towards the ERC Synergy project REVEAL.Stellar magnetic activity induces both distortions and Doppler-shifts in the absorption line profiles of Sun-like stars. Those effects produce apparent radial velocity (RV) signals which greatly hamper the search for potentially habitable, Earth-like planets. In this work, we investigate these distortions in the Sun using cross-correlation functions (CCFs), derived from intensive monitoring with the high-precision spectrograph HARPS-N. We show that the RV signal arising from line-shape variations on time-scales associated with the Sun’s rotation and activity cycle can be robustly extracted from the data, reducing the RV dispersion by half. Once these have been corrected, activity-induced Doppler-shifts remain, that are modulated at the solar rotation period, and that are most effectively modelled in the time domain, using Gaussian processes (GPs). Planet signatures are still best retrieved with multidimensonal GPs, when activity is jointly modelled from the raw RVs and indicators of the line width or of the Ca ii H & K emission. After GP modelling, the residual RVs exhibit a dispersion of 0.6–0.8 m s−1, likely to be dominated by signals induced by supergranulation. Finally, we find that the statistical properties of the RVs evolve significantly over time, and that this evolution is primarily driven by sunspots, which control the smoothness of the signal. Such evolution, which reduces the sensitivity to long-period planet signatures, is no longer seen in the activity-induced Doppler-shifts, which is promising for long term RV monitoring surveys such as the Terra Hunting Experiment or the PLATO follow-up campaign.Peer reviewe

Investigating stellar activity through eight years of Sun-as-a-star observations

Stellar magnetic activity induces both distortions and Doppler-shiftsin the absorption line profiles of Sun-like stars. Those effects produce apparent radial velocity (RV) signals which greatly hamper the search for potentially habitable, Earth-like planets. In this work, we investigate these distortions in the Sun using cross-correlation functions (CCFs), derived from intensive monitoring with the high-precision spectrograph HARPS-N. We show that the RV signal arising from line-shape variations on time-scales associated with the Sun’s rotation and activity cycle can be robustly extracted from the data, reducing the RV dispersion by half. Once these have been corrected, activity-induced Doppler-shifts remain, that are modulated at the solar rotation period, and that are most effectively modelled in the time domain, using Gaussian processes (GPs). Planet signatures are still best retrieved with multidimensonal GPs, when activity is jointly modelled from the raw RVs and indicators of the line width or of the Ca II H & K emission. After GP modelling, the residual RVs exhibit a dispersion of 0.6–0.8 m s−1, likely to be dominated by signals induced by supergranulation. Finally, we find that the statistical properties of the RVs evolve significantly over time, and that this evolution is primarily driven by sunspots, which control the smoothness of the signal. Such evolution, which reduces the sensitivity to long-period planet signatures, is no longer seen in the activity-induced Doppler-shifts, which is promising for long term RV monitoring surveys such as the Terra Hunting Experiment or the PLATO follow-up campaign

The EXPRES Stellar Signals Project II. State of the Field in Disentangling Photospheric Velocities

Measured spectral shifts due to intrinsic stellar variability (e.g., pulsations, granulation) and activity (e.g., spots, plages) are the largest source of error for extreme-precision radial-velocity (EPRV) exoplanet detection. Several methods are designed to disentangle stellar signals from true center-of-mass shifts due to planets. The Extreme-precision Spectrograph (EXPRES) Stellar Signals Project (ESSP) presents a self-consistent comparison of 22 different methods tested on the same extreme-precision spectroscopic data from EXPRES. Methods derived new activity indicators, constructed models for mapping an indicator to the needed radial-velocity (RV) correction, or separated out shape- and shift-driven RV components. Since no ground truth is known when using real data, relative method performance is assessed using the total and nightly scatter of returned RVs and agreement between the results of different methods. Nearly all submitted methods return a lower RV rms than classic linear decorrelation, but no method is yet consistently reducing the RV rms to sub-meter-per-second levels. There is a concerning lack of agreement between the RVs returned by different methods. These results suggest that continued progress in this field necessitates increased interpretability of methods, high-cadence data to capture stellar signals at all timescales, and continued tests like the ESSP using consistent data sets with more advanced metrics for method performance. Future comparisons should make use of various well-characterized data sets—such as solar data or data with known injected planetary and/or stellar signals—to better understand method performance and whether planetary signals are preserved