Search

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

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 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 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

TOI-2525 is a K-type star with an estimated mass of M = 0.849$_{-0.033}^{+0.024}$ M$_\odot$ and radius of R = 0.785$_{-0.007}^{+0.007}$ R$_\odot$ observed by the TESS mission in 22 sectors (within sectors 1 and 39). The TESS light curves yield significant transit events of two companions, which show strong transit timing variations (TTVs) with a semi-amplitude of a $\sim$6 hours. We performed TTV dynamical, and photo-dynamical light curve analysis of the TESS data, combined with radial velocity (RV) measurements from FEROS and PFS, and we confirmed the planetary nature of these companions. The TOI-2525 system consists of a transiting pair of planets comparable to Neptune and Jupiter with estimated dynamical masses of $m_{\rm b}$ = 0.088$_{-0.004}^{+0.005}$ M$_{\rm Jup.}$, and $m_{\rm c}$ = 0.709$_{-0.033}^{+0.034}$ M$_{\rm Jup.}$, radius of $r_b$ = 0.88$_{-0.02}^{+0.02}$ R$_{\rm Jup.}$ and $r_c$ = 0.98$_{-0.02}^{+0.02}$ R$_{\rm Jup.}$, and with orbital periods of $P_{\rm b}$ = 23.288$_{-0.002}^{+0.001}$ days and $P_{\rm c}$ = 49.260$_{-0.001}^{+0.001}$ days for the inner and the outer planet, respectively. The period ratio is close to the 2:1 period commensurability, but the dynamical simulations of the system suggest that it is outside the mean motion resonance (MMR) dynamical configuration. TOI-2525 b is among the lowest density Neptune-mass planets known to date, with an estimated median density of $\rho_{\rm b}$ = 0.174$_{-0.015}^{+0.016}$ g\,cm$^{-3}$. The TOI-2525 system is very similar to the other K-dwarf systems discovered by TESS, TOI-2202 and TOI-216, which are composed of almost identical K-dwarf primary and two warm giant planets near the 2:1 MMR., Comment: Accepted for publication in AJ

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

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

Aiming to detect planetary companions to young stars with debris disks via the radial velocity method, we observed HD114082 during April 2018 - August 2022 as one of the targets of our RVSPY program (Radial Velocity Survey for Planets around Young stars). We used the FEROS spectrograph, mounted to the MPG/ESO 2.2 m telescope in Chile, to obtain high signal-to-noise spectra and time series of precise radial velocities (RVs). Additionally, we analyzed archival HARPS spectra and TESS photometric data. We used the CERES, CERES++ and SERVAL pipelines to derive RVs and activity indicators and ExoStriker for the independent and combined analysis of the RVs and TESS photometry. We report the discovery of a warm super-Jovian companion around HD114082 based on a 109.8$\pm$0.4 day signal in the combined RV data from FEROS and HARPS, and on one transit event in the TESS photometry. The best-fit model indicates a 8.0$\pm$1.0 Mjup companion with a radius of 1.00$\pm$0.03 Rjup in an orbit with a semi-major axis of 0.51$\pm$0.01 au and an eccentricity of 0.4$\pm$0.04. The companions orbit is in agreement with the known near edge-on debris disk located at about 28 au. HD114082b is possibly the youngest (15$\pm$6 Myr), and one of only three younger than 100 Myr giant planetary companions for which both their mass and radius have been determined observationally. It is probably the first properly model-constraining giant planet that allows distinguishing between hot and cold-start models. It is significantly more compatible with the cold-start model., Comment: 10 pages, 9 figures, 5 tables; Accepted for publication in A&A Letters

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 introduce our Radial Velocity Survey for Planets around Young stars (RVSPY), characterise our target stars, and search for substellar companions at orbital separations smaller than a few au from the host star. We use the FEROS spectrograph to obtain high signal-to-noise spectra and time series of precise radial velocities (RVs) of 111 stars most of which are surrounded by debris discs. Our target stars have spectral types between early F and late K, a median age of 400 Myr, and a median distance of 45 pc. We determine for all target stars their basic stellar parameters and present the results of the high-cadence RV survey and activity characterization. We achieve a median single-measurement RV precision of 6 m/s and derive the short-term intrinsic RV scatter of our targets (median 22 m/s), which is mostly caused by stellar activity and decays with age from >100 m/s at <20 Myr to <20 m/s at >500 Myr. We discover six previously unknown close companions with orbital periods between 10 and 100 days, three of which are low-mass stars, and three are in the brown dwarf mass regime. We detect no hot companion with an orbital period <10 days down to a median mass limit of ~1 M_Jup for stars younger than 500 Myr, which is still compatible with the established occurrence rate of such companions around main-sequence stars. We find significant RV periodicities between 1.3 and 4.5 days for 14 stars, which are, however, all caused by rotational modulation due to starspots. We also analyse the TESS photometric time series data and find significant periodicities for most of the stars. For 11 stars, the photometric periods are also clearly detected in the RV data. We also derive stellar rotation periods ranging from 1 to 10 days for 91 stars, mostly from TESS data. From the intrinsic activity-related short-term RV jitter, we derive the expected mass-detection thresholds for longer-period companions., Comment: 24 pages, 14 figures, 4 tables; Accepted for publication in A&A

Based on recently-taken and archival HARPS, FEROS and HIRES radial velocities (RVs), we present evidence for a new planet orbiting the first ascent red giant star HD33142 (with an improved mass estimate of 1.52$\pm$0.03 M$_\odot$), already known to host two planets. We confirm the Jovian mass planets HD33142 b and c with periods of $P_{\rm b}$ = 330.0$_{-0.4}^{+0.4}$ d and $P_{\rm c}$ = 810.2$_{-4.2}^{+3.8}$ d and minimum dynamical masses of $m_{\rm b}\sin{i}$ = 1.26$_{-0.05}^{+0.05}$ M$_{\rm Jup}$ and $m_{\rm c}\sin{i}$ = 0.89$_{-0.05}^{+0.06}$ M$_{\rm Jup}$. Furthermore, our periodogram analysis of the precise RVs shows strong evidence for a short-period Doppler signal in the residuals of a two-planet Keplerian fit, which we interpret as a third, Saturn-mass planet with $m_\mathrm{d}\sin{i}$ = 0.20$_{-0.03}^{+0.02}$ M$_{\rm Jup}$ on a close-in orbit with an orbital period of $P_{\rm d}$ =89.9$_{-0.1}^{+0.1}$ d. We study the dynamical behavior of the three-planet system configurations with an N-body integration scheme, finding it long-term stable with the planets alternating between low and moderate eccentricities episodes. We also performed N-body simulations, including stellar evolution and second-order dynamical effects such as planet-stellar tides and stellar mass-loss on the way to the white dwarf phase. We find that planets HD33142 b, c and d are likely to be engulfed near the tip of the red giant branch phase due to tidal migration. These results make the HD33142 system an essential benchmark for the planet population statistics of the multiple-planet systems found around evolved stars., Comment: Accepted for publication in AJ

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

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

We present an independent Doppler validation and dynamical orbital analysis of the two-planet system HD 107148, which was recently announced in Rosenthal et al. (2021). Our detailed analyses are based on literature HIRES data and newly obtained HARPS and CARMENES radial velocity (RV) measurements as part of our survey in search for additional planets around single planet systems. We perform a periodogram analysis of the available HIRES and HARPS precise RVs and stellar activity indicators. We do not find any apparent correlation between the RV measurements and the stellar activity indicators, thus linking the two strong periodicities to a moderately compact multiple-planet system. We carry out orbital fitting analysis by testing various one- and two-planet orbital configurations and studying the posterior probability distribution of the fitted parameters. Our results solidify the existence of a Saturn-mass planet (HD 107148 b, discovered first) with a period $P_b\sim77.2$ d, and a second, eccentric ($e_c \sim$ 0.4), Neptune-mass exoplanet (HD 107148 c), with an orbital period of $P_c\sim18.3$ d. Finally, we investigate the two-planet system's long-term stability and overall orbital dynamics with the posterior distribution of our preferred orbital configuration. Our N-body stability simulations show that the system is long-term stable and exhibits large secular osculations in eccentricity but in no particular mean-motion resonance configuration. The HD 107148 system, consisting of a Solar-type main sequence star with two giant planets in a rare configuration, features a common proper motion white dwarf companion and is, therefore, a valuable target for understanding the formation and evolution of planetary systems., Comment: Accepted for publication in AJ

We report the discovery and characterization of TOI-1759~b, a temperate (400 K) sub-Neptune-sized exoplanet orbiting the M~dwarf TOI-1759 (TIC 408636441). TOI-1759 b was observed by TESS to transit on sectors 16, 17 and 24, with only one transit observed per sector, creating an ambiguity on the orbital period of the planet candidate. Ground-based photometric observations, combined with radial-velocity measurements obtained with the CARMENES spectrograph, confirm an actual period of $18.85019 \pm 0.00014$ d. A joint analysis of all available photometry and radial velocities reveal a radius of $3.17 \pm 0.10\,R_\oplus$ and a mass of $10.8 \pm 1.5\,M_\oplus$. Combining this with the stellar properties derived for TOI-1759 ($R_\star = 0.597 \pm 0.015\,R_\odot$; $M_\star = 0.606 \pm 0.020\,M_\odot$; $T_{\textrm{eff}} = 4065 \pm 51$ K), we compute a transmission spectroscopic metric (TSM) value of over 80 for the planet, making it a good target for transmission spectroscopy studies. TOI-1759 b is among the top five temperate, small exoplanets ($T_\textrm{eq} < 500$ K, $R_p < 4 \,R_\oplus$) with the highest TSM discovered to date. Two additional signals with periods of 80 d and $>$ 200 d seem to be present in our radial velocities. While our data suggest both could arise from stellar activity, the later signal's source and periodicity are hard to pinpoint given the $\sim 200$ d baseline of our radial-velocity campaign with CARMENES. Longer baseline radial-velocity campaigns should be performed in order to unveil the true nature of this long period signal., Comment: 22 pages, 10 figures, 7 tables. AJ in press

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

There are several strong K I lines found in the spectra of M dwarfs, among them the doublet near 7700 AA and another doublet near 12 500 AA. We study these optical and near-infrared doublets in a sample of 324 M dwarfs, observed with CARMENES, the high-resolution optical and near-infrared spectrograph at Calar Alto, and investigate how well the lines can be used as photospheric and chromospheric diagnostics. Both doublets have a dominant photospheric component in inactive stars and can be used as tracers of effective temperature and gravity. For variability studies using the optical doublet, we concentrate on the red line component because this is less prone to artefacts from telluric correction in individual spectra. The optical doublet lines are sensitive to activity, especially for M dwarfs later than M5.0 V where the lines develop an emission core. For earlier type M dwarfs, the red component of the optical doublet lines is also correlated with H$\alpha$ activity. We usually find positive correlation for stars with H$\alpha$ in emission, while early-type M stars with H$\alpha$ in absorption show anti-correlation. During flares, the optical doublet lines can exhibit strong fill-in or emission cores for our latest spectral types. On the other hand, the near-infrared doublet lines very rarely show correlation or anti-correlation to H$\alpha$ and do not change line shape significantly even during the strongest observed flares. Nevertheless, the near-infrared doublet lines show notable resolved Zeeman splitting for about 20 active stars which allows to estimate the magnetic fields B., Comment: accepted to A&A

We present the discovery of a transiting mini-Neptune around TOI-1201, a relatively bright and moderately young early M dwarf ($J \approx$ 9.5 mag, $\sim$600-800 Myr) in an equal-mass $\sim$8 arcsecond-wide binary system, using data from the Transiting Exoplanet Survey Satellite (TESS), along with follow-up transit observations. With an orbital period of 2.49 d, TOI-1201 b is a warm mini-Neptune with a radius of $R_\mathrm{b} = 2.415\pm0.090 R_\oplus$. This signal is also present in the precise radial velocity measurements from CARMENES, confirming the existence of the planet and providing a planetary mass of $M_\mathrm{b} = 6.28\pm0.88 M_\oplus$ and, thus, an estimated bulk density of $2.45^{+0.48}_{-0.42}$ g cm$^{-3}$. The spectroscopic observations additionally show evidence of a signal with a period of 19 d and a long periodic variation of undetermined origin. In combination with ground-based photometric monitoring from WASP-South and ASAS-SN, we attribute the 19 d signal to the stellar rotation period ($P_{rot}=$ 19-23 d), although we cannot rule out that the variation seen in photometry belongs to the visually close binary companion. We calculate precise stellar parameters for both TOI-1201 and its companion. The transiting planet is an excellent target for atmosphere characterization (the transmission spectroscopy metric is $97^{+21}_{-16}$) with the upcoming James Webb Space Telescope. It is also feasible to measure its spin-orbit alignment via the Rossiter-McLaughlin effect using current state-of-the-art spectrographs with submeter per second radial velocity precision., Comment: 33 pages; 18 figures; accepted for publication in A&A

The CARMENES exoplanet survey of M dwarfs has obtained more than 18 000 spectra of 329 nearby M dwarfs over the past five years as part of its guaranteed time observations (GTO) program. We determine planet occurrence rates with the 71 stars from the GTO program for which we have more than 50 observations. We use injection-and-retrieval experiments on the radial-velocity (RV) time series to measure detection probabilities. We include 27 planets in 21 planetary systems in our analysis. We find 0.06+0.04-0.03 giant planets (100 M_Earth < M_pl sin i < 1000 M_Earth) per star in periods of up to 1000 d, but due to a selection bias this number could be up to a factor of five lower in the whole 329-star sample. The upper limit for hot Jupiters (orbital period of less than 10 d) is 0.03 planets per star, while the occurrence rate of planets with intermediate masses (10 M_Earth < M_pl sin i < 100 M_Earth) is 0.18+0.07-0.05 planets per star. Less massive planets with 1 M_Earth < M_pl sin i < 10 M_Earth are very abundant, with an estimated rate of 1.32+0.33-0.31 planets per star for periods of up to 100 d. When considering only late M dwarfs with masses M_star < 0.34 M_sol, planets more massive than 10 M_Earth become rare. Instead, low-mass planets with periods shorter than 10 d are significantly overabundant. For orbital periods shorter than 100 d, our results confirm the known stellar mass dependences from the Kepler survey: M dwarfs host fewer giant planets and at least two times more planets with M_pl sin i < 10 M_Earth than G-type stars. In contrast to previous results, planets around our sample of very low-mass stars have a higher occurrence rate in short-period orbits of less than 10 d. Our results demonstrate the need to take into account host star masses in planet formation models., Comment: 15 pages, 12 figures. Accepted for publication in Astronomy & Astrophysics

M dwarfs are ideal targets for the search of Earth-size planets in the habitable zone using the radial velocity method, attracting the attention of many ongoing surveys. As a by-product of these surveys, new multiple stellar systems are also found. This is the case also for the CARMENES survey, from which nine new SB2 systems have already been announced. Throughout the five years of the survey, the accumulation of new observations has resulted in the detection of several new multiple stellar systems with long periods and low radial-velocity amplitudes. Here, we newly characterise the spectroscopic orbits and constrain the masses of eight systems and update the properties of a system that we reported earlier. We derive the radial velocities of the stars using two-dimensional cross correlation techniques and template matching. The measurements are modelled to determine the orbital parameters of the systems. We combine CARMENES spectroscopic observations with archival high-resolution spectra from other instruments to increase the time-span of the observations and improve our analysis. When available, we also added archival photometric, astrometric, and adaptive optics imaging data to constrain the rotation periods and absolute masses of the components. We determine the spectroscopic orbits of nine multiple systems, eight of which are presented for the first time. The sample is composed of five SB1s, two SB2s, and two ST3s. The companions of two of the single-line binaries, GJ 3626 and GJ 912, have minimum masses below the stellar boundary and, thus, could be brown dwarfs. We find a new white dwarf in a close binary orbit around the M star GJ 207.1. From a global fit to radial velocities and astrometric measurements, we are able to determine the absolute masses of the components of GJ 282C, which is one of the youngest systems with measured dynamical masses., Comment: Accepted for publication in A&A

We report the discovery of two planetary systems, namely G 264-012, an M4.0 dwarf with two terrestrial planets ($M_{\rm b}\sin{i} = 2.50^{+0.29}_{-0.30}$ M$_{\oplus}$ and $M_{\rm c}\sin{i} = 3.75^{+0.48}_{-0.47}$ M$_{\oplus}$), and Gl 393, a bright M2.0 dwarf with one terrestrial planet ($M_{\rm b}\sin{i} = 1.71 \pm 0.24$ M$_{\oplus}$). Although both stars were proposed to belong to young stellar kinematic groups, we estimate their ages to be older than about 700 Ma. The two planets around G 264-012 were discovered using only radial-velocity (RV) data from the CARMENES exoplanet survey, with estimated orbital periods of $2.30$ d and $8.05$ d, respectively. Photometric monitoring and analysis of activity indicators reveal a third signal present in the RV measurements, at about 100 d, caused by stellar rotation. The planet Gl 393 b was discovered in the RV data from the HARPS, CARMENES, and HIRES instruments. Its identification was only possible after modelling, with a Gaussian process (GP), the variability produced by the magnetic activity of the star. For the earliest observations, this variability produced a forest of peaks in the periodogram of the RVs at around the 34 d rotation period determined from {\em Kepler} data, which disappeared in the latest epochs. After correcting for them with this GP model, a significant signal showed at a period of $7.03$ d. No significant signals in any of our spectral activity indicators or contemporaneous photometry were found at any of the planetary periods. Given the orbital and stellar properties, the equilibrium temperatures of the three planets are all higher than that for Earth. Current planet formation theories suggest that these two systems represent a common type of architecture. This is consistent with formation following the core accretion paradigm., Comment: 24 pages, 18 figures. Paper in press in A&A

Stellar activity poses one of the main obstacles for the detection and characterisation of small exoplanets around cool stars, as it can induce radial velocity (RV) signals that can hide or mimic the presence of planetary companions. Several indicators of stellar activity are routinely used to identify activity-related signals in RVs, but not all indicators trace exactly the same activity effects, nor are any of them always effective in all stars. We evaluate the performance of a set of spectroscopic activity indicators for M dwarf stars with different masses and activity levels with the aim of finding a relation between the indicators and stellar properties. In a sample of 98 M dwarfs observed with CARMENES, we analyse the temporal behaviour of RVs and nine spectroscopic activity indicators: cross-correlation function (CCF) full-width-at-half-maximum (FWHM), contrast, and bisector inverse slope (BIS), chromatic index (CRX), differential line width (dLW), and indices of the chromospheric lines H$\alpha$ and calcium infrared triplet. A total of 56 stars of the initial sample show periodic signals related to activity in at least one of these ten parameters. RV is the parameter for which most of the targets show an activity-related signal. CRX and BIS are effective activity tracers for the most active stars in the sample, especially stars with a relatively high mass, while for less active stars, chromospheric lines perform best. FWHM and dLW show a similar behaviour in all mass and activity regimes, with the highest number of activity detections in the low-mass, high-activity regime. Most of the targets for which we cannot identify any activity-related signals are stars at the low-mass end of the sample. These low-mass stars also show the lowest RV scatter, which indicates that ultracool M dwarfs could be better candidates for planet searches than earlier types, which show larger RV jitter., Comment: 16 pages, 8 figures, 1 table (37 pages, 14 figures, 12 tables with appendices). Accepted for publication in A&A

The interaction between Earth-like exoplanets and the magnetic field of low-mass host stars are considered to produce weak emission signals at radio frequencies. A study using LOFAR data announced the detection of radio emission from the mid M-type dwarf GJ 1151 that could potentially arise from a close-in terrestrial planet. Recently, the presence of a 2.5-Me planet orbiting GJ 1151 with a 2-day period has been claimed using 69 radial velocities (RVs) from the HARPS-N and HPF instruments. We have obtained 70 new high-precision RV measurements in the framework of the CARMENES M-dwarf survey and use these data to confirm the presence of the claimed planet and to place limits on possible planetary companions in the GJ 1151 system. We analyse the periodicities present in the combined RV data sets from all three instruments and calculate the detection limits for potential planets in short-period orbits. We cannot confirm the recently-announced candidate planet and conclude that the 2-day signal in the HARPS-N and HPF data sets is most probably produced by a long-term RV variability possibly arising from an outer planetary companion yet unconstrained. We calculate a 99.9% significance detection limit of 1.50 ms-1 in the RV semi-amplitude, which places upper limits of 0.7 Me and 1.2 Me to the minimum masses of any potential exoplanets with orbital periods of 1 and 5 days, respectively.

Spectroscopy of transiting exoplanets can be used to investigate their atmospheric properties and habitability. Combining radial velocity (RV) and transit data provides additional information on exoplanet physical properties. We detect a transiting rocky planet with an orbital period of 1.467 days around the nearby red dwarf star Gliese 486. The planet Gliese 486 b is 2.81 Earth masses and 1.31 Earth radii, with uncertainties of 5%, as determined from RV data and photometric light curves. The host star is at a distance of ~8.1 parsecs, has a J-band magnitude of ~7.2, and is observable from both hemispheres of Earth. On the basis of these properties and the planet's short orbital period and high equilibrium temperature, we show that this terrestrial planet is suitable for emission and transit spectroscopy., Comment: Published in Science

Dynamical histories of planetary systems, as well as atmospheric evolution of highly irradiated planets, can be studied by characterizing the ultra-short-period planet population, which the TESS mission is particularly well suited to discover. Here, we report on the follow-up of a transit signal detected in the TESS sector 19 photometric time series of the M3.0 V star TOI-1685 (2MASS J04342248+4302148). We confirm the planetary nature of the transit signal, which has a period of P_b=0.6691403+0.0000023-0.0000021 d, using precise radial velocity measurements taken with the CARMENES spectrograph. From the joint photometry and radial velocity analysis, we estimate the following parameters for TOI-1685 b: a mass of M_b=3.78+/-0.63 M_Earth, a radius of R_b=1.70+/-0.07 R_Earth, which together result in a bulk density of rho_b=4.21+0.95-0.82 g/cm3, and an equilibrium temperature of Teq_b=1069+/-16 K. TOI-1685 b is the least dense ultra-short period planet around an M dwarf known to date. TOI-1685 b is also one of the hottest transiting Earth-size planets with accurate dynamical mass measurements, which makes it a particularly attractive target for thermal emission spectroscopy. Additionally, we report a further non-transiting planet candidate in the system, TOI-1685[c], with an orbital period of P_[c]=9.02+0.10-0.12 d., Comment: 18 pages, 15 figures

We report the discovery of a Neptune-like planet (LP 714-47 b, P = 4.05204 d, m_b = 30.8 +/- 1.5 M_earth , R_b = 4.7 +/- 0.3 R_earth ) located in the 'hot Neptune desert'. Confirmation of the TESS Object of Interest (TOI 442.01) was achieved with radial-velocity follow-up using CARMENES, ESPRESSO, HIRES, iSHELL, and PFS, as well as from photometric data using TESS, Spitzer, and ground-based photometry from MuSCAT2, TRAPPIST- South, MONET-South, the George Mason University telescope, the Las Cumbres Observatory Global Telescope network, the El Sauce telescope, the TUBITAK National Observatory, the University of Louisville Manner Telescope, and WASP-South. We also present high-spatial resolution adaptive optics imaging with the Gemini Near-Infrared Imager. The low uncertainties in the mass and radius determination place LP 714-47 b among physically well-characterised planets, allowing for a meaningful comparison with planet structure models. The host star LP 714-47 is a slowly rotating early M dwarf (T_eff = 3950 +/- 51 K) with a mass of 0.59 +/- 0.02 M_sun and a radius of 0.58 +/- 0.02 R_sun. From long-term photometric monitoring and spectroscopic activity indicators, we determine a stellar rotation period of about 33 d. The stellar activity is also manifested as correlated noise in the radial-velocity data. In the power spectrum of the radial-velocity data, we detect a second signal with a period of 16 days in addition to the four-day signal of the planet. This could be shown to be a harmonic of the stellar rotation period or the signal of a second planet. It may be possible to tell the difference once more TESS data and radial-velocity data are obtained., Comment: Accepted for publication in A&A. 24 pages, 21 figures, 5 tables

We announce the discovery of two planets orbiting the M dwarfs GJ 251 ($0.360\pm0.015$ M$_\odot$) and HD 238090 ($0.578\pm0.021$ M$_\odot$) based on CARMENES radial velocity (RV) data. In addition, we independently confirm with CARMENES data the existence of Lalande 21185 b, a planet that has recently been discovered with the SOPHIE spectrograph. All three planets belong to the class of warm or temperate super-Earths and share similar properties. The orbital periods are 14.24 d, 13.67 d, and 12.95 d and the minimum masses are $4.0\pm0.4$ $M_\oplus$, $6.9\pm0.9$ $M_\oplus$, and $2.7\pm0.3$ $M_\oplus$ for GJ 251 b, HD 238090 b, and Lalande 21185 b, respectively. Based on the orbital and stellar properties, we estimate equilibrium temperatures of $351.0\pm1.4$ K for GJ 251 b, $469.6\pm2.6$ K for HD 238090 b, and $370.1\pm6.8$ K for Lalande 21185 b. For the latter we resolve the daily aliases that were present in the SOPHIE data and that hindered an unambiguous determination of the orbital period. We find no significant signals in any of our spectral activity indicators at the planetary periods. The RV observations were accompanied by contemporaneous photometric observations. We derive stellar rotation periods of $122.1\pm2.2$ d and $96.7\pm3.7$ d for GJ 251 and HD 238090, respectively. The RV data of all three stars exhibit significant signals at the rotational period or its first harmonic. For GJ 251 and Lalande 21185, we also find long-period signals around 600 d, and 2900 d, respectively, which we tentatively attribute to long-term magnetic cycles. We apply a Bayesian approach to carefully model the Keplerian signals simultaneously with the stellar activity using Gaussian process regression models and extensively search for additional significant planetary signals hidden behind the stellar activity., Comment: Accepted for publication in A&A

We present the confirmation and characterisation of GJ 3473 b (G 50--16, TOI-488.01), a hot Earth-sized planet orbiting an M4 dwarf star, whose transiting signal ($P=1.1980035\pm0.0000018\mathrm{\,d}$) was first detected by the Transiting Exoplanet Survey Satellite (TESS). Through a joint modelling of follow-up radial velocity observations with CARMENES, IRD, and HARPS together with extensive ground-based photometric follow-up observations with LCOGT, MuSCAT, and MuSCAT2, we determined a precise planetary mass, $M_b = 1.86\pm0.30\,\mathrm{M_\oplus},$ and radius, $R_b = {1.264\pm0.050}\,\mathrm{R_\oplus}$. Additionally, we report the discovery of a second, temperate, non-transiting planet in the system, GJ 3473 c, which has a minimum mass, $M_c \sin{i} = {7.41\pm0.91}\,\mathrm{M_\oplus,}$ and orbital period, $P_c={15.509\pm0.033}\,\mathrm{d}$. The inner planet of the system, GJ 3473 b, is one of the hottest transiting Earth-sized planets known thus far, accompanied by a dynamical mass measurement, which makes it a particularly attractive target for thermal emission spectroscopy., Comment: 21 pages; accepted for publication in A&A

In this study, abundances of the neutron-capture elements Rb, Sr, and Zr are derived, for the first time, in a sample of nearby M dwarfs. We focus on stars in the metallicity range -0.5<[Fe/H]<+0.3, an interval poorly explored for Rb abundances in previous analyses. To do this we use high-resolution, high-signal-to-noise-ratio, optical and near-infrared spectra of 57 M dwarfs observed with CARMENES. The resulting [Sr/Fe] and [Zr/Fe] ratios for most M dwarfs are almost constant at about the solar value, and are identical to those found in GK dwarfs of the same metallicity. However, for Rb we find systematic underabundances ([Rb/Fe]<0.0) by a factor two on average. Furthermore, a tendency is found for Rb-but not for other heavy elements (Sr, Zr) -to increase with increasing metallicity such that [Rb/Fe]>0.0 is attained at metallicities higher than solar. These are surprising results, never seen for any other heavy element, and are difficult to understand within the formulation of the s- and r-processes, both contributing sources to the Galactic Rb abundance. We discuss the reliability of these findings for Rb in terms of non-LTE effects, stellar activity, or an anomalous Rb abundance in the Solar System, but no explanation is found. We then interpret the full observed [Rb/Fe] versus [Fe/H] trend within the framework of theoretical predictions from state-of-the-art chemical evolution models for heavy elements, but a simple interpretation is not found either. In particular, the possible secondary behaviour of the [Rb/Fe] ratio at super-solar metallicities would require a much larger production of Rb than currently predicted in AGB stars through the s-process without overproducing Sr and Zr., Comment: 16 pag, 8 Figs, 1 table. Accepted in Astronomy & Astrophysics

Existing and upcoming instrumentation is collecting large amounts of astrophysical data, which require efficient and fast analysis techniques. We present a deep neural network architecture to analyze high-resolution stellar spectra and predict stellar parameters such as effective temperature, surface gravity, metallicity, and rotational velocity. With this study, we firstly demonstrate the capability of deep neural networks to precisely recover stellar parameters from a synthetic training set. Secondly, we analyze the application of this method to observed spectra and the impact of the synthetic gap (i.e., the difference between observed and synthetic spectra) on the estimation of stellar parameters, their errors, and their precision. Our convolutional network is trained on synthetic PHOENIX-ACES spectra in different optical and near-infrared wavelength regions. For each of the four stellar parameters, $T_{\rm eff}$, $\log{g}$, [M/H], and $v \sin{i}$, we constructed a neural network model to estimate each parameter independently. We then applied this method to 50 M dwarfs with high-resolution spectra taken with CARMENES (Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs), which operates in the visible (520-960 nm) and near-infrared wavelength range (960-1710 nm) simultaneously. Our results are compared with literature values for these stars. They show mostly good agreement within the errors, but also exhibit large deviations in some cases, especially for [M/H], pointing out the importance of a better understanding of the synthetic gap.

Context. Variability caused by stellar activity represents a challenge to the discovery and characterization of terrestrial exoplanets and complicates the interpretation of atmospheric planetary signals. Aims. We aim to use a detailed modeling tool to reproduce the effect of active regions on radial velocity measurements, which aids the identification of the key parameters that have an impact on the induced variability. Methods. We analyzed the effect of stellar activity on radial velocities as a function of wavelength by simulating the impact of the properties of spots, shifts induced by convective motions, and rotation. We focused our modeling effort on the active star YZ CMi (GJ 285), which was photometrically and spectroscopically monitored with CARMENES and the Telescopi Joan Or\'o. Results. We demonstrate that radial velocity curves at different wavelengths yield determinations of key properties of active regions, including spot filling factor, temperature contrast, and location, thus solving the degeneracy between them. Most notably, our model is also sensitive to convective motions. Results indicate a reduced convective shift for M dwarfs when compared to solar-type stars (in agreement with theoretical extrapolations) and points to a small global convective redshift instead of blueshift. Conclusions. Using a novel approach based on simultaneous chromatic radial velocities and light curves, we can set strong constraints on stellar activity, including an elusive parameter such as the net convective motion effect., Comment: A&A, in press

The He I infrared (IR) triplet at 10830 \AA is known as an activity indicator in solar-type stars and has become a primary diagnostic in exoplanetary transmission spectroscopy. He I lines are a tracer of the stellar extreme-ultraviolet irradiation from the transition region and corona. We study the variability of the He I IR triplet lines in a spectral time series of 319 M~dwarf stars that was obtained with the CARMENES high-resolution optical and near-infrared spectrograph at Calar Alto. We detect He I IR line variability in 18% of our sample stars, all of which show H$\alpha$ in emission. Therefore, we find detectable He I IR variability in 78% of the sub-sample of stars with H$\alpha$ emission. Detectable variability is strongly concentrated in the latest spectral sub-types, where the He I IR lines during quiescence are typically weak. The fraction of stars with detectable He I IR variation remains lower than 10% for stars earlier than M3.0 V, while it exceeds 30% for the later spectral sub-types. Flares are accompanied by particularly pronounced line variations, including strongly broadened lines with red and blue asymmetries. However, we also find evidence for enhanced He I IR absorption, which is potentially associated with increased high-energy irradiation levels at flare onset. Generally, He I IR and H$\alpha$ line variations tend to be correlated, with H$\alpha$ being the most sensitive indicator in terms of pseudo-equivalent width variation. This makes the He I IR triplet a favourable target for planetary transmission spectroscopy., Comment: 14 pages, 12 figures, accepted by A&A, full Table 2 only available electronically at CDS

The high-resolution, dual channel, visible and near-infrared spectrograph CARMENES offers exciting opportunities for stellar and exoplanetary research on M dwarfs. In this work we address the challenge of reaching the highest radial velocity precision possible with a complex, actively cooled, cryogenic instrument, such as the near-infrared channel. We describe the performance of the instrument and the work flow used to derive precise Doppler measurements from the spectra. The capability of both CARMENES channels to detect small exoplanets is demonstrated with the example of the nearby M5.0 V star CD Cet (GJ 1057), around which we announce a super-Earth ($4.0\pm0.4\,M_\oplus$) companion on a 2.29 d orbit.

The He I infrared (IR) line at a vacuum wavelength of 10833 A is a diagnostic for the investigation of atmospheres of stars and planets orbiting them. For the first time, we study the behavior of the He I IR line in a set of chromospheric models for M-dwarf stars, whose much denser chromospheres may favor collisions for the level population over photoionization and recombination, which are believed to be dominant in solar-type stars. For this purpose, we use published PHOENIX models for stars of spectral types M2 V and M3 V and also compute new series of models with different levels of activity following an ansatz developed for the case of the Sun. We perform a detailed analysis of the behavior of the He I IR line within these models. We evaluate the line in relation to other chromospheric lines and also the influence of the extreme ultraviolet (EUV) radiation field. The analysis of the He I IR line strengths as a function of the respective EUV radiation field strengths suggests that the mechanism of photoionization and recombination is necessary to form the line for inactive models, while collisions start to play a role in our most active models. Moreover, the published model set, which is optimized in the ranges of the Na I D2, H$\alpha$, and the bluest Ca II IR triplet line, gives an adequate prediction of the He I IR line for most stars of the stellar sample. Because especially the most inactive stars with weak He I IR lines are fit worst by our models, it seems that our assumption of a 100% filling factor of a single inactive component no longer holds for these stars., Comment: Accepted for publication in A&A, 13 pages, 9 figures, 4 tables

We report the confirmation of a transiting planet around the bright, inactive M0.5 V star TOI-1235 (TYC 4384-1735-1, V = 11.5 mag), whose transit signal was detected in the photometric time series of Sectors 14, 20, and 21 of the TESS space mission. We confirm the planetary nature of the transit signal, which has a period of 3.44 d, by using precise radial velocity measurements with CARMENES and HARPS-N spectrographs. A comparison of the properties derived for TOI-1235 b's with theoretical models reveals that the planet has a rocky composition, with a bulk density slightly higher than Earth's. In particular, we measure a mass of M_p = 5.9+/-0.6 M_Earth and a radius of R_p = 1.69+/-0.08 R_Earth, which together result in a density of rho_p = 6.7+1.3-1.1 g/cm3. When compared with other well-characterized exoplanetary systems, the particular combination of planetary radius and mass puts our discovery in the radius gap, a transition region between rocky planets and planets with significant atmospheric envelopes, with few known members. While the exact location of the radius gap for M dwarfs is still a matter of debate, our results constrain it to be located at around 1.7 R_Earth or larger at the insolation levels received by TOI-1235 b (~60 S_Earth), which makes it an extremely interesting object for further studies of planet formation and atmospheric evolution., Comment: A&A, in press

We report on radial velocity time series for two M0.0V stars, GJ338B and GJ338A, using the CARMENES spectrograph, complemented by ground-telescope photometry from Las Cumbres and Sierra Nevada observatories. We aim to explore the presence of small planets in tight orbits using the spectroscopic radial velocity technique. We obtained 159 and 70 radial velocity measurements of GJ338B and A, respectively, with the CARMENES visible channel. We also compiled additional relative radial velocity measurements from the literature and a collection of astrometric data that cover 200 a of observations to solve for the binary orbit. We found dynamical masses of 0.64$\pm$0.07M$_\odot$ for GJ338B and 0.69$\pm$0.07M$_\odot$ for GJ338A. The CARMENES radial velocity periodograms show significant peaks at 16.61$\pm$0.04 d (GJ338B) and 16.3$^{+3.5}_{-1.3}$ d (GJ338A), which have counterparts at the same frequencies in CARMENES activity indicators and photometric light curves. We attribute these to stellar rotation. GJ338B shows two additional, significant signals at 8.27$\pm$0.01 and 24.45$\pm$0.02 d, with no obvious counterparts in the stellar activity indices. The former is likely the first harmonic of the star's rotation, while we ascribe the latter to the existence of a super-Earth planet with a minimum mass of 10.27$^{+1.47}_{-1.38}$$M_{\oplus}$ orbiting GJ338B. GJ338B b lies inside the inner boundary of the habitable zone around its parent star. It is one of the least massive planets ever found around any member of stellar binaries. The masses, spectral types, brightnesses, and even the rotational periods are very similar for both stars, which are likely coeval and formed from the same molecular cloud, yet they differ in the architecture of their planetary systems., Comment: All the figures were size reduced for arXiv preprint archive. astro-ph.EP astro-ph.SR

For years, the standard procedure to measure radial velocities (RVs) of spectral observations consisted in cross-correlating the spectra with a binary mask, that is, a simple stellar template that contains information on the position and strength of stellar absorption lines. The cross-correlation function (CCF) profiles also provide several indicators of stellar activity. We present a methodology to first build weighted binary masks and, second, to compute the CCF of spectral observations with these masks from which we derive radial velocities and activity indicators. These methods are implemented in a python code that is publicly available. To build the masks, we selected a large number of sharp absorption lines based on the profile of the minima present in high signal-to-noise ratio (S/N) spectrum templates built from observations of reference stars. We computed the CCFs of observed spectra and derived RVs and the following three standard activity indicators: full-width-at-half-maximum as well as contrast and bisector inverse slope.We applied our methodology to CARMENES high-resolution spectra and obtain RV and activity indicator time series of more than 300 M dwarf stars observed for the main CARMENES survey. Compared with the standard CARMENES template matching pipeline, in general we obtain more precise RVs in the cases where the template used in the standard pipeline did not have enough S/N. We also show the behaviour of the three activity indicators for the active star YZ CMi and estimate the absolute RV of the M dwarfs analysed using the CCF RVs., Comment: 24 pages, 14 figures. Accepted for publication in A&A. Code available at https://github.com/mlafarga/raccoon

We present the discovery and characterisation of two transiting planets observed by the Transiting Exoplanet Survey Satellite (TESS) orbiting the nearby (d ~ 22 pc), bright (J ~ 9 mag) M3.5 dwarf LTT 3780 (TOI-732). We confirm both planets and their association with LTT 3780 via ground-based photometry and determine their masses using precise radial velocities measured with the CARMENES spectrograph. Precise stellar parameters determined from CARMENES high resolution spectra confirm that LTT 3780 is a mid-M dwarf with an effective temperature of T_eff = 3360 +\- 51 K, a surface gravity of log(g) = 4.81 +/- 0.04 (cgs), and an iron abundance of [Fe/H] = 0.09 +/- 0.16 dex, with an inferred mass of M_star = 0.379 +/- 0.016 M_sun and a radius of R_star = 0.382 +/- 0.012 R_sun. The ultra-short-period planet LTT 3780 b (P_b = 0.77 d) with a radius of 1.35^{+0.06}_{-0.06} R_earth, a mass of 2.34^{+0.24}_{-0.23} M_earth, and a bulk density of 5.24^{+0.94}_{-0.81} g cm^{-3} joins the population of Earth-size planets with rocky, terrestrial composition. The outer planet, LTT 3780 c, with an orbital period of 12.25 d, radius of 2.42^{+0.10}_{-0.10} R_earth, mass of 6.29^{+0.63}_{-0.61} M_earth, and mean density of 2.45^{+0.44}_{-0.37} g cm^{-3} belongs to the population of dense sub-Neptunes. With the two planets located on opposite sides of the radius gap, this planetary system is an excellent target for testing planetary formation, evolution and atmospheric models. In particular, LTT 3780 c is an ideal object for atmospheric studies with the James Webb Space Telescope., Comment: 21 pages, 15 figures, 4 tables, submitted to A&A

The nearby ultra-compact multiplanetary system YZ Ceti consists of at least three planets. The orbital period of each planet is the subject of discussion in the literature due to strong aliasing in the radial velocity data. The stellar activity of this M dwarf also hampers significantly the derivation of the planetary parameters. With an additional 229 radial velocity measurements obtained since the discovery publication, we reanalyze the YZ Ceti system and resolve the alias issues. We use model comparison in the framework of Bayesian statistics and periodogram simulations based on a method by Dawson and Fabrycky to resolve the aliases. We discuss additional signals in the RV data, and derive the planetary parameters by simultaneously modeling the stellar activity with a Gaussian process regression model. To constrain the planetary parameters further we apply a stability analysis on our ensemble of Keplerian fits. We resolve the aliases: the three planets orbit the star with periods of $2.02$ d, $3.06$ d, and $4.66$ d. We also investigate an effect of the stellar rotational signal on the derivation of the planetary parameters, in particular the eccentricity of the innermost planet. Using photometry we determine the stellar rotational period to be close to $68$ d. From the absence of a transit event with TESS, we derive an upper limit of the inclination of $i_\mathrm{max} = 87.43$ deg. YZ Ceti is a prime example of a system where strong aliasing hindered the determination of the orbital periods of exoplanets. Additionally, stellar activity influences the derivation of planetary parameters and modeling them correctly is important for the reliable estimation of the orbital parameters in this specific compact system. Stability considerations then allow additional constraints to be placed on the planetary parameters., Comment: Accepted for publication in A&A

Context. GJ 1148 is an M-dwarf star hosting a planetary system composed of two Saturn-mass planets in eccentric orbits with periods of 41.38 and 532.02 days. Aims. We reanalyze the orbital configuration and dynamics of the GJ 1148 multi-planetary system based on new precise radial velocity (RV) measurements taken with CARMENES. Methods. We combined new and archival precise Doppler measurements from CARMENES with those available from HIRES for GJ 1148 and modeled these data with a self-consistent dynamical model. We studied the orbital dynamics of the system using the secular theory and direct N-body integrations. The prospects of potentially habitable moons around GJ 1148 b were examined. Results. The refined dynamical analyses show that the GJ 1148 system is long-term stable in a large phase-space of orbital parameters with an orbital configuration suggesting apsidal alignment, but not in any particular high-order mean-motion resonant commensurability. GJ 1148 b orbits inside the optimistic habitable zone (HZ). We find only a narrow stability region around the planet where exomoons can exist. However, in this stable region exomoons exhibit quick orbital decay due to tidal interaction with the planet. Conclusions. The GJ 1148 planetary system is a very rare M-dwarf planetary system consisting of a pair of gas giants, the inner of which resides in the HZ. We conclude that habitable exomoons around GJ 1148 b are very unlikely to exist., Comment: Accepted for publication in A&A

The He I infrared (IR) triplet at 10830 AA is an important activity indicator for the Sun and in solar-type stars, however, it has rarely been studied in relation to M dwarfs to date. In this study, we use the time-averaged spectra of 319 single stars with spectral types ranging from M0.0 V to M9.0 V obtained with the CARMENES high resolution optical and near-infrared spectrograph at Calar Alto to study the properties of the He I IR triplet lines. In quiescence, we find the triplet in absorption with a decrease of the measured pseudo equivalent width (pEW) towards later sub-types. For stars later than M5.0 V, the He I triplet becomes undetectable in our study. This dependence on effective temperature may be related to a change in chromospheric conditions along the M dwarf sequence. When an emission in the triplet is observed, we attribute it to flaring. The absence of emission during quiescence is consistent with line formation by photo-ionisation and recombination, while flare emission may be caused by collisions within dense material. The He I triplet tends to increase in depth according to increasing activity levels, ultimately becoming filled in; however, we do not find a correlation between the pEW(He IR) and X-ray properties. This behaviour may be attributed to the absence of very inactive stars ($L_{\rm X}/L_{\rm bol}$ < -5.5) in our sample or to the complex behaviour with regard to increasing depth and filling in., Comment: 18 pages, 15 figures, A&A accepted

Stars show various amounts of radial velocity (RV) jitter due to varying stellar activity levels. The typical amount of RV jitter as a function of stellar age and observational timescale has not yet been systematically quantified, although it is often larger than the instrumental precision of modern high-resolution spectrographs used for Doppler planet detection and characterization. We aim to empirically determine the intrinsic stellar RV variation for mostly G and K dwarf stars on different timescales and for different stellar ages independently of stellar models. We also focus on young stars ($\lesssim$ 30 Myr), where the RV variation is known to be large. We use archival FEROS and HARPS RV data of stars which were observed at least 30 times spread over at least two years. We then apply the pooled variance (PV) technique to these data sets to identify the periods and amplitudes of underlying, quasiperiodic signals. We show that the PV is a powerful tool to identify quasiperiodic signals in highly irregularly sampled data sets. We derive activity-lag functions for 20 putative single stars, where lag is the timescale on which the stellar jitter is measured. Since the ages of all stars are known, we also use this to formulate an activity--age--lag relation which can be used to predict the expected RV jitter of a star given its age and the timescale to be probed. The maximum RV jitter on timescales of decades decreases from over 500 m/s for 5 Myr-old stars to 2.3 m/s for stars with ages of around 5 Gyr. The decrease in RV jitter when considering a timescale of only 1 d instead of 1 yr is smaller by roughly a factor of 4 for 5 Myr old stars, and a factor of 1.5 for stars with an age of 5 Gyr. The rate at which the RV jitter increases with lag strongly depends on stellar age and ranges from a few days for a few 10 Myr old stars to presumably decades for stars with an age of a few gigayears., Comment: 15 pages, 7 Figures; Changelog v2: Updated link to CDS for table E.1; rearranged Fig. 2 to match journal layout

Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an eccentric 204-day orbit around the very low-mass star GJ 3512. Dynamical models show that the high eccentricity of the orbit is most likely explained from planet-planet interactions. The reported planetary system challenges current formation theories and puts stringent constraints on the accretion and migration rates of planet formation and evolution models, indicating that disc instability may be more efficient in forming planets than previously thought., Comment: Manuscript author version. 41 pages, 11 figures