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.

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.

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.

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.

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.

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

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

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

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

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

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

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

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

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

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.

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

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

We present the analysis of four M dwarf stars -plus one M giant that seeped past our selection criteria- observed in Cycle 3 of Kepler Guest Observer program (GO3) in a search for intrinsic pulsations. Stellar oscillations in M dwarfs were theoretically predicted by Rodr\'iguez-L\'opez et al. (2012) to be in the range ~20-40 min and ~4-8 h, depending on the age and the excitation mechanism. We requested Kepler short cadence observations to have an adequate sampling of the oscillations. The targets were chosen on the basis of detectable rotation in the initial Kepler results, biasing towards youth.The analysis reveals no oscillations attributable to pulsations at a detection limit of several parts per million, showing that either the driving mechanisms are not efficient in developing the oscillations to observable amplitudes, or that if pulsations are driven, the amplitudes are very low. The size of the sample, and the possibility that the instability strip is not pure, allowing the coexistence of pulsators and non-pulsators, prevent us from deriving definite conclusions. Inmediate plans include more M dwarfs photometric observations of similar precision with Kepler K2 mission and spectroscopic searches already underway within the Cool Tiny Beats Project (Anglada-Escud\'e et al. 2014, Berdi\~nas et al. 2014) with the high-resolution spectrographs HARPS and HARPS-N.

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.

The overstability of the fundamental radial mode in M dwarf models was theoretically predicted by Rodr\'iguez-L\'opez et al. (2012). The periods were found to be in the ranges ~25-40 min and ~4-8 h, depending on stellar age and excitation mechanism. We have extended our initial M dwarf model grid in mass, metallicity, and mixing length parameter. We have also considered models with boundary conditions from PHOENIX NextGen atmospheres to test their influence on the pulsation spectra. We find instability of non-radial modes with radial orders up to k=3, degree l=0-3, including p and g modes, with the period range extending from 20 min up to 11 h. Furthermore, we find theoretical evidence of the potential of M dwarfs as solar-like oscillators.

(Abridged) The origin of hot subdwarf O-type stars (sdOs) remains unclear since their discovery in 1947. Among others, a post-Asymptotic Giant Branch (post-AGB) origin is possible for a fraction of sdOs. We are involved in a comprehensive ongoing study to search for and to analyze planetary nebulae (PNe) around sdOs with the aim of establishing the fraction and properties of sdOs with a post-AGB origin. We use deep Halpha and [OIII] images of sdOs to detect nebular emission and intermediate resolution, long-slit optical spectroscopy of the detected nebulae and their sdO central stars. These data are complemented with other observations for further analysis of the detected nebulae. We report the detection of an extremely faint, complex PN around 2MASS J19310888+4324577 (2M1931+4324), a star classified as sdO in a binary system. The PN shows a bipolar and an elliptical shell, whose major axes are oriented perpendicular to each other, and high-excitation structures outside the two shells. WISE archive images show faint, extended emission at 12 and 22 microns in the inner nebular regions. The internal nebular kinematics is consistent with a bipolar and a cylindrical/ellipsoidal shell, in both cases with the main axis mainly perpendicular to the line of sight. The nebular spectrum only exhibits Halpha, Hbeta and [OIII]4959,5007 emission lines, but suggests a very low-excitation ([OIII]/Hbeta = 1.5), in strong contrast with the absence of low-excitation emission lines. The spectrum of 2M1931+4324 presents narrow, ionized helium absorptions that confirm the previous sdO classification and suggest an effective temperature >= 60000 K. The binary nature of 2M1931+4324, its association with a complex PN, and several properties of the system provide strong support for the idea that binary central stars are a crucial ingredient in the formation of complex PNe.

In this work, we present the first results of an ongoing survey to search for planetary nebulae (PNe) around hot subdwarf stars (sdOs). Deep images and intermediate-resolution long-slit spectra of RWT 152, the only confirmed PN+sdO system in the northern hemisphere, as well as preliminary results for other sdO+PN candidate are presented.

We present the results of the first theoretical non-radial non-adiabatic pulsational study of M dwarf stellar models with masses in the range 0.1 to 0.5M_solar. We find the fundamental radial mode to be unstable due to an \epsilon mechanism caused by deuterium (D-) burning for the young 0.1 and 0.2M_solar models, by non-equilibrium He^3 burning for the 0.2 and 0.25M_solar models of 10^4Myr, and by a flux blocking mechanism for the partially convective 0.4 and 0.5M_solar models once they reach the age of 500Myr. The periods of the overstable modes excited by the D-burning are in the range 4.2 to 5.2h for the 0.1M_solar models and is of order 8.4h for the 0.2M_solar models. The periods of the modes excited by He^3 burning and flux blocking are in the range 23 to 40min. The more massive and oldest models are more promising for the observational detection of pulsations, as their ratio of instability e-folding time to age is more favourable.

Hot subdwarf stars are faint, blue objects, and are the main contributors to the far-UV excess observed in elliptical galaxies. They offer an excellent laboratory to study close and wide binary systems, and to scrutinize their interiors through asteroseismology, as some of them undergo stellar oscillations. However, their origins are still uncertain, and increasing the number of detections is crucial to undertake statistical studies. In this work, we aim at defining a strategy to find new, uncatalogued hot subdwarfs. Making use of Virtual Observatory tools we thoroughly search stellar catalogues to retrieve multi-colour photometry and astrometric information of a known sample of blue objects, including hot subdwarfs, white dwarfs, cataclysmic variables and main sequence OB stars. We define a procedure to discriminate among these spectral classes, particularly designed to obtain a hot subdwarf sample with a low contamination factor. In order to check the validity of the method, this procedure is then applied to two test sky regions: the Kepler FoV and to a test region of around (RA:225, DEC:5) deg. As a result, we obtained 38 hot subdwarf candidates, 23 of which had already a spectral classification. We have acquired spectroscopy for three other targets, and four additional ones have an available SDSS spectrum, which we used to determine their spectral type. A temperature estimate is provided for the candidates based on their spectral energy distribution, considering two-atmospheres fit for objects with clear infrared excess. Eventually, out of 30 candidates with spectral classification, 26 objects were confirmed to be hot subdwarfs, yielding a contamination factor of only 13%. The high rate of success demonstrates the validity of the proposed strategy to find new uncatalogued hot subdwarfs. An application of this method to the entire sky will be presented in a forthcoming work.

We present results from the final six months of a survey to search for pulsations in white dwarfs and hot subdwarf stars with the Kepler spacecraft. Spectroscopic observations are used to separate the objects into accurate classes, and we explore the physical parameters of the subdwarf B (sdB) stars and white dwarfs in the sample. From the Kepler photometry and our spectroscopic data, we find that the sample contains 5 new pulsators of the V1093 Her type, one AM CVn type cataclysmic variable, and a number of other binary systems. This completes the survey for compact pulsators with Kepler. No V361 Hya type of short-period pulsating sdB stars were found in this half, leaving us with a total of one single multiperiodic V361 Hya and 13 V1093 Her pulsators for the full survey. Except for the sdB pulsators, no other clearly pulsating hot subdwarfs or white dwarfs were found, although a few low-amplitude candidates still remain. The most interesting targets discovered in this survey will be observed throughout the remainder of the Kepler Mission, providing the most long-term photometric datasets ever made on such compact, evolved stars. Asteroseismic investigations of these datasets will be invaluable in revealing the interior structure of these stars, and will boost our understanding of their evolutionary history.

We present results from the first two quarters of a survey to search for pulsations in compact stellar objects with the Kepler spacecraft. The survey sample and the various methods applied in its compilation are described, and spectroscopic observations are presented to separate the objects into accurate classes. From the Kepler photometry we clearly identify nine compact pulsators, and a number of interesting binary stars. Of the pulsators, one shows the strong, rapid pulsations typical for a V361 Hya type sdB variable (sdBV), seven show long-period pulsations characteristic of V1093 Her type sdBVs, and one shows low-amplitude pulsations with both short and long periods. We derive effective temperatures and surface gravities for all the subdwarf B stars in the sample and demonstrate that below the boundary region where hybrid sdB pulsators are found, all our targets are pulsating. For the stars hotter than this boundary temperature a low fraction of strong pulsators (<10 per cent) is confirmed. Interestingly, the short-period pulsator also shows a low-amplitude mode in the long-period region, and several of the V1093 Her pulsators show low amplitude modes in the short-period region, indicating that hybrid behaviour may be common in these stars, also outside the boundary temperature region where hybrid pulsators have hitherto been found.

In the last years, space missions such as COROT, Kepler or MOST have provided very accurate photometric observational data. In the particular case of $\delta$ Scuti stars, the observed frequency spectra have hundreds (if not thousands) of modes and a clear amplitude distribution. In this work we present new techniques for modelling these observations and the results obtained. We searched for regular patterns in the observational data, which yields something resembling the large separation. This allows to reduce the possible positions of the star in the HR diagram, yielding a value of the mean density with an accuracy never reached before for isolated stars of this type. Finally, we answer whether a $\delta$ Scuti star is stable despite all of the observed frequencies are simultaneously excited.

The recent analyses of the light curves provided by CoRoT have revealed pulsation spectra of unprecedented richness and precision, in particular, thousands of pulsating modes, and a clear distribution of amplitudes with frequency. In the community, some scientists have started doubting about the validity of the classical tools to analyze these very accurate light curves. This work provides the asteroseismic community with answers to this question showing that (1) it is physically possible for a star to excite at a time and with the observed amplitudes such a large number of modes; and (2) that the kinetic energy accumulated in all those modes does not destroy the equilibrium of the star. Consequently, mathematical tools presently applied in the analyses of light curves can a priori be trusted. This conclusion is even more important now, when a large amount of space data coming from Kepler are currently being analyzed. The power spectrum of different stellar cases, and the non-adiabatic code GraCo have been used to estimate the upper limit of the energy per second required to excite all the observed modes, and their total kinetic energy. A necessary previous step for this study is to infer the relative radial pulsational amplitude from the observed photometric amplitude, scaling our linear pulsational solutions to absolute values. The derived upper limits for the required pulsational energy were compared with 1) the luminosity of the star; and 2) the gravitational energy. We obtained that both upper energy limits are orders of magnitude smaller.

We present the first description of mode trapping for sdO models. Mode trapping of gravity modes caused by the He/H chemical transition is found for a particular model, providing a selection effect for high radial order trapped modes. Low- and intermediate-radial order {\em p}-modes (mixed modes with a majority of nodes in the P-mode region) are found to be trapped by the C-O/He transition, but with no significant effects on the driving. This region seems to have also a subtle effect on the trapping of low radial order {\em g}-modes (mixed modes with a majority of nodes in the G-mode region), but again with no effect on the driving. We found that for mode trapping to have an influence on the driving of sdO modes (1) the mode should be trapped in a way that the amplitude of the eigenfunctions is lower in a damping region and (2) in this damping region significant energy interchange has to be produced.

We present new photometric and spectroscopic observations for 2M 1533+3759 (= NSVS 07826147). It has an orbital period of 0.16177042 day, significantly longer than the 2.3--3.0 hour periods of the other known eclipsing sdB+dM systems. Spectroscopic analysis of the hot primary yields Teff = 29230 +/- 125 K, log g = 5.58 +/- 0.03 and log N(He)/N(H) = -2.37 +/- 0.05. The sdB velocity amplitude is K1 = 71.1 +/- 1.0 km/s. The only detectable light contribution from the secondary is due to the surprisingly strong reflection effect. Light curve modeling produced several solutions corresponding to different values of the system mass ratio, q(M2/M1), but only one is consistent with a core helium burning star, q=0.301. The orbital inclination is 86.6 degree. The sdB primary mass is M1 = 0.376 +/- 0.055 Msun and its radius is R1 = 0.166 +/- 0.007 Rsun. 2M1533+3759 joins PG0911+456 (and possibly also HS2333+3927) in having an unusually low mass for an sdB star. SdB stars with masses significantly lower than the canonical value of 0.48 Msun, down to as low as 0.30 Msun, were theoretically predicted by Han et al. (2002, 2003), but observational evidence has only recently begun to confirm the existence of such stars. The existence of core helium burning stars with masses lower than 0.40--0.43 Msun implies that at least some sdB progenitors have initial main sequence masses of 1.8--2.0 Msun or more, i.e. they are at least main sequence A stars. The secondary is a main sequence M5 star.

We have explored the possibility of driving pulsation modes in models of sdO stars in which the effects of element diffusion, gravitational settling and radiative levitation have been neglected so that the distribution of iron-peak elements remains uniform throughout the evolution. The stability of these models was determined using a non-adiabatic oscillations code. We analysed 27 sdO models from 16 different evolutionary sequences and discovered the first ever sdO models capable of driving high-radial order g-modes. In one model, the driving is by a classical kappa-mechanism due to the opacity bump from iron-peak elements at temperature ~200,000 K. In a second model, the driving result from the combined action of kappa-mechanisms operating in three distinct regions of the star: (i) a carbon-oxygen partial ionization zone at temperature ~2 10^6 K, (ii) a deeper region at temperature ~2 10^7 K, which we attribute to ionization of argon, and (iii) at the transition from radiative to conductive opacity in the core of the star.

HD172189 is a spectroscopic eclipsing binary system with a rapidly-rotating pulsating delta Scuti component. It is also a member of the open cluster IC4756. These combined characteristics make it an excellent laboratory for asteroseismic studies. To date, HD172189 has been analysed in detail photometrically but not spectroscopically. For this reason we have compiled a set of spectroscopic data to determine the absolute and atmospheric parameters of the components. We determined the radial velocities (RV) of both components using four different techniques. We disentangled the binary spectra using KOREL, and performed the first abundance analysis on both disentangled spectra. By combining the spectroscopic results and the photometric data, we obtained the component masses, 1.8 and 1.7 solar masses, and radii, 4.0 and 2.4 solar radii, for inclination i = 73.2 degrees, eccentricity e = 0.28, and orbital period 5.70198 days. Effective temperatures of 7600 K and 8100 K were also determined. The measured vsini are 78 and 74 km/s, respectively, giving rotational periods of 2.50 and 1.55 days for the components. The abundance analysis shows [Fe/H] = -0.28 for the primary (pulsating) star, consistent with observations of IC4756. We also present an assessment of the different analysis techniques used to obtain the RVs and the global parameters.

A spectroscopic analysis of SDSS J160043.6+074802.9, a binary system containing a pulsating subdwarf-O (sdO) star with a late-type companion, yields Teff = 70 000 +/- 5000 K and log g = 5.25 +/- 0.30, together with a most likely type of K3V for the secondary star. We compare our results with atmospheric parameters derived by Fontaine et al. (2008) and in the context of existing evolution models for sdO stars. New and more extensive photometry is also presented which recovers most, but not all, frequencies found in an earlier paper. It therefore seems probable that some pulsation modes have variable amplitudes. A non-adiabatic pulsation analysis of uniform metallicity sdO models show those having log g > 5.3 to be more likely to be unstable and capable of driving pulsation in the observed frequency range.