No Arabic abstract
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.
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 Oro. 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.
Stellar activity poses one of the main obstacles for the detection and characterisation of small exoplanets around cool stars, as it can induce radial velocity (RV) signals that can hide or mimic the presence of planetary companions. Several indicators of stellar activity are routinely used to identify activity-related signals in RVs, but not all indicators trace exactly the same activity effects, nor are any of them always effective in all stars. We evaluate the performance of a set of spectroscopic activity indicators for M dwarf stars with different masses and activity levels with the aim of finding a relation between the indicators and stellar properties. In a sample of 98 M dwarfs observed with CARMENES, we analyse the temporal behaviour of RVs and nine spectroscopic activity indicators: cross-correlation function (CCF) full-width-at-half-maximum (FWHM), contrast, and bisector inverse slope (BIS), chromatic index (CRX), differential line width (dLW), and indices of the chromospheric lines H$alpha$ and calcium infrared triplet. A total of 56 stars of the initial sample show periodic signals related to activity in at least one of these ten parameters. RV is the parameter for which most of the targets show an activity-related signal. CRX and BIS are effective activity tracers for the most active stars in the sample, especially stars with a relatively high mass, while for less active stars, chromospheric lines perform best. FWHM and dLW show a similar behaviour in all mass and activity regimes, with the highest number of activity detections in the low-mass, high-activity regime. Most of the targets for which we cannot identify any activity-related signals are stars at the low-mass end of the sample. These low-mass stars also show the lowest RV scatter, which indicates that ultracool M dwarfs could be better candidates for planet searches than earlier types, which show larger RV jitter.
We announce the discovery of two planets orbiting the M dwarfs GJ 251 ($0.360pm0.015$ M$_odot$) and HD 238090 ($0.578pm0.021$ M$_odot$) based on CARMENES radial velocity (RV) data. In addition, we independently confirm with CARMENES data the existence of Lalande 21185 b, a planet that has recently been discovered with the SOPHIE spectrograph. All three planets belong to the class of warm or temperate super-Earths and share similar properties. The orbital periods are 14.24 d, 13.67 d, and 12.95 d and the minimum masses are $4.0pm0.4$ $M_oplus$, $6.9pm0.9$ $M_oplus$, and $2.7pm0.3$ $M_oplus$ for GJ 251 b, HD 238090 b, and Lalande 21185 b, respectively. Based on the orbital and stellar properties, we estimate equilibrium temperatures of $351.0pm1.4$ K for GJ 251 b, $469.6pm2.6$ K for HD 238090 b, and $370.1pm6.8$ K for Lalande 21185 b. For the latter we resolve the daily aliases that were present in the SOPHIE data and that hindered an unambiguous determination of the orbital period. We find no significant signals in any of our spectral activity indicators at the planetary periods. The RV observations were accompanied by contemporaneous photometric observations. We derive stellar rotation periods of $122.1pm2.2$ d and $96.7pm3.7$ d for GJ 251 and HD 238090, respectively. The RV data of all three stars exhibit significant signals at the rotational period or its first harmonic. For GJ 251 and Lalande 21185, we also find long-period signals around 600 d, and 2900 d, respectively, which we tentatively attribute to long-term magnetic cycles. We apply a Bayesian approach to carefully model the Keplerian signals simultaneously with the stellar activity using Gaussian process regression models and extensively search for additional significant planetary signals hidden behind the stellar activity.
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.
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.