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.

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

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

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

The nearby ultra-compact multiplanetary system YZ Ceti consists of at least three planets. The orbital period of each planet is the subject of discussion in the literature due to strong aliasing in the radial velocity data. The stellar activity of this M dwarf also hampers significantly the derivation of the planetary parameters. With an additional 229 radial velocity measurements obtained since the discovery publication, we reanalyze the YZ Ceti system and resolve the alias issues. We use model comparison in the framework of Bayesian statistics and periodogram simulations based on a method by Dawson and Fabrycky to resolve the aliases. We discuss additional signals in the RV data, and derive the planetary parameters by simultaneously modeling the stellar activity with a Gaussian process regression model. To constrain the planetary parameters further we apply a stability analysis on our ensemble of Keplerian fits. We resolve the aliases: the three planets orbit the star with periods of $2.02$ d, $3.06$ d, and $4.66$ d. We also investigate an effect of the stellar rotational signal on the derivation of the planetary parameters, in particular the eccentricity of the innermost planet. Using photometry we determine the stellar rotational period to be close to $68$ d. From the absence of a transit event with TESS, we derive an upper limit of the inclination of $i_\mathrm{max} = 87.43$ deg. YZ Ceti is a prime example of a system where strong aliasing hindered the determination of the orbital periods of exoplanets. Additionally, stellar activity influences the derivation of planetary parameters and modeling them correctly is important for the reliable estimation of the orbital parameters in this specific compact system. Stability considerations then allow additional constraints to be placed on the planetary parameters.

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

Helicity is a fundamental property of a magnetic field but to date it has only been possible to observe its evolution in one star - the Sun. In this paper we provide a simple technique for mapping the large-scale helicity density across the surface of any star using only observable quantities: the poloidal and toroidal magnetic field components (which can be determined from Zeeman-Doppler imaging) and the stellar radius. We use a sample of 51 stars across a mass range of 0.1-1.34 M$_\odot$ to show how the helicity density relates to stellar mass, Rossby number, magnetic energy and age. We find that the large-scale helicity density increases with decreasing Rossby number $R_o$, peaking at $R_o \simeq 0.1$, with a saturation or decrease below that. For both fully- and partially-convective stars we find that the mean absolute helicity density scales with the mean squared toroidal magnetic flux density according to the power law: $|\langle{h\,}\rangle|$ $\propto$ $\langle{\rm{B_{tor}}^2_{}\,\rangle}^{0.86\,\pm\,0.04}$. The scatter in this relation is consistent with the variation across a solar cycle, which we compute using simulations and observations across solar cycles 23 and 24 respectively. We find a significant decrease in helicity density with age.

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

Main sequence low-mass stars are known to spin-down as a consequence of their magnetised stellar winds. However, estimating the precise rate of this spin-down is an open problem. The mass-loss rate, angular momentum-loss rate and the magnetic field properties of low-mass stars are fundamentally linked making this a challenging task. Of particular interest is the stellar magnetic field geometry. In this work, we consider whether non-dipolar field modes contribute significantly to the spin-down of low-mass stars. We do this using a sample of stars that have all been previously mapped with Zeeman-Doppler imaging. For a given star, as long as its mass-loss rate is below some critical mass-loss rate, only the dipolar fields contribute to its spin-down torque. However, if it has a larger mass-loss rate, higher order modes need to be considered. For each star, we calculate this critical mass-loss rate, which is a simple function of the field geometry. Additionally, we use two methods of estimating mass-loss rates for our sample of stars. In the majority of cases, we find that the estimated mass-loss rates do not exceed the critical mass-loss rate and hence, the dipolar magnetic field alone is sufficient to determine the spin-down torque. However, we find some evidence that, at large Rossby numbers, non-dipolar modes may start to contribute.

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.

Previous studies have related surface temperature maps, obtained with the Doppler imaging (DI) technique, of LQ Hya with long-term photometry. We compare surface magnetic field maps, obtained with the Zeeman Doppler imaging (ZDI) technique, with contemporaneous photometry, with the aim of quantifying the star's magnetic cycle characteristics. We inverted Stokes IV spectropolarimetry into magnetic field and surface brightness maps using a tomographic inversion code that models high signal-to-noise ratio mean line profiles produced by the least squares deconvolution (LSD) technique. The magnetic field and surface brightness maps reveal similar patterns to previous DI and ZDI studies: non-axisymmetric polar magnetic field structure, void of fields at mid-latitudes, and a complex structure in the equatorial regions. There is a weak but clear tendency of the polar structures to be linked with strong radial field and the equatorial ones with the azimuthal. We find a polarity reversal in the radial field between 2016 and 2017 coincident with an activity minimum seen in the long-term photometry. The inverted field strengths cannot easily be related with the observed spottedness, but we find that they are partially connected with the retrieved field complexity. Comparing to global magnetoconvection models for rapidly rotating young Suns, this field topology and dominance of the poloidal field component could be explained by a turbulent dynamo, where differential rotation does not play a major role (so called alpha^2 Omega or alpha^2 dynamos), and axi- and non-axisymmetric modes are excited simultaneously. The complex equatorial magnetic field structure could arise from the twisted (helical) wreaths often seen in these simulations, while the polar feature would be connected to the mostly poloidal non-axisymmetric component having a smooth spatial structure.

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.

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.

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.

In this work, we simulate the evolution of the solar wind along its main sequence lifetime and compute its thermal radio emission. To study the evolution of the solar wind, we use a sample of solar mass stars at different ages. All these stars have observationally-reconstructed magnetic maps, which are incorporated in our 3D magnetohydrodynamic simulations of their winds. We show that angular-momentum loss and mass-loss rates decrease steadily on evolutionary timescales, although they can vary in a magnetic cycle timescale. Stellar winds are known to emit radiation in the form of thermal bremsstrahlung in the radio spectrum. To calculate the expected radio fluxes from these winds, we solve the radiative transfer equation numerically from first principles. We compute continuum spectra across the frequency range 100 MHz - 100 GHz and find maximum radio flux densities ranging from 0.05 - 8.3 $\mu$Jy. At a frequency of 1 GHz and a normalised distance of d = 10 pc, the radio flux density follows 0.24 $(\Omega/\Omega_{\odot})^{0.9}$ (d/[10pc])$^2$ $\mu$Jy, where $\Omega$ is the rotation rate. This means that the best candidates for stellar wind observations in the radio regime are faster rotators within distances of 10 pc, such as $\kappa^1$ Ceti (2.83 $\mu$Jy) and $\chi^1$ Ori (8.3 $\mu$Jy). These flux predictions provide a guide to observing solar-type stars across the frequency range 0.1 - 100 GHz in the future using the next generation of radio telescopes, such as ngVLA and SKA.

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.

Low-mass stars are known to have magnetic fields that are believed to be of dynamo origin. Two complementary techniques are principally used to characterise them. Zeeman-Doppler imaging (ZDI) can determine the geometry of the large-scale magnetic field while Zeeman broadening can assess the total unsigned flux including that associated with small-scale structures such as spots. In this work, we study a sample of stars that have been previously mapped with ZDI. We show that the average unsigned magnetic flux follows an activity-rotation relation separating into saturated and unsaturated regimes. We also compare the average photospheric magnetic flux recovered by ZDI, $\langle B_V\rangle$, with that recovered by Zeeman broadening studies, $\langle B_I\rangle$. In line with previous studies, $\langle B_V\rangle$ ranges from a few % to $\sim$20% of $\langle B_I\rangle$. We show that a power law relationship between $\langle B_V\rangle$ and $\langle B_I\rangle$ exists and that ZDI recovers a larger fraction of the magnetic flux in more active stars. Using this relation, we improve on previous attempts to estimate filling factors, i.e. the fraction of the stellar surface covered with magnetic field, for stars mapped only with ZDI. Our estimated filling factors follow the well-known activity-rotation relation which is in agreement with filling factors obtained directly from Zeeman broadening studies. We discuss the possible implications of these results for flux tube expansion above the stellar surface and stellar wind models.

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.

We have observed the active star $\xi$ Boo A (HD 131156A) with high precision broadband linear polarimetry contemporaneously with circular spectropolarimetry. We find both signals are modulated by the 6.43 day rotation period of $\xi$ Boo A. The signals from the two techniques are 0.25 out of phase, consistent with the broadband linear polarization resulting from differential saturation of spectral lines in the global transverse magnetic field. The mean magnitude of the linear polarization signal is ~4 ppm/G but its structure is complex and the amplitude of the variations suppressed relative to the longitudinal magnetic field. The result has important implications for current attempts to detect polarized light from hot Jupiters orbiting active stars in the combined light of the star and planet. In such work stellar activity will manifest as noise, both on the time scale of stellar rotation, and on longer time scales - where changes in activity level will manifest as a baseline shift between observing runs.

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

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

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

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

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.

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.

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.

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.

One of the aims of the BCool programme is to search for cycles in other stars and to understand how similar they are to the Sun. In this paper we aim to monitor the evolution of $\tau$ Boo's large-scale magnetic field using high-cadence observations covering its chromospheric activity maximum. For the first time, we detect a polarity switch that is in phase with $\tau$ Boo's 120 day chromospheric activity maximum and its inferred X-ray activity cycle maximum. This means that $\tau$ Boo has a very fast magnetic cycle of only 240 days. At activity maximum $\tau$ Boo's large-scale field geometry is very similar to the Sun at activity maximum: it is complex and there is a weak dipolar component. In contrast, we also see the emergence of a strong toroidal component which has not been observed on the Sun, and a potentially overlapping butterfly pattern where the next cycle begins before the previous one has finished.

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

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.

Magnetic activity is known to be correlated to the rotation period for moderately active main sequence solar-like stars. In turn, the stellar rotation period evolves as a result of magnetised stellar winds that carry away angular momentum. Understanding the interplay between magnetic activity and stellar rotation is therefore a central task for stellar astrophysics. Angular momentum evolution models typically employ spin-down torques that are formulated in terms of the surface magnetic field strength. However, these formulations fail to account for the magnetic field geometry, unlike those that are expressed in terms of the open flux, i.e. the magnetic flux along which stellar winds flow. In this work, we model the angular momentum evolution of main sequence solar-mass stars using a torque law formulated in terms of the open flux. This is done using a potential field source surface model in conjunction with the Zeeman-Doppler magnetograms of a sample of roughly solar-mass stars. We explore how the open flux of these stars varies with stellar rotation and choice of source surface radii. We also explore the effect of field geometry by using two methods of determining the open flux. The first method only accounts for the dipole component while the second accounts for the full set of spherical harmonics available in the Zeeman-Doppler magnetogram. We find only a small difference between the two methods, demonstrating that the open flux, and indeed the spin-down, of main sequence solar-mass stars is likely dominated by the dipolar component of the magnetic field.

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.

The young and magnetically active K dwarf Epsilon Eridani exhibits a chromospheric activity cycle of about 3 years. Previous reconstructions of its large-scale magnetic field show strong variations at yearly epochs. To understand how Epsilon Eridani's large-scale magnetic field geometry evolves over its activity cycle we focus on high cadence observations spanning 5 months at its activity minimum. Over this timespan we reconstruct 3 maps of Epsilon Eridani's large-scale magnetic field using the tomographic technique of Zeeman Doppler Imaging. The results show that at the minimum of its cycle, Epsilon Eridani's large-scale field is more complex than the simple dipolar structure of the Sun and 61 Cyg A at minimum. Additionally we observe a surprisingly rapid regeneration of a strong axisymmetric toroidal field as Epsilon Eridani emerges from its S-index activity minimum. Our results show that all stars do not exhibit the same field geometry as the Sun and this will be an important constraint for the dynamo models of active solar-type stars.

We aim to understand how stellar parameters such as mass and rotation impact the distribution of starspots on the stellar surface. To this purpose, we have used Doppler imaging to reconstruct the surface brightness distributions of three fully convective M dwarfs with similar rotation rates. We secured high cadence spectral time series observations of the 5.5 AU separation binary GJ 65, comprising GJ 65A (M5.5V, Prot = 0.24 d) and GJ 65B (M6V, Prot = 0.23 d). We also present new observations of GJ 791.2A (M4.5V, Prot = 0.31 d). Observations of each star were made on two nights with UVES, covering a wavelength range from 0.64 - 1.03 microns. The time series spectra reveal multiple line distortions, which we interpret as cool starspots and which are persistent on both nights suggesting stability on the timescale of 3 days. Spots are recovered with resolutions down to 8.3 degs at the equator. The global spot distributions for GJ 791.2A are similar to observations made a year earlier. Similar high latitude and circumpolar spot structure is seen on GJ791.2A and GJ 65A. However, they are surprisingly absent on GJ 65B, which instead reveals more extensive, larger, spots concentrated at intermediate latitudes. All three stars show small amplitude latitude-dependent rotation that is consistent with solid body rotation. We compare our measurements of differential rotation with previous Doppler imaging studies and discuss the results in the wider context of other observational estimates and recent theoretical predictions.

Magnetic activity and rotation are known to be intimately linked for low-mass stars. Understanding rotation evolution over the stellar lifetime is therefore an important goal within stellar astrophysics. In recent years, there has been increased focus on how the complexity of the stellar magnetic field affects the rate of angular momentum-loss from a star. This is a topic that Zeeman-Doppler imaging (ZDI), a technique that is capable of reconstructing the large-scale magnetic field topology of a star, can uniquely address. Using a potential field source surface model, we estimate the open flux, mass loss-rate and angular momentum-loss rates for a sample of 66 stars that have been mapped with ZDI. We show that the open flux of a star is predominantly determined by the dipolar component of its magnetic field for our choice of source surface radius. We also show that, on the main sequence, the open flux, mass- and angular momentum-loss rates increase with decreasing Rossby number. The exception to this rule is stars less massive than $0.3M_{\odot}$. Previous work suggests that low mass M dwarfs may possess either strong, ordered and dipolar fields or weak and complex fields. This range of field strengths results in a large spread of angular momentum-loss rates for these stars and has important consequences for their spin down behaviour. Additionally, our models do not predict a transition in the mass-loss rates at the so called wind dividing line noted from Ly$\alpha$ studies.

Magnetic fields are essential ingredients of many physical processes in the interiors and envelopes of cool stars. Yet their direct detection and characterisation is notoriously difficult, requiring high-quality observations and advanced analysis techniques. Significant progress has been recently achieved by several types of direct magnetic field studies on the surfaces of cool active stars. In particular, complementary techniques of the field topology mapping with polarisation data and total magnetic flux measurements from intensity spectra have been systematically applied to different classes of active stars leading to interesting and occasionally controversial results. In this paper we summarise the current status of direct magnetic field studies of cool stars, and investigations of surface inhomogeneities caused by the field, based on the material presented at the Cool Stars 19 splinter session.

Accounting for stellar activity is a crucial component of the search for ever-smaller planets orbiting stars of all spectral types. We use Doppler imaging methods to demonstrate that starspot induced radial velocity variability can be effectively reduced for moderately rotating, fully convective stars. Using starspot distributions extrapolated from sunspot observations, we adopt typical M dwarf starspot distributions with low contrast spots to synthesise line profile distortions. The distortions are recovered using maximum entropy regularised fitting and the corresponding stellar radial velocities are measured. The procedure is demonstrated for a late-M star harbouring an orbiting planet in the habitable zone. The technique is effective for stars with vsini = 1-10 km/s, reducing the stellar noise contribution by factors of nearly an order of magnitude. With a carefully chosen observing strategy, the technique can be used to determine the stellar rotation period and is robust to uncertainties such as unknown stellar inclination. While demonstrated for late-type M stars, the procedure is applicable to all spectral types.

The magnetic fields, activity and dynamos of young solar-type stars can be empirically studied using time-series of spectropolarimetric observations and tomographic imaging techniques such as Doppler imaging and Zeeman Doppler imaging. In this paper we use these techniques to study the young Sun-like star EK Draconis (Sp-Type: G1.5V, HD 129333) using ESPaDOnS at the Canada-France-Hawaii Telescope and NARVAL at the T\`elescope Bernard Lyot. This multi-epoch study runs from late 2006 until early 2012. We measure high levels of chromospheric activity indicating an active, and varying, chromosphere. Surface brightness features were constructed for all available epochs. The 2006/7 and 2008 data show large spot features appearing at intermediate-latitudes. However, the 2012 data indicate a distinctive polar spot. We observe a strong, almost unipolar, azimuthal field during all epochs that is similar to that observed on other Sun-like stars. Using magnetic features, we determined an average equatorial rotational velocity, \Omega_eq, of 2.50 +/- 0.08 rad/d. High levels of surface differential rotation were measured with an average rotational shear, \Delta\Omega, of 0.27 +0.24-0.26 rad/d. During an intensively observed 3-month period from December 2006 until February 2007, the magnetic field went from predominantly toroidal ( approx. 80%) to a more balanced poloidal-toroidal (approx. 40-60%) field. Although the large-scale magnetic field evolved over the epochs of our observations, no polarity reversals were found in our data.

We present a spectropolarimetric snapshot survey of solar-type planet hosting stars. In addition to 14 planet-hosting stars observed as part of the BCool magnetic snapshot survey, we obtained magnetic observations of a further 19 planet-hosting solar-type stars in order to see if the presence of close-in planets had an effect on the measured surface magnetic field (|B$_{\ell}$|). Our results indicate that the magnetic activity of this sample is congruent with that of the overall BCool sample. The effects of the planetary systems on the magnetic activity of the parent star, if any, are too subtle to detect compared to the intrinsic dispersion and correlations with rotation, age and stellar activity proxies in our sample. Four of the 19 newly observed stars, two of which are subgiants, have unambiguously detected magnetic fields and are future targets for Zeeman Doppler Mapping.

The magnetic geometry of the surface magnetic fields of more than 55 cool stars have now been mapped using spectropolarimetry. In order to better understand these observations, we compare the magnetic field topology at different surface scale sizes of observed and simulated cool stars. For ease of comparison between the high-resolution non-potential magnetofrictional simulations and the relatively low-resolution observations, we filter out the small-scale field in the simulations using a spherical harmonics decomposition. We show that the large-scale field topologies of the solar-based simulations produce values of poloidal/toroidal fields and fractions of energy in axisymmetric modes that are similar to the observations. These global non-potential evolution model simulations capture key magnetic features of the observed solar-like stars through the processes of surface flux transport and magnetic flux emergence. They do not, however, reproduce the magnetic field of M-dwarfs or stars with dominantly toroidal field. Furthermore, we analyse the magnetic field topologies of individual spherical harmonics for the simulations and discover that the dipole is predominately poloidal, while the quadrupole shows the highest fraction of toroidal fields. Magnetic field structures smaller than a quadrupole display a fixed ratio between the poloidal and toroidal magnetic energies.

Whether BP Piscium (BP Psc) is either a pre-main sequence T Tauri star at d ~ 80 pc, or a post-main sequence G giant at d ~ 300 pc is still not clear. As a first-ascent giant, it is the first to be observed with a molecular and dust disk. Alternatively, BP Psc would be among the nearest T Tauri stars with a protoplanetary disk (PPD). We investigate whether the disk geometry resembles typical PPDs, by comparing polarimetric images with radiative transfer models. Our VLT/SPHERE/ZIMPOL observations allow us to perform Polarimetric Differential Imaging; Reference Star Differential Imaging; and Richardson-Lucy deconvolution. We present the first visible light polarization and intensity images of the disk of BP Psc. Our deconvolution confirms the disk shape as detected before, mainly showing the southern side of the disk. In polarized intensity the disk is imaged at larger detail and also shows the northern side, giving it the typical shape of high inclination flared disks. We explain the observed disk features by retrieving the large-scale geometry with MCMax radiative transfer modeling, which yields a strongly flared model, atypical for disks of T Tauri stars.

Zeeman Doppler imaging has successfully mapped the large-scale magnetic fields of stars over a large range of spectral types, rotation periods and ages. When observed over multiple epochs, some stars show polarity reversals in their global magnetic fields. On the Sun, polarity reversals are a feature of its activity cycle. In this paper, we examine the magnetic properties of stars with existing chromospherically determined cycle periods. Previous authors have suggested that cycle periods lie on multiple branches, either in the cycle period-Rossby number plane or the cycle period-rotation period plane. We find some evidence that stars along the active branch show significant average toroidal fields that exhibit large temporal variations while stars exclusively on the inactive branch remain dominantly poloidal throughout their entire cycle. This lends credence to the idea that different shear layers are in operation along each branch. There is also evidence that the short magnetic polarity switches observed on some stars are characteristic of the inactive branch while the longer chromospherically determined periods are characteristic of the active branch. This may explain the discrepancy between the magnetic and chromospheric cycle periods found on some stars. These results represent a first attempt at linking global magnetic field properties obtained form ZDI and activity cycles.

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.

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.

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.

We present six epochs of spectropolarimetric observations of the hot-Jupiter-hosting star $\tau$ Bo\"otis that extend the exceptional previous multi-year data set of its large-scale magnetic field. Our results confirm that the large-scale magnetic field of $\tau$ Bo\"otis varies cyclicly, with the observation of two further magnetic reversals; between December 2013 and May 2014 and between January and March 2015. We also show that the field evolves in a broadly solar-type manner in contrast to other F-type stars. We further present new results which indicate that the chromospheric activity cycle and the magnetic activity cycles are related, which would indicate a very rapid magnetic cycle. As an exemplar of long-term magnetic field evolution, $\tau$ Bo\"otis and this longterm monitoring campaign presents a unique opportunity for studying stellar magnetic cycles.

We present new wind models for {\tau} Bo\"otis ({\tau} Boo), a hot-Jupiter-host-star whose observable magnetic cycles makes it a uniquely useful target for our goal of monitoring the temporal variability of stellar winds and their exoplanetary impacts. Using spectropolarimetric observations from May 2009 to January 2015, the most extensive information of this type yet available, to reconstruct the stellar magnetic field, we produce multiple 3D magnetohydrodynamic stellar wind models. Our results show that characteristic changes in the large-scale magnetic field as the star undergoes magnetic cycles produce changes in the wind properties, both globally and locally at the position of the orbiting planet. Whilst the mass loss rate of the star varies by only a minimal amount ($\sim$ 4 percent), the rates of angular momentum loss and associated spin-down timescales are seen to vary widely (up to $\sim$ 140 percent), findings consistent with and extending previous research. In addition, we find that temporal variation in the global wind is governed mainly by changes in total magnetic flux rather than changes in wind plasma properties. The magnetic pressure varies with time and location and dominates the stellar wind pressure at the planetary orbit. By assuming a Jovian planetary magnetic field for {\tau} Boo b, we nevertheless conclude that the planetary magnetosphere can remain stable in size for all observed stellar cycle epochs, despite significant changes in the stellar field and the resulting local space weather environment.

We report magnetic field measurements for Kappa1~Cet, a proxy of the young Sun when life arose on Earth. We carry out an analysis of the magnetic properties determined from spectropolarimetric observations and reconstruct its large-scale surface magnetic field to derive the magnetic environment, stellar winds and particle flux permeating the interplanetary medium around Kappa1~Cet. Our results show a closer magnetosphere and mass-loss rate of Mdot = 9.7 x 10^{-13} Msol/yr, i.e., a factor 50 times larger than the current solar wind mass-loss rate, resulting in a larger interaction via space weather disturbances between the stellar wind and a hypothetical young-Earth analogue, potentially affecting the planet's habitability. Interaction of the wind from the young Sun with the planetary ancient magnetic field may have affected the young Earth and its life conditions

The surface rotation rates of young solar-type stars vary rapidly with age from the end of the pre-main sequence through the early main sequence. Important changes in the dynamos operating in these stars may result from this evolution, which should be observable in their surface magnetic fields. Here we present a study aimed at observing the evolution of these magnetic fields through this critical time period. We observed stars in open clusters and stellar associations of known ages, and used Zeeman Doppler Imaging to characterize their complex magnetic large-scale fields. Presented here are results for 15 stars, from 5 associations, with ages from 20 to 250 Myr, masses from 0.7 to 1.2 solar masses, and rotation periods from 0.4 to 6 days. We find complex large-scale magnetic field geometries, with global average strengths from 14 to 140 G. There is a clear trend towards decreasing average large-scale magnetic field strength with age, and a tight correlation between magnetic field strength and Rossby number. Comparing the magnetic properties of our zero-age main sequence sample to those of both younger and older stars, it appears that the magnetic evolution of solar-type stars during the pre-main sequence is primarily driven by structural changes, while it closely follows the stars' rotational evolution on the main sequence.

(shortened for arXiv) We aim to progress towards more efficient exoplanet detection around active stars by optimizing the use of Doppler Imaging in radial velocity measurements. We propose a simple method to simultaneously extract a brightness map and a set of orbital parameters through a tomographic inversion technique derived from classical Doppler mapping. Based on the maximum entropy principle, the underlying idea is to determine the set of orbital parameters that minimizes the information content of the resulting Doppler map. We carry out a set of numerical simulations to perform a preliminary assessment of the robustness of our method, using an actual Doppler map of the very active star HR 1099 to produce a realistic synthetic data set for various sets of orbital parameters of a single planet in a circular orbit. Using a simulated time-series of 50 line profiles affected by a peak-to-peak activity jitter of 2.5 km/s, we are able in most cases to recover the radial velocity amplitude, orbital phase and orbital period of an artificial planet down to a radial velocity semi-amplitude of the order of the radial velocity scatter due to the photon noise alone (about 50 m/s in our case). One noticeable exception occurs when the planetary orbit is close to co-rotation, in which case significant biases are observed in the reconstructed radial velocity amplitude, while the orbital period and phase remain robustly recovered. The present method constitutes a very simple way to extract orbital parameters from heavily distorted line profiles of active stars, when more classical radial velocity detection methods generally fail. It is easily adaptable to most existing Doppler Imaging codes, paving the way towards a systematic search for close-in planets orbiting young, rapidly-rotating stars.

Wood et al suggested that mass-loss rate is a function of X-ray flux ($\dot{M} \propto F_x^{1.34}$) for dwarf stars with $F_x \lesssim F_{x,6} \equiv 10^6$ erg cm$^{-2}$ s$^{-1}$. However, more active stars do not obey this relation. These authors suggested that the break at $F_{x,6}$ could be caused by significant changes in magnetic field topology that would inhibit stellar wind generation. Here, we investigate this hypothesis by analysing the stars in Wood et al's sample that had their surface magnetic fields reconstructed through Zeeman-Doppler Imaging (ZDI). Although the solar-like outliers in the $\dot{M}$ -- $F_x$ relation have higher fractional toroidal magnetic energy, we do not find evidence of a sharp transition in magnetic topology at $F_{x,6}$. To confirm this, further wind measurements and ZDI observations at both sides of the break are required. As active stars can jump between states with highly toroidal to highly poloidal fields, we expect significant scatter in magnetic field topology to exist for stars with $F_x \gtrsim F_{x,6}$. This strengthens the importance of multi-epoch ZDI observations. Finally, we show that there is a correlation between $F_x$ and magnetic energy, which implies that $\dot{M}$ -- magnetic energy relation has the same qualitative behaviour as the original $\dot{M}$ -- $F_x$ relation. No break is seen in any of the $F_x$ -- magnetic energy

Since M4.5 - M9 dwarfs exhibit equatorial rotation velocities of order 10 km/s on average, radial velocity surveys targeting this stellar population will likely need to find methods to effectively remove starspot jitter. We present the first high resolution Doppler images of the M4.5 dwarf, GJ 791.2A, and the M9 dwarf, LP 944-20. The time series spectra of both objects reveal numerous line profile distortions over the rotation period of each star which we interpret as starspots. The transient distortions are modelled with spot/photosphere contrast ratios that correspond to model atmosphere temperature differences of Tphot-Tspot = 200 and 300 K. GJ 791.2A is a fully convective star with vsini = 35.1 km/s. Although we find more starspot structure at high latitudes, we reconstruct spots at a range of phases and latitudes with a mean spot filling of ~3%. LP 944-20 is one of the brightest known late-M dwarfs, with spectral type M9V and vsini = 30.8 km/s. Its spectral time series exhibits two dominant transient line distortions that are reconstructed as high latitude spots, while a mean spot filling factor of only 1.5% is found. The occurrence of low-contrast spots at predominantly high latitudes, which we see in both targets here, is in general likely to be responsible for the low amplitude photometric variability seen in late-M dwarfs. For GJ 791.2A, the radial velocities induced by the starspot features yield an rms velocity variability of 138 m/s, which can be reduced by a factor of 1.9 using our reconstructed surface brightness distributions.

Spectropolarimetric observations have been used to map stellar magnetic fields, many of which display strong bands of azimuthal fields that are toroidal. A number of explanations have been proposed to explain how such fields might be generated though none are definitive. In this paper, we examine the toroidal fields of a sample of 55 stars with magnetic maps, with masses in the range 0.1-1.5$\,{\rm M}_\odot$. We find that the energy contained in toroidal fields has a power law dependence on the energy contained in poloidal fields. However the power index is not constant across our sample, with stars less and more massive than 0.5$\,{\rm M}_\odot$ having power indices of 0.72$\pm$0.08 and 1.25$\pm$0.06 respectively. There is some evidence that these two power laws correspond to stars in the saturated and unsaturated regimes of the rotation-activity relation. Additionally, our sample shows that strong toroidal fields must be generated axisymmetrically. The latitudes at which these bands appear depend on the stellar rotation period with fast rotators displaying higher latitude bands than slow rotators. The results in this paper present new constraints for future dynamo studies.

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.

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 the magnetic fields of young solar-type stars provide a way to investigate the signatures of their magnetic activity and dynamos. Spectropolarimetry enables the study of these stellar magnetic fields and was thus employed at the T\'{e}lescope Bernard Lyot and the Anglo-Australian Telescope to investigate two moderately rotating young Sun-like stars, namely HD 35296 (V119 Tau, HIP 25278) and HD 29615 (HIP 21632). The results indicate that both stars display rotational variation in chromospheric indices consistent with their spot activity, with variations indicating a probable long-term cyclic period for HD 35296. Additionally, both stars have complex, and evolving, large-scale surface magnetic fields with a significant toroidal component. High levels of surface differential rotation were measured for both stars. For the F8V star HD 35296 a rotational shear of $\Delta\Omega$ = 0.22$^{+0.04}_{-0.02}$ rad/d was derived from the observed magnetic profiles. For the G3V star HD 29615 the magnetic features indicate a rotational shear of $\Delta\Omega$ = 0.48$_{-0.12}^{+0.11}$ rad/d, while the spot features, with a distinctive polar spot, provide a much lower value of $\Delta\Omega$ of 0.07$_{-0.03}^{+0.10}$ rad/d. Such a significant discrepancy in shear values between spot and magnetic features for HD 29615 is an extreme example of the variation observed for other lower-mass stars. From the extensive and persistent azimuthal field observed for both targets it is concluded that a distributed dynamo operates in these moderately rotating Sun-like stars, in marked contrast to the Sun's interface-layer dynamo.

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

Understanding the formation of planetary nebulae remains elusive because in the preceding asymtotic giant branch (AGB) phase these stars are heavily enshrouded in an optically thick dusty envelope. To further understand the morphology of the circumstellar environments of AGB stars we observe the closest carbon-rich AGB star IRC+10216 in scattered light. When imaged in scattered light at optical wavelengths, IRC+10216 surprisingly shows a narrow equatorial density enhancement, in contrast to the large-scale spherical rings that have been imaged much further out. We use radiative transfer models to interpret this structure in terms of two models: firstly, an equatorial density enhancement, commonly observed in the more evolved post-AGB stars, and secondly, in terms of a dust rings model, where a local enhancement of mass-loss creates a spiral ring as the star rotates. We conclude that both models can be used to reproduce the dark lane in the scattered light images, which is caused by an equatorially density enhancement formed by dense dust rather than a bipolar outflow as previously thought. We are unable to place constraints on the formation of the equatorial density enhancement by a binary system.

The large-scale magnetic field of solar-type stars reconstructed from their spectropolarimetric observations provide important insight into their underlying dynamo processes.We aim to investigate the temporal variability of the large-scale surface magnetic field and chromospheric activity of a young solar analogue, the G0 dwarf HN Peg.The large-scale surface magnetic field topology is reconstructed using Zeeman Doppler Imaging at six observational epochs covering seven years.We also investigated the chromospheric activity variations by measuring the flux in the line cores of the three chromospheric activity indicators: Ca II H&K, H alpha, and the Ca II IRT lines.The magnetic topology of HN Peg shows a complex and variable geometry. While the radial field exhibits a stable positive polarity magnetic region at the poles at each observational epoch, the azimuthal field is strongly variable in strength, where a strong band of positive polarity magnetic field is present at equatorial latitudes. This field disappears during the middle of our time span, reappearing again during the last two epochs of observations. The mean magnetic field derived from the magnetic maps also follow a similar trend to the toroidal field, with the field strength at a minimum in epoch 2009.54. Summing the line of sight magnetic field over the visible surface at each observation, HN Peg exhibits a weak longitudinal magnetic field ranging from -14 G to 13 G, with no significant long-term trend, although there is significant rotational variability within each epoch. Those chromospheric activity indicators exhibit more long-term variations over the time span of observations, where the minimal is observed in Epoch 2008.71.

Context: The principle definition of habitability for exoplanets is whether they can sustain liquid water on their surfaces, i.e. that they orbit within the habitable zone. However, the planet's magnetosphere should also be considered, since without it, an exoplanet's atmosphere may be eroded away by stellar winds. Aims: The aim of this paper is to investigate magnetospheric protection of a planet from the effects of stellar winds from solar-mass stars. Methods: We study hypothetical Earth-like exoplanets orbiting in the host star's habitable zone for a sample of 124 solar-mass stars. These are targets that have been observed by the Bcool collaboration. Using two wind models, we calculate the magnetospheric extent of each exoplanet. These wind models are computationally inexpensive and allow the community to quickly estimate the magnetospheric size of magnetised Earth-analogues orbiting cool stars. Results: Most of the simulated planets in our sample can maintain a magnetosphere of ~5 Earth radii or larger. This suggests that magnetised Earth analogues in the habitable zones of solar analogues are able to protect their atmospheres and is in contrast to planets around young active M dwarfs. In general, we find that Earth-analogues around solar-type stars, of age 1.5 Gyr or older, can maintain at least a Paleoarchean Earth sized magnetosphere. Our results indicate that planets around 0.6 - 0.8 solar-mass stars on the low activity side of the Vaughan-Preston gap are the optimum observing targets for habitable Earth analogues.

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.

PolarBase is an evolving data base that contains all stellar data collected with the ESPaDOnS and NARVAL high-resolution spectropolarimeters, in their reduced form, as soon as they become public. As of early 2014, observations of 2,000 stellar objects throughout the Hertzsprung-Russell diagram are available. Intensity spectra are available for all targets, and the majority of the observations also include simultaneous spectra in circular or linear polarization, with the majority of the polarimetric measurements being performed only in circularly polarized light (Stokes V). Observations are associated with a cross-correlation pseudo-line profile in all available Stokes parameters, greatly increasing the detectability of weak polarized signatures. Stokes V signatures are detected for more than 300 stars of all masses and evolutionary stages, and linear polarization is detected in 35 targets. The detection rate in Stokes V is found to be anti-correlated with the stellar effective temperature. This unique set of Zeeman detections offers the first opportunity to run homogeneous magnetometry studies throughout the H-R diagram. The web interface of PolarBase is available at http://polarbase.irap.omp.eu.

We investigate how the observed large-scale surface magnetic fields of low-mass stars (~0.1 -- 2 Msun), reconstructed through Zeeman-Doppler imaging (ZDI), vary with age t, rotation and X-ray emission. Our sample consists of 104 magnetic maps of 73 stars, from accreting pre-main sequence to main-sequence objects (1 Myr < t < 10 Gyr). For non-accreting dwarfs we empirically find that the unsigned average large-scale surface field <|Bv|> is related to age as $t^{-0.655 \pm 0.045}$. This relation has a similar dependence to that identified by Skumanich (1972), used as the basis for gyrochronology. Likewise, our relation could be used as an age-dating method ("magnetochronology"). The trends with rotation we find for the large-scale stellar magnetism are consistent with the trends found from Zeeman broadening measurements (sensitive to large- and small-scale fields). These similarities indicate that the fields recovered from both techniques are coupled to each other, suggesting that small- and large-scale fields could share the same dynamo field generation processes. For the accreting objects, fewer statistically significant relations are found, with one being a correlation between the unsigned magnetic flux and rotation period. We attribute this to a signature of star-disc interaction, rather than being driven by the dynamo.

We present radial velocity measurements of a sample of M5V-M9V stars from our Red-Optical Planet Survey, ROPS, operating at 0.65-1.025 micron. Radial velocities for 15 stars, with r.m.s. precision down to 2.5 m/s over a week long time scale are achieved using Thorium-Argon reference spectra. We are sensitive to planets with m_psin(i) >= 1.5 MEarth (3 MEarth at 2-sigma) in the classical habitable zone and our observations currently rule out planets with m_psin(i) > 0.5 MJup at 0.03 AU for all our targets. A total of 9 of the 15 targets exhibit r.m.s. < 16 m/s, which enables us to rule out the presence of planets with m_psin(i) > 10 MEarth in 0.03 AU orbits. Since the mean rotation velocity is of order 8 km/s for an M6V star and 15 km/s by M9V, we avoid observing only slow rotators that would introduce a bias towards low axial inclination i << 90 deg systems, which are unfavourable for planet detection. Our targets with the highest vsini values exhibit radial velocities significantly above the photon-noise limited precision, even after accounting for vsini. We therefore monitored stellar activity via chromospheric emission from the Halpha and Ca II infrared triplet lines. A clear trend of log_10(Halpha/Lbol) with radial velocity r.m.s. is seen, implying that significant starspot activity is responsible for the observed radial velocity precision floor. The implication that most late M dwarfs are significantly spotted, and hence exhibit time varying line distortions, indicates that observations to detect orbiting planets need strategies to reliably mitigate against the effects of activity induced radial velocity variations.

Zeeman-Doppler imaging is a spectropolarimetric technique that is used to map the large-scale surface magnetic fields of stars. These maps in turn are used to study the structure of the stars' coronae and winds. This method, however, misses any small-scale magnetic flux whose polarisation signatures cancel out. Measurements of Zeeman broadening show that a large percentage of the surface magnetic flux may be neglected in this way. In this paper we assess the impact of this 'missing flux' on the predicted coronal structure and the possible rates of spin down due to the stellar wind. To do this we create a model for the small-scale field and add this to the Zeeman-Doppler maps of the magnetic fields of a sample of 12 M dwarfs. We extrapolate this combined field and determine the structure of a hydrostatic, isothermal corona. The addition of small-scale surface field produces a carpet of low-lying magnetic loops that covers most of the surface, including the stellar equivalent of solar 'coronal holes' where the large-scale field is opened up by the stellar wind and hence would be X-ray dark. We show that the trend of the X-ray emission measure with rotation rate (the so-called 'activity-rotation relation') is unaffected by the addition of small-scale field, when scaled with respect to the large-scale field of each star. The addition of small-scale field increases the surface flux; however, the large-scale open flux that governs the loss of mass and angular momentum in the wind remains unaffected. We conclude that spin-down times and mass loss rates calculated from surface magnetograms are unlikely to be significantly influenced by the neglect of small-scale field.

The circumstellar environments of classical T Tauri stars are challenging to directly image because of their high star-to-disk contrast ratio. One method to overcome this is by using imaging polarimetry where scattered and consequently polarised starlight from the star's circumstellar disk can be separated from the unpolarised light of the central star. We present images of the circumstellar environment of SU Aur, a classical T Tauri star at the transition of T Tauri to Herbig stars. The images directly show that the disk extends out to ~500 au with an inclination angle of $\sim$ 50$^\circ$. Using interpretive models, we derived very small grains in the surface layers of its disk, with a very steep size- and surface-density distribution. Additionally, we resolved a large and extended nebulosity in our images that is most likely a remnant of the prenatal molecular cloud. The position angle of the disk, determined directly from our images, rules out a polar outflow or jet as the cause of this large-scale nebulosity.

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.

Among the solar proxies studied in the Sun in Time k Cet (HD 20630) stands out as potentially having a mass very close to solar and a young age. On this study, we monitored the magnetic field and the chromospheric activity from the Ca II H & K lines of k Cet. We used the Least-Square-Deconvolution (LSD, Donati et al. 1997) by simultaneously extracting the information contained in all 8,000 photospheric lines of the echelogram (K1 type star). To reconstruct a reliable magnetic map and characterize the surface differential rotation of k Cet we used 14 exposures spread over 2 months, in order to cover at least two rotational cycles (Prot ~ 9.2 days). The LSD technique was applied to detect the Zeeman signature of the magnetic field in each of our 14 observations and to measure its longitudinal component. In order to reconstruct the magnetic field geometry of k Cet, we applied the Zeeman Doppler Imaging (ZDI) inversion method. ZDI revealed a structure in the radial magnetic field consisting of a polar magnetic spot. On this study, we present the fisrt look results of a high-resolution spectropolarimetric campaign to characterize the activity and the magnetic fields of this young solar proxy.

As part of the Bcool project, over 150 solar-type stars chosen mainly from planet search databases have been observed between 2006 and 2013 using the NARVAL and ESPaDOnS spectropolarimeters on the Telescope Bernard Lyot (Pic du Midi, France) and the Canada France Hawaii Telescope (Mauna Kea, USA), respectively. These single 'snapshot' observations have been used to detect the presence of magnetic fields on 40% of our sample, with the highest detection rates occurring for the youngest stars. From our observations we have determined the mean surface longitudinal field (or an upper limit for stars without detections) and the chromospheric surface fluxes, and find that the upper envelope of the absolute value of the mean surface longitudinal field is directly correlated to the chromospheric emission from the star and increases with rotation rate and decreases with age.

High-resolution spectropolarimetric observations provide simultaneous information about stellar magnetic field topologies and three-dimensional distributions of chemical elements. Here we present analysis of a unique full Stokes vector spectropolarimetric data set, acquired for the cool magnetic Ap star HD 24712. The goal of our work is to examine circular and linear polarization signatures inside spectral lines and to study variation of the stellar spectrum and magnetic observables as a function of rotational phase. HD 24712 was observed with the HARPSpol instrument at the 3.6-m ESO telescope over a period of 2010-2011. The resulting spectra have S/N ratio of 300-600 and resolving power exceeding 100000. The multiline technique of least-squares deconvolution (LSD) was applied to combine information from the spectral lines of Fe-peak and rare-earth elements. We used the HARPSPol spectra of HD 24712 to study the morphology of the Stokes profile shapes in individual spectral lines and in LSD Stokes profiles corresponding to different line masks. From the LSD Stokes V profiles we measured the longitudinal component of the magnetic field, <Bz>, with an accuracy of 5-10 G. We also determined the net linear polarization from the LSD Stokes Q and U profiles. We determined an improved rotational period of the star, P_rot = 12.45812 +/- 0.00019d. We measured <Bz> from the cores of Halpha and Hbeta lines. The analysis of <Bz> measurements showed no evidence for a significant radial magnetic field gradient in the atmosphere of HD 24712. We used our <Bz> and net linear polarization measurements to determine parameters of the dipolar magnetic field topology. We found that magnetic observables can be reasonably well reproduced by the dipolar model. We discovered rotational modulation of the Halpha core and related it a non-uniform surface distribution of rare-earth elements.

Several recent spectrophotometric studies failed to detect significant global magnetic fields in late-B HgMn chemically peculiar stars, but some investigations have suggested the presence of strong unstructured or tangled fields in these objects. We used detailed spectrum synthesis analysis to search for evidence of tangled magnetic fields in high-quality observed spectra of 8 slowly rotating HgMn stars and one normal late-B star. We also evaluated recent sporadic detections of weak longitudinal magnetic fields in HgMn stars based on the moment technique. Our analysis of the Zeeman broadening of magnetically sensitive spectral lines reveals no evidence of tangled magnetic fields in any of the studied HgMn or normal stars. We infer upper limits of 200-700 G for the mean magnetic field modulus -- much smaller than the field strengths implied by studies based on differential magnetic line intensification and quadratic field diagnostics. The new HARPSpol longitudinal field measurements for the extreme HgMn star HD 65949 and the normal late-B star 21 Peg are consistent with zero at a precision of 3-6 G. Re-analysis of our Stokes V spectra of the spotted HgMn star HD 11753 shows that the recent moment technique measurements retrieved from the same data are incompatible with the lack of circular polarization signatures in the spectrum of this star. We conclude that there is no evidence for substantial tangled magnetic fields on the surfaces of studied HgMn stars. We cannot independently confirm the presence of very strong quadratic or marginal longitudinal fields for these stars, so results from the moment technique are likely to be spurious.

Context. Many protoplanetary nebulae show strong asymmetries in their surrounding shell, pointing to asymmetries during the mass loss phase. Questions concerning the origin and the onset of deviations from spherical symmetry are important for our understanding of the evolution of these objects. Here we focus on the circumstellar shell of the post-AGB star HD 161796. Aims. We aim at detecting signatures of an aspherical outflow, as well as to derive the properties of it. Methods. We use the imaging polarimeter ExPo (the extreme polarimeter), a visitor instrument at the William Herschel Telescope, to accurately image the dust shell surrounding HD 161796 in various wavelength filters. Imaging polarimetry allows us to separate the faint, polarized, light from circumstellar material from the bright, unpolarized, light from the central star. Results. The shell around HD 161796 is highly aspherical. A clear signature of an equatorial density enhancement can be seen. This structure is optically thick at short wavelengths and changes its appearance to optically thin at longer wavelengths. In the classification of the two different appearances of planetary nebulae from HST images it changes from being classified as DUPLEX at short wavelengths to SOLE at longer wavelengths. This strengthens the interpretation that these two appearances are manifestations of the same physical structure. Furthermore, we find that the central star is hotter than often assumed and the relatively high observed reddening is due to circumstellar rather than interstellar extinction.

We present high spectral resolution ($R\approx108,000$) Stokes $V$ polarimetry of the Classical T Tauri stars (CTTSs) GQ Lup and TW Hya obtained with the polarimetric upgrade to the HARPS spectrometer on the ESO 3.6 m telescope. We present data on both photospheric lines and emission lines, concentrating our discussion on the polarization properties of the \ion{He}{1} emission lines at 5876 \AA\ and 6678 \AA. The \ion{He}{1} lines in these CTTSs contain both narrow emission cores, believed to come from near the accretion shock region on these stars, and broad emission components which may come from either a wind or the large scale magnetospheric accretion flow. We detect strong polarization in the narrow component of the two \ion{He}{1} emission lines in both stars. We observe a maximum implied field strength of $6.05 \pm 0.24$ kG in the 5876 \AA\ line of GQ Lup, making it the star with the highest field strength measured in this line for a CTTS. We find field strengths in the two \ion{He}{1} lines that are consistent with each other, in contrast to what has been reported in the literature on at least one star. We do not detect any polarization in the broad component of the \ion{He}{1} lines on these stars, strengthening the conclusion that they form over a substantially different volume relative the formation region of the narrow component of the \ion{He}{1} lines.

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.

We present an analysis of the photometric variability of M dwarfs in the WFCAM Transit Survey. Although periodic lightcurve variability in low mass stars is generally dominated by photospheric star spot activity, M dwarf variability in the J band has not been as thoroughly investigated as at visible wavelengths. Spectral type estimates for a sample of over 200,000 objects are made using spectral type-colour relations, and over 9600 dwarfs (J<17) with spectral types later than K7 were found. The light curves of the late-type sample are searched for periodicity using a Lomb-Scargle periodogram analysis. A total of 68 periodic variable M dwarfs are found in the sample with periods ranging from 0.16 days to 90.33 days, with amplitudes in the range of ~0.009 to ~0.115 in the J band. We simulate active M dwarfs with a range of latitude-independent spot coverages and estimate a periodically variable fractions of 1-3 per cent for stars where spots cover more than 10 per cent of the star's surface. Our simulated spot distributions indicate that operating in the J band, where spot contrast ratios are minimised, enables variability in only the most active of stars to be detected. These findings affirm the benefits of using the $J$ band for planetary transit searches compared to visible bands. We also serendipitously find a \Delta J>0.2 mag flaring event from an M4V star in our sample.

Well over 700 exoplanets have been detected to date. Only a handful of these have been observed directly. Direct observation is extremely challenging due to the small separation and very large contrast involved. Imaging polarimetry offers a way to decrease the contrast between the unpolarized starlight and the light that has become linearly polarized after scattering by circumstellar material. This material can be the dust and debris found in circumstellar disks, but also the atmosphere or surface of an exoplanet. We present the design, calibration approach, polarimetric performance and sample observation results of the Extreme Polarimeter, an imaging polarimeter for the study of circumstellar environments in scattered light at visible wavelengths. The polarimeter uses the beam-exchange technique, in which the two orthogonal polarization states are imaged simultaneously and a polarization modulator swaps the polarization states of the two beams before the next image is taken. The instrument currently operates without the aid of Adaptive Optics. To reduce the effects of atmospheric seeing on the polarimetry, the images are taken at a frame rate of 35 fps, and large numbers of frames are combined to obtain the polarization images. Four successful observing runs have been performed using this instrument at the 4.2 m William Herschel Telescope on La Palma, targeting young stars with protoplanetary disks as well as evolved stars surrounded by dusty envelopes. In terms of fractional polarization, the instrument sensitivity is better than 10^-4. The contrast achieved between the central star and the circumstellar source is of the order 10^-6. We show that our calibration approach yields absolute polarization errors below 1%.

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.

Z CMa is a complex binary system, composed of a Herbig Be and an FU Ori star. The Herbig star is surrounded by a dust cocoon of variable geometry, and the whole system is surrounded by an infalling envelope. Previous spectropolarimetric observations have reported a preferred orientation of the polarization angle, perpendicular to the direction of a large, parsec-sized jet associated with the Herbig star. The variability in the amount of polarized light has been associated to changes in the geometry of the dust cocoon that surrounds the Herbig star. We aim to constrain the properties of Z CMa by means of imaging polarimetry at optical wavelengths. Using ExPo, a dual-beam imaging polarimeter which operates at optical wavelengths, we have obtained imaging (linear) polarimetric data of Z CMa. Our observations were secured during the return to quiescence after the 2008 outburst. We detect three polarized features over Z CMa. Two of these features are related to the two jets reported in this system: the large jet associated to the Herbig star, and the micro-jet associated to the FU Ori star. Our results suggest that the micro-jet extends to a distance ten times larger than reported in previous studies. The third feature suggests the presence of a hole in the dust cocoon that surrounds the Herbig star of this system. According to our simulations, this hole can produce a pencil beam of light that we see scattered off the low-density envelope surrounding the system.

We present the first results from our Red Optical Planet Survey (ROPS) to search for low mass planets orbiting late type dwarfs (M5.5V - M9V) in their habitable zones (HZ). Our observations, with the red arm of the MIKE spectrograph (0.5 - 0.9 microns) at the 6.5 m Magellan Clay telescope at Las Campanas Observatory indicate that >= 92 per cent of the flux lies beyond 0.7 microns. We use a novel approach that is essentially a hybrid of the simultaneous iodine and ThAr methods for determining precision radial velocities. We apply least squares deconvolution to obtain a single high S/N ratio stellar line for each spectrum and cross correlate against the simultaneously observed telluric line profile, which we derive in the same way. Utilising the 0.62 - 0.90 micron region, we have achieved an r.m.s. precision of 10 m/s for an M5.5V spectral type star with spectral S/N ~160 on 5 minute timescales. By M8V spectral type, a precision of ~30 m/s at S/N = 25 is suggested, although more observations are needed. An assessment of our errors and scatter in the radial velocity points hints at the presence of stellar radial velocity variations. Of our sample of 7 stars, 2 show radial velocity signals at 6-sigma and 10-sigma of the cross correlation uncertainties. If the signals are planetary in origin, our findings are consistent with estimates of Neptune mass planets that predict a frequency of 13 - 27 per cent for early M dwarfs.Our current analysis indicates the we can achieve a sensitivity that is equivalent to the amplitude induced by a 6 M_Earth planet orbiting in the habitable zone. Based on simulations, we estimate that <10 M_Earth habitable zone planets will be detected in a new stellar mass regime, with <=20 epochs of observations.

We present recent polarimetric images of the highly variable star R CrB using ExPo and archival WFPC2 images from the HST. We observed R CrB during its current dramatic minimum where it decreased more than 9 mag due to the formation of an obscuring dust cloud. Since the dust cloud is only in the line-of-sight, it mimics a coronograph allowing the imaging of the star's circumstellar environment. Our polarimetric observations surprisingly show another scattering dust cloud at approximately 1.3" or 2000 AU from the star. We find that to obtain a decrease in the stellar light of 9 mag and with 30% of the light being reemitted at infrared wavelengths (from R CrB's SED) the grains in R CrB's circumstellar environment must have a very low albedo of approximately 0.07%. We show that the properties of the dust clouds formed around R CrB are best fitted using a combination of two distinct populations of grains size. The first are the extremely small 5 nm grains, formed in the low density continuous wind, and the second population of large grains (~0.14 {\mu}m) which are found in the ejected dust clouds. The observed scattering cloud, not only contains such large grains, but is exceptionally massive compared to the average cloud.

Mercury-manganese (HgMn) stars have been considered as non-magnetic and non-variable chemically peculiar (CP) stars for a long time. However, recent discoveries of the variability in spectral line profiles have suggested an inhomogeneous surface distribution of chemical elements in some HgMn stars. From the studies of other CP stars it is known that magnetic field plays a key role in the formation of surface spots. All attempts to find magnetic fields in HgMn stars have yielded negative results. In this study, we investigate the possible presence of a magnetic field in phi Phe (HD 11753) and reconstruct surface distribution of chemical elements that show variability in spectral lines.We also test a hypothesis that a magnetic field is concentrated in chemical spots and look into the possibility that some chemical elements are stratified with depth in the stellar atmosphere. Our analysis is based on high-quality spectropolarimetric time-series observations, covering a full rotational period of the star. Spectra were obtained with the HARPSpol at the ESO 3.6-m telescope. Combining information from all suitable spectral lines, we set an upper limit of 4 G on the mean longitudinal magnetic field. For chemical spots, an upper limit on the longitudinal field varies between 8 and 15 G. We confirmed the variability of Y, Sr, and Ti and detected variability in Cr lines. Stratification analysis showed that Y and Ti are not concentrated in the uppermost atmospheric layers. Our spectropolarimetric observations rule out the presence of a strong, globally-organised magnetic field in phi Phe. This implies an alternative mechanism of spot formation, which could be related to a non-equilibrium atomic diffusion. However, the typical time scales of the variation in stratification predicted by the recent time-dependent diffusion models exceed significantly the spot evolution time-scale reported for phi Phe.

Chemically peculiar stars of the mercury-manganese (HgMn) type represent a new class of spotted late-B stars, in which evolving surface chemical inhomogeneities are apparently unrelated to the presence of strong magnetic fields but are produced by some hitherto unknown astrophysical mechanism. The goal of this study is to perform a detailed line profile variability analysis and carry out a sensitive magnetic field search for one of the brightest HgMn stars - mu Lep. We acquired a set of very high-quality intensity and polarization spectra of mu Lep with the HARPSpol polarimeter. These data were analyzed with the multiline technique of least-squares deconvolution in order to extract information on the magnetic field and line profile variability. Our spectra show very weak but definite variability in the lines of Sc, all Fe-peak elements represented in the spectrum of mu Lep, as well as Y, Sr, and Hg. Variability might also be present in the lines of Si and Mg. Anomalous profile shapes of Ti II and Y II lines suggest a dominant axisymmetric distribution of these elements. At the same time, we found no evidence of the magnetic field in mu Lep, with the 3 sigma upper limit of only 3 G for the mean longitudinal magnetic field. This is the most stringent upper limit on the possible magnetic field derived for a spotted HgMn star. The very weak variability detected for many elements in the spectrum mu Lep suggests that low-contrast chemical inhomogeneities may be common in HgMn stars and that they have not been recognized until now due to the limited precision of previous spectroscopic observations and a lack of time-series data. The null result of the magnetic field search reinforces the conclusion that formation of chemical spots in HgMn stars is not magnetically driven.

Imaging polarimetry is a powerful tool for detecting and characterizing exoplanets and circumstellar environments. Polarimetry allows a separation of the light coming from an unpolarized source such as a star and the polarized source such as a planet or a protoplanetary disk. Future facilities like SPHERE at the VLT or EPICS at the E-ELT will incorporate imaging polarimetry to detect exoplanets. The Extreme Polarimeter (ExPo) is a dual-beam imaging polarimeter that currently can reach contrast ratios of 10^5, enough to characterize circumstellar environments. We present the data reduction steps for a dual-beam imaging polarimeter that can reach contrast ratios of 10^5. The data obtained with ExPo at the William Herschel Telescope (WHT) are analyzed. Instrumental artifacts and noise sources are discussed for an unpolarized star and for a protoplanetary disk (AB Aurigae). The combination of fast modulation and dual-beam techniques allow us to minimize instrumental artifacts. A proper data processing and alignment of the images is fundamental when dealing with large contrasts. Imaging polarimetry proves to be a powerful method to resolve circumstellar environments even without a coronagraph mask or an Adaptive Optics system.

The application of high-resolution spectropolarimetry has led to major progress in understanding the magnetism and activity of late-type stars. During the last decade, magnetic fields have been discovered and mapped for many types of active cool stars using spectropolarimetric data. However, these observations and modeling attempts are fundamentally incomplete since they are based on the interpretation of the circular polarization alone. Taking advantage of the newly built HARPS polarimeter, we have obtained the first systematic observations of several cool active stars in all four Stokes parameters. Here we report the detection of magnetically-induced linear polarization for the primary component of the very active RS CVn binary HR 1099 and for the moderately active K dwarf epsilon Eri. For both stars the amplitude of linear polarization signatures is measured to be ~10^{-4} of the unpolarized continuum, which is approximately a factor of ten lower than for circular polarization. This is the first detection of the linear polarization in line profiles of cool active stars. Our observations of the inactive solar-like star alpha Cen A show neither circular nor linear polarization above the level of 10^{-5}, indicating the absence of a net longitudinal magnetic field stronger than 0.2 G.

According to our current understanding, a subclass of the upper main sequence chemically peculiar stars, called mercury-manganese (HgMn), is non-magnetic. Nevertheless, chemical inhomogeneities were recently discovered on their surfaces. At the same time, no global magnetic fields stronger than 1-100 G are detected by modern studies. The goals of our study are to search for magnetic field in the HgMn binary system 66 Eri and to investigate chemical spots on the stellar surfaces of both components. Our analysis is based on high quality spectropolarimetric time-series observations obtained during 10 consecutive nights with the HARPSpol instrument at the ESO 3.6-m telescope. To increase the sensitivity of the magnetic field search we employed a least-squares deconvolution (LSD). We used spectral disentangling to measure radial velocities and study line profile variability. Chemical spot geometry was reconstructed using multi-line Doppler imaging. We report a non-detection of magnetic field in 66 Eri, with error bars 10-24 G for the longitudinal field. Circular polarization profiles also do not indicate any signatures of complex surface magnetic fields. For a simple dipolar field configuration we estimated an upper limit of the polar field strength to be 60-70 G. For the HgMn component we found variability in spectral lines of Ti, Ba, Y, and Sr with the rotational period equal to the orbital one. The surface maps of these elements reconstructed with the Doppler imaging technique, show relative underabundance on the hemisphere facing the secondary component. The contrast of chemical inhomogeneities ranges from 0.4 for Ti to 0.8 for Ba.

In light of the growing interest in searching for low mass, rocky planets, we investigate the impact of starspots on radial velocity searches for earth-mass planets in orbit about M dwarf stars. Since new surveys targeting M dwarfs will likely be carried out at infrared wavelengths, a comparison between V and Y band starspot induced jitter is made, indicating a reduction of up to an order of magnitude when observing in the Y band. The exact reduction in jitter is dependent on the photosphere to spot contrast ratio, with greater improvements at smaller contrasts. We extrapolate a model used to describe solar spot distributions to simulate the spot patterns that we expect to find on M dwarfs. Under the assumption that M dwarfs are near or fully convective, we randomly place starspots on the stellar surface, simulating different levels of spot coverage. Line profiles, distorted by spots are derived and are used to investigate the starspot induced jitter. By making assumptions about the degree of spot activity, detection limits for earth-mass planets in habitable zones are simulated for between 10 and 500 observation epochs. We find that <= 50 epochs are required to detect 1 - 2 MEarth planets (with < 1 per cent false alarm probability) orbiting slowly rotating 0.1 and 0.2 MSun stars. This sensitivity decreases when typical rotation velocities and activity levels for each stellar mass/spectral type are considered. No detections of below 20 MEarth planets are expected for <= 500 observations for the most active stars with vsini >= 20 km/s and dark spots.

We performed a highly sensitive search for magnetic fields on a large set of HgMn stars. With the aid of a new polarimeter attached to the HARPS spectrometer at the ESO 3.6m-telescope, we obtained high-quality circular polarization spectra of 41 single and double HgMn stars. Using a multi-line analysis technique on each star, we co-added information from hundreds of spectral lines resulting in significantly greater sensitivity to the presence of magnetic fields, including very weak fields. For the 47 individual objects studied, including 6 components of SB2 systems, we do not detect any magnetic fields at greater than the 3 sigma level. The lack of detection in the circular polarization profiles indicates that if strong fields are present on these stars, they must have complex surface topologies. For simple global fields, our detection limits imply upper limits to the fields present of 2-10 Gauss in the best cases. We conclude that HgMn stars lack large-scale magnetic fields, typical for spotted magnetic Ap stars, sufficient to form and sustain the chemical spots observed on HgMn stars. Our study confirms that in addition to magnetically altered atomic diffusion, there exists another differentiation mechanism operating in the atmospheres of late-B main sequence stars which can compositional inhomogeneities on their surfaces.

We recently commissioned the polarimetric upgrade of the HARPS spectrograph at ESO's 3.6-m telescope at La Silla, Chile. The HARPS polarimeter is capable of full Stokes spectropolarimetry with large sensitivity and accuracy, taking advantage of the large spectral resolution and stability of HARPS. In this paper we present the instrument design and its polarimetric performance. The first HARPSpol observations show that it can attain a polarimetric sensitivity of ~10^-5 (after addition of many lines) and that no significant instrumental polarization effects are present.

A small fraction of intermediate-mass main sequence (A and B type) stars have strong, organised magnetic fields. The large majority of such stars, however, show no evidence for magnetic fields, even when observed with very high precision. In this paper we describe a simple model, motivated by qualitatively new observational results, that provides a natural physical explanation for the small fraction of observed magnetic stars.

We have investigated a sample of 28 well-known spectroscopically-identified magnetic Ap/Bp stars, with weak, poorly-determined or previously undetected magnetic fields, with the aim of exploring the weak part of the magnetic field distribution of Ap/Bp stars. Using the MuSiCoS and NARVAL spectropolarimeters we have obtained 282 LSD Stokes V signatures of our 28 sample stars. All stars were detected, showing clearly that when observed with sufficient precision, all firmly classified Ap/Bp stars show detectable surface magnetic fields. To better characterise the surface magnetic field intensities and geometries of the sample, we have inferred the dipolar field intensity and the magnetic obliquity. The distribution of derived dipole strengths for these stars exhibits a plateau at about 1 kG, falling off to larger and smaller field strengths. Remarkably, in this sample of stars selected for their presumably weak magnetic fields, we find only 2 stars for which the derived dipole strength is weaker than 300 G. We interpret this "magnetic threshold" as a critical value necessary for the stability of large-scale magnetic fields, and develop a simple quantitative model that is able to approximately reproduce the observed threshold characteristics. This scenario leads to a natural explanation of the small fraction of intermediate-mass magnetic stars. It may also explain the near-absence of magnetic fields in more massive B and O-type stars.

From observations collected with the ESPaDOnS spectropolarimeter, we report the discovery of magnetic fields at the surface of the mildly accreting classical T Tauri star V2129 Oph. Zeeman signatures are detected, both in photospheric lines and in the emission lines formed at the base of the accretion funnels linking the disc to the protostar, and monitored over the whole rotation cycle of V2129 Oph. We observe that rotational modulation dominates the temporal variations of both unpolarized and circularly polarized line profiles. We reconstruct the large-scale magnetic topology at the surface of V2129 Oph from both sets of Zeeman signatures simultaneously. We find it to be rather complex, with a dominant octupolar component and a weak dipole of strengths 1.2 and 0.35 kG, respectively, both slightly tilted with respect to the rotation axis. The large-scale field is anchored in a pair of 2-kG unipolar radial field spots located at high latitudes and coinciding with cool dark polar spots at photospheric level. This large-scale field geometry is unusually complex compared to those of non-accreting cool active subgiants with moderate rotation rates. As an illustration, we provide a first attempt at modelling the magnetospheric topology and accretion funnels of V2129 Oph using field extrapolation. We find that the magnetosphere of V2129 Oph must extend to about 7R* to ensure that the footpoints of accretion funnels coincide with the high-latitude accretion spots on the stellar surface. It suggests that the stellar magnetic field succeeds in coupling to the accretion disc as far out as the corotation radius, and could possibly explain the slow rotation of V2129 Oph. The magnetospheric geometry we derive produces X-ray coronal fluxes typical of those observed in cTTSs.

Surface brightness maps for the young K0 dwarf AB Doradus are reconstructed from archival data sets for epochs spanning 1988 to 1994. By using the signal-to-noise enhancement technique of Least-Squares Deconvolution, our results show a greatly increased resolution of spot features than obtained in previously published surface brightness reconstructions. These images show that for the exception of epoch 1988.96, the starspot distributions are dominated by a long-lived polar cap, and short-lived low to high latitude features. The fragmented polar cap at epoch 1988.96 could indicate a change in the nature of the dynamo in the star. For the first time we measure differential rotation for epochs with sufficient phase coverage (1992.05, 1993.89, 1994.87). These measurements show variations on a timescale of at least one year, with the strongest surface differential rotation ever measured for AB Dor occurring in 1994.86. In conjunction with previous investigations, our results represent the first long-term analysis of the temporal evolution of differential rotation on active stars.

We have used maximum entropy eclipse mapping to recover images of the visual surface brightness distribution of the primary component of the RS CVn eclipsing binary SV Cam, using high-precision photometry data obtained during three primary eclipses with STIS aboard the Hubble Space Telescope. These were augmented by contemporaneous ground-based photometry secured around the rest of the orbit. The goal of these observations was to determine the filling factor and size distribution of starspots too small to be resolved by Doppler imaging. The information content of the final image and the fit to the data were optimised with respect to various system parameters using the chi^2 landscape method, using an eclipse mapping code that solves for large-scale spot coverage. It is only with the unprecedented photometric precision of the HST data (0.00015 mag) that it is possible to see strong discontinuities at the four contact points in the residuals of the fit to the lightcurve. These features can only be removed from the residual lightcurve by the reduction of the photospheric temperature, to synthesise high unresolvable spot coverage, and the inclusion of a polar spot. We show that this spottedness of the stellar surface can have a significant impact on the determination of the stellar binary parameters and the fit to the lightcurve by reducing the secondary radius from 0.794 +/- 0.009 Rsun to 0.727 +/- 0.009 Rsun. This new technique can also be applied to other binary systems with high precision spectrophotometric observations.

The variation of the specific intensity across the stellar disc is essential input parameter in surface brightness reconstruction techniques such as Doppler imaging, where the relative intensity contributions of different surface elements are important in detecting starspots. We use PHOENIX and ATLAS model atmospheres to model lightcurves derived from high precision (S/N ~ 5000) HST data of the eclipsing binary SV Cam (F9V + K4V), where the variation of specific intensity across the stellar disc will determine the contact points of the binary system lightcurve. For the first time we use chi^2 comparison fits to the first derivative profiles to determine the best-fitting model atmosphere. We show the wavelength dependence of the limb darkening and that the first derivative profile is sensitive to the limb-darkening profile very close to the limb of the primary star. It is concluded that there is only a marginal difference (< 1sigma) between the chi^2 comparison fits of the two model atmospheres to the HST lightcurve at all wavelengths. The usefulness of the second derivative of the light-curve for measuring the sharpness of the primary's limb is investigated, but we find that the data are too noisy to permit a quantitative analysis.

In this paper, size distributions of starspots extrapolated from the case of the Sun, are modelled on the eclipsing binary SV Cam to synthesise images of stellar photospheres with high spot filling factors. These spot distributions pepper the primary's surface with spots, many of which are below the resolution capabilities of eclipse mapping and Doppler imaging techniques. The lightcurves resulting from these modelled distributions are used to determine the limitations of image reconstruction from photometric data. Surface brightness distributions reconstructed from these lightcurves show distinctive spots on the primary star at its quadrature points. It is concluded that two-spot modelling or chi-squared minimisation techniques are more susceptible to spurious structures being generated by systematic errors, arising from incorrect assumptions about photospheric surface brightness, than simple Fourier analysis of the light-curves.

We have used spectrophotometric data from the Hubble Space Telescope to eclipse-map the primary component of the RS CVn binary SV Cam over 9 HST orbits. We find from these observations and the HIPPARCOS parallax that the surface flux in the eclipsed low-latitude region of the primary is about 30% lower than computed from a PHOENIX model atmosphere at the effective temperature that best fits the spectral energy distribution of the eclipsed flux. This can only be accounted for if about a third of the primary's surface is covered with unresolved dark star-spots. Extending this to the full surface of the primary, we find that even taking into account this spot filling factor there is an additional flux deficit on the primary star. This can only be explained if there is a large polar spot on the primary star extending from the pole to latitude 48 degrees.

Self-gravitating protostellar discs are unstable to fragmentation if the gas can cool on a time scale that is short compared to the orbital period. We use a combination of hydrodynamic simulations and N-body orbit integrations to study the long term evolution of a fragmenting disc with an initial mass ratio to the star of M_disc/M_star = 0.1. For a disc which is initially unstable across a range of radii, a combination of collapse and subsequent accretion yields substellar objects with a spectrum of masses extending (for a Solar mass star) up to ~0.01 M_sun. Subsequent gravitational evolution ejects most of the lower mass objects within a few million years, leaving a small number of very massive planets or brown dwarfs in eccentric orbits at moderately small radii. Based on these results, systems such as HD 168443 -- in which the companions are close to or beyond the deuterium burning limit -- appear to be the best candidates to have formed via gravitational instability. If massive substellar companions originate from disc fragmentation, while lower-mass planetary companions originate from core accretion, the metallicity distribution of stars which host massive substellar companions at radii of ~1 au should differ from that of stars with lower mass planetary companions.