The properties of inhomogeneities on the surface of active stars (i.e. dark spots and bright faculae) significantly influence the determination of the parameters of an exoplanet. The chromatic effect they have on transmission spectroscopy could affect the analysis of data from future space missions such as JWST and Ariel. To quantify and mitigate the effects of those surface phenomena, we developed a modelling approach to derive the surface distribution and properties of active regions by modelling simultaneous multi-wavelength time-series observables. By using the StarSim code, now featuring the capability to solve the inverse problem, we analysed $\sim$ 600 days of BVRI multiband photometry from TJO and STELLA telescopes exoplanet host star WASP-52. From the results, we simulated the chromatic contribution of surface phenomena on the observables of its transits. We are able to determine the relevant activity parameters of WASP-52 and reconstruct the time-evolving longitudinal map of active regions. The star shows a heterogeneous surface composed of dark spots with a mean temperature contrast of $575\pm150$ K with filling factors ranging from 3 to 14 %. We studied the chromatic effects on the depths of transits obtained at different epochs with different stellar spot distributions. For WASP-52, with peak-to-peak photometric variations of $\sim$7 % in the visible, the residual effects of dark spots on the measured transit depth, after applying the calculated corrections, are about $10^{-4}$ at 550 nm and $3\times10^{-5}$ at 6$\mu$m. We demonstrate that by using contemporaneous ground-based multiband photometry of an active star, it is possible to reconstruct the parameters and distribution of active regions over time, and thus, quantify the chromatic effects on the planetary radii measured with transit spectroscopy and mitigate them by about an order of magnitude.

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.

AU Microscopii (AU Mic) is the second closest pre main sequence star, at a distance of 9.79 parsecs and with an age of 22 million years. AU Mic possesses a relatively rare and spatially resolved3 edge-on debris disk extending from about 35 to 210 astronomical units from the star, and with clumps exhibiting non-Keplerian motion. Detection of newly formed planets around such a star is challenged by the presence of spots, plage, flares and other manifestations of magnetic activity on the star. Here we report observations of a planet transiting AU Mic. The transiting planet, AU Mic b, has an orbital period of 8.46 days, an orbital distance of 0.07 astronomical units, a radius of 0.4 Jupiter radii, and a mass of less than 0.18 Jupiter masses at 3 sigma confidence. Our observations of a planet co-existing with a debris disk offer the opportunity to test the predictions of current models of planet formation and evolution.

We describe a flexible data reduction package for high resolution cross-dispersed echelle data. This open-source package is developed in Python and includes optional GUIs for most of the steps. It does not require any pre-knowledge about the form or position of the echelle-orders. It has been tested on cross-dispersed echelle spectrographs between 13k and 115k resolution (bifurcated fiber-fed spectrogaph ESO-HARPS and single fiber-fed spectrograph TNT-MRES). HiFLEx can be used to determine radial velocities and is designed to use the TERRA package but can also control the radial velocity packages such as CERES and SERVAL to perform the radial velocity analysis. Tests on HARPS data indicates radial velocities results within 3m/s of the literature pipelines without any fine tuning of extraction parameters.

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

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

The Dispersed Matter Planet Project targets stars with anomalously low Ca II H&K chromospheric emission. High precision, high cadence radial velocity measurements of the F8V star HD 38677 / DMPP-1 reveal four short period planets. DMPP-1 has log(R'HK) = -5.16 which probably indicates the presence of circumstellar absorbing gas arising from an ablating hot planet. The planets have P$_{\rm orb}$ ~ 2.9 - 19 d, i.e., a compact planetary system with super-Earth (~3 M$_\oplus$) to Neptune-mass (~24 M$_\oplus$) planets. These irradiated planets may be chthonian: remnant cores of giant planets after mass-loss while crossing the Neptune desert. Modelling the possible long-term activity indicators while searching for Keplerian signals modifies the recovered planetary signals. A priori inferences about the presence of short period planets allowed the efficient discovery of the DMPP-1 planets. They have great potential for novel and informative follow-up characterisation studies.

We report the first discovery of a thick-disk planet, LHS 1815b (TOI-704b, TIC 260004324), detected in the TESS survey. LHS 1815b transits a bright (V = 12.19 mag, K = 7.99 mag) and quiet M dwarf located $ 29.87\pm0.02 pc$ away with a mass of $0.502\pm0.015 M_{\odot}$ and a radius of $0.501\pm0.030 R_{\odot}$. We validate the planet by combining space and ground-based photometry, spectroscopy, and imaging. The planet has a radius of $1.088\pm 0.064 R_{\oplus}$ with a $3 \sigma$ mass upper-limit of $8.7 M_{\oplus}$. We analyze the galactic kinematics and orbit of the host star LHS1815 and find that it has a large probability ($P_{thick}/P_{thin} = 6482$) to be in the thick disk with a much higher expected maximal height ($Z_{max} = 1.8 kpc$) above the Galactic plane compared with other TESS planet host stars. Future studies of the interior structure and atmospheric properties of planets in such systems using for example the upcoming James Webb Space Telescope (JWST), can investigate the differences in formation efficiency and evolution for planetary systems between different Galactic components (thick and thin disks, and halo).

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.

Some highly irradiated close-in exoplanets orbit stars showing anomalously low stellar chromospheric emission. We attribute this to absorption by circumstellar gas replenished by mass loss from ablating planets. Here we report statistics validating this hypothesis. Among ~3000 nearby, bright, main sequence stars ~40 show depressed chromospheric emission indicative of undiscovered mass-losing planets. The Dispersed Matter Planet Project uses high precision, high cadence radial velocity measurements to detect these planets. We summarise results for two planetary systems (DMPP-1 and DMPP-3) and fully present observations revealing a Mp sin i = 0.469 M$_{\rm J}$ planet in a 5.207 d orbit around the $\gamma$-Doradus pulsator HD 11231 (DMPP-2). We have detected short period planets wherever we have made more than 60 RV measurements, demonstrating that we have originated a very efficient method for detecting nearby compact planetary systems. Our shrouded, ablating planetary systems may be a short-lived phase related to the Neptunian desert: i.e. the dearth of intermediate-mass planets at short orbital periods. The circumstellar gas facilitates compositional analysis; allowing empirical exogeology in the cases of sublimating rocky planets. Dispersed Matter Planet Project discoveries will be important for establishing the empirical mass-radius-composition relationship(s) for low mass planets.

Earth mass exoplanets are difficult to detect. The Dispersed Matter Planet Project (DMPP) identifies stars which are likely to host the most detectable low mass exoplanets. The star DMPP-3 (HD 42936) shows signs of circumstellar absorption, indicative of mass loss from ablating planets. Here we report the radial velocity (RV) discovery of a highly eccentric 507 d binary companion and a hot super-Earth planet in a 6.67 d orbit around the primary star. DMPP-3A is a solar type star while DMPP-3B is just massive enough to fuse hydrogen. The binary, with semi-major axis 1.22 $\pm$ 0.02 AU, is considerably tighter than others known to host planets orbiting only one of the component stars. The configuration of the DMPP-3 planetary system is rare and indicates dynamical interactions, though the evolutionary history is not entirely clear. DMPP-3Ab is possibly the residual core of a giant planet precursor, consistent with the inferred circumstellar gas shroud.

The presence of Jupiter is crucial to the architecture of the Solar System and models underline this to be a generic feature of planetary systems. We find the detection of the difference between the position and motion recorded by the contemporary astrometric satellite Gaia and its precursor Hipparcos can be used to discover Jupiter-like planets. We illustrate how observations of the nearby star $\varepsilon$ Indi A giving astrometric and radial velocity data can be used to independently find the orbit of its suspected companion. The radial velocity and astrometric data provide complementary detections which allow for a much stronger solution than either technique would provide individually. We quantify $\varepsilon$ Indi A b as the closest Jupiter-like exoplanet with a mass of 3 $M_{Jup}$ on a slightly eccentric orbit with an orbital period of 45 yr. While other long-period exoplanets have been discovered, $\varepsilon$ Indi A b provides a well constrained mass and along with the well-studied brown dwarf binary in orbit around $\varepsilon$ Indi A means that the system provides a benchmark case for our understanding of the formation of gas giant planets and brown dwarfs.

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.

The search for Earth-like planets around late-type stars using ultra-stable spectrographs requires a very precise characterization of the stellar activity and the magnetic cycle of the star, since these phenomena induce radial velocity (RV) signals that can be misinterpreted as planetary signals. Among the nearby stars, we have selected Barnard's Star (Gl 699) to carry out a characterization of these phenomena using a set of spectroscopic data that covers about 14.5 years and comes from seven different spectrographs: HARPS, HARPS-N, CARMENES, HIRES, UVES, APF, and PFS; and a set of photometric data that covers about 15.1 years and comes from four different photometric sources: ASAS, FCAPT-RCT, AAVSO, and SNO. We have measured different chromospheric activity indicators (H$\alpha$, Ca~{\sc II}~HK and Na I D), as well as the FWHM of the cross-correlation function computed for a sub-set of the spectroscopic data. The analysis of Generalized Lomb-Scargle periodograms of the time series of different activity indicators reveals that the rotation period of the star is 145 $\pm$ 15 days, consistent with the expected rotation period according to the low activity level of the star and previous claims. The upper limit of the predicted activity-induced RV signal corresponding to this rotation period is about 1 m/s. We also find evidence of a long-term cycle of 10 $\pm$ 2 years that is consistent with previous estimates of magnetic cycles from photometric time series in other M stars of similar activity levels. The available photometric data of the star also support the detection of both the long-term and the rotation signals.

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

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

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

We report Spitzer Space Telescope observations during predicted transits of the exoplanet Proxima Centauri b. As the nearest terrestrial habitable-zone planet we will ever discover, any potential transit of Proxima b would place strong constraints on its radius, bulk density, and atmosphere. Subsequent transmission spectroscopy and secondary-eclipse measurements could then probe the atmospheric chemistry, physical processes, and orbit, including a search for biosignatures. However, our photometric results rule out planetary transits at the 200~ppm level at 4.5$~{\mu}m$, yielding a 3$\sigma$ upper radius limit of 0.4~$R_\rm{\oplus}$ (Earth radii). Previous claims of possible transits from optical ground- and space-based photometry were likely correlated noise in the data from Proxima Centauri's frequent flaring. Follow-up observations should focus on planetary radio emission, phase curves, and direct imaging. Our study indicates dramatically reduced stellar activity at near-to-mid infrared wavelengths, compared to the optical. Proxima b is an ideal target for space-based infrared telescopes, if their instruments can be configured to handle Proxima's brightness.

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

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

A low-amplitude periodic signal in the radial velocity (RV) time series of Barnard's Star was recently attributed to a planetary companion with a minimum mass of ${\sim}$3.2 $M_\oplus$ at an orbital period of $\sim$233 days. The relatively long orbital period and the proximity of Barnard's Star to the Sun raises the question whether the true mass of the planet can be constrained by accurate astrometric measurements. By combining the assumption of an isotropic probability distribution of the orbital orientation with the RV-analysis results, we calculated the probability density function of the astrometric signature of the planet. In addition, we reviewed the astrometric capabilities and limitations of current and upcoming astrometric instruments. We conclude that Gaia and the Hubble Space Telescope (HST) are currently the best-suited instruments to perform the astrometric follow-up observations. Taking the optimistic estimate of their single-epoch accuracy to be $\sim$30 $\mu$as, we find a probability of $\sim$10% to detect the astrometric signature of Barnard's Star b with $\sim$50 individual-epoch observations. In case of no detection, the implied mass upper limit would be $\sim$8 $M_\oplus$, which would place the planet in the super-Earth mass range. In the next decade, observations with the Wide-Field Infrared Space Telescope (WFIRST) may increase the prospects of measuring the true mass of the planet to $\sim$99%.

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

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.

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

AD Leonis is a nearby magnetically active M dwarf. We find Doppler variability with a period of 2.23 days as well as photometric signals: (1) a short period signal which is similar to the radial velocity signal albeit with considerable variability; and (2) a long term activity cycle of 4070$\pm$120 days. We examine the short-term photometric signal in the available ASAS and MOST photometry and find that the signal is not consistently present and varies considerably as a function of time. This signal undergoes a phase change of roughly 0.8 rad when considering the first and second halves of the MOST data set which are separated in median time by 3.38 days. In contrast, the Doppler signal is stable in the combined HARPS and HIRES radial velocities for over 4700 days and does not appear to vary in time in amplitude, phase, period or as a function of extracted wavelength. We consider a variety of star-spot scenarios and find it challenging to simultaneously explain the rapidly varying photometric signal and the stable radial velocity signal as being caused by starspots co-rotating on the stellar surface. This suggests that the origin of the Doppler periodicity might be the gravitational tug of a planet orbiting the star in spin-orbit resonance. For such a scenario and no spin-orbit misalignment, the measured $v \sin i$ indicates an inclination angle of 15.5$\pm$2.5 deg and a planetary companion mass of 0.237$\pm$0.047 M$_{\rm Jup}$.

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

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

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

Proxima Centauri, the star closest to our Sun, is known to host at least one terrestrial planet candidate in a temperate orbit. Here we report the ALMA detection of the star at 1.3 mm wavelength and the discovery of a belt of dust orbiting around it at distances ranging between 1 and 4 au, approximately. Given the low luminosity of the Proxima Centauri star, we estimate a characteristic temperature of about 40 K for this dust, which might constitute the dust component of a small-scale analog to our solar system Kuiper belt. The estimated total mass, including dust and bodies up to 50 km in size, is of the order of 0.01 Earth masses, which is similar to that of the solar Kuiper belt. Our data also show a hint of warmer dust closer to the star. We also find signs of two additional features that might be associated with the Proxima Centauri system, which, however, still require further observations to be confirmed: an outer extremely cold (about 10 K) belt around the star at about 30 au, whose orbital plane is tilted about 45 degrees with respect to the plane of the sky; and additionally, we marginally detect a compact 1.3 mm emission source at a projected distance of about 1.2 arcsec from the star, whose nature is still unknown.

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.

Parallaxes for 331 classical Cepheids, 31 Type II Cepheids and 364 RR Lyrae stars in common between Gaia and the Hipparcos and Tycho-2 catalogues are published in Gaia Data Release 1 (DR1) as part of the Tycho-Gaia Astrometric Solution (TGAS). In order to test these first parallax measurements of the primary standard candles of the cosmological distance ladder, that involve astrometry collected by Gaia during the initial 14 months of science operation, we compared them with literature estimates and derived new period-luminosity ($PL$), period-Wesenheit ($PW$) relations for classical and Type II Cepheids and infrared $PL$, $PL$-metallicity ($PLZ$) and optical luminosity-metallicity ($M_V$-[Fe/H]) relations for the RR Lyrae stars, with zero points based on TGAS. The new relations were computed using multi-band ($V,I,J,K_{\mathrm{s}},W_{1}$) photometry and spectroscopic metal abundances available in the literature, and applying three alternative approaches: (i) by linear least squares fitting the absolute magnitudes inferred from direct transformation of the TGAS parallaxes, (ii) by adopting astrometric-based luminosities, and (iii) using a Bayesian fitting approach. TGAS parallaxes bring a significant added value to the previous Hipparcos estimates. The relations presented in this paper represent first Gaia-calibrated relations and form a "work-in-progress" milestone report in the wait for Gaia-only parallaxes of which a first solution will become available with Gaia's Data Release 2 (DR2) in 2018.

Context. The first Gaia Data Release contains the Tycho-Gaia Astrometric Solution (TGAS). This is a subset of about 2 million stars for which, besides the position and photometry, the proper motion and parallax are calculated using Hipparcos and Tycho-2 positions in 1991.25 as prior information. Aims. We investigate the scientific potential and limitations of the TGAS component by means of the astrometric data for open clusters. Methods. Mean cluster parallax and proper motion values are derived taking into account the error correlations within the astrometric solutions for individual stars, an estimate of the internal velocity dispersion in the cluster, and, where relevant, the effects of the depth of the cluster along the line of sight. Internal consistency of the TGAS data is assessed. Results. Values given for standard uncertainties are still inaccurate and may lead to unrealistic unit-weight standard deviations of least squares solutions for cluster parameters. Reconstructed mean cluster parallax and proper motion values are generally in very good agreement with earlier Hipparcos-based determination, although the Gaia mean parallax for the Pleiades is a significant exception. We have no current explanation for that discrepancy. Most clusters are observed to extend to nearly 15 pc from the cluster centre, and it will be up to future Gaia releases to establish whether those potential cluster-member stars are still dynamically bound to the clusters. Conclusions. The Gaia DR1 provides the means to examine open clusters far beyond their more easily visible cores, and can provide membership assessments based on proper motions and parallaxes. A combined HR diagram shows the same features as observed before using the Hipparcos data, with clearly increased luminosities for older A and F dwarfs.

Radial velocity monitoring has found the signature of a $M \sin i = 1.3$~M$_\oplus$ planet located within the Habitable Zone (HZ) of Proxima Centauri \citep{Anglada16}. Despite a hotter past and an active host star the planet Proxima~b could have retained enough volatiles to sustain surface habitability \citep{Ribas2016}. Here we use a 3D Global Climate Model (GCM) to simulate Proxima b's atmosphere and water cycle for its two likely rotation modes (1:1 and 3:2 spin-orbit resonances) while varying the unconstrained surface water inventory and atmospheric greenhouse effect. We find that a broad range of atmospheric compositions allow surface liquid water. On a tidally-locked planet with sufficient surface water inventory, liquid water is always present, at least in the substellar region. With a non-synchronous rotation, this requires a minimum greenhouse warming ($\sim$10~mbar of CO$_2$ and 1~bar of N$_2$). If the planet is dryer, $\sim$0.5~bar/1.5~bars of CO$_2$ (respectively for asynchronous/synchronous rotation) suffice to prevent the trapping of any arbitrary small water inventory into polar/nightside ice caps. We produce reflection/emission spectra and phase curves for the simulated climates. We find that atmospheric characterization will be possible by direct imaging with forthcoming large telescopes. The angular separation of $7 \lambda/D$ at 1~$\mu$m (with the E-ELT) and a contrast of $\sim$10$^{-7}$ will enable high-resolution spectroscopy and the search for molecular signatures, including H$_2$O, O$_2$, and CO$_2$. The observation of thermal phase curves can be attempted with JWST, thanks to a contrast of $2\times10^{-5}$ at 10~$\mu$m. Proxima~b will also be an exceptional target for future IR interferometers. Within a decade it will be possible to image Proxima~b and possibly determine whether this exoplanet's surface is habitable.

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 detected excess absorption in the emission cores of Ca II H & K during transits of HD 189733b for the first time. Using observations of three transits we investigate the origin of the absorption, which is also seen in H{\alpha} and the Na I D lines. Applying differential spectrophotometry methods to the Ca II H and Ca II K lines combined, using respective passband widths of $\Delta\lambda$ = 0.4 & 0.6 \AA\ yields excess absorption of t$_d$ = 0.0074 $\pm$ 0.0044 (1.7$\sigma$; Transit 1) and 0.0214 +/- 0.0022 (9.8$\sigma$; Transit 2). Similarly, we detect excess H{\alpha} absorption in a passband of width $\Delta\lambda$ = 0.7 \AA, with t$_d$ = 0.0084 $\pm$ 0.0016 (5.2$\sigma$) and 0.0121 $\pm$ 0.0012 (9.9$\sigma$). For both lines, Transit 2 is thus significantly deeper. Combining all three transits for the Na I D lines yields excess absorption of t$_d$ = 0.0041 $\pm$ 0.0006 (6.5$\sigma$). By considering the time series observations of each line, we find that the excess apparent absorption is best recovered in the stellar reference frame. These findings lead us to postulate that the main contribution to the excess transit absorption in the differential light curves arises because the normalising continuum bands form in the photosphere, whereas the line cores contain a chromospheric component. We can not rule out that part of the excess absorption signature arises from the planetary atmosphere, but we present evidence which casts doubt on recent claims to have detected wind motions in the planet's atmosphere in these data.

We report trigonometric parallaxes for 134 low mass stars and brown dwarfs, of which 38 have no previously published measurement and 79 more have improved uncertainties. Our survey targeted nearby targets, so 119 are closer than 30 pc. Of the 38 stars with new parallaxes, 14 are within 20 pc and seven are likely brown dwarfs (spectral types later than L0). These parallaxes are useful for studies of kinematics, multiplicity, and spectrophotometric calibration. Two objects with new parallaxes are confirmed as young stars with membership in nearby young moving groups: LP 870-65 in AB Doradus and G 161-71 in Argus. We also report the first parallax for the planet-hosting star GJ 3470; this allows us to refine the density of its Neptune-mass planet. One T-dwarf, 2MASS J12590470-4336243, previously thought to lie within 4 pc, is found to be at 7.8 pc, and the M-type star 2MASS J01392170-3936088 joins the ranks of nearby stars as it is found to be within 10 pc. Five stars that are over-luminous and/or too red for their spectral types are identified and deserve further study as possible young stars.

Understanding the sources of instrumental systematic noise is a must to improve the design of future spectrographs. In this study, we alternated observations of the well-suited pair of M-stars GJ 725A+B to delve into the sub-night HARPS-N response. Besides the possible presence of a low-mass planet orbiting GJ 725B, our observations reveal changes in the spectral energy distribution (SED) correlated with measurements of the width of the instrumental line profile and, to a lower degree, with the Doppler measurements. To study the origin of these effects, we searched for correlations among several quantities defined and measured on the spectra and on the acquisition images. We find that the changes in apparent SED are very likely related to flux losses at the fibre input. Further tests indicate that such flux losses do not seriously affect the shape of the instrumental point spread function of HARPS-N, but identify an inefficient fitting of the continuum as the most likely source of the systematic variability observed in the FWHM. This index, accounting for the HARPS-N cross-correlation profiles width, is often used to decorrelate Doppler time-series. We show that the Doppler measurement obtained by a parametric least-squares fitting of the spectrum accounting for continuum variability is insensitive to changes in the slope of the SED, suggesting that forward modeling techniques to measure moments of the line profile are the optimal way to achieve higher accuracy. Remaining residual variability at ~1 m/s suggests that for M-stars Doppler surveys the current noise floor still has an instrumental origin.

Given that low-mass stars have intrinsically low luminosities at optical wavelengths and a propensity for stellar activity, it is advantageous for radial velocity (RV) surveys of these objects to use near-infrared (NIR) wavelengths. In this work we describe and test a novel RV extraction pipeline dedicated to retrieving RVs from low mass stars using NIR spectra taken by the CSHELL spectrograph at the NASA Infrared Telescope Facility, where a methane isotopologue gas cell is used for wavelength calibration. The pipeline minimizes the residuals between the observations and a spectral model composed of templates for the target star, the gas cell, and atmospheric telluric absorption; models of the line spread function, continuum curvature, and sinusoidal fringing; and a parameterization of the wavelength solution. The stellar template is derived iteratively from the science observations themselves without a need for separate observations dedicated to retrieving it. Despite limitations from CSHELL's narrow wavelength range and instrumental systematics, we are able to (1) obtain an RV precision of 35 m/s for the RV standard star GJ 15 A over a time baseline of 817 days, reaching the photon noise limit for our attained SNR, (2) achieve ~3 m/s RV precision for the M giant SV Peg over a baseline of several days and confirm its long-term RV trend due to stellar pulsations, as well as obtain nightly noise floors of ~2 - 6 m/s, and (3) show that our data are consistent with the known masses, periods, and orbital eccentricities of the two most massive planets orbiting GJ 876. Future applications of our pipeline to RV surveys using the next generation of NIR spectrographs, such as iSHELL, will enable the potential detection of Super-Earths and Mini-Neptunes in the habitable zones of M dwarfs.

The Kapteyn moving group has been postulated as tidal debris from $\omega$ Centauri. If true, members of the group should show some of the chemical abundance patterns known for stars in the cluster. We present an optical and near-infrared high-resolution, high-S/N spectroscopic study of 14 stars of the Kapteyn group, plus 10 additional stars (the $\omega$ Cen-group) that, while not listed as members of the Kapteyn group as originally defined, have been nevertheless associated dynamically with $\omega$ Centauri. Abundances for Na, O, Mg, Al, Ca and Ba were derived from the optical spectra, while the strength of the chromospheric He I 10830 {\AA} line is studied as a possible helium abundance indicator. The resulting Na-O and Mg-Al patterns for stars of the combined Kapteyn and $\omega$ Cen-group samples do not resemble those of $\omega$ Centauri, and are not different from those of field stars of the Galactic halo. The distribution of equivalent widths of the He I 10830 {\AA} line is consistent with that found among non-active field stars. Therefore, no evidence is found for second-generation stars within our samples, which most likely rules out a globular-cluster origin. Moreover, no hint of the unique Ba-overabundance at the metal-rich end, well-established for $\omega$ Centauri stars, is seen among stars of the combined samples. Because this specific Ba pattern is present in $\omega$ Centauri irrespective of stellar generation, this would rule out the possibility that our entire sample might be composed of only first generation stars from the cluster. Finally, for the stars of the Kapteyn group, the possibility of an origin in the hypothetical $\omega$ Centauri's parent galaxy is disfavored by the different run of $\alpha$-elements with metallicity between our targets and stars from present-day dwarf galaxies.

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.

We perform an analysis of ~80000 photometric measurements for the following 10 stars hosting transiting planets: WASP-2, -4, -5, -52, Kelt-1, CoRoT-2, XO-2, TrES-1, HD 189733, GJ 436. Our analysis includes mainly transit lightcurves from the Exoplanet Transit Database, public photometry from the literature, and some proprietary photometry privately supplied by other authors. Half of these lightcurves were obtained by amateurs. From this photometry we derive 306 transit timing measurements, as well as improved planetary transit parameters. Additionally, for 6 of these 10 stars we present a set of radial velocity measurements obtained from the spectra stored in the HARPS, HARPS-N, and SOPHIE archives using the HARPS-TERRA pipeline. Our analysis of these TTV and RV data did not reveal significant hints of additional orbiting bodies in almost all of the cases. In the WASP-4 case, we found hints of marginally significant TTV signals having amplitude 10-20 sec, although their parameters are model-dependent and uncertain, while radial velocities did not reveal statistically significant Doppler signals.

Robertson et al.(Reports, July 25 2014, p440-444)(1) claimed that activity-induced variability is responsible for the Doppler signal of the proposed planet candidate GJ 581d. We point out that their analysis using periodograms of residual data is incorrect, further promoting inadequate tools. Since the claim challenges the viability of the method to detect exo-Earths, we urge for more appropriate analyses (see appendix).

Characterizing the local space density of double degenerate binary systems is a complementary approach to broad sky surveys of double degenerates to determine the expected rates of white dwarf binary mergers, in particular those that may evolve into other observable phenomena such as extreme helium stars, Am CVn systems, and supernovae Ia. However, there have been few such systems detected in local space. We report here the discovery that WD 1242$-$105, a nearby bright WD, is a double-line spectroscopic binary consisting of two degenerate DA white dwarfs of similar mass and temperature, despite it previously having been spectroscopically characterized as a single degenerate. Follow-up photometry, spectroscopy, and trigonometric parallax have been obtained in an effort to determine the fundamental parameters of each component of this system. The binary has a mass ratio of 0.7 and a trigonometric parallax of 25.5 mas, placing it at a distance of 39 pc. The system's total mass is 0.95 M$_\odot$ and has an orbital period of 2.85 hours, making it the strongest known gravitational wave source ($\log h = -20.78$) in the mHz regime. Because of its orbital period and total mass, WD 1242$-$105 is predicted to merge via gravitational radiation on a timescale of 740 Myr, which will most likely not result in a catastrophic explosion.

This paper presents the simultaneous Bayesian analysis of the radial velocity and activity indices for the K5Ve active star BD+20 1790, that was proposed to host a close-in massive planet by Hern\'an-Obispo et al. (2010). The Bayesian analysis supports the hypothesis of a planetary companion. We present a new orbital solution after removing the two main contributions to the stellar jitter, one that varies with the photometric period and another that varies with the synodic period of the star+planet system. We use a new method to determine these jitter components, considering them as a second and third signals in the system. A discussion on possible magnetic star-planet-interaction is included, based on the Bayesian analysis of the activity indices. We propose two possible sources for flare events of this star: one related to the geometry of the system and the relative movement of the star+planet, and a second one purely stochastic source that is related to the evolution of stellar active regions. Also, we observe for the first time the magnetic field of the star, from spectropolarimetric data.

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.

Due to their higher planet-star mass-ratios, M dwarfs are the easiest targets for detection of low-mass planets orbiting nearby stars using Doppler spectroscopy. Furthermore, because of their low masses and luminosities, Doppler measurements enable the detection of low-mass planets in their habitable zones that correspond to closer orbits than for Solar-type stars. We re-analyse literature UVES radial velocities of 41 nearby M dwarfs in a combination with new velocities obtained from publicly available spectra from the HARPS-ESO spectrograph of these stars in an attempt to constrain any low-amplitude Keplerian signals. We apply Bayesian signal detection criteria, together with posterior sampling techniques, in combination with noise models that take into account correlations in the data and obtain estimates for the number of planet candidates in the sample. More generally, we use the estimated detection probability function to calculate the occurrence rate of low-mass planets around nearby M dwarfs. We report eight new planet candidates in the sample (orbiting GJ 27.1, GJ 160.2, GJ 180, GJ 229, GJ 422, and GJ 682), including two new multiplanet systems, and confirm two previously known candidates in the GJ 433 system based on detections of Keplerian signals in the combined UVES and HARPS radial velocity data that cannot be explained by periodic and/or quasiperiodic phenomena related to stellar activities. Finally, we use the estimated detection probability function to calculate the occurrence rate of low-mass planets around nearby M dwarfs. According to our results, M dwarfs are hosts to an abundance of low-mass planets and the occurrence rate of planets less massive than 10 M$_{\oplus}$ is of the order of one planet per star, possibly even greater. ...

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.

Since low-mass stars have low luminosities, orbits at which liquid water can exist on Earth-sized planets are relatively close-in, which produces Doppler signals that are detectable using state-of-the-art Doppler spectroscopy. GJ 667C is already known to be orbited by two super-Earth candidates. We investigate whether the data supports the presence of additional companions. We obtain new Doppler measurements from HARPS extracted spectra and combined them with those obtained from the PFS and HIRES spectrographs. We used Bayesian and periodogram-based methods to re-assess the number of candidates and evaluated the confidence of each detection. Among other tests, we validated the planet candidates by analyzing correlations of each Doppler signal activity indices and investigate quasi-periodicity. Doppler measurements of GJ 667C are described better by six Keplerian-like signals: the two known candidates (b and c); three additional few-Earth mass candidates with periods of 92, 62 and 39 days (d, e and f); a cold super-Earth in a 260-day orbit (g) and tantalizing evidence of a $\sim$ 1 M$_\oplus$ object in a close-in orbit of 17 days (h). We explore whether long-term stable orbits are compatible with the data by integrating 8$\times 10^4$ solutions derived from the Bayesian samplings. The system consisting of six planets is compatible with dynamically stable configurations. As for the solar system, the most stable solutions do not contain mean-motion resonances and are described well by analytic Laplace-Lagrange solutions. The presence of a seventh planet (h) is supported by the fact that it appears squarely centered on the only island of stability left in the six-planet solution. Habitability assessments accounting for the stellar flux, as well as tidal dissipation effects, indicate that three (maybe four) planets are potentially habitable...

We report two low mass companions orbiting the nearby K7 dwarf GJ 221 that have emerged from re-analyzing 4.4 years of publicly available HARPS spectra complemented with 2 years of high precision Doppler measurements with Magellan/PFS. The HARPS measurements alone contain the clear signal of a low mass companion with a period of 125 days and a minimum mass of 53.2 \mearth (GJ 221b), falling in a mass range where very few planet candidates have been found (sub-Saturn desert). The addition of 17 PFS observations allow the confident detection of a second low mass companion (6.5 \mearth) in a hot orbit (3.87 days period, GJ 221c). Spectrocopic and photometric calibrations suggest that GJ 221 is slightly depleted ([Fe/H]$\sim$ -0.1) compared to the Sun so the presence of two low mass companions in the system confirms the trend that slightly reduced stellar metallicity does not prevent the formation of planets in the super-Earth to sub-Saturn mass regime.

GJ 1214 is orbited by a transiting super-Earth-mass planet. It is a primary target for ongoing efforts to understand the emerging population of super-Earth-mass planets around M dwarfs. We present new precision astrometric measurements, a re-analysis of HARPS radial velocity measurements, and medium-resolution infrared spectroscopy of GJ 1214. We combine these measurements with recent transit follow-up observations and new catalog photometry to provide a comprehensive update of the star-planet properties. The distance is obtained with 0.6% relative uncertainty using CAPScam astrometry. The new value increases the nominal distance to the star by ~10% and is significantly more precise than previous measurements. Updated Doppler measurements combined with published transit observations significantly refine the constraints on the orbital solution. The analysis of the infrared spectrum and photometry confirm that the star is enriched in metals compared to the Sun. Using all this information, combined with empirical mass-luminosity relations for low mass stars, we derive updated values for the bulk properties of the star-planet system. We also use infrared absolute fluxes to estimate the stellar radius and to re-derive the star-planet properties. Both approaches provide very consistent values for the system. Our analysis shows indicates that the favoured mean density of GJ 1214b is 1.6 +/-0.6 g cm^{-3}. We illustrate how fundamental properties of M dwarfs are of paramount importance in the proper characterization of the low mass planetary candidates orbiting them. Given that the distance is now known to better than 1%, interferometric measurements of the stellar radius, additional high precision Doppler observations, and/or or detection of the secondary transit (occultation), are necessary to further improve the constraints on the GJ 1214 star-planet properties.

From common proper motion and signatures of youth, researchers have identified about 30 members of a putative TW Hydrae Association. Only four of these had parallactic distances from Hipparcos. We have measured parallaxes and proper motions for 14 primary members. We combine these with literature values of radial velocities to show that the Galactic space motions of the stars, with the exception of TWA 9 and 22, are parallel and do not indicate convergence at a common formation point sometime in the last few million years. The space motions of TWA 9 and 22 do not agree with the others and indicate that they are not TWA members. The median parallax is 18 mas or 56 pc. We further analyze the stars' absolute magnitudes on pre-main sequence evolutionary tracks and find a range of ages with a median of 10.1 Myr and no correlation between age and Galactic location. The TWA stars may have formed from an extended and filamentary molecular cloud but are not necessarily precisely coeval.

Several M dwarfs are targets of systematical monitoring in searches for Doppler signals caused by low-mass exoplanet companions. As a result, an emerging population of high-multiplicity planetary systems around low-mass stars are being detected as well. We optimize classic data analysis methods and develop new ones to enhance sensitivity towards low-amplitude planets in high-multiplicity systems. We apply these methods to the public HARPS observations of GJ 676A, a nearby M dwarf with one reported gas giant companion. We re-derived Doppler measurements using the template matching method (HARPS-TERRA). We used refined versions of periodograms to assess the presence of additional low-mass companions. We analyse the same dataset with Bayesian tools and compared the performance of both approaches. We confirm the reported massive gas giant candidate and a long period trend, whose curvature is now well detected. We also find very secure evidence of two new candidates in close-in orbits and masses in the super-Earth mass regime. Despite the increased sensitivity of the new periodogram tools, we find that Bayesian methods are more sensitive in the early detection of candidate signals. While hardware development is important, development of data analysis techniques can help to reveal new results from existing data sets with significantly fewer resources. This new system holds the record of minimum-mass range (from Msin i = 4.5 M_Earth to 5 M_Jup) and period range (from 3.6 days to more than 10 years). Although all the planet candidates are substantially more massive, it is the first exoplanetary system with a general architecture similar to our solar system. GJ 676A can be happily added to the family of high-multiplicity planetary systems around M dwarfs.

We present a method to optimize absorption cells for precise wavelength calibration in the near-infrared. We apply it to design and optimize methane isotopologue cells for precision radial velocity measurements in the K band. We also describe the construction and installation of two such cells for the CSHELL spectrograph at NASA's IRTF. We have obtained their high-resolution laboratory spectra, which we can then use in precision radial velocity measurements and which can also have other applications. In terms of obtainable RV precision methane should out-perform other proposed cells, such as the ammonia cell ($^{14}$NH$_{3}$) recently demonstrated on CRIRES/VLT. The laboratory spectra of Ammonia and the Methane cells show strong absorption features in the H band that could also be exploited for precision Doppler measurements. We present spectra and preliminary radial velocity measurements obtained during our first-light run. These initial results show that a precision down to 20-30 m s$^{-1}$ can be obtained using a wavelength interval of only 5 nm in the K band and S/N$\sim$150. This supports the prediction that a precision down to a few m s$^{-1}$ can be achieved on late M dwarfs using the new generation of NIR spectrographs, thus enabling the detection of terrestrial planets in their habitable zones. Doppler measurements in the NIR can also be used to mitigate the radial velocity jitter due to stellar activity enabling more efficient surveys on young active stars.

Doppler spectroscopy has uncovered or confirmed all the known planets orbiting nearby stars. Two main techniques are used to obtain precision Doppler measurements at optical wavelengths. The first approach is the gas cell method, which consists on the least-squares matching of the spectrum of Iodine imprinted on the spectrum of the star. The second method relies on the construction of a stabilized spectrograph externally calibrated in wavelength. The most precise stabilized spectrometer in operation is HARPS, operated by ESO in La Silla Observatory/Chile. The Doppler measurements obtained with HARPS are typically obtained using the Cross-Correlation Function technique (CCF). It consists of multiplying the stellar spectrum with a weighted binary mask and finding the minimum of such product as a function of the Doppler shift. It is known that CCF is suboptimal in exploiting the Doppler information in the stellar spectrum. Here, we describe an algorithm to obtain precision RV measurements using least-squares matching of each observed spectrum to a high signal-to-noise ratio template derived from the same observations. Such algorithm is implemented in our software called HARPS-TERRA (Template Enhanced Radial velocity Re-analysis Application). New radial velocity measurements on a representative sample of stars observed by HARPS is used to illustrate the benefits of the proposed method. We show that, compared to CCF, template matching provides a significant improvement in accuracy, specially when applied to M dwarfs.

We re-analyze 4 years of HARPS spectra of the nearby M1.5 dwarf GJ 667C available through the ESO public archive. The new radial velocity (RV) measurements were obtained using a new data analysis technique that derives the Doppler measurement and other instrumental effects using a least-squares approach. Combining these new 143 measurements with 41 additional RVs from the Magellan/PFS and Keck/HIRES spectrometers, reveals 3 additional signals beyond the previously reported 7.2-day candidate, with periods of 28 days, 75 days, and a secular trend consistent with the presence of a gas giant (Period sim 10 years). The 28-day signal implies a planet candidate with a minimum mass of 4.5 Mearth orbiting well within the canonical definition of the star's liquid water habitable zone, this is, the region around the star at which an Earth-like planet could sustain liquid water on its surface. Still, the ultimate water supporting capability of this candidate depends on properties that are unknown such as its albedo, atmospheric composition and interior dynamics. The 75-day signal is less certain, being significantly affected by aliasing interactions among a potential 91-day signal, and the likely rotation period of the star at 105 days detected in two activity indices. GJ 667C is the common proper motion companion to the GJ 667AB binary, which is metal poor compared to the Sun. The presence of a super-Earth in the habitable zone of a metal poor M dwarf in a triple star system, supports the evidence that such worlds should be ubiquitous in the Galaxy.

We have obtained precision astrometry of the planet hosting M dwarf GJ 317 in the framework of the Carnegie Astrometric Planet Search project. The new astrometric measurements give a distance determination of 15.3 pc, 65% further than previous estimates. The resulting absolute magnitudes suggest it is metal rich and more massive than previously assumed. This result strengthens the correlation between high metallicity and the presence of gas giants around low mass stars. At 15.3 pc, the minimal astrometric amplitude for planet candidate GJ 317b is 0.3 milliarcseconds (edge-on orbit), just below our astrometric sensitivity. However, given the relatively large number of observations and good astrometric precision, a Bayesian Monte Carlo Markov Chain analysis indicates that the mass of planet b has to be smaller than twice the minimum mass with a 99% confidence level, with a most likely value of 2.5 Mjup. Additional RV measurements obtained with Keck by the Lick-Carnegie Planet search program confirm the presence of an additional very long period planet candidate, with a period of 20 years or more. Even though such an object will imprint a large astrometric wobble on the star, its curvature is yet not evident in the astrometry. Given high metallicity, and the trend indicating that multiple systems are rich in low mass companions, this system is likely to host additional low mass planets in its habitable zone that can be readily detected with state-of-the-art optical and near infrared RV measurements.

The radial velocity (RV) method for detecting extrasolar planets has been the most successful to date. The RV signal imprinted by a few Earth-mass planet around a cool star is at the limit of the typical single measurement uncertainty obtained using state-of-the-art spectrographs. This requires relying on statistics in order to unearth signals buried below noise. Artifacts introduced by observing cadences can produce spurious signals or mask genuine planets that should be easily detected otherwise. Here we discuss a particularly confusing statistical degeneracy resulting from the yearly aliasing of the first eccentric harmonic of an already-detected planet. This problem came sharply into focus after the recent announcement of the detection of a 3.1 Earth mass planet candidate in the habitable zone of the nearby low mass star GJ 581. The orbital period of the new candidate planet (GJ 581g) corresponds to an alias of the first eccentric harmonic of a previously reported planet, GJ 581d. Although the star is stable, the combination of the observing cadence and the presence of multiple planets can cause period misinterpretations. In this work, we determine whether the detection of GJ 581g is justified given this degeneracy. We also discuss the implications of our analysis for the recent Bayesian studies of the same data set, which failed to confirm the existence of the new planet. Performing a number of statistical tests, we show that, despite some caveats, the existence of GJ 581g remains the most likely orbital solution to the currently available RV data.

BD+20 1790 is a young active, metal-rich, late-type K5Ve star. We have undertaken a study of stellar activity and kinematics for this star over the past few years. Previous results show a high level of stellar activity, with the presence of prominence-like structures, spots on the surface and strong flare events, despite the moderate rotational velocity of the star. In addition, radial velocity variations with a semi-amplitude of up to 1 km/s were detected. We investigated the nature of these radial velocity variations, in order to determine whether they are due to stellar activity or the reflex motion of the star induced by a companion. We have analysed high-resolution echelle spectra and also two-band photometry was obtained to produce the light curve and determine the photometric period. Based upon the analysis of the bisector velocity span, as well as spectroscopic indices of chromospheric indicators and taking into account the photometric analysis, we report that the best explanation for the RV variation is the presence of a sub-stellar companion. The Keplerian fit of the RV data yields a solution for a close-in massive planet with an orbital period of 7.78 days. The presence of the close-in massive planet could also be an interpretation for the high level of stellar activity detected. Since the RV data are not part of a planet search program, we can consider our results as a serendipitous evidence of a planetary companion. To date, this is the youngest main sequence star for which a planetary candidate has been reported.

Spectroscopic and eclipsing binary systems offer the best means for determining accurate physical properties of stars, including their masses and radii. The data available for low-mass stars have yielded firm evidence that stellar structure models predict smaller radii and higher effective temperatures than observed, but the number of systems with detailed analyses is still small. In this paper we present a complete reanalysis of one of such eclipsing systems, CM Dra, composed of two dM4.5 stars. New and existing light curves as well as a radial velocity curve are modeled to measure the physical properties of both components. The masses and radii determined for the components of CM Dra are M1=0.2310+/-0.0009 Msun, M2=0.2141+/-0.0010 Msun, R1=0.2534+/-0.0019 Rsun, and R2=0.2396+/-0.0015 Rsun. With relative uncertainties well below the 1% level, these values constitute the most accurate properties to date for fully convective stars. This makes CM Dra a valuable benchmark for testing theoretical models. In comparing our measurements with theory, we confirm the discrepancies reported previously for other low-mass eclipsing binaries. These discrepancies seem likely to be due to the effects of magnetic activity. We find that the orbit of this system is slightly eccentric, and we have made use of eclipse timings spanning three decades to infer the apsidal motion and other related properties.