Search

The sub-Jovian desert is a region in the mass-period and radius-period parameter space, typically encompassing short-period ranges between super-Earths and hot Jupiters, that exhibits an intrinsic dearth of planets. This scarcity is likely shaped by photoevaporation caused by the stellar irradiation received by giant planets that have migrated inward. We report the detection and characterization of TOI-3568 b, a transiting super-Neptune with a mass of $26.4\pm1.0$ M$_\oplus$, a radius of $5.30\pm0.27$ R$_\oplus$, a bulk density of $0.98\pm0.15$ g cm$^{-3}$, and an orbital period of 4.417965(5) d situated in the vicinity of the sub-Jovian desert. This planet orbiting a K dwarf star with solar metallicity, was identified photometrically by TESS. It was characterized as a planet by our high-precision radial velocity monitoring program using MAROON-X at Gemini North, supplemented by additional observations from the SPICE large program with SPIRou at CFHT. We performed a Bayesian MCMC joint analysis of the TESS and ground-based photometry, MAROON-X and SPIRou radial velocities, to measure the orbit, radius, and mass of the planet, as well as a detailed analysis of the high-resolution flux and polarimetric spectra to determine the physical parameters and elemental abundances of the host star. Our results reveal TOI-3568 b as a hot super-Neptune, rich in hydrogen and helium with a core of heavier elements with a mass between 10 and 25 M$_\oplus$. We analyzed the photoevaporation status of TOI-3568 b and found that it experiences one of the highest EUV luminosities among planets with a mass M$_{\rm p}$ $<2$ M$_{\rm Nep}$, yet it has an evaporation lifetime exceeding 5 Gyr. Positioned in the transition between two significant populations of exoplanets on the mass-period and energy diagrams, this planet presents an opportunity to test theories concerning the origin of the sub-Jovian desert., Comment: Accepted for publication in A&A on September 4, 2024

This paper analyses spectropolarimetric observations of the classical T Tauri star (CTTS) GM Aurigae collected with SPIRou, the near-infrared spectropolarimeter at the Canada-France-Hawaii Telescope, as part of the SLS and SPICE Large Programs. We report for the first time results on the large-scale magnetic field at the surface of GM Aur using Zeeman Doppler imaging. Its large-scale magnetic field energy is almost entirely stored in an axisymmetric poloidal field, which places GM Aur close to other CTTSs with similar internal structures. A dipole of about 730 G dominates the large-scale field topology, while higher-order harmonics account for less than 30 per-cent of the total magnetic energy. Overall, we find that the main difference between our three reconstructed maps (corresponding to sequential epochs) comes from the evolving tilt of the magnetic dipole, likely generated by non-stationary dynamo processes operating in this largely convective star rotating with a period of about 6 d. Finally, we report a 5.5$\sigma$ detection of a signal in the activity-filtered radial velocity data of semi-amplitude 110 $\pm$ 20 m/s at a period of 8.745 $\pm$ 0.009 d. If attributed to a close-in planet in the inner accretion disc of GM Aur, it would imply that this planet candidate has a minimum mass of 1.10 $\pm$ 0.30 Mjup and orbits at a distance of 0.082 $\pm$ 0.002 au., Comment: Published in MNRAS. 24 pages, 15 figures

We report the discovery of two new exoplanet systems around fully convective stars, found from the radial-velocity (RV) variations of their host stars measured with the nIR spectropolarimeter CFHT/SPIRou over multiple years. GJ 3378 b is a planet with minimum mass of $5.26^{+0.94}_{-0.97}$ Mearth in an eccentric 24.73-day orbit around an M4V star of 0.26 Msun. GJ 1289 b has a minimum mass of $6.27\pm1.25$ Mearth in a 111.74-day orbit, in a circular orbit around an M4.5V star of mass 0.21 Msun. Both stars are in the solar neighbourhood, at respectively 7.73 and 8.86 pc. The low-amplitude RV signals are detected after line-by-line post-processing treatment. These potential sub-Neptune class planets around cool stars may have temperate atmospheres and be interesting nearby systems for further studies. We also recovered the large-scale magnetic field of both stars, found to be mostly axisymmetric and dipolar, and with a polar strength of 20-30 G and 200-240 G for GJ 3378 (in 2019-21) and GJ 1289 (in 2022-23), respectively. The rotation periods measured with the magnetic field differ from the orbital periods, and in general, stellar activity is not seen in the studied nIR RV time series of both stars. GJ 3378 b detection is not confirmed by optical RVs and is therefore considered a candidate at this point., Comment: accepted in A&A

For several years, the metastable helium triplet line has been successfully used as a tracer to probe atmospheric escape in transiting exoplanets. This absorption in the near-infrared (1083.3 nm) can be observed from the ground using high-resolution spectroscopy, providing new constraints on the mass-loss rate and the temperature characterizing the upper atmosphere of close-in exoplanets. The aim of this work is to search for the He triplet signature in 15 transiting exoplanets -- ranging from super-Earths to ultrahot Jupiters -- observed with SPIRou, a high-resolution (R~70 000) near-infrared spectropolarimeter at the CFHT, in order to bring new constraints or to improve existing ones regarding atmospheric escape through a homogeneous study. We developed a full data processing and analysis pipeline to correct for the residual telluric and stellar contributions. We then used two different 1D models based on the Parker-wind equations and nonlocal thermodynamic equilibrium (NLTE) radiative transfer to interpret the observational results. We confirm published He triplet detections for HAT-P-11 b, HD 189733 b, and WASP-69 b. We tentatively detect the signature of escaping He in HD 209458 b, GJ 3470 b, and WASP-76 b. We report new constraints on the mass-loss rate and temperature for our three detections and set upper limits for the tentative and nondetections. We notably report improved constraints on the mass-loss rate and temperature of the escaping gas for TOI-1807 b, and report a nondetection for the debated atmospheric escape in GJ 1214 b. We also conducted the first search for the He signature in GJ 486 b since its discovery and report a nondetection of the He triplet. Finally, we studied the impact of important model assumptions on our retrieved parameters, notably the limitations of 1D models and the influence of the H/He ratio on the derived constraints., Comment: 33 pages, 14 figures

In this paper we report near-infrared observations of the classical T Tauri star TW Hya with the SPIRou high-resolution spectropolarimeter and velocimeter at the 3.6-m Canada-France-Hawaii Telescope in 2019, 2020, 2021 and 2022. By applying Least-Squares Deconvolution (LSD) to our circularly polarized spectra, we derived longitudinal fields that vary from year to year from -200 to +100 G, and exhibit low-level modulation on the 3.6 d rotation period of TW Hya, despite the star being viewed almost pole-on. We then used Zeeman-Doppler Imaging to invert our sets of unpolarized and circularly-polarized LSD profiles into brightness and magnetic maps of TW Hya in all 4 seasons, and obtain that the large-scale field of this T Tauri star mainly consists of a 1.0-1.2 kG dipole tilted at about 20{\deg} to the rotation axis, whereas the small-scale field reaches strengths of up to 3-4 kG. We find that the large-scale field is strong enough to allow TW Hya to accrete material from the disc on the polar regions at the stellar surface in a more or less geometrically stable accretion pattern, but not to succeed in spinning down the star. We also report the discovery of a radial velocity signal of semi-amplitude $11.1^{+3.3}_{-2.6}$ m/s (detected at 4.3$\sigma$ at a period of 8.3 d in the spectrum of TW Hya, whose origin may be attributed to either a non-axisymmetric density structure in the inner accretion disc, or to a $0.55^{+0.17}_{-0.13}$ Jupiter mass candidate close-in planet (if orbiting in the disc plane), at an orbital distance of $0.075\pm0.001$ au., Comment: MNRAS, in press (23 pages, 16 figures, 6 tables)

Dynamo models of stellar magnetic fields for partly and fully convective stars are guided by observational constraints. Zeeman-Doppler imaging has revealed a variety of magnetic field geometries and, for fully convective stars in particular, a dichotomy: either strong, mostly axisymmetric, and dipole-dominated or weak, non-axisymmetric, and multipole-dominated. This dichotomy is explained by dynamo bistability or by long-term magnetic cycles, but there is no definite conclusion on the matter. We analysed optical spectropolarimetric data sets collected with ESPaDOnS and Narval between 2005 and 2016, and near-infrared SPIRou data obtained between 2019 and 2022 for three active M dwarfs with masses between 0.1 and 0.6 MSun: EV Lac, DS Leo, and CN Leo. We looked for changes in time series of longitudinal magnetic field, width of unpolarised mean-line profiles, and large-scale field topology as retrieved with principal component analysis and Zeeman-Doppler imaging. We retrieved pulsating (EV Lac), stable (DS Leo), and sine-like (CN Leo) long-term trends in longitudinal field. The width of near-infrared mean-line profiles exhibits rotational modulation only for DS Leo, whereas in the optical it is evident for both EV Lac and DS Leo. The line width variations are not necessarily correlated to those of the longitudinal field, suggesting complex relations between small- and large-scale field. We also recorded topological changes: a reduced axisymmetry for EV Lac and a transition from toroidal- to poloidal-dominated regime for DS Leo. For CN Leo, the topology remained dipolar and axisymmetric, with only an oscillation in field strength. Our results show a peculiar evolution of the magnetic field for each M dwarf, confirming that M dwarfs with distinct masses and rotation periods can undergo magnetic long-term variations, and suggesting a variety of cyclic behaviours of their magnetic fields., Comment: 50 pages, 26 figures, 12 tables, accepted in Astronomy & Astrophysics

Detecting planets within protoplanetary disks around young stars is essential for understanding planet formation and evolution. However, planet detection using the radial velocity method faces challenges due to strong stellar activity in these early stages. We aim to detect long-term periodicities in photometric and spectroscopic time series of the classical T Tauri star (CTTS) CI Tau, and retrieve evidence for inner embedded planets in its disk. The study conducted photometric and spectroscopic analyses using K2 and Las Cumbres Observatory Global Network light curves, and high-resolution spectra from ESPaDOnS and SPIRou. We focus our radial velocity analysis on a wavelength domain less affected by spot activity. To account for spot effects, a quasi-periodic Gaussian process model was applied to K2 light curve, ESPaDOnS, and SPIRou radial velocity data. Additionally, a detailed bisector analysis on cross-correlation functions was carried out to understand the cause of long-term periodicity. We detect coherent periods at $\sim$ 6.6 d, 9 d, $\sim$ 11.5 d, $\sim$ 14.2 d and $\sim$ 25.2 d, the latter is seen consistently across all datasets. Bisector analysis of the cross-correlation functions provides strong hints for combined activity-induced and Doppler reflex signal in the radial velocities at a period of 25.2 d. Our analysis suggests that this periodicity is best explained by the presence of a 3.6$\pm$0.3 M$_{Jup}$, eccentric (e$\sim$0.58) planet at a semi-major axis of 0.17 au. Our study outlines the difficulty of searching for disk-embedded planets in the inner 0.1 au's of young and active systems. We demonstrate that, when searching for planets in actively accreting stars such as CI Tau, the primary limitation is stellar activity rather than the precision of RV measurements provided by the instrument., Comment: 19 pages, 16 figures, Accepted to be published in Astronomy and Astrophysics

We report new observations of the classical T~Tauri star CI~Tau with the SPIRou near-infrared spectropolarimeter and velocimeter at the Canada-France-Hawaii Telescope (CFHT) in late 2019, 2020 and 2022, complemented with observations obtained with the ESPaDOnS optical spectropolarimeter at CFHT in late 2020. From our SPIRou and ESPaDOnS spectra, to which we applied Least-Squares Deconvolution, we infer longitudinal fields clearly modulated with the 9-d rotation period of CI~Tau. Using Zeeman-Doppler imaging, we reconstruct the large-scale magnetic topology, first from SPIRou data only in all three seasons, then from our 2020 SPIRou and ESPaDOnS data simultaneously. We find that CI~Tau hosts a mainly axisymmetric poloidal field, with a 1~kG dipole slightly tilted to the rotation axis and dark spots close to the pole that coincide with the footpoints of accretion funnels linking the star to the inner disc. Our results also suggest that CI~Tau accretes mass from the disc in a stable fashion. We further find that radial velocities (RV) derived from atomic and CO lines in SPIRou spectra are both rotationally modulated, but with a much lower amplitude than that expected from the putative candidate planet CI~Tau~b. We confirm the presence of a RV signal at a period of 23.86~d reported in a separate analysis, but detect it clearly in CO lines only and not in atomic lines, suggesting that it likely traces a non-axisymmetric structure in the inner disc of CI~Tau rather than a massive close-in planet., Comment: MNRAS, in press (22 pages, 11 figures, 6 tables)

We report the characterization of two planet candidates detected by the Transiting Exoplanet Survey Satellite (TESS), TOI-1199 b and TOI-1273 b, with periods of 3.7 and 4.6 days, respectively. Follow-up observations for both targets, which include several ground-based light curves, confirmed the transit events. High-precision radial velocities from the SOPHIE spectrograph revealed signals at the expected frequencies and phases of the transiting candidates and allowed mass determinations with a precision of $8.4\%$ and $6.7\%$ for TOI-1199 b and TOI-1273 b, respectively. The planetary and orbital parameters were derived from a joint analysis of the radial velocities and photometric data. We find that the planets have masses of $0.239\,\pm\,0.020\,M_{\mathrm{J}}$ and $0.222\,\pm\,0.015\,M_{\mathrm{J}}$ and radii of $0.938\,\pm\,0.025\,R_{\mathrm{J}}$ and $0.99\,\pm\,0.22\,R_{\mathrm{J}}$, respectively. The grazing transit of TOI-1273 b translates to a larger uncertainty in its radius, and hence also in its bulk density, compared to TOI-1199 b. The inferred bulk densities of $0.358\,\pm\,0.041\,\mathrm{g}\,\mathrm{cm}^{-3}$ and $0.28\,\pm\,0.11\,\mathrm{g}\,\mathrm{cm}^{-3}$ are among the lowest known for exoplanets in this mass range, which, considering the brightness of the host stars ($V \approx 11\,\mathrm{mag}$), render them particularly amenable to atmospheric characterization via the transit spectroscopy technique. The better constraints on the parameters of TOI-1199 b provide a transmission spectroscopy metric of $134\,\pm\,17$, making it the better suited of the two planets for atmospheric studies., Comment: Accepted for publication in A&A. 20 pages, 7 tables, and 15 figures

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

Transiting planets with orbital periods longer than 40 d are extremely rare among the 5000+ planets discovered so far. The lack of discoveries of this population poses a challenge to research into planetary demographics, formation, and evolution. Here, we present the detection and characterization of HD88986b, a potentially transiting sub-Neptune, possessing the longest orbital period among known transiting small planets (< 4 R$_{\oplus}$) with a precise mass measurement ($\sigma_M/M$ > 25%). Additionally, we identified the presence of a massive companion in a wider orbit around HD88986. Our analysis reveals that HD88986b, based on two potential single transits on sector 21 and sector 48 which are both consistent with the predicted transit time from the RV model, is potentially transiting. The joint analysis of RV and photometric data show that HD88986b has a radius of 2.49$\pm$0.18 R$_{\oplus}$, a mass of 17.2$^{+4.0}_{-3.8}$ M$_{\oplus}$, and it orbits every 146.05$^{+0.43}_{-0.40}$ d around a subgiant HD88986 which is one of the closest and brightest exoplanet host stars (G2V type, R=1.543 $\pm$0.065 R$_{\odot}$, V=$6.47\pm 0.01$ mag, distance=33.37$\pm$0.04 pc). The nature of the outer, massive companion is still to be confirmed; a joint analysis of RVs, Hipparcos, and Gaia astrometric data shows that with a 3$\sigma$ confidence interval, its semi-major axis is between 16.7 and 38.8 au and its mass is between 68 and 284 M$_{Jup}$. HD88986b's wide orbit suggests the planet did not undergo significant mass loss due to extreme-ultraviolet radiation from its host star. Therefore, it probably maintained its primordial composition, allowing us to probe its formation scenario. Furthermore, the cold nature of HD88986b (460$\pm$8 K), thanks to its long orbital period, will open up exciting opportunities for future studies of cold atmosphere composition characterization., Comment: 37 pages, accepted to be published in A&A

Planetary systems around solar analogs inform us about how planets form and evolve in Solar System-like environments. We report the detection and characterization of two planetary systems around the solar analogs TOI-1736 and TOI-2141 using TESS photometry data and spectroscopic data obtained with the SOPHIE instrument on the 1.93 m telescope at the Observatoire de Haute-Provence (OHP). We performed a detailed spectroscopic analysis of these systems to obtain the precise radial velocities (RV) and physical properties of their host stars. TOI-1736 and TOI-2141 each host a transiting sub-Neptune with radii of $2.44\pm0.18$ R$_{\oplus}$ and $3.05\pm0.23$ R$_{\oplus}$, orbital periods of $7.073088(7)$ d and $18.26157(6)$ d, and masses of $12.8\pm1.8$ M$_{\oplus}$ and $24\pm4$ M$_{\oplus}$, respectively. TOI-1736 shows long-term RV variations that are consistent with a two-planet solution plus a linear trend of $-0.177$ ms$^{-1}$d$^{-1}$. We measured an RV semi-amplitude of $201.1\pm0.7$ ms$^{-1}$ for the outer companion, TOI-1736 c, implying a projected mass of $m_{c}\sin{i}=8.09\pm0.20$ M$_{\rm Jup}$. From the GAIA DR3 astrometric excess noise, we constrained the mass of TOI-1736 c at $8.7^{+1.5}_{-0.6}$ M$_{\rm Jup}$. This planet is in an orbit of $570.2\pm0.6$ d with an eccentricity of $0.362\pm0.003$ and a semi-major axis of $1.381\pm0.017$ au, where it receives a flux of $0.71\pm0.08$ times the bolometric flux incident on Earth, making it an interesting case of a supergiant planet that has settled into an eccentric orbit in the habitable zone of a solar analog. Our analysis of the mass-radius relation for the transiting sub-Neptunes shows that both TOI-1736 b and TOI-2141 b likely have an Earth-like dense rocky core and a water-rich envelope., Comment: Accepted for publication in A&A on October 6, 2023

In this paper, we study six slowly rotating mid-to-late M~dwarfs (rotation period $P_{\mathrm{rot}} \approx 40-190\,\mathrm{dy}$) by analysing spectropolarimetric data collected with SPIRou at the Canada-France-Hawaii Telescope as part of the SPIRou Legacy Survey from 2019 to 2022. From $\approx$100--200 Least-Squares-Deconvolved (LSD) profiles of circularly polarised spectra of each star, we confirm the stellar rotation periods of the six M~dwarfs and explore their large-scale magnetic field topology and its evolution with time using both the method based on Principal Component Analysis (PCA) proposed recently and Zeeman-Doppler Imaging. All M~dwarfs show large-scale field variations on the time-scale of their rotation periods, directly seen from the circularly polarised LSD profiles using the PCA method. We detect a magnetic polarity reversal for the fully-convective M~dwarf GJ~1151, and a possible inversion in progress for Gl~905. The four fully-convective M~dwarfs of our small sample (Gl~905, GJ~1289, GJ~1151, GJ~1286) show a larger amount of temporal variations (mainly in field strength and axisymmetry) than the two partly-convective ones (Gl~617B, Gl~408). Surprisingly, the six M~dwarfs show large-scale field strengths in the range between 20 to 200\,G similar to those of M~dwarfs rotating significantly faster. Our findings imply that the large-scale fields of very slowly rotating M~dwarfs are likely generated through dynamo processes operating in a different regime than those of the faster rotators that have been magnetically characterized so far., Comment: 16 pages, 35 figures, accepted for publication in MNRAS

We present the results of an analysis aimed at probing the small-scale magnetic fields of M dwarfs observed with SPIRou, the nIR high-resolution spectro-polarimeter installed at the Canada-France-Hawaii Telescope, in the context of the SPIRou Legacy Survey. Our analysis relies on high-resolution median spectra built from several tens of spectra recorded between 2019 and 2022, and on synthetic spectra computed with the ZeeTurbo code for various combination of atmospheric parameters and magnetic field strengths. We pursue the efforts undertaken in a previous study and focus on 44 weakly to moderately active M dwarfs. We derive average magnetic field strengths (<$B$>) ranging from 0.05 to 1.15 kG, in good agreement with activity estimates and rotation periods. We found that including magnetic fields in our models has virtually no impact on our derived atmospheric parameters, and that a priori assumptions on the stellar surface gravity can affect our estimated <$B$>. Our results suggest that small-scale magnetic fields account for more than 70% of the overall average magnetic field for most targets whose large-scale fields were previously measured. We derived low magnetic fluxes for several targets in our sample, and found no clear evidence that <$B$> decreases with increasing Rossby number in the unsaturated dynamo regime. We even identified counterexamples (GJ 1289 and GJ 1286) where the small-scale field is unusually strong despite the long rotation period. Along with similar results on the large-scale fields, our findings further suggest that dynamo processes may operate in a non-conventional mode in these strongly magnetic, slowly-rotating stars., Comment: 27 pages, 25 figures, including appendix. Accepted for publication in MNRAS

Molecular bands of metal oxides and hydrides dominate the optical and near-infrared spectra of M dwarfs. High-resolution spectra of these bands have immense potential for determining many properties of these stars, such as effective temperature, surface gravity, elemental abundances, radial velocity, or surface magnetic fields. Techniques are being developed to do this but remain limited by the current availability and accuracy of molecular data and spectral line lists. This paper reports metal monohydride line lists selected from near-infrared and visible laboratory data to show that specific bands in several electronic transitions can be used to identify CrH, NiH, and FeH in M stars and to determine radial velocities from Doppler shifts. The possibility of measuring magnetic fields is also investigated for FeH and CrH. We used systematic cross-correlation analysis between unpolarised spectra from a selection of M stars and state-specific laboratory line lists. These lists were generated from a combination of existing data and new laboratory laser-excitation spectra recorded at Doppler-limited resolution, in zero-field conditions or in magnetic fields up to 0.6 tesla. Results. We show that transitions at visible wavelengths in FeH and NiH, usually neglected in the analysis of the spectra of M-type stars, do in fact contribute to the spectra, and we demonstrate the influence of magnetic sensitivity on selected transitions in CrH and FeH. Although the new line lists focus on transitions recorded at temperatures significantly lower than those of stellar objects, they remain pertinent because they cover some band-head regions of high spectral density. FeH bands can provide a useful supplement to atomic lines for the analysis of high-resolution optical and near-infrared spectra of M dwarfs. We demonstrate the influence of a magnetic field on CrH signatures around 862 nm., Comment: Accepted for publication in A&A, 10 pages, 12 figures, 2 tables, additional material will be made available at CDS/Vizier

We report results of a spectropolarimetric monitoring of the young Sun-like star V1298~Tau based on data collected with the near-infrared spectropolarimeter SPIRou at the Canada-France-Hawaii Telescope between late 2019 and early 2023. Using Zeeman-Doppler Imaging and the Time-dependent Imaging of Magnetic Stars methods on circularly polarized spectra, we reconstructed the large-scale magnetic topology of the star (and its temporal evolution), found to be mainly poloidal and axisymmetric with an average strength varying from 90 to 170 G over the ~3.5 years of monitoring. The magnetic field features a dipole whose strength evolves from 85 to 245 G, and whose inclination with respect to the stellar rotation axis remains stable until 2023 where we observe a sudden change, suggesting that the field may undergo a polarity reversal, potentially similar to those periodically experienced by the Sun. Our data suggest that the differential rotation shearing the surface of V1298 Tau is about 1.5 times stronger than that of the Sun. When coupling our data with previous photometric results from K2 and TESS and assuming circular orbits for all four planets, we report a $3.9\sigma$ detection of the radial velocity signature of the outermost planet (e), associated with a most probable mass, density and orbital period of $M_e=0.95^{+0.33}_{-0.24} \ \rm M_{\rm jup}$, $\rho_e=1.66^{+0.61}_{-0.48}$ $\rm g\,cm^{-3}$ and $P_e=53.0039\pm0.0001 \ \rm d$, respectively. For the 3 inner planets, we only derive 99\% confidence upper limits on their mass of $0.44\ \rm M_{\rm jup}$, $0.22\ \rm M_{\rm jup}$ and $0.25\ \rm M_{\rm jup}$, for b, c and d, respectively., Comment: 45 pages, 28 figures, accepted for publication in MNRAS

In a companion paper, we introduced a publicly-available pipeline to characterise exoplanet atmospheres through high-resolution spectroscopy. In this paper, we use this pipeline to study the biases and degeneracies that arise in atmospheric characterisation of exoplanets in near-infrared ground-based transmission spectroscopy. We inject synthetic planetary transits into sequences of SPIRou spectra of the well known M dwarf star Gl 15 A, and study the effects of different assumptions on the retrieval. We focus on (i) mass and radius uncertainties, (ii) non isothermal vertical profiles and (iii) identification and retrieval of multiple species. We show that the uncertainties on mass and radius should be accounted for in retrievals and that depth-dependent temperature information can be derived from high-resolution transmission spectroscopy data. Finally, we discuss the impact of selecting wavelength orders in the retrieval and the issues that arise when trying to identify a single species in a multi-species atmospheric model. This analysis allows us to understand better the results obtained through transmission spectroscopy and their limitations in preparation to the analysis of actual SPIRou data., Comment: Accepted in MNRAS

Atmospheric characterisation of exoplanets from the ground is an actively growing field of research. In this context we have created the ATMOSPHERIX consortium: a research project aimed at characterizing exoplanets atmospheres using ground-based high resolution spectroscopy. This paper presents the publicly-available data analysis pipeline and demonstrates the robustness of the recovered planetary parameters from synthetic data. Simulating planetary transits using synthetic transmission spectra of a hot Jupiter that were injected into real SPIRou observations of the non-transiting system Gl 15 A, we show that our pipeline is successful at recovering the planetary signal and input atmospheric parameters. We also introduce a deep learning algorithm to optimise data reduction which proves to be a reliable, alternative tool to the commonly used principal component analysis. We estimate the level of uncertainties and possible biases when retrieving parameters such as temperature and composition and hence the level of confidence in the case of retrieval from real data. Finally, we apply our pipeline onto two real transits of HD~189733 b observed with SPIRou and obtain similar results than in the literature. In summary, we have developed a publicly available and robust pipeline for the forthcoming studies of the targets to be observed in the framework of the ATMOSPHERIX consortium, which can easily be adapted to other high resolution instruments than SPIRou (e.g. VLT-CRIRES, MAROON-X, ELT-ANDES), Comment: Accepted in MNRAS

We report the discovery and characterisation of a giant transiting planet orbiting a nearby M3.5V dwarf (d = 80.4 pc, $G$ = 15.1 mag, $K$=11.2 mag, R$_\star$ = 0.358 $\pm$ 0.015 R$_\odot$, M$_\star$ = 0.340 $\pm$ 0.009 M$_\odot$). Using the photometric time series from TESS sectors 10, 36, 46, and 63 and near-infrared spectrophotometry from ExTrA, we measured a planetary radius of 0.77 $\pm$ 0.03 R$_J$ and an orbital period of 1.52 days. With high-resolution spectroscopy taken by the CFHT/SPIRou and ESO/ESPRESSO spectrographs, we refined the host star parameters ([Fe/H] = 0.27 $\pm$ 0.12) and measured the mass of the planet (0.273 $\pm$ 0.006 M$_J$). Based on these measurements, TOI-4860 b joins the small set of massive planets ($>$80 M$_E$) found around mid to late M dwarfs ($<$0.4 R$_\odot$), providing both an interesting challenge to planet formation theory and a favourable target for further atmospheric studies with transmission spectroscopy. We identified an additional signal in the radial velocity data that we attribute to an eccentric planet candidate ($e=0.66\pm0.09$) with an orbital period of $427\pm7$~days and a minimum mass of $1.66\pm 0.26$ M$_J$, but additional data would be needed to confirm this., Comment: 16 pages, 14 figures, accepted for publication in A&A

We present near-infrared spectropolarimetric observations of a sample of 43 weakly- to moderately-active M dwarfs, carried with SPIRou at the Canada-France-Hawaii Telescope in the framework of the SPIRou Legacy Survey from early 2019 to mid 2022. We use the 6700 circularly polarised spectra collected for this sample to investigate the longitudinal magnetic field and its temporal variations for all sample stars, from which we diagnose, through quasi-periodic Gaussian process regression, the periodic modulation and longer-term fluctuations of the longitudinal field. We detect the large-scale field for 40 of our 43 sample stars, and infer a reliable or tentative rotation period for 38 of them, using a Bayesian framework to diagnose the confidence level at which each rotation period is detected. We find rotation periods ranging from 14 to over 60d for the early-M dwarfs, and from 70 to 200d for most mid- and late-M dwarfs (potentially up to 430d for one of them). We also find that the strength of the detected large-scale fields does not decrease with increasing period or Rossby number for the slowly rotating dwarfs of our sample as it does for higher-mass, more active stars, suggesting that these magnetic fields may be generated through a different dynamo regime than those of more rapidly rotating stars. We also show that the large-scale fields of most sample stars evolve on long timescales, with some of them globally switching sign as stars progress on their putative magnetic cycles., Comment: MNRAS, in press (25 pages, 15 figures, 3 tables)

SPIRou is a near-infrared spectropolarimeter and a high-precision velocimeter. The SPIRou Legacy Survey collected data from February 2019 to June 2022, half of the time devoted to a blind search for exoplanets around nearby cool stars. The aim of this paper is to present this program and an overview of its properties, and to revisit the radial velocity (RV) data of two multiplanet systems, including new visits with SPIRou. From SPIRou data, we can extract precise RVs using efficient telluric correction and line-by-line measurement techniques, and we can reconstruct stellar magnetic fields from the collection of polarized spectra using the Zeeman-Doppler imaging method. The stellar sample of our blind search in the solar neighborhood, the observing strategy, the RV noise estimates, chromatic behavior, and current limitations of SPIRou RV measurements on bright M dwarfs are described. In addition, SPIRou data over a 2.5-year time span allow us to revisit the known multiplanet systems GJ~876 and GJ~1148. For GJ~876, the new dynamical analysis including the four planets is consistent with previous models and confirms that this system is deep in the Laplace resonance and likely chaotic. The large-scale magnetic map of GJ~876 over two consecutive observing seasons is obtained and shows a dominant dipolar field with a polar strength of 30~G, which defines the magnetic environment in which the inner planet with a period of 1.94~d is embedded. For GJ~1148, we refine the known two-planet model., Comment: accepted in A&A

One manifestation of dynamo action on the Sun is the 22-yr magnetic cycle, exhibiting a polarity reversal and a periodic conversion between poloidal and toroidal fields. For M dwarfs, several authors claim evidence of activity cycles from photometry and analyses of spectroscopic indices, but no clear polarity reversal has been identified from spectropolarimetric observations. Our aim is to monitor the evolution of the large-scale field of AD Leo, which has shown hints of a secular evolution from past dedicated spectropolarimetric campaigns. We analysed near-infrared spectropolarimetric observations of the active M dwarf AD Leo taken with SPIRou between 2019 and 2020 and archival optical data collected with ESPaDOnS and Narval between 2006 and 2019. We searched for long-term variability in the longitudinal field, the width of unpolarised Stokes profiles, the unsigned magnetic flux derived from Zeeman broadening, and the geometry of the large-scale magnetic field using both Zeeman-Doppler Imaging and Principal Component Analysis. We found evidence of a long-term evolution of the magnetic field, featuring a decrease in axisymmetry (from 99% to 60%). This is accompanied by a weakening of the longitudinal field (-300 to -50 G) and a correlated increase in the unsigned magnetic flux (2.8 to 3.6 kG). Likewise, the width of the mean profile computed with selected near-infrared lines manifests a long-term evolution corresponding to field strength changes over the full time series, but does not exhibit modulation with the stellar rotation of AD Leo in individual epochs. The large-scale magnetic field of AD Leo manifested first hints of a polarity reversal in late 2020 in the form of a substantially increased dipole obliquity, while the topology remained predominantly poloidal and dipolar. This suggests that low-mass M dwarfs with a dipole-dominated magnetic field can undergo magnetic cycles., Comment: 26 pages, 18 figures, 8 tables

The space environment in which planets are embedded depends mainly on the host star and impacts the evolution of the planetary atmosphere. The quiet M dwarf GJ 436 hosts a close-in hot Neptune which is known to feature a comet-like tail of hydrogen atoms escaped from its atmosphere due to energetic stellar irradiation. Understanding such star-planet interactions is essential to shed more light on planet formation and evolution theories, in particular the scarcity of Neptune-size planets below 3 d orbital period, also known as ``Neptune desert''. We aimed at characterising the stellar environment around GJ 436, which requires an accurate knowledge of the stellar magnetic field. The latter is studied efficiently with spectropolarimetry, since it is possible to recover the geometry of the large-scale magnetic field by applying tomographic inversion on time series of circularly polarised spectra. We used spectropolarimetric data collected in the optical domain with Narval in 2016 to compute the longitudinal magnetic field, examine its periodic content via Lomb-Scargle periodogram and Gaussian Process Regression analysis, and finally reconstruct the large-scale field configuration by means of Zeeman-Doppler Imaging. We found an average longitudinal field of -12 G and a stellar rotation period of 46.6 d using a Gaussian Process model and 40.1 d using Zeeman-Doppler Imaging, both consistent with the literature. The Lomb-Scargle analysis did not reveal any significant periodicity. The reconstructed large-scale magnetic field is predominantly poloidal, dipolar and axisymmetric, with a mean strength of 16 G. This is in agreement with magnetic topologies seen for other stars of similar spectral type and rotation rate., Comment: 9 pages, 8 figures, 3 tables

Our aim is to detect and characterise long-period companions around main sequence stars (spectral types late F to early M). We use the RV method to search for exoplanets around stars. The RV variations are measured with HARPS at the ESO 3.6 metre telescope. The true mass and inclination of our heavier companions are provided by astrometry, for which we use proper motions from Hipparcos and Gaia. Five Jupiter-mass exoplanets are reported to orbit HIP54597, BD-210397 (x2), HD74698, and HD94771 with 8.9 yr, 5.2 yr, 17.4 yr, 9.4 yr, and 5.9 yr orbits, and to have minimum masses of $2.01 \pm 0.03$, $0.7 \pm 0.1$, $2.4^{+1.5}_{-0.2}$, $0.40 \pm 0.06$, and $0.53 \pm 0.03 M_J$ respectively. HD74698 also hosts a highly irradiated Neptune in a 15 day orbit with a minimum mass of $0.07 \pm 0.01 M_J$. The mass of HIP54597 b can maximally increase by 10% - 30%, the minimum mass of HD74698 c is likely equal to its true mass, and BD-210397 c has a mass of $2.66^{+0.63}_{-0.32} M_J$. HD62364 hosts a brown dwarf with a true mass of $18.77^{+0.66}_{-0.63} M_J$ in an orbit of 14 yr. HD56380B, HD221638B, and HD33473C have minimum masses within the brown dwarf limits, in orbits of 8.9 yr, 16.6 yr, and 50 yr respectively; however, astrometric measurements reveal them to be stellar binaries, with masses of $375.3^{+8.6}_{-8.4}$, $110.0^{+3.9}_{-3.7}$, and $271.0^{+3.9}_{-3.8} M_J$. The orbits of the stellar binaries HD11938 and HD61383 are incomplete. The preliminary result for HD61383 is a 0.190 $M_{\odot}$ binary in a 39 yr orbit. The secondary of the binary system HD11938 has a mass of 0.33 $M_{\odot}$ - which is confirmed by a secondary peak in the CCF - and a preliminary period of 35 yr. The origin of the 3.0 yr RV signal of HD3964 is uncertain as it shows entanglement with the magnetic cycle of the star. We finally report one more star, HD11608, with a magnetic cycle that mimics a planetary signal., Comment: 22 pages, 32 figures, accepted for publication in A&A

Context: Recent advances in the development of precise radial velocity (RV) instruments in the near-infrared (nIR) domain, such as SPIRou, have facilitated the study of M-type stars to more effectively characterize planetary systems. However, the nIR presents unique challenges in exoplanet detection due to various sources of planet-independent signals which can result in systematic errors in the RV data. Aims: In order to address the challenges posed by the detection of exoplanetary systems around M-type stars using nIR observations, we introduce a new data-driven approach for correcting systematic errors in RV data. The effectiveness of this method is demonstrated through its application to the star GJ 251. Methods: Our proposed method, referred to as $\texttt{Wapiti}$ (Weighted principAl comPonent analysIs reconsTructIon), uses a dataset of per-line RV time-series generated by the line-by-line (LBL) algorithm and employs a weighted principal component analysis (wPCA) to reconstruct the original RV time-series. A multi-step process is employed to determine the appropriate number of components, with the ultimate goal of subtracting the wPCA reconstruction of the per-line RV time-series from the original data in order to correct systematic errors. Results: The application of $\texttt{Wapiti}$ to GJ 251 successfully eliminates spurious signals from the RV time-series and enables the first detection in the nIR of GJ 251b, a known temperate super-Earth with an orbital period of 14.2 days. This demonstrates that, even when systematics in SPIRou data are unidentified, it is still possible to effectively address them and fully realize the instrument's capability for exoplanet detection. Additionally, in contrast to the use of optical RVs, this detection did not require to filter out stellar activity, highlighting a key advantage of nIR RV measurements., Comment: Submitted to A&A. For the publicly available Wapiti code, see https://github.com/HkmMerwan/wapiti

In this paper we present an analysis of near-infrared spectropolarimetric and velocimetric data of the young M dwarf AU Mic, collected with SPIRou at the Canada-France-Hawaii telescope from 2019 to 2022, mostly within the SPIRou Legacy Survey. With these data, we study the large- and small-scale magnetic field of AU Mic, detected through the unpolarized and circularly-polarized Zeeman signatures of spectral lines. We find that both are modulated with the stellar rotation period (4.86 d), and evolve on a timescale of months under differential rotation and intrinsic variability. The small-scale field, estimated from the broadening of spectral lines, reaches $2.61\pm0.05$ kG. The large-scale field, inferred with Zeeman-Doppler imaging from Least-Squares Deconvolved profiles of circularly-polarized and unpolarized spectral lines, is mostly poloidal and axisymmetric, with an average intensity of $550\pm30$ G. We also find that surface differential rotation, as derived from the large-scale field, is $\simeq$30% weaker than that of the Sun. We detect the radial velocity (RV) signatures of transiting planets b and c, although dwarfed by activity, and put an upper limit on that of candidate planet d, putatively causing the transit-timing variations of b and c. We also report the detection of the RV signature of a new candidate planet (e) orbiting further out with a period of $33.39\pm0.10$ d, i.e., near the 4:1 resonance with b. The RV signature of e is detected at 6.5$\sigma$ while those of b and c show up at $\simeq$4$\sigma$, yielding masses of $10.2^{+3.9}_{-2.7}$ and $14.2^{+4.8}_{-3.5}$ Earth masses for b and c, and a minimum mass of $35.2^{+6.7}_{-5.4}$ Earth masses for e., Comment: MNRAS, in press (20 pages and 12 figures + 9 pages of supplementary material)

Context: The search for extrasolar planets around the nearest M dwarfs is a crucial step towards identifying the nearest Earth-like planets. One of the main challenges in this search is that M dwarfs can be magnetically active and stellar activity can produce radial velocity (RV) signals that could mimic those of a planet. Aims: We aim to investigate whether the 2.2 day period observed in optical RVs of the nearby active M dwarf star Gl 388 (AD Leo) is due to stellar activity or to a planet that corotates with the star as suggested in the past. Methods: We obtained quasi-simultaneous optical RVs of Gl 388 from 2019 to 2021 with SOPHIE (R~75k) at the OHP in France, and near-IR RV and Stokes V measurements with SPIRou at the CFHT (R~70k). Results: The SOPHIE RV time series displays a periodic signal with a 2.23+-0.01 day period and 23.6+-0.5 m/s amplitude, which is consistent with previous HARPS observations obtained in 2005-2006. The SPIRou RV time series is flat at 5 m/s rms and displays no periodic signals. RV signals of amplitude higher than 5.3 m/s at a period of 2.23 days can be excluded with a confidence level higher than 99%. Using the modulation of the longitudinal magnetic field (Bl) measured with SPIRou, we derive a stellar rotation period of 2.2305+-0.0016 days. Conclusions: SPIRou RV measurements provide solid evidence that the periodic variability of the optical RVs of Gl 388 is due to stellar activity rather than to a corotating planet. The magnetic activity nature of the optical RV signal is further confirmed by the modulation of Bl with the same period. The SPIRou campaign on Gl 388 demonstrates the power of near-IR RV to confirm or infirm planet candidates discovered in the optical around active stars. SPIRou observations reiterate how effective spectropolarimetry is at determining the stellar rotation period., Comment: 25 pages, 23 figures, Accepted by Astronomy and Astrophysics; [v2] version implementing A&A language editor suggestions; [v3] improved ascii characters for ADS metadata visualization; [v4] minor corrections added in proof

We report first results on a method aimed at simultaneously characterising atmospheric parameters and magnetic properties of M dwarfs from high-resolution nIR spectra recorded with SPIRou in the framework of the SPIRou Legacy Survey. Our analysis relies on fitting synthetic spectra computed from MARCS model atmospheres to selected spectral lines, both sensitive and insensitive to magnetic fields. We introduce a new code, $\texttt{ZeeTurbo}$, obtained by including the Zeeman effect and polarised radiative transfer capabilities to $\texttt{Turbospectrum}$. We compute a grid of synthetic spectra with $\texttt{ZeeTurbo}$ for different magnetic field strengths and develop a process to simultaneously constrain $T_{\rm eff}$, $\log{g}$, [M/H], [$\alpha$/Fe] and the average surface magnetic flux. In this paper, we present our approach and assess its performance using simulations, before applying it to six targets observed in the context of the SPIRou Legacy Survey (SLS), namely AU Mic, EV Lac, AD Leo, CN Leo, PM J18482+0741, and DS Leo. Our method allows us to retrieve atmospheric parameters in good agreement with the literature, and simultaneously yields surface magnetic fluxes in the range 2-4 kG with a typical precision of 0.05 kG, in agreement with literature estimates, and consistent with the saturated dynamo regime in which most of these stars are., Comment: 17 pages plus supplementary material. Accepted for publication in MNRAS

Context. The rotation period of stars is an important parameter along with mass, radius, effective temperature. It is an essential parameter for any radial velocity monitoring, as stellar activity can mimic the presence of a planet at the stellar rotation period. Several methods exist to measure it, including long sequences of photometric measurements or temporal series of stellar activity indicators. Aims. Here, we use the circular polarization in near-infrared spectral lines for a sample of 43 quiet M dwarfs and compare the measured rotation periods to those obtained with other methods. Methods. From Stokes V spectropolarimetric sequences observed with SPIRou at CFHT and the data processed with the APERO pipeline, we compute the least squares deconvolution profiles using different masks of atomic stellar lines with known Land\'e factor appropriate to the effective temperature of the star. We derive the longitudinal magnetic field to examine its possible variation along the 50 to 200 observations of each star. For determining the stellar rotation period, we apply a Gaussian process regression enabling us to determine the rotation period of stars with evolving longitudinal field. Results. Among the 43 stars of our sample, we were able to measure a rotation period for 27 stars. For 8 stars, the rotation period was previously unknown. We find a good agreement of our rotation periods with periods found in the literature based on photometry and activity indicators and confirm that near-infrared spectropolarimetry is an important tool to measure rotation periods, even for magnetically quiet stars. Furthermore, we compute ages for 20 stars of our sample using gyrochronology.

The stellar activity of M dwarfs is the main limitation for discovering and characterizing exoplanets orbiting them since it induces quasi-periodic RV variations. We aim to characterize the magnetic field and stellar activity of the early, moderately active, M dwarf Gl205 in the optical and nIR domains. We obtained high-precision quasi-simultaneous spectra in the optical and nIR with the SOPHIE spectrograph and SPIRou spectropolarimeter between 2019 and 2022. We computed the RVs from both instruments and the SPIRou Stokes V profiles. We used ZDI to map the large-scale magnetic field over the time span of the observations. We studied the temporal behavior of optical and nIR RVs and activity indicators with the Lomb-Scargle periodogram and a quasi-periodic GP regression. In the nIR, we studied the equivalent width of Al I, Ti I, K I, Fe I, and He I. We modeled the activity-induced RV jitter using a multi-dimensional GP regression with activity indicators as ancillary time series. The optical and nIR RVs have similar scatter but nIR shows a more complex temporal evolution. We observe an evolution of the magnetic field topology from a poloidal dipolar field in 2019 to a dominantly toroidal field in 2022. We measured a stellar rotation period of Prot=34.4$\pm$0.5 d in the longitudinal magnetic field. Using ZDI we measure the amount of latitudinal differential rotation (DR) shearing the stellar surface yielding rotation periods of Peq=32.0$\pm$1.8 d at the stellar equator and Ppol=45.5$\pm$0.3 d at the poles. We observed inconsistencies in the activity indicators' periodicities that could be explained by these DR values. The multi-dimensional GP modeling yields an RMS of the RV residuals down to the noise level of 3 m/s for both instruments, using as ancillary time series H$\alpha$ and the BIS in the optical, and the FWHM in the nIR., Comment: 41 pages, 24 figures. Accepted for publication in A&A. Improved quality of figures and reduced size of Appendix

We describe advances on a method designed to derive accurate parameters of M dwarfs. Our analysis consists in comparing high-resolution infrared spectra acquired with the near-infrared spectro-polarimeter SPIRou to synthetic spectra computed from MARCS model atmospheres, in order to derive the effective temperature ($T_{\rm eff}$), surface gravity ($\rm \log{g}$), metallicity ([M/H]) and alpha-enhancement ($\rm [\alpha/Fe]$) of 44 M dwarfs monitored within the SPIRou Legacy Survey (SLS). Relying on 12 of these stars, we calibrated our method by refining our selection of well modelled stellar lines, and adjusted the line list parameters to improve the fit when necessary. Our retrieved $T_{\rm eff}$, $\rm \log{g}$ and [M/H] are in good agreement with literature values, with dispersions of the order of 50 K in $T_{\rm eff}$ and 0.1 dex in $\rm \log{g}$ and [M/H]. We report that fitting $\rm [\alpha/Fe]$ has an impact on the derivation of the other stellar parameters, motivating us to extend our fitting procedure to this additional parameter. We find that our retrieved $\rm [\alpha/Fe]$ are compatible with those expected from empirical relations derived in other studies., Comment: 13 pages, 9 figures + appendix and supplementary material; accepted for publication in MNRAS

We investigate the exoplanet candidate TOI-1422b, which was discovered by the TESS space telescope around the high proper-motion G2V star TOI-1422 ($V=10.6$ mag), 155pc away, with the primary goal of confirming its planetary nature and characterising its properties. We monitored TOI-1422 with the HARPS-N spectrograph for 1.5 years to precisely quantify its radial velocity variation. The radial velocity measurements are analyzed jointly with TESS photometry and we also check for blended companions through high-spatial resolution images using the AstraLux instrument. We estimate that the parent star has a radius and a mass of $R^*=1.019_{-0.013}^{+0.014} R_{\odot}$, $M^*=0.981_{-0.065}^{+0.062} M_{\odot}$, respectively. Our analysis confirms the planetary nature of TOI-1422b and also suggests the presence of a Neptune-mass planet on a more distant orbit, the candidate TOI-1422c, which is not detected in TESS light curves. The inner planet, TOI-1422b, orbits on a period $P_{\rm b}=12.9972\pm0.0006$ days and has an equilibrium temperature $T_{\rm eq, b}=867\pm17$ K. With a radius of $R_{\rm b}=3.96^{+0.13}_{-0.11} R_{\oplus}$, a mass of $M_{\rm b}=9.0^{+2.3}_{-2.0} M_{\oplus}$ and, consequently, a density of $\rho_{\rm b}=0.795^{+0.290}_{-0.235}$ g cm$^{-3}$, it can be considered a warm Neptune-size planet. Compared to other exoplanets of similar mass range, TOI-1422b is among the most inflated ones and we expect this planet to have an extensive gaseous envelope that surrounds a core with a mass fraction around $10\%-25\%$ of the total mass of the planet. The outer non-transiting planet candidate, TOI-1422c, has an orbital period of $P_{\rm c}=29.29^{+0.21}_{-0.20}$ days, a minimum mass, $M_{\rm c}\sin{i}$, of $11.1^{+2.6}_{-2.3} M_{\oplus}$, an equilibrium temperature of $T_{\rm eq, c}=661\pm13$ K and, therefore, if confirmed, it could be considered as another warm Neptune., Comment: 22 pages, 20 figures, 6 tables. Accepted for publication in Astronomy & Astrophysics on July 7, 2022. Abstract abridged

We report the discovery and characterization of the transiting extrasolar planet TOI-1710$\:$b. It was first identified as a promising candidate by the Transiting Exoplanet Survey Satellite (TESS). Its planetary nature was then established with SOPHIE and HARPS-N spectroscopic observations via the radial-velocity method. The stellar parameters for the host star are derived from the spectra and a joint Markov chain Monte-Carlo (MCMC) adjustment of the spectral energy distribution and evolutionary tracks of TOI-1710. A joint MCMC analysis of the TESS light curve and the radial-velocity evolution allows us to determine the planetary system properties. From our analysis, TOI-1710$\:$b is found to be a massive warm super-Neptune ($M_{\rm p}=28.3\:\pm\:4.7\:{\rm M}_{\rm Earth}$ and $R_{\rm p}=5.34\:\pm\:0.11\:{\rm R}_{\rm Earth}$) orbiting a G5V dwarf star ($T_{\rm eff}=5665\pm~55\mathrm{K}$) on a nearly circular 24.3-day orbit ($e=0.16\:\pm\:0.08$). The orbital period of this planet is close to the estimated rotation period of its host star $P_{\rm rot}=22.5\pm2.0~\mathrm{days}$ and it has a low Keplerian semi-amplitude $K=6.4\pm1.0~\mathrm{m\:s^{-1}}$; we thus performed additional analyses to show the robustness of the retrieved planetary parameters. With a low bulk density of $1.03\pm0.23~\mathrm{g\:cm^{-3}}$ and orbiting a bright host star ($J=8.3$, $V=9.6$), TOI-1710$\:$b is one of the best targets in this mass-radius range (near the Neptunian desert) for atmospheric characterization via transmission spectroscopy, a key measurement in constraining planet formation and evolutionary models of sub-Jovian planets., Comment: 18 pages, 16 figures, 5 tables, A&A in press

Magnetic interactions between stars and close-in planets may lead to a detectable signal on the stellar disk. HD 189733 is one of the key exosystems thought to harbor magnetic interactions, which may have been detected in August 2013. We present a set of twelve wind models at that period, covering the possible coronal states and coronal topologies of HD 189733 at that time. We assess the power available for the magnetic interaction and predict its temporal modulation. By comparing the predicted signal with the observed signal, we find that some models could be compatible with an interpretation based on star-planet magnetic interactions. We also find that the observed signal can be explained only with a stretch-and-break interaction mechanism, while that the Alfv\'en wings scenario cannot deliver enough power. We finally demonstrate that the past observational cadence of HD 189733 leads to a detection rate of only between 12 to 23%, which could explain why star-planet interactions have been hard to detect in past campaigns. We conclude that the firm confirmation of their detection will require dedicated spectroscopic observations covering densely the orbital and rotation period, combined with scarcer spectropolarimetric observations to assess the concomitant large-scale magnetic topology of the star., Comment: 16 pages, 8 figures, 4 tables, accepted for publication in MNRAS

We report the detection and characterization of the transiting sub-Neptune TOI-1759 b, using photometric time-series from TESS and near infrared spectropolarimetric data from SPIRou on the CFHT. TOI-1759 b orbits a moderately active M0V star with an orbital period of $18.849975\pm0.000006$ d, and we measure a planetary radius and mass of $3.06\pm0.22$ R$_\oplus$ and $6.8\pm2.0$ M$_\oplus$. Radial velocities were extracted from the SPIRou spectra using both the CCF and the LBL methods, optimizing the velocity measurements in the near infrared domain. We analyzed the broadband SED of the star and the high-resolution SPIRou spectra to constrain the stellar parameters and thus improve the accuracy of the derived planet parameters. A LSD analysis of the SPIRou Stokes $V$ polarized spectra detects Zeeman signatures in TOI-1759. We model the rotational modulation of the magnetic stellar activity using a GP regression with a quasi-periodic covariance function, and find a rotation period of $35.65^{+0.17}_{-0.15}$ d. We reconstruct the large-scale surface magnetic field of the star using ZDI, which gives a predominantly poloidal field with a mean strength of $18\pm4$ G. Finally, we perform a joint Bayesian MCMC analysis of the TESS photometry and SPIRou RVs to optimally constrain the system parameters. At $0.1176\pm0.0013$ au from the star, the planet receives $6.4$ times the bolometric flux incident on Earth, and its equilibrium temperature is estimated at $433\pm14$ K. TOI-1759 b is a likely gas-dominated sub-Neptune with an expected high rate of photoevaporation. Therefore, it is an interesting target to search for neutral hydrogen escape, which may provide important constraints on the planetary formation mechanisms responsible for the observed sub-Neptune radius desert., Comment: Accepted for publication in the 10. Planets and planetary systems section of Astronomy & Astrophysics

We present the results of a study aiming at retrieving the fundamental parameters of M dwarfs from spectra secured with SPIRou, the near-infrared high-resolution spectropolarimeter installed at the Canada-France-Hawaii Telescope (CFHT), in the framework of the SPIRou Legacy Survey (SLS). Our study relies on comparing observed spectra with two grids of synthetic spectra, respectively computed from PHOENIX and MARCS model atmospheres, with the ultimate goal of optimizing the precision at which fundamental parameters can be determined. In this first step, we applied our technique to 12 inactive M dwarfs with effective temperatures ($T_{\rm eff}$) ranging from 3000 to 4000 K. We implemented a benchmark to carry out a comparison of the two models used in this study. We report that the choice of model has a significant impact on the results and may lead to discrepancies in the derived parameters of 30 K in $T_{\rm eff}$ and 0.05 dex to 0.10 dex in surface gravity ($\log{g}$) and metallicity ([M/H]), as well as systematic shifts of up to 50 K in $T_{\rm eff}$ and 0.4 dex $\log{g}$ and [M/H]. The analysis is performed on high signal-to-noise ratio template SPIRou spectra, averaged over multiple observations corrected from telluric absorption features and sky lines, using both a synthetic telluric transmission model and principal component analysis. With both models, we retrieve $T_{\rm eff}$ , $\log{g}$ and [M/H] estimates in good agreement with reference literature studies, with internal error bars of about 30 K, 0.05 dex and 0.1 dex, respectively., Comment: 20 pages, plus supplementary material, accepted for publication in MNRAS

We present the discovery and characterization of a transiting sub-Neptune orbiting with a 16.20 day period around a nearby (28 pc) and bright(V=8.37) K0V star HD207897 (TOI-1611). This discovery is based on photometric measurements from the Transiting Exoplanet Survey Satellite(TESS) mission and radial velocity (RV) observations from the SOPHIE, Automated Planet Finder (APF) and HIRES high precision spectrographs. We used EXOFASTv2 for simultaneously modeling the parameters of the planet and its host star, combining photometric and RV data to determine the planetary system parameters. We show that the planet has a radius of 2.50+/-0.08 RE and a mass of either 14.4+/-1.6 ME or 15.9+/-1.6 ME with nearly equal probability; the two solutions correspond to two possibilities for the stellar activity period. Hence, the density is either 5.1+/-0.7 g cm^-3 or 5.5^{+0.8}_{-0.7} g cm^-3, making it one of the relatively rare dense sub-Neptunes. The existence of such a dense planet at only 0.12 AU from its host star is unusual in the currently observed sub-Neptune (2

We present the discovery of two new transiting extrasolar planet candidates identified as TOI-1296.01 and TOI-1298.01 by the Transiting Exoplanet Survey Satellite (TESS). The planetary nature of these candidates has been secured with the SOPHIE high-precision spectrograph through the measurement of the companion's mass with the radial velocity method. Both planets are similar to Saturn in mass and have similar orbital periods of a few days. They, however, show discrepant radii and therefore different densities. The radius discrepancy might be explained by the different levels of irradiation by the host stars. The subgiant star TOI-1296 hosts a low-density planet with 1.2 RJup while the less luminous, lower-size star TOI-1298 hosts a much denser planet with a 0.84 RJup radius, resulting in bulk densities of 0.198 and 0.743 g/cm3, respectively. In addition, both stars are strongly enriched in heavy elements, having metallicities of 0.44 and 0.49 dex, respectively. The planet masses and orbital periods are 0.298 (0.039) MJup and 3.9443715 days for TOI-1296b, and 0.356 (0.032) MJup and 4.537164 days for TOI-1298b.

Context. We present precise radial-velocity measurements of five solar-type stars observed with the HARPS Echelle spectrograph mounted on the 3.6-m telescope in La Silla (ESO, Chile). With a time span of more than 10 years and a fairly dense sampling, the survey is sensitive to low mass planets down to super-Earths on orbital periods up to 100 days. Aims. Our goal was to search for planetary companions around the stars HD39194, HD93385, HD96700, HD154088, and HD189567 and use Bayesian model comparison to make an informed choice on the number of planets present in the systems based on the radial velocity observations. These findings will contribute to the pool of known exoplanets and better constrain their orbital parameters. Methods. A first analysis was performed using the DACE (Data & Analysis Center for Exoplanets) online tools to assess the activity level of the star and the potential planetary content of each system. We then used Bayesian model comparison on all targets to get a robust estimate of the number of planets per star. We did this using the nested sampling algorithm PolyChord. For some targets, we also compared different noise models to disentangle planetary signatures from stellar activity. Lastly, we ran an efficient MCMC (Markov chain Monte Carlo) algorithm for each target to get reliable estimates for the planets' orbital parameters. Results. We identify 12 planets within several multiplanet systems. These planets are all in the super-Earth and sub-Neptune mass regime with minimum masses ranging between 4 and 13 M$_\oplus$ and orbital periods between 5 and 103 days. Three of these planets are new, namely HD 93385 b, HD 96700 c, and HD 189567 c., Comment: 19 pages, 17 figures

Context. Due to their low transit probability, the long-period planets are, as a population, only partially probed by transit surveys. Radial velocity surveys thus have a key role to play, in particular for giant planets. Cold Jupiters induce a typical radial velocity semi-amplitude of 10m.s^{-1}, which is well within the reach of multiple instruments that have now been in operation for more than a decade. Aims. We take advantage of the ongoing radial velocity survey with the sophie high-resolution spectrograph, which continues the search started by its predecessor elodie to further characterize the cold Jupiter population. Methods. Analyzing the radial velocity data from six bright solar-like stars taken over a period of up to 15 years, we attempt the detection and confirmation of Keplerian signals. Results. We announce the discovery of six planets, one per system, with minimum masses in the range 2.99-8.3 Mjup and orbital periods between 200 days and 10 years. The data do not provide enough evidence to support the presence of additional planets in any of these systems. The analysis of stellar activity indicators confirms the planetary nature of the detected signals. Conclusions. These six planets belong to the cold and massive Jupiter population, and four of them populate its eccentric tail. In this respect, HD 80869 b stands out as having one of the most eccentric orbits, with an eccentricity of 0.862^{+0.028}_{-0.018}. These planets can thus help to better constrain the migration and evolution processes at play in the gas giant population. Furthermore, recent works presenting the correlation between small planets and cold Jupiters indicate that these systems are good candidates to search for small inner planets., Comment: 24 pages, 12 figures, Accepted for publication in Astronomy & Astrophysics

Distinguishing classes within substellar objects and understanding their formation and evolution need larger samples of substellar companions such as exoplanets, brown dwarfs, and low-mass stars. In this paper, we look for substellar companions using radial velocity surveys of FGK stars with the SOPHIE spectrograph at the Observatoire de Haute-Provence. We assign here the radial velocity variations of 27 stars to their orbital motion induced by low-mass companions. We also constrained their plane-of-the-sky motion using HIPPARCOS and Gaia Data Release 1 measurements, which constrain the true masses of some of these companions. We report the detection and characterization of six cool Jupiters, three brown dwarf candidates, and 16 low-mass stellar companions. We additionally update the orbital parameters of the low-mass star HD 8291 B, and we conclude that the radial velocity variations of HD 204277 are likely due to stellar activity despite resembling the signal of a giant planet. One of the new giant planets, BD+631405 b, adds to the population of highly eccentric cool Jupiters, and it is presently the most massive member. Two of the cool Jupiter systems also exhibit signatures of an additional outer companion. The orbital periods of the new companions span 30 days to 11.5 years, their masses 0.72 Jupiter mass to 0.61 Solar mass, and their eccentricities 0.04 to 0.88. These discoveries probe the diversity of substellar objects and low-mass stars, which will help constrain the models of their formation and evolution., Comment: 27 pages, 14 figures, 13 tables, Accepted in A&A

Classical T Tauri stars are young low-mass systems still accreting material from their disks. These systems are dynamic on timescales of hours to years. The observed variability can help us infer the physical processes that occur in the circumstellar environment. We aim at understanding the dynamics of the magnetic interaction between the star and the inner accretion disk in young stellar objects. We present the case of the young stellar system V2129 Oph, which is a well-known T Tauri star. We performed a time series analysis of this star using high-resolution spectroscopic data at optical and infrared wavelengths from CFHT/ESPaDOnS, ESO/HARPS and CFHT/SPIRou. The new data sets allowed us to characterize the accretion-ejection structure in this system and to investigate its evolution over a timescale of a decade via comparisons to previous observational data. We measure radial velocity variations and recover a stellar rotation period of 6.53d. However, we do not recover the stellar rotation period in the variability of various circumstellar lines, such as H$\alpha$ and H$\beta$ in the optical or HeI 1083nm and Pa$\beta$ in the infrared. Instead, we show that the optical and infrared line profile variations are consistent with a magnetospheric accretion scenario that shows variability with a period of about 6.0d, shorter than the stellar rotation period. Additionally, we find a period of 8.5d in H$\alpha$ and H$\beta$ lines, probably due to a structure located beyond the corotation radius, at a distance of 0.09au. We investigate whether this could be accounted for by a wind component, twisted or multiple accretion funnel flows, or an external disturbance in the inner disk. We conclude that the dynamics of the accretion-ejection process can vary significantly on a timescale of just a few years, presumably reflecting the evolving magnetic field topology at the stellar surface.

We present the multiple stellar systems observed within the SpHere INfrared survey for Exoplanet (SHINE). SHINE searched for substellar companions to young stars using high contrast imaging. Although stars with known stellar companions within SPHERE field of view (<5.5 arcsec) were removed from the original target list, we detected additional stellar companions to 78 of the 463 SHINE targets observed so far. 27% of the systems have three or more components. Given the heterogeneity of the sample in terms of observing conditions and strategy, tailored routines were used for data reduction and analysis, some of which were specifically designed for these data sets. We then combined SPHERE data with literature and archival ones, TESS light curves and Gaia parallaxes and proper motions, to characterise these systems as completely as possible. Combining all data, we were able to constrain the orbits of 25 systems. We carefully assessed the completeness of our sample for the separation range 50-500 mas (period range a few years - a few tens of years), taking into account the initial selection biases and recovering part of the systems excluded from the original list due to their multiplicity. This allowed us to compare the binary frequency for our sample with previous studies and highlight some interesting trends in the mass ratio and period distribution. We also found that, for the few objects for which such estimate was possible, the values of the masses derived from dynamical arguments were in good agreement with the model predictions. Stellar and orbital spins appear fairly well aligned for the 12 stars having enough data, which favour a disk fragmentation origin. Our results highlight the importance of combining different techniques when tackling complex problems such as the formation of binaries and show how large samples can be useful for more than one purpose., Comment: 59 pages, 11 tables, 18 figures

SPIRou is a near-infrared (nIR) spectropolarimeter at the CFHT, covering the YJHK nIR spectral bands ($980-2350\,\mathrm{nm}$). We describe the development and current status of the SPIRou wavelength calibration in order to obtain precise radial velocities (RVs) in the nIR. We make use of a UNe hollow-cathode lamp and a Fabry-P\'erot \'etalon to calibrate the pixel-wavelength correspondence for SPIRou. Different methods are developed for identifying the hollow-cathode lines, for calibrating the wavelength dependence of the Fabry-P\'erot cavity width, and for combining the two calibrators. The hollow-cathode spectra alone do not provide a sufficiently accurate wavelength solution to meet the design requirements of an internal error of $\mathrm{<0.45\,m\,s^{-1}}$, for an overall RV precision of $\mathrm{1\,m\,s^{-1}}$. However, the combination with the Fabry-P\'erot spectra allows for significant improvements, leading to an internal error of $\mathrm{\sim 0.15\,m\,s^{-1}}$. We examine the inter-night stability, intra-night stability, and impact on the stellar RVs of the wavelength solution., Comment: 12 pages, 9 figures. Accepted to A&A

Young stars interact with their accretion disk through their strong magnetosphere. We investigate the magnetospheric accretion process in the young stellar system DoAr 44. We monitored the system over several rotational cycles, combining high-resolution optical and near-IR spectropolarimetry with long-baseline near-IR interferometry and multicolor photometry. DoAr 44 is a young 1.2 solar mass star, moderately accreting from its disk, and seen at a low inclination. We derive a rotational period of 2.96 d from the system's light curve. Several optical and near-IR line profiles probing the accretion funnel flows and the accretion shock are modulated at the stellar rotation period. The most variable line profile, HeI 1083 nm, exhibits modulated redshifted wings a signature of accretion funnel flows, as well as deep blueshifted absorptions indicative of transient outflows. The Zeeman-Doppler analysis suggests the star hosts a mainly dipolar magnetic field, inclined by about 20 deg. onto the spin axis, with an intensity reaching about 800 G at the photosphere, and up to 2 +/- 0.8 kG close to the accretion shock. The magnetic field appears strong enough to disrupt the inner disk close to the corotation radius, at a distance of about 4.6 stellar radii (0.043 au). This supports the upper limit of 5 stellar radii (0.047 au) we derived for the size of the magnetosphere from long baseline interferometry. DoAr 44 is a pre-transitional disk system, exhibiting a 25-30 au gap in its circumstellar disk, with the inner and outer disks being misaligned. On a scale of 0.1 au or less, our results indicate that the system steadily accretes from its inner disk through its tilted dipolar magnetosphere. We conclude that in spite of a highly structured outer disk, perhaps the signature of ongoing planetary formation, the magnetospheric accretion process proceeds unimpeded at the star-disk interaction level., Comment: 18 pages, 22 figures, Astronomy & Astrophysics, in press

This paper presents an overview of SPIRou, the new-generation near-infrared spectropolarimeter / precision velocimeter recently installed on the 3.6-m Canada-France-Hawaii Telescope (CFHT). Starting from the two main science goals, namely the quest for planetary systems around nearby M dwarfs and the study of magnetized star / planet formation, we outline the instrument concept that was designed to efficiently address these forefront topics, and detail the in-lab and on-sky instrument performances measured throughout the intensive testing phase that SPIRou was submitted to before passing the final acceptance review in early 2019 and initiating science observations. With a central position among the newly started programmes, the SPIRou Legacy Survey (SLS) Large Programme was allocated 300 CFHT nights until at least mid 2022. We also briefly describe a few of the first results obtained in the various science topics that SPIRou started investigating, focusing in particular on planetary systems of nearby M dwarfs, transiting exoplanets and their atmospheres, magnetic fields of young stars, but also on alternate science goals like the atmospheres of M dwarfs and the Earth's atmosphere. We finally conclude on the essential role that SPIRou and the CFHT can play in coordination with forthcoming major facilities like the JWST, the ELTs, PLATO and ARIEL over the decade., Comment: MNRAS, in press (20 pages, 18 figures, 1 table)

The SHINE project is a 500-star survey performed with SPHERE on the VLT for the purpose of directly detecting new substellar companions and understanding their formation and early evolution. Here we present an initial statistical analysis for a subsample of 150 stars that are representative of the full SHINE sample. Our goal is to constrain the frequency of substellar companions with masses between 1 and 75 MJup and semimajor axes between 5 and 300 au. We adopt detection limits as a function of angular separation from the survey data for all stars converted into mass and projected orbital separation using the BEX-COND-hot evolutionary tracks and known distance to each system. Based on the results obtained for each star and on the 13 detections in the sample, we use a MCMC tool to compare our observations to two different types of models. The first is a parametric model based on observational constraints, and the second type are numerical models that combine advanced core accretion and gravitational instability planet population synthesis. Using the parametric model, we show that the frequencies of systems with at least one substellar companion are $23.0_{-9.7}^{+13.5}\%$, $5.8_{-2.8}^{+4.7}\%$, and $12.6_{-7.1}^{+12.9}\%$ for BA, FGK, and M stars, respectively. We also demonstrate that a planet-like formation pathway probably dominates the mass range from 1-75 MJup for companions around BA stars, while for M dwarfs, brown dwarf binaries dominate detections. In contrast, a combination of binary star-like and planet-like formation is required to best fit the observations for FGK stars. Using our population model and restricting our sample to FGK stars, we derive a frequency of $5.7_{-2.8}^{+3.8}\%$, consistent with predictions from the parametric model. More generally, the frequency values that we derive are in excellent agreement with values obtained in previous studies., Comment: 24 pages, 14 figures, 3 tables. Accepted for publication in A&A

We present high resolution near-infrared spectropolarimetric observations using the SPIRou instrument at CFHT during a transit of the recently detected young planet AU Mic b, with supporting spectroscopic data from iSHELL at IRTF. We detect Zeeman signatures in the Stokes V profiles, and measure a mean longitudinal magnetic field of $\overline{B}_\ell=46.3\pm0.7$~G. Rotationally modulated magnetic spots likely cause long-term variations of the field with a slope of $d{B_\ell}/dt=-108.7\pm7.7$~G/d. We apply the cross-correlation technique to measure line profiles and obtain radial velocities through CCF template matching. We find an empirical linear relationship between radial velocity and $B_\ell$, which allows us to estimate the radial velocity variations which stellar activity induces through rotational modulation of spots for the five hours of continuous monitoring of AU Mic with SPIRou. We model the corrected radial velocities for the classical Rossiter-McLaughlin effect, using MCMC to sample the posterior distribution of the model parameters. This analysis shows that the orbit of AU Mic b is prograde and aligned with the stellar rotation axis with a sky-projected spin-orbit obliquity of $\lambda=0^{+18}_{-15}$ degrees. The aligned orbit of AU Mic b indicates that it formed in the protoplanetary disk that evolved to the current debris disk around AU Mic., Comment: Manuscript accepted for publication on 07/08/2020 as a Letter in Astronomy and Astrophysics, section 1

Context. 55 Cancri hosts five known exoplanets, most notably the hot super-Earth 55 Cnc e, which is one of the hottest known transiting super-Earths. Aims. Due to the short orbital separation and host star brightness, 55 Cnc e provides one of the best opportunities for studying star-planet interactions (SPIs). We aim to understand possible SPIs in this system, which requires a detailed understanding of the stellar magnetic field and wind impinging on the planet. Methods. Using spectropolarimetric observations, and Zeeman Doppler Imaging, we derive a map of the large-scale stellar magnetic field. We then simulate the stellar wind starting from the magnetic field map, using a 3D MHD model. Results. The map of the large-scale stellar magnetic field we derive has an average strength of 3.4 G. The field has a mostly dipolar geometry, with the dipole tilted by 90 degrees with respect to the rotation axis, and dipolar strength of 5.8 G at the magnetic pole. The wind simulations based on this magnetic geometry lead us to conclude that 55 Cnc e orbits inside the Alfv\'en surface of the stellar wind, implying that effects from the planet on the wind can propagate back to the stellar surface and result in SPI., Comment: Accepted for publication in A&A, 6 pages, 5 figures

Since 2011, the SOPHIE spectrograph has been used to search for Neptunes and super-Earths in the Northern Hemisphere. As part of this observational program, 290 radial velocity measurements of the 6.4 V magnitude star HD 158259 were obtained. Additionally, TESS photometric measurements of this target are available. We present an analysis of the SOPHIE data and compare our results with the output of the TESS pipeline. The radial velocity data, ancillary spectroscopic indices, and ground-based photometric measurements were analyzed with classical and $\ell_1$ periodograms. The stellar activity was modeled as a correlated Gaussian noise and its impact on the planet detection was measured with a new technique. The SOPHIE data support the detection of five planets, each with $m \sin i \approx 6 M_\oplus$, orbiting HD 158259 in 3.4, 5.2, 7.9, 12, and 17.4 days. Though a planetary origin is strongly favored, the 17.4 d signal is classified as a planet candidate due to a slightly lower statistical significance and to its proximity to the expected stellar rotation period. The data also present low frequency variations, most likely originating from a magnetic cycle and instrument systematics. Furthermore, the TESS pipeline reports a significant signal at 2.17 days corresponding to a planet of radius $\approx 1.2 R_\oplus$. A compatible signal is seen in the radial velocities, which confirms the detection of an additional planet and yields a $\approx 2 M_\oplus$ mass estimate. We find a system of five planets and a strong candidate near a 3:2 mean motion resonance chain orbiting HD 158259. The planets are found to be outside of the two and three body resonances., Comment: Accepted for publication in Astronomy & Astrophysics