With the purpose of assessing classic spectroscopic methods on high-resolution and high signal-to-noise ratio spectra in the near-infrared wavelength region, we selected a sample of 65 F-, G-, and K-type stars observed with CARMENES, the new, ultra-stable, double-channel spectrograph at the 3.5 m Calar Alto telescope. We computed their stellar atmospheric parameters ($T_{\rm eff}$, $\log{g}$, $\xi$, and [Fe/H]) by means of the StePar code, a Python implementation of the equivalent width method that employs the 2017 version of the MOOG code and a grid of MARCS model atmospheres. We compiled four Fe I and Fe II line lists suited to metal-rich dwarfs, metal-poor dwarfs, metal-rich giants, and metal-poor giants that cover the wavelength range from 5300 to 17100 angstroms, thus substantially increasing the number of identified Fe I and Fe II lines up to 653 and 23, respectively. We examined the impact of the near-infrared Fe I and Fe II lines upon our parameter determinations after an exhaustive literature search, placing special emphasis on the 14 $Gaia$ benchmark stars contained in our sample. Even though our parameter determinations remain in good agreement with the literature values, the increase in the number of Fe I and Fe II lines when the near-infrared region is taken into account reveals a deeper $T_{\rm eff}$ scale that might stem from a higher sensitivity of the near-infrared lines to $T_{\rm eff}$.

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.

Context. Teegarden's Star is the brightest and one of the nearest ultra-cool dwarfs in the solar neighbourhood. For its late spectral type (M7.0V), the star shows relatively little activity and is a prime target for near-infrared radial velocity surveys such as CARMENES. Aims. As part of the CARMENES search for exoplanets around M dwarfs, we obtained more than 200 radial-velocity measurements of Teegarden's Star and analysed them for planetary signals. Methods. We find periodic variability in the radial velocities of Teegarden's Star. We also studied photometric measurements to rule out stellar brightness variations mimicking planetary signals. Results. We find evidence for two planet candidates, each with $1.1M_\oplus$ minimum mass, orbiting at periods of 4.91 and 11.4 d, respectively. No evidence for planetary transits could be found in archival and follow-up photometry. Small photometric variability is suggestive of slow rotation and old age. Conclusions. The two planets are among the lowest-mass planets discovered so far, and they are the first Earth-mass planets around an ultra-cool dwarf for which the masses have been determined using radial velocities.

Context: StePar is an automatic code written in Python 3.X designed to compute the stellar atmospheric parameters Teff, log(g), [Fe/H], and of FGK-type stars by means of the EW method. This code has already been extensively tested in different spectroscopic studies of FGK-type stars with several spectrographs and against myriads of Gaia-ESO Survey UVES U580 spectra of late-type, low-mass stars as one of its thirteen pipelines. Aims: We describe the code and test it against a library of well characterised Gaia benchmark stars. We also release the code to the community and provide the link for download. Methods: We carried out the required EW determination of Fe I and Fe II spectral lines using the automatic tool TAME. StePar implements a grid of MARCS model atmospheres and the MOOG radiative transfer code to compute stellar atmospheric parameters by means of a Downhill Simplex minimisation algorithm. Results: We show the results of the benchmark star test and also discuss the limitations of the EW method, and hence the code. In addition, we found a small internal scatter for the benchmark stars of 9 +- 32 K in Teff, 0.00 +- 0.07 dex in log(g), and 0.00 +- 0.03 dex in [Fe/H]. Finally, we advise against using StePar on double-lined spectroscopic binaries or spectra with R < 30,000, SNR < 20 or vsini >15 km/s as well as stars later than K4 or earlier than F6.

We determine the radii and masses of 293 nearby, bright M dwarfs of the CARMENES survey. This is the first time that such a large and homogeneous high-resolution (R>80 000) spectroscopic survey has been used to derive these fundamental stellar parameters. We derived the radii using Stefan-Boltzmann's law. We obtained the required effective temperatures $T_{\rm eff}$ from a spectral analysis and we obtained the required luminosities L from integrated broadband photometry together with the Gaia DR2 parallaxes. The mass was then determined using a mass-radius relation that we derived from eclipsing binaries known in the literature. We compared this method with three other methods: (1) We calculated the mass from the radius and the surface gravity log g, which was obtained from the same spectral analysis as $T_{\rm eff}$. (2) We used a widely used infrared mass-magnitude relation. (3) We used a Bayesian approach to infer stellar parameters from the comparison of the absolute magnitudes and colors of our targets with evolutionary models. Between spectral types M0V and M7V our radii cover the range $0.1\,R_{\normalsize\odot}<R<0.6\,R_{\normalsize\odot}$ with an error of 2-3% and our masses cover $0.09\,{\mathcal M}_{\normalsize\odot}<{\mathcal M}<0.6\,{\mathcal M}_{\normalsize\odot}$ with an error of 3-5%. We find good agreement between the masses determined with these different methods for most of our targets. Only the masses of very young objects show discrepancies. This can be well explained with the assumptions that we used for our methods.

We investigated almost 500 stars distributed among 193 binary or multiple systems made of late-F, G-, or early-K primaries and late-K or M dwarf companion candidates. For all of them, we compiled or measured coordinates, J-band magnitudes, spectral types, distances, and proper motions.With these data, we established a sample of 192 physically bound systems. In parallel, we carried out observations with HERMES/Mercator and obtained high resolution spectra for the 192 primaries and five secondaries. We used these spectra and the automatic StePar code for deriving precise stellar atmospheric parameters: Teff, log g, Vmicro, and chemical abundances for 13 atomic species, including [Fe/H]. After computing Galactocentric space velocities for all the primary stars, we performed a kinematic analysis and classified them in different Galactic populations and stellar kinematic groups of very different ages, which match our own metallicity determinations and isochronal age estimations. In particular, we identified three systems in the halo and 33 systems in the young Local Association, Ursa Major and Castor moving groups, and IC 2391 and Hyades Superclusters. We finally studied the exoplanet-metallicity relation in our 193 primaries and made a list 13 M-dwarf companions with very high metallicity that can be the targets of new dedicated exoplanet surveys. All in all, our dataset will be of great help for future works on the accurate determination of metallicity of M dwarfs.

The CARMENES radial velocity (RV) survey is observing 324 M dwarfs to search for any orbiting planets. In this paper, we present the survey sample by publishing one CARMENES spectrum for each M dwarf. These spectra cover the wavelength range 520--1710nm at a resolution of at least $R > 80,000$, and we measure its RV, H$\alpha$ emission, and projected rotation velocity. We present an atlas of high-resolution M-dwarf spectra and compare the spectra to atmospheric models. To quantify the RV precision that can be achieved in low-mass stars over the CARMENES wavelength range, we analyze our empirical information on the RV precision from more than 6500 observations. We compare our high-resolution M-dwarf spectra to atmospheric models where we determine the spectroscopic RV information content, $Q$, and signal-to-noise ratio. We find that for all M-type dwarfs, the highest RV precision can be reached in the wavelength range 700--900nm. Observations at longer wavelengths are equally precise only at the very latest spectral types (M8 and M9). We demonstrate that in this spectroscopic range, the large amount of absorption features compensates for the intrinsic faintness of an M7 star. To reach an RV precision of 1ms$^{-1}$ in very low mass M dwarfs at longer wavelengths likely requires the use of a 10m class telescope. For spectral types M6 and earlier, the combination of a red visual and a near-infrared spectrograph is ideal to search for low-mass planets and to distinguish between planets and stellar variability. At a 4m class telescope, an instrument like CARMENES has the potential to push the RV precision well below the typical jitter level of 3-4ms$^{-1}$.

Context: The main goal of the CARMENES survey is to find Earth-mass planets around nearby M-dwarf stars. Seven M-dwarfs included in the CARMENES sample had been observed before with HIRES and HARPS and either were reported to have one short period planetary companion (GJ15A, GJ176, GJ436, GJ536 and GJ1148) or are multiple planetary systems (GJ581 and GJ876). Aims: We aim to report new precise optical radial velocity measurements for these planet hosts and test the overall capabilities of CARMENES. Methods: We combined our CARMENES precise Doppler measurements with those available from HIRES and HARPS and derived new orbital parameters for the systems. Bona-fide single planet systems are fitted with a Keplerian model. The multiple planet systems were analyzed using a self-consistent dynamical model and their best fit orbits were tested for long-term stability. Results: We confirm or provide supportive arguments for planets around all the investigated stars except for GJ15A, for which we find that the post-discovery HIRES data and our CARMENES data do not show a signal at 11.4 days. Although we cannot confirm the super-Earth planet GJ15Ab, we show evidence for a possible long-period ($P_{\rm c}$ = 7025$_{-629}^{+972}$ d) Saturn-mass ($m_{\rm c} \sin i$ = 51.8$_{-5.8}^{+5.5}M_\oplus$) planet around GJ15A. In addition, based on our CARMENES and HIRES data we discover a second planet around GJ1148, for which we estimate a period $P_{\rm c}$ = 532.6$_{-2.5}^{+4.1}$ d, eccentricity $e_{\rm c}$ = 0.34$_{-0.06}^{+0.05}$ and minimum mass $m_{\rm c} \sin i$ = 68.1$_{-2.2}^{+4.9}M_\oplus$. Conclusions: The CARMENES optical radial velocities have similar precision and overall scatter when compared to the Doppler measurements conducted with HARPS and HIRES. We conclude that CARMENES is an instrument that is up to the challenge of discovering rocky planets around low-mass stars.

We report on the first star discovered to host a planet detected by radial velocity (RV) observations obtained within the CARMENES survey for exoplanets around M dwarfs. HD 147379 ($V = 8.9$ mag, $M = 0.58 \pm 0.08$ M$_{\odot}$), a bright M0.0V star at a distance of 10.7 pc, is found to undergo periodic RV variations with a semi-amplitude of $K = 5.1\pm0.4$ m s$^{-1}$ and a period of $P = 86.54\pm0.06$ d. The RV signal is found in our CARMENES data, which were taken between 2016 and 2017, and is supported by HIRES/Keck observations that were obtained since 2000. The RV variations are interpreted as resulting from a planet of minimum mass $m_{\rm p}\sin{i} = 25 \pm 2$ M$_{\oplus}$, 1.5 times the mass of Neptune, with an orbital semi-major axis $a = 0.32$ au and low eccentricity ($e < 0.13$). HD 147379b is orbiting inside the temperate zone around the star, where water could exist in liquid form. The RV time-series and various spectroscopic indicators show additional hints of variations at an approximate period of 21.1d (and its first harmonic), which we attribute to the rotation period of the star.