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

This paper presents near-infrared spectropolarimetric and velocimetric observations of the young planet-hosting T Tauri star PDS 70, collected with SPIRou at the 3.6m Canada-France-Hawaii Telescope from 2020 to 2024. Clear Zeeman signatures from magnetic fields at the surface of PDS 70 are detected in our data set of 40 circularly polarized spectra. Longitudinal fields inferred from Zeeman signatures, ranging from -116 to 176 G, are modulated on a timescale of 3.008$\pm$0.006 d, confirming that this is the rotation period of PDS 70. Applying Zeeman-Doppler imaging to subsets of unpolarized and circularly polarised line profiles, we show that PDS 70 hosts low-contrast brightness spots and a large-scale magnetic field in its photosphere, featuring in particular a dipole component of strength 200-420 G that evolves on a timescale of months. From the broadening of spectral lines, we also infer that PDS 70 hosts a small-scale field of 2.51$\pm$0.12 kG. Radial velocities derived from unpolarized line profiles are rotationally modulated as well, and exhibit additional longer-term chromatic variability, most likely attributable to magnetic activity rather than to a close-in giant planet (with a 3sigma upper limit on its minimum mass of ~4 Mjup at a distance of ~0.2 au). We finally confirm that accretion occurs at the surface of PDS 70, generating modulated red-shifted absorption in the 1083.3-nm He i triplet, and show that the large-scale magnetic field, often strong enough to disrupt the inner accretion disc up to the corotation radius, weakens as the star gets fainter and redder (as in 2022), suggesting that dust from the disc more easily penetrates the stellar magnetosphere in such phases., Comment: MNRAS in press (20 pages, 12 figures, 5 tables)

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 present in this paper an analysis of near-infrared observations of the 0.3-Msun protostar V347 Aur collected with the SPIRou high-resolution spectropolarimeter and velocimeter at the 3.6-m Canada-France-Hawaii Telescope from October 2019 to April 2023. From a set of 79 unpolarized and circularly polarized spectra of V347 Aur to which we applied Least-Squares Deconvolution (LSD), we derived radial velocities and longitudinal fields, along with their temporal variations over our monitoring campaign of 1258 d. Our data show that V347 Aur is an eccentric binary system with an orbital period 154.6$\pm$0.7 d, experiencing strong to extreme accretion events near periastron. The companion is a 29.0$\pm$1.6-Mjup brown dwarf, a rare member of the brown dwarf desert of close companions around M dwarfs. We detect weak longitudinal fields (<100 G) at the surface of V347 Aur, significantly weaker than those of more evolved prototypical T Tauri stars. These fields show small-amplitude rotational modulation, indicating a mainly axisymmetric parent large-scale magnetic topology, and larger fluctuations at half the orbital period, suggesting that what we dub a ``pulsed dynamo'' triggered by orbital motion and pulsed accretion operates in V347 Aur. Applying Zeeman-Doppler imaging to our circularly polarised LSD profiles, we find that the large-scale field of V347 Aur is mainly toroidal for most of our observations, with the toroidal component switching sign near periastron and apoastron. The weak large-scale dipole (~30 G) is not able to disrupt the disc beyond 1.3 Rstar even at the lowest accretion rates, implying longitudinally distributed (rather than localized) accretion at the surface of the protostar., Comment: MNRAS in press (18 pages, 12 figures, 5 tables)

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

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)

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

The M dwarf star GJ 436 hosts a warm-Neptune that is losing substantial amount of atmosphere, which is then shaped by the interactions with the wind of the host star. The stellar wind is formed by particles and magnetic fields that shape the exo-space weather around the exoplanet GJ 436 b. Here, we use the recently published magnetic map of GJ 436 to model its 3D Alfv\'en-wave driven wind. By comparing our results with previous transmission spectroscopic models and measurements of non-thermal velocities at the transition region of GJ 436, our models indicate that the wind of GJ 436 is powered by a smaller flux of Alfv\'en waves than that powering the wind of the Sun. This suggests that the canonical flux of Alfv\'en waves assumed in solar wind models might not be applicable to the winds of old M dwarf stars. Compared to the solar wind, GJ 436's wind has a weaker acceleration and an extended sub-Alfv\'enic region. This is important because it places the orbit of GJ 436 b inside the region dominated by the stellar magnetic field (i.e., inside the Alfv\'en surface). Due to the sub-Alfv\'enic motion of the planet through the stellar wind, magnetohydrodynamic waves and particles released in reconnection events can travel along the magnetic field lines towards the star, which could power the anomalous ultraviolet flare distribution recently observed in the system. For an assumed planetary magnetic field of $B_p \simeq 2$ G, we derive the power released by stellar wind-planet interactions as $\mathcal{P} \sim 10^{22}$ -- $10^{23}$ erg s$^{-1}$, which is consistent with the upper limit of $10^{26}$ erg s$^{-1}$ derived from ultraviolet lines. We further highlight that, because star-planet interactions depend on stellar wind properties, observations that probe these interactions and the magnetic map used in 3D stellar wind simulations should be contemporaneous for deriving realistic results., Comment: 12 pages, 7 figures, 1 table; accepted for publication in A&A

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

The metastable helium triplet in the near-infrared (10833{\AA}) is among the most important probes of exoplanet atmospheres. It can trace their extended outer layers and constrain mass-loss. We use the near-infrared high-resolution spectropolarimeter SPIRou on the CFHT to search for the spectrally resolved helium triplet in the atmospheres of eleven exoplanets, ranging from warm mini-Neptunes to hot Jupiters and orbiting G, K, and M dwarfs. Observations were obtained as part of the SPIRou Legacy Survey and complementary open-time programs. We apply a homogeneous data reduction to all datasets and set constraints on the presence of metastable helium, despite the presence of systematics in the data. We confirm published detections for HAT-P-11b, HD189733b, and WASP-69b and set upper limits for the other planets. We apply the p-winds open source code to set upper limits on the mass-loss rate for the non-detections and to constrain the thermosphere temperature, mass-loss rate, line-of-sight velocity, and the altitude of the thermosphere for the detections. We confirm that the presence of metastable helium correlates with the stellar mass and the XUV flux received by the planets. We investigated the correlation between the mass-loss rate and the presence of metastable helium, but it remains difficult to draw definitive conclusions. Finally, some of our results are in contradiction with previous results in the literature, therefore we stress the importance of repeatable, homogeneous, and larger-scale analyses of the helium triplet to obtain robust statistics, study temporal variability, and better understand how the helium triplet can be used to explore the evolution of exoplanets., Comment: 28 pages, 13 figures, Accepted in A&A for publication

A comprehensive multi-wavelength campaign has been carried out to probe stellar activity and variability in the young Sun-like star $\iota$-Horologii. We present the results from long-term spectropolarimetric monitoring of the system by using the ultra-stable spectropolarimeter/velocimeter HARPS at the ESO 3.6-m telescope. Additionally, we included high-precision photometry from the NASA Transiting Exoplanet Survey Satellite (TESS) and observations in the far- and near-ultraviolet spectral regions using the STIS instrument on the NASA/ESA Hubble Space Telescope (HST). The high-quality dataset allows a robust characterisation of the star's rotation period, as well as a probe of the variability using a range of spectroscopic and photometric activity proxies. By analyzing the gradient of the power spectra (GPS) of the TESS lightcurves we constrained the faculae-to-spot driver ratio ($\rm S_{fac}/S_{spot}$) to 0.510$\pm$0.023, which indicates that the stellar surface is spot dominated during the time of the observations. We compared the photospheric activity properties derived from the GPS method with a magnetic field map of the star derived using Zeeman-Doppler imaging (ZDI) from simultaneous spectropolarimetric data for the first time. Different stellar activity proxies enable a more complete interpretation of the observed variability. For example, we observed enhanced emission in the HST transition line diagnostics C IV and C III, suggesting a flaring event. From the analysis of TESS data acquired simultaneously with the HST data, we investigate the photometric variability at the precise moment that the emission increased and derive correlations between different observables, probing the star from its photosphere to its corona.

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

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.

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 investigate the structure and dynamics of the magnetospheric accretion region and associated outflows on a scale smaller than 0.1 au around the young transitional disk system GM Aur. We monitored the variability of the system on timescales ranging from days to months, using high-resolution optical and near-infrared spectroscopy, multiwavelength photometry, and low-resolution near-infrared spectroscopy, over a total duration of six months (30 rotational cycles). We analyzed the photometric and line profile variability to characterize the accretion and ejection processes. The luminosity of the system is modulated by surface spots at the stellar rotation period of 6.04 days. The Balmer, Paschen, and Brackett hydrogen lines as well as the HeI 5876 A and HeI 10830 A line profiles are modulated on the same period. The PaB line flux correlates with the photometric excess in the u' band, which suggests that most of the line emission originates from the accretion process. High-velocity redshifted absorptions reaching below the continuum periodically appear in the near-infrared line profiles at the rotational phase in which the veiling and line fluxes are the largest. These are signatures of a stable accretion funnel flow and associated accretion shock at the stellar surface. This large-scale magnetospheric accretion structure appears fairly stable over at least 15 and possibly up to 30 rotational periods. In contrast, outflow signatures randomly appear as blueshifted absorption components in the Balmer and HeI 10830 A line profiles and disappear on a timescale of a few days. The coexistence of a stable, large-scale accretion pattern and episodic outflows supports magnetospheric ejections as the main process occurring at the star-disk interface. Stable magnetospheric accretion and episodic outflows appear to be physically linked on a scale of a few stellar radii in this system., Comment: 30 pages, 28 figures, 12 tables, accepted for publication in Astronomy & Astrophysics

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

Veiling is ubiquitous at different wavelength ranges in accreting stars. However, the origin of the veiling in the IR domain is not well understood. The accretion spot alone is not enough to explain the shallow photospheric IR lines in accreting systems, suggesting that another source is contributing to the veiling in the NIR. The inner disk is often quoted as the additional emitting source meant to explain the IR veiling. In this work, we aim to measure and discuss the NIR veiling to understand its origins and variability timescale, using a sample of 14 accreting stars observed with the CFHT/SPIRou spectrograph, within the framework of the SPIRou Legacy Survey. We compared the veiling measurements with accretion and inner disk diagnostics. The measured veiling grows from the Y to the K band for most of the targets in our sample. The IR veiling agrees with NIR emission excess obtained using photometric data. However, we also find a linear correlation between the veiling and the accretion properties of the system, showing that accretion contributes to the inner disk heating and, consequently, to the inner disk emission excess. We also show a connection between the NIR veiling and the system's inclination with respect to our line of sight. This is probably due to the reduction of the visible part of the inner disk edge, where the NIR emission excess is expected to arise, as the inclination of the system increases. The NIR veiling appears variable on a timescale of a day, showing the night-by-night dynamics of the optical veiling variability. In the long term, the mean NIR veiling seems to be stable for most of the targets on timescales of a month to a few years. However, during occasional episodes of high accretion, which affect the system's dynamic, the veiling also seems to be much more prominent at such times, as we found in the case of the target RU Lup.

Misalignments between a forming star's rotation axis and its outer disk axis, although not predicted by standard theories of stellar formation, have been observed in several classical T Tauri stars (cTTs). The low-mass cTTs DK Tau is suspected of being among them. It is also an excellent subject to investigate the interaction between stellar magnetic fields and material accreting from the circumstellar disk, as it presents clear signatures of accretion. The goal of this paper is to study DK Tau's average line-of-sight magnetic field (Blos) in both photospheric absorption lines and emission lines linked to accretion, using spectropolarimetric observations, as well as to examine inconsistencies regarding its rotation axis. We used data collected with the ESPaDOnS and NARVAL spectropolarimeters, probing two distinct epochs (2010 and 2012). We first determined the stellar parameters, such as effective temperature and v sin i. Next, we removed the effect of veiling from the spectra, then obtained least-squares deconvolution profiles of the absorption lines, before determining the Blos. We also investigated emission lines, the 587.6 nm HeI line and the CaII infrared triplet, as tracers of the magnetic fields present in the accretion shocks. We find that DK Tau experiences accretion onto a magnetic pole at an angle of about 30 degrees from the pole of its rotation axis, with a positive field at the base of the accretion funnels. In 2010 we find a magnetic field of up to 1.77kG, and in 2012 up to 1.99kG. Additionally, using our derived values of period, v sin i and stellar radius, we find a value of 58 degrees (+18)(-11) for the inclination of the stellar rotation axis, which is significantly different from the outer disk axis inclination of 21 degrees given in the literature. We find that DK Tau's outer disk axis is likely misaligned compared to its rotation axis by 37 degrees.

TOI-1695 is a V-mag=13 M-dwarf star from the northern hemisphere at 45$\,$pc from the Sun, around which a 3.134-day periodic transit signal from a super-Earth candidate was identified in TESS photometry. With a transit depth of 1.3$\,$mmag, the radius of candidate TOI-1695.01 was estimated by the TESS pipeline to be 1.82$\,$R$_\oplus$ with an equilibrium temperature of $\sim 620\,$K. We successfully detect a reflex motion of the star and establish it is due to a planetary companion at an orbital period consistent with the photometric transit period thanks to a year-long radial-velocity monitoring of TOI-1695 by the SPIRou infrared spectropolarimeter. We use and compare different methods to reduce and analyse those data. We report a 5.5-$\sigma$ detection of the planetary signal, giving a mass of $5.5 \pm 1.0\,$M$_\oplus$ and a radius of $2.03 \pm 0.18\,$R$_\oplus$. We derive a mean equilibrium planet temperature of $590 \pm 90\,$K. The mean density of this small planet of $3.6 \pm 1.1\,$g$\,$cm$^{-3}$ is similar (1.7-$\sigma$ lower) than that of the Earth. It leads to a non-negligible fraction of volatiles in its atmosphere with $f_{H,He}=0.28^{+0.46}_{-0.23}$% or $f_\text{water}=23 \pm 12$%. TOI-1695$\,$b is a new sub-Neptune planet at the border of the M-dwarf radius valley that can help test formation scenarios for super-Earth/sub-Neptune-like planets., Comment: 15 pages (+9 appendix pages), 19 figures, 8 tables, accepted for publication in A&A

Observations of surface magnetic fields of cool stars reveal a large diversity of configurations. Although there is now a consensus that these fields are generated through dynamo processes occurring within the convective zone, the physical mechanism driving such a variety of field topologies is still debated. This paper discusses the possible origins of dipole and multipole-dominated morphologies using three-dimensional numerical simulations of stratified systems where the magnetic feedback on the fluid motion is significant. Our main result is that dipolar solutions are found at Rossby numbers up to 0.4 in strongly stratified simulations, where previous works suggested that only multipolar fields should exist. We argue that these simulations are reminiscent of the outlier stars observed at Rossby numbers larger than 0.1, whose large-scale magnetic field is dominated by their axisymmetric poloidal component. As suggested in previous Boussinesq calculations, the relative importance of inertial over Lorentz forces is again controlling the dipolar to multipolar transition. Alternatively, we find that the ratio of kinetic to magnetic energies can equally well capture the transition in the field morphology. We test the ability of this new proxy to predict the magnetic morphology of a few M-dwarf stars whose internal structure matches that of our simulations and for which homogeneous magnetic field characterization is available. Finally, the magnitude of the differential rotation obtained in our simulations is compared to actual measurements reported in the literature for M-dwarfs. In our simulations, we find a clear relationship between anti-solar differential rotation and the emergence of dipolar fields., Comment: 15 pages, 15 figures, accepted for publication in MNRAS

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

Insights on stellar surface large-scale magnetic field topologies are usually drawn by applying Zeeman-Doppler-Imaging (ZDI) to the observed spectropolarimetric time series. However, ZDI requires experience for reliable results to be reached and is based on a number of prior assumptions that may not be valid, e.g., when the magnetic topology is evolving on timescales comparable to or shorter than the time span over which observations are collected. In this paper, we present a method based on Principal Component Analysis (PCA) applied to circularly polarised (Stokes~$V$) line profiles of magnetic stars to retrieve the main characteristics of the parent large-scale magnetic topologies, like for instance, the relative strength of the poloidal and toroidal components, and the degree of axisymmetry of the dominant field component and its complexity (dipolar or more complex). We show that this method can also be used to diagnose the temporal variability of the large-scale magnetic field. Performing best for stars with moderate projected equatorial velocities hosting relatively simple magnetic field topologies, this new method is simpler than ZDI, making it convenient to rapidly diagnose the main characteristics of the large-scale fields of non-degenerate stars and to provide insights into the temporal evolution of the field topology., Comment: 13 pages, 14 figures, accepted for publication in MNRAS

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

A detailed investigation of the magnetic properties of young Sun-like stars can provide valuable information on our Sun's magnetic past and its impact on the early Earth. We determine the properties of the moderately rotating young Sun-like star kappa Ceti's magnetic and activity cycles using 50 years of chromospheric activity data and six epochs of spectropolarimetric observations. The chromospheric activity was determined by measuring the flux in the Ca II H and K lines. A generalised Lomb-Scargle periodogram and a wavelet decomposition were used on the chromospheric activity data to establish the associated periodicities. The vector magnetic field of the star was reconstructed using the technique of Zeeman Doppler imaging on the spectropolarimetric observations. Our period analysis algorithms detect a 3.1 year chromospheric cycle in addition to the star's well-known ~6 year cycle period. Although the two cycle periods have an approximate 1:2 ratio, they exhibit an unusual temporal evolution. Additionally, the spectropolarimetric data analysis shows polarity reversals of the star's large-scale magnetic field, suggesting a ~10 year magnetic or Hale cycle. The unusual evolution of the star's chromospheric cycles and their lack of a direct correlation with the magnetic cycle establishes kappa Ceti as a curious young Sun. Such complex evolution of magnetic activity could be synonymous with moderately active young Suns, which is an evolutionary path that our own Sun could have taken., Comment: 15 pages including appendices, accepted for publication in A&A

Aims: We investigate magnetic tracers in the photosphere and the chromosphere of the ultra-rapid rotator ($P\sim0.32d$) V530 Per, a cool member of the open cluster $\alpha$ Persei, to characterize the short-term variability of the magnetic activity and large-scale magnetic field of this prototypical young, rapidly rotating solar-like star. Methods: With time-resolved spectropolarimetric observations spread over four close-by nights, we reconstructed the brightness distribution and large-scale magnetic field geometry of V530 Per through Zeeman-Doppler imaging. Simultaneously, we estimated the short-term variability of the surface through latitudinal differential rotation. Using the same data set, we also mapped the spatial distribution of prominences through tomography of H$\alpha$ emission. Results: As in our previous study, a large dark spot occupies the polar region of V530 Per with smaller, dark, and bright spots at lower latitudes. The large-scale magnetic field is dominated by a toroidal, mostly axisymmetric component. The maximal radial field strength is equal to $\sim1$ kG. The surface differential rotation is consistent with a smooth Sun-like shear d$\Omega = 0.053 \pm 0.004$ rad.d$^{-1}$, close to the solar shear level. The prominence pattern displays a stable component that is confined close to the corotation radius. We also observe rapidly evolving H$\alpha$ emitting structures, over timescales ranging from minutes to days. The fast H$\alpha$ evolution was not linked to any detected photospheric changes in the spot or magnetic coverage., Comment: Accepted by A&A, 17 pages, 13 figures

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 an analysis of spectropolarimetric observations of the low-mass weak-line T Tauri stars TWA 25 and TWA 7. The large-scale surface magnetic fields have been reconstructed for both stars using the technique of Zeeman Doppler imaging. Our surface maps reveal predominantly toroidal and non-axisymmetric fields for both stars. These maps reinforce the wide range of surface magnetic fields that have been recovered, particularly in pre-main sequence stars that have stopped accreting from the (now depleted) central regions of their discs. We reconstruct the large scale surface brightness distributions for both stars, and use these reconstructions to filter out the activity-induced radial velocity jitter, reducing the RMS of the radial velocity variations from 495 m/s to 32 m/s for TWA 25, and from 127 m/s to 36 m/s for TWA 7, ruling out the presence of close-in giant planets for both stars. The TWA 7 radial velocities provide an example of a case where the activity-induced radial velocity variations mimic a Keplerian signal that is uncorrelated with the spectral activity indices. This shows the usefulness of longitudinal magnetic field measurements in identifying activity-induced radial velocity variations., Comment: 14 pages, 19 figures, Accepted for publication in the Monthly Notices of the Royal Astronomical Society

The helicity of a magnetic field is a fundamental property that is conserved in ideal MHD. It can be explored in the stellar context by mapping large-scale magnetic fields across stellar surfaces using Zeeman-Doppler imaging. A recent study of 51 stars in the mass range 0.1-1.34 M$_\odot$ showed that the photospheric magnetic helicity density follows a single power law when plotted against the toroidal field energy, but splits into two branches when plotted against the poloidal field energy. These two branches divide stars above and below $\sim$ 0.5 M$_\odot$. We present here a novel method of visualising the helicity density in terms of the linkage of the toroidal and poloidal fields that are mapped across the stellar surface. This approach allows us to classify the field linkages that provide the helicity density for stars of different masses and rotation rates. We find that stars on the lower-mass branch tend to have toroidal fields that are non-axisymmetric and so link through regions of positive and negative poloidal field. A lower-mass star may have the same helicity density as a higher-mass star, despite having a stronger poloidal field. Lower-mass stars are therefore less efficient at generating large-scale helicity., Comment: 8 pages, 7 figures, 1 table

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

We report on observations of the active K2 dwarf $\epsilon$ Eridani based on contemporaneous SPIRou, NARVAL, and TESS data obtained over two months in late 2018, when the activity of the star was reported to be in a non-cyclic phase. We first recover the fundamental parameters of the target from both visible and nIR spectral fitting. The large-scale magnetic field is investigated from polarimetric data. From unpolarized spectra, we estimate the total magnetic flux through Zeeman broadening of magnetically sensitive nIR lines and the chromospheric emission using the CaII H & K lines. The TESS photometric monitoring is modeled with pseudo-periodic Gaussian Process Regression. Fundamental parameters of $\epsilon$ Eridani derived from visible and near-infrared wavelengths provide us with consistent results, also in agreement with published values. We report a progressive increase of macroturbulence towards larger nIR wavelengths. Zeeman broadening of individual lines highlights an unsigned surface magnetic field $B_{\rm mono} = 1.90 \pm 0.13$ kG, with a filling factor $f = 12.5 \pm 1.7$% (unsigned magnetic flux $Bf = 237 \pm 36$ G). The large-scale magnetic field geometry, chromospheric emission, and broadband photometry display clear signs of non-rotational evolution over the course of data collection. Characteristic decay times deduced from the light curve and longitudinal field measurements fall in the range 30-40 d, while the characteristic timescale of surface differential rotation, as derived through the evolution of the magnetic geometry, is equal to $57 \pm 5$ d. The large-scale magnetic field exhibits a combination of properties not observed previously for $\epsilon$ Eridani, with a surface field among the weakest previously reported, but also mostly axisymmetric, and dominated by a toroidal component., Comment: A&A, in press, abridged abstract

Context: Pollux is considered as an archetype of a giant star hosting a planet. We then discovered a weak magnetic field at its surface using spectropolarimetry. Aims and Methods: We followed up our investigations on Pollux first using ESPaDOnS at CFHT and then Narval at TBL to obtain Stokes I and Stokes V spectra to study their variations for a duration of 4.25 years, that is, for more than two periods of about 590 d of the RV variations. We used the least-squares deconvolution (LSD) profiles to measure the longitudinal magnetic field and to perform a Zeeman Doppler imaging (ZDI) investigation. Results: The longitudinal magnetic field of Pollux is found to vary with a sinusoidal behavior and a period similar to that of the RV variations. From the ZDI investigation a rotation period of Pollux is determined to be equal to 660+/-15 days and possibly different than the period of variations of the RV. As to the magnetic topology, the poloidal component is dominant and almost purely dipolar with an inclination of 10.5{\deg} of the dipole with respect to the rotation axis. The mean strength of the surface magnetic field is 0.44 G. Conclusions: As to the origin of the magnetic field of Pollux, we favor the hypothesis that it is maintained through contemporaneous dynamo action. Pollux appears as the representative of a class of slowly rotating and weakly magnetic G-K red giants. To explain the sinusoidal RV variations of Pollux, two scenarios are proposed. If the RV period is different from the rotation period, the observed periodic RV variations are due to the hosted planet and the contribution of Pollux magnetic activity is not significantly detected. In the peculiar case in which the two periods are equal, we cannot discard the possibility that the activity of Pollux could explain the total RV variations and that the planet hypothesis would appear unnecessary., Comment: 9 pages, 7 figures, accepted for publication in Astronomy & Astrophysics

The distribution of eccentricities of warm giant exoplanets is commonly explained through planet--planet interactions, although no physically sound argument favours the ubiquity of such interactions. No simple, generic explanation has been put forward to explain the high mean eccentricity of these planets. In this paper, we revisit a simple, plausible explanation to account for the eccentricities of warm Jupiters: migration inside a cavity in the protoplanetary disc. Such a scenario allows to excite the outer eccentric resonances, a working mechanism for higher mass planets, leading to a growth in the eccentricity while preventing other, closer resonances to damp eccentricity. We test this idea with diverse numerical simulations, which show that the eccentricity of a Jupiter-mass planet around a Sun-like star can increase up to 0.4, a value never reached before with solely planet--disc interactions. This high eccentricity is comparable to, if not larger than, the median eccentricity of warm Saturn- to Jupiter-mass exoplanets. We also discuss the effects such a mechanism would have on exoplanet observations. This scenario could have strong consequences on the discs lifetime and the physics of inner disc dispersal, which could be constrained by the eccentricity distribution of gas giants.

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

(Abridged) We use optical spectroscopy to investigate the disk, wind, and accretion during the 2008 ZCMa NW outburst. Over 1000 optical emission lines reveal accretion, a variable, multi-component wind, and double-peaked lines of disk origin. The variable, non-axisymmetric, accretion-powered wind has slow ($\sim $0 km s$^{-1}$), intermediate ($\sim -$100 km s$^{-1}$) and fast ($\geq -$400 km s$^{-1}$) components. The fast components are of stellar origin and disappear in quiescence, while the slow component is less variable and could be related to a disk wind. The changes in the optical depth of the lines between outburst and quiescence are consistent with increased accretion being responsible for the observed outburst. We derive an accretion rate of 10$^{-4}$ M$_\odot$/yr in outburst. The Fe I and weak Fe II lines arise from an irradiated, flared disk at $\sim$0.5-3 $\times$M$_*$/16M$_\odot$ au with asymmetric upper layers, revealing that the energy from the accretion burst is deposited at scales below 0.5 au. Some line profiles have redshifted asymmetries, but the system is unlikely sustained by magnetospheric accretion, especially in outburst. The accretion-related structures extend over several stellar radii and, like the wind, are likely non-axisymmetric. The stellar mass may be $\sim$6-8 M$_\odot$, lower than previously thought ($\sim$16 M$_\odot$). Emission line analysis is found to be a powerful tool to study the innermost regions and accretion in stars within a very large range of effective temperatures. The density ranges in the disk and accretion structures are higher than in late-type stars, but the overall behavior, including the innermost disk emission and variable wind, is very similar independently of the spectral type. Our work suggests a common outburst behavior for stars with spectral types ranging from M-type to intermediate-mass stars., Comment: Accepted by A&A, 21 pages plus online appendix. Version 2: Minor typos corrected. Version 3: Note on DIBs added, typos corrected