Search

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 perform an analysis of kinematic, rotation, and activity indicators (rotational periods and projected velocities, Halpha, X-rays, and UV emission) for 1546 M dwarfs. We defined five rotation-activity-colour relations satisfied by young ({\tau} <= 800 Ma) stars. Results. We identified 191 young M dwarf candidates (~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.

As part of our comprehensive, ongoing characterisation of the low-mass end of the main sequence in the Solar neighbourhood, we used the OSIRIS instrument at the 10.4 m Gran Telescopio Canarias to acquire low- and mid-resolution (R$\approx$300 and R$\approx$2500) optical spectroscopy of 53 late-M and L ultracool dwarfs. Most of these objects are known but poorly investigated and lacking complete kinematics. We measured spectral indices, determined spectral types (six of which are new) and inferred effective temperature and surface gravity from BT-Settl synthetic spectra fits for all objects. We were able to measure radial velocities via line centre fitting and cross correlation for 46 objects, 29 of which lacked previous radial velocity measurements. Using these radial velocities in combination with the latest Gaia DR3 data, we also calculated Galactocentric space velocities. From their kinematics, we identified two candidates outside of the thin disc and four in young stellar kinematic groups. Two further ultracool dwarfs are apparently young field objects: 2MASSW J1246467$+$402715 (L4$\beta$), which has a potential, weak lithium absorption line, and G 196$-$3B (L3$\beta$), which was already known as young due to its well-studied primary companion., Comment: 33 pages, 14 figures, Accepted by MNRAS

We present the abundances of magnesium (Mg) and silicon (Si) for 314 dwarf stars with spectral types in the interval K7.0-M5.5 (Teff range ~4200-3050 K) observed with the CARMENES high-resolution spectrograph at the 3.5 m telescope at the Calar Alto Observatory. Our analysis employs the BT-Settl model atmospheres, the radiative transfer code Turbospectrum, and a state-of-the-art selection of atomic and molecular data. These Mg and Si abundances are critical for understanding both the chemical evolution and assembly of the Milky Way and the formation and composition of rocky planets. Our chemical abundances show a line-to-line scatter at the level of 0.1 dex for all studied spectral types. The typical error bar of our chemical abundance measurements is +- 0.11 dex (Mg) and +- 0.16 dex (Si) for all spectral types based on the comparison of the results obtained for stellar components of multiple systems. The derived abundances are consistent with the galactic evolution trends and observed chemical abundance distribution of earlier FGK-type stars in the solar neighbourhood. Besides, our analysis provides compatible abundances for stars in multiple systems. In addition, we studied the abundances of different galactic stellar populations. In this paper, we also explore the relation of the Mg and Si abundances of stars with and without known planets., Comment: 11 Pages, 5 figures. Accepted for publication in A&A

Understanding planet formation is important in the context of the origin of planetary systems in general and of the Solar System in particular, as well as to predict the likelihood of finding Jupiter, Neptune, and Earth analogues around other stars. We aim to precisely determine the radii and dynamical masses of transiting planets orbiting the young M star AU Mic using public photometric and spectroscopic datasets. We characterise the stellar activity and physical properties (radius, mass, density) of the transiting planets in the young AU Mic system through joint transit and radial velocity fits with Gaussian processes. We determine a radius of $R^{b}$= 4.79 +/- 0.29 R$_\oplus$, a mass of $M^{b}$= 9.0 +/- 2.7 M$_\oplus$, and a bulk density of $\rho^{b}$ = 0.49 +/- 0.16 g cm$^{-3}$ for the innermost transiting planet AU Mic b. For the second known transiting planet, AU Mic c, we infer a radius of $R^{c}$= 2.79 +/- 0.18 R$_\oplus$, a mass of $M^{c}$= 14.5 +/- 3.4 M$_\oplus$, and a bulk density of $\rho^{c}$ = 3.90 +/- 1.17 g cm$^{-3}$. According to theoretical models, AU Mic b may harbour an H2 envelope larger than 5\% by mass, with a fraction of rock and a fraction of water. AU Mic c could be made of rock and/or water and may have an H2 atmosphere comprising at most 5\% of its mass. AU Mic b has retained most of its atmosphere but might lose it over tens of millions of years due to the strong stellar radiation, while AU Mic c likely suffers much less photo-evaporation because it lies at a larger separation from its host. Using all the datasets in hand, we determine a 3$\sigma$ upper mass limit of $M^{[d]}\sin{i}$ = 8.6 M$_{\oplus}$ for the AU Mic 'd' TTV-candidate. In addition, we do not confirm the recently proposed existence of the planet candidate AU Mic 'e' with an orbital period of 33.4 days., Comment: Accepted in A&A. 28 pages, 15 figures

GJ 581 is a nearby M dwarf known to host a packed multiple planet system with 2 super-Earths and a Neptune-mass planet. We present new orbital analyses of the system, utilizing recent RV data obtained from the CARMENES spectrograph combined with newly reprocessed archival data from the HARPS and HIRES spectrographs. Our aim was to analyze the post-discovery spectroscopic data of GJ 581, 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. 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 explored the long-term stability and overall orbital dynamics of the system. We refined the orbital parameters of the system using the most precise and complete set of Doppler data available. Consistent with the existing literature, we confirm that the system is unequivocally composed of only 3 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 deg, which implies that the planets could be up to 30% more massive than their previously reported minimum masses. Furthermore, we report that the system exhibits long-term stability, as indicated by the posterior probability distribution, characterized by secular dynamical interactions without the involvement of mean motion resonances., Comment: 21 pages, 8 Figures, 5 Tables, 1 Appendix, Accepted for publication in A&A (20th March 2024), Data available at the CDS

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). 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. 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. 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 Teff, logg, and [Fe/H], respectively. Our results are broadly consistent with those of recent studies using CARMENES data, with a systematic deviation in our Teff 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. 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., Comment: Accepted in A&A

Ultrahot Jupiters are a type of gaseous exoplanet that orbit extremely close to their host star, resulting in significantly high equilibrium temperatures. In recent years, high-resolution emission spectroscopy has been broadly employed in observing the atmospheres of ultrahot Jupiters. We used the CARMENES spectrograph to observe the high-resolution spectra of the dayside hemisphere of MASCARA-1b in both visible and near-infrared. Through cross-correlation analysis, we detected signals of \ion{Fe}{i} and \ion{Ti}{i}. Based on these detections, we conducted an atmospheric retrieval and discovered the presence of a strong inversion layer in the planet's atmosphere. The retrieved Ti and Fe abundances are broadly consistent with solar abundances. In particular, we obtained a relative abundance of [Ti/Fe] as $-1.0 \pm 0.8$ under the free retrieval and $-0.4^{+0.5}_{-0.8}$ under the chemical equilibrium retrieval, suggesting the absence of significant titanium depletion on this planet. Furthermore, we considered the influence of planetary rotation on spectral line profiles. The resulting equatorial rotation speed was determined to be $4.4^{+1.6}_{-2.0}\,\mathrm{km\,s^{-1}}$, which agrees with the rotation speed induced by tidal locking., Comment: 18 pages, 11 figures; accepted for publication in A&A

During the first Gyr of their life, exoplanet atmospheres suffer from different atmospheric escape phenomena that can strongly affect the shape and morphology of the exoplanet itself. These processes can be studied with Ly$\alpha$, H$\alpha$ and/or He I triplet observations. We present high-resolution spectroscopy observations from CARMENES and GIARPS checking for He I and H$\alpha$ signals in 20 exoplanetary atmospheres: V1298Tau c, K2-100b, HD63433b, HD63433c, HD73583b, HD73583c, K2-77b, TOI-2076b, TOI-2048b, HD235088b, TOI-1807b, TOI-1136d, TOI-1268b, TOI-1683b, TOI-2018b, MASCARA-2b, WASP-189b, TOI-2046b, TOI-1431b, and HAT-P-57b. We report two new high-resolution spectroscopy He I detections for TOI-1268b and TOI-2018b, and an H$\alpha$ detection for TOI-1136d. The MOPYS (Measuring Out-flows in Planets orbiting Young Stars) project aims to understand the evaporating phenomena and test their predictions from the current observations. We compiled a list of 70 exoplanets with He I and/or H$\alpha$ observations, from this work and the literature, and we considered the He I and H$\alpha$ results as proxy for atmospheric escape. Our principal results are that 0.1-1Gyr-old planets do not exhibit more He I or H$\alpha$ detections than older planets, and evaporation signals are more frequent for planets orbiting $\sim$1-3Gyr-old stars. We provide new constrains to the cosmic shoreline, the empirical division between rocky planets and planets with atmosphere, by using the evaporation detections and explore the capabilities of a new dimensionless parameter, $R_{\rm He}/R_{\rm Hill}$, to explain the He I triplet detections. Furthermore, we present a statistically significant upper boundary for the He I triplet detections in the $T_{\rm eq}$ vs $\rho_{\rm p}$ parameter space. Planets located above that boundary are unlikely to show He I absorption signals., Comment: Accepted in A&A. 64 pages, many figures. Supplementary material in Zenodo

The Li abundance observed in pre-main sequence and main sequence late-type stars is strongly age-dependent, but also shows a complex pattern depending on several parameters, such as rotation, chromospheric activity and metallicity. The best way to calibrate these effects, with the aim of studying Li as an age indicator for FGK stars, is to calibrate coeval groups of stars, such as open clusters (OCs) and associations. We present a considerable target sample of 42 OCs and associations, ranging from 1 Myr to 5 Gyr, observed within the Gaia-ESO survey (GES), and using the latest data provided by GES iDR6 and the most recent release of Gaia that was then available, EDR3. As part of this study, we update and improve the membership analysis for all 20 OCs presented in our previous article. We perform detailed membership analyses for all target clusters to identify likely candidates, using all available parameters provided by GES and based on numerous criteria: from radial velocity distributions, to the astrometry and photometry provided by Gaia, to gravity indicators, [Fe/H] metallicity, and Li content. We obtain updated lists of cluster members for the whole target sample, as well as a selection of Li-rich giant contaminants obtained as an additional result of the membership process. Each selection of cluster candidates was thoroughly contrasted with numerous existing membership studies using data from Gaia to ensure the most robust results. These final cluster selections will be used in the third and last paper of this series, which reports the results of a comparative study characterising the observable Li dispersion in each cluster and analysing its dependence on several parameters, allowing us to calibrate a Li-age relation and obtain a series of empirical Li envelopes for key ages in our sample., Comment: 83 pages, 246 figures. To be published in A&A, accepted January 27, 2024

The Gaia-ESO Survey is an European Southern Observatory (ESO) public spectroscopic survey that targeted $10^5$ stars in the Milky Way covering the major populations of the disk, bulge and halo. The observations were made using FLAMES on the VLT obtaining both UVES high ($R\sim47,000$) and GIRAFFE medium ($R\sim20,000$) resolution spectra. The analysis of the Gaia-ESO spectra was the work of multiple analysis teams (nodes) within five working groups (WG). The homogenisation of the stellar parameters within WG11 (high resolution observations of FGK stars) and the homogenisation of the stellar parameters within WG10 (medium resolution observations of FGK stars) is described here. In both cases, the homogenisation was carried out using a bayesian Inference method developed specifically for the Gaia-ESO Survey by WG11. The WG10 homogenisation primarily used the cross-match of stars with WG11 as the reference set in both the stellar parameter and chemical abundance homogenisation. In this way the WG10 homogenised results have been placed directly onto the WG11 stellar parameter and chemical abundance scales. The reference set for the metal-poor end was sparse which limited the effectiveness of the homogenisation in that regime. For WG11, the total number of stars for which stellar parameters were derived was 6,231 with typical uncertainties for Teff, log g and [Fe/H] of 32~K, 0.05 and 0.05 respectively. One or more chemical abundances out of a possible 39 elements were derived for 6,188 of the stars. For WG10, the total number of stars for which stellar parameters were derived was 76,675 with typical uncertainties for Teff, log g and [Fe/H] of 64~K, 0.15 and 0.07 respectively. One or more chemical abundances out of a possible 30 elements were derived for 64,177 of the stars., Comment: 31 pages, 19 figures

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., Comment: Astronomy & Astrophysics, accepted; 21 pages, 18 figures

Planets with orbital periods shorter than 1 day are rare and have formation histories that are not completely understood. Small ($R_\mathrm{p} < 2\; R_\oplus$) ultra-short-period (USP) planets are highly irradiated, probably have rocky compositions with high bulk densities, and are often found in multi-planet systems. Additionally, USP planets found around small stars are excellent candidates for characterization using present-day instrumentation. Of the current full sample of approximately 5500 confirmed exoplanets, only 130 are USP planets and around 40 have mass and radius measurements. Wolf 327 (TOI-5747) is an M dwarf ($R_\star = 0.406 \pm 0.015 \; R_\odot$, $M_\star = 0.405 \pm 0.019 \; M_\odot$, $T_{\mathrm{eff}}=3542 \pm 70$ K, and $V = 13$ mag) located at a distance $d = 28.5$ pc. NASA's planet hunter satellite, TESS, detected transits in this star with a period of 0.573 d (13.7 h) and with a transit depth of 818 ppm. Ground-based follow-up photometry, high resolution imaging, and radial velocity (RV) measurements taken with the CARMENES spectrograph confirm the presence of this new USP planet. Wolf 327b is a super-Earth with a radius of $R_\mathrm{p} = 1.24 \pm 0.06 \; R_\oplus$ and a mass of $M_\mathrm{p} = 2.53 \pm 0.46 \; M_\oplus$, yielding a bulk density of $7.24 \pm 1.66 $\,g cm$^{-3}$ and thus suggesting a rocky composition. Owing to its close proximity to its host star ($a = 0.01$ au), Wolf 327b has an equilibrium temperature of $996 \pm 22$ K. This planet has a mass and radius similar to K2-229b, a planet with an inferred Mercury-like internal composition. Planet interior models suggest that Wolf 327b has a large iron core, a small rocky mantle, and a negligible (if any) H/He atmosphere., Comment: 26 pages, 19 figures. Accepted for publication in A&A

To date, the hot Jupiter WASP-12 b has been the only planet with confirmed orbital decay. The late F-type host star has been hypothesized to be surrounded by a large structure of circumstellar material evaporated from the planet. We obtained two high-resolution spectral transit time series with CARMENES and extensively searched for absorption signals by the atomic species Na, H, Ca, and He using transmission spectroscopy, thereby covering the He I triplet with high resolution for the first time. We apply SYSREM for atomic line transmission spectroscopy, introduce the technique of signal protection to improve the results for individual absorption lines, and compare the outcomes to those of established methods. No transmission signals were detected and the most stringent upper limits as of yet were derived for the individual indicators. Nonetheless, we found variation in the stellar Halpha and He I lines, the origin of which remains uncertain but is unlikely to be activity. To constrain the enigmatic activity state of WASP-12, we analyzed XMM-Newton X-ray data and found the star to be moderately active at most. We deduced an upper limit for the X-ray luminosity and the irradiating X-ray and extreme ultraviolet (XUV) flux of WASP-12 b. Based on the XUV flux upper limit and the lack of the He I signal, our hydrodynamic models slightly favor a moderately irradiated planet with a thermospheric temperature of <= 12,000 K, and a conservative upper limit of <= 4e12 g/s on the mass-loss rate. Our study does not provide evidence for an extended planetary atmosphere or absorption by circumstellar material close to the planetary orbit., Comment: Accepted for publication in A&A

Light from celestial objects interacts with the molecules of the Earth's atmosphere, resulting in the production of telluric absorption lines in ground-based spectral data. Correcting for these lines, which strongly affect red and infrared wavelengths, is often needed in a wide variety of scientific applications. Here, we present the template division telluric modeling (TDTM) technique, a method for accurately removing telluric absorption lines in stars that exhibit numerous intrinsic features. Based on the Earth's barycentric motion throughout the year, our approach is suited for disentangling telluric and stellar spectral components. By fitting a synthetic transmission model, telluric-free spectra are derived. We demonstrate the performance of the TDTM technique in correcting telluric contamination using a high-resolution optical spectral time series of the feature-rich M3.0 dwarf star Wolf 294 that was obtained with the CARMENES spectrograph. We apply the TDTM approach to the CARMENES survey sample, which consists of 382 targets encompassing 22357 optical and 20314 near-infrared spectra, to correct for telluric absorption. The corrected spectra are coadded to construct template spectra for each of our targets. This library of telluric-free, high signal-to-noise ratio, high-resolution (R>80000) templates comprises the most comprehensive collection of spectral M-dwarf data available to date, both in terms of quantity and quality, and is available at the project website (http://carmenes.cab.inta-csic.es)., Comment: 31 pages, 24 figures, 3 tables, accepted for publication by A&A

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 $\pm$ 1.46 $M_\oplus$, which together with a radius of 2.08 $\pm$ 0.12 $R_\oplus$, 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., Comment: Accepted in A&A. 29 pages, 21 figures

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., Comment: A&A in press

The hydrogen Paschen lines are known activity indicators, but studies of them in M~dwarfs during quiescence are as rare as their reports in flare studies. This situation is mostly caused by a lack of observations, owing to their location in the near-infrared regime, which is covered by few high-resolution spectrographs. We study the Pa$\beta$ line, using a sample of 360 M~dwarfs observed by the CARMENES spectrograph. Descending the spectral sequence of inactive M~stars in quiescence, we find the Pa$\beta$ line to get shallower until about spectral type M3.5 V, after which a slight re-deepening is observed. Looking at the whole sample, for stars with H$\alpha$ in absorption, we find a loose anti-correlation between the (median) pseudo-equivalent widths (pEWs) of H$\alpha$ and Pa$\beta$ for stars of similar effective temperature. Looking instead at time series of individual stars, we often find correlation between pEW(H$\alpha$) and pEW(Pa$\beta$) for stars with H$\alpha$ in emission and an anti-correlation for stars with H$\alpha$ in absorption. Regarding flaring activity, we report the automatic detection of 35 Paschen line flares in 20 stars. Additionally we found visually six faint Paschen line flares in these stars plus 16 faint Paschen line flares in another 12 stars. In strong flares, Paschen lines can be observed up to Pa 14. Moreover, we find that Paschen line emission is almost always coupled to symmetric H$\alpha$ line broadening, which we ascribe to Stark broadening, indicating high pressure in the chromosphere. Finally we report a few Pa$\beta$ line asymmetries for flares that also exhibit strong H$\alpha$ line asymmetries., Comment: 21 pages, 22 figures, 7 tables, accepted to A&A

HD235088 (TOI-1430) is a young star known to host a sub-Neptune-sized planet candidate. We validated the planetary nature of HD235088 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 HD235088 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$\pm$0.11%, which is slightly deeper (2$\sigma$) than the previous measurement. Furthermore, we simulated the He I signal with a spherically symmetric 1D hydrodynamic model, finding that the upper atmosphere of HD235088 b escapes hydrodynamically with a significant mass loss rate of (1.5-5) $\times$10$^{10}$g s$^{-1}$, in a relatively cold outflow, with $T$=3125$\pm$375 K, in the photon-limited escape regime. HD235088 b ($R_{p}$ = 2.045$\pm$0.075 R$_{\oplus}$) 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., Comment: Accepted for publication in A&A. 17 pages, 18 figures

Aims. A transiting planet candidate with a sub-Neptune radius orbiting the nearby ($d$ = 51.9$\pm$0.07 pc) M1.5 V star TOI-1470 with a period of $\sim$2.5 d was announced by the NASA Transiting Exoplanet Survey Satellite (TESS), which observed the field of TOI-1470 in four different sectors. We aim to validate its planetary nature using precise radial velocities (RVs) taken with the CARMENES spectrograph. Methods. We obtained 44 RV measurements with CARMENES spanning eight months between 3 June 2020 and 17 January 2021. For a better characterization of the parent star activity, we also collected contemporaneous optical photometric observations at the Joan Or\'o and Sierra Nevada Observatories, and we retrieved archival photometry from the literature. We used ground-based photometric observations from MuSCAT and also from MuSCAT2 and MuSCAT3 to confirm the planetary transit signals. We performed a combined photometric and spectroscopic analysis by including Gaussian processes and Keplerian orbits to simultaneously account for the stellar activity and planetary signals. Results. We estimate that TOI-1470 has a rotation period of 29$\pm$3 d based on photometric and spectroscopic data. The combined analysis confirms the discovery of the announced transiting planet, TOI-1470 b, with an orbital period of 2.527093$\pm$0.000003 d, a mass of $7.32^{+1.21}_{-1.24}$ M$_{\oplus}$, and a radius of $2.18^{+0.04}_{-0.04}$ R$_{\oplus}$. We also discover a second transiting planet that was not announced previously by TESS, TOI-1470 c, with an orbital period of 18.08816$\pm$0.00006 d, a mass of $7.24^{+2.87}_{-2.77}$ M$_{\oplus}$, and a radius of $2.47^{+0.02}_{-0.02}$ R$_{\oplus}$. The two planets are placed on the same side of the radius valley of M dwarfs and lie between TOI-1470 and the inner border of its habitable zone., Comment: arXiv admin note: substantial text overlap with arXiv:2111.14602

We report the discovery of HN Lib b, a sub-Neptunian mass planet orbiting the nearby ($d \approx$ = 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_\text{b}\sin i = $ 5.46 $\pm$ 0.75 $\text{M}_\oplus$ and an orbital period of $P_\text{b} = $ 36.116 $\pm$ 0.029 d, using $\sim$5 yr of CARMENES data, as well as archival RVs from HARPS and HIRES spanning more than 13 years. 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_\text{[c]}\sin i = $ 9.7 $\pm$ 1.9 $\text{M}_\oplus$ and $P_\text{[c]} = $ 113.46 $\pm$ 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.

The Gaia-ESO Survey is a public spectroscopic survey that has targeted $\gtrsim10^5$ stars covering all major components of the Milky Way from the end of 2011 to 2018, delivering its public final release in May 2022. Unlike other spectroscopic surveys, Gaia-ESO is the only survey that observed stars across all spectral types with dedicated, specialised analyses: from O ($T_\mathrm{eff} \sim 30,000-52,000$~K) all the way to K-M ($\gtrsim$3,500~K). The physics throughout these stellar regimes varies significantly, which has previously prohibited any detailed comparisons between stars of significantly different type. In the final data release (internal data release 6) of the Gaia-ESO Survey, we provide the final database containing a large number of products such as radial velocities, stellar parameters and elemental abundances, rotational velocity, and also, e.g., activity and accretion indicators in young stars and membership probability in star clusters for more than 114,000 stars. The spectral analysis is coordinated by a number of Working Groups (WGs) within the Survey, which specialise in the various stellar samples. Common targets are analysed across WGs to allow for comparisons (and calibrations) amongst instrumental setups and spectral types. Here we describe the procedures employed to ensure all Survey results are placed on a common scale to arrive at a single set of recommended results for all Survey collaborators to use. We also present some general quality and consistency checks performed over all Survey results., Comment: A&A accepted, minor revision, 36 pages, 38 figures

The large amounts of astrophysical data being provided by existing and future instrumentation require efficient and fast analysis tools. Transfer learning is a new technique promising higher accuracy in the derived data products, with information from one domain being transferred to improve the accuracy of a neural network model in another domain. In this work, we demonstrate the feasibility of applying the deep transfer learning (DTL) approach to high-resolution spectra in the framework of photospheric stellar parameter determination. To this end, we used 14 stars of the CARMENES survey sample with interferometric angular diameters to calculate the effective temperature, as well as six M dwarfs that are common proper motion companions to FGK-type primaries with known metallicity. After training a deep learning (DL) neural network model on synthetic PHOENIX-ACES spectra, we used the internal feature representations together with those 14+6 stars with independent parameter measurements as a new input for the transfer process. We compare the derived stellar parameters of a small sample of M dwarfs kept out of the training phase with results from other methods in the literature. Assuming that temperatures from bolometric luminosities and interferometric radii and metallicities from FGK+M binaries are sufficiently accurate, DTL provides a higher accuracy than our previous state-of-the-art DL method (mean absolute differences improve by 20 K for temperature and 0.2 dex for metallicity from DL to DTL when compared with reference values from interferometry and FGK+M binaries). Furthermore, the machine learning (internal) precision of DTL also improves as uncertainties are five times smaller on average. These results indicate that DTL is a robust tool for obtaining M-dwarf stellar parameters comparable to those obtained from independent estimations for well-known stars.

Although the number of exoplanets reported in the literature exceeds 5000 so far, only a few dozen of them are young planets ($\le$900 Myr). However, a complete characterization of these young planets is key to understanding the current properties of the entire population. Hence, it is necessary to constrain the planetary formation processes and the timescales of dynamical evolution by measuring the masses of exoplanets transiting young stars. We characterize and measure the masses of two transiting planets orbiting the 400 Myr old solar-type star HD\,63433, which is a member of the Ursa Major moving group. We analysed precise photometric light curves of five sectors of the TESS mission with a baseline of $\sim$750 days and obtained $\sim$150 precise radial velocity measurements with the visible and infrared arms of the CARMENES instrument at the Calar Alto 3.5 m telescope in two different campaigns of $\sim$500 days. We performed a combined photometric and spectroscopic analysis to retrieve the planetary properties of two young planets. The strong stellar activity signal was modelled by Gaussian regression processes. We have updated the transit parameters of HD\,63433\,b and c and obtained planet radii of R$_p^b$\,=\,2.140\,$\pm$\,0.087 R$_\oplus$ and R$_p^c$\,=\,2.692\,$\pm$\,0.108 R$_\oplus$. Our analysis allowed us to determine the dynamical mass of the outer planet with a 4$\sigma$ significance ($M_p^c$\,=\,15.54\,$\pm$\,3.86 M$_\oplus$) and set an upper limit on the mass of the inner planet at 3$\sigma$ ($M_p^b$\,$<$\,21.76 M$_\oplus$). According to theoretical models, both planets are expected to be sub-Neptunes, whose interiors mostly consist of silicates and water with no dominant composition of iron, and whose gas envelopes are lower than 2\% in the case of HD\,63433\,c. The envelope is unconstrained in HD\,63433\,b.

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., Comment: Published in A&A (https://www.aanda.org/10.1051/0004-6361/202244879), 25 pages, 12 figures, Tables 1 and 2 only available online

Radial velocities (RVs) measured from high-resolution stellar spectra are routinely used to detect and characterise orbiting exoplanet companions. The different lines present in stellar spectra are created by several species, which are non-uniformly affected by stellar variability features such as spots or faculae. Stellar variability distorts the shape of the spectral absorption lines from which precise RVs are measured, posing one of the main problems in the study of exoplanets. In this work we aim to study how the spectral lines present in M dwarfs are independently impacted by stellar activity. We used CARMENES optical spectra of six active early- and mid-type M dwarfs to compute line-by-line RVs and study their correlation with several well-studied proxies of stellar activity. We are able to classify spectral lines based on their sensitivity to activity in five M dwarfs displaying high levels of stellar activity. We further used this line classification to compute RVs with activity-sensitive lines and less sensitive lines, enhancing or mitigating stellar activity effects in the RV time series. For specific sets of the least activity-sensitive lines, the RV scatter decreases by ~ 2 to 5 times the initial one, depending on the star. Finally, we compare these lines in the different stars analysed, finding the sensitivity to activity to vary from star to star. Despite the high density of lines and blends present in M dwarf stellar spectra, we find that a line-by-line approach is able to deliver precise RVs. Line-by-line RVs are also sensitive to stellar activity effects, and they allow for an accurate selection of activity-insensitive lines to mitigate activity effects in RV. However, we find stellar activity effects to vary in the same insensitive lines from star to star., Comment: Accepted for publication in A&A

Knowledge of stellar atmospheric parameters ($T_{\rm eff}$, $\log{g}$, [Fe/H]) of M dwarfs can be used to constrain both theoretical stellar models and Galactic chemical evolutionary models, and guide exoplanet searches, but their determination is difficult due to the complexity of the spectra of their cool atmospheres. In our ongoing effort to characterize M dwarfs, and in particular their chemical composition, we carried out multiband photometric calibrations of metallicity for early- and intermediate-type M dwarfs. The third Gaia data release provides high-precision astrometry and three-band photometry. This information, combined with the 2MASS and CatWISE2020 infrared photometric surveys and a sample of 4919 M dwarfs with metallicity values determined with high-resolution spectroscopy by The Cannon and APOGEE spectra, allowed us to study the effect of the metallicity in color-color and color-magnitude diagrams. We divided this sample into two subsamples: we used 1000 stars to train the calibrations with Bayesian statistics and Markov Chain Monte Carlo techniques, and the remaining 3919 stars to check the accuracy of the estimations. We derived several photometric calibrations of metallicity applicable to M dwarfs in the range of $-0.45\leq\text{[Fe/H]}\leq +0.45$ dex and spectral types down to M5.0 V that yield uncertainties down to the $0.10$ dex level. Lastly, we compared our results with other photometric estimations published in the literature for an additional sample of 46 M dwarfs in wide binary systems with FGK-type primary stars, and found a great predictive performance., Comment: Accepted for publication in ApJ, 24 pages, 14 figures and 7 tables

One of the main scientific goals of the TESS mission is the discovery of transiting small planets around the closest and brightest stars in the sky. Here, using data from the CARMENES, MAROON-X, and HIRES spectrographs, together with TESS, we report the discovery and mass determination of a planetary system around the M1.5 V star GJ 806 (TOI-4481). GJ 806 is a bright (V=10.8 mag, J=7.3 mag) and nearby (d=12 pc) M dwarf that hosts at least two planets. The innermost planet, GJ 806 b, is transiting and has an ultra-short orbital period of 0.93 d, a radius of 1.331+-0.023 Re, a mass of 1.90+-0.17 Me, a mean density of 4.40+-0.45 g/cm3, and an equilibrium temperature of 940+-10 K. We detect a second, non-transiting, super-Earth planet in the system, GJ 806c, with an orbital period of 6.6 d, a minimum mass of 5.80+-0.30 Me, and an equilibrium temperature of 490+-5 K. The radial velocity data also shows evidence for a third periodicity at 13.6 d, although the current dataset does not provide sufficient evidence to unambiguously distinguish between a third super-Earth mass (Msin(i)=8.50+-0.45 Me) planet or stellar activity. Additionally, we report one transit observation of GJ 806 b taken with CARMENES in search for a possible extended atmosphere of H or He, but we can only place upper limits to its existence. This is not surprising as our evolutionary models support the idea that any possible primordial H/He atmosphere that GJ 806 b might have had, would long have been lost. However, GJ 806b's bulk density makes it likely that the planet hosts some type of volatile atmosphere. In fact, with a transmission spectroscopy metrics (TSM) of 44 and an emission spectroscopy metrics (ESM) of 24, GJ 806 b the third-ranked terrestrial planet around an M dwarf suitable for transmission spectroscopy studies, and the most promising terrestrial planet for emission spectroscopy studies., Comment: Under second review in A&A. This paper is NOT yet accepted, but it is made openly available to the community due to the approaching JWST deadline

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. 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. 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. We found a M$>$10.6 M$_{\oplus}$ 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 b, 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$_{\oplus}$ for the minimum mass., Comment: 18 pages, 11 figures. Accepted version, A&A (2023)

We present the discovery of an Earth-mass planet ($M_b\sin i = 1.26\pm0.21M_\oplus$) on a 15.6d orbit of a relatively nearby ($d\sim$9.6pc) and low-mass ($0.167\pm0.011 M_\odot$) M5.0V star, Wolf 1069. Sitting at a separation of $0.0672\pm0.0014$au away from the host star puts Wolf 1069b in the habitable zone (HZ), receiving an incident flux of $S=0.652\pm0.029S_\oplus$. The planetary signal was detected using telluric-corrected radial-velocity (RV) data from the CARMENES spectrograph, amounting to a total of 262 spectroscopic observations covering almost four years. There are additional long-period signals in the RVs, one of which we attribute to the stellar rotation period. This is possible thanks to our photometric analysis including new, well-sampled monitoring campaigns undergone with the OSN and TJO facilities that supplement archival photometry (i.e., from MEarth and SuperWASP), and this yielded an updated rotational period range of $P_{rot}=150-170$d, with a likely value at $169.3^{+3.7}_{-3.6}$d. The stellar activity indicators provided by the CARMENES spectra likewise demonstrate evidence for the slow rotation period, though not as accurately due to possible factors such as signal aliasing or spot evolution. Our detectability limits indicate that additional planets more massive than one Earth mass with orbital periods of less than 10 days can be ruled out, suggesting that perhaps Wolf 1069 b had a violent formation history. This planet is also the 6th closest Earth-mass planet situated in the conservative HZ, after Proxima Centauri b, GJ 1061d, Teegarden's Star c, and GJ 1002 b and c. Despite not transiting, Wolf 1069b is nonetheless a very promising target for future three-dimensional climate models to investigate various habitability cases as well as for sub-ms$^{-1}$ RV campaigns to search for potential inner sub-Earth-mass planets in order to test planet formation theories., Comment: 26 pages, 15 figures

We report the discovery and characterisation of two Earth-mass planets orbiting in the habitable zone of the nearby M-dwarf GJ~1002 based on the analysis of the radial-velocity (RV) time series from the ESPRESSO and CARMENES spectrographs. The host star is the quiet M5.5~V star GJ~1002 (relatively faint in the optical, $V \sim 13.8$ mag, but brighter in the infrared, $J \sim 8.3$ mag), located at 4.84 pc from the Sun. We analyse 139 spectroscopic observations taken between 2017 and 2021. We performed a joint analysis of the time series of the RV and full-width half maximum (FWHM) of the cross-correlation function (CCF) to model the planetary and stellar signals present in the data, applying Gaussian process regression to deal with the stellar activity. We detect the signal of two planets orbiting GJ~1002. GJ~1002~b is a planet with a minimum mass $m_p \sin i $ of 1.08 $\pm$ 0.13 M$_{\oplus}$ with an orbital period of 10.3465 $\pm$ 0.0027 days at a distance of 0.0457 $\pm$ 0.0013 au from its parent star, receiving an estimated stellar flux of 0.67 $F_{\oplus}$. GJ~1002 c is a planet with a minimum mass $m_p \sin i $ of 1.36 $\pm$ 0.17 M$_{\oplus}$ with an orbital period of 21.202 $\pm$ 0.013 days at a distance of 0.0738 $\pm$ 0.0021 au from its parent star, receiving an estimated stellar flux of 0.257 $F_{\oplus}$. We also detect the rotation signature of the star, with a period of 126 $\pm$ 15 days. GJ~1002 is one of the few known nearby systems with planets that could potentially host habitable environments. The closeness of the host star to the Sun makes the angular sizes of the orbits of both planets ($\sim$ 9.7 mas and $\sim$ 15.7 mas, respectively) large enough for their atmosphere to be studied via high-contrast high-resolution spectroscopy with instruments such as the future spectrograph ANDES for the ELT or the LIFE mission., Comment: 23 pages, 20 figures, 3 tables ; Replacement to fix some typos

It is clearly established that the Sun has an 11-year cycle that is caused by its internal magnetic field. This cycle is also observed in a sample of M dwarfs. In the framework of exoplanet detection or atmospheric characterisation of exoplanets, the activity status of the host star plays a crucial role, and inactive states are preferable for such studies. This means that it is important to know the activity cycles of these stars. We study systematic long-term variability in a sample of 211 M dwarfs observed with CARMENES, the high-resolution optical and near-infrared spectrograph at Calar Alto Observatory. In an automatic search using time series of different activity indicators, we identified 26 stars with linear or quadratic trends or with potentially cyclic behaviour. Additionally, we performed an independent search in archival R$^{\prime}_{\rm HK}$ data collected from different instruments whose time baselines were usually much longer. These data are available for a subset of 186 of our sample stars. Our search revealed 22 cycle candidates in the data. We found that the percentage of stars showing long-term variations drops dramatically to the latest M dwarfs. Moreover, we found that the pseudo-equivalent width (pEW) of the H$\alpha$ and Ca ii infrared triplet more often triggers automatic detections of long-term variations than the TiO index, differential line width, chromatic index, or radial velocity. This is in line with our comparison of the median relative amplitudes of the different indicators. For stars that trigger our automatic detection, this leads to the highest amplitude variation in R$^{\prime}_{\rm HK}$, followed by pEW(H$\alpha$), pEW(Ca ii IRT), and the TiO index., Comment: Accepted to A&A, 16 pages, 12 figures

Ultra-hot Jupiters are highly irradiated gas giant exoplanets on close-in orbits around their host stars. We analyzed high-resolution spectra from CARMENES, HARPS-N, and ESPaDOnS taken over eight observation nights to study the emission spectrum of WASP-33b and draw conclusions about its atmosphere. By applying the cross-correlation technique, we detected the spectral signatures of Ti I, V I, and a tentative signal of Ti II for the first time via emission spectroscopy. These detections are an important finding because of the fundamental role of Ti- and V-bearing species in the planetary energy balance. Moreover, we assessed and confirm the presence of OH, Fe I, and Si I from previous studies. The spectral lines are all detected in emission, which unambiguously proves the presence of an inverted temperature profile in the planetary atmosphere. By performing retrievals on the emission lines of all the detected species, we determined a relatively weak atmospheric thermal inversion extending from approximately 3400 K to 4000 K. We infer a supersolar metallicity close to 1.5 dex in the planetary atmosphere, and find that its emission signature undergoes significant line broadening with a Gaussian FWHM of about 4.5 km/s. Also, we find that the atmospheric temperature profile retrieved at orbital phases far from the secondary eclipse is about 300 K to 700 K cooler than that measured close to the secondary eclipse, which is consistent with different day- and nightside temperatures. Moreover, retrievals performed on the emission lines of the individual chemical species lead to consistent results, which gives additional confidence to our retrieval method. Increasing the number of species included in the retrieval and expanding the set of retrieved atmospheric parameters will further advance our understanding of exoplanet atmospheres., Comment: Accepted for publication in A&A

Context: A challenge with radial-velocity (RV) data is disentangling the origin of signals either due to a planetary companion or to stellar activity. In fact, the existence of a planetary companion has been proposed, as well as contested, around the relatively bright, nearby M3.0V star AD Leo at the same period as the stellar rotation of 2.23d. Aims: We further investigate the nature of this signal. We introduce new CARMENES optical and near-IR RV data and an analysis in combination with archival data taken by HIRES and HARPS, along with more recent data from HARPS-N, GIANO-B, and HPF. Also, we address the confusion concerning the binarity of AD Leo. Methods: We consider possible correlations between the RVs and various stellar activity indicators accessible with CARMENES. We applied models within a Bayesian framework to determine whether a Keplerian model, a red-noise quasi-periodic model using a Gaussian process, or a mixed model would explain the observed data best. We also exclusively focus on spectral lines potentially associated with stellar activity. Results: The CARMENES RV data agree with the previously reported periodicity of 2.23d, correlate with some activity indicators, and exhibit chromaticity. However, when considering the entire RV data set, we find that a mixed model composed of a stable and a variable component performs best. Moreover, when recomputing the RVs using only spectral lines insensitive to activity, there appears to be some residual power at the period of interest. We therefore conclude that it is not possible to determinedly prove that there is no planet orbiting in synchronization with the stellar rotation given our data, current tools, machinery, and knowledge of how stellar activity affects RVs. We do rule out planets more massive than 27M_E (=0.084M_J). We also exclude any binary companion around AD Leo with Msini > 3-6M_J on orbital periods <14yr., Comment: 27 pages, 11 figures

We report the discovery and characterization of two small transiting planets orbiting the bright M3.0V star TOI-1468 (LSPM J0106+1913), whose transit signals were detected in the photometric time series in three sectors of the TESS mission. We confirm the e planetary nature of both of them using precise radial velocity measurements from the CARMENES and MAROON-X spectrographs, and supplement them with ground-based transit photometry. A joint analysis of all these data reveals that the shorter-period planet, TOI-1468 b ($P_{\rm b}$ = 1.88 d), has a planetary mass of $M_{\rm b} = 3.21\pm0.24$ $M_{\oplus}$ and a radius of $R_{\rm b} =1.280^{+0.038}_{-0.039} R_{\oplus}$, resulting in a density of $\rho_{\rm b} = 8.39^{+ 1.05}_{- 0.92}$ g cm$^{-3}$, which is consistent with a mostly rocky composition. For the outer planet, TOI-1468 c ($P_{\rm c} = 15.53$ d), we derive a mass of $M_{\rm c} = 6.64^{+ 0.67}_{- 0.68}$ $M_{\oplus}$, a radius of $R_{\rm c} = 2.06\pm0.04\,R_{\oplus}$, and a bulk density of $\rho_{c} = 2.00^{+ 0.21}_{- 0.19}$ g cm$^{-3}$, which corresponds to a rocky core composition with a H/He gas envelope. These planets are located on opposite sides of the radius valley, making our system an interesting discovery as there are only a handful of other systems with the same properties. This discovery can further help determine a more precise location of the radius valley for small planets around M dwarfs and, therefore, shed more light on planet formation and evolution scenarios., Comment: 26 pages, 15 figures

The Gaia-ESO Public Spectroscopic Survey is an ambitious project designed to obtain astrophysical parameters and elemental abundances for 100,000 stars, including large representative samples of the stellar populations in the Galaxy, and a well-defined sample of 60 (plus 20 archive) open clusters. We provide internally consistent results calibrated on benchmark stars and star clusters, extending across a very wide range of abundances and ages. This provides a legacy data set of intrinsic value, and equally a large wide-ranging dataset that is of value for homogenisation of other and future stellar surveys and Gaia's astrophysical parameters. This article provides an overview of the survey methodology, the scientific aims, and the implementation, including a description of the data processing for the GIRAFFE spectra. A companion paper (arXiv:2206.02901) introduces the survey results. Gaia-ESO aspires to quantify both random and systematic contributions to measurement uncertainties. Thus all available spectroscopic analysis techniques are utilised, each spectrum being analysed by up to several different analysis pipelines, with considerable effort being made to homogenise and calibrate the resulting parameters. We describe here the sequence of activities up to delivery of processed data products to the ESO Science Archive Facility for open use. The Gaia-ESO Survey obtained 202,000 spectra of 115,000 stars using 340 allocated VLT nights between December 2011 and January 2018 from GIRAFFE and UVES. The full consistently reduced final data set of spectra was released through the ESO Science Archive Facility in late 2020, with the full astrophysical parameters sets following in 2022., Comment: 38 pages. A&A in press

We report on the high-resolution spectroscopic observations of two planetary transits of the hot Jupiter KELT-7b ($M_{\rm p}$ $=$ 1.28 $\pm$ 0.17 M$_{\rm Jup}$, $T_{\rm eq}$ $=$ 2028 K) observed with the High Optical Resolution Spectrograph (HORuS) mounted on the 10.4-m Gran Telescopio Canarias (GTC). A new set of stellar parameters are obtained for the rapidly rotating parent star from the analysis of the spectra. Using the newly derived stellar mass and radius, and the planetary transit data of the Transiting Exoplanet Survey Satellite (TESS) together with the HORuS velocities and the photometric and spectroscopic data available in the literature, we update and improve the ephemeris of KELT-7b. Our results indicate that KELT-7 has an angle $\lambda$ = $-$10.55 $\pm$ 0.27 deg between the sky projections of the star's spin axis and the planet's orbital axis. By combining this angle and our newly derived stellar rotation period of 1.38 $\pm$ 0.05 d, we obtained a 3D obliquity $\psi$ = 12.4 $\pm$ 11.7 deg (or 167.6 deg), thus reinforcing that KELT-7 is a well-aligned planetary system. We search for the presence of H$\alpha$, Li I, Na I, Mg I, and Ca II features in the transmission spectrum of KELT-7b but we are only able to determine upper limits of 0.08-1.4 % on their presence after accounting for the contribution of the stellar variability to the extracted planetary spectrum. We also discuss the impact of stellar variability in the planetary data. Our results reinforce the importance of monitoring the parent star when performing high-resolution transmission spectroscopy of the planetary atmosphere in the presence of stellar activity., Comment: Accepted for publication in MNRAS

The Gl 486 system consists of a very nearby, relatively bright, weakly active M3.5 V star at just 8 pc with a warm transiting rocky planet of about 1.3 R_Terra and 3.0 M_Terra that is ideal for both transmission and emission spectroscopy and for testing interior models of telluric planets. To prepare for future studies, we collected light curves of seven new transits observed with the CHEOPS space mission and new radial velocities obtained with MAROON-X/Gemini North and CARMENES/Calar Alto telescopes, together with previously published spectroscopic and photometric data from the two spectrographs and TESS. We also performed interferometric observations with the CHARA Array and new photometric monitoring with a suite of smaller telescopes. From interferometry, we measure a limb-darkened disc angular size of the star Gl 486. Together with a corrected Gaia EDR3 parallax, we obtain a stellar radius. We also measure a stellar rotation period at P_rot ~ 49.9 d, an upper limit to its XUV (5-920 AA) flux with new Hubble/STIS data, and, for the first time, a variety of element abundances (Fe, Mg, Si, V, Sr, Zr, Rb) and C/O ratio. Besides, we impose restrictive constraints on the presence of additional components, either stellar or substellar, in the system. With the input stellar parameters and the radial-velocity and transit data, we determine the radius and mass of the planet Gl 486 b at R_p = 1.343+/0.063 R_Terra and M_p = 3.00+/-0.13 M_Terra. From the planet parameters and the stellar element abundances, we infer the most probable models of planet internal structure and composition, which are consistent with a relatively small metallic core with respect to the Earth, a deep silicate mantle, and a thin volatile upper layer. With all these ingredients, we outline prospects for Gl 486 b atmospheric studies, especially with forthcoming James Webb Space Telescope observations (abridged)., Comment: A&A, in press. See https://carmenes.caha.es/ext/pressreleases/GJ486/Exoearth.Video.mp4

In the last 15 years different ground-based spectroscopic surveys have been started (and completed) with the general aim of delivering stellar parameters and elemental abundances for large samples of Galactic stars, complementing Gaia astrometry. Among those surveys, the Gaia-ESO Public Spectroscopic Survey (GES), the only one performed on a 8m class telescope, was designed to target 100,000 stars using FLAMES on the ESO VLT (both Giraffe and UVES spectrographs), covering all the Milky Way populations, with a special focus on open star clusters. This article provides an overview of the survey implementation (observations, data quality, analysis and its success, data products, and releases), of the open cluster survey, of the science results and potential, and of the survey legacy. A companion article (Gilmore et al.) reviews the overall survey motivation, strategy, Giraffe pipeline data reduction, organisation, and workflow. The GES has determined homogeneous good-quality radial velocities and stellar parameters for a large fraction of its more than 110,000 unique target stars. Elemental abundances were derived for up to 31 elements for targets observed with UVES. Lithium abundances are delivered for about 1/3 of the sample. The analysis and homogenisation strategies have proven to be successful; several science topics have been addressed by the Gaia-ESO consortium and the community, with many highlight results achieved. The final catalogue has been released through the ESO archive at the end of May 2022, including the complete set of advanced data products. In addition to these results, the Gaia-ESO Survey will leave a very important legacy, for several aspects and for many years to come., Comment: Accepted for publication in Astronomy and Astrophysics. 30 pages, 30 figures, 4 tables

The CARMENES instrument is searching for periodic radial-velocity (RV) variations of M dwarfs, which are induced by orbiting planets. However, there are other potential sources of such variations, including rotational modulation caused by stellar activity. We aim to investigate four M dwarfs (Ross 318, YZ CMi, TYC 3529-1437-1, and EV Lac) with different activity levels and spectral sub-types. Our goal is to compare the periodicities seen in 22 activity indicators and the stellar RVs, and to examine their stability over time. For each star, we calculated GLS periodograms of pseudo-equivalent widths of chromospheric lines, indices of photospheric bands, the differential line width as a measure of the width of the average photospheric absorption line, the RV, the chromatic index that describes the wavelength dependence of the RV, and CCF parameters. We also calculated periodograms for subsets of the data and compared our results to TESS photometry. We find the rotation periods of all four stars to manifest themselves in the RV and photospheric indicators, particularly the TiO~7050 index, whereas the chromospheric lines show clear signals only at lower activity levels. For EV Lac and TYC 3529-1437-1, we find episodes during which indicators vary with the rotation period, and episodes during which they vary with half the rotation period, similarly to photometric light curves. The changing periodicities reflect the evolution of stellar activity features on the stellar surface. We therefore conclude that our results not only emphasise the importance of carefully analysing indicators complementary to the RV in RV surveys, but they also suggest that it is also useful to search for signals in activity indicators in subsets of the dataset, because an activity signal that is present in the RV may not be visible in the activity indicators all the time, in particular for the most active stars., Comment: 8 pages + 9 pages appendix, 5+4 figures, 1+4 tables, accepted for publication in A&A

We report the discovery of a multi-planetary system transiting the M0 V dwarf HD 260655 (GJ 239, TOI-4599). The system consists of at least two transiting planets, namely HD 260655 b, with a period of 2.77 d, a radius of R$_b$ = 1.240$\pm$0.023 R$_\oplus$, a mass of M$_b$ = 2.14$\pm$0.34 M$_\oplus$, and a bulk density of $\rho_b$ = 6.2$\pm$1.0 g cm$^{-3}$, and HD 260655 c, with a period of 5.71 d, a radius of R$_c$ = 1.533$^{+0.051}_{-0.046}$ R$_\oplus$, a mass of M$_c$ = 3.09$\pm$0.48 M$_\oplus$, and a bulk density of $\rho_c$ = 4.7$^{+0.9}_{-0.8}$ g cm$^{-3}$. The planets were detected in transit by the TESS mission and confirmed independently with archival and new precise radial velocities obtained with the HIRES and CARMENES instruments since 1998 and 2016, respectively. At a distance of 10 pc, HD 260655 becomes the fourth closest known multi-transiting planet system after HD 219134, LTT 1445 A, and AU Mic. Due to the apparent brightness of the host star (J = 6.7 mag), both planets are among the most suitable rocky worlds known today for atmospheric studies with the JWST, both in transmission and emission., Comment: 24 pages, 17 figures. Accepted for publication in A&A

We investigated the presence of planetary companions around the nearby (7.6 pc) and bright ($V=9$ mag) early-type M dwarf Gl 514, analysing 540 radial velocities collected over nearly 25 years with the HIRES, HARPS, and CARMENES spectrographs. The data are affected by time-correlated signals at the level of 2-3 ms$^{-1}$ due to stellar activity, that we filtered out testing three different models based on Gaussian process regression. As a sanity cross-check, we repeated the analyses using HARPS radial velocities extracted with three different algorithms. We used HIRES radial velocities and Hipparcos-Gaia astrometry to put constraints on the presence of long-period companions, and we analysed TESS photometric data. We found strong evidence that Gl 514 hosts a super-Earth on a likely eccentric orbit, residing in the conservative habitable zone for nearly $34\%$ of its orbital period. The planet Gl 514 b has minimum mass $m_b\sin i_b=5.2\pm0.9$ $M_{\rm Earth}$, orbital period $P_b=140.43\pm0.41$ days, and eccentricity $e_b=0.45^{+0.15}_{-0.14}$. No evidence for transits is found in the TESS light curve. There is no evidence for a longer period companion in the radial velocities and, based on astrometry, we can rule out a $\sim0.2$ $M_{\rm Jup}$ planet at a distance of $\sim3-10$ au, and massive giant planets/brown dwarfs out to several tens of au. We discuss the possible presence of a second low-mass companion at a shorter distance from the host than Gl 514 b. Gl 514 b represents an interesting science case to study the habitability of planets on eccentric orbits. We advocate for additional spectroscopic follow-up to get more accurate and precise planetary parameters. Further follow-up is also needed to investigate sub \ms and shorter period signals., Comment: 34 pages, accepted for publication in Astronomy & Astrophysics

Stellar dynamos generate magnetic fields that are of fundamental importance to the variability and evolution of Sun-like and low-mass stars, and for the development of their planetary systems. We report measurements of surface-average magnetic fields in 292 M dwarfs from a comparison with radiative transfer calculations; for 260 of them, this is the first measurement of this kind. Our data were obtained from more than 15,000 high-resolution spectra taken during the CARMENES project. They reveal a relation between average field strength, , and Rossby number, $Ro$, resembling the well-studied rotation-activity relation. Among the slowly rotating stars, we find that magnetic flux, $\Phi_\textrm{B}$, is proportional to rotation period, $P$, and among the rapidly rotating stars that average surface fields do not grow significantly beyond the level set by the available kinetic energy. Furthermore, we find close relations between nonthermal coronal X-ray emission, chromospheric H$\alpha$ and Ca H&K emission, and magnetic flux. Taken together, these relations demonstrate empirically that the rotation-activity relation can be traced back to a dependence of the magnetic dynamo on rotation. We advocate the picture that the magnetic dynamo generates magnetic flux on the stellar surface proportional to rotation rate with a saturation limit set by the available kinetic energy, and we provide relations for average field strengths and nonthermal emission that are independent of the choice of the convective turnover time. We also find that Ca H&K emission saturates at average field strengths of $\langle B \rangle \approx 800$ G while H$\alpha$ and X-ray emission grow further with stronger fields in the more rapidly rotating stars. This is in conflict with the coronal stripping scenario predicting that in the most rapidly rotating stars coronal plasma would be cooled to chromospheric temperatures., Comment: A&A accepted, abstract abbreviated, 20 pages, Table B.1 included in arXiv version, full MCMC posteriors and linefit plots are available at http://carmenes.cab.inta-csic.es/

The CARMENES radial-velocity survey is currently searching for planets in a sample of 387 M dwarfs. Here we report on two Saturn-mass planets orbiting TYC 2187-512-1 ($M_\star = 0.50 M_\odot$) and TZ Ari ($M_\star = 0.15 M_\odot$), respectively. We obtained supplementary photometric time series, which we use along with spectroscopic information to determine the rotation periods of the two stars. In both cases, the radial velocities also show strong modulations at the respective rotation period. We thus modeled the radial velocities as a Keplerian orbit plus a Gaussian process representing the stellar variability. TYC 2187-512-1 is found to harbor a planet with a minimum mass of 0.33 $M_{\rm Jup}$ in a near-circular 692-day orbit. The companion of TZ Ari has a minimum mass of 0.21 $M_{\rm Jup}$, orbital period of 771 d, and orbital eccentricity of 0.46. We provide an overview of all known giant planets in the CARMENES sample, from which we infer an occurrence rate of giant planets orbiting M dwarfs with periods up to 2 years in the range between 2% to 6%. TZ Ari b is only the second giant planet discovered orbiting a host with mass less than 0.3 $M_\odot$. These objects occupy an extreme location in the planet mass versus host mass plane. It is difficult to explain their formation in core-accretion scenarios, so they may possibly have been formed through a disk fragmentation process., Comment: Accepted for publication in Astronomy and Astrophysics

Current exoplanet surveys using the radial velocity (RV) technique are targeting M dwarfs because any habitable zone terrestrial-mass planets will induce a high RV and orbit on shorter periods than for more massive stars. One of the main caveats is that M dwarfs show a wide range of activity levels from inactive to very active, which can induce an asymmetry in the line profiles and, consequently, a spurious RV measurement. We aim to benchmark the impact of stellar activity on high-precision RV measurements using regular-cadence CARMENES visible and near-infrared observations of the active M3.5 dwarf EV Lac. We used the newly developed technique of low-resolution Doppler imaging to determine the centre-of-light, or spot-induced RV component, for eight observational epochs. We confirm a statistically significant and strong correlation between the independently measured centre-of-light and the chromatic index, which is a measure of the amplitude variation with wavelength of the RVs. We also find circular closed-loop relations of several activity indices with RV for a subset of data that covers only several rotation periods. We also investigate the implications of large phase gaps in the periodograms of activity indicators. Finally, by removing the spot-induced RV component we improve the planet-mass sensitivity by a factor of at least three. We conclude that for active M stars, a regular-cadence observing strategy is the most efficient way to identify and eliminate sources of correlated noise., Comment: accepted by A&A

The last decade has seen a revolution in our knowledge of the Galaxy thanks to the Gaia and asteroseismic space missions and the ground-based spectroscopic surveys. To complete this picture, it is necessary to map the ages of its stellar populations. During recent years, the dependence on time of abundance ratios involving slow (s) neutron-capture and $\alpha$ elements (called chemical-clocks) has been used to provide estimates of stellar ages, usually in a limited volume close to the Sun. We aim to analyse the relations of chemical clocks in the Galactic disc extending the range to R$_{\rm GC}\sim$6-20~kpc. Using the sixth internal data release of the Gaia-ESO survey, we calibrated several relations between stellar ages and abundance ratios [s/$\alpha$] using a sample of open clusters, the largest one so far used with this aim. Thanks to their wide galactocentric coverage, we investigated the radial variations of the shape of these relations, confirming their non-universality. We estimated our accuracy and precision in recovering the global ages of open clusters, and the ages of their individual members. We applied the multi-variate relations with the highest correlation coefficients to the field star population. We confirm that there is no single age-chemical clock relationship valid for the whole disc, but that there is a dependence on the galactocentric position, which is related to the radial variation of the star formation history combined with the non-monotonic dependence on metallicity of the yields of the s-process elements from low- and intermediate-mass stars. Finally, the abundance ratios [Ba/$\alpha$] are more sensitive to age than those with [Y/$\alpha$] for young disc stars, and their slopes vary less with galactocentric distance., Comment: 14 pages, 10 figures + Appendix (3 tables and 2 figures)

We report the discovery of GJ 3929 b, a hot Earth-sized planet orbiting the nearby M3.5 V dwarf star, GJ 3929 (G 180--18, TOI-2013). Joint modelling of photometric observations from TESS sectors 24 and 25 together with 73 spectroscopic observations from CARMENES and follow-up transit observations from SAINT-EX, LCOGT, and OSN yields a planet radius of $R_b = 1.150 +/- 0.040$ R$_{earth}$, a mass of $M_b = 1.21 +/- 0.42$ M$_{earth}$, and an orbital period of $P_b = 2.6162745 +/- 0.0000030$ d. The resulting density of $\rho_b= 4.4 +/- 1.6$ g/cm$^{-3}$ is compatible with the Earth's mean density of about 5.5 g/cm$^{-3}$. Due to the apparent brightness of the host star (J=8.7 mag) and its small size, GJ 3929 b is a promising target for atmospheric characterisation with the JWST. Additionally, the radial velocity data show evidence for another planet candidate with $P_{[c]} = 14.303 +/- 0.035$ d, which is likely unrelated to the stellar rotation period, $P_{rot} = 122+/-13$ d, which we determined from archival HATNet and ASAS-SN photometry combined with newly obtained TJO data., Comment: 24 pages; accepted for publication in A&A

The $\ion{He}{I}$ 10833\,$\AA$ triplet is a powerful tool for characterising the upper atmosphere of exoplanets and tracing possible mass loss. Here, we analysed one transit of GJ\,1214\,b observed with the CARMENES high-resolution spectrograph to study its atmosphere via transmission spectroscopy around the $\ion{He}{I}$ triplet. Although previous studies using lower resolution instruments have reported non-detections of $\ion{He}{I}$ in the atmosphere of GJ\,1214\,b, we report here the first potential detection. We reconcile the conflicting results arguing that previous transit observations did not present good opportunities for the detection of $\ion{He}{I}$, due to telluric H$_2$O absorption and OH emission contamination.We simulated those earlier observations, and show evidence that the planetary signal was contaminated. From our single non-telluric-contaminated transit, we determined an excess absorption of 2.10$^{+0.45}_{-0.50}$\,\% (4.6\,$\sigma$) with a full width at half maximum (FWHM) of 1.30$^{+0.30}_{-0.25}$\,\AA. The detection of \ion{He}{I} is statistically significant at the 4.6\,$\sigma$ level, but repeatability of the detection could not be confirmed due to the availability of only one transit. By applying a hydrodynamical model and assuming an H/He composition of 98/2, we found that GJ\,1214\,b would undergo hydrodynamic escape in the photon-limited regime, losing its primary atmosphere with a mass-loss rate of (1.5--18)\,$\times$\,10$^{10}$\,g\,s$^{-1}$ and an outflow temperature in the range of 2900--4400\,K. Our $\ion{He}{I}$ excess absorption is the first tentative detection of a chemical species in the atmosphere of this benchmark sub-Neptune planet., Comment: Accepted for publication in A&A. 12 pages, 7 figures

Ultra-hot Jupiters (UHJs) are gas giants with very high equilibrium temperatures. In recent years, multiple chemical species, including various atoms and ions, have been discovered in their atmospheres. Most of these observations have been performed with transmission spectroscopy, although UHJs are also ideal targets for emission spectroscopy due to their strong thermal radiation. We present high-resolution thermal emission spectroscopy of the transiting UHJ KELT-20b/MASCARA-2b. The observation was performed with the CARMENES spectrograph at orbital phases before and after the secondary eclipse. We detected atomic Fe using the cross-correlation technique. The detected Fe lines are in emission, which unambiguously indicates a temperature inversion on the dayside hemisphere. We furthermore retrieved the temperature structure with the detected Fe lines. The result shows that the atmosphere has a strong temperature inversion with a temperature of $4900\pm{700}$ K and a pressure of $10^{-4.8_{-1.1}^{+1.0}}$ bar at the upper layer of the inversion. A joint retrieval of the CARMENES data and the TESS secondary eclipse data returns a temperature of $2550_{-250}^{+150}$ K and a pressure of $10^{-1.5_{-0.6}^{+0.7}}$ bar at the lower layer of the temperature inversion. The detection of such a strong temperature inversion is consistent with theoretical simulations that predict an inversion layer on the dayside of UHJs. The joint retrieval of the CARMENES and TESS data demonstrates the power of combing high-resolution emission spectroscopy with secondary eclipse photometry in characterizing atmospheric temperature structures., Comment: 14 pages, 13 figures; accepted for publication in A&A

Atmospheres of highly irradiated gas giant planets host a large variety of atomic and ionic species. Here we observe the thermal emission spectra of the two ultra-hot Jupiters WASP-33b and KELT-20b/MASCARA-2b in the near-infrared wavelength range with CARMENES. Via high-resolution Doppler spectroscopy, we searched for neutral silicon (Si) in their dayside atmospheres. We detect the Si spectral signature of both planets via cross-correlation with model spectra. Detection levels of 4.8$\mathrm{\sigma}$ and 5.4$\mathrm{\sigma}$, respectively, are observed when assuming a solar atmospheric composition. This is the first detection of Si in exoplanet atmospheres. The presence of Si is an important finding due to its fundamental role in cloud formation and, hence, for the planetary energy balance. Since the spectral lines are detected in emission, our results also confirm the presence of an inverted temperature profile in the dayside atmospheres of both planets., Comment: Accepted for publication in A&A

Two transiting planet candidates with super-Earth radii around the nearby K7--M0 dwarf star TOI-1238 were announced by TESS. We aim to validate their planetary nature using precise radial velocities (RV) taken with the CARMENES spectrograph. We obtained 55 CARMENES RV data that span 11 months. For a better characterization of the parent star's activity, we also collected contemporaneous optical photometric observations and retrieved archival photometry from the literature. We performed a combined TESS+CARMENES photometric and spectroscopic analysis by including Gaussian processes and Keplerian orbits to account for the stellar activity and planetary signals simultaneously. We estimate that TOI-1238 has a rotation period of 40 $\pm$ 5 d based on photometric and spectroscopic data. The combined analysis confirms the discovery of two transiting planets, TOI-1238 b and c, with orbital periods of $0.764597^{+0.000013}_{-0.000011}$ d and $3.294736^{+0.000034}_{-0.000036}$ d, masses of 3.76$^{+1.15}_{-1.07}$ M$_{\oplus}$ and 8.32$^{+1.90}_{-1.88}$ M$_{\oplus}$, and radii of $1.21^{+0.11}_{-0.10}$ R$_{\oplus}$ and $2.11^{+0.14}_{-0.14}$ R$_{\oplus}$. They orbit their parent star at semimajor axes of 0.0137$\pm$0.0004 au and 0.036$\pm$0.001 au, respectively. The two planets are placed on opposite sides of the radius valley for M dwarfs and lie between the star and the inner border of TOI-1238's habitable zone. The inner super-Earth TOI-1238 b is one of the densest ultra-short-period planets ever discovered ($\rho=11.7^{+4.2}_{-3.4}$ g $\rm cm^{-3}$). The CARMENES data also reveal the presence of an outer, non-transiting, more massive companion with an orbital period and radial velocity amplitude of $\geq$600 d and $\geq$70 m s$^{-1}$, which implies a likely mass of $M \geq 2 \sqrt{1-e^2}$ M$_{\rm Jup}$ and a separation $\geq$1.1 au from its parent star.