Stellar variability is driven by a multitude of internal physical processes that depend on fundamental stellar properties. These properties are our bridge to reconciling stellar observations with stellar physics, and for understanding the distribution of stellar populations within the context of galaxy formation. Numerous ongoing and upcoming missions are charting brightness fluctuations of stars over time, which encode information about physical processes such as rotation period, evolutionary state (such as effective temperature and surface gravity), and mass (via asteroseismic parameters). Here, we explore how well we can predict these stellar properties, across different evolutionary states, using only photometric time series data. To do this, we implement a convolutional neural network, and with data-driven modeling we predict stellar properties from light curves of various baselines and cadences. Based on a single quarter of \textit{Kepler} data, we recover stellar properties, including surface gravity for red giant stars (with an uncertainty of $\lesssim$ 0.06 dex), and rotation period for main sequence stars (with an uncertainty of $\lesssim$ 5.2 days, and unbiased from $\approx$5 to 40 days). Shortening the \textit{Kepler} data to a 27-day TESS-like baseline, we recover stellar properties with a small decrease in precision, $\sim$0.07 dex for log $g$ and $\sim$5.5 days for $P_{\rm rot}$, unbiased from $\approx$5 to 35 days. Our flexible data-driven approach leverages the full information content of the data, requires minimal feature engineering, and can be generalized to other surveys and datasets. This has the potential to provide stellar property estimates for many millions of stars in current and future surveys.

The rotational evolution of cool dwarfs is poorly constrained after around 1-2 Gyr due to a lack of precise ages and rotation periods for old main-sequence stars. In this work we use velocity dispersion as an age proxy to reveal the temperature-dependent rotational evolution of low-mass Kepler dwarfs, and demonstrate that kinematic ages could be a useful tool for calibrating gyrochronology in the future. We find that a linear gyrochronology model, calibrated to fit the period-Teff relationship of the Praesepe cluster, does not apply to stars older than around 1 Gyr. Although late-K dwarfs spin more slowly than early-K dwarfs when they are young, at old ages we find that late-K dwarfs rotate at the same rate or faster than early-K dwarfs of the same age. This result agrees qualitatively with semi-empirical models that vary the rate of surface-to-core angular momentum transport as a function of time and mass. It also aligns with recent observations of stars in the NGC 6811 cluster, which indicate that the surface rotation rates of K dwarfs go through an epoch of inhibited evolution. We find that the oldest Kepler stars with measured rotation periods are late-K and early-M dwarfs, indicating that these stars maintain spotted surfaces and stay magnetically active longer than more massive stars. Finally, based on their kinematics, we confirm that many rapidly rotating GKM dwarfs are likely to be synchronized binaries.

The claim of an exomoon candidate in the Kepler-1625b system has generated substantial discussion regarding possible alternative explanations for the purported signal. In this work we examine in detail these possibilities. First, the effect of more flexible trend models is explored and we show that sufficiently flexible models are capable of attenuating the signal, although this is an expected byproduct of invoking such models. We also explore trend models using X and Y centroid positions and show that there is no data-driven impetus to adopt such models over temporal ones. We quantify the probability that the 500 ppm moon-like dip could be caused by a Neptune-sized transiting planet to be < 0.75%. We show that neither autocorrelation, Gaussian processes nor a Lomb-Scargle periodogram are able to recover a stellar rotation period, demonstrating that K1625 is a quiet star with periodic behavior < 200 ppm. Through injection and recovery tests, we find that the star does not exhibit a tendency to introduce false-positive dip-like features above that of pure Gaussian noise. Finally, we address a recent re-analysis by Kreidberg et al (2019) and show that the difference in conclusions is not from differing systematics models but rather the reduction itself. We show that their reduction exhibits i) slightly higher intra-orbit and post-fit residual scatter, ii) $\simeq$ 900 ppm larger flux offset at the visit change, iii) $\simeq$ 2 times larger Y-centroid variations, and iv) $\simeq$ 3.5 times stronger flux-centroid correlation coefficient than the original analysis. These points could be explained by larger systematics in their reduction, potentially impacting their conclusions.

We present a new age-dating technique that combines gyrochronology with isochrone fitting to infer ages for FGKM main-sequence and subgiant field stars. Gyrochronology and isochrone fitting are each capable of providing relatively precise ages for field stars in certain areas of the Hertzsprung-Russell diagram: gyrochronology works optimally for cool main-sequence stars, and isochrone fitting can provide precise ages for stars near the main-sequence turnoff. Combined, these two age-dating techniques can provide precise and accurate ages for a broader range of stellar masses and evolutionary stages than either method used in isolation. We demonstrate that the position of a star on the Hertzsprung- Russell or color-magnitude diagram can be combined with its rotation period to infer a precise age via both isochrone fitting and gyrochronology simultaneously. We show that incorporating rotation periods with 5% uncertainties into stellar evolution models improves age precision for FGK stars on the main sequence, and can, on average, provide age estimates up to three times more precise than isochrone fitting alone. In addition, we provide a new gyrochronology relation, calibrated to the Praesepe cluster and the Sun, that includes a variance model to capture the rotational behavior of stars whose rotation periods do not lengthen with the square-root of time, and parts of the Hertzsprung-Russell diagram where gyrochronology has not been calibrated. This publication is accompanied by an open source Python package, stardate, for inferring the ages of main-sequence and subgiant FGKM stars from rotation periods, spectroscopic parameters and/or apparent magnitudes and parallaxes.

We present a sample of 74,216 M and L dwarfs constructed from two existing catalogs of cool dwarfs spectroscopically identified in the Sloan Digital Sky Survey (SDSS). We cross-matched the SDSS catalog with Gaia DR2 to obtain parallaxes and proper motions and modified the quality cuts suggested by the Gaia Collaboration to make them suitable for late-M and L dwarfs. We also provide relations between Gaia colors and absolute magnitudes with spectral type and conclude that (G-RP) has the tightest relation to spectral type for M and L dwarfs. In addition, we study magnetic activity as a function of position on the color-magnitude diagram, finding that Halpha magnetically active stars have, on average, redder colors and/or brighter magnitudes than inactive stars. This effect cannot be explained by youth alone and might indicate that active stars are magnetically inflated, binaries and/or high metallicity. Moreover, we find that vertical velocity and vertical action dispersion are correlated with Halpha emission, confirming that these two parameters are age indicators. We also find that stars below the main sequence have high tangential velocity which is consistent with a low metallicity and old population of stars that belong to the halo or thick disk.

We use archival daily spot coverage measurements from Howard et al. (1984) to study the rotational modulation of the Sun as though it were a distant star. A quasi-periodic Gaussian process measures the solar rotation period $P_\mathrm{rot} = 26.3 \pm 0.1$ days, and activity cycle period $P_\mathrm{cyc} = 10.7 \pm 0.3$ years. We attempt to search for evidence of differential rotation in variations of the apparent rotation period throughout the activity cycle and do not detect a clear signal of differential rotation, consistent with the null results of the hare-and-hounds exercise of Aigrain et al. (2015). The full reconstructed solar light curve is available online.

It has been demonstrated that the time variability of a star's brightness at different frequencies can be used to infer its surface gravity, radius, mass, and age. With large samples of light curves now available from Kepler and K2, and upcoming surveys like TESS, we wish to quantify the overall information content of this data and identify where the information resides. As a first look into this question we ask which stellar parameters we can predict from the brightness variations in red-giant stars data and to what precision, using a data-driven model. We demonstrate that the long-cadence (30-minute) Kepler light curves for 2000 red-giant stars can be used to predict their $T_{\rm eff}$ and $\log g$. Our inference makes use of a data-driven model of a part of the autocorrelation function (ACF) of the light curve, where we posit a polynomial relationship between stellar parameters and the ACF pixel values. We find that this model, trained using 1000 stars, can be used to recover the temperature $T_{\rm eff}$ to $<$100 K, the surface gravity to $<$ 0.1 dex, and the asteroseismic power-spectrum parameters $\rm \nu_{max}$ and $\rm \Delta{\nu}$ to $<11$ $\mu$Hz and $<0.9$ $\mu$Hz ($\lesssim$ 15\%). We recover $T_{\rm eff}$ from range of time-lags 0.045 $<$ $T_{\rm lag}$ $<$ 370 days and the $\log g$, $\rm \nu_{max}$ and $\rm \Delta{\nu}$ from the range 0.045 $<$ $T_{\rm lag}$ $<$ 35 days. We do not discover any information about stellar metallicity. The information content of the data about each parameter is empirically quantified using this method, enabling comparisons to theoretical expectations about convective granulation.

The Large Synoptic Survey Telescope is designed to provide an unprecedented optical imaging dataset that will support investigations of our Solar System, Galaxy and Universe, across half the sky and over ten years of repeated observation. However, exactly how the LSST observations will be taken (the observing strategy or "cadence") is not yet finalized. In this dynamically-evolving community white paper, we explore how the detailed performance of the anticipated science investigations is expected to depend on small changes to the LSST observing strategy. Using realistic simulations of the LSST schedule and observation properties, we design and compute diagnostic metrics and Figures of Merit that provide quantitative evaluations of different observing strategies, analyzing their impact on a wide range of proposed science projects. This is work in progress: we are using this white paper to communicate to each other the relative merits of the observing strategy choices that could be made, in an effort to maximize the scientific value of the survey. The investigation of some science cases leads to suggestions for new strategies that could be simulated and potentially adopted. Notably, we find motivation for exploring departures from a spatially uniform annual tiling of the sky: focusing instead on different parts of the survey area in different years in a "rolling cadence" is likely to have significant benefits for a number of time domain and moving object astronomy projects. The communal assembly of a suite of quantified and homogeneously coded metrics is the vital first step towards an automated, systematic, science-based assessment of any given cadence simulation, that will enable the scheduling of the LSST to be as well-informed as possible.

The growing field of large-scale time domain astronomy requires methods for probabilistic data analysis that are computationally tractable, even with large datasets. Gaussian Processes are a popular class of models used for this purpose but, since the computational cost scales, in general, as the cube of the number of data points, their application has been limited to small datasets. In this paper, we present a novel method for Gaussian Process modeling in one-dimension where the computational requirements scale linearly with the size of the dataset. We demonstrate the method by applying it to simulated and real astronomical time series datasets. These demonstrations are examples of probabilistic inference of stellar rotation periods, asteroseismic oscillation spectra, and transiting planet parameters. The method exploits structure in the problem when the covariance function is expressed as a mixture of complex exponentials, without requiring evenly spaced observations or uniform noise. This form of covariance arises naturally when the process is a mixture of stochastically-driven damped harmonic oscillators -- providing a physical motivation for and interpretation of this choice -- but we also demonstrate that it can be a useful effective model in some other cases. We present a mathematical description of the method and compare it to existing scalable Gaussian Process methods. The method is fast and interpretable, with a range of potential applications within astronomical data analysis and beyond. We provide well-tested and documented open-source implementations of this method in C++, Python, and Julia.

Variability in the light curves of spotted, rotating stars is often non-sinusoidal and quasi-periodic --- spots move on the stellar surface and have finite lifetimes, causing stellar flux variations to slowly shift in phase. A strictly periodic sinusoid therefore cannot accurately model a rotationally modulated stellar light curve. Physical models of stellar surfaces have many drawbacks preventing effective inference, such as highly degenerate or high-dimensional parameter spaces. In this work, we test an appropriate effective model: a Gaussian Process with a quasi-periodic covariance kernel function. This highly flexible model allows sampling of the posterior probability density function of the periodic parameter, marginalising over the other kernel hyperparameters using a Markov Chain Monte Carlo approach. To test the effectiveness of this method, we infer rotation periods from 333 simulated stellar light curves, demonstrating that the Gaussian process method produces periods that are more accurate than both a sine-fitting periodogram and an autocorrelation function method. We also demonstrate that it works well on real data, by inferring rotation periods for 275 Kepler stars with previously measured periods. We provide a table of rotation periods for these 1132 Kepler objects of interest and their posterior probability density function samples. Because this method delivers posterior probability density functions, it will enable hierarchical studies involving stellar rotation, particularly those involving population modelling, such as inferring stellar ages, obliquities in exoplanet systems, or characterising star-planet interactions. The code used to implement this method is available online.

The Large Synoptic Survey Telescope (LSST) will be a discovery machine for the astronomy and physics communities, revealing astrophysical phenomena from the Solar System to the outer reaches of the observable Universe. While many discoveries will be made using LSST data alone, taking full scientific advantage of LSST will require ground-based optical-infrared (OIR) supporting capabilities, e.g., observing time on telescopes, instrumentation, computing resources, and other infrastructure. This community-based study identifies, from a science-driven perspective, capabilities that are needed to maximize LSST science. Expanding on the initial steps taken in the 2015 OIR System Report, the study takes a detailed, quantitative look at the capabilities needed to accomplish six representative LSST-enabled science programs that connect closely with scientific priorities from the 2010 decadal surveys. The study prioritizes the resources needed to accomplish the science programs and highlights ways that existing, planned, and future resources could be positioned to accomplish the science goals.

The relations between observable stellar parameters are usually assumed to be deterministic. That is, given an infinitely precise measurement of independent variable, `$x$', and some model, the value of dependent variable, `$y$' can be known exactly. In practise this assumption is rarely valid and intrinsic stochasticity means that two stars with exactly the same `$x$', will have slightly different `$y$'s. The relation between short-timescale brightness fluctuations (flicker) of stars and both surface gravity and stellar density are two such stochastic relations that have, until now, been treated as deterministic ones. We recalibrate these relations in a probabilistic framework, using Hierarchical Bayesian Modelling (HBM) to constrain the intrinsic scatter in the relations. We find evidence for additional scatter in the relationships, that cannot be accounted for by the observational uncertainties alone. The scatter in surface gravity and stellar density does not depend on flicker, suggesting that using flicker as a proxy for $\log g$ and $\rho_\star$ is equally valid for dwarf and giant stars, despite the fact that the observational uncertainties tend to be larger for dwarfs.

This is the stars chapter of the Kavli workshop report, which resulted from the community-based study of needed US OIR capabilities in the LSST era. The full report, which will include this chapter, is anticipated to be available in Fall 2016. See NOAO website (http://www.noao.edu/meetings/lsst-oir-study/) for more details.

White dwarfs are the end state of most stars, including the Sun, after they exhaust their nuclear fuel. Between 1/4 and 1/2 of white dwarfs have elements heavier than helium in their atmospheres, even though these elements should rapidly settle into the stellar interiors unless they are occasionally replenished. The abundance ratios of heavy elements in white dwarf atmospheres are similar to rocky bodies in the Solar system. This and the existence of warm dusty debris disks around about 4% of white dwarfs suggest that rocky debris from white dwarf progenitors' planetary systems occasionally pollute the stars' atmospheres. The total accreted mass can be comparable to that of large asteroids in the solar system. However, the process of disrupting planetary material has not yet been observed. Here, we report observations of a white dwarf being transited by at least one and likely multiple disintegrating planetesimals with periods ranging from 4.5 hours to 4.9 hours. The strongest transit signals occur every 4.5 hours and exhibit varying depths up to 40% and asymmetric profiles, indicative of a small object with a cometary tail of dusty effluent material. The star hosts a dusty debris disk and the star's spectrum shows prominent lines from heavy elements like magnesium, aluminium, silicon, calcium, iron, and nickel. This system provides evidence that heavy element pollution of white dwarfs can originate from disrupted rocky bodies such as asteroids and minor planets.

From pulsating stars to transiting exoplanets, the search for periodic signals in K2 data, Kepler's 2-wheeled extension, is relevant to a long list of scientific goals. Systematics affecting K2 light curves due to the decreased spacecraft pointing precision inhibit the easy extraction of periodic signals from the data. We here develop a method for producing periodograms of K2 light curves that are insensitive to pointing-induced systematics; the Systematics-Insensitive Periodogram (SIP). Traditional sine-fitting periodograms use a generative model to find the frequency of a sinusoid that best describes the data. We extend this principle by including systematic trends, based on a set of 'Eigen light curves', following Foreman-Mackey et al. (2015), in our generative model as well as a sum of sine and cosine functions over a grid of frequencies. Using this method we are able to produce periodograms with vastly reduced systematic features. The quality of the resulting periodograms are such that we can recover acoustic oscillations in giant stars and measure stellar rotation periods without the need for any detrending. The algorithm is also applicable to the detection of other periodic phenomena such as variable stars, eclipsing binaries and short-period exoplanet candidates. The SIP code is available at https://github.com/RuthAngus/SIPK2.

Among the available methods for dating stars, gyrochronology is a powerful one because it requires knowledge of only the star's mass and rotation period. Gyrochronology relations have previously been calibrated using young clusters, with the Sun providing the only age dependence, and are therefore poorly calibrated at late ages. We used rotation period measurements of 310 Kepler stars with asteroseismic ages, 50 stars from the Hyades and Coma Berenices clusters and 6 field stars (including the Sun) with precise age measurements to calibrate the gyrochronology relation, whilst fully accounting for measurement uncertainties in all observable quantities. We calibrated a relation of the form $P=A^n\times(B-V-c)^b$, where $P$ is rotation period in days, $A$ is age in Myr, $B$ and $V$ are magnitudes and $a$, $b$ and $n$ are the free parameters of our model. We found $a = 0.40^{+0.3}_{-0.05}$, $b = 0.31^{+0.05}_{-0.02}$ and $n = 0.55^{+0.02}_{-0.09}$. Markov Chain Monte Carlo methods were used to explore the posterior probability distribution functions of the gyrochronology parameters and we carefully checked the effects of leaving out parts of our sample, leading us to find that no single relation beween rotation period, colour and age can adequately describe all the subsets of our data. The Kepler asteroseismic stars, cluster stars and local field stars cannot all be described by the same gyrochronology relation. The Kepler asteroseismic stars may be subject to observational biases, however the clusters show unexpected deviations from the predicted behaviour, providing concerns for the overall reliability of gyrochronology as a dating method.

We report the first planet discovery from the two-wheeled Kepler (K2) mission: HIP 116454 b. The host star HIP 116454 is a bright (V = 10.1, K = 8.0) K1-dwarf with high proper motion, and a parallax-based distance of 55.2 +/- 5.4 pc. Based on high-resolution optical spectroscopy, we find that the host star is metal-poor with [Fe/H] = -.16 +/- .18, and has a radius R = 0.716 +/- .0024 R_sun and mass M = .775 +/- .027 Msun. The star was observed by the Kepler spacecraft during its Two-Wheeled Concept Engineering Test in February 2014. During the 9 days of observations, K2 observed a single transit event. Using a new K2 photometric analysis technique we are able to correct small telescope drifts and recover the observed transit at high confidence, corresponding to a planetary radius of Rp = 2.53 +/- 0.18 Rearth. Radial velocity observations with the HARPS-N spectrograph reveal a 11.82 +/- 1.33 Mearth planet in a 9.1 day orbit, consistent with the transit depth, duration, and ephemeris. Follow-up photometric measurements from the MOST satellite confirm the transit observed in the K2 photometry and provide a refined ephemeris, making HIP 116454 b amenable for future follow-up observations of this latest addition to the growing population of transiting super-Earths around nearby, bright stars.

We report the discovery of a massive and dense transiting planet CoRoT-27b on a 3.58 day orbit around a 4.2 Gyr-old G2 star. The planet candidate was identified from the CoRoT photometry, and was confirmed as a planet with ground-based spectroscopy. The confirmation of the planet candidate is based on radial velocity observations combined with imaging to rule out blends. The characterisation of the planet and its host star is carried out using a Bayesian approach where all the data (CoRoT photometry, radial velocities, and spectroscopic characterisation of the star) are used jointly. The Bayesian analysis includes a study whether the assumption of white normally distributed noise holds for the CoRoT photometry, and whether the use of a non-normal noise distribution offers advantages in parameter estimation and model selection. CoRoT-27b has a mass of $10.39 \pm 0.55$ $\mathrm{M}_{\rm Jup}$, a radius of $1.01 \pm 0.04$ $\mathrm{R}_{\rm Jup}$, a mean density of $12.6_{-1.67}^{+1.92}$ $\mathrm{g\,cm^{-3}}$, and an effective temperature of $1500 \pm 130$ K. The planet orbits around its host star, a 4.2 Gyr-old G2-star with a mass $M_{\star}=1.06$ $M_{\odot}$, and a radius $R_{\star}=1.05$ $R_{\odot}$, on a $0.048 \pm 0.007$ AU orbit every 3.58 days. The radial velocity observations allow us to exclude highly eccentric orbits, namely, $e<0.065$ with a 99% confidence. Given its high mass and density, theoretical modelling of CoRoT-27b is demanding. We identify two solutions with heavy element mass fractions of $0.11\pm0.08$ $\mathrm{M_{\oplus}}$ and $0.07\pm0.06$ $\mathrm{M_{\oplus}}$, but even solutions void of heavy elements cannot be excluded. We carry out a secondary eclipse search from the CoRoT photometry using a method based on Bayesian model selection, but conclude that the noise level is too high to detect eclipses shallower than 9% of the transit depth.

We outline a proposal to use the Kepler spacecraft in two-wheel mode to monitor a handful of young associations and open clusters, for a few weeks each. Judging from the experience of similar projects using ground-based telescopes and the CoRoT spacecraft, this program would transform our understanding of early stellar evolution through the study of pulsations, rotation, activity, the detection and characterisation of eclipsing binaries, and the possible detection of transiting exoplanets. Importantly, Kepler's wide field-of-view would enable key spatially extended, nearby regions to be monitored in their entirety for the first time, and the proposed observations would exploit unique synergies with the GAIA ESO spectroscopic survey and, in the longer term, the GAIA mission itself. We also outline possible strategies for optimising the photometric performance of Kepler in two-wheel mode by modelling pixel sensitivity variations and other systematics.

In today's mailing, Hogg et al. propose image modeling techniques to maintain 10-ppm-level precision photometry in Kepler data with only two working reaction wheels. While these results are relevant to many scientific goals for the repurposed mission, all modeling efforts so far have used a toy model of the Kepler telescope. Because the two-wheel performance of Kepler remains to be determined, we advocate for the consideration of an alternate strategy for a >1 year program that maximizes the science return from the "low-torque" fields across the ecliptic plane. Assuming we can reach the precision of the original Kepler mission, we expect to detect 800 new planet candidates in the first year of such a mission. Our proposed strategy has benefits for transit timing variation and transit duration variation studies, especially when considered in concert with the future TESS mission. We also expect to help address the first key science goal of Kepler: the frequency of planets in the habitable zone as a function of spectral type.

Kepler's immense photometric precision to date was maintained through satellite stability and precise pointing. In this white paper, we argue that image modeling--fitting the Kepler-downlinked raw pixel data--can vastly improve the precision of Kepler in pointing-degraded two-wheel mode. We argue that a non-trivial modeling effort may permit continuance of photometry at 10-ppm-level precision. We demonstrate some baby steps towards precise models in both data-driven (flexible) and physics-driven (interpretably parameterized) modes. We demonstrate that the expected drift or jitter in positions in the two-weel era will help with constraining calibration parameters. In particular, we show that we can infer the device flat-field at higher than pixel resolution; that is, we can infer pixel-to-pixel variations in intra-pixel sensitivity. These results are relevant to almost any scientific goal for the repurposed mission; image modeling ought to be a part of any two-wheel repurpose for the satellite. We make other recommendations for Kepler operations, but fundamentally advocate that the project stick with its core mission of finding and characterizing Earth analogs. [abridged]

Eight years of optical photometry from OGLE-III are presented of the optical counterpart to the High Mass X-ray Binary system, SMC X-1. The optical data provide the best view to date of the modulation of the light from the system at the binary period of 3.9d. In addition, it is shown for the first time that the light is also modulated at the superorbital period - a period associated with the precessing of a warped accretion disk around the neutron star partner. The implications for the sources of the various components of the optical light in this system are discussed.