Search

We report the discovery of a super-Earth candidate orbiting the nearby mid M dwarf Gl\,725A using the radial velocity (RV) method. The planetary signal has been independently identified using high-precision RVs from the SOPHIE and SPIRou spectrographs, in the optical and near-infrared domains, respectively. We modelled the stellar activity signal jointly with the planet using two Gaussian Processes, one for each instrument to account for the chromaticity of the stellar activity and instrumental systematics, along with a Keplerian model. The signal is significantly detected with a RV semi-amplitude of $1.67\pm0.20$ m/s. The planet Gl 725A b is found to be in an orbit compatible with circular with a period of $11.2201\pm0.0051$ days. We analysed 27 sectors of TESS photometry on which no transit event was found. We determined a minimum mass of $M_{p}\sin{i}=2.78\pm0.35\,M_{\oplus}$ which places the planet in the super-Earth regime. Using Mass-Radius relationships we predict a planetary radius to be between 1.2 and $2.0\,R_{\oplus}$. The proximity of Gl 725A, of only 3.5 pc, makes this new exoplanet one of the closest to Earth and joins the group of S-type low-mass planets in short orbits ($P<15$ d) around close M dwarfs., Comment: 26 pages, 17 figures. Accepted in A&A

Hot Jupiters with close-by planetary companions are rare, with only a handful of them having been discovered so far. This could be due to their suggested dynamical histories, leading to the possible ejection of other planets. TOI-2109 b is special in this regard because it is the hot Jupiter with the closest relative separation from its host star, being separated by less than 2.3 stellar radii. Unexpectedly, transit timing measurements from recently obtained CHEOPS observations show low amplitude transit-timing variations (TTVs). We aim to search for signs of orbital decay and to characterise the apparent TTVs, trying to gain information about a possible companion. We fit the newly obtained CHEOPS light curves using TLCM and extract the resulting mid-transit timings. Successively, we use these measurements in combination with TESS and archival photometric data and radial velocity data to estimate the rate of tidal orbital decay of TOI-2109 b, as well as characterise the TTVs using the N-body code TRADES and the photodynamical approach of PyTTV. We find tentative evidence at $3\sigma$ for orbital decay in the TOI-2109 system, when we correct the mid-transit timings using the best-fitting sinusoidal model of the TTVs. We do not detect additional transits in the available photometric data, but find evidence towards the authenticity of the apparent TTVs, indicating a close-by, outer companion with $P_\mathrm{c} > 1.125\,$d. Due to the fast rotation of the star, the new planetary candidate cannot be detected in the available radial velocity (RV) measurements, and its parameters can only be loosely constrained by our joint TTV and RV modelling. TOI-2109 could join a small group of rare hot Jupiter systems that host close-by planetary companions, only one of which (WASP-47 b) has an outer companion. More high-precision photometric measurements are necessary to confirm the planetary companion., Comment: Accepted for publication in A&A. 21 pages, 21 figures

We present new observations from CHEOPS and TESS to clarify the architecture of the planetary system hosted by the old Galactic thick disk star TOI-561. Our global analysis, which also includes previously published photometric and radial velocity data, incontrovertibly proves that TOI-561 is hosting at least four transiting planets with periods of 0.44 days (TOI-561 b), 10.8 days (TOI-561 c), 25.7 days (TOI-561 d), and 77.1 days (TOI-561 e) and a fifth non-transiting candidate, TOI-561f with a period of 433 days. The precise characterisation of TOI-561's orbital architecture is interesting since old and metal-poor thick disk stars are less likely to host ultra-short period Super-Earths like TOI-561 b. The new period of planet -e is consistent with the value obtained using radial velocity alone and is now known to be $77.14399\pm0.00025$ days, thanks to the new CHEOPS and TESS transits. The new data allowed us to improve its radius ($R_p = 2.517 \pm 0.045 R_{\oplus}$ from 5$\%$ to 2$\%$ precision) and mass ($M_p = 12.4 \pm 1.4 M_{\oplus}$) estimates, implying a density of $\rho_p = 0.778 \pm 0.097 \rho_{\oplus}$. Thanks to recent TESS observations and the focused CHEOPS visit of the transit of TOI-561 e, a good candidate for exomoon searches, the planet's period is finally constrained, allowing us to predict transit times through 2030 with 20-minute accuracy. We present an updated version of the internal structure of the four transiting planets. We finally performed a detailed stability analysis, which confirmed the long-term stability of the outer planet TOI-561 f., Comment: 13 pages, 10 Figures. Accepted on MNRAS. Updated the author list

We present a series of 9 Doppler images of the magnetically active K component of the RS CVn-type binary SZ Psc, based on the high-resolution spectroscopic data collected from 2014 to 2018. We apply least-squares deconvolution to all spectra to extract the average profiles with high signal-to-noise ratios (SNRs) for Doppler imaging. The surface maps of the K subgiant show starspots widely distributed along latitude and longitude. A prominent, non-axisymmetric polar spot around phase 0 is revealed by all images with sufficient phase coverage, which may be a stable feature on the K component. The starspots evolve in a time scale of one month. We have determined the surface shear rate of the K component from the starspot maps reconstructed 10 days apart in 2017 Nov--Dec, through the cross-correlation method. The surface differential rotation parameters are $\Omega_{eq} = 1.591 \pm 0.002$ rad d$^{-1}$ and $\Delta \Omega = 0.035 \pm 0.003$ rad d$^{-1}$. The absorption lines contributed from the tertiary component are detected in all LSD profiles of SZ Psc, and we measure the radial velocity of the binary system and the tertiary component to derive an elliptical orbit with a period of $1530 \pm 3$ days and a mass of $0.75 \pm 0.06$ M$\odot$ for the tertiary component., Comment: 15 pages, 8 figures, accepted for publication in ApJ

Hot Jupiters are giant planets subject to intense stellar radiation. The physical and chemical properties of their atmosphere makes them the most amenable targets for the atmospheric characterization. In this paper we analyze the photometry collected during the secondary eclipses of the hot Jupiter WASP-3 b by CHEOPS, TESS and Spitzer. Our aim is to characterize the atmosphere of the planet by measuring the secondary eclipse depth in several passbands and constrain the planetary dayside spectrum. Our update of the stellar and planetary properties is consistent with previous works. The analysis of the occultations returns an eclipse depth of 92+-21 ppm in the CHEOPS passband, 83+-27 ppm for TESS and >2000 ppm in the IRAC 1-2-4 Spitzer passbands. Using the eclipse depths in the Spitzer bands we propose a set of likely emission spectra which constrain the emission contribution in the \cheops and TESS passbands to approximately a few dozens of parts per million. This allowed us to measure a geometric albedo of 0.21+-0.07 in the CHEOPS passband, while the TESS data lead to a 95\% upper limit of $\sim$0.2. WASP-3 b belongs to the group of ultra-hot Jupiters which are characterized by low Bond albedo (<0.3+-0.1), as predicted by different atmospheric models. On the other hand, it unexpectedly seems to efficiently recirculate the absorbed stellar energy, unlike similar highly irradiated planets. To explain this inconsistency, we propose that other energy recirculation mechanisms may be at play other than advection (for example, dissociation and recombination of H_2). Another possibility is that the observations in different bandpasses probe different atmospheric layers, making the atmospheric analysis difficult without an appropriate modeling of the thermal emission spectrum of WASP-3 b, which is not feasible with the limited spectroscopic data available to date.

Recent observations of giant planets have revealed unexpected bulk densities. Hot Jupiters, in particular, appear larger than expected for their masses compared to planetary evolution models, while warm Jupiters seem denser than expected. These differences are often attributed to the influence of the stellar incident flux, but could they also result from different planet formation processes? Is there a trend linking the planetary density to the chemical composition of the host star? In this work we present the confirmation of three giant planets in orbit around solar analogue stars. TOI-2714 b ($P \simeq 2.5$ d, $R_{\rm p} \simeq 1.22 R_{\rm J}$, $M_{\rm p} = 0.72 M_{\rm J}$) and TOI-2981 b ($P \simeq 3.6$ d, $R_{\rm p} \simeq 1.2 R_{\rm J}$, $M_{\rm p} = 2 M_{\rm J}$) are hot Jupiters on nearly circular orbits, while TOI-4914 b ($P \simeq 10.6$ d, $R_{\rm p} \simeq 1.15 R_{\rm J}$, $M_{\rm p} = 0.72 M_{\rm J}$) is a warm Jupiter with a significant eccentricity ($e = 0.41 \pm 0.02$) that orbits a star more metal-poor ([Fe/H]$~= -0.13$) than most of the stars known to host giant planets. Our radial velocity (RV) follow-up with the HARPS spectrograph allows us to detect their Keplerian signals at high significance (7, 30, and 23$\sigma$, respectively) and to place a strong constraint on the eccentricity of TOI-4914 b (18$\sigma$). TOI-4914 b, with its large radius and low insolation flux ($F_\star < 2 \times 10^8~{\rm erg~s^{-1}~cm^{-2}}$), appears to be more inflated than what is supported by current theoretical models for giant planets. Moreover, it does not conform to the previously noted trend that warm giant planets orbiting metal-poor stars have low eccentricities. This study thus provides insights into the diverse orbital characteristics and formation processes of giant exoplanets, in particular the role of stellar metallicity in the evolution of planetary systems., Comment: Accepted for publication in Astronomy & Astrophysics. 21 pages, 26 figures, and 8 tables. Abstract abridged

K2-24 is a planetary system composed of two transiting low-density Neptunians locked in an almost perfect 2:1 resonance and showing large TTVs, i.e., an excellent laboratory to search for signatures of planetary migration. Previous studies performed with K2, Spitzer and RV data tentatively claimed a significant non-zero eccentricity for one or both planets, possibly high enough to challenge the scenario of pure disk migration through resonant capture. With 13 new CHEOPS light curves (seven of planet -b, six of planet -c), we carried out a global photometric and dynamical re-analysis by including all the available literature data as well. We got the most accurate set of planetary parameters to date for the K2-24 system, including radii and masses at 1% and 5% precision (now essentially limited by the uncertainty on stellar parameters) and non-zero eccentricities $e_b=0.0498_{-0.0018}^{+0.0011}$, $e_c=0.0282_{-0.0007}^{+0.0003}$ detected at very high significance for both planets. Such relatively large values imply the need for an additional physical mechanism of eccentricity excitation during or after the migration stage. Also, while the accuracy of the previous TTV model had drifted by up to 0.5 days at the current time, we constrained the orbital solution firmly enough to predict the forthcoming transits for the next ~15 years, thus enabling an efficient follow-up with top-level facilities such as JWST or ESPRESSO., Comment: 19 pages, 11 figures, 8 tables. Accepted for publication in A&A on September 4, 2024. Typos corrected

The Wide Angle Search for Planets (WASP) survey used transit photometry to discover nearly 200 gas-giant exoplanets and derive their planetary and stellar parameters. Reliable determination of the planetary density depends on accurate measurement of the planet's radius, obtained from the transit depth and photodynamical determination of the stellar radius. The stellar density, and hence the stellar radius are typically determined in a model-independent way from the star's reflex orbital acceleration and the transit profile. Additional flux coming from the system due to a bright, undetected stellar binary companion can, however, potentially dilute the transit curve and radial velocity signal, leading to under-estimation of the planet's mass and radius, and to overestimation of the planet's density. In this study, we cross-check the published radii of all the WASP planet host stars, determined from their transit profiles and radial-velocity curves, against radiometric measurements of stellar radii derived from their angular diameters (via the Infrared Flux method) and trigonometric parallaxes. We identify eight systems showing radiometric stellar radii significantly greater than their published photodynamical values: WASPs 20, 85, 86, 103, 105, 129, 144 and 171. We investigate these systems in more detail to establish plausible ranges of angular and radial-velocity separations within which such "stealth binaries" could evade detection, and deduce their likely orbital periods, mass ratios, and flux ratios. After accounting for the dilution of transit depth and radial velocity amplitude, we find that on average, the planetary densities for the identified stealth binary systems should be reduced by a factor of 1.3.

We report the spectroscopic confirmation and fundamental properties of TOI-757 b, a mini-Neptune on a 17.5-day orbit transiting a bright star ($V = 9.7$ mag) discovered by the TESS mission. We acquired high-precision radial velocity measurements with the HARPS, ESPRESSO, and PFS spectrographs to confirm the planet detection and determine its mass. We also acquired space-borne transit photometry with the CHEOPS space telescope to place stronger constraints on the planet radius, supported with ground-based LCOGT photometry. WASP and KELT photometry were used to help constrain the stellar rotation period. We also determined the fundamental parameters of the host star. We find that TOI-757 b has a radius of $R_{\mathrm{p}} = 2.5 \pm 0.1 R_{\oplus}$ and a mass of $M_{\mathrm{p}} = 10.5^{+2.2}_{-2.1} M_{\oplus}$, implying a bulk density of $\rho_{\text{p}} = 3.6 \pm 0.8$ g cm$^{-3}$. Our internal composition modeling was unable to constrain the composition of TOI-757 b, highlighting the importance of atmospheric observations for the system. We also find the planet to be highly eccentric with $e$ = 0.39$^{+0.08}_{-0.07}$, making it one of the very few highly eccentric planets among precisely characterized mini-Neptunes. Based on comparisons to other similar eccentric systems, we find a likely scenario for TOI-757 b's formation to be high eccentricity migration due to a distant outer companion. We additionally propose the possibility of a more intrinsic explanation for the high eccentricity due to star-star interactions during the earlier epoch of the Galactic disk formation, given the low metallicity and older age of TOI-757., Comment: Accepted for publication in MNRAS; 26 pages, 14 figures, 6 tables

Among the thousands of exoplanets discovered to date, approximately a few hundred gas giants on short-period orbits are classified as "lonely" and only a few are in a multi-planet system with a smaller companion on a close orbit. The processes that formed multi-planet systems hosting gas giants on close orbits are poorly understood, and only a few examples of this kind of system have been observed and well characterised. Within the contest of multi-planet system hosting gas-giant on short orbits, we characterise TOI-1130 system by measuring masses and orbital parameters. This is a 2-transiting planet system with a Jupiter-like planet (c) on a 8.35 days orbit and a Neptune-like planet (b) on an inner (4.07 days) orbit. Both planets show strong anti-correlated transit timing variations (TTVs). Furthermore, radial velocity (RV) analysis showed an additional linear trend, a possible hint of a non-transiting candidate planet on a far outer orbit. Since 2019, extensive transit and radial velocity observations of the TOI-1130 have been acquired using TESS and various ground-based facilities. We present a new photo-dynamical analysis of all available transit and RV data, with the addition of new CHEOPS and ASTEP+ data that achieve the best precision to date on the planetary radii and masses and on the timings of each transit. We were able to model interior structure of planet b constraining the presence of a gaseous envelope of H/He, while it was not possible to assess the possible water content. Furthermore, we analysed the resonant state of the two transiting planets, and we found that they lie just outside the resonant region. This could be the result of the tidal evolution that the system underwent. We obtained both masses of the planets with a precision less than 1.5%, and radii with a precision of about 1% and 3% for planet b and c, respectively., Comment: 22 pages, 10 figures, Accepted for publication by Astronomy and Astrophysics

This work focuses on the analysis of the mean longitudinal magnetic field as a stellar activity tracer in the context of small exoplanet detection and characterisation in radial-velocity (RV) surveys. We use SDO/HMI filtergrams to derive Sun-as-a-star magnetic field measurements, and show that the mean longitudinal magnetic field is an excellent rotational period detector and a useful tracer of the solar magnetic cycle. To put these results into context, we compare the mean longitudinal magnetic field to three common activity proxies derived from HARPS-N Sun-as-a-star data: the full-width at half-maximum, the bisector span and the S-index. The mean longitudinal magnetic field does not correlate with the RVs and therefore cannot be used as a one-to-one proxy. However, with high cadence and a long baseline, the mean longitudinal magnetic field outperforms all other considered proxies as a solar rotational period detector, and can be used to inform our understanding of the physical processes happening on the surface of the Sun. We also test the mean longitudinal magnetic field as a "stellar proxy" on a reduced solar dataset to simulate stellar-like observational sampling. With a Gaussian Process regression analysis, we confirm that the solar mean longitudinal magnetic field is the most effective of the considered indicators, and is the most efficient rotational period indicator over different levels of stellar activity. This work highlights the need for polarimetric time series observations of stars., Comment: 23 pages, 27 figures

Multiplanetary systems spanning the radius valley are ideal testing grounds for exploring the proposed explanations for the observed bimodality in the radius distribution of close-in exoplanets. One such system is HIP 29442 (TOI-469), an evolved K0V star hosting two super-Earths and a sub-Neptune. We observe HIP 29442 with CHEOPS for a total of 9.6 days, which we model jointly with 2 sectors of TESS data to derive planetary radii of $3.410\pm0.046$, $1.551\pm0.045$ and $1.538\pm0.049$ R$_\oplus$ for planets b, c and d, which orbit HIP 29442 with periods of 13.6, 3.5 and 6.4 days. For planet d, this value deviates by more than 3 sigma from the median value reported in the discovery paper, leading us to conclude that caution is required when using TESS photometry to determine the radii of small planets with low per-transit S/N and large gaps between observations. Given the high precision of these new radii, combining them with published RVs from ESPRESSO and HIRES provides us with ideal conditions to investigate the internal structure and formation pathways of the planets in the system. We introduce the publicly available code plaNETic, a fast and robust neural network-based Bayesian internal structure modelling framework. We then apply hydrodynamic models to explore the upper atmospheric properties of these inferred structures. Finally, we identify planetary system analogues in a synthetic population generated with the Bern model for planet formation and evolution. Based on this analysis, we find that the planets likely formed on opposing sides of the water iceline from a protoplanetary disk with an intermediate solid mass. We finally report that the observed parameters of the HIP 29442 system are compatible with both a scenario where the second peak in the bimodal radius distribution corresponds to sub-Neptunes with a pure H/He envelope as well as a scenario with water-rich sub-Neptunes., Comment: 30 pages, 17 figures, accepted for publication in A&A

We identify the bright Am-type star HD 181793 to be a previously-unknown eclipsing, chemically peculiar heartbeat binary, the second of its kind known. The system carries an orbital period of $P = 11.47578275 \pm 0.00000055$ days. We use TESS photometry and LCOGT NRES radial velocity data to build a self-consistent orbital model and determine the fundamental stellar characteristics of the primary. We use a spectral separation method to unveil the secondary and measure the masses of both stars. The radial velocity amplitude of the primary, $K_1 = 47.41+0.13-0.12 km s^{-1}$, gives a mass $M_1 = 1.57 \pm 0.01 $ Msun. The secondary radial velocity amplitude $K_2 = 84.95+0.12-0.09 km s^{-1}$ yields a mass ratio $q = 0.558 \pm 0.002$ and a secondary mass $M_2 = 0.87 \pm 0.01 $ Msun. From the spectral energy distribution and Gaia parallax we find a radius $R_1 = 2.04 \pm 0.05$ Rsun. The grazing transit profile and spectroscopic luminosity ratio indicate $R_2 = 1.04+0.15-0.10$ Rsun, suggesting an early-K spectral type. We show that the heartbeat feature in the TESS light curve can be explained by time-varying ellipsoidal variation, driven by the orbital eccentricity of $e = 0.3056+0.0024-0.0026$, and relativistic beaming of the light of the primary. We find no evidence of tidally-excited oscillations., Comment: 14 pages, 15 figures

Planets orbiting binary systems are relatively unexplored compared to those around single stars. Detections of circumbinary planets and planetary systems offer a first detailed view into our understanding of circumbinary planet formation and dynamical evolution. The BEBOP (Binaries Escorted by Orbiting Planets) radial velocity survey plays a special role in this adventure as it focuses on eclipsing single-lined binaries with an FGK dwarf primary and M dwarf secondary allowing for the highest-radial velocity precision using the HARPS and SOPHIE spectrographs. We obtained 4512 high-resolution spectra for the 179 targets in the BEBOP survey which we used to derive the stellar atmospheric parameters using both equivalent widths and spectral synthesis. We furthermore derive stellar masses, radii, and ages for all targets. With this work, we present the first homogeneous catalogue of precise stellar parameters for these eclipsing single-lined binaries., Comment: accepted for publication in MNRAS

CHEOPS is a space telescope specifically designed to monitor transiting exoplanets orbiting bright stars. In September 2023, CHEOPS completed its nominal mission and remains in excellent operational conditions. The mission has been extended until the end of 2026. Scientific and instrumental data have been collected throughout in-orbit commissioning and nominal operations, enabling a comprehensive analysis of the mission's performance. In this article, we present the results of this analysis with a twofold goal. First, we aim to inform the scientific community about the present status of the mission and what can be expected as the instrument ages. Secondly, we intend for this publication to serve as a legacy document for future missions, providing insights and lessons learned from the successful operation of CHEOPS. To evaluate the instrument performance in flight, we developed a comprehensive monitoring and characterisation programme. It consists of dedicated observations that allow us to characterise the instrument's response. In addition to the standard collection of nominal science and housekeeping data, these observations provide input for detecting, modelling, and correcting instrument systematics, discovering and addressing anomalies, and comparing the instrument's actual performance with expectations. The precision of the CHEOPS measurements has enabled the mission objectives to be met and exceeded. Careful modelling of the instrumental systematics allows the data quality to be significantly improved during the light curve analysis phase, resulting in more precise scientific measurements. CHEOPS is compliant with the driving scientific requirements of the mission. Although visible, the ageing of the instrument has not affected the mission's performance., Comment: Accepted for publication in Astronomy and Astrophysics

HIP 41378 d is a long-period planet that has only been observed to transit twice, three years apart, with K2. According to stability considerations and a partial detection of the Rossiter-McLaughlin effect, $P_\mathrm{d} = 278.36$ d has been determined to be the most likely orbital period. We targeted HIP 41378 d with CHEOPS at the predicted transit timing based on $P_\mathrm{d}= 278.36$ d, but the observations show no transit. We find that large ($>22.4$ hours) transit timing variations (TTVs) could explain this non-detection during the CHEOPS observation window. We also investigated the possibility of an incorrect orbital solution, which would have major implications for our knowledge of this system. If $P_\mathrm{d} \neq 278.36$ d, the periods that minimize the eccentricity would be $101.22$ d and $371.14$ d. The shortest orbital period will be tested by TESS, which will observe HIP 41378 in Sector 88 starting in January 2025. Our study shows the importance of a mission like CHEOPS, which today is the only mission able to make long observations (i.e., from space) to track the ephemeris of long-period planets possibly affected by large TTVs., Comment: Accepted for publication in Astronomy and Astrophysics

The TOI-178 system consists of a nearby late K-dwarf transited by six planets in the super-Earth to mini-Neptune regime, with radii ranging from 1.2 to 2.9 earth radius and orbital periods between 1.9 and 20.7 days. All planets but the innermost one form a chain of Laplace resonances. The fine-tuning and fragility of such orbital configurations ensure that no significant scattering or collision event has taken place since the formation and migration of the planets in the protoplanetary disc, hence providing important anchors for planet formation models. We aim to improve the characterisation of the architecture of this key system, and in particular the masses and radii of its planets. In addition, since this system is one of the few resonant chains that can be characterised by both photometry and radial velocities, we aim to use it as a test bench for the robustness of the planetary mass determination with each technique. We perform a global analysis of all available photometry and radial velocity. We also try different sets of priors on the masses and eccentricity, as well as different stellar activity models, to study their effects on the masses estimated by each method. We show how stellar activity is preventing us from obtaining a robust mass estimation for the three outer planets using radial velocity data alone. We also show that our joint photo-dynamical and radial velocity analysis resulted in a robust mass determination for planets c to g, with precision of 12% for the mass of planet c, and better than 10% for planets d to g. The new precisions on the radii range from 2 to 3%. The understanding of this synergy between photometric and radial velocity measurements will be valuable during the PLATO mission. We also show that TOI-178 is indeed currently locked in the resonant configuration, librating around an equilibrium of the chain.

Many novel methods have been proposed to mitigate stellar activity for exoplanet detection as the presence of stellar activity in radial velocity (RV) measurements is the current major limitation. Unlike traditional methods that model stellar activity in the RV domain, more methods are moving in the direction of disentangling stellar activity at the spectral level. The goal of this paper is to present a novel convolutional neural network-based algorithm that efficiently models stellar activity signals at the spectral level, enhancing the detection of Earth-like planets. We trained a convolutional neural network to build the correlation between the change in the spectral line profile and the corresponding RV, full width at half maximum (FWHM) and bisector span (BIS) values derived from the classical cross-correlation function. This algorithm has been tested on three intensively observed stars: Alpha Centauri B (HD128621), Tau ceti (HD10700), and the Sun. By injecting simulated planetary signals at the spectral level, we demonstrate that our machine learning algorithm can achieve, for HD128621 and HD10700, a detection threshold of 0.5 m/s in semi-amplitude for planets with periods ranging from 10 to 300 days. This threshold would correspond to the detection of a $\sim$4$\mathrm{M}_{\oplus}$ in the habitable zone of those stars. On the HARPS-N solar dataset, our algorithm is even more efficient at mitigating stellar activity signals and can reach a threshold of 0.2 m/s, which would correspond to a 2.2$\mathrm{M}_{\oplus}$ planet on the orbit of the Earth. To the best of our knowledge, it is the first time that such low detection thresholds are reported for the Sun, but also for other stars, and therefore this highlights the efficiency of our convolutional neural network-based algorithm at mitigating stellar activity in RV measurements., Comment: Accepted for publication in A&A

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 solar 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 multi-dimensonal GPs, when activity is jointly modelled from the raw RVs and indicators of the line width or of the Ca II H and K emission. After GP modelling, the residual RVs exhibit a dispersion of 0.6-0.8 m/s, likely to be dominated by signals induced by super-granulation. 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., Comment: Accepted by MNRAS on the 2024 May 20. 17 pages, 14 figures (plus 8 pages of Appendix, 9 Appendix figures)

The TOI-421 planetary system contains two sub-Neptune-type planets and is a prime target to study the formation and evolution of planets and their atmospheres. The inner planet is especially interesting as the existence of a hydrogen-dominated atmosphere at its orbital separation cannot be explained by current formation models without previous orbital migration. We jointly analysed photometric data of three TESS sectors and six CHEOPS visits as well as 156 radial velocity data points to retrieve improved planetary parameters. We also searched for TTVs and modelled the interior structure of the planets. Finally, we simulated the evolution of the primordial H-He atmospheres of the planets using two different modelling frameworks. We determine the planetary radii and masses of TOI-421 b and c to be $R_{\rm b} = 2.64 \pm 0.08 \, R_{\oplus}$, $M_{\rm b} = 6.7 \pm 0.6 \, M_{\oplus}$, $R_{\rm c} = 5.09 \pm 0.07 \, R_{\oplus}$, and $M_{\rm c} = 14.1 \pm 1.4 \, M_{\oplus}$. We do not detect any statistically significant TTV signals. Assuming the presence of a hydrogen-dominated atmosphere, the interior structure modelling results in both planets having extensive envelopes. While the modelling of the atmospheric evolution predicts for TOI-421 b to have lost any primordial atmosphere that it could have accreted at its current orbital position, TOI-421 c could have started out with an initial atmospheric mass fraction somewhere between 10 and 35%. We conclude that the low observed mean density of TOI-421 b can only be explained by either a bias in the measured planetary parameters (e.g. driven by high-altitude clouds) and/or in the context of orbital migration. We also find that the results of atmospheric evolution models are strongly dependent on the employed planetary structure model., Comment: 21 pages, 12 figures, accepted for publication in A&A

Context. White-light stellar flares are proxies for some of the most energetic types of flares, but their triggering mechanism is still poorly understood. As they are associated with strong X and UV emission, their study is particularly relevant to estimate the amount of high-energy irradiation onto the atmospheres of exoplanets, especially those in their stars' habitable zone. Aims. We used the high-cadence, high-photometric capabilities of the CHEOPS and TESS space telescopes to study the detailed morphology of white-light flares occurring in a sample of 130 late-K and M stars, and compared our findings with results obtained at a lower cadence. We developed dedicated software for this purpose. Results. Multi-peak flares represent a significant percentage ($\gtrsim 30$\%) of the detected outburst events. Our findings suggest that high-impulse flares are more frequent than suspected from lower-cadence data, so that the most impactful flux levels that hit close-in exoplanets might be more time-limited than expected. We found significant differences in the duration distributions of single-peak and complex flare components, but not in their peak luminosity. A statistical analysis of the flare parameter distributions provides marginal support for their description with a log-normal instead of a power-law function, leaving the door open to several flare formation scenarios. We tentatively confirmed previous results about quasi-periodic pulsations in high-cadence photometry, report the possible detection of a pre-flare dip, and did not find hints of photometric variability due to an undetected flare background. Conclusions. The high-cadence study of stellar hosts might be crucial to evaluate the impact of their flares on close-in exoplanets, as their impulsive phase emission might otherwise be incorrectly estimated. Future telescopes such as PLATO and Ariel will help in this respect., Comment: 28 pages, 25 figures, 4 tables, to be published in Astronomy & Astrophysics

We aim to constrain the internal structure and composition of HD 15337 b and c, two short-period planets situated on opposite sides of the radius valley, using new transit photometry and radial velocity data. We acquire 6 new transit visits with the CHaracterising ExOPlanet Satellite (CHEOPS) and 32 new radial velocity measurements from the High Accuracy Radial Velocity Planet Searcher (HARPS) to improve the accuracy of the mass and radius estimates for both planets. We reanalyse light curves from TESS sectors 3 and 4 and analyse new data from sector 30, correcting for long-term stellar activity. Subsequently, we perform a joint fit of the TESS and CHEOPS light curves, and all available RV data from HARPS and the Planet Finder Spectrograph (PFS). Our model fits the planetary signals, the stellar activity signal and the instrumental decorrelation model for the CHEOPS data simultaneously. The stellar activity was modelled using a Gaussian-process regression on both the RV and activity indicators. We finally employ a Bayesian retrieval code to determine the internal composition and structure of the planets. We derive updated and highly precise parameters for the HD 15337 system. Our improved precision on the planetary parameters makes HD 15337 b one of the most precisely characterised rocky exoplanets, with radius and mass measurements achieving a precision better than 2\% and 7\%, respectively. We are able to improve the precision of the radius measurement of HD 15337 c to 3\%. Our results imply that the composition of HD 15337 b is predominantly rocky, while HD 15337 c exhibits a gas envelope with a mass of at least $0.01\ M_\oplus$.Our results lay the groundwork for future studies, which can further unravel the atmospheric evolution of these exoplanets and give new insights into their composition and formation history and the causes behind the radius gap., Comment: 17 pages, including appendix

Ultra-hot Jupiters present a unique opportunity to understand the physics and chemistry of planets at extreme conditions. WASP-12b stands out as an archetype of this class of exoplanets. We performed comprehensive analyses of the transits, occultations, and phase curves of WASP-12b by combining new CHEOPS observations with previous TESS and Spitzer data to measure the planet's tidal deformation, atmospheric properties, and orbital decay rate. The planet was modeled as a triaxial ellipsoid parameterized by the second-order fluid Love number, $h_2$, which quantifies its radial deformation and provides insight into the interior structure. We measured the tidal deformation of WASP-12b and estimated a Love number of $h_2=1.55_{-0.49}^{+0.45}$ (at 3.2$\sigma$) from its phase curve. We measured occultation depths of $333\pm24$ppm and $493\pm29$ppm in the CHEOPS and TESS bands, respectively, while the dayside emission spectrum indicates that CHEOPS and TESS probe similar pressure levels in the atmosphere at a temperature of 2900K. We also estimated low geometric albedos of $0.086\pm0.017$ and $0.01\pm0.023$ in the CHEOPS and TESS passbands, respectively, suggesting the absence of reflective clouds in the dayside of the WASP-12b. The CHEOPS occultations do not show strong evidence for variability in the dayside atmosphere of the planet. Finally, we refine the orbital decay rate by 12% to a value of -30.23$\pm$0.82 ms/yr. WASP-12b becomes the second exoplanet, after WASP-103b, for which the Love number has been measured (at 3$sigma$) from the effect of tidal deformation in the light curve. However, constraining the core mass fraction of the planet requires measuring $h_2$ with a higher precision. This can be achieved with high signal-to-noise observations with JWST since the phase curve amplitude, and consequently the induced tidal deformation effect, is higher in the infrared., Comment: accepted for publication in A&A

We present the discovery and characterization of two warm mini-Neptunes transiting the K3V star TOI-815 in a K-M binary system. Analysis of the spectra and rotation period reveal it to be a young star with an age of $200^{+400}_{-200}$Myr. TOI-815b has a 11.2-day period and a radius of 2.94$\pm$0.05$\it{R_{\rm\mathrm{\oplus}}}$ with transits observed by TESS, CHEOPS, ASTEP, and LCOGT. The outer planet, TOI-815c, has a radius of 2.62$\pm$0.10$\it{R_{\rm\mathrm{\oplus}}}$, based on observations of three non-consecutive transits with TESS, while targeted CHEOPS photometry and radial velocity follow-up with ESPRESSO were required to confirm the 35-day period. ESPRESSO confirmed the planetary nature of both planets and measured masses of 7.6$\pm$1.5 $\it{M_{\rm \mathrm{\oplus}}}$ ($\rho_\mathrm{P}$=1.64$^{+0.33}_{-0.31}$gcm$^{-3}$) and 23.5$\pm$2.4$\it{M_{\rm\mathrm{\oplus}}}$ ($\rho_\mathrm{P}$=7.2$^{+1.1}_{-1.0}$gcm$^{-3}$) respectively. Thus, the planets have very different masses, unlike the usual similarity of masses in compact multi-planet systems. Moreover, our statistical analysis of mini-Neptunes orbiting FGK stars suggests that weakly irradiated planets tend to have higher bulk densities compared to those suffering strong irradiation. This could be ascribed to their cooler atmospheres, which are more compressed and denser. Internal structure modeling of TOI-815b suggests it likely has a H-He atmosphere constituting a few percent of the total planet mass, or higher if the planet is assumed to have no water. In contrast, the measured mass and radius of TOI-815c can be explained without invoking any atmosphere, challenging planetary formation theories. Finally, we infer from our measurements that the star is viewed close to pole-on, which implies a spin-orbit misalignment at the 3$\sigma$ level., Comment: 24 pages, 27 figures, 6 tables

Observations of low-mass stars have frequently shown a disagreement between observed stellar radii and radii predicted by theoretical stellar structure models. This ``radius inflation'' problem could have an impact on both stellar and exoplanetary science. We present the final results of our observation programme with the CHEOPS satellite to obtain high-precision light curves of eclipsing binaries with low mass stellar companions (EBLMs). Combined with the spectroscopic orbits of the solar-type companion, we can derive the masses, radii and effective temperatures of 23 M-dwarf stars. We use the PYCHEOPS data analysis software to analyse their primary and secondary occultations. For all but one target, we also perform analyses with TESS light curves for comparison. We have assessed the impact of starspot-induced variation on our derived parameters and account for this in our radius and effective temperature uncertainties using simulated light curves. We observe trends for inflation with both metallicity and orbital separation. We also observe a strong trend in the difference between theoretical and observational effective temperatures with metallicity. There is no such trend with orbital separation. These results are not consistent with the idea that observed inflation in stellar radius combines with lower effective temperature to preserve the luminosity predicted by low-mass stellar models. Our EBLM systems are high-quality and homogeneous measurements that can be used in further studies into radius inflation., Comment: 21 pages, 10 figures, accepted for publication in MNRAS, Supplementary material provided as ancillary files

TOI-732 is an M dwarf hosting two transiting planets, which are located on the two opposite sides of the radius valley. By doubling the number of available space-based observations and increasing the number of radial velocity (RV) measurements, we aim at refining the parameters of TOI-732 b and c. We also aim at using the results to study the slope of both the radius valley and the density valley for a well characterised sample of M dwarf exoplanets. We performed a global MCMC analysis by jointly modelling ground-based light curves, CHEOPS and TESS observations, along with RV time series both taken from the literature and obtained with the MAROON-X spectrograph. The slopes of the M dwarf valleys were quantified via a Support Vector Machine (SVM) procedure. TOI-732 b is an ultra short period planet ($P\sim0.77$ d) with a radius $R_b=1.325_{-0.058}^{+0.057}$ $R_{\oplus}$ and a mass $M_b=2.46\pm0.19$ $M_{\oplus}$ (mean density $\rho_b=5.8_{-0.8}^{+1.0}$ g cm$^{-3}$), while the outer planet at $P\sim12.25$ d has $R_c=2.39_{-0.11}^{+0.10}$ $R_{\oplus}$, $M_c=8.04_{-0.48}^{+0.50}$ $M_{\oplus}$, and thus $\rho_c=3.24_{-0.43}^{+0.55}$ g cm$^{-3}$. Also considering our interior structure calculations TOI-732 b is a super-Earth and TOI-732 c is a mini-Neptune. Following the SVM approach, we quantified $\mathrm{d}\log{R_{p,{\mathrm{valley}}}}/\mathrm{d}\log{P}=-0.065_{-0.013}^{+0.024}$, which is flatter than for Sun-like stars. In line with former analyses, we noted a more filled radius valley for M-planets and we further quantified the density valley slope as $\mathrm{d}\log{\hat{\rho}_{\mathrm{valley}}}/\mathrm{d}\log{P}=-0.02_{-0.04}^{+0.12}$. Compared to FGK stars, the weaker dependence of the position of the radius valley with orbital period might indicate a heavier influence of formation relative to evolution mechanisms in shaping the radius valley around M-dwarfs., Comment: 28 pages (17 in the main text), 18 figures (9 in the main text), 11 tables (7 in the main text). Accepted for publication in A&A

Occultations are windows of opportunity to indirectly peek into the dayside atmosphere of exoplanets. High-precision transit events provide information on the spin-orbit alignment of exoplanets around fast-rotating hosts. We aim to precisely measure the planetary radius and geometric albedo of the ultra-hot Jupiter (UHJ) KELT-20 b as well as the system's spin-orbit alignment. We obtained optical high-precision transits and occultations of KELT-20 b using CHEOPS observations in conjunction with the simultaneous TESS observations. We interpret the occultation measurements together with archival infrared observations to measure the planetary geometric albedo and dayside temperatures. We further use the host star's gravity-darkened nature to measure the system's obliquity. We present time-averaged precise occultation depth of 82(6) ppm measured with seven CHEOPS visits and 131(+8/-7) ppm from the analysis of all available TESS photometry. Using these measurements, we precisely constrain the geometric albedo of KELT-20 b to 0.26(0.04) and the brightness temperature of the dayside hemisphere to 2566(+77/-80) K. Assuming Lambertian scattering law, we constrain the Bond albedo to 0.36(+0.04/-0.05) along with a minimal heat transfer to the nightside. Furthermore, using five transit observations we provide stricter constraints of 3.9(1.1) degrees on the sky-projected obliquity of the system. The aligned orbit of KELT-20 b is in contrast to previous CHEOPS studies that have found strongly inclined orbits for planets orbiting other A-type stars. KELT-20 b's comparably high planetary geometric albedo corroborates a known trend of strongly irradiated planets being more reflective. Finally, we tentatively detect signs of temporal variability in the occultation depths, which might indicate variable cloud cover advecting onto the planetary dayside., Comment: 27 pages, 15 figures, Accepted for publication in Astronomy & Astrophysics

We present the characterisation of the two transiting planets around HD 152843 (TOI 2319, TIC 349488688) using an intensive campaign of HARPS-N radial velocities, and two sectors of TESS data. These data reveal a unique and fascinating system: HD 152843 b and c have near equal masses of around 9 M$_\oplus$ but differing radii of $3.05 \pm 0.11$ R$_\oplus$ and $5.94^{+0.18}_{-0.16}$ R$_\oplus$ , respectively, and orbital periods of $11.62071^{+9.6e-05}_{-0.000106}$ days and $19.502104^{+7.4e-05}_{-8.5e-05}$ days. This indicates that HD 152843 c is in the lowest fifth percentile in density of the known exoplanet population, and has the longest orbital period among these low density planets. Further, HD 152843 c's radius places it in the Saturn valley, the observed lack of planets larger than Neptune, but smaller than Saturn. The orbital periods of these planets indicate they are near a 5:3 mean motion resonance, indicating the possibility of transit timing variations, and hints at the possibility of interaction with a third planet at some point in the evolution of this system. Further, the brightness of the host star and the low density of HD 152843 c make it a key target for atmospheric characterisation., Comment: Submitted for review to the Monthly Notices of the Royal Astronomical Society. 13 pages, 16 figures

We present the characterisation of an inner mini-Neptune in a 9.2292005$\pm$0.0000063 day orbit and an outer mono-transiting sub-Saturn planet in a 95.50$^{+0.36}_{-0.25}$ day orbit around the moderately active, bright (mv=8.9 mag) K5V star TOI-2134. Based on our analysis of five sectors of TESS data, we determine the radii of TOI-2134b and c to be 2.69$\pm$0.16 R$_{e}$ for the inner planet and 7.27$\pm$0.42 R$_{e}$ for the outer one. We acquired 111 radial-velocity spectra with HARPS-N and 108 radial-velocity spectra with SOPHIE. After careful periodogram analysis, we derive masses for both planets via Gaussian Process regression: 9.13$^{+0.78}_{-0.76}$ M$_{e}$ for TOI-2134b and 41.86$^{+7.69}_{-7.83}$ M$_{e}$ for TOI-2134c. We analysed the photometric and radial-velocity data first separately, then jointly. The inner planet is a mini-Neptune with density consistent with either a water-world or a rocky core planet with a low-mass H/He envelope. The outer planet has a bulk density similar to Saturn's. The outer planet is derived to have a significant eccentricity of 0.67$^{+0.05}_{-0.06}$ from a combination of photometry and RVs. We compute the irradiation of TOI-2134c as 1.45$\pm$0.10 times the bolometric flux received by Earth, positioning it for part of its orbit in the habitable sone of its system. We recommend further RV observations to fully constrain the orbit of TOI-2134c. With an expected Rossiter-McLaughlin (RM) effect amplitude of 7.2$\pm$1.3 m/s, we recommend TOI-2134c for follow-up RM analysis to study the spin-orbit architecture of the system. We calculate the Transmission Spectroscopy Metric, and both planets are suitable for bright-mode NIRCam atmospheric characterisation., Comment: 24 pages, 7 tables, 14 figures

Is the population of close-in planets orbiting M dwarfs sculpted by thermally driven escape or is it a direct outcome of the planet formation process? A number of recent empirical results strongly suggest the latter. However, the unique architecture of the TOI-1266 system presents a challenge to models of planet formation and atmospheric escape given its seemingly "inverted" architecture of a large sub-Neptune ($P_b=10.9$ days, $R_{p,b}=2.62\pm 0.11\, \mathrm{R}_{\oplus}$) orbiting interior to that of the system's smaller planet ($P_c=18.8$ days, $R_{p,c}=2.13\pm 0.12\, \mathrm{R}_{\oplus}$). Here we present revised planetary radii based on new TESS and diffuser-assisted ground-based transit observations, and characterize both planetary masses using a set of 145 radial velocity measurements from HARPS-N ($M_{p,b}=4.23\pm 0.69\, \mathrm{M}_{\oplus}, M_{p,c}=2.88\pm 0.80\, \mathrm{M}_{\oplus}$). Our analysis also reveals a third planet candidate ($P_d=32.3$ days, $M_{p,d}\sin{i} = 4.59^{+0.96}_{-0.94}\, \mathrm{M}_{\oplus}$), which if real, would form a chain of near 5:3 period ratios, although the system is likely not in a mean motion resonance. Our results indicate that TOI-1266 b and c are among the lowest density sub-Neptunes around M dwarfs and likely exhibit distinct bulk compositions of a gas-enveloped terrestrial ($X_{\mathrm{env},b}=5.5\pm 0.7$%) and a water-rich world (WMF$_c=59\pm 14$%), which is supported by hydrodynamic escape models. If distinct bulk compositions are confirmed through atmospheric characterization, the system's unique architecture would represent an interesting test case of inside-out sub-Neptune formation at pebble traps., Comment: Accepted for publication in MNRAS. 21 pages. Our spectroscopic time series are included in the arXiv source files as table6.csv

HD 139139 (a.k.a. 'The Random Transiter') is a star that exhibited enigmatic transit-like features with no apparent periodicity in K2 data. The shallow depth of the events ($\sim$200 ppm -- equivalent to transiting objects with radii of $\sim$1.5 R$_\oplus$ in front of a Sun-like star), and their non-periodicity, constitutes a challenge for the photometric follow-up of this star. The goal of this study is to confirm with independent measurements the presence of shallow, non-periodic transit-like features on this object. We performed observations with CHEOPS, for a total accumulated time of 12.75 d, distributed in visits of roughly 20 h in two observing campaigns in years 2021 and 2022. The precision of the data is sufficient to detect 150 ppm features with durations longer than 1.5 h. We use the duration and times of the events seen in the K2 curve to estimate how many should have been detected in our campaigns, under the assumption that their behaviour during the CHEOPS observations would be the same as in the K2 data of 2017. We do not detect events with depths larger than 150 ppm in our data set. If the frequency, depth, and duration of the events were the same as in the K2 campaign, we estimate the probability of having missed all events due to our limited observing window would be 4.8 %. We suggest three different scenarios to explain our results: 1) Our observing window was not long enough, and the events were missed with the estimated 4.8 % probability. 2) The events recorded in the K2 observations were time critical, and the mechanism producing them was either not active in the 2021 and 2022 campaigns or created shallower events under our detectability level. 3) The enigmatic events in the K2 data are the result of an unidentified and infrequent instrumental noise in the original data set or its data treatment., Comment: 7 pages, 3 figures, accepted for publication in A&A. Language-corrected version

Multiwavelength photometry of the secondary eclipses of extrasolar planets is able to disentangle the reflected and thermally emitted light radiated from the planetary dayside. This leads to the measurement of the planetary geometric albedo $A_g$, which is an indicator of the presence of clouds in the atmosphere, and the recirculation efficiency $\epsilon$, which quantifies the energy transport within the atmosphere. In this work we aim to measure $A_g$ and $\epsilon$ for the planet WASP-178 b, a highly irradiated giant planet with an estimated equilibrium temperature of 2450 K.} We analyzed archival spectra and the light curves collected by Cheops and Tess to characterize the host WASP-178, refine the ephemeris of the system and measure the eclipse depth in the passbands of the two respective telescopes. We measured a marginally significant eclipse depth of 70$\pm$40 ppm in the Tess passband and statistically significant depth of 70$\pm$20 ppm in the Cheops passband. Combining the eclipse depth measurement in the Cheops (lambda_eff=6300 AA) and Tess (lambda_eff=8000 AA) passbands we constrained the dayside brightness temperature of WASP-178 b in the 2250-2800 K interval. The geometric albedo 0.1<$\rm A_g$<0.35 is in general agreement with the picture of poorly reflective giant planets, while the recirculation efficiency $\epsilon>$0.7 makes WASP-178 b an interesting laboratory to test the current heat recirculation models., Comment: Accepted by Astronomy and Astrophysics on 31/08/2023

We present an analysis of Sun-as-a-star observations from four different high-resolution, stabilized spectrographs -- HARPS, HARPS-N, EXPRES, and NEID. With simultaneous observations of the Sun from four different instruments, we are able to gain insight into the radial velocity precision and accuracy delivered by each of these instruments and isolate instrumental systematics that differ from true astrophysical signals. With solar observations, we can completely characterize the expected Doppler shift contributed by orbiting Solar System bodies and remove them. This results in a data set with measured velocity variations that purely trace flows on the solar surface. Direct comparisons of the radial velocities measured by each instrument show remarkable agreement with residual intra-day scatter of only 15-30 cm/s. This shows that current ultra-stabilized instruments have broken through to a new level of measurement precision that reveals stellar variability with high fidelity and detail. We end by discussing how radial velocities from different instruments can be combined to provide powerful leverage for testing techniques to mitigate stellar signals., Comment: 17 pages, 9 figures, accepted for publication

The TOI-178 system consists of a nearby late K-dwarf transited by six planets in the super-Earth to mini-Neptune regime, with orbital periods between 1.9 and 20.7 days. All planets but the innermost one form a chain of Laplace resonances. Mass estimates derived from a preliminary radial velocity (RV) dataset suggest that the planetary densities do not decrease in a monotonic way with the orbital distance to the star, contrary to what one would expect based on simple formation and evolution models. To improve the characterisation of this key system and prepare for future studies (in particular with JWST), we perform a detailed photometric study based on 40 new CHEOPS visits, one new TESS sector, as well as previously published CHEOPS, TESS, and NGTS data. First we perform a global analysis of the 100 transits contained in our data to refine the transit parameters of the six planets and study their transit timing variations (TTVs). We then use our extensive dataset to place constraints on the radii and orbital periods of potential additional transiting planets in the system. Our analysis significantly refines the transit parameters of the six planets, most notably their radii, for which we now obtain relative precisions $\lesssim$3%, with the exception of the smallest planet $b$ for which the precision is 5.1%. Combined with the RV mass estimates, the measured TTVs allow us to constrain the eccentricities of planets $c$ to $g$, which are found to be all below 0.02, as expected from stability requirements. Taken alone, the TTVs also suggest a higher mass for planet $d$ than the one estimated from the RVs, which had been found to yield a surprisingly low density for this planet. However, the masses derived from the current TTV dataset are very prior-dependent and further observations, over a longer temporal baseline, are needed to deepen our understanding of this iconic planetary system., Comment: 20 pages, 5 figures, 9 tables. Accepted for publication in A&A

Ultra-short period planets (USPs) are a unique class of super-Earths with an orbital period of less than a day and hence subject to intense radiation from their host star. While most of them are consistent with bare rocks, some show evidence of a heavyweight envelope, which could be a water layer or a secondary metal-rich atmosphere sustained by an outgassing surface. Much remains to be learned about the nature of USPs. The prime goal of the present work is to study the bulk planetary properties and atmosphere of TOI-561b, through the study of its transits and occultations. We obtained ultra-precise transit photometry of TOI-561b with CHEOPS and performed a joint analysis of this data with four TESS sectors. Our analysis of TOI-561b transit photometry put strong constraints on its properties, especially on its radius, Rp=1.42 +/- 0.02 R_Earth (at ~2% error). The internal structure modelling of the planet shows that the observations are consistent with negligible H/He atmosphere, however requiring other lighter materials, in addition to pure iron core and silicate mantle to explain the observed density. We find that this can be explained by the inclusion of a water layer in our model. We searched for variability in the measured Rp/R* over time to trace changes in the structure of the planetary envelope but none found within the data precision. In addition to the transit event, we tentatively detect occultation signal in the TESS data with an eclipse depth of ~27 +/- 11 ppm. Using the models of outgassed atmospheres from the literature we find that the thermal emission from the planet can mostly explain the observation. Based on this, we predict that NIR/MIR observations with JWST should be able to detect silicate species in the atmosphere of the planet. This could also reveal important clues about the planetary interior and help disentangle planet formation and evolution models., Comment: 17 pages, 10 + 3 figures, 4 tables, accepted for publication in A&A (abstract abbreviated)

Minimising the impact of stellar variability in Radial Velocity (RV) measurements is a critical challenge in achieving the 10 cm s$^{-1}$ precision needed to hunt for Earth twins. Since 2012, a dedicated programme has been underway with HARPS-N, to conduct a blind RV Rocky Planets Search (RPS) around bright stars in the Northern Hemisphere. Here we describe the results of a comprehensive search for planetary systems in two RPS targets, HD 166620 and HD 144579. Using wavelength-domain line-profile decorrelation vectors to mitigate the stellar activity and performing a deep search for planetary reflex motions using a trans-dimensional nested sampler, we found no significant planetary signals in the data sets of either of the stars. We validated the results via data-splitting and injection recovery tests. Additionally, we obtained the 95th percentile detection limits on the HARPS-N RVs. We found that the likelihood of finding a low-mass planet increases noticeably across a wide period range when the inherent stellar variability is corrected for using scalpels U-vectors. We are able to detect planet signals with $M\sin i \leq 1$ M$_\oplus$ for orbital periods shorter than 10 days. We demonstrate that with our decorrelation technique, we are able to detect signals as low as 54 cm s$^{-1}$, which brings us closer to the calibration limit of 50 cm s$^{-1}$ demonstrated by HARPS-N. Therefore, we show that we can push down towards the RV precision required to find Earth analogues using high-precision radial velocity data with novel data-analysis techniques., Comment: 7 tables, 24 figures (including those in appendix)

55 Cnc e is an ultra-short period super-Earth transiting a Sun-like star. Previous observations in the optical range detected a time-variable flux modulation that is phased with the planetary orbital period, whose amplitude is too large to be explained by reflected light and thermal emission alone. The goal of the study is to investigate the origin of the variability and timescale of the phase-curve modulation in 55 Cnc e. To this end, we used the CHaracterising ExOPlanet Satellite (CHEOPS), whose exquisite photometric precision provides an opportunity to characterise minute changes in the phase curve from one orbit to the next. CHEOPS observed 29 individual visits of 55 Cnc e between March 2020 and February 2022. Based on these observations, we investigated the different processes that could be at the origin of the observed modulation. In particular, we built a toy model to assess whether a circumstellar torus of dust driven by radiation pressure and gravity might match the observed flux variability timescale. We find that the phase-curve amplitude and peak offset of 55 Cnc e do vary between visits. The sublimation timescales of selected dust species reveal that silicates expected in an Earth-like mantle would not survive long enough to explain the observed phase-curve modulation. We find that silicon carbide, quartz, and graphite are plausible candidates for the circumstellar torus composition because their sublimation timescales are long. The extensive CHEOPS observations confirm that the phase-curve amplitude and offset vary in time.We find that dust could provide the grey opacity source required to match the observations. However, the data at hand do not provide evidence that circumstellar material with a variable grain mass per unit area causes the observed variability. Future observations with the James Webb Space Telescope promise exciting insights into this iconic super-Earth., Comment: 27 pages, 22 figures. Accepted for publication on A&A

We report the confirmation of three exoplanets transiting TOI-4010 (TIC-352682207), a metal-rich K dwarf observed by TESS in Sectors 24, 25, 52, and 58. We confirm these planets with HARPS-N radial velocity observations and measure their masses with 8 - 12% precision. TOI-4010 b is a sub-Neptune ($P = 1.3$ days, $R_{p} = 3.02_{-0.08}^{+0.08}~R_{\oplus}$, $M_{p} = 11.00_{-1.27}^{+1.29}~M_{\oplus}$) in the hot Neptune desert, and is one of the few such planets with known companions. Meanwhile, TOI-4010 c ($P = 5.4$ days, $R_{p} = 5.93_{-0.12}^{+0.11}~R_{\oplus}$, $M_{p} = 20.31_{-2.11}^{+2.13}~M_{\oplus}$) and TOI-4010 d ($P = 14.7$ days, $R_{p} = 6.18_{-0.14}^{+0.15}~R_{\oplus}$, $M_{p} = 38.15_{-3.22}^{+3.27}~M_{\oplus}$) are similarly-sized sub-Saturns on short-period orbits. Radial velocity observations also reveal a super-Jupiter-mass companion called TOI-4010 e in a long-period, eccentric orbit ($P \sim 762$ days and $e \sim 0.26$ based on available observations). TOI-4010 is one of the few systems with multiple short-period sub-Saturns to be discovered so far., Comment: 26 pages, 16 figures, published in AJ; (v3) added missing citation

We report the discovery of two warm sub-Neptunes transiting the bright (G = 9.5 mag) K-dwarf HD 15906 (TOI 461, TIC 4646810). This star was observed by the Transiting Exoplanet Survey Satellite (TESS) in sectors 4 and 31, revealing two small transiting planets. The inner planet, HD 15906 b, was detected with an unambiguous period but the outer planet, HD 15906 c, showed only two transits separated by $\sim$ 734 days, leading to 36 possible values of its period. We performed follow-up observations with the CHaracterising ExOPlanet Satellite (CHEOPS) to confirm the true period of HD 15906 c and improve the radius precision of the two planets. From TESS, CHEOPS and additional ground-based photometry, we find that HD 15906 b has a radius of 2.24 $\pm$ 0.08 R$_\oplus$ and a period of 10.924709 $\pm$ 0.000032 days, whilst HD 15906 c has a radius of 2.93$^{+0.07}_{-0.06}$ R$_\oplus$ and a period of 21.583298$^{+0.000052}_{-0.000055}$ days. Assuming zero bond albedo and full day-night heat redistribution, the inner and outer planet have equilibrium temperatures of 668 $\pm$ 13 K and 532 $\pm$ 10 K, respectively. The HD 15906 system has become one of only six multiplanet systems with two warm ($\lesssim$ 700 K) sub-Neptune sized planets transiting a bright star (G $\leq$ 10 mag). It is an excellent target for detailed characterisation studies to constrain the composition of sub-Neptune planets and test theories of planet formation and evolution., Comment: 30 pages, 20 figures, 11 tables (including appendix). Published in MNRAS

Multi-planet systems are important sources of information regarding the evolution of planets. However, the long-period planets in these systems often escape detection. HD 22946 is a bright star around which 3 transiting planets were identified via TESS photometry, but the true orbital period of the outermost planet d was unknown until now. We aim to use CHEOPS to uncover the true orbital period of HD 22946d and to refine the orbital and planetary properties of the system, especially the radii of the planets. We used the available TESS photometry of HD 22946 and observed several transits of the planets b, c, and d using CHEOPS. We identified 2 transits of planet d in the TESS photometry, calculated the most probable period aliases based on these data, and then scheduled CHEOPS observations. The photometric data were supplemented with ESPRESSO radial velocity data. Finally, a combined model was fitted to the entire dataset. We successfully determined the true orbital period of the planet d to be 47.42489 $\pm$ 0.00011 d, and derived precise radii of the planets in the system, namely 1.362 $\pm$ 0.040 R$_\oplus$, 2.328 $\pm$ 0.039 R$_\oplus$, and 2.607 $\pm$ 0.060 R$_\oplus$ for planets b, c, and d, respectively. Due to the low number of radial velocities, we were only able to determine 3$\sigma$ upper limits for these respective planet masses, which are 13.71 M$_\oplus$, 9.72 M$_\oplus$, and 26.57 M$_\oplus$. We estimated that another 48 ESPRESSO radial velocities are needed to measure the predicted masses of all planets in HD 22946. Planet c appears to be a promising target for future atmospheric characterisation. We can also conclude that planet d, as a warm sub-Neptune, is very interesting because there are only a few similar confirmed exoplanets to date. Such objects are worth investigating in the near future, for example in terms of their composition and internal structure.

HIP 9618 (HD 12572, TOI-1471, TIC 306263608) is a bright ($G=9.0$ mag) solar analogue. TESS photometry revealed the star to have two candidate planets with radii of $3.9 \pm 0.044$ $R_\oplus$ (HIP 9618 b) and $3.343 \pm 0.039$ $R_\oplus$ (HIP 9618 c). While the 20.77291 day period of HIP 9618 b was measured unambiguously, HIP 9618 c showed only two transits separated by a 680-day gap in the time series, leaving many possibilities for the period. To solve this issue, CHEOPS performed targeted photometry of period aliases to attempt to recover the true period of planet c, and successfully determined the true period to be 52.56349 d. High-resolution spectroscopy with HARPS-N, SOPHIE and CAFE revealed a mass of $10.0 \pm 3.1 M_\oplus$ for HIP 9618 b, which, according to our interior structure models, corresponds to a $6.8\pm1.4\%$ gas fraction. HIP 9618 c appears to have a lower mass than HIP 9618 b, with a 3-sigma upper limit of $< 18M_\oplus$. Follow-up and archival RV measurements also reveal a clear long-term trend which, when combined with imaging and astrometric information, reveal a low-mass companion ($0.08^{+0.12}_{-0.05} M_\odot$) orbiting at $26^{+19}_{-11}$ au. This detection makes HIP 9618 one of only five bright ($K<8$ mag) transiting multi-planet systems known to host a planet with $P>50$ d, opening the door for the atmospheric characterisation of warm ($T_{\rm eq}<750$ K) sub-Neptunes., Comment: 19 pages, 16 figures, 9 tables. Accepted at MNRAS. CHEOPS, RV and ground-based photometric data is available on CDS at https://cdsarc.cds.unistra.fr/viz-bin/cat/J/MNRAS/523/3069

A large sample of long-period giant planets has been discovered thanks to long-term radial velocity surveys, but only a few dozen of these planets have a precise radius measurement. Transiting gas giants are crucial targets for the study of atmospheric composition across a wide range of equilibrium temperatures and for shedding light on the formation and evolution of planetary systems. Indeed, compared to hot Jupiters, the atmospheric properties and orbital parameters of cooler gas giants are unaltered by intense stellar irradiation and tidal effects. We identify long-period planets in the Transiting Exoplanet Survey Satellite (TESS) data as duo-transit events. To solve the orbital periods of TESS duo-transit candidates, we use the CHaracterising ExOPlanet Satellite (CHEOPS) to observe the highest-probability period aliases in order to discard or confirm a transit event at a given period. We also collect spectroscopic observations with CORALIE and HARPS in order to confirm the planetary nature and measure the mass of the candidates. We report the discovery of a warm transiting Neptune-mass planet orbiting TOI-5678. After four non-detections corresponding to possible periods, CHEOPS detected a transit event matching a unique period alias. Joint modeling reveals that TOI-5678 hosts a 47.73 day period planet. TOI-5678 b has a mass of 20 (+-4) Me and a radius of 4.91 (+-0.08 Re) . Using interior structure modeling, we find that TOI-5678 b is composed of a low-mass core surrounded by a large H/He layer with a mass of 3.2 (+1.7, -1.3) Me. TOI-5678 b is part of a growing sample of well-characterized transiting gas giants receiving moderate amounts of stellar insolation (11 Se). Precise density measurement gives us insight into their interior composition, and the objects orbiting bright stars are suitable targets to study the atmospheric composition of cooler gas giants., Comment: 17 pages, 10 figures, accepted to A&A

The photospheric unsigned magnetic flux has been shown to be highly correlated with radial velocity (RV) variations caused by solar surface activity. This activity indicator is therefore a prime candidate to unlock the potential of RV surveys to discover Earth twins orbiting Sun-like stars. We show for the first time how a precise proxy of the unsigned magnetic flux ($\Delta\alpha B^2$) can be obtained from Sun-as-a-star intensity spectra by harnessing the magnetic information contained in over 4000 absorption lines in the wavelength range from 380 to 690 nm. This novel activity proxy can thus be obtained from the same spectra from which RVs are routinely extracted. We derived $\Delta\alpha B^2$ from 500 randomly selected spectra from the HARPS-N public solar data set, which spans from 2015 to 2018. We compared our estimates with the unsigned magnetic flux values from the Solar Dynamics Observatory (SDO) finding excellent agreement (median absolute deviation: 4.9 per cent). The extracted indicator $\Delta\alpha B^2$ correlates with SDO's unsigned magnetic flux estimates on the solar rotational timescale (Pearson correlation coefficient 0.67) and on the three-year timescale of our data set (correlation coefficient 0.91). We find correlations of $\Delta\alpha B^2$ with the HARPS-N solar RV variations of 0.49 on the rotational timescale and 0.78 on the three-year timescale. The Pearson correlation of $\Delta\alpha B^2$ with the RVs is found to be greater than the correlation of the classical activity indicators with the RVs. For solar-type stars, $\Delta\alpha B^2$ therefore represents the best simultaneous activity proxy known to date., Comment: Typo corrected in Eq. 6. Note added to Eq. 16 highlighting its similarity to the second derivative. Results unchanged. 17 pages, 10 figures, published in MNRAS

Moons orbiting exoplanets ("exomoons") may hold clues about planet formation, migration, and habitability. In this work, we investigate the plausibility of exomoons orbiting the temperate ($T_\text{eq}=294$ K) giant ($R = 9.2$ R$_\oplus$) planet HIP 41378 f, which has been shown to have a low apparent bulk density of $0.09\,\text{g}\,\text{cm}^{-3}$ and a flat near-infrared transmission spectrum, hinting that it may possess circumplanetary rings. Given this planet's long orbital period ($P\approx1.5$ yr), it has been suggested that it may also host a large exomoon. Here, we analyze the orbital stability of a hypothetical exomoon with a satellite-to-planet mass ratio of 0.0123 orbiting HIP 41378 f. Combining a new software package, astroQTpy, with REBOUND and EqTide, we conduct a series of N-body and tidal migration simulations, demonstrating that satellites up to this size are largely stable against dynamical escape and collisions. We simulate the expected transit signal from this hypothetical exomoon and show that current transit observations likely cannot constrain the presence of exomoons orbiting HIP 41378 f, though future observations may be capable of detecting exomoons in other systems. Finally, we model the combined transmission spectrum of HIP 41378 f and a hypothetical moon with a low-metallicity atmosphere, and show that the total effective spectrum would be contaminated at the $\sim$10 ppm level. Our work not only demonstrates the feasibility of exomoons orbiting HIP 41378 f, but also shows that large exomoons may be a source of uncertainty in future high-precision measurements of exoplanet systems., Comment: 28 pages, 10 figures, 2 tables; accepted to AJ

TOI-1055 is a Sun-like star known to host a transiting Neptune-sized planet on a 17.5-day orbit (TOI-1055 b). Radial velocity (RV) analyses carried out by two independent groups using nearly the same set of HARPS spectra have provided measurements of planetary masses that differ by $\sim$ 2$\sigma$. Our aim in this work is to solve the inconsistency in the published planetary masses by significantly extending the set of HARPS RV measurements and employing a new analysis tool that is able to account and correct for stellar activity. Our further aim was to improve the precision on measurements of the planetary radius by observing two transits of the planet with the CHEOPS space telescope. We fit a skew normal (SN) function to each cross correlation function extracted from the HARPS spectra to obtain RV measurements and hyperparameters to be used for the detrending. We evaluated the correlation changes of the hyperparameters along the RV time series using the breakpoint technique. We performed a joint photometric and RV analysis using a Markov chain Monte Carlo (MCMC) scheme to simultaneously detrend the light curves and the RV time series. We firmly detected the Keplerian signal of TOI-1055 b, deriving a planetary mass of $M_b=20.4_{-2.5}^{+2.6} M_{\oplus}$ ($\sim$12%). This value is in agreement with one of the two estimates in the literature, but it is significantly more precise. Thanks to the TESS transit light curves combined with exquisite CHEOPS photometry, we also derived a planetary radius of $R_b=3.490_{-0.064}^{+0.070} R_{\oplus}$ ($\sim$1.9%). Our mass and radius measurements imply a mean density of $\rho_b=2.65_{-0.35}^{+0.37}$ g cm$^{-3}$ ($\sim$14%). We further inferred the planetary structure and found that TOI-1055 b is very likely to host a substantial gas envelope with a mass of $0.41^{+0.34}_{-0.20}$ M$_\oplus$ and a thickness of $1.05^{+0.30}_{-0.29}$ R$_\oplus$., Comment: 13 pages, 6 figures, 5 tables. Accepted for publication in A&A

DE Boo is a unique system, with an edge-on view through the debris disk around the star. The disk, which is analogous to the Kuiper belt in the Solar System, was reported to extend from 74 to 84 AU from the central star. The high photometric precision of the Characterising Exoplanet Satellite (CHEOPS) provided an exceptional opportunity to observe small variations in the light curve due to transiting material in the disk. This is a unique chance to investigate processes in the debris disk. Photometric observations of DE Boo of a total of four days were carried out with CHEOPS. Photometric variations due to spots on the stellar surface were subtracted from the light curves by applying a two-spot model and a fourth-order polynomial. The photometric observations were accompanied by spectroscopic measurements with the 1m RCC telescope at Piszk\'estet\H{o} and with the SOPHIE spectrograph in order to refine the astrophysical parameters of DE Boo. We present a detailed analysis of the photometric observation of DE Boo. We report the presence of nonperiodic transient features in the residual light curves with a transit duration of 0.3-0.8 days. We calculated the maximum distance of the material responsible for these variations to be 2.47 AU from the central star, much closer than most of the mass of the debris disk. Furthermore, we report the first observation of flaring events in this system. We interpreted the transient features as the result of scattering in an inner debris disk around DE Boo. The processes responsible for these variations were investigated in the context of interactions between planetesimals in the system., Comment: 10 pages, 4 figures, 5 tables; accepted by Astronomy and Astrophysics

The planetary system around the naked-eye star $\nu^2$ Lupi (HD 136352; TOI-2011) is composed of three exoplanets with masses of 4.7, 11.2, and 8.6 Earth masses. The TESS and CHEOPS missions revealed that all three planets are transiting and have radii straddling the radius gap separating volatile-rich and volatile-poor super-earths. Only a partial transit of planet d had been covered so we re-observed an inferior conjunction of the long-period 8.6 Earth-mass exoplanet $\nu^2$ Lup d with the CHEOPS space telescope. We confirmed its transiting nature by covering its whole 9.1 h transit for the first time. We refined the planet transit ephemeris to P = 107.1361 (+0.0019/-0.0022) days and Tc = 2,459,009.7759 (+0.0101/-0.0096) BJD_TDB, improving by ~40 times on the previously reported transit timing uncertainty. This refined ephemeris will enable further follow-up of this outstanding long-period transiting planet to search for atmospheric signatures or explore the planet's Hill sphere in search for an exomoon. In fact, the CHEOPS observations also cover the transit of a large fraction of the planet's Hill sphere, which is as large as the Earth's, opening the tantalising possibility of catching transiting exomoons. We conducted a search for exomoon signals in this single-epoch light curve but found no conclusive photometric signature of additional transiting bodies larger than Mars. Yet, only a sustained follow-up of $\nu^2$ Lup d transits will warrant a comprehensive search for a moon around this outstanding exoplanet., Comment: 16 pages, 9 figures. Accepted for publication in Astronomy & Astrophysics

Among the hundreds of known hot Jupiters (HJs), only five have been found to have companions on short-period orbits. Within this rare class of multiple planetary systems, the architecture of WASP-47 is unique, hosting an HJ (planet -b) with both an inner and an outer sub-Neptunian mass companion (-e and -d, respectively) as well as an additional non-transiting, long-period giant (-c). The small period ratio between planets -b and -d boosts the transit time variation (TTV) signal, making it possible to reliably measure the masses of these planets in synergy with the radial velocity (RV) technique. In this paper, we present new space- and ground-based photometric data of WASP-47b and WASP-47-d, including 11 unpublished light curves from the ESA mission CHEOPS. We analyzed the light curves in a homogeneous way together with all the publicly available data to carry out a global $N$-body dynamical modeling of the TTV and RV signals. We retrieved, among other parameters, a mass and density for planet -d of $M_\mathrm{d}=15.5\pm 0.8$ $M_\oplus$ and $\rho_\mathrm{d}=1.69\pm 0.22$ g\,cm$^{-3}$, which is in good agreement with the literature and consistent with a Neptune-like composition. For the inner planet (-e), we found a mass and density of $M_\mathrm{e}=9.0\pm 0.5$ $M_\oplus$ and $\rho_\mathrm{e}=8.1\pm 0.5$ g\,cm$^{-3}$, suggesting an Earth-like composition close to other ultra-hot planets at similar irradiation levels. Though this result is in agreement with previous RV+TTV studies, it is not in agreement with the most recent RV analysis (at 2.8$\sigma$), which yielded a lower density compatible with a pure silicate composition. This discrepancy highlights the still unresolved issue of suspected systematic offsets between RV and TTV measurements. In this paper, we also significantly improve the orbital ephemerides of all transiting planets, which will be crucial for any future follow-up., Comment: 17 pages, 8 figures, 10 tables, A&A in press. Typos corrected

We present a Zeeman-Doppler imaging study of two weak-line T Tauri stars TAP 4 and TAP 40, based on the high-resolution spectropolarimetric observations with ESPaDOnS at the Canada-France-Hawaii Telescope in November 2013, in the framework of the MaTYSSE large programme. We apply two Zeeman-Doppler imaging codes to the Stokes I and V profiles to reconstruct their brightness and large-scale magnetic field images. The results given by the two imaging codes are in good agreement with each other. TAP 4 shows a large polar cool spot and several intermediate-latitude warm spots on its surface, whereas TAP 40 exhibits very weak variations in its Stokes I profiles suggesting a mostly unspotted photosphere. We detect Zeeman signatures in the Stokes V profiles of both stars. The reconstructed magnetic maps reveal dominantly toroidal fields, which enclose about 60 per cent of the total magnetic energy for both of TAP 4 and TAP 40. Both stars show prominent circular ring features of the azimuthal magnetic field. We derive a solar-like surface differential rotation on TAP 4 from the tomographic modelling. The brightness image of TAP 4 is used to predict the radial velocity jitters induced by its activity. After filtering out the activity jitter, the RMS of its RVs is reduced from 1.7 km s$^{-1}$ to 0.2 km s$^{-1}$, but we do not detect any periodic signals in the filtered RVs of TAP 4, implying that it is unlikely to host a close-in exoplanet more massive than $\sim$3.5 M$_{\rm Jup}$ at 0.1 au., Comment: 10 pages, 11 figures, accepted for publication in MNRAS

We report the detection of a gas-giant planet in orbit around both stars of an eclipsing binary star system that also contains the smaller, inner transiting planet TOI-1338b. The new planet, called TOI-1338/BEBOP-1c, was discovered using radial-velocity data collected with the HARPS and ESPRESSO spectrographs. Our analysis reveals it is a $65.2~\rm{M_{\oplus}}$ circumbinary planet with a period of $215.5~$days. This is the first detection of a circumbinary planet using radial-velocity observations alone, and makes TOI-1338/BEBOP-1 only the second confirmed multiplanet circumbinary system to date. We do not detect the smaller inner transiting planet with radial-velocity data, and can place an upper limit on the inner planet's mass at $21.8~\mathrm{M}_\oplus$ with $99\%$ confidence. The inner planet is the first circumbinary planet amenable for atmospheric characterisation, using the James Webb Space Telescope., Comment: The Version of Record of this article is published in Nature Astronomy, and is available online at https://doi.org/10.1038/s41550-023-01948-4 Main: 14-pages, 4-Figures, 1-table. Supporting material: 48-pages, 12-figures, 4-tables