Kepler Presearch Data Conditioning I - Architecture and Algorithms for Error Correction in Kepler Light Curves

Kepler provides light curves of 156,000 stars with unprecedented precision. However, the raw data as they come from the spacecraft contain significant systematic and stochastic errors. These errors, which include discontinuities, systematic trends, and outliers, obscure the astrophysical signals in the light curves. To correct these errors is the task of the Presearch Data Conditioning (PDC) module of the Kepler data analysis pipeline. The original version of PDC in Kepler did not meet the extremely high performance requirements for the detection of miniscule planet transits or highly accurate analysis of stellar activity and rotation. One particular deficiency was that astrophysical features were often removed as a side-effect to removal of errors. In this paper we introduce the completely new and significantly improved version of PDC which was implemented in Kepler SOC 8.0. This new PDC version, which utilizes a Bayesian approach for removal of systematics, reliably corrects errors in the light curves while at the same time preserving planet transits and other astrophysically interesting signals. We describe the architecture and the algorithms of this new PDC module, show typical errors encountered in Kepler data, and illustrate the corrections using real light curve examples.Comment: Submitted to PASP. Also see companion paper "Kepler Presearch Data Conditioning II - A Bayesian Approach to Systematic Error Correction" by Jeff C. Smith et a

Initial Characteristics of Kepler Short Cadence Data

The Kepler Mission offers two options for observations -- either Long Cadence (LC) used for the bulk of core mission science, or Short Cadence (SC) which is used for applications such as asteroseismology of solar-like stars and transit timing measurements of exoplanets where the 1-minute sampling is critical. We discuss the characteristics of SC data obtained in the 33.5-day long Quarter 1 (Q1) observations with Kepler which completed on 15 June 2009. The truly excellent time series precisions are nearly Poisson limited at 11th magnitude providing per-point measurement errors of 200 parts-per-million per minute. For extremely saturated stars near 7th magnitude precisions of 40 ppm are reached, while for background limited measurements at 17th magnitude precisions of 7 mmag are maintained. We note the presence of two additive artifacts, one that generates regularly spaced peaks in frequency, and one that involves additive offsets in the time domain inversely proportional to stellar brightness. The difference between LC and SC sampling is illustrated for transit observations of TrES-2.Comment: 5 pages, 4 figures, ApJ Letters in pres

Kepler Presearch Data Conditioning II - A Bayesian Approach to Systematic Error Correction

With the unprecedented photometric precision of the Kepler Spacecraft, significant systematic and stochastic errors on transit signal levels are observable in the Kepler photometric data. These errors, which include discontinuities, outliers, systematic trends and other instrumental signatures, obscure astrophysical signals. The Presearch Data Conditioning (PDC) module of the Kepler data analysis pipeline tries to remove these errors while preserving planet transits and other astrophysically interesting signals. The completely new noise and stellar variability regime observed in Kepler data poses a significant problem to standard cotrending methods such as SYSREM and TFA. Variable stars are often of particular astrophysical interest so the preservation of their signals is of significant importance to the astrophysical community. We present a Bayesian Maximum A Posteriori (MAP) approach where a subset of highly correlated and quiet stars is used to generate a cotrending basis vector set which is in turn used to establish a range of "reasonable" robust fit parameters. These robust fit parameters are then used to generate a Bayesian Prior and a Bayesian Posterior Probability Distribution Function (PDF) which when maximized finds the best fit that simultaneously removes systematic effects while reducing the signal distortion and noise injection which commonly afflicts simple least-squares (LS) fitting. A numerical and empirical approach is taken where the Bayesian Prior PDFs are generated from fits to the light curve distributions themselves.Comment: 43 pages, 21 figures, Submitted for publication in PASP. Also see companion paper "Kepler Presearch Data Conditioning I - Architecture and Algorithms for Error Correction in Kepler Light Curves" by Martin C. Stumpe, et a

Validation of Twelve Small Kepler Transiting Planets in the Habitable Zone

We present an investigation of twelve candidate transiting planets from Kepler with orbital periods ranging from 34 to 207 days, selected from initial indications that they are small and potentially in the habitable zone (HZ) of their parent stars. Few of these objects are known. The expected Doppler signals are too small to confirm them by demonstrating that their masses are in the planetary regime. Here we verify their planetary nature by validating them statistically using the BLENDER technique, which simulates large numbers of false positives and compares the resulting light curves with the Kepler photometry. This analysis was supplemented with new follow-up observations (high-resolution optical and near-infrared spectroscopy, adaptive optics imaging, and speckle interferometry), as well as an analysis of the flux centroids. For eleven of them (KOI-0571.05, 1422.04, 1422.05, 2529.02, 3255.01, 3284.01, 4005.01, 4087.01, 4622.01, 4742.01, and 4745.01) we show that the likelihood they are true planets is far greater than that of a false positive, to a confidence level of 99.73% (3 sigma) or higher. For KOI-4427.01 the confidence level is about 99.2% (2.6 sigma). With our accurate characterization of the GKM host stars, the derived planetary radii range from 1.1 to 2.7 R_Earth. All twelve objects are confirmed to be in the HZ, and nine are small enough to be rocky. Excluding three of them that have been previously validated by others, our study doubles the number of known rocky planets in the HZ. KOI-3284.01 (Kepler-438b) and KOI-4742.01 (Kepler-442b) are the planets most similar to the Earth discovered to date when considering their size and incident flux jointly.Comment: 27 pages in emulateapj format, including tables and figures. To appear in The Astrophysical Journa

Preparation of Kepler lightcurves for asteroseismic analyses

The Kepler mission is providing photometric data of exquisite quality for the asteroseismic study of different classes of pulsating stars. These analyses place particular demands on the pre-processing of the data, over a range of timescales from minutes to months. Here, we describe processing procedures developed by the Kepler Asteroseismic Science Consortium (KASC) to prepare light curves that are optimized for the asteroseismic study of solar-like oscillating stars in which outliers, jumps and drifts are corrected.Comment: Accepted for publication in MNRAS. 5 pages, 2 figure

Overview of the Kepler Science Processing Pipeline

The Kepler Mission Science Operations Center (SOC) performs several critical functions including managing the ~156,000 target stars, associated target tables, science data compression tables and parameters, as well as processing the raw photometric data downlinked from the spacecraft each month. The raw data are first calibrated at the pixel level to correct for bias, smear induced by a shutterless readout, and other detector and electronic effects. A background sky flux is estimated from ~4500 pixels on each of the 84 CCD readout channels, and simple aperture photometry is performed on an optimal aperture for each star. Ancillary engineering data and diagnostic information extracted from the science data are used to remove systematic errors in the flux time series that are correlated with these data prior to searching for signatures of transiting planets with a wavelet-based, adaptive matched filter. Stars with signatures exceeding 7.1 sigma are subjected to a suite of statistical tests including an examination of each star's centroid motion to reject false positives caused by background eclipsing binaries. Physical parameters for each planetary candidate are fitted to the transit signature, and signatures of additional transiting planets are sought in the residual light curve. The pipeline is operational, finding planetary signatures and providing robust eliminations of false positives.Comment: 8 pages, 3 figure

Asteroseismology of old open clusters with Kepler: direct estimate of the integrated RGB mass loss in NGC6791 and NGC6819

Mass loss of red giant branch (RGB) stars is still poorly determined, despite its crucial role in the chemical enrichment of galaxies. Thanks to the recent detection of solar-like oscillations in G-K giants in open clusters with Kepler, we can now directly determine stellar masses for a statistically significant sample of stars in the old open clusters NGC6791 and NGC6819. The aim of this work is to constrain the integrated RGB mass loss by comparing the average mass of stars in the red clump (RC) with that of stars in the low-luminosity portion of the RGB (i.e. stars with L <~ L(RC)). Stellar masses were determined by combining the available seismic parameters numax and Dnu with additional photometric constraints and with independent distance estimates. We measured the masses of 40 stars on the RGB and 19 in the RC of the old metal-rich cluster NGC6791. We find that the difference between the average mass of RGB and RC stars is small, but significant (Delta M=0.09 +- 0.03 (random) +- 0.04 (systematic) Msun). Interestingly, such a small DeltaM does not support scenarios of an extreme mass loss for this metal-rich cluster. If we describe the mass-loss rate with Reimers' prescription, a first comparison with isochrones suggests that the observed DeltaM is compatible with a mass-loss efficiency parameter in the range 0.1 <~ eta <~ 0.3. Less stringent constraints on the RGB mass-loss rate are set by the analysis of the ~ 2 Gyr-old NGC6819, largely due to the lower mass loss expected for this cluster, and to the lack of an independent and accurate distance determination. In the near future, additional constraints from frequencies of individual pulsation modes and spectroscopic effective temperatures, will allow further stringent tests of the Dnu and numax scaling relations, which provide a novel, and potentially very accurate, means of determining stellar radii and masses.Comment: 13 pages, 7 figures, accepted for publication in MNRA

Discovery and Rossiter-McLaughlin Effect of Exoplanet Kepler-8b

We report the discovery and the Rossiter-McLaughlin effect of Kepler-8b, a transiting planet identified by the NASA Kepler Mission. Kepler photometry and Keck-HIRES radial velocities yield the radius and mass of the planet around this F8IV subgiant host star. The planet has a radius RP = 1.419 RJ and a mass, MP = 0.60 MJ, yielding a density of 0.26 g cm^-3, among the lowest density planets known. The orbital period is P = 3.523 days and orbital semima jor axis is 0.0483+0.0006/-0.0012 AU. The star has a large rotational v sin i of 10.5 +/- 0.7 km s^-1 and is relatively faint (V = 13.89 mag), both properties deleterious to precise Doppler measurements. The velocities are indeed noisy, with scatter of 30 m s^-1, but exhibit a period and phase consistent with the planet implied by the photometry. We securely detect the Rossiter-McLaughlin effect, confirming the planet's existence and establishing its orbit as prograde. We measure an inclination between the projected planetary orbital axis and the projected stellar rotation axis of lambda = -26.9 +/- 4.6 deg, indicating a moderate inclination of the planetary orbit. Rossiter-McLaughlin measurements of a large sample of transiting planets from Kepler will provide a statistically robust measure of the true distribution of spin-orbit orientations for hot jupiters in general.Comment: 26 pages, 8 figures, 2 tables; In preparation for submission to the Astrophysical Journa