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 nearest exoplanets to the Sun are our best possibilities for detailed characterization. We report the discovery of a compact multi-planet system of super-Earths orbiting the nearby red dwarf GJ 887, using radial velocity measurements. The planets have orbital periods of 9.3 and 21.8~days. Assuming an Earth-like albedo, the equilibrium temperature of the 21.8 day planet is approx 350 K; which is interior, but close to the inner edge, of the liquid-water habitable zone. We also detect a further unconfirmed signal with a period of 50 days which could correspond to a third super-Earth in a more temperate orbit. GJ 887 is an unusually magnetically quiet red dwarf with a photometric variability below 500 parts-per-million, making its planets amenable to phase-resolved photometric characterization.

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

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.

High resolution transit spectroscopy has proven to be a reliable technique for the characterization of the chemical composition of exoplanet atmospheres. Taking advantage of the broad spectral coverage of the CARMENES spectrograph, we initiated a survey aimed at characterizing a broad range of planetary systems. Here, we report our observations of three transits of \tplanet with CARMENES in search of \het\ absorption. On one of the nights, the He~{\sc i} region was heavily contaminated by OH$^-$ telluric emission and, thus, it was not useful for our purposes. The remaining two nights had a very different signal-to-noise ratio (S/N) due to weather. They both indicate the presence of \het\ absorption in the transmission spectrum of \tplanet, although a statistically valid detection can only be claimed for the night with higher S/N. For that night, we retrieved a 1.5$\pm$0.3\% absorption depth, translating into a $R_p(\lambda)/R_p = 1.15\pm 0.14$ at this wavelength. Spectro-photometric light curves for this same night also indicate the presence of extra absorption during the planetary transit with a consistent absorption depth. The \het\ absorption is modeled in detail using a radiative transfer code, and the results of our modeling efforts are compared to the observations. We find that the mass-loss rate, \mlr, is confined to a range of 3\,$\times\,10^{10}$\,\gs\ for $T$ = 6000\,K to 10\,$\times\,10^{10}$\,\gs\ for $T$ = 9000\,K. We discuss the physical mechanisms and implications of the He~{\sc i} detection in \tplanet and put it in context as compared to similar detections and non-detections in other Neptune-size planets. We also present improved stellar and planetary parameter determinations based on our visible and near-infrared observations.

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.

HD 209458b was the first transiting planet discovered, and the first for which an atmosphere, in particular NaI, was detected. With time, it has become one of the most frequently studied planets, with a large diversity of atmospheric studies using low- and high-resolution spectroscopy. Here, we present transit spectroscopy observations of HD 209458b using the HARPS-N and CARMENES spectrographs. We fit the Rossiter-McLaughlin effect by combining radial velocity data from both instruments (nine transits in total), measuring a projected spin-orbit angle of $-1.6\pm0.3~{\rm deg}$. We also present the analysis of high-resolution transmission spectroscopy around the NaI region at $590~{\rm nm}$, using a total of five transit observations. In contrast to previous studies where atmospheric NaI absorption is detected, we find that for all of the nights, whether individually or combined, the transmission spectra can be explained by the combination of the centre-to-limb variation and the Rossiter-McLaughlin effect. This is also observed in the time-evolution maps and transmission light curves, but at lower signal-to-noise ratio. Other strong lines such as H$\alpha$, CaII IRT, the MgI triplet region, and KI D1 are analysed, and are also consistent with the modelled effects, without considering any contribution from the exoplanet atmosphere. Thus, the transmission spectrum reveals no detectableNaI absorption in HD 209458b. We discuss how previous pioneering studies of this benchmark object may have overlooked these effects. While for some star-planet systems these effects are small, for other planetary atmospheres the results reported in the literature may require revision.

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.

The star Kepler-1625 recently attracted considerable attention when an analysis of the stellar photometric time series from the Kepler mission was interpreted as showing evidence of a large exomoon around the transiting Jupiter-sized planet candidate Kepler-1625b. We aim to detect the radial velocity (RV) signal imposed by Kepler-1625b (and its putative moon) on the host star or, as the case may be, determine an upper limit on the mass of the transiting object. We took a total of 22 spectra of Kepler-1625 using CARMENES, 20 of which were useful. Observations were spread over a total of seven nights between October 2017 and October 2018, covering $125\%$ of one full orbit of Kepler-1625b. We used the automatic Spectral Radial Velocity Analyser (SERVAL) pipeline to deduce the stellar RVs and uncertainties. Then we fitted the RV curve model of a single planet on a Keplerian orbit to the observed RVs using a $\chi^2$ minimisation procedure. We derive upper limits on the mass of Kepler-1625b under the assumption of a single planet on a circular orbit. In this scenario, the $1\,\sigma$, $2\,\sigma$, and $3\,\sigma$ confidence upper limits for the mass of Kepler-1625b are $2.90\,M_{\rm J}$, $7.15\,M_{\rm J}$, and $11.60\,M_{\rm J}$, respectively. We present strong evidence for the planetary nature of Kepler-1625b, making it the 10th most long-period confirmed planet known today. Our data does not answer the question about a second, possibly more short-period planet that could be responsible for the observed transit timing variation of Kepler-1625b.

Context. The HARPS spectrograph provides state-of-the-art stellar radial velocity (RV) measurements with a precision down to 1 m/s. The spectra are extracted with a dedicated data-reduction software (DRS) and the RVs are computed by CCF with a numerical mask. Aims. The aim of this study is three-fold: (i) Create easy access to the public HARPS RV data set. (ii) Apply the new public SERVAL pipeline to the spectra, and produce a more precise RV data set. (iii) Check whether the precision of the RVs can be further improved by correcting for small nightly systematic effects. Methods. For each star observed with HARPS, we downloaded the publicly available spectra from the ESO archive and recomputed the RVs with SERVAL. We then computed nightly zero points (NZPs) by averaging the RVs of quiet stars. Results. Analysing the RVs of the most RV-quiet stars, whose RV scatter is < 5 m/s, we find that SERVAL RVs are on average more precise than DRS RVs by a few percent. We find three significant systematic effects, whose magnitude is independent of the software used for the RV derivation: (i) stochastic variations with a magnitude of 1 m/s; (ii) long-term variations, with a magnitude of 1 m/s and a typical timescale of a few weeks; and (iii) 20-30 NZPs significantly deviating by a few m/s. In addition, we find small (< 1 m/s) but significant intra-night drifts in DRS RVs before the 2015 intervention, and in SERVAL RVs after it. We confirm that the fibre exchange in 2015 caused a discontinuous RV jump, which strongly depends on the spectral type of the observed star: from 14 m/s for late F-type stars, to -3 m/s for M dwarfs. Conclusions. Our NZP-corrected SERVAL RVs can be retrieved from a user-friendly, public database. It provides more than 212 000 RVs for about 3000 stars along with many auxiliary information, NZP corrections, various activity indices, and DRS-CCF products.

The determination of extrasolar planet masses with the radial velocity (RV) technique requires spectroscopic Doppler information from the planet's host star, which varies with stellar brightness and temperature. We analyze Doppler information in spectra of F--M dwarfs utilizing empirical information from HARPS and CARMENES, and from model spectra. We come to the conclusions that an optical setup ($BVR$-bands) is more efficient that a near-infrared one ($YJHK$) in dwarf stars hotter than 3200\,K. We publish a catalogue of 46,480 well-studied F--M dwarfs in the solar neighborhood and compare their distribution to more than one million stars from Gaia DR2. For all stars, we estimate the RV photon noise achievable in typical observations assuming no activity jitter and slow rotation. We find that with an ESPRESSO-like instrument at an 8m-telescope, a photon noise limit of 10\,cm\,s$^{-1}$ or lower can be reached in more than 280 stars in a 5\,min observation. At 4m-telescopes, a photon noise limit of 1\,m\,s$^{-1}$ can be reached in a 10\,min exposure in approx.\ 10,000 predominantly sun-like stars with a HARPS-like (optical) instrument. The same applies to $\sim$3000 stars for a red-optical setup covering the $RIz$-bands, and to $\sim$700 stars for a near-infrared instrument. For the latter two, many of the targets are nearby M dwarfs. Finally, we identify targets in which Earth-mass planets within the liquid water habitable zone can cause RV amplitudes comparable to the RV photon noise. Assuming the same exposure times, we find that an ESPRESSO-like instrument can reach this limit for 1\,M$_\Earth$ planets in more than 1000 stars. The optical, red-optical, and near-infrared configurations reach the limit for 2\,M$_\Earth$ planets in approximately 500, 700, and 200 stars, respectively.

Ultra-hot Jupiters are emerging as a new class of exoplanets. Studying their chemical compositions and temperature structures will improve the understanding of their mass loss rate as well as their formation and evolution. We present the detection of ionized calcium in the two hottest giant exoplanets - KELT-9b and WASP-33b. By utilizing transit datasets from CARMENES and HARPS-N observations, we achieved high confidence level detections of Ca II using the cross-correlation method. We further obtain the transmission spectra around the individual lines of the Ca II H&K doublet and the near-infrared triplet, and measure their line profiles. The Ca II H&K lines have an average line depth of 2.02 $\pm$ 0.17 % (effective radius of 1.56 Rp) for WASP-33b and an average line depth of 0.78 $\pm$ 0.04 % (effective radius of 1.47 Rp) for KELT-9b, which indicates that the absorptions are from very high upper atmosphere layers close to the planetary Roche lobes. The observed Ca II lines are significantly deeper than the predicted values from the hydrostatic models. Such a discrepancy is probably a result of hydrodynamic outflow that transports a significant amount of Ca II into the upper atmosphere. The prominent Ca II detection with the lack of significant Ca I detection implies that calcium is mostly ionized in the upper atmospheres of the two planets.

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.

Small low-mass stars are favourable targets for the detection of rocky habitable planets. In particular, planetary systems in the solar neighbourhood are interesting and suitable for precise characterisation. The Red Dots campaigns seek to discover rocky planets orbiting nearby low-mass stars. The 2018 campaign targeted GJ 1061, which is the 20$^{\rm th}$ nearest star to the Sun. For three consecutive months we obtained nightly, high-precision radial velocity measurements with the HARPS spectrograph. We analysed these data together with archival HARPS data. We report the detection of three planet candidates with periods of $3.204\pm 0.001$, $6.689\pm 0.005$ and $13.03\pm 0.03$ days, which is close to 1:2:4 period commensurability. After several considerations related to the properties of the noise and sampling, we conclude that a 4$^{\rm th}$ signal is most likely explained by stellar rotation, although it may be due to a planet. The proposed three-planet system (and the potential four-planet solution) is long-term dynamically stable. Planet-planet gravitational interactions are below our current detection threshold. The minimum masses of the three planets range from $1.4\pm 0.2$ to $1.8\pm 0.3$ M$_\oplus$. Planet d, with $m \sin i = 1.68\pm0.25$ M$_\oplus$, receives a similar amount of energy as Earth receives from the Sun. Consequently it lies within the liquid-water habitable zone of the star and has a similar equilibrium temperature to Earth. GJ 1061 has very similar properties to Proxima Centauri but activity indices point to lower levels of stellar activity.

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.

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.

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.

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.

The HIRES spectrograph, mounted on the $10$-m Keck-I telescope, belongs to a small group of radial-velocity (RV) instruments that produce stellar RVs with long-term precision down to $\sim1$ ms$^{-1}$. In $2017$, the HIRES team published $64,480$ RVs of $1,699$ stars, collected between $1996$ and $2014$. In this bank of RVs, we identify a sample of RV-quiet stars, whose RV scatter is $<10$ ms$^{-1}$, and use them to reveal two small but significant nightly zero-point effects: a discontinuous jump, caused by major modifications of the instrument in August $2004$, and a long-term drift. The size of the $2004$ jump is $1.5\pm0.1$ ms$^{-1}$, and the slow zero-point variations have a typical magnitude of $\lesssim1$ ms$^{-1}$. In addition, we find a small but significant correlation between stellar RVs and the time relative to local midnight, indicative of an average intra-night drift of $0.051\pm0.004$ ms$^{-1}$hr$^{-1}$. We correct the $64,480$ HIRES RVs for the systematic effects we find, and make the corrected RVs publicly available. Our findings demonstrate the importance of observing RV-quiet stars, even in the era of simultaneously-calibrated RV spectrographs. We hope that the corrected HIRES RVs will facilitate the search for new planet candidates around the observed stars.

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.

Hot gas giant exoplanets can lose part of their atmosphere due to strong stellar irradiation, affecting their physical and chemical evolution. Studies of atmospheric escape from exoplanets have mostly relied on space-based observations of the hydrogen Lyman-{\alpha} line in the far ultraviolet which is strongly affected by interstellar absorption. Using ground-based high-resolution spectroscopy we detect excess absorption in the helium triplet at 1083 nm during the transit of the Saturn-mass exoplanet WASP-69b, at a signal-to-noise ratio of 18. We measure line blue shifts of several km/s and post transit absorption, which we interpret as the escape of part of the atmosphere trailing behind the planet in comet-like form. [Additional notes by authors: Furthermore, we provide upper limits for helium signals in the atmospheres of the exoplanets HD 209458b, KELT-9b, and GJ 436b. We investigate the host stars of all planets with detected helium signals and those of the three planets we derive upper limits for. In each case we calculate the X-ray and extreme ultraviolet flux received by these planets. We find that helium is detected in the atmospheres of planets (orbiting the more active stars and) receiving the larger amount of irradiation from their host stars.]

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.

Context. Many algorithms to solve Kepler's equations require the evaluation of trigonometric or root functions. Aims. We present an algorithm to compute the eccentric anomaly and even its cosine and sine terms without usage of other transcendental functions at run-time. With slight modifications it is applicable for the hyperbolic case, too. Methods. Based on the idea of CORDIC, it requires only additions and multiplications and a short table. The table is independent of eccentricity and can be hardcoded. Its length depends on the desired precision. Results. The code is short. The convergence is linear for all mean anomalies and eccentricities e (including e = 1). As a stand-alone algorithm, single and double precision is obtained with 29 and 55 iterations, respectively. One half or two third of the iterations can be saved in combination with Newton's or Halley's method at the cost of one division.

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

CARMENES is a spectrograph for radial velocity surveys of M dwarfs with the aim of detecting Earth-mass planets orbiting in the habitable zones of their host stars. To ensure an optimal use of the CARMENES Guaranteed Time Observations, in this paper we investigate the correlation of activity and rotation for approximately 2200 M dwarfs, ranging in spectral type from M0.0 V to M9.0 V. We present new high-resolution spectroscopic observations with FEROS, CAFE, and HRS of approximately 500 M dwarfs. For each new observation, we determined its radial velocity and measured its Halpha activity index and its rotation velocity. Additionally, we have multiple observations of many stars to investigate if there are any radial velocity variations due to multiplicity. The results of our survey confirm that early-M dwarfs are Halpha inactive with low rotational velocities and that late-M dwarfs are Halpha active with very high rotational velocities. The results of this high-resolution analysis comprise the most extensive catalogue of rotation and activity in M dwarfs currently available.

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.

Stellar oscillations appear all across the Hertzsprung-Russell diagram. Recent theoretical studies support their existence also in the atmospheres of M dwarfs. These studies predict for them short periodicities ranging from 20~min to 3~h. Our Cool Tiny Beats (CTB) programme aims at finding these oscillations for the very first time. With this goal, CTB explores the short time domain of M dwarfs using radial velocity data from the HARPS-ESO and HARPS-N high-precision spectrographs. Here we present the results for the two most long-term stable targets observed to date with CTB, GJ~588 and GJ~699 (i.e. Barnard's star). In the first part of this work we detail the correction of several instrumental effects. These corrections are specially relevant when searching for sub-night signals. Results show no significant signals in the range where M dwarfs pulsations were predicted. However, we estimate that stellar pulsations with amplitudes larger than $\sim0.5\,\mathrm{m\,s}^{-1}$ can be detected with a 90% completeness with our observations. This result, along with the excess of power regions detected in the periodograms, open the possibility of non-resolved very low amplitude pulsation signals. Next generation more precise instrumentation would be required to detect such oscillations. However, the possibility of detecting pulsating M-dwarf stars with larger amplitudes is feasible due to the short size of the analysed sample. This motivates the need for completeness of the CTB survey.

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.

We have been monitoring 373 very bright (V < 6 mag) G and K giants with high precision optical Doppler spectroscopy for more than a decade at Lick Observatory. Our goal was to discover planetary companions around those stars and to better understand planet formation and evolution around intermediate-mass stars. However, in principle, long-term, g-mode nonradial stellar pulsations or rotating stellar features, such as spots, could effectively mimic a planetary signal in the radial velocity data. Our goal is to compare optical and infrared radial velocities for those stars with periodic radial velocity patterns and to test for consistency of their fitted radial velocity semiamplitudes. Thereby, we distinguish processes intrinsic to the star from orbiting companions as reason for the radial velocity periodicity observed in the optical. Stellar spectra with high spectral resolution have been taken in the H-band with the CRIRES near-infrared spectrograph at ESO's VLT for 20 stars of our Lick survey. Radial velocities are derived using many deep and stable telluric CO2 lines for precise wavelength calibration. We find that the optical and near-infrared radial velocities of the giant stars in our sample are consistent. We present detailed results for eight stars in our sample previously reported to have planets or brown dwarf companions. All eight stars passed the infrared test. We conclude that the planet hypothesis provides the best explanation for the periodic radial velocity patterns observed for these giant stars.

Stellar activity may induce Doppler variability at the level of a few m/s which can then be confused by the Doppler signal of an exoplanet orbiting the star. To first order, linear correlations between radial velocity measurements and activity indices have been proposed to account for any such correlation. The likely presence of two super-Earths orbiting Kapteyn's star was reported in Anglada et al. (2014, MNRAS 443L, 89A), but this claim was recently challenged by Robertson et al. (2015, ApJ 805L, 22R) arguing evidence of a rotation period (143 days) at three times the orbital period of one of the proposed planets (Kapteyn's b, P=48.6 days), and the existence of strong linear correlations between its Doppler signal and activity data. By re-analyzing the data using global optimization methods and model comparison, we show that such claim is incorrect given that; 1) the choice of a rotation period at 143 days is unjustified, and 2) the presence of linear correlations is not supported by the data. We conclude that the radial velocity signals of Kapteyn's star remain more simply explained by the presence of two super-Earth candidates orbiting it. We also advocate for the use of global optimization procedures and objective arguments, instead of claims lacking of a minimal statistical support.

Over the past decades hollow-cathode lamps have been calibration standards for spectroscopic measurements. Advancing to cm/s radial velocity precisions with the next generation of instruments requires more suitable calibration sources with more lines and less dynamic range problems. Fabry-Perot interferometers provide a regular and dense grid of lines and homogeneous amplitudes making them good candidates for next generation calibrators. We investigate the usefulness of Fabry-Perot etalons in wavelength calibration, present an algorithm to incorporate the etalon spectrum in the wavelength solution and examine potential problems. The quasi periodic pattern of Fabry-Perot lines is used along with a hollow-cathode lamp to anchor the numerous spectral features on an absolute scale. We test our method with the HARPS spectrograph and compare our wavelength solution to the one derived from a laser frequency comb. The combined hollow-cathode lamp/etalon calibration overcomes large distortion (50 m/s) in the wavelength solution of the HARPS data reduction software. Direct comparison to the laser frequency comb bears differences of only maximum 10 m/s. Combining hollow-cathode lamps with Fabry-Perot Interferometers can lead to substantial improvements in the wavelength calibration of echelle spectrographs. Etalons can provide economical alternatives to the laser frequency comb, especially for smaller projects.

We investigate the nature of the long-period radial velocity variations in Alpha Tau first reported over 20 years ago. We analyzed precise stellar radial velocity measurements for Alpha Tau spanning over 30 years. An examination of the Halpha and Ca II 8662 spectral lines, and Hipparcos photometry was also done to help discern the nature of the long-period radial velocity variations. Our radial velocity data show that the long-period, low amplitude radial velocity variations are long-lived and coherent. Furthermore, Halpha equivalent width measurements and Hipparcos photometry show no significant variations with this period. Another investigation of this star established that there was no variability in the spectral line shapes with the radial velocity period. An orbital solution results in a period of P = 628.96 +/- 0.90 d, eccentricity, e = 0.10 +/- 0.05, and a radial velocity amplitude, K = 142.1 +/- 7.2 m/s. Evolutionary tracks yield a stellar mass of 1.13 +/- 0.11 M_sun, which corresponds to a minimum companion mass of 6.47 +/- 0.53 M_Jup with an orbital semi-major axis of a = 1.46 +/- 0.27 AU. After removing the orbital motion of the companion, an additional period of ~ 520 d is found in the radial velocity data, but only in some time spans. A similar period is found in the variations in the equivalent width of Halpha and Ca II. Variations at one-third of this period are also found in the spectral line bisector measurements. The 520 d period is interpreted as the rotation modulation by stellar surface structure. Its presence, however, may not be long-lived, and it only appears in epochs of the radial velocity data separated by $\sim$ 10 years. This might be due to an activity cycle. The data presented here provide further evidence of a planetary companion to Alpha Tau, as well as activity-related radial velocity variations.

Exoplanets of a few Earth masses can be now detected around nearby low-mass stars using Doppler spectroscopy. In this paper, we investigate the radial velocity variations of Kapteyn's star, which is both a sub-dwarf M-star and the nearest halo object to the Sun. The observations comprise archival and new HARPS, HIRES and PFS Doppler measurements. Two Doppler signals are detected at periods of 48 and 120 days using likelihood periodograms and a Bayesian analysis of the data. Using the same techniques, the activity indicies and archival ASAS-3 photometry show evidence for low-level activity periodicities of the order of several hundred days. However, there are no significant correlations with the radial velocity variations on the same time-scales. The inclusion of planetary Keplerian signals in the model results in levels of correlated and excess white noise that are remarkably low compared to younger G, K and M dwarfs. We conclude that Kapteyn's star is most probably orbited by two super-Earth mass planets, one of which is orbiting in its circumstellar habitable zone, becoming the oldest potentially habitable planet known to date. The presence and long-term survival of a planetary system seems a remarkable feat given the peculiar origin and kinematic history of Kapteyn's star. The detection of super-Earth mass planets around halo stars provides important insights into planet-formation processes in the early days of the Milky Way.

Optimal extraction is a key step in processing the raw images of spectra as registered by two-dimensional detector arrays to a one-dimensional format. Previously reported algorithms reconstruct models for a mean one-dimensional spatial profile to assist a properly weighted extraction. We outline a simple optimal extraction algorithm including error propagation, which is very suitable for stabilised, fibre-fed spectrographs and does not model the spatial profile shape. A high signal-to-noise, master-flat image serves as reference image and is directly used as an extraction profile mask. Each extracted spectral value is the scaling factor relative to the cross-section of the unnormalised master-flat which contains all information about the spatial profile as well as pixel-to-pixel variations, fringing, and blaze. The extracted spectrum is measured relative to the flat spectrum. Using echelle spectra of the HARPS spectrograph we demonstrate a competitive extraction performance in terms of signal-to-noise and show that extracted spectra can be used for high precision radial velocity measurement. Pre- or post-flat-fielding of the data is not necessary, since all spectrograph inefficiencies inherent to the extraction mask are automatically accounted for. Also the reconstruction of the mean spatial profile by models is not needed, thereby reducing the number of operations to extract spectra. Flat-relative optimal extraction is a simple, efficient, and robust method that can be applied easily to stabilised, fibre-fed spectrographs.

[ABRIDGED]: Aims: two candidates in the KOI 1574 system are relatively long-period (about 114d and 191d) and in 5:3 resonance. We therefore search for TTVs in this particularly promising system. Methods: The full Kepler data was used, allowing to search for TTVs as well as for additional transit-like signals. Results: We detect strong anti-correlated TTVs of the 114d and 191d signals, dynamically confirming them as members of the same system. Dynamical simulations reproducing the observed TTVs allow us to also determine the masses of the planets. KOI 1574.01 (hereafter Kepler-87 b) was found to have a radius of 13.49 +/- 0.55 R_earth and a mass of 324.2 +/- 8.8M_earth, and KOI 1574.02 (Kepler-87 c) was found to have a radius of 6.14 +/- 0.29R_earth and a mass of 6.4 +/- 0.8M_earth. Both planets have low densities of 0.729 and 0.152 g cm^-3, respectively, which is non-trivial for such cold and old (7-8 Gyr) planets. Specifically, Kepler-87 c is the lowest- density planet in the super-Earth mass range. Both planets are thus particularly amenable to modeling and planetary structure studies, and also present an interesting case were ground-based photometric follow-up of Kepler planets is very desirable. Finally, we also detect two more short period super-Earth sized planetary (< 2R_earth) candidates in the system, making the relatively high multiplicity of this system notable against the general paucity of multiple systems in the presence of giant planets like Kepler-87 b.

Stellar activity signatures such as spots and plage can significantly limit the search for extrasolar planets. Current models of activity-induced radial velocity (RV) signals focused on the impact of temperature contrast in spots predicting the signal to diminish toward longer wavelengths. On the other hand, the relative importance of the Zeeman effect on RV measurements should grow with wavelength, because the Zeeman displacement itself grows with \lambda, and because a magnetic and cool spot contributes more to the total flux at longer wavelengths. We model the impact of active regions on stellar RV measurements including both temperature contrast in spots and Zeeman line broadening. We calculate stellar line profiles using polarized radiative transfer models including atomic and molecular Zeeman splitting from 0.5 to 2.3\mum. Our results show that the amplitude of the RV signal caused by the Zeeman effect alone can be comparable to that caused by temperature contrast. Furthermore, the RV signal caused by cool and magnetic spots increases with wavelength. We also calculate the RV signal due to variations in average magnetic field strength from one observation to the next, but find it unlikely that this can significantly influence the search for extrasolar planets. As an example, we derive the RV amplitude of the active M dwarf AD Leo as a function of wavelength using data from the HARPS spectrograph. The RV signal does not diminish at longer wavelengths but shows evidence for the opposite behavior. We conclude that the RV signal of active stars does not vanish at longer wavelength but sensitively depends on the combination of spot temperature and magnetic field; in active low-mass stars, it is even likely to grow with wavelength.

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.

We use precise radial velocity measurements and photometric data to derive the frequency spacing of the p-mode oscillation spectrum of the planet-hosting star Beta Gem. This spacing along with the interferometric radius for this star is used to derive an accurate stellar mass. A long time series of over 60 hours of precise stellar radial velocity measurements of Beta Gem were taken with an iodine absorption cell and the echelle spectrograph mounted on the 2m Alfred Jensch Telescope. Complementary photometric data for this star were also taken with the MOST microsatellite spanning 3.6 d. A Fourier analysis of the radial velocity data reveals the presence of up to 17 significant pulsation modes in the frequency interval 10-250 micro-Hz. Most of these fall on a grid of equally-spaced frequencies having a separation of 7.14 +/- 0.12 micro-Hz. An analysis of 3.6 days of high precision photometry taken with the MOST space telescope shows the presence of up to 16 modes, six of which are consistent with modes found in the spectral (radial velocity) data. This frequency spacing is consistent with high overtone radial pulsations; however, until the pulsation modes are identified we cannot be sure if some of these are nonradial modes or even mixed modes. The radial velocity frequency spacing along with angular diameter measurements of Beta Gem via interferometry results in a stellar mass of M = 1.91 +/- 0.09 solar masses. This value confirms the intermediate mass of the star determined using stellar evolutionary tracks. Beta Gem is confirmed to be an intermediate mass star. Stellar pulsations in giant stars along with interferometric radius measurements can provide accurate determinations of the stellar mass of planet hosting giant stars. These can also be used to calibrate stellar evolutionary tracks.

KIC 4247791 is an eclipsing binary observed by the Kepler satellite mission. We wish to determine the nature of its components and in particular the origin of a shallow dip in its Kepler light curve that previous investigations have been unable to explain in a unique way. We analyze newly obtained high-resolution spectra of the star using synthetic spectra based on atmosphere models, derive the radial velocities of the stellar components from cross-correlation with a synthetic template, and calculate the orbital solution. We use the JKTEBOP program to model the Kepler light curve of KIC 4247791. We find KIC 4247791 to be a SB4 star. The radial velocity variations of its four components can be explained by two separate eclipsing binaries. In contradiction to previous photometric findings, we show that the observed composite spectrum as well as the derived masses of all four of its components correspond to spectral type F. The observed small dip in the light curve is not caused by a transit-like phenomenon but by the eclipses of the second binary system. We find evidence that KIC 4247791 might belong to the very rare hierarchical SB4 systems with two eclipsing binaries.

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.

Eps Ind A is one of the nearest sun-like stars, located only 3.6 pc away. It is known to host a binary brown dwarf companion, Eps Ind Ba/Bb, at a large projected separation of 6.7", but radial velocity measurements imply that an additional, yet unseen component is present in the system, much closer to Eps Ind A. Previous direct imaging has excluded the presence of any stellar or high-mass brown dwarf companion at small separations, indicating that the unseen companion may be a low-mass brown dwarf or high-mass planet. We present the results of a deep high-contrast imaging search for the companion, using active angular differential imaging (aADI) at 4 micron, a particularly powerful technique for planet searches around nearby and relatively old stars. We also develop an additional PSF reference subtraction scheme based on locally optimized combination of images (LOCI) to further enhance the detection limits. No companion is seen in the images, although we are sensitive to significantly lower masses than previously achieved. Combining the imaging data with the known radial velocity trend, we constrain the properties of the companion to within approximately 5-20 Mjup, 10-20 AU, and i > 20 deg, unless it is an exotic stellar remnant. The results also imply that the system is probably older than the frequently assumed age of ~1 Gyr.

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.

$\iota$ Dra (HIP 75458) is a well-known example for a K giant hosting a substellar companion since its discovery by Frink et al. (2002). We present radial velocity measurements of this star from observations taken with three different instruments spanning nearly 8 years. They show more clearly that the RV period is long-lived and coherent thus supporting the companion hypothesis. The longer time baseline now allows for a more accurate determination of the orbit with a revised period of P=511 d and an additional small linear trend, indicative of another companion in a wide orbit. Moreover we show that the star exhibits low amplitude, solar like oscillations with frequencies around 3-4 d$^{-1}$ (34.7-46.3 $\mu$Hz).

We present the results of a long time series of precise stellar radial velocity measurements of the planet hosting K giant star Beta Geminorum. A total of 20 hours of observations spanning three nights were obtained and the radial velocity variations show the presence of solar-like stellar oscillations. Our period analysis yields six significant pulsation modes that have frequencies in the range of 30 - 150 microHz. The dominant mode is at a frequency of 86.9 microHz and has an amplitude of 5.3 m/s. These values are consistent with stellar oscillations for a giant star with a stellar mass of approximately 2 solar masses. This stellar mass implies a companion minimum mass of 2.6 Jupiter masses. Beta Gem is the first planet hosting giant star in which multi-periodic stellar oscillations have been detected. The study of stellar oscillations in planet hosting giant stars may provide an independent, and more accurate determination of the stellar mass.