Temperature inversion layers are predicted to be present in ultra-hot giant planet atmospheres. Although such inversion layers have been observed in several ultra-hot Jupiters recently, the chemical species responsible for creating the inversion remain unidentified. Here, we present observation of the thermal emission spectrum of an ultra-hot Jupiter, WASP-189b, at high spectral resolution using the HARPS-N spectrograph. Using the cross-correlation technique, we detect a strong Fe I signal. The detected Fe I spectral lines are found in emission, which is a direct evidence of a temperature inversion in the planetary atmosphere. We further performed a retrieval on the observed spectrum using a forward model with an MCMC approach. When assuming a solar metallicity, the best-fit result returns a temperature of $4320_{-100}^{+120}$ K at the top of the inversion, which is significantly hotter than the planetary equilibrium temperature (2641 K). The temperature at the bottom of the inversion is determined as $2200_{-800}^{+1000}$ K. Such a strong temperature inversion is probably created by the absorption of atomic species like Fe I.

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

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

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

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

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

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

Exoplanets orbiting very close to their host star are strongly irradiated. This can lead the upper atmospheric layers to expand and evaporate into space. The metastable helium (HeI) triplet at 1083.3nm has recently been shown to be a powerful diagnostic to probe extended and escaping exoplanetary atmosphere. We perform high-resolution transmission spectroscopy of the transiting hot Jupiter HD189733b with the GIARPS (GIANO-B + HARPS-N) observing mode of the Telescopio Nazionale Galileo, taking advantage of the simultaneous optical+near infrared spectral coverage to detect HeI in the planet's extended atmosphere and to gauge the impact of stellar magnetic activity on the planetary absorption signal. Observations were performed during five transit events of HD189733b. By comparison of the in- and out-of-transit GIANO-B observations we compute high-resolution transmission spectra, on which we perform equivalent width measurements and light-curves analyses to gauge the excess in-transit absorption in the HeI triplet. We detect an absorption signal during all five transits. The mean in-transit absorption depth amounts to 0.75+/-0.03%. We detect night-to-night variations in the HeI absorption signal likely due to the transit events occurring in presence of stellar surface inhomogeneities. We evaluate the impact of stellar-activity pseudo-signals on the true planetary absorption using a comparative analysis of the HeI and the H$\alpha$ lines. We interpret the time-series of the HeI absorption lines in the three nights not affected by stellar contamination -exhibiting a mean in-transit absorption depth of 0.77+/-0.04%- using a 3-d atmospheric code. Our simulations suggest that the helium layers only fill part of the Roche lobe. Observations can be explained with a thermosphere heated to $\sim$12000 K, expanding up to $\sim$1.2 planetary radii, and losing $\sim$1 g/s of metastable helium.

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

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

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

HD 209458b is an exoplanet with an upper atmosphere undergoing blow-off escape that has mainly been studied using measurements of the Ly-alpha absorption. Recently, high-resolution measurements of absorption in the He I triplet line at 10830 angstroms of several exoplanets (including HD 209458b) have been reported, creating a new opportunity to probe escaping atmospheres. We aim to better understand the atmospheric regions of HD 209458b from where the escape originates. We developed a 1D hydrodynamic model with spherical symmetry for the HD 209458 b thermosphere coupled with a non-local thermodynamic model for the population of the He triplet state. In addition, we performed high-resolution radiative transfer calculations of synthetic spectra for the He triplet lines and compared them with the measured absorption spectrum in order to retrieve information about the atmospheric parameters. We find that the measured spectrum constrains the [H]/[H$^{+}$] transition altitude occurring in the range of 1.2 to 1.9Rp. H is almost fully ionised at altitudes above 2.9Rp. We also find that the X-ray and EUV absorption takes place at effective radii from 1.16 to 1.30Rp, and that the He triplet peak density occurs at altitudes from 1.04 to 1.60Rp. Additionally, the averaged mmw is confined to the 0.61-0.73 g/mole interval, and the thermospheric H/He ratio should be larger than 90/10, and most likely approximately 98/2. We also provide a one-to-one relationship between mass-loss rate and temperature. Based on the energy-limited escape approach and assuming heating efficiencies of 0.1-0.2, we find a mass-loss rate in the range of (0.42-1.00)$\times 10^{11}$ g/s and a corresponding temperature range of 7125 to 8125K. The analysis of the measured He triplet absorption spectrum significantly constrains the thermospheric structure of HD 209458b and advances our knowledge of its escaping atmosphere.

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

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

Mostly multiband photometric transit observations have been used so far to retrieve broadband transmission spectra of transiting exoplanets in order to study their atmosphere. An alternative method has been proposed and has only been used once to recover transmission spectra using chromatic Rossiter-McLaughlin observations. Stellar activity has been shown to potentially imitate narrow and broadband features in the transmission spectra retrieved from multiband photometric observations; however, there has been no study regarding the influence of stellar activity on the retrieved transmission spectra through chromatic Rossiter-McLaughlin. In this study with the modified SOAP3.0 tool, we consider different types of stellar activity features (spots and plages), and we generated a large number of realistic chromatic Rossiter-McLaughlin curves for different types of planets and stars. We were then able to retrieve their transmission spectra to evaluate the impact of stellar activity on them. We find that chromatic Rossiter-McLaughlin observations are also not immune to stellar activity, which can mimic broadband features, such as Rayleigh scattering slope, in their retrieved transmission spectra. We also find that the influence is independent of the planet radius, orbital orientations, orbital period, and stellar rotation rate. However, more general simulations demonstrate that the probability of mimicking strong broadband features is lower than 25% and that can be mitigated by combining several Rossiter-McLaughlin observations obtained during several transits.

Young and active stars generally have regular, almost sinusoidal, patterns of variability attributed to their rotation, while the majority of older and less active stars, including the Sun, have more complex and non-regular light-curves which do not have clear rotational-modulation signals. Consequently, the rotation periods have been successfully determined only for a small fraction of the Sun-like stars observed by transit-based planet-hunting missions, such as CoRoT, Kepler, and TESS. This suggests that only a small fraction of such systems have been properly identified as solar-like analogs. We apply a new method for determining rotation periods of low-activity stars, like the Sun. The method is based on calculating the gradient of the power spectrum (GPS) of stellar brightness variations and identifying a tell-tale inflection point in the spectrum. The rotation frequency is then proportional to the frequency of that inflection point. In this paper test this GPS method against available photometric records of the Sun. We apply GPS, autocorrelation functions, Lomb-Scargle periodograms, and wavelet analyses to the total solar irradiance (TSI) time series obtained from the Total Irradiance Monitor (TIM) on the Solar Radiation and Climate Experiment (SORCE) and the Variability of solar IRradiance and Gravity Oscillations (VIRGO) experiment on the SOlar and Heliospheric Observatory (SoHO) missions. We analyse the performance of all methods at various levels of solar activity. We show that the GPS method returns accurate values of solar rotation independently of the level of solar activity. In particular, it performs well during periods of high solar activity, when TSI variability displays an irregular pattern and other methods fail. Our results suggest that the GPS method can successfully determine the rotational periods of stars with both regular and non-regular light-curves.

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.

We report the discovery of a new planetary system with three transiting planets, one super-Earth and two sub-Neptunes, that orbit EPIC\,249893012, a G8\,IV-V evolved star ($M_\star$\,=\,1.05\,$\pm$\,0.05\,$M_\odot$, $R_\star$\,=\,1.71\,$\pm$\,0.04\,$R_\odot$, $T_\mathrm{eff}$\,=5430\,$\pm$\,85\,K). The star is just leaving the main sequence. We combined \ktwo \ photometry with IRCS adaptive-optics imaging and HARPS, HARPS-N, and CARMENES high-precision radial velocity measurements to confirm the planetary system, determine the stellar parameters, and measure radii, masses, and densities of the three planets. With an orbital period of $3.5949^{+0.0007}_{-0.0007}$ days, a mass of $8.75^{+1.09}_{-1.08}\ M_{\oplus}$ , and a radius of $1.95^{+0.09}_{-0.08}\ R_{\oplus}$, the inner planet b is compatible with nickel-iron core and a silicate mantle ($\rho_b= 6.39^{+1.19}_{-1.04}$ g cm$^{-3}$). Planets c and d with orbital periods of $15.624^{+0.001}_{-0.001}$ and $35.747^{+0.005}_{-0.005}$ days, respectively, have masses and radii of $14.67^{+1,84}_{-1.89}\ M_{\oplus}$ and $3.67^{+0.17}_{-0.14}\ R_{\oplus}$ and $10.18^{+2.46}_{-2.42}\ M_{\oplus}$ and $3.94^{+0.13}_{-0.12}\ R_{\oplus}$, respectively, yielding a mean density of $1.62^{+0.30}_{-0.29}$ and $0.91^{+0.25}_{-0.23}$ g cm$^{-3}$, respectively. The radius of planet b lies in the transition region between rocky and gaseous planets, but its density is consistent with a rocky composition. Its semimajor axis and the corresponding photoevaporation levels to which the planet has been exposed might explain its measured density today. In contrast, the densities and semimajor axes of planets c and d suggest a very thick atmosphere. The singularity of this system, which orbits a slightly evolved star that is just leaving the main sequence, makes it a good candidate for a deeper study from a dynamical point of view.

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

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

Observations of the Sun as a star have been key to guiding models of stellar atmospheres and additionally provide useful insights on the effects of granulation and stellar activity on radial velocity measurements. Most high resolution solar atlases contain telluric lines that span the optical and limit the spectral regions useful for analysis. We present here a telluric-corrected solar atlas covering 0.5-1.0 micron derived from solar spectra taken with a Fourier Transform Spectrograph (FTS) at the Institut f\"ur Astrophysik, G\"ottingen. This atlas is the highest resolution spectrum with a wavelength calibration precise to +/-10 m/s across this 500 nm spectral window. We find that the atlas matches to within 3% of the telluric-corrected Kitt Peak atlas in regions containing telluric absorption weaker than 50% in transmission. The telluric component of the spectral data is fit with a semi-empirical model composed of Lorentz profiles initialized to the HITRAN parameters for each absorption feature. Comparisons between the best-fit telluric parameters describing the Lorentz profile for each absorption feature to the original HITRAN values in general show excellent agreement considering the effects atmospheric pressure and temperature have on our final parameters. However, we identify a small subset of absorption features with larger offsets relative to the catalogued line parameters. We make our final solar atlas available online. We additionally make available the telluric spectra extracted from the data that, given the high resolution of the spectrum, would be useful for studying the time evolution of telluric line shapes and their impact on Doppler measurements.

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

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

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

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

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.

Thanks to the advances in modern instrumentation we have learned about many exoplanets that span a wide range of masses and composition. Studying their atmospheres provides insight into planetary origin, evolution, dynamics, and habitability. Present and future observing facilities will address these important topics in great detail by using more precise observations, high-resolution spectroscopy, and improved analysis methods. We investigate the feasibility of retrieving the vertical temperature distribution and molecular number densities from expected exoplanet spectra in the near-infrared. We use the test case of the CRIRES+, instrument at the Very Large Telescope which will operate in the near-infrared between 1 and 5 micron and resolving powers of R=100000 and R=50000. We also determine the optimal wavelength coverage and observational strategies for increasing accuracy in the retrievals. We used the optimal estimation approach to retrieve the atmospheric parameters from the simulated emission observations of the hot Jupiter HD~189733b. The radiative transfer forward model is calculated using a public version of the tauREx software package. Our simulations show that we can retrieve accurate temperature distribution in a very wide range of atmospheric pressures between 1 bar and $10^{-6}$ bar depending on the chosen spectral region. Retrieving molecular mixing ratios is very challenging, but a simultaneous observations in two separate infrared regions around 1.6 micron and 2.3 micron helps to obtain accurate estimates; the exoplanetary spectra must be of relatively high signal-to-noise ratio S/N>10, while the temperature can already be derived accurately with the lowest value that we considered in this study (S/N=5).

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.

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

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

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

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

High resolution spectroscopy of the lowest-mass stars and brown dwarfs reveals their origins, multiplicity, compositions and physical properties, with implications for the star formation and chemical evolution history of the Milky Way. We motivate the need for high-resolution, infrared spectroscopic surveys to reach these faint sources.

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

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

We present the discovery and characterization of a new transiting planet from Campaign 17 of the Kepler extended mission K2. HD 119130 b is a warm sub-Neptune on a 17-d orbit around a bright (V = 9.9 mag) solar-like G3 V star with a mass and radius of $M_\star = 1.00\pm0.03\,\mathrm{M_\odot}$ and $R_\star = 1.09\pm0.03\,\mathrm{R_\odot}$, respectively. We model simultaneously the K2 photometry and CARMENES spectroscopic data and derive a radius of $R_\mathrm{p} = 2.63_{-0.10}^{+0.12}\,\mathrm{R_\oplus}$ and mass of $M_\mathrm{p} = 24.5_{-4.4}^{+4.4}\,\mathrm{M_\oplus}$, yielding a mean density of $\rho_\mathrm{p} = 7.4_{-1.5}^{+1.6}\,\mathrm{g\,cm^{-3}}$, which makes it one of the densest sub-Neptune planets known to date. We also detect a linear trend in radial velocities of HD 119130 ($\dot{\gamma}_{\rm RV}= -0.40^{+0.07}_{-0.07}\,\mathrm{m\,s^{-1}\,d^{-1}}$) that suggests a long-period companion with a minimum mass on the order of $33\,\mathrm{M_\oplus}$. If confirmed, it would support a formation scenario of HD 119130 b by migration caused by Kozai-Lidov oscillations.

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

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

The behaviour of the large-scale dipolar field, during a star's magnetic cycle, can provide valuable insight into the stellar dynamo and associated magnetic field manifestations such as stellar winds. We investigate the temporal evolution of the dipolar field of the K dwarf 61 Cyg A using spectropolarimetric observations covering nearly one magnetic cycle equivalent to two chromospheric activity cycles. The large-scale magnetic field geometry is reconstructed using Zeeman Doppler imaging, a tomographic inversion technique. Additionally, the chromospheric activity is also monitored. The observations provide an unprecedented sampling of the large-scale field over a single magnetic cycle of a star other than the Sun. Our results show that 61 Cyg A has a dominant dipolar geometry except at chromospheric activity maximum. The dipole axis migrates from the southern to the northern hemisphere during the magnetic cycle. It is located at higher latitudes at chromospheric activity cycle minimum and at middle latitudes during cycle maximum. The dipole is strongest at activity cycle minimum and much weaker at activity cycle maximum. The behaviour of the large-scale dipolar field during the magnetic cycle resembles the solar magnetic cycle. Our results are further confirmation that 61 Cyg A indeed has a large-scale magnetic geometry that is comparable to the Sun's, despite being a slightly older and cooler K dwarf.

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

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

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

Aims. The main goal of this work is to measure rotation periods of the M-type dwarf stars being observed by the CARMENES exoplanet survey to help distinguish radial-velocity signals produced by magnetic activity from those produced by exoplanets. Rotation periods are also fundamental for a detailed study of the relation between activity and rotation in late-type stars. Methods. We look for significant periodic signals in 622 photometric time series of 337 bright, nearby M dwarfs obtained by long-time baseline, automated surveys (MEarth, ASAS, SuperWASP, NSVS, Catalina, ASAS-SN, K2, and HATNet) and for 20 stars which we obtained with four 0.2-0.8 m telescopes at high geographical latitudes. Results. We present 142 rotation periods (73 new) from 0.12 d to 133 d and ten long-term activity cycles (six new) from 3.0 a to 11.5 a. We compare our determinations with those in the existing literature; we investigate the distribution of P rot in the CARMENES input catalogue,the amplitude of photometric variability, and their relation to vsin i and pEW(Halfa); and we identify three very active stars with new rotation periods between 0.34 d and 23.6 d.

One of the most powerful methods used to estimate sky-projected spin-orbit angles of exoplanetary systems is through a spectroscopic transit observation known as the Rossiter-McLaughlin (RM) effect. So far mostly single RM observations have been used to estimate the spin-orbit angle, and thus there have been no studies regarding the variation of estimated spin-orbit angle from transit to transit. Stellar activity can alter the shape of photometric transit light curves and in a similar way they can deform the RM signal. In this paper we discuss several RM observations, obtained using the HARPS spectrograph, of known transiting planets that all transit extremely active stars, and by analyzing them individually we assess the variation in the estimated spin-orbit angle. Our results reveal that the estimated spin-orbit angle can vary significantly (up to 42 degrees) from transit to transit, due to variation in the configuration of stellar active regions over different nights. This finding is almost two times larger than the expected variation predicted from simulations. We could not identify any meaningful correlation between the variation of estimated spin-orbit angles and the stellar magnetic activity indicators. We also investigated two possible approaches to mitigate the stellar activity influence on RM observations. The first strategy was based on obtaining several RM observations and folding them to reduce the stellar activity noise. Our results demonstrated that this is a feasible and robust way to overcome this issue. The second approach is based on acquiring simultaneous high-precision short-cadence photometric transit light curves using TRAPPIST/SPECULOOS telescopes, which provide more information about the stellar active region's properties and allow a better RM modeling.

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

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

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

Atmospheric composition provides essential markers of the most fundamental properties of giant exoplanets, such as their formation mechanism or internal structure. New-generation exoplanet imagers, like VLT/SPHERE or Gemini/GPI, have been designed to achieve very high contrast (>15 mag) at small angular separations ($<$0.5\as) for the detection of young giant planets in the near-infrared, but they only provide very low spectral resolutions ($R<100$) for their characterization. High-dispersion spectroscopy at resolutions up to $10^5$ is one of the most promising pathways for the detailed characterization of exoplanets, but it is currently out of reach for most directly imaged exoplanets because current high-dispersion spectrographs in the near-infrared lack coronagraphs to attenuate the stellar signal and the spatial resolution necessary to resolve the planet. Project HiRISE (High-Resolution Imaging and Spectroscopy of Exoplanets) ambitions to develop a demonstrator that will combine the capabilities of two flagship instruments installed on the ESO Very Large Telescope, the high-contrast exoplanet imager SPHERE and the high-resolution spectrograph CRIRES+, with the goal of answering fundamental questions on the formation, composition and evolution of young planets. In this work, we will present the project, the first set of realistic simulations and the preliminary design of the fiber injection unit that will be implemented in SPHERE.

Thorium hollow cathode lamps (HCLs) are used as frequency calibrators for many high resolution astronomical spectrographs, some of which aim for Doppler precision at the 1 m/s level. We aim to determine the most suitable combination of elements (Th or U, Ar or Ne) for wavelength calibration of astronomical spectrographs, to characterize differences between similar HCLs, and to provide a new U line-list. We record high resolution spectra of different HCLs using a Fourier transform spectrograph: (i) U-Ne, U-Ar, Th-Ne, and Th-Ar lamps in the spectral range from 500 to 1000 nm and U-Ne and U-Ar from 1000 to 1700 nm; (ii) we systematically compare the number of emission lines and the line intensity ratio for a set of 12 U-Ne HCLs; and (iii) we record a master spectrum of U-Ne to create a new U line-list. Uranium lamps show more lines suitable for calibration than Th lamps from 500 to 1000 nm. The filling gas of the lamps significantly affects their performance because Ar and Ne lines contaminate different spectral regions. We find differences (up to 88 %) in the line intensity of U lines in different lamps from the same batch. We find 8239 isolated lines between 500 and 1700 nm that we attribute to U, 3379 of which were not contained in earlier line-lists. The U line-list is available at the http://www.astro.physik.uni-goettingen.de/research/U_atlas . We suggest using a combination of U-Ne and U-Ar lamps to wavelength-calibrate astronomical spectrographs up to 1000 nm. From 1000 to 1700 nm, U-Ne shows better properties. The differences in line strength between different HCLs underline the importance of characterizing HCLs in the laboratory. The new 3379 U lines can significantly improve the radial velocity precision of astronomical spectrographs.

Targeted spectroscopic exoplanet surveys face the challenge of maximizing their planet detection rates by means of careful planning. The number of possible observation combinations for a large exoplanet survey, i.e., the sequence of observations night after night, both in total time and amount of targets, is enormous. Sophisticated scheduling tools and the improved understanding of the exoplanet population are employed to investigate an efficient and optimal way to plan the execution of observations. This is applied to the CARMENES instrument, which is an optical and infrared high-resolution spectrograph that has started a survey of about 300 M-dwarf stars in search for terrestrial exoplanets. We use evolutionary computation techniques to create an automatic scheduler that minimizes the idle periods of the telescope and that distributes the observations among all the targets using configurable criteria. We simulate the case of the CARMENES survey with a realistic sample of targets, and we estimate the efficiency of the planning tool both in terms of telescope operations and planet detection. Our scheduling simulations produce plans that use about 99$\%$ of the available telescope time (including overheads) and optimally distribute the observations among the different targets. Under such conditions, and using current planet statistics, the optimized plan using this tool should allow the CARMENES survey to discover about 65$\%$ of the planets with radial-velocity semi-amplitudes greater than 1$~m\thinspace s^{-1}$ when considering only photon noise. The simulations using our scheduling tool show that it is possible to optimize the survey planning by minimizing idle instrument periods and fulfilling the science objectives in an efficient manner to maximize the scientific return.

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

An interesting question about ultracool dwarfs recently raised in the literature is whether their emission is purely internally driven or partially powered by external processes similar to planetary aurora known from the solar system. In this work we present Hubble Space Telescope observations of the energy fluxes of the M8.5 ultracool dwarf LSR J1835+3259 throughout the UV. The obtained spectra reveal that the object is generally UV-fainter compared to other earlier-type dwarfs. We detect the \ion{Mg}{2} doublet at 2800 \A and constrain an average flux throughout the Near-UV. In the Far-UV without Lyman alpha, the ultracool dwarf is extremely faint with an energy output at least a factor of 1000 smaller as expected from auroral emission physically similar to that on Jupiter. We also detect the red wing of the Lyman alpha emission. Our overall finding is that the observed UV spectrum of LSR J1835+3259 resembles the spectrum of mid/late-type M-dwarf stars relatively well, but it is distinct from a spectrum expected from Jupiter-like auroral processes.

Chromospheric activity monitoring of a wide range of cool stars can provide valuable information on stellar magnetic activity and its dependence on fundamental stellar parameters such as effective temperature and rotation.We compile a chromospheric activity catalogue of 4454 cool stars from a combination of archival HARPS spectra and multiple other surveys, including the Mount Wilson data that have recently been released by the NSO. We explore the variation in chromospheric activity of cool stars along the main sequence for stars with different effective temperatures. Additionally, we also perform an activity-cycle period search and investigate its relation with rotation. The chromospheric activity index, S-index, was measured for 304 main-sequence stars from archived high-resolution HARPS spectra. Additionally, the measured and archived S-indices were converted into the chromospheric flux ratio log R'HK. The activity-cycle periods were determined using the generalised Lomb-Scargle periodogram to study the active and inactive branches on the rotation-activity-cycle period plane. The global sample shows that the bimodality of chromospheric activity, known as the Vaughan-Preston gap, is not prominent, with a significant percentage of the stars at an intermediate-activity level around log R'HK = -4.75. Independently, the cycle period search shows that stars can lie in the region intermediate between the active and inactive branch, which means that the active branch is not as clearly distinct as previously thought. The weakening of the Vaughan-Preston gap indicates that cool stars spin down from a higher activity level and settle at a lower activity level without a sudden break at intermediate activity. Some cycle periods are close to the solar value between the active and inactive branch, which suggests that the solar dynamo is most likely a common case of the stellar dynamo.

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

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

The new CARMENES instrument comprises two high-resolution and high-stability spectrographs that are used to search for habitable planets around M dwarfs in the visible and near-infrared regime via the Doppler technique. Characterising our target sample is important for constraining the physical properties of any planetary systems that are detected. The aim of this paper is to determine the fundamental stellar parameters of the CARMENES M-dwarf target sample from high-resolution spectra observed with CARMENES. We also include several M-dwarf spectra observed with other high-resolution spectrographs, that is CAFE, FEROS, and HRS, for completeness. We used a {chi}^2 method to derive the stellar parameters effective temperature T_eff, surface gravity log g, and metallicity [Fe/H] of the target stars by fitting the most recent version of the PHOENIX-ACES models to high-resolution spectroscopic data. These stellar atmosphere models incorporate a new equation of state to describe spectral features of low-temperature stellar atmospheres. Since T_eff, log g, and [Fe/H] show degeneracies, the surface gravity is determined independently using stellar evolutionary models. We derive the stellar parameters for a total of 300 stars. The fits achieve very good agreement between the PHOENIX models and observed spectra. We estimate that our method provides parameters with uncertainties of {sigma} T_eff = 51 K, {sigma} log g = 0.07, and {sigma} [Fe/H] = 0.16, and show that atmosphere models for low-mass stars have significantly improved in the last years. Our work also provides an independent test of the new PHOENIX-ACES models, and a comparison for other methods using low-resolution spectra. In particular, our effective temperatures agree well with literature values, while metallicities determined with our method exhibit a larger spread when compared to literature results.

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

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

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

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

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

Simultaneous, high-quality measurements of the neutral helium triplet features at 5876~\AA\ and 10830~\AA, respectively, in a sample of solar-type stars are presented. The observations were made with ESO telescopes at the La Silla Paranal Observatory under program ID 088.D-0028(A) and MPG Utility Run for FEROS 088.A-9029(A). The equivalent widths of these features combined with chromospheric models are utilized to infer the fractional area coverage, or filling factor, of magnetic regions outside of spots. We find that the majority of the sample is characterized by filling factors less than unity. However, discrepancies occur among the coolest K-type and warmest and most rapidly rotating F-type dwarf stars. We discuss these apparently anomalous results and find that in the case of K-type stars they are an artifact of the application of chromospheric models best suited to the Sun than to stars with significantly lower $T_\mathrm{eff}$. The case of the F-type rapid rotators can be explained with the measurement uncertainties of the equivalent widths, but they may also be due to a non-magnetic heating component in their atmospheres. With the exceptions noted above, preliminary results suggest that the average heating rates in the active regions are the same from one star to the other, differing in the spatially integrated, observed level of activity due to the area coverage. Hence, differences in activity in this sample are mainly due to the filling factor of active regions.

We propose a dual cavity Fabry-Perot interferometer as a wavelength calibrator and a stability tracking device for astronomical spectrograph. The FPI consists of two adjoining cavities engraved on a low expansion monoblock spacer. A low-finesse astro-cavity is intended for generating a uniform grid of reference lines to calibrate the spectrograph and a high-finesse lock-cavity is meant for tracking the stability of the reference lines using optical frequency standards. The differential length changes in two cavities due to temperature and vibration perturbations are quantitatively analyzed using finite element method. An optimized mounting geometry with fractional length changes $\Delta L/L \approx 1.5\times 10^{-12}$ is suggested. We also identify conditions necessary to suppress relative length variations between two cavities well below 10$^{-10}$~m, thus facilitating accurate dimension tracking and generation of stable reference spectra for Doppler measurement at 10 cms$^{-1}$ level.

Proxima b is a planet with a minimum mass of 1.3 MEarth orbiting within the habitable zone (HZ) of Proxima Centauri, a very low-mass, active star and the Sun's closest neighbor. Here we investigate a number of factors related to the potential habitability of Proxima b and its ability to maintain liquid water on its surface. We set the stage by estimating the current high-energy irradiance of the planet and show that the planet currently receives 30 times more EUV radiation than Earth and 250 times more X-rays. We compute the time evolution of the star's spectrum, which is essential for modeling the flux received over Proxima b's lifetime. We also show that Proxima b's obliquity is likely null and its spin is either synchronous or in a 3:2 spin-orbit resonance, depending on the planet's eccentricity and level of triaxiality. Next we consider the evolution of Proxima b's water inventory. We use our spectral energy distribution to compute the hydrogen loss from the planet with an improved energy-limited escape formalism. Despite the high level of stellar activity we find that Proxima b is likely to have lost less than an Earth ocean's worth of hydrogen before it reached the HZ 100-200 Myr after its formation. The largest uncertainty in our work is the initial water budget, which is not constrained by planet formation models. We conclude that Proxima b is a viable candidate habitable planet.

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.

On March 20, 2015, we obtained 159 spectra of the Sun as a star with the solar telescope and the FTS at the Institut f\"ur Astrophysik G\"ottingen, 76 spectra were taken during partial solar eclipse. We obtained RVs using $I_2$ as wavelength reference and determined the RM curve with a peak-to-peak amplitude of almost 1.4 km s$^{-1}$ at typical RV precision better than 1 m s$^{-1}$. We modeled disk-integrated solar RVs using surface velocities, limb darkening, and information about convective blueshift from 3D magneto-hydrodynamic simulations. We confirm that convective blueshift is crucial to understand solar RVs during eclipse. Our best model reproduced the observations to within a relative precision of 10% with residuals less than 30 m s$^{-1}$. We cross-checked parameterizations of velocity fields using a Dopplergram from the Solar Dynamics Observatory and conclude that disk-integration of the Dopplergram does not provide correct information about convective blueshift necessary for m s$^{-1}$ RV work. As main limitation for modeling RVs during eclipses, we identified limited knowledge about convective blueshift and line shape as functions of solar limb angle. We suspect that our model line profiles are too shallow at limb angles larger than $\mu = 0.6$ resulting in incorrect weighting of the velocities across the solar disk. Alternative explanations cannot be excluded like suppression of convection in magnetic areas and undiscovered systematics during eclipse observations. Accurate observations of solar line profiles across the solar disk are suggested. We publish our RVs taken during solar eclipse as a benchmark curve for codes calculating the RM effect and for models of solar surface velocities and line profiles.

The first generation of E-ELT instruments will include an optical-infrared High Resolution Spectrograph, conventionally indicated as EELT-HIRES, which will be capable of providing unique breakthroughs in the fields of exoplanets, star and planet formation, physics and evolution of stars and galaxies, cosmology and fundamental physics. A 2-year long phase A study for EELT-HIRES has just started and will be performed by a consortium composed of institutes and organisations from Brazil, Chile, Denmark, France, Germany, Italy, Poland, Portugal, Spain, Sweden, Switzerland and United Kingdom. In this paper we describe the science goals and the preliminary technical concept for EELT-HIRES which will be developed during the phase A, as well as its planned development and consortium organisation during the study.

Aims: We search for low-mass companions of M dwarfs and characterize their multiplicity fraction with the purpose of helping in the selection of the most appropriate targets for the CARMENES exoplanet survey. Methods: We obtained high-resolution images in the I band with the lucky imaging instrument FastCam at the 1.5 m Telescopio Carlos Sanchez for 490 mid- to late-M dwarfs. For all the detected binaries, we measured angular separations, position angles, and magnitude differences in the I band. We also calculated the masses of each individual component and estimated orbital periods, using the available magnitude and colour relations for M dwarfs and our own MJ-spectral type and mass-MI relations. To avoid biases in our sample selection, we built a volume-limited sample of M0.0-M5.0 dwarfs that is complete up to 86% within 14 pc. Results: From the 490 observed stars, we detected 80 companions in 76 systems, of which 30 are new discoveries. The multiplicity fraction in our observed sample is 16.7+-2.0% . In our volume-limited sample it is 19.5+-2.3% for angular separations of 0.2 to 5.0 arcsec (1.4-65.6 au), The distribution of the projected physical separations peaks at 2.5-7.5 au. For M0.0-M3.5 V primaries, our search is sensitive to mass ratios higher than 0.3. Binaries with projected physical separations shorter than 50 au tend to be of equal mass. For 26 of our systems, we estimated orbital periods shorter than 50 a, 10 of which are presented here for the first time. We measured variations in angular separation and position angle that are due to orbital motions in 17 of these systems. The contribution of binaries and multiples with angular separations shorter than 0.2 arcsec, longer than 5.0 arcsec, and of spectroscopic binaries identified from previous searches, although not complete, may increase the multiplicity fraction of M dwarfs in our volume-limited sample to at least 36%.

CARMENES (Calar Alto high-Resolution search for M dwarfs with Exoearths with Near-infrared and optical Echelle Spectrographs) started a new planet survey on M-dwarfs in January this year. The new high-resolution spectrographs are operating in the visible and near-infrared at Calar Alto Observatory. They will perform high-accuracy radial-velocity measurements (goal 1 m s-1) of about 300 M-dwarfs with the aim to detect low-mass planets within habitable zones. We characterised the candidate sample for CARMENES and provide fundamental parameters for these stars in order to constrain planetary properties and understand star-planet systems. Using state-of-the-art model atmospheres (PHOENIX-ACES) and chi2-minimization with a downhill-simplex method we determine effective temperature, surface gravity and metallicity [Fe/H] for high-resolution spectra of around 480 stars of spectral types M0.0-6.5V taken with FEROS, CAFE and HRS. We find good agreement between the models and our observed high-resolution spectra. We show the performance of the algorithm, as well as results, parameter and spectral type distributions for the CARMENES candidate sample, which is used to define the CARMENES target sample. We also present first preliminary results obtained from CARMENES spectra.

Several studies have shown that the occultation of stellar active regions by the transiting planet can generate anomalies in the high-precision transit light curves, and these anomalies may lead to an inaccurate estimate of the planetary parameters (e.g., the planet radius). Since the physics and geometry behind the transit light curve and the Rossiter- McLaughlin effect (spectroscopic transit) are the same, the Rossiter-McLaughlin observations are expected to be affected by the occultation of stellar active regions in a similar way. In this paper we perform a fundamental test on the spin-orbit angles as derived by Rossiter-McLaughlin measurements, and we examine the impact of the occultation of stellar active regions by the transiting planet on the spin-orbit angle estimations. Our results show that the inaccurate estimation on the spin-orbit angle due to stellar activity can be quite significant (up to 30 degrees), particularly for the edge-on, aligned, and small transiting planets. Therefore, our results suggest that the aligned transiting planets are the ones that can be easily misinterpreted as misaligned owing to the stellar activity. In other words, the biases introduced by ignoring stellar activity are unlikely to be the culprit for the highly misaligned systems.

The long-term monitoring of magnetic cycles in cool stars is a key diagnostic in understanding how dynamo generation and amplification of magnetic fields occur in stars similar in structure to the Sun. We investigate the temporal evolution of a possible magnetic cycle, determined from its large-scale field, of 61 Cyg A over its activity cycle using spectropolarimetric observations and compare it to the solar case. We use Zeeman Doppler imaging (ZDI) to reconstruct the large-scale magnetic geometry over multiple observational epochs spread over a time span of nine years. We investigate the time evolution of the different components of the large-scale field and compare it with the evolution of its chromospheric activity by measuring the flux in three different chromospheric indicators: Ca II H&K, H-alpha and Ca II infra red triplet lines. We also compare our results with the star's coronal activity using XMM-Newton observations. The large-scale magnetic geometry of 61 Cyg A exhibits polarity reversals in both poloidal and toroidal field components, in phase with the chromospheric activity cycle. We also detect weak solar-like differential rotation with a shear level similar to the Sun. During our observational time span of nine years, 61 Cyg A exhibits solar-like variations in its large-scale field geometry as it evolves from minimum activity to maximum activity and vice versa. During its activity minimum in epoch 2007.59, ZDI reconstructs a simple dipolar geometry which becomes more complex close to activity maximum in epoch 2010.55. The radial field flips polarity and reverts back to a simple geometry in epoch 2013.61. The field is strongly dipolar and the evolution of the dipole component of the field is reminiscent of the solar behaviour. The polarity reversal of the large-scale field indicates a magnetic cycle that is in phase with the chromospheric and coronal cycle.

Close-in extrasolar giant planets are expected to cool their thermospheres by producing H3+ emission in the near-infrared (NIR), but simulations predict H3+ emission intensities that differ in the resulting intensity by several orders of magnitude. We want to test the observability of H3+ emission with CRIRES at the Very Large Telescope (VLT), providing adequate spectral resolution for planetary atmospheric lines in NIR spectra. We search for signatures of planetary H3+ emission in the L` band, using spectra of HD 209458 obtained during and after secondary eclipse of its transiting planet HD 209458 b. We searched for H3+ emission signatures in spectra containing the combined light of the star and, possibly, the planet. With the information on the ephemeris of the transiting planet, we derive the radial velocities at the time of observation and search for the emission at the expected line positions and search for planetary signals and use a shift and add technique combining all observed spectra taken after sec. eclipse to calculate an upper emission limit. We do not find signatures of atmospheric H3+ emission in the spectra containing the combined light of HD 209458 and planet b. We calculate the emission limit for the H3+ line at 3953.0 nm (Q(1, 0)) to be 8.32 E18W and a limit of 5.34E18 W for the line at 3985.5 nm (Q(3, 0)). Comparing our emission limits to the theoretical predictions suggests that we lack 1 to 3 magnitudes of sensitivity to measure H3+ emission in our target object. We show that under more favorable weather conditions the data quality can be improved significantly, reaching 5 E16W for star-planet systems that are close to Earth. We estimate that pushing the detection limit down to 1E15W will be possible with ground-based observations with future instrumentation, for example, the E-ELT.

Radial velocity observations of stars are entering the sub-m/s domain revealing fundamental barriers for Doppler precision experiments. Observations of the Sun as a star can easily overcome the m/s photon limit but face other obstacles. We introduce the G\"ottingen Solar Radial Velocity Project with the goal to obtain high precision (cm/s) radial velocity measurements of the Sun as a star with a Fourier Transform Spectrograph. In this first paper, we present the project and first results. The photon limit of our 2 min observations is at the 2 cm/s level but currently limited by strong instrumental systematics. A drift of a few m/s per h is visible in all observing days probably caused by vignetting of the solar disk in our fiber coupled setup, and imperfections of our guiding system adds further offsets in our data. Binning the data into 30 min groups shows m/s stability after correcting for a daily and linear instrumental trend. Our results show the potential of Sun-as-a-star radial velocity measurements that can possibly be achieved after a substantial upgrade of our spectrograph coupling strategy. Sun-as-a-star observations can provide crucial empirical information about the radial velocity signal of convective motion and stellar activity, and on the wavelength dependence of radial velocity signals caused by stellar line profile variations.

The Second Workshop on Extreme Precision Radial Velocities defined circa 2015 the state of the art Doppler precision and identified the critical path challenges for reaching 10 cm/s measurement precision. The presentations and discussion of key issues for instrumentation and data analysis and the workshop recommendations for achieving this precision are summarized here. Beginning with the HARPS spectrograph, technological advances for precision radial velocity measurements have focused on building extremely stable instruments. To reach still higher precision, future spectrometers will need to produce even higher fidelity spectra. This should be possible with improved environmental control, greater stability in the illumination of the spectrometer optics, better detectors, more precise wavelength calibration, and broader bandwidth spectra. Key data analysis challenges for the precision radial velocity community include distinguishing center of mass Keplerian motion from photospheric velocities, and the proper treatment of telluric contamination. Success here is coupled to the instrument design, but also requires the implementation of robust statistical and modeling techniques. Center of mass velocities produce Doppler shifts that affect every line identically, while photospheric velocities produce line profile asymmetries with wavelength and temporal dependencies that are different from Keplerian signals. Exoplanets are an important subfield of astronomy and there has been an impressive rate of discovery over the past two decades. Higher precision radial velocity measurements are required to serve as a discovery technique for potentially habitable worlds and to characterize detections from transit missions. The future of exoplanet science has very different trajectories depending on the precision that can ultimately be achieved with Doppler measurements.

M-dwarf stars are the most numerous stars in the Universe; they span a wide range in mass and are in the focus of ongoing and planned exoplanet surveys. To investigate and understand their physical nature, detailed spectral information and accurate stellar models are needed. We use a new synthetic atmosphere model generation and compare model spectra to observations. To test the model accuracy, we compared the models to four benchmark stars with atmospheric parameters for which independent information from interferometric radius measurements is available. We used $\chi^2$ -based methods to determine parameters from high-resolution spectroscopic observations. Our synthetic spectra are based on the new PHOENIX grid that uses the ACES description for the equation of state. This is a model generation expected to be especially suitable for the low-temperature atmospheres. We identified suitable spectral tracers of atmospheric parameters and determined the uncertainties in $T_{\rm eff}$, $\log{g}$, and [Fe/H] resulting from degeneracies between parameters and from shortcomings of the model atmospheres. The inherent uncertainties we find are {\sigma}$T_{\rm eff}$= 35 K, {\sigma}$\log{g}$ = 0.14, and {\sigma}[Fe/H] = 0.11. The new model spectra achieve a reliable match to our observed data; our results for $T_{\rm eff}$ and $\log{g}$ are consistent with literature values to within 1{\sigma}. However, metallicities reported from earlier photometric and spectroscopic calibrations in some cases disagree with our results by more than 3 {\sigma}. A possible explanation are systematic errors in earlier metallicity determinations that were based on insufficient descriptions of the cool atmospheres. At this point, however, we cannot definitely identify the reason for this discrepancy, but our analysis indicates that there is a large uncertainty in the accuracy of M-dwarf parameter estimates.

We study rotating thermal convection in spherical shells. We base our analysis on a set of about 450 direct numerical simulations of the (magneto)hydrodynamic equations under the Boussinesq approximation. The Ekman number ranges from $10^{-3}$ to $10^{-5}$. The supercriticality of the convection reaches about 1000 in some models. Four sets of simulations are considered: non-magnetic simulations and dynamo simulations with either free-slip or no-slip flow boundary conditions. The non-magnetic setup with free-slip boundaries generates the strongest zonal flows. Both non-magnetic simulations with no-slip flow boundary conditions and self-consistent dynamos with free-slip boundaries have drastically reduced zonal-flows. Suppression of shear leads to a substantial gain in heat-transfer efficiency, increasing by a factor of 3 in some cases. Such efficiency enhancement occurs as long as the convection is significantly influenced by rotation. At higher convective driving the heat-transfer efficiency tends towards that of the classical non-rotating Rayleigh-B\'enard system. Analysis of the latitudinal distribution of heat flow at the outer boundary reveals that the shear is most effective at suppressing heat-transfer in the equatorial regions. Furthermore, we explore the influence of the magnetic field on the {\em non-zonal} flow components of the convection. For this we compare the heat-transfer efficiency of no-slip non-magnetic cases with that of the no-slip dynamo simulations. We find that at $E=10^{-5}$ magnetic field significantly affects the convection and a maximum gain of about 30\% (as compared to the non-magnetic case) in heat-transfer efficiency is obtained for an Elsasser number of about 3. Our analysis motivates us to speculate that convection in the polar regions in dynamos at $E=10^{-5}$ is probably in a `magnetostrophic' regime.

We present a new solar flux atlas with the aim to understand wavelength precision and accuracy in solar benchmark data. The atlas covers the wavelength range 405--2300 nm and was observed at the Institut f\"ur Astrophysik, G\"ottingen (IAG) with a Fourier Transform Spectrograph. In contrast to other FTS atlases, the entire visible wavelength range was observed simultaneously using only one spectrograph setting. We compare the wavelength solution of the new atlas to the Kitt Peak solar flux atlases and to the HARPS frequency-comb calibrated solar atlas. Comparison reveals systematics in the two Kitt Peak FTS atlases resulting from their wavelength scale construction, and shows consistency between the IAG and the HARPS atlas. We conclude that the IAG atlas is precise and accurate on the order of $\pm 10$ m s$^{-1}$ in the wavelength range 405--1065 nm while the Kitt Peak atlases show deviations as large as several ten to 100 m s$^{-1}$. We determine absolute convective blueshift across the spectrum from the IAG atlas and report slight differences relative to results from the Kitt Peak atlas that we attribute to the differences between wavelength scales. We conclude that benchmark solar data with accurate wavelength solution are crucial to better understand the effect of convection on stellar RV measurements, which is one of the main limitations of Doppler spectroscopy at m s$^{-1}$ precision.

Despite the lack of a shear-rich tachocline region low-mass fully convective stars are capable of generating strong magnetic fields, indicating that a dynamo mechanism fundamentally different from the solar dynamo is at work in these objects. We present a self-consistent three dimensional model of magnetic field generation in low-mass fully convective stars. The model utilizes the anelastic magnetohydrodynamic equations to simulate compressible convection in a rotating sphere. A distributed dynamo working in the model spontaneously produces a dipole-dominated surface magnetic field of the observed strength. The interaction of this field with the turbulent convection in outer layers shreds it, producing small-scale fields that carry most of the magnetic flux. The Zeeman-Doppler-Imaging technique applied to synthetic spectropolarimetric data based on our model recovers most of the large-scale field. Our model simultaneously reproduces the morphology and magnitude of the large-scale field as well as the magnitude of the small-scale field observed on low-mass fully convective stars.

A discharge of nitrogen gas, as created in a microwave-induced plasma, exhibits a very dense molecular emission line spectrum. Emission spectra of this kind could serve as wavelength calibrators for high-resolution astrophysical spectrographs in the near-infrared, where only very few calibration sources are currently available. The compilation of a spectral line list and the characterization of line intensities and line density belong to the initial steps when investigating the feasibility of potential wavelength calibration sources. Although the molecular nitrogen spectrum was extensively studied in the past, to our knowledge, no line list exists that covers a continuous range of several thousand wavenumbers in the near-infrared. We recorded three high-resolution ($\Delta \tilde{\nu} = 0.018$cm$^{-1}$) spectra of a nitrogen gas discharge operated at different microwave powers. The nitrogen gas is kept inside a sealed glass cell at a pressure of 2mbar. The emission lines in the spectra were fitted by a superposition of Gaussian profiles to determine their position, relative intensity, and width. The line parameters were corrected for an absolute wavelength scale, instrumental line broadening, and intensity modulation. Molecular and atomic transitions of nitrogen were identified with available line positions from the literature. We report line lists with more than 40000 emission lines in the spectral range $4500-11000$cm$^{-1}$ ($0.9-2.2$$\mu$m). The spectra exhibit emission lines over the complete spectral range under investigation with about $350-1300$ lines per 100cm$^{-1}$. Depending on the microwave power, a fraction of $35\% - 55\%$ of all lines are blended. The total dynamic range of the detected lines covers about four orders of magnitude.

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

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

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

The convective envelopes of cool main-sequence stars harbour magnetic fields with a complex global and local structure. These fields affect the near-surface convection and the outer stellar atmospheres in many ways and are responsible for the observable magnetic activity of stars. Our aim is to understand the local structure in unipolar regions with moderate average magnetic flux density. These correspond to plage regions covering a substantial fraction of the surface of the Sun (and likely also the surface of other Sun-like stars) during periods of high magnetic activity. We analyse the results of 18 local-box magnetohydrodynamics simulations covering the upper layers of the convection zones and the photospheres of cool main-sequence stars of spectral types F to early M. The average vertical field in these simulations ranges from 20 to 500G. We find a substantial variation of the properties of the surface magnetoconvection between main-sequence stars of different spectral types. As a consequence of a reduced efficiency of the convective collapse of flux tubes, M dwarfs lack bright magnetic structures in unipolar regions of moderate field strength. The spatial correlation between velocity and the magnetic field as well as the lifetime of magnetic structures and their sizes relative to the granules vary significantly along the model sequence of stellar types.

Magnetic fields affect the local structure of the photosphere of stars. They can considerably influence the radiative properties near the optical surface, flow velocities, and the temperature and pressure profiles. We aim at understanding qualitatively the influence of small magnetic flux concentrations in unipolar plage regions on the centre-to-limb variation of the intensity and its contrast and on the shape of spectral line profiles in cool main-sequence stars. We analyse the bolometric and continuum intensity and its angular dependence of 24 radiative magnetohydrodynamic simulations of the near-surface layers of main-sequence stars with six different sets of stellar parameters (spectral types F to early M) and four different average magnetic field strengths (including the non-magnetic case). We also calculated disc-integrated profiles of three spectral lines. The small magnetic flux concentrations formed in the magnetic runs of simulations have a considerable impact on the intensity and its centre-to-limb variation. Spectral lines are not only broadened owing to the Zeeman effect, but are also strongly affected by the modified thermodynamical structure and flow patterns. This indirect magnetic impact on the line profiles is often bigger than that of the Zeeman effect. The effects of the magnetic field on the radiation leaving the star can be considerable and is not restricted to spectral line broadening and polarisation by the Zeeman effect. The inhomogeneous structure of the magnetic field on small length scales and its impact on (and spatial correlation with) the local thermodynamical structure and the flow field near the surface influence the measurement of the global field properties and stellar parameters. These effects need to be taken into account in the interpretation of observations.

Exoplanet science is now in its full expansion, particularly after the CoRoT and Kepler space missions that led us to the discovery of thousands of extra-solar planets. The last decade has taught us that UV observations play a major role in advancing our understanding of planets and of their host stars, but the necessary UV observations can be carried out only by HST, and this is going to be the case for many years to come. It is therefore crucial to build a treasury data archive of UV exoplanet observations formed by a dozen "golden systems" for which observations will be available from the UV to the infrared. Only in this way we will be able to fully exploit JWST observations for exoplanet science, one of the key JWST science case.

Context. CARMENES is a stabilised, high-resolution, double-channel spectrograph at the 3.5 m Calar Alto telescope. It is optimally designed for radial-velocity surveys of M dwarfs with potentially habitable Earth-mass planets. Aims. We prepare a list of the brightest, single M dwarfs in each spectral subtype observable from the northern hemisphere, from which we will select the best planet-hunting targets for CARMENES. Methods. In this first paper on the preparation of our input catalogue, we compiled a large amount of public data and collected low-resolution optical spectroscopy with CAFOS at the 2.2 m Calar Alto telescope for 753 stars. We derived accurate spectral types using a dense grid of standard stars, a double least-squares minimisation technique, and 31 spectral indices previously defined by other authors. Additionally, we quantified surface gravity, metallicity, and chromospheric activity for all the stars in our sample. Results. We calculated spectral types for all 753 stars, of which 305 are new and 448 are revised. We measured pseudo-equivalent widths of Halpha for all the stars in our sample, concluded that chromospheric activity does not affect spectral typing from our indices, and tabulated 49 stars that had been reported to be young stars in open clusters, moving groups, and stellar associations. Of the 753 stars, two are new subdwarf candidates, three are T Tauri stars, 25 are giants, 44 are K dwarfs, and 679 are M dwarfs. Many of the 261 investigated dwarfs in the range M4.0-8.0 V are among the brightest stars known in their spectral subtype. Conclusions. This collection of low-resolution spectroscopic data serves as a candidate target list for the CARMENES survey and can be highly valuable for other radial-velocity surveys of M dwarfs and for studies of cool dwarfs in the solar neighbourhood.

Observations of cool stars reveal dark spot-like features on their surfaces. Compared to sunspots, starspots can be bigger or cover a larger fraction of the stellar surface. While sunspots appear only at low latitudes, starspots are also found in polar regions, in particular on rapidly rotating stars. Sunspots are believed to result from the eruption of magnetic flux-tubes rising from the deep interior of the Sun. The strong magnetic field locally reduces convective heat transport to the solar surface. Such flux-tube models have also been invoked to explain starspot properties. However, these models use several simplifications and so far the generation of either sunspots or starspots has not been demonstrated in a self-consistent simulation of stellar magnetic convection. Here we show that direct numerical simulations of a distributed dynamo operating in a density-stratified rotating spherical shell can spontaneously generate cool spots. Convection in the interior of the model produces a large scale magnetic field which interacts with near surface granular convection leading to strong concentrations of magnetic flux and formation of starspots. Prerequisites for the formation of sizeable high-latitude spots in the model are sufficiently strong density stratification and rapid rotation. Our model presents an alternate mechanism for starspot formation by distributed dynamo action.

State-of-the-art Doppler experiments require wavelength calibration with precision at the cm/s level. A low-finesse Fabry-Perot interferometer (FPI) can provide a wavelength comb with a very large bandwidth as required for astronomical experiments, but unavoidable spectral drifts are difficult to control. Instead of actively controlling the FPI cavity, we propose to passively stabilize the interferometer and track the time-dependent cavity length drift externally. A dual-finesse cavity allows drift tracking during observation. The drift of the cavity length is monitored in the high-finesse range relative to an external standard: a single narrow transmission peak is locked to an external cavity diode laser and compared to an atomic frequency. Following standard locking schemes, tracking at sub-mm/s precision can be achieved. This is several orders of magnitude better than currently planned high-precision Doppler experiments. It allows freedom for relaxed designs rendering this approach particularly interesting for upcoming Doppler experiments. We also show that the large number of interference modes used in an FPI allows us to unambiguously identify the interference mode of each FPI transmission peak defining its absolute wavelength solution. The accuracy reached in each resonance with the laser concept is then defined by the cavity length that is determined from the one locked peak and by the group velocity dispersion. The latter can vary by several 100m/s over the relevant frequency range and severely limits the accuracy of individual peak locations. A potential way to determine the absolute peak positions is to externally measure the frequency of each individual peak with a laser frequency comb (LFC). Thus, the concept of laser-locked FPIs may be useful for applying the absolute accuracy of an LFC to astronomical spectrographs without the need for an LFC at the observatory.

Magnetic activity in Sun-like and low-mass stars causes X-ray coronal emission, which is stronger for more rapidly rotating stars. This relation is often interpreted in terms of the Rossby number, i.e., the ratio of rotation period to convective overturn time. We reconsider this interpretation on the basis of the observed X-ray emission and rotation periods of 821 stars with masses below 1.4 Msun. A generalized analysis of the relation between X-ray luminosity normalized by bolometric luminosity, L_X/L_bol, and combinations of rotational period, P, and stellar radius, R, shows that the Rossby formulation does not provide the solution with minimal scatter. Instead, we find that the relation L_X/L_bol ~ P^{-2}R^{-4} optimally describes the non-saturated fraction of the stars. This relation is equivalent to L_X ~ P^{-2}, indicating that the rotation period alone determines the total X-ray emission. Since L_X is directly related to the magnetic flux at the stellar surface, this means that the surface flux is determined solely by the star's rotation and is independent of other stellar parameters. While a formulation in terms of a Rossby number would be consistent with these results if the convective overturn time scales exactly as L_bol^{-1/2}, our generalized approach emphasizes the need to test a broader range of mechanisms for dynamo action in cool stars.

The solar photospheric magnetic field is highly structured owing to its interaction with the convective flows. Its local structure has a strong influence on the profiles of spectral lines not only by virtue of the Zeeman effect, but also through the modification of the thermodynamical structure (e.g. line weakening in hot small-scale magnetic structures). Many stars harbor surface magnetic fields comparable to or larger than the Sun at solar maximum. Therefore, a strong influence of the field on the surface convection and on spectral line profiles can be expected. We carried out 3D local-box MHD simulations of unipolar magnetized regions (average fields of 20, 100, and 500G) with parameters corresponding to six main-sequence stars (spectral types F3V to M2V). The influence of the magnetic field on the convection and the local thermodynamical structure were analyzed in detail. For three spectral lines, we determined the impact of the magnetic field on the disc-integrated Stokes-I profiles. Line weakening has in many cases a stronger impact on the spectral line profiles than the Zeeman effect. Moreover, for some stars, the correlation between the magnetic field and the vertical velocity strongly influences the line shapes. These effects can impair determinations of stellar magnetic fields since currently used methods neglect the local structure of the magnetic field and its interaction with the convective flows. The MHD simulations presented can be used to quantify these effects and thus help to improve magnetic field measurements of cool main-sequence stars.

CRIRES, the ESO high resolution infrared spectrometer, is a unique instrument which allows astronomers to access a parameter space which up to now was largely uncharted. In its current setup, it consists of a single-order spectrograph providing long-slit, single-order spectroscopy with resolving power up to R=100,000 over a quite narrow spectral range. This has resulted in sub-optimal efficiency and use of telescope time for all the scientific programs requiring broad spectral coverage of compact objects (e.g. chemical abundances of stars and intergalactic medium, search and characterization of extra-solar planets). To overcome these limitations, a consortium was set-up for upgrading CRIRES to a cross-dispersed spectrometer, called CRIRES+. This paper presents the updated optical design of the crossdispersion module for CRIRES+. This new module can be mounted in place of the current pre-disperser unit. The new system yields a factor of >10 increase in simultaneous spectral coverage and maintains a quite long slit (10"), ideal for observations of extended sources and for precise sky-background subtraction.

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.

Based on detailed spectral synthesis we carry out quantitative measurements of the strength and complexity of surface magnetic fields in the four well-known M-dwarfs GJ 388, GJ 729, GJ 285, and GJ 406 populating the mass regime around the boundary between partially and fully convective stars. Very high resolution R=100000, high signal-to-noise (up to 400) near-infrared Stokes I spectra were obtained with CRIRES at ESO's Very Large Telescope covering regions of the FeH Wing-Ford transitions at 1mum. The field distributions in all four stars are characterized by three distinct groups of field components, the data are neither consistent with a smooth distribution of different field strengths, nor with one average field strength covering the full star. We find evidence of a subtle difference in the field distribution of GJ 285 compared to the other three targets. GJ 285 also has the highest average field of 3.5kG and the strongest maximum field component of 7-7.5kG. The maximum local field strengths in our sample seem to be correlated with rotation rate. While the average field strength is saturated, the maximum local field strengths in our sample show no evidence for saturation. We find no difference between the field distributions of partially and fully convective stars. The one star with evidence for a field distribution different to the other three is the most active star (i.e. with largest x-ray luminosity and mean surface magnetic field) rotating relatively fast. A possible explanation is that rotation determines the distribution of surface magnetic fields, and that local field strengths grow with rotation even in stars in which the average field is already saturated.

Context. The study of rotational velocity distributions for normal stars requires an accurate spectral characterization of the objects in order to avoid polluting the results with undetected binary or peculiar stars. This piece of information is a key issue in the understanding of the link between rotation and the presence of chemical peculiarities. Aims. A sample of 47 low v sin i A0-A1 stars (v sin i < 65km/s), initially selected as main-sequence normal stars, are investigated with high-resolution and high signal-to-noise spectroscopic data. The aim is to detect spectroscopic binaries and chemically peculiar stars, and eventually establish a list of confirmed normal stars. Methods. A detailed abundance analysis and spectral synthesis is performed to derive abundances for 14 chemical species. A hierarchical classification, taking measurement errors into account, is applied to the abundance space and splits the sample into two different groups, identified as the chemically peculiar stars and the normal stars. Results. We show that about one third of the sample is actually composed of spectroscopic binaries (12 double-lined and five single-lined spectroscopic binaries). The hierarchical classification breaks down the remaining sample into 13 chemically peculiar stars (or uncertain) and 17 normal stars.

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.

Theoretical predictions suggest that the distribution of planets in very young stars could be very different to that typically observed in Gyr old systems that are the current focus of radial velocity surveys. However, the detection of planets around young stars is hampered by the increased stellar activity associated with young stars, the signatures of which can bias the detection of planets. In this paper we place realistic limitations on the possibilities for detecting planets around young active G and K dwarfs. The models of stellar activity based on tomographic imaging of the G dwarf HD 141943 and the K1 dwarf AB Dor and also include contributions from plage and many small random starspots. Our results show that the increased stellar activity levels present on young Solar-type stars strongly impacts the detection of Earth-mass and Jupiter mass planets and that the degree of activity jitter is directly correlated with stellar \vsinis. We also show that for G and K dwarfs, the distribution of activity in individual stars is more important than the differences in induced radial velocities as a function of spectral type. We conclude that Jupiter mass planets can be detected close-in around fast-rotating young active stars, Neptune-mass planets around moderate rotators and that Super-Earths are only detectable around very slowly rotating stars. The effects of an increase in stellar activity jitter by observing younger stars can be compensated for by extending the observational base-line to at least 100 epochs.

Stellar magnetic field measurements obtained from spectropolarimetry offer key data for activity and dynamo studies, and we present the results of a major high-resolution spectropolarimetric Bcool project magnetic snapshot survey of 170 solar-type stars from observations with the Telescope Bernard Lyot and the Canada-France-Hawaii Telescope. For each target star a high signal-to-noise circularly polarised Stokes V profile has been obtained using Least-Squares Deconvolution, and used to detect surface magnetic fields and measure the corresponding mean surface longitudinal magnetic field ($B_{l}$). Chromospheric activity indicators were also measured. Surface magnetic fields were detected for 67 stars, with 21 of these stars classified as mature solar-type stars, a result that increases by a factor of four the number of mature solar-type stars on which magnetic fields have been observed. In addition, a magnetic field was detected for 3 out of 18 of the subgiant stars surveyed. For the population of K-dwarfs the mean value of $B_{l}$ ($|B_{l}|_{mean}$) was also found to be higher (5.7 G) than $|B_{l}|_{mean}$ measured for the G-dwarfs (3.2 G) and the F-dwarfs (3.3 G). For the sample as a whole $|B_{l}|_{mean}$ increases with rotation rate and decreases with age, and the upper envelope for $|B_{l}|$ correlates well with the observed chromospheric emission. Stars with a chromospheric S-index greater than about 0.2 show a high magnetic field detection rate and so offer optimal targets for future studies. This survey constitutes the most extensive spectropolarimetric survey of cool stars undertaken to date, and suggests that it is feasible to pursue magnetic mapping of a wide range of moderately active solar-type stars to improve understanding of their surface fields and dynamos.

Stellar differential rotation can be separated into two main regimes: solar-like when the equator rotates faster than the poles and anti-solar when the polar regions rotate faster than the equator. We investigate the transition between these two regimes with 3-D numerical simulations of rotating spherical shells. We conduct a systematic parameter study which also includes models from different research groups. We find that the direction of the differential rotation is governed by the contribution of the Coriolis force in the force balance, independently of the model setup (presence of a magnetic field, thickness of the convective layer, density stratification). Rapidly-rotating cases with a small Rossby number yield solar-like differential rotation, while weakly-rotating models sustain anti-solar differential rotation. Close to the transition, the two kinds of differential rotation are two possible bistable states. This study provides theoretical support for the existence of anti-solar differential rotation in cool stars with large Rossby numbers.

Stellar magnetic dynamos are driven by rotation, rapidly rotating stars produce stronger magnetic fields than slowly rotating stars do. The Zeeman effect is the most important indicator of magnetic fields, but Zeeman broadening must be disentangled from other broadening mechanisms, mainly rotation. The relations between rotation and magnetic field generation, between Doppler and Zeeman line broadening, and between rotation, stellar radius, and angular momentum evolution introduce several observational biases that affect our picture of stellar magnetism. In this overview, a few of these relations are explicitly shown, and the currently known distribution of field measurements is presented.

Building on the experience of the high-resolution community with the suite of VLT high-resolution spectrographs, which has been tremendously successful, we outline here the (science) case for a high-fidelity, high-resolution spectrograph with wide wavelength coverage at the E-ELT. Flagship science drivers include: the study of exo-planetary atmospheres with the prospect of the detection of signatures of life on rocky planets; the chemical composition of planetary debris on the surface of white dwarfs; the spectroscopic study of protoplanetary and proto-stellar disks; the extension of Galactic archaeology to the Local Group and beyond; spectroscopic studies of the evolution of galaxies with samples that, unlike now, are no longer restricted to strongly star forming and/or very massive galaxies; the unraveling of the complex roles of stellar and AGN feedback; the study of the chemical signatures imprinted by population III stars on the IGM during the epoch of reionization; the exciting possibility of paradigm-changing contributions to fundamental physics. The requirements of these science cases can be met by a stable instrument with a spectral resolution of R~100,000 and broad, simultaneous spectral coverage extending from 370nm to 2500nm. Most science cases do not require spatially resolved information, and can be pursued in seeing-limited mode, although some of them would benefit by the E-ELT diffraction limited resolution. Some multiplexing would also be beneficial for some of the science cases. (Abridged)

Observations of active M dwarfs show a broad variety of large-scale magnetic fields encompassing dipole-dominated and multipolar geometries. We detail the analogy between some anelastic dynamo simulations and spectropolarimetric observations of 23 M stars. In numerical models, the relative contribution of inertia and Coriolis force in the global force balance -estimated by the so-called local Rossby number- is known to have a strong impact on the magnetic field geometry. We discuss the relevance of this parameter in setting the large-scale magnetic field of M dwarfs.

The near-surface layers of cool main-sequence stars are structured by convective flows, which are overshooting into the atmosphere. The flows and the associated spatio-temporal variations of density and temperature affect spectral line profiles and thus have an impact on estimates of stellar properties such as effective temperature, gravitational acceleration, and abundances. We aim at identifying distinctive properties of the thermodynamic structure of the atmospheres of different stars and understand their causes. We ran comprehensive 3D radiation hydrodynamics simulations of the near-surface layers of six simulated stars of spectral type F3V to M2V with the MURaM code. We carry out a systematic parameter study of the mean stratifications, flow structures, and the energy flux in these stars.\par Results: We find monotonic trends along the lower main sequence in granule size, flow velocity, and intensity contrast. The convection in the M-star models differs substantially from that of the hotter stars, mainly owing to the more gradual transition from convective to radiative energy transport. While the basic mechanisms driving surface convection in cool stars are the same, the properties of the convection vary along the main sequence. Apart from monotonic trends in rms velocity, intensity contrast, granule size, etc., there is a transition between "naked" and "hidden" granulation around spectral type K5V caused by the (highly non-linear) temperature dependence of the opacity. These variations have to be taken into account when stellar parameters are derived from spectra.

The atmospheres of cool main-sequence stars are structured by convective flows from the convective envelope that penetrate the optically thin layers and lead to structuring of the stellar atmospheres analogous to solar granulation. The flows have considerable influence on the 3D structure of temperature and pressure and affect the profiles of spectral lines formed in the photosphere. For the set of six 3D radiative (M)HD simulations of cool main-sequence stars described in the first paper of this series, we analyse the near-surface layers. We aim at describing the properties of granulation of different stars and at quantifying the effects on spectral lines of the thermodynamic structure and flows of 3D convective atmospheres. We detected and tracked granules in brightness images from the simulations to analyse their statistical properties, as well as their evolution and lifetime. We calculated spatially resolved spectral line profiles using the line synthesis code SPINOR. To enable a comparison to stellar observations, we implemented a numerical disc-integration, which includes (differential) rotation. Although the stellar parameters change considerably along the model sequence, the properties of the granules are very similar. The impact of the 3D structure of the atmospheres on line profiles is measurable in disc-integrated spectra. Line asymmetries caused by convection are modulated by stellar rotation. The 3D structure of cool stellar atmospheres as shaped by convective flows has to be taken into account when using photospheric lines to determine stellar parameters.

We present rotation periods for thousands of active stars in the Kepler field derived from Q3 data. In most cases a second period close to the rotation period was detected, which we interpreted as surface differential rotation (DR). Active stars were selected from the whole sample using the range of the variability amplitude. To detect different periods in the light curves we used the Lomb-Scargle periodogram in a pre-whitening approach to achieve parameters for a global sine fit. The most dominant periods from the fit were ascribed to different surface rotation periods, but spot evolution could also play a role. Due to the large number of stars the period errors were estimated in a statistical way. We thus cannot exclude the existence of false positives among our periods. In our sample of 40.661 active stars we found 24.124 rotation periods $P_1$ between 0.5-45 days. The distribution of stars with 0.5 < B-V < 1.0 and ages derived from angular momentum evolution that are younger than 300 Myr is consistent with a constant star-formation rate. A second period $P_2$ within $\pm30$% of the rotation period $P_1$ was found in 18.619 stars (77.2%). Attributing these two periods to DR we found that the relative shear $\alpha=\Delta\Omega/\Omega$ increases with rotation period, and slightly decreases with effective temperature. The absolute shear $\Delta\Omega$ slightly increases between $T_{eff}=3500-6000$ K. Above 6000 K $\Delta\Omega$ shows much larger scatter. We found weak dependence of $\Delta\Omega$ on rotation period. Latitudinal differential rotation measured for the first time in more than 18.000 stars provides a comprehensive picture of stellar surface shear, consistent with major predictions from mean-field theory. To what extent our observations are prone to false positives and selection bias is not fully explored, and needs to be addressed using more Kepler data.

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

Co-rotating spots at different latitudes on the stellar surface generate periodic photometric variability and can be useful proxies to detect Differential Rotation (DR). DR is a major ingredient of the solar dynamo but observations of stellar DR are rather sparse. In view of the Kepler space telescope we are interested in the detection of DR using photometric information of the star, and to develop a fast method to determine stellar DR from photometric data. We ran a large Monte-Carlo simulation of differentially rotating spotted stars with very different properties to investigate the detectability of DR. For different noise levels the resulting light curves are prewhitened using Lomb-Scargle periodograms to derive parameters for a global sine fit to detect periodicities. We show under what conditions DR can successfully be detected from photometric data, and in which cases the light curve provides insufficient or even misleading information on the stellar rotation law. In our simulations, the most significant period P1_{out} could be detected in 96.2% of all light curves. Detection of a second period close to P1_{out} is the signature of DR in our model. For the noise-free case, in 64.2% of all stars such a period was found. Calculating the measured latitudinal shear of two distinct spots \alpha_{out}, and comparing it to the known original spot rotation rates shows that the real value is on average 3.2% lower. Comparing the total equator-to-pole shear $\alpha$ to $\alpha_{out}$ we find that $\alpha$ is underestimated by 8.8%, esp. the detection of DR for stars with $\alpha$ < 6% is challenging. Finally, we apply our method to four differentially rotating Kepler stars and find close agreement with results from detailed modeling. Our method is capable of measuring stellar rotation periods and detecting DR with relatively high accuracy and is suitable for large data sets.

Transit-spectroscopy of exoplanets is one of the key observational techniques to characterize the extrasolar planet and its atmosphere. The observational challenges of these measurements require dedicated instrumentation and only the space environment allows an undisturbed access to earth-like atmospheric features such as water or carbon-dioxide. Therefore, several exoplanet-specific space missions are currently being studied. One of them is EChO, the Exoplanet Characterization Observatory, which is part of ESA's Cosmic Vision 2015-2025 program, and which is one of four candidates for the M3 launch slot in 2024. In this paper we present the results of our assessment study of the EChO spectrometer, the only science instrument onboard this spacecraft. The instrument is a multi-channel all-reflective dispersive spectrometer, covering the wavelength range from 400 nm to 16 microns simultaneously with a moderately low spectral resolution. We illustrate how the key technical challenge of the EChO mission - the high photometric stability - influences the choice of spectrometer concept and drives fundamentally the instrument design. First performance evaluations underline the fitness of the elaborated design solution for the needs of the EChO mission.

We present a new library of high-resolution synthetic spectra based on the stellar atmosphere code PHOENIX that can be used for a wide range of applications of spectral analysis and stellar parameter synthesis. The spherical mode of PHOENIX was used to create model atmospheres and to derive detailed synthetic stellar spectra from them. We present a new self-consistent way of describing micro-turbulence for our model atmospheres. The synthetic spectra cover the wavelength range from 500AA to 50.000AA with resolutions of R=500.000 in the optical and near IR, R=100.000 in the IR and a step size of 0.1AA in the UV. The parameter space covers 2.300K<=Teff<=12.000K, 0.0<=log(g)<=+6.0, -4.0<=[Fe/H]<=+1.0, and -0.2<=[alpha/Fe]<=+1.2. The library is a work in progress and we expect to extend it up to Teff=25.000 K.

We utilize Kepler data to study the precision differential photometric variability of solar-type and cooler stars at different timescales, ranging from half an hour to 3 months. We define a diagnostic that characterizes the median differential intensity change between data bins of a given timescale. We apply the same diagnostics to SOHO data that has been rendered comparable to Kepler. The Sun exhibits similar photometric variability on all timescales as comparable solar-type stars in the Kepler field (it is not unusually quiet). The previously-defined photometric "range" serves as our activity proxy (driven by starspot coverage). We revisit the fraction of comparable stars in the Kepler field that are more active than the Sun. The exact active fraction depends on what is meant by "more active than the Sun", and on the magnitude limit of the sample of stars considered. This active fraction is between a quarter and a third (depending on the timescale). We argue that a reliable result requires timescales of half a day or longer and stars brighter than Kepler magnitude of 14, otherwise non-stellar noise distorts it. We also analyze main sequence stars grouped by temperature from 6500-3500K. As one moves to cooler stars, the active fraction of stars becomes steadily larger (greater than 90% for early M dwarfs). The Sun is a good photometric model at all timescales for those cooler stars that have long-term variability within the span of solar variability.

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

Context: observations of rapidly rotating M dwarfs show a broad variety of large-scale magnetic fields encompassing dipole-dominated and multipolar geometries. In dynamo models, the relative importance of inertia in the force balance -- quantified by the local Rossby number -- is known to have a strong impact on the magnetic field geometry. Aims: we aim to assess the relevance of the local Rossby number in controlling the large-scale magnetic field geometry of M dwarfs. Methods: we explore the similarities between anelastic dynamo models in spherical shells and observations of active M-dwarfs, focusing on field geometries derived from spectropolarimetric studies. To do so, we construct observation-based quantities aimed to reflect the diagnostic parameters employed in numerical models. Results: the transition between dipole-dominated and multipolar large-scale fields in early to mid M dwarfs is tentatively attributed to a Rossby number threshold. We interpret late M dwarfs magnetism to result from a dynamo bistability occurring at low Rossby number. By analogy with numerical models, we expect different amplitudes of differential rotation on the two dynamo branches.

The K2.5 dwarf HD 40307 has been reported to host three super-Earths. The system lacks massive planets and is therefore a potential candidate for having additional low-mass planetary companions. We re-derive Doppler measurements from public HARPS spectra of HD 40307 to confirm the significance of the reported signals using independent data analysis methods. We also investigate these measurements for additional low-amplitude signals. We used Bayesian analysis of our radial velocities to estimate the probability densities of different model parameters. We also estimated the relative probabilities of models with differing numbers of Keplerian signals and verified their significance using periodogram analyses. We investigated the relation of the detected signals with the chromospheric emission of the star. As previously reported for other objects, we found that radial velocity signals correlated with the S-index are strongly wavelength dependent. We identify two additional clear signals with periods of 34 and 51 days, both corresponding to planet candidates with minimum masses a few times that of the Earth. An additional sixth candidate is initially found at a period of 320 days. However, this signal correlates strongly with the chromospheric emission from the star and is also strongly wavelength dependent. When analysing the red half of the spectra only, the five putative planetary signals are recovered together with a very significant periodicity at about 200 days. This signal has a similar amplitude as the other new signals reported in the current work and corresponds to a planet candidate with M sin i = 7 Me (HD 40307 g). ...

CARMENES (Calar Alto high-Resolution search for M dwarfs with Exo-earths with Near-infrared and optical Echelle Spectrographs) is a next generation instrument being built for the 3.5-m telescope at the Calar Alto Observatory by a consortium of eleven Spanish and German institutions. Conducting a five-year exoplanet survey targeting 300 M dwarfs with the completed instrument is an integral part of the project. The CARMENES instrument consists of two separate echelle spectrographs covering the wavelength range from 550 to 1700 nm at a spectral resolution of R=82,000, fed by fibers from the Cassegrain focus of the telescope. The spectrographs are housed in vacuum tanks providing the temperature-stabilized environments necessary to enable a 1 m/s radial velocity precision employing a simultaneous calibration with emission-line lamps.

Magnetic fields are regarded as a crucial element for our understanding of stellar physics. They can be studied with a variety of methods which provide complementary - and sometimes contradictory - information about the structure, strength and dynamics of the magnetic field and its role in the evolution of stars. Stellar magnetic fields can be investigated either with direct methods based on the Zeeman effect or through the observation of activity phenomena resulting from the interaction of the field with the stellar atmosphere. In this Cool Stars XVII Splinter Session we discussed the results obtained by the many ongoing studies of stellar activity and direct studies of surface magnetic fields, as well as the state- of-the-art techniques on which they are based. We show the strengths and limitations of the various approaches currently used and to point out their evolution as well as the interest of coupling various magnetism and activity proxies.

We present projected rotational velocities and new measurements of the rotational profile of some 180 nearby stars with spectral types A-F. The overall broadening profile is derived analysing spectral line shape from hundreds of spectral lines by the method of least-squares deconvolution. Rigid and differential rotation can be distinguished in 56 cases. Ten stars with significant differential rotation rates are identified. As of now, 33 differential rotators detected by line profile analysis have been confirmed. The frequency of differential rotators decreases towards high effective temperature and rapid rotation. There is evidence for two populations of differential rotators with a gap in between at spectral type early-F. The gap can only partly be explained by an upper bound found for the horizontal shear of F stars. Apparently, the physical conditions of differential rotation change at early-F spectral types.

Cool stars like the Sun harbor convection zones capable of producing substantial surface magnetic fields leading to stellar magnetic activity. The influence of stellar parameters like rotation, radius, and age on cool-star magnetism, and the importance of the shear layer between a radiative core and the convective envelope for the generation of magnetic fields are keys for our understanding of low-mass stellar dynamos, the solar dynamo, and also for other large-scale and planetary dynamos. Our observational picture of cool-star magnetic fields has improved tremendously over the last years. Sophisticated methods were developed to search for the subtle effects of magnetism, which are difficult to detect particularly in cool stars. With an emphasis on the assumptions and capabilities of modern methods used to measure magnetism in cool stars, I review the different techniques available for magnetic field measurements. I collect the analyses on cool-star magnetic fields and try to compare results from different methods, and I review empirical evidence that led to our current picture of magnetic fields and their generation in cool stars and brown dwarfs.

We present a catalogue of rotation and chromospheric activity in a sample of 334 M dwarfs of spectral types M0--M4.5 populating the parameter space around the boundary to full convection. We obtained high-resolution optical spectra for 206 targets and determined projected rotational velocity, vsini, and Halpha emission. The data are combined with measurements of vsini in field stars of the same spectral type from the literature. Our sample adds 157 new rotation measurements to the existing literature and almost doubles the sample of available vsini. The final sample provides a statistically meaningful picture of rotation and activity at the transition to full convection in the solar neighborhood. We confirm the steep rise in the fraction of active stars at the transition to full convection known from earlier work. In addition, we see a clear rise in rotational velocity in the same stars. In very few stars, no chromospheric activity but a detection of rotational broadening was reported. We argue that all of them are probably spurious detections; we conclude that in our sample all significantly rotating stars are active, and all active stars are significantly rotating. The rotation-activity relation is valid in partially and in fully convective stars. Thus, we do not observe any evidence for a transition from a rotationally dominated dynamo in partially convective stars to a rotation-independent turbulent dynamo in fully convective stars; turbulent dynamos in fully convective stars of spectral types around M4 are still driven by rotation. Finally, we compare projected rotational velocities of 33 stars to rotational periods derived from photometry in the literature and determine inclinations for a few of them.

A dedicated mission to investigate exoplanetary atmospheres represents a major milestone in our quest to understand our place in the universe by placing our Solar System in context and by addressing the suitability of planets for the presence of life. EChO -the Exoplanet Characterisation Observatory- is a mission concept specifically geared for this purpose. EChO will provide simultaneous, multi-wavelength spectroscopic observations on a stable platform that will allow very long exposures. EChO will build on observations by Hubble, Spitzer and groundbased telescopes, which discovered the first molecules and atoms in exoplanetary atmospheres. EChO will simultaneously observe a broad enough spectral region -from the visible to the mid-IR- to constrain from one single spectrum the temperature structure of the atmosphere and the abundances of the major molecular species. The spectral range and resolution are tailored to separate bands belonging to up to 30 molecules to retrieve the composition and temperature structure of planetary atmospheres. The target list for EChO includes planets ranging from Jupiter-sized with equilibrium temperatures Teq up to 2000 K, to those of a few Earth masses, with Teq ~300 K. We have baselined a dispersive spectrograph design covering continuously the 0.4-16 micron spectral range in 6 channels (1 in the VIS, 5 in the IR), which allows the spectral resolution to be adapted from several tens to several hundreds, depending on the target brightness. The instrument will be mounted behind a 1.5 m class telescope, passively cooled to 50 K, with the instrument structure and optics passively cooled to ~45 K. EChO will be placed in a grand halo orbit around L2. We have also undertaken a first-order cost and development plan analysis and find that EChO is easily compatible with the ESA M-class mission framework.

The relativley large spread in the derived metallicities ([Fe/H]) of M dwarfs shows that various approaches have not yet converged to consistency. The presence of strong molecular features, and incomplete line lists for the corresponding molecules have made metallicity determinations of M dwarfs difficult. Furthermore, the faint M dwarfs require long exposure times for a signal-to-noise ratio sufficient for a detailed spectroscopic abundance analysis. We present a high-resolution (R~50,000) spectroscopic study of a sample of eight single M dwarfs and three wide-binary systems observed in the infrared J-band. The absence of large molecular contributions allow for a precise continuum placement. We derive metallicities based on the best fit synthetic spectra to the observed spectra. To verify the accuracy of the applied atmospheric models and test our synthetic spectrum approach, three binary systems with a K-dwarf primary and an M-dwarf companion were observed and analysed along with the single M dwarfs. We obtain a good agreement between the metallicities derived for the primaries and secondaries of our test binaries and thereby confirm the reliability of our method of analysing M dwarfs. Our metallicities agree well with certain earlier determinations, and deviate from others. We conclude that spectroscopic abundance analysis in the J band is a reliable method for establishing the metallicity scale for M dwarfs. We recommend its application to a larger sample covering lower as well as higher metallicities. Further prospects of the method include abundance determinations for individual elements.

Angular momentum evolution in low-mass stars is determined by initial conditions during star formation, stellar structure evolution, and the behaviour of stellar magnetic fields. Here we show that the empirical picture of angular momentum evolution arises naturally if rotation is related to magnetic field strength instead of to magnetic flux, and formulate a corrected braking law based on this. Angular momentum evolution then becomes a strong function of stellar radius, explaining the main trends observed in open clusters and field stars at a few Gyr: the steep transition in rotation at the boundary to full convection arises primarily from the large change in radius across this boundary, and does not require changes in dynamo mode or field topology. Additionally, the data suggest transient core-envelope decoupling among solar-type stars, and field saturation at longer periods in very low-mass stars. For solar-type stars, our model is also in good agreement with the empirical Skumanich law. Finally, in further support of the theory, we show that the predicted age at which low-mass stars spin down from the saturated to unsaturated field regimes in our model corresponds remarkably well to the observed lifetime of magnetic activity in these stars.

We report simultaneous observations of the nearby flare star Proxima Centauri with VLT/UVES and XMM-Newton over three nights in March 2009. Our optical and X-ray observations cover the star's quiescent state, as well as its flaring activity and allow us to probe the stellar atmospheric conditions from the photosphere into the chromosphere, and then the corona during its different activity stages. Using the X-ray data, we investigate variations in coronal densities and abundances and infer loop properties for an intermediate-sized flare. The optical data are used to investigate the magnetic field and its possible variability, to construct an emission line list for the chromosphere, and use certain emission lines to construct physical models of Proxima Centauri's chromosphere. We report the discovery of a weak optical forbidden Fe xiii line at 3388 AA during the more active states of Proxima Centauri. For the intermediate flare, we find two secondary flare events that may originate in neighbouring loops, and discuss the line asymmetries observed during this flare in H i, He i, and Ca ii lines. The high time-resolution in the H alpha line highlights strong temporal variations in the observed line asymmetries, which re-appear during a secondary flare event. We also present theoretical modelling with the stellar atmosphere code PHOENIX to construct flaring chromospheric models.

Close M-dwarf binaries and higher multiples allow the investigation of rotational evolution and mean magnetic flux unbiased from scatter in inclination angle and age since the orientation of the spin axis of the components is most likely parallel and the individual systems are coeval. Systems composed of an early (M0.0 -- M4.0) and a late (M4.0 -- M8.0) type component offer the possibility to study differences in rotation and magnetism between partially and fully convective stars. We have selected 10 of the closest dM systems to determine the rotation velocities and the mean magnetic field strengths based on spectroscopic analysis of FeH lines of Wing-Ford transitions at 1 $\mu$m observed with VLT/CRIRES. We also studied the quality of our spectroscopic model regarding atmospheric parameters including metallicity. A modified version of the Molecular Zeeman Library (MZL) was used to compute Land\'e g-factors for FeH lines. Magnetic spectral synthesis was performed with the Synmast code. We confirmed previously reported findings that less massive M-dwarfs are braked less effectively than objects of earlier types. Strong surface magnetic fields were detected in primaries of four systems (GJ 852, GJ 234, LP 717-36, GJ 3322), and in the secondary of the triple system GJ 852. We also confirm strong 2 kG magnetic field in the primary of the triple system GJ 2005. No fields could be accurately determined in rapidly rotating stars with $\vsini>10$ \kms. For slow and moderately rotating stars we find the surface magnetic field strength to increase with the rotational velocity $\vsini$ which is consistent with other results from studying field stars.

[Abridged] We present a new radio survey of about 100 late-M and L dwarfs undertaken with the VLA. The sample was chosen to explore the role of rotation in the radio activity of ultracool dwarfs. Combining the new sample with results from our previous studies and from the literature, we compile the largest sample to date of ultracool dwarfs with radio observations and measured rotation velocities (167 objects). In the spectral type range M0-M6 we find a radio activity-rotation relation, with saturation at log(L_rad/L_bol) 10^(-7.5) above vsini~5 km/s, similar to the relation in H-alpha and X-rays. However, at spectral types >M7 the ratio of radio to bolometric luminosity increases regardless of rotation velocity, and the scatter in radio luminosity increases. In particular, while the most rapid rotators (vsini>20 km/s) exhibit "super-saturation" in X-rays and H-alpha, this effect is not seen in the radio. We also find that ultracool dwarfs with vsini>20 km/s have a higher radio detection fraction by about a factor of 3 compared to objects with vsini<10 km/s. When measured in terms of the Rossby number (Ro), the radio activity-rotation relation follows a single trend and with no apparent saturation from G to L dwarfs and down to Ro~10^-3; in X-rays and H-alpha there is clear saturation at Ro<0.1, with super-saturation beyond M7. A similar trend is observed for the radio surface flux (L_rad/R^2) as a function of Ro. The continued role of rotation in the overall level of radio activity and in the fraction of active sources, and the single trend of L_rad/L_bol and L_rad/R^2 as a function of Ro from G to L dwarfs indicates that rotation effects are important in regulating the topology or strength of magnetic fields in at least some fully-convective dwarfs. The fact that not all rapid rotators are detected in the radio provides additional support to the idea of dual dynamo states.

We present multi-epoch radio and optical observations of the M7 dwarf 2MASS J13142039+1320011. We detect a ~1 mJy source at 1.43, 4.86, 8.46 and 22.5 GHz, making it the most luminous radio emission over the widest frequency range detected from an ultracool dwarf to date. A 10 hr VLA observation reveals that the radio emission varies sinusoidally with a period of 3.89+/-0.05 hr, and an amplitude of ~30% at 4.86 GHz and ~20% at 8.46 GHz. The periodicity is also seen in circular polarization, where at 4.86 GHz the polarization reverses helicity from left- to right-handed in phase with the total intensity. An archival detection in the FIRST survey indicates that the radio emission has been stable for at least a decade. We also detect periodic photometric variability in several optical filters with a period of 3.79 hr, and measure a rotation velocity of vsini=45+/-5 km/s, in good agreement with the radio and optical periods. The period and rotation velocity allow us to place a lower limit on the radius of the source of >0.12 R_sun, about 30% larger than theoretical expectations. The properties of the radio emission can be explained with a simple model of a magnetic dipole mis-aligned relative to the stellar rotation axis, with the sinusoidal variations and helicity reversal due to the rotation of the magnetic poles relative to our line of sight. The long-term stability of the radio emission indicates that the magnetic field (and hence the dynamo) is stable on a much longer timescale than the convective turn-over time of ~0.2 yr. If the radio emission is due to the electron cyclotron maser process, the inferred magnetic field strength reaches at least 8 kG.

Close-in, giant planets are expected to influence their host stars via tidal or magnetic interaction. But are these effects strong enough in suitable targets known so far to be observed with today's instrumentation? The upsilon And system, an F8V star with a Hot Jupiter, was claimed to undergo cyclic changes in chromospheric activity indicators with its innermost planet's period. We want to investigate the stellar chromospheric and coronal activity over several months. We therefore monitored the star in X-rays as well as at optical wavelengths to test coronal and chromospheric activity indicators for planet-induced variability, making use of the Chandra X-ray Observatory as well as the echelle spectrographs FOCES and HRS at Calar Alto (Spain) and the Hobby-Eberly Telescope (Texas, US). The stellar activity level is low, as seen both in X-rays as in Ca II line fluxes; the chromospheric data show variability with the stellar rotation period. We do not find activity variations in X-rays or in the optical which can be traced back to the planet. Gaining observational evidence for star-planet interactions in X-rays remains challenging.

We carried out 3D radiative magnetohydrodynamic simulations of the convective and magnetic structure in the surface layers (uppermost part of the convection zone and photosphere) of main-sequence stars of spectral types F3 to M2. The simulation results were analyzed in terms of sizes and properties of the convection cells (granules) and magnetic flux concentrations as well as velocity, pressure, density, and temperature profiles. Our numerical simulations show for the first time a qualitative difference in the magneto-convection between solar-like stars and M dwarfs. Owing to higher surface gravity, lower opacity (resulting in higher density at optical depth unity), and more stable downflows, small-scale magnetic structures concentrate into pore-like configurations of reduced intensity. This implies that in very cool stars magnetic surface structures like plage regions and starspots significantly differ from the solar example. Such a difference would have major impact on the interpretation of Doppler imaging data and the analysis of M dwarf spectra.

We present our analysis aimed at inferring average magnetic fields in slowly-rotating solar-like stars. Using the spectral line inversion code SPINOR, we perform high-accuracy line profile fitting and investigate whether Zeeman broadening can be reliably detected in optical data of unprecedented quality. We argue that our usage of both high- and low-g_{eff} lines does provide a certain sensitivity to magnetic fields that may, indeed, be detected. However, the measurement is subject to a model dependence and prone to ambiguities, e.g. due to spectral blends. Hence, while a field may be successfully recovered, the quantification of this field is subject to large uncertainties, even for the highest-quality optical data.

The Kepler space telescope monitors over 156.000 stars with an unprecedented photometric precision. We are interested in stellar rotational periods which we find using Lomb-Scargle periodograms. This work focuses on the 306 exoplanet candidate host stars released on June 15, 2010. We present statistics on how many of them show periodic photometric variability, providing preliminary periods and estimates of stellar activity. In the future, our work will focus on spot evolution and differential rotation.

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

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

Strong surface magnetic fields are ubiquitously found in M-dwarfs with mean intensities on the order of few thousand Gauss-three orders of magnitude higher than the mean surface magnetic field of the Sun. These fields and their interaction with photospheric convection are the main source of stellar activity, which is of big interest to study links between parent stars and their planets. Moreover, the understanding of stellar magnetism, as well as the role of different dynamo-actions in particular, is impossible without explaining magnetic fields in M-dwarfs. Measuring magnetic field intensities and geometries in such cool objects, however, is strongly limited to our ability to simulate the Zeeman effect in molecular lines. In this work, we present quantitative results of modelling and analysis of the magnetic fields in selected M-dwarfs in FeH Wing-Ford lines and strong atomic lines. Some particular FeH lines are found to be the excellent probes of the magnetic field.

Magnetic fields are an important ingredient to cool star physics, and there is great interest in measuring fields and their geometry in order to understand stellar dynamos and their influence on star formation and stellar evolution. During the last few years, a large number of magnetic field measurements became available. Two main approaches are being followed to measure the Zeeman effect in cool stars; 1) the measurement of polarized light, for example to produce magnetic maps, and 2) the measurement of integrated Zeeman broadening to measure the average magnetic field strength on the stellar surface. This article briefly reviews the two methods and compares results between them that are now available for about a dozen M-type stars. It seems that we see a great variety of magnetic geometries and field strengths with typical average fields of a few kG in active M-type stars. The interpretation of geometries, however, has not yet led to a clear picture of magnetic dynamos and field configuration, and work is needed on more observational data but also on the fundamental understanding of our measurements.

M-type stars exhibit strong magnetic fields towards decreasing effective temperatures. The measurement of these fields is complicated due to missing indicators. Molecular FeH lines provide an excellent means to determine magnetic field strengths from the Zeeman broadening of magnetically sensitive lines. Our aim is the investigation of possible dependencies of the amount of sensitivity to magnetic fields from rotational quantum number, branch, and the projection of the total angular momentum onto the internuclear axis ({\Omega}). We also compare results from com-putations with those from observations. We use high resolution CRIRES spectra of the two M dwarfs GJ1002 (M5.5 inactive) and GJ1224 (M4.5 active). Individual lines are fitted by Gaussians and the obtained line depths and widths from the active and inactive star can be compared with each other. In this way, magnetically sensitive lines can be detected. For test purposes, we do the same with computed spectra of FeH. One with zero magnetic field and the other with a 2 kG magnetic field vector used in the disc integration (i.e. pure radial at the disc center). We found, in agreement with theory, that lines with high {\Omega} show strong sensitivity to magnetic fields. No obvious correlation with branch or J was found, which was also expected for lines formed in intermediate Hund's case. The computations agreed in general well with the observations, but in many cases the individual splitting of certain lines can be very different to observations.

We continue our studies on stellar latitudinal differential rotation. The presented work is a sequel of the work of Reiners et al. who studied the spectral line broadening profile of hundreds of stars of spectral types A through G at high rotational speed (vsini > 12 km/s). While most stars were found to be rigid rotators, only a few tens show the signatures of differential rotation. The present work comprises the rotational study of some 180 additional stars. The overall broadening profile is derived according to Reiners et al. from hundreds of spectral lines by least-squares deconvolution, reducing spectral noise to a minimum. Projected rotational velocities vsini are measured for about 120 of the sample stars. Differential rotation produces a cuspy line shape which is best measured in inverse wavelength space by the first two zeros of its Fourier transform. Rigid and differential rotation can be distinguished for more than 50 rapid rotators (vsini > 12 km/s) among the sample stars from the available spectra. Ten stars with significant differential rotation rates of 10-54 % are identified, which add to the few known rapid differential rotators. Differential rotation measurements of 6 % and less for four of our targets are probably spurious and below the detection limit. Including these objects, the line shapes of more than 40 stars are consistent with rigid rotation.

We present first quantitative results of the surface magnetic field measurements in selected M-dwarfs based on detailed spectra synthesis conducted simultaneously in atomic and molecular lines of the FeH Wing-Ford $F^4\,\Delta-X^4\,\Delta$ transitions. A modified version of the Molecular Zeeman Library (MZL) was used to compute Land\'e g-factors for FeH lines in different Hund's cases. Magnetic spectra synthesis was performed with the Synmast code. We show that the implementation of different Hund's case for FeH states depending on their quantum numbers allows us to achieve a good fit to the majority of lines in a sunspot spectrum in an automatic regime. Strong magnetic fields are confirmed via the modelling of atomic and FeH lines for three M-dwarfs YZ~CMi, EV~Lac, and AD~Leo, but their mean intensities are found to be systematically lower than previously reported. A much weaker field ($1.7-2$~kG against $2.7$~kG) is required to fit FeH lines in the spectra of GJ~1224. Our method allows us to measure average magnetic fields in very low-mass stars from polarized radiative transfer. The obtained results indicate that the fields reported in earlier works were probably overestimated by about $15-30$\%. Higher quality observations are needed for more definite results.

We investigate the detectability of Zeeman broadening in optical Stokes I spectra of slowly rotating sun-like stars. To this end, we apply the LTE spectral line inversion package SPINOR to very-high quality CES data and explore how fit quality depends on the average magnetic field, Bf . One-component (OC) and two-component (TC) models are adopted. In OC models, the entire surface is assumed to be magnetic. Under this assumption, we determine formal 3{\sigma} upper limits on the average magnetic field of 200 G for the Sun, and 150 G for 61 Vir (G6V). Evidence for an average magnetic field of ~ 500 G is found for 59 Vir (G0V), and of ~ 1000 G for HD 68456 (F6V). A distinction between magnetic and non-magnetic regions is made in TC models, while assuming a homogeneous distribution of both components. In our TC inversions of 59 Vir, we investigate three cases: both components have equal temperatures; warm magnetic regions; cool magnetic regions. Our TC model with equal temperatures does not yield significant improvement over OC inversions for 59 Vir. The resulting Bf values are consistent for both. Fit quality is significantly improved, however, by using two components of different temperatures. The inversions for 59 Vir that assume different temperatures for the two components yield results consistent with 0 - 450 G at the formal 3{\sigma} confidence level. We thus find a model dependence of our analysis and demonstrate that the influence of an additional temperature component can dominate over the Zeeman broadening signature, at least in optical data. Previous comparable analyses that neglected effects due to multiple temperature components may be prone to the same ambiguities.

Molecular FeH provides a large number of sharp and isolated absorption lines that can be used to measure radial velocity, rotation, or magnetic field strength with high accuracy. Our aim is to provide an FeH atlas for M-type stars in the spectral region from 986 nm to 1077 nm (Wing-Ford band). To identify these lines in CRIRES spectra of the magnetically inactive, slowly rotating, M5.5 dwarf GJ1002, we calculated model spectra for the selected spectral region with theoretical FeH line data. In general this line list agrees with the observed data, but several individual lines differ significantly in position or in line strength. After identification of as many as possible FeH lines, we correct the line data for position and line strength to provide an accurate atlas of FeH absorption lines for use in high precision spectroscopy of low mass stars. For all lines, we use a Voigt function to obtain their positions and equivalent widths. Identification with theoretical lines is done by hand. For confirmation of the identified lines, we use statistical methods, cross- correlation techniques, and line intensities. Eventually, we were able to identify FeH lines from the (0, 0), (1, 0), (1, 1), (2, 1), (2, 2), (3, 2), and (4, 3) vibrational bands in the observed spectra and correct the positions of the lines if necessary. The deviations between theoretical and observed positions follow a normal distribution approximately around zero. In order to empirically correct the line strength, we determined Teff, instrumental broadening (rotational broadening) and a van der Waals enhancement factor for FeH lines in GJ1002. We also give scaling factors for the Einstein A values to correct the line strengths. With the identified lines, we derived rotational temperatures from line intensities for GJ1002. ... .

Very little is known about magnetic fields of extrasolar planets and brown dwarfs. We use the energy flux scaling law presented by Christensen et al. (2009) to calculate the evolution of average magnetic fields in extrasolar planets and brown dwarfs under the assumption of fast rotation, which is probably the case for most of them. We find that massive brown dwarfs of about 70 M_Jup can have fields of a few kilo-Gauss during the first few hundred Million years. These fields can grow by a factor of two before they weaken after deuterium burning has stopped. Brown dwarfs with weak deuterium burning and extrasolar giant planets start with magnetic fields between ~100G and ~1kG at the age of a few Myr, depending on their mass. Their magnetic field weakens steadily until after 10Gyr it has shrunk by about a factor of 10. We use observed X-ray luminosities to estimate the age of the known extrasolar giant planets that are more massive than 0.3M_Jup and closer than 20pc. Taking into account the age estimate, and assuming sun-like wind-properties and radio emission processes similar to those at Jupiter, we calculate their radio flux and its frequency. The highest radio flux we predict comes out as 700mJy at a frequency around 150MHz for $\tau$Boob, but the flux is below 60mJy for the rest. Most planets are expected to emit radiation between a few Mhz and up to 100MHz, well above the ionospheric cutoff frequency.

In near-infrared NaCo observations of the young brown dwarf 2MASS J0041353-562112, we discovered a companion a little less than a magnitude fainter than the primary. The binary candidate has a separation of 143 mas, the spectral types are M6.5 and M9.0 for the two components. Colors and flux ratios are consistent with the components being located at the same distance minimizing the probability of the secondary being a background object. The brown dwarf is known to show Li absorption constraining the age to less than ~200 Myr, and it was suspected to show ongoing accretion, indicating an age as low as ~10 Myr. We estimate distance and orbital parameters of the binary as a function of age. For an age of 10 Myr, the distance to the system is 50 pc, the orbital period is 126 yr, and the masses of the components are ~30 and ~15 MJup. The binary brown dwarf fills a so far unoccupied region in the parameters mass and age; it is a valuable new benchmark object for brown dwarf atmospheric and evolutionary models.

Radial velocities measured from near-infrared spectra are a potentially powerful tool to search for planets around cool stars and sub-stellar objects. However, no technique currently exists that yields near-infrared radial velocity precision comparable to that routinely obtained in the visible. We describe a method for measuring high-precision relative radial velocities of these stars from K-band spectra. The method makes use of a glass cell filled with ammonia gas to calibrate the spectrograph response similar to the "iodine cell" technique that has been used very successfully in the visible. Stellar spectra are obtained through the ammonia cell and modeled as the product of a Doppler-shifted template spectrum of the object and a spectrum of the cell, convolved with a variable instrumental profile model. A complicating factor is that a significant number of telluric absorption lines are present in the spectral regions containing useful stellar and ammonia lines. The telluric lines are modeled simultaneously as well using spectrum synthesis with a time-resolved model of the atmosphere over the observatory. The free parameters in the complete model are the wavelength scale of the spectrum, the instrumental profile, adjustments to the water and methane abundances in the atmospheric model, telluric spectrum Doppler shift, and stellar Doppler shift. Tests of the method based on the analysis of hundreds of spectra obtained for late M dwarfs over six months demonstrate that precisions of ~5 m/s are obtainable over long timescales, and precisions of better than 3 m/s can be obtained over timescales up to a week. The obtained precision is comparable to the predicted photon-limited errors, but primarily limited over long timescales by the imperfect modeling of the telluric lines.

We present high-precision relative radial velocities of the very low-mass star VB10 that were obtained over a time span of 0.61 yr as part of an ongoing search for planets around stars at the end of the main sequence. The radial velocities were measured from high-resolution near-infrared spectra obtained using the CRIRES instrument on the VLT with an ammonia gas cell. The typical internal precision of the measurements is 10 m/s. These data do not exhibit significant variability and are essentially constant at a level consistent with the measurement uncertainties. Therefore, we do not detect the radial velocity variations of VB10 expected due to the presence of an orbiting giant planet similar to that recently proposed by Pravdo and Shaklan based on apparent astrometric perturbations. In addition, we do not confirm the ~1 km/s radial velocity variability of the star tentatively detected by Zapatero Osorio and colleagues with lower precision measurements. Our measurements rule out planets with M_p > 3 M_Jup and the orbital period and inclination suggested by Pravdo and Shaklan at better than 5 sigma confidence. We conclude that the planet detection claimed by Pravdo and Shaklan is spurious on the basis of this result. Although the outcome of this work is a non-detection, it illustrates the potential of using ammonia cell radial velocities to detect planets around very low-mass stars.

We present radial velocity measurements of a sample of L0-L8 dwarfs observed with VLT/UVES and Keck/HIRES. We combine these measurements with distance and proper motion from the literature to determine space motions for 43 of our targets. We identify nine candidate members of young moving groups, which have ages of 50-600 Myr according to their space motion. From the total velocity dispersion of the 43 L dwarfs, we calculate a kinematic age of ~5 Gyr for our sample. This age is significantly higher than the ~3 Gyr age known for late M dwarfs in the solar neighbourhood. We find that the distributions of the U and V velocity components of our sample are clearly non-Gaussian, placing the age estimate inferred from the full space motion vector into question. The W-component exhibits a distribution more consistent with a normal distribution, and from W alone we derive an age of ~3 Gyr, which is the same age found for late-M dwarf samples. Our brightness-limited sample is probably contaminated by a number of outliers that predominantly bias the U and V velocity components. The origin of the outliers remain unclear, but we suggest that these brown dwarfs may have gained their high velocities by means of ejection from multiple systems during their formation.

In a volume-limited sample of 63 ultracool dwarfs of spectral type M7-M9.5, we have obtained high-resolution spectroscopy with UVES at the Very Large Telescope and HIRES at Keck Observatory. In this second paper, we present projected rotation velocities, average magnetic field strengths, and chromospheric emission from the Halpha line. We confirm earlier results that the mean level of normalized Halpha luminosity decreases with lower temperature, and we find that the scatter among Halpha luminosities is larger at lower temperature. We measure average magnetic fields between 0 and 4kG with no indication for a dependence on temperature between M7 and M9.5. For a given temperature, Halpha luminosity is related to magnetic field strength, consistent with results in earlier stars. A few very slowly rotating stars show very weak magnetic fields and Halpha emission, all stars rotating faster than our detection limit show magnetic fields of at least a few hundred Gauss. In contrast to earlier-type stars, we observe magnetic fields weaker than 1kG in stars rotating faster than ~3km/s, but we find no correlation between rotation and magnetic flux generation among them. We interpret this as a fundamental change in the dynamo mechanism; in ultracool dwarfs, magnetic field generation is predominantly achieved by a turbulent dynamo, while other mechanisms can operate more efficiently at earlier spectral type.

The detection of planets around very low-mass stars with the radial velocity method is hampered by the fact that these stars are very faint at optical wavelengths where the most high-precision spectrometers operate. We investigate the precision that can be achieved in radial velocity measurements of low mass stars in the near infrared (nIR) Y-, J-, and H-bands, and we compare it to the precision achievable in the optical. For early-M stars, radial velocity measurements in the nIR offer no or only marginal advantage in comparison to optical measurements. Although they emit more flux in the nIR, the richness of spectral features in the optical outweighs the flux difference. We find that nIR measurement can be as precise than optical measurements in stars of spectral type ~M4, and from there the nIR gains in precision towards cooler objects. We studied potential calibration strategies in the nIR finding that a stable spectrograph with a ThAr calibration can offer enough wavelength stability for m/s precision. Furthermore, we simulate the wavelength-dependent influence of activity (cool spots) on radial velocity measurements from optical to nIR wavelengths. Our spot simulations reveal that the radial velocity jitter does not decrease as dramatically towards longer wavelengths as often thought. The jitter strongly depends on the details of the spots, i.e., on spot temperature and the spectral appearance of the spot. Forthcoming nIR spectrographs will allow the search for planets with a particular advantage in mid- and late-M stars. Activity will remain an issue, but simultaneous observations at optical and nIR wavelengths can provide strong constraints on spot properties in active stars.

CARMENES, Calar Alto high-Resolution search for M dwarfs with Exo-earths with a Near-infrared Echelle Spectrograph, is a study for a next-generation instrument for the 3.5m Calar Alto Telescope to be designed, built, integrated, and operated by a consortium of nine German and Spanish institutions. Our main objective is finding habitable exoplanets around M dwarfs, which will be achieved by radial velocity measurements on the m/s level in the near-infrared, where low-mass stars emit the bulk of their radiation.

We report the results of the analysis of high resolution photospheric line spectra obtained with the UVES instrument on the VLT for a sample of 15 solar-type stars selected from a recent survey of the distribution of H and K chromospheric line strengths in the solar-age open cluster M67. We find upper limits to the projected rotation velocities that are consistent with solar-like rotation (i.e., v sini ~< 2-3 km/s) for objects with Ca II chromospheric activity within the range of the contemporary solar cycle. Two solar-type stars in our sample exhibit chromospheric emission well in excess of even solar maximum values. In one case, Sanders 1452, we measure a minimum rotational velocity of vsini = 4 +/- 0.5 km/s, or over twice the solar equatorial rotational velocity. The other star with enhanced activity, Sanders 747, is a spectroscopic binary. We conclude that high activity in solar-type stars in M67 that exceeds solar levels is likely due to more rapid rotation rather than an excursion in solar-like activity cycles to unusually high levels. We estimate an upper limit of 0.2% for the range of brightness changes occurring as a result of chromospheric activity in solar-type stars and, by inference, in the Sun itself. We discuss possible implications for our understanding of angular momentum evolution in solar-type stars, and we tentatively attribute the rapid rotation in Sanders 1452 to a reduced braking efficiency.

We present an investigation of the velocity fields in early to late M-type star hydrodynamic models, and we simulate their influence on FeH molecular line shapes. The M star model parameters range between log g of 3.0 - 5.0 and Teff of 2500 K and 4000 K. Our aim is to characterize the Teff- and log g -dependence of the velocity fields and express them in terms of micro- and macro-turbulent velocities in the one dimensional sense. We present also a direct comparison between 3D hydrodynamical velocity fields and 1D turbulent velocities. The velocity fields strongly affect the line shapes of FeH, and it is our goal to give a rough estimate for the log g and Teff parameter range in which 3D spectral synthesis is necessary and where 1D synthesis suffices. In order to calculate M-star structure models we employ the 3D radiative-hydrodynamics (RHD) code CO5BOLD. The spectral synthesis on these models is performed with the line synthesis code LINFOR3D. We describe the 3D velocity fields in terms of a Gaussian standard deviation and project them onto the line of sight to include geometrical and limb-darkening effects. The micro- and macro-turbulent velocities are determined with the "Curve of Growth" method and convolution with a Gaussian velocity profile, respectively. To characterize the log g and Teff dependence of FeH lines, the equivalent width, line width, and line depth are regarded. The velocity fields in M-stars strongly depend on log g and Teff. They become stronger with decreasing log g and increasing Teff.

In a volume-limited sample of 63 ultracool dwarfs of spectral type M7--M9.5, we have obtained high-resolution spectroscopy with UVES at the Very Large Telescope and HIRES at Keck Observatory. In this first paper we introduce our volume-complete sample from DENIS and 2MASS targets, and we derive radial velocities and space motion. Kinematics of our sample are consistent with the stars being predominantly members of the young disk. The kinematic age of the sample is 3.1 Gyr. We find that six of our targets show strong Li lines implying that they are brown dwarfs younger than several hundred million years. Five of the young brown dwarfs were unrecognized before. Comparing the fraction of Li detections to later spectral types, we see a hint of an unexpected local maximum of this fraction at spectral type M9. It is not yet clear whether this maximum is due to insufficient statistics, or to a combination of physical effects including spectral appearance of young brown dwarfs, Li line formation, and the star formation rate at low masses.

We report on the discovery of the young, nearby, brown dwarf 2MASS J0041353$-$562112. The object has a spectral type of M7.5, it shows Li absorption and signatures of accretion, which implies that it still has a disk and suggests an age below 10 Myr. The space motion vector and position on the sky indicate that the brown dwarf is probably a member of the $\sim$20 Myr old Tuc-Hor association, or that it may be an ejected member of the $\sim$12 Myr old $\beta$ Pic association, both would imply that 2MASS J0041353$-$562112 may in fact be older than 10 Myr. No accreting star or brown dwarf was previously known in these associations. Assuming an age of 10 Myr, the brown dwarf has a mass of about 30 M$_{\rm Jup}$ and is located at 35 pc distance. The newly discovered object is the closest accreting brown dwarf known. Its membership to an association older than 10 Myr implies that either disks in brown dwarfs can survive as long as in more massive stars, perhaps even longer, or that star formation in Tuc-Hor or $\beta$ Pic occured more recently than previously thought. The history and evolution of this object can provide new fundamental insight into the formation process of stars, brown dwarfs, and planets.

We present an investigation of velocity fields in early to late M-type hydrodynamic stellar atmosphere models. These velocities will be expressed in classical terms of micro- and macro-turbulent velocities for usage in 1D spectral synthesis. The M-star model parameters range between log g of 3.0 - 5.0 and Teff of 2500 K - 4000 K. We characterize the Teff- and log g-dependence of the hydrodynamical velocity fields in these models with a binning method, and for the determination of micro-turbulent velocities, the Curve of Growth method is used. The macro-turbulent velocities are obtained by convolutions with Gaussian profiles. Velocity fields in M-stars strongly depend on log g and Teff. Their velocity amplitudes increase with decreasing log g and increasing Teff. The 3D hydrodynamical and 1D macro-turbulent velocities range from ~100 m/s for cool high gravity models to ~ 800 m/s - 1000 m/s for hot models or models with low log g. The micro-turbulent velocities range in the order of ~100 m/s for cool models, to ~600 m/s for hot or low log g models. Our M-star structure models are calculated with the 3D radiative-hydrodynamics (RHD) code CO5BOLD. The spectral synthesis on these models is performed with the line synthesis code LINFOR3D.

This White Paper, submitted to the National Academy of Sciences' Astro2010 Decadal Review Committee, focuses on 2 central themes in the study of young brown dwarfs -- their formation mechanism and disk characteristics -- which are of direct relevance to fundamental questions of stellar and planetary origins and properties.

We investigate the effect of stellar activity and flares on short-term radial velocity measurements in the mid-M flare star CN Leo. Radial velocity variations are calculated from 181 UVES spectra obtained during three nights. We searched for spectral orders that contain very few atmospheric absorption lines and calibrated them against the telluric A-band from O$_2$ in the Earth's atmosphere. One giant flare occurred during our observations, which has a very strong effect on radial velocity. The apparent radial velocity shift due to the flare is several hundred m s$^{-1}$ and clearly correlated with H$\alpha$ emission. Outside the flare, only spectral orders containing the most prominent emission lines of H, He, and Ca show a correlation to chromospheric activity together with a radial velocity jitter exceeding a few 10 m s$^{-1}$. We identify a number of spectral orders that are free of strong emission lines and show no flaring-related radial velocity jitter, although flares occurred as strong as 0.4 dex in normalized H$\alpha$ luminosity. The mean radial velocity jitter due to moderate flaring is less than 10 m s$^{-1}$. Strong flares are easily recognized directly in the spectra and should be neglected for planet searches.

The detection of extrasolar planets' atmospheres requires very demanding observations. For planets that can not be spatially separated from their host stars, i.e. the vast majority of planets, the transiting planets are the only ones allowing to probe their atmospheres. This is possible from transmission spectroscopy or from measurements taken during secondary eclipse. An alternative is the measurement of the Rossiter-McLaughlin-effect, which is sensitive to the size of the planetary radius. Since the radius is wavelength-dependent due to contributions of strong planetary absorption lines, this opens a path to probe planetary atmospheres also with ground-based high-resolution spectroscopy. The major goal of our numerical simulations is to provide a reliable estimate of the amplitude of the wavelength-dependent RM-effect. Our numerical simulations provide phase resolved synthetic spectra modeling the partly eclipsed stellar surface during the transit in detail. Using these spectra we can obtain RM-curves for different wavelength regions and for a wavelength-dependent planetary radius. Curves from regions with high and low contributions of absorption lines within the planetary atmosphere can be compared. From these differential effects observable quantities are derived. We applied our simulations to HD209458. Our numerical simulations show that a detailed treatment of the limb darkening for the synthetic spectra is important for a precise analysis. Compared to a parameterized limb darkening law, systematic errors of 6 m/sec occur. The wavelength dependency of the planetary atmospheres over the NaD-doublet produce a differential effect in the RM-curve of 1.5 m/sec for a star with a rotation velocity of 4.5 km/sec which increases to 4 m/sec for twice the rotation velocity.

We have measured the surface magnetic flux on four accreting young brown dwarfs and one non-accreting young very low-mass star utilizing high resolution spectra of absorption lines of the FeH molecule. A magnetic field of 1-2 kG had been proposed for one of the brown dwarfs, 2MASS J1207334$-$393254, because of its similarities to higher mass T Tauri stars as manifested in accretion and the presence of a jet. We do not find clear evidence for a kilo-Gauss field in any of our young brown dwarfs but do find a 2 kG field on the young VLM star. Our 3-$\sigma$ upper limit for the magnetic flux in 2MASS J1207334$-$393254 just reaches 1 kG. We estimate the magnetic field required for accretion in young brown dwarfs given the observed rotations, and find that fields of only a few hundred Gauss are sufficient for magnetospheric accretion. This predicted value is less than our observed upper limit. We conclude that magnetic fields in young brown dwarfs are a factor of five or more lower than in young stars of about one solar mass, and in older stars with spectral types similar to our young brown dwarfs. It is interesting that, during the first few million years, the fields scale down with mass in line with what is needed for magnetospheric accretion, yet no such scaling is observed at later ages within the same effective temperature range. This scaling is opposite to the trend in rotation, with shorter rotation periods for very young accreting brown dwarfs compared with accreting solar-mass objects (and very low Rossby numbers in all cases). We speculate that in young objects a deeper intrinsic connection may exist between magnetospheric accretion and magnetic field strength, or that magnetic field generation in brown dwarfs may be less efficient than in stars. Neither of these currently have an easy physical explanation.

Direct measurements of magnetic fields in low-mass stars of spectral class M have become available during the last years. This contribution summarizes the data available on direct magnetic measurements in M dwarfs from Zeeman analysis in integrated and polarized light. Strong magnetic fields at kilo-Gauss strength are found throughout the whole M spectral range, and so far all field M dwarfs of spectral type M6 and later show strong magnetic fields. Zeeman Doppler images from polarized light find weaker fields, which may carry important information on magnetic field generation in partially and fully convective stars.

We compare the amount of magnetic flux measured in Stokes V and Stokes I in a sample of early- and mid-M stars around the boundary to full convection (~M3.5). Early-M stars possess a radiative core, mid-M stars are fully convective. While Stokes V is sensitive to the net polarity of magnetic flux arising mainly from large-scale configurations, Stokes I measurements can see the total mean flux. We find that in early-M dwarfs, only ~6% of the total magnetic flux is detected in Stokes V. This ratio is more than twice as large, ~14%, in fully convective mid-M dwarfs. The bulk of the magnetic flux on M-dwarfs is not seen in Stokes V. This is presumably because magnetic flux is mainly stored in small scale components. There is also more to learn about the effect of the weak-field approximation on the accuracy of strong field detections. In our limited sample, we see evidence for a change in magnetic topology at the boundary to full convection. Fully convective stars store a 2-3 times higher fraction of their flux in fields visible to Stokes V. We estimate the total magnetic energy detected in Stokes I and compare it to results from Stokes V. We find that in early-M dwarfs only ~0.5% of the total magnetic energy is detected in Stokes V while this fraction is ~2.5% in mid-M dwarfs.

We present synthetic FeH band spectra in the z-filter range for several M-dwarf models with logg=3.0-5.0 [cgs] and Teff=2800K -3450K. Our aim is to characterize convective velocities in M-dwarfs and to give a rough estimate of the range in which 3D-atmosphere treatment is necessary and where 1D-atmosphere models suffice for the interpretation of molecular spectral features. This is also important in order to distinguish between the velocity-broadening and the rotational- or Zeeman-broadening. The synthetic spectra were calculated using 3D CO5BOLD radiative-hydrodynamic (RHD) models and the line synthesis code LINFOR3D. We used complete 3D-models and high resolution 3D spectral synthesis for the detailed study of some well isolated FeH lines. The FeH line strength shows a dependence on surface gravity and effective temperature and could be employed to measure both quantities in M-type objects. The line width is related to the velocity-field in the model stars, which depends strongly on surface gravity. Furthermore, we investigate the velocity-field in the 3D M-dwarf models together with the related micro- and macro-turbulent velocities in the 1D case. We also search for effects on the lineshapes.

We present first results from a multi-object spectroscopy campaign in IC2602, the Hyades, the Pleiades, and the Coma cluster using VLT/FLAMES. We analysed the data for radial velocity, rotational velocity, and H-alpha activity. Here, we highlight three aspects of this study in the context of rotational braking and the rotation-activity relationship among low mass stars. Finally we discuss the cluster membership of sources in IC2602.

We present magnetic flux measurements in seven rapidly rotating M dwarfs. Our sample stars have X-ray and H-alpha emission indicative of saturated emission, i.e., emission at a high level independent of rotation rate. Our measurements are made using near-infrared FeH molecular spectra observed with HIRES at Keck. Because of their large convective overturn times, the rotation velocity of M stars with small Rossby numbers is relatively slow and does not hamper the measurement of Zeeman splitting. The Rossby numbers of our sample stars are as small as 0.01. All our sample stars exhibit magnetic flux of kilo-Gauss strength. We find that the magnetic flux saturates in the same regime as saturation of coronal and chromospheric emission, at a critical Rossby number of around 0.1. The filling factors of both field and emission are near unity by then. We conclude that the strength of surface magnetic fields remains independent of rotation rate below that; making the Rossby number yet smaller by a factor of ten has little effect. These saturated M-star dynamos generate an integrated magnetic flux of roughly 3 kG, with a scatter of about 1 kG. The relation between emission and flux also has substantial scatter.

Today, the generation of magnetic fields in solar-type stars and its relation to activity and rotation can coherently be explained, although it is certainly not understood in its entirety. Rotation facilitates the generation of magnetic flux that couples to the stellar wind, slowing down the star. There are still many open questions, particularly at early phases (young age), and at very low mass. It is vexing that rotational braking becomes inefficient at the threshold to fully convective interiors, although no threshold in magnetic activity is seen, and the generation of large scale magnetic fields is still possible for fully convective stars. This article briefly outlines our current understanding of the rotation-magnetic field relation.

We report moderate magnetic flux of 450G < Bf < 750G (3sigma) on the nearby M5.5 flare star Proxima Centauri. A high resolution UVES spectrum was used to measure magnetic flux from Zeeman broadening in absorption lines of molecular FeH around 1mu. The magnetic flux we find is relatively weak compared with classical strong flare stars, but so are Proxima's flaring rates and actual emission levels. We compare what is known about the rotation rate, Rossby number, and activity levels in this star to relations between these quantities that have been recently being developed more generally for M dwarfs. We conclude that the magnetic flux is higher than the best estimates of the Rossby number from period measurements. On the other hand, the activity levels on Proxima Centauri are at the high end of what could be expected based on the measured field, but not so high as to exceed the natural scatter in these relations (other stars lie along this high envelope as well).

Flares on dM stars contain plasmas at very different temperatures and thus affect a wide wavelength range in the electromagnetic spectrum. While the coronal properties of flares are studied best in X-rays, the chromosphere of the star is observed best in the optical and ultraviolet ranges. Therefore, multiwavelength observations are essential to study flare properties throughout the atmosphere of a star. We analysed simultaneous observations with UVES/VLT and XMM-Newton of the active M5.5 dwarf CN Leo (Gl 406) exhibiting a major flare. The optical data cover the wavelength range from 3000 to 10000 Angstrom. From our optical data, we find an enormous wealth of chromospheric emission lines occurring throughout the spectrum. We identify a total of 1143 emission lines, out of which 154 are located in the red arm, increasing the number of observed emission lines in this red wavelength range by about a factor of 10. Here we present an emission line list and a spectral atlas. We also find line asymmetries for H I, He I, and Ca II lines. For the last, this is the first observation of asymmetries due to a stellar flare. During the flare onset, there is additional flux found in the blue wing, while in the decay phase, additional flux is found in the red wing. We interpret both features as caused by mass motions. In addition to the lines, the flare manifests itself in the enhancement of the continuum throughout the whole spectrum, inverting the normal slope for the net flare spectrum.

We present time series photometry for six partial transits of GJ436b obtained with the Fine Guidance Sensor instrument on the Hubble Space Telescope (HST). Our analysis of these data yields independent estimates of the host star's radius R_star = 0.505 +0.029/-0.020 R_sun, and the planet's orbital period P = 2.643882 +0.000060/-0.000058 d, orbital inclination i = 85.80 +0.21/-0.25 deg, mean central transit time T_c = 2454455.279241 +0.00026/-0.00025 HJD, and radius R_p = 4.90 +0.45/-0.33 R_earth. The radius we determine for the planet is larger than the previous findings from analyses of an infrared light curve obtained with the Spitzer Space Telescope. Although this discrepancy has a 92% formal significance (1.7 sigma), it might be indicative of systematic errors that still influence the analyses of even the highest-precision transit light curves. Comparisons of all the measured radii to theoretical models suggest that GJ436b has a H/He envelope of ~10% by mass. We also find that the transit times for GJ436b are constant to within 10 s over the 11 planetary orbits that the HST data span. However, the ensemble of published values exhibits a long-term drift and our mean transit time is 128 s later than that expected from the Spitzer ephemeris. The sparseness of the currently available data hinders distinguishing between an error in the orbital period or perturbations arising from an additional object in the system as the cause of the apparent trend. Assuming the drift is due to an error in the orbital period we obtain an improved estimate for it of P = 2.643904 +/- 0.000005 d. This value and our measured transit times will serve as important benchmarks in future studies of the GJ436 system. (abridged)

We present results from a high-resolution spectroscopic survey of 45 L dwarfs, which includes both very low-mass stars and brown dwarfs. Our spectra allow us to derive a significant number of new rotational velocities, and discover a slowly rotating (in projected velocity) L dwarf that allows more accurate measurement of spectroscopic rotations for these objects. We measure chromospheric activity (and often its variability) through the H$\alpha$ emission line. Our primary new result is good evidence that magnetic braking dominates the angular momentum evolution of even brown dwarfs, although spindown times appear to increase as mass decreases. We confirm that activity decreases as effective temperature decreases, though a larger fraction of L dwarfs are active than has previously been reported. Essentially all active objects are also variable. We confirm the lack of a rotation-activity connection for L dwarfs. We find a minimum limit for rotational velocities that increases with later spectral types, rising from near zero in older mid-M stars to more than 20 km s$^{-1}$ for mid-L objects. There is strong evidence that all L dwarfs are rapid rotators. We derive a braking law that can depend on either temperature or mass which can explain all the rotational results and provides an age dependence for the angular momentum evolution. It is clear that angular momentum loss mechanisms in smaller and cooler objects become more inefficient starting at the fully convective boundary.

We present a refined orbital solution for the components A, B, and C of the nearby late-M type multiple system LHS 1070. By combining astrometric datapoints from NACO/VLT, CIAO/SUBARU, and PUEO/CFHT, as well as a radial velocity measurement from the newly commissioned near infrared high-resolution spectrograph CRIRES/VLT, we achieve a very precise orbital solution for the B and C components and a first realistic constraint on the much longer orbit of the A-BC system. Both orbits appear to be co-planar. Masses for the B and C components calculated from the new orbital solution (M_(B+C) = 0.157 +/- 0.009 M_sun) are in excellent agreement with theoretical models, but do not match empirical mass-luminosity tracks. The preliminary orbit of the A-BC system reveals no mass excess for the A component, giving no indication for a previously proposed fourth (D) component in LHS 1070.

We present simultaneous high-temporal and high-spectral resolution observations at optical and soft X-ray wavelengths of the nearby flare star CN Leo. During our observing campaign a major flare occurred, raising the star's instantaneous energy output by almost three orders of magnitude. The flare shows the often observed impulsive behavior, with a rapid rise and slow decay in the optical and a broad soft X-ray maximum about 200 seconds after the optical flare peak. However, in addition to this usually encountered flare phenomenology we find an extremely short (~2 sec) soft X-ray peak, which is very likely of thermal, rather than non-thermal nature and temporally coincides with the optical flare peak. While at hard X-ray energies non-thermal bursts are routinely observed on the Sun at flare onset, thermal soft X-ray bursts on time scales of seconds have never been observed in a solar nor stellar context. Time-dependent, one-dimensional hydrodynamic modeling of this event requires an extremely short energy deposition time scale of a few seconds to reconcile theory with observations, thus suggesting that we are witnessing the results of a coronal explosion on CN Leo. Thus the flare on CN Leo provides the opportunity to observationally study the physics of the long-sought "micro-flares" thought to be responsible for coronal heating.

The bottom of the main sequence hosts objects with fundamentally different properties. At masses of about 0.3 M$_{\odot}$, stars become fully convective and at about 0.08 M$_{\odot}$ the hydrogen-burning main sequence ends; less massive objects are brown dwarfs. While stars and brown dwarfs experience very different evolutions, their inner structure has relatively little impact on the atmospheres. The generation of magnetic fields and activity is obviously connected to the threshold between partial and complete convection, because dynamo mechanisms involving a layer of shear like the solar $\alpha\Omega$-dynamo must cease. Hence a change in stellar activity can be expected there. Observations of stellar activity do not confirm a rapid break in activity at the convection boundary, but the fraction of active stars and rapid rotators is higher on the fully convective side. I summarize the current picture of stellar activity and magnetic field measurements at the bottom of the main sequence and present recent results on rotational braking beyond.

We show high resolution spectra of the eclipsing brown dwarf binary 2MASSJ05352184-0546085 taken at the two opposite radial velocity maxima. Comparisons of the TiO bands to model and template spectra are fully consistent with the temperatures previously derived for this system. In particular, the reversal of temperatures with mass - in which the higher-mass primary is cooler than its companion - is confirmed. We measure the projected rotation velocities of the compononents; the primary is rotating at least twice as rapidly as the secondary. At the two radial velocity maxima, Halpha emission lines of both components stick out to either sides of the Halpha central wavelength, which is dominated by nebula emission. This enables us to model the individual Halpha lines of the primary and the secondary. We find that the Halpha emission from the primary is at least 7 times stronger than the emission from the secondary. We conclude that the temperature reversal is very likely due to strong magnetic fields inhibiting convection on the primary.

The study of rotation and activity in low-mass stars or brown dwarfs of spectral classes M and L has seen enormous progress during the last years. I summarize the results from different works that measured activity, rotation, and sometimes magnetic fields. The generation of magnetic activity seems to be unchanged at the threshold to completely convective stars, i.e. no change in the efficiency of the magnetic dynamos is observed. On the other hand, a sudden change in the strength of rotational braking appears at the threshold mass to full convection, and strong evidence exists for rotational braking weakening with lower mass. A probable explanation is that the field topology changes from dipolar to small scale structure as the objects become fully convective.

Differential rotation can be detected in single line profiles of stars rotating more rapidly than about $v \sin{i} = 10$ km s$^{-1}$ with the Fourier transform technique. This allows to search for differential rotation in large samples to look for correlations between differential rotation and other stellar parameters. I analyze the fraction of differentially rotating stars as a function of color, rotation, and activity in a large sample of F-type stars. Color and rotation exhibit a correlation with differential rotation in the sense that more stars are rotating differentially in the cooler, less rapidly rotating stars. Effects of rotation and color, however, cannot be disentangled in the underlying sample. No trend with activity is found.

We present a UVES/VLT high resolution atlas of three L dwarfs and one T dwarf system, spectral classes at which most of the objects are brown dwarfs. Our atlas covers the optical region from H$\alpha$ up to the near infrared at 1 $\mu$m. We present spectral details of ultra-cool atmospheres at very high resolution ($R \sim 33 000$) and compare the spectra to model calculations. Our comparison shows that molecular features from VO and CaH, and atomic features from Cs and Rb are reasonably well fit by current models. On the other hand, features due to TiO, CrH, and water, and atomic Na and K reveal large discrepancies between model calculations and our observations.

We show individual high resolution spectra of components A, B, and C of the nearby late-M type multiple system LHS 1070. Component A is a mid-M star, B and C are known to have masses at the threshold to brown dwarfs. From our spectra we measure rotation velocities and the mean magnetic field for all three components individually. We find magnetic flux on the order of several kilo-Gauss in all components. The rotation velocities of the two late-M objects B and C are similar (vsini = 16km/s), the earlier A component is spinning only at about half that rate. This suggests weakening of net rotational braking at late-M spectral type, and that the lack of slowly rotating late-M and L dwarfs is real. Furthermore, we found that magnetic flux in the B component is about twice as strong as in component C at similar rotation rate. This indicates that rotational braking is not proportional to magnetic field strength in fully convective objects, and that a different field topology is the reason for the weak braking in low mass objects.

Context: Planets outside our solar system transiting their host star, i. e. those with an orbital inclination near 90 degree, are of special interest to derive physical properties of extrasolar planets. With the knowledge of the host star's physical parameters, the planetary radius can be determined. Combined with spectroscopic observations the mass and therefore the density can be derived from Doppler-measurements. Depending on the brightness of the host star, additional information, e. g. about the spin-orbit alignment between the host star and planetary orbit, can be obtained. Aims: The last few years have witnessed a growing success of transit surveys. Among other surveys, the MACHO project provided nine potential transiting planets, several of them with relatively bright parent stars. The photometric signature of a transit event is, however, insufficient to confirm the planetary nature of the faint companion. The aim of this paper therefore is a determination of the spectroscopic parameters of the host stars as well as a dynamical mass determination through Doppler-measurements. Methods: We have obtained follow-up high-resolution spectra for five stars selected from the MACHO sample, which are consistent with transits of low-luminosity objects. Radial velocities have been determined by means of cross-correlation with model spectra. The MACHO light curves have been compared to simulations based on the physical parameters of the system derived from the radial velocities and spectral analyses. Aims: We show that all transit light curves of the exoplanet candidates analysed in this work can be explained by eclipses of stellar objects, hence none of the five transiting objects is a planet.

We present UVES/VLT observations of the nearby flare star CN Leo covering the Wing-Ford FeH band near 1mu with high spectral resolution. Some of the FeH absorption lines in this band are magnetically sensitive and allow a measurement of the mean magnetic flux on CN Leo. Our observations, covering three nights separated by 48 hours each, allow a clear detection of a mean magnetic field of Bf~2.2kG. The differential flux measurements show a night-to-night variability with extremely high significance. Finally, our data strongly suggest magnetic flux variability on time scales as low as 6 hours in line with chromospheric variability.

The rotation-activity connection explains stellar activity in terms of rotation and convective overturn time. It is well established in stars of spectral types F--K and in M-type stars of young clusters. It is not established in field M-dwarfs, because they rotate very slowly, and detecting rotation periods or rotational line broadening is a challenge. In field M-dwarfs, saturation sets in below v_rot = 5 km/s, hence they are expected to populate the non-saturated part of the rotation-activity connection. This work for the first time shows intrinsically resolved spectral lines of slowly rotating M-dwarfs and makes a first comparison to estimates of convective velocities. By measuring rotation velocities in a sample of mostly inactive M-dwarfs, the unsaturated part of the rotation-activity connection is followed into the regime of very low activity. Spectra of ten M-dwarfs are taken at a resolving power of R = 200 000 in the molecular FeH absorption band. The intrinsically narrow lines are compared to model calculations and rotational broadening is measured. In one star, an upper limit of vsini = 1 km/s was found, significant rotation was detected in the other nine objects. All inactive stars show rotation below or equal to 2 km/s. In the two active stars AD Leo and YZ CMi, rotation velocities are found to be 40-50% below the results from earlier studies. The rotation activity connection holds in field early-M stars, too. Activity and rotation velocities of the sample stars are well in agreement with the relation found in earlier and younger stars. The intrinsic absorption profiles of molecular FeH lines are consistent with calculations from atomic Fe lines. Investigation of FeH line profiles is a very promising tool to measure convection patterns at the surfaces of M-stars.

We present the first direct magnetic field measurements on M dwarfs cooler than spectral class M4.5. Utilizing a new method based on the effects of a field on the FeH band near 1 micron, we obtain information on whether the integrated surface magnetic flux (Bf) is low (well under 1 kilogauss), intermediate (between 1 and about 2.5 kG), or strong (greater than about 3 kG) on a set of stars ranging from M2 down to M9. We also measure the rotational broadening (vsini) and Halpha emission for more than 20 stars. Our goal is to advance the understanding of how dynamo field production varies with stellar parameters for very low-mass stars, how the field and emission activity are related, and whether there is a connection between the rotation and magnetic flux. We find that fields are produced throughout the M-dwarfs. Among the early M stars we have too few targets to yield conclusive results. In the mid-M stars, there is a clear connection between slow rotation and weak fields. In the late-M stars, rotation is always measureable, and the strongest fields go with the most rapid rotators. These very cool rapid rotators have the largest magnetic flux in the whole sample. Halpha emission is found to be a good general proxy for magnetic fields. The drop-off in fractional emission near the bottom of the main sequence is not accompanied by a drop-off in magnetic flux, lending credence to the hypothesis that it is due to atmospheric coupling to the field rather than changes in the field itself. It is clear that the methodology we have developed can be further applied to discover more about the behavior of magnetic dynamos and magnetic activity in cool and fully convective objects.

We report on the results of a survey for radial velocity variability in a heterogeneous sample of very low-mass stars and brown dwarfs. One distinguishing characteristic of the survey is its timespan, which allows an overlap between spectroscopic binaries and those which can be found by high angular-resolution imaging. We are able to place a new constraint on the total binary fraction in these objects, which suggests that they are more likely the result of extending the same processes at work at higher masses into this mass range, rather than a distinct mode of formation. Our basic result is that there are $6 \pm 2$ out of 53, or $11^{+0.07}_{-0.04}$% spectroscopic binaries in the separation range 0-6 AU, nearly as many as resolved binaries. This leads to an estimate of an upper limit of $26 \pm 10$% for the binary fraction of VLM objects (it is an upper limit because of the possible overlap between the spectroscopic and resolved populations). A reasonable estimate for the very low-mass binary fraction is $20 - 25$%. We consider several possible separation and frequency distributions, including the same one as found for GK stars, a compressed version of that, a version of the compressed distribution truncated at 15 AU, and a theoretical distribution which considers the evaporation of small-N clusters. We conclude that the latter two bracket the observations, which may mean that these systems form with intrinsically smaller separations due to their smaller mass, and then are truncated due to their smaller binding energy. We do not find support for the ``ejection hypothesis'' as their dominant mode of formation, particularly in view of the similarity in the total binary fraction compared with slightly more massive stars, and the difficulty this mechanism has in producing numerous binary systems.

We present a new method for direct measurement of magnetic fields on ultracool stars and brown dwarfs. It takes advantage of the Wing-Ford band of FeH, which are seen throughout the M and L spectral types. These molecular features are not as blended as other optical molecular bands, are reasonably strong through most of the spectral range, and exhibit a response to magnetic fields which is easier to detect than other magnetic diagnostics, including the usual optical and near-infrared atomic spectral lines that have heretofore been employed. The FeH bands show a systematic growth as the star gets cooler. We do not find any contamination by CrH in the relevant spectral region. We are able to model cool and rapidly-rotating spectra from warmer, slowly-rotating spectra utilizing an interpolation scheme based on optical depth scaling. We show that the FeH features can distinguish between negligible, moderate, and high magnetic fluxes on low-mass dwarfs, with a current accuracy of about one kilogauss. Two different approaches to extracting the information from the spectra are developed and compared. Which one is superior depends on a number of factors. We demostrate the validity of our new procedures by comparing the spectra of three M stars whose magnetic fluxes are already known from atomic line analysis. The low and high field stars are used to produce interpolated moderate-strength spectra which closely resemble the moderate-field star. The assumption of linear behavior for the magnetic effects appears to be reasonable, but until the molecular constants are better understood the method is subject to that assumption, and rather approximate. Nonetheless, it opens a new regime of very low-mass objects to direct confirmation and testing of their magnetic dynamos.

We present the first high-resolution (R ~ 31,000) spectra of the cool sdL 2MASS0532, and what was originally identified as an early-type L subdwarf (sdL) LSR1610-0040. Our work, in combination with contemporaneous work by Cushing and Vacca, makes it clear that the latter object is more probably a mid-M dwarf with an unusual composition that gives it some sub-dwarf spectral features. We use the data to derive precise radial velocities for both objects and to estimate space motion; both are consistent with halo kinematics. We measure the projected rotational velocities, revealing very slow rotation for the old sd?M6 object \lsr. \twom exhibits rapid rotation of vsini = 65 +- 15km/s, consistent with the behavior of L dwarfs. This means that the braking time for L dwarfs is extremely long, or that perhaps they never slow down. A detailed comparison of the atomic Rb and Cs lines to spectra of field L dwarfs shows the spectral type \twom is consistent with being mid- to late-L. The Rb I and K I lines of \lsr\ are like an early-L dwarf, but the Cs I line is like a mid-M dwarf. The appearance of the Ca II triplet in absorption in this object is very hard to understand if it is not as least as warm as M6. We explain these effects in a consistent way using a mildly metal-poor mid-M model. M subdwarfs have weak metal-oxides and enhanced metal-hydrides relative to normal M dwarfs. \lsr\ exhibits metal-hydrides like an M dwarf but metal-oxides like a subdwarf. The same explanation that resolves the atomic line discrepancy explains this as well. We identify atomic lines of Ti around 9600 \AA and a small contribution of FeH, but we cannot confirm a detection of TiH in the spectra of cool L subdwarfs. High resolution spectroscopy has aided in beginning to understand the complex molecular chemistry in metal-deficient and ultralow-mass objects.

Evolutionary and atmospheric models have become available for young ultralow-mass objects. These models are being used to determine fundamental parameters from observational properties. TiO bands are used to determine effective temperatures in ultralow-mass objects, and together with Na- and K-lines to derive gravities at the substellar boundary. Unfortunately, model calibrations in (young) ultralow-mass objects are rare. As a first step towards a calibration of synthetic spectral features, I show molecular bands of TiO, which is a main opacity source in late M-dwarfs. The TiO $\epsilon$-band at 8450\AA is systematically too weak. This implies that temperatures determined from that band are underestimated, and I discuss implications for determining fundamental parameters from high resolution spectra.

More than 600 high resolution spectra of stars with spectral type F and later were obtained in order to search for signatures of differential rotation in line profiles. In 147 stars, the rotation law could be measured, 28 of them are found to be differentially rotating. Comparison to rotation laws in stars of spectral type A reveals that differential rotation sets in at the convection boundary in the HR-diagram; no star that is significantly hotter than the convection boundary exhibits the signatures of differential rotation. Four late A-/early F-type stars close to the convection boundary and at vsini~100 km/s show extraordinarily strong absolute shear at short rotation periods around one day. It is suggested that this is due to their small convection zone depth and that it is connected to a narrow range in surface velocity. Detection frequencies of differential rotation were analyzed in stars with varying temperature and rotation velocity. Measurable differential rotation is more frequent in late-type stars and slow rotators. The strength of absolute shear and differential rotation are examined as functions of the stellar effective temperature and rotation period. The strongest shear is found at rotation periods between two and three days. In slower rotators, the strongest shear at a given rotation rate is given approximately by DOmega_max ~ P^{-1}. In faster rotators, alpha_max and DOmega_max diminish less rapidly. A comparison with differential rotation measurements in stars of later spectral type shows that F-stars exhibit stronger shear than cooler stars do, the upper boundary in absolute shear DOmega with temperature is consistent with the temperature scaling law found in Doppler Imaging measurements.

We report the discovery of a new low-mass spectroscopic (SB2) stellar binary system in the star-forming region of Upper Scorpius. This object, UScoCTIO5, was discovered by Ardila (2000), who assigned it a spectral class of M4. A KeckI HIRES spectrum revealed it to be double-lined, and we then carried out a program at several observatories to determine its orbit. The orbital period is 34 days, and the eccentricity is nearly 0.3. The importance of such a discovery is that it can be used to help calibrate evolutionary models at low masses and young ages. This is one of the outstanding problems in the study of formation mechanisms and initial mass functions at low masses. The orbit allows us to place a lower limit of 0.64 +- 0.02 M_sol on the total system mass. The components appear to be of almost equal mass. We are able to show that this mass is significantly higher than predicted by evolutionary models for an object of this luminosity and age, in agreement with other recent results. More precise determination of the temperature and surface gravity of the components would be helpful in further solidifying this conclusion.

We have derived the mean broadening profile of the star V102 in the region of the open cluster IC4665 from high resolution spectroscopy. At a projected equatorial rotation velocity of vsini = (105 +- 12)km/s we find strong deviation from classical rotation. We discuss several scenarios, the most plausible being strong differential rotation in latitudinal direction. For this scenario we find a difference in angular velocity of DeltaOmega = 3.6 +- 0.8 rad/d (DeltaOmega/Omega = 0.42 +- 0.09). From the Halpha line we derive a spectral type of A9 and support photometric measurements classifying IC4665 V102 as a non-member of IC4665. At such early spectral type this is the strongest case of differential rotation observed so far. Together with three similar stars, IC4665 V102 seems to form a new class of objects that exhibit extreme latitudinal shear in a very shallow convective envelope.

We present a detailed spectroscopic study of line broadening in the A7IV-V star Altair. In a wavelength region covering 690AA we reconstruct the overall broadening profile taking into account more than 650 spectral lines. From the broadening profile we determine the projected rotational velocity vsini, derive an upper limit for the equatorial velocity v from the shape of the profile and search for signatures of differential rotation. Our redetermined value of vsini is (227 +-11)km/s. Measuring the first two zeros of the Fourier transformed broadening profile yield no signatures of differential rotation. We derive that Altair is seen under an inclination angle higher than i=68deg and it rotates at v<245km/s or slower than 53% of breakup velocity on a 1sigma level.

We reanalyzed high quality spectra of 158 stars of spectral types A0--F1 and $v \sin{i}$ between 60 and 150 km s$^{-1}$. Using a Least Squares Deconvolution technique we extracted high $S/N$ broadening profiles and determined the loci of the Fourier transform zeros $q_{1}$ and $q_{2}$ where the $S/N$-ratio was high enough. The values of $v \sin{i}$ were redetermined and found to be consistent with the values derived by \cite{Royer02a}. For 78 stars $q_{2}$ could be determined and the ratio $q_{2}/q_{1}$ was used as a shape parameter sensitive for solar-like differential rotation (Equator faster than Pole). 74 of the 78 stars have values of $q_{2}/q_{1}$ consistent with solid body rotation; in four of the 78 cases, values of $q_{2}/q_{1}$ are not consistent with rigid rotation. Although these stars may be binaries, none of the profiles shows signatures of a companion. The Fourier transforms do not indicate any distortions and the broadening profiles can be considered due to single objects. One of those candidates may be an extremely rapid rotator seen pole-on, but for the other three stars of spectral types as early as A6, differential rotation seems to be the most plausible explanation for the peculiar profiles.

We obtained high quality spectra of 135 stars of spectral types F and later and derived ``overall'' broadening functions in selected wavelength regions utilizing a Least Squares Deconvolution (LSD) procedure. Precision values of the projected rotational velocity $v \sin{i}$ were derived from the first zero of the Fourier transformed profiles and the shapes of the profiles were analyzed for effects of differential rotation. The broadening profiles of 70 stars rotating faster than $v \sin{i} = 45$ km s$^{-1}$ show no indications of multiplicity nor of spottedness. In those profiles we used the ratio of the first two zeros of the Fourier transform $q_2/q_1$ to search for deviations from rigid rotation. In the vast majority the profiles were found to be consistent with rigid rotation. Five stars were found to have flat profiles probably due to cool polar caps, in three stars cuspy profiles were found. Two out of those three cases may be due to extremely rapid rotation seen pole on, only in one case ($v \sin{i} = 52$ km s$^{-1}$) solar-like differential rotation is the most plausible explanation for the observed profile. These results indicate that the strength of differential rotation diminishes in stars rotating as rapidly as $v \sin{i} \ga 50$ km s$^{-1}$.

Mechanisms influencing absorption line profiles of fast rotating stars can be sorted into two groups; (i) intrinsic variations sensitive to temperature and pressure, and (ii) global effects common to all spectral lines. I present a detailed study on the latter effects focusing on gravity darkening and inclination for various rotational velocities and spectral types. It is shown that the line shapes of rapidly and rigidly rotating stars mainly depend on the equatorial velocity $v_{\rm e}$, not on the projected rotational velocity $v \sin{i}$ which determines the lines width. The influence of gravity darkening and spectral type on the line profiles is shown. The results demonstrate the possibility of determining the inclination angle $i$ of single fast rotators, and they show that constraints on gravity darkening can be drawn for stellar samples. While significant line profile deformation occurs in stars rotating as fast as $v_{\rm e} \ga 200$ km s$^{-1}$, for slower rotators profile distortion are marginal. In these cases spectral signatures induced by, e.g., differential rotation are not affected by gravity darkening and the methods applicable to slow rotators can be applied to these faster rotators, too.

We report detection of differential rotation on the F5 dwarf psi Cap using line profile analysis. The Fourier transform of both FeI lambda 5775 and SiI lambda 5772 are used to obtain a projected rotational velocity of v sini = (42+-1)km/s. Modelling of the Fourier transformed profiles shows that the combined effects of equatorial velocity, inclination and differential rotation dominate the line profile while limb darkening and turbulence velocities have only minor effects. Rigid rotation is shown to be inconsistent with the measured profiles. Modelling the line profiles analogous to solar differential rotation we find a differential rotation parameter of alpha = 0.15+-0.1 (15+-10%) comparable to the solar case. To our knowledge this is the first successful measurement of differential rotation through line profile analysis.