The 4-year light curves of 156,717 stars observed with NASA's Kepler mission are analyzed using the AutoRegressive Planet Search (ARPS) methodology described by Caceres et al. (2019). The three stages of processing are: maximum likelihood ARIMA modeling of the light curves to reduce stellar brightness variations; constructing the Transit Comb Filter periodogram to identify transit-like periodic dips in the ARIMA residuals; Random Forest classification trained on Kepler Team confirmed planets using several dozen features from the analysis. Orbital periods between 0.2 and 100 days are examined. The result is a recovery of 76% of confirmed planets, 97% when period and transit depth constraints are added. The classifier is then applied to the full Kepler dataset; 1,004 previously noticed and 97 new stars have light curve criteria consistent with the confirmed planets, after subjective vetting removes clear False Alarms and False Positive cases. The 97 Kepler ARPS Candidate Transits mostly have periods $P<10$ days; many are UltraShort Period hot planets with radii $<1$% of the host star. Extensive tabular and graphical output from the ARPS time series analysis is provided to assist in other research relating to the Kepler sample.

The detection of periodic signals from transiting exoplanets is often impeded by extraneous aperiodic photometric variability, either intrinsic to the star or arising from the measurement process. Frequently, these variations are autocorrelated wherein later flux values are correlated with previous ones. In this work, we present the methodology of the Autoregessive Planet Search (ARPS) project which uses Autoregressive Integrated Moving Average (ARIMA) and related statistical models that treat a wide variety of stochastic processes, as well as nonstationarity, to improve detection of new planetary transits. Providing a time series is evenly spaced or can be placed on an evenly spaced grid with missing values, these low-dimensional parametric models can prove very effective. We introduce a planet-search algorithm to detect periodic transits in the residuals after the application of ARIMA models. Our matched-filter algorithm, the Transit Comb Filter (TCF), is closely related to the traditional Box-fitting Least Squares and provides an analogous periodogram. Finally, if a previously identified or simulated sample of planets is available, selected scalar features from different stages of the analysis -- the original light curves, ARIMA fits, TCF periodograms, and folded light curves -- can be collectively used with a multivariate classifier to identify promising candidates while efficiently rejecting false alarms. We use Random Forests for this task, in conjunction with Receiver Operating Characteristic (ROC) curves, to define discovery criteria for new, high fidelity planetary candidates. The ARPS methodology can be applied to both evenly spaced satellite light curves and densely cadenced ground-based photometric surveys.

The study of accurate methods to estimate the distribution of stellar rotational velocities is important for understanding many aspects of stellar evolution. From such observations we obtain the projected rotational speed v sin(i) in order to recover the true distribution of the rotational velocity. To that end, we need to solve a difficult inverse problem that can be posed as a Fredholm integral of the first kind. n this work we have used a novel approach based on Maximum likelihood (ML) estimation to obtain an approximation of the true rotational velocity probability density function expressed as a sum of known distribution families. In our proposal, the measurements have been treated as random variables drawn from the projected rotational velocity probability density function. We analyzed the case of Maxwellian sum approximation, where we estimated the parameters that define the sum of distributions. The performance of the proposed method is analyzed using Monte Carlo simulations considering two theoretical cases for the probability density function of the true rotational stellar velocities: i) an unimodal Maxwellian probability density distribution and ii) a bimodal Maxwellian probability density distribution. The results show that the proposed method yielded more accurate estimates in comparison with the Tikhonov regularization method, especially for small sample length N=50. Our proposal was evaluated using real data from three sets of measurements, and our findings were validated using three statistical tests. The ML approach with Maxwellian sum approximation is a accurate method to deconvolve the rotational velocity probability density function, even when the sample length is small (N= 50)

Knowing the distribution of stellar rotational velocities is essential for the understanding stellar evolution. Because we measure the projected rotational speed vsini, we need to solve an ill-posed problem given by a Fredholm integral of the first kind to recover the true rotational velocity distribution. After discretization of the Fredholm integral, we apply the Tikhonov regularization method to obtain directly the probability distribution function for stellar rotational velocities. We propose a simple and straightforward procedure to determine the Tikhonov parameter. We applied Monte Carlo simulations to prove that Tikhonov method is a consistent estimator and asymptotically unbiased. This method is applied to a sample of cluster stars. We obtain confidences intervals using bootstrap method. Our results are in good agreement with the one obtained using the Lucy method, in recovering the probability density distribution of rotational velocities. Furthermore, Lucy estimation lies inside our confidence interval. Tikhonov regularization is a very robust method that deconvolve the rotational velocity probability density function from a sample of vsini data straightforward without needing any convergence criteria.

To better understand the evolution of stars in binary systems as well as to constrain the formation of binary stars, it is important to know the binary mass-ratio distribution. However, in most cases, i.e. for single-lined spectroscopic binaries, the mass ratio cannot be measured directly but only derived as the convolution of a function that depends on the mass ratio and the unknown inclination angle of the orbit on the plane of the sky. We extend our previous method to deconvolve this inverse problem (Cure et al. 2014), i.e., we obtain as an integral the cumulative distribution function (CDF) for the mass ratio distribution. After a suitable transformation of variables it turns out that this problem is the same as the one for rotational velocities $v \sin i$, allowing a close analytic formulation for the CDF. We then apply our method to two real datasets: a sample of Am stars binary systems, and a sample of massive spectroscopic binaries in the Cyg OB2 Association.} {We are able to reproduce the previous results of Boffin (2010) for the sample of Am stars, while we show that the mass ratio distribution of massive stars shows an excess of small mass ratio systems, contrarily to what was claimed by Kobulnicky et al. (2014). Our method proves very robust and deconvolves the distribution from a sample in just a single step.

Rotational speed is an important physical parameter of stars and knowing the distribution of stellar rotational velocities is essential for the understanding stellar evolution. However, it cannot be measured directly but the convolution of the rotational speed and the sine of the inclination angle, $v \sin i$. We developed a method to deconvolve this inverse problem and obtain the cumulative distribution function (CDF) for stellar rotational velocities extending the work of Chandrasekhar & M\"unch (1950). This method is applied a) to theoretical synthetic data recovering the original velocity distribution with very small error; b) to a sample of about 12.000 field main--sequence stars, corroborating that the velocity distribution function is non--Maxwellian, but is better described by distributions based on the concept of maximum entropy, such as Tsallis or Kaniadakis distribution functions. This is a very robust and novel method that deconvolve the rotational velocity cumulative distribution function from a sample of $v \sin i$ data in just one single step without needing any convergence criteria.