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

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

EarthFinder is a NASA Astrophysics Probe mission concept selected for study as input to the 2020 Astrophysics National Academies Decadal Survey. The EarthFinder concept is based on a dramatic shift in our understanding of how PRV measurements should be made. We propose a new paradigm which brings the high precision, high cadence domain of transit photometry as demonstrated by Kepler and TESS to the challenges of PRV measurements at the cm/s level. This new paradigm takes advantage of: 1) broad wavelength coverage from the UV to NIR which is only possible from space to minimize the effects of stellar activity; 2) extremely compact, highly stable, highly efficient spectrometers (R>150,000) which require the diffraction-limited imaging possible only from space over a broad wavelength range; 3) the revolution in laser-based wavelength standards to ensure cm/s precision over many years; 4) a high cadence observing program which minimizes sampling-induced period aliases; 5) exploiting the absolute flux stability from space for continuum normalization for unprecedented line-by-line analysis not possible from the ground; and 6) focusing on the bright stars which will be the targets of future imaging missions so that EarthFinder can use a ~1.5 m telescope.

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.

Statistical analyses from exoplanet surveys around low-mass stars indicate that super-Earth and Neptune-mass planets are more frequent than gas giants around such stars, in agreement with core accretion theory of planet formation. Using precise radial velocities derived from visual and near-infrared spectra, we report the discovery of a giant planet with a minimum mass of 0.46 Jupiter masses in an eccentric 204-day orbit around the very low-mass star GJ 3512. Dynamical models show that the high eccentricity of the orbit is most likely explained from planet-planet interactions. The reported planetary system challenges current formation theories and puts stringent constraints on the accretion and migration rates of planet formation and evolution models, indicating that disc instability may be more efficient in forming planets than previously thought.

Context. Teegarden's Star is the brightest and one of the nearest ultra-cool dwarfs in the solar neighbourhood. For its late spectral type (M7.0V), the star shows relatively little activity and is a prime target for near-infrared radial velocity surveys such as CARMENES. Aims. As part of the CARMENES search for exoplanets around M dwarfs, we obtained more than 200 radial-velocity measurements of Teegarden's Star and analysed them for planetary signals. Methods. We find periodic variability in the radial velocities of Teegarden's Star. We also studied photometric measurements to rule out stellar brightness variations mimicking planetary signals. Results. We find evidence for two planet candidates, each with $1.1M_\oplus$ minimum mass, orbiting at periods of 4.91 and 11.4 d, respectively. No evidence for planetary transits could be found in archival and follow-up photometry. Small photometric variability is suggestive of slow rotation and old age. Conclusions. The two planets are among the lowest-mass planets discovered so far, and they are the first Earth-mass planets around an ultra-cool dwarf for which the masses have been determined using radial velocities.

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.

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.

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.

Brightness inhomogeneities in the stellar photosphere (dark spots or bright regions) affect the measurements of the planetary transmission spectrum. To investigate the star spots of the M dwarf GJ 1214, we conducted a multicolor photometric monitoring from 2012 to 2016. The measured variability shows a periodicity of 125 +- 5 days, which we interpret as the signature of the stellar rotation period. This value overrules previous suggestions of a significantly shorter stellar rotation period. A light curve inversion of the monitoring data yields an estimation of the flux dimming of a permanent spot filling factor not contributing to the photometric variability, a temperature contrast of the spots of about 370 K and persistent active longitudes. The derived surface maps over all five seasons were used to estimate the influence of the star spots on the transmission spectrum of the planet from 400 nm to 2000 nm. The monitoring data presented here do not support a recent interpretation of a measured transmission spectrum of GJ 1214b as to be caused by bright regions in the stellar photosphere. Instead, we list arguments as to why the effect of dark spots likely dominated over bright regions in the period of our monitoring. Furthermore, our photometry proves an increase in variability over at least four years, indicative for a cyclic activity behavior. The age of GJ 1214 is likely between 6 and 10 Gyr. The long-term photometry allows for a correction of unocculted spots. For an active star such as GJ 1214, there remains a degeneracy between occulted spots and the transit parameters used to build the transmission spectrum. This degeneracy can only be broken by high-precision transit photometry resolving the spot crossing signature in the transit light curve.

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 aim to investigate the presence of signatures of magnetic cycles and rotation on a sample of 71 early M-dwarfs from the HADES RV programme using high-resolution time-series spectroscopy of the Ca II H & K and Halpha chromospheric activity indicators, the radial velocity series, the parameters of the cross correlation function and the V-band photometry. We used mainly HARPS-N spectra, acquired over four years, and add HARPS spectra from the public ESO database and ASAS photometry light-curves as support data, extending the baseline of the observations of some stars up to 12 years. We provide log(R'hk) measurements for all the stars in the sample, cycle length measurements for 13 stars, rotation periods for 33 stars and we are able to measure the semi-amplitude of the radial velocity signal induced by rotation in 16 stars. We complement our work with previous results and confirm and refine the previously reported relationships between the mean level of chromospheric emission, measured by the log(R'hk), with the rotation period, and with the measured semi-amplitude of the activity induced radial velocity signal for early M-dwarfs. We searched for a possible relation between the measured rotation periods and the lengths of the magnetic cycle, finding a weak correlation between both quantities. Using previous v sin i measurements we estimated the inclinations of the star's poles to the line of sight for all the stars in the sample, and estimate the range of masses of the planets GJ 3998 b and c (2.5 - 4.9 Mearth and 6.3 - 12.5 Mearth), GJ 625 b (2.82 Mearth), GJ 3942 b (7.1 - 10.0 Mearth) and GJ 15A b (3.1 - 3.3 Mearth), assuming their orbits are coplanar with the stellar rotation.

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.

Short- to mid-term magnetic phenomena on the stellar surface of M-type stars cannot only resemble the effects of planets in radial velocity data, but also may hide them. We analyze 145 spectroscopic HARPS-N observations of GJ 3942 taken over the past five years and additional photometry to disentangle stellar activity effects from genuine Doppler signals as a result of the orbital motion of the star around the common barycenter with its planet. To achieve this, we use the common methods of pre-whitening, and treat the correlated red noise by a first-order moving average term and by Gaussian-process regression following an MCMC analysis. We identify the rotational period of the star at 16.3 days and discover a new super-Earth, GJ 3942 b, with an orbital period of 6.9 days and a minimum mass of 7.1 Me. An additional signal in the periodogram of the residuals is present but we cannot claim it to be related to a second planet with sufficient significance at this point. If confirmed, such planet candidate would have a minimum mass of 6.3 Me and a period of 10.4 days, which might indicate a 3:2 mean-motion resonance with the inner planet.

We report the discovery of a super-Earth orbiting at the inner edge of the habitable zone of the star GJ 625 based on the analysis of the radial-velocity (RV) time series from the HARPS-N spectrograph, consisting in 151 HARPS-N measurements taken over 3.5 yr. GJ 625 b is a planet with a minimum mass M sin $i$ of 2.82 $\pm$ 0.51 M$_{\oplus}$ with an orbital period of 14.628 $\pm$ 0.013 days at a distance of 0.078 AU of its parent star. The host star is the quiet M2 V star GJ 625, located at 6.5 pc from the Sun. We find the presence of a second radial velocity signal in the range 74-85 days that we relate to stellar rotation after analysing the time series of Ca II H\&K and H${\alpha}$ spectroscopic indicators, the variations of the FWHM of the CCF and and the APT2 photometric light curves. We find no evidence linking the short period radial velocity signal to any activity proxy.

The distribution of exoplanets around low-mass stars is still not well understood. Such stars, however, present an excellent opportunity of reaching down to the rocky and habitable planet domains. The number of current detections used for statistical purposes is still quite modest and different surveys, using both photometry and precise radial velocities, are searching for planets around M dwarfs. Our HARPS-N red dwarf exoplanet survey is aimed at the detection of new planets around a sample of 78 selected stars, together with the subsequent characterization of their activity properties. Here we investigate the survey performance and strategy. From 2700 observed spectra, we compare the radial velocity determinations of the HARPS-N DRS pipeline and the HARPS-TERRA code, we calculate the mean activity jitter level, we evaluate the planet detection expectations, and we address the general question of how to define the strategy of spectroscopic surveys in order to be most efficient in the detection of planets. We find that the HARPS-TERRA radial velocities show less scatter and we calculate a mean activity jitter of 2.3 m/s for our sample. For a general radial velocity survey with limited observing time, the number of observations per star is key for the detection efficiency. In the case of an early M-type target sample, we conclude that approximately 50 observations per star with exposure times of 900 s and precisions of about 1 m/s maximizes the number of planet detections.

M dwarfs are prime targets for planet search programs, particularly of those focused on the detection and characterization of rocky planets in the habitable zone. Understanding their magnetic activity is important because it affects our ability to detect small planets, and it plays a key role in the characterization of the stellar environment. We analyze observations of the Ca II H&K and H{\alpha} lines as diagnostics of chromospheric activity for low-activity early-type M dwarfs. We analyze the time series of spectra of 71 early-type M dwarfs collected for the HADES project for planet search purposes. The HARPS-N spectra provide simultaneously the H&K doublet and the H{\alpha} line. We develop a reduction scheme able to correct the HARPS-N spectra for instrumental and atmospheric effects, and to provide flux-calibrated spectra in units of flux at the stellar surface. The H&K and H{\alpha} fluxes are compared with each other, and their variability is analyzed. We find that the H and K flux excesses are strongly correlated with each other, while the H{\alpha} flux excess is generally less correlated with the H&K doublet. We also find that H{\alpha} emission does not increase monotonically with the H&K line flux, showing some absorption before being filled in by chromospheric emission when H&K activity increases. Analyzing the time variability of the emission fluxes, we derive a tentative estimate of the rotation period (of the order of a few tens of days) for some of the program stars, and the typical lifetime of chromospheric active regions (a few stellar rotations). Our results are in good agreement with previous studies. In particular, we find evidence that the chromospheres of early-type M dwarfs could be characterized by different filaments coverage, affecting the formation mechanism of the H{\alpha} line. We also show that chromospheric structure is likely related to spectral type.

Context. M dwarfs are considered ideal targets for Doppler radial velocity searches. Nonetheless, the statistics of frequency of low-mass planets hosted by low mass stars remains poorly constrained. Aims. Our M-dwarf radial velocity monitoring with HARPS-N can provide a major contribution to the widening of the current statistics through the in-depth analysis of accurate radial velocity observations in a narrow range of spectral sub-types (79 stars, between dM0 to dM3). Spectral accuracy will enable us to reach the precision needed to detect small planets with a few earth masses. Our survey will bring a contribute to the surveys devoted to the search for planets around M-dwarfs, mainly focused on the M-dwarf population of the northern emisphere, for which we will provide an estimate of the planet occurence. Methods. We present here a long duration radial velocity monitoring of the M1 dwarf star GJ 3998 with HARPS-N to identify periodic signals in the data. Almost simultaneous photometric observations were carried out within the APACHE and EXORAP programs to characterize the stellar activity and to distinguish from the periodic signals those due to activity and to the presence of planetary companions. Results. The radial velocities have a dispersion in excess of their internal errors due to at least four superimposed signals, with periods of 30.7, 13.7, 42.5 and 2.65 days. The analysis of spectral indices based on Ca II H & K and Halpha lines demonstrates that the periods of 30.7 and 42.5 days are due to chromospheric inhomogeneities modulated by stellar rotation and differential rotation. The shorter periods of 13.74 +/- 0.02 d and 2.6498 +/- 0.0008 d are well explained with the presence of two planets, with minimum masses of 6.26 +/- 0.79 MEarth and 2.47 +/- 0.27 MEarth and distances of 0.089 AU and 0.029 AU from the host, respectively.

Stellar activity patterns are responsible for jitter effects that are observed at different timescales and amplitudes. These effects are currently in the focus of many exoplanet search projects, since the lack of a well-defined characterization and correction strategy hampers the detection of the signals associated with small exoplanets. Accurate simulations of the stellar photosphere can provide synthetic time series data. These may help to investigate the relation between activity jitter and stellar parameters when considering different active region patterns. Moreover, jitters can be analysed at different wavelength scales in order to design strategies to remove or minimize them. In this work we present the StarSim tool, which is based on a model for a spotted rotating photosphere built from the integration of the spectral contribution of a fine grid of surface elements. The model includes all significant effects affecting the flux intensities and the wavelength of spectral features produced by active regions and planets. A specific application for the characterization and modelling of the spectral signature of active regions is considered, showing that the chromatic effects of faculae are dominant for low temperature contrasts of spots. Synthetic time series are modelled for HD 189733. Our algorithm reproduces both the photometry and the RVs to good precision, generally better than the studies published to date. We evaluate the RV signature of the activity in HD 189733 by exploring a grid of solutions from the photometry. We find that the use of RV data in the inverse problem could break degeneracies and allow for a better determination of some stellar and activity parameters. In addition, the effects of spots are studied for a set of simulated transit photometry, showing that these can introduce variations which are very similar to the signal of an atmosphere dominated by dust.

Certain instrumental effects and data reduction anomalies introduce systematic errors in photometric time-series. Detrending algorithms such as the Trend Filtering Algorithm (TFA) (Kov\'{a}cs et al. 2004) have played a key role in minimizing the effects caused by these systematics. Here we present the results obtained after applying the TFA, Savitszky-Golay (Savitzky & Golay 1964) detrending algorithms and the Box Least Square phase folding algorithm (Kov\'{a}cs et al. 2002) to the TFRM-PSES data (Fors et al. 2013). Tests performed on this data show that by applying these two filtering methods together, the photometric RMS is on average improved by a factor of 3-4, with better efficiency towards brighter magnitudes, while applying TFA alone yields an improvement of a factor 1-2. As a result of this improvement, we are able to detect and analyze a large number of stars per TFRM-PSES field which present some kind of variability. Also, after porting these algorithms to Python and parallelizing them, we have improved, even for large data samples, the computing performance of the overall detrending+BLS algorithm by a factor of $\sim$10 with respect to Kov\'{a}cs et al. (2004).