Context. Multiplicity studies greatly benefit from focusing on M dwarfs because they are often paired in a variety of configurations with both stellar and substellar objects, including exoplanets. Aims. We aim to address the observed multiplicity of M dwarfs by conducting a systematic analysis using the latest available astropho-tometric data. Methods. For every star in a sample of 2214 M dwarfs from the CARMENES catalogue, we investigated the existence of resolved and unresolved physical companions in the literature and in all-sky surveys, especially in Gaia DR3 data products. We covered a very wide range of separations, from known spectroscopic binaries in tight arrangements (similar to 0.01 au) to remarkably separated ultra-wide pairs (similar to 10(5) au). Results. We identified 835 M dwarfs in 720 multiple systems, predominantly binaries. Thus, we propose 327 new binary candidates based on Gaia data. If these candidates are finally confirmed, we expect the multiplicity fraction of M dwarfs to be 40.3(-2.0)(+2.1)%. When only considering the systems already identified, the multiplicity fraction is reduced to 27.8(-1.8)(+1.9)%. This result is in line with most of the values published in the literature. We also identified M-dwarf multiple systems with FGK, white dwarf, ultra-cool dwarf, and exoplanet companions, as well as those in young stellar kinematic groups. We studied their physical separations, orbital periods, binding energies, and mass ratios. Conclusions. We argue that based on reliable astrometric data and spectroscopic investigations from the literature (even when considering detection biases), the multiplicity fraction of M dwarfs could still be significantly underestimated. This calls for further high-resolution follow-up studies to validate these findings.

Context. K-dwarf stars are promising targets in the exploration of potentially habitable planets. Their properties, falling between G and M dwarfs, provide an optimal trade-off between the prospect of habitability and ease of detection. The KOBE experiment is a blind-search survey exploiting this niche, monitoring the radial velocity of 50 late-type K-dwarf stars. It employs the CARMENES spectrograph, with an observational strategy designed to detect planets in the habitable zone of their system. Aims. In this work, we exploit the KOBE data set to characterize planetary signals in the K7 V star HIP 5957 (KOBE-1) and to constrain the planetary population within its habitable zone. Methods. We used 82 CARMENES spectra over a time span of three years. We employed a generalized Lomb-Scargle periodogram to search for significant periodic signals that would be compatible with Keplerian motion on KOBE-1. We carried out a model comparison within a Bayesian framework to ensure the significance of the planetary model over alternative configurations of lower complexity. We also inspected two available TESS sectors in search of planetary signals. Results. We identified two signals: at P-b = 8.5 d and P-c = 29.7 d. We confirmed their planetary nature through ruling out other non-planetary configurations. Their minimum masses are 8.80 +/- 0.76 M-circle plus (KOBE-1 b), and 12.4 +/- 1.1 M-circle plus (KOBE-1 c), corresponding to absolute masses within the planetary regime at a high certainty (>99.7%). By analyzing the sensitivity of the CARMENES time series to additional signals, we discarded planets above 8.5 M-circle plus within the habitable zone. We identified a single transit-like feature in TESS, whose origin is still uncertain, but still compatible within 1 sigma with a transit from planet c. Conclusions. The KOBE-1 multi-planetary system, consisting of a relatively quiet K7-dwarf hosting two sub-Neptune-minimum- mass planets, establishes the first discovery from the KOBE experiment. We have explored future prospects for characterizing this system, concluding that Gaia DR4 will be insensitive to their astrometric signature. Meanwhile, nulling interferometry with the Large Interferometer For Exoplanets (LIFE) mission could be capable of directly imaging both planets and characterizing their atmospheres in future studies.

Context. GJ 581 is a nearby M dwarf known to host a packed multiple planet system composed of two super-Earths and a Neptune-mass planet. We present new orbital analyses of the GJ 581 system, utilizing recent radial velocity (RV) data obtained from the CARMENES spectrograph combined with newly reprocessed archival data from the HARPS and HIRES spectrographs. Aims. Our aim was to analyze the post-discovery spectroscopic data of GJ581, which were obtained with CARMENES. In addition, we used publicly available HIRES and HARPS spectroscopic data to seek evidence of the known and disputed exoplanets in this system. We aimed to investigate the stellar activity of GJ 581 and update the planetary system's orbital parameters using state-of-the-art numerical models and techniques. Methods. We performed a periodogram analysis of the available precise CARMENES, HIRES, and HARPS RVs and of stellar activity indicators. We conducted detailed orbital analyses by testing various orbital configurations consistent with the RV data. We studied the posterior probability distribution of the parameters fit to the data and we explored the long-term stability and overall orbital dynamics of the GJ 581 system. Results. We refined the orbital parameters of the GJ 581 system using the most precise and complete set of Doppler data available. Consistent with the existing literature, our analysis confirms that the system is unequivocally composed of only three planets detectable in the present data, dismissing the putative planet GJ 581 d as an artifact of stellar activity. Our N-body fit reveals that the system's inclination is i = 47.0(-13.0)(+14.6) deg, which implies that the planets could be up to 30% more massive than their previously reported minimum masses. Furthermore, we report that the GJ 581 system exhibits long-term stability, as indicated by the posterior probability distribution, characterized by secular dynamical interactions without the involvement of mean motion resonances.

The two known planets in the planetary system of Teegarden's Star are among the most Earth-like exoplanets currently known. Revisiting this nearby planetary system with two planets in the habitable zone aims at a more complete census of planets around very low-mass stars. A significant number of new radial velocity measurements from CARMENES, ESPRESSO, MAROON-X, and HPF, as well as photometry from TESS motivated a deeper search for additional planets. We confirm and refine the orbital parameters of the two know planets Teegarden's Star b and c. We also report the detection of a third planet d with an orbital period of 26.13+-0.04 d and a minimum mass of 0.82+-0.17 M_Earth. A signal at 96 d is attributed to the stellar rotation period. The interpretation of a signal at 172 d remains open. The TESS data exclude transiting short-period planets down to about half an Earth radius. We compare the planetary system architecture of very low-mass stars. In the currently known configuration, the planetary system of Teegarden's star is dynamically quite different from that of TRAPPIST-1, which is more compact, but dynamically similar to others such as GJ 1002.

Radial velocity (RV) jitter represents an intrinsic limitation on the precision of Doppler searches for exoplanets that can originate from both instrumental and astrophysical sources. We aim to determine the RV jitter floor in M dwarfs and investigate the stellar properties that lead to RV jitter induced by stellar activity. We determined the RV jitter in 239 M dwarfs from the CARMENES survey that are predominantly of mid to late spectral type and solar metallicity. We also investigated the correlation between stellar rotation and magnetic fields with RV jitter. The median jitter in the CARMENES sample is 3.1 m/s, and it is 2.3 m/s for stars with an upper limit of 2 km/s on their projected rotation velocities. We provide a relation between the stellar equatorial rotation velocity and RV jitter in M dwarfs based on a subsample of 129 well-characterized CARMENES stars. RV jitter induced by stellar rotation dominates for stars with equatorial rotation velocities greater than 1 km/s. A jitter floor of 2 m/s dominates in stars with equatorial rotation velocities below 1 km/s. This jitter floor likely contains contributions from stellar jitter, instrumental jitter, and undetected companions. We study the impact of the average magnetic field and the distributions of magnetic filling factors on the RV jitter. We find a series of stars with excess RV jitter and distinctive distributions of magnetic filling factors. These stars are characterized by a dominant magnetic field component between 2-4 kG. An RV jitter floor can be distinguished from RV jitter induced by activity and rotation based on the stellar equatorial rotation velocity. RV jitter induced by activity and rotation primarily depends on the equatorial rotation velocity. This RV jitter is also related to the distribution of magnetic filling factors, and this emphasizes the role of the magnetic field in the generation of RV jitter.

Aims. We intended to quantify the impact of stellar multiplicity on the presence and properties of exoplanets. Methods. We investigated all exoplanet host stars at less than 100 pc using the latest astrometric data from Gaia DR3 and advanced statistical methodologies. We complemented our search for common proper motion and parallax companions with data from the Washington Double Star catalogue and the literature. After excluding a number of systems based on radial velocity data, and membership in clusters and open associations, or with resolved ultracool companions, we kept 215 exoplanet host stars in 212 multiple-star systems. Results. We found 17 new companions in the systems of 15 known exoplanet host stars, and we measured precise angular and projected physical separations and position angles for 236 pairs of stars, compiled key parameters for 276 planets in multiple systems, and established a comparison sample comprising 687 single stars with exoplanets. With all of this, we statistically analysed a series of hypotheses regarding planets in multiple stellar systems. Although they are only statistically significant at a 2 sigma level, our analysis pointed to several interesting results on the comparison in the mean number of planets in multiple versus single stellar systems and the tendency of high-mass planets to be located in closer orbits in multiple systems. We confirm that planets in multiple systems tend to have orbits with larger eccentricities than those in single systems. In particular, we found a significant (>4 sigma) preference for planets to exhibit high orbital eccentricities at small ratios between star-star projected physical separations and star-planet semi-major axes.

Context. GJ 581 is a nearby M dwarf known to host a packed multiple planet system composed of two super-Earths and a Neptune-mass planet. We present new orbital analyses of the GJ 581 system, utilizing recent radial velocity (RV) data obtained from the CARMENES spectrograph combined with newly reprocessed archival data from the HARPS and HIRES spectrographs. Aims. Our aim was to analyze the post-discovery spectroscopic data of GJ581, which were obtained with CARMENES. In addition, we used publicly available HIRES and HARPS spectroscopic data to seek evidence of the known and disputed exoplanets in this system. We aimed to investigate the stellar activity of GJ 581 and update the planetary system’s orbital parameters using state-of-the-art numerical models and techniques. Methods. We performed a periodogram analysis of the available precise CARMENES, HIRES, and HARPS RVs and of stellar activity indicators. We conducted detailed orbital analyses by testing various orbital configurations consistent with the RV data. We studied the posterior probability distribution of the parameters fit to the data and we explored the long-term stability and overall orbital dynamics of the GJ 581 system. Results. We refined the orbital parameters of the GJ 581 system using the most precise and complete set of Doppler data available. Consistent with the existing literature, our analysis confirms that the system is unequivocally composed of only three planets detectable in the present data, dismissing the putative planet GJ 581 d as an artifact of stellar activity. Our N-body fit reveals that the system’s inclination is i = 47.0−13.0+14.6 deg, which implies that the planets could be up to 30% more massive than their previously reported minimum masses. Furthermore, we report that the GJ 581 system exhibits long-term stability, as indicated by the posterior probability distribution, characterized by secular dynamical interactions without the involvement of mean motion resonances.

Aims. Knowledge of rotation periods (P-rot) is important for understanding the magnetic activity and angular momentum evolution of late-type stars, as well as for evaluating radial velocity signals of potential exoplanets and identifying false positives. We measured photometric and spectroscopic P-rot for a large sample of nearby bright M dwarfs with spectral types from M0 to M9, as part of our continual effort to fully characterize the Guaranteed Time Observation programme stars of the CARMENES survey. Methods. We analyse light curves chiefly from the SuperWASP survey and TESS data. We supplemented these with our own follow-up photometric monitoring programme from ground-based facilities, as well as spectroscopic indicator time series derived directly from the CARMENES spectra. Results. From our own analysis, we determined P-rot for 129 stars. Combined with the literature, we tabulated P-rot for 261 stars, or 75% of our sample. We developed a framework to evaluate the plausibility of all periods available for this sample by comparing them with activity signatures and checking for consistency between multiple measurements. We find that 166 of these stars have independent evidence that confirmed their P-rot. There are inconsistencies in 27 periods, which we classify as debated. A further 68 periods are identified as provisional detections that could benefit from independent verification. We provide an empirical relation for the P-rot uncertainty as a function of the P-rot value, based on the dispersion of the measurements. We show that published formal errors seem to be often underestimated for periods longwards of similar to 10 d. We examined rotation-activity relations with emission in X-rays, H alpha, Ca II H&K, and surface magnetic field strengths for this sample of M dwarfs. We find overall agreement with previous works, as well as tentative differences in the partially versus fully convective subsamples. We show P-rot as a function of stellar mass, age, and galactic kinematics. With the notable exception of three transiting planet systems and TZ Ari, all known planet hosts in this sample have P-rot greater than or similar to 15 d. Conclusions. Inherent challenges in determining accurate and precise stellar P-rot means independent verification is important, especially for inactive M dwarfs. Evidence of potential mass dependence in activity-rotation relations would suggest physical changes in the magnetic dynamo that warrants further investigation using larger samples of M dwarfs on both sides of the fully convective boundary. Important limitations need to be overcome before the radial velocity technique can be routinely used to detect and study planets around young and active stars.

Aims. Our goals are to characterise the kinematic properties and to identify young and old stars among the M dwarfs of the CARMENES input catalogue. Methods. We compiled the spectral types, proper motions, distances, and radial velocities for 2187 M dwarfs. We used the public code SteParKin to derive their galactic space velocities and identify members in the different galactic populations. We also identified candidate members in young stellar kinematic groups, with ages ranging from 1 Ma to 800 Ma with SteParKin, LACEwING, and BANYAN Sigma. We removed known close binaries and performed an analysis of kinematic, rotation, and activity indicators (rotational periods and projected velocities, H alpha, X-rays, and UV emission) for 1546 M dwarfs. We defined five rotation-activity-colour relations satisfied by young (tau less than or similar to 800 Ma) stars. Results. We identified 191 young M dwarf candidates (similar to 12%), 113 of which are newly recognised in this work. In this young sample, there are 118 very active stars based on H alpha emission, fast rotation, and X-ray and UV emission excess. Of them, 27 have also strong magnetic fields, 9 of which are likely younger than 50 Ma. Additionally, there are 87 potentially young stars and 99 stars with a dubious youth classification, which may increase the fraction of young stars to an astounding 24%. Only one star out of the 2187 exhibits kinematics typical of the old Galactic halo. Conclusions. A combined analysis of kinematic and rotation-activity properties provides a robust method for identifying young M dwarfs from archival data. However, more observational efforts are needed to ascertain the true nature of numerous young star candidates in the field and, perhaps more importantly, to precisely quantify their age.

This work aims to investigate the behaviour of the lithium abundance in stars with and without detected planets. Our study is based on a sample of 1332 FGK main-sequence stars with measured lithium abundances, for 257 of which planets were detected. Our method reviews the sample statistics and is addressed specifically to the influence of tides and orbital decay, with special attention to planets on close orbits, whose stellar rotational velocity is higher than the orbital period of the planet. In this case, tidal effects are much more pronounced. The analysis also covers the orbital decay on a short timescale, with planets spiralling into their parent star. Furthermore, the sample allows us to study the relation between the presence of planets and the physical properties of their host stars, such as the chromospheric activity, metallicity, and lithium abundance. In the case of a strong tidal influence, we cannot infer from any of the studies described that the behaviour of Li differs between stars that host planets and those that do not. Our sample includes stars with super-solar metallicity ([Fe/H]>0.15 dex) and a low lithium abundance (A(Li) <1.0 dex). This enabled us to analyse scenarios of the origin and existence of these stars. Considering the possible explanation of the F dip, we show that it is not a plausible scenario. Our analysis is based on a kinematic study and concludes that the possible time that elapsed in the travel from their birth places in the central regions of the Galaxy to their current positions in the solar neighbourhood is not enough to explain the high lithium depletion. It is remarkable that those of our high-metallicity low-lithium stars with the greatest eccentricity (e>0.2) are closest to the Galactic centre. A dedicated study of a set of high-metallicity low-Li stars is needed to test the migration-depletion scenario.

This work aims to investigate the behaviour of the lithium abundance in stars with and without detected planets. Our study is based on a sample of 1332 FGK main-sequence stars with measured lithium abundances, for 257 of which planets were detected. Our method reviews the sample statistics and is addressed specifically to the influence of tides and orbital decay, with special attention to planets on close orbits, whose stellar rotational velocity is higher than the orbital period of the planet. In this case, tidal effects are much more pronounced. The analysis also covers the orbital decay on a short timescale, with planets spiralling into their parent star. Furthermore, the sample allows us to study the relation between the presence of planets and the physical properties of their host stars, such as the chromospheric activity, metallicity, and lithium abundance. In the case of a strong tidal influence, we cannot infer from any of the studies described that the behaviour of Li differs between stars that host planets and those that do not. Our sample includes stars with super-solar metallicity ([Fe/H] > 0.15 dex) and a low lithium abundance (A(Li) < 1.0 dex). This enabled us to analyse scenarios of the origin and existence of these stars. Considering the possible explanation of the F dip, we show that it is not a plausible scenario. Our analysis is based on a kinematic study and concludes that the possible time that elapsed in the travel from their birth places in the central regions of the Galaxy to their current positions in the solar neighbourhood is not enough to explain the high lithium depletion. It is remarkable that those of our high-metallicity low-lithium stars with the greatest eccentricity (e > 0.2) are closest to the Galactic centre. A dedicated study of a set of high-metallicity low-Li stars is needed to test the migration-depletion scenario.

Context. Deep learning (DL) techniques are a promising approach among the set of methods used in the ever-challenging determination of stellar parameters in M dwarfs. In this context, transfer learning could play an important role in mitigating uncertainties in the results due to the synthetic gap (i.e. difference in feature distributions between observed and synthetic data). Aims. We propose a feature-based deep transfer learning (DTL) approach based on autoencoders to determine stellar parameters from high-resolution spectra. Using this methodology, we provide new estimations for the effective temperature, surface gravity, metallicity, and projected rotational velocity for 286 M dwarfs observed by the CARMENES survey. Methods. Using autoencoder architectures, we projected synthetic PHOENIX-ACES spectra and observed CARMENES spectra onto a new feature space of lower dimensionality in which the differences between the two domains are reduced. We used this low-dimensional new feature space as input for a convolutional neural network to obtain the stellar parameter determinations. Results. We performed an extensive analysis of our estimated stellar parameters, ranging from 3050 to 4300 K, 4.7 to 5.1 dex, and -0.53 to 0.25 dex for T-eff, log g, and [Fe/H], respectively. Our results are broadly consistent with those of recent studies using CARMENES data, with a systematic deviation in our T-eff scale towards hotter values for estimations above 3750 K. Furthermore, our methodology mitigates the deviations in metallicity found in previous DL techniques due to the synthetic gap. Conclusions. We consolidated a DTL-based methodology to determine stellar parameters in M dwarfs from synthetic spectra, with no need for high-quality measurements involved in the knowledge transfer. These results suggest the great potential of DTL to mitigate the differences in feature distributions between the observations and the PHOENIX-ACES spectra.

We report the confirmation and mass determination of a mini-Neptune transiting the M3.5 V star TOI-4438 (G 182-34) every 7.44 days. A transit signal was detected with NASA's TESS space mission in the sectors 40, 52, and 53. In order to validate the planet TOI-4438 b and to determine the system properties, we combined TESS data with high-precision radial velocity measurements from the CARMENES spectrograph, spanning almost one year, and ground-based transit photometry. We found that TOI-4438 b has a radius of Rb = 2.52 +/- 0.13 R_Earth (5% precision), which together with a mass of Mb=5.4 +/- 1.1 M_Earth (20% precision), results in a bulk density of rho = 1.85+0.51-0.44 g cm-3 (28% precision), aligning the discovery with a volatile-rich planet. Our interior structure retrieval with a pure water envelope yields a minimum water mass fraction of 46% (1-sigma). TOI-4438 b is a volatile-rich mini-Neptune with likely H/He mixed with molecules, such as water, CO_2, and CH_4. The primary star has a J-band magnitude of 9.7, and the planet has a high transmission spectroscopy metric (TSM) of 136 +/- 13. Taking into account the relatively warm equilibrium temperature of T_eq = 435 +/- 15 K, and the low activity level of its host star, TOI-4438 b is one of the most promising mini-Neptunes around an M dwarf for transmission spectroscopy studies.

HD 235088 (TOI-1430) is a young star known to host a sub-Neptune-sized planet candidate. We validated the planetary nature of HD 235088 b with multiband photometry, refined its planetary parameters, and obtained a new age estimate of the host star, placing it at 600–800 Myr. Previous spectroscopic observations of a single transit detected an excess absorption of He I coincident in time with the planet candidate transit. Here, we confirm the presence of He I in the atmosphere of HD 235088 b with one transit observed with CARMENES. We also detected hints of variability in the strength of the helium signal, with an absorption of −0.91 ± 0.11%, which is slightly deeper (2σ) than the previous measurement. Furthermore, we simulated the He I signal with a spherically symmetric 1D hydrodynamic model, finding that the upper atmosphere of HD 235088 b escapes hydrodynamically with a significant mass loss rate of (1.5−5) × 1010 g s−1 in a relatively cold outflow, with T = 3125 ±375 K, in the photon-limited escape regime. HD 235088 b (Rp = 2.045 ± 0.075 R⊕) is the smallest planet found to date with a solid atmospheric detection – not just of He I but any other atom or molecule. This positions it a benchmark planet for further analyses of evolving young sub-Neptune atmospheres.

The nature of the few known solar-mass stars simultaneously containing debris disks and planets remains an open question. A number of works have shown that this property appears to be independent of planetary masses as well as of stellar age, but possible correlations with stellar kinematics and metallicity have not been investigated. In this paper, we show that the majority of known stars containing both debris disks and planets belong to the metal-enriched Galactic thin disk. The few exceptions are stars that seem to be born in the star formation peak occurring in times of thick disk formation (i.e., HD 10700, HD 20794, and HD 40307), that is, between 11 and 8 Gyr. The mass of the dusty disk of these three old stars measured at 70 μm is very small - in fact, it is lower than that of the Kuiper belt of our Solar system by several orders of magnitude. These results are not surprising, as they remain within the values expected for the stellar disk evolution of such primitive stars. In parallel, we found another six thick-disk stars containing only debris disks or planets. These results enable us to establish a correlation between stellar metallicity and the mass of the dust disk modulated by the different formation epochs of the thick and thin Galactic disks.

Context. Detecting a planetary companion in a short-period orbit through radio emission from the interaction with its host star is a new prospect in exoplanet science. Recently, a tantalising signal was found close to the low-mass stellar system GJ 1151 using LOFAR observations. Aims. We studied spectroscopic time-series data of GJ 1151 in order to search for planetary companions, investigate possible signatures of stellar magnetic activity, and to find possible explanations for the radio signal. Methods. We used the combined radial velocities measured from spectra acquired with the CARMENES, HARPS-N, and HPF instruments, extracted activity indices from those spectra in order to mitigate the impact of stellar magnetic activity on the data, and performed a detailed analysis of Gaia astrometry and all available photometric time series coming from the MEarth and ASAS-SN surveys. Results. We found a M > 10.6 M⊕ companion to GJ 1151 in a 390d orbit at a separation of 0.57 au. Evidence for a second modulation is also present; this could be due to long-term magnetic variability or a second (substellar) companion. The star shows episodes of elevated magnetic activity, one of which could be linked to the observed LOFAR radio emission. We show that it is highly unlikely that the detected GJ 1151 c, or any additional outer companion is the source of the detected signal. We cannot firmly rule out the suggested explanation of an undetected short-period planet that could be related to the radio emission, as we establish an upper limit of 1.2 M⊕ for the minimum mass.

The CARMENES instrument was conceived to deliver high-accuracy radial velocity (RV) measurements with long-term stability to search for temperate rocky planets around a sample of nearby cool stars. The broad wavelength coverage was designed to provide a range of stellar activity indicators to assess the nature of potential RV signals and to provide valuable spectral information to help characterise the stellar targets. The CARMENES Data Release 1 (DR1) makes public all observations obtained during the CARMENES guaranteed time observations, which ran from 2016 to 2020 and collected 19,633 spectra for a sample of 362 targets. The CARMENES survey target selection was aimed at minimising biases, and about 70% of all known M dwarfs within 10 pc and accessible from Calar Alto were included. The data were pipeline-processed, and high-level data products, including 18,642 precise RVs for 345 targets, were derived. Time series data of spectroscopic activity indicators were also obtained. We discuss the characteristics of the CARMENES data, the statistical properties of the stellar sample, and the spectroscopic measurements. We show examples of the use of CARMENES data and provide a contextual view of the exoplanet population revealed by the survey, including 33 new planets, 17 re-analysed planets, and 26 confirmed planets from transiting candidate follow-up. A subsample of 238 targets was used to derive updated planet occurrence rates, yielding an overall average of 1.44+/-0.20 planets with 1 M_Earth < M sin i < 1000 M_Earth and 1 d < P_orb < 1000 d per star, and indicating that nearly every M dwarf hosts at least one planet. CARMENES data have proven very useful for identifying and measuring planetary companions as well as for additional applications, such as the determination of stellar properties, the characterisation of stellar activity, and the study of exoplanet atmospheres.

We report the discovery of two exoplanets around the M dwarfs GJ 724 and GJ 3988 using the radial velocity (RV) method. We obtained a total of 153 3.5 m Calar Alto/CARMENES spectra for both targets and measured their RVs and activity indicators. We also added archival ESO/HARPS data for GJ 724 and infrared RV measurements from Subaru/IRD for GJ 3988. We searched for periodic and stable signals to subsequently construct Keplerian models, considering different numbers of planets, and we selected the best models based on their Bayesian evidence. Gaussian process (GP) regression was included in some models to account for activity signals. For both systems, the best model corresponds to one single planet. The minimum masses are $10.75^{+0.96}_{-0.87}$ and $3.69^{+0.42}_{-0.41}$ Earth-masses for GJ 724 b and GJ 3988 b, respectively. Both planets have short periods (P < 10 d) and, therefore, they orbit their star closely (a < 0.05 au). GJ 724 b has an eccentric orbit (e = $0.577^{+0.055}_{-0.052}$), whereas the orbit of GJ 3988 b is circular. The high eccentricity of GJ 724 b makes it the most eccentric single exoplanet (to this date) around an M dwarf. Thus, we suggest a further analysis to understand its configuration in the context of planetary formation and architecture. In contrast, GJ 3988 b is an example of a common type of planet around mid-M dwarfs.

We report the discovery, mass, and radius determination of TOI-1801 b, a temperate mini-Neptune around a young M dwarf. TOI-1801 b was observed in TESS sectors 22 and 49, and the alert that this was a TESS planet candidate with a period of 21.3 days went out in April 2020. However, ground-based follow-up observations, including seeing-limited photometry in and outside transit together with precise radial velocity (RV) measurements with CARMENES and HIRES revealed that the true period of the planet is 10.6 days. These observations also allowed us to retrieve a mass of 5.74 ± 1.46 M_⊕, which together with a radius of 2.08 ± 0.12 R_⊕, means that TOI-1801 b is most probably composed of water and rock, with an upper limit of 2% by mass of H_2 in its atmosphere. The stellar rotation period of 16 days is readily detectable in our RV time series and in the ground-based photometry. We derived a likely age of 600–800 Myr for the parent star TOI-1801, which means that TOI-1801 b is the least massive young mini-Neptune with precise mass and radius determinations. Our results suggest that if TOI-1801 b had a larger atmosphere in the past, it must have been removed by some evolutionary mechanism on timescales shorter than 1 Gyr.

We report the discovery of HN Lib b, a sub-Neptunian mass planet orbiting the nearby ( d ≈ = 6.25 pc) M4.0 V star HN Lib detected by our CARMENES radial-velocity (RV) survey. We determined a planetary minimum mass of M b sin i = 5.46 ± 0.75 M ⊕ and an orbital period of P b = 36.116 ± 0.029 d, using ~5 yr of CARMENES data, as well as archival RVs from HARPS and HIRES spanning more than 13 yr. The flux received by the planet equals half the instellation on Earth, which places it in the middle of the conservative habitable zone (HZ) of its host star. The RV data show evidence for another planet candidate with M [c] sin i = 9.7 ± 1.9 M ⊕ and P [c] = 113.46 ± 0.20 d. The long-term stability of the signal and the fact that the best model for our data is a two-planet model with an independent activity component stand as strong arguments for establishing a planetary origin. However, we cannot rule out stellar activity due to its proximity to the rotation period of HN Lib, which we measured using CARMENES activity indicators and photometric data from a ground-based multi-site campaign as well as archival data. The discovery adds HN Lib b to the shortlist of super-Earth planets in the habitable zone of M dwarfs, but HN Lib [c] probably cannot be inhabited because, if confirmed, it would most likely be an icy giant.