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 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.

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.