Kepler Quarters Research Articles

A Search for Temporal Atmospheric Variability of Kepler Hot Jupiters

Abstract We perform a systematic search for atmospheric variability in short-period gas-giant planets (hot Jupiters) observed by the Kepler mission, by looking for temporal variability of their secondary eclipse depths. This is motivated by a recent detection of a decrease in the dayside brightness of KELT-1 b between TESS Sectors 17 and 57, separated by about 3 yr. We fit the Kepler light curves of 53 hot Jupiters and measure their secondary eclipse depths during individual Kepler quarters and four-quarter windows. We detect the secondary eclipses in individual quarters or four-quarter windows for 17 out of the 53 systems. In those 17 systems we do not detect statistically significant astrophysical variation in the secondary eclipse depths. We show that the data is sensitive to the variability seen for KELT-1 b in TESS data. Therefore, the absence of detected secondary eclipse variability in Kepler data suggests that the atmospheric variability in KELT-1 b is not common. In addition, several of the 53 targets we investigated display variability in their transit depths with a period of four quarters (1 yr). This instrumental signal is likely present in the light curves of other transiting planets we did not analyze and other variable stars observed by Kepler. Finally, we find that Kepler-488 b has a secondary eclipse depth that is unphysically large for a planet, and thus is likely a misclassified red dwarf.

Kepler-411: a four-planet system with an active host star

We present a detailed characterization of the Kepler-411 system (KOI 1781). This system was previously known to host two transiting planets: one with a period of 3 days (R = 2.4 R⊕; Kepler-411b) and one with a period of 7.8 days (R = 4.4 R⊕; Kepler-411c), as well as a transiting planetary candidate with a 58-day period (R = 3.3 R⊕; KOI 1781.03) from Kepler photometry. Here, we combine Kepler photometry data and new transit timing variation (TTV) measurements from all the Kepler quarters with previous adaptive-optics imaging results, and dynamical simulations, in order to constrain the properties of the Kepler-411 system. From our analysis, we obtain masses of 25.6 ± 2.6 M⊕ for Kepler-411b and 26.4 ± 5.9 M⊕ for Kepler-411c, and we confirm the planetary nature of KOI 1781.03 with a mass of 15.2 ± 5.1 M⊕, hence the name Kepler-411d. Furthermore, by assuming near-coplanarity of the system (mutual inclination below 30°), we discover a nontransiting planet, Kepler-411e, with a mass of 10.8 ± 1.1 M⊕ on a 31.5-day orbit, which has a strong dynamical interaction with Kepler-411d. With densities of 1.71 ± 0.39 g cm−3 and 2.32 ± 0.83 g cm−3, both Kepler-411c and Kepler-411d belong to the group of planets with a massive core and a significant fraction of volatiles. Although Kepler-411b has a sub-Neptune size, it belongs to the group of rocky planets.

Masses of the Kepler-419 planets from transit timing variations analysis

We perform dynamical fits to the Transit Timing Variations (TTVs) of Kepler-419b. The TTVs from 17 Kepler quarters are obtained from Holczer et al (2016, ApJS 225,9). The dynamical fits are performed using the MultiNest Bayesian inference tool, coupled to an efficient symplectic N-body integrator. We find that the existing TTV data alone are able to uniquely constrain the planetary masses of Kepler-419b and c. Our estimates are in a good agreement with previous mass determinations that combined different techniques and observations, such as TTVs, radial velocity measurements, and the photoeccentric effect. As expected, however, our mass estimates have larger uncertainty. We study the global stability of the system within the parameters uncertainties to discard possible unstable solutions. We conclude that our method applied to the Holczer et al. data can provide reliable determinations of the planetary parameters. In the fore-coming work, we will use Holczer et al. (2016b) data to determine planetary parameters for a large set of Kepler systems.

Spin–Orbit Misalignment and Precession in the Kepler-13Ab Planetary System

Abstract Gravity darkening induced by rapid stellar rotation provides us with a unique opportunity to characterize the spin–orbit misalignment of a planetary system through analysis of its photometric transit. We use the gravity-darkened transit modeling code simuTrans to reproduce the transit light curve of Kepler-13Ab by separately analyzing phase-folded transits for 12 short-cadence Kepler quarters. We verify the temporal change in impact parameter indicative of spin–orbit precession identified by Szabó et al. and Masuda, reporting a rate of change day−1. We further investigate the effect of light dilution on the fitted impact parameter and find that less than 1% of additional light is sufficient to explain the seasonal variation seen in the Kepler quarter data. We then extend our precession analysis to the phase curve data from which we report a rate of change day−1. This value is consistent with that of the transit data at a lower significance and provides the first evidence of spin–orbit precession based solely on the temporal variation of the secondary eclipse.

The Kepler-19 System: A Thick-envelope Super-Earth with Two Neptune-mass Companions Characterized Using Radial Velocities and Transit Timing Variations

Abstract We report a detailed characterization of the Kepler-19 system. This star was previously known to host a transiting planet with a period of 9.29 days, a radius of 2.2 R ⊕, and an upper limit on the mass of 20 M ⊕. The presence of a second, non-transiting planet was inferred from the transit time variations (TTVs) of Kepler-19b over eight quarters of Kepler photometry, although neither the mass nor period could be determined. By combining new TTVs measurements from all the Kepler quarters and 91 high-precision radial velocities obtained with the HARPS-N spectrograph, using dynamical simulations we obtained a mass of 8.4 ± 1.6 M ⊕ for Kepler-19b. From the same data, assuming system coplanarity, we determined an orbital period of 28.7 days and a mass of 13.1 ± 2.7 M ⊕ for Kepler-19c and discovered a Neptune-like planet with a mass of 20.3 ± 3.4 M ⊕ on a 63-day orbit. By comparing dynamical simulations with non-interacting Keplerian orbits, we concluded that neglecting interactions between planets may lead to systematic errors that can hamper the precision in the orbital parameters when the data set spans several years. With a density of 4.32 ± 0.87 g cm−3 (0.78 ± 0.16 ρ ⊕) Kepler-19b belongs to the group of planets with a rocky core and a significant fraction of volatiles, in opposition to low-density planets characterized only by transit time variations and an increasing number of rocky planets with Earth-like density. Kepler-19 joins the small number of systems that reconcile transit timing variation and radial velocity measurements.

Periodicities of the RV Tauri-type pulsating star DF Cygni: A combination ofKeplerdata with ground-based observations

The RV Tauri stars constitute a small group of classical pulsating stars with some dozen known members in the Milky Way. The light variation is caused predominantly by pulsations, but these alone do not explain the full complexity of the light curves. High quality photometry of RV Tau-type stars is very rare. DF Cygni is the only member of this class of stars in the original Kepler field, hence allowing the most accurate photometric investigation of an RV Tauri star to date. The main goal is to analyse the periodicities of the RV Tauri-type star DF Cygni by combining four years of high-quality Kepler photometry with almost half a century of visual data collected by the American Association of Variable Star Observers. Kepler quarters of data have been stitched together to minimize the systematic effects of the space data. The mean levels have been matched with the AAVSO visual data. Both datasets have been submitted to Fourier and wavelet analyses, while the stability of the main pulsations has been studied with the O-C method and the analysis of the time-dependent amplitudes. DF Cygni shows a very rich behaviour on all time-scales. The slow variation has a period of 779.606 d and it has been remarkably coherent during the whole time-span of the combined data. On top of the long-term cycles the pulsations appear with a period of 24.925 d. Both types of light variation significantly fluctuate in time, with a constantly changing interplay of amplitude and phase modulations. The long-period change (i.e. the RVb signature) somewhat resembles the Long Secondary Period (LSP) phenomenon of the pulsating red giants, whereas the short-period pulsations are very similar to those of the Cepheid variables. Comparing the pulsation patterns with the latest models of Type-II Cepheids, we found evidence of strong non-linear effects directly observable in the Kepler light curve.

Photometry of very bright stars with Kepler and K2 smear data

Abstract High-precision time series photometry with the Kepler satellite has been crucial to our understanding both of exoplanets, and via asteroseismology, of stellar physics. After the failure of two reaction wheels, the Kepler satellite has been repurposed as Kepler-2 (K2), observing fields close to the ecliptic plane. As these fields contain many more bright stars than the original Kepler field, K2 provides an unprecedented opportunity to study nearby objects amenable to detailed follow-up with ground-based instruments. Due to bandwidth constraints, only a small fraction of pixels can be downloaded, with the result that most bright stars which saturate the detector are not observed. We show that engineering data acquired for photometric calibration, consisting of collateral ‘smear’ measurements, can be used to reconstruct light curves for bright targets not otherwise observable with Kepler/K2. Here we present some examples from Kepler Quarter 6 and K2 Campaign 3, including the δ Scuti variables HD 178875 and 70 Aqr, and the red giant HR 8500 displaying solar-like oscillations. We compare aperture and smear photometry where possible, and also study targets not previously observed. These encouraging results suggest this new method can be applied to most Kepler and K2 fields.

A Bayesian method for detecting stellar flares

We present a Bayesian-odds-ratio-based algorithm for detecting stellar flares in light curve data. We assume flares are described by a model in which there is a rapid rise with a half-Gaussian profile, followed by an exponential decay. Our signal model also contains a polynomial background model. This is required to fit underlying light curve variations that are expected in the data, which could otherwise partially mimic a flare. We characterise the false alarm probability and efficiency of this method and compare it with a simpler thresholding method based on that used in Walkowicz et al (2011). We find our method has a significant increase in detection efficiency for low signal-to-noise ratio (S/N) flares. For a conservative false alarm probability our method can detect 95% of flares with S/N less than ~20, as compared to S/N of ~25 for the simpler method. As an example we have applied our method to a selection of stars in Kepler Quarter 1 data. The method finds 687 flaring stars with a total of 1873 flares after vetos have been applied. For these flares we have characterised their durations and and signal-to-noise ratios.

Mapping a star with transits: orbit precession effects in the Kepler-13 system

Kepler-13b is a most intriguing exoplanet system due to the rapid precession rate, exhibiting several exotic phenomena. We analysed Kepler short cadence data up to Quarter 14, with a total time-span of 928 d, to reveal changes in transit duration, depth, asymmetry and identify the possible signals of stellar rotation and low-level activity. We investigated long-term variations of transit light curves testing for duration, peak depth and asymmetry. We also performed cluster analysis on Kepler quarters. We computed the autocorrelation function of the out-of-transit light variations. Transit duration, peak depth and asymmetry evolve slowly, due to the slowly drifting transit path through the stellar disc. The detected transit shapes will map the stellar surface on the time-scale of decades. We found a very significant clustering pattern with 3-orbit period. Its source is very probably the rotating stellar surface, in the 5:3 spin–orbit resonance reported in a previous study. The autocorrelation function of the out-of-transit light variations, filtered to 25.4 h and harmonics, shows slow variations and a peak around 300–360 d period, which could be related to the activity cycle of the host star.

Rotation periods of 12 000 main-sequenceKeplerstars: Dependence on stellar spectral type and comparison withvsiniobservations

Aims: We aim to measure the starspot rotation periods of active stars in the Kepler field as a function of spectral type and to extend reliable rotation measurements from F-, G-, and K-type to M-type stars. Methods: Using the Lomb-Scargle periodogram we searched more than 150 000 stellar light curves for periodic brightness variations. We analyzed periods between 1 and 30 days in eight consecutive Kepler quarters, where 30 days is an estimated maximum for the validity of the PDC_MAP data correction pipeline. We selected stable rotation periods, i.e., periods that do not vary from the median by more than one day in at least six of the eight quarters. We averaged the periods for each stellar spectral class according to B - V color and compared the results to archival vsini data, using stellar radii estimates from the Kepler Input Catalog. Results: We report on the stable starspot rotation periods of 12 151 Kepler stars. We find good agreement between starspot velocities and vsini data for all F-, G- and early K-type stars. The 795 M-type stars in our sample have a median rotation period of 15.4 days. We find an excess of M-type stars with periods less than 7.5 days that are potentially fast-rotating and fully convective. Measuring photometric variability in multiple Kepler quarters appears to be a straightforward and reliable way to determine the rotation periods of a large sample of active stars, including late-type stars.

KEPLERPHOTOMETRY OF FOUR RADIO-LOUD ACTIVE GALACTIC NUCLEI IN 2010-2012

We have used Kepler photometry to characterize variability in four radio-loud active galactic nuclei (three quasars and one object tentatively identified as a Seyfert 1.5 galaxy) on timescales from minutes to months, comparable to the light crossing time of the accretion disk around the central supermassive black hole or the base of the relativistic jet. Kepler's almost continuous observations provide much better temporal coverage than is possible from ground-based observations. We report the first such data analyzed for quasars. We have constructed power spectral densities using 8 Kepler quarters of long-cadence (30-minute) data for three AGN, 6 quarters for one AGN and 2 quarters of short-cadence (1-minute) data for all four AGN. On timescales longer than about 0.2-0.6 day, we find red noise with mean power-law slopes ranging from -1.8 to -1.2, consistent with the variability originating in turbulence either behind a shock or within an accretion disk. Each AGN has a range of red noise slopes which vary slightly by month and quarter of observation. No quasi-periodic oscillations of astrophysical origin were detected. We detected several days-long flares when brightness increased by 3% - 7% in two objects. No flares on timescales of minutes to hours were detected. Our observations imply that the duty cycle for enhanced activity in these radio-loud AGN is small. These well-sampled AGN light curves provide an impetus to develop more detailed models of turbulence in jets and instabilities in accretion disks.

PERIOD ERROR ESTIMATION FOR THE KEPLER ECLIPSING BINARY CATALOG

The Kepler Eclipsing Binary Catalog (KEBC)describes 2165 eclipsing binaries identified in the 115 deg^2 Kepler Field based on observations from Kepler quarters Q0, Q1, and Q2. The periods in the KEBC are given in units of days out to six decimal places but no period errors are provided. We present the PEC (Period Error Calculator) algorithm which can be used to estimate the period errors of strictly periodic variables observed by the Kepler Mission. The PEC algorithm is based on propagation of error theory and assumes that observation of every light curve peak/minimum in a long time-series observation can be unambiguously identified. The PEC algorithm can be efficiently programmed using just a few lines of C computer language code. The PEC algorithm was used to develop a simple model which provides period error estimates for eclipsing binaries in the KEBC with periods less than 62.5 days. KEBC systems with periods >=62.5 days have KEBC period errors of about 0.0144 days. Periods and period errors of 7 eclipsing binary systems in the KEBC were measured using the NASA Exoplanet Archive Periodogram Service and compared to period errors estimated using the PEC algorithm.

THE NEPTUNE-SIZED CIRCUMBINARY PLANET KEPLER-38b

We discuss the discovery and characterization of the circumbinary planet Kepler-38b. The stellar binary is single-lined, with a period of 18.8 days, and consists of a moderately evolved main-sequence star (M_A = 0.949 +/- 0.059 solar masses and R_A = 1.757 +/- 0.034 solar radii) paired with a low-mass star (M_B = 0.249 +/- 0.010 solar masses and R_B = 0.2724 +/- 0.0053 solar radii) in a mildly eccentric (e=0.103) orbit. A total of eight transits due to a circumbinary planet crossing the primary star were identified in the Kepler light curve (using Kepler Quarters 1 through 11), from which a planetary period of 105.595 +/- 0.053 days can be established. A photometric dynamical model fit to the radial velocity curve and Kepler light curve yields a planetary radius of 4.35 +/- 0.11 Earth radii, or 1.12 +/- 0.03 Neptune radii. Since the planet is not sufficiently massive to observably alter the orbit of the binary from Keplerian motion, we can only place an upper limit on the mass of the planet of 122 Earth masses (7.11 Neptune masses or 0.384 Jupiter masses) at 95% confidence. This upper limit should decrease as more Kepler data become available.

PLANET HUNTERS: ASSESSING THEKEPLERINVENTORY OF SHORT-PERIOD PLANETS

We present the results from a search of data from the first 33.5 days of the Kepler science mission (Quarter 1) for exoplanet transits by the Planet Hunters citizen science project. Planet Hunters enlists members of the general public to visually identify transits in the publicly released Kepler light curves via the World Wide Web. Over 24,000 volunteers reviewed the Kepler Quarter 1 data set. We examine the abundance of \geq 2 R\oplus planets on short period (< 15 days) orbits based on Planet Hunters detections. We present these results along with an analysis of the detection efficiency of human classifiers to identify planetary transits including a comparison to the Kepler inventory of planet candidates. Although performance drops rapidly for smaller radii, \geq 4 R\oplus Planet Hunters \geq 85% efficient at identifying transit signals for planets with periods less than 15 days for the Kepler sample of target stars. Our high efficiency rate for simulated transits along with recovery of the majority of Kepler \geq 4 R\oplus planets suggest suggests the Kepler inventory of \geq 4 R\oplus short period planets is nearly complete.

THE EVIL-MC MODEL FOR ELLIPSOIDAL VARIATIONS OF PLANET-HOSTING STARS AND APPLICATIONS TO THE HAT-P-7 SYSTEM

We present a new model for Ellipsoidal Variations Induced by a Low-Mass Companion, the EVIL-MC model. We employ several approximations appropriate for planetary systems to substantially increase the computational efficiency of our model relative to more general ellipsoidal variation models and improve upon the accuracy of simpler models. This new approach gives us a unique ability to rapidly and accurately determine planetary system parameters. We use the EVIL-MC model to analyze Kepler Quarter 0-2 (Q0-2) observations of the HAT-P-7 system, an F-type star orbited by a nearly Jupiter-mass companion. Our analysis corroborates previous estimates of the planet-star mass ratio q = (1.10 +/- 0.06) x 10^(-3), and we have revised the planet's dayside brightness temperature to 2680 +10/-20 K. We also find a large difference between the day- and nightside planetary flux, with little nightside emission. Preliminary dynamical+radiative modeling of the atmosphere indicates this result is qualitatively consistent with high altitude absorption of stellar heating. Similar analyses of Kepler and CoRoT photometry of other planets using EVIL-MC will play a key role in providing constraints on the properties of many extrasolar systems, especially given the limited resources for follow-up and characterization of these systems. However, as we highlight, there are important degeneracies between the contributions from ellipsoidal variations and planetary emission and reflection. Consequently, for many of the hottest and brightest Kepler and CoRoT planets, accurate estimates of the planetary emission and reflection, diagnostic of atmospheric heat budgets, will require accurate modeling of the photometric contribution from the stellar ellipsoidal variation.

DETECTION OF KOI-13.01 USING THE PHOTOMETRIC ORBIT

We use the KOI-13 transiting star-planet system as a test case for the recently developed BEER algorithm (Faigler & Mazeh 2011), aimed at identifying non-transiting low-mass companions by detecting the photometric variability induced by the companion along its orbit. Such photometric variability is generated by three mechanisms, including the beaming effect, tidal ellipsoidal distortion, and reflection/heating. We use data from three Kepler quarters, from the first year of the mission, while ignoring measurements within the transit and occultation, and show that the planet's ephemeris is clearly detected. We fit for the amplitude of each of the three effects and use the beaming effect amplitude to estimate the planet's minimum mass, which results in M_p sin i = 9.2 +/- 1.1 M_J (assuming the host star parameters derived by Szabo et al. 2011). Our results show that non-transiting star-planet systems similar to KOI-13.01 can be detected in Kepler data, including a measurement of the orbital ephemeris and the planet's minimum mass. Moreover, we derive a realistic estimate of the amplitudes uncertainties, and use it to show that data obtained during the entire lifetime of the Kepler mission, of 3.5 years, will allow detecting non-transiting close-in low-mass companions orbiting bright stars, down to the few Jupiter mass level. Data from the Kepler Extended Mission, if funded by NASA, will further improve the detection capabilities.

WHITE-LIGHT FLARES ON COOL STARS IN THEKEPLERQUARTER 1 DATA

We present the results of a search for white light flares on the ~23,000 cool dwarfs in the Kepler Quarter 1 long cadence data. We have identified 373 flaring stars, some of which flare multiple times during the observation period. We calculate relative flare energies, flare rates and durations, and compare these with the quiescent photometric variability of our sample. We find that M dwarfs tend to flare more frequently but for shorter durations than K dwarfs, and that they emit more energy relative to their quiescent luminosity in a given flare than K dwarfs. Stars that are more photometrically variable in quiescence tend to emit relatively more energy during flares, but variability is only weakly correlated with flare frequency. We estimate distances for our sample of flare stars and find that the flaring fraction agrees well with other observations of flare statistics for stars within 300 pc above the Galactic Plane. These observations provide a more rounded view of stellar flares by sampling stars that have not been pre-selected by their activity, and are informative for understanding the influence of these flares on planetary habitability.