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.

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.

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.

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.

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.

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 LP 729-54 (TOI-732). We confirm both planets and their association with LP 729-54 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 LP 729-54 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 LP 729-54 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, LP 729-54 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, LP 729-54 c is an ideal object for atmospheric studies with the James Webb Space Telescope.

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.

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.

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.

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.

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.

Aims: We aim at detecting H$_2$O in the atmosphere of the hot Jupiter HD 209458 b and perform a multi-band study in the near infrared with CARMENES. Methods: The H$_2$O absorption lines from the planet's atmosphere are Doppler-shifted due to the large change in its radial velocity during transit. This shift is of the order of tens of km s$^{-1}$, whilst the Earth's telluric and the stellar lines can be considered quasi-static. We took advantage of this to remove the telluric and stellar lines using SYSREM, a principal component analysis algorithm. The residual spectra contain the signal from thousands of planetary molecular lines well below the noise level. We retrieve this information by cross-correlating the spectra with models of the atmospheric absorption. Results: We find evidence of H$_2$O in HD 209458 b with a signal-to-noise ratio (S/N) of 6.4. The signal is blueshifted by --5.2 $^{+2.6}_{-1.3}$ km s$^{-1}$, which, despite the error bars, is a firm indication of day-to-night winds at the terminator of this hot Jupiter. Additionally, we performed a multi-band study for the detection of H$_2$O individually from the three NIR bands covered by CARMENES. We detect H$_2$O from its 1.0 $\mu$m band with a S/N of 5.8, and also find hints from the 1.15 $\mu$m band, with a low S/N of 2.8. No clear planetary signal is found from the 1.4 $\mu$m band. Conclusions: Our significant signal from the 1.0 $\mu$m band in HD 209458 b represents the first detection of H$_2$O from this band, the bluest one to date. The unfavorable observational conditions might be the reason for the inconclusive detection from the stronger 1.15 and 1.4 $\mu$m bands. H$_2$O is detected from the 1.0 $\mu$m band in HD 209458 b, but hardly in HD 189733 b, which supports a stronger aerosol extinction in the latter.

Context: Recently, the He I triplet at 10830 \r{A} has been rediscovered as an excellent probe of the extended and possibly evaporating atmospheres of close-in transiting planets. This has already resulted in detections of this triplet in the atmospheres of a handful of planets, both from space and from the ground. However, while a strong signal is expected for the hot Jupiter HD 209458 b, only upper limits have been obtained so far. Aims: Our goal is to measure the helium excess absorption from HD 209458 b and assess the extended atmosphere of the planet and possible evaporation. Methods: We obtained new high-resolution spectral transit time-series of HD 209458 b using CARMENES at the 3.5 m Calar Alto telescope, targeting the He I triplet at 10830 \r{A} at a spectral resolving power of 80 400. The observed spectra were corrected for stellar absorption lines using out of transit data, for telluric absorption using the molecfit software, and for the sky emission lines using simultaneous sky measurements through a second fibre. Results: We detect He I absorption at a level of 0.91 $\pm$ 0.10 % (9 $\sigma$) at mid-transit. The absorption follows the radial velocity change of the planet during transit, unambiguously identifying the planet as the source of the absorption. The core of the absorption exhibits a net blueshift of 1.8 $\pm$ 1.3 km s$^{-1}$. Possible low-level excess absorption is seen further blueward from the main absorption near the centre of the transit, which could be caused by an extended tail. However, this needs to be confirmed. Conclusions: Our results further support a close relationship between the strength of planetary absorption in the helium triplet lines and the level of ionising, stellar X-ray and extreme-UV irradiation.

We report the detection of a transiting Earth-size planet around GJ 357, a nearby M2.5V star, using data from the Transiting Exoplanet Survey Satellite (TESS). GJ 357 b (TOI-562.01) is a transiting, hot, Earth-sized planet (Teq=525+-11 K) with a radius of Rb=1.217+-0.084 Re and an orbital period of Pb=3.93 d. Precise stellar radial velocities from CARMENES and PFS, as well as archival data from HIRES, UVES, and HARPS also display a 3.93-day periodicity, confirming the planetary nature and leading to a planetary mass of Mb=1.84+-0.31 Me. In addition to the radial velocity signal for GJ 357 b, more periodicities are present in the data indicating the presence of two further planets in the system: GJ 357 c, with a minimum mass of Mc=3.40+-0.46 Me in a 9.12 d orbit, and GJ 357 d, with a minimum mass of Md=6.1+-1.0 Me in a 55.7 d orbit inside the habitable zone. The host is relatively inactive and exhibits a photometric rotation period of Prot=78+-2 d. GJ 357 b is to date the second closest transiting planet to the Sun, making it a prime target for further investigations such as transmission spectroscopy. Therefore, GJ 357 b represents one of the best terrestrial planets suitable for atmospheric characterization with the upcoming JWST and ground-based ELTs.

We present precise photospheric parameters of 282 M dwarfs determined from fitting the most recent version of PHOENIX models to high-resolution CARMENES spectra in the visible (0.52 - 0.96 $\mu$m) and near-infrared wavelength range (0.96 - 1.71 $\mu$m). With its aim to search for habitable planets around M dwarfs, several planets of different masses have been detected. The characterization of the target sample is important for the ability to derive and constrain the physical properties of any planetary systems that are detected. As a continuation of previous work in this context, we derived the fundamental stellar parameters effective temperature, surface gravity, and metallicity of the CARMENES M-dwarf targets from PHOENIX model fits using a $\chi^2$ method. We calculated updated PHOENIX stellar atmosphere models that include a new equation of state to especially account for spectral features of low-temperature stellar atmospheres as well as new atomic and molecular line lists. We show the importance of selecting magnetically insensitive lines for fitting to avoid effects of stellar activity in the line profiles. For the first time, we directly compare stellar parameters derived from multiwavelength range spectra, simultaneously observed for the same star. In comparison with literature values we show that fundamental parameters derived from visible spectra and visible and near-infrared spectra combined are in better agreement than those derived from the same spectra in the near-infrared alone.

We present the discovery of TYC9191-519-1b (TOI-150b, TIC 271893367) and HD271181b (TOI-163b, TIC 179317684), two hot Jupiters initially detected using 30-minute cadence Transiting Exoplanet Survey Satellite TESS photometry from Sector 1 and thoroughly characterized through follow-up photometry (CHAT, Hazelwood, LCO/CTIO, El Sauce, TRAPPIST-S), high-resolution spectroscopy (FEROS, CORALIE) and speckle imaging (Gemini/DSSI), confirming the planetary nature of the two signals. A simultaneous joint fit of photometry and radial velocity using a new fitting package juliet reveals that TOI-150b is a $1.254\pm0.016\ R_J$, massive ($2.61^{+0.19}_{-0.12}\ M_J$) hot Jupiter in a $5.857$-day orbit, while TOI-163b is an inflated ($R_P$ = $1.478^{+0.022}_{-0.029} R_J$, $M_P$ = $1.219\pm0.11 M_J$) hot Jupiter on a $P$ = $4.231$-day orbit; both planets orbit F-type stars. A particularly interesting result is that TOI-150b shows an eccentric orbit ($e=0.262^{+0.045}_{-0.037}$), which is quite uncommon among hot Jupiters. We estimate that this is consistent, however, with the circularization timescale which is slightly larger than the age of the system. These two hot Jupiters are both prime candidates for further characterization --- in particular, both are excellent candidates for determining spin-orbit alignments via the Rossiter-McLaughlin (RM) effect and for characterizing atmospheric thermal structures using secondary eclipse observations considering they are both located closely to the James Webb Space Telescope (JWST) Continuous Viewing Zone (CVZ).

Although M dwarfs are known for high levels of stellar activity, they are ideal targets for the search of low-mass exoplanets with the radial velocity (RV) method. We report the discovery of a planetary-mass companion around LSPM J2116+0234 (M3.0 V) and confirm the existence of a planet orbiting GJ 686 (BD+18 3421; M1.0 V). The discovery of the planet around LSPM J2116+0234 is based on CARMENES RV observations in the visual and near-infrared channels. We confirm the planet orbiting around GJ 686 by analyzing the RV data spanning over two decades of observations from CARMENES VIS, HARPS-N, HARPS, and HIRES. We find planetary signals at 14.44 and 15.53 d in the RV data for LSPM J2116+0234 and GJ 686, respectively. Additionally, the RV, photometric time series, and various spectroscopic indicators show hints of variations of 42 d for LSPM J2116+0234 and 37 d for GJ 686, which we attribute to the stellar rotation periods. The orbital parameters of the planets are modeled with Keplerian fits together with correlated noise from the stellar activity. A mini-Neptune with a minimum mass of 11.8 Me orbits LSPM J2116+0234 producing an RV semi-amplitude of 6.19 m/s, while a super-Earth of mass 6.6 Me orbits GJ 686 and produces an RV semi-amplitude of 3.0 m/s. Both LSPM J2116+0234 and GJ 686 have planetary companions populating the regime of exoplanets with masses lower than 15 Me and orbital periods <20 d.

M dwarfs are known to generate the strongest magnetic fields among main-sequence stars with convective envelopes, but the link between the magnetic fields and underlying dynamo mechanisms, rotation, and activity still lacks a consistent picture. In this work we measure magnetic fields from the high-resolution near-infrared spectra taken with the CARMENES radial-velocity planet survey in a sample of 29 active M dwarfs and compare our results against stellar parameters. We use the state-of-the-art radiative transfer code to measure total magnetic flux densities from the Zeeman broadening of spectral lines and filling factors. We detect strong kG magnetic fields in all our targets. In 16 stars the magnetic fields were measured for the first time. Our measurements are consistent with the magnetic field saturation in stars with rotation periods P<4d. The analysis of the magnetic filling factors reveal two different patterns of either very smooth distribution or a more patchy one, which can be connected to the dynamo state of the stars and/or stellar mass. Our measurements extend the list of M dwarfs with strong surface magnetic fields. They also allow us to better constrain the interplay between the magnetic energy, stellar rotation, and underlying dynamo action. The high spectral resolution and observations at near-infrared wavelengths are the beneficial capabilities of the CARMENES instrument that allow us to address important questions about the stellar magnetism.

Small planets around low-mass stars often show orbital periods in a range that corresponds to the temperate zones of their host stars which are therefore of prime interest for planet searches. Surface phenomena such as spots and faculae create periodic signals in radial velocities and in observational activity tracers in the same range, so they can mimic or hide true planetary signals. We aim to detect Doppler signals corresponding to planetary companions, determine their most probable orbital configurations, and understand the stellar activity and its impact on different datasets. We analyze 22 years of data of the M1.5V-type star Gl49 (BD+61 195) including HARPS-N and CARMENES spectrographs, complemented by APT2 and SNO photometry. Activity indices are calculated from the observed spectra, and all datasets are analyzed with periodograms and noise models. We investigate how the variation of stellar activity imprints on our datasets. We further test the origin of the signals and investigate phase shifts between the different sets. To search for the best-fit model we maximize the likelihood function in a Markov Chain Monte Carlo approach. As a result of this study, we are able to detect the super-Earth Gl49b with a minimum mass of 5.6 Ms. It orbits its host star with a period of 13.85d at a semi-major axis of 0.090 au and we calculate an equilibrium temperature of 350 K and a transit probability of 2.0%. The contribution from the spot-dominated host star to the different datasets is complex, and includes signals from the stellar rotation at 18.86d, evolutionary time-scales of activity phenomena at 40-80d, and a long-term variation of at least four years.

We aim at the discovery of new planetary systems by exploiting the transit light curve results from TESS orbital observatory's Sector 1 and 2 observations and validating them with precise Doppler measurements obtained from archival HARPS data. Taking advantage of the reported TESS transit events around GJ143 (TOI 186) and HD23472 (TOI 174) we model their HARPS precise Doppler measurements and derive orbital parameters for these two systems. For the GJ143 system TESS has reported only a single transit, and thus its period is unconstrained from photometry. Our RV analysis of GJ143 reveal the full Keplerian solution of the system, which is consistent with an eccentric planet with a mass almost twice that of Neptune and a period of $P_{\rm b}$ = $35.59_{-0.01}^{+0.01}$ days. Our estimates of the GJ143 b planet are fully consistent with the transit timing from TESS. We confirm the two-planet system around HD23472, which according to our analysis is composed of two Neptune mass planets in a possible 5:3 MMR.

Chromospheric modeling of observed differences in stellar activity lines is imperative to fully understand the upper atmospheres of late-type stars. We present one-dimensional parametrized chromosphere models computed with the atmosphere code PHOENIX using an underlying photosphere of 3500 K. The aim of this work is to model chromospheric lines of a sample of 50 M2-3 dwarfs observed in the framework of the CARMENES, the Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs, exoplanet survey. The spectral comparison between observed data and models is performed in the chromospheric lines of Na I D2, H$\alpha$, and the bluest Ca II infrared triplet line to obtain best-fit models for each star in the sample. We find that for inactive stars a single model with a VAL C-like temperature structure is sufficient to describe simultaneously all three lines adequately. Active stars are rather modeled by a combination of an inactive and an active model, also giving the filling factors of inactive and active regions. Moreover, the fitting of linear combinations on variable stars yields relationships between filling factors and activity states, indicating that more active phases are coupled to a larger portion of active regions on the surface of the star.

The CARMENES survey is searching for Earth-like planets orbiting M dwarfs using the radial velocity method. Studying the stellar activity of the target stars is important to avoid false planet detections and to improve our understanding of the atmospheres of late-type stars. In this work we present measurements of activity indicators at visible and near-infrared wavelengths for 331 M dwarfs observed with CARMENES. Our aim is to identify the activity indicators that are most sensitive and easiest to measure, and the correlations among these indicators. We also wish to characterise their variability. Using a spectral subtraction technique, we measured pseudo-equivalent widths of the He I D3, H$\alpha$, He I $\lambda$10833 {\AA}, and Pa$\beta$ lines, the Na I D doublet, and the Ca II infrared triplet, which have a chromospheric component in active M dwarfs. In addition, we measured an index of the strength of two TiO and two VO bands, which are formed in the photosphere. We also searched for periodicities in these activity indicators for all sample stars using generalised Lomb-Scargle periodograms. We find that the most slowly rotating stars of each spectral subtype have the strongest H$\alpha$ absorption. H$\alpha$ is correlated most strongly with He I D3, whereas Na I D and the Ca II infrared triplet are also correlated with H$\alpha$. He I $\lambda$10833 {\AA} and Pa$\beta$ show no clear correlations with the other indicators. The TiO bands show an activity effect that does not appear in the VO bands. We find that the relative variations of H$\alpha$ and He I D3 are smaller for stars with higher activity levels, while this anti-correlation is weaker for Na I D and the Ca II infrared triplet, and is absent for He I $\lambda$10833 {\AA} and Pa$\beta$. Periodic variation with the rotation period most commonly appears in the TiO bands, H$\alpha$, and in the Ca II infrared triplet.

We use spectra from CARMENES, the Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs, to search for periods in chromospheric indices in 16 M0 to M2 dwarfs. We measure spectral indices in the H$\alpha$, the Ca II infrared triplet (IRT), and the Na I D lines to study which of these indices are best-suited to find rotation periods in these stars. Moreover, we test a number of different period-search algorithms, namely the string length method, the phase dispersion minimisation, the generalized Lomb-Scargle periodogram, and the Gaussian process regression with quasi-periodic kernel. We find periods in four stars using H$\alpha$ and in five stars using the Ca II IRT, two of which have not been found before. Our results show that both H$\alpha$ and the Ca II IRT lines are well suited for period searches, with the Ca II IRT index performing slightly better than H$\alpha$. Unfortunately, the Na I D lines are strongly affected by telluric airglow, and we could not find any rotation period using this index. Further, different definitions of the line indices have no major impact on the results. Comparing the different search methods, the string length method and the phase dispersion minimisation perform worst, while Gaussian process models produce the smallest numbers of false positives and non-detections.

We report the detection of a Neptune-mass exoplanet around the M4.0 dwarf GJ 4276 (G 232-070) based on radial velocity (RV) observations obtained with the CARMENES spectrograph. The RV variations of GJ 4276 are best explained by the presence of a planetary companion that has a minimum mass of $m_{\rm b}\sin i \approx 16\, M_\oplus$ on a $P_{\rm b}=13.35$ day orbit. The analysis of the activity indicators and spectral diagnostics exclude stellar induced RV perturbations and prove the planetary interpretation of the RV signal. We show that a circular single-planet solution can be excluded by means of a likelihood ratio test. Instead, we find that the RV variations can be explained either by an eccentric orbit or interpreted as a pair of planets on circular orbits near a period ratio of 2:1. Although the eccentric single-planet solution is slightly preferred, our statistical analysis indicates that none of these two scenarios can be rejected with high confidence using the RV time series obtained so far. Based on the eccentric interpretation, we find that GJ 4276 b is the most eccentric ($e_{\rm b} = 0.37$) exoplanet around an M dwarf with such a short orbital period known today.

We present three transit observations of HD 189733 b obtained with the high-resolution spectrograph CARMENES at Calar Alto. A strong absorption signal is detected in the near-infrared He I triplet at 10830 \AA{} in all three transits. During mid-transit, the mean absorption level is $0.88\pm0.04$ % measured in a $\pm$10 km s$^{-1}$ range at a net blueshift of $-3.5\pm0.4$ km s$^{-1}$ (10829.84--10830.57 \AA{}). The absorption signal exhibits radial velocities of $+6.5\pm3.1$ km s$^{-1}$ and $-12.6\pm1.0$ km s$^{-1}$ during ingress and egress, respectively; measured in the planetary rest frame. We show that stellar activity related pseudo-signals interfere with the planetary atmospheric absorption signal. They could contribute as much as 80% of the observed signal and might also affect the radial velocity signature, but pseudo-signals are very unlikely to explain the entire signal. The observed line ratio between the two unresolved and the third line of the He I triplet is $2.8\pm0.2$, which strongly deviates from the value expected for an optically thin atmospheres. When interpreted in terms of absorption in the planetary atmosphere, this favors a compact helium atmosphere with an extent of only 0.2 planetary radii and a substantial column density on the order of $4\times 10^{12}$ cm$^{-2}$. The observed radial velocities can be understood either in terms of atmospheric circulation with equatorial superrotation or as a sign of an asymmetric atmospheric component of evaporating material. We detect no clear signature of ongoing evaporation, like pre- or post-transit absorption, which could indicate material beyond the planetary Roche lobe, or radial velocities in excess of the escape velocity. These findings do not contradict planetary evaporation, but only show that the detected helium absorption in HD 189733 b does not trace the atmospheric layers that show pronounced escape signatures.

At a distance of 1.8 parsecs, Barnard's star (Gl 699) is a red dwarf with the largest apparent motion of any known stellar object. It is the closest single star to the Sun, second only to the alpha Centauri triple stellar system. Barnard's star is also among the least magnetically active red dwarfs known and has an estimated age older than our Solar System. Its properties have made it a prime target for planet searches employing techniques such as radial velocity, astrometry, and direct imaging, all with different sensitivity limits but ultimately leading to disproved or null results. Here we report that the combination of numerous measurements from high-precision radial velocity instruments reveals the presence of a low-amplitude but significant periodic signal at 233 days. Independent photometric and spectroscopic monitoring, as well as the analysis of instrumental systematic effects, show that this signal is best explained as arising from a planetary companion. The candidate planet around Barnard's star is a cold super-Earth with a minimum mass of 3.2 Earth masses orbiting near its snow-line. The combination of all radial velocity datasets spanning 20 years additionally reveals a long-term modulation that could arise from a magnetic activity cycle or from a more distant planetary object. Because of its proximity to the Sun, the proposed planet has a maximum angular separation of 220 milli-arcseconds from Barnard's star, making it an excellent target for complementary direct imaging and astrometric observations.

Aims: We explore the capabilities of CARMENES for characterizing hot-Jupiter atmospheres by targeting multiple water bands, in particular, those at 1.15 and 1.4 $\mu$m. Hubble Space Telescope observations suggest that this wavelength region is relevant for distinguishing between hazy/cloudy and clear atmospheres. Methods: We observed one transit of the hot Jupiter HD 189733 b with CARMENES. Telluric and stellar absorption lines were removed using Sysrem, which performs a principal component analysis including proper error propagation. The residual spectra were analysed for water absorption with cross-correlation techniques using synthetic atmospheric absorption models. Results: We report a cross-correlation peak at a signal-to-noise ratio (SNR) of 6.6, revealing the presence of water in the transmission spectrum of HD 189733 b. The absorption signal appeared slightly blueshifted at -3.9 $\pm$ 1.3 kms$^{-1}$. We measured the individual cross-correlation signals of the water bands at 1.15 and 1.4 $\mu$m, finding cross-correlation peaks at SNRs of 4.9 and 4.4, respectively. The 1.4 $\mu$m feature is consistent with that observed with the Hubble Space Telescope. Conclusions: The water bands studied in this work have been mainly observed in a handful of planets from space. The ability of also detecting them individually from the ground at higher spectral resolution can provide insightful information to constrain the properties of exoplanet atmospheres. Although the current multiband detections can not yet constrain atmospheric haze models for HD 189733 b, future observations at higher signal-to-noise ratio could provide an alternative way to achieve this aim.

We announce the discovery of two planetary companions orbiting around the low mass stars Ross 1020 (GJ 3779, M4.0V) and LP 819-052 (GJ 1265, M4.5V). The discovery is based on the analysis of CARMENES radial velocity observations in the visual channel as part of its survey for exoplanets around M dwarfs. In the case of GJ 1265, CARMENES observations were complemented with publicly available Doppler measurements from HARPS. The datasets reveal one planetary companion for each star that share very similar properties: minimum masses of $8.0\pm0.5$ M$_{\oplus}$ and $7.4\pm0.5$ M$_{\oplus}$ in low-eccentricity orbits with periods of $3.023\pm0.001$ d and $3.651\pm0.001$ d for GJ 3779 b and GJ 1265 b, respectively. The periodic signals around three days found in the radial velocity data have no counterpart in any spectral activity indicator. Besides, we collected available photometric data for the two host stars, which confirm that the additional Doppler variations found at periods around 95 d can be attributed to the rotation of the stars. The addition of these planets in a mass-period diagram of known planets around M dwarfs suggests a bimodal distribution with a lack of short-period low-mass planets in the range of 2-5 M$_{\oplus}$. It also indicates that super-Earths (> 5 M$_{\oplus}$) currently detected by radial velocity and transit techniques around M stars are usually found in systems dominated by a single planet.

In this paper we present new precise Doppler data of GJ1046 taken between November 2005 and July 2018 with the HARPS and the FEROS high-resolution spectographs. In addition, we provide a new stellar mass estimate of GJ1046 and we update the orbital parameters of the GJ1046 system. These new data and analysis could be used together with the GAIA epoch astrometry, when available, for braking the $\sin i$ degeneracy and revealing the true mass of the GJ1046 system.

Context. The CARMENES spectrograph is surveying ~300 M dwarf stars in search for exoplanets. Among the target stars, spectroscopic binary systems have been discovered, which can be used to measure fundamental properties of stars. Aims. Using spectroscopic observations, we determine the orbital and physical properties of nine new double-line spectroscopic binary systems by analysing their radial velocity curves. Methods. We use two-dimensional cross-correlation techniques to derive the radial velocities of the targets, which are then employed to determine the orbital properties. Photometric data from the literature are also analysed to search for possible eclipses and to measure stellar variability, which can yield rotation periods. Results. Out of the 342 stars selected for the CARMENES survey, 9 have been found to be double-line spectroscopic binaries, with periods ranging from 1.13 to ~8000 days and orbits with eccentricities up to 0.54. We provide empirical orbital properties and minimum masses for the sample of spectroscopic binaries. Absolute masses are also estimated from mass-luminosity calibrations, ranging between ~0.1 and ~0.6 Msol . Conclusions. These new binary systems increase the number of double-line M dwarf binary systems with known orbital parameters by 15%, and they have lower mass ratios on average.

Despite their activity, low-mass stars are of particular importance for the search of exoplanets by the means of Doppler spectroscopy, as planets with lower masses become detectable. We report on the discovery of a planetary companion around HD 180617, a bright J = 5.58 mag, low-mass M = 0.45 M_{sun} star of spectral type M2.5 V. The star, located at a distance of 5.9 pc, is the primary of the high proper motion binary system containing vB 10, a star with one of the lowest masses known in most of the twentieth century. Our analysis is based on new radial velocity (RV) measurements made at red-optical wavelengths provided by the high-precision spectrograph CARMENES, which was designed to carry out a survey for Earth-like planets around M dwarfs. The available CARMENES data are augmented by archival Doppler measurements from HIRES and HARPS. Altogether, the RVs span more than 16 years. The modeling of the RV variations, with a semi-amplitude of K = 2.85-0.25/+0.16m/s yields a Neptune-like planet with a minimum mass of 12.2-1.4/+1.0 M_{Earth} on a 105.90-0.10/+0.09d circumprimary orbit, which is partly located in the host star's habitable zone. The analysis of time series of common activity indicators does not show any dependence on the detected RV signal. The discovery of HD 180617 b not only adds information to a currently hardly filled region of the mass-period diagram of exoplanets around M dwarfs, but the investigated system becomes the third known binary consisting of M dwarfs and hosting an exoplanet in an S-type configuration. Its proximity makes it an attractive candidate for future studies.

The Kepler extended mission, also known as K2, has provided the community with a wealth of planetary candidates that orbit stars typically much brighter than the targets of the original mission. These planet candidates are suitable for further spectroscopic follow-up and precise mass determinations, leading ultimately to the construction of empirical mass-radius diagrams. Particularly interesting is to constrain the properties of planets between the Earth and Neptune in size, the most abundant type of planets orbiting Sun-like stars with periods less than a few years. Among many other K2 candidates, we discovered a multi-planetary system around EPIC246471491, with four planets ranging in size from twice the size of Earth, to nearly the size of Neptune. We measure the mass of the planets of the EPIC246471491 system by means of precise radial velocity measurements using the CARMENES spectrograph and the HARPS-N spectrograph. With our data we are able to determine the mass of the two inner planets of the system with a precision better than 15%, and place upper limits on the masses of the two outer planets. We find that EPIC246471491b has a mass of 9.68 Me, and a radius of 2.59 Re, yielding a mean density of 3.07 g/cm3, while EPIC246471491c has a mass of 15.68 Me, radius of 3.53 Re, and a mean density of 19.5 g/cm3. For EPIC246471491d (R=2.48Re) and EPIC246471491e (R=1.95Re) the upper limits for the masses are 6.5 and 10.7 Me, respectively. The system is thus composed of a nearly Neptune-twin planet (in mass and radius), two sub-Neptunes with very different densities and presumably bulk composition, and a fourth planet in the outermost orbit that resides right in the middle of the super-Earth/sub-Neptune radius gap. Future comparative planetology studies of this system can provide useful insights into planetary formation, and also a good test of atmospheric escape and evolution theories.

K2-18 is a nearby M2.5 dwarf, located at 34 pc and hosting a transiting planet which was first discovered by the {\it K2} mission and later confirmed with {\it Spitzer Space Telescope} observations. With a radius of $\sim 2 \, R_{\oplus}$ and an orbital period of $\sim 33$ days, the planet lies in the temperate zone of its host star and receives stellar irradiation similar to Earth. Here we perform radial velocity follow-up observations with the visual channel of CARMENES with the goal of determining the mass and density of the planet. We measure a planetary semi-amplitude of $K_b \sim 3.5$ \mpersec\ and a mass of $M_b \sim 9 \, M_{\oplus}$, yielding a bulk density around $\rho_b \sim 4 \, \mathrm{g \,cm^{-3}}$. This indicates a low-mass planet with a composition consistent with a solid core and a volatile-rich envelope. A signal at 9 days was recently reported using radial velocity measurements taken with the HARPS spectrograph. This was interpreted as being due to a second planet. We see a weaker, time and wavelength dependent signal in the CARMENES data set and thus favor stellar activity for its origin. \ktwo\ joins the growing group of low-mass planets detected in the temperate zone of M dwarfs. The brightness of the host star in the near-infrared makes the system a good target for detailed atmospheric studies with the {\it James Webb Space Telescope}.

Previous simulations predicted the activity-induced radial-velocity (RV) variations of M dwarfs to range from $\sim1$ cm/s to $\sim1$ km/s, depending on various stellar and activity parameters. We investigate the observed relations between RVs, stellar activity, and stellar parameters of M dwarfs by analyzing CARMENES high-resolution visual-channel spectra ($0.5$$-$$1$$\mu$m), which were taken within the CARMENES RV planet survey during its first $20$ months of operation. During this time, $287$ of the CARMENES-sample stars were observed at least five times. From each spectrum we derived a relative RV and a measure of chromospheric H$\alpha$ emission. In addition, we estimated the chromatic index (CRX) of each spectrum, which is a measure of the RV wavelength dependence. Despite having a median number of only $11$ measurements per star, we show that the RV variations of the stars with RV scatter of $>10$ m/s and a projected rotation velocity $v \sin{i}>2$ km/s are caused mainly by activity. We name these stars `active RV-loud stars' and find their occurrence to increase with spectral type: from $\sim3\%$ for early-type M dwarfs (M$0.0$$-$$2.5$V) through $\sim30\%$ for mid-type M dwarfs (M$3.0$$-$$5.5$V) to $>50\%$ for late-type M dwarfs (M$6.0$$-$$9.0$V). Their RV-scatter amplitude is found to be correlated mainly with $v \sin{i}$. For about half of the stars, we also find a linear RV$-$CRX anticorrelation, which indicates that their activity-induced RV scatter is lower at longer wavelengths. For most of them we can exclude a linear correlation between RV and H$\alpha$ emission. Our results are in agreement with simulated activity-induced RV variations in M dwarfs. The RV variations of most active RV-loud M dwarfs are likely to be caused by dark spots on their surfaces, which move in and out of view as the stars rotate.

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}$.

Stellar activity is ubiquitously encountered in M dwarfs and often characterised by the H$\alpha$ line. In the most active M dwarfs, H$\alpha$ is found in emission, sometimes with a complex line profile. Previous studies have reported extended wings and asymmetries in the H$\alpha$ line during flares. We used a total of 473 high-resolution spectra of 28 active M dwarfs obtained by the CARMENES (Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs) spectrograph to study the occurrence of broadened and asymmetric H$\alpha$ line profiles and their association with flares, and examine possible physical explanations. We detected a total of 41 flares and 67 broad, potentially asymmetric, wings in H$\alpha$. The broadened H$\alpha$ lines display a variety of profiles with symmetric cases and both red and blue asymmetries. Although some of these line profiles are found during flares, the majority are at least not obviously associated with flaring. We propose a mechanism similar to coronal rain or chromospheric downward condensations as a cause for the observed red asymmetries; the symmetric cases may also be caused by Stark broadening. We suggest that blue asymmetries are associated with rising material, and our results are consistent with a prevalence of blue asymmetries during the flare onset. Besides the H$\alpha$ asymmetries, we find some cases of additional line asymmetries in \ion{He}{i} D$_{3}$, \ion{Na}{i}~D lines, and the \ion{He}{i} line at 10830\,\AA\, taken all simultaneously thanks to the large wavelength coverage of CARMENES. Our study shows that asymmetric H$\alpha$ lines are a rather common phenomenon in M~dwarfs and need to be studied in more detail to obtain a better understanding of the atmospheric dynamics in these objects.

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.

Context: The CARMENES survey is a high-precision radial velocity (RV) programme that aims to detect Earth-like planets orbiting low-mass stars. Aims: We develop least-squares fitting algorithms to derive the RVs and additional spectral diagnostics implemented in the SpEctrum Radial Velocity Analyser (SERVAL), a publicly available python code. Methods: We measured the RVs using high signal-to-noise templates created by coadding all available spectra of each star.We define the chromatic index as the RV gradient as a function of wavelength with the RVs measured in the echelle orders. Additionally, we computed the differential line width by correlating the fit residuals with the second derivative of the template to track variations in the stellar line width. Results: Using HARPS data, our SERVAL code achieves a RV precision at the level of 1m/s. Applying the chromatic index to CARMENES data of the active star YZ CMi, we identify apparent RV variations induced by stellar activity. The differential line width is found to be an alternative indicator to the commonly used full width half maximum. Conclusions: We find that at the red optical wavelengths (700--900 nm) obtained by the visual channel of CARMENES, the chromatic index is an excellent tool to investigate stellar active regions and to identify and perhaps even correct for activity-induced RV variations.

Aims. Our new program with HARPS aims to detect mean motion resonant planetary systems around stars which were previously reported to have a single bona fide planet, often based only on sparse radial velocity data. Methods. Archival and new HARPS radial velocities for the K2V star HD 27894 were combined and fitted with a three-planet self-consistent dynamical model. The best-fit orbit was tested for long-term stability. Results. We find clear evidence that HD 27894 is hosting at least three massive planets. In addition to the already known Jovian planet with a period $P_{\rm b}$ $\approx$ 18 days we discover a Saturn-mass planet with $P_{\rm c}$ $\approx$ 36 days, likely in a 2:1 mean motion resonance with the first planet, and a cold massive planet ($\approx$ 5.3 $M_{\mathrm{Jup}}$) with a period $P_{\rm d}$ $\approx$ 5170 days on a moderately eccentric orbit ($e_{\rm d}$ = 0.39). Conclusions. HD 27894 is hosting a massive, eccentric giant planet orbiting around a tightly packed inner pair of massive planets likely involved in an asymmetric 2:1 mean motion resonance. HD 27894 may be an important milestone for probing planetary formation and evolution scenarios.

At a distance of 1.295 parsecs, the red-dwarf Proxima Centauri ($\alpha$ Centauri C, GL 551, HIP 70890, or simply Proxima) is the Sun's closest stellar neighbor and one of the best studied low-mass stars. It has an effective temperature of only $\sim$ 3050 K, a luminosity of $\sim$0.1 per cent solar, a measured radius of 0.14 R$_\odot$ and a mass of about 12 per cent the mass of the Sun. Although Proxima is considered a moderately active star, its rotation period is $\sim$ 83 days, and its quiescent activity levels and X-ray luminosity are comparable to the Sun's. New observations reveal the presence of a small planet orbiting Proxima with a minimum mass of 1.3~Earth masses and an orbital period of $\sim$11.2 days. Its orbital semi-major axis is $\sim0.05$ AU, with an equilibrium temperature in the range where water could be liquid on its surface.

In radial velocity observations, a pair of extrasolar planets near a 2:1 orbital resonance can be misinterpreted as a single eccentric planet, if data are sparse and measurement precision insufficient to distinguish between these models. We determine the fraction of alleged single-planet RV detected systems for which a 2:1 resonant pair of planets is also a viable model and address the question of how the models can be disentangled. By simulation we quantified the mismatch arising from applying the wrong model. Model alternatives are illustrated using the supposed single-planet system HD 27894 for which we also study the dynamical stability of near-2:1 resonant solutions. From the data scatter around the fitted single-planet Keplerians, we find that for $74\% $ of the $254$ putative single-planet systems, a 2:1 resonant pair cannot be excluded as a viable model, since the error due to the wrong model is smaller than the scatter. For $187$ stars $\chi ^2$-probabilities can be used to reject the Keplerian models with a confidence of $95\% $ for $54\% $ of the stars and with $99.9\% $ for $39\% $ of the stars. For HD 27894 a considerable fit improvement is obtained when adding a low-mass planet near half the orbital period of the known Jovian planet. Dynamical analysis demonstrates that this system is stable when both planets are initially placed on circular orbits. For fully Keplerian orbits a stable system is only obtained if the eccentricity of the inner planet is constrained to $<0.3$. A large part of the allegedly RV detected single-planet systems should be scrutinized in order to determine the fraction of systems containing near-2:1 resonant pairs of planets. Knowing the abundance of such systems will allow us to revise the eccentricity distribution for extrasolar planets and provide direct constraints for planetary system formation.

The full LINC-NIRVANA instrument will be one of the most complex ground-based astronomical systems ever built. It will consist of multiple subsystems, including two multi-conjugate ground layer AO systems (MCAO) that drive the LBT adaptive secondaries, two mid-high layer AO systems with their own Xynetics 349 actuator DM's , a fringe tracker, a beam combiner, and the NIR science camera. In order to mitigate risk, we take a modular approach to instrument testing and commissioning by decoupling these subsystems individually. The first subsystem tested on-sky will be one of the ground-layer AO systems, part of a test-bed known as the Pathfinder. The Pathfinder consists of a 12-star pyramid wavefront sensor (PWFS) that drives one of the LBT's adaptive secondaries, a support structure known as "The Foot," and the infrared test camera (IRTC), which is used for acquisition and alignment. The 12 natural guide stars are acquired by moveable arms called "star enlargers," each of which contains its own optical path. The Pathfinder was shipped from MPIA in Heidelberg, Germany to the LBT mountain lab on Mt. Graham, Arizona in February 2013. The system was unpacked, assembled in the LBT clean room, and internally optically aligned. We present the results of our system tests, including star enlarger alignment and system alignment. We also present our immediate plans for on-sky closed loop tests on the LBT scheduled for late Fall. Because plans for all ELTs call for ground layer correction, the Pathfinder provides valuable preliminary information not only for the full LINC-NIRVANA system, but also for future advanced MCAO systems.

We present descriptions of the alignment and calibration tests of the Pathfinder, which achieved first light during our 2013 commissioning campaign at the LBT. The full LINC-NIRVANA instrument is a Fizeau interferometric imager with fringe tracking and 2-layer natural guide star multi-conjugate adaptive optics (MCAO) systems on each eye of the LBT. The MCAO correction for each side is achieved using a ground layer wavefront sensor that drives the LBT adaptive secondary mirror and a mid-high layer wavefront sensor that drives a Xinetics 349 actuator DM conjugated to an altitude of 7.1 km. When the LINC-NIRVANA MCAO system is commissioned, it will be one of only two such systems on an 8-meter telescope and the only such system in the northern hemisphere. In order to mitigate risk, we take a modular approach to commissioning by decoupling and testing the LINC-NIRVANA subsystems individually. The Pathfinder is the ground-layer wavefront sensor for the DX eye of the LBT. It uses 12 pyramid wavefront sensors to optically co-add light from natural guide stars in order to make four pupil images that sense ground layer turbulence. Pathfinder is now the first LINC-NIRVANA subsystem to be fully integrated with the telescope and commissioned on sky. Our 2013 commissioning campaign consisted of 7 runs at the LBT with the tasks of assembly, integration and communication with the LBT telescope control system, alignment to the telescope optical axis, off-sky closed loop AO calibration, and finally closed loop on-sky AO. We present the programmatics of this campaign, along with the novel designs of our alignment scheme and our off-sky calibration test, which lead to the Pathfinder's first on-sky closed loop images.

With time-series spectroscopic observations taken with the Near Infrared Spectrometer (NIRSPEC) at Keck II, we investigated the atmosphere of the close orbiting transiting extrasolar giant planet, HD 189733b. In particular, we intended to measure the dense absorption line forest around 2.3 micron, which is produced by carbon monoxide (CO). CO is expected to be present in the planetary atmosphere, although no detection of this molecule has been claimed yet. To identify the best suited data analysis method, we created artificial spectra of planetary atmospheres and analyzed them by three approaches found in the literature, the deconvolution method, data modeling via chi^2-minimization, and cross-correlation. As a result, we found that cross-correlation and chi^2-data modeling show systematically a higher sensitivity than the deconvolution method. We analyzed the NIRSPEC data with cross-correlation and detect CO absorption in the day-side spectrum of HD 189733b at the known planetary radial velocity semi-amplitude with 3.4 sigma confidence.

In 1992 we began a precision radial velocity (RV) survey for planets around solar-like stars with the Coude Echelle Spectrograph and the Long Camera (CES LC) at the 1.4 m telescope in La Silla (Chile). We have continued the survey with the upgraded CES Very Long Camera (VLC) and HARPS, both at the 3.6 m telescope, until 2007. The observations for 31 stars cover a time span of up to 15 years and the RV precision permit a search for Jupiter analogues. We perform a joint analysis for variability, trends, periodicities, and Keplerian orbits and compute detection limits. Moreover, the HARPS RVs are analysed for correlations with activity indicators (CaII H&K and CCF shape). We achieve a long-term RV precision of 15 m/s (CES+LC, 1992-1998), 9 m/s (CES+VLC, 1999-2006), and 2.8 m/s (HARPS, 2003-2009, including archive data), resp. This enables us to confirm the known planets around Iota Hor, HR 506, and HR 3259. A steady RV trend for Eps Ind A can be explained by a planetary companion. On the other hand, we find previously reported trends to be smaller for Beta Hyi and not present for Alp Men. The candidate planet Eps Eri b was not detected despite our better precision. Also the planet announced for HR 4523 cannot be confirmed. Long-term trends in several of our stars are compatible with known stellar companions. We provide a spectroscopic orbital solution for the binary HR 2400 and refined solutions for the planets around HR 506 and Iota Hor. For some other stars the variations could be attributed to stellar activity. The occurrence of two Jupiter-mass planets in our sample is in line with the estimate of 10% for the frequency of giant planets with periods smaller than 10 yr around solar-like stars. We have not detected a Jupiter analogue, while the detections limits for circular orbits indicate at 5 AU a sensitivity for minimum mass of at least 1 M_Jup (2 M_Jup) for 13% (61%) of the stars.

Massive planets in very close orbits around their central stars can induce so-called star-planet interactions (SPI), which may be of magnetic or gravitational nature. In both cases, SPI can potentially cause recurring chromospheric emission on the host star visible in Ca II H & K and/or H$\alpha$. The emission would be bound to the planetary orbit, not to the rotation period of the star. We searched for SPI in a sample of 7 stars with massive close-in planets using high-resolution spectroscopic data taken at HRS (HET) and FEROS (La Silla). We find no periodically recurring emission in the planet-hosting stars. In the case of HD 41004 AB, a binary system consisting of a K dwarf and an M dwarf, where the M dwarf is orbited by a brown dwarf companion, we find signs of cyclic variation in the Ca II K and H$\alpha$ emission lines that could be associated to interactions between the M dwarf and its companion. We present our first results of this interesting system that may become an important system for the understanding of SPI.

We conduct a study of K to M type stars to investigate the activity and the rotation limit in the Hyades. We use a sample of 40 stars in this intermediate-age cluster (~625 Myr) to probe stellar rotation in the threshold region where stellar activity becomes prevalent. Here we present projected equatorial velocities (vsin i) and chromospheric activity measurements (H{\alpha}) that indicate the existence of fast rotators in the Hyades at spectral types where also the fraction of stars with H{\alpha} emission shows a rapid increase ("H{\alpha} limit"). The locus of enhanced rotation (and activity) thus seems to be shifted to earlier types in contrast to what is seen as the rotation limit in field stars. The relation between activity and rotation appears to be similar to the one observed in fields stars.

We present radial velocity (RV) measurements of our sample of 40 M dwarfs from our planet search programme with VLT+UVES begun in 2000. Although with our RV precision down to 2 - 2.5 m/s and timebase line of up to 7 years, we are capable of finding planets of a few Earth masses in the close-in habitable zones of M dwarfs, there is no detection of a planetary companion. To demonstrate this we present mass detection limits allowing us to exclude Jupiter-mass planets up to 1 AU for most of our sample stars. We identified 6 M dwarfs that host a brown dwarf or low-mass stellar companion. With the exception of these, all other sample stars show low RV variability with an rms < 20 m/s. Some high proper motion stars exhibit a linear RV trend consistent with their secular acceleration. Furthermore, we examine our data sets for a possible correlation between RVs and stellar activity as seen in variations of the Halpha line strength. For Barnard's star we found a significant anticorrelation, but most of the sample stars do not show such a correlation.

The Lomb-Scargle periodogram is a common tool in the frequency analysis of unequally spaced data equivalent to least-squares fitting of sine waves. We give an analytic solution for the generalisation to a full sine wave fit, including an offset and weights ($\chi^{2}$ fitting). Compared to the Lomb-Scargle periodogram, the generalisation is superior as it provides more accurate frequencies, is less susceptible to aliasing, and gives a much better determination of the spectral intensity. Only a few modifications are required for the computation and the computational effort is similar. Our approach brings together several related methods that can be found in the literature, viz. the date-compensated discrete Fourier transform, the floating-mean periodogram, and the "spectral significance" estimator used in the SigSpec program, for which we point out some equivalences. Furthermore, we present an algorithm that implements this generalisation for the evaluation of the Keplerian periodogram that searches for the period of the best-fitting Keplerian orbit to radial velocity data. The systematic and non-random algorithm is capable of detecting eccentric orbits, which is demonstrated by two examples and can be a useful tool in searches for the orbital periods of exoplanets.

Context. Brown dwarf companions to stars at separations of a few AU or less are rare objects, and none have been found so far around early-type M dwarfs M0V-M5V). With GJ 1046 (M2.5V), a strong candidate for such a system with a separation of 0.42 AU is presented. Aims. We aim at constraining the mass of the companion in order to decide whether it is a brown dwarf or a low-mass star. Methods. We employed precision RV measurements to determine the orbital parameters and the minimum companion mass. We then derived an upper limit to the companion mass from the lack of disturbances of the RV measurements by a secondary spectrum. An even tighter upper limit is subsequently established by combining the RV-derived orbital parameters with the recent new version of the Hipparcos Intermediate Astrometric Data. Results. For the mass of the companion, we derive m>26.9 MJup from the RV data. Based on the RV data alone, the probability that the companion exceeds the stellar mass threshold is just 6.2%. The absence of effects from the secondary spectrum lets us constrain the companion mass to m <229 MJup. The combination of RV and Hipparcos data yields a 3sigma upper mass limit to the companion mass of 112 MJup with a formal optimum value at m=47.2 MJup. From the combination of RV and astrometric data, the chance probability that the companion is a star is 2.9%. Conclusions. We have found a low-mass, close companion to an early-type M dwarf. While the most likely interpretation of this object is that it is a brown dwarf, a low-mass stellar companion is not fully excluded.

We monitored the chromospheric activity in the Ca II H & K lines of 13 solar-type stars (including the Sun); 8 of them over three years at the CFHT and 5 in a single run at the VLT. Ten of the 13 targets have close planetary companions. All of the stars observed at the CFHT show long-term (months to years) changes in H & K intensity levels. Four stars display short-term (days) cyclical activity. For two, HD 73256 and kappa^1 Ceti, the activity is likely associated with an active region rotating with the star, however, the flaring in excess of the rotational modulation may be associated with a hot jupiter. A planetary companion remains a possibility for kappa^1 Ceti. For the other two, HD 179949 and upsilon And, the cyclic variation is synchronized to the hot jupiter's orbit. For both stars this synchronicity with the orbit is clearly seen in two out of three epochs. The effect is only marginal in the third epoch at which the seasonal level of chromospheric activity had changed for both stars. Short-term chromospheric activity appears weakly dependent on the mean K-line reversal intensities for the sample of 13 stars. Also, a suggestive correlation exists between this activity and the M_p sin(i) of the star's hot jupiter. Because of their small separation (<= 0.1 AU), many of the hot jupiters lie within the Alfv\'en radius of their host stars which allows a direct magnetic interaction with the stellar surface. We discuss the conditions under which a planet's magnetic field might induce activity on the stellar surface and why no such effect was seen for the prime candidate, tau Boo. This work opens up the possibility of characterizing planet-star interactions, with implications for extrasolar planet magnetic fields and the energy contribution to stellar atmospheres.

High precision radial velocity (RV) measurements spanning the years 1980.8--2000.0 are presented for the nearby (3.22 pc) K2 V star $\epsilon$ Eri. These data, which represent a combination of six independent data sets taken with four different telescopes, show convincing variations with a period of $\approx$ 7 yrs. A least squares orbital solution using robust estimation yields orbital parameters of period, $P$ = 6.9 yrs, velocity $K$-amplitude $=$ 19 {\ms}, eccentricity $e$ $=$ 0.6, projected companion mass $M$ sin $i$ = 0.86 $M_{Jupiter}$, and semi-major axis $a_2$ $=$ 3.3 AU. Ca II H&K S-index measurements spanning the same time interval show significant variations with periods of 3 and 20 yrs, yet none at the RV period. If magnetic activity were responsible for the RV variations then it produces a significantly different period than is seen in the Ca II data. Given the lack of Ca II variation with the same period as that found in the RV measurements, the long-lived and coherent nature of these variations, and the high eccentricity of the implied orbit, Keplerian motion due to a planetary companion seems to be the most likely explanation for the observed RV variations. The wide angular separation of the planet from the star (approximately 1 arc-second) and the long orbital period make this planet a prime candidate for both direct imaging and space-based astrometric measurements.

We present differential radial velocity measurements of Proxima Centauri collected over 4 years with the ESO CES with a mean precision of $54 ms^{-1}$. We find no evidence of a periodic signal that could corroborate the existence of a sub-stellar companion. We put upper limits (97% confidence) to the companion mass ranging from 1.1 to 22 M_{Jup} at orbital periods of 0.75 to 3000 d, i.e. separations 0.008-2 AU from Prox Cen. Our mass limits concur with limits found by precise astrometry (Benedict et al. 1998a and priv. comm.) which strongly constrain the period range 50-600 d to 1.1-0.22 M_{Jup}. Combining both results we exclude a brown dwarf or supermassive planet at separations 0.008-0.69 AU from Prox Cen. We also find that, at the level of our precision, the RV data are not affected by stellar activity.